Improving the Performance of MM/PBSA in Protein-Protein Interactions via the Screening Electrostatic Energy

Accurate calculation of protein-protein binding free energy is of great importance in biological and medical science, yet it remains a hugely challenging problem. In this work, we develop a new strategy in which a screened electrostatic energy (i.e., adding an exponential damping factor to the Coulombic interaction energy) is used within the framework of the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method. Our results show that the Pearson correlation coefficient in the modified MM/PBSA is over 0.70, which is much better than that in the standard MM/PBSA, especially in the Amber14SB force field. In particular, the performance of the standard MM/PBSA is very poor in a system where the proteins carry like charges. Moreover, we also calculated the mean absolute error (MAE) between the calculated and experimental ΔG values and found that the MAE in the modified MM/PBSA was indeed much smaller than that in the standard MM/PBSA. Furthermore, the effect of the dielectric constant of the proteins and the salt conditions on the results was also investigated. The present study highlights the potential power of the modified MM/PBSA for accurately predicting the binding energy in highly charged biosystems.

[1]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[2]  W. Sherman,et al.  Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. , 2010, Journal of chemical theory and computation.

[3]  Huan‐Xiang Zhou,et al.  Fundamental aspects of protein-protein association kinetics. , 2009, Chemical reviews.

[4]  Gail J. Bartlett,et al.  Predicting and Experimentally Validating Hot-Spot Residues at Protein–Protein Interfaces , 2019, ACS chemical biology.

[5]  Samuel Genheden,et al.  Effect of explicit water molecules on ligand-binding affinities calculated with the MM/GBSA approach , 2014, Journal of Molecular Modeling.

[6]  J. Homola.  Surface plasmon resonance sensors for detection of chemical and biological species. , 2008, Chemical reviews.

[7]  Bing Xie,et al.  On Restraints in End‐Point Protein–Ligand Binding Free Energy Calculations , 2019, J. Comput. Chem..

[8]  Ray Luo,et al.  Virtual screening using molecular simulations , 2011, Proteins.

[9]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  W. Rocchia,et al.  dMM-PBSA: A New HADDOCK Scoring Function for Protein-Peptide Docking , 2016, Front. Mol. Biosci..

[12]  Mathias Jucker,et al.  The Amyloid State of Proteins in Human Diseases , 2012, Cell.

[13]  H. Gohlke,et al.  Free Energy Calculations by the Molecular Mechanics Poisson−Boltzmann Surface Area Method , 2012, Molecular informatics.

[14]  T. Straatsma,et al.  Inversion of receptor binding preferences by mutagenesis: free energy thermodynamic integration studies on sugar binding to L-arabinose binding proteins. , 1993, Biochemistry.

[15]  Ye Mei,et al.  Electrostatic polarization makes a substantial contribution to the free energy of avidin-biotin binding. , 2010, Journal of the American Chemical Society.

[16]  Junmei Wang,et al.  End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. , 2019, Chemical reviews.

[17]  Yu Kang,et al.  Development and Evaluation of MM/GBSA Based on a Variable Dielectric GB Model for Predicting Protein-Ligand Binding Affinities , 2020, J. Chem. Inf. Model..

[18]  John Z. H. Zhang,et al.  Interaction Entropy: A New Paradigm for Highly Efficient and Reliable Computation of Protein-Ligand Binding Free Energy. , 2016, Journal of the American Chemical Society.

[19]  Alessandro Contini,et al.  Explicit Ligand Hydration Shells Improve the Correlation between MM-PB/GBSA Binding Energies and Experimental Activities. , 2013, Journal of chemical theory and computation.

[20]  J. H. Zhang,et al.  Interaction entropy for protein-protein binding. , 2017, The Journal of chemical physics.

[21]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[22]  David D. L. Minh,et al.  Simple Entropy Terms for End-Point Binding Free Energy Calculations. , 2018, Journal of chemical theory and computation.

[23]  Edward D Harder,et al.  Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems , 2001 .

[24]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

[25]  Weiliang Zhu,et al.  Improving the accuracy of predicting protein-ligand binding-free energy with semiempirical quantum chemistry charge. , 2019, Future medicinal chemistry.

[26]  T. Straatsma,et al.  Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration , 1994 .

[27]  Giulio Rastelli,et al.  Fast and accurate predictions of binding free energies using MM‐PBSA and MM‐GBSA , 2009, J. Comput. Chem..

[28]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[29]  Alessandro Contini,et al.  Improved Computation of Protein-Protein Relative Binding Energies with the Nwat-MMGBSA Method , 2016, J. Chem. Inf. Model..

[30]  W. Stites,et al.  Protein−Protein Interactions: Interface Structure, Binding Thermodynamics, and Mutational Analysis , 1997 .

[31]  Ray Luo,et al.  Calculating protein–ligand binding affinities with MMPBSA: Method and error analysis , 2016, J. Comput. Chem..

[32]  Harel Weinstein,et al.  Thermodynamic basis for promiscuity and selectivity in protein-protein interactions: PDZ domains, a case study. , 2006, Journal of the American Chemical Society.

[33]  William L. Jorgensen,et al.  Quantum and statistical mechanical studies of liquids. 25. Solvation and conformation of methanol in water , 1983 .

[34]  Xiao Liu,et al.  Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants. , 2018, Journal of chemical information and modeling.

[35]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[36]  Xiao Hu,et al.  Rescoring Virtual Screening Results with the MM-PBSA Methods: Beware of Internal Dielectric Constants , 2019, J. Chem. Inf. Model..

[37]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[38]  F. Kirchhoff,et al.  An optimized MM/PBSA virtual screening approach applied to an HIV‐1 gp41 fusion peptide inhibitor , 2011, Proteins.

[39]  S. Jones,et al.  Principles of protein-protein interactions. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[40]  H Videler,et al.  Protein-protein interactions in colicin E9 DNase-immunity protein complexes. 2. Cognate and noncognate interactions that span the millimolar to femtomolar affinity range. , 1995, Biochemistry.

[41]  Christian Kramer,et al.  MM/GBSA Binding Energy Prediction on the PDBbind Data Set: Successes, Failures, and Directions for Further Improvement , 2013, J. Chem. Inf. Model..

[42]  Y. Levin,et al.  Electrostatic correlations: from plasma to biology , 2002 .

[43]  A. Hopkins,et al.  Emerging role of surface plasmon resonance in fragment-based drug discovery. , 2011, Future medicinal chemistry.

[44]  William L Jorgensen,et al.  Perspective on Free-Energy Perturbation Calculations for Chemical Equilibria. , 2008, Journal of chemical theory and computation.

[45]  Samuel Genheden,et al.  How to obtain statistically converged MM/GBSA results , 2009, J. Comput. Chem..

[46]  Ray Luo,et al.  Recent Developments and Applications of the MMPBSA Method , 2018, Front. Mol. Biosci..

[47]  B. L. de Groot,et al.  CHARMM36m: an improved force field for folded and intrinsically disordered proteins , 2016, Nature Methods.

[48]  Chris-Kriton Skylaris,et al.  Protein-protein interactions from linear-scaling first-principles quantum-mechanical calculations , 2010 .

[49]  F. Khuri,et al.  Targeting protein-protein interactions as an anticancer strategy. , 2013, Trends in pharmacological sciences.

[50]  Holger Gohlke,et al.  Efficient Approximation of Ligand Rotational and Translational Entropy Changes upon Binding for Use in MM-PBSA Calculations , 2017, J. Chem. Inf. Model..

[51]  Carlos Simmerling,et al.  Improved Generalized Born Solvent Model Parameters for Protein Simulations. , 2013, Journal of chemical theory and computation.

[52]  Samuel Genheden,et al.  A semiempirical approach to ligand‐binding affinities: Dependence on the Hamiltonian and corrections , 2012, J. Comput. Chem..

[53]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[54]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[55]  Jianing Song,et al.  Calculation of hot spots for protein–protein interaction in p53/PMI‐MDM2/MDMX complexes , 2018, J. Comput. Chem..

[56]  Luigi F Agnati,et al.  Fluorescence resonance energy transfer-based technologies in the study of protein-protein interactions at the cell surface. , 2012, Methods.