Unrestrained Molecular Dynamics Simulations of [d(AT)5]2 Duplex in Aqueous Solution: Hydration and Binding of Sodium Ions in the Minor Groove

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.

An investigation of the ionic and solvation patterns of dsDNG versus dsDNA by use of molecular dynamics simulations

Molecular dynamics has been employed to analyze the counterion and water atmospheres around the deoxynucleic guanidine (DNG) duplex G(12)-C(12). These features are compared to those of DNA G(12)-C(12). The chloride counterions of cationic DNG demonstrate fewer penetrations and only fleeting residence times in the minor groove, as opposed to the multi-nanosecond visits seen by sodium ions in DNA minor grooves. The 10 ns simulations also show the differences in hydration patterns between the DNG and DNA duplexes.

The dynamics of the B-A transition of natural DNA double helices

The dynamics of the B-A transition of DNA double helices with different GC contents and various chain lengths has been characterized by an electric field pulse technique. The field-induced B-A reaction is separated from orientation effects using the magic angle technique. Amplitudes reflecting the B-A reaction are observed selectively in the limited range of ethanol contents, where CD spectra demonstrate the B-A transition. The maximum amplitude appears at 1-2% higher ethanol content than the center of the B-A transition observed by CD because electric field pulses induce a relatively large perturbation from the A- toward the B-form. The relaxation curves measured after pulse termination reflect a spectrum of up to three relaxation processes. For DNA's with approximately 50% GC, the main part of the amplitude ( approximately 75%) is associated with time constants of approximately 2 micros, and another major component appears with time constants of 50-100 micros. These relaxation effects have been observed for DNA samples with 859, 2629, 7160, and 48501 bp. The time constant associated with the main amplitude increases with decreasing GC content from approximately 2 micros at 50% GC to approximately 3 mus at 41% GC and approximately 10 micros at 0% GC at the center of the B-A transition. Model calculations on the kinetics of cooperative linear Ising lattices predict the appearance of a distinct maximum of the mean relaxation time at the center of the transition. The absence of such maximum in our experimental data indicates a low cooperativity of the B-A transition with a nucleation parameter of approximately 0.1. The rate of the B-A transition is lower by approximately 3 orders of magnitude than that predicted by molecular dynamics simulations.

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.

Ion motions in molecular dynamics simulations on DNA

Counterions play a significant role in DNA structure and function, and molecular dynamics (MD) simulations offer the prospect of detailed description of the dynamical structure of ions at the molecular level. However, the motions of mobile counterions are notably slow to converge in MD on DNA. Obtaining accurate and reliable MD simulations requires knowing just how much sampling is required for convergence of each of the properties of interest. To address this issue, MD on a d(CGCGAATTCGCG) duplex in a dilute aqueous solution of water and 22 Na+ counterions was performed until convergence was achieved. The calculated first shell ion occupancies and DNA-Na+ radial distribution functions were computed as a function of time to assess convergence, and compared with relaxation times of the DNA internal parameters shift, slide, rise, tilt, roll, and twist. The sequence dependence of fractional occupancies of ions in the major and minor grooves of the DNA is examined, and the possibility of correlation between ion proximity and DNA minor groove widths is investigated.

DNA A-tract bending in three dimensions: solving the dA4T4 vs. dT4A4 conundrum

DNA A-tracts have been defined as four or more consecutive A.T base pairs without a TpA step. When inserted in phase with the DNA helical repeat, bending is manifested macroscopically as anomalous migration on polyacrylamide gels, first observed >20 years ago. An unsolved conundrum is why DNA containing in-phase A-tract repeats of A(4)T(4) are bent, whereas T(4)A(4) is straight. We have determined the solution structures of the DNA duplexes formed by d(GCAAAATTTTGC) [A4T4] and d(CGTTTTAAAACG) [T4A4] with NH(4)(+) counterions by using NMR spectroscopy, including refinement with residual dipolar couplings. Analysis of the structures shows that the ApT step has a large negative roll, resulting in a local bend toward the minor groove, whereas the TpA step has a positive roll and locally bends toward the major groove. For A4T4, this bend is nearly in phase with bends at the two A-tract junctions, resulting in an overall bend toward the minor groove of the A-tract, whereas for T4A4, the bends oppose each other, resulting in a relatively straight helix. NMR-based structural modeling of d(CAAAATTTTG)(15) and d(GTTTTAAAAC)(15) reveals that the former forms a left-handed superhelix with a diameter of approximately 110 A and pitch of 80 A, similar to DNA in the nucleosome, whereas the latter has a gentle writhe with a pitch of >250 A and diameter of approximately 50 A. Results of gel electrophoretic mobility studies are consistent with the higher-order structure of the DNA and furthermore depend on the nature of the monovalent cation present in the running buffer.

Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies

Ab initio quantum-chemical calculations with inclusion of electron correlation significantly contributed to our understanding of molecular interactions of DNA and RNA bases. Some of the most important findings are introduced in the present overview: structures and energies of hydrogen bonded base pairs, nature of base stacking, interactions between metal cations and nucleobases, nonplanarity of isolated nucleobases and other monomer properties, tautomeric equilibria of nucleobases, out-of-plane hydrogen bonds and amino acceptor interactions. The role of selected molecular interactions in nucleic acids is discussed and representative examples where these interactions occur are given. Also, accuracy of density functional theory, semiempirical methods, distributed multipole analysis and empirical potentials is commented on. Special attention is given to our very recent reference calculations on base stacking and H-bonding. Finally, we briefly comment on the relationship between advanced ab initio quantum-chemical methods and large-scale explicit solvent molecular dynamics simulations of nucleic acids.

DNA dynamics in a water drop without counterions

Because of its polyionic character, the DNA double helix is stable and biologically active only in salty aqueous media where its charge is compensated by solvent counterions. Monovalent metal ions are ubiquitous in DNA environment, and they are usually considered as the possible driving force of sequence-dependent modulations of DNA structure that make it recognizable by proteins. In an effort to directly examine this hypothesis, MD simulations of DNA in a water drop surrounded by vacuum were carried out, which relieves the requirement of charge neutrality. Surprisingly, with zero concentration of counterions, a dodecamer DNA duplex appears metastable, and its structure remains similar to that observed in experiment, including the minor groove narrowing in the dodecamer d(CGCGAATTCGCG)(2) often considered as the most evident cation effect. It is suggested that the same computational approach will allow one to simulate dynamics of long DNA chains more efficiently than with periodical boundary conditions.

Simulations of nucleic acids and their complexes

Recent years have seen considerable progress in simulations of nucleic acids. Improvements in force fields, simulation techniques and protocols, and increasing computer power have all contributed to making nanosecond-scale simulations of both DNA and RNA commonplace. The results are already helping to explain how nucleic acids respond to their environment and to their base sequence and to reveal the factors underlying recognition processes by probing biologically important nucleic acid-protein interactions and medically important nucleic acid-drug complexation. This Account summarizes methodological progress and applications of molecular dynamics to nucleic acids over the past few years and tries to identify remaining challenges.

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

Electronic properties, hydrogen bonding, stacking, and cation binding of DNA and RNA bases

This review summarizes results concerning molecular interactions of nucleic acid bases as revealed by advanced ab initio quantum chemical (QM) calculations published in last few years. We first explain advantages and limitations of modern QM calculations of nucleobases and provide a brief history of this still rather new field. Then we provide an overview of key electronic properties of standard and selected modified nucleobases, such as their charge distributions, dipole moments, polarizabilities, proton affinities, tautomeric equilibria, and amino group hybridization. Then we continue with hydrogen bonding of nucleobases, by analyzing energetics of standard base pairs, mismatched base pairs, thio-base pairs, and others. After this, the nature of aromatic stacking interactions is explained. Also, nonclassical interactions in nucleic acids such as interstrand bifurcated hydrogen bonds, interstrand close amino group contacts, C [bond] H...O interbase contacts, sugar-base stacking, intrinsically nonplanar base pairs, out-of-plane hydrogen bonds, and amino-acceptor interactions are commented on. Finally, we overview recent calculations on interactions between nucleic acid bases and metal cations. These studies deal with effects of cation binding on the strength of base pairs, analysis of specific differences among cations, such as the difference between zinc and magnesium, the influence of metalation on protonation and tautomeric equlibria of bases, and cation-pi interactions involving nucleobases. In this review, we do not provide methodological details, as these can be found in our preceding reviews. The interrelation between advanced QM approaches and classical molecular dynamics simulations is briefly discussed.

Molecular dynamics of the frame-shifting pseudoknot from beet western yellows virus: the role of non-Watson-Crick base-pairing, ordered hydration, cation binding and base mutations on stability and unfolding

Molecular dynamics simulations of the frame-shifting pseudoknot from beet western yellows virus (BWYV, NDB file UR0004) were performed with explicit inclusion of solvent and counterions. In all, 33 ns of simulation were carried out, including 10 ns of the native structure with protonation of the crucial cytosine residue, C8(N3+). The native structure exhibited stable trajectories retaining all Watson-Crick and tertiary base-pairs, except for fluctuations or transient disruptions at specific sites. The most significant fluctuations involved the change or disruption of hydrogen-bonding between C8(N3+) and bases G12, A25, and C26, as well as disruption of the water bridges linking C8(N3+) with A25 and C26. To increase sampling of rare events, the native simulation was continued at 400 K. A partial, irreversible unfolding of the molecule was initiated by slippage of C8(N3+) relative to G12 and continued by sudden concerted changes in hydrogen-bonding involving A23, A24, and A25. These events were followed by a gradual loss of stacking interactions in loop 2. Of the Watson-Crick base-pairs, only the 5'-terminal pair of stem 1 dissociated at 400 K, while the trans sugar-edge/sugar-edge A20.G4 interaction remained surprisingly stable. Four additional room-temperature simulations were carried out to obtain insights into the structural and dynamic effects of selected mutations. In two of these, C8 was left unprotonated. Considerable local rearrangements occurred that were not observed in the crystal structure, thus confirming N3-protonation of C8 in the native molecule. We also investigated the effect of mutating C8(N3+) to U8, to correlate with experimental and phylogenetic studies, and of changing the G4 x C17 base-pair to A4 x U17 to weaken the trans sugar-edge interaction between positions 4 and 20 and to test models of unfolding. The simulations indicate that the C8 x G12 x C26 base-triple at the junction is the most labile region of the frame-shifting pseudoknot. They provide insights into the roles of the other non-Watson-Crick base-pairs in the early stages of unfolding of the pseudoknot, which must occur to allow readthrough of the message by the ribosome. The simulations revealed several critical, highly ordered hydration sites with close to 100 % occupancies and residency times of individual water molecules of up to 5 ns. Sodium cation coordination sites with occupancies above 50 % were also observed.

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.

Structural dynamics and cation interactions of DNA quadruplex molecules containing mixed guanine/cytosine quartets revealed by large-scale MD simulations

Large-scale molecular dynamics (MD) simulations have been utilized to study G-DNA quadruplex molecules containing mixed GCGC and all-guanine GGGG quartet layers. Incorporation of mixed GCGC quartets into G-DNA stems substantially enhances their sequence variability. The mixed quadruplexes form rigid assemblies that require integral monovalent cations for their stabilization. The interaction of cations with the all-guanine quartets is the leading contribution for the stability of the four-stranded assemblies, while the mixed quartets are rather tolerated within the structure. The simulations predict that two cations are preferred to stabilize a four-layer quadruplex stem composed of two GCGC and two all-guanine quartets. The distribution of cations in the structure is influenced by the position of the GCGC quartets within the quadruplex, the presence and arrangement of thymidine loops connecting the guanine/cytosine stretches forming the stems, and the cation type present (Na(+) or K(+)). The simulations identify multiple nanosecond-scale stable arrangements of the thymidine loops present in the molecules investigated. In these thymidine loops, several structured pockets are identified capable of temporarily coordinating cations. However, no stable association of cations to a loop has been observed. The simulations reveal several paths through the thymidine loop regions that can be followed by the cations when exchanging between the central ion channel in the quadruplex stem and the surrounding solvent. We have carried out 20 independent simulations while the length of simulations reaches a total of 90 ns, rendering this study one of the most extensive MD investigations carried out on nucleic acids so far. The trajectories provide a largely converged characterization of the structural dynamics of these four-stranded G-DNA molecules.

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.

Molecular dynamics of DNA quadruplex molecules containing inosine, 6-thioguanine and 6-thiopurine

The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.