Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA

J. Wang, P. Morin, W. Wang, P. A. Kollman

Research output: Contribution to journalArticle

569 Citations (Scopus)

Abstract

In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA™ another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 Å resolution crystal (structure the root-mean-square deviation of the 54 Cα around the binding site and the inhibitor is 1.1 Å). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.

Original languageEnglish (US)
Pages (from-to)5221-5230
Number of pages10
JournalJournal of the American Chemical Society
Volume123
Issue number22
DOIs
StatePublished - 2001

Fingerprint

efavirenz
Molecular mechanics
Mechanics
Free energy
HIV-1
Derivatives
Molecular Dynamics Simulation
Molecular dynamics
Crystal structure
Ligands
Proteins
Computer simulation
Binding sites
Mean square error
Binding Sites
Databases

ASJC Scopus subject areas

  • Chemistry(all)

Cite this

@article{8eeb167453fe49bc8e49fa96a1354303,
title = "Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA",
abstract = "In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA™ another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 {\AA} resolution crystal (structure the root-mean-square deviation of the 54 Cα around the binding site and the inhibitor is 1.1 {\AA}). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.",
author = "J. Wang and P. Morin and W. Wang and Kollman, {P. A.}",
year = "2001",
doi = "10.1021/ja003834q",
language = "English (US)",
volume = "123",
pages = "5221--5230",
journal = "Journal of the American Chemical Society",
issn = "0002-7863",
publisher = "American Chemical Society",
number = "22",

}

TY - JOUR

T1 - Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA

AU - Wang, J.

AU - Morin, P.

AU - Wang, W.

AU - Kollman, P. A.

PY - 2001

Y1 - 2001

N2 - In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA™ another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 Å resolution crystal (structure the root-mean-square deviation of the 54 Cα around the binding site and the inhibitor is 1.1 Å). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.

AB - In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA™ another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 Å resolution crystal (structure the root-mean-square deviation of the 54 Cα around the binding site and the inhibitor is 1.1 Å). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.

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

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

U2 - 10.1021/ja003834q

DO - 10.1021/ja003834q

M3 - Article

C2 - 11457384

AN - SCOPUS:0034811498

VL - 123

SP - 5221

EP - 5230

JO - Journal of the American Chemical Society

JF - Journal of the American Chemical Society

SN - 0002-7863

IS - 22

ER -