1
|
Schnee P, Pleiss J, Jeltsch A. Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques. Crit Rev Biochem Mol Biol 2024; 59:20-68. [PMID: 38449437 DOI: 10.1080/10409238.2024.2318547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Accepted: 02/10/2024] [Indexed: 03/08/2024]
Abstract
Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.
Collapse
Affiliation(s)
- Philipp Schnee
- Institute of Biochemistry and Technical Biochemistry, University of Stuttgart, Stuttgart, Germany
| | - Jürgen Pleiss
- Institute of Biochemistry and Technical Biochemistry, University of Stuttgart, Stuttgart, Germany
| | - Albert Jeltsch
- Institute of Biochemistry and Technical Biochemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
2
|
Li Z, Zhang L, Xu K, Jiang Y, Du J, Zhang X, Meng LH, Wu Q, Du L, Li X, Hu Y, Xie Z, Jiang X, Tang YJ, Wu R, Guo RT, Li S. Molecular insights into the catalytic promiscuity of a bacterial diterpene synthase. Nat Commun 2023; 14:4001. [PMID: 37414771 PMCID: PMC10325987 DOI: 10.1038/s41467-023-39706-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2023] [Accepted: 06/19/2023] [Indexed: 07/08/2023] Open
Abstract
Diterpene synthase VenA is responsible for assembling venezuelaene A with a unique 5-5-6-7 tetracyclic skeleton from geranylgeranyl pyrophosphate. VenA also demonstrates substrate promiscuity by accepting geranyl pyrophosphate and farnesyl pyrophosphate as alternative substrates. Herein, we report the crystal structures of VenA in both apo form and holo form in complex with a trinuclear magnesium cluster and pyrophosphate group. Functional and structural investigations on the atypical 115DSFVSD120 motif of VenA, versus the canonical Asp-rich motif of DDXX(X)D/E, reveal that the absent second Asp of canonical motif is functionally replaced by Ser116 and Gln83, together with bioinformatics analysis identifying a hidden subclass of type I microbial terpene synthases. Further structural analysis, multiscale computational simulations, and structure-directed mutagenesis provide significant mechanistic insights into the substrate selectivity and catalytic promiscuity of VenA. Finally, VenA is semi-rationally engineered into a sesterterpene synthase to recognize the larger substrate geranylfarnesyl pyrophosphate.
Collapse
Affiliation(s)
- Zhong Li
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Lilan Zhang
- State Key Laboratory of Biocatalysis and Enzyme Engineering, Hubei Hongshan Laboratory, Hubei Collaborative Innovation Center for Green Transformation of Bio-Resources, Hubei Key Laboratory of Industrial Biotechnology, School of Life Sciences, Hubei University, Wuhan, Hubei, 430062, China
| | - Kangwei Xu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, Guangdong, 510006, China
| | - Yuanyuan Jiang
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Jieke Du
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Xingwang Zhang
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Ling-Hong Meng
- Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Nanhai Road 7, Qingdao, Shandong, 266071, China
- Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266237, China
| | - Qile Wu
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Lei Du
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Xiaoju Li
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Yuechan Hu
- State Key Laboratory of Biocatalysis and Enzyme Engineering, Hubei Hongshan Laboratory, Hubei Collaborative Innovation Center for Green Transformation of Bio-Resources, Hubei Key Laboratory of Industrial Biotechnology, School of Life Sciences, Hubei University, Wuhan, Hubei, 430062, China
| | - Zhenzhen Xie
- State Key Laboratory of Biocatalysis and Enzyme Engineering, Hubei Hongshan Laboratory, Hubei Collaborative Innovation Center for Green Transformation of Bio-Resources, Hubei Key Laboratory of Industrial Biotechnology, School of Life Sciences, Hubei University, Wuhan, Hubei, 430062, China
| | - Xukai Jiang
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Ya-Jie Tang
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China
| | - Ruibo Wu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, Guangdong, 510006, China.
| | - Rey-Ting Guo
- State Key Laboratory of Biocatalysis and Enzyme Engineering, Hubei Hongshan Laboratory, Hubei Collaborative Innovation Center for Green Transformation of Bio-Resources, Hubei Key Laboratory of Industrial Biotechnology, School of Life Sciences, Hubei University, Wuhan, Hubei, 430062, China.
| | - Shengying Li
- State Key Laboratory of Microbial Technology, Shandong University, No. 72 Binhai Road, Qingdao, Shandong, 266237, China.
- Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266237, China.
| |
Collapse
|
3
|
Kamerzell TJ, Mikell B, Chen L, Elias H, Dawn B, MacRae C, Middaugh CR. The structural basis of histone modifying enzyme specificity and promiscuity: Implications for metabolic regulation and drug design. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2022; 130:189-243. [PMID: 35534108 DOI: 10.1016/bs.apcsb.2022.02.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Histone modifying enzymes regulate chromatin architecture through covalent modifications and ultimately control multiple aspects of cellular function. Disruption of histone modification leads to changes in gene expression profiles and may lead to disease. Both small molecule inhibitors and intermediary metabolites have been shown to modulate histone modifying enzyme activity although our ability to identify successful drug candidates or novel metabolic regulators of these enzymes has been limited. Using a combination of large scale in silico screens and in vivo phenotypic analysis, we identified several small molecules and intermediary metabolites with distinctive HME activity. Our approach using unsupervised learning identifies the chemical fingerprints of both small molecules and metabolites that facilitate recognition by the enzymes active sites which can be used as a blueprint to design novel inhibitors. Furthermore, this work supports the idea that histone modifying enzymes sense intermediary metabolites integrating genes, environment and cellular physiology.
Collapse
Affiliation(s)
- Tim J Kamerzell
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, KS, United States; Division of Internal Medicine, HCA MidWest Health, Overland Park, KS, United States; Cardiovascular Division, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States; Division of Cardiovascular Diseases, Cardiovascular Research Institute, University of Kansas Medical Center, Kansas City, KS, United States; Applied AI Technologies, LLC, Overland Park, KS, United States.
| | - Brittney Mikell
- Cardiovascular Division, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States
| | - Lei Chen
- Division of Cardiovascular Diseases, Cardiovascular Research Institute, University of Kansas Medical Center, Kansas City, KS, United States
| | - Harold Elias
- Division of Cardiovascular Diseases, Cardiovascular Research Institute, University of Kansas Medical Center, Kansas City, KS, United States
| | - Buddhadeb Dawn
- Division of Cardiovascular Diseases, Cardiovascular Research Institute, University of Kansas Medical Center, Kansas City, KS, United States
| | - Calum MacRae
- Cardiovascular Division, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States
| | - C Russell Middaugh
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, KS, United States
| |
Collapse
|
4
|
The Impact of Electron Correlation on Describing QM/MM Interactions in the Attendant Molecular Dynamics Simulations of CO in Myoglobin. Sci Rep 2020; 10:8539. [PMID: 32444817 PMCID: PMC7244521 DOI: 10.1038/s41598-020-65475-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 05/05/2020] [Indexed: 01/10/2023] Open
Abstract
The impact of the dispersion and electron correlation effects on describing quantum mechanics/molecular mechanics (QM/MM) interactions in QM/MM molecular dynamics (MD) simulations was explored by performing a series of up to 2 ns QM/MM MD simulations on the B states of the myoglobin-carbon monoxide (MbCO) system. The results indicate that both dispersion and electron correlations play significant roles in the simulation of the ratios of two B states (B1/B2), which suggests that the inclusion of the electron correlation effects is essential for accurately modeling the interactions between QM and MM subsystems. We found that the QM/MM interaction energies between the CO and the surroundings statistically present a linear correlation with the electric fields along the CO bond. This indicates that QM/MM interactions can be described by a simple physical model of a dipole with constant moment under the action of the electric fields. The treatment provides us with an accurate and effective approach to account for the electron correlation effects in QM/MM MD simulations.
Collapse
|
5
|
Czarnota S, Johannissen LO, Baxter NJ, Rummel F, Wilson AL, Cliff MJ, Levy CW, Scrutton NS, Waltho JP, Hay S. Equatorial Active Site Compaction and Electrostatic Reorganization in Catechol- O-methyltransferase. ACS Catal 2019; 9:4394-4401. [PMID: 31080692 PMCID: PMC6503465 DOI: 10.1021/acscatal.9b00174] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Revised: 03/26/2019] [Indexed: 12/18/2022]
Abstract
Catechol-O-methyltransferase (COMT) is a model S-adenosyl-l-methionine (SAM) dependent methyl transferase, which catalyzes the methylation of catecholamine neurotransmitters such as dopamine in the primary pathway of neurotransmitter deactivation in animals. Despite extensive study, there is no consensus view of the physical basis of catalysis in COMT. Further progress requires experimental data that directly probes active site geometry, protein dynamics and electrostatics, ideally in a range of positions along the reaction coordinate. Here we establish that sinefungin, a fungal-derived inhibitor of SAM-dependent enzymes that possess transition state-like charge on the transferring group, can be used as a transition state analog of COMT when combined with a catechol. X-ray crystal structures and NMR backbone assignments of the ternary complexes of the soluble form of human COMT containing dinitrocatechol, Mg2+ and SAM or sinefungin were determined. Comparison and further analysis with the aid of density functional theory calculations and molecular dynamics simulations provides evidence for active site "compaction", which is driven by electrostatic stabilization between the transferring methyl group and "equatorial" active site residues that are orthogonal to the donor-acceptor (pseudo reaction) coordinate. We propose that upon catecholamine binding and subsequent proton transfer to Lys 144, the enzyme becomes geometrically preorganized, with little further movement along the donor-acceptor coordinate required for methyl transfer. Catalysis is then largely facilitated through stabilization of the developing charge on the transferring methyl group via "equatorial" H-bonding and electrostatic interactions orthogonal to the donor-acceptor coordinate.
Collapse
Affiliation(s)
- Sylwia Czarnota
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Linus O. Johannissen
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Nicola J. Baxter
- Krebs
Institute for Biomolecular Research, Department of Molecular Biology
and Biotechnology, The University of Sheffield, Firth Court, Western Bank, Sheffield, S10 2TN, United Kingdom
| | - Felix Rummel
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Alex L. Wilson
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Matthew J. Cliff
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Colin W. Levy
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Nigel S. Scrutton
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Jonathan P. Waltho
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
- Krebs
Institute for Biomolecular Research, Department of Molecular Biology
and Biotechnology, The University of Sheffield, Firth Court, Western Bank, Sheffield, S10 2TN, United Kingdom
| | - Sam Hay
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| |
Collapse
|
6
|
Abstract
Protein lysine methylation is a distinct posttranslational modification that causes minimal changes in the size and electrostatic status of lysine residues. Lysine methylation plays essential roles in regulating fates and functions of target proteins in an epigenetic manner. As a result, substrates and degrees (free versus mono/di/tri) of protein lysine methylation are orchestrated within cells by balanced activities of protein lysine methyltransferases (PKMTs) and demethylases (KDMs). Their dysregulation is often associated with neurological disorders, developmental abnormalities, or cancer. Methyllysine-containing proteins can be recognized by downstream effector proteins, which contain methyllysine reader domains, to relay their biological functions. While numerous efforts have been made to annotate biological roles of protein lysine methylation, limited work has been done to uncover mechanisms associated with this modification at a molecular or atomic level. Given distinct biophysical and biochemical properties of methyllysine, this review will focus on chemical and biochemical aspects in addition, recognition, and removal of this posttranslational mark. Chemical and biophysical methods to profile PKMT substrates will be discussed along with classification of PKMT inhibitors for accurate perturbation of methyltransferase activities. Semisynthesis of methyllysine-containing proteins will also be covered given the critical need for these reagents to unambiguously define functional roles of protein lysine methylation.
Collapse
Affiliation(s)
- Minkui Luo
- Chemical Biology Program , Memorial Sloan Kettering Cancer Center , New York , New York 10065 , United States.,Program of Pharmacology, Weill Graduate School of Medical Science , Cornell University , New York , New York 10021 , United States
| |
Collapse
|
7
|
Fan F, Chen N, Wang Y, Wu R, Cao Z. QM/MM and MM MD Simulations on the Pyrimidine-Specific Nucleoside Hydrolase: A Comprehensive Understanding of Enzymatic Hydrolysis of Uridine. J Phys Chem B 2018; 122:1121-1131. [PMID: 29285933 DOI: 10.1021/acs.jpcb.7b10524] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The pyrimidine-specific nucleoside hydrolase Yeik (CU-NH) from Escherichia coli cleaves the N-glycosidic bond of uridine and cytidine with a 102-104-fold faster rate than that of purine nucleoside substrates, such as inosine. Such a remarkable substrate specificity and the plausible hydrolytic mechanisms of uridine have been explored by using QM/MM and MM MD simulations. The present calculations show that the relatively stronger hydrogen-bond interactions between uridine and the active-site residues Gln227 and Tyr231 in CU-NH play an important role in enhancing the substrate binding and thus promoting the N-glycosidic bond cleavage, in comparison with inosine. The estimated energy barrier of 30 kcal/mol for the hydrolysis of inosine is much higher than 22 kcal/mol for uridine. Extensive MM MD simulations on the transportation of substrates to the active site of CU-NH indicate that the uridine binding is exothermic by ∼23 kcal/mol, more remarkable than inosine (∼12 kcal/mol). All of these arise from the noncovalent interactions between the substrate and the active site featured in CU-NH, which account for the substrate specificity. Quite differing from other nucleoside hydrolases, here the enzymatic N-glycosidic bond cleavage of uridine is less influenced by its protonation.
Collapse
Affiliation(s)
- Fangfang Fan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University , Xiamen 360015, China
| | - Nanhao Chen
- Department of Chemistry, University of California , Davis, California 95616, United States
| | - Yongheng Wang
- School of Pharmaceutical Sciences, Sun Yat-sen University , Guangzhou 510006, China
| | - Ruibo Wu
- School of Pharmaceutical Sciences, Sun Yat-sen University , Guangzhou 510006, China
| | - Zexing Cao
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University , Xiamen 360015, China
| |
Collapse
|
8
|
Linscott JA, Kapilashrami K, Wang Z, Senevirathne C, Bothwell IR, Blum G, Luo M. Kinetic isotope effects reveal early transition state of protein lysine methyltransferase SET8. Proc Natl Acad Sci U S A 2016; 113:E8369-E8378. [PMID: 27940912 PMCID: PMC5206543 DOI: 10.1073/pnas.1609032114] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Protein lysine methyltransferases (PKMTs) catalyze the methylation of protein substrates, and their dysregulation has been linked to many diseases, including cancer. Accumulated evidence suggests that the reaction path of PKMT-catalyzed methylation consists of the formation of a cofactor(cosubstrate)-PKMT-substrate complex, lysine deprotonation through dynamic water channels, and a nucleophilic substitution (SN2) transition state for transmethylation. However, the molecular characters of the proposed process remain to be elucidated experimentally. Here we developed a matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS) method and corresponding mathematic matrix to determine precisely the ratios of isotopically methylated peptides. This approach may be generally applicable for examining the kinetic isotope effects (KIEs) of posttranslational modifying enzymes. Protein lysine methyltransferase SET8 is the sole PKMT to monomethylate histone 4 lysine 20 (H4K20) and its function has been implicated in normal cell cycle progression and cancer metastasis. We therefore implemented the MS-based method to measure KIEs and binding isotope effects (BIEs) of the cofactor S-adenosyl-l-methionine (SAM) for SET8-catalyzed H4K20 monomethylation. A primary intrinsic 13C KIE of 1.04, an inverse intrinsic α-secondary CD3 KIE of 0.90, and a small but statistically significant inverse CD3 BIE of 0.96, in combination with computational modeling, revealed that SET8-catalyzed methylation proceeds through an early, asymmetrical SN2 transition state with the C-N and C-S distances of 2.35-2.40 Å and 2.00-2.05 Å, respectively. This transition state is further supported by the KIEs, BIEs, and steady-state kinetics with the SAM analog Se-adenosyl-l-selenomethionine (SeAM) as a cofactor surrogate. The distinct transition states between protein methyltransferases present the opportunity to design selective transition-state analog inhibitors.
Collapse
Affiliation(s)
- Joshua A Linscott
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065
- Program of Pharmacology, Weill Graduate School of Medical Science, Cornell University, New York, NY 10021
| | - Kanishk Kapilashrami
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Zhen Wang
- Department of Biochemistry, Albert Einstein College of Medicine, Bronx, NY 10461
| | - Chamara Senevirathne
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Ian R Bothwell
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065
- Tri-Institutional PhD Program in Chemical Biology, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Gil Blum
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065
- Tri-Institutional PhD Program in Chemical Biology, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Minkui Luo
- Chemical Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY 10065;
- Program of Pharmacology, Weill Graduate School of Medical Science, Cornell University, New York, NY 10021
| |
Collapse
|
9
|
Zhou Y, Wang S, Li Y, Zhang Y. Born-Oppenheimer Ab Initio QM/MM Molecular Dynamics Simulations of Enzyme Reactions. Methods Enzymol 2016; 577:105-18. [PMID: 27498636 PMCID: PMC4986621 DOI: 10.1016/bs.mie.2016.05.013] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
There are two key requirements for reliably simulating enzyme reactions: one is a reasonably accurate potential energy surface to describe the bond-forming/breaking process as well as to adequately model the heterogeneous enzyme environment; the other is to perform extensive sampling since an enzyme system consists of at least thousands of atoms and its energy landscape is very complex. One attractive approach to meet both daunting tasks is Born-Oppenheimer ab initio QM/MM molecular dynamics (aiQM/MM-MD) simulation with umbrella sampling. In this chapter, we describe our recently developed pseudobond Q-Chem-Amber interface, which employs a combined electrostatic-mechanical embedding scheme with periodic boundary condition and the particle mesh Ewald method for long-range electrostatics interactions. In our implementation, Q-Chem and the sander module of Amber are combined at the source code level without using system calls, and all necessary data communications between QM and MM calculations are achieved via computer memory. We demonstrate the applicability of this pseudobond Q-Chem-Amber interface by presenting two examples, one reaction in aqueous solution and one enzyme reaction. Finally, we describe our established aiQM/MM-MD enzyme simulation protocol, which has been successfully applied to study more than a dozen enzymes.
Collapse
Affiliation(s)
- Y Zhou
- Laboratory of Mesoscopic Chemistry, Collaborative Innovation Center of Chemistry for Life Sciences, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - S Wang
- New York University, New York, NY, United States
| | - Y Li
- International Center of Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai, China
| | - Y Zhang
- New York University, New York, NY, United States; NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China.
| |
Collapse
|
10
|
Mechanisms of histone lysine-modifying enzymes: A computational perspective on the role of the protein environment. J Mol Graph Model 2016; 67:69-84. [DOI: 10.1016/j.jmgm.2016.04.011] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2016] [Revised: 04/28/2016] [Accepted: 04/29/2016] [Indexed: 12/13/2022]
|
11
|
Yao Y, Liu J, Zheng F, Zhan CG. Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine. Theor Chem Acc 2016; 135. [PMID: 28250715 DOI: 10.1007/s00214-015-1788-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.
Collapse
Affiliation(s)
- Yuan Yao
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536; The Academy of Fundamental and Interdisciplinary Sciences, Harbin Institute of Technology, Harbin 150080, P.R. China
| | - Junjun Liu
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536; Tongji School of Pharmacy, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei 430030, P.R. China
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| |
Collapse
|
12
|
Abstract
SAHA (vorinostat, Merck) is a famous clinical drug for zinc-containing histone deacetylase (HDAC) targets against cancer and several other human disorders, whose inhibition mechanism (namely the protonation mechanism) upon binding to HDAC has been debated for more than ten years. It is very challenging to verify experimentally and is still controversial theoretically. The popular "Class-dependent" (namely "Tyr-dependent") hypothesis is that the deprotonation of SAHA is mostly regulated by the conserved Tyr308 in class I HDAC while it is replaced by the His843 in class IIa HDAC. Herein, by elaborate QM(DFT)/MM MD simulations, we exclude the prevalent "Class-dependent" mechanism and advance a novel "Metal-dependent" mechanism, where the remote second metal site (K(+) in most HDAC and Ca(2+) in HDAC2) determines the protonation of SAHA. This proof-of-principle "Metal-dependent" mechanism opens up a new avenue to utilize the second metal site for isoform-selective inhibitor design.
Collapse
Affiliation(s)
- Jingwei Zhou
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, P. R. China.
| | | | | |
Collapse
|
13
|
Gong W, Wu R, Zhang Y. Thiol versus hydroxamate as zinc binding group in HDAC inhibition: An ab initio QM/MM molecular dynamics study. J Comput Chem 2015; 36:2228-35. [PMID: 26452222 DOI: 10.1002/jcc.24203] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Revised: 07/17/2015] [Accepted: 08/30/2015] [Indexed: 12/21/2022]
Abstract
Zinc-dependent histone deacetylases (HDACs) play a critical role in transcriptional repression and gene silencing, and are among the most attractive targets for the development of new therapeutics against cancer and various other diseases. Two HDAC inhibitors have been approved by FDA as anti-cancer drugs: one is SAHA whose hydroxamate is directly bound to zinc, the other is FK228 whose active form may use thiol as the zinc binding group. In spite of extensive studies, it remains to be ambiguous regarding how thiol and hydroxamate are bound to the zinc active site of HDACs. In this work, our computational approaches center on Born-Oppenheimer ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics with umbrella sampling, which allow for modeling of the zinc active site with reasonable accuracy while properly including dynamics and effects of protein environment. Meanwhile, an improved short-long effective function (SLEF2) to describe non-bonded interactions between zinc and other atoms has been employed in initial MM equilibrations. Our ab initio QM/MM MD simulations have confirmed that hydroxamate is neutral when it is bound to HDAC8, and found that thiol is deprotonated when directly bound to zinc in the HDAC active site. By comparing thiol and hydroxamate, our results elucidated the differences in their binding environment in the HDAC active sites, and emphasized the importance of the linker design to achieve more specific binding toward class IIa HDACs.
Collapse
Affiliation(s)
- Wenjing Gong
- Department of Chemistry, New York University, New York, New York, 10003
| | - Ruibo Wu
- School of Pharmaceutical Sciences, Sun Yat-Sen University, Guangzhou, 510006, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York, 10003.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| |
Collapse
|
14
|
Chu Y, Guo H. QM/MM MD and Free Energy Simulation Study of Methyl Transfer Processes Catalyzed by PKMTs and PRMTs. Interdiscip Sci 2015; 7:309-18. [PMID: 26267708 DOI: 10.1007/s12539-015-0280-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2014] [Revised: 09/19/2014] [Accepted: 10/17/2014] [Indexed: 11/27/2022]
Abstract
Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here, we review the results of quantum mechanics/molecular mechanics molecular dynamics and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.
Collapse
Affiliation(s)
- Yuzhuo Chu
- School of Life Science and Biotechnology, Dalian University of Technology, Dalian, 116024, China.
| | - Hong Guo
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN, 37996, USA
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, TN, 37831-6164, USA
| |
Collapse
|
15
|
Chu Y, Guo H. QM/MM MD and free energy simulation study of methyl transfer processes catalyzed by PKMTs and PRMTs. Interdiscip Sci 2015. [PMID: 25595588 DOI: 10.1007/s12539-014-0228-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2014] [Revised: 09/19/2014] [Accepted: 10/17/2014] [Indexed: 09/29/2022]
Abstract
Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here we review the results of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.
Collapse
Affiliation(s)
- Yuzhuo Chu
- School of Life Science and Biotechnology, Dalian University of Technology, Dalian, 116024, China,
| | | |
Collapse
|
16
|
Li Z, Wu Y, Feng LJ, Wu R, Luo HB. Ab Initio QM/MM Study Shows a Highly Dissociated SN2 Hydrolysis Mechanism for the cGMP-Specific Phosphodiesterase-5. J Chem Theory Comput 2014; 10:5448-57. [DOI: 10.1021/ct500761d] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Zhe Li
- School of Pharmaceutical
Sciences, Sun Yat-Sen University, Guangzhou 510006, P. R. China
| | - Yinuo Wu
- School of Pharmaceutical
Sciences, Sun Yat-Sen University, Guangzhou 510006, P. R. China
| | - Ling-Jun Feng
- School of Pharmaceutical
Sciences, Sun Yat-Sen University, Guangzhou 510006, P. R. China
| | - Ruibo Wu
- School of Pharmaceutical
Sciences, Sun Yat-Sen University, Guangzhou 510006, P. R. China
| | - Hai-Bin Luo
- School of Pharmaceutical
Sciences, Sun Yat-Sen University, Guangzhou 510006, P. R. China
| |
Collapse
|
17
|
Zhou J, Xie H, Liu Z, Luo HB, Wu R. Structure–Function Analysis of the Conserved Tyrosine and Diverse π-Stacking among Class I Histone Deacetylases: A QM (DFT)/MM MD Study. J Chem Inf Model 2014; 54:3162-71. [DOI: 10.1021/ci500513n] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Jingwei Zhou
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, 510006 Guangdong, P.R. China
| | - Hujun Xie
- School
of Food Science and Biotechnology, Zhejiang Gongshang University, Hangzhou, 310035 Zhejiang, P.R. China
| | - Zhihong Liu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, 510006 Guangdong, P.R. China
| | - Hai-Bin Luo
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, 510006 Guangdong, P.R. China
| | - Ruibo Wu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, 510006 Guangdong, P.R. China
| |
Collapse
|
18
|
Mlýnský V, Banáš P, Šponer J, van der Kamp MW, Mulholland AJ, Otyepka M. Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme. J Chem Theory Comput 2014; 10:1608-22. [PMID: 26580373 DOI: 10.1021/ct401015e] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics , Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
- CEITEC-Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, U.K
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, U.K
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| |
Collapse
|
19
|
Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration. ENTROPY 2013. [DOI: 10.3390/e16010163] [Citation(s) in RCA: 291] [Impact Index Per Article: 26.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
20
|
Navrátil V, Klusák V, Rulíšek L. Theoretical Aspects of Hydrolysis of Peptide Bonds by Zinc Metalloenzymes. Chemistry 2013; 19:16634-45. [DOI: 10.1002/chem.201302663] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2013] [Revised: 08/08/2013] [Indexed: 12/11/2022]
Affiliation(s)
- Václav Navrátil
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
- Charles University in Prague, Faculty of Science, Department of Biochemistry, Hlavova 2030, 128 43 Praha 2 (Czech Republic)
| | - Vojtěch Klusák
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry, Gilead Sciences & IOCB Research Center, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 166 10 Praha 6 (Czech Republic)
| |
Collapse
|
21
|
Shi Y, Zhou Y, Wang S, Zhang Y. Sirtuin Deacetylation Mechanism and Catalytic Role of the Dynamic Cofactor Binding Loop. J Phys Chem Lett 2013; 4:491-495. [PMID: 23585919 PMCID: PMC3621114 DOI: 10.1021/jz302015s] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Sirtuins constitute a novel family of protein deacetylases and play critical roles in epigenetics, cell death, and metabolism. In spite of numerous experimental studies, the key and most complicated stage of its NAD+-dependent catalytic mechanism remains to be elusive. Herein by employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, a state-of-the-art computational approach to study enzyme reactions, we have characterized the complete deacetylation mechanism for a sirtuin enzyme, determined its multistep free-energy reaction profile, and elucidated essential catalytic roles of the conserved dynamic cofactor binding loop. These new detailed mechanistic insights would facilitate the design of novel mechanism-based sirtuin modulators.
Collapse
Affiliation(s)
- Yawei Shi
- Department of Chemistry, New York University, New York, New York 10003
| | - Yanzi Zhou
- Department of Chemistry, New York University, New York, New York 10003
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, P.R. China
| | - Shenglong Wang
- Department of Chemistry, New York University, New York, New York 10003
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York 10003
| |
Collapse
|
22
|
Recent Progress in Density Functional Methodology for Biomolecular Modeling. STRUCTURE AND BONDING 2013. [DOI: 10.1007/978-3-642-32750-6_1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
|
23
|
Yao Y, Liu J, Zhan CG. Why does the G117H mutation considerably improve the activity of human butyrylcholinesterase against sarin? Insights from quantum mechanical/molecular mechanical free energy calculations. Biochemistry 2012; 51:8980-92. [PMID: 23092211 DOI: 10.1021/bi3009246] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Human butyrylcholinesterase (BChE) is recognized as the most promising bioscavenger for organophosphorus (OP) warfare nerve agents. The G117H mutant of human BChE has been identified as a potential catalytic bioscavenger with a remarkably improved activity against OP nerve agents such as sarin, but it still does not satisfy the clinical use. For further design of the higher-activity mutants against OP nerve agents, it is essential to understand how the G117H mutation improves the activity. The reaction mechanisms and the free energy profiles for spontaneous reactivation of wild-type BChE and its G117H mutant phosphorylated by sarin have been explored, in this study, by performing first-principles quantum mechanical/molecular mechanical free energy calculations, and the remarkable role of the G117H mutation on the activity has been elucidated. For both the wild-type and G117H mutant enzymes, H438 acts as a general base to initiate the spontaneous reactivation that consists of two reaction steps: the nucleophilic attack at the phosphorus by a water molecule and decomposition of the pentacoordinated phosphorus intermediate. The calculated overall free energy barriers, i.e., 30.2 and 23.9 kcal/mol for the wild type and G117H mutant, respectively, are in good agreement with available experimental kinetic data. On the basis of the calculated results, the mutated residue (H117 in the G117H mutant) cannot initiate the spontaneous reactivation as a general base. Instead, it skews the oxyanion hole and makes the phosphorus more open to the nucleophilic water molecule, resulting in a remarkable change in the rate-determining step and significantly improved catalytic activity of human BChE.
Collapse
Affiliation(s)
- Yuan Yao
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, Lexington, Kentucky 40536, United States
| | | | | |
Collapse
|
24
|
Sirin GS, Zhou Y, Lior-Hoffmann L, Wang S, Zhang Y. Aging mechanism of soman inhibited acetylcholinesterase. J Phys Chem B 2012; 116:12199-207. [PMID: 22984913 PMCID: PMC3475498 DOI: 10.1021/jp307790v] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nervous system that hydrolyzes neurotransmitter acetylcholine (ACh) and terminates synaptic signals. The catalytic serine of AChE can be phosphonylated by soman, one of the most potent nerve agents, and subsequently undergo an aging reaction. This phosphonylation and aging process leads to irreversible AChE inhibition, results in accumulation of excess ACh at the synaptic clefts, and causes neuromuscular paralysis. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the aging mechanism of soman phosphonylated AChE and determined its free energy profile. This aging reaction starts with the scission of the O2-Cα bond, which is followed by methyl migration, and results in a tertiary carbenium intermediate. At the transition state, the scissile O2-Cα bond is already cleaved with an average O-C distance of 3.2 ± 0.3 Å and the migrating methyl group is shared between Cα and Cβ carbons with C-C distances of 1.9 ± 0.1 and 1.8 ± 0.1 Å, respectively. The negatively charged phosphonate group is stabilized by a salt bridge with the imidazole ring of the catalytic histidine. A major product of aging, 2,3-dimethyl-2-butanol can be formed swiftly by the reaction of a water molecule. Our characterized mechanism and simulation results provide new detailed insights into this important biochemical process.
Collapse
Affiliation(s)
- Gulseher Sarah Sirin
- Sackler Institute of Graduate Biomedical Sciences, New York University School of Medicine, New York, New York, 10016
- Department of Chemistry, New York University, New York, New York, 10003
| | - Yanzi Zhou
- Department of Chemistry, New York University, New York, New York, 10003
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, P.R. China
| | - Lee Lior-Hoffmann
- Department of Chemistry, New York University, New York, New York, 10003
| | - Shenglong Wang
- Department of Chemistry, New York University, New York, New York, 10003
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York, 10003
| |
Collapse
|
25
|
Rooklin DW, Lu M, Zhang Y. Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase. J Am Chem Soc 2012; 134:15595-603. [PMID: 22928549 PMCID: PMC3461190 DOI: 10.1021/ja307267y] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.
Collapse
Affiliation(s)
- David W. Rooklin
- Department of Chemistry, New York University, New York, NY 10003
| | - Min Lu
- Public Health Research Institute Center, Department of Microbiology and Molecular Genetics, UMDNJ – New Jersey Medical School, Newark, NJ 07103
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003
| |
Collapse
|
26
|
Chu Y, Yao J, Guo H. QM/MM MD and free energy simulations of G9a-like protein (GLP) and its mutants: understanding the factors that determine the product specificity. PLoS One 2012; 7:e37674. [PMID: 22624060 PMCID: PMC3356298 DOI: 10.1371/journal.pone.0037674] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2012] [Accepted: 04/25/2012] [Indexed: 12/02/2022] Open
Abstract
Certain lysine residues on histone tails could be methylated by protein lysine methyltransferases (PKMTs) using S-adenosyl-L-methionine (AdoMet) as the methyl donor. Since the methylation states of the target lysines play a fundamental role in the regulation of chromatin structure and gene expression, it is important to study the property of PKMTs that allows a specific number of methyl groups (one, two or three) to be added (termed as product specificity). It has been shown that the product specificity of PKMTs may be controlled in part by the existence of specific residues at the active site. One of the best examples is a Phe/Tyr switch found in many PKMTs. Here quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed on wild type G9a-like protein (GLP) and its F1209Y and Y1124F mutants for understanding the energetic origin of the product specificity and the reasons for the change of product specificity as a result of single-residue mutations at the Phe/Tyr switch as well as other positions. The free energy barriers of the methyl transfer processes calculated from our simulations are consistent with experimental data, supporting the suggestion that the relative free energy barriers may determine, at least in part, the product specificity of PKMTs. The changes of the free energy barriers as a result of the mutations are also discussed based on the structural information obtained from the simulations. The results suggest that the space and active-site interactions around the ε-amino group of the target lysine available for methyl addition appear to among the key structural factors in controlling the product specificity and activity of PKMTs.
Collapse
Affiliation(s)
- Yuzhuo Chu
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee, United States of America
- National Institute for Mathematical and Biological Synthesis, University of Tennessee, Knoxville, Tennessee, United States of America
| | - Jianzhuang Yao
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee, United States of America
| | - Hong Guo
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee, United States of America
- University of Tennessee/Oak Ridge National Laboratory Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee, United States of America
| |
Collapse
|
27
|
Liu J, Zhan CG. Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine. J Chem Theory Comput 2012; 8:1426-1435. [PMID: 23066354 PMCID: PMC3469271 DOI: 10.1021/ct200810d] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.
Collapse
Affiliation(s)
- Junjun Liu
- Tongji School of Pharmacy, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei 430030, P.R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
28
|
Yao J, Chu Y, An R, Guo H. Understanding Product Specificity of Protein Lysine Methyltransferases from QM/MM Molecular Dynamics and Free Energy Simulations: The Effects of Mutation on SET7/9 beyond the Tyr/Phe Switch. J Chem Inf Model 2012; 52:449-56. [DOI: 10.1021/ci200364m] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Jianzhuang Yao
- Department
of Biochemistry and Cellular and Molecular Biology and ‡National Institute for Mathematical
and Biological Synthesis, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Yuzhuo Chu
- Department
of Biochemistry and Cellular and Molecular Biology and ‡National Institute for Mathematical
and Biological Synthesis, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ran An
- Department
of Biochemistry and Cellular and Molecular Biology and ‡National Institute for Mathematical
and Biological Synthesis, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Hong Guo
- Department
of Biochemistry and Cellular and Molecular Biology and ‡National Institute for Mathematical
and Biological Synthesis, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
29
|
Wu R, Gong W, Liu T, Zhang Y, Cao Z. QM/MM Molecular Dynamics Study of Purine-Specific Nucleoside Hydrolase. J Phys Chem B 2012; 116:1984-91. [DOI: 10.1021/jp211403j] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Ruibo Wu
- School of
Pharmaceutical Sciences,
East Campus, Sun Yat-sen University, Guangzhou
510006, China
- State Key
Laboratory of Physical
Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of
Theoretical and Computational Chemistry, College of Chemistry and
Chemical Engineering, Xiamen University, Xiamen 361005, China
- Department
of Chemistry, New York University, New
York, New York 10003, United
States
| | - Wengjin Gong
- Department
of Chemistry, New York University, New
York, New York 10003, United
States
| | - Ting, Liu
- State Key
Laboratory of Physical
Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of
Theoretical and Computational Chemistry, College of Chemistry and
Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yingkai Zhang
- Department
of Chemistry, New York University, New
York, New York 10003, United
States
| | - Zexing Cao
- State Key
Laboratory of Physical
Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of
Theoretical and Computational Chemistry, College of Chemistry and
Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
30
|
In Silico Strategies Toward Enzyme Function and Dynamics. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2012. [DOI: 10.1016/b978-0-12-398312-1.00009-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
31
|
Smith GK, Ke Z, Guo H, Hengge AC. Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies. J Phys Chem B 2011; 115:13713-22. [PMID: 21999515 DOI: 10.1021/jp207532s] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Phosphorylation reactions catalyzed by kinases and phosphatases play an indispensible role in cellular signaling, and their malfunctioning is implicated in many diseases. A better understanding of the catalytic mechanism will help design novel and effective mechanism-based inhibitors of these enzymes. In this work, ab initio quantum mechanical/molecular mechanical studies are reported for the phosphoryl transfer reaction catalyzed by a cyclin-dependent kinase, CDK2. Our results suggest that an active-site Asp residue, rather than ATP as previously proposed, serves as the general base to activate the Ser nucleophile. The corresponding transition state features a dissociative, metaphosphate-like structure, stabilized by the Mg(2+) ion and several hydrogen bonds. The calculated free-energy barrier is consistent with experimental values. Implications of our results in this and other protein kinases are discussed.
Collapse
Affiliation(s)
- Gregory K Smith
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | | | | | | |
Collapse
|
32
|
Bai Q, Shen Y, Yao X, Wang F, Du Y, Wang Q, Jin N, Hai J, Hu T, Yang J. Modeling a new water channel that allows SET9 to dimethylate p53. PLoS One 2011; 6:e19856. [PMID: 21625555 PMCID: PMC3098259 DOI: 10.1371/journal.pone.0019856] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2011] [Accepted: 04/05/2011] [Indexed: 01/25/2023] Open
Abstract
SET9, a protein lysine methyltransferase, has been thought to be capable of transferring only one methyl group to target lysine residues. However, some reports have pointed out that SET9 can dimethylate Lys372 of p53 (p53-K372) and Lys4 of histone H3 (H3-K4). In order to understand how p53 can be dimethylated by SET9, we measured the radius of the channel that surrounds p53-K372, first on the basis of the crystal structure of SET9, and we show that the channel is not suitable for water movement. Second, molecular dynamic (MD) simulations were carried out for 204 ns on the crystal structure of SET9. The results show that water leaves the active site of SET9 through a new channel, which is made of G292, A295, Y305 and Y335. In addition, the results of molecular docking and MD simulations indicate that the new water channel continues to remain open when S-adenosyl-L-methionine (AdoMet) or S-adenosyl-L-homocysteine (AdoHcy) is bound to SET9. The changes in the radii of these two channels were measured in the equilibrium phase at the constant temperature of 300 K. The results indicate that the first channel still does not allow water to get into or out of the active site, but the new channel is large enough to allow this water to circulate. Our results indicate that water can be removed from the active site, an essential process for allowing the dimethylation reaction to occur.
Collapse
Affiliation(s)
- Qifeng Bai
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
| | - Yulin Shen
- Gansu Computing Center, Lanzhou, Gansu, China
| | - Xiaojun Yao
- School of Chemistry, Lanzhou University, Lanzhou, Gansu, China
| | - Fang Wang
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
| | - Yuping Du
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
| | - Qin Wang
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
| | - Nengzhi Jin
- Gansu Computing Center, Lanzhou, Gansu, China
| | - Jun Hai
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
| | - Tiejun Hu
- Gansu Computing Center, Lanzhou, Gansu, China
- * E-mail: (TH); (JY)
| | - Jinbo Yang
- School of Life Science, Lanzhou University, Lanzhou, Gansu, China
- Department of Molecular Genetics, Lerner Research Institute, The Cleveland Clinic, Cleveland, Ohio, United States of America
- * E-mail: (TH); (JY)
| |
Collapse
|
33
|
Pan K, Deem MW. Predicting Fixation Tendencies of the H3N2 Influenza Virus by Free Energy Calculation. J Chem Theory Comput 2011; 7:1259-1272. [PMID: 21691431 PMCID: PMC3117619 DOI: 10.1021/ct100540p] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Influenza virus evolves to escape from immune system antibodies that bind to it. We used free energy calculations with Einstein crystals as reference states to calculate the difference of antibody binding free energy (ΔΔG) induced by amino acid substitution at each position in epitope B of the H3N2 influenza hemagglutinin, the key target for antibody. A substitution with positive ΔΔG value decreases the antibody binding constant. On average an uncharged to charged amino acid substitution generates the highest ΔΔG values. Also on average, substitutions between small amino acids generate ΔΔG values near to zero. The 21 sites in epitope B have varying expected free energy differences for a random substitution. Historical amino acid substitutions in epitope B for the A/Aichi/2/1968 strain of influenza A show that most fixed and temporarily circulating substitutions generate positive ΔΔG values. We propose that the observed pattern of H3N2 virus evolution is affected by the free energy landscape, the mapping from the free energy landscape to virus fitness landscape, and random genetic drift of the virus. Monte Carlo simulations of virus evolution are presented to support this view.
Collapse
Affiliation(s)
- Keyao Pan
- Department of Bioengineering, Rice University, Houston, Texas 77005
| | - Michael W. Deem
- Department of Bioengineering, Rice University, Houston, Texas 77005
- Department of Physics & Astronomy, Rice University, Houston, Texas 77005
| |
Collapse
|
34
|
Liu J, Zhao X, Yang W, Zhan CG. Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step. J Phys Chem B 2011; 115:5017-25. [PMID: 21486046 PMCID: PMC3087188 DOI: 10.1021/jp200975v] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to examine the catalytic mechanism for cocaine esterase (CocE)-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. It has been shown that the acylation of (+)-cocaine consists of nucleophilic attack of the hydroxyl group of Ser117 on the carbonyl carbon of (+)-cocaine benzoyl ester and the dissociation of (+)-cocaine benzoyl ester. The first reaction step of deacylation of (+)-cocaine, which is identical to that of (-)-cocaine, is rate-determining, indicating that CocE-catalyzed hydrolyses of (+)- and (-)-cocaine have a common rate-determining step. The computational results predict that the catalytic rate constant of CocE against (+)-cocaine should be the same as that of CocE against (-)-cocaine, in contrast with the remarkable difference between human butyrylcholinesterase-catalyzed hydrolyses of (+)- and (-)-cocaine. The prediction has been confirmed by experimental kinetic analysis on CocE-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. The determined common rate-determining step indicates that rational design of a high-activity mutant of CocE should be focused on the first reaction step of the deacylation. Furthermore, the obtained mechanistic insights into the detailed differences in the acylation between the (+)- and (-)-cocaine hydrolyses provide indirect clues for rational design of amino acid mutations that could more favorably stabilize the rate-determining transition state in the deacylation and, thus, improve the catalytic activity of CocE. This study provides a valuable mechanistic base for rational design of an improved esterase for therapeutic treatment of cocaine abuse.
Collapse
Affiliation(s)
| | | | - Wenchao Yang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
35
|
Wu R, Lu Z, Cao Z, Zhang Y. Zinc chelation with hydroxamate in histone deacetylases modulated by water access to the linker binding channel. J Am Chem Soc 2011; 133:6110-3. [PMID: 21456530 PMCID: PMC3086135 DOI: 10.1021/ja111104p] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
It is of significant biological interest and medical importance to develop class- and isoform-selective histone deacetylase (HDAC) modulators. The impact of the linker component on HDAC inhibition specificity has been revealed but is not understood. Using Born-Oppenheimer ab initio QM/MM MD simulations, a state-of-the-art approach to simulating metallo-enzymes, we have found that the hydroxamic acid remains to be protonated upon its binding to HDAC8, and thus disproved the mechanistic hypothesis that the distinct zinc-hydroxamate chelation modes between two HDAC subclasses come from different protonation states of the hydroxamic acid. Instead, our simulations suggest a novel mechanism in which the chelation mode of hydroxamate with the zinc ion in HDACs is modulated by water access to the linker binding channel. This new insight into the interplay between the linker binding and the zinc chelation emphasizes its importance and gives guidance regarding linker design for the development of new class-IIa-specific HDAC inhibitors.
Collapse
Affiliation(s)
- Ruibo Wu
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen 361005, China
| | | | | | | |
Collapse
|
36
|
Yan S, Li T, Yao L. Mutational Effects on the Catalytic Mechanism of Cellobiohydrolase I from Trichoderma reesei. J Phys Chem B 2011; 115:4982-9. [DOI: 10.1021/jp200384m] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Shihai Yan
- Key Lab of Biofuels, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, PR China
| | - Tong Li
- Key Lab of Biofuels, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, PR China
| | - Lishan Yao
- Key Lab of Biofuels, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, PR China
| |
Collapse
|
37
|
Ke Z, Guo H. Ab initio QM/MM free-energy studies of arginine deiminase catalysis: the protonation state of the Cys nucleophile. J Phys Chem B 2011; 115:3725-33. [PMID: 21395290 PMCID: PMC3070061 DOI: 10.1021/jp200843s] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The first step of the hydrolytic deimination of L-arginine catalyzed by arginine deiminase is examined using ab initio quantum mechanical/molecular mechanical molecular dynamics simulations. Two possible protonation states of the nucleophilic Cys406 residue were investigated, and the corresponding activation free energies were obtained via umbrella sampling. Our calculations indicated a reaction free-energy barrier of 21.3 kcal/mol for the neutral cysteine, which is in reasonably good agreement with the experimental k(cat) value of 6.3 s(-1), i.e., a barrier of 16.7 kcal/mol. On the other hand, the deprotonated Cys nucleophile yields a free-energy barrier of 6.7 kcal/mol, much lower than the experimental result. The reaction free-energy barriers along with other data suggest that the Cys nucleophile is dominated by its protonated state in the Michaelis complex, and the reaction barrier corresponds largely to its deprotonation.
Collapse
Affiliation(s)
- Zhihong Ke
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131
| |
Collapse
|
38
|
Siegbahn PE, Himo F. The quantum chemical cluster approach for modeling enzyme reactions. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.13] [Citation(s) in RCA: 166] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
39
|
Wu R, Lu Z, Cao Z, Zhang Y. A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins. J Chem Theory Comput 2011; 7:433-443. [PMID: 21552372 PMCID: PMC3087386 DOI: 10.1021/ct100525r] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Herein we introduce a novel practical strategy to overcome the well-known challenge of modeling the divalent zinc cation in metalloproteins. The main idea is to design short-long effective functions (SLEF) to describe charge interactions between the zinc ion and all other atoms. This SLEF approach has the following desired features: (1). It is pairwise, additive and compatible with widely used atomic pair-wise force fields for modeling biomolecules; (2). It only changes interactions between the zinc ion and other atoms, and does not affect force field parameters that model other interactions in the system; (3). It is a non-bonded model that is inherently capable to describe different zinc ligands and coordination modes. By optimizing two SLEF parameters as well as zinc vdW parameters through force matching based on Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have successfully developed the first SLEF force field (SLEF1) to describe zinc interactions. Extensive molecular dynamics simulations of seven zinc enzyme systems with different coordination ligands and distinct chelation modes (4-,5-,6-fold), including the binuclear zinc active site, yielded zinc coordination numbers and binding distances in good agreement with the corresponding crystal structures as well as ab initio QM/MM MD results. This not only demonstrates the transferability and adequacy of the new SLEF1 force field in describing a variety of zinc proteins, but also indicates that this novel SLEF approach is a promising direction to explore for improving force field description of metal ion interactions.
Collapse
Affiliation(s)
- Ruibo Wu
- Department of Chemistry, New York University, New York, NY 10003 USA
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Zhenyu Lu
- Department of Chemistry, New York University, New York, NY 10003 USA
| | - Zexing Cao
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003 USA
| |
Collapse
|
40
|
Zhou Y, Zhang Y. Serine protease acylation proceeds with a subtle re-orientation of the histidine ring at the tetrahedral intermediate. Chem Commun (Camb) 2010; 47:1577-9. [PMID: 21116528 DOI: 10.1039/c0cc04112b] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The acylation mechanism of a prototypical serine protease trypsin and its complete free energy reaction profile have been determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling.
Collapse
Affiliation(s)
- Yanzi Zhou
- Department of Chemistry, New York University, New York, NY 10003, USA
| | | |
Collapse
|
41
|
Xie HB, Zhou Y, Zhang Y, Johnson JK. Reaction Mechanism of Monoethanolamine with CO2 in Aqueous Solution from Molecular Modeling. J Phys Chem A 2010; 114:11844-52. [DOI: 10.1021/jp107516k] [Citation(s) in RCA: 126] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Hong-Bin Xie
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States, National Energy Technology Laboratory, Pittsburgh, Pennsylvania 15236, United States, and Department of Chemistry, New York University, New York, New York 10003, United States
| | - Yanzi Zhou
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States, National Energy Technology Laboratory, Pittsburgh, Pennsylvania 15236, United States, and Department of Chemistry, New York University, New York, New York 10003, United States
| | - Yingkai Zhang
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States, National Energy Technology Laboratory, Pittsburgh, Pennsylvania 15236, United States, and Department of Chemistry, New York University, New York, New York 10003, United States
| | - J. Karl Johnson
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States, National Energy Technology Laboratory, Pittsburgh, Pennsylvania 15236, United States, and Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
42
|
Zhou Y, Wang S, Zhang Y. Catalytic reaction mechanism of acetylcholinesterase determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations. J Phys Chem B 2010; 114:8817-25. [PMID: 20550161 DOI: 10.1021/jp104258d] [Citation(s) in RCA: 96] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1-2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.
Collapse
Affiliation(s)
- Yanzi Zhou
- Department of Chemistry, New York University, New York, New York 10003, USA
| | | | | |
Collapse
|
43
|
Georgieva P, Himo F. Quantum chemical modeling of enzymatic reactions: the case of histone lysine methyltransferase. J Comput Chem 2010; 31:1707-14. [PMID: 20082388 DOI: 10.1002/jcc.21458] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N-terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods.
Collapse
Affiliation(s)
- Polina Georgieva
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden
| | | |
Collapse
|
44
|
Wu R, Wang S, Zhou N, Cao Z, Zhang Y. A proton-shuttle reaction mechanism for histone deacetylase 8 and the catalytic role of metal ions. J Am Chem Soc 2010; 132:9471-9. [PMID: 20568751 PMCID: PMC2908479 DOI: 10.1021/ja103932d] [Citation(s) in RCA: 101] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Zinc-dependent histone deacetylase 8 (HDAC8) catalyzes the removal of acetyl moieties from histone tails, and is critically involved in regulating chromatin structure and gene expression. The detailed knowledge of its catalytic process is of high importance since it has been established as a most promising target for the development of new antitumor drugs. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations and umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have provided further evidence against the originally proposed general acid-base catalytic pair mechanism for Zinc-dependent histone deacetylases. Instead, our results indicated that HDAC8 employs a proton-shuttle catalytic mechanism, in which a neutral His143 first serves as the general base to accept a proton from the zinc-bound water molecule in the initial rate-determining nucleophilic attack step, and then shuttles it to the amide nitrogen atom to facilitate the cleavage of the amide bond. During the deacetylation process, the Zn(2+) ion changes its coordination mode and plays multiple catalytic roles. For the K(+) ion, which is located about 7 A from the catalytic Zn(2+) ion and conserved in class I and II HDACs, our simulations indicated that its removal would lead to the different transition state structure and a higher free energy reaction barrier for the rate-determining step. It is found that the existence of this conserved K(+) ion would enhance the substrate binding, increase the basicity of His143, strengthen the catalytic role of zinc ion, and improve the transition state stabilization by the enzyme environment.
Collapse
Affiliation(s)
- Ruibo Wu
- Department of Chemistry, New York University, New York, NY 10003 USA
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Shenglong Wang
- Department of Chemistry, New York University, New York, NY 10003 USA
| | - Nengjie Zhou
- Department of Chemistry, New York University, New York, NY 10003 USA
| | - Zexing Cao
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003 USA
| |
Collapse
|
45
|
Zhang X, Bruice TC. Product specificity and mechanism of protein lysine methyltransferases: insights from the histone lysine methyltransferase SET8. Biochemistry 2010; 47:6671-7. [PMID: 18512960 DOI: 10.1021/bi800244s] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations employing a molecular mechanics (MM) force field and hybrid quantum mechanics (QM) and MM (QM/MM) have been carried out to investigate the product specificity and mechanism of the histone H4 lysine 20 (H4-K20) methylation by human histone lysine methyltransferase SET8. At neutral pH, the target lysine is available to only the enzyme in the protonated state. The first step in the methylation reaction must be deprotonation of the lysine target which is followed by the (+)AdoMet methylation of the neutral lysine [Enz.Lys-CH(2)-NH(3)(+).(+)AdoMet --> H(+) + Enz.Lys-CH(2)-NH(2).(+)AdoMet -->--> Enz.Lys-CH(2)-N(Me)H(2)(+).AdoHcy]. The electrostatic interactions between two positive charges on (+)AdoMet and Lys20-NH(3)(+) decrease the pK(a) of Lys20-NH(3)(+). Upon formation of Enz.Lys-NH(3)(+).(+)AdoMet, a water channel by which the proton escapes to the outer solvent phase is formed. The formation of a water channel for the escape of a proton from Lys20-N(Me)H(2)(+) in Enz.Lys20-N(Me)H(2)(+).(+)AdoMet is not formed because the methyl substituent blocks the starting of the water channel. Thus, a second methylation does not take place. The dependence of the occurrence of methyl transfer on the formation of a water channel in SET8 is in accord with our previous reports on product specificity by histone lysine monomethyltransferase SET7/9, large subunit lysine dimethyltransferase (LSMT), and viral histone lysine trimethyltransferase (vSET). The average value of the experimental DeltaG(E)() for the six lysine methyl transfer reactions catalyzed by vSET, LSMT, and SET7/9 with p53 as a substrate is 22.1 +/- 1.0 kcal/mol, and the computed average (DeltaG(C)()) is 22.2 +/- 0.8 kcal/mol. In this study, the computed free energy barrier of the methyl transfer reaction [Lys20-NH(2) + (+)AdoMet --> Lys20-N(Me)H(2)(+) + AdoHcy] catalyzed by SET8 is 20.8 kcal/mol. This is in agreement with the value of 20.6 kcal/mol calculated from the experimental rate constant (0.43 +/- 0.02 min(-1)). Our bond-order computations establish that the H4-K20 monomethylation in SET8 is a concerted linear S(N)2 displacement reaction.
Collapse
Affiliation(s)
- Xiaodong Zhang
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| | | |
Collapse
|
46
|
Ke Z, Wang S, Xie D, Zhang Y. Born-Oppenheimer ab initio QM/MM molecular dynamics simulations of the hydrolysis reaction catalyzed by protein arginine deiminase 4. J Phys Chem B 2009; 113:16705-10. [PMID: 20028143 PMCID: PMC2801900 DOI: 10.1021/jp9080614] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein arginine deiminase 4 (PAD4) catalyzes the citrullination of the peptidylarginine via two successive stages: deimination and hydrolysis. Herein, by employing state-of-the-art Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with the umbrella sampling method, we characterized the catalytic mechanism of the hydrolysis reaction: first, the nucleophilic attack of a water molecule to the C(zeta) of the thiouronium intermediate yields a stable tetrahedral intermediate, and then the S-C(zeta) bond breaks to generate the final product, citrulline. Throughout the hydrolysis reaction, His471 and Asp473 play pivotal catalytic roles by first enhancing the nucleophilic ability of the active water through forming shorter and low-barrier hydrogen bonds and then by serving as proton-accepting groups to deprotonate the water molecule, which is consistent with experimental findings. At the transition state, the spontaneous proton transfer among the reactive water, His471 and Asp473 have been observed. The determined overall free energy barrier for this hydrolysis stage is 16.5 kcal x mol(-1), which is lower than the barrier of 20.9 kcal x mol(-1) for the deimination stage determined previously with the same computational approach [J. Phys. Chem. B 2009, 113, 12750-12758]. Thus, the rate-determining step of the PAD4-catalyzed citrullination is the first step of the deimination. Our current work further demonstrates the strength and applicability of the ab initio QM/MM MD approach in simulating enzyme reactions.
Collapse
Affiliation(s)
- Zhihong Ke
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
- Department of Chemistry, New York University, New York, NY 10003
| | - Shenglong Wang
- Department of Chemistry, New York University, New York, NY 10003
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003
| |
Collapse
|
47
|
Liu J, Zhang Y, Zhan CG. Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase. J Phys Chem B 2009. [PMID: 19924840 DOI: 10.1021/jp903600r] [Citation(s) in RCA: 111] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.
Collapse
Affiliation(s)
- Junjun Liu
- Key Laboratory of Pesticide and Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, People's Republic of China
| | | | | |
Collapse
|
48
|
Liu J, Zhang Y, Zhan CG. Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase. J Phys Chem B 2009; 113:16226-36. [PMID: 19924840 PMCID: PMC2795121 DOI: 10.1021/jp9055335] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.
Collapse
Affiliation(s)
- Junjun Liu
- Key Laboratory of Pesticide and Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P. R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| | - Yingkai Zhang
- Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| |
Collapse
|
49
|
Wu R, Hu P, Wang S, Cao Z, Zhang Y. Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born-Oppenheimer ab initio QM/MM Molecular Dynamics Study. J Chem Theory Comput 2009; 6:337. [PMID: 20161624 PMCID: PMC2812930 DOI: 10.1021/ct9005322] [Citation(s) in RCA: 72] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
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.
Collapse
Affiliation(s)
- Ruibo Wu
- Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003 USA
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Po Hu
- Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003 USA
| | - Shenglong Wang
- Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003 USA
| | - Zexing Cao
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003 USA
| |
Collapse
|
50
|
Zheng H, Wang S, Zhang Y. Increasing the time step with mass scaling in Born-Oppenheimer ab initio QM/MM molecular dynamics simulations. J Comput Chem 2009; 30:2706-11. [PMID: 19399770 PMCID: PMC2773245 DOI: 10.1002/jcc.21296] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Born-Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state-of-the-art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass-scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass-scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles.
Collapse
Affiliation(s)
- Han Zheng
- Department of Chemistry, New York University, New York, New York 10003, USA
| | | | | |
Collapse
|