Temperature-dependent elasticity of DNA, RNA, and hybrid double helices.

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

RNA kink-turns are highly anisotropic with respect to lateral displacement of the flanking stems.

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

The temperature dependence of the helical twist of DNA.

DNA is the carrier of all cellular genetic information and increasingly used in nanotechnology. Quantitative understanding and optimization of its functions requires precise experimental characterization and accurate modeling of DNA properties. A defining feature of DNA is its helicity. DNA unwinds with increasing temperature, even for temperatures well below the melting temperature. However, accurate quantitation of DNA unwinding under external forces and a microscopic understanding of the corresponding structural changes are currently lacking. Here we combine single-molecule magnetic tweezers measurements with atomistic molecular dynamics and coarse-grained simulations to obtain a comprehensive view of the temperature dependence of DNA twist. Experimentally, we find that DNA twist changes by ΔTw(T) = (-11.0 ± 1.2)°/(°C·kbp), independent of applied force, in the range of forces where torque-induced melting is negligible. Our atomistic simulations predict ΔTw(T) = (-11.1 ± 0.3)°/(°C·kbp), in quantitative agreement with experiments, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates, while the effects of changes in substate populations are minor. Coarse-grained simulations using the oxDNA framework yield a value of ΔTw(T) = (-6.4 ± 0.2)°/(°C·kbp) in semi-quantitative agreement with experiments.

Interstrand cross-linking implies contrasting structural consequences for DNA: insights from molecular dynamics.

Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.

Long-timescale dynamics of the Drew-Dickerson dodecamer.

We present a systematic study of the long-timescale dynamics of the Drew-Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na(+)Cl(-) or K(+)Cl(-) The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 µs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 µs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.

Influence of BII Backbone Substates on DNA Twist: A Unified View and Comparison of Simulation and Experiment for All 136 Distinct Tetranucleotide Sequences.

Reliable representation of the B-DNA base-pair step twist is one of the crucial requirements for theoretical modeling of DNA supercoiling and other biologically relevant phenomena in B-DNA. It has long been suspected that the twist is inaccurately described by current empirical force fields. Unfortunately, comparison of simulation results with experiments is not straightforward because of the presence of BII backbone substates, whose populations may differ in experimental and simulation ensembles. In this work, we provide a comprehensive view of the effect of BII substates on the overall B-DNA helix twist and show how to reliably compare twist values from experiment and simulation in two scenarios. First, for longer DNA segments freely moving in solution, we show that sequence-averaged twists of different BI/BII ensembles can be compared directly because of approximate cancellation of the opposing BII effects. Second, for sequence-specific data, such as a particular base-pair step or tetranucleotide twist, can be compared only for a clearly defined BI/BII backbone conformation. For the purpose of force field testing, we designed a compact set of fourteen 22-base-pair B-DNA duplexes (Set 14) containing all 136 distinct tetranucleotide sequences and carried out a total of 84 µs of molecular dynamics simulations, primarily with the OL15 force field. Our results show that the ff99bsc0εζOL1χOL4, parmbsc1, and OL15 force fields model the B-DNA helical twist in good agreement with X-ray and minicircle ligation experiments. The comprehensive understanding obtained regarding the effect of BII substates on the base-pair step geometry should aid meaningful comparisons of various conformational ensembles in future research.

Explaining the striking difference in twist-stretch coupling between DNA and RNA: A comparative molecular dynamics analysis.

Double stranded helical DNA and RNA are flexible molecules that can undergo global conformational fluctuations. Their bending, twisting and stretching deformabilities are of similar magnitude. However, recent single-molecule experiments revealed a striking qualitative difference indicating an opposite sign for the twist-stretch couplings of dsDNA and dsRNA [Lipfert et al. 2014. Proc. Natl. Acad. Sci. U.S.A. 111, 15408] that is not explained by existing models. Employing unconstrained Molecular Dynamics (MD) simulations we are able to reproduce the qualitatively different twist-stretch coupling for dsDNA and dsRNA in semi-quantitative agreement with experiment. Similar results are also found in simulations that include an external torque to induce over- or unwinding of DNA and RNA. Detailed analysis of the helical deformations coupled to twist indicate that the interplay of helical rise, base pair inclination and displacement from the helix axis upon twist changes are responsible for the different twist-stretch correlations. Overwinding of RNA results in more compact conformations with a narrower major groove and consequently reduced helical extension. Overwinding of DNA decreases the size of the minor groove and the resulting positive base pair inclination leads to a slender and more extended helical structure.

On the Use of Molecular Dynamics Simulations for Probing Allostery through DNA.

A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning.

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

µABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

Insights into the structure of intrastrand cross-link DNA lesion-containing oligonucleotides: G[8-5m]T and G[8-5]C from molecular dynamics simulations.

Oxidatively generated complex DNA lesions occur more rarely than single-nucleotide defects, yet they play an important role in carcinogenesis and aging diseases because they have proved to be more mutagenic than simple lesions. Whereas their formation pathways are rather well understood, the field suffers from the absence of structural data that are crucial for interpreting the lack of repair. No experimental structures are available for oligonucleotides featuring such a lesion. Hence, the detailed structural basis of such damaged duplexes has remained elusive. We propose the use of explicit solvent molecular dynamics simulations to build up damaged oligonucleotides containing two intrastrand cross-link defects, namely, the guanine-thymine and guanine-cytosine defects. Each of these lesions, G[8-5m]T and G[8-5]C, is placed in the middle of a dodecameric sequence, which undergoes an important structural rearrangement that we monitor and analyze. In both duplexes, the structural evolution is dictated by the more favorable stacking of guanine G6, which aims to restore π-stacking with the 3' purine nucleobase. Subsequently, transient formation of hydrogen bonds with a strand shifting is observed. Our simulations are combined with density functional theory to rationalize the structural evolution. We report converging computational evidence that the G[8-5m]T- and G[8-5]C-containing structures evolve toward "abasic-like" duplexes, with a stabilization of the interstrand pairing noncovalent interactions. Meanwhile, both lesions restore B-helicity within tens of nanoseconds. The identification of plausible structures characterizes the last hydrogen abstraction step toward the formation of such defects as a non-innocent chemical reaction.

Effect O6-guanine alkylation on DNA flexibility studied by comparative molecular dynamics simulations.

Alkylation of guanine at the O6 atom is a highly mutagenic DNA lesion because it alters the coding specificity of the base causing G:C to A:T transversion mutations. Specific DNA repair enzymes, e.g. O(6)-alkylguanin-DNA-Transferases (AGT), recognize and repair such damage after looping out the damaged base to transfer it into the enzyme active site. The exact mechanism how the repair enzyme identifies a damaged site within a large surplus of undamaged DNA is not fully understood. The O(6)-alkylation of guanine may change the deformability of DNA which may facilitate the initial binding of a repair enzyme at the damaged site. In order to characterize the effect of O(6)-methyl-guanine (O(6)-MeG) containing base pairs on the DNA deformability extensive comparative molecular dynamics (MD) simulations on duplex DNA with central G:C, O(6)-MeG:C or O(6)-MeG:T base pairs were performed. The simulations indicate significant differences in the helical deformability due to the presence of O(6)-MeG compared to regular undamaged DNA. This includes enhanced base pair opening, shear and stagger motions and alterations in the backbone fine structure caused in part by transient rupture of the base pairing at the damaged site and transient insertion of water molecules. It is likely that the increased opening motions of O(6)-MeG:C or O(6)-MeG:T base pairs play a decisive role for the induced fit recognition or for the looping out of the damaged base by repair enzymes.

Structure and mechanical properties of the ribosomal L1 stalk three-way junction.

The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations.

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

Structure, dynamics, and interactions of a C4'-oxidized abasic site in DNA: a concomitant strand scission reverses affinities.

Apurinic/apyrimidinic (AP) sites constitute the most frequent form of DNA damage. They have proven to produce oxidative interstrand cross-links, but the structural mechanism of cross-link formation within a DNA duplex is poorly understood. In this work, we study three AP-containing d[GCGCGCXCGCGCG]·d[CGCGCGKGCGCGC] duplexes, where X = C, A, or G and K denotes an α,ß-unsaturated ketoaldehyde derived from elimination of a C4'-oxidized AP site featuring a 3' single-strand break. We use explicit solvent molecular dynamics simulations, complemented by quantum chemical density functional theory calculations on isolated X:K pairs. When X = C, the K moiety in the duplex flips around its glycosidic bond to form a stable C:K pair in a near-optimal geometry with two hydrogen bonds. The X = A duplex shows no stable interaction between K and A, which contrasts with AP sites lacking a strand scission that present a preferential affinity for adenine. Only one, transient G:K hydrogen bond is formed in the X = G duplex, although the isolated G:K pair is the most stable one. In the duplex, the stable C:K pair induces unwinding and sharp bending into the major groove at the lesion site, while the internal structure of the flanking DNA remains unperturbed. Our simulations also unravel transient hydrogen bonding between K and the cytosine 5' to the orphan base X = A. Taken together, our results provide a mechanistic explanation for the experimentally proven high affinity of C:K sites in forming cross-links in DNA duplexes and support experimental hints that interstrand cross-links can be formed with a strand offset.

Strikingly different effects of hydrogen bonding on the photodynamics of individual nucleobases in DNA: comparison of guanine and cytosine.

Ab initio surface hopping dynamics calculations were performed to study the photophysical behavior of cytosine and guanine embedded in DNA using a quantum mechanical/molecular mechanics (QM/MM) approach. It was found that the decay rates of photo excited cytosine and guanine were affected in a completely different way by the hydrogen bonding to the DNA environment. In case of cytosine, the geometrical restrictions exerted by the hydrogen bonds did not influence the relaxation time of cytosine significantly due to the generally small cytosine ring puckering required to access the crossing region between excited and ground state. On the contrary, the presence of hydrogen bonds significantly altered the photodynamics of guanine. The analysis of the dynamics indicates that the major contribution to the lifetime changes comes from the interstrand hydrogen bonds. These bonds considerably restricted the out-of-plane motions of the NH(2) group of guanine which are necessary for the ultrafast decay to the ground state. As a result, only a negligible amount of trajectories decayed into the ground state for guanine embedded in DNA within the simulation time of 0.5 ps, while for comparison, the isolated guanine relaxed to the ground state with a lifetime of about 0.22 ps. These examples show that, in addition to phenomena related to electronic interactions between nucleobases, there also exist relatively simple mechanisms in DNA by which the lifetime of a nucleobase is significantly enhanced as compared to the gas phase.

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA.

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein-DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50-100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

A comparative study of the binding of QSY 21 and Rhodamine 6G fluorescence probes to DNA: structure and dynamics.

Molecular dynamics (MD) simulations and ab initio quantum chemical calculations were employed to investigate the structure, dynamics and interactions of the QSY 21 nonfluorescent quencher and the fluorescence dye Rhodamine 6G bound to a B-DNA decamer. For QSY 21, two binding motifs were observed. In the first motif, the central xanthene ring is stacked on one base of the adjacent cytosine-guanine DNA base pair, whereas one of the 2,3-dihydro-1-indolyl aromatic side rings is stacked on the other base. In the second motif, the QSY 21 stacking interaction with the DNA base pair is mediated only by one of the side rings. Several transitions between the motifs are observed during a MD simulation. The ab initio calculations show that none of these motifs is energetically preferred. Two binding motifs were found also for Rhodamine 6G, with the xanthene ring stacked predominantly either on the cytosine or on the guanine. These results suggest that the side rings of QSY 21 play a crucial role in its stacking on the DNA and indicate novel binding mode absent in the case of Rhodamine 6G, which lacks aromatic side rings.

Towards a molecular dynamics consensus view of B-DNA flexibility.

We present a systematic study of B-DNA flexibility in aqueous solution using long-scale molecular dynamics simulations with the two more recent versions of nucleic acids force fields (CHARMM27 and parmbsc0) using four long duplexes designed to contain several copies of each individual base pair step. Our study highlights some differences between pambsc0 and CHARMM27 families of simulations, but also extensive agreement in the representation of DNA flexibility. We also performed additional simulations with the older AMBER force fields parm94 and parm99, corrected for non-canonical backbone flips. Taken together, the results allow us to draw for the first time a consensus molecular dynamics picture of B-DNA flexibility.