A measure of bending in nucleic acids structures applied to A-tract DNA.

A method is proposed to measure global bending in DNA and RNA structures. It relies on a properly defined averaging of base-fixed coordinate frames, computes mean frames of suitably chosen groups of bases and uses these mean frames to evaluate bending. The method is applied to DNA A-tracts, known to induce considerable bend to the double helix. We performed atomistic molecular dynamics simulations of sequences containing the A(4)T(4) and T(4)A(4) tracts, in a single copy and in two copies phased with the helical repeat. Various temperature and salt conditions were investigated. Our simulations indicate bending by roughly 10 degrees per A(4)T(4) tract into the minor groove, and an essentially straight structure containing T(4)A(4), in agreement with electrophoretic mobility data. In contrast, we show that the published NMR structures of analogous sequences containing A(4)T(4) and T(4)A(4) tracts are significantly bent into the minor groove for both sequences, although bending is less pronounced for the T(4)A(4) containing sequence. The bending magnitudes obtained by frame averaging are confirmed by the analysis of superhelices composed of repeated tract monomers.

Different intrastrand and interstrand contributions to stacking account for roll variations at the alternating purine-pyrimidine sequences in A-DNA and A-RNA.

An explanation is suggested for the roll alternation between low and high values in A-type nucleic acid duplexes containing alternating sequences of purine and pyrimidine residues. The explanation combines two points. (1) Roll inevitably occurs in A-type duplexes due to geometrical reasons. (2) Intrastrand base stacking is much more impaired by roll than interstrand base stacking in A-type duplexes. Therefore purine-pyrimidine steps, whose bases mainly exhibit an intrastrand stacking, resist roll and decrease it. By contrast, bases at pyrimidine-purine steps exhibit a significant interstrand stacking that is tolerant to roll in A-type nucleic acid duplexes. In consequence, it is favourable if the purine-pyrimidine and pyrimidine-purine steps adopt low and high rolls, respectively in A-conformations of DNA and RNA molecules containing alternating purine-pyrimidine sequences. This is actually observed in the relevant molecular crystal structures.

Sequence-dependent elastic properties of DNA.

Harmonic elastic constants of 3-11 bp duplex DNA fragments were evaluated using four 5 ns unrestrained molecular dynamics simulation trajectories of 17 bp duplexes with explicit inclusion of solvent and counterions. All simulations were carried out with the Cornell et al. force-field and particle mesh Ewald method for long-range electrostatic interactions. The elastic constants including anisotropic bending and all coupling terms were derived by analyzing the correlations of fluctuations of structural properties along the trajectories. The following sequences have been considered: homopolymer d(ApA)(n) and d(GpG)(n), and alternating d(GPC)(n) and d(APT)(n). The calculated values of elastic constants are in very good overall agreement with experimental values for random sequences. The atomic-resolution molecular dynamics approach, however, reveals a pronounced sequence-dependence of the stretching and torsional rigidity of DNA, while sequence-dependence of the bending rigidity is smaller for the sequences considered. The earlier predicted twist-bend coupling emerged as the most important cross-term for fragments shorter than one helical turn. The calculated hydrodynamic relaxation times suggest that damping of bending motions may play a role in molecular dynamics simulations of long DNA fragments. A comparison of elasticity calculations using global and local helicoidal analyses is reported. The calculations reveal the importance of the fragment length definition. The present work shows that large-scale molecular dynamics simulations represent a unique source of data to study various aspects of DNA elasticity including its sequence-dependence.

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.

On the potential role of the amino nitrogen atom as a hydrogen bond acceptor in macromolecules.

Crystallographic studies of duplex DNA have indicated that opposing exocyclic amino groups may form close NH⋯:N contacts. To study the nature of such interactions, we have examined the database of small molecule, high-resolution crystal structures for more accurate examples of this type of unconventional interaction. We have found cases where the amino groups in guanine and adenine bases accept hydrogen bonds from conventional donors, such as amino or hydroxyl groups. More frequently, the purine amino group was found to contact closely electropositive C-H groups. Searches of the nucleic acid structural databases also yielded several examples where the purine amino group is contacted by hydrogen bond donors in macromolecules. Ab initio calculations indicate that the hydrogen-amino contact is improved energetically when the amino group moves from the conventional geometry, where all atoms are co-planar with the base, to one in which the hydrogen atoms lie out of the plane and the nitrogen is at the apex of a pyramid, resulting in polarization of the amino group. The combined structural and theoretical data suggest that the amino group is flexible, and can accommodate close contacts, because the resulting polarization permits electropositive atoms to approach the amino group nitrogen more closely than expected for their conventional van der Waals radii. The flexibility of the amino group may permit particular DNA conformations that enforce hydrogen-amino contacts to optimize favorable stacking interactions, and it may play a role in the recognition of nucleosides. We speculate that the amino group can accept hydrogen bonds under special circumstances in macromolecules, and that this ability might play a mechanistic role in catalytic processes such as deamination or amino transfer.

Can We Execute Stable Microsecond-Scale Atomistic Simulations of Protein-RNA Complexes?

We report over 30 µs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.

Theoretical analysis of the base stacking in DNA: choice of the force field and a comparison with the oligonucleotide crystal structures.

It follows from previous studies that changes in the base pair vertical separation (BPVS) influence the architecture of DNA much more than any other conformational parameter. This inspired us to compare BPVS in the available oligonucleotide crystal structures with the optimum values provided by nine different empirical potentials employed in the theoretical studies of DNA conformation. This comparison shows that BPVS is reproduced by three fields in all steps of the highly resolved oligonucleotide crystal structures while the remaining six empirical potentials, including AMBER, GROMOS and CHARMM, provide systematic deviations. We further find that the base pairs are poorly stacked (mostly compressed) in some other refined DNA crystal structures. Our analysis indicates that this poor stacking originates from improperly determined positions of the bases. The approach described in the present communication can be used to identify DNA structures which are not accurate enough for studies of the relationships between the base sequence and DNA conformation.

G.C base pair in parallel-stranded DNA--a novel type of base pairing: an ab initio quantum chemical study.

It was shown recently that parallel stranded DNA also incorporates the Reverse Watson-Crick Guanine. Cytosine (RWC G.C) pair. It contains only one hydrogen bond surrounded by mutual carbonyl group and amino group contacts. Here we describe an ab initio quantum chemical analysis of this unusual base pair. The conclusions can be summarized as follows. (1) The RWC G.C base pair exhibits a previously unknown type of DNA base pairing. (2) Its stabilization energy is estimated to be above 20 kJ/mol. A significant part of the stabilization originates in attractive interactions between the amino group hydrogen atom on one base and lone electron pair of the amino group nitrogen atom of the other (opposite) base. (3) In contrast to any other base pairing patterns, observed in the DNA molecules, the two DNA base amino groups in the RWC G.C base pair are highly nonplanar. They can form additional bifurcated hydrogen bonds with neighboring base pairs. Coplanarity of the two DNA base rings is energetically acceptable. (4) The unique RWC G.C pair cannot be analyzed by the empirical potential calculations.

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting.

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.

Relationships among rise, cup, roll and stagger in DNA suggested by empirical potential studies of base stacking.

Empirical analyses of experimental data have recently revealed a strong correlation in B- and A-DNA crystal structures between rise of the base pair steps and their cup, or the difference between roll and cup. We show here using empirical potentials that a major part of this correlation can be explained by the base stacking forces. Our calculations further demonstrate that the correlation depends on the base sequence while the dependence is strongest with the C-G step. We also show that small values (which lie beyond resolution of the X-ray diffraction data obtained with the DNA fragment single crystals) of base pair stagger can completely substitute for the effects of roll in the correlation. The present and our previous studies demonstrate that the base pair buckle and stagger can substantially affect base stacking in DNA so that variability of these parameters cannot be neglected in the theoretical analysis of the base sequence effects on DNA conformation.

Cation-pi and amino-acceptor interactions between hydrated metal cations and DNA bases. A quantum-chemical view.

Cation-pi interactions between cytosine and hexahydrated cations have been characterized using ab initio method with inclusion of electron correlation effects, assuming idealized and crystal geometries of the interacting species. Hydrated metal cations can interact with nucleobases in a cation-pi manner. The stabilization energy of such complexes would be large and comparable to the one for cation-pi complex with benzene. Further, polarized water molecules belonging to the hydration shell of the cation are capable to form a strong hydrogen bond interaction with the nitrogen lone electron pair of the amino groups of bases and enforce a pronounced sp3 pyramidalization of the nucleobase amino groups. However, in contrast to the benzene-cation complexes, the cation-pi configurations are highly unstable for a nucleobase since the conventional in plane binding of hydrated cations to the acceptor sites on the nucleobase is strongly preferred. Thus, a cation-pi interaction with a nucleobase can occur only if the position of the cation is locked above the nucleobase plane by another strong interaction. This indeed can occur in biopolymers and may have an effect on the local DNA architecture. Nevertheless, nucleobases have no intrinsic propensity to form cation-pi interactions.

Interactions between amino groups in DNA. An Ab initio study and a comparison with empirical potentials.

An ab initio quantum chemical analysis of the close amino group contacts, existing in many DNA crystal structures, is presented. The calculations are made at the Hartree-Fock (HF) level with medium 6-31G* and 6-31G(NH2*) basis sets as well as with inclusion of correlation energy using the second order Møller-Plesset theory (MP2) with the 6-31G* basis set. We demonstrate that the model system (methylamine dimer, cytosine dimer) amino groups are forced to adopt significantly non-planar geometry to stabilize their mutual interaction. Comparison is made with a representative set of empirical potentials including AMBER, CHARMM and GROMOS. The empirical potentials are not reliable enough to analyze the amino group contacts occurring in the DNA double helices. We propose that the mutual amino group interactions contribute to the conformational variability of the CpG and ApT B-DNA steps.

Hydrogen bonding and stacking of DNA bases: a review of quantum-chemical ab initio studies.

Ab initio quantum-chemical calulations with inclusion of electron correlation made since 1994 (such reliable calculations were not feasible before) significantly modified our view on interactions of nucleic acid bases. These calculations allowed to perform the first reliable comparison of the strength of stacked and hydrogen bonded pairs of nucleic acid bases, and to characterize the nature of the base-base interactions. Although hydrogen-bonded complexes of nucleobases are primarily stabilized by the electrostatic interaction, the dispersion attraction is also important. The stacked pairs are stabilized by dispersion attraction, however, the mutual orientation of stacked bases is determined rather by the electrostatic energy. Some popular theories of stacking were ruled out: The theory based on attractive interactions of polar exocyclic groups of bases with delocalized electrons of the aromatic rings (Bugg et al., Biopolymers 10, 175 (1971), and the pi-pi interactions model (C.A. Hunter, J. Mol. Biol. 230, 1025 (1993)). The calculations demonstrated that amino groups of nucleobases are very flexible and intrinsically nonplanar, allowing hydrogen-bond-like interactions which are oriented out of the plane of the nucleobase. Many H-bonded DNA base pairs are intrinsically nonplanar. Higher-level ab initio calculations provide a unique set of reliable and consistent data for parametrization and verification of empirical potentials. In this article, we present a short survey of the recent calculations, and discuss their significance and limitations. This summary is written for readers which are not experts in computational quantum chemistry.

Nonplanar DNA base pairs.

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp3 hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G-A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)...A(anti) pair of the ApG step of the d(CCAAGATTGG)2 crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.

Structure, energetics, vibrational frequencies and charge transfer of base pairs, nucleoside pairs, nucleotide pairs and B-DNA pairs of trinucleotides: ab initio HF/MINI-1 and empirical force field study.

Geometries, interaction energies and vibrational frequencies of base pairs, nucleoside pairs and nucleotide pairs were studied by ab initio Hartree-Fock (HF) method using MINI-1 basis set and empirical Cornell et al. force field (AMBER 4.1). A good agreement was found between HF/MINI-1 and AMBER results. In addition, both methods provide reasonable agreement with available high-level ab initio data. Finally, AMBER potential was used to determine the structure, energetics and vibrational frequencies of B-DNA pairs of trinucleotides. Stabilization energies of clusters are lowered when passing from base pairs to nucleoside pairs, nucleotide pairs and to pairs of trinucleotides. The lowest vibrations of base pairs and nucleoside pairs correspond to intermolecular motions of bases, specifically to buckle and propeller motions. In the case of pairs of larger subunits the lowest vibrations are of intramolecular nature (rotation around glycosidic bond, sugar and phosphate vibration). The spectra of these clusters became more complicated and quasi-degenerate. Intermolecular charge transfer between bases in H-bonded and stacked pairs is negligible, while a significant intramolecular charge transfer was observed.

Base stacking and hydrogen bonding in protonated cytosine dimer: the role of molecular ion-dipole and induction interactions.

An ab initio quantum-chemical study of stacked and hydrogen-bonded protonated cytosine dimer has been carried out. The calculations were made using the second-order Moller-Plesset perturbational method (MP2) with a medium-sized polarized set of atomic orbitals. H-bonded as well as stacked protonated base pairs are more stable than the neutral base pairs. Two energy contributions not present in the neutral base pairs stabilize the protonated base pairs: the molecular ion - dipole interaction, and the induction interaction. The molecular ion - dipole stabilization dominates in base pairs with highly polar neutral monomers, such as the C...CH+ base pair. The induction interaction is not included in the commonly used empirical potentials, which do not reproduce the changes in intermolecular stabilization due to protonation. We demonstrate that the base stacking of several consecutive C...CH + pairs, as proposed for polycytidylic acid and i-DNA, is strongly repulsive. We also show that the intermolecular interactions strongly prefer protonation of adenine in protonated adenine-cytosine pairs.

Hydrogen-bonded trimers of DNA bases and their interaction with metal cations: ab initio quantum-chemical and empirical potential study.

Neutral (G.GC, A.AT, G.AT, T.AT, and C(imino).GC) and protonated (CH+.GC and AH+.GC) hydrogen-bonded trimers of nucleic acid bases were characterized by ab initio methods with the inclusion of electron correlation. In addition, the influence of metal cations on the third-strand binding in Purine-Purine-Pyrimidine (Pu.PuPy) reverse-Hoogsteen triplets has been studied. The ab initio calculations were compared with those from recently introduced force fields (AMBER4.1, CHARMM23, and CFF95). The three-body term in neutral trimers is mostly negligible, and the use of empirical potentials is justified. The only exception is the neutral G.GC Hoogsteen trimer with a three-body term of -4 kcal/mol. Protonated trimers are stabilized by molecular ion-molecular dipole attraction and the interaction within the complex is nonadditive, with the three-body term on the order of -3 kcal/mol. There is a significant induction interaction between the third-strand protonated base and guanine. The calculations indicate an enhancement of the third-strand binding in the G.GC reverse-Hoogsteen trimer due to-metal cation coordination to the N7/O6 position of the third-strand guanine. Interactions between metal cations and complexes of DNA bases are in general highly non-additive; the three-body term is above-10 kcal/mol in a complex of a divalent cation (Ca2+) with the GG reverse-Hoogsteen pair. The pairwise additive empirical potentials qualitatively underestimate the binding energy between cation and base.

Interactions of hydrated IIa and IIb group metal cations with thioguanine-cytosine DNA base pair: Ab initio and density functional theory investigation of polarization effects, differences among cations, and flexibility of the cation hydration shell.

The structures and energies of the thioguanine-cytosine Watson-Crick (thioGC WC) base pair interacting with hydrated IIa (Mg2+, Ca2+, Ba2+) and IIb group (Zn2+, Cd2+, Hg2+) cations have been studied using ab initio techniques. Furthermore, complexes between guanine and thioguanine with hydrated cations have been characterized assuming various structures of the hydration shells. The complexes of the thioGC WC base pair with hydrated cations have similar properties as the previously studied GC WC base pair. There is substantial polarization stabilization of the base pairing due to cation binding which amounts to 7 - 11 kcal/mol. Soft Cd2+ and Hg2+ cations have a uniquely strong interaction with the thiogroup and induce substantial nonplanarity of the pairing. The thiogroup tends to reduce the number of water molecules in the first hydration shell of the cation. All complexes were optimized within the Hartree-Fock (HF) approximation while their energetics has been evaluated using the second-order Moller-Plesset perturbational method (MP2). All interaction energy evaluations and a substantial portion of the optimizations of the hydrated cation-(thio)guanine complexes have been repeated using Becke-3LYP Density Functional Theory method. All three approximations used (HF, Becke-3LYP, and MP2) give qualitatively the same results for the present cationic complexes. The results demonstrate specific differences among the cations and provide a set of reference structures and energies for verification and/or parametrization of empirical potentials and other theoretical methods.

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.

Close mutual contacts of the amino groups in DNA.

We analysed close contacts between the neighbouring base pairs in available oligonucleotide crystal structures and found that they mostly involved the amino groups of bases. Mutual amino group close contacts are abundant at the ApT steps in B-DNA but they also occur on the minor and major groove sides of the CpG steps in B-DNA. The close amino group contacts are much more frequent than the known bifurcated hydrogen bonds between the amino and carbonyl groups. The present empirical analysis indicates a previously unknown attractive interaction between the amino groups of bases in DNA which may significantly contribute to the base sequence-dependent variations of DNA architecture.