Dataset Information

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 due to 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 QM/MM molecular dynamics 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 non-amino-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 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.

SUBMITTER: Wu R

PROVIDER: S-EPMC2812930 | biostudies-literature | 2009 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study.

Wu Ruibo R Hu Po P Wang Shenglong S Cao Zexing Z Zhang Yingkai Y

Journal of chemical theory and computation 20091201 1

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 due to 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 QM/MM molecular dynamics simulation has been employe ...[more]

PMID: 20161624

Similar Datasets

Project description:Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.

Dataset Information

Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study.

Publications

Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets