Flexibility of catalytic zinc coordination in thermolysin and HDAC8

A Born-Oppenheimer ab initio QM/MM molecular dynamics study

Ruibo Wu, Po Hu, Shenglong Wang, Zexing Cao, Yingkai Zhang

Research output: Contribution to journalArticle

Abstract

The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly because of the well-known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio Quantum Mechanical/Molecular Mechanical (QM/MM) molecular dynamics (MD) simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in thermolysin (TLN) and histone deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time scale. For TLN, the coordination of the carboxylate group of Glu166 to zinc is found to continuously change between monodentate and bidentate manner dynamically, while for HDAC8, the flexibility mainly comes from the coordination to a nonamino acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of the amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.

Original languageEnglish (US)
Pages (from-to)337-343
Number of pages7
JournalJournal of Chemical Theory and Computation
Volume6
Issue number1
DOIs
StatePublished - Jan 12 2010

Fingerprint

Thermolysin
Histone Deacetylases
Molecular dynamics
Zinc
flexibility
zinc
molecular dynamics
Ligands
ligands
Metals
simulation
Computer simulation
Amides
viability
metals
carboxylates
amides
Catalysis
catalysis
hydrolysis

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Computer Science Applications

Cite this

Flexibility of catalytic zinc coordination in thermolysin and HDAC8 : A Born-Oppenheimer ab initio QM/MM molecular dynamics study. / Wu, Ruibo; Hu, Po; Wang, Shenglong; Cao, Zexing; Zhang, Yingkai.

In: Journal of Chemical Theory and Computation, Vol. 6, No. 1, 12.01.2010, p. 337-343.

Research output: Contribution to journalArticle

@article{e5bb14b4d9794567b1b9841756918d2a,
title = "Flexibility of catalytic zinc coordination in thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM molecular dynamics study",
abstract = "The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly because of the well-known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio Quantum Mechanical/Molecular Mechanical (QM/MM) molecular dynamics (MD) simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in thermolysin (TLN) and histone deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time scale. For TLN, the coordination of the carboxylate group of Glu166 to zinc is found to continuously change between monodentate and bidentate manner dynamically, while for HDAC8, the flexibility mainly comes from the coordination to a nonamino acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of the amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.",
author = "Ruibo Wu and Po Hu and Shenglong Wang and Zexing Cao and Yingkai Zhang",
year = "2010",
month = "1",
day = "12",
doi = "10.1021/ct9005322",
language = "English (US)",
volume = "6",
pages = "337--343",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "1",

}

TY - JOUR

T1 - Flexibility of catalytic zinc coordination in thermolysin and HDAC8

T2 - A Born-Oppenheimer ab initio QM/MM molecular dynamics study

AU - Wu, Ruibo

AU - Hu, Po

AU - Wang, Shenglong

AU - Cao, Zexing

AU - Zhang, Yingkai

PY - 2010/1/12

Y1 - 2010/1/12

N2 - The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly because of the well-known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio Quantum Mechanical/Molecular Mechanical (QM/MM) molecular dynamics (MD) simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in thermolysin (TLN) and histone deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time scale. For TLN, the coordination of the carboxylate group of Glu166 to zinc is found to continuously change between monodentate and bidentate manner dynamically, while for HDAC8, the flexibility mainly comes from the coordination to a nonamino acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of the amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.

AB - The different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly because of the well-known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio Quantum Mechanical/Molecular Mechanical (QM/MM) molecular dynamics (MD) simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in thermolysin (TLN) and histone deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's time scale. For TLN, the coordination of the carboxylate group of Glu166 to zinc is found to continuously change between monodentate and bidentate manner dynamically, while for HDAC8, the flexibility mainly comes from the coordination to a nonamino acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of the amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.

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

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

U2 - 10.1021/ct9005322

DO - 10.1021/ct9005322

M3 - Article

VL - 6

SP - 337

EP - 343

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 1

ER -