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.

Visualizing the protons in a metalloenzyme electron proton transfer pathway.

In redox metalloenzymes, the process of electron transfer often involves the concerted movement of a proton. These processes are referred to as proton-coupled electron transfer, and they underpin a wide variety of biological processes, including respiration, energy conversion, photosynthesis, and metalloenzyme catalysis. The mechanisms of proton delivery are incompletely understood, in part due to an absence of information on exact proton locations and hydrogen bonding structures in a bona fide metalloenzyme proton pathway. Here, we present a 2.1-Å neutron crystal structure of the complex formed between a redox metalloenzyme (ascorbate peroxidase) and its reducing substrate (ascorbate). In the neutron structure of the complex, the protonation states of the electron/proton donor (ascorbate) and all of the residues involved in the electron/proton transfer pathway are directly observed. This information sheds light on possible proton movements during heme-catalyzed oxygen activation, as well as on ascorbate oxidation.

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.

Conformation control of the histidine kinase BceS of Bacillus subtilis by its cognate ABC-transporter facilitates need-based activation of antibiotic resistance.

Bacteria closely control gene expression to ensure optimal physiological responses to their environment. Such careful gene expression can minimize the fitness cost associated with antibiotic resistance. We previously described a novel regulatory logic in Bacillus subtilis enabling the cell to directly monitor its need for detoxification. This cost-effective strategy is achieved via a two-component regulatory system (BceRS) working in a sensory complex with an ABC-transporter (BceAB), together acting as a flux-sensor where signaling is proportional to transport activity. How this is realized at the molecular level has remained unknown. Using experimentation and computation we here show that the histidine kinase is activated by piston-like displacements in the membrane, which are converted to helical rotations in the catalytic core via an intervening HAMP-like domain. Intriguingly, the transporter was not only required for kinase activation, but also to actively maintain the kinase in its inactive state in the absence of antibiotics. Such coupling of kinase activity to that of the transporter ensures the complete control required for transport flux-dependent signaling. Moreover, we show that the transporter likely conserves energy by signaling with sub-maximal sensitivity. These results provide the first mechanistic insights into transport flux-dependent signaling, a unique strategy for energy-efficient decision making.

Mitochondrial import, health and mtDNA copy number variability seen when using type II and type V CRISPR effectors.

Current methodologies for targeting the mitochondrial genome for research and/or therapy development in mitochondrial diseases are restricted by practical limitations and technical inflexibility. A molecular toolbox for CRISPR-mediated mitochondrial genome editing is desirable, as this could enable targeting of mtDNA haplotypes using the precision and tuneability of CRISPR enzymes. Such 'MitoCRISPR' systems described to date lack reproducibility and independent corroboration. We have explored the requirements for MitoCRISPR in human cells by CRISPR nuclease engineering, including the use of alternative mitochondrial protein targeting sequences and smaller paralogues, and the application of guide (g)RNA modifications for mitochondrial import. We demonstrate varied mitochondrial targeting efficiencies and effects on mitochondrial dynamics/function of different CRISPR nucleases, with Lachnospiraceae bacterium ND2006 (Lb) Cas12a being better targeted and tolerated than Cas9 variants. We also provide evidence of Cas9 gRNA association with mitochondria in HeLa cells and isolated yeast mitochondria, even in the absence of a targeting RNA aptamer. Our data link mitochondrial-targeted LbCas12a/crRNA with increased mtDNA copy number dependent upon DNA binding and cleavage activity. We discuss reproducibility issues and the future steps necessary for MitoCRISPR.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases.

Monoterpene synthases are often promiscuous enzymes, yielding product mixtures rather than pure compounds due to the nature of the branched reaction mechanism involving reactive carbocations. Two previously identified bacterial monoterpene synthases, a linalool synthase (bLinS) and a cineole synthase (bCinS), produce nearly pure linalool and cineole from geranyl diphosphate, respectively. We used a combined experimental and computational approach to identify critical residues involved in bacterial monoterpenoid synthesis. Phe77 is essential for bCinS activity, guiding the linear carbocation intermediate towards the formation of the cyclic α-terpinyl intermediate; removal of the aromatic ring results in variants that produce acyclic products only. Computational chemistry confirmed the importance of Phe77 in carbocation stabilisation. Phe74, Phe78 and Phe179 are involved in maintaining the active site shape in bCinS without a specific role for the aromatic ring. Phe295 in bLinS, and the equivalent Ala301 in bCinS, are essential for linalool and cineole formation, respectively. Where Phe295 places steric constraints on the carbocation intermediates, Ala301 is essential for bCinS initial cyclisation and activity. Our multidisciplinary approach gives unique insights into how carefully placed amino acid residues in the active site can direct carbocations down specific paths, by placing steric constraints or offering stabilisation via cation-π interactions.

Reliable In Silico Ranking of Engineered Therapeutic TCR Binding Affinities with MMPB/GBSA.

Accurate and efficient in silico ranking of protein-protein binding affinities is useful for protein design with applications in biological therapeutics. One popular approach to rank binding affinities is to apply the molecular mechanics Poisson-Boltzmann/generalized Born surface area (MMPB/GBSA) method to molecular dynamics (MD) trajectories. Here, we identify protocols that enable the reliable evaluation of T-cell receptor (TCR) variants binding to their target, peptide-human leukocyte antigens (pHLAs). We suggest different protocols for variant sets with a few (≤4) or many mutations, with entropy corrections important for the latter. We demonstrate how potential outliers could be identified in advance and that just 5-10 replicas of short (4 ns) MD simulations may be sufficient for the reproducible and accurate ranking of TCR variants. The protocols developed here can be applied toward in silico screening during the optimization of therapeutic TCRs, potentially reducing both the cost and time taken for biologic development.

Loop Dynamics and Enzyme Catalysis in Protein Tyrosine Phosphatases.

Protein tyrosine phosphatases (PTPs) play an important role in cellular signaling and have been implicated in human cancers, diabetes, and obesity. Despite shared catalytic mechanisms and transition states for the chemical steps of catalysis, catalytic rates within the PTP family vary over several orders of magnitude. These rate differences have been implied to arise from differing conformational dynamics of the closure of a protein loop, the WPD-loop, which carries a catalytically critical residue. The present work reports computational studies of the human protein tyrosine phosphatase 1B (PTP1B) and YopH from Yersinia pestis, for which NMR has demonstrated a link between their respective rates of WPD-loop motion and catalysis rates, which differ by an order of magnitude. We have performed detailed structural analysis, both conventional and enhanced sampling simulations of their loop dynamics, as well as empirical valence bond simulations of the chemical step of catalysis. These analyses revealed the key residues and structural features responsible for these differences, as well as the residues and pathways that facilitate allosteric communication in these enzymes. Curiously, our wild-type YopH simulations also identify a catalytically incompetent hyper-open conformation of its WPD-loop, sampled as a rare event, previously only experimentally observed in YopH-based chimeras. The effect of differences within the WPD-loop and its neighboring loops on the modulation of loop dynamics, as revealed in this work, may provide a facile means for the family of PTP enzymes to respond to environmental changes and regulate their catalytic activities.

Antimicrobial Resistance Conferred by OXA-48 ß-Lactamases: Towards a Detailed Mechanistic Understanding.

OXA-48-type ß-lactamases are now routinely encountered in bacterial infections caused by carbapenem-resistant Enterobacterales These enzymes are of high and growing clinical significance due to the importance of carbapenems in treatment of health care-associated infections by Gram-negative bacteria, the wide and increasing dissemination of OXA-48 enzymes on plasmids, and the challenges posed by their detection. OXA-48 confers resistance to penicillin (which is efficiently hydrolyzed) and carbapenem antibiotics (which is more slowly broken down). In addition to the parent enzyme, a growing array of variants of OXA-48 is now emerging. The spectrum of activity of these variants varies, with some hydrolyzing expanded-spectrum oxyimino-cephalosporins. The growth in importance and diversity of the OXA-48 group has motivated increasing numbers of studies that aim to elucidate the relationship between structure and specificity and establish the mechanistic basis for ß-lactam turnover in this enzyme family. In this review, we collate recently published structural, kinetic, and mechanistic information on the interactions between clinically relevant ß-lactam antibiotics and inhibitors and OXA-48 ß-lactamases. Collectively, these studies are starting to form a detailed picture of the underlying bases for the differences in ß-lactam specificity between OXA-48 variants and the consequent differences in resistance phenotype. We focus specifically on aspects of carbapenemase and cephalosporinase activities of OXA-48 ß-lactamases and discuss ß-lactamase inhibitor development in this context. Throughout the review, we also outline key open research questions for future investigation.

Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible.

MOTIVATION: Experimental structural data can allow detailed insight into protein structure and protein-ligand interactions, which is crucial for many areas of bioscience, including drug design and enzyme engineering. Typically, however, little more than a static picture of protein-ligand interactions is obtained, whereas dynamical information is often required for deeper understanding and to assess the effect of mutations. Molecular dynamics (MD) simulations can provide such information, but setting up and running these simulations is not straightforward and requires expert knowledge. There is thus a need for a tool that makes protein-ligand simulation easily accessible to non-expert users. RESULTS: We present Enlighten2: efficient simulation protocols for protein-ligand systems alongside a user-friendly plugin to the popular visualization program PyMOL. With Enlighten2, non-expert users can straightforwardly run and visualize MD simulations on protein-ligand models of interest. There is no need to learn new programs and all underlying tools are free and open source. AVAILABILITY AND IMPLEMENTATION: The Enlighten2 Python package and PyMOL plugin are free to use under the GPL3.0 licence and can be found at https://enlighten2.github.io. We also provide a lightweight Docker image via DockerHub that includes Enlighten2 with all the required utilities.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins.

Evaluating the availability of molecular oxygen (O2 ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O2 accessibility to the ß-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O2 pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O2 accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Taming the Reactivity of Monoterpene Synthases To Guide Regioselective Product Hydroxylation.

Monoterpenoids are industrially important natural products with applications in the flavours, fragrances, fuels and pharmaceutical industries. Most monoterpenoids are produced by plants, but recently two bacterial monoterpene synthases have been identified, including a cineole synthase (bCinS). Unlike plant cineole synthases, bCinS is capable of producing nearly pure cineole from geranyl diphosphate in a complex cyclisation cascade that is tightly controlled. Here we have used a multidisciplinary approach to show that Asn305 controls water attack on the α-terpinyl cation and subsequent cyclisation and deprotonation of the α-terpineol intermediate, key steps in the cyclisation cascade which direct product formation towards cineole. Mutation of Asn305 results in variants that no longer produce α-terpineol or cineole. Molecular dynamics simulations revealed that water coordination is disrupted in all variants tested. Quantum mechanics calculations indicate that Asn305 is most likely a (transient) proton acceptor for the final deprotonation step. Our synergistic approach gives unique insight into how a single residue, Asn305, tames the promiscuous chemistry of monoterpene synthase cyclisation cascades. It does this by tightly controlling the final steps in cineole formation catalysed by bCinS to form a single hydroxylated monoterpene product.

Exposing the Interplay Between Enzyme Turnover, Protein Dynamics, and the Membrane Environment in Monoamine Oxidase B.

There is an increasing realization that structure-based drug design may show improved success by understanding the ensemble of conformations accessible to an enzyme and how the environment affects this ensemble. Human monoamine oxidase B (MAO-B) catalyzes the oxidation of amines and is inhibited for the treatment of both Parkinson's disease and depression. Despite its clinical importance, its catalytic mechanism remains unclear, and routes to drugging this target would be valuable. Evidence of a radical in either the transition state or the resting state of MAO-B is present throughout the literature and is suggested to be a flavin semiquinone, a tyrosyl radical, or both. Here we see evidence of a resting-state flavin semiquinone, via absorption redox studies and electron paramagnetic resonance, suggesting that the anionic semiquinone is biologically relevant. On the basis of enzyme kinetic studies, enzyme variants, and molecular dynamics simulations, we find evidence for the importance of the membrane environment in mediating the activity of MAO-B and that this mediation is related to the protein dynamics of MAO-B. Further, our MD simulations identify a hitherto undescribed entrance for substrate binding, membrane modulated substrate access, and indications for half-site reactivity: only one active site is accessible to binding at a time. Our study combines both experimental and computational evidence to illustrate the subtle interplay between enzyme activity and protein dynamics and the immediate membrane environment. Understanding key biomedical enzymes to this level of detail will be crucial to inform strategies (and binding sites) for rational drug design for these targets.

An Efficient Computational Assay for ß-Lactam Antibiotic Breakdown by Class A ß-Lactamases.

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.

Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site.

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.

An Esterase-like Lyase Catalyzes Acetate Elimination in Spirotetronate/Spirotetramate Biosynthesis.

Spirotetronate and spirotetramate natural products include a multitude of compounds with potent antimicrobial and antitumor activities. Their biosynthesis incorporates many unusual biocatalytic steps, including regio- and stereo-specific modifications, cyclizations promoted by Diels-Alderases, and acetylation-elimination reactions. Here we focus on the acetate elimination catalyzed by AbyA5, implicated in the formation of the key Diels-Alder substrate to give the spirocyclic system of the antibiotic abyssomicin C. Using synthetic substrate analogues, it is shown that AbyA5 catalyzes stereospecific acetate elimination, establishing the (R)-tetronate acetate as a biosynthetic intermediate. The X-ray crystal structure of AbyA5, the first of an acetate-eliminating enzyme, reveals a deviant acetyl esterase fold. Molecular dynamics simulations and enzyme assays show the use of a His-Ser dyad to catalyze either elimination or hydrolysis, via disparate mechanisms, under substrate control.

Multiscale Simulations of Clavulanate Inhibition Identify the Reactive Complex in Class A ß-Lactamases and Predict the Efficiency of Inhibition.

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).

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case.

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.

QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes.

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.

Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase.

Creatininase catalyzes the conversion of creatinine (a biosensor for kidney function) to creatine via a two-step mechanism: water addition followed by ring opening. Water addition is common to other known cyclic amidohydrolases, but the precise mechanism for ring opening is still under debate. The proton donor in this step is either His178 or a water molecule bound to one of the metal ions, and the roles of His178 and Glu122 are unclear. Here, the two possible reaction pathways have been fully examined by means of combined quantum mechanics/molecular mechanics simulations at the SCC-DFTB/CHARMM22 level of theory. The results indicate that His178 is the main catalytic residue for the whole reaction and explain its role as proton shuttle during the ring-opening step. In the first step, His178 provides electrostatic stabilization to the gem-diolate tetrahedral intermediate. In the second step, His178 abstracts the hydroxyl proton of the intermediate and delivers it to the cyclic amide nitrogen, leading to ring opening. The latter is the rate-limiting step with a free energy barrier of 18.5 kcal/mol, in agreement with the experiment. We find that Glu122 must be protonated during the enzyme reaction, so that it can form a stable hydrogen bond with its neighboring water molecule. Simulations of the E122Q mutant showed that this replacement disrupts the H-bond network formed by three conserved residues (Glu34, Ser78, and Glu122) and water, increasing the energy barrier. Our computational studies provide a comprehensive explanation for previous structural and kinetic observations, including why the H178A mutation causes a complete loss of activity but the E122Q mutation does not.