Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution.

Molecular dynamics simulation in explicit solvent and continuum solvent models are applied to investigate the relative stability of A- and B-form helices for two DNA sequences, dA10-dT10 and dG10-dC10 in three structural forms. One structural form is based on an unrestrained molecular dynamics (MD) trajectory starting from a canonical B-DNA structure, the second is based on a MD trajectory starting in a canonical B-DNA structure with the sugars constrained to be C2'-endo and the third simulation started from a canonical A-DNA structure with the sugars constrained to C3'-endo puckers. For the energetic analysis, structures were taken as snapshots from nanosecond length molecular dynamics simulations computed in a consistent fashion in explicit solvent, applying the particle mesh Ewald method and the Cornell et al. force field. The electrostatic contributions to solvation free energies are computed using both a finite-difference Poisson-Boltzmann model and a pairwise Generalized Born model. The non-electrostatic contributions to the solvation free energies are estimated with a solvent accessible surface area dependent term. To estimate the gas phase component of the relative free energy between the various structures, the mean solute internal energies (determined with the Cornell et al. molecular mechanics potential including all pairwise interactions within the solute) and estimates of the solute entropy (using a harmonic approximation) were used. Consistent with experiment, the polyG-polyC (GC) structures are found to be much more A-phillic than the polyA-polyT (AT) structures, the latter being quite A-phobic. The dominant energy components responsible for this difference comes from the internal and van der Waal energies. A perhaps less appreciated difference between the GC and AT rich sequences is suggested by the calculated salt dependence which demonstrates a significantly enhanced ability to drive GC rich sequences towards an A-form structure compared to AT rich sequences. In addition to being A-phobic, the AT structure also has a noticably larger helical repeat than GC and other mixed sequence duplexes, consistent with experiment. Analysis of the average solvent density from the trajectories shows hydration patterns in qualitative agreement with experiment and previous theoretical treatments.

[1]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[2]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[3]  P. Kollman,et al.  Molecular Dynamics Simulations on Solvated Biomolecular Systems: The Particle Mesh Ewald Method Leads to Stable Trajectories of DNA, RNA, and Proteins , 1995 .

[4]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[5]  D. Beveridge,et al.  A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation. , 1997, Biophysical journal.

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

[7]  D. Beveridge,et al.  Molecular dynamics studies of the conformational preferences of a DNA double helix in water and an ethanol/water mixture: Theoretical considerations of the a ⇔ B transition , 1998 .

[8]  P. Kollman,et al.  A molecular level picture of the stabilization of A-DNA in mixed ethanol-water solutions. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Alexander D. MacKerell.  Observations on the A versus B Equilibrium in Molecular Dynamics Simulations of Duplex DNA and RNA , 1997 .

[10]  S J Weiner,et al.  Molecular dynamics simulations of d(C-G-C-G-A) X d(T-C-G-C-G) with and without "hydrated" counterions. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Theoretical account of the 'spine of hydration' in the minor groove of duplex d(CGCGAATTCGCG). , 1990, Journal of biomolecular structure & dynamics.

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

[13]  H R Drew,et al.  Structure of a B-DNA dodecamer. III. Geometry of hydration. , 1981, Journal of molecular biology.

[14]  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 .

[15]  B. Honig,et al.  Evaluation of the conformational free energies of loops in proteins , 1994, Proteins.

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

[17]  S. Swaminathan,et al.  Molecular Dynamics Simulations on the Hydration, Structure and Motions of DNA Oligomers , 1993 .

[18]  H M Berman,et al.  A systematic method for studying the spatial distribution of water molecules around nucleic acid bases. , 1993, Biophysical journal.

[19]  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 .

[20]  R Lavery,et al.  The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids. , 1988, Journal of biomolecular structure & dynamics.

[21]  M. Eisenstein,et al.  Hydration patterns and intermolecular interactions in A-DNA crystal structures. Implications for DNA recognition. , 1995, Journal of molecular biology.

[22]  B. Pettitt,et al.  Experiment vs force fields: DNA conformation from molecular dynamics simulations , 1997 .

[23]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[24]  A. Klug,et al.  Sequence-dependent helical periodicity of DNA , 1981, Nature.

[25]  R. Sarma,et al.  Structure, Motion, Interaction and Expression of Biological Macromolecules. , 1998, Journal of biomolecular structure & dynamics.

[26]  J. Ponder,et al.  Protein structure prediction using a combination of sequence homology and global energy minimization: II. Energy functions , 1998 .

[27]  M. Sundaralingam,et al.  Comparison of major groove hydration in isomorphous A-DNA octamers and dependence on base sequence and local helix geometry. , 1997, Biochemistry.

[28]  Peter A. Kollman,et al.  Molecular dynamics simulations highlight the structural differences among DNA: DNA, RNA:RNA, and DNA:RNA hybrid duplexes , 1997 .

[29]  R. Dickerson,et al.  Definitions and nomenclature of nucleic acid structure components. , 1989, Nucleic acids research.

[30]  Wolfram Saenger,et al.  Principles of Nucleic Acid Structure , 1983 .

[31]  A. Wang,et al.  Neomycin, spermine and hexaamminecobalt (III) share common structural motifs in converting B- to A-DNA. , 1996, Nucleic acids research.

[32]  G G Hu,et al.  The B-DNA dodecamer at high resolution reveals a spine of water on sodium. , 1998, Biochemistry.

[33]  Gregory D. Hawkins,et al.  Pairwise solute descreening of solute charges from a dielectric medium , 1995 .

[34]  K. Sharp,et al.  Macroscopic models of aqueous solutions : biological and chemical applications , 1993 .

[35]  R Lavery,et al.  Conformational sub-states in B-DNA. , 1992, Journal of molecular biology.

[36]  西村 善文.  W. Saenger: Principles of Nucleic Acid Structure, Springer-Verlag, New York and Berlin, 1984, xx+556ページ, 24.5×16.5cm, 14,160円 (Springer Advanced Texts in Chemistry). , 1985 .

[37]  R. Ornstein,et al.  Effect of warmup protocol and sampling time on convergence of molecular dynamics simulations of a DNA dodecamer using AMBER 4.1 and particle-mesh Ewald method. , 1997, Journal of biomolecular structure & dynamics.

[38]  David A. Case,et al.  A computational study of the role of solvation effects in reverse turn formation in the tetrapeptides APGD and APGN , 1997 .

[39]  P A Kollman,et al.  Observation of the A-DNA to B-DNA transition during unrestrained molecular dynamics in aqueous solution. , 1996, Journal of molecular biology.

[40]  Charles L. Brooks,et al.  Dielectric Continuum Models for Hydration Effects on Peptide Conformational Transitions , 1996 .

[41]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[42]  D. Case,et al.  Thermodynamics of a reverse turn motif. Solvent effects and side-chain packing. , 1997, Journal of molecular biology.

[43]  Tami I. Spector,et al.  Unrestrained Molecular Dynamics of Photodamaged DNA in Aqueous Solution , 1997 .

[44]  J. Piškur,et al.  Aggregated DNA in ethanol solution , 1995, FEBS letters.

[45]  L. J. Peck,et al.  Sequence dependence of the helical repeat of DNA in solution , 1981, Nature.

[46]  R Lavery,et al.  Conformational and helicoidal analysis of 30 PS of molecular dynamics on the d(CGCGAATTCGCG) double helix: "curves", dials and windows. , 1989, Journal of biomolecular structure & dynamics.

[47]  G. Hummer,et al.  Theoretical description of biomolecular hydration. Application to A-DNA. , 1994, Basic life sciences.

[48]  B. Honig,et al.  Free energy determinants of secondary structure formation: III. beta-turns and their role in protein folding. , 1996, Journal of molecular biology.

[49]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[50]  D W Hukins,et al.  Optimised parameters for A-DNA and B-DNA. , 1972, Biochemical and biophysical research communications.

[51]  M. Sundaralingam,et al.  Conformational analysis of the sugar ring in nucleosides and nucleotides. A new description using the concept of pseudorotation. , 1972, Journal of the American Chemical Society.

[52]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[53]  P A Kollman,et al.  Insight into the stabilization of A-DNA by specific ion association: spontaneous B-DNA to A-DNA transitions observed in molecular dynamics simulations of d[ACCCGCGGGT]2 in the presence of hexaamminecobalt(III). , 1997, Structure.

[54]  Neocles B. Leontis,et al.  Molecular modeling of nucleic acids , 1997 .

[55]  H M Berman,et al.  Hydration of the DNA bases is local. , 1995, Biophysical journal.

[56]  Benno P. Schoenborn,et al.  Neutrons in Biology , 1996, Basic Life Sciences.

[57]  P. Kollman,et al.  Molecular Dynamics Simulations Find That 3‘ Phosphoramidate Modified DNA Duplexes Undergo a B to A Transition and Normal DNA Duplexes an A to B Transition , 1997 .

[58]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[59]  Conrad C. Huang,et al.  The MIDAS display system , 1988 .

[60]  D. Beveridge,et al.  A theoretical study of the aqueous hydration of canonical B d(CGCGAATTCGCG): Monte Carlo simulation and comparison with crystallographic ordered water sites. , 1989, Journal of biomolecular structure & dynamics.

[61]  M Feig,et al.  Structural equilibrium of DNA represented with different force fields. , 1998, Biophysical journal.

[62]  Theoretical considerations on the "spine of hydration" in the minor groove of d(CGCGAATTCGCG).d(GCGCTTAAGCGC): Monte Carlo computer simulation. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[63]  B. Honig,et al.  Free energy determinants of secondary structure formation: I. alpha-Helices. , 1995, Journal of molecular biology.

[64]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[65]  S. Bram.  Secondary structure of DNA depends on base composition. , 1971, Nature: New biology.

[66]  C. Gaillard,et al.  Helical periodicity of DNA, Poly(dA) . poly(dT) and poly(dA-dT). poly(dA-dT) in solution. , 2005, European journal of biochemistry.