76
|
Kamsri P, Punkvang A, Hannongbua S, Suttisintong K, Kittakoop P, Spencer J, Mulholland AJ, Pungpo P. In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR, CoMSIA and molecular dynamics simulations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2019; 30:775-800. [PMID: 31607177 DOI: 10.1080/1062936x.2019.1658218] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 08/17/2019] [Indexed: 06/10/2023]
Abstract
Mycobacterium tuberculosis DNA gyrase subunit B (GyrB) has been identified as a promising target for rational drug design against fluoroquinolone drug-resistant tuberculosis. In this study, we attempted to identify the key structural feature for highly potent GyrB inhibitors through 2D-QSAR using HQSAR, 3D-QSAR using CoMSIA and molecular dynamics (MD) simulations approaches on a series of thiazole urea core derivatives. The best HQSAR and CoMSIA models based on IC50 and MIC displayed the structural basis required for good activity against both GyrB enzyme and mycobacterial cell. MD simulations and binding free energy analysis using MM-GBSA and waterswap calculations revealed that the urea core of inhibitors has the strongest interaction with Asp79 via hydrogen bond interactions. In addition, cation-pi interaction and hydrophobic interactions of the R2 substituent with Arg82 and Arg141 help to enhance the binding affinity in the GyrB ATPase binding site. Thus, the present study provides crucial structural features and a structural concept for rational design of novel DNA gyrase inhibitors with improved biological activities against both enzyme and mycobacterial cell, and with good pharmacokinetic properties and drug safety profiles.
Collapse
|
77
|
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: 41] [Impact Index Per Article: 8.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
|
78
|
Kiani YS, Ranaghan KE, Jabeen I, Mulholland AJ. Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors. Int J Mol Sci 2019; 20:ijms20184468. [PMID: 31510073 PMCID: PMC6769491 DOI: 10.3390/ijms20184468] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 08/22/2019] [Accepted: 09/01/2019] [Indexed: 12/20/2022] Open
Abstract
The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.
Collapse
|
79
|
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
|
80
|
Bunzel HA, Kries H, Marchetti L, Zeymer C, Mittl PRE, Mulholland AJ, Hilvert D. Emergence of a Negative Activation Heat Capacity during Evolution of a Designed Enzyme. J Am Chem Soc 2019; 141:11745-11748. [PMID: 31282667 DOI: 10.1021/jacs.9b02731] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Temperature influences the reaction kinetics and evolvability of all enzymes. To understand how evolution shapes the thermodynamic drivers of catalysis, we optimized the modest activity of a computationally designed enzyme for an elementary proton-transfer reaction by nearly 4 orders of magnitude over 9 rounds of mutagenesis and screening. As theorized for primordial enzymes, the catalytic effects of the original design were almost entirely enthalpic in origin, as were the rate enhancements achieved by laboratory evolution. However, the large reductions in ΔH⧧ were partially offset by a decrease in TΔS⧧ and unexpectedly accompanied by a negative activation heat capacity, signaling strong adaptation to the operating temperature. These findings echo reports of temperature-dependent activation parameters for highly evolved natural enzymes and are relevant to explanations of enzymatic catalysis and adaptation to changing thermal environments.
Collapse
|
81
|
Nutho B, Mulholland AJ, Rungrotmongkol T. The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study. Phys Chem Chem Phys 2019; 21:14945-14956. [PMID: 31236554 DOI: 10.1039/c9cp02377a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.
Collapse
|
82
|
Oliveira ASF, Shoemark DK, Campello HR, Wonnacott S, Gallagher T, Sessions RB, Mulholland AJ. Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations. Structure 2019; 27:1171-1183.e3. [PMID: 31130483 DOI: 10.1016/j.str.2019.04.008] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 01/28/2019] [Accepted: 04/10/2019] [Indexed: 02/02/2023]
Abstract
Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.
Collapse
|
83
|
O'Connor MB, Bennie SJ, Deeks HM, Jamieson-Binnie A, Jones AJ, Shannon RJ, Walters R, Mitchell TJ, Mulholland AJ, Glowacki DR. Interactive molecular dynamics in virtual reality from quantum chemistry to drug binding: An open-source multi-person framework. J Chem Phys 2019; 150:220901. [PMID: 31202243 DOI: 10.1063/1.5092590] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
As molecular scientists have made progress in their ability to engineer nanoscale molecular structure, we face new challenges in our ability to engineer molecular dynamics (MD) and flexibility. Dynamics at the molecular scale differs from the familiar mechanics of everyday objects because it involves a complicated, highly correlated, and three-dimensional many-body dynamical choreography which is often nonintuitive even for highly trained researchers. We recently described how interactive molecular dynamics in virtual reality (iMD-VR) can help to meet this challenge, enabling researchers to manipulate real-time MD simulations of flexible structures in 3D. In this article, we outline various efforts to extend immersive technologies to the molecular sciences, and we introduce "Narupa," a flexible, open-source, multiperson iMD-VR software framework which enables groups of researchers to simultaneously cohabit real-time simulation environments to interactively visualize and manipulate the dynamics of molecular structures with atomic-level precision. We outline several application domains where iMD-VR is facilitating research, communication, and creative approaches within the molecular sciences, including training machines to learn potential energy functions, biomolecular conformational sampling, protein-ligand binding, reaction discovery using "on-the-fly" quantum chemistry, and transport dynamics in materials. We touch on iMD-VR's various cognitive and perceptual affordances and outline how these provide research insight for molecular systems. By synergistically combining human spatial reasoning and design insight with computational automation, technologies such as iMD-VR have the potential to improve our ability to understand, engineer, and communicate microscopic dynamical behavior, offering the potential to usher in a new paradigm for engineering molecules and nano-architectures.
Collapse
|
84
|
Ranaghan KE, Shchepanovska D, Bennie SJ, Lawan N, Macrae SJ, Zurek J, Manby FR, Mulholland AJ. Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution. J Chem Inf Model 2019; 59:2063-2078. [PMID: 30794388 DOI: 10.1021/acs.jcim.8b00940] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.
Collapse
|
85
|
Leferink NGH, Ranaghan KE, Karuppiah V, Currin A, van der Kamp MW, Mulholland AJ, Scrutton NS. Experiment and Simulation Reveal How Mutations in Functional Plasticity Regions Guide Plant Monoterpene Synthase Product Outcome. ACS Catal 2019; 8:3780-3791. [PMID: 31157124 PMCID: PMC6542672 DOI: 10.1021/acscatal.8b00692] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Monoterpenes (C10 isoprenoids) are a structurally diverse group of natural compounds that are attractive to industry as flavours and fragrances. Monoterpenes are produced from a single linear substrate, geranyl diphosphate, by a group of enzymes called the monoterpene cyclases/synthases (mTC/Ss) that catalyse high-energy cyclisation reactions involving unstable carbocation intermediates. Efforts towards producing monoterpenes via biocatalysis or metabolic engineering often result in the formation of multiple products due to the nature of the highly branched reaction mechanism of mTC/Ss. Rational engineering of mTC/Ss is hampered by the lack of correlation between the active site sequence and cyclisation type. We used available mutagenesis data to show that amino acids involved in product outcome are clustered and spatially conserved within the mTC/S family. Consensus sequences for three such plasticity regions were introduced in different mTC/S with increasingly complex cyclisation cascades, including the model enzyme limonene synthase (LimS). In all three mTC/S studied, mutations in the first two regions mostly give rise to products that result from premature quenching of the linalyl or α-terpinyl cations, suggesting that both plasticity regions are involved in the formation and stabilisation of cations early in the reaction cascade. A LimS variant with mutations in the second region (S454G, C457V, M458I), produced mainly more complex bicyclic products. QM/MM MD simulations reveal that the second cyclisation is not due to compression of the C2-C7 distance in the α-terpinyl cation, but is the result of an increased distance between C8 of the α-terpinyl cation and two putative bases (W324, H579) located on the other side of the active site, preventing early termination by deprotonation. Such insights into the impact of mutations can only be obtained using integrated experimental and computational approaches, and will aid the design of altered mTC/S activities towards clean monoterpenoid products.
Collapse
|
86
|
Raza S, Ranaghan KE, van der Kamp MW, Woods CJ, Mulholland AJ, Azam SS. Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site. J Comput Aided Mol Des 2019; 33:461-475. [PMID: 30989572 DOI: 10.1007/s10822-019-00200-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Accepted: 04/05/2019] [Indexed: 02/08/2023]
Abstract
Kallikrein-8, a serine protease, is a target for structure-based drug design due to its therapeutic potential in treating Alzheimer's disease and is also useful as a biomarker in ovarian cancer. We present a binding assessment of ligands to kallikrein-8 using a residue-wise decomposition of the binding energy. Binding of four putative inhibitors of kallikrein-8 is investigated through molecular dynamics simulation and ligand binding energy evaluation with two methods (MM/PBSA and WaterSwap). For visualization of the residue-wise decomposition of binding energies, chemical energy-wise decomposition or CHEWD is introduced as a plugin to UCSF Chimera and Pymol. CHEWD allows easy comparison between ligands using individual residue contributions to the binding energy. Molecular dynamics simulations indicate one ligand binds stably to the kallikrein-8 S1 binding site. Comparison with other members of the kallikrein family shows that residues responsible for binding are specific to kallikrein-8. Thus, ZINC02927490 is a promising lead for development of novel kallikrein-8 inhibitors.
Collapse
|
87
|
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
|
88
|
O'Hagan MP, Haldar S, Duchi M, Oliver TAA, Mulholland AJ, Morales JC, Galan MC. A Photoresponsive Stiff-Stilbene Ligand Fuels the Reversible Unfolding of G-Quadruplex DNA. Angew Chem Int Ed Engl 2019; 58:4334-4338. [PMID: 30682233 PMCID: PMC6563076 DOI: 10.1002/anie.201900740] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Indexed: 11/14/2022]
Abstract
The polymorphic nature of G-quadruplex (G4) DNA structures points to a range of potential applications in nanodevices and an opportunity to control G4 in biological settings. Light is an attractive means for the regulation of oligonucleotide structure as it can be delivered with high spatiotemporal precision. However, surprisingly little attention has been devoted towards the development of ligands for G4 that allow photoregulation of G4 folding. We report a novel G4-binding chemotype derived from stiff-stilbene. Surprisingly however, whilst the ligand induces high stabilization in the potassium form of human telomeric DNA, it causes the unfolding of the same G4 sequence in sodium buffer. This effect can be reversed on demand by irradiation with 400 nm light through deactivation of the ligand by photo-oxidation. By fuelling the system with the photolabile ligand, the conformation of G4 DNA was switched five times.
Collapse
|
89
|
O'Hagan MP, Haldar S, Duchi M, Oliver TAA, Mulholland AJ, Morales JC, Galan MC. A Photoresponsive Stiff‐Stilbene Ligand Fuels the Reversible Unfolding of G‐Quadruplex DNA. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201900740] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
90
|
Limb MAL, Suardíaz R, Grant IM, Mulholland AJ. Quantum Mechanics/Molecular Mechanics Simulations Show Saccharide Distortion is Required for Reaction in Hen Egg-White Lysozyme. Chemistry 2019; 25:764-768. [PMID: 30347479 DOI: 10.1002/chem.201805250] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Indexed: 01/07/2023]
Abstract
Hybrid quantum mechanics/molecular mechanics (QM/MM) calculations on lysozyme show significant distortion of the bound saccharide is required to facilitate the catalytic reaction.
Collapse
|
91
|
Haldar S, Comitani F, Saladino G, Woods C, van der Kamp MW, Mulholland AJ, Gervasio FL. A Multiscale Simulation Approach to Modeling Drug-Protein Binding Kinetics. J Chem Theory Comput 2018; 14:6093-6101. [PMID: 30208708 DOI: 10.1021/acs.jctc.8b00687] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Drug-target binding kinetics has recently emerged as a sometimes critical determinant of in vivo efficacy and toxicity. Its rational optimization to improve potency or reduce side effects of drugs is, however, extremely difficult. Molecular simulations can play a crucial role in identifying features and properties of small ligands and their protein targets affecting the binding kinetics, but significant challenges include the long time scales involved in (un)binding events and the limited accuracy of empirical atomistic force fields (lacking, e.g., changes in electronic polarization). In an effort to overcome these hurdles, we propose a method that combines state-of-the-art enhanced sampling simulations and quantum mechanics/molecular mechanics (QM/MM) calculations at the BLYP/VDZ level to compute association free energy profiles and characterize the binding kinetics in terms of structure and dynamics of the transition state ensemble. We test our combined approach on the binding of the anticancer drug Imatinib to Src kinase, a well-characterized target for cancer therapy with a complex binding mechanism involving significant conformational changes. The results indicate significant changes in polarization along the binding pathways, which affect the predicted binding kinetics. This is likely to be of widespread importance in binding of ligands to protein targets.
Collapse
|
92
|
Liao Q, Kulkarni Y, Sengupta U, Petrović D, Mulholland AJ, van der Kamp MW, Strodel B, Kamerlin SCL. Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case. J Am Chem Soc 2018; 140:15889-15903. [PMID: 30362343 DOI: 10.1021/jacs.8b09378] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.
Collapse
|
93
|
Zhu B, Tiller BP, Walker AJ, Mulholland AJ, Windmill JFC. "Pipe Organ" Inspired Air-Coupled Ultrasonic Transducers With Broader Bandwidth. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2018; 65:1873-1881. [PMID: 30072320 DOI: 10.1109/tuffc.2018.2861575] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Piezoelectric micromachined ultrasonic transducers (PMUTs) are used to receive and transmit ultrasonic signals in industrial and biomedical applications. This type of transducer can be miniaturized and integrated with electronic systems since each element is small and the power requirements are low. The bandwidth of the PMUT may be narrow in some conventional designs; however, it is possible to apply modified structures to enhance this. This paper presents a methodology for improving the bandwidth of air-coupled PMUTs without sensitivity loss by connecting a number of resonating pipes of various lengths to a cavity. A prototype piezoelectric diaphragm ultrasonic transducer is presented to prove the theory. This novel device was fabricated by additive manufacturing (3-D printing), and consists of a polyvinylidene fluoride thin film over a stereolithography designed backplate. The backplate design is inspired by a pipe organ musical instrument, where the resonant frequency (pitch) of each pipe is mainly determined by its length. The -6-dB bandwidth of the "pipe organ" air-coupled transducer is 55.7% and 58.5% in transmitting and receiving modes, respectively, which is ~5 times wider than a custom-built standard device.
Collapse
|
94
|
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]
|
95
|
Thomas F, Dawson WM, Lang EJM, Burton AJ, Bartlett GJ, Rhys GG, Mulholland AJ, Woolfson DN. De Novo-Designed α-Helical Barrels as Receptors for Small Molecules. ACS Synth Biol 2018; 7:1808-1816. [PMID: 29944338 DOI: 10.1021/acssynbio.8b00225] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We describe de novo-designed α-helical barrels (αHBs) that bind and discriminate between lipophilic biologically active molecules. αHBs have five or more α-helices arranged around central hydrophobic channels the diameters of which scale with oligomer state. We show that pentameric, hexameric, and heptameric αHBs bind the environmentally sensitive dye 1,6-diphenylhexatriene (DPH) in the micromolar range and fluoresce. Displacement of the dye is used to report the binding of nonfluorescent molecules: palmitic acid and retinol bind to all three αHBs with submicromolar inhibitor constants; farnesol binds the hexamer and heptamer; but β-carotene binds only the heptamer. A co-crystal structure of the hexamer with farnesol reveals oriented binding in the center of the hydrophobic channel. Charged side chains engineered into the lumen of the heptamer facilitate binding of polar ligands: a glutamate variant binds a cationic variant of DPH, and introducing lysine allows binding of the biosynthetically important farnesol diphosphate.
Collapse
|
96
|
Fritz RA, Alzate-Morales JH, Spencer J, Mulholland AJ, van der Kamp MW. Multiscale Simulations of Clavulanate Inhibition Identify the Reactive Complex in Class A β-Lactamases and Predict the Efficiency of Inhibition. Biochemistry 2018; 57:3560-3563. [PMID: 29812917 DOI: 10.1021/acs.biochem.8b00480] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Clavulanate is used as an effective drug in combination with β-lactam antibiotics to treat infections of some antibiotic resistant bacteria. Here, we perform combined quantum mechanics/molecular mechanics simulations of several covalent complexes of clavulanate with class A β-lactamases KPC-2 and TEM-1. Simulations of the deacylation reactions identify the decarboxylated trans-enamine complex as being responsible for inhibition. Further, the obtained free energy barriers discriminate clinically relevant inhibition (TEM-1) from less effective inhibition (KPC-2).
Collapse
|
97
|
Lence E, van der Kamp MW, González-Bello C, Mulholland AJ. QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes. Org Biomol Chem 2018; 16:4443-4455. [PMID: 29767194 PMCID: PMC6011038 DOI: 10.1039/c8ob00066b] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Accepted: 04/18/2018] [Indexed: 11/29/2022]
Abstract
Type II dehydroquinase enzymes (DHQ2), recognized targets for antibiotic drug discovery, show significantly different activities dependent on the species: DHQ2 from Mycobacterium tuberculosis (MtDHQ2) and Helicobacter pylori (HpDHQ2) show a 50-fold difference in catalytic efficiency. Revealing the determinants of this activity difference is important for our understanding of biological catalysis and further offers the potential to contribute to tailoring specificity in drug design. Molecular dynamics simulations using a quantum mechanics/molecular mechanics potential, with correlated ab initio single point corrections, identify and quantify the subtle determinants of the experimentally observed difference in efficiency. The rate-determining step involves the formation of an enolate intermediate: more efficient stabilization of the enolate and transition state of the key step in MtDHQ2, mainly by the essential residues Tyr24 and Arg19, makes it more efficient than HpDHQ2. Further, a water molecule, which is absent in MtDHQ2 but involved in generation of the catalytic Tyr22 tyrosinate in HpDHQ2, was found to destabilize both the transition state and the enolate intermediate. The quantification of the contribution of key residues and water molecules in the rate-determining step of the mechanism also leads to improved understanding of higher potencies and specificity of known inhibitors, which should aid ongoing inhibitor design.
Collapse
|
98
|
O’Connor M, Deeks HM, Dawn E, Metatla O, Roudaut A, Sutton M, Thomas LM, Glowacki BR, Sage R, Tew P, Wonnacott M, Bates P, Mulholland AJ, Glowacki DR. Sampling molecular conformations and dynamics in a multiuser virtual reality framework. SCIENCE ADVANCES 2018; 4:eaat2731. [PMID: 29963636 PMCID: PMC6025904 DOI: 10.1126/sciadv.aat2731] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Accepted: 05/18/2018] [Indexed: 05/28/2023]
Abstract
We describe a framework for interactive molecular dynamics in a multiuser virtual reality (VR) environment, combining rigorous cloud-mounted atomistic physics simulations with commodity VR hardware, which we have made accessible to readers (see isci.itch.io/nsb-imd). It allows users to visualize and sample, with atomic-level precision, the structures and dynamics of complex molecular structures "on the fly" and to interact with other users in the same virtual environment. A series of controlled studies, in which participants were tasked with a range of molecular manipulation goals (threading methane through a nanotube, changing helical screw sense, and tying a protein knot), quantitatively demonstrate that users within the interactive VR environment can complete sophisticated molecular modeling tasks more quickly than they can using conventional interfaces, especially for molecular pathways and structural transitions whose conformational choreographies are intrinsically three-dimensional. This framework should accelerate progress in nanoscale molecular engineering areas including conformational mapping, drug development, synthetic biology, and catalyst design. More broadly, our findings highlight the potential of VR in scientific domains where three-dimensional dynamics matter, spanning research and education.
Collapse
|
99
|
Ainsley J, Mulholland AJ, Black GW, Sparagano O, Christov CZ, Karabencheva-Christova TG. Structural Insights from Molecular Dynamics Simulations of Tryptophan 7-Halogenase and Tryptophan 5-Halogenase. ACS OMEGA 2018; 3:4847-4859. [PMID: 31458701 PMCID: PMC6641897 DOI: 10.1021/acsomega.8b00385] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2018] [Accepted: 04/24/2018] [Indexed: 05/08/2023]
Abstract
Many natural organic compounds with pharmaceutical applications, including antibiotics (chlortetracycline and vancomycin), antifungal compounds (pyrrolnitrin), and chemotherapeutics (salinosporamide A and rebeccamycin) are chlorinated. Halogenating enzymes like tryptophan 7-halogenase (PrnA) and tryptophan 5-halogenase (PyrH) perform regioselective halogenation of tryptophan. In this study, the conformational dynamics of two flavin-dependent tryptophan halogenases-PrnA and PyrH-was investigated through molecular dynamics simulations, which are in agreement with the crystallographic and kinetic experimental studies of both enzymes and provide further explanation of the experimental data at an atomistic level of accuracy. They show that the binding sites of the cofactor-flavin adenine dinucleotide and the substrate do not come into close proximity during the simulations, thus supporting an enzymatic mechanism without a direct contact between them. Two catalytically important active site residues, glutamate (E346/E354) and lysine (K79/K75) in PrnA and PyrH, respectively, were found to play a key role in positioning the proposed chlorinating agent, hypochlorous acid. The changes in the regioselectivity between PrnA and PyrH arise as a consequence of differences in the orientation of substrate in its binding site.
Collapse
|
100
|
Rinaldi S, Van der Kamp MW, Ranaghan KE, Mulholland AJ, Colombo G. Understanding Complex Mechanisms of Enzyme Reactivity: The Case of Limonene-1,2-Epoxide Hydrolases. ACS Catal 2018. [DOI: 10.1021/acscatal.8b00863] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|