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

Yuchen Li, Yalong Cong, Guoqiang Feng, Susu Zhong, John Zhang, Huiyong Sun, Lili Duan

Research output: Contribution to journalArticle

Abstract

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.

Original languageEnglish (US)
Article number064101
JournalStructural Dynamics
Volume5
Issue number6
DOIs
StatePublished - Nov 1 2018

Fingerprint

human immunodeficiency virus
Free energy
Permittivity
free energy
permittivity
predictions
protease
Entropy
entropy
correlation coefficients
inhibitors
interactions
Molecular mechanics
Protease Inhibitors
Viruses
Polarization
Molecules
Water
polarization
electronics

ASJC Scopus subject areas

  • Radiation
  • Instrumentation
  • Condensed Matter Physics
  • Spectroscopy

Cite this

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction. / Li, Yuchen; Cong, Yalong; Feng, Guoqiang; Zhong, Susu; Zhang, John; Sun, Huiyong; Duan, Lili.

In: Structural Dynamics, Vol. 5, No. 6, 064101, 01.11.2018.

Research output: Contribution to journalArticle

Li, Yuchen ; Cong, Yalong ; Feng, Guoqiang ; Zhong, Susu ; Zhang, John ; Sun, Huiyong ; Duan, Lili. / The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction. In: Structural Dynamics. 2018 ; Vol. 5, No. 6.
@article{aa750d4fe9fa461997e49cf23305995e,
title = "The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction",
abstract = "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.",
author = "Yuchen Li and Yalong Cong and Guoqiang Feng and Susu Zhong and John Zhang and Huiyong Sun and Lili Duan",
year = "2018",
month = "11",
day = "1",
doi = "10.1063/1.5058172",
language = "English (US)",
volume = "5",
journal = "Structural Dynamics",
issn = "2329-7778",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "6",

}

TY - JOUR

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

AU - Li, Yuchen

AU - Cong, Yalong

AU - Feng, Guoqiang

AU - Zhong, Susu

AU - Zhang, John

AU - Sun, Huiyong

AU - Duan, Lili

PY - 2018/11/1

Y1 - 2018/11/1

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=85058447538&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85058447538&partnerID=8YFLogxK

U2 - 10.1063/1.5058172

DO - 10.1063/1.5058172

M3 - Article

VL - 5

JO - Structural Dynamics

JF - Structural Dynamics

SN - 2329-7778

IS - 6

M1 - 064101

ER -