Sodium and chlorine ions as part of the DNA solvation shell

Accurate small and wide angle x-ray scattering profiles from atomic models of proteins and nucleic acids

A new method is introduced to compute X-ray solution scattering profiles from atomic models of macromolecules. The three-dimensional version of the Reference Interaction Site Model (RISM) from liquid-state statistical mechanics is employed to compute the solvent distribution around the solute, including both water and ions. X-ray scattering profiles are computed from this distribution together with the solute geometry. We describe an efficient procedure for performing this calculation employing a Lebedev grid for the angular averaging. The intensity profiles (which involve no adjustable parameters) match experiment and molecular dynamics simulations up to wide angle for two proteins (lysozyme and myoglobin) in water, as well as the small-angle profiles for a dozen biomolecules taken from the BioIsis.net database. The RISM model is especially well-suited for studies of nucleic acids in salt solution. Use of fiber-diffraction models for the structure of duplex DNA in solution yields close agreement with the observed scattering profiles in both the small and wide angle scattering (SAXS and WAXS) regimes. In addition, computed profiles of anomalous SAXS signals (for Rb(+) and Sr(2+)) emphasize the ionic contribution to scattering and are in reasonable agreement with experiment. In cases where an absolute calibration of the experimental data at q = 0 is available, one can extract a count of the excess number of waters and ions; computed values depend on the closure that is assumed in the solution of the Ornstein-Zernike equations, with results from the Kovalenko-Hirata closure being closest to experiment for the cases studied here.

Free energetics of carbon nanotube association in aqueous inorganic NaI salt solutions: Temperature effects using all-atom molecular dynamics simulations

In this study, we examine the temperature dependence of free energetics of nanotube association using graphical processing unit-enabled all-atom molecular dynamics simulations (FEN ZI) with two (10,10) single-walled carbon nanotubes in 3 m NaI aqueous salt solution. Results suggest that the free energy, enthalpy and entropy changes for the association process are all reduced at the high temperature, in agreement with previous investigations using other hydrophobes. Via the decomposition of free energy into individual components, we found that solvent contribution (including water, anion, and cation contributions) is correlated with the spatial distribution of the corresponding species and is influenced distinctly by the temperature. We studied the spatial distribution and the structure of the solvent in different regions: intertube, intratube and the bulk solvent. By calculating the fluctuation of coarse-grained tube-solvent surfaces, we found that tube-water interfacial fluctuation exhibits the strongest temperature dependence. By taking ions to be a solvent-like medium in the absence of water, tube-anion interfacial fluctuation shows similar but weaker dependence on temperature, while tube-cation interfacial fluctuation shows no dependence in general. These characteristics are discussed via the malleability of their corresponding solvation shells relative to the nanotube surface. Hydrogen bonding profiles and tetrahedrality of water arrangement are also computed to compare the structure of solvent in the solvent bulk and intertube region. The hydrophobic confinement induces a relatively lower concentration environment in the intertube region, therefore causing different intertube solvent structures which depend on the tube separation. This study is relevant in the continuing discourse on hydrophobic interactions (as they impact generally a broad class of phenomena in biology, biochemistry, and materials science and soft condensed matter research), and interpretations of hydrophobicity in terms of alternative but parallel signatures such as interfacial fluctuations, dewetting transitions, and enhanced fluctuation probabilities at interfaces.

Analyzing ion distributions around DNA: sequence-dependence of potassium ion distributions from microsecond molecular dynamics

Microsecond molecular dynamics simulations of B-DNA oligomers carried out in an aqueous environment with a physiological salt concentration enable us to perform a detailed analysis of how potassium ions interact with the double helix. The oligomers studied contain all 136 distinct tetranucleotides and we are thus able to make a comprehensive analysis of base sequence effects. Using a recently developed curvilinear helicoidal coordinate method we are able to analyze the details of ion populations and densities within the major and minor grooves and in the space surrounding DNA. The results show higher ion populations than have typically been observed in earlier studies and sequence effects that go beyond the nature of individual base pairs or base pair steps. We also show that, in some special cases, ion distributions converge very slowly and, on a microsecond timescale, do not reflect the symmetry of the corresponding base sequence.

DNA exit ramps are revealed in the binding landscapes obtained from simulations in helical coordinates

DNA molecules are highly charged semi-flexible polymers that are involved in a wide variety of dynamical processes such as transcription and replication. Characterizing the binding landscapes around DNA molecules is essential to understanding the energetics and kinetics of various biological processes. We present a curvilinear coordinate system that fully takes into account the helical symmetry of a DNA segment. The latter naturally allows to characterize the spatial organization and motions of ligands tracking the minor or major grooves, in a motion reminiscent of sliding. Using this approach, we performed umbrella sampling (US) molecular dynamics (MD) simulations to calculate the three-dimensional potentials of mean force (3D-PMFs) for a Na+ cation and for methyl guanidinium, an arginine analog. The computed PMFs show that, even for small ligands, the free energy landscapes are complex. In general, energy barriers of up to ~5 kcal/mol were measured for removing the ligands from the minor groove, and of ~1.5 kcal/mol for sliding along the minor groove. We shed light on the way the minor groove geometry, defined mainly by the DNA sequence, shapes the binding landscape around DNA, providing heterogeneous environments for recognition by various ligands. For example, we identified the presence of dissociation points or "exit ramps" that naturally would terminate sliding. We discuss how our findings have important implications for understanding how proteins and ligands associate and slide along DNA.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

Close encounters with DNA

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

Structure and mechanical characterization of DNA i-motif nanowires by molecular dynamics simulation

We studied the structure and mechanical properties of DNA i-motif nanowires by means of molecular dynamics computer simulations. We built up to 230 nm-long nanowires, based on a repeated TC5 sequence from crystallographic data, fully relaxed and equilibrated in water. The unusual C⋅C(+) stacked structure, formed by four ssDNA strands arranged in an intercalated tetramer, is here fully characterized both statically and dynamically. By applying stretching, compression, and bending deformations with the steered molecular dynamics and umbrella sampling methods, we extract the apparent Young's and bending moduli of the nanowire, as well as estimates for the tensile strength and persistence length. According to our results, the i-motif nanowire shares similarities with structural proteins, as far as its tensile stiffness, but is closer to nucleic acids and flexible proteins, as far as its bending rigidity is concerned. Furthermore, thanks to its very thin cross section, the apparent tensile toughness is close to that of a metal. Besides their yet to be clarified biological significance, i-motif nanowires may qualify as interesting candidates for nanotechnology templates, due to such outstanding mechanical properties.

Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamics simulation and a Poisson-Boltzmann approach

The ion atmosphere created by monovalent (Na(+) ) or divalent (Mg(2+) ) cations surrounding a B-form DNA duplex were examined using atomistic molecular dynamics (MD) simulations and the nonlinear Poisson-Boltzmann (PB) equation. The ion distributions predicted by the two methods were compared using plots of radial and two-dimensional cation concentrations and by calculating the total number of cations and net solution charge surrounding the DNA. Na(+) ion distributions near the DNA were more diffuse in PB calculations than in corresponding MD simulations, with PB calculations predicting lower concentrations near DNA groove sites and phosphate groups and a higher concentration in the region between these locations. Other than this difference, the Na(+) distributions generated by the two methods largely agreed, as both predicted similar locations of high Na(+) concentration and nearly identical values of the number of cations and the net solution charge at all distances from the DNA. In contrast, there was greater disagreement between the two methods for Mg(2+) cation concentration profiles, as both the locations and magnitudes of peaks in Mg(2+) concentration were different. Despite experimental and simulation observations that Mg(2+) typically maintains its first solvation shell when interacting with nucleic acids, modeling Mg(2+) as an unsolvated ion during PB calculations improved the agreement of the Mg(2+) ion atmosphere predicted by the two methods and allowed for values of the number of bound ions and net solution charge surrounding the DNA from PB calculations that approached the values observed in MD simulations.

Structure, stability and behaviour of nucleic acids in ionic liquids

Nucleic acids have become a powerful tool in nanotechnology because of their conformational polymorphism. However, lack of a medium in which nucleic acid structures exhibit long-term stability has been a bottleneck. Ionic liquids (ILs) are potential solvents in the nanotechnology field. Hydrated ILs, such as choline dihydrogen phosphate (choline dhp) and deep eutectic solvent (DES) prepared from choline chloride and urea, are 'green' solvents that ensure long-term stability of biomolecules. An understanding of the behaviour of nucleic acids in hydrated ILs is necessary for developing DNA materials. We here review current knowledge about the structures and stabilities of nucleic acids in choline dhp and DES. Interestingly, in choline dhp, A-T base pairs are more stable than G-C base pairs, the reverse of the situation in buffered NaCl solution. Moreover, DNA triplex formation is markedly stabilized in hydrated ILs compared with aqueous solution. In choline dhp, the stability of Hoogsteen base pairs is comparable to that of Watson-Crick base pairs. Moreover, the parallel form of the G-quadruplex is stabilized in DES compared with aqueous solution. The behaviours of various DNA molecules in ILs detailed here should be useful for designing oligonucleotides for the development of nanomaterials and nanodevices.

Balancing the interactions of ions, water, and DNA in the Drude polarizable force field

Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems.

DNA conformational flexibility study using phosphate backbone neutralization model

Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.

Origin of current blockades in nanopore translocation experiments

We present a detailed investigation of the ionic current in a cylindrical model nanopore in the absence and the presence of a double stranded DNA homopolymer. Our atomistic simulations are capable of reproducing almost exactly the experimental data obtained by Smeets et al., including notably the crossover salt concentration that yields equal current measurements in both situations. We can rule out that the observed current blockade is due to the steric exclusion of charge carriers from the DNA, since for all investigated salt concentrations the charge carrier density is higher when the DNA is present. Calculations using a mean-field electrokinetic model proposed by van Dorp et al. fail quantitatively in predicting this effect. We can relate the shortcomings of the mean-field model to a surface related molecular drag that the ions feel in the presence of the DNA. This drag is independent of the salt concentration and originates from electrostatic, hydrodynamic, and excluded volume interactions.

Control of stability and structure of nucleic acids using cosolutes

The stabilities, structures, and functions of nucleic acids are responsive to surrounding conditions. Living cells contain biomolecules, including nucleic acids, proteins, polysaccharides, and other soluble and insoluble low-molecular weight components, that occupy a significant fraction of the cellular volume (up to 40%), resulting in a highly crowded intracellular environment. We have proven that conditions that mimic features of this intra-cellular environment alter the physical properties affect the stability, structure, and function of nucleic acids. The ability to control structure of nucleic acids by mimicking intra-cellular conditions will be useful in nanotechnology applications of nucleic acids. This paper describes methods that can be used to analyze quantitatively the intra-cellular environment effects caused by cosolutes on nucleic acid structures and to regulate properties of nucleic acids using cosolutes.

Probing the electrostatics and pharmacological modulation of sequence-specific binding by the DNA-binding domain of the ETS family transcription factor PU.1: a binding affinity and kinetics investigation

Members of the ETS family of transcription factors regulate a functionally diverse array of genes. All ETS proteins share a structurally conserved but sequence-divergent DNA-binding domain, known as the ETS domain. Although the structure and thermodynamics of the ETS-DNA complexes are well known, little is known about the kinetics of sequence recognition, a facet that offers potential insight into its molecular mechanism. We have characterized DNA binding by the ETS domain of PU.1 by biosensor-surface plasmon resonance (SPR). SPR analysis revealed a striking kinetic profile for DNA binding by the PU.1 ETS domain. At low salt concentrations, it binds high-affinity cognate DNA with a very slow association rate constant (≤10(5)M(-)(1)s(-)(1)), compensated by a correspondingly small dissociation rate constant. The kinetics are strongly salt dependent but mutually balance to produce a relatively weak dependence in the equilibrium constant. This profile contrasts sharply with reported data for other ETS domains (e.g., Ets-1, TEL) for which high-affinity binding is driven by rapid association (>10(7)M(-)(1)s(-)(1)). We interpret this difference in terms of the hydration properties of ETS-DNA binding and propose that at least two mechanisms of sequence recognition are employed by this family of DNA-binding domain. Additionally, we use SPR to demonstrate the potential for pharmacological inhibition of sequence-specific ETS-DNA binding, using the minor groove-binding distamycin as a model compound. Our work establishes SPR as a valuable technique for extending our understanding of the molecular mechanisms of ETS-DNA interactions as well as developing potential small-molecule agents for biotechnological and therapeutic purposes.

K(+) , Na(+) , and Mg(2+) on DNA translocation in silicon nitride nanopores

In this work, we report on how salt concentration and cation species affect DNA translocation in voltage-biased silicon nitride nanopores. The translocation of dsDNA in linear, circular, and supercoiled forms was measured in salt solutions containing KCl, NaCl, and MgCl(2) . As the KCl concentrations were decreased from 1 to 0.1 M, the time taken by a DNA molecule to pass through a nanopore was shorter and the frequency of the translocation in a folded configuration was reduced, suggesting an increase in DNA electrophoretic mobility and DNA persistence length. When the salt concentration was kept at 1 M, but replacing K(+) with Na(+) , longer DNA translocation times (t(d) ) were observed. The addition of low concentrations of MgCl(2) with 1.6 M KCl resulted in longer t(d) and an increased frequency of supercoiled DNA molecules in a branched form. These observations were consistent with the greater counterion charge screening ability of Na(+) and Mg(2+) as compared to K(+) . In addition, we demonstrated that dsDNA molecules indeed translocated through a ∼10 nm nanopore by PCR amplification and gel electrophoresis. We also compared the dependence of DNA mobility and conformation on KCl concentration and cation species measured at single molecule level by silicon nitride nanopores with existing bulk-based experimental results and theoretical predictions.

Effect of initial ion positions on the interactions of monovalent and divalent ions with a DNA duplex as revealed with atomistic molecular dynamics simulations

Monovalent (Na(+)) and divalent (Mg(2+)) ion distributions around the Dickerson-Drew dodecamer were studied by atomistic molecular dynamics (MD) simulations with AMBER molecular modeling software. Different initial placements of ions were tried and the resulting effects on the ion distributions around DNA were investigated. For monovalent ions, results were found to be nearly independent of initial cation coordinates. However, Mg(2+) ions demonstrated a strong initial coordinate dependent behavior. While some divalent ions initially placed near the DNA formed essentially permanent direct coordination complexes with electronegative DNA atoms, Mg(2+) ions initially placed further away from the duplex formed a full, nonexchanging, octahedral first solvation shell. These fully solvated cations were still capable of binding with DNA with events lasting up to 20 ns, and in comparison were bound much longer than Na(+) ions. Force field parameters were also investigated with modest and little differences arising from ion (ions94 and ions08) and nucleic acid description (ff99, ff99bsc0, and ff10), respectively. Based on known Mg(2+) ion solvation structure, we conclude that in most cases Mg(2+) ions retain their first solvation shell, making only solvent-mediated contacts with DNA duplex. The proper way to simulate Mg(2+) ions around DNA duplex, therefore, should begin with ions placed in the bulk water.

The ionic atmosphere around A-RNA: Poisson-Boltzmann and molecular dynamics simulations

The distributions of different cations around A-RNA are computed by Poisson-Boltzmann (PB) equation and replica exchange molecular dynamics (MD). Both the nonlinear PB and size-modified PB theories are considered. The number of ions bound to A-RNA, which can be measured experimentally, is well reproduced in all methods. On the other hand, the radial ion distribution profiles show differences between MD and PB. We showed that PB results are sensitive to ion size and functional form of the solvent dielectric region but not the solvent dielectric boundary definition. Size-modified PB agrees with replica exchange molecular dynamics much better than nonlinear PB when the ion sizes are chosen from atomistic simulations. The distribution of ions 14 Å away from the RNA central axis are reasonably well reproduced by size-modified PB for all ion types with a uniform solvent dielectric model and a sharp dielectric boundary between solvent and RNA. However, this model does not agree with MD for shorter distances from the A-RNA. A distance-dependent solvent dielectric function proposed by another research group improves the agreement for sodium and strontium ions, even for shorter distances from the A-RNA. However, Mg(2+) distributions are still at significant variances for shorter distances.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Dynamics of ion migration in nanopores and the effect of DNA-ion interaction

We have carried out Brownian dynamics calculations to investigate the effect of DNA-ion interaction on ion transport in a nanopore. We calculated the self-diffusion coefficient of monovalent ions in the presence of DNA in a nanopore and compared the result with that through an open pore, that is, without the presence of DNA. We find that the self-diffusion coefficient of the co-ions is essentially unaffected by the DNA. The self-diffusion coefficient of the counterions, on the other hand, is significantly reduced depending on the ion concentration. At high ion concentration, around 1 M, the effect of DNA on counterion diffusion is relatively small, causing a slight reduction in counterion diffusion coefficient. At low concentrations of a few millimolar, the effect is much larger, resulting in a reduction in counterion diffusion coefficient by a factor of about 2.5. The variation in the self-diffusion of the counterion is well described by accounting for the contributions from two components: the adsorbed counterions and the free counterions. Detailed dynamics of the DNA-counterion interaction is characterized by the varying length of the transient adsorption time of the counterions to the DNA charge sites and the exchange rate with the environment. This variation in counterion adsorption time is attributed to the ionic electric screening effect, which is in turn determined by the ion concentration.

Infrared multiple photon dissociation action spectroscopy and theoretical studies of diethyl phosphate complexes: effects of protonation and sodium cationization on structure

The gas-phase structures of deprotonated, protonated, and sodium-cationized complexes of diethyl phosphate (DEP) including [DEP - H](-), [DEP + H](+), [DEP + Na](+), and [DEP - H + 2Na](+) are examined via infrared multiple photon dissociation (IRMPD) action spectroscopy using tunable IR radiation generated by a free electron laser, a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) with an electrospray ionization (ESI) source, and theoretical electronic structure calculations. Measured IRMPD spectra are compared to linear IR spectra calculated at the B3LYP/6-31G(d,p) level of theory to identify the structures accessed in the experimental studies. For comparison, theoretical studies of neutral complexes are also performed. These experiments and calculations suggest that specific geometric changes occur upon the binding of protons and/or sodium cations, including changes correlating to nucleic acid backbone geometry, specifically P-O bond lengths and ∠OPO bond angles. Information from these observations may be used to gain insight into the structures of more complex systems, such as nucleotides and solvated nucleic acids.

Ion and solvent density distributions around canonical B-DNA from integral equations

We calculate the water and ion spatial distributions around charged oligonucleotides using a renormalized three-dimensional reference interaction site theory coupled with the HNC closure. Our goal is to understand the balance between inter-DNA strand forces and solvation forces as a function of oligonucleotide length in the short strand limit. The DNA is considered in aqueous electrolyte solutions of 1 M KCl, 0.1 M KCl, or 0.1 M NaCl. The current theoretical results are compared to molecular dynamics (MD) simulations and experiments. It is found that the integral equation (IE) theory replicates the MD and the experimental results for the base-specific hydration patterns in both the major and the minor grooves. We are also able to discern characteristic structural pattern differences between Na(+) and K(+) ions. When compared to Poisson-Boltzmann methods, the IE theory, like simulation, predicts a richly structured ion environment, which is better described as multilayer rather than double layer.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.

A systematically coarse-grained model for DNA and its predictions for persistence length, stacking, twist, and chirality

We introduce a coarse-grained model of DNA with bases modeled as rigid-body ellipsoids to capture their anisotropic stereochemistry. Interaction potentials are all physicochemical and generated from all-atom simulation/parameterization with minimal phenomenology. Persistence length, degree of stacking, and twist are studied by molecular dynamics simulation as functions of temperature, salt concentration, sequence, interaction potential strength, and local position along the chain for both single- and double-stranded DNA where appropriate. The model of DNA shows several phase transitions and crossover regimes in addition to dehybridization, including unstacking, untwisting, and collapse, which affect mechanical properties such as rigidity and persistence length. The model also exhibits chirality with a stable right-handed and metastable left-handed helix.

Mediating role of multivalent cations in DNA electrostatics: an epsilon-modified Poisson-Boltzmann study of B-DNA-B-DNA interactions in mixture of NaCl and MgCl2 solutions

Potentials of mean force acting between two ions in SPC/E water have been determined via molecular dynamics simulations using the spherical cavity approach ( J. Phys. Chem. B 2006 , 110 , 10878 ). The potentials were obtained for Me(2+)-Me(+) pairs, where Me(2+) means cations Mg(2+) and Ca(2+) and Me(+) denotes monovalent ions Li(+), Na(+), and K(+). The hard-core interaction distance for effective Me(2+)-Me(+) potentials appears to be of about 5 A that looks like a sum of the effective radii of a Me(2+) ion (3 A) and of an alkali metal ion Me(+) (about 2 A). These ion-ion interaction parameters were used in the epsilon-Modified Poisson-Boltzmann (epsilon-MPB) calculations ( J. Phys. Chem. B 2007 , 111 , 5264 ) of ionic distributions around DNA generalized for the arbitrary mixture of different ion species. Ionic distributions around an all-atom geometry model of B-DNA in solution of a mixture of NaCl and MgCl(2) were obtained. It was found that even a small fraction of ions Mg(2+) led to sharp condensation of Mg(2+) near the phosphate groups of DNA due to polarization deficiency of cluster [Mg(H(2)O)(6)](2+) in an external field. The epsilon-MPB calculations of the B-DNA-B-DNA interaction energies suggest that adding 1 mM of Mg(2+) to 50 mM solution of NaCl notably affects the force acting between the two macromolecules. Being compared to Poisson-Boltzmann results and to MPB calculations for the primitive model of ions, the epsilon-MPB results also indicate an important contribution of dielectric saturation effects to the mediating role of divalent cations in the DNA-DNA interaction energies.

Monovalent cation binding in the minor groove of DNA A-tracts

The binding of five different monovalent cations to DNA oligomers containing A-tracts, runs of four or more contiguous adenine residues, has been assessed by capillary electrophoresis, using the Replacement Ion method. In this method, a nonbinding cation in the background electrolyte is gradually replaced by a binding cation, keeping the ionic strength of the solution constant. Monovalent cation binding reduces the effective charge of an A-tract-containing oligomer, decreasing its free solution mobility. The cations bind in the A-tract minor groove, because the binding site can be blocked by the minor groove binding drug netropsin. Li(+), NH(4)(+), and Tris(+) ions bind to A-tracts with similar affinities; the binding of Na(+) ions is weaker, and K(+) ion binding is highly variable. Each A-tract appears to bind one monovalent cation upon saturation of the binding site(s). For a given cation, the apparent dissociation constants depend on A-tract sequence and orientation, but not on the phasing of the A-tracts with respect to the helix repeat. Differences in the cooperativity of binding of the various cations to A-tracts with different sequences suggest that monovalent cation binding may be coupled with a conformational transition leading to the formation of the characteristic narrow minor groove A-tract structure.

Assessing the performance of implicit solvation models at a nucleic acid surface

Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.

Electrostatics of B-DNA in NaCl and CaCl2 solutions: ion size, interionic correlation, and solvent dielectric saturation effects

The epsilon-modified Poisson-Boltzmann (-MPB) equations ( J. Phys. Chem. B, 2007, 111, 5264) have been solved on a three-dimensional grid for an all-atom geometry model of B-DNA. The approach is based on the implicit solvent model including finite sizes of hydrated ions and a dielectric approximation of the ion hydration shell. Results were obtained for the detailed geometry model of B-DNA in dilute and moderately concentrated solutions of NaCl and CaCl(2). All -MPB parameters of ions and dielectric medium were extracted from published results of all-atom molecular dynamics simulations. The study allows evaluations of the ion size, interionic correlation, and the solvent dielectric saturation effects on the ion distributions around DNA. It unambiguously suggests that the difference between the -MPB and Poisson-Boltzmann distributions of ions is low for Na(+) counterions. Such a difference in the case of divalent counterions Ca(2+) is dramatic: the dielectric saturation of the ion hydration shell leads to point-like adsorption of Ca(2+) on the phosphate groups of DNA. The -MPB equations were also applied to calculate the energy of interaction between two B-DNA molecules. Results agree with previously published simulations and experimental data. Some aspects of ion specificity of polyelectrolyte properties are discussed.

Packing structures of single-stranded DNA and double-stranded DNA thiolates on Au(111): a molecular simulation study

The packing structures of single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) thiolates on implicit gold surfaces were studied in explicit aqueous solutions of 1M NaCl using molecular dynamics simulations. The simulations were based on individual DNA chains placed in hexagonal simulation boxes of different sizes, representing various packing densities. The total potential energy per DNA chain was compared. The optimal packing structures were determined based on the minimal potential energy within the limits of the conditions that were evaluated in this study. The optimal packing density of ssDNA was found to be 0.19 DNA chains/nm(2), which is consistent with that determined experimentally. Furthermore, the optimal packing density of dsDNA was shown to be approximately 58% of the packing density for ssDNA, indicating that the packing of ssDNA should be approximately 58% of its optimal packing in order to achieve the best hybridization.

Free Energy Calculations for DNA Near Surfaces Using an Ellipsoidal Geometry

The change in some thermodynamic quantities such as Gibbs' free energy, entropy and enthalpy of the binding of two DNA strands (forming a double helix), while one is tethered to a surface and are analytically calculated. These particles are submerged in an electrolytic solution; the ionic strength of the media allows the linearized version of the Poisson-Boltzmann equation (from the theory of the double layer interaction) to properly describe the interactions [13]. There is experimental and computational evidence that an ion penetrable ellipsoid is an adequate model for the single strand and the double helix [22-25]. The analytic solution provides simple calculations useful for DNA chip design. The predicted electrostatic effects suggest the feasibility of electronic control and detection of DNA hybridization in the fast growing area of DNA recognition.

Molecular dynamics studies of ion distributions for DNA duplexes and DNA clusters: salt effects and connection to DNA melting

We present extensive molecular dynamics simulations of the ion distributions for DNA duplexes and DNA clusters using the Amber force field with implicit water. The distribution of ions and the electrostatic energy of ions around an isolated DNA duplex and clusters of DNA duplexes in different salt (NaCl) concentrations over the range 0.2-1.0 mol/L are determined on the basis of the simulation results. Using the electrostatic energy profile, we determine a local net charge fraction phi, which is found to increase with increasing of salt concentration. For DNA clusters containing two DNA duplexes (DNA pair) or four DNA duplexes, phi increases as the distance between the duplexes decreases. Combining this result with experimental results for the dependence of the DNA melting temperature on bulk salt concentration, we conclude that for a pair of DNA duplexes the melting temperature increases by 5-10 K for interaxis separations of 25-40 A. For a cluster of four DNA duplexes, an even larger melting temperature increase should occur. We argue that this melting temperature increase in dense DNA clusters is responsible for the cooperative melting mechanism in DNA-linked nanoparticle aggregates and DNA-linked polymer aggregates.

Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment

Molecular dynamics simulations of biomolecules with implicit solvent reduce the computational cost and complexity of such simulations so that longer time scales and larger system sizes can be reached. While implicit solvent simulations of proteins have become well established, the success of implicit solvent in the simulation of nucleic acids has not been fully established to date. Results obtained in this study demonstrate that stable and efficient simulations of DNA and a protein-DNA complex can be achieved with an implicit solvent model based on continuum dielectric electrostatics. Differences in conformational sampling of DNA with two sets of atomic radii that are used to define the dielectric interface between the solute and the continuum dielectric model of the solvent are investigated. Results suggest that depending on the choice of atomic radii agreement is either closer to experimental data or to explicit solvent simulations. Furthermore, partial conformational transitions toward A-DNA conformations when salt is added within the implicit solvent framework are observed.

Ion atmosphere around nucleic acid

The poyion-ion preferential interaction coefficient Gamma describes the exclusion of coions and accumulations of counterions in the vicinity of a polyion in an aqueous solution. We give tight upper and lower bounds for Gamma when the polyion can be modeled by a cylinder of infinite length but of arbitrary charge density. This case can be used as a model for long strands of DNA or RNA in an aqueous solution containing univalent cations. The salt dependence of Gamma is predicted from low to intermediate and high salt concentrations. We also indicate how the bounds for the infinite polyion can be exploited to place bounds for polyions of length greater than a constant on the order of the inverse Debye screening length.

Nucleic acid solvation: from outside to insight

Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.

Removal of polycyclic aromatic hydrocarbons from contaminated soil by aqueous DNA solution

An aqueous DNA solution was applied for the extraction of polycyclic aromatic hydrocarbons (PAHs) from a spiked soil. Solubilities of 0.56, 11.78, and 11.24 mg L(-1) for anthracene, phenanthrene, and pyrene, respectively, were achieved after 1 day equilibration in 1% DNA. Using a spiked soil that contained 72 mg kg(-1) anthracene, 102 mg kg(-1) phenanthrene, and 99 mg kg(-1) pyrene, extractions of close to 88, 78, and 94%, respectively, were attained with 5% DNA at a 1:50 soil/extractant ratio. Maximum PAH dissolution occurred after 4-6 h. Comparative tests showed the main advantage of DNA over methyl-beta- and gamma-cyclodextrins and Tween 80 for pyrene removal. An ionic strength of 0.1 M NaCl was found necessary for maximum PAH dissolution and extraction. The performance of hexane regenerated DNA also remained stable after three stages of recycling. These results show the huge potential of DNA as an aqueous washing agent for PAH-contaminated soil.

Visualization of Different Pathways of DNA Release from Interpolyelectrolyte Complex

The release of double-stranded DNA from its interpolyelectrolyte complex with positively charged poly(allylamine hydrochloride) via exchange reaction with added polyanion, poly(sodium styrenesulfonate), is directly observed by fluorescence microscopy. It is shown that the pathways of DNA release depend essentially on the amount of added low-molecular-weight salt. At low salt content, the DNA release proceeds via the formation of an intermediate "beads-on-string" structure, whereas at high salt content the release goes directly from globule to coil states without any intermediate structures. The reasons for different character of DNA release are discussed.

Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site

The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.