CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model

A fluctuating charge (FQ) force field is applied to molecular dynamics simulations for six small proteins in explicit polarizable solvent represented by the TIP4P‐FQ potential. The proteins include 1FSV, 1ENH, 1PGB, 1VII, 1H8K, and 1CRN, representing both helical and β‐sheet secondary structural elements. Constant pressure and temperature (NPT) molecular dynamics simulations are performed on time scales of several nanoseconds, the longest simulations yet reported using explicitly polarizable all‐atom empirical potentials (for both solvent and protein) in the condensed phase. In terms of structure, the FQ force field allows deviations from native structure up to 2.5 Å (with a range of 1.0 to 2.5 Å). This is commensurate to the performance of the CHARMM22 nonpolarizable model and other currently existing polarizable models. Importantly, secondary structural elements maintain native structure in general to within 1 Å (both helix and β‐strands), again in good agreement with the nonpolarizable case. In qualitative agreement with QM/MM ab initio dynamics on crambin (Liu et al. Proteins 2001, 44, 484), there is a sequence dependence of average condensed phase atomic charge for all proteins, a dependence one would anticipate considering the differing chemical environments around individual atoms; this is a subtle quantum mechanical feature captured in the FQ model but absent in current state‐of‐the‐art nonpolarizable models. Furthermore, there is a mutual polarization of solvent and protein in the condensed phase. Solvent dipole moment distributions within the first and second solvation shells around the protein display a shift towards higher dipole moments (increases on the order of 0.2–0.3 Debye) relative to the bulk; protein polarization is manifested via the enhanced condensed phase charges of typical polar atoms such as backbone carbonyl oxygens, amide nitrogens, and amide hydrogens. Finally, to enlarge the sample set of proteins, gas‐phase minimizations and 1 ps constant temperature simulations are performed on various‐sized proteins to compare to earlier work by Kaminsky et al. (J Comp Chem 2002, 23, 1515). The present work establishes the feasibility of applying a fully polarizable force field for protein simulations and demonstrates the approach employed in extending the CHARMM force field to include these effects. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1504–1514, 2004

[1]  Jacek Korchowiec,et al.  Molecular hardness and softness parameters and their use in chemistry , 1988 .

[2]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[3]  Steven J. Stuart,et al.  Fluctuating charge force fields for aqueous solutions , 1995 .

[4]  Nohad Gresh,et al.  Energetics of Zn2+ binding to a series of biologically relevant ligands: A molecular mechanics investigation grounded on ab initio self‐consistent field supermolecular computations , 1995, J. Comput. Chem..

[5]  Douglas J. Tobias,et al.  Ions at the Air/Water Interface , 2002 .

[6]  Gerhard Hummer,et al.  Simulation and Theory of Electrostatic Interactions in Solution: Computational Chemistry, Biophysics and Aqueous Solutions, Santa Fe, New Mexico, U. S. A., 23-25 June 1999 , 1999 .

[7]  Darrin M. York,et al.  A chemical potential equalization method for molecular simulations , 1996 .

[8]  Steven J. Stuart,et al.  Effects of Polarizability on the Hydration of the Chloride Ion , 1996 .

[9]  Akihiro Morita.  Water polarizability in condensed phase: Ab initio evaluation by cluster approach , 2002, J. Comput. Chem..

[10]  Pastore,et al.  Theory of ab initio molecular-dynamics calculations. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[11]  Noriyuki Yoshii,et al.  A molecular dynamics study of sub- and supercritical water using a polarizable potential model , 1998 .

[12]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

[13]  Marimuthu Krishnan,et al.  Computer simulation study of water using a fluctuating charge model , 2001 .

[14]  R. T. Sanderson,et al.  An Interpretation of Bond Lengths and a Classification of Bonds. , 1951, Science.

[15]  C. Brooks,et al.  A molecular dynamics simulation study of segment B1 of protein G , 1997, Proteins.

[16]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[17]  Koji Ando,et al.  A stable fluctuating-charge polarizable model for molecular dynamics simulations: Application to aqueous electron transfers , 2001 .

[18]  Harry A. Stern,et al.  Development of a polarizable force field for proteins via ab initio quantum chemistry: First generation model and gas phase tests , 2002, J. Comput. Chem..

[19]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[20]  Serdar Kuyucak,et al.  Gramicidin A channel as a test ground for molecular dynamics force fields. , 2003, Biophysical journal.

[21]  Nohad Gresh,et al.  Comparative binding energetics of Mg2+, Ca2+, Zn2+, and Cd2+ to biologically relevant ligands: Combined ab initio SCF supermolecule and molecular mechanics investigation , 1996, J. Comput. Chem..

[22]  Charles L. Brooks,et al.  An abinitio study of hydrated chloride ion complexes: Evidence of polarization effects and nonadditivity , 1987 .

[23]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[24]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[25]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[26]  L. Dang,et al.  The nonadditive intermolecular potential for water revised , 1992 .

[27]  D. R. Garmer,et al.  Modeling of inhibitor–metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry , 1998, Proteins.

[28]  L. E. Chirlian,et al.  Atomic charges derived from electrostatic potentials: A detailed study , 1987 .

[29]  Noriyuki Yoshii,et al.  A molecular dynamics study of dielectric constant of water from ambient to sub- and supercritical conditions using a fluctuating-charge potential model , 2001 .

[30]  Alexander D. MacKerell,et al.  Computational Biochemistry and Biophysics , 2001 .

[31]  Charles L. Brooks,et al.  CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations , 2004, J. Comput. Chem..

[32]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[33]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[34]  Harry A. Stern,et al.  Fluctuating Charge, Polarizable Dipole, and Combined Models: Parameterization from ab Initio Quantum Chemistry , 1999 .

[35]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[36]  L. Dang,et al.  A Mechanism for Ion Transport across the Water/Dichloromethane Interface: A Molecular Dynamics Study Using Polarizable Potential Models , 2001 .

[37]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[38]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[39]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[40]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[41]  L. Dang,et al.  MOLECULAR DYNAMICS STUDY OF WATER CLUSTERS, LIQUID, AND LIQUID-VAPOR INTERFACE OF WATER WITH MANY-BODY POTENTIALS , 1997 .

[42]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[43]  C. Brooks.  Computer simulation of liquids , 1989 .

[44]  Steven J. Stuart,et al.  Surface Curvature Effects in the Aqueous Ionic Solvation of the Chloride Ion , 1999 .

[45]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[46]  Michiel Sprik,et al.  Solvent polarization and hydration of the chlorine anion , 1990 .

[47]  Kirk A. Peterson,et al.  Computer Simulation of Chloroform with a Polarizable Potential Model , 1997 .

[48]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[49]  S. L. Mayo,et al.  De novo protein design: fully automated sequence selection. , 1997, Science.

[50]  L. Dang,et al.  Many-body interactions in liquid methanol and its liquid/vapor interface: A molecular dynamics study , 2003 .

[51]  David van der Spoel,et al.  Molecular Dynamics Simulations of Water with Novel Shell-Model Potentials , 2001 .

[52]  B. Berne,et al.  Free Energy of the Hydrophobic Interaction from Molecular Dynamics Simulations: The Effects of Solute and Solvent Polarizability , 1997 .

[53]  P. Kollman,et al.  Structure and Properties of Neat Liquids Using Nonadditive Molecular Dynamics: Water, Methanol, and N-Methylacetamide , 1995 .

[54]  S. Rick.  Simulations of ice and liquid water over a range of temperatures using the fluctuating charge model , 2001 .

[55]  R. T. Sanderson.  Chemical Bonds and Bond Energy , 1976 .

[56]  Donald G Truhlar,et al.  Importance of substrate and cofactor polarization in the active site of dihydrofolate reductase. , 2003, Journal of molecular biology.

[57]  Ruhong Zhou,et al.  Parametrizing a polarizable force field from ab initio data. I. The fluctuating point charge model , 1999 .

[58]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[59]  Mauro C. C. Ribeiro,et al.  Fluctuating charge model for polyatomic ionic systems: A test case with diatomic anions , 1999 .

[60]  Alexander D. MacKerell,et al.  Force field influence on the observation of π-helical protein structures in molecular dynamics simulations , 2003 .

[61]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[62]  J. Ilja Siepmann,et al.  Development of Polarizable Water Force Fields for Phase Equilibrium Calculations , 2000 .

[63]  E. Baker,et al.  Hydrogen bonding in globular proteins. , 1984, Progress in biophysics and molecular biology.

[64]  F. Artzner,et al.  THERMOTROPIC PHASE BEHAVIOR OF CATIONIC LIPID: DNA COMPLEXES COMPARED TO BINARY LIPID MIXTURES , 1999 .

[65]  Koji Ando,et al.  Fluctuating Charge Study of Polarization Effects in Chlorinated Organic Liquids , 2001 .

[66]  O. Olsen,et al.  A simple and realistic model system for studying hydrogen bonds in β-sheets , 2003 .

[67]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[68]  M Elstner,et al.  Quantum mechanics simulation of protein dynamics on long timescale , 2001, Proteins.

[69]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[70]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[71]  M. Parrinello,et al.  Crystal structure and pair potentials: A molecular-dynamics study , 1980 .

[72]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

[73]  B. Berne,et al.  Combined fluctuating charge and polarizable dipole models: Application to a five-site water potential function , 2001 .

[74]  Richard A. Friesner,et al.  Accurate ab Initio Quantum Chemical Determination of the Relative Energetics of Peptide Conformations and Assessment of Empirical Force Fields , 1997 .

[75]  Bruce J. Berne,et al.  Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides , 1996 .

[76]  Nathaniel O. J. Malcolm,et al.  Cooperative effects in the structuring of fluoride water clusters: Ab initio hybrid quantum mechanical/molecular mechanical model incorporating polarizable fluctuating charge solvent , 1998 .

[77]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[78]  Michiel Sprik,et al.  Hydrogen bonding and the static dielectric constant in liquid water , 1991 .

[79]  Michiel Sprik,et al.  A polarizable model for water using distributed charge sites , 1988 .

[80]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .