Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations

A point‐charge force field for molecular mechanics simulations of proteins based on... 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 ε = 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 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of Computational Chemistry Wiley

A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations

Download PDF
14 pages

Loading next page...
 
/lp/wiley/a-point-charge-force-field-for-molecular-mechanics-simulations-of-c6iGmjornO

References (77)

Publisher
Wiley
Copyright
Copyright © 2003 Wiley Periodicals, Inc.
ISSN
0192-8651
eISSN
1096-987X
DOI
10.1002/jcc.10349
pmid
14531054
Publisher site
See Article on Publisher Site

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 ε = 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

Journal

Journal of Computational Chemistry – Wiley

Published: Dec 1, 2003

There are no references for this article.