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
Nature of Glycine and Its α-Carbon Radical in Aqueous Solution: A Theoretical Investigation.
TL;DR: Analysis of the microsolvation environments of each species provides evidence that the change in preference between glycine and glycyl radical is due to the inherent gas-phase stability of the neutral α-carbon radical rather than to any significant difference in the solvation behavior of the constituent species.
26
Long range recognition and selection in IDPs: the interactions of the C-terminus of p53
TL;DR: The C-terminal domain of p53 is an extensively studied IDP, interacting with different partners through multiple distinct conformations, and the free peptide segment rapidly interconverts between ordered and disordered states with significant populations of the conformations that are seen in the complexed states.
HIV-1 Protease and Substrate Coevolution Validates the Substrate Envelope As the Substrate Recognition Pattern.
TL;DR: This study on the conformational and mutational ensembles of protease-substrate complexes validates the substrate envelope as the substrate recognition motif for HIV-1 protease, and allows for the elucidation of possible drug resistance mutation patterns in the polyprotein cleavage sites.
How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics.
TL;DR: The tunnel region at triosephosphate isomerase (TIM)'s dimer interface is revealed as a key binding site that can alter global dynamics of the enzyme and allosteric changes in the catalytic loop's dynamics, flexibility, and correlations are observed.
26
A Candidate Ion-Retaining State in the Inward-Facing Conformation of Sodium/Galactose Symporter: Clues from Atomistic Simulations
TL;DR: This work has identified a more stable Na(+) binding site corresponding to a putative ion-retaining state of the transporter and highlights the importance of D189 that, without being one of the Na(+)-coordinating residues, regulates its binding/release.
26
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).