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 条
  • [31] All-atom molecular dynamics simulation of HPMA polymers
    Meleshko, Glib
    Kulhavy, Jiri
    Paul, Alison
    Willock, David J.
    Platts, James A.
    RSC ADVANCES, 2014, 4 (14) : 7003 - 7012
  • [32] Understanding the Function of the Cyclic Antifungal Lipopeptide Fengycin using All-Atom Md Simulation
    Sur, Sreyoshi
    Grossfield, Alan
    BIOPHYSICAL JOURNAL, 2015, 108 (02) : 84A - 84A
  • [34] All-atom simulation approaches to probe dynamics in proteins
    Vashisth, Harish
    ABSTRACTS OF PAPERS OF THE AMERICAN CHEMICAL SOCIETY, 2016, 251
  • [35] Asymmetric Periodic Boundary Conditions for All-Atom Molecular Dynamics and Coarse-Grained Simulations of Nucleic Acids
    Erban, Radek
    Togashi, Yuichi
    JOURNAL OF PHYSICAL CHEMISTRY B, 2023, 127 (38): : 8257 - 8267
  • [36] MOBILITY AND FUNCTION IN PROTEINS AND NUCLEIC-ACIDS - CLOSING REMARKS
    RICHARDS, FM
    CIBA FOUNDATION SYMPOSIA, 1983, 93 : 344 - 345
  • [37] MOBILITY AND FUNCTION IN PROTEINS AND NUCLEIC-ACIDS - SUMMARY AND OUTLOOK
    FRAUENFELDER, H
    CIBA FOUNDATION SYMPOSIA, 1983, 93 : 329 - 339
  • [38] NUCLEIC-ACIDS AND/OR THEIR COMPONENTS - A POSSIBLE ROLE IN IMMUNE FUNCTION
    ADJEI, AA
    YAMAMOTO, S
    KULKARNI, A
    JOURNAL OF NUTRITIONAL SCIENCE AND VITAMINOLOGY, 1995, 41 (01) : 1 - 16
  • [39] HOMO AND LUMO ENERGY VALUES FOR BASES OF NUCLEIC-ACIDS
    LAL, KB
    AGRAWAL, VK
    SRIVASTAVA, AK
    KRISHNA, B
    INDIAN JOURNAL OF PURE & APPLIED PHYSICS, 1983, 21 (11) : 670 - 671
  • [40] POLAROGRAPHIC INVESTIGATION OF THE HALOGEN ATOM MOBILITY IN THE HALOGENOSUBSTITUTED COMPONENTS OF NUCLEIC-ACIDS
    KADYSH, VP
    KAMINSKII, YL
    RUMYANTSEVA, LN
    EFIMOVA, VL
    STRADINSH, JP
    KHIMIYA GETEROTSIKLICHESKIKH SOEDINENII, 1992, (10): : 1404 - 1408