Dataset Information

Active site cysteine is protonated in the PAD4 Michaelis complex: evidence from Born-Oppenheimer ab initio QM/MM molecular dynamics simulations.

ABSTRACT: The protein arginine deiminase 4 (PAD4) catalyzes the citrullination of the peptidylarginine and plays a critical role in rheumatoid arthritis (RA) and gene regulation. Understanding its catalytic mechanism is not only of fundamental importance but also of significant medical interest for the rational design of new inhibitors. By employing on-the-fly Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have demonstrated that it is unlikely for the active site cysteine and histidine to exist as a thiolate-imidazolium ion pair in the PAD4 Michaelis reactant complex. Instead, a substrate-assisted proton transfer mechanism for the deimination reaction step has been characterized: both Cys645 and His471 in the PAD4 active site are neutral prior to the reaction; the deprotonation of Cys645 by the substrate arginine occurs in concert with the nucleophilic addition of the Cys thiolate to Czeta of the substrate, and leads to a covalent tetrahedral intermediate; then, the Czeta-Neta1 bond cleaves and the resulted ammonia is displaced by a solvent water molecule. The initial deprotonation and nucleophilic attack step is found to be rate-determining. The computed free energy barrier with B3LYP(6-31G*) QM/MM MD simulations and umbrella sampling is 20.9 kcal.mol(-1), consistent with the experimental kinetic data. During the deimination, His471 plays an important role in stabilizing the transition state through the formation of the hydrogen bond with the guanidinium group. Our current studies further demonstrated the viability and strength of the ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

SUBMITTER: Ke Z

PROVIDER: S-EPMC2773242 | biostudies-literature | 2009 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Active site cysteine is protonated in the PAD4 Michaelis complex: evidence from Born-Oppenheimer ab initio QM/MM molecular dynamics simulations.

Ke Zhihong Z Zhou Yanzi Y Hu Po P Wang Shenglong S Xie Daiqian D Zhang Yingkai Y

The journal of physical chemistry. B 20090901 38

The protein arginine deiminase 4 (PAD4) catalyzes the citrullination of the peptidylarginine and plays a critical role in rheumatoid arthritis (RA) and gene regulation. Understanding its catalytic mechanism is not only of fundamental importance but also of significant medical interest for the rational design of new inhibitors. By employing on-the-fly Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have demonstrated that it is unlikely for the active site cysteine and histidin ...[more]

PMID: 19507815

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

Active site cysteine is protonated in the PAD4 Michaelis complex: evidence from Born-Oppenheimer ab initio QM/MM molecular dynamics simulations.

Publications

Active site cysteine is protonated in the PAD4 Michaelis complex: evidence from Born-Oppenheimer ab initio QM/MM molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets