Sodium and chlorine ions as part of the DNA solvation shell

Counterion Condensation on Charged Spheres, Cylinders, and Planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.

Quantitative analysis of monovalent counterion binding to random-sequence, double-stranded DNA using the replacement ion method

A variation of affinity capillary electrophoresis, called the replacement ion (RI) method, has been developed to measure the binding of monovalent cations to random sequence, double-stranded (ds) DNA. In this method, the ionic strength is kept constant by gradually replacing a non-binding ion in the solution with a binding ion and measuring the mobility of binding and non-binding analytes as a function of binding ion concentration. The method was validated by measuring the binding of Li+ ions to adenosine nucleotides; the apparent dissociation constants obtained by the RI method are comparable to literature values obtained by other methods. The binding of Tris+, NH4+, Li+, Na+, and K+ to dsDNA was then investigated. The apparent dissociation constants observed for counterion binding to a random-sequence 26-base pair (bp) oligomer ranged from 71 mM for Tris+ to 173 mM for Na+ and K+. Hence, positively charged Tris buffer ions will compete with other monovalent cations in Tris-buffered solutions. The bound cations identified in this study may correspond to the strongly correlated, tightly bound ions recently postulated to exist as a class of ions near the surface of dsDNA (Tan, Z.-J., and Chen, S.-J. (2006) Biophys. J. 91, 518-536). Monovalent cation binding to random-sequence dsDNA would be expected to occur in addition to any site-specific binding of cations to A-tracts or other DNA sequence motifs. Single-stranded DNA oligomers do not bind the five tested cations under the conditions investigated here.

Electrostatic, Steric, and Hydration Interactions Favor Na+ Condensation around DNA Compared with K+

Condensation of monovalent counterions around DNA influences polymer properties of the DNA chain. For example, the Na(+) ions show markedly stronger propensity to induce multiple DNA chains to assemble into compact structures compared with the K(+) ions. To investigate the similarities and differences in the sodium and potassium ion condensation around DNA, we carried out a number of extensive all-atom molecular dynamics simulations of a DNA oligomer consisting of 16 base pairs, [d(CGAGGTTTAAACCTCG)](2), in explicit water. We found that the Na(+) ions penetrate the DNA interior and condense around the DNA exterior to a significantly larger degree compared with the K(+) ions. We have provided a microscopic explanation for the larger Na(+) affinity toward DNA that is based on a combination of steric, electrostatic, and hydration effects. Unexpectedly, we found that the Cl(-) co-ions provide more efficient electrostatic screening for the K(+) ions than for the Na(+) ions, contributing to the larger Na(+) condensation around DNA. To examine the importance of the discrete nature of water and ions, we also computed the counterion distributions from the mean-field electrostatic theory, demonstrating significant disagreements with the all-atom simulations. Prior experimental results on the relative extent of the Na(+) and K(+) condensation around DNA were somewhat contradictory. Recent DNA compaction experiments may be interpreted to suggest stronger Na(+) condensation around DNA compared to K(+), which is consistent with our simulations. We also provide a simple interpretation for the experimentally observed increase in DNA electrophoretic mobility in the alkali metal series, Li(+) < Na(+) < K(+) < Rb(+). We compare the DNA segment conformational preferences in various buffers with the proposed NMR models.

Competition between cation-π interactions and intermolecular hydrogen bonds in alkali metal ion-phenol clusters. II. Phenol trimer

The competition between ion-molecule and molecule-molecule interactions was investigated in M+(phenol)3 cluster ions for M=Li, Na, K, and Cs. Infrared predissociation spectroscopy in the O-H stretch region was used to characterize the structure of the cluster ions. By adjusting the experimental conditions, it was possible to generate species where argon was additionally bound in order to investigate cold cluster ions. From a comparison of the M+(phenol)3 spectra with the M+(phenol)3Ar spectra, it is clear that the relative populations of hydrogen-bonded configurations are significantly higher in the colder (argon-bearing) species. For the cold species, the IR spectra were compared with minimum energy ab initio calculations to elucidate the hydrogen-bonded structures. The experimental spectra are most consistent with a cyclic hydrogen-bonded configuration for Cs+(phenol)3 in which the ion binds to the phenol molecules via cation-pi interactions, and noncyclic configurations for Li+, Na+, and K+.

Electrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte

A statistical mechanical model is presented which explicitly accounts for the fluctuations, the electrostatic, and the excluded volume correlations for ions bound to a polyelectrolyte such as DNA. The method can be employed to treat a wide range of ionic conditions including multivalent ions. The microscopic framework of the theory permits the use of realistic finite length and grooved structural model for the polyelectrolyte and modeling of the finite size of the bound ions. Test against Monte Carlo simulations suggests that the theory can give accurate predictions for the ion distribution and the thermodynamic properties. For multivalent ions, the theory makes improved predictions as compared with the mean-field approach. Moreover, for long polyelectrolyte and dilute salt concentration, the theory predicts ion binding properties that agree with the counterion condensation theory.

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

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

Effects of counterion valency on the damping of phonons propagating along the axial direction of liquid-crystalline DNA

The phonon propagation and damping along the axial direction of films of aligned 40 wt % calf-thymus DNA rods are studied by inelastic x-ray scattering (IXS). The IXS spectra are analyzed with the generalized three effective eigenmode theory, from which we extract the dynamic structure factor S(Q,E) as a function of transferred energy E=variant Planck's over 2piomega, and the magnitude of the transferred wave vector Q. S(Q,E) of a DNA sample typically consists of three peaks, one central Rayleigh scattering peak, and two symmetric Stokes and anti-Stokes Brillouin side peaks. By analyzing the Brillouin peaks, the phonon excitation energy and damping can be extracted at different Q values from about 4 to 30 nm(-1). A high-frequency sound speed is obtained from the initial slope of the linear portion of the dispersion relation below Q=4 nm(-1). The high-frequency sound speed obtained in this Q range is 3100 ms, which is about twice faster than the ultrasound speed of 1800 ms, measured by Brillouin light scattering at Q approximately 0.01 nm(-1) at the similar hydration level. Our observations provide further evidence of the strong coupling between the internal dynamics of a DNA molecule and the dynamics of the solvent. The effect on damping and propagation of phonons along the axial direction of DNA rods due to divalent and trivalent counterions has been studied. It is found that the added multivalent counterions introduce stronger phonon damping. The phonons at the range between approximately 12.5 and approximately 22.5 nm(-1) are overdamped by the added counterions according to our model analyses. The intermediate scattering function is extracted and it shows a clear two-step relaxation with the fast relaxation time ranging from 0.1 to 4 ps.

The effect of salt concentration on DNA conformation transition: a molecular-dynamics study

We performed three 3-ns molecular dynamics simulations of d(CGCGAATTCGCG)2 using the AMBER 8 package to determine the effect of salt concentration on DNA conformational transitions. All the simulations were started with A-DNA, with different salt concentrations, and converged with B-DNA with similar conformational parameters. However, the dynamic processes of the three MD simulations were very different. We found that the conformation transition was slow in the solution with higher salt concentration. To determine the cause of this retardation, we performed three additional 1.5-ns simulations starting with B-DNA and with the salt concentrations corresponding to the simulations mentioned above. However, astonishingly, there was no delayed conformation evolution found in any of the three simulations. Thus, our simulation conclusion is that higher salt concentrations slows the A --> B conformation transition, but have no effect on the final stable structure. [Figure: see text].

Differential Solvation and Tautomer Stability of a Model Base Pair within the Minor and Major Grooves of DNA

2-(2'-Hydroxyphenyl)benzoxazole (HBO) may be used as a model base pair to study solvation, duplex environment, and tautomerization within the major and minor groves of DNA duplexes. In its ground state, HBO possesses an enol moiety which may be oriented syn or anti relative to the imino nitrogen of the benzoxazole ring. In the absence of external hydrogen-bond donors and acceptors HBO exists as the internally hydrogen-bonded syn-enol, a mimic of the rare base pair tautomer found in DNA, which may be photoinduced to tautomerize and form the keto tautomer, a mimic of the dominant base pair tautomer. Previously, we demonstrated that when incorporated into DNA such that the enol moiety is positioned in the major groove, HBO is not solvated, exists exclusively as the internally hydrogen-bonded syn-enol which is efficiently photoinduced to tautomerize, and the corresponding keto tautomer is preferentially stabilized. In stark contrast, we now show that when HBO is incorporated in DNA such that the enol moiety is positioned in the minor groove, the enol tautomer is preferentially stabilized. Molecular dynamics simulations suggest that this results from the formation of a stable hydrogen-bond between the HBO enol and the O4' atom of an adjacent nucleotide, an H-bond acceptor that is only available in the minor groove. The differential stabilization of the enol and keto tautomers in the major and minor grooves may reflect the functions for which these environments evolved, including duplex replication, stability, and recognition.

Competition between cation-π interactions and intermolecular hydrogen bonds in alkali metal ion-phenol clusters. I. Phenol dimer

The competition between ion-molecule and molecule-molecule interactions was investigated in M+(phenol)2 cluster ions for M=Li, Na, K, and Cs. Infrared predissociation spectroscopy in the O-H stretch region was used to characterize the structure of the cluster ions. By adjusting the experimental conditions, it was possible to generate species where argon was additionally bound in order to investigate cold cluster ions. The spectra showed the presence of hydrogen bonding in the colder M+(phenol)2Ar cluster ions but the absence of hydrogen bonding in the warmer M+(phenol)2 species. For the cold species, the IR spectra were compared with minimum-energy ab initio calculations to elucidate the hydrogen-bonded structures. In the dominant hydrogen-bonded configurations observed experimentally, the phenol molecules form hydrogen-bonded dimers and the alkali-metal ions bind to the phenol via a cation-pi interaction with the aromatic ring. Increasing the strength of the cation-pi interaction by decreasing the ion size forces the distance between the phenol O-H groups to increase, thus weakening the intermolecular hydrogen bond. Free-energy differences of different configurations relative to the ground state demonstrate that hydrogen-bonded structures are enthalpically favored, while non-hydrogen-bonded structures are entropically favored and are thus observed in the warm cluster ions.

Ions and RNA Folding

The problem of how ions influence the folding of RNA into specific tertiary structures is being addressed from both thermodynamic (by how much do different salts affect the free energy change of folding) and structural (how are ions arranged on or near an RNA and what kinds of environments do they occupy) points of view. The challenge is to link these different approaches in a theoretical framework that relates the energetics of ion-RNA interactions to the spatial distribution of ions. This review distinguishes three different kinds of ion environments that differ in the extent of direct ion-RNA contacts and the degree to which the ion hydration is perturbed, and summarizes the current understanding of the way each environment relates to the overall energetics of RNA folding.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ∼30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

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

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

Electrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte

A statistical mechanical model is presented which explicitly accounts for the fluctuations, the electrostatic, and the excluded volume correlations for ions bound to a polyelectrolyte such as DNA. The method can be employed to treat a wide range of ionic conditions including multivalent ions. The microscopic framework of the theory permits the use of realistic finite length and grooved structural model for the polyelectrolyte and modeling of the finite size of the bound ions. Test against Monte Carlo simulations suggests that the theory can give accurate predictions for the ion distribution and the thermodynamic properties. For multivalent ions, the theory makes improved predictions as compared with the mean-field approach. Moreover, for long polyelectrolyte and dilute salt concentration, the theory predicts ion binding properties that agree with the counterion condensation theory.

Anion binding to nucleic acids

Nucleic acids are generally considered as efficient cation binders. Therefore, the likelihood that negatively charged ions might intrude their first hydration shell is rarely considered. Here, we show on the basis of (i) a survey of the Nucleic Acid Database, (ii) several structures extracted from the Cambridge Structural Database, and (iii) molecular dynamics simulations, that the nucleotide electropositive edges involving mainly amino, imino, and hydroxyl groups can cast specific anion binding sites. These binding sites constitute also good locations for the binding of the negatively charged groups of the Asp and Glu residues or the nucleic acid phosphate groups. Furthermore, it is observed in several instances that anions, like water molecules and cations, do mediate protein/nucleic acid interactions. Thus, anions as well as negatively charged groups are directly involved in specific recognition and folding phenomena involving polyanionic nucleic acids.

Ion motions in molecular dynamics simulations on DNA

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

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.

DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na+ and K+, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na+ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K+ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive alpha/gamma transitions in the phosphate backbone during the simulation with Na+. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.

Attraction between DNA molecules mediated by multivalent ions

The effective force between two parallel DNA molecules is calculated as a function of their mutual separation for different valencies of counterion and salt ions and different salt concentrations. Computer simulations of the primitive model are used and the shape of the DNA molecules is accurately modeled using different geometrical shapes. We find that multivalent ions induce a significant attraction between the DNA molecules whose strength can be tuned by the averaged valency of the ions. The physical origin of the attraction is traced back either to electrostatics or to entropic contributions. For multivalent counterions and monovalent salt ions, we find a salt-enhanced repulsion effect: the force is first attractive but gets repulsive with increasing salt concentration. Furthermore, we show that the multivalent-ion-induced attraction does not necessarily correlate with DNA overcharging.

Brownian dynamics simulations of ion atmospheres around polyalanine and B-DNA: effects of biomolecular dielectric

We have extended an earlier Brownian dynamics simulation algorithm for simulating the structural dynamics of ions around biomolecules to accommodate dielectric inhomogeneity. The electrostatic environment of a biomolecule immersed in water was obtained by numerically solving the Poisson equation with the biomolecule treated as a low dielectric region and the solvent treated as a high dielectric region. Instead of using the mean-field type approximations of ion interactions as in the Poisson-Boltzmann model, the ions were treated explicitly by allowing them to evolve dynamically under the electrostatic field of the biomolecule. This model thus accounts for ion-ion correlations and the finite-size effects of the ions. For a 13-residue alpha-helical polyalanine and a 12-base-pair bp B-form DNA, we found that the choice of the dielectric constant of the biomolecule has much larger effects on the mean ionic structure around the biomolecule than on the fluctuational and dynamical properties of the ions surrounding the biomolecule.

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.

Symmetric K+ and Mg2+ Ion-binding Sites in the 5S rRNA Loop E Inferred from Molecular Dynamics Simulations

Potassium binding to the 5S rRNA loop E motif has been studied by molecular dynamics at high (1.0 M) and low (0.2 M) concentration of added KCl in the presence and absence of Mg2+. A clear pattern of seven deep groove K+ binding sites or regions, in all cases connected with guanine N7/O6 atoms belonging to GpG, GpA, and GpU steps, was identified, indicating that the LE deep groove is significantly more ionophilic than the equivalent groove of regular RNA duplexes. Among all, two symmetry-related sites (with respect to the central G.A pair) were found to accommodate K+ ions with particularly long residence times. In a preceding molecular dynamics study by Auffinger et al. in the year 2003, these two sites were described as constituting important Mg2+ binding locations. Altogether, the data suggest that these symmetric sites correspond to the loop E main ion binding regions. Indeed, they are located in the deep groove of an important ribosomal protein binding motif associated with a fragile pattern of non-Watson-Crick pairs that has certainly to be stabilized by specific Mg2+ ions in order to be efficiently recognized by the protein. Besides, the other sites accommodate monovalent ions in a more diffuse way pointing out their lesser significance for the structure and function of this motif. Ion binding to the shallow groove and backbone atoms was generally found to be of minor importance since, at the low concentration, no well defined binding site could be characterized while high K+ concentration promoted mostly unspecific potassium binding to the RNA backbone. In addition, several K+ binding sites were located in positions equivalent to water molecules from the first hydration shell of divalent ions in simulations performed with magnesium, indicating that ion binding regions are able to accommodate both mono- and divalent ionic species. Overall, the simulations provide a more precise but, at the same time, a more intricate view of the relations of this motif with its ionic surrounding.

Adsorption of monovalent and multivalent cations and anions on DNA molecules

Adsorption of monovalent and multivalent cations and anions on a deoxyribose nucleic acid (DNA) molecule from a salt solution is investigated by computer simulation. The ions are modeled as charged hard spheres, the DNA molecule as a point charge pattern following the double-helical phosphate strands. The geometrical shape of the DNA molecules is modeled on different levels ranging from a simple cylindrical shape to structured models which include the major and minor grooves between the phosphate strands. The densities of the ions adsorbed on the phosphate strands in the major and in the minor grooves are calculated. First, we find that the adsorption pattern on the DNA surface depends strongly on its geometrical shape: counterions adsorb preferentially along the phosphate strands for a cylindrical model shape, but in the minor groove for a geometrically structured model. Second, we find that an addition of monovalent salt ions results in an increase of the charge density in the minor groove while the total charge density of ions adsorbed in the major groove stays unchanged. The adsorbed ion densities are highly structured along the minor groove while they are almost smeared along the major groove. Furthermore, for a fixed amount of added salt, the major-groove cationic charge is independent of the counterion valency. For increasing salt concentration the major groove is neutralized while the total charge adsorbed in the minor groove is constant. DNA overcharging is detected for multivalent salts. Simulations for larger ion radii, which mimic the effect of ion hydration, indicate an increased adsorbtion of cations in the major groove.

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.

The Mg2+ binding sites of the 5S rRNA loop E motif as investigated by molecular dynamics simulations

Molecular dynamics simulations have been used to investigate the binding of Mg(2+) ions to the deep groove of the eubacterial 5S rRNA loop E. The simulations suggest that long-lived and specific water-mediated interactions established between the hydrated ions and the RNA atoms lining up the binding sites contribute to the stabilization of this motif. The Mg(2+) binding specificity is modulated by two factors: (i) a required electrostatic complementarity and (ii) a structural correspondence between the hydrated ion and its binding pocket that can be estimated by its degree of dehydration and the resulting number and lifetime of the intervening water-mediated contacts. Two distinct binding modes for pentahydrated Mg(2+) ions that result in a significant freezing of the tumbling motions of the ions are described, and mechanistic details related to the stabilization of nucleic acids by divalent ions are provided.

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.

Comments on selected aspects of nucleic acid electrostatics

Recent experimental, theoretical, and computational developments in the field of nucleic acid electrostatics have brought interesting concepts to the fore. The phosphate charge on the double helix apparently influences its structure. When the charge is neutralized asymmetrically, the resulting force imbalance drives bending toward the neutralized side. When the charge is uniformly neutralized, the force imbalance acts to buckle the helical axis, resulting in a compact tertiary conformation. Sharing of condensed counterions by single strands is a stabilizing factor for formation of the double helix. Sharing of condensed counterions by two double helices causes clustering of DNA and may be a factor in RNA folding. Support for these statements is reviewed.

The role of minor groove functional groups in DNA hydration

Here we describe the crystal structure of modified [d(CGCGAATTCGCG)]2 refined to 2.04 A. The modification, which affects only the two thymines at the central ApT step, involves isosteric removal of the 2-keto oxygen atoms and substitution of the N1 nitrogen with carbon. The crystal structure reveals the ability of this modified thymine to effectively base pair with adenine in [d(CGCGAAtTCGCG)]2. The structure also suggests that the minor groove 'spine of hydration' is destabilized but essentially intact.

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.

DNA dynamics in a water drop without counterions

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

Electrostatic free energy of the DNA double helix in counterion condensation theory

Polyelectrolyte theory based on counterion condensation is extended from the standard line charge model to helical and double helical charge arrays. The number of condensed counterions turns out to be the same as for a line charge with charge density equal to the axial charge density of the helix. Also, the logarithmic salt dependence of the electrostatic free energy is the same in the range of lower salt concentration, so that the limiting laws remain unchanged. However, the internal free energy of the condensed layer of counterions and the overall electrostatic free energy depend on the helical parameters. At higher salt, the free energies of both single and double helix are negative, indicating electrostatic stabilization of the helical charge lattices due to the mixing entropy of the condensed counterions. Except at very low salt, the free energy of a single helix is higher than the free energy of a double helix with twice the charge density. With B-DNA parameters and single strands modeled as single helices, the predicted salt dependence of the free energy of transition from double helix to separated single strands has a maximum at approximately 0.2 M salt, close to the location in the laboratory of this well-known feature of the DNA strand separation transition. We also calculate the electrostatic free energy for the transition of the DNA double helix from the B to the A conformation. The B form is electrostatically stable over most of the salt range, but there is a spontaneous electrostatic transition to A near 1 M salt. The electrostatic free energy values are close to the experimental values of the overall (electrostatic plus non-electrostatic) transition free energies for A-philic base pair sequences. We are led to suggest that the experimentally observed B-to-A transition for A-philic sequences near 1 M salt in water is governed by the polyelectrolyte properties of these two conformations of the DNA double helix. The effect of ethanol, however, cannot be attributed to lowering of the bulk dielectric constant.

Solvent interactions of halophilic malate dehydrogenase

Malate dehydrogenase from the extreme halophilic Haloarcula marismortui (Hm MalDH) is an acidic protein that is unstable below molar salt concentrations. The solvated folded protein was studied by small-angle neutron scattering in solvents containing salt: NaCl, NaCH(3)CO(2), KF, NH(4)Cl, NH(4)CH(3)CO(2), (NH(4))(2)SO(4), MgCl(2), and MgSO(4). It was found that the global solvent interactions depend mainly on the nature of the cation. Complementary mass density measurements in MgCl(2), NaCl, NaCH(3)CO(2), and (NH(4))(2)SO(4) allowed determining the partial molal volumes of the protein, which were found to increase slightly with the salt, and the preferential salt binding parameters for each solvent condition. These are strongly dependent on the cation type and salt concentration. Hm MalDH can be modeled as an invariant particle binding 4100 water molecules in MgCl(2) and 2000 +/- 200 in NaCl, NaCH(3)CO(2), or (NH(4))(2)SO(4). The number of salt molecules associated to the particle decreases from about 85 to 0 in the order MgCl(2) > NaCl = NaCH(3)CO(2) > (NH(4))(2)SO(4). Alternatively, we considered exchangeable sites for water and salt with the effects of solvent nonideality. It does not change the description of the solvent interactions. Solvent anions act on Hm MalDH stability through a limited number of strong binding sites, as those seen at the interfaces of Hm MalDH by crystallography. Cations would act through some strong and numerous weak binding sites defined on the folded protein, in possible addition to nonspecific hydration effects.

Influence of netropsin's charges on the minor groove width of d(CGCGAATTCGCG)2

The exact understanding of the interaction of minor groove binding drugs with DNA is of interest due to their importance as transcription controlling drugs. In this study we performed four molecular dynamics simulations, one of the uncomplexed d(CGCGAATTCGCG)(2) dodecamer and three simulations of the DNA complexed with the minor groove binder netropsin. The charged guanidinium and amidinium ends of the small ligand were in one simulation formally uncharged, in the second one normally charged, and in the third simulation we doubled the charges of the two ends. So we are able to filter out the influence the charges exert on the DNA structure. The positive charges reduce the width of the minor groove showing that charges are able to modify the groove width by charge neutralization of the negative phosphate groups. The quality of the used force field was successfully tested by comparing the results of the uncomplexed dodecamer with already reported NMR and x-ray studies. Thus our simulations should be able to describe the minor groove width of DNA in a correct manner underlying the validity of the results.

Characterization of divalent cation localization in the minor groove of the A(n)T(n) and T(n)A(n) DNA sequence elements by (1)H NMR spectroscopy and manganese(II)

The localization of Mn(2+) in A-tract DNA has been studied by (1)H NMR spectroscopy using a series of self-complementary dodecamer oligonucleotides that contain the sequence motifs A(n)(n) and T(n)A(n), where n = 2, 3, or 4. Mn(2+) localization in the minor groove is observed for all the sequences that have been studied, with the position and degree of localization being highly sequence-dependent. The site most favored for Mn(2+) localization in the minor groove is near the 5'-most ApA step for both the T(n)A(n) and the A(n)T(n) series. For the T(n)A(n) series, this results in two closely spaced symmetry-related Mn(2+) localization sites near the center of each duplex, while for the A(n)T(n) series, the two symmetry-related sites are separated by as much as one half-helical turn. The degree of Mn(2+) localization in the minor groove of the T(n)A(n) series decreases substantially as the AT sequence element is shortened from T(4)A(4) to T(2)A(2). The A(n)T(n) series also exhibits length-dependent Mn(2+) localization; however, the degree of minor groove occupancy by Mn(2+) is significantly less than that observed for the T(n)A(n) series. For both A(n)T(n) and T(n)A(n) sequences, the 3'-most AH2 resonance is the least broadened of the AH2 resonances. This is consistent with the observation that the minor groove of A-tract DNA narrows in the 5' to 3' direction, apparently becoming too narrow after two base pairs for the entry of a fully hydrated divalent cation. The results that are reported illustrate the delicate interplay that exists between DNA nucleotide sequence, minor groove width, and divalent cation localization. The proposed role of cation localization in helical axis bending by A-tracts is also discussed.

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.

Free energy and structural pathways of base flipping in a DNA GCGC containing sequence

Structural distortions of DNA are essential for its biological function due to the genetic information of DNA not being physically accessible in the duplex state. Base flipping is one of the simplest structural distortions of DNA and may represent an initial event in strand separation required to access the genetic code. Flipping is also utilized by DNA-modifying and repair enzymes to access specific bases. It is typically thought that base flipping (or base-pair opening) occurs via the major groove whereas minor groove flipping is only possible when mediated by DNA-binding proteins. Here, umbrella sampling with a novel center-of-mass pseudodihedral reaction coordinate was used to calculate the individual potentials of mean force (PMF) for flipping of the Watson-Crick (WC) paired C and G bases in the CCATGCGCTGAC DNA dodecamer. The novel reaction coordinate allowed explicit investigation of the complete flipping process via both the minor and major groove pathways. The minor and major groove barriers to flipping are similar for C base flipping while the major groove barrier is slightly lower for G base flipping. Minor groove flipping requires distortion of the WC partner while the flipping base pulls away from its partner during major groove flipping. The flipped states are represented by relatively flat free energy surfaces, with a small, local minimum observed for the flipped G base. Conserved patterns of phosphodiester backbone dihedral distortions during flipping indicate their essential role in the flipping process. During flipping, the target base tracks along the respective grooves, leading to hydrogen-bonding interactions with neighboring base-pairs. Such hydrogen-bonding interactions with the neighboring sequence suggest a novel mechanism of sequence dependence in DNA dynamics.

Bending of DNA by asymmetric charge neutralization: all-atom energy simulations

DNA dodecamers of the alternating d(CG).d(CG) sequence with six phosphate groups either charge-neutralized or substituted by neutral methylphosphonates across the major or minor groove have been subjected to energy minimization to determine the conformational effect of the asymmetric elimination of phosphate charge. We report bending angles, directions of bending, and detailed structural characteristics such as groove widths and local base-pair parameters. Our principal results are that charge neutralization on one face of the DNA induces significant bending toward the neutralized face, in agreement with theoretical predictions on a simplified model and experimental data on a similar base-pair sequence, and that the DNA conformation averaged over all stereospecific methylphosphonate substitutions is nearly the same as the conformation produced by charge neutralization of the phosphates. Individual isomers, however, cover a wide range of structures, with the magnitude and direction of overall bending sensitive to the precise stereochemical pattern of neutralization. Our simulation does not explicitly contain counterions, and the results therefore suggest that counterions can influence DNA structure by neutralizing the phosphate charge. These data provide new hints into the molecular mechanisms which underlie the deformations of DNA structure induced by the binding of positively charged proteins and other tightly associated cationic species.

Explicit solvent molecular dynamics simulation of duplex formed by the modified oligonucleotide with alternating phosphate/phosphonate internucleoside linkages and its natural counterpart

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT(11).dA(11) duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5'-O-PO(2)-CH(2)-O-3' (5'PC3') or 5'-O-CH(2)-PO(2)-O-3' (5'CP3') type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5'PC3' and ggg, tgg for the 5'CP3', respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.

Molecular dynamics simulation of nucleic acids: successes, limitations, and promise

In the last five years we have witnessed a significant increase in the number publications describing accurate and reliable all-atom molecular dynamics simulations of nucleic acids. This increase has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions, the availability of faster parallel computers, and the development of well-validated empirical molecular mechanical force fields. With these technologies, it has been demonstrated that simulation is not only capable of consistently reproducing experimental observations of sequence specific fine structure of DNA, but also can give detailed insight into prevalent problems in nucleic acid structure, ion association and specific hydration of nucleic acids, polyadenine tract bending, and the subtle environmental dependence of the A-DNA-B-DNA duplex equilibrium. Despite the advances, there are still issues with the methods that need to be resolved through rigorous controlled testing. In general, these relate to deficiencies of the underlying molecular mechanical potentials or applied methods (such as the imposition of true periodicity in Ewald simulations and the need for energy conservation), and significant limits in effective conformational sampling. In this perspective, we provide an overview of our experiences, provide some cautionary notes, and provide recommendations for further study in molecular dynamics simulation of nucleic acids.

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.