Collective‐variable Monte Carlo simulation of DNA

Monte Carlo simulations have been carried out on DNA oligomers using an internal coordinate model associated with a pseudorotational representation of sugar repuckering. It is shown that when this model is combined with the scaled collective variable approach of Noguti and Go, much more efficient simulations are obtained than with simple single variable steps. Application of this method to a DNA oligomer containing a recognition site for the TATA‐box binding protein leads to striking similarities with results recently obtained from a 1‐ns molecular dynamics simulation using explicit solvent and counterions. In particular, large amplitude bending fluctuations are observed directed toward the major groove. Conformational analysis of the Monte Carlo simulation shows clear base sequence effects on conformational fluctuations and also that the DNA energy hypersurface, like that of proteins, is complex with many local, conformational substates. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 2001–2011, 1997

[1]  Twists and Turns in DNA: Predicting Base Sequence Effects on the Conformation of the Double Helix , 1995 .

[2]  R. H. Ritchie,et al.  Dielectric effects in biopolymers: The theory of ionic saturation revisited , 1985 .

[3]  K. J. Miller.  Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B‐DNA structure and a counterclockwise helix , 1979, Biopolymers.

[4]  Bhyravabhotla Jayaram,et al.  Monte Carlo Simulation Studies on the Structure of the Counterion Atmosphere of B-DNA. Variations on the Primitive Dielectric Model , 1990 .

[5]  J. Mccammon,et al.  Simulation methods for protein structure fluctuations , 1980, Biopolymers.

[6]  Henry A. Gabb,et al.  Efficient conformational space sampling for nucleosides using internal coordinate Monte Carlo simulations and a modified furanose description , 1995, J. Comput. Chem..

[7]  K. Zakrzewska,et al.  Calculation and analysis of low frequency normal modes for DNA , 1997 .

[8]  R L Jernigan,et al.  Static and statistical bending of DNA evaluated by Monte Carlo simulations. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[9]  R Lavery,et al.  Local DNA stretching mimics the distortion caused by the TATA box-binding protein. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[10]  S C Harvey,et al.  Molecular modelling of (A4T4NN)n and (T4A4NN)n: sequence elements responsible for curvature. , 1996, Nucleic acids research.

[11]  R Lavery,et al.  Modelling DNA conformational mechanics. , 1994, Biophysical chemistry.

[12]  Heinz Sklenar,et al.  JUMNA (junction minimisation of nucleic acids) , 1995 .

[13]  Analytic Theory of Pseudorotation in Five-Membered Rings. Cyclopentane, Tetrahydrofuran, Ribose, and Deoxyribose , 1995 .

[14]  E. Westhof,et al.  Three-center hydrogen bonds in DNA: molecular dynamics of poly(dA).cntdot.poly(dT) , 1991 .

[15]  D. Beveridge,et al.  A Nanosecond Molecular Dynamics Trajectory for a B DNA Double Helix: Evidence for Substates , 1994 .

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

[17]  R. Lavery,et al.  Conformational properties of the TATA-box binding sequence of DNA. , 1997, Journal of biomolecular structure & dynamics.

[18]  R. Lavery,et al.  Defining the structure of irregular nucleic acids: conventions and principles. , 1989, Journal of biomolecular structure & dynamics.

[19]  V. Zhurkin,et al.  Different families of double‐stranded conformations of DNA as revealed by computer calculations , 1978 .

[20]  Steven Hahn,et al.  Crystal structure of a yeast TBP/TATA-box complex , 1993, Nature.

[21]  P. Wolynes,et al.  The energy landscapes and motions of proteins. , 1991, Science.

[22]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[23]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[24]  Nidhi Arora,et al.  Strength of hydrogen bonds in α helices , 1997 .

[25]  N Go,et al.  Efficient monte carlo method for simulation of fluctuating conformations of native proteins , 1985, Biopolymers.

[26]  M. Karplus,et al.  Multiple conformational states of proteins: a molecular dynamics analysis of myoglobin. , 1987, Science.

[27]  V. Zhurkin,et al.  Anisotropic flexibility of DNA and the nucleosomal structure. , 1979, Nucleic acids research.

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

[29]  Mahmoud Ghomi,et al.  Monte Carlo simulations on short single‐stranded oligonucleotides. I. Application to RNA trimers , 1994, J. Comput. Chem..

[30]  P. Sharp,et al.  Pre-bending of a promoter sequence enhances affinity for the TATA-binding factor , 1995, Nature.

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