Computation of relative binding free energy for an inhibitor and its analogs binding with erk kinase using thermodynamic integration md simulation

Kuan Wei Wu, Po Chin Chen, Jun Wang, Ying-Chieh Sun

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

In the present study, we carried out thermodynamic integration molecular dynamics simulation for a pair of analogous inhibitors binding with Erk kinase to investigate how computation performs in reproducing the relative binding free energy. The computation with BCC-AM1 charges for ligands gave -1.1 kcal/mol, deviated from experimental value of -2.3 kcal/mol by 1.2 kcal/mol, in good agreement with experimental result. The error of computed value was estimated to be 0.5 kcal/mol. To obtain convergence, switching vdw interaction on and off required approximately 10 times more CPU time than switching charges. Residue-based contributions and hydrogen bonding were analyzed and discussed. Furthermore, subsequent simulation using RESP charge for ligand gave ΔΔG of -1.6 kcal/mol. The computed results are better than the result of -5.6 kcal/mol estimated using PBSA method in a previous study. Based on these results, we further carried out computations to predict ΔΔG for five new analogs, focusing on placing polar and nonpolar functional groups at the meta site of benzene ring shown in the Fig. 1, to see if these ligands have better binding affinity than the above ligands. The computations resulted that a ligand with polar -OH group has better binding affinity than the previous examined ligand by ~ 2.0 kcal/ mol and two other ligands have better affinity by ~ 1.0 kcal/mol. The predicted better inhibitors of this kind should be of interest to experimentalist for future experimental enzyme and/or cell assays.

Original languageEnglish
Pages (from-to)1159-1169
Number of pages11
JournalJournal of Computer-Aided Molecular Design
Volume26
Issue number10
DOIs
Publication statusPublished - 2012 Oct 1

Fingerprint

Thermodynamics
inhibitors
Free energy
Phosphotransferases
free energy
Ligands
analogs
thermodynamics
ligands
simulation
affinity
Molecular Dynamics Simulation
Hydrogen Bonding
Benzene
Functional groups
Program processors
Molecular dynamics
enzymes
Assays
Hydrogen bonds

Keywords

  • Erk kinase
  • Inhibitor
  • MD simulation
  • Relative binding free energy
  • Thermodynamic integration

ASJC Scopus subject areas

  • Drug Discovery
  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Computation of relative binding free energy for an inhibitor and its analogs binding with erk kinase using thermodynamic integration md simulation. / Wu, Kuan Wei; Chen, Po Chin; Wang, Jun; Sun, Ying-Chieh.

In: Journal of Computer-Aided Molecular Design, Vol. 26, No. 10, 01.10.2012, p. 1159-1169.

Research output: Contribution to journalArticle

@article{eb54b3eb1d1748679876deabf24ae94d,
title = "Computation of relative binding free energy for an inhibitor and its analogs binding with erk kinase using thermodynamic integration md simulation",
abstract = "In the present study, we carried out thermodynamic integration molecular dynamics simulation for a pair of analogous inhibitors binding with Erk kinase to investigate how computation performs in reproducing the relative binding free energy. The computation with BCC-AM1 charges for ligands gave -1.1 kcal/mol, deviated from experimental value of -2.3 kcal/mol by 1.2 kcal/mol, in good agreement with experimental result. The error of computed value was estimated to be 0.5 kcal/mol. To obtain convergence, switching vdw interaction on and off required approximately 10 times more CPU time than switching charges. Residue-based contributions and hydrogen bonding were analyzed and discussed. Furthermore, subsequent simulation using RESP charge for ligand gave ΔΔG of -1.6 kcal/mol. The computed results are better than the result of -5.6 kcal/mol estimated using PBSA method in a previous study. Based on these results, we further carried out computations to predict ΔΔG for five new analogs, focusing on placing polar and nonpolar functional groups at the meta site of benzene ring shown in the Fig. 1, to see if these ligands have better binding affinity than the above ligands. The computations resulted that a ligand with polar -OH group has better binding affinity than the previous examined ligand by ~ 2.0 kcal/ mol and two other ligands have better affinity by ~ 1.0 kcal/mol. The predicted better inhibitors of this kind should be of interest to experimentalist for future experimental enzyme and/or cell assays.",
keywords = "Erk kinase, Inhibitor, MD simulation, Relative binding free energy, Thermodynamic integration",
author = "Wu, {Kuan Wei} and Chen, {Po Chin} and Jun Wang and Ying-Chieh Sun",
year = "2012",
month = "10",
day = "1",
doi = "10.1007/s10822-012-9606-6",
language = "English",
volume = "26",
pages = "1159--1169",
journal = "Journal of Computer-Aided Molecular Design",
issn = "0920-654X",
publisher = "Springer Netherlands",
number = "10",

}

TY - JOUR

T1 - Computation of relative binding free energy for an inhibitor and its analogs binding with erk kinase using thermodynamic integration md simulation

AU - Wu, Kuan Wei

AU - Chen, Po Chin

AU - Wang, Jun

AU - Sun, Ying-Chieh

PY - 2012/10/1

Y1 - 2012/10/1

N2 - In the present study, we carried out thermodynamic integration molecular dynamics simulation for a pair of analogous inhibitors binding with Erk kinase to investigate how computation performs in reproducing the relative binding free energy. The computation with BCC-AM1 charges for ligands gave -1.1 kcal/mol, deviated from experimental value of -2.3 kcal/mol by 1.2 kcal/mol, in good agreement with experimental result. The error of computed value was estimated to be 0.5 kcal/mol. To obtain convergence, switching vdw interaction on and off required approximately 10 times more CPU time than switching charges. Residue-based contributions and hydrogen bonding were analyzed and discussed. Furthermore, subsequent simulation using RESP charge for ligand gave ΔΔG of -1.6 kcal/mol. The computed results are better than the result of -5.6 kcal/mol estimated using PBSA method in a previous study. Based on these results, we further carried out computations to predict ΔΔG for five new analogs, focusing on placing polar and nonpolar functional groups at the meta site of benzene ring shown in the Fig. 1, to see if these ligands have better binding affinity than the above ligands. The computations resulted that a ligand with polar -OH group has better binding affinity than the previous examined ligand by ~ 2.0 kcal/ mol and two other ligands have better affinity by ~ 1.0 kcal/mol. The predicted better inhibitors of this kind should be of interest to experimentalist for future experimental enzyme and/or cell assays.

AB - In the present study, we carried out thermodynamic integration molecular dynamics simulation for a pair of analogous inhibitors binding with Erk kinase to investigate how computation performs in reproducing the relative binding free energy. The computation with BCC-AM1 charges for ligands gave -1.1 kcal/mol, deviated from experimental value of -2.3 kcal/mol by 1.2 kcal/mol, in good agreement with experimental result. The error of computed value was estimated to be 0.5 kcal/mol. To obtain convergence, switching vdw interaction on and off required approximately 10 times more CPU time than switching charges. Residue-based contributions and hydrogen bonding were analyzed and discussed. Furthermore, subsequent simulation using RESP charge for ligand gave ΔΔG of -1.6 kcal/mol. The computed results are better than the result of -5.6 kcal/mol estimated using PBSA method in a previous study. Based on these results, we further carried out computations to predict ΔΔG for five new analogs, focusing on placing polar and nonpolar functional groups at the meta site of benzene ring shown in the Fig. 1, to see if these ligands have better binding affinity than the above ligands. The computations resulted that a ligand with polar -OH group has better binding affinity than the previous examined ligand by ~ 2.0 kcal/ mol and two other ligands have better affinity by ~ 1.0 kcal/mol. The predicted better inhibitors of this kind should be of interest to experimentalist for future experimental enzyme and/or cell assays.

KW - Erk kinase

KW - Inhibitor

KW - MD simulation

KW - Relative binding free energy

KW - Thermodynamic integration

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

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

U2 - 10.1007/s10822-012-9606-6

DO - 10.1007/s10822-012-9606-6

M3 - Article

C2 - 22986633

AN - SCOPUS:84869499151

VL - 26

SP - 1159

EP - 1169

JO - Journal of Computer-Aided Molecular Design

JF - Journal of Computer-Aided Molecular Design

SN - 0920-654X

IS - 10

ER -