1
|
Kirilin EM, Švedas VK. Analysis of Glycosyl-Enzyme Intermediate Formation in the Catalytic Mechanism of Influenza Virus Neuraminidase Using Molecular Modeling. BIOCHEMISTRY. BIOKHIMIIA 2020; 85:490-498. [PMID: 32569556 DOI: 10.1134/s0006297920040094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2020] [Revised: 03/02/2020] [Accepted: 03/02/2020] [Indexed: 06/11/2023]
Abstract
Using classical molecular dynamics, constant-pH molecular dynamics simulation, metadynamics, and combined quantum mechanical and molecular mechanical approach, we identified an alternative pathway of glycosyl-enzyme intermediate formation during oligosaccharide substrate conversion by the influenza H5N1 neuraminidase. The Asp151 residue located in the enzyme mobile loop plays a key role in catalysis within a wide pH range due to the formation of a network of interactions with water molecules. Considering that propagation of influenza virus takes place in the digestive tract of birds at low pH values and in the human respiratory tract at pH values close to neutral, the existence of alternative reaction pathways functioning at different medium pH can explain the dual tropism of the virus and circulation of H5N1 viral strains capable of transmission from birds to humans.
Collapse
Affiliation(s)
- E M Kirilin
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, 119991, Russia.
| | - V K Švedas
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, 119991, Russia.
- Lomonosov Moscow State University, Faculty of Bioengineering and Bioinformatics, Moscow, 119991, Russia
| |
Collapse
|
2
|
Coines J, Raich L, Rovira C. Modeling catalytic reaction mechanisms in glycoside hydrolases. Curr Opin Chem Biol 2019; 53:183-191. [PMID: 31731209 DOI: 10.1016/j.cbpa.2019.09.007] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Revised: 09/22/2019] [Accepted: 09/24/2019] [Indexed: 01/21/2023]
Abstract
Modeling catalysis in carbohydrate-active enzymes is a daunting challenge because of the high flexibility and diversity of both enzymes and carbohydrates. Glycoside hydrolases (GHs) are an illustrative example, where conformational changes and subtle interactions have been shown to be critical for catalysis. GHs have pivotal roles in industry (e.g. biofuel or detergent production) and biomedicine (e.g. targets for cancer and diabetes), and thus, a huge effort is devoted to unveil their molecular mechanisms. Besides experimental techniques, computational methods have served to provide an in-depth understanding of GH mechanisms, capturing complex reaction coordinates and the conformational itineraries that substrates follow during the whole catalytic pathway, providing a framework that ultimately may assist the engineering of these enzymes and the design of new inhibitors.
Collapse
Affiliation(s)
- Joan Coines
- Departament de Química Inorgànica I Orgànica (secció de Química Orgànica) and Institut de Química Teòrica I Computacional (IQTCUB), Universitat de Barcelona, Martí I Franquès 1 08028 Barcelona, Spain
| | - Lluís Raich
- Departament de Química Inorgànica I Orgànica (secció de Química Orgànica) and Institut de Química Teòrica I Computacional (IQTCUB), Universitat de Barcelona, Martí I Franquès 1 08028 Barcelona, Spain
| | - Carme Rovira
- Departament de Química Inorgànica I Orgànica (secció de Química Orgànica) and Institut de Química Teòrica I Computacional (IQTCUB), Universitat de Barcelona, Martí I Franquès 1 08028 Barcelona, Spain; Institució Catalana de Recerca I Estudis Avançats (ICREA) Passeig Lluís Companys 23, 08010 Barcelona, Spain.
| |
Collapse
|
3
|
Xiao K, Wang X, Yu H. Comparative studies of catalytic pathways for Streptococcus pneumoniae sialidases NanA, NanB and NanC. Sci Rep 2019; 9:2157. [PMID: 30770840 PMCID: PMC6377674 DOI: 10.1038/s41598-018-38131-z] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Accepted: 12/14/2018] [Indexed: 11/09/2022] Open
Abstract
Streptococcus pneumoniae (S. pneumoniae) is a leading human pathogen, which takes large responsibility for severe otitis media, acute meningitis and septicaemia. It encodes up to three distinct sialidases: NanA, NanB and NanC, which are promising drug targets. Recent experimental studies have shown that these three sialidases might work together up to the ultimate step, where NanA and NanB produce N-acetylneuraminic acid (Neu5Ac) and 2,7-anhydro-Neu5Ac following the functions of sialidase and intramolecular trans-sialidase, whilst NanC carries on a ping-pong mechanism that produces or removes 2-deoxy-2,3-didehydro-Neu5AC. It is intriguing that these sialidases have similar active sites but operate via three distinct reaction pathways. To clarify this issue, herein we present the first systematic computational investigation on the catalytic pathways for S. pneumoniae NanA, NanB and NanC based on combined quantum mechanics/molecular mechanics simulations, and propose the most preferred routes for the three S. pneumoniae sialidases. Our findings support the mechanisms of NanA and NanC that were proposed by previous experimental studies, whereas the role of water in NanB was found to differ slightly from our current understandings. The mechanistic insights obtained from this work are expected to assist in the design of potent inhibitors targeting these key enzymes for therapeutic applications.
Collapse
Affiliation(s)
- Kela Xiao
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong, NSW, 2500, Australia.,Molecular Horizons, University of Wollongong, Wollongong, NSW, 2500, Australia
| | - Xingyong Wang
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong, NSW, 2500, Australia.,Molecular Horizons, University of Wollongong, Wollongong, NSW, 2500, Australia
| | - Haibo Yu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong, NSW, 2500, Australia. .,Molecular Horizons, University of Wollongong, Wollongong, NSW, 2500, Australia. .,Illawarra Health and Medical Research Institute, Wollongong, NSW, 2500, Australia.
| |
Collapse
|
4
|
Shi N, Zheng QC, Zhang HX. QM/MM calculations and MD simulations on acetate CoA-transferase to reveal its catalytic mechanism and illuminate the role of residue Asn347. MOLECULAR CATALYSIS 2018. [DOI: 10.1016/j.mcat.2018.01.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
5
|
Rogers IL, Naidoo KJ. Producing DFT/MM enzyme reaction trajectories from SCC-DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction. J Comput Chem 2017; 38:1789-1798. [PMID: 28488320 DOI: 10.1002/jcc.24820] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2016] [Revised: 03/14/2017] [Accepted: 04/03/2017] [Indexed: 02/05/2023]
Abstract
The SCC-DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC-DFTB/MIO/CHARMM reaction surface to drive B3LYP/6-31G/MM and B3LYP/6-31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6-31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6-31G/MM-computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi-empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ian L Rogers
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch, 7701, South Africa
| | - Kevin J Naidoo
- Scientific Computing Research Unit and Department of Chemistry, University of Cape Town, Rondebosch, 7701, South Africa
| |
Collapse
|
6
|
Sánchez-Murcia PA, Bueren-Calabuig JA, Camacho-Artacho M, Cortés-Cabrera Á, Gago F. Stepwise Simulation of 3,5-Dihydro-5-methylidene-4H-imidazol-4-one (MIO) Biogenesis in Histidine Ammonia-lyase. Biochemistry 2016; 55:5854-5864. [DOI: 10.1021/acs.biochem.6b00744] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pedro A. Sánchez-Murcia
- Área
de Farmacología, Departamento de Ciencias Biomédicas,
Unidad Asociada al IQM-CSIC, Universidad de Alcalá, E-28805 Alcalá de Henares, Spain
| | - Juan A. Bueren-Calabuig
- Área
de Farmacología, Departamento de Ciencias Biomédicas,
Unidad Asociada al IQM-CSIC, Universidad de Alcalá, E-28805 Alcalá de Henares, Spain
| | - Marta Camacho-Artacho
- Structural
Biology Department, Centro Nacional de Investigaciones Oncológicas (CNIO), E-28029 Madrid, Spain
| | - Álvaro Cortés-Cabrera
- Área
de Farmacología, Departamento de Ciencias Biomédicas,
Unidad Asociada al IQM-CSIC, Universidad de Alcalá, E-28805 Alcalá de Henares, Spain
| | - Federico Gago
- Área
de Farmacología, Departamento de Ciencias Biomédicas,
Unidad Asociada al IQM-CSIC, Universidad de Alcalá, E-28805 Alcalá de Henares, Spain
| |
Collapse
|
7
|
In silico structural characterization of protein targets for drug development against Trypanosoma cruzi. J Mol Model 2016; 22:244. [DOI: 10.1007/s00894-016-3115-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2016] [Accepted: 09/02/2016] [Indexed: 10/21/2022]
|
8
|
Rogers IL, Naidoo KJ. Multidimensional Reaction Dynamics Reveal How the Enzyme TcTS Suppresses Competing Side Reactions and Their Side Products. ACS Catal 2016. [DOI: 10.1021/acscatal.6b01522] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ian L. Rogers
- Scientific Computing Research Unit and ‡Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| | - Kevin J. Naidoo
- Scientific Computing Research Unit and ‡Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| |
Collapse
|
9
|
Raich L, Borodkin V, Fang W, Castro-López J, van Aalten DMF, Hurtado-Guerrero R, Rovira C. A Trapped Covalent Intermediate of a Glycoside Hydrolase on the Pathway to Transglycosylation. Insights from Experiments and Quantum Mechanics/Molecular Mechanics Simulations. J Am Chem Soc 2016; 138:3325-32. [DOI: 10.1021/jacs.5b10092] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Lluís Raich
- Departament
de Química Inorgànica i Orgànica and Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
| | | | | | - Jorge Castro-López
- Institute
of Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, BIFI-IQFR (CSIC) Joint Unit, Mariano Esquillor s/n, Campus Rio
Ebro, Edificio I+D, 50018 Zaragoza, Spain
| | | | - Ramón Hurtado-Guerrero
- Fundación ARAID, Edificio CEEI
Aragón, 50018 Zaragoza, Spain
- Institute
of Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, BIFI-IQFR (CSIC) Joint Unit, Mariano Esquillor s/n, Campus Rio
Ebro, Edificio I+D, 50018 Zaragoza, Spain
| | - Carme Rovira
- Departament
de Química Inorgànica i Orgànica and Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys, 23, 08020 Barcelona, Spain
| |
Collapse
|
10
|
Silva JRA, Roitberg AE, Alves CN. A QM/MM free energy study of the oxidation mechanism of dihydroorotate dehydrogenase (class 1A) from Lactococcus lactis. J Phys Chem B 2015; 119:1468-73. [PMID: 25564307 DOI: 10.1021/jp512860r] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The dihydroorotate dehydrogenase (DHOD) class 1A enzyme catalyzes is the only redox enzyme in the biosynthetic pathway (de novo) of pyrimidines where dihydroorotate (DHO) is oxidized to orotate (OA) coupled to reduction of a flavin mononucleotide (FMN) cofactor. The rupture of two DHO C-H bonds can proceed in a concerted or stepwise way. Herein, the catalytic mechanism of DHOD from Lactococcus lactis involving DHO oxidation (first half-reaction) was described using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach and molecular dynamics simulations. The free energy profile obtained from self-consistent charge-density functional tight binding/molecular mechanics calculations (corrected by DFT/MM) reveals that this occurs with the proton abstraction from DHO C5 to Cys130 deprotonated and DHO H6 is transferred to FMN N5 in a concerted mechanism with a very low barrier of 5.64 kcal/mol. Finally, through a residual decomposition analysis, the residues that have the main influence on the redox reaction were identified.
Collapse
Affiliation(s)
- José Rogério A Silva
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará , Belém, PA 66075-110, Brazil
| | | | | |
Collapse
|
11
|
Rogers IL, Naidoo KJ. Profiling Transition-State Configurations on the Trypanosoma cruzi trans-Sialidase Free-Energy Reaction Surfaces. J Phys Chem B 2014; 119:1192-201. [DOI: 10.1021/jp506824r] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ian L. Rogers
- Scientific
Computing Research Unit and ‡Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| | - Kevin J. Naidoo
- Scientific
Computing Research Unit and ‡Department of Chemistry, University of Cape Town, Rondebosch 7701, South Africa
| |
Collapse
|