Calculation of protein-ligand binding affinities based on a fragment quantum mechanical method

Jinfeng Liu, Xianwei Wang, John Zhang, Xiao He

Research output: Contribution to journalArticle

Abstract

An electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method (J. Phys. Chem. A, 2013, 117, 7149) has been successfully used for efficient linear-scaling quantum mechanical (QM) calculations of protein energies. Furthermore, an efficient approach that combined the EE-GMFCC method with a conductor-like polarizable continuum model (CPCM), termed as EE-GMFCC-CPCM (J. Chem. Phys., 2013, 139, 214104), was developed for ab initio calculations of the electrostatic solvation energy of proteins. In this study, we applied the EE-GMFCC-CPCM method for the calculation of binding affinities of 14 avidin-biotin analogues. The calculation delineated the contributions of interaction energy and electrostatic solvation energy to the binding affinity. The binding affinity of each ligand bound to avidin was calculated at the HF/6-31G∗ and B3LYP/6-31G∗ levels with empirical dispersion corrections, respectively. The correlation coefficient (R) between the calculated binding energies and experimental values is 0.75 at the HF/6-31G∗-D level based on single complex structure calculations, as compared to 0.73 of the force field result. On the other hand, the correlation coefficient between the calculated binding energies and the experimental values is 0.85 at the B3LYP/6-31G∗-D level based on single complex structure calculations, and this correlation can be further improved to 0.88 when multiple snapshots are considered. Our study demonstrates that the EE-GMFCC-CPCM method is capable of providing reliable predictions of binding affinities for various ligands bound to the same target.

Original languageEnglish (US)
Pages (from-to)107020-107030
Number of pages11
JournalRSC Advances
Volume5
Issue number129
DOIs
StatePublished - Dec 3 2015

Fingerprint

Fractionation
Ligands
Proteins
Avidin
Solvation
Binding energy
Electrostatics
Biotin

ASJC Scopus subject areas

  • Chemical Engineering(all)
  • Chemistry(all)

Cite this

Calculation of protein-ligand binding affinities based on a fragment quantum mechanical method. / Liu, Jinfeng; Wang, Xianwei; Zhang, John; He, Xiao.

In: RSC Advances, Vol. 5, No. 129, 03.12.2015, p. 107020-107030.

Research output: Contribution to journalArticle

Liu, Jinfeng ; Wang, Xianwei ; Zhang, John ; He, Xiao. / Calculation of protein-ligand binding affinities based on a fragment quantum mechanical method. In: RSC Advances. 2015 ; Vol. 5, No. 129. pp. 107020-107030.
@article{37c5b747cdbe4a78a1165afbe93babb4,
title = "Calculation of protein-ligand binding affinities based on a fragment quantum mechanical method",
abstract = "An electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method (J. Phys. Chem. A, 2013, 117, 7149) has been successfully used for efficient linear-scaling quantum mechanical (QM) calculations of protein energies. Furthermore, an efficient approach that combined the EE-GMFCC method with a conductor-like polarizable continuum model (CPCM), termed as EE-GMFCC-CPCM (J. Chem. Phys., 2013, 139, 214104), was developed for ab initio calculations of the electrostatic solvation energy of proteins. In this study, we applied the EE-GMFCC-CPCM method for the calculation of binding affinities of 14 avidin-biotin analogues. The calculation delineated the contributions of interaction energy and electrostatic solvation energy to the binding affinity. The binding affinity of each ligand bound to avidin was calculated at the HF/6-31G∗ and B3LYP/6-31G∗ levels with empirical dispersion corrections, respectively. The correlation coefficient (R) between the calculated binding energies and experimental values is 0.75 at the HF/6-31G∗-D level based on single complex structure calculations, as compared to 0.73 of the force field result. On the other hand, the correlation coefficient between the calculated binding energies and the experimental values is 0.85 at the B3LYP/6-31G∗-D level based on single complex structure calculations, and this correlation can be further improved to 0.88 when multiple snapshots are considered. Our study demonstrates that the EE-GMFCC-CPCM method is capable of providing reliable predictions of binding affinities for various ligands bound to the same target.",
author = "Jinfeng Liu and Xianwei Wang and John Zhang and Xiao He",
year = "2015",
month = "12",
day = "3",
doi = "10.1039/c5ra20185c",
language = "English (US)",
volume = "5",
pages = "107020--107030",
journal = "RSC Advances",
issn = "2046-2069",
publisher = "Royal Society of Chemistry",
number = "129",

}

TY - JOUR

T1 - Calculation of protein-ligand binding affinities based on a fragment quantum mechanical method

AU - Liu, Jinfeng

AU - Wang, Xianwei

AU - Zhang, John

AU - He, Xiao

PY - 2015/12/3

Y1 - 2015/12/3

N2 - An electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method (J. Phys. Chem. A, 2013, 117, 7149) has been successfully used for efficient linear-scaling quantum mechanical (QM) calculations of protein energies. Furthermore, an efficient approach that combined the EE-GMFCC method with a conductor-like polarizable continuum model (CPCM), termed as EE-GMFCC-CPCM (J. Chem. Phys., 2013, 139, 214104), was developed for ab initio calculations of the electrostatic solvation energy of proteins. In this study, we applied the EE-GMFCC-CPCM method for the calculation of binding affinities of 14 avidin-biotin analogues. The calculation delineated the contributions of interaction energy and electrostatic solvation energy to the binding affinity. The binding affinity of each ligand bound to avidin was calculated at the HF/6-31G∗ and B3LYP/6-31G∗ levels with empirical dispersion corrections, respectively. The correlation coefficient (R) between the calculated binding energies and experimental values is 0.75 at the HF/6-31G∗-D level based on single complex structure calculations, as compared to 0.73 of the force field result. On the other hand, the correlation coefficient between the calculated binding energies and the experimental values is 0.85 at the B3LYP/6-31G∗-D level based on single complex structure calculations, and this correlation can be further improved to 0.88 when multiple snapshots are considered. Our study demonstrates that the EE-GMFCC-CPCM method is capable of providing reliable predictions of binding affinities for various ligands bound to the same target.

AB - An electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method (J. Phys. Chem. A, 2013, 117, 7149) has been successfully used for efficient linear-scaling quantum mechanical (QM) calculations of protein energies. Furthermore, an efficient approach that combined the EE-GMFCC method with a conductor-like polarizable continuum model (CPCM), termed as EE-GMFCC-CPCM (J. Chem. Phys., 2013, 139, 214104), was developed for ab initio calculations of the electrostatic solvation energy of proteins. In this study, we applied the EE-GMFCC-CPCM method for the calculation of binding affinities of 14 avidin-biotin analogues. The calculation delineated the contributions of interaction energy and electrostatic solvation energy to the binding affinity. The binding affinity of each ligand bound to avidin was calculated at the HF/6-31G∗ and B3LYP/6-31G∗ levels with empirical dispersion corrections, respectively. The correlation coefficient (R) between the calculated binding energies and experimental values is 0.75 at the HF/6-31G∗-D level based on single complex structure calculations, as compared to 0.73 of the force field result. On the other hand, the correlation coefficient between the calculated binding energies and the experimental values is 0.85 at the B3LYP/6-31G∗-D level based on single complex structure calculations, and this correlation can be further improved to 0.88 when multiple snapshots are considered. Our study demonstrates that the EE-GMFCC-CPCM method is capable of providing reliable predictions of binding affinities for various ligands bound to the same target.

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

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

U2 - 10.1039/c5ra20185c

DO - 10.1039/c5ra20185c

M3 - Article

VL - 5

SP - 107020

EP - 107030

JO - RSC Advances

JF - RSC Advances

SN - 2046-2069

IS - 129

ER -