1
|
Movilla S, Roca M, Moliner V, Magistrato A. Molecular Basis of RNA-Driven ATP Hydrolysis in DExH-Box Helicases. J Am Chem Soc 2023; 145:6691-6701. [PMID: 36926902 DOI: 10.1021/jacs.2c11980] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Abstract
The spliceosome machinery catalyzes precursor messenger (pre-m)RNA splicing. In each cycle, the spliceosome experiences massive compositional and conformational remodeling fueled by the concerted action of specific RNA-dependent ATPases/helicases. Intriguingly, these enzymes are allosterically activated to perform ATP hydrolysis and trigger helicase activity only upon pre-mRNA binding. Yet, the molecular mechanism underlying the RNA-driven regulation of their ATPase function remains elusive. Here, we focus on the Prp2 ATPase/helicase which contributes to reshaping the spliceosome into its catalytic competent state. By performing classical and quantum-classical molecular dynamics simulations, we unprecedentedly unlock the molecular terms governing the Prp2 ATPase/helicase function. Namely, we dissect the molecular mechanism of ATP hydrolysis, and we disclose that RNA binding allosterically triggers the formation of a set of interactions linking the RNA binding tunnel to the catalytic site. This activates the Prp2's ATPase function by optimally placing the nucleophilic water and the general base of the enzymatic process to perform ATP hydrolysis. The key structural motifs, mechanically coupling RNA gripping and the ATPase/helicase functions, are conserved across all DExH-box helicases. This mechanism could thus be broadly applicable to all DExH-box helicase family.
Collapse
Affiliation(s)
- Santiago Movilla
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| | - Maite Roca
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| | - Vicent Moliner
- BioComp Group, Institute of Advanced Materials (INAM), Universitat Jaume I, 12071 Castellón, Spain
| | - Alessandra Magistrato
- Department National Research Council of Italy (CNR), Institute of Material (IOM) c/o International School for Advanced Studies (SISSA), 34136 Trieste, Italy
| |
Collapse
|
2
|
Recabarren R, Zinovjev K, Tuñón I, Alzate-Morales J. How a Second Mg 2+ Ion Affects the Phosphoryl-Transfer Mechanism in a Protein Kinase: A Computational Study. ACS Catal 2020. [DOI: 10.1021/acscatal.0c03304] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Rodrigo Recabarren
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, 1 Poniente, 1141 Talca, Chile
| | - Kirill Zinovjev
- School of Biochemistry, University of Bristol, Biomedical Sciences Building, University Walk, Bristol BS8 1TD, U.K
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, Valencia 46010, Spain
| | - Jans Alzate-Morales
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, 1 Poniente, 1141 Talca, Chile
| |
Collapse
|
3
|
Erxleben A. Mechanistic Studies of Homo- and Heterodinuclear Zinc Phosphoesterase Mimics: What Has Been Learned? Front Chem 2019; 7:82. [PMID: 30847339 PMCID: PMC6393734 DOI: 10.3389/fchem.2019.00082] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2018] [Accepted: 01/30/2019] [Indexed: 11/13/2022] Open
Abstract
Phosphoesterases hydrolyze the phosphorus oxygen bond of phosphomono-, di- or triesters and are involved in various important biological processes. Carboxylate and/or hydroxido-bridged dizinc(II) sites are a widespread structural motif in this enzyme class. Much effort has been invested to unravel the mechanistic features that provide the enormous rate accelerations observed for enzymatic phosphate ester hydrolysis and much has been learned by using simple low-molecular-weight model systems for the biological dizinc(II) sites. This review summarizes the knowledge and mechanistic understanding of phosphoesterases that has been gained from biomimetic dizinc(II) complexes, showing the power as well as the limitations of model studies.
Collapse
Affiliation(s)
- Andrea Erxleben
- School of Chemistry, National University of Ireland Galway, Galway, Ireland
| |
Collapse
|
4
|
Deazapurine Analogues Bearing a 1 H-Pyrazolo[3,4- b]pyridin-3(2 H)-one Core: Synthesis and Biological Activity. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
5
|
Chhetri N, Bhoite SA, Singh AK. Micellar catalyzed hydrolysis of mono-2,3-dichloroaniline phosphate. J DISPER SCI TECHNOL 2018. [DOI: 10.1080/01932691.2017.1379022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Nisha Chhetri
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur, Chhattisgarh, India
| | - S. A. Bhoite
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur, Chhattisgarh, India
| | - A. K. Singh
- Department of Chemistry, Government V. Y. T. P.G. Autonomous College, Durg, Chhattisgarh, India
| |
Collapse
|
6
|
Purg M, Kamerlin SCL. Empirical Valence Bond Simulations of Organophosphate Hydrolysis: Theory and Practice. Methods Enzymol 2018; 607:3-51. [DOI: 10.1016/bs.mie.2018.06.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
7
|
Zhang H, Yang L, Ding W, Ma Y. The pH-dependent activation mechanism of Ser102 in Escherichia coli alkaline phosphatase: a theoretical study. J Biol Inorg Chem 2017; 23:277-284. [DOI: 10.1007/s00775-017-1529-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2017] [Accepted: 12/17/2017] [Indexed: 12/22/2022]
|
8
|
Roston D, Cui Q. Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects. J Am Chem Soc 2016; 138:11946-57. [PMID: 27541005 PMCID: PMC5705178 DOI: 10.1021/jacs.6b07347] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.
Collapse
Affiliation(s)
- Daniel Roston
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin , Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin , Madison, Wisconsin 53706, United States
| |
Collapse
|
9
|
Roston D, Demapan D, Cui Q. Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site. J Am Chem Soc 2016; 138:7386-94. [PMID: 27186960 DOI: 10.1021/jacs.6b03156] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A reaction's transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.
Collapse
Affiliation(s)
- Daniel Roston
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin , Madison, Wisconsin 53706, United States
| | - Darren Demapan
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin , Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin , Madison, Wisconsin 53706, United States
| |
Collapse
|
10
|
Pabis A, Duarte F, Kamerlin SCL. Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer. Biochemistry 2016; 55:3061-81. [PMID: 27187273 PMCID: PMC4899807 DOI: 10.1021/acs.biochem.6b00297] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
The
enzymes that facilitate phosphate and sulfate hydrolysis are
among the most proficient natural catalysts known to date. Interestingly,
a large number of these enzymes are promiscuous catalysts that exhibit
both phosphatase and sulfatase activities in the same active site
and, on top of that, have also been demonstrated to efficiently catalyze
the hydrolysis of other additional substrates with varying degrees
of efficiency. Understanding the factors that underlie such multifunctionality
is crucial both for understanding functional evolution in enzyme superfamilies
and for the development of artificial enzymes. In this Current Topic,
we have primarily focused on the structural and mechanistic basis
for catalytic promiscuity among enzymes that facilitate both phosphoryl
and sulfuryl transfer in the same active site, while comparing this
to how catalytic promiscuity manifests in other promiscuous phosphatases.
We have also drawn on the large number of experimental and computational
studies of selected model systems in the literature to explore the
different features driving the catalytic promiscuity of such enzymes.
Finally, on the basis of this comparative analysis, we probe the plausible
origins and determinants of catalytic promiscuity in enzymes that
catalyze phosphoryl and sulfuryl transfer.
Collapse
Affiliation(s)
- Anna Pabis
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University , BMC Box 596, S-751 24 Uppsala, Sweden
| | - Fernanda Duarte
- Chemistry Research Laboratory, University of Oxford , 12 Mansfield Road, Oxford OX1 3TA, U.K.,Physical and Theoretical Chemistry Laboratory, University of Oxford , South Parks Road, Oxford OX1 3QZ, U.K
| | - Shina C L Kamerlin
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University , BMC Box 596, S-751 24 Uppsala, Sweden
| |
Collapse
|
11
|
|
12
|
Pabis A, Kamerlin SCL. Promiscuity and electrostatic flexibility in the alkaline phosphatase superfamily. Curr Opin Struct Biol 2015; 37:14-21. [PMID: 26716576 DOI: 10.1016/j.sbi.2015.11.008] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2015] [Revised: 11/09/2015] [Accepted: 11/17/2015] [Indexed: 12/17/2022]
Abstract
Catalytic promiscuity, that is, the ability of single enzymes to facilitate the turnover of multiple, chemically distinct substrates, is a widespread phenomenon that plays an important role in the evolution of enzyme function. Additionally, such pre-existing multifunctionality can be harnessed in artificial enzyme design. The members of the alkaline phosphatase superfamily have served extensively as both experimental and computational model systems for enhancing our understanding of catalytic promiscuity. In this Opinion, we present key recent computational studies into the catalytic activity of these highly promiscuous enzymes, highlighting the valuable insight they have provided into both the molecular basis for catalytic promiscuity in general, and its implications for the evolution of phosphatase activity.
Collapse
Affiliation(s)
- Anna Pabis
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Shina Caroline Lynn Kamerlin
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden.
| |
Collapse
|
13
|
Liao RZ, Siegbahn PEM. Phosphate Hydrolysis by the Fe2–Ca3-Dependent Alkaline Phosphatase PhoX: Mechanistic Insights from DFT calculations. Inorg Chem 2015; 54:11941-7. [DOI: 10.1021/acs.inorgchem.5b02268] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rong-Zhen Liao
- Key Laboratory of Material Chemistry for
Energy Conversion and Storage, Ministry of Education, School of Chemistry
and Chemical Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Per E. M. Siegbahn
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691 Stockholm, Sweden
| |
Collapse
|
14
|
Bordes I, Ruiz-Pernía JJ, Castillo R, Moliner V. A computational study of the phosphoryl transfer reaction between ATP and Dha in aqueous solution. Org Biomol Chem 2015; 13:10179-90. [PMID: 26303076 DOI: 10.1039/c5ob01079a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Phosphoryl transfer reactions are ubiquitous in biology, being involved in processes ranging from energy and signal transduction to the replication genetic material. Dihydroxyacetone phosphate (Dha-P), an intermediate of the synthesis of pyruvate and a very important building block in nature, can be generated by converting free dihydroxyacetone (Dha) through the action of the dihydroxyacetone kinase enzyme. In this paper the reference uncatalyzed reaction in solution has been studied in order to define the foundations of the chemical reaction and to determine the most adequate computational method to describe this electronically complex reaction. In particular, the phosphorylation reaction mechanism between adenosine triphosphate (ATP) and Dha in aqueous solution has been studied by means of quantum mechanics/molecular mechanics (QM/MM) Molecular Dynamics (MD) simulations with the QM subset of atoms described with semi-empirical and DFT methods. The results appear to be strongly dependent on the level of calculation, which will have to be taken into account for future studies of the reaction catalyzed by enzymes. In particular, PM3/MM renders lower free energy barriers and a less endergonic process than AM1d/MM and PM6/MM methods. Nevertheless, the concerted pathway was not located with the former combination of potentials.
Collapse
Affiliation(s)
- I Bordes
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | | | | | | |
Collapse
|
15
|
Sgrignani J, Magistrato A. QM/MM MD Simulations on the Enzymatic Pathway of the Human Flap Endonuclease (hFEN1) Elucidating Common Cleavage Pathways to RNase H Enzymes. ACS Catal 2015. [DOI: 10.1021/acscatal.5b00178] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Jacopo Sgrignani
- Institute of Research in Biomedicine (IRB), Via Vincenzo Vela, 6500 Bellinzona, Switzerland
| | - Alessandra Magistrato
- CNR-IOM-Democritos
National Simulation Center c/o International School for Advanced Studies (SISSA/ISAS), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
16
|
Theoretical studies of the hydrolysis of antibiotics catalyzed by a metallo-β-lactamase. Arch Biochem Biophys 2015; 582:116-26. [PMID: 25622886 DOI: 10.1016/j.abb.2015.01.013] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2014] [Revised: 01/13/2015] [Accepted: 01/14/2015] [Indexed: 12/27/2022]
Abstract
In this paper, hybrid QM/MM molecular dynamics (MD) simulations have been performed to explore the mechanisms of hydrolysis of two antibiotics, Imipenen (IMI), an antibiotic belonging to the subgroup of carbapenems, and the Cefotaxime (CEF), a third-generation cephalosporin antibiotic, in the active site of a mono-nuclear β-lactamase, CphA from Aeromonas hydrophila. Significant different transition state structures are obtained for the hydrolysis of both antibiotics: while the TS of the CEF is an ionic species with negative charge on nitrogen, the IMI TS presents a tetrahedral-like character with negative charge on oxygen atom of the carbonyl group of the lactam ring. Thus, dramatic conformational changes can take place in the cavity of CphA to accommodate different substrates, which would be the origin of its substrate promiscuity. Since CphA shows only activity against carbapenem antibiotic, this study sheds some light into the origin of the selectivity of the different MbL and, as a consequence, into the discovery of specific and potent MβL inhibitors against a broad spectrum of bacterial pathogens. We have finally probed that a re-parametrization of semiempirical methods should be done to properly describe the behavior the metal cation in active site, Zn(2+), when used in QM/MM calculations.
Collapse
|
17
|
Sánchez-Tarín M, Swiderek K, Roca M, Tuñón I. Enzyme Promiscuity in Enolase Superfamily. Theoretical Study of o-Succinylbenzoate Synthase Using QM/MM Methods. J Phys Chem B 2015; 119:1899-911. [DOI: 10.1021/jp511147b] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- María Sánchez-Tarín
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Katarzyna Swiderek
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
- Institute
of Applied Radiation Chemistry, Lodz University of Technology, 90-924, Lodz, Poland
| | - Maite Roca
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
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
|
Chen SL, Liao RZ. Phosphate monoester hydrolysis by trinuclear alkaline phosphatase; DFT study of transition States and reaction mechanism. Chemphyschem 2014; 15:2321-30. [PMID: 24683174 DOI: 10.1002/cphc.201402016] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2014] [Revised: 02/25/2014] [Indexed: 12/23/2022]
Abstract
Alkaline phosphatase (AP) is a trinuclear metalloenzyme that catalyzes the hydrolysis of a broad range of phosphate monoesters to form inorganic phosphate and alcohol (or phenol). In this paper, by using density functional theory with a model based on a crystal structure, the AP-catalyzed hydrolysis of phosphate monoesters is investigated by calculating two substrates, that is, methyl and p-nitrophenyl phosphates, which represent alkyl and aryl phosphates, respectively. The calculations confirm that the AP reaction employs a "ping-pong" mechanism involving two chemical displacement steps, that is, the displacement of the substrate leaving group by a Ser102 alkoxide and the hydrolysis of the phosphoseryl intermediate by a Zn2-bound hydroxide. Both displacement steps proceed via a concerted associative pathway no matter which substrate is used. Other mechanistic aspects are also studied. Comparison of our calculations with linear free energy relationships experiments shows good agreement.
Collapse
Affiliation(s)
- Shi-Lu Chen
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic, Conversion Materials, School of Chemistry, Beijing Institute of Technology, Beijing 100081 (China), Fax: (+86) 01-6891-1354.
| | | |
Collapse
|
20
|
Blomberg MRA, Borowski T, Himo F, Liao RZ, Siegbahn PEM. Quantum chemical studies of mechanisms for metalloenzymes. Chem Rev 2014; 114:3601-58. [PMID: 24410477 DOI: 10.1021/cr400388t] [Citation(s) in RCA: 436] [Impact Index Per Article: 43.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Margareta R A Blomberg
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University , SE-106 91 Stockholm, Sweden
| | | | | | | | | |
Collapse
|
21
|
Density functional calculations on alcoholysis and thiolysis of phosphate triesters: Stepwise or concerted? COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2013.05.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
22
|
Retegan M, Neese F, Pantazis DA. Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II. J Chem Theory Comput 2013; 9:3832-42. [PMID: 26584129 DOI: 10.1021/ct400477j] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.
Collapse
Affiliation(s)
- Marius Retegan
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-38, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-38, 45470 Mülheim an der Ruhr, Germany
| | - Dimitrios A Pantazis
- Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-38, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
23
|
Rudbeck ME, Blomberg MRA, Barth A. Hydrolysis of the E2P Phosphoenzyme of the Ca2+-ATPase: A Theoretical Study. J Phys Chem B 2013; 117:9224-32. [DOI: 10.1021/jp4049814] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Maria E. Rudbeck
- Department of Biochemistry and
Biophysics, Arrhenius Laboratories, Stockholm University, 10691 Stockholm, Sweden
| | - Margareta R. A. Blomberg
- Department of Organic Chemistry,
Arrhenius Laboratories, Stockholm University, 10691 Stockholm, Sweden
| | - Andreas Barth
- Department of Biochemistry and
Biophysics, Arrhenius Laboratories, Stockholm University, 10691 Stockholm, Sweden
| |
Collapse
|
24
|
Duarte F, Amrein BA, Kamerlin SCL. Modeling catalytic promiscuity in the alkaline phosphatase superfamily. Phys Chem Chem Phys 2013; 15:11160-77. [PMID: 23728154 PMCID: PMC3693508 DOI: 10.1039/c3cp51179k] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2013] [Accepted: 05/02/2013] [Indexed: 12/19/2022]
Abstract
In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.
Collapse
Affiliation(s)
- Fernanda Duarte
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | - Beat Anton Amrein
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | | |
Collapse
|
25
|
Hou G, Cui Q. Stabilization of different types of transition states in a single enzyme active site: QM/MM analysis of enzymes in the alkaline phosphatase superfamily. J Am Chem Soc 2013; 135:10457-69. [PMID: 23786365 PMCID: PMC3759165 DOI: 10.1021/ja403293d] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The first step for the hydrolysis of a phosphate monoester (pNPP(2-)) in enzymes of the alkaline phosphatase (AP) superfamily, R166S AP and wild-type NPP, is studied using QM/MM simulations based on an approximate density functional theory (SCC-DFTBPR) and a recently introduced QM/MM interaction Hamiltonian. The calculations suggest that similar loose transition states are involved in both enzymes, despite the fact that phosphate monoesters are the cognate substrates for AP but promiscuous substrates for NPP. The computed loose transition states are clearly different from the more synchronous ones previously calculated for diester reactions in the same AP enzymes. Therefore, our results explicitly support the proposal that AP enzymes are able to recognize and stabilize different types of transition states in a single active site. Analysis of the structural features of computed transition states indicates that the plastic nature of the bimetallic site plays a minor role in accommodating multiple types of transition states and that the high degree of solvent accessibility of the AP active site also contributes to its ability to stabilize diverse transition-state structures without the need of causing large structural distortions of the bimetallic motif. The binding mode of the leaving group in the transition state highlights that vanadate may not always be an ideal transition state analog for loose phosphoryl transfer transition states.
Collapse
Affiliation(s)
- Guanhua Hou
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706
| |
Collapse
|
26
|
Radak BK, Harris ME, York DM. Molecular simulations of RNA 2'-O-transesterification reaction models in solution. J Phys Chem B 2013; 117:94-103. [PMID: 23214417 PMCID: PMC3574632 DOI: 10.1021/jp3084277] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
We employ quantum mechanical/molecular mechanical umbrella sampling simulations to probe the free energy surfaces of a series of increasingly complex reaction models of RNA 2'-O-transesterification in aqueous solution under alkaline conditions. Such models are valuable for understanding the uncatalyzed processes underlying catalytic cleavage of the phosphodiester backbone of RNA, a reaction of fundamental importance in biology. The chemically reactive atoms are modeled by the AM1/d-PhoT quantum model for phosphoryl transfer, whereas the aqueous solvation environment is modeled with a molecular mechanics force field. Several simulation protocols were compared that used different ionic conditions and force field models. The results provide insight into how variation of the structural environment of the nucleophile and leaving group affects the free energy profile for the transesterification reaction. Results for a simple RNA backbone model are compared with recent experiments by Harris et al. on the specific base-catalyzed cleavage of a UpG dinucleotide. The calculated and measured free energies of activation match extremely well (ΔF(‡) = 19.9-20.8 vs 19.9 kcal/mol). Solvation is seen to play a crucial role and is characterized by a network of hydrogen bonds that envelopes the pentacoordinate dianionic phosphorane transition state and provides preferential stabilization relative to the reactant state.
Collapse
Affiliation(s)
- Brian K. Radak
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431 USA
| | - Michael E. Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, OH 44106 USA
| | - Darrin M. York
- BioMaPS Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
| |
Collapse
|
27
|
Navarro-Whyte L, Kellie JL, Lenz SAP, Wetmore SD. Hydrolysis of the damaged deoxythymidine glycol nucleoside and comparison to canonical DNA. Phys Chem Chem Phys 2013; 15:19343-52. [DOI: 10.1039/c3cp53217h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Lee TS, Radak BK, Pabis A, York DM. A New Maximum Likelihood Approach for Free Energy Profile Construction from Molecular Simulations. J Chem Theory Comput 2012; 9:153-164. [PMID: 23457427 DOI: 10.1021/ct300703z] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A novel variational method for construction of free energy profiles from molecular simulation data is presented. The variational free energy profile (VFEP) method uses the maximum likelihood principle applied to the global free energy profile based on the entire set of simulation data (e.g from multiple biased simulations) that spans the free energy surface. The new method addresses common obstacles in two major problems usually observed in traditional methods for estimating free energy surfaces: the need for overlap in the re-weighting procedure and the problem of data representation. Test cases demonstrate that VFEP outperforms other methods in terms of the amount and sparsity of the data needed to construct the overall free energy profiles. For typical chemical reactions, only ~5 windows and ~20-35 independent data points per window are sufficient to obtain an overall qualitatively correct free energy profile with sampling errors an order of magnitude smaller than the free energy barrier. The proposed approach thus provides a feasible mechanism to quickly construct the global free energy profile and identify free energy barriers and basins in free energy simulations via a robust, variational procedure that determines an analytic representation of the free energy profile without the requirement of numerically unstable histograms or binning procedures. It can serve as a new framework for biased simulations and is suitable to be used together with other methods to tackle with the free energy estimation problem.
Collapse
Affiliation(s)
- Tai-Sung Lee
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | | | | | | |
Collapse
|
29
|
Zienau J, Cui Q. Implementation of the Solvent Macromolecule Boundary Potential and Application to Model and Realistic Enzyme Systems. J Phys Chem B 2012; 116:12522-34. [DOI: 10.1021/jp308218m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Jan Zienau
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, 1101 University Avenue,
Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, 1101 University Avenue,
Madison, Wisconsin 53706, United States
| |
Collapse
|
30
|
López-Canut V, Ruiz-Pernía JJ, Castillo R, Moliner V, Tuñón I. Hydrolysis of Phosphotriesters: A Theoretical Analysis of the Enzymatic and Solution Mechanisms. Chemistry 2012; 18:9612-21. [DOI: 10.1002/chem.201103615] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2011] [Revised: 04/26/2012] [Indexed: 11/06/2022]
|
31
|
Liao RZ, Thiel W. Comparison of QM-Only and QM/MM Models for the Mechanism of Tungsten-Dependent Acetylene Hydratase. J Chem Theory Comput 2012; 8:3793-803. [DOI: 10.1021/ct3000684] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rong-Zhen Liao
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470, Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470, Mülheim an der Ruhr, Germany
| |
Collapse
|
32
|
Hou G, Cui Q. QM/MM analysis suggests that Alkaline Phosphatase (AP) and nucleotide pyrophosphatase/phosphodiesterase slightly tighten the transition state for phosphate diester hydrolysis relative to solution: implication for catalytic promiscuity in the AP superfamily. J Am Chem Soc 2012; 134:229-46. [PMID: 22097879 PMCID: PMC3257412 DOI: 10.1021/ja205226d] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Several members of the Alkaline Phosphatase (AP) superfamily exhibit a high level of catalytic proffciency and promiscuity in structurally similar active sites. A thorough characterization of the nature of transition state for different substrates in these enzymes is crucial for understanding the molecular mechanisms that govern those remarkable catalytic properties. In this work, we study the hydrolysis of a phosphate diester, MpNPP(-), in solution, two experimentally well-characterized variants of AP (R166S AP, R166S/E322Y AP) and wild type Nucleotide pyrophosphatase/phosphodiesterase (NPP) by QM/MM calculations in which the QM method is an approximate density functional theory previously parametrized for phosphate hydrolysis (SCC-DFTBPR). The general agreements found between these calculations and available experimental data for both solution and enzymes support the use of SCC-DFTBPR/MM for a semiquantitative analysis of the catalytic mechanism and nature of transition state in AP and NPP. Although phosphate diesters are cognate substrates for NPP but promiscuous substrates for AP, the calculations suggest that their hydrolysis reactions catalyzed by AP and NPP feature similar synchronous transition states that are slightly tighter in nature compared to that in solution, due in part to the geometry of the bimetallic zinc motif. Therefore, this study provides the first direct computational support to the hypothesis that enzymes in the AP superfamily catalyze cognate and promiscuous substrates via similar transition states to those in solution. Our calculations do not support the finding of recent QM/MM studies by López-Canut and co-workers, who suggested that the same diester substrate goes through a much looser transition state in NPP/AP than in solution, a result likely biased by the large structural distortion of the bimetallic zinc site in their simulations. Finally, our calculations for different phosphate diester orientations and phosphorothioate diesters highlight that the interpretation of thio-substitution experiments is not always straightforward.
Collapse
Affiliation(s)
- Guanhua Hou
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| | | |
Collapse
|
33
|
Bobyr E, Lassila JK, Wiersma-Koch HI, Fenn TD, Lee JJ, Nikolic-Hughes I, Hodgson KO, Rees DC, Hedman B, Herschlag D. High-resolution analysis of Zn(2+) coordination in the alkaline phosphatase superfamily by EXAFS and x-ray crystallography. J Mol Biol 2012; 415:102-17. [PMID: 22056344 PMCID: PMC3249517 DOI: 10.1016/j.jmb.2011.10.040] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2011] [Revised: 10/21/2011] [Accepted: 10/24/2011] [Indexed: 10/15/2022]
Abstract
Comparisons among evolutionarily related enzymes offer opportunities to reveal how structural differences produce different catalytic activities. Two structurally related enzymes, Escherichia coli alkaline phosphatase (AP) and Xanthomonas axonopodis nucleotide pyrophosphatase/phosphodiesterase (NPP), have nearly identical binuclear Zn(2+) catalytic centers but show tremendous differential specificity for hydrolysis of phosphate monoesters or phosphate diesters. To determine if there are differences in Zn(2+) coordination in the two enzymes that might contribute to catalytic specificity, we analyzed both x-ray absorption spectroscopic and x-ray crystallographic data. We report a 1.29-Å crystal structure of AP with bound phosphate, allowing evaluation of interactions at the AP metal site with high resolution. To make systematic comparisons between AP and NPP, we measured zinc extended x-ray absorption fine structure for AP and NPP in the free-enzyme forms, with AMP and inorganic phosphate ground-state analogs and with vanadate transition-state analogs. These studies yielded average zinc-ligand distances in AP and NPP free-enzyme forms and ground-state analog forms that were identical within error, suggesting little difference in metal ion coordination among these forms. Upon binding of vanadate to both enzymes, small increases in average metal-ligand distances were observed, consistent with an increased coordination number. Slightly longer increases were observed in NPP relative to AP, which could arise from subtle rearrangements of the active site or differences in the geometry of the bound vanadyl species. Overall, the results suggest that the binuclear Zn(2+) catalytic site remains very similar between AP and NPP during the course of a reaction cycle.
Collapse
Affiliation(s)
- Elena Bobyr
- Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | | | | | - Timothy D. Fenn
- Department of Molecular and Cellular Physiology, Stanford University, Stanford, CA 94305, USA
| | - Jason J. Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
| | - Ivana Nikolic-Hughes
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
- Department of Chemical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Keith O. Hodgson
- Department of Chemistry, Stanford University, Stanford, CA 94305, USA
- Stanford Synchrotron Radiation Lightsource, SLAC, Stanford University, Menlo Park, CA 94025, USA
| | - Douglas C. Rees
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
- Howard Hughes Medical Institute, California Institute of Technology, Pasadena, CA 91125, USA
| | - Britt Hedman
- Stanford Synchrotron Radiation Lightsource, SLAC, Stanford University, Menlo Park, CA 94025, USA
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
34
|
Meier K, Thiel W, van Gunsteren WF. On the effect of a variation of the force field, spatial boundary condition and size of the QM region in QM/MM MD simulations. J Comput Chem 2011; 33:363-78. [DOI: 10.1002/jcc.21962] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2011] [Revised: 09/10/2011] [Accepted: 09/11/2011] [Indexed: 12/12/2022]
|
35
|
López-Canut V, Roca M, Bertrán J, Moliner V, Tuñón I. Promiscuity in Alkaline Phosphatase Superfamily. Unraveling Evolution through Molecular Simulations. J Am Chem Soc 2011; 133:12050-62. [DOI: 10.1021/ja2017575] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Violeta López-Canut
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Maite Roca
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Juan Bertrán
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
36
|
Murgia S, Lampis S, Zucca P, Sanjust E, Monduzzi M. Nucleotide recognition and phosphate linkage hydrolysis at a lipid cubic interface. J Am Chem Soc 2010; 132:16176-84. [PMID: 20977215 DOI: 10.1021/ja1069745] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Mononucleotides, when entrapped within a mono-olein-based cubic Ia3d liquid crystalline phase, have been found to undergo hydrolysis at the sugar-phosphate ester bond in spite of their natural inertness toward hydrolysis. Here, kinetics of the hydrolysis reaction and interactions between the lipid matrix and the mononucleotide adenosine 5'-monophosphate disodium salt (AMP) and its 2'-deoxy derivative (dAMP) are thoroughly investigated in order to shed some light on the mechanism of the nucleotide recognition and phosphate ester hydrolysis. Experiments evidenced that molecular recognition occurs essentially through the sn-2 and the sn-3 alcoholic OH groups of mono-olein. As deduced from the apparent activation energies, the mechanism underlying the hydrolysis reaction is the same for AMP and dAMP. Nevertheless, the reaction proceeds slower for the latter, highlighting a substantial difference in the chemical behavior of the two nucleotides. A model that explains the hydrolysis reaction is presented. Remarkably, the hydrolysis mechanism appears to be highly specific for the Ia3d phase.
Collapse
Affiliation(s)
- Sergio Murgia
- Department of Chemical Science, Cagliari University, CNBS and CSGI, ss 554, bivio Sestu, 09042 Monserrato (CA), Italy.
| | | | | | | | | |
Collapse
|
37
|
López-Canut V, Roca M, Bertrán J, Moliner V, Tuñón I. Theoretical study of phosphodiester hydrolysis in nucleotide pyrophosphatase/phosphodiesterase. Environmental effects on the reaction mechanism. J Am Chem Soc 2010; 132:6955-63. [PMID: 20429564 DOI: 10.1021/ja908391v] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We here present a theoretical study of the alkaline hydrolysis of methyl p-nitrophenyl phosphate (MpNPP(-)) in aqueous solution and in the active site of nucleotide pyrophosphatase/phosphodiesterase (NPP). The analysis of our simulations, carried out by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods, shows that the reaction takes place through different reaction mechanisms depending on the environment. Thus, while in aqueous solution the reaction occurs by means of an A(N)D(N) mechanism, the enzymatic process takes place through a D(N)A(N) mechanism. In the first case, we found associative transition-state (TS) structures, while in the enzyme TS structures have dissociative character. The reason for this change is rationalized in terms of the very different nature of the electrostatic interactions established in each of the environments: while the aqueous solution reduces the repulsion between the negatively charged reacting fragments, assisting their approach, the NPP active site stabilizes the charge distribution of dissociative TS structures, allowing the reaction to proceed with a significantly reduced free energy cost. Interestingly, the NPP active site is able to accommodate different substrates, and it seems that the nature of the TSs depends on their electronic characteristics. So, in the case of the MpNPP(-) substrate, the nitro group establishes hydrogen-bond interactions with water molecules and residues found in the outer part of the catalytic site, while the leaving group oxygen atom does not coordinate directly with any of the zinc atoms of the active site. If methyl phenyl phosphate is used as substrate, then the charge on the leaving group is supported to larger extent by the oxygen atom and the phenolate anion can be then coordinated to one of the two zinc atoms present in the active site.
Collapse
Affiliation(s)
- Violeta López-Canut
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | | | | | | | | |
Collapse
|
38
|
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
|
39
|
Ferreira DEC, Florentino BPD, Rocha WR, Nome F. Quantum Mechanical/Effective Fragment Potential (QM/EFP) Study of Phosphate Monoester Aminolysis in Aqueous Solution. J Phys Chem B 2009; 113:14831-6. [DOI: 10.1021/jp907014r] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Dalva E. C. Ferreira
- LQC-MM: Laboratório de Química Computacional e Modelagem Molecular, Departamento de Química, ICEx, Universidade Federal de Minas Gerais, Campus Universitário, Pampulha, Belo Horizonte, MG 31270-901, Brazil, and LACFI: Laboratório de Catálise e Fenômenos Interfaciais Departamento de Química, Universidade Federal de Santa Catarina, Campus Trindade, Florianópolis, Santa Catarina 88040-900, Brazil
| | - Bruno P. D. Florentino
- LQC-MM: Laboratório de Química Computacional e Modelagem Molecular, Departamento de Química, ICEx, Universidade Federal de Minas Gerais, Campus Universitário, Pampulha, Belo Horizonte, MG 31270-901, Brazil, and LACFI: Laboratório de Catálise e Fenômenos Interfaciais Departamento de Química, Universidade Federal de Santa Catarina, Campus Trindade, Florianópolis, Santa Catarina 88040-900, Brazil
| | - Willian R. Rocha
- LQC-MM: Laboratório de Química Computacional e Modelagem Molecular, Departamento de Química, ICEx, Universidade Federal de Minas Gerais, Campus Universitário, Pampulha, Belo Horizonte, MG 31270-901, Brazil, and LACFI: Laboratório de Catálise e Fenômenos Interfaciais Departamento de Química, Universidade Federal de Santa Catarina, Campus Trindade, Florianópolis, Santa Catarina 88040-900, Brazil
| | - Faruk Nome
- LQC-MM: Laboratório de Química Computacional e Modelagem Molecular, Departamento de Química, ICEx, Universidade Federal de Minas Gerais, Campus Universitário, Pampulha, Belo Horizonte, MG 31270-901, Brazil, and LACFI: Laboratório de Catálise e Fenômenos Interfaciais Departamento de Química, Universidade Federal de Santa Catarina, Campus Trindade, Florianópolis, Santa Catarina 88040-900, Brazil
| |
Collapse
|