Calculation of adsorption free energy for solute-surface interactions using biased replica-exchange molecular dynamics

The adsorption behavior of a biomolecule, such as a peptide or protein, to a functionalized surface is of fundamental importance for a broad range of applications in biotechnology. The adsorption free energy for these types of interactions can be determined from a molecular dynamics simulation using the partitioning between adsorbed and nonadsorbed states, provided that sufficient sampling of both states is obtained. However, if interactions between the solute and the surface are strong, the solute will tend to be trapped near the surface during the simulation, thus preventing the adsorption free energy from being calculated by this method. This situation occurs even when using an advanced sampling algorithm such as replica-exchange molecular dynamics (REMD). In this paper, the authors demonstrate the fundamental basis of this problem using a model system consisting of one sodium ion (Na+) as the solute positioned over a surface functionalized with one negatively charged group (COO-) in explicit water. With this simple system, the authors show that sufficient sampling in the coordinate normal to the surface cannot be obtained by conventional REMD alone. The authors then present a method to overcome this problem through the use of an adaptive windowed-umbrella sampling technique to develop a biased-energy function that is combined with REMD. This approach provides an effective method for the calculation of adsorption free energy for solute-surface interactions.

[1]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[2]  Richard Lavery,et al.  Base pair opening within B-DNA: free energy pathways for GC and AT pairs from umbrella sampling simulations. , 2003, Nucleic acids research.

[3]  M. Dobrovolskaia,et al.  Immunological properties of engineered nanomaterials , 2007, Nature Nanotechnology.

[4]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[5]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[6]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[7]  Cezary Czaplewski,et al.  Replica Exchange and Multicanonical Algorithms with the coarse-grained UNRES force field. , 2006, Journal of chemical theory and computation.

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

[9]  J. Skolnick,et al.  Ab initio protein structure prediction via a combination of threading, lattice folding, clustering, and structure refinement , 2001, Proteins.

[10]  Clive D McFarland,et al.  Platelet adhesion to polystyrene-based surfaces preadsorbed with plasmas selectively depleted in fibrinogen, fibronectin, vitronectin, or von Willebrand's factor. , 2002, Journal of biomedical materials research.

[11]  Kerstin Länge,et al.  Chemical modification of parylene C coatings for SAW biosensors , 2007 .

[12]  T. Horbett,et al.  Passivating protein coatings for implantable glucose sensors: evaluation of protein retention. , 2007, Journal of biomedical materials research. Part B, Applied biomaterials.

[13]  J. Skolnick,et al.  A new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics , 2001 .

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

[15]  J. Overbeek,et al.  Recent developments in the understanding of colloid stability , 1977 .

[16]  V. Vacek,et al.  REMC computer simulations of the thermodynamic properties of argon and air plasmas , 2002 .

[17]  Umbrella sampling along linear combinations of generalized coordinates. Theory and application to a glycine dipeptide , 1995 .

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

[19]  Robert A Latour,et al.  Molecular dynamics simulations of peptide-surface interactions. , 2005, Langmuir : the ACS journal of surfaces and colloids.

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

[21]  Joshua A. Kritzer,et al.  Relationship between side chain structure and 14-helix stability of beta3-peptides in water. , 2005, Journal of the American Chemical Society.

[22]  F. Hartwig,et al.  Exploratory Data Analysis , 2008, Using Science in Cybersecurity.

[23]  Edward Lyman,et al.  A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation. , 2006, Journal of chemical theory and computation.

[24]  Christian Bartels,et al.  Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations , 1997 .

[25]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[26]  John W. Tukey,et al.  Exploratory Data Analysis. , 1979 .

[27]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[28]  Jeffrey Skolnick,et al.  All-atom ab initio folding of a diverse set of proteins. , 2006, Structure.

[29]  Y. Sugita,et al.  Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. I. Formulation and benchmark test , 2003 .

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

[31]  Martin Karplus,et al.  Probability Distributions for Complex Systems: Adaptive Umbrella Sampling of the Potential Energy , 1998 .

[32]  Christian Bartels,et al.  Determination of equilibrium properties of biomolecular systems using multidimensional adaptive umbrella sampling , 1999 .

[33]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[34]  B. Dominy,et al.  Development of a generalized Born model parameterization for proteins and nucleic acids , 1999 .

[35]  Eugene I Shakhnovich,et al.  A knowledge‐based move set for protein folding , 2007, Proteins.

[36]  V. Hlady,et al.  Protein adsorption on solid surfaces. , 1996, Current opinion in biotechnology.

[37]  M. Karplus,et al.  Adaptive umbrella sampling of the potential energy: modified updating procedure of the umbrella potential and application to peptide folding , 1999 .

[38]  Fernando A Escobedo,et al.  A general framework for non-Boltzmann Monte Carlo sampling. , 2006, The Journal of chemical physics.

[39]  K. Sanbonmatsu,et al.  Exploring the energy landscape of a β hairpin in explicit solvent , 2001 .

[40]  M. Ricci,et al.  Solid lipid nanoparticles for targeted brain drug delivery. , 2007, Advanced drug delivery reviews.

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

[42]  I. Tavernelli,et al.  A Novel Hamiltonian Replica Exchange MD Protocol to Enhance Protein Conformational Space Sampling. , 2006, Journal of chemical theory and computation.

[43]  Y. Okamoto,et al.  Prediction of membrane protein structures by replica-exchange Monte Carlo simulations: case of two helices. , 2004, The Journal of chemical physics.

[44]  A. Amanda.  Semicrystalline poly(vinyl alcohol) ultrafiltration membranes for bioseparations , 2000 .

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

[46]  Holger H. Hoos,et al.  A replica exchange Monte Carlo algorithm for protein folding in the HP model , 2007, BMC Bioinformatics.

[47]  Anthony K. Felts,et al.  Temperature weighted histogram analysis method, replica exchange, and transition paths. , 2005, The journal of physical chemistry. B.

[48]  Hartmut Löwen,et al.  Umbrella sampling in non-equilibrium computer simulations , 2005, Comput. Phys. Commun..

[49]  D. Castner,et al.  Biomedical surface science: Foundations to frontiers , 2002 .

[50]  Mihaly Mezei,et al.  The Potentials of Mean Force of Sodium Chloride and Sodium Dimethylphosphate in Water : An Application of Adaptive Umbrella Sampling , 1995 .

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

[52]  Robert A. Latour,et al.  Adsorption Thermodynamics Of A Mid-Chain Peptide Residue On Functionalized SAM Surfaces Using SPR , 2005 .