Transition-state structural refinement with GRACE and CHARMM: Flexible QM/MM modelling for lactate dehydrogenase

Alexander J. Turner, Vicente Moliner, Ian H. Williams

Research output: Contribution to journalArticle

122 Citations (Scopus)

Abstract

Realistic simulations of chemical reactions require the use not only of methods capable of describing accurately the energy of molecules undergoing bonding changes within a particular chemical environment, but also of methods capable of exploring topographical features of significance on energy hypersurfaces spanning perhaps several thousand degrees of freedom. Hybrid quantum-mechanical/molecular-mechanical techniques show much promise for the first task, but existing computer codes are inadequate for the second. Application of these methods to real chemical problems demands new tools for location and characterisation of saddle-points, intrinsic reaction coordinates, hessians and vibrational frequencies for very large flexible systems. Algorithms capable of performing these tasks have been incorporated in a new software package, GRACE, which provides a non-invasive interface between popular codes for quantum chemistry and molecular dynamics and modelling. Transition structures (TSs) have been refined by this novel procedure, using a combined AM1/CHARMM24/TIP3P potential, involving full gradient relaxation of the positions of 1900-2000 atoms of a solvated enzyme- substrate complex (lactate dehydrogenase/NADH/pyruvate/water). Six different starting structures (arbitrarily selected from a molecular dynamics trajectory for the enzyme-substrate complex) lead to six different TSs. Although the essential features of these TSs are invariant, the relative dispositions of active-site residues differ quite significantly. The transition state for the enzymic reaction would represent an average of the properties of many, nearly degenerate TSs. This insight emerges only as a consequence of the flexible model of the active site employed in this study.

Original languageEnglish
Pages (from-to)1323-1331
Number of pages9
JournalPhysical Chemistry Chemical Physics
Volume1
Issue number6
DOIs
Publication statusPublished - 15 Mar 1999

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Transition-state structural refinement with GRACE and CHARMM : Flexible QM/MM modelling for lactate dehydrogenase. / Turner, Alexander J.; Moliner, Vicente; Williams, Ian H.

In: Physical Chemistry Chemical Physics, Vol. 1, No. 6, 15.03.1999, p. 1323-1331.

Research output: Contribution to journalArticle

@article{d4bc2b9039c54e82ac489782bef61271,
title = "Transition-state structural refinement with GRACE and CHARMM: Flexible QM/MM modelling for lactate dehydrogenase",
abstract = "Realistic simulations of chemical reactions require the use not only of methods capable of describing accurately the energy of molecules undergoing bonding changes within a particular chemical environment, but also of methods capable of exploring topographical features of significance on energy hypersurfaces spanning perhaps several thousand degrees of freedom. Hybrid quantum-mechanical/molecular-mechanical techniques show much promise for the first task, but existing computer codes are inadequate for the second. Application of these methods to real chemical problems demands new tools for location and characterisation of saddle-points, intrinsic reaction coordinates, hessians and vibrational frequencies for very large flexible systems. Algorithms capable of performing these tasks have been incorporated in a new software package, GRACE, which provides a non-invasive interface between popular codes for quantum chemistry and molecular dynamics and modelling. Transition structures (TSs) have been refined by this novel procedure, using a combined AM1/CHARMM24/TIP3P potential, involving full gradient relaxation of the positions of 1900-2000 atoms of a solvated enzyme- substrate complex (lactate dehydrogenase/NADH/pyruvate/water). Six different starting structures (arbitrarily selected from a molecular dynamics trajectory for the enzyme-substrate complex) lead to six different TSs. Although the essential features of these TSs are invariant, the relative dispositions of active-site residues differ quite significantly. The transition state for the enzymic reaction would represent an average of the properties of many, nearly degenerate TSs. This insight emerges only as a consequence of the flexible model of the active site employed in this study.",
author = "Turner, {Alexander J.} and Vicente Moliner and Williams, {Ian H.}",
year = "1999",
month = "3",
day = "15",
doi = "10.1039/a808364i",
language = "English",
volume = "1",
pages = "1323--1331",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "6",

}

TY - JOUR

T1 - Transition-state structural refinement with GRACE and CHARMM

T2 - Flexible QM/MM modelling for lactate dehydrogenase

AU - Turner, Alexander J.

AU - Moliner, Vicente

AU - Williams, Ian H.

PY - 1999/3/15

Y1 - 1999/3/15

N2 - Realistic simulations of chemical reactions require the use not only of methods capable of describing accurately the energy of molecules undergoing bonding changes within a particular chemical environment, but also of methods capable of exploring topographical features of significance on energy hypersurfaces spanning perhaps several thousand degrees of freedom. Hybrid quantum-mechanical/molecular-mechanical techniques show much promise for the first task, but existing computer codes are inadequate for the second. Application of these methods to real chemical problems demands new tools for location and characterisation of saddle-points, intrinsic reaction coordinates, hessians and vibrational frequencies for very large flexible systems. Algorithms capable of performing these tasks have been incorporated in a new software package, GRACE, which provides a non-invasive interface between popular codes for quantum chemistry and molecular dynamics and modelling. Transition structures (TSs) have been refined by this novel procedure, using a combined AM1/CHARMM24/TIP3P potential, involving full gradient relaxation of the positions of 1900-2000 atoms of a solvated enzyme- substrate complex (lactate dehydrogenase/NADH/pyruvate/water). Six different starting structures (arbitrarily selected from a molecular dynamics trajectory for the enzyme-substrate complex) lead to six different TSs. Although the essential features of these TSs are invariant, the relative dispositions of active-site residues differ quite significantly. The transition state for the enzymic reaction would represent an average of the properties of many, nearly degenerate TSs. This insight emerges only as a consequence of the flexible model of the active site employed in this study.

AB - Realistic simulations of chemical reactions require the use not only of methods capable of describing accurately the energy of molecules undergoing bonding changes within a particular chemical environment, but also of methods capable of exploring topographical features of significance on energy hypersurfaces spanning perhaps several thousand degrees of freedom. Hybrid quantum-mechanical/molecular-mechanical techniques show much promise for the first task, but existing computer codes are inadequate for the second. Application of these methods to real chemical problems demands new tools for location and characterisation of saddle-points, intrinsic reaction coordinates, hessians and vibrational frequencies for very large flexible systems. Algorithms capable of performing these tasks have been incorporated in a new software package, GRACE, which provides a non-invasive interface between popular codes for quantum chemistry and molecular dynamics and modelling. Transition structures (TSs) have been refined by this novel procedure, using a combined AM1/CHARMM24/TIP3P potential, involving full gradient relaxation of the positions of 1900-2000 atoms of a solvated enzyme- substrate complex (lactate dehydrogenase/NADH/pyruvate/water). Six different starting structures (arbitrarily selected from a molecular dynamics trajectory for the enzyme-substrate complex) lead to six different TSs. Although the essential features of these TSs are invariant, the relative dispositions of active-site residues differ quite significantly. The transition state for the enzymic reaction would represent an average of the properties of many, nearly degenerate TSs. This insight emerges only as a consequence of the flexible model of the active site employed in this study.

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

U2 - 10.1039/a808364i

DO - 10.1039/a808364i

M3 - Article

AN - SCOPUS:0033560128

VL - 1

SP - 1323

EP - 1331

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 6

ER -