Journal Article10.1002/JCC.10349
A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations.
Yong Duan,Chun Wu,Shibasish Chowdhury,Mathew C. Lee,Guoming Xiong,Wei Zhang,Rong Yang,Piotr Cieplak,Piotr Cieplak,Ray Luo,Tai-Sung Lee,Tai-Sung Lee,James W. Caldwell,Junmei Wang,Peter A. Kollman +14 more
TL;DR: A third‐generation point‐charge all‐atom force field for proteins is developed and initial tests on peptides demonstrated a high‐degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace‐Gly‐nme and Ace‐Ala‐Nme di‐peptides.
read more
Abstract: Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc-pVTZ//HF/6-31G** quantum mechanical methods. The main-chain torsion parameters were obtained by fitting to the energy profiles of Ace-Ala-Nme and Ace-Gly-Nme di-peptides calculated using MP2/cc-pVTZ//HF/6-31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of e = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed-phase simulations of proteins. Initial tests on peptides demonstrated a high-degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace-Gly-Nme and Ace-Ala-Nme di-peptides. Some highlights of our results include (1) well-preserved balance between the extended and helical region distributions, and (2) favorable type-II poly-proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand-binding calculations. Test simulations on a large set of proteins are also discussed. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1999–2012, 2003
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Citations
Energetics of bacterial photosynthesis.
TL;DR: The results of extensive numerical simulations and theoretical calculations of electronic transitions in the reaction center of Rhodobacter sphaeroides photosynthetic bacterium suggest that the experimentally observed break in the Arrhenius slope of the primary recombination rate can be traced back to a significant drop of the solvent reorganization energy close to that temperature.
Molecular Dynamics Simulations and Free Energy Analyses on the Dimer Formation of an Amyloidogenic Heptapeptide from Human β2-Microglobulin: Implication for the Protofibril Structure
TL;DR: It is proposed that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds, and a double-layer protofibril structure is also proposed in which two antipar parallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.
33
Modeling biominerals formed by apatites and DNA
Guillermo Revilla-López,Jordi Casanovas,Oscar Bertran,Pau Turon,Jordi Puiggalí,Carlos Alemán +5 more
TL;DR: Results indicate that duplexes of DNA adopting a B double helix can be encapsulated inside nanopores of hydroxyapatite without undergoing significant distortions in the inter-strand hydrogen bonds and the intra-Strand stacking.
Synthesis, biological evaluation and computational study of novel isoniazid containing 4H-Pyrimido[2,1-b]benzothiazoles derivatives
Manoj N. Bhoi,Mayuri A. Borad,Divya J. Jethava,Prachi T. Acharya,Edwin A. Pithawala,Chirag N. Patel,Himanshu A. Pandya,Hitesh D. Patel +7 more
TL;DR: The calculated ADMET descriptors for the synthesized compounds validated good pharmacokinetic properties, confirming that these compounds could be used as templates for the development of new Anti-mycobacterial agents.
32
An Integrated Pharmacophore/Docking/3D-QSAR Approach to Screening a Large Library of Products in Search of Future Botulinum Neurotoxin A Inhibitors.
Davide Gentile,Giuseppe Floresta,Vincenzo Patamia,Rita Chiaramonte,Giulia Letizia Mauro,Antonio Rescifina,Michele Vecchio +6 more
TL;DR: In this paper, 13,592,287 compounds were screened through a pharmacophore filter, a 3D-QSAR model, and a virtual screening; then, the compounds with the best affinity were selected and identified compounds with a calculated inhibitory activity in the range of 316–500 nM.
32
References
Density‐functional thermochemistry. III. The role of exact exchange
TL;DR: In this article, a semi-empirical exchange correlation functional with local spin density, gradient, and exact exchange terms was proposed. But this functional performed significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density
TL;DR: Numerical calculations on a number of atoms, positive ions, and molecules, of both open- and closed-shell type, show that density-functional formulas for the correlation energy and correlation potential give correlation energies within a few percent.
94.9K
Comparison of simple potential functions for simulating liquid water
TL;DR: In this article, the authors compared the Bernal Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P potential functions for liquid water in the NPT ensemble at 25°C and 1 atm.
39.4K
A smooth particle mesh Ewald method
TL;DR: It is demonstrated that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N), which is comparable to that of a simple truncation method of 10 A or less.
21.3K
Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions
TL;DR: In this paper, a reliable procedure for calculating the electron affinity of an atom and present results for hydrogen, boron, carbon, oxygen, and fluorine (hydrogen is included for completeness).