The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4–2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4–1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8–2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.

[1]  Julian Tirado-Rives,et al.  Improving MM-GB/SA Scoring through the Application of the Variable Dielectric Model. , 2011, Journal of chemical theory and computation.

[2]  Philip C Biggin,et al.  Statistical Analysis on the Performance of Molecular Mechanics Poisson–Boltzmann Surface Area versus Absolute Binding Free Energy Calculations: Bromodomains as a Case Study , 2017, J. Chem. Inf. Model..

[3]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. , 2014, Physical chemistry chemical physics : PCCP.

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

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

[6]  J. H. Zhang,et al.  Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations , 2017, Scientific Reports.

[7]  P A Kollman,et al.  Free energy calculations on dimer stability of the HIV protease using molecular dynamics and a continuum solvent model. , 2000, Journal of molecular biology.

[8]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 7. Entropy effects on the performance of end-point binding free energy calculation approaches. , 2018, Physical chemistry chemical physics : PCCP.

[9]  Paolo Carloni,et al.  Conformational flexibility of the catalytic Asp dyad in HIV‐1 protease: An ab initio study on the free enzyme , 2000, Proteins.

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

[11]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

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

[13]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[14]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[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]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[17]  Martin J. Field,et al.  Free energy perturbation method for chemical reactions in the condensed phase: a dynamic approach based on a combined quantum and molecular mechanics potential , 1987 .

[18]  Robert T. C. Brownlee,et al.  Effect of atomic charge, solvation, entropy, and ligand protonation state on MM‐PB(GB)SA binding energies of HIV protease , 2012, J. Comput. Chem..

[19]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[20]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

[21]  Harry A. Stern,et al.  Development of a polarizable force field for proteins via ab initio quantum chemistry: First generation model and gas phase tests , 2002, J. Comput. Chem..

[22]  Mai Suan Li,et al.  Protocol for fast screening of multi-target drug candidates: Application to Alzheimer's disease. , 2017, Journal of molecular graphics & modelling.

[23]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[24]  J. H. Zhang,et al.  Energetics of protein backbone hydrogen bonds and their local electrostatic environment , 2014, Science China Chemistry.

[25]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

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

[27]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[28]  Xiao Zhu,et al.  Fast and accurate computation schemes for evaluating vibrational entropy of proteins , 2011, J. Comput. Chem..

[29]  Ye Mei,et al.  Folding of a helix is critically stabilized by polarization of backbone hydrogen bonds: study in explicit water. , 2012, The journal of physical chemistry. B.

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

[31]  U. Singh,et al.  A free energy perturbation study of solvation in methanol and dimethyl sulfoxide , 1990 .

[32]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

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

[34]  Jianing Song,et al.  An efficient method for computing excess free energy of liquid , 2017, Science China Chemistry.

[35]  Mengxin Li,et al.  Trypsin-Ligand binding affinities calculated using an effective interaction entropy method under polarized force field , 2017, Scientific Reports.

[36]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[37]  Emppu Salonen,et al.  Polarizable force fields. , 2013, Methods in molecular biology.

[38]  Cristiano Ruch Werneck Guimarães,et al.  Addressing Limitations with the MM-GB/SA Scoring Procedure using the WaterMap Method and Free Energy Perturbation Calculations , 2010, J. Chem. Inf. Model..

[39]  Ye Mei,et al.  A numerically stable restrained electrostatic potential charge fitting method , 2013, J. Comput. Chem..

[40]  C. Ji,et al.  Protein polarization is critical to stabilizing AF-2 and helix-2' domains in ligand binding to PPAR-gamma. , 2008, Journal of the American Chemical Society.

[41]  Parimal Kar,et al.  Origin of decrease in potency of darunavir and two related antiviral inhibitors against HIV-2 compared to HIV-1 protease. , 2012, The journal of physical chemistry. B.

[42]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

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

[44]  John Z. H. Zhang,et al.  Interaction Entropy for Computational Alanine Scanning , 2017, J. Chem. Inf. Model..

[45]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[46]  Zhihai Liu,et al.  Comparative Assessment of Scoring Functions on an Updated Benchmark: 2. Evaluation Methods and General Results , 2014, J. Chem. Inf. Model..

[47]  John Z H Zhang,et al.  Quantifying the stabilizing energy of the intraprotein hydrogen bond due to local mutation. , 2011, The journal of physical chemistry. B.

[48]  William A. McLaughlin,et al.  Evaluating the potency of HIV‐1 protease drugs to combat resistance , 2008, Proteins.

[49]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.

[50]  P A Kollman,et al.  An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis. , 2001, Journal of the American Chemical Society.

[51]  Jianing Song,et al.  Interaction entropy for computational alanine scanning in protein–protein binding , 2018 .

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

[53]  Michael E. Johnson,et al.  The Flavonoid Apigenin Is a Progesterone Receptor Modulator with In Vivo Activity in the Uterus , 2018, Hormones and Cancer.

[54]  Peter V. Coveney,et al.  Accurate Ensemble Molecular Dynamics Binding Free Energy Ranking of Multidrug-Resistant HIV-1 Proteases , 2010, J. Chem. Inf. Model..

[55]  Guanghong Wei,et al.  Orcein-Related Small Molecule O4 Destabilizes hIAPP Protofibrils by Interacting Mostly with the Amyloidogenic Core Region. , 2017, The journal of physical chemistry. B.

[56]  Marta A. S. Perez,et al.  Computational Alanine Scanning Mutagenesis: MM-PBSA vs TI. , 2013, Journal of chemical theory and computation.

[57]  Qinggang Zhang,et al.  Effect of electrostatic polarization and bridging water on CDK2-ligand binding affinities calculated using a highly efficient interaction entropy method. , 2017, Physical chemistry chemical physics : PCCP.

[58]  Chao-Yie Yang,et al.  Binding free energy contributions of interfacial waters in HIV-1 protease/inhibitor complexes. , 2006, Journal of the American Chemical Society.

[59]  Dan Li,et al.  Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking. , 2016, Physical chemistry chemical physics : PCCP.

[60]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

[61]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[62]  Peter A. Kollman,et al.  Computational alanine scanning of the 1:1 human growth hormone–receptor complex , 2002, J. Comput. Chem..

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

[65]  Ye Mei,et al.  Folding of a helix at room temperature is critically aided by electrostatic polarization of intraprotein hydrogen bonds. , 2010, Journal of the American Chemical Society.

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

[67]  Wei Zhang,et al.  Predicting drug resistance of the HIV‐1 protease using molecular interaction energy components , 2009, Proteins.

[68]  Lili Duan,et al.  Large-scale molecular dynamics simulation: Effect of polarization on thrombin-ligand binding energy , 2016, Scientific Reports.

[69]  J. H. Zhang,et al.  Protein's native structure is dynamically stabilized by electronic polarization , 2014 .

[70]  Xiaojie Xu,et al.  Recent Advances in Free Energy Calculations with a Combination of Molecular Mechanics and Continuum Models , 2006 .

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

[72]  Àngels González-Lafont,et al.  Inhibition of Mammalian 15-Lipoxygenase by Three Ebselen-like Drugs. A QM/MM and MM/PBSA Comparative Study. , 2017, The journal of physical chemistry. A.

[73]  Jing Zhang,et al.  Detecting and understanding combinatorial mutation patterns responsible for HIV drug resistance , 2010, Proceedings of the National Academy of Sciences.

[74]  T. Bhat,et al.  Structure of HIV-1 protease with KNI-272, a tight-binding transition-state analog containing allophenylnorstatine. , 1995, Structure.

[75]  Ye Mei,et al.  Developing polarized protein-specific charges for protein dynamics: MD free energy calculation of pKa shifts for Asp26/Asp20 in thioredoxin. , 2008, Biophysical journal.

[76]  Ruhong Zhou,et al.  Parametrizing a polarizable force field from ab initio data. I. The fluctuating point charge model , 1999 .

[77]  Ye Mei,et al.  Some practical approaches to treating electrostatic polarization of proteins. , 2014, Accounts of chemical research.

[78]  David A. Case,et al.  On finding stationary states on large-molecule potential energy surfaces , 1985 .

[79]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 3. The impact of force fields and ligand charge models. , 2013, The journal of physical chemistry. B.

[80]  Pedro Alexandrino Fernandes,et al.  New Parameters for Higher Accuracy in the Computation of Binding Free Energy Differences upon Alanine Scanning Mutagenesis on Protein-Protein Interfaces , 2017, J. Chem. Inf. Model..

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

[82]  Shantenu Jha,et al.  Computing Clinically Relevant Binding Free Energies of HIV-1 Protease Inhibitors , 2014, Journal of chemical theory and computation.

[83]  U. Singh,et al.  Free energy perturbation calculations on binding and catalysis after mutating Asn 155 in subtilisin , 1987, Nature.

[84]  J. H. Zhang,et al.  Simulation of NMR data reveals that proteins' local structures are stabilized by electronic polarization. , 2009, Journal of the American Chemical Society.

[85]  T. Bhat,et al.  Bound Water Molecules at the Interface between the HIV-1 Protease and a Potent Inhibitor, KNI-272, Determined by NMR , 1996 .

[86]  T. Lin,et al.  Contribution of the surface free energy perturbation to protein-solvent interactions. , 1994, Biochemistry.

[87]  C. R. W. Guimarães.  A Direct Comparison of the MM-GB/SA Scoring Procedure and Free-Energy Perturbation Calculations Using Carbonic Anhydrase as a Test Case: Strengths and Pitfalls of Each Approach. , 2011, Journal of chemical theory and computation.

[88]  J. H. Zhang,et al.  Computational Alanine Scanning with Interaction Entropy for Protein-Ligand Binding Free Energies. , 2018, Journal of chemical theory and computation.

[89]  J. H. Zhang,et al.  Computational search for aflatoxin binding proteins , 2017 .