Using thermodynamic integration MD simulation to compute relative protein-ligand binding free energy of a GSK3β kinase inhibitor and its analogs

Hsing Chou Lee, Wen Chi Hsu, An Lun Liu, Chia Jen Hsu, Ying Chieh Sun

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

Thermodynamic integration molecular dynamics simulation was used to investigate how TI-MD simulation preforms in reproducing relative protein-ligand binding free energy of a pair of analogous GSK3β kinase inhibitors of available experimental data (see Fig. 1), and to predict the affinity for other analogs. The computation for the pair gave a ΔΔG of 1.0 kcal/mol, which was in reasonably good agreement with the experimental value of -0.1 kcal/mol. The error bar was estimated at 0.5 kcal/mol. Subsequently, we employed the same protocol to proceed with simulations to find analogous inhibitors with a stronger affinity. Four analogs with a substitution at one site inside the binding pocket were the first to be tried, but no significant enhancement in affinity was found. Subsequent simulations for another 7 analogs was focused on substitutions at the benzene ring of another site, which gave two analogs (analogs 9 and 10) with ΔΔG values of -0.6 and -0.8 kcal/mol, respectively. Both analogs had a OH group at the meta position and another OH group at the ortho position at the other side of the benzene ring, as shown in Table 3. To explore further, another 4 analogs with this characteristic were investigated. Three analogs with ΔΔG values of -2.2, -1.7 and -1.2 kcal/mol, respectively, were found. Hydrogen bond analysis suggested that the additional hydrogen bonds of the added OH groups with Gln185 and/or Asn64, which did not appear in the reference inhibitor or as an analog with one substitution only in the examined cases, were the main contributors to an enhanced affinity. A prediction for better inhibitors should interest experimentalists of enzyme and/or cell assays. Analysis of the interactions between GSK3β kinase and the investigated analogs will be useful in the design of GSK3β kinase inhibitors for compounds of this class.

Original languageEnglish
Pages (from-to)37-49
Number of pages13
JournalJournal of Molecular Graphics and Modelling
Volume51
DOIs
Publication statusPublished - 2014 Jun

Fingerprint

inhibitors
Free energy
Substitution reactions
Phosphotransferases
free energy
Ligands
Thermodynamics
analogs
Benzene
proteins
Proteins
thermodynamics
ligands
Hydrogen bonds
simulation
affinity
Molecular dynamics
Assays
Enzymes
Binding Sites

Keywords

  • GSK3β kinase
  • Inhibitor
  • MD simulation
  • Relative binding free energy
  • Thermodynamic integration

ASJC Scopus subject areas

  • Spectroscopy
  • Physical and Theoretical Chemistry
  • Computer Graphics and Computer-Aided Design
  • Materials Chemistry

Cite this

Using thermodynamic integration MD simulation to compute relative protein-ligand binding free energy of a GSK3β kinase inhibitor and its analogs. / Lee, Hsing Chou; Hsu, Wen Chi; Liu, An Lun; Hsu, Chia Jen; Sun, Ying Chieh.

In: Journal of Molecular Graphics and Modelling, Vol. 51, 06.2014, p. 37-49.

Research output: Contribution to journalArticle

@article{f5032ca920a741919eb6a6f64c72febd,
title = "Using thermodynamic integration MD simulation to compute relative protein-ligand binding free energy of a GSK3β kinase inhibitor and its analogs",
abstract = "Thermodynamic integration molecular dynamics simulation was used to investigate how TI-MD simulation preforms in reproducing relative protein-ligand binding free energy of a pair of analogous GSK3β kinase inhibitors of available experimental data (see Fig. 1), and to predict the affinity for other analogs. The computation for the pair gave a ΔΔG of 1.0 kcal/mol, which was in reasonably good agreement with the experimental value of -0.1 kcal/mol. The error bar was estimated at 0.5 kcal/mol. Subsequently, we employed the same protocol to proceed with simulations to find analogous inhibitors with a stronger affinity. Four analogs with a substitution at one site inside the binding pocket were the first to be tried, but no significant enhancement in affinity was found. Subsequent simulations for another 7 analogs was focused on substitutions at the benzene ring of another site, which gave two analogs (analogs 9 and 10) with ΔΔG values of -0.6 and -0.8 kcal/mol, respectively. Both analogs had a OH group at the meta position and another OH group at the ortho position at the other side of the benzene ring, as shown in Table 3. To explore further, another 4 analogs with this characteristic were investigated. Three analogs with ΔΔG values of -2.2, -1.7 and -1.2 kcal/mol, respectively, were found. Hydrogen bond analysis suggested that the additional hydrogen bonds of the added OH groups with Gln185 and/or Asn64, which did not appear in the reference inhibitor or as an analog with one substitution only in the examined cases, were the main contributors to an enhanced affinity. A prediction for better inhibitors should interest experimentalists of enzyme and/or cell assays. Analysis of the interactions between GSK3β kinase and the investigated analogs will be useful in the design of GSK3β kinase inhibitors for compounds of this class.",
keywords = "GSK3β kinase, Inhibitor, MD simulation, Relative binding free energy, Thermodynamic integration",
author = "Lee, {Hsing Chou} and Hsu, {Wen Chi} and Liu, {An Lun} and Hsu, {Chia Jen} and Sun, {Ying Chieh}",
year = "2014",
month = "6",
doi = "10.1016/j.jmgm.2014.04.010",
language = "English",
volume = "51",
pages = "37--49",
journal = "Journal of Molecular Graphics and Modelling",
issn = "1093-3263",
publisher = "Elsevier Inc.",

}

TY - JOUR

T1 - Using thermodynamic integration MD simulation to compute relative protein-ligand binding free energy of a GSK3β kinase inhibitor and its analogs

AU - Lee, Hsing Chou

AU - Hsu, Wen Chi

AU - Liu, An Lun

AU - Hsu, Chia Jen

AU - Sun, Ying Chieh

PY - 2014/6

Y1 - 2014/6

N2 - Thermodynamic integration molecular dynamics simulation was used to investigate how TI-MD simulation preforms in reproducing relative protein-ligand binding free energy of a pair of analogous GSK3β kinase inhibitors of available experimental data (see Fig. 1), and to predict the affinity for other analogs. The computation for the pair gave a ΔΔG of 1.0 kcal/mol, which was in reasonably good agreement with the experimental value of -0.1 kcal/mol. The error bar was estimated at 0.5 kcal/mol. Subsequently, we employed the same protocol to proceed with simulations to find analogous inhibitors with a stronger affinity. Four analogs with a substitution at one site inside the binding pocket were the first to be tried, but no significant enhancement in affinity was found. Subsequent simulations for another 7 analogs was focused on substitutions at the benzene ring of another site, which gave two analogs (analogs 9 and 10) with ΔΔG values of -0.6 and -0.8 kcal/mol, respectively. Both analogs had a OH group at the meta position and another OH group at the ortho position at the other side of the benzene ring, as shown in Table 3. To explore further, another 4 analogs with this characteristic were investigated. Three analogs with ΔΔG values of -2.2, -1.7 and -1.2 kcal/mol, respectively, were found. Hydrogen bond analysis suggested that the additional hydrogen bonds of the added OH groups with Gln185 and/or Asn64, which did not appear in the reference inhibitor or as an analog with one substitution only in the examined cases, were the main contributors to an enhanced affinity. A prediction for better inhibitors should interest experimentalists of enzyme and/or cell assays. Analysis of the interactions between GSK3β kinase and the investigated analogs will be useful in the design of GSK3β kinase inhibitors for compounds of this class.

AB - Thermodynamic integration molecular dynamics simulation was used to investigate how TI-MD simulation preforms in reproducing relative protein-ligand binding free energy of a pair of analogous GSK3β kinase inhibitors of available experimental data (see Fig. 1), and to predict the affinity for other analogs. The computation for the pair gave a ΔΔG of 1.0 kcal/mol, which was in reasonably good agreement with the experimental value of -0.1 kcal/mol. The error bar was estimated at 0.5 kcal/mol. Subsequently, we employed the same protocol to proceed with simulations to find analogous inhibitors with a stronger affinity. Four analogs with a substitution at one site inside the binding pocket were the first to be tried, but no significant enhancement in affinity was found. Subsequent simulations for another 7 analogs was focused on substitutions at the benzene ring of another site, which gave two analogs (analogs 9 and 10) with ΔΔG values of -0.6 and -0.8 kcal/mol, respectively. Both analogs had a OH group at the meta position and another OH group at the ortho position at the other side of the benzene ring, as shown in Table 3. To explore further, another 4 analogs with this characteristic were investigated. Three analogs with ΔΔG values of -2.2, -1.7 and -1.2 kcal/mol, respectively, were found. Hydrogen bond analysis suggested that the additional hydrogen bonds of the added OH groups with Gln185 and/or Asn64, which did not appear in the reference inhibitor or as an analog with one substitution only in the examined cases, were the main contributors to an enhanced affinity. A prediction for better inhibitors should interest experimentalists of enzyme and/or cell assays. Analysis of the interactions between GSK3β kinase and the investigated analogs will be useful in the design of GSK3β kinase inhibitors for compounds of this class.

KW - GSK3β kinase

KW - Inhibitor

KW - MD simulation

KW - Relative binding free energy

KW - Thermodynamic integration

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

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

U2 - 10.1016/j.jmgm.2014.04.010

DO - 10.1016/j.jmgm.2014.04.010

M3 - Article

C2 - 24858254

AN - SCOPUS:84901271322

VL - 51

SP - 37

EP - 49

JO - Journal of Molecular Graphics and Modelling

JF - Journal of Molecular Graphics and Modelling

SN - 1093-3263

ER -