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
Stabilization of antibody structure upon association to a human carbonic anhydrase IX epitope studied by X-ray crystallography, microcalorimetry, and molecular dynamics simulations.
Kral,P. Mader,P. Mader,R Collard,Milan Fábry,M. Horejsi,Pavlina Rezacova,Milan Kozisek,Milan Kozisek,J. Zavada,Juraj Sedlacek,L Rulisek,Jiří Brynda +12 more
TL;DR: A correlation between structural elements and thermodynamic parameters of the association of antibody fragment Fab M75 to a peptide corresponding to its epitope in the proteoglycan‐like domain of CA IX, and it was observed that one proton is transferred from the environment to the protein‐ligand complex upon binding.
29
Molecular modeling and molecular dynamics simulation studies of the GSK3β/ATP/substrate complex: understanding the unique P+4 primed phosphorylation specificity for GSK3β substrates.
TL;DR: 3D model of the ternary complex of GSK3β, ATP, and the phosphorylated glycogen synthase (pGS) is constructed using a hierarchical approach and by integrating molecular modeling and molecular dynamics simulations, showing that for an optimal phosphorylation reaction to occur, the pGS priming phosphate in the GSK2β/ATP/pGS system optimizes the proper orientation of the G SK3β N- and C-terminal domains.
29
Functional and Structural Characterization of the Arylamine N-Acetyltransferase from the Opportunistic Pathogen Nocardia farcinica
Marta Martins,Benjamin Pluvinage,Inès Li de la Sierra-Gallay,Florent Barbault,Julien Dairou,Jean-Marie Dupret,Fernando Rodrigues-Lima +6 more
TL;DR: It is shown that this novel prokaryotic isoform is structurally and functionally related to the mycobacterial NAT enzymes, and was found to display high N-acetyltransferase activity towards several known NAT substrates including the antitubercular drug isoniazid.
29
The fast-folding HP35 double mutant has a substantially reduced primary folding free energy barrier
TL;DR: The LYS24/29NLE double mutant of villin headpiece subdomain (HP35) is the fastest folding protein known so far with a folding time constant of 0.6 micros, and unlike the proposed downhill folding, this work observed a well-defined folding transition that was consistent with experiments.
29
Molecular Dynamics Simulations and Structure-Based Rational Design Lead to Allosteric HCV NS5B Polymerase Thumb Pocket 2 Inhibitor with Picomolar Cellular Replicon Potency
Oliver Hucke,René Coulombe,Pierre R. Bonneau,Megan Bertrand-Laperle,Christian Brochu,James Gillard,Marc-André Joly,Serge Landry,Olivier Lepage,Montse Llinas-Brunet,Marc Pesant,Martin Poirier,Maude Poirier,Ginette McKercher,Martin Marquis,George Kukolj,Pierre L. Beaulieu,Timothy A. Stammers +17 more
TL;DR: This approach predicted a small but critical ligand-binding induced movement of a protein backbone region which increases the pocket size and improves access to the backbone carbonyl groups of Val 494 and Pro 495, and was subsequently confirmed by X-ray crystallography.
29
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).