Molecular basis of hyper-thermostability in the thermophilic archaeal aldolase MfnB.

Methanogenic archaea are chemolithotrophic prokaryotes that can reduce carbon dioxide with hydrogen gas to form methane. These microorganisms make a significant contribution to the global carbon cycle, with methanogenic archaea from anoxic environments estimated to contribute > 500 million tons of global methane annually. Archaeal methanogenesis is dependent on the methanofurans; aminomethylfuran containing coenzymes that act as the primary C1 acceptor molecule during carbon dioxide fixation. Although the biosynthetic pathway to the methanofurans has been elucidated, structural adaptations which confer thermotolerance to Mfn enzymes from extremophilic archaea are yet to be investigated. Here we focus on the methanofuran biosynthetic enzyme MfnB, which catalyses the condensation of two molecules of glyceralde-3-phosphate to form 4­(hydroxymethyl)-2-furancarboxaldehyde-phosphate. In this study, MfnB enzymes from the hyperthermophile Methanocaldococcus jannaschii and the mesophile Methanococcus maripaludis have been recombinantly overexpressed and purified to homogeneity. Thermal unfolding studies, together with steady-state kinetic assays, demonstrate thermoadaptation in the M. jannaschii enzyme. Molecular dynamics simulations have been used to provide a structural explanation for the observed properties. These reveal a greater number of side chain interactions in the M. jannaschii enzyme, which may confer protection from heating effects by enforcing spatial residue constraints.

Antibiotic origami: selective formation of spirotetronates in abyssomicin biosynthesis.

The abyssomicins are a structurally intriguing family of bioactive natural products that include compounds with potent antibacterial, antitumour and antiviral activities. The biosynthesis of the characteristic abyssomicin spirotetronate core occurs via an enzyme-catalysed intramolecular Diels-Alder reaction, which proceeds via one of two distinct stereochemical pathways to generate products differing in configuration at the C15 spirocentre. Using the purified spirotetronate cyclases AbyU (from abyssomicin C/atrop-abyssomicin C biosynthesis) and AbmU (from abyssomicin 2/neoabyssomicin biosynthesis), in combination with synthetic substrate analogues, here we show that stereoselectivity in the spirotetronate-forming [4 + 2]-cycloaddition is controlled by a combination of factors attributable to both the enzyme and substrate. Furthermore, an achiral substrate was enzymatically cyclised to a single enantiomer of a spirocyclic product. X-ray crystal structures, molecular dynamics simulations, and assessment of substrate binding affinity and reactivity in both AbyU and AbmU establish the molecular determinants of stereochemical control in this important class of biocatalysts.

Simulation-Guided Engineering Enables a Functional Switch in Selinadiene Synthase toward Hydroxylation.

Engineering sesquiterpene synthases to form predefined alternative products is a major challenge due to their diversity in cyclization mechanisms and our limited understanding of how amino acid changes affect the steering of these mechanisms. Here, we use a combination of atomistic simulation and site-directed mutagenesis to engineer a selina-4(15),7(11)-diene synthase (SdS) such that its final reactive carbocation is quenched by trapped active site water, resulting in the formation of a complex hydroxylated sesquiterpene (selin-7(11)-en-4-ol). Initially, the SdS G305E variant produced 20% selin-7(11)-en-4-ol. As suggested by modeling of the enzyme-carbocation complex, selin-7(11)-en-4-ol production could be further improved by varying the pH, resulting in selin-7(11)-en-4-ol becoming the major product (48%) at pH 6.0. We incorporated the SdS G305E variant along with genes from the mevalonate pathway into bacterial BL21(DE3) cells and demonstrated the production of selin-7(11)-en-4-ol at a scale of 10 mg/L in batch fermentation. These results highlight opportunities for the simulation-guided engineering of terpene synthases to produce predefined complex hydroxylated sesquiterpenes.

Delineation of the complete reaction cycle of a natural Diels-Alderase.

The Diels-Alder reaction is one of the most effective methods for the synthesis of substituted cyclohexenes. The development of protein catalysts for this reaction remains a major priority, affording new sustainable routes to high value target molecules. Whilst a small number of natural enzymes have been shown capable of catalysing [4 + 2] cycloadditions, there is a need for significant mechanistic understanding of how these prospective Diels-Alderases promote catalysis to underpin their development as biocatalysts for use in synthesis. Here we present a molecular description of the complete reaction cycle of the bona fide natural Diels-Alderase AbyU, which catalyses formation of the spirotetronate skeleton of the antibiotic abyssomicin C. This description is derived from X-ray crystallographic studies of AbyU in complex with a non-transformable synthetic substrate analogue, together with transient kinetic analyses of the AbyU catalysed reaction and computational reaction simulations. These studies reveal the mechanistic intricacies of this enzyme system and establish a foundation for the informed reengineering of AbyU and related biocatalysts.

Epistasis arises from shifting the rate-limiting step during enzyme evolution of a ß-lactamase.

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a ß-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity.

Electric Fields Are a Key Determinant of Carbapenemase Activity in Class A ß-Lactamases.

Resistance to antibiotics is a public health crisis. Although carbapenems are less susceptible to resistance than other ß-lactam antibiotics, ß-lactamases mediating resistance against these drugs are spreading. Here, we dissect the contributions of electric fields to carbapenemase activity in class A ß-lactamases. We perform QM/MM molecular dynamics simulations of meropenem acyl-enzyme hydrolysis that correctly discriminate carbapenemases. Electric field analysis shows that active-site fields in the deacylation transition state and tetrahedral intermediate are important determinants of activity. The active-site fields identify several residues, some distal, that distinguish efficient carbapenemases. Our field analysis script (www.github.com/bunzela/FieldTools) may help in understanding and combating antibiotic resistance.

emle-engine: A Flexible Electrostatic Machine Learning Embedding Package for Multiscale Molecular Dynamics Simulations.

We present in this work the emle-engine package (https://github.com/chemle/emle-engine)─the implementation of a new machine learning embedding scheme for hybrid machine learning potential/molecular-mechanics (ML/MM) dynamics simulations. The package is based on an embedding scheme that uses a physics-based model of the electronic density and induction with a handful of tunable parameters derived from in vacuo properties of the subsystem to be embedded. This scheme is completely independent of the in vacuo potential and requires only the positions of the atoms of the machine learning subsystem and the positions and partial charges of the molecular mechanics environment. These characteristics allow emle-engine to be employed in existing QM/MM software. We demonstrate that the implemented electrostatic machine learning embedding scheme (named EMLE) is stable in enhanced sampling molecular dynamics simulations. Through the calculation of free energy surfaces of alanine dipeptide in water with two different ML options for the in vacuo potential and three embedding models, we test the performance of EMLE. When compared to the reference DFT/MM surface, the EMLE embedding is clearly superior to the MM one based on fixed partial charges. The configurational dependence of the electronic density and the inclusion of the induction energy introduced by the EMLE model leads to a systematic reduction in the average error of the free energy surface when compared to MM embedding. By enabling the usage of EMLE embedding in practical ML/MM simulations, emle-engine will make it possible to accurately model systems and processes that feature significant variations in the charge distribution of the ML subsystem and/or the interacting environment.

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis.

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis.

Active Site Loop Engineering Abolishes Water Capture in Hydroxylating Sesquiterpene Synthases.

Terpene synthases (TS) catalyze complex reactions to produce a diverse array of terpene skeletons from linear isoprenyl diphosphates. Patchoulol synthase (PTS) from Pogostemon cablin converts farnesyl diphosphate into patchoulol. Using simulation-guided engineering, we obtained PTS variants that eliminate water capture. Further, we demonstrate that modifying the structurally conserved Hα-1 loop also reduces hydroxylation in PTS, as well as in germacradiene-11-ol synthase (Gd11olS), leading to cyclic neutral intermediates as products, including α-bulnesene (PTS) and isolepidozene (Gd11olS). Hα-1 loop modification could be a general strategy for engineering sesquiterpene synthases to produce complex cyclic hydrocarbons without the need for structure determination or modeling.

Free energy along drug-protein binding pathways interactively sampled in virtual reality.

We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.

Promiscuous recognition of MR1 drives self-reactive mucosal-associated invariant T cell responses.

Mucosal-associated invariant T (MAIT) cells use canonical semi-invariant T cell receptors (TCR) to recognize microbial riboflavin precursors displayed by the antigen-presenting molecule MR1. The extent of MAIT TCR crossreactivity toward physiological, microbially unrelated antigens remains underexplored. We describe MAIT TCRs endowed with MR1-dependent reactivity to tumor and healthy cells in the absence of microbial metabolites. MAIT cells bearing TCRs crossreactive toward self are rare but commonly found within healthy donors and display T-helper-like functions in vitro. Experiments with MR1-tetramers loaded with distinct ligands revealed significant crossreactivity among MAIT TCRs both ex vivo and upon in vitro expansion. A canonical MAIT TCR was selected on the basis of extremely promiscuous MR1 recognition. Structural and molecular dynamic analyses associated promiscuity to unique TCRß-chain features that were enriched within self-reactive MAIT cells of healthy individuals. Thus, self-reactive recognition of MR1 represents a functionally relevant indication of MAIT TCR crossreactivity, suggesting a potentially broader role of MAIT cells in immune homeostasis and diseases, beyond microbial immunosurveillance.

Interrogation of an Enzyme Library Reveals the Catalytic Plasticity of Naturally Evolved [4+2] Cyclases.

Stereoselective carbon-carbon bond forming reactions are quintessential transformations in organic synthesis. One example is the Diels-Alder reaction, a [4+2] cycloaddition between a conjugated diene and a dienophile to form cyclohexenes. The development of biocatalysts for this reaction is paramount for unlocking sustainable routes to a plethora of important molecules. To obtain a comprehensive understanding of naturally evolved [4+2] cyclases, and to identify hitherto uncharacterised biocatalysts for this reaction, we constructed a library comprising forty-five enzymes with reported or predicted [4+2] cycloaddition activity. Thirty-one library members were successfully produced in recombinant form. In vitro assays employing a synthetic substrate incorporating a diene and a dienophile revealed broad-ranging cycloaddition activity amongst these polypeptides. The hypothetical protein Cyc15 was found to catalyse an intramolecular cycloaddition to generate a novel spirotetronate. The crystal structure of this enzyme, along with docking studies, establishes the basis for stereoselectivity in Cyc15, as compared to other spirotetronate cyclases.

Synthetic cannabinoid receptor agonists are monoamine oxidase-A selective inhibitors.

Synthetic cannabinoid receptor agonists (SCRAs) are one of the fastest growing classes of recreational drugs. Despite their growth in use, their vast chemical diversity and rapidly changing landscape of structures make understanding their effects challenging. In particular, the side effects for SCRA use are extremely diverse, but notably include severe outcomes such as cardiac arrest. These side effects appear at odds with the main putative mode of action, as full agonists of cannabinoid receptors. We have hypothesized that SCRAs may act as MAO inhibitors, owing to their structural similarity to known monoamine oxidase inhibitors (MAOI's) as well as matching clinical outcomes (hypertensive crisis) of 'monoaminergic toxicity' for users of MAOIs and some SCRA use. We have studied the potential for SCRA-mediated inhibition of MAO-A and MAO-B via a range of SCRAs used commonly in the UK, as well as structural analogues to prove the atomistic determinants of inhibition. By combining in silico and experimental kinetic studies we demonstrate that SCRAs are MAO-A-specific inhibitors and their affinity can vary significantly between SCRAs, most notably affected by the nature of the SCRA 'head' group. Our data allow us to posit a putative mechanism of inhibition. Crucially our data demonstrate that SCRA activity is not limited to just cannabinoid receptor agonism and that alternative interactions might account for some of the diversity of the observed side effects and that these effects can be SCRA-specific.

The Role of Cytochrome†P450 AbyV in the Final Stages of Abyssomicin†C Biosynthesis.

Abyssomicin C and its atropisomer are potent inhibitors of bacterial folate metabolism. They possess complex polycyclic structures, and their biosynthesis has been shown to involve several unusual enzymatic transformations. Using a combination of synthesis and in vitro assays we reveal that AbyV, a cytochrome P450 enzyme from the aby gene cluster, catalyses a key late-stage epoxidation required for the installation of the characteristic ether-bridged core of abyssomicin C. The X-ray crystal structure of AbyV has been determined, which in combination with molecular dynamics simulations provides a structural framework for our functional data. This work demonstrates the power of combining selective carbon-13 labelling with NMR spectroscopy as a sensitive tool to interrogate enzyme-catalysed reactions in vitro with no need for purification.

The Role of Cytochrome†P450 AbyV in the Final Stages of Abyssomicin†C Biosynthesis.

Abyssomicin C and its atropisomer are potent inhibitors of bacterial folate metabolism. They possess complex polycyclic structures, and their biosynthesis has been shown to involve several unusual enzymatic transformations. Using a combination of synthesis and in vitro assays we reveal that AbyV, a cytochrome P450 enzyme from the aby gene cluster, catalyses a key late-stage epoxidation required for the installation of the characteristic ether-bridged core of abyssomicin C. The X-ray crystal structure of AbyV has been determined, which in combination with molecular dynamics simulations provides a structural framework for our functional data. This work demonstrates the power of combining selective carbon-13 labelling with NMR spectroscopy as a sensitive tool to interrogate enzyme-catalysed reactions in vitro with no need for purification.

Structurally Informed Mutagenesis of a Stereochemically Promiscuous Aldolase Produces Mutants That Catalyze the Diastereoselective Syntheses of All Four Stereoisomers of 3-Deoxy-hexulosonic Acid.

A 2-keto-3-deoxygluconate aldolase from the hyperthermophile Sulfolobus solfataricus catalyzes the nonstereoselective aldol reaction of pyruvate and d-glyceraldehyde to produce 2-keto-3-deoxygluconate (d-KDGlc) and 2-keto-3-deoxy-d-galactonate (d-KDGal). Previous investigations into curing the stereochemical promiscuity of this hyperstable aldolase used high-resolution structures of the aldolase bound to d-KDGlc or d-KDGal to identify critical amino acids involved in substrate binding for mutation. This structure-guided approach enabled mutant variants to be created that could stereoselectively catalyze the aldol reaction of pyruvate and natural d-glyceraldehyde to selectively afford d-KDGlc or d-KDGal. Here we describe the creation of two further mutants of this Sulfolobus aldolase that can be used to catalyze aldol reactions between pyruvate and non-natural l-glyceraldehyde to enable the diastereoselective synthesis of l-KDGlc and l-KDGal. High-resolution crystal structures of all four variant aldolases have been determined (both unliganded and liganded), including Variant 1 with d-KDGlc, Variant 2 with pyruvate, Variant 3 with l-KDGlc, and Variant 4 with l-KDGal. These structures have enabled us to rationalize the observed changes in diastereoselectivities in these variant-catalyzed aldol reactions at a molecular level. Interestingly, the active site of Variant 4 was found to be sufficiently flexible to enable catalytically important amino acids to be replaced while still retaining sufficient enzymic activity to enable production of l-KDGal.

QM/MM Simulations Reveal the Determinants of Carbapenemase Activity in Class A ß-Lactamases.

ß-lactam antibiotic resistance in Gram-negative bacteria, primarily caused by ß-lactamase enzymes that hydrolyze the ß-lactam ring, has become a serious clinical problem. Carbapenems were formerly considered "last resort" antibiotics because they escaped breakdown by most ß-lactamases, due to slow deacylation of the acyl-enzyme intermediate. However, an increasing number of Gram-negative bacteria now produce ß-lactamases with carbapenemase activity: these efficiently hydrolyze the carbapenem ß-lactam ring, severely limiting the treatment of some bacterial infections. Here, we use quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reactions of acyl-enzyme complexes of eight ß-lactamases of class A (the most widely distributed ß-lactamase group) with the carbapenem meropenem to investigate differences between those inhibited by carbapenems (TEM-1, SHV-1, BlaC, and CTX-M-16) and those that hydrolyze them (SFC-1, KPC-2, NMC-A, and SME-1). QM/MM molecular dynamics simulations confirm the two enzyme groups to differ in the preferred acyl-enzyme orientation: carbapenem-inhibited enzymes favor hydrogen bonding of the carbapenem hydroxyethyl group to deacylating water (DW). QM/MM simulations of deacylation give activation free energies in good agreement with experimental hydrolysis rates, correctly distinguishing carbapenemases. For the carbapenem-inhibited enzymes, free energies for deacylation are significantly higher than for the carbapenemases, even when the hydroxyethyl group was restrained to prevent interaction with the DW. Analysis of these simulations, and additional simulations of mutant enzymes, shows how factors including the hydroxyethyl orientation, the active site volume, and architecture (conformations of Asn170 and Asn132; organization of the oxyanion hole; and the Cys69-Cys238 disulfide bond) collectively determine catalytic efficiency toward carbapenems.

Path to Actinorhodin: Regio- and Stereoselective Ketone Reduction by a Type II Polyketide Ketoreductase Revealed in Atomistic Detail.

In type II polyketide synthases (PKSs), which typically biosynthesize several antibiotic and antitumor compounds, the substrate is a growing polyketide chain, shuttled between individual PKS enzymes, while covalently tethered to an acyl carrier protein (ACP): this requires the ACP interacting with a series of different enzymes in succession. During biosynthesis of the antibiotic actinorhodin, produced by Streptomyces coelicolor, one such key binding event is between an ACP carrying a 16-carbon octaketide chain (actACP) and a ketoreductase (actKR). Once the octaketide is bound inside actKR, it is likely cyclized between C7 and C12 and regioselective reduction of the ketone at C9 occurs: how these elegant chemical and conformational changes are controlled is not yet known. Here, we perform protein-protein docking, protein NMR, and extensive molecular dynamics simulations to reveal a probable mode of association between actACP and actKR; we obtain and analyze a detailed model of the C7-C12-cyclized octaketide within the actKR active site; and we confirm this model through multiscale (QM/MM) reaction simulations of the key ketoreduction step. Molecular dynamics simulations show that the most thermodynamically stable cyclized octaketide isomer (7R,12R) also gives rise to the most reaction competent conformations for ketoreduction. Subsequent reaction simulations show that ketoreduction is stereoselective as well as regioselective, resulting in an S-alcohol. Our simulations further indicate several conserved residues that may be involved in selectivity of C7-12 cyclization and C9 ketoreduction. Detailed insights obtained on ACP-based substrate presentation in type II PKSs can help design ACP-ketoreductase systems with altered regio- or stereoselectivity.

Multiscale Simulations Identify Origins of Differential Carbapenem Hydrolysis by the OXA-48 ß-Lactamase.

OXA-48 ß-lactamases are frequently encountered in bacterial infections caused by carbapenem-resistant Gram-negative bacteria. Due to the importance of carbapenems in the treatment of healthcare-associated infections and the increasingly wide dissemination of OXA-48-like enzymes on plasmids, these ß-lactamases are of high clinical significance. Notably, OXA-48 hydrolyzes imipenem more efficiently than other commonly used carbapenems, such as meropenem. Here, we use extensive multiscale simulations of imipenem and meropenem hydrolysis by OXA-48 to dissect the dynamics and to explore differences in the reactivity of the possible conformational substates of the respective acylenzymes. Quantum mechanics/molecular mechanics (QM/MM) simulations of the deacylation reaction for both substrates demonstrate that deacylation is favored when the 6α-hydroxyethyl group is able to hydrogen bond to the water molecule responsible for deacylation but disfavored by the increasing hydration of either oxygen of the carboxylated Lys73 general base. Differences in free energy barriers calculated from the QM/MM simulations correlate well with the experimentally observed differences in hydrolytic efficiency between meropenem and imipenem. We conclude that the impaired breakdown of meropenem, compared to imipenem, which arises from a subtle change in the hydrogen bonding pattern between the deacylating water molecule and the antibiotic, is most likely induced by the meropenem 1ß-methyl group. In addition to increased insights into carbapenem breakdown by OXA ß-lactamases, which may aid in future efforts to design antibiotics or inhibitors, our approach exemplifies the combined use of atomistic simulations in determining the possible different enzyme-substrate substates and their influence on enzyme reaction kinetics.