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.

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.

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.

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.

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.

Molecular Rules Underpinning Enhanced Affinity Binding of Human T Cell Receptors Engineered for Immunotherapy.

Immuno-oncology approaches that utilize T cell receptors (TCRs) are becoming highly attractive because of their potential to target virtually all cellular proteins, including cancer-specific epitopes, via the recognition of peptide-human leukocyte antigen (pHLA) complexes presented at the cell surface. However, because natural TCRs generally recognize cancer-derived pHLAs with very weak affinities, efforts have been made to enhance their binding strength, in some cases by several million-fold. In this study, we investigated the mechanisms underpinning human TCR affinity enhancement by comparing the crystal structures of engineered enhanced affinity TCRs with those of their wild-type progenitors. Additionally, we performed molecular dynamics simulations to better understand the energetic mechanisms driving the affinity enhancements. These data demonstrate that supra-physiological binding affinities can be achieved without altering native TCR-pHLA binding modes via relatively subtle modifications to the interface contacts, often driven through the addition of buried hydrophobic residues. Individual energetic components of the TCR-pHLA interaction governing affinity enhancements were distinct and highly variable for each TCR, often resulting from additive, or knock-on, effects beyond the mutated residues. This comprehensive analysis of affinity-enhanced TCRs has important implications for the future rational design of engineered TCRs as efficacious and safe drugs for cancer treatment.

Specificity of bispecific T cell receptors and antibodies targeting peptide-HLA.

Tumor-associated peptide-human leukocyte antigen complexes (pHLAs) represent the largest pool of cell surface-expressed cancer-specific epitopes, making them attractive targets for cancer therapies. Soluble bispecific molecules that incorporate an anti-CD3 effector function are being developed to redirect T cells against these targets using 2 different approaches. The first achieves pHLA recognition via affinity-enhanced versions of natural TCRs (e.g., immune-mobilizing monoclonal T cell receptors against cancer [ImmTAC] molecules), whereas the second harnesses an antibody-based format (TCR-mimic antibodies). For both classes of reagent, target specificity is vital, considering the vast universe of potential pHLA molecules that can be presented on healthy cells. Here, we made use of structural, biochemical, and computational approaches to investigate the molecular rules underpinning the reactivity patterns of pHLA-targeting bispecifics. We demonstrate that affinity-enhanced TCRs engage pHLA using a comparatively broad and balanced energetic footprint, with interactions distributed over several HLA and peptide side chains. As ImmTAC molecules, these TCRs also retained a greater degree of pHLA selectivity, with less off-target activity in cellular assays. Conversely, TCR-mimic antibodies tended to exhibit binding modes focused more toward hot spots on the HLA surface and exhibited a greater degree of crossreactivity. Our findings extend our understanding of the basic principles that underpin pHLA selectivity and exemplify a number of molecular approaches that can be used to probe the specificity of pHLA-targeting molecules, aiding the development of future reagents.

Peptide cargo tunes a network of correlated motions in human leucocyte antigens.

Most biomolecular interactions are typically thought to increase the (local) rigidity of a complex, for example, in drug-target binding. However, detailed analysis of specific biomolecular complexes can reveal a more subtle interplay between binding and rigidity. Here, we focussed on the human leucocyte antigen (HLA), which plays a crucial role in the adaptive immune system by presenting peptides for recognition by the αß T-cell receptor (TCR). The role that the peptide plays in tuning HLA flexibility during TCR recognition is potentially crucial in determining the functional outcome of an immune response, with obvious relevance to the growing list of immunotherapies that target the T-cell compartment. We have applied high-pressure/temperature perturbation experiments, combined with molecular dynamics simulations, to explore the drivers that affect molecular flexibility for a series of different peptide-HLA complexes. We find that different peptide sequences affect peptide-HLA flexibility in different ways, with the peptide cargo tuning a network of correlated motions throughout the pHLA complex, including in areas remote from the peptide-binding interface, in a manner that could influence T-cell antigen discrimination.

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.

A Multiscale Simulation Approach to Modeling Drug-Protein Binding Kinetics.

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.

Conformational change and ligand binding in the aristolochene synthase catalytic cycle.

Terpene synthases are potentially useful biocatalysts for the synthesis of valuable compounds, such as anticancer drugs and antibiotics. The design of altered activities requires better knowledge of their mechanisms, for example, an understanding of the complex conformational changes that are part of their catalytic cycle, how they are coordinated, and what drives them. Crystallographic studies of the sesquiterpene synthase artistolochene synthase have led to a proposed sequence of ligand binding and conformational change but have provided only indirect insight. Here, we have performed extensive molecular dynamics simulations of multiple enzyme-ligand complexes (over 2 µs in total). The simulations provide clear evidence of what drives the conformational changes required for reaction. They support a picture in which the substrate farnesyl diphosphate binds first, followed by three magnesium ions in sequence, and, after reaction, the release of aristolochene and two magnesium ions followed by the final magnesium ion and diphosphate. Binding of farnesyl diphosphate leads to an increased level of sampling of open conformations, allowing the first two magnesium ions to bind. The closed enzyme conformation is maintained with a diphosphate moiety and two magnesium ions bound. The open-to-closed transition reduces flexibility around the active site entrance, partly through a lid closing over it. The simulations with all three magnesium ions and farnesyl diphosphate bound provide, for the first time, a realistic model of the Michaelis complex involved in reaction, which is inaccessible to experimental structural studies. These insights could help with the design of altered activities in a range of terpene synthases.

Affinity of Avr2 for tomato cysteine protease Rcr3 correlates with the Avr2-triggered Cf-2-mediated hypersensitive response.

The Cladosporium fulvum Avr2 effector is a novel type of cysteine protease inhibitor with eight cysteine residues that are all involved in disulphide bonds. We have produced wild-type Avr2 protein in Pichia pastoris and determined its disulphide bond pattern. By site-directed mutagenesis of all eight cysteine residues, we show that three of the four disulphide bonds are required for Avr2 stability. The six C-terminal amino acid residues of Avr2 contain one disulphide bond that is not embedded in its overall structure. Avr2 is not processed by the tomato cysteine protease Rcr3 and is an uncompetitive inhibitor of Rcr3. We also produced mutant Avr2 proteins in which selected amino acid residues were individually replaced by alanine, and, in one mutant, all six C-terminal amino acid residues were deleted. We determined the inhibitory constant (K(i) ) of these mutants for Rcr3 and their ability to trigger a Cf-2-mediated hypersensitive response (HR) in tomato. We found that the two C-terminal cysteine residues and the six amino acid C-terminal tail of Avr2 are required for both Rcr3 inhibitory activity and the ability to trigger a Cf-2-mediated HR. Individual replacement of the lysine-17, lysine-20 or tyrosine-21 residue by alanine did not affect significantly the biological activity of Avr2. Overall, our data suggest that the affinity of the Avr2 mutants for Rcr3 correlates with their ability to trigger a Cf-2-mediated HR.