Dna Dynamics in Aqueous Solution: Opening the Double Helix

The opening of a DNA base pair is a simple reaction that is a prerequisite for replication, transcription, and other vital biological functions. Understanding the mo lecular mechanisms of biological reactions is crucial for predicting and, ultimately, controlling them. Realistic computer simulations of the reactions can provide the needed understanding. To model even the simplest re action in aqueous solution requires hundreds of hours of supercomputing time. We have used molecular dy namics techniques to simulate fraying of the ends of a six base pair double strand of DNA, [TCGCGA]2, where the four bases of DNA are denoted by T (thymine), C (cytosine), G (guanine), and A (adenine), and to estimate the free energy barrier to this process. The calculations, in which the DNA was surrounded by 2,594 water mol ecules, required 50 hours of CRAY-2 CPU time for every simulated 100 picoseconds. A free energy barrier to fraying, which is mainly characterized by the movement of adenine away from thymine into aqueous environ ment, was estimated to be 4 kcal/mol. Another fraying pathway, which leads to stacking between terminal adenine and thymine, was also observed. These detailed pictures of the motions and energetics of DNA base pair opening in water are a first step toward understanding how DNA will interact with any molecule.

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

[2]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[3]  I. Tinoco,et al.  Kinetics for exchange of the imino protons of the d(C-G-C-G-A-A-T-T-C-G-C-G) double helix in complexes with the antibiotics netropsin and/or actinomycin. , 1983, Biochemistry.

[4]  J. Pitha,et al.  THE SPECIFICITY OF HYDROGEN BOND FORMATION BETWEEN DERIVATIVES OF NUCLEIC ACID BASES AND SOME ANALOGUES , 1966 .

[5]  G L Seibel,et al.  A molecular dynamics simulation of double-helical B-DNA including counterions and water. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[6]  H. Blöcker,et al.  Predicting DNA duplex stability from the base sequence. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A. Pohorille,et al.  Monte Carlo simulation of the influence of solvent on nucleic acid base associations. , 1984, Journal of the American Chemical Society.

[8]  Jacob Katzenelson Computational structure of the N-body problem , 1989 .

[9]  S. Chan,et al.  Interaction and association of bases and nucleosides in aqueous solutions. IV. Proton magnetic resonance studies of the association of pyrimidine nucleosides and their interactions with purine. , 1963, Journal of the American Chemical Society.

[10]  H. Berendsen,et al.  A Molecular Dynamics Computer Simulation of an Eight‐Base‐Pair DNA Fragment in Aqueous Solution: Comparison with Experimental Two‐Dimensional NMR Data a , 1986, Annals of the New York Academy of Sciences.

[11]  I. Tinoco,et al.  Kinetics for exchange of imino protons in the d(C-G-C-G-A-A-T-T-C-G-C-G) double helix and in two similar helices that contain a G . T base pair, d(C-G-T-G-A-A-T-T-C-G-C-G), and an extra adenine, d(C-G-C-A-G-A-A-T-T-C-G-C-G). , 1982, Biochemistry.

[12]  P. Kollman,et al.  Molecular mechanical studies of base‐pair opening in d(CGCGC):d(GCGCG), dG5·dC5, d(TATAT):d(ATATA), and dA5·dT5 in the B and Z forms of DNA , 1984 .

[13]  M. P. Schweizer,et al.  Studies of inter- and intramolecular interaction in mononucleotides by proton magnetic resonance. , 1968, Journal of the American Chemical Society.

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

[15]  L. Greengard The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[16]  M Karplus,et al.  Dynamics of DNA oligomers. , 1983, Journal of biomolecular structure & dynamics.

[17]  Peter A. Kollman,et al.  Molecular dynamics simulations study of the free energy of association of 9-methyladenine and 1-methylthymine bases in water , 1990 .

[18]  Ron Elber,et al.  Calculations of a list of neighbors in Molecular Dynamics simulations , 1989 .

[19]  L. Sukhodub,et al.  Experimental studies of molecular interactions between nitrogen bases of nucleic acids , 1979, Biopolymers.

[20]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[21]  Bernard Pettitt,et al.  Chloride ion pairs in water , 1987 .

[22]  S. Penman,et al.  Association by hydrogen bonding of free nucleosides in non-aqueous solution. , 1966, Journal of molecular biology.

[23]  A. Pohorille,et al.  Solution influence on biomolecular equilibria: nucleic acid base associations. , 1984, Journal of biomolecular structure & dynamics.

[24]  P. Rossky,et al.  Free energies of association for the sodium-dimethyl phosphate ion pair in aqueous solution , 1989 .

[25]  A. Rich,et al.  An infrared study of hydrogen bonding between adenine and uracil derivatives in chloroform solution. , 1967, Journal of the American Chemical Society.

[26]  W. L. Jorgensen,et al.  Energetics and hydration of the constituent ion pairs of tetramethylammonium chloride , 1989 .

[27]  P. Weiner,et al.  Temperature-dependent molecular dynamics and restrained X-ray refinement simulations of a Z-DNA hexamer. , 1986, Journal of molecular biology.

[28]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[29]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[30]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[31]  B. Montgomery Pettitt,et al.  Alkali halides in water: Ion–solvent correlations and ion–ion potentials of mean force at infinite dilution , 1986 .

[32]  A. Pohorille,et al.  Theoretical methods for obtaining free energies of biomolecular equilibria in aqueous solutions. , 1986, Methods in enzymology.

[33]  A. Laaksonen,et al.  Molecular dynamics simulation of double helix Z-DNA in solution , 1989 .