No Arabic abstract
Far from being a passive information store, the genome is a mechanically dynamic and diverse system in which torsion and tension fluctuate and combine to determine structure and help regulate gene expression. Much of this mechanical perturbation is due to molecular machines such as topoisomerases which must stretch and twist DNA as part of various functions including DNA repair and replication. While the broad-scale mechanical response of nucleic acids to tension and torsion is well characterized, detail at the single base pair level is beyond the limits of even super-resolution imaging. Here, we present a straightforward, flexible, and extensible umbrella-sampling protocol to twist and stretch nucleic acids in silico using the popular biomolecular simulation package Amber -- though the principles we describe are applicable also to other packages such as GROMACS. We discuss how to set up the simulation system, decide forcefields and solvation models, and equilibrate. We then introduce the torsionally-constrained stretching protocol, and finally we present some analysis techniques we have used to characterize structural motif formation. Rather than define forces or fictional pseudoatoms, we instead define a fixed translation of specified atoms between each umbrella sampling step, which allows comparison with experiment without needing to estimate applied forces by simply using the fractional end-to-end displacement as a comparison metric. We hope that this easy to implement solution will be valuable for interrogating optical and magnetic tweezers data on nucleic acids at base pair resolution.
Loops are essential secondary structure elements in folded DNA and RNA molecules and proliferate close to the melting transition. Using a theory for nucleic acid secondary structures that accounts for the logarithmic entropy c ln m for a loop of length m, we study homopolymeric single-stranded nucleic acid chains under external force and varying temperature. In the thermodynamic limit of a long strand, the chain displays a phase transition between a low temperature / low force compact (folded) structure and a high temperature / high force molten (unfolded) structure. The influence of c on phase diagrams, critical exponents, melting, and force extension curves is derived analytically. For vanishing pulling force, only for the limited range of loop exponents 2 < c < 2.479 a melting transition is possible; for c <= 2 the chain is always in the folded phase and for 2.479 < c always in the unfolded phase. A force induced melting transition with singular behavior is possible for all loop exponents c < 2.479 and can be observed experimentally by single molecule force spectroscopy. These findings have implications for the hybridization or denaturation of double stranded nucleic acids. The Poland-Scheraga model for nucleic acid duplex melting does not allow base pairing between nucleotides on the same strand in denatured regions of the double strand. If the sequence allows these intra-strand base pairs, we show that for a realistic loop exponent c ~ 2.1 pronounced secondary structures appear inside the single strands. This leads to a lower melting temperature of the duplex than predicted by the Poland-Scheraga model. Further, these secondary structures renormalize the effective loop exponent c^, which characterizes the weight of a denatured region of the double strand, and thus affect universal aspects of the duplex melting transition.
EXSY, TOCSY and NOESY lie at the foundation of homonuclear NMR experiments in organic and pharmaceutical chemistry, as well as in structural biology. Limited magnetization transfer efficiency is an intrinsic downside of these methods, particularly when targeting rapidly exchanging species such as labile protons ubiquitous in polysaccharides, sidechains and backbones of proteins, and in bases and sugars of nucleic acids: the fast decoherence imparted on these protons through solvent exchanges, greatly reduces their involvement in homonuclear correlation experiments. We have recently discussed how these decoherences can be visualized as an Anti-Zeno Effect, that can be harnessed to enhance the efficiency of homonuclear transfers within Looped PROjected SpectroscopY (L-PROSY) leading to 200-300% enhancements in NOESY and TOCSY cross-peaks for amide groups in biomolecules. This study demonstrates that even larger sensitivity gains per unit time, equivalent to reductions by several hundred-folds in the duration of experiments, can be achieved by looping inversion or using saturation procedures. In the ensuing experiments a priori selected frequencies are encoded according to Hadamard recipes, and subsequently resolved along the indirect dimension via linear combinations. Magnetization-transfer (MT) processes reminiscent of those occurring in CEST provide significant enhancements in the resulting cross-peaks, in only a fraction of acquisition time of a normal 2D experiment. The effectiveness of the ensuing three-way polarization transfer interplay between water, labile and non-labile protons was corroborated experimentally for proteins, homo-oligosaccharides and nucleic acids. In all cases, cross-peaks barely detectable in conventional 2D NMR counterparts, were measured ca. 10-fold faster and with 200-600% signal enhancements by the Hadamard MT counterparts.
In several previous works, I presented the mirror symmetry in the set of protein amino acids, expressed through the number of atoms. Here, however, the same thing is shown but over the number of nucleons and molecules mass. Compared to the previous version of the paper, minimal changes have been made, and Display 2 as well as Figures 3 and 4 have been added.
Helicases are molecular motors that unwind double-stranded nucleic acids (dsNA), such as DNA and RNA). Typically a helicase translocates along one of the NA single strands while unwinding and uses adenosine triphosphate (ATP) hydrolysis as an energy source. Here we model of a helicase motor that can switch between two states, which could represent two different points in the ATP hydrolysis cycle. Our model is an extension of the earlier Betterton-Julicher model of helicases to incorporate switching between two states. The main predictions of the model are the speed of unwinding of the dsNA and fluctuations around the average unwinding velocity. Motivated by a recent claim that the NS3 helicase of Hepatitis C virus follows a flashing ratchet mechanism, we have compared the experimental results for the NS3 helicase with a special limit of our model which corresponds to the flashing ratchet scenario. Our model accounts for one key feature of the experimental data on NS3 helicase. However, contradictory observations in experiments carried out under different conditions limit the ability to compare the model to experiments.
Optimal hematocrit $H_o$ maximizes oxygen transport. In healthy humans, the average hematocrit $H$ is in the range of 40-45$%$, but it can significantly change in blood pathologies such as severe anemia (low $H$) and polycythemia (high $H$). Whether the hematocrit level in humans corresponds to the optimal one is a long standing physiological question. Here, using numerical simulations with the Lattice Boltzmann method and two mechanical models of the red blood cell (RBC) we predict the optimal hematocrit, and explore how altering the mechanical properties of RBCs affects $H_o$. We develop a simplified analytical theory that accounts for results obtained from numerical simulations and provides insight into the physical mechanisms determining $H_o$. Our numerical and analytical models can easily be modified to incorporate a wide range of mechanical properties of RBCs as well as other soft particles thereby providing means for the rational design of blood substitutes. Our work lays the foundations for systematic theoretical study of the optimal hematocrit and its link with pathological RBCs associated with various diseases (e.g. sickle cell anemia, diabetes mellitus, malaria, elliptocytosis).