Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

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.

Knowledge-based elastic potentials for docking drugs or proteins with nucleic acids

Elastic ellipsoidal functions defined by the observed hydration patterns around the DNA bases provide a new basis for measuring the recognition of ligands in the grooves of double-helical structures. Here a set of knowledge-based potentials suitable for quantitative description of such behavior is extracted from the observed positions of water molecules and amino acid atoms that form hydrogen bonds with the nitrogenous bases in high resolution crystal structures. Energies based on the displacement of hydrogen-bonding sites on drugs in DNA-crystal complexes relative to the preferred locations of water binding around the heterocyclic bases are low, pointing to the reliability of the potentials and the apparent displacement of water molecules by drug atoms in these structures. The validity of the energy functions has been further examined in a series of sequence substitution studies based on the structures of DNA bound to polyamides that have been designed to recognize the minor-groove edges of Watson-Crick basepairs. The higher energies of binding to incorrect sequences superimposed (without conformational adjustment or displacement of polyamide ligands) on observed high resolution structures confirm the hypothesis that the drug subunits associate with specific DNA bases. The knowledge-based functions also account satisfactorily for the measured free energies of DNA-polyamide association in solution and the observed sites of polyamide binding on nucleosomal DNA. The computations are generally consistent with mechanisms by which minor-groove binding ligands are thought to recognize DNA basepairs. The calculations suggest that the asymmetric distributions of hydrogen-bond-forming atoms on the minor-groove edge of the basepairs may underlie ligand discrimination of G.C from C.G pairs, in addition to the commonly believed role of steric hindrance. The analysis of polyamide-bound nucleosomal structures reveals other discrepancies in the expected chemical design, including unexpected contacts to DNA and modified basepair targets of some ligands. The ellipsoidal potentials thus appear promising as a mathematical tool for the study of drug- and protein-DNA interactions and for gaining new insights into DNA-binding mechanisms.

Theoretical evaluation of dielectric absorption of microwave energy at the scale of nucleic acids

A theoretical model is proposed for the evaluation of dielectric properties of the cell nucleus between 0.3 and 3 GHz, as a function of its nucleic acids (NA) concentration (CNA). It is based on literature data on dielectric properties of DNA solutions and nucleoplasm. In skeletal muscle cells, the specific absorption rate (SAR) ratio between nucleoplasm and cytoplasm is found to be larger than one for CNA above 30 mg/ml. A nearly linear relationship is found between CNA and this nucleocytoplasmic SAR ratio. Considering the nanoscale of the layer of condensed counterions and bound water molecules at the NA-solution interface, the power absorption per unit volume is evaluated at this precise location. It is found to be between one and two orders of magnitude above that in muscle tissue as a whole. Under realistic microwave (MW) exposure conditions, however, these SAR inhomogeneities do not generate any significant thermal gradient at the scale considered here. Nevertheless, the question arises of a possible biological relevance of nonnegligible and preferential heat production at the location of the cell nucleus and of the NA molecules.

Probing hydration of monovalent cations condensed around polymeric nucleic acids

We report the relative molar sound velocity increments, [U], partial molar volumes, V(o), and partial molar adiabatic compressibilities, K(S)(o), of the Li(+), Na(+), K(+), Rb(+), Cs(+), NH(4)(+), and N(CH(3))(4)(+) salts of poly(dAdT)poly(dAdT), poly(dGdC)poly(dGdC), poly(dIdC)poly(dIdC), poly(rA)poly(rU), poly(rG)poly(rC), poly(rI)poly(rC), and poly(rU) at 25 degrees C. When analyzing these data, we take into account the Donnan membrane equilibrium effect. Comparison between the values of [U], V(o), and K(S)(o) exhibited by the nucleic acid salts and respective chlorides (LiCl, NaCl, KCl, RbCl, CsCl, NH(4)Cl, and N(CH(3))(4)Cl) yields information about the state of counterion hydration in the vicinity of each nucleic acid structure studied here. Our analysis reveals that the poly(dGdC)poly(dGdC), poly(dIdC)poly(dIdC), and poly(rI)poly(rC) duplexes and single-stranded poly(rU) do not significantly influence the hydration properties of their condensed counterions. In the vicinity of these polymers, counterions retain their full hydration shells (within +/-15%). By contrast, counterions condensed around the poly(dAdT)poly(dAdT), poly(rA)poly(rU), and poly(rG)poly(rC) duplexes are significantly dehydrated and retain, respectively, only 65(+/-18)%, 34(+/-21)%, and 33(+/-9)% of their original hydration shells. Taken together, the volumetric data reported here provide important new information that ultimately may help us understand the central role that hydration and counterions play in modulating the conformational preferences of nucleic acids and the energetics of DNA recognition events.

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

Influence of Spermine on DNA Conformation in a Molecular Dynamics Trajectory of d(CGCGAATTCGCG)2: Major Groove Binding by One Spermine Molecule Delays the A→B Transition

The effect of spermine on the A-DNA to B-DNA transition in d(CGCGAATTCGCG)(2) has been investigated by five A-start molecular dynamics simulations, using the Cornell et al. potential. In the absence of spermine an A-->B transition is initiated immediately and the DNA becomes equidistant from the A- and B-forms at 200ps. In three DNA-spermine simulations, when a spermine is located across the major groove of A-DNA in one of three different initial locations, the time taken to reach equidistance from the A- and B-forms is delayed until 800, 950 or 1000ps. In each case the A-form appears to be temporarily stabilized by spermine's electrostatic interactions with phosphates on both sides of the major groove. The onset of the A-->B transition can be correlated with the spermine losing contact with phosphates on one side of the groove and with A-like --> B-like sugar pucker transitions in the vicinity of the spermine bridge. However in the fifth trajectory, in which the spermine initially threads from the major groove via the backbone into the minor groove, the B-->A transition occurs rapidly once again and the DNA is equidistant between the A- and B-forms within 300ps. This indicates that the mere presence of spermine is insufficient to delay the transition and that major groove binding stabilizes A-DNA.

Water-mediated contacts in thetrp-repressor operator complex recognition process

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.

Atomic-resolution STM structure of DNA and localization of the retinoic acid binding site

Single-molecule imaging by scanning tunnelling microscopy (STM) yields the atomic-resolution (0.6A) structure of individual B-type DNA molecules. The strong correlation between these STM structures and those predicted from the known base sequence indicates that sequencing of single DNA molecules using STM may be feasible. There is excellent agreement between the STM and X-ray structures, but subtle differences exist due to radial distortions. We show that the interactions of other molecules with DNA, their binding configurations, and the structure of these complexes can be studied at the single-molecule level. The anti-cancer drug retinoic acid (RA) binds selectively to the minor groove of DNA with up to 6 RA molecules per DNA turn and with the plane of the RA molecule approximately parallel to the DNA symmetry axis. Similar studies for other drug molecules will be valuable in the a priori evaluation of the effectiveness of anti-cancer drugs.

DNA polymorphism: a comparison of force fields for nucleic acids

The improvements of the force fields and the more accurate treatment of long-range interactions are providing more reliable molecular dynamics simulations of nucleic acids. The abilities of certain nucleic acid force fields to represent the structural and conformational properties of nucleic acids in solution are compared. The force fields are AMBER 4.1, BMS, CHARMM22, and CHARMM27; the comparison of the latter two is the primary focus of this paper. The performance of each force field is evaluated first on its ability to reproduce the B-DNA decamer d(CGATTAATCG)(2) in solution with simulations in which the long-range electrostatics were treated by the particle mesh Ewald method; the crystal structure determined by Quintana et al. (1992) is used as the starting point for all simulations. A detailed analysis of the structural and solvation properties shows how well the different force fields can reproduce sequence-specific features. The results are compared with data from experimental and previous theoretical studies.

Surface electrostatic effects in oligonucleotide microarrays: control and optimization of binding thermodynamics

We present a theoretical thermodynamic framework for the design of more efficient oligonucleotide microarrays. A general thermodynamic relation is derived to describe the electrostatic surface effects on the binding of the assayed biomolecule to a surface-tethered molecular probe. The relation is applied to analyze how the nucleic acid target, the oligonuleotide probe, and their DNA duplex electrostatic interactions with the surface affect the hybridization on DNA arrays. Taking advantage of a closed form exact solution of the linear Poisson-Boltzmann equation for a charged ion-penetrable sphere in electrolyte solution interacting with a plane wall, we study the effects of the surface and solution conditions. Binding free energy is found as a function of the surface material, dielectric or metal, the surface charge density, linker molecule length, temperature, and added salt content. The charge or electric potential of the dielectric or metal surface, respectively, is shown to dominate the hybridization, especially at low added salt or short linker length. We predict that substantial enhancement of sensitivity, selectivity, and reliability of microarrays can be achieved by control of the surface conditions. As examples, we discuss how to overcome two limitations of current technologies: nonequal sensitivity of the probes with different GC and AT bases content, and poor match/mismatch discrimination. In addition, we suggest the design of microarray conditions where the tested nucleic acid is unfolded, thus making possible the screening of a larger sequence with single nucleotide resolution. These promising findings are discussed and further experimental tests suggested.

Theoretical methods for the simulation of nucleic acids

Different theoretical methods for the description of nucleic acid structures are reviewed. Firstly, we introduce the concept of classical force-field in the context of nucleic acid structures, discussing their accuracy. We then examine theoretical approaches to the description of nucleic acids based on: i) a rigid or quasi-rigid description of the molecule, ii) molecular mechanics optimization, and iii) molecular dynamics. Special emphasis is made ion current state of the art molecular dynamics simulations of nucleic acids structures.

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A<-->B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a 'slide first, roll later' mechanism for the B-->A helix transition. In the distribution of helical parameters in protein-DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.

Bias of purine stretches in sequenced chromosomes

We examined more than 700 DNA sequences (full length chromosomes and plasmids) for stretches of purines (R) or pyrimidines (Y) and alternating YR stretches; such regions will likely adopt structures which are different from the canonical B-form. Since one turn of the DNA helix is roughly 10 bp, we measured the fraction of each genome which contains purine (or pyrimidine) tracts of lengths of 10 bp or longer (hereafter referred to as 'purine tracts'), as well as stretches of alternating pyrimidines/purine (pyr/pur tracts') of the same length. Using this criteria, a random sequence would be expected to contain 1.0% of purine tracts and also 1.0% of the alternating pyr/pur tracts. In the vast majority of cases, there are more purine tracts than would be expected from a random sequence, with an average of 3.5%, significantly larger than the expectation value. The fraction of the chromosomes containing pyr/pur tracts was slightly less than expected, with an average of 0.8%. One of the most surprising findings is a clear difference in the length distributions of the regions studied between prokaryotes and eukaryotes. Whereas short-range correlations can explain the length distributions in prokaryotes, in eukaryotes there is an abundance of long stretches of purines or alternating purine/pyrimidine tracts, which cannot be explained in this way; these sequences are likely to play an important role in eukaryotic chromosome organisation.

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.

On the competition between water, sodium ions, and spermine in binding to DNA: a molecular dynamics computer simulation study

The interaction of DNA with the polyamine spermine(4+) (Spm(4+)), sodium ions, and water molecules has been studied using molecular dynamics computer simulations in a system modeling a DNA crystal. The simulation model consisted of three B-DNA decamers in a periodic hexagonal cell, containing 1200 water molecules, 8 Spm(4+), 32 Na(+), and 4 Cl(-) ions. The present paper gives a more detailed account of a recently published report of this system and compares results on this mixed Spm(4+)/Na(+)-cation system with an molecular dynamics simulation carried out for the same DNA decamer under similar conditions with only sodium counterions (Korolev et al., J. Mol. Biol. 308:907). The presence of Spm(4+) makes significant influence on the DNA hydration and on the interaction of the sodium ions with DNA. Spermine pushes water molecules out of the minor groove, whereas Na(+) attracts and organizes water around DNA. The major binding site of the Spm(4+) amino groups and the Na(+) ions is the phosphate group of DNA. The flexible polyamine spermine displays a high presence in the minor groove but does not form long-lived and structurally defined complexes. Sodium ions compete with Spm(4+) for binding to the DNA bases in the minor groove. Sodium ions also have several strong binding sites in the major groove. The ability of water molecules, Spm(4+), and Na(+) to modulate the local structure of the DNA double helix is discussed.

Solvation and hydration of proteins and nucleic acids: a theoretical view of simulation and experiment

Many theoretical, computational, and experimental techniques recently have been successfully used for description of the solvent distribution around macromolecules. In this Account, we consider recent developments in the areas of protein and nucleic acid solvation and hydration as seen by experiment, theory, and simulations. We find that in most cases not only the general phenomena of solvation but even local hydration patterns are more accurately discussed in the context of water distributions rather than individual molecules of water. While a few localized or high-residency waters are often associated with macromolecules in solution (or crystals from aqueous liquors), these are readily and accurately included in this more general description. The goal of this Account is to review the theoretical models used for the description of the interfacial solvent structure on the border near DNA and protein molecules. In particular, we hope to highlight the progress in this field over the past five years with a focus on comparison of simulation and experimental results.

RNA solvation: a molecular dynamics simulation perspective

With the availability of accurate methods to treat the electrostatic long-range interactions, molecular dynamics simulations have resulted in refined dynamical models of the structure of the hydration shell around RNA motifs. The models reviewed here range from basic Watson-Crick to more specific noncanonical base pairs, from "simple" double helices to RNA molecules displaying more complex tertiary folds, and from DNA/RNA hybrid double helices to RNA hybrids formed with a chemically modified strand.

Steered molecular dynamics investigations of protein function

Molecular recognition and mechanical properties of proteins govern molecular processes in the cell that can cause disease and can be targeted for drug design. Single molecule measurement techniques have greatly advanced knowledge but cannot resolve enough detail to be interpreted in terms of protein structure. We seek to complement the observations through so-called Steered Molecular Dynamics (SMD) simulations that link directly to experiments and provide atomic-level descriptions of the underlying events. Such a research program has been initiated in our group and has involved, for example, studies of elastic properties of immunoglobulin and fibronectin domains as well as the binding of biotin and avidin. In this article we explain the SMD method and suggest how it can be applied to the function of three systems that are the focus of modern molecular biology research: force transduction by the muscle protein titin and extracellular matrix protein fibronectin, recognition of antibody-antigene pairs, and ion selective conductivity of the K+ channel.

Sequence-dependent B<-->A transition in DNA evaluated with dimeric and trimeric scales

Experimental data on the sequence-dependent B<-->A conformational transition in 24 oligo- and polymeric duplexes yield optimal dimeric and trimeric scales for this transition. The 10 sequence dimers and the 32 trimers of the DNA duplex were characterized by the free energy differences between the B and A forms in water solution. In general, the trimeric scale describes the sequence-dependent DNA conformational propensities more accurately than the dimeric scale, which is likely related to the trimeric model accounting for the two interfaces between adjacent base pairs on both sides (rather than only one interface in the dimeric model). The exceptional preference of the B form for the AA:TT dimers and AAN:N'TT trimers is consistent with the cooperative interactions in both grooves. In the minor groove, this is the hydration spine that stabilizes adenine runs in B form. In the major groove, these are hydrophobic interactions between the thymine methyls and the sugar methylene groups from the preceding nucleotides, occurring in B form. This interpretation is in accord with the key role played by hydration in the B<-->A transition in DNA. Importantly, our trimeric scale is consistent with the relative occurrences of the DNA trimers in A form in protein-DNA cocrystals. Thus, we suggest that the B/A scales developed here can be used for analyzing genome sequences in search for A-philic motifs, putatively operative in the protein-DNA recognition.

Molecular dynamics simulations of B '-DNA: sequence effects on A-tract-induced bending and flexibility

Molecular dynamics (MD) simulations including water and counterions are reported on five examples of A-tract DNA oligonucleotide dodecamer duplexes for which crystal structures are available, the homopolymeric duplex sequences poly(dA) and poly(dG), and two related sequences that serve as controls. MD was performed using the AMBER suite of programs for 3 ns on each sequence. These results, combined with previously reported MDs on 25-mer and 30-mer oligonucleotides on sequences with phased A-tracts carried out under a similar simulation protocol, are used to examine salient issues in the structural chemistry of ApA steps and A-tract induced axis bending. MD modeling successfully describes the distinctive B' structure of A-tracts in solution as essentially straight (wedge angles of <1 degrees ), more rigid than generic B-form DNA, with slight base-pair inclination, high propeller twist and a minor groove narrowing 5' to 3'. The MD structures in solution agree closely with corresponding crystal structures, supporting the idea that crystal structures provide a good model for A-tract DNA structure in solution. From the collective MD results, bending and flexibility are calculated by step. Pyrimidine-purine steps are predicted to be most intrinsically bent and also most bendable, i.e. susceptible to bending. Pyrimidine-pyrimidine ( approximately purine-purine) and purine-pyrimidine steps show less intrinsic deformation and deformability. The MD calculated flexibility correlates well with the protein-induced bendability derived independently from the protein DNA crystal structures. The MD results indicate that bending and flexibility of base-pair steps in DNA are highly correlated, i.e. steps that exhibit the most intrinsic deformation from B-form DNA turn are also the most dynamically deformable. The MD description of A-tract-induced axis bending shows most consistency with the non A-tract, general-sequence model, in which the sequence curvature originates primarily in base-pair roll towards the major groove in non-A-tract regions of the sequence, particularly pyrimidine-purine steps. The direction of curvature is towards the minor groove viewed from opposite the A-tracts, but the A-tracts per se exhibit only minor deformation. The MD results are found to be consistent with the directionality of bending inferred for DNA sequences from gel retardation and cyclization experiments.

The 1.19 A X-ray structure of 2'-O-Me(CGCGCG)(2) duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove

The crystal and molecular structure of 2'-O-Me(CGCGCG)(2) has been determined at 1.19 A resolution, at 100 K, using synchrotron radiation. The structure in space group P3(2)12 is a half-turn right-handed helix that includes two 2-methyl-2,4-pentanediol (MPD) molecules bound in the minor groove. The structure deviates from A-form RNA. The duplex is overwound with an average value of 9.7 bp per turn, characterised as having a C3'-endo sugar pucker, very low base pair rise and high helical twist and inclination angles. The structure includes 65 ordered water molecules. Only a single row of water molecules is observed in the minor groove due to the presence of hydrophobic 2'-O-methyl groups. As many as five magnesium ions are located in the structure. Two are in the major groove and interact with O(6) and N(7) of guanosine and N(4) of cytidine residues through their hydration spheres. This work provides the first example of molecular interactions of nucleic acids with MPD, which was used as a precipitant, cryo-solvent and resolution enhancing agent. The two MPD molecules intrude into the hydration network in the minor groove, each forming hydrogen bonds between their secondary hydroxyl group and exo-amino functions of guanosine residues. Comparison of the 2'-O-Me(CGCGCG)(2) structure in the P3(2)12 and P6(1)22 crystals delineates stability of the water network within the minor groove to dehydration by MPD and is of interest for evaluating factors governing small molecule binding to RNA. Intrusion of MPD into the minor groove of 2'-O-Me(CGCGCG)(2) is discussed with respect to RNA dehydration, a prerequisite of Z-RNA formation.

Volumetric characterization of the hydration properties of heterocyclic bases and nucleosides

We have determined the partial molar volumes, expansibilities, and adiabatic compressibilities of six heterocyclic nucleic acid bases, five ribonucleosides, and six 2'-deoxyribonucleosides within the temperature range 18-55 degrees C. We interpret the resulting data in terms of the hydration of the component hydrophobic and polar atomic groups. From our temperature-dependent volumetric studies, we found that the total contraction of water caused by polar groups of each individual heterocyclic base and nucleoside depends on the proximity and chemical nature of other functional groups of the solute. In addition, the compressibility contributions of polar groups vary greatly in sign and magnitude depending on the surrounding functional groups. In agreement with previous studies, our results are suggestive of little or no interaction between the sugar and base moieties of a nucleoside. In general, our data shed light into the hydration properties of individual heterocyclic bases and nucleosides, which may have significant implications for the sequence-dependent hydration of nucleic acids. We discuss the potential importance of our results for developing an understanding of the role that solvent plays in the stabilization/destabilization of nucleic acid structures.

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.

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.

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.

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA<RNA (+3 H(2)O) and B-DNA<RNA (+2 H(2)O). The calculated residence times of water molecules in the first hydration shell of the helices range from 0.5 to 1 ns, in good agreement with available experimental data. Such water molecules are essentially found in the vicinity of the phosphate groups and in the DNA minor groove. The calculated number of ions that break into the first hydration shell of the nucleic acids is close to 0.5 per base-pair for both RNA and DNA. These ions form contacts essentially with the oxygen atoms of the phosphate groups and with the guanine N7 and O6 atoms; they display residence times in the deep/major groove approaching 500 ps. Further, a significant sequence-dependent effect on ion binding has been noted. Despite slight structural differences, K(+) binds essentially to GpC and not to CpG steps. These results may be of importance for understanding various sequence-dependent binding affinities. Additionally, the data help to rationalize the experimentally observed differences in gel electrophoretic mobility between RNA and DNA as due to the difference in hydration (two water molecules in favor of RNA) rather than to strong ion-binding features, which are largely similar for both nucleic acid structures.

A new crystal form for the dodecamer C-G-C-G-A-A-T-T-C-G-C-G: symmetry effects on sequence-dependent DNA structure

The dodecanucleotide d(CGCGAATTCGCG) has been crystallised in the space group P3(2)12, representing a new crystal form for this sequence. The structure has been solved by molecular replacement and refined at 1.8 A resolution. The present structure contrasts with previous ones for this sequence since it is situated on a crystallographic 2-fold axis, and the crystal symmetry reflects the palindromic nature of this sequence. Some features accord with previous observations, notably that the minor groove is hydrated with a continuous spine of solvent. There is no evidence of alkali metal ions within this spine. The minor groove retains its narrow width, although it is now symmetric and extends over the A/T tract. Various base and base-pair morphological parameters have been examined. Their values do not show significant correlations with earlier reports, suggesting that crystal packing effects on them are more dominant than has been hitherto realised.

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.

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.

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.

Sequence-specific binding of counterions to B-DNA

Recent studies by x-ray crystallography, NMR, and molecular simulations have suggested that monovalent counterions can penetrate deeply into the minor groove of B form DNA. Such groove-bound ions potentially could play an important role in AT-tract bending and groove narrowing, thereby modulating DNA function in vivo. To address this issue, we report here (23)Na magnetic relaxation dispersion measurements on oligonucleotides, including difference experiments with the groove-binding drug netropsin. The exquisite sensitivity of this method to ions in long-lived and intimate association with DNA allows us to detect sequence-specific sodium ion binding in the minor groove AT tract of three B-DNA dodecamers. The sodium ion occupancy is only a few percent, however, and therefore is not likely to contribute importantly to the ensemble of B-DNA structures. We also report results of ion competition experiments, indicating that potassium, rubidium, and cesium ions bind to the minor groove with similarly weak affinity as sodium ions, whereas ammonium ion binding is somewhat stronger. The present findings are discussed in the light of previous NMR and diffraction studies of sequence-specific counterion binding to DNA.

Sodium and chlorine ions as part of the DNA solvation shell

The distribution of sodium and chlorine ions around DNA is presented from two molecular dynamics simulations of the DNA fragment d(C(5)T(5)). (A(5)G(5)) in explicit solvent with 0.8 M additional NaCl salt. One simulation was carried out for 10 ns with the CHARMM force field that keeps the DNA structure close to A-DNA, the other for 12 ns with the AMBER force field that preferentially stabilizes B-DNA conformations (, Biophys. J. 75:134-149). From radial distributions of sodium and chlorine ions a primary ion shell is defined. The ion counts and residence times of ions within this shell are compared between conformations and with experiment. Ordered sodium ion sites were found in minor and major grooves around both A and B-DNA conformations. Changes in the surrounding hydration structure are analyzed and implications for the stabilization of A-DNA and B-DNA conformations are discussed.

Probing the role of structural water in a duplex oligodeoxyribonucleotide containing a water-mimicking base analog

Molecular dynamics simulations were performed on models of the dodecamer DNA double-stranded segment, [d(CGCGAATTCGCG)](2), in which each of the adenine residues, individually or jointly, was replaced by the water-mimicking analog 2'-deoxy-7-(hydroxy-methyl)-7-deazaadenosine (hm(7)c(7)dA) [Rockhill, J.K., Wilson,S.R. and Gumport,R.I. (1996) J. Am. Chem. Soc.,118, 10065-10068]. The simulations, when compared with those of the dodecamer itself, show that incorporation of the analog affects neither the overall DNA structure nor its hydrogen-bonding and stacking interactions when it replaces a single individual base. Furthermore, the water molecules near the bases in the singly-substituted oligonucleotides are similarly unaffected. Double substitutions lead to differences in all the aforementioned parameters with respect to the reference sequence. The results suggest that the analog provides a good mimic of specific 'ordered' water molecules observed in contact with DNA itself and at the interface between protein and DNA in specific complexes.