+ Site Statistics
+ Search Articles
+ PDF Full Text Service
How our service works
Request PDF Full Text
+ Follow Us
Follow on Facebook
Follow on Twitter
Follow on LinkedIn
+ Subscribe to Site Feeds
Most Shared
PDF Full Text
+ Translate
+ Recently Requested

Improving Prediction Accuracy of Binding Free Energies and Poses of HIV Integrase Complexes Using the Binding Energy Distribution Analysis Method with Flattening Potentials



Improving Prediction Accuracy of Binding Free Energies and Poses of HIV Integrase Complexes Using the Binding Energy Distribution Analysis Method with Flattening Potentials



Journal of Chemical Information and Modeling 58(7): 1356-1371



To accelerate conformation sampling of slow dynamics from receptor or ligand, we introduced flattening potentials on selected bonded and nonbonded intramolecular interactions to the binding energy distribution analysis method (BEDAM) for calculating absolute binding free energies of protein-ligand complexes using an implicit solvent model and implemented flattening BEDAM using the asynchronous replica exchange (AsyncRE) framework for performing large scale replica exchange molecular dynamics (REMD) simulations. The advantage of using the flattening feature to reduce high energy barriers was exhibited first by the p-xylene-T4 lysozyme complex, where the intramolecular interactions of a protein side chain on the binding site were flattened to accelerate the conformational transition of the side chain from the trans to the gauche state when the p-xylene ligand is present in the binding site. Much more extensive flattening BEDAM simulations were performed for 53 experimental binders and 248 nonbinders of HIV-1 integrase which formed the SAMPL4 challenge, with the total simulation time of 24.3 μs. We demonstrated that the flattening BEDAM simulations not only substantially increase the number of true positives (and reduce false negatives) but also improve the prediction accuracy of binding poses of experimental binders. Furthermore, the values of area under the curve (AUC) of receiver operating characteristic (ROC) and the enrichment factors at 20% cutoff calculated from the flattening BEDAM simulations were improved significantly in comparison with that of simulations without flattening as we previously reported for the whole SAMPL4 database. Detailed analysis found that the improved ability to discriminate the binding free energies between the binders and nonbinders is due to the fact that the flattening simulations reduce the reorganization free energy penalties of binders and decrease the overlap of binding free energy distributions of binders relative to that of nonbinders. This happens because the conformational ensemble distributions for both the ligand and protein in solution match those at the fully coupled (complex) state more closely when the systems are more fully sampled after the flattening potentials are applied to the intermediate states.

Please choose payment method:






(PDF emailed within 0-6 h: $19.90)

Accession: 065370730

Download citation: RISBibTeXText

PMID: 29927237

DOI: 10.1021/acs.jcim.8b00194


Related references

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation. Journal of Chemical Theory and Computation 12(5): 2459-2470, 2016

Assessing the performance of MM/PBSA and MM/GBSA methods. 8. Predicting binding free energies and poses of protein-RNA complexes. Rna 24(9): 1183-1194, 2018

D3R Grand Challenge 2: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies. Journal of Computer-Aided Molecular Design 32(1): 1-20, 2018

Improving the LIE Method for Binding Free Energy Calculations of Protein-Ligand Complexes. Journal of Chemical Information and Modeling 55(9): 1867-1877, 2015

Improving the Accuracy of the Linear Interaction Energy Method for Solvation Free Energies. Journal of Chemical Theory and Computation 3(6): 2162-2175, 2007

Assessing the performance of MM/PBSA and MM/GBSA methods. 9. Prediction reliability of binding affinities and binding poses for protein-peptide complexes. Physical Chemistry Chemical Physics 21(19): 10135-10145, 2019

Prediction of the binding free energies of new TIBO-like HIV-1 reverse transcriptase inhibitors using a combination of PROFEC, PB/SA, CMC/MD, and free energy calculations. Journal of Medicinal Chemistry 42(5): 868-881, 1999

How to deal with multiple binding poses in alchemical relative protein-ligand binding free energy calculations. Journal of Chemical Theory and Computation 11(6): 2670-2679, 2015

Absolute binding free energies for octa-acids and guests in SAMPL5 : Evaluating binding free energies for octa-acid and guest complexes in the SAMPL5 blind challenge. Journal of Computer-Aided Molecular Design 31(1): 107-118, 2017

Free energies underlying ion binding and transport in protein channels free energy perturbation simulations of ion binding and selectivity for valinomycin. Journal of the Chemical Society Faraday Transactions 87(13): 2099-2109, 1991

Computational prediction of monosaccharide binding free energies to lectins with linear interaction energy models. Journal of Computational Chemistry 33(29): 2340-2350, 2012

Structure and thermodynamics of RNA-protein binding: using molecular dynamics and free energy analyses to calculate the free energies of binding and conformational change. Journal of Molecular Biology 297(5): 1145-1158, 2000

Advances in binding free energies calculations: QM/MM-based free energy perturbation method for drug design. Current Pharmaceutical Design 19(26): 4674-4686, 2013

Advances in Binding Free Energies Calculations: Qm/Mm-Based Free Energy Perturbation Method for Drug Design. Current Pharmaceutical Design 19(26): 4674-4686, 2013

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. Physical Chemistry Chemical Physics 18(32): 22129-22139, 2016