AN ALL-ATOM EMPIRICAL ENERGY FUNCTION FOR THE SIMULATION OF NUCLEIC-ACIDS

被引:637
|
作者
MACKERELL, AD [1 ]
WIORKIEWICZKUCZERA, J [1 ]
KARPLUS, M [1 ]
机构
[1] UNIV MARYLAND, SCH PHARM, DEPT PHARMACEUT SCI, BALTIMORE, MD 21201 USA
关键词
D O I
10.1021/ja00153a017
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
Nucleic acid parameters are developed for the all-atom empirical energy function used in the CHARMM program. The parameters were determined by use of results for model compounds, including the nucleic acid bases, dimethyl phosphate and anionic and dianionic methyl phosphate, ribose, and deoxyribose. Internal parameters (bond length, bond angle, Urey-Bradley, dihedral, and improper dihedral terms) were chosen to reproduce geometries and vibrational spectra from experimental crystal structures, infrared and Raman spectroscopic data, and ab initio calculations. Interaction parameters (electrostatic and van der Waals terms) were derived from 6-31G* ab initio interaction energies and geometries for water molecules bonded to polar sites of the model compounds and from the experimentally measured gas phase Watson-Crick base pair energies and geometries, base heats of sublimation, and experimental and 6-31G* ab initio dipole moments. Emphasis was placed on a proper balance between solvent-solvent, solvent-solute, and solute-solute interactions with reference to the TIP3P water model. Tests on nucleic acid base crystals showed satisfactory agreement between calculated and experimental values for the lattice parameters, nonbonded interactions, and heats of sublimation. Base pair variations in stacking energies are consistent with experiment and ab initio calculations. Further testing was performed on GpC and B-DNA dodecamer crystal structures, including water molecules and counterions. Simulations of these systems revealed the parameters to accurately reproduce Watson-Crick base pairing, internal geometries including the backbone dihedrals, sugar puckering and glycosidic linkages, and the hydration of the nucleic acids. The present parameters should be useful for modeling and simulation studies of nucleic acids, including both structural and energetic analysis. Further, they provide a parameter set that is consistent with the protein parameters in CHARMM so that complexes between proteins and nucleic acids can be modeled. A list of the parameter values is included in an Appendix (supporting information). In the test simulations, a number of interesting results were obtained. Significant anharmonicity is present due to nonbonded terms in certain bond angles (P-O-C) of the phosphate group that may play a role in the observed variation of these angles. There is a wide range of stacking energies of the bases of B-DNA due to repulsive electrostatic interactions. The motion of the two strands appears to be correlated for bases involved in Watson-Crick interactions while there is a lack of correlation for the sugar and phosphate moieties. The importance of conformational substates in the B-DNA dodecamer is pointed out. Comparisons with nucleic acid simulations made with other energy functions show the importance of the parameters in determining the internal geometry and the interactions with the surroundings. A summary of these results is given in the concluding section.
引用
收藏
页码:11946 / 11975
页数:30
相关论文
共 50 条
  • [21] CHROMATOGRAPHY AND FUNCTION OF MODIFIED NUCLEOSIDES IN NUCLEIC-ACIDS
    GEHRKE, CW
    ZUMWALT, RW
    AGRIS, PF
    EHRLICH, M
    KUO, KC
    FEDERATION PROCEEDINGS, 1983, 42 (07) : 1883 - 1883
  • [22] MOBILITY AND FUNCTION IN PROTEINS AND NUCLEIC-ACIDS - INTRODUCTION
    RICHARDS, FM
    CIBA FOUNDATION SYMPOSIA, 1983, 93 : 1 - 3
  • [23] MODELING OF STRUCTURE AND FUNCTION OF PROTEINS AND NUCLEIC-ACIDS
    PLOCHOCKA, D
    RABCZENKO, A
    ACTA BIOCHIMICA POLONICA, 1991, 38 (03) : 281 - 293
  • [24] DISSIPATION OF UV ENERGY IN NUCLEIC-ACIDS AND NUCLEOPROTEINS
    LONGWORTH, JW
    RADIATION RESEARCH, 1974, 59 (01) : 175 - 176
  • [25] Single-atom test of all-atom empirical potentials: Fe in myoglobin
    Rai, BK
    Prohofsky, EW
    Durbin, SM
    JOURNAL OF PHYSICAL CHEMISTRY B, 2005, 109 (40): : 18983 - 18987
  • [26] All-atom level direct folding simulation of a ββα miniprotein
    Jang, Soonmin
    Kim, Eunae
    Pak, Youngshang
    JOURNAL OF CHEMICAL PHYSICS, 2008, 128 (10):
  • [27] All-Atom Molecular Dynamics Simulation of Stealth Liposomes
    Tabari, Seyed Hamid
    Hettige, Jeevapani
    Moradi, Mahmoud
    BIOPHYSICAL JOURNAL, 2017, 112 (03) : 75A - 75A
  • [28] All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data
    Foloppe, N
    MacKerell, AD
    JOURNAL OF COMPUTATIONAL CHEMISTRY, 2000, 21 (02) : 86 - 104
  • [29] Fracture of silica aerogels: An all-atom simulation study
    Patil, Sandeep P.
    Rege, Ameya
    Itskov, Mikhail
    Markert, Bernd
    JOURNAL OF NON-CRYSTALLINE SOLIDS, 2018, 498 : 125 - 129
  • [30] Protein folding simulation by all-atom CSAW method
    Sun, Weitao
    2007 IEEE INTERNATIONAL CONFERENCE ON BIOINFORMATICS AND BIOMEDICINE WORKSHOPS, PROCEEDINGS, 2007, : 45 - 52