1
|
Song Z, Trozzi F, Tian H, Yin C, Tao P. Mechanistic Insights into Enzyme Catalysis from Explaining Machine-Learned Quantum Mechanical and Molecular Mechanical Minimum Energy Pathways. ACS PHYSICAL CHEMISTRY AU 2022; 2:316-330. [PMID: 35936506 PMCID: PMC9344433 DOI: 10.1021/acsphyschemau.2c00005] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
With the increasing popularity of machine learning (ML) applications, the demand for explainable artificial intelligence techniques to explain ML models developed for computational chemistry has also emerged. In this study, we present the development of the Boltzmann-weighted cumulative integrated gradients (BCIG) approach for effective explanation of mechanistic insights into ML models trained on high-level quantum mechanical and molecular mechanical (QM/MM) minimum energy pathways. Using the acylation reactions of the Toho-1 β-lactamase and two antibiotics (ampicillin and cefalexin) as the model systems, we show that the BCIG approach could quantitatively attribute the energetic contribution in one system and the relative reactivity of individual steps across different systems to specific chemical processes such as the bond making/breaking and proton transfers. The proposed BCIG contribution attribution method quantifies chemistry-interpretable insights in terms of contributions from each elementary chemical process, which is in agreement with the validating QM/MM calculations and our intuitive mechanistic understandings of the model reactions.
Collapse
|
2
|
Gervasoni S, Spencer J, Hinchliffe P, Pedretti A, Vairoletti F, Mahler G, Mulholland AJ. A multiscale approach to predict the binding mode of metallo beta-lactamase inhibitors. Proteins 2022; 90:372-384. [PMID: 34455628 PMCID: PMC8944931 DOI: 10.1002/prot.26227] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2021] [Revised: 06/09/2021] [Accepted: 08/18/2021] [Indexed: 02/03/2023]
Abstract
Antibiotic resistance is a major threat to global public health. β-lactamases, which catalyze breakdown of β-lactam antibiotics, are a principal cause. Metallo β-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all β-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems.
Collapse
Affiliation(s)
- Silvia Gervasoni
- Department of Pharmaceutical Sciences, University of Milan, Milan, Italy
| | - James Spencer
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, UK
| | - Philip Hinchliffe
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, UK
| | | | - Franco Vairoletti
- Laboratorio de Química Farmacéutica, Departamento de Química Orgánica, Facultad de Química, Universidad de la República (UdelaR), Avda. General Flores 2124, Montevideo, Uruguay
| | - Graciela Mahler
- Laboratorio de Química Farmacéutica, Departamento de Química Orgánica, Facultad de Química, Universidad de la República (UdelaR), Avda. General Flores 2124, Montevideo, Uruguay
| | | |
Collapse
|
3
|
Galdadas I, Qu S, Oliveira ASF, Olehnovics E, Mack AR, Mojica MF, Agarwal PK, Tooke CL, Gervasio FL, Spencer J, Bonomo RA, Mulholland AJ, Haider S. Allosteric communication in class A β-lactamases occurs via cooperative coupling of loop dynamics. eLife 2021; 10:e66567. [PMID: 33755013 PMCID: PMC8060031 DOI: 10.7554/elife.66567] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 03/19/2021] [Indexed: 12/16/2022] Open
Abstract
Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A β-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.
Collapse
Affiliation(s)
- Ioannis Galdadas
- University College London, Department of ChemistryLondonUnited Kingdom
| | - Shen Qu
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| | - Ana Sofia F Oliveira
- University of Bristol, Centre for Computational Chemistry, School of ChemistryBristolUnited Kingdom
| | - Edgar Olehnovics
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| | - Andrew R Mack
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Molecular Biology and MicrobiologyClevelandUnited States
| | - Maria F Mojica
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Infectious Diseases, School of MedicineClevelandUnited States
| | - Pratul K Agarwal
- Department of Physiological Sciences and High-Performance Computing Center, Oklahoma State UniversityStillwaterUnited States
| | - Catherine L Tooke
- University of Bristol, School of Cellular and Molecular MedicineBristolUnited Kingdom
| | - Francesco Luigi Gervasio
- University College London, Department of ChemistryLondonUnited Kingdom
- University College London, Institute of Structural and Molecular BiologyLondonUnited Kingdom
- University of Geneva, Pharmaceutical SciencesGenevaSwitzerland
| | - James Spencer
- University of Bristol, School of Cellular and Molecular MedicineBristolUnited Kingdom
| | - Robert A Bonomo
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Molecular Biology and MicrobiologyClevelandUnited States
- Case Western Reserve University, Department of Infectious Diseases, School of MedicineClevelandUnited States
- Case Western Reserve University, Department of BiochemistryClevelandUnited States
- Case Western Reserve University, Department of PharmacologyClevelandUnited States
- Case Western Reserve University, Department of Proteomics and BioinformaticsClevelandUnited States
- CWRU-Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES)ClevelandUnited States
| | - Adrian J Mulholland
- University of Bristol, Centre for Computational Chemistry, School of ChemistryBristolUnited Kingdom
| | - Shozeb Haider
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| |
Collapse
|
4
|
Oliveira ASF, Ciccotti G, Haider S, Mulholland AJ. Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:144. [PMID: 34720710 PMCID: PMC8549953 DOI: 10.1140/epjb/s10051-021-00157-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/05/2021] [Indexed: 05/03/2023]
Abstract
ABSTRACT A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ UK
| | - Giovanni Ciccotti
- Institute for Applied Computing “Mauro Picone” (IAC), CNR, Via dei Taurini 19, 00185 Rome, Italy
- School of Physics, University College of Dublin, UCD-Belfield, Dublin 4, Ireland
- Università di Roma La Sapienza, Ple. A. Moro 5, 00185 Rome, Italy
| | - Shozeb Haider
- School of Pharmacy, University College London, London, WC1N 1AX UK
| | - Adrian J. Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
5
|
Prunotto A, Bahr G, González LJ, Vila AJ, Dal Peraro M. Molecular Bases of the Membrane Association Mechanism Potentiating Antibiotic Resistance by New Delhi Metallo-β-lactamase 1. ACS Infect Dis 2020; 6:2719-2731. [PMID: 32865963 DOI: 10.1021/acsinfecdis.0c00341] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Resistance to last-resort carbapenem antibiotics is an increasing threat to human health, as it critically limits therapeutic options. Metallo-β-lactamases (MBLs) are the largest family of carbapenemases, enzymes that inactivate these drugs. Among MBLs, New Delhi metallo-β-lactamase 1 (NDM-1) has experienced the fastest and largest worldwide dissemination. This success has been attributed to the fact that NDM-1 is a lipidated protein anchored to the outer membrane of bacteria, while all other MBLs are soluble periplasmic enzymes. By means of a combined experimental and computational approach, we show that NDM-1 interacts with the surface of bacterial membranes in a stable, defined conformation, in which the active site is not occluded by the bilayer. Although the lipidation is required for a long-lasting interaction, the globular domain of NDM-1 is tuned to interact specifically with the outer bacterial membrane. In contrast, this affinity is not observed for VIM-2, a natively soluble MBL. Finally, we identify key residues involved in the membrane interaction with NDM-1, which constitute potential targets for developing therapeutic strategies able to combat resistance granted by this enzyme.
Collapse
Affiliation(s)
- Alessio Prunotto
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Guillermo Bahr
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Lisandro J. González
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Alejandro J. Vila
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Matteo Dal Peraro
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
6
|
Marcos-Alcalde I, López-Viñas E, Gómez-Puertas P. MEPSAnd: minimum energy path surface analysis over n-dimensional surfaces. Bioinformatics 2020; 36:956-958. [PMID: 31418769 DOI: 10.1093/bioinformatics/btz649] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2019] [Revised: 07/30/2019] [Accepted: 08/14/2019] [Indexed: 01/01/2023] Open
Abstract
SUMMARY n-dimensional energy surfaces are becoming computationally accessible, yet interpreting their information is not straightforward. We present minimum energy path surface analysis over n-dimensional surfaces (MEPSAnd), an open source GUI-based program that natively calculates minimum energy paths across energy surfaces of any number of dimensions. Among other features, MEPSAnd can compute the path through lowest barriers and automatically provide a set of alternative paths. MEPSAnd offers distinct plotting solutions as well as direct python scripting. AVAILABILITY AND IMPLEMENTATION MEPSAnd is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSAnd/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Iñigo Marcos-Alcalde
- Molecular Modelling Group, Centro de Biología Molecular Severo Ochoa, CBMSO (Consejo Superior de Investigaciones Científicas-Universidad Autónoma de Madrid, CSIC-UAM), Cabrera, CL Nicolás 1, E-28049 Madrid, Spain.,Biosciences Research Institute, School of Experimental Sciences, Universidad Francisco de Vitoria, UFV, Building E, Ctra. M-515 Pozuelo-Majadahonda Km 1,800, Pozuelo de Alarcón, Madrid 28223, Spain
| | - Eduardo López-Viñas
- Biosciences Research Institute, School of Experimental Sciences, Universidad Francisco de Vitoria, UFV, Building E, Ctra. M-515 Pozuelo-Majadahonda Km 1,800, Pozuelo de Alarcón, Madrid 28223, Spain
| | - Paulino Gómez-Puertas
- Molecular Modelling Group, Centro de Biología Molecular Severo Ochoa, CBMSO (Consejo Superior de Investigaciones Científicas-Universidad Autónoma de Madrid, CSIC-UAM), Cabrera, CL Nicolás 1, E-28049 Madrid, Spain
| |
Collapse
|
7
|
Hirvonen VHA, Mulholland AJ, Spencer J, van der Kamp MW. Small Changes in Hydration Determine Cephalosporinase Activity of OXA-48 β-Lactamases. ACS Catal 2020. [DOI: 10.1021/acscatal.0c00596] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Viivi H. A. Hirvonen
- School of Biochemistry, University of Bristol, University Walk, Bristol, BS8 1TD United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol, BS8 1TS United Kingdom
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol, BS8 1TS United Kingdom
| | - James Spencer
- School of Cellular and Molecular Medicine, University of Bristol, University Walk, Bristol, BS8 1TD United Kingdom
| | - Marc W. van der Kamp
- School of Biochemistry, University of Bristol, University Walk, Bristol, BS8 1TD United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol, BS8 1TS United Kingdom
| |
Collapse
|
8
|
Tooke CL, Hinchliffe P, Lang PA, Mulholland AJ, Brem J, Schofield CJ, Spencer J. Molecular Basis of Class A β-Lactamase Inhibition by Relebactam. Antimicrob Agents Chemother 2019; 63:AAC.00564-19. [PMID: 31383664 PMCID: PMC6761529 DOI: 10.1128/aac.00564-19] [Citation(s) in RCA: 51] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2019] [Accepted: 07/04/2019] [Indexed: 12/23/2022] Open
Abstract
β-Lactamase production is the major β-lactam resistance mechanism in Gram-negative bacteria. β-Lactamase inhibitors (BLIs) efficacious against serine β-lactamase (SBL) producers, especially strains carrying the widely disseminated class A enzymes, are required. Relebactam, a diazabicyclooctane (DBO) BLI, is in phase 3 clinical trials in combination with imipenem for the treatment of infections by multidrug-resistant Enterobacteriaceae We show that relebactam inhibits five clinically important class A SBLs (despite their differing spectra of activity), representing both chromosomal and plasmid-borne enzymes, i.e., the extended-spectrum β-lactamases L2 (inhibition constant 3 μM) and CTX-M-15 (21 μM) and the carbapenemases KPC-2, -3, and -4 (1 to 5 μM). Against purified class A SBLs, relebactam is an inferior inhibitor compared with the clinically approved DBO avibactam (9- to 120-fold differences in half maximal inhibitory concentration [IC50]). MIC assays indicate relebactam potentiates β-lactam (imipenem) activity against KPC-producing Klebsiella pneumoniae, with similar potency to avibactam (with ceftazidime). Relebactam is less effective than avibactam in combination with aztreonam against Stenotrophomonas maltophilia K279a. X-ray crystal structures of relebactam bound to CTX-M-15, L2, KPC-2, KPC-3, and KPC-4 reveal its C2-linked piperidine ring can sterically clash with Asn104 (CTX-M-15) or His/Trp105 (L2 and KPCs), rationalizing its poorer inhibition activity than that of avibactam, which has a smaller C2 carboxyamide group. Mass spectrometry and crystallographic data show slow, pH-dependent relebactam desulfation by KPC-2, -3, and -4. This comprehensive comparison of relebactam binding across five clinically important class A SBLs will inform the design of future DBOs, with the aim of improving clinical efficacy of BLI-β-lactam combinations.
Collapse
Affiliation(s)
- Catherine L Tooke
- School of Cellular and Molecular Medicine, Biomedical Sciences Building, University of Bristol, Bristol, United Kingdom
| | - Philip Hinchliffe
- School of Cellular and Molecular Medicine, Biomedical Sciences Building, University of Bristol, Bristol, United Kingdom
| | - Pauline A Lang
- Department of Chemistry, University of Oxford, Oxford, United Kingdom
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | - Jürgen Brem
- Department of Chemistry, University of Oxford, Oxford, United Kingdom
| | | | - James Spencer
- School of Cellular and Molecular Medicine, Biomedical Sciences Building, University of Bristol, Bristol, United Kingdom
| |
Collapse
|
9
|
Hirvonen VHA, Hammond K, Chudyk EI, Limb MAL, Spencer J, Mulholland AJ, van der Kamp MW. An Efficient Computational Assay for β-Lactam Antibiotic Breakdown by Class A β-Lactamases. J Chem Inf Model 2019; 59:3365-3369. [PMID: 31361944 DOI: 10.1021/acs.jcim.9b00442] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Class A β-lactamases cause clinically relevant resistance to β-lactam antibiotics. Carbapenem degradation is a particular concern. We present an efficient QM/MM molecular simulation protocol that accurately predicts the activity of β-lactamases against carbapenems. Simulations take less than 24 CPU hours, a greater than 99% reduction, and do not require fitting against experimental data or significant parametrization. This computational assay also reveals mechanistic details of β-lactam breakdown and should assist in evaluating emerging β-lactamase variants and developing new antibiotics.
Collapse
Affiliation(s)
- Viivi H A Hirvonen
- School of Biochemistry , University of Bristol , University Walk, Bristol BS8 1TD , United Kingdom
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| | - Katharine Hammond
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| | - Ewa I Chudyk
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| | - Michael A L Limb
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| | - James Spencer
- School of Cellular and Molecular Medicine , University of Bristol , University Walk, Bristol BS8 1TD , United Kingdom
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| | - Marc W van der Kamp
- School of Biochemistry , University of Bristol , University Walk, Bristol BS8 1TD , United Kingdom
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Cantock's Close, Bristol BS8 1TS , United Kingdom
| |
Collapse
|
10
|
Nutho B, Mulholland AJ, Rungrotmongkol T. Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate. J Phys Chem B 2019; 123:2889-2903. [PMID: 30845796 DOI: 10.1021/acs.jpcb.9b02157] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.
Collapse
Affiliation(s)
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol BS8 1TS , U.K
| | | |
Collapse
|
11
|
Cortina GA, Kasson PM. Predicting allostery and microbial drug resistance with molecular simulations. Curr Opin Struct Biol 2018; 52:80-86. [PMID: 30243041 PMCID: PMC6296865 DOI: 10.1016/j.sbi.2018.09.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2018] [Revised: 08/31/2018] [Accepted: 09/07/2018] [Indexed: 11/30/2022]
Abstract
Beta-lactamase enzymes mediate the most common forms of gram-negative antibiotic resistance affecting clinical treatment. They also constitute an excellent model system for the difficult problem of understanding how allosteric mutations can augment catalytic activity of already-competent enzymes. Multiple allosteric mutations have been identified that alter catalytic activity or drug-resistance spectrum in class A beta lactamases, but predicting these in advance continues to be challenging. Here, we review computational techniques based on structure and/or molecular simulation to predict such mutations. Structure-based techniques have been particularly helpful in developing graph algorithms for analyzing critical residues in beta-lactamase function, while classical molecular simulation has recently shown the ability to prospectively predict allosteric mutations increasing beta-lactamase activity and drug resistance. These will ultimately achieve the greatest power when combined with simulation methods that model reactive chemistry to calculate activation free energies directly.
Collapse
Affiliation(s)
- George A Cortina
- Departments of Molecular Physiology and of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, United States
| | - Peter M Kasson
- Departments of Molecular Physiology and of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, United States; Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala 75146, Sweden.
| |
Collapse
|
12
|
Huggins DJ, Biggin PC, Dämgen MA, Essex JW, Harris SA, Henchman RH, Khalid S, Kuzmanic A, Laughton CA, Michel J, Mulholland AJ, Rosta E, Sansom MSP, van der Kamp MW. Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1393] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- David J. Huggins
- TCM Group, Cavendish Laboratory University of Cambridge Cambridge UK
- Unilever Centre, Department of Chemistry University of Cambridge Cambridge UK
- Department of Physiology and Biophysics Weill Cornell Medical College New York NY
| | | | - Marc A. Dämgen
- Department of Biochemistry University of Oxford Oxford UK
| | - Jonathan W. Essex
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | - Sarah A. Harris
- School of Physics and Astronomy University of Leeds Leeds UK
- Astbury Centre for Structural and Molecular Biology University of Leeds Leeds UK
| | - Richard H. Henchman
- Manchester Institute of Biotechnology The University of Manchester Manchester UK
- School of Chemistry The University of Manchester Oxford UK
| | - Syma Khalid
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | | | - Charles A. Laughton
- School of Pharmacy University of Nottingham Nottingham UK
- Centre for Biomolecular Sciences University of Nottingham Nottingham UK
| | - Julien Michel
- EaStCHEM school of Chemistry University of Edinburgh Edinburgh UK
| | - Adrian J. Mulholland
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
| | - Edina Rosta
- Department of Chemistry King's College London London UK
| | | | - Marc W. van der Kamp
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
- School of Biochemistry, Biomedical Sciences Building University of Bristol Bristol UK
| |
Collapse
|