Base pair openings and temperature dependence of DNA flexibility

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Temperature-dependent elastic properties of DNA

Knowledge of the elastic properties, e.g., the persistence length or interphosphate distance, of single-stranded (ss) and double-stranded (ds) DNA under different experimental conditions is critical to characterizing molecular reactions studied with single-molecule techniques. While previous experiments have addressed the dependence of the elastic parameters upon varying ionic strength and contour length, temperature-dependent effects are less studied. Here, we examine the temperature-dependent elasticity of ssDNA and dsDNA in the range 5°C-50°C using a temperature-jump optical trap. We find a temperature softening for dsDNA and a temperature stiffening for ssDNA. Our results highlight the need for a general theory explaining the phenomenology observed.

Improving the mesoscopic modeling of DNA denaturation dynamics

Although previously developed mesoscopic DNA models have successfully reproduced thermodynamic denaturation data, recent studies show that these overestimate the rate of denaturation by orders of magnitude. Using adapted Peyrard-Bishop-Dauxois (PBD) models, we have calculated the denaturation rates of several DNA hairpins and made comparison with experimental data. We show that the addition of a barrier at the onsite potential of the PBD model gives a more accurate description of the unzipping dynamics of short DNA sequences. The new models provide a refined theoretical insight on the dynamical mechanisms of unzipping which can have implications for the understanding of transcription and replication processes. Still, this class of adapted PBD models seems to have a fundamental limitation which implies that it is not possible to get agreement with available experimental results on the dynamics of long DNA sequences and at the same time maintain the good agreement regarding its thermodynamics. The reason for this is that the denaturation rate of long DNA chains is not dramatically lowered by the additional barrier-as the base-pairs that open are more likely to remain open, facilitating the opening of the full DNA molecule. Some care has to be taken, since experimental techniques suitable to the study of denaturation rates of long sequences seem not to agree with other experimental data on short DNA sequences. Further research, both theoretical and experimental, is therefore needed to resolve these inconsistencies-which will be a starting point for new minimalistic models that are able to describe both thermodynamics and dynamics at a predictive level.

Tensile elasticity of semiflexible polymers with hinge defects

It has become clear in recent years that the simple uniform wormlike chain model needs to be modified in order to account for more complex behavior which has been observed experimentally in some important biopolymers. For example, the large flexibility of short ds-DNA has been attributed to kink or hinge defects. In this paper, we calculate analytically, within the weak bending approximation, the force-extension relation of a wormlike chain with a permanent hinge defect along its contour. The defect is characterized by its bending energy (which can be zero, in the completely flexible case) and its position along the polymer contour. Besides the bending rigidity of the chain, these are the only parameters which describe our model. We show that a hinge defect causes a significant increase in the differential tensile compliance of a prestressed chain. In the small force limit, a hinge defect significantly increases the entropic elasticity. Our results apply to any pair of semiflexible segments connected by a hinge. As such, they may also be relevant to cytoskeletal filaments (F-actin, microtubules), where one may treat the cross-link connecting two filaments as a hinge defect.

On the mechanical analogue of DNA

The creation of mechanical analogues of biological systems is known as a useful instrument that helps to understand better the dynamical mechanisms of the functioning of living organisms. Mechanical analogues of biomolecules are usually constructed for imitation of their internal mobility, which is one of the most important properties of the molecules. Among the different types of internal motions, angular oscillations of nitrous bases are of special interest because they make a substantial contribution to the base pairs opening that in turn is an important element of the process of the DNA-protein recognition. In this paper, we investigate the possibility to construct a mechanical analogue for imitation of angular oscillations of nitrous bases in inhomogeneous DNA. It is shown that the analogue has the form of a mechanical chain of non-identical pendulums that oscillate in the gravitational field of the Earth and coupled by identical springs. The masses and lengths of pendulums, as well as the distances between neighboring pendulums and the rigidity of springs are calculated. To illustrate the approach, we present the result of construction of the mechanical analogue of the fragment of the sequence of bacteriophage T7D.

Observation of coherent delocalized phonon-like modes in DNA under physiological conditions

Underdamped terahertz-frequency delocalized phonon-like modes have long been suggested to play a role in the biological function of DNA. Such phonon modes involve the collective motion of many atoms and are prerequisite to understanding the molecular nature of macroscopic conformational changes and related biochemical phenomena. Initial predictions were based on simple theoretical models of DNA. However, such models do not take into account strong interactions with the surrounding water, which is likely to cause phonon modes to be heavily damped and localized. Here we apply state-of-the-art femtosecond optical Kerr effect spectroscopy, which is currently the only technique capable of taking low-frequency (GHz to THz) vibrational spectra in solution. We are able to demonstrate that phonon modes involving the hydrogen bond network between the strands exist in DNA at physiologically relevant conditions. In addition, the dynamics of the solvating water molecules is slowed down by about a factor of 20 compared with the bulk.

Rigidity of melting DNA

The temperature dependence of DNA flexibility is studied in the presence of stretching and unzipping forces. Two classes of models are considered. In one case the origin of elasticity is entropic due to the polymeric correlations, and in the other the double-stranded DNA is taken to have an intrinsic rigidity for bending. In both cases single strands are completely flexible. The change in the elastic constant for the flexible case due to thermally generated bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic constant is found to be proportional to the square root of the bubble number fluctuation.

Flexibility of short DNA helices with finite-length effect: From base pairs to tens of base pairs

Flexibility of short DNA helices is important for the biological functions such as nucleosome formation and DNA-protein recognition. Recent experiments suggest that short DNAs of tens of base pairs (bps) may have apparently higher flexibility than those of kilo bps, while there is still the debate on such high flexibility. In the present work, we have studied the flexibility of short DNAs with finite-length of 5-50 bps by the all-atomistic molecular dynamics simulations and Monte Carlo simulations with the worm-like chain model. Our microscopic analyses reveal that short DNAs have apparently high flexibility which is attributed to the significantly strong bending and stretching flexibilities of ∼6 bps at each helix end. Correspondingly, the apparent persistence length lp of short DNAs increases gradually from ∼29 nm to ∼45 nm as DNA length increases from 10 to 50 bps, in accordance with the available experimental data. Our further analyses show that the short DNAs with excluding ∼6 bps at each helix end have the similar flexibility with those of kilo bps and can be described by the worm-like chain model with lp ∼ 50 nm.

Shape Persistence of Polyproline II Helical Oligoprolines

Oligoprolines are commonly used as molecular scaffolds. Past studies on the persistence length of their secondary structure, the polyproline II (PPII) helix, and on the fraction of backbone cis amide bonds have provided conflicting results. We resolved this debate by studying a series of spin-labeled proline octadecamers with EPR spectroscopy. Distance distributions between an N-terminal Gd(III) -DOTA (DOTA=1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) label and a nitroxide label at one of five evenly spaced backbone sites allowed us to discriminate between the flexibility of the PPII helix and the cis amide contributions. An upper limit of 2 % cis amide bonds per residue was found in a 7:3 (v/v) water/glycerol mixture, whereas cis amides were not observed in trifluoroethanol. Extrapolation of Monte Carlo models from the glass transition to ambient temperature predicts a persistence length of ≈3-3.5 nm in both solvents. The method is generally applicable to any type of oligomer for which the persistence length is of interest.

Temperature dependence of the DNA double helix at the nanoscale: structure, elasticity, and fluctuations

Biological organisms exist over a broad temperature range of -15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ ≃ +0.03) in the interval 0-100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.

Elasticity of cross-linked semiflexible biopolymers under tension

Aiming at the mechanical properties of cross-linked biopolymers, we set up and analyze a model of two weakly bending wormlike chains subjected to a tensile force, with regularly spaced inter-chain bonds (cross-links) represented by harmonic springs. Within this model, we compute the force-extension curve and the differential stiffness exactly and discuss several limiting cases. Cross-links effectively stiffen the chain pair by reducing thermal fluctuations transverse to the force and alignment direction. The extra alignment due to cross-links increases both with growing number and with growing strength of the cross-links, and is most prominent for small force f. For large f, the additional, cross-link-induced extension is subdominant except for the case of linking the chains rigidly and continuously along their contour. In this combined limit, we recover asymptotically the elasticity of a weakly bending wormlike chain without constraints, stiffened by a factor of 4. The increase in differential stiffness can be as large as 100% for small f or large numbers of cross-links.

Stretching and bending fluctuations of short DNA molecules

Recent measurements of the distribution of end-to-end distance in short DNA molecules infer cooperative stretching fluctuations. The assumptions underlying the analysis can be questioned if transient, thermally induced defects producing a localized decrease in bending stiffness are present in thermal equilibrium, such as regions in which DNA melts locally (bubbles), sustains large-angle bends (kinks), or can locally transform into an alternative (S-DNA) state. We study a generalized discrete worm-like chain model for DNA, capable of describing these experiments, showing that the model yields accurate fits to available experimental data. Our results indicate that DNA bending arising from such localized defects, rather than solely stretching, can be an equal contributor to end-to-end distance fluctuations for 35-bp DNA and contributes nontrivially to such fluctuations at all scales below the persistence length. The analysis suggests that such fluctuations should exhibit a scale-dependent cooperativity, specifically relevant in determining the behavior of short chains, but which saturates rapidly to a length-independent value for longer DNA, to ensure a consistent physical description of DNA across multiple scales. Our approach provides a minimal, yet accurate, coarse-grained description of DNA at the subpersistence length scales of current experimental interest.

Biophysical and ultrastructural characterization of adeno-associated virus capsid uncoating and genome release

We describe biophysical and ultrastructural differences in genome release from adeno-associated virus (AAV) capsids packaging wild-type DNA, recombinant single-stranded DNA (ssDNA), or dimeric, self-complementary DNA (scDNA) genomes. Atomic force microscopy and electron microscopy (EM) revealed that AAV particles release packaged genomes and undergo marked changes in capsid morphology upon heating in physiological buffer (pH 7.2). When different AAV capsids packaging ss/scDNA varying in length from 72 to 123% of wild-type DNA (3.4 to 5.8 kb) were incrementally heated, the proportion of uncoated AAV capsids decreased with genome length as observed by EM. Genome release was further characterized by a fluorimetric assay, which demonstrated that acidic pH and high osmotic pressure suppress genome release from AAV particles. In addition, fluorimetric analysis corroborated an inverse correlation between packaged genome length and the temperature needed to induce uncoating. Surprisingly, scAAV vectors required significantly higher temperatures to uncoat than their ssDNA-packaging counterparts. However, externalization of VP1 N termini appears to be unaffected by packaged genome length or self-complementarity. Further analysis by tungsten-shadowing EM revealed striking differences in the morphologies of ssDNA and scDNA genomes upon release from intact capsids. Computational modeling and molecular dynamics simulations suggest that the unusual thermal stability of scAAV vectors might arise from partial base pairing and optimal organization of packaged scDNA. Our work further defines the biophysical mechanisms underlying adeno-associated virus uncoating and genome release.

Magnetic tweezers measurements of the nanomechanical stability of DNA against denaturation at various conditions of pH and ionic strength

The opening of DNA double strands is extremely relevant to several biological functions, such as replication and transcription or binding of specific proteins. Such opening phenomenon is particularly sensitive to the aqueous solvent conditions in which the DNA molecule is dispersed, as it can be observed by considering the classical dependence of DNA melting temperature on pH and salt concentration. In the present work, we report a single-molecule study of the stability of DNA against denaturation when subjected to changes in solvent. We investigated the appearance of DNA instability under specific external applied force and imposed twist values, which was revealed by an increase in the temporal fluctuations in the DNA extension. These fluctuations occur in the presence of a continuous interval of equilibrium states, ranging from a plectonemic state to a state characterized by denaturation bubbles. In particular, we observe the fluctuations only around a characteristic force value. Moreover, this characteristic force is demonstrated to be notably sensitive to variations in the pH and ionic strength. Finally, an extension of a theoretical model of plectoneme formation is used to estimate the average denaturation energy, which is found to be linearly correlated to the melting temperature of the DNA double strands.

Single-molecule study of the DNA denaturation phase transition in the force-torsion space

We use the "magnetic tweezers" technique to show the structural transitions that the DNA undergoes in the force-torsion space. In particular, we focus on the regions corresponding to negative supercoiling. These regions are characterized by the formation of the so-called denaturation bubbles, which play an essential role in the replication and transcription of DNA. We experimentally map the region of the force-torsion space where the denaturation takes place. We observe that large fluctuations in DNA extension occur at one of the boundaries of this region, i.e., when the formation of denaturation bubbles and of plectonemes compete. To describe the experiments, we introduce a suitable extension of the classical model. The model correctly describes the position of the denaturation regions, the transition boundaries, and the measured values of the DNA extension fluctuations.