Calculation of Standard Binding Free Energies:  Aromatic Molecules in the T4 Lysozyme L99A Mutant.

Calculations of the binding free energy of various nonpolar aromatic ligands with the L99A mutant of T4 lysozyme using molecular dynamics (MD) simulation are presented. To ensure better convergence, biasing potentials are used to restrain the ligand orientation and center-of-mass movement relative to the binding site when the ligand is decoupled from its environment in the binding pocket. The bias introduced by the restraint potentials is removed once the ligand fully interacts with the rest of the system and the calculated binding free energy is independent of the applied restraints. To decrease the computational cost, the simulations are generated with a reduced system in which protein and water atoms within a 15 Å-radius sphere around the ligand are included explicitly, while the rest of the system is treated with the generalized solvent boundary potential (GSBP). For all the ligands, the precision of the calculated free energy is less than 0.5 kcal/mol. For small nonpolar ligands such as benzene, toluene, and ethylbenzene, the calculated binding free energies are within 1.1 kcal/mol of the experimental values. For larger ligands, the computed binding free energies are slightly more favorable than the experimental values. The nonbinding polar molecule, phenol, has a calculated binding free energy of -0.88 kcal/mol. The simulation protocol presented here provides a way to calculate the binding free energy of small molecules to receptors at moderate computational cost.

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

[2]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[3]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[4]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[5]  A. Cheng,et al.  Design of RNA-binding proteins and ligands. , 2001, Current opinion in structural biology.

[6]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

[7]  B. Matthews,et al.  A model binding site for testing scoring functions in molecular docking. , 2002, Journal of molecular biology.

[8]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[9]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[10]  B. Montgomery Pettitt,et al.  MODELING SOLVENT IN BIOMOLECULAR SYSTEMS , 1994 .

[11]  Brian K Shoichet,et al.  Testing a flexible-receptor docking algorithm in a model binding site. , 2004, Journal of molecular biology.

[12]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[13]  B. Matthews,et al.  Specificity of ligand binding in a buried nonpolar cavity of T4 lysozyme: linkage of dynamics and structural plasticity. , 1995, Biochemistry.

[14]  K. Sharp,et al.  On the calculation of absolute macromolecular binding free energies , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Donald G. Truhlar,et al.  New Class IV Charge Model for Extracting Accurate Partial Charges from Wave Functions , 1998 .

[16]  B. G. Rao,et al.  Free Energy Calculations: Methods and Applications , 2002 .

[17]  M. Karplus,et al.  On the Treatment of Electrostatic Interactions in Biomolecular Simulations , 1991 .

[18]  W. Im,et al.  Generalized solvent boundary potential for computer simulations , 2001 .

[19]  Lu Wang,et al.  Inclusion of Loss of Translational and Rotational Freedom in Theoretical Estimates of Free Energies of Binding. Application to a Complex of Benzene and Mutant T4 Lysozyme , 1997 .

[20]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[21]  M Karplus,et al.  The energetics of off-rotamer protein side-chain conformations. , 2001, Journal of molecular biology.

[22]  John Marelius,et al.  Calculation of Ligand Binding Free Energies from Molecular Dynamics Simulations , 1998 .

[23]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

[24]  KumarShankar,et al.  The weighted histogram analysis method for free-energy calculations on biomolecules. I , 1992 .

[25]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[26]  J Hermans,et al.  Modeling protein-small molecule interactions: structure and thermodynamics of noble gases binding in a cavity in mutant phage T4 lysozyme L99A. , 2000, Journal of molecular biology.

[27]  A. Vindigni Energetic dissection of specificity in serine proteases. , 1999, Combinatorial chemistry & high throughput screening.

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

[29]  Alexander D. MacKerell,et al.  CH/π interactions involving aromatic amino acids: Refinement of the CHARMM tryptophan force field , 2005, J. Comput. Chem..

[30]  Michael K Gilson,et al.  On the theory of noncovalent binding. , 2004, Biophysical journal.

[31]  D. Beglov,et al.  Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations , 1994 .

[32]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[33]  B. Roux,et al.  Energetics of ion conduction through the gramicidin channel , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[34]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[35]  Roland L. Dunbrack,et al.  Conformational analysis of the backbone-dependent rotamer preferences of protein sidechains , 1994, Nature Structural Biology.

[36]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

[37]  Matthew Clark,et al.  Grand Canonical Monte Carlo Simulation of Ligand-Protein Binding , 2006, J. Chem. Inf. Model..

[38]  Emilio Gallicchio,et al.  On the nonpolar hydration free energy of proteins: surface area and continuum solvent models for the solute-solvent interaction energy. , 2003, Journal of the American Chemical Society.

[39]  B. Matthews,et al.  Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4 lysozyme. , 1995, Biochemistry.

[40]  J. Hermans,et al.  The Free Energy of Xenon Binding to Myoglobin from Molecular Dynamics Simulation , 1986 .

[41]  Saher Afshan Shaikh,et al.  A molecular thermodynamic view of DNA-drug interactions: a case study of 25 minor-groove binders. , 2004, Archives of biochemistry and biophysics.

[42]  C. Garvie,et al.  Recognition of specific DNA sequences. , 2001, Molecular cell.

[43]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[44]  P. Charifson,et al.  Are free energy calculations useful in practice? A comparison with rapid scoring functions for the p38 MAP kinase protein system. , 2001, Journal of medicinal chemistry.

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

[46]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[47]  G. Klebe,et al.  Approaches to the Description and Prediction of the Binding Affinity of Small-Molecule Ligands to Macromolecular Receptors , 2002 .

[48]  Benoît Roux,et al.  Electrostatic free energy calculations using the generalized solvent boundary potential method , 2002 .

[49]  Benoît Roux,et al.  Grand canonical Monte Carlo simulations of water in protein environments. , 2004, The Journal of chemical physics.

[50]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[51]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[52]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[53]  Jeremy C. Smith,et al.  Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study. , 1996, Biophysical journal.

[54]  Frederick W. Dahlquist,et al.  Studying excited states of proteins by NMR spectroscopy , 2001, Nature Structural Biology.

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

[56]  Benoît Roux,et al.  Hydration of Amino Acid Side Chains: Nonpolar and Electrostatic Contributions Calculated from Staged Molecular Dynamics Free Energy Simulations with Explicit Water Molecules , 2004 .

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

[58]  William L. Jorgensen Interactions between amides in solution and the thermodynamics of weak binding , 1989 .

[59]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[60]  K. Dill,et al.  On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. , 2006, The Journal of chemical physics.

[61]  Christophe Chipot,et al.  Can Absolute Free Energies of Association Be Estimated from Molecular Mechanical Simulations? The Biotin−Streptavidin System Revisited , 2001 .

[62]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[63]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.