No Arabic abstract
Many biological electron transfer (ET) reactions are mediated by metal centres in proteins. NADH:ubiquinone oxidoreductase (complex I) contains an intramolecular chain of seven iron-sulphur (FeS) clusters, one of the longest chains of metal centres in biology and a test case for physical models of intramolecular ET. In biology, intramolecular ET is commonly described as a diffusive hopping process, according to the semi-classical theories of Marcus and Hopfield. However, recent studies have raised the possibility that non-trivial quantum mechanical effects play a functioning role in certain biomolecular processes. Here, we extend the semi-classical model for biological ET to incorporate both semi-classical and coherent quantum phenomena using a quantum master equation based on the Holstein Hamiltonian. We test our model on the structurally-defined chain of FeS clusters in complex I. By exploring a wide range of realistic parameters we find that, when the energy profile for ET along the chain is relatively flat, just a small coherent contribution can provide a robust and significant increase in ET rate (above the semi-classical diffusive-hopping rate), even at physiologically-relevant temperatures. Conversely, when the on-site energies vary significantly along the chain the coherent contribution is negligible. For complex I, a crucial respiratory enzyme that is linked to many neuromuscular and degenerative diseases, our results suggest a new contribution towards ensuring that intramolecular ET does not limit the rate of catalysis. For the emerging field of quantum biology, our model is intended as a basis for elucidating the general role of coherent ET in biological ET reactions.
A major challenge in molecular simulations is to describe denaturant-dependent folding of proteins order to make direct comparisons with {it in vitro} experiments. We use the molecular transfer model, which is currently the only method that accomplishes this goal albeit phenomenologically, to quantitatively describe urea-dependent folding of PDZ domain, which plays a significant role in molecular recognition and signaling. Experiments show that urea-dependent unfolding rates of the PDZ2 domain exhibit a downward curvature at high urea concentrations, which has been interpreted by invoking the presence of a sparsely populated high energy intermediate. Simulations using the MTM and a coarse-grained model of PDZ2 are used to show that the intermediate, which has some native-like character, is present in equilibrium both in the presence and absence of urea. The free energy profiles show that there are two barriers separating the folded and unfolded states. Structures of the transition state ensembles, ($TSE1$ separating the unfolded and $I_{EQ}$ and $TSE2$ separating $I_{EQ}$ and the native state), determined using the $P_{fold}$ method, show that $TSE1$ is expanded; $TSE2$ and native-like. Folding trajectories reveal that PDZ2 folds by parallel routes. In one pathway folding occurs exclusively through $I_1$, which resembles $I_{EQ}$. In a fraction of trajectories, constituting the second pathway, folding occurs through a combination of $I_{1}$ and a kinetic intermediate. The radius of gyration ($R_g^{U}$) of the unfolded state is more compact (by $sim$ 9%) under native conditions. Decrease in $R_g^{U}$ occurs on the time scale on the order of utmost $sim$ 20 $mu s$. The modest decrease in $R_g^{U}$ and the rapid collapse suggest that high spatial and temporal resolution are needed to detect compaction in finite-sized proteins.
The conformational change of biological macromolecule is investigated from the point of quantum transition. A quantum theory on protein folding is proposed. Compared with other dynamical variables such as mobile electrons, chemical bonds and stretching-bending vibrations the molecular torsion has the lowest energy and can be looked as the slow variable of the system. Simultaneously, from the multi-minima property of torsion potential the local conformational states are well defined. Following the idea that the slow variables slave the fast ones and using the nonadiabaticity operator method we deduce the Hamiltonian describing conformational change. It is proved that the influence of fast variables on the macromolecule can fully be taken into account through a phase transformation of slow variable wave function. Starting from the conformation- transition Hamiltonian the nonradiative matrix element is calculated in two important cases: A, only electrons are fast variables and the electronic state does not change in the transition process; B, fast variables are not limited to electrons but the perturbation approximation can be used. Then, the general formulas for protein folding rate are deduced. The analytical form of the formula is utilized to study the temperature dependence of protein folding rate and the curious non-Arrhenius temperature relation is interpreted. The decoherence time of quantum torsion state is estimated and the quantum coherence degree of torsional angles in the protein folding is studied by using temperature dependence data. The proposed folding rate formula gives a unifying approach for the study of a large class problems of biological conformational change.
We analyze the problem of the helix-coil transition in explicit solvents analytically by using spin-based models incorporating two different mechanisms of solvent action: explicit solvent action through the formation of solvent-polymer hydrogen bonds that can compete with the intrinsic intra-polymer hydrogen bonded configurations (competing interactions) and implicit solvent action, where the solvent-polymer interactions tune biopolymer configurations by changing the activity of the solvent (non-competing interactions). The overall spin Hamiltonian is comprised of three terms: the background emph{in vacuo} Hamiltonian of the Generalized Model of Polypeptide Chain type and two additive terms that account for the two above mechanisms of solvent action. We show that on this level the solvent degrees of freedom can be {sl explicitly} and {sl exactly} traced over, the ensuing effective partition function combining all the solvent effects in a unified framework. In this way we are able to address helix-coil transitions for polypeptides, proteins, and DNA, with different buffers and different external constraints. Our spin-based effective Hamiltonian is applicable for treatment of such diverse phenomena as cold denaturation, effects of osmotic pressure on the cold and warm denaturation, complicated temperature dependence of the hydrophobic effect as well as providing a conceptual base for understanding the behavior of Intrinsically Disordered Proteins and their analogues.
Solutions of manually purified gastric mucins have been shown to be promising lubricants for biomedical purposes, where they can efficiently reduce friction and wear. However, so far, such mucin solutions have been mostly tested in specific settings, and variations in the composition of the lubricating fluid have not been systematically explored. We here fill this gap and determine the viscosity, adsorption behavior, and lubricity of porcine gastric mucin solutions on hydrophobic surfaces at different pH levels, mucin and salt concentrations and in the presence of other proteins. We demonstrate that mucin solutions provide excellent lubricity even at very low concentrations of 0.01 % (w/v), over a broad range of pH levels and even at elevated ionic strength. Furthermore, we provide mechanistic insights into mucin lubricity, which help explain how certain variations in physiologically relevant parameters can limit the lubricating potential of mucin solutions. Our results motivate that solutions of manually purified mucin solutions can be powerful biomedical lubricants, e.g. serving as eye drops, mouth sprays or as a personal lubricant for intercourse.
The ongoing effort to detect and characterize physical entanglement in biopolymers has so far established that knots are present in many globular proteins and also abound in viral DNA packaged inside bacteriophages. RNA molecules, on the other hand, have not yet been systematically screened for the occurrence of physical knots. We have accordingly undertaken the systematic profiling of the ~6,000 RNA structures present in the protein data bank. The search identified no more than three deeply-knotted RNA molecules. These are ribosomal RNAs solved by cryo-em and consist of about 3,000 nucleotides. Compared to the case of proteins and viral DNA, the observed incidence of RNA knots is therefore practically negligible. This suggests that either evolutionary selection, or thermodynamic and kinetic folding mechanisms act towards minimizing the entanglement of RNA to an extent that is unparalleled by other types of biomolecules. The properties of the three observed RNA knotting patterns provide valuable clues for designing RNA sequences capable of self-tying in a twist-knot fold.