Large-scale conformational sampling of proteins using temperature- accelerated molecular dynamics

Cameron F. Abrams, Eric Vanden Eijnden

Research output: Contribution to journalArticle

Abstract

We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t → r″ transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 Å, with a twist of 90° relative to the equatorial domain, and the root-mean-squared deviation relative to the r″ conformer is reduced from 13 to 5 Å, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 Å from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.

Original languageEnglish (US)
Pages (from-to)4961-4966
Number of pages6
JournalProceedings of the National Academy of Sciences of the United States of America
Volume107
Issue number11
DOIs
StatePublished - Mar 16 2010

Fingerprint

Molecular Dynamics Simulation
Temperature
Proteins
HIV-1
Antibodies

Keywords

  • Biophysical simulation
  • Collective variables
  • Domain motion
  • Free energy

ASJC Scopus subject areas

  • General

Cite this

@article{09bde1ee383d4afea9307d4f3fb6cd2e,
title = "Large-scale conformational sampling of proteins using temperature- accelerated molecular dynamics",
abstract = "We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t → r″ transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 {\AA}, with a twist of 90° relative to the equatorial domain, and the root-mean-squared deviation relative to the r″ conformer is reduced from 13 to 5 {\AA}, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 {\AA} from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.",
keywords = "Biophysical simulation, Collective variables, Domain motion, Free energy",
author = "Abrams, {Cameron F.} and {Vanden Eijnden}, Eric",
year = "2010",
month = "3",
day = "16",
doi = "10.1073/pnas.0914540107",
language = "English (US)",
volume = "107",
pages = "4961--4966",
journal = "Proceedings of the National Academy of Sciences of the United States of America",
issn = "0027-8424",
number = "11",

}

TY - JOUR

T1 - Large-scale conformational sampling of proteins using temperature- accelerated molecular dynamics

AU - Abrams, Cameron F.

AU - Vanden Eijnden, Eric

PY - 2010/3/16

Y1 - 2010/3/16

N2 - We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t → r″ transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 Å, with a twist of 90° relative to the equatorial domain, and the root-mean-squared deviation relative to the r″ conformer is reduced from 13 to 5 Å, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 Å from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.

AB - We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t → r″ transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 Å, with a twist of 90° relative to the equatorial domain, and the root-mean-squared deviation relative to the r″ conformer is reduced from 13 to 5 Å, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 Å from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.

KW - Biophysical simulation

KW - Collective variables

KW - Domain motion

KW - Free energy

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

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

U2 - 10.1073/pnas.0914540107

DO - 10.1073/pnas.0914540107

M3 - Article

C2 - 20194785

AN - SCOPUS:77950448255

VL - 107

SP - 4961

EP - 4966

JO - Proceedings of the National Academy of Sciences of the United States of America

JF - Proceedings of the National Academy of Sciences of the United States of America

SN - 0027-8424

IS - 11

ER -