Development of molecular simulation methods to accurately represent protein-surface interactions: The effect of pressure and its determination for a system with constrained atoms

When performing molecular dynamics simulations for a system with constrained (fixed) atoms, traditional isobaric algorithms (e.g., NPT simulation) often cannot be used. In addition, the calculation of the internal pressure of a system with fixed atoms may be highly inaccurate due to the nonphysical nature of the atomic constraints and difficulties in accurately defining the volume occupied by the unconstrained atoms in the system. The inability to properly set and control pressure can result in substantial problems for the accurate simulation of condensed-phase systems if the behavior of the system (e.g., peptide/protein adsorption) is sensitive to pressure. To address this issue, the authors have developed an approach to accurately determine the internal pressure for a system with constrained atoms. As the first step in this method, a periodically extendable portion of the mobile phase of the constrained system (e.g., the solvent atoms) is used to create a separate unconstrained system for which the pressure can be accurately calculated. This model system is then used to create a pressure calibration plot for an intensive local effective virial parameter for a small volume cross section or “slab” of the system. Using this calibration plot, the pressure of the constrained system can then be determined by calculating the virial parameter for a similarly sized slab of mobile atoms. In this article, the authors present the development of this method and demonstrate its application using the CHARMM molecular simulation program to characterize the adsorption behavior of a peptide in explicit water on a hydrophobic surface whose lattice spacing is maintained with atomic constraints. The free energy of adsorption for this system is shown to be dramatically influenced by pressure, thus emphasizing the importance of properly maintaining the pressure of the system for the accurate simulation of protein-surface interactions.

[1]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[2]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[3]  R. Winkler.  Virial pressure of periodic systems with long range forces , 2002 .

[4]  S J Tendler,et al.  Surface plasmon resonance analysis of dynamic biological interactions with biomaterials. , 2000, Biomaterials.

[5]  M. Klein,et al.  Molecular dynamics simulations of a hydrated protein vectorially oriented on polar and nonpolar soft surfaces. , 2002, Biophysical journal.

[6]  R. Netz.  Water and ions at interfaces , 2004 .

[7]  Robert A Latour,et al.  Molecular simulation to characterize the adsorption behavior of a fibrinogen gamma-chain fragment. , 2005, Langmuir : the ACS journal of surfaces and colloids.

[8]  D A Bruce,et al.  Modeling of peptide adsorption interactions with a poly(lactic acid) surface. , 2008, Langmuir : the ACS journal of surfaces and colloids.

[9]  Nadeem A. Vellore,et al.  An improved replica-exchange sampling method: temperature intervals with global energy reassignment. , 2007, The Journal of chemical physics.

[10]  D. H. Tsai.  The virial theorem and stress calculation in molecular dynamics , 1979 .

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

[12]  M. Grunze,et al.  Grand Canonical Monte Carlo Simulations of the Hydration Interaction between Oligo(ethylene glycol)-Terminated Alkanethiol Self-Assembled Monolayers , 2002 .

[13]  Roberto C. Salvarezza,et al.  Surface characterization of sulfur and alkanethiol self-assembled monolayers on Au(111) , 2006 .

[14]  E. Paci.  High pressure simulations of biomolecules. , 2002, Biochimica et biophysica acta.

[15]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[16]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[17]  Robert A Latour,et al.  A molecular modeling study of the effect of surface chemistry on the adsorption of a fibronectin fragment spanning the 7-10th type III repeats. , 2004, Journal of biomedical materials research. Part A.

[18]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[19]  A. Leach.  Molecular Modelling: Principles and Applications , 1996 .

[20]  Yu Sun,et al.  Comparison of solvation‐effect methods for the simulation of peptide interactions with a hydrophobic surface , 2007, J. Comput. Chem..

[21]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .

[22]  Tomohiro Hayashi,et al.  Interaction of water with self-assembled monolayers: Neutron reflectivity measurements of the water density in the interface region , 2003 .

[23]  Robert A Latour,et al.  Benchmark experimental data set and assessment of adsorption free energy for peptide-surface interactions. , 2009, Langmuir : the ACS journal of surfaces and colloids.

[24]  W. G. McMillan,et al.  The Virial Theorem , 2007 .

[25]  H. Heinz,et al.  Calculation of local and average pressure tensors in molecular simulations , 2007 .

[26]  A. Wallqvist.  Pressure dependence of methane solvation in aqueous mixtures and the relation to the structure of liquid water , 1992 .

[27]  Free Energy, Entropy and Heat Capacity of the Hydrophobic Interaction as a Function of Pressure , 2000 .

[28]  On the Equivalence of Atomic and Molecular Pressure , 2004 .

[29]  Kurt Binder,et al.  Calculation of local pressure tensors in systems with many-body interactions. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  M. Grunze,et al.  Computer Simulation of Water near the Surface of Oligo(ethylene glycol)-Terminated Alkanethiol Self-Assembled Monolayers† , 2000 .

[31]  D. Y. Yoon,et al.  Novel molecular dynamics simulations at constant pressure , 1992 .

[32]  Robert A. Latour,et al.  Molecular modeling of biomaterial surfaces , 1999 .

[33]  G. Hummer,et al.  The pressure dependence of hydrophobic interactions is consistent with the observed pressure denaturation of proteins. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

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

[35]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[36]  Baldomero Oliva,et al.  Calculation of the group-based pressure in molecular simulations. II. Numerical tests and application to liquid water , 2002 .

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

[38]  R. Winkler,et al.  Liquid benzene confined between graphite surfaces. A constant pressure molecular dynamics study , 1993 .

[39]  Aidan P Thompson,et al.  General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions. , 2009, The Journal of chemical physics.

[40]  Evert Jan Baerends,et al.  Calculation of pressure in case of periodic boundary conditions , 2006 .

[41]  G. Wnek,et al.  Encyclopedia of biomaterials and biomedical engineering , 2008 .

[42]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[43]  Feng Wang,et al.  Calculation of adsorption free energy for solute-surface interactions using biased replica-exchange molecular dynamics , 2008, Biointerphases.

[44]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[45]  Nadeem A. Vellore,et al.  Assessment of the transferability of a protein force field for the simulation of peptide-surface interactions. , 2010, Langmuir : the ACS journal of surfaces and colloids.

[46]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[47]  H. Berendsen,et al.  AN EFFICIENT, BOX SHAPE INDEPENDENT NONBONDED FORCE AND VIRIAL ALGORITHM FOR MOLECULAR-DYNAMICS , 1995 .

[48]  Nadeem A. Vellore,et al.  Development of molecular simulation methods to accurately represent protein-surface interactions: Method assessment for the calculation of electrostatic effects , 2009, Biointerphases.

[49]  Philippe H. Hünenberger,et al.  Calculation of the group-based pressure in molecular simulations. I. A general formulation including Ewald and particle-particle-particle-mesh electrostatics , 2002 .

[50]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

[51]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[52]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[53]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[54]  Larry L Hench,et al.  A theoretical analysis of the thermodynamic contributions for the adsorption of individual protein residues on functionalized surfaces. , 2002, Biomaterials.

[55]  Robert A Latour,et al.  Determination of the adsorption free energy for peptide-surface interactions by SPR spectroscopy. , 2008, Langmuir : the ACS journal of surfaces and colloids.

[56]  W. L. Jorgensen.  Revised TIPS for simulations of liquid water and aqueous solutions , 1982 .

[57]  S. Garde,et al.  Molecular dynamics simulations of pressure effects on hydrophobic interactions. , 2001, Journal of the American Chemical Society.

[58]  J. Erpenbeck,et al.  Molecular Dynamics Techniques for Hard-Core Systems , 1977 .

[59]  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..

[60]  S. Neyertz,et al.  A GENERAL PRESSURE TENSOR CALCULATION FOR MOLECULAR DYNAMICS SIMULATIONS , 1995 .

[61]  Yu Sun,et al.  Comparison of implicit solvent models for the simulation of protein–surface interactions , 2006, J. Comput. Chem..