Role of sugar re-puckering in the transition of A and B forms of DNA in solution. A molecular dynamics study

Base pair compositional variability influences DNA structural stability and tunes hydration thermodynamics and dynamics

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes.

Sequence-Dependent Properties of the RNA Duplex

Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

The core structure of phi29 prohead RNA (pRNA) is composed of three major helices organized into three-way junction pRNA (3WJ-pRNA) and has stout structural rigidity along the coaxial helices. Prohead RNAs of the other Bacillus subtilis bacteriophages such as GA1 and SF5 share similar secondary structure and function with phi29; whether these pRNAs have similar mechanical rigidity remains to be elucidated. In this study, we constructed the tertiary structures of GA1 and SF5 3WJ-pRNAs by comparative modeling. Both GA1 and SF5 3WJ-pRNAs adopt a similar structure, in which three helices are organized as the three-way junction and two of the three helices are stacked coaxially. Moreover, detailed structural features of GA1 and SF5 3WJ-pRNAs are also similar to those of phi29 3WJ-pRNA: all of the bases of the coaxial helices are paired, and all of the adenines in the junction region are paired, which eliminates the interference of A-minor tertiary interactions. Hence, the coaxial helices tightly join to each other, and the major groove between them is very narrow. Two Mg²⁺ ions can thus fit into this major groove and form double Mg clamps. A steered molecular dynamics simulation was used to study the mechanical properties of these 3WJ-pRNAs. Both GA1 and SF5 3WJ-pRNAs show strong resistance to applied force in the direction of their coaxial helices. Such mechanical stability can be attributed to the Mg clamps.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.

Force field parameters for rotation around chi torsion axis in nucleic acids

To raise the accuracy of the force field for nucleic acids, several parameters were elaborated, focusing on the rotation around chi torsion axis. The reliability of molecular dynamics (MD) simulation was significantly increased by improving the torsion parameters at C8--N9--C1'--X (X = H1', C2', O4') in A, G and those at C6--N1--C1'--X in C, T, and U. In this work, we constructed small models representing the chemical structure of A, G, C, T, and U, and estimated energy profile for chi-axis rotation by executing numerous quantum mechanical (QM) calculations. The parameters were derived by discrete Fourier transformation of the calculated QM data. A comparison in energy profile between molecular mechanical (MM) calculation and QM one shows that our presently derived parameters well reproduce the energy surface of QM calculation for all the above torsion terms. Furthermore, our parameters show a good performance in MD simulations of some nucleic acids. Hence, the present refinement of parameters will enable us to perform more accurate simulations for various types of nucleic acids.

Collective headgroup conformational transition in twisted micellar superstructures

Predictions on amphiphilic self-assemblies traditionally rely on considerations on molecular shape and charge of the surfactant. In the case of functional surfactants a more sophisticated toolbox becomes necessary to design amphiphiles encoding chemical functionalities that provide additional responsive properties to their self-assemblies. Here we report on a comprehensive and combined structural-spectroscopic characterization of 1,2-dilauroyl-phosphatidyl-adenosine (DLPA) micelles in phosphate buffer. The temperature dependence, more precisely the thermal history of the sample, is explicitly taken into account. The experimental data, supplemented with MD simulations, indicate the presence of two possible states at room temperature, characterized by distinctly different structural properties that depend on the thermal history of the sample. The twisted superstructures, produced by aging DLPA micelles through intermicellar assembly of locally cylindrical aggregates at room temperature, collapse upon warming at 35 °C, yielding aligned filaments and/or wormlike structures. The initial superstructures cannot be recovered by thermal inversion. The reason for this behaviour is that the thermal activation causes a redistribution of syn-anti conformations of adenosine headgroups, as indicated by spectroscopic results (NMR, CD, FTIR), which is then collectively frozen thanks to molecular constraints present in the aggregate.

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.

Theoretical study of large conformational transitions in DNA: the B<-->A conformational change in water and ethanol/water

We explore here the possibility of determining theoretically the free energy change associated with large conformational transitions in DNA, like the solvent-induced B<-->A conformational change. We find that a combination of targeted molecular dynamics (tMD) and the weighted histogram analysis method (WHAM) can be used to trace this transition in both water and ethanol/water mixture. The pathway of the transition in the A-->B direction mirrors the B-->A pathway, and is dominated by two processes that occur somewhat independently: local changes in sugar puckering and global rearrangements (particularly twist and roll) in the structure. The B-->A transition is found to be a quasi-harmonic process, which follows closely the first spontaneous deformation mode of B-DNA, showing that a physiologically-relevant deformation is in coded in the flexibility pattern of DNA.

Are Isolated Nucleic Acid Bases Really Planar? A Car−Parrinello Molecular Dynamics Study

Car-Parrinello molecular dynamics simulations of the flexibility of isolated DNA bases have been carried out. The comparison of lowest ring out-of-plane vibrations calculated by using MP2/cc-pvdz and BLYP/PW methods reveals that the DFT method with the plane wave basis set reasonably reproduces out-of-plane deformability of the pyrimidine ring in nucleic acid bases and could be used for reliable modeling of conformational flexibility of nucleobases. The conformational phase space of pyrimidine rings in thymine, cytosine, guanine, and adenine has been investigated by using the ab initio Car-Parrinello molecular dynamics method. It is demonstrated that all nucleic acid bases are highly flexible molecules and possess a nonplanar effective conformation of the pyrimidine ring despite the fact that the planar geometry corresponds to a minimum on the potential energy surface. The population of the planar geometry of the pyrimidine ring does not exceed 30%. Among the nonplanar conformations of the pyrimidine rings, the boat-like and half-chair conformations are the most populated.

Photochemical probing of the B--a conformational transition in a linearized pUC19 DNA and its polylinker region

We induced the B-to-A conformational transition by ethanol in a linearized pUC19 DNA. A primer extension method was used in combination with UV light irradiation to follow the transition, based on pausing of DNA synthesis due to the presence of damaged bases in the template. Primer extension data highly correlated with the results of another method monitoring the B-A transition, i.e. inhibition of restriction endonuclease cleavage of UV light-irradiated DNA. Primer extension enabled us to locate damaged nucleotides within the region of interest. Most damaged nucleotides were located in B-form trimers, exclusively containing both pyrimidine bases (TTC, TCT, CTC, and CTT), and in a cytosine tetramer. The amount of damaged bases decreased in the course of B-A transition. Some of the damage even disappeared in the A-form, which mainly concerns the C(4) and C(3) blocks. The cleavage was nearly restored in the A-form within this region (Eco88I). On the contrary the decrease of damage was less significant with thymine dimers, only dropping to 50-60% of the B-form level. Consequently, the cleavage with EcoRI and HindIII remained mostly as before the transition (75% and 60% of uncleaved DNA preserved). We found significant differences in the B- and A-form pattern of UV light-damaged bases within the same region (polylinker) of DNA embedded within long (plasmid) or short (127 bp fragment) DNA molecules. The B-A transition of the fragment was found less cooperative than with linearized plasmid, which was confirmed by both CD spectroscopy and restriction cleavage inhibition.

Relative Flexibility of DNA and RNA: a Molecular Dynamics Study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.

Titration in silico of reversible B <=> A transitions in DNA

Reversible transitions between the A- and B-forms of DNA are obtained in free molecular dynamics simulations of a single double helix immersed in a water drop with Na(+) counterions. The dynamics of the transitions agrees with their supposed cooperative character. In silico titration of the transitions was carried out by smooth variation of the drop size. The estimated range of hydration numbers corresponding to the transition roughly agrees with experimental data. The chain length dependence was studied for double helices from 6 to 16 base pairs. It appeared that the B --> A transition is hindered for DNA shorter than one helical turn. With increased NaCl concentration in the drop, stabilization of the B-form is observed accompanied by the salt crystallization. The results strongly suggest that the B --> A transition at low hydration is caused by Na(+) ions sandwiched between phosphate strands in the major groove and is driven by direct medium range electrostatic interactions. The role of the reduced water shell apparently consists of increasing the counterion concentration in the opening of the major groove. Analysis of the available experimental data suggests that this mechanism is perhaps generally responsible for the A/B polymorphism in DNA.

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.

Solution structure of a DNA duplex with a chiral alkyl phosphonate moiety

The solution structures of two DNA decamers of sequence d(CCACCpxGGAAC).(GTTCCGGTGG) with a chiral alkyl phosphonate moiety (px) have been determined using NMR and restrained molecular dynamics simulations and compared with the solution structure of the unmodified duplex. The (1)H NMR spectra of two samples with pure stereochemistry in the modified phosphate have been assigned. The structures of both diastereoisomers, as well as the unmodified control duplex, have been determined from NMR-derived distance and torsion angle constraints. Accurate distance constraints were obtained from a complete relaxation matrix analysis of the NOE intensities. The structures have been refined with state of the art molecular dynamics methods, including explicit solvent and applying the particle mesh Ewald method to properly evaluate the long-range electrostatic interactions. In both cases, the calculations converge to well-defined structures, with RMSDs of approximately 1 A. The resulting structures belong to the general B family of DNA structures, even though the presence of the alkyl phosphonate moiety induces some slight displacement to the A-form in the neighborhood of the modified phosphate. Partial neutralization of this phosphate and the steric effect of the alkyl moiety provoke moderate bending in the DNA. This effect is more pronounced in the S diastereoisomer, where the alkyl group points inwards to the double helix.

The effect of amino groups on the stability of DNA duplexes and triplexes based on purines derived from inosine

The effect of amino groups attached at positions 2 and 8 of the hypoxanthine moiety in the structure, reactivity and stability of DNA duplexes and triplexes is studied by means of quantum mechanical calculations, as well as extended molecular dynamics (MD) and thermodynamic integration (MD/TI) simulations. Theoretical estimates of the change in stability related to 2'-deoxyguanosine (G) --> 2'-deoxyinosine (I) --> 8-amino-2'-deoxyinosine (8AI) mutations have been experimentally verified, after synthesis of the corresponding compounds. An amino group placed at position 2 stabilizes the duplex, as expected, and surprisingly also the triplex. The presence of an amino group at position 8 of the hypoxanthine moiety stabilizes the triplex but, surprisingly, destabilizes the duplex. The subtle electronic redistribution occurring upon the introduction of an amino group on the purine seems to be responsible for this surprising behavior. Interesting 'universal base' properties are found for 8AI.

Molecular dynamics simulation of nucleic acids

We review molecular dynamics simulations of nucleic acids, including those completed from 1995 to 2000, with a focus on the applications and results rather than the methods. After the introduction, which discusses recent advances in the simulation of nucleic acids in solution, we describe force fields for nucleic acids and then provide a detailed summary of the published literature. We emphasize simulations of small nucleic acids ( approximately 6 to 24 mer) in explicit solvent with counterions, using reliable force fields and modern simulation protocols that properly represent the long-range electrostatic interactions. We also provide some limited discussion of simulation in the absence of explicit solvent. Absent from this discussion are results from simulations of protein-nucleic acid complexes and modified DNA analogs. Highlights from the molecular dynamics simulation are the spontaneous observation of A B transitions in duplex DNA in response to the environment, specific ion binding and hydration, and reliable representation of protein-nucleic acid interactions. We close by examining major issues and the future promise for these methods.

Misincorporation of 2'-deoxyoxanosine into DNA: a molecular basis for NO-induced mutagenesis derived from theoretical calculations

A wide range of theoretical methods, including high level ab initio, density functional, self-consistent reaction field, molecular dynamics and thermodynamic integration calculations, have been used to analyze the mutagenic properties of oxanosine. The major tautomeric forms in the gas phase and aqueous solution have been determined. The ability of oxanosine to recognize thymine and cytosine in the gas phase and in the DNA environment has been compared with that of guanine. A physicochemical explanation for the mutagenic properties of oxanosine is suggested.

Modeling DNA deformations

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.