Analyzing ion distributions around DNA

We present a new method for analyzing ion, or molecule, distributions around helical nucleic acids and illustrate the approach by analyzing data derived from molecular dynamics simulations. The analysis is based on the use of curvilinear helicoidal coordinates and leads to highly localized ion densities compared to those obtained by simply superposing molecular dynamics snapshots in Cartesian space. The results identify highly populated and sequence-dependent regions where ions strongly interact with the nucleic and are coupled to its conformational fluctuations. The data from this approach is presented as ion populations or ion densities (in units of molarity) and can be analyzed in radial, angular and longitudinal coordinates using 1D or 2D graphics. It is also possible to regenerate 3D densities in Cartesian space. This approach makes it easy to understand and compare ion distributions and also allows the calculation of average ion populations in any desired zone surrounding a nucleic acid without requiring references to its constituent atoms. The method is illustrated using microsecond molecular dynamics simulations for two different DNA oligomers in the presence of 0.15 M potassium chloride. We discuss the results in terms of convergence, sequence-specific ion binding and coupling with DNA conformation.

[1]  G. Hong,et al.  Nucleic Acids Research , 2015, Nucleic Acids Research.

[2]  Sebastian Doniach,et al.  Understanding nucleic acid-ion interactions. , 2014, Annual review of biochemistry.

[3]  D. Case,et al.  Ion counting from explicit-solvent simulations and 3D-RISM. , 2014, Biophysical journal.

[4]  Aatto Laaksonen,et al.  Insight into nucleic acid counterion interactions from inside molecular dynamics simulations is “worth its salt” , 2012 .

[5]  M. Mezei,et al.  Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations , 2012, Journal of Biosciences.

[6]  N. Hud,et al.  B-DNA structure is intrinsically polymorphic: even at the level of base pair positions , 2011, Nucleic acids research.

[7]  Christophe Blanchet,et al.  CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures , 2011, Nucleic Acids Res..

[8]  L. Pollack.  SAXS studies of ion-nucleic acid interactions. , 2011, Annual review of biophysics.

[9]  Gillian C. Lynch,et al.  Ion and solvent density distributions around canonical B-DNA from integral equations. , 2011, The journal of physical chemistry. B.

[10]  D. Case,et al.  A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA , 2009, Nucleic acids research.

[11]  J. H. Maddocks,et al.  Conformational analysis of nucleic acids revisited: Curves+ , 2009, Nucleic acids research.

[12]  Wei Yang,et al.  Efficient sampling of ion motions in molecular dynamics simulations on DNA: Variant Hamiltonian replica exchange method , 2008 .

[13]  D. Herschlag,et al.  Quantitative and comprehensive decomposition of the ion atmosphere around nucleic acids. , 2007, Journal of the American Chemical Society.

[14]  F. J. Luque,et al.  Dynamics of B-DNA on the microsecond time scale. , 2007, Journal of the American Chemical Society.

[15]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[16]  N. Korolev,et al.  Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study , 2006, Nucleic acids research.

[17]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[18]  Heinz Sklenar,et al.  Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps. , 2005, Biophysical journal.

[19]  Davenport Rj.  Worth its salt. Elusive enzyme generates hormone that lowers blood pressure. , 2005 .

[20]  Heinz Sklenar,et al.  Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps. , 2004, Biophysical journal.

[21]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[22]  Charles A Laughton,et al.  Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations? , 2004, Biophysical journal.

[23]  D. VanDerveer,et al.  Locating monovalent cations in the grooves of B-DNA. , 2001, Biochemistry.

[24]  D. Beveridge,et al.  DNA structure: what's in charge? , 2000, Journal of molecular biology.

[25]  Donald Hamelberg,et al.  Flexible Structure of DNA: Ion Dependence of Minor-Groove Structure and Dynamics , 2000 .

[26]  E Westhof,et al.  Water and ion binding around RNA and DNA (C,G) oligomers. , 2000, Journal of molecular biology.

[27]  R. Dickerson,et al.  Absence of minor groove monovalent cations in the crosslinked dodecamer C-G-C-G-A-A-T-T-C-G-C-G. , 1999, Journal of molecular biology.

[28]  L. Williams,et al.  DNA structure: cations in charge? , 1999, Current opinion in structural biology.

[29]  J. Feigon,et al.  Localization of ammonium ions in the minor groove of DNA duplexes in solution and the origin of DNA A-tract bending. , 1999, Journal of molecular biology.

[30]  P. Kollman,et al.  A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat. , 1999, Journal of biomolecular structure & dynamics.

[31]  Martin Egli,et al.  The Dickerson-Drew B-DNA Dodecamer Revisited-At Atomic Resolution , 1998 .

[32]  S. Harvey,et al.  The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition , 1998, J. Comput. Chem..

[33]  Bhyravabhotla Jayaram,et al.  Intrusion of Counterions into the Spine of Hydration in the Minor Groove of B-DNA: Fractional Occupancy of Electronegative Pockets , 1997 .

[34]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[35]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[36]  L. Dang,et al.  Mechanism and Thermodynamics of Ion Selectivity in Aqueous Solutions of 18-Crown-6 Ether: A Molecular Dynamics Study , 1995 .

[37]  J. Bond,et al.  Importance of coulombic end effects on cation accumulation near oligoelectrolyte B-DNA: a demonstration using 23Na NMR. , 1995, Biophysical journal.

[38]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[39]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[40]  M. Mezei,et al.  Monte Carlo studies of the structure of dilute aqueous sclutions of Li+, Na+, K+, F−, and Cl− , 1981 .

[41]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[42]  Conrad C. Huang,et al.  Visualizing density maps with UCSF Chimera. , 2007, Journal of structural biology.

[43]  Daniel Svozil,et al.  Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers. , 2007, Biophysical journal.

[44]  Krystyna Zakrzewska,et al.  DNA and its counterions: a molecular dynamics study. , 2004, Nucleic acids research.