Dependency of ligand free energy landscapes on charge parameters and solvent models

We performed replica-exchange molecular dynamics (REMD) simulations of six ligands to examine the dependency of their free energy landscapes on charge parameters and solvent models. Six different charge parameter sets for each ligand were first generated by RESP and AM1-BCC methods using three different conformations independently. RESP charges showed some conformational dependency. On the other hand, AM1-BCC charges did not show conformational dependency and well reproduced the overall trend of RESP charges. The free energy landscapes obtained from the REMD simulations of ligands in vacuum, Generalized-Born (GB), and TIP3P solutions were then analyzed. We found that even small charge differences can produce qualitatively different landscapes in vacuum condition, but the differences tend to be much smaller under GB and TIP3P conditions. The simulations in the GB model well reproduced the landscapes in the TIP3P model using only a fraction of the computational cost. The protein-bound ligand conformations were rarely the global minimum states, but similar conformations were found to exist in aqueous solution without proteins in regions close to the global minimum, local minimum or intermediate states.

[1]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[2]  David L Mobley,et al.  Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent. , 2007, The journal of physical chemistry. B.

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

[4]  David A. Case,et al.  Harmonic and quasiharmonic descriptions of crambin , 1990 .

[5]  B Montgomery Pettitt,et al.  Molecular basis of the apparent near ideality of urea solutions. , 2007, Biophysical journal.

[6]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[7]  Kenneth M. Merz,et al.  Can we separate active from inactive conformations? , 2002, J. Comput. Aided Mol. Des..

[8]  P Argos,et al.  Optimal protocol and trajectory visualization for conformational searches of peptides and proteins. , 1992, Journal of molecular biology.

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

[10]  David L Mobley,et al.  Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. , 2009, Journal of chemical theory and computation.

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

[12]  Eric J. Martin,et al.  Conformational Sampling of Bioactive Molecules: A Comparative Study , 2007, J. Chem. Inf. Model..

[13]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[14]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[15]  Jörg Weiser,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999, J. Comput. Chem..

[16]  H. Fujitani,et al.  Massively parallel computation of absolute binding free energy with well-equilibrated states. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  P. Charifson,et al.  Conformational analysis of drug-like molecules bound to proteins: an extensive study of ligand reorganization upon binding. , 2004, Journal of medicinal chemistry.

[18]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[19]  B Montgomery Pettitt,et al.  Preferential solvation in urea solutions at different concentrations: properties from simulation studies. , 2007, The journal of physical chemistry. B.

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

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

[22]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

[23]  N. Go,et al.  Investigating protein dynamics in collective coordinate space. , 1999, Current opinion in structural biology.

[24]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[25]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[26]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. , 2001, Advanced drug delivery reviews.

[27]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[28]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[29]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[30]  Fumio Hirata,et al.  The effects of solvent on the conformation and the collective motions of protein: normal mode analysis and molecular dynamics simulations of melittin in water and in vacuum , 1991 .

[31]  David L Mobley,et al.  Treating entropy and conformational changes in implicit solvent simulations of small molecules. , 2008, The journal of physical chemistry. B.

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

[33]  David L. Mobley,et al.  Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations , 2010, J. Comput. Aided Mol. Des..

[34]  David L Mobley,et al.  Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site. , 2009, Journal of molecular biology.

[35]  Jonas Boström,et al.  Reproducing the conformations of protein-bound ligands: A critical evaluation of several popular conformational searching tools , 2001, J. Comput. Aided Mol. Des..

[36]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[37]  Nicolas Foloppe,et al.  Conformational Sampling of Druglike Molecules with MOE and Catalyst: Implications for Pharmacophore Modeling and Virtual Screening , 2008, J. Chem. Inf. Model..

[38]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

[39]  Jonas Boström,et al.  Assessing the performance of OMEGA with respect to retrieving bioactive conformations. , 2003, Journal of molecular graphics & modelling.

[40]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[41]  Benoît Roux,et al.  Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations. , 2008, The Journal of chemical physics.

[42]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.