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.

Impact of the Nucleosome Histone Core on the Structure and Dynamics of DNA-Containing Pyrimidine-Pyrimidone (6-4) Photoproduct.

The pyrimidine-pyrimidone (6-4) photoproduct (64-PP) is an important photoinduced DNA lesion constituting a mutational signature for melanoma. The structural impact of 64-PP on DNA complexed with histones affects the lesion mutagenicity and repair but remains poorly understood. Here we investigate the conformational dynamics of DNA-containing 64-PP within the nucleosome core particle by atomic-resolution molecular dynamics simulations and multiscale data analysis. We demonstrate that the histone core exerts important mechanical restraints that largely decrease global DNA structural fluctuations. However, the local DNA flexibility at the damaged site is enhanced due to imperfect structural adaptation to restraints imposed by the histone core. If 64-PP faces the histone core and is therefore not directly accessible by the repair protein, the complementary strand facing the solvent is deformed and exhibits higher flexibility than the corresponding strand in a naked, undamaged DNA. This may serve as an initial recognition signal for repair. Our simulations also pinpoint the structural role of proximal residues from the truncated histone tails.

Compensatory Mechanisms in Temperature Dependence of DNA Double Helical Structure: Bending and Elongation.

Changes in the structure of double-stranded (ds) DNA with temperature affect processes in thermophilic organisms and are important for nanotechnological applications. Here we investigate temperature-dependent conformational changes of dsDNA at the scale of several helical turns and at the base pair step level, inferred from extensive all-atom molecular dynamics simulations of DNA at temperatures from 7 to 47 °C. Our results suggest that, contrary to twist, the overall bending of dsDNA without A-tracts depends only very weakly on temperature, due to the mutual compensation of directional local bends. Investigating DNA length as a function of temperature, we find that the sum of distances between base pair centers (the wire length) exhibits a large expansion coefficient of ∼2 × 10-4 °C-1, similar to values reported for thermoplastic materials. However, the wire length increase with temperature is absorbed by expanding helix radius, so the length measured along the helical axis (the spring length) seems to suggest a very small negative thermal expansion coefficient. These compensatory mechanisms contribute to thermal stability of DNA structure on the biologically relevant scale of tens of base pairs and longer.

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.

rRNA C-Loops: Mechanical Properties of a Recurrent Structural Motif.

C-loop is an internal loop motif found in the ribosome and used in artificial nanostructures. While its geometry has been partially characterized, its mechanical properties remain elusive. Here we propose a method to evaluate global shape and stiffness of an internal loop. The loop is flanked by short A-RNA helices modeled as rigid bodies. Their relative rotation and displacement are fully described by six interhelical coordinates. The deformation energy of the loop is assumed to be a general quadratic function of the interhelical coordinates. The model parameters for isolated C-loops are inferred from unrestrained all-atom molecular dynamics simulations. C-loops exhibit high twist as reported earlier, but also a bend and a lateral displacement of the flanking helices. Their bending stiffness and lateral displacement stiffness are nearly isotropic and similar to the control A-RNA duplexes. Nevertheless, we found systematic variations with the C-loop position in the ribosome and the organism of origin. The results characterize global properties of C-loops in the full six-dimensional interhelical space and enable one to choose an optimally stiff C-loop for use in a nanostructure. Our approach can be readily applied to other internal loops and extended to more complex structural motifs.

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.

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.

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.

Multiscale modelling of DNA mechanics.

Mechanical properties of DNA are important not only in a wide range of biological processes but also in the emerging field of DNA nanotechnology. We review some of the recent developments in modeling these properties, emphasizing the multiscale nature of the problem. Modern atomic resolution, explicit solvent molecular dynamics simulations have contributed to our understanding of DNA fine structure and conformational polymorphism. These simulations may serve as data sources to parameterize rigid base models which themselves have undergone major development. A consistent buildup of larger entities involving multiple rigid bases enables us to describe DNA at more global scales. Free energy methods to impose large strains on DNA, as well as bead models and other approaches, are also briefly discussed.

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.

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.

Mechanical Model of DNA Allostery.

The importance of allosteric effects in DNA is becoming increasingly appreciated, but the underlying mechanisms remain poorly understood. In this work, we propose a general modeling framework to study DNA allostery. We describe DNA in a coarse-grained manner by intra-base pair and base pair step coordinates, complemented by groove widths. Quadratic deformation energy is assumed, yielding linear relations between the constraints and their effect. Model parameters are inferred from standard unrestrained, explicit-solvent molecular dynamics simulations of naked DNA. We applied the approach to study minor groove binding of diamidines and pyrrole-imidazole polyamides. The predicted DNA bending is in quantitative agreement with experiment and suggests that diamidine binding to the alternating TA sequence brings the DNA closer to the A-tract conformation, with potentially important functional consequences. The approach can be readily applied to other allosteric effects in DNA and generalized to model allostery in various molecular systems.

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.

Effect of 8-oxoguanine on DNA structure and deformability.

8-Oxoguanine (oxoG) is an abundant product of oxidative DNA damage. It is removed by repair glycosylases, but exactly how the enzymes recognize oxoG in the large surplus of undamaged bases is not fully understood. The lesion may induce changes in the properties of naked DNA that facilitate the recognition. In this work, we assess the effect of oxoG on DNA structure and mechanical deformability. We performed extensive unrestrained, atomic resolution molecular dynamics simulations to parametrize a nonlocal, rigid base mechanical model of DNA. Our data indicate that oxoG induces unwinding of the base pair step at the 5'-side of the lesion. This brings the damaged DNA closer to its conformation in the initial complex with bacterial glycosylase MutM. The untwisting is partially caused by different BII substate populations and is further enhanced by the base-sugar repulsion within oxoG. On the other hand, our analysis shows that damaged and undamaged DNA have very similar harmonic stiffness. These results suggest an indirect readout component of the MutM-DNA initial complex formation. They also help one to understand the effect of oxoG on the formation of nucleosomes and looped gene regulatory complexes.

Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer.

The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]2 is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.

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.