Conformational dynamics and internal friction in homo-polymer globules: equilibrium vs. non-equilibrium simulations

We study the conformational dynamics within homo-polymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength {epsilon} and the globule size NG is observed. We find two distinct dynamical regimes: a liquid- like regime (for {epsilon} < {epsilon}s) with fast internal dynamics and a solid-like regime (for {epsilon} > {epsilon}s) with slow internal dynamics. The cohesion strength {epsilon}s of this freezing transition depends on NG. Equilibrium simulations, where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompassing low enough velocities so that the linear-response, viscous regime is reached). From both simulation protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal friction increases continuously with {epsilon} and scales extensive in NG. This suggests an internal friction scenario where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the globular structure.

Secondary structure formation of homopolymeric single-stranded nucleic acids including force and loop entropy: implications for DNA hybridization

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.

Theory for RNA folding, stretching, and melting including loops and salt

Secondary structure formation of nucleic acids strongly depends on salt concentration and temperature. We develop a theory for RNA folding that correctly accounts for sequence effects, the entropic contributions associated with loop formation, and salt effects. Using an iterative expression for the partition function that neglects pseudoknots, we calculate folding free energies and minimum free energy configurations based on the experimentally derived base pairing free energies. The configurational entropy of loop formation is modeled by the asymptotic expression -c ln m, where m is the length of the loop and c the loop exponent, which is an adjustable constant. Salt effects enter in two ways: first, we derive salt induced modifications of the free energy parameters for describing base pairing and, second, we include the electrostatic free energy for loop formation. Both effects are modeled on the Debye-Hueckel level including counterion condensation. We validate our theory for two different RNA sequences: For tRNA-phe, the resultant heat capacity curves for thermal denaturation at various salt concentrations accurately reproduce experimental results. For the P5ab RNA hairpin, we derive the global phase diagram in the three-dimensional space spanned by temperature, stretching force, and salt concentration and obtain good agreement with the experimentally determined critical unfolding force. We show that for a proper description of RNA melting and stretching, both salt and loop entropy effects are needed.

Impact of loop statistics on the thermodynamics of RNA folding

Loops are abundant in native RNA structures and proliferate close to the unfolding transition. By including a statistical weight ~ l^{-c} for loops of length l in the recursion relation for the partition function, we show that the calculated heat capacity depends sensitively on the presence and value of the exponent c, even of short t-RNA. For homo-RNA we analytically calculate the critical temperature and critical exponents which exhibit a non-universal dependence on c.