High-resolution Crystal Structures Research Articles

A Comparative Study for the Accuracy of Three Molecular Docking Programs Using HIV-1 Protease Inhibitors as a Model

Flexible molecular docking is a computational method of structure-based drug design to evaluate binding interactions between receptor and ligand and identify the ligand conformation within the receptor pocket. Currently, various molecular docking programs are extensively applied; therefore, realizing accuracy and performance of the various docking programs could have a significant value. In this comparative study, the performance and accuracy of three widely used non-commercial docking software (AutoDock Vina, 1-Click Docking, and UCSF DOCK) was evaluated through investigations of the predicted binding affinity and binding conformation of the same set of small molecules (HIV-1 protease inhibitors) and a protein target HIV-1 protease enzyme. The tested sets are composed of eight receptor-ligand complexes with high resolution crystal structures downloaded from Protein Data Bank website. Molecular dockings were applied between approved HIV-1 protease inhibitors and the HIV-1 protease using AutoDock Vina, 1-Click Docking, and DOCK6. Then, docking poses of the top-ranked solution was realized using UCSF Chimera. Furthermore, Pearson correlation coefficient (r) and coefficient of determination (r2) between the experimental results and the top scored docking results of each program were calculated using Graphpad prism V9.2. After comparing saquinavir top scored binding poses of each docking program with the crystal structure, various conformational changes were observed. Moreover, according to the relative comparison between the top ranked calculated ?Gbinding values against the experimental results, r2 value of AutoDock Vina, 1-Click Docking, and DOCK6 were 0.65, 0.41, and 0.005, respectively. The outcome of this study shows that the top scored binding free energy could not produce the best pose prediction. In addition, AutoDock Vina results have the highest correlation with the experimental results.

Insights into the molecular mechanismof translation inhibition by the ribosome-targeting antibiotic thermorubin.

Thermorubin (THR) is an aromatic anthracenopyranone antibiotic active against both Gram-positive and Gram-negative bacteria. It is known to bind to the 70S ribosome at the intersubunit bridge B2a and was thought to inhibit factor-dependent initiation of translation and obstruct the accommodation of tRNAs into the A site. Here, we show that thermorubin causes ribosomes to stall in vivo and in vitro at internal and termination codons, thereby allowing the ribosome to initiate protein synthesis and translate at least a few codons before stalling. Our biochemical data show that THR affects multiple steps of translation elongation with a significant impact on the binding stability of the tRNA in the A site, explaining premature cessation of translation. Our high-resolution crystal and cryo-EM structures of the 70S-THR complex show that THR can co-exist with P- and A-site tRNAs, explaining how ribosomes can elongate in the presence of the drug. Remarkable is the ability of THR to arrest ribosomes at the stop codons. Our data suggest that by causing structural re-arrangements in the decoding center, THR interferes with the accommodation of tRNAs or release factors into the ribosomal A site.

Effects of Active-Center Reduction of Plant-Type Ferredoxin on Its Structure and Dynamics: Computational Analysis Using Molecular Dynamics Simulations.

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1.

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

Calmodulin variant E140G associated with long QT syndrome impairs CaMKIIδ autophosphorylation and L-type calcium channel inactivation

Long QT syndrome (LQTS) is a human inherited heart condition that can cause life-threatening arrhythmia including sudden cardiac death. Mutations in the ubiquitous Ca2+-sensingprotein calmodulin (CaM) are associated with LQTS, but the molecular mechanism by which these mutations lead to irregular heartbeats is not fully understood. Here, we use a multidisciplinary approach including protein biophysics, structural biology, confocal imaging, and patch-clamp electrophysiology to determine the effect of the disease-associated CaM mutation E140G on CaM structure and function. We present novel data showing that mutant-regulated CaMKIIδ kinase activity is impaired with a significant reduction in enzyme autophosphorylation rate. We report the first high-resolution crystal structure of a LQTS-associated CaM variant in complex with the CaMKIIδ peptide, which shows significant structural differences, compared to the WT complex. Furthermore, we demonstrate that the E140G mutation significantly disrupted Cav1.2 Ca2+/CaM-dependent inactivation, while cardiac ryanodine receptor (RyR2) activity remained unaffected. In addition, we show that the LQTS-associated mutation alters CaM's Ca2+-binding characteristics, secondary structure content, and interaction with key partners involved in excitation-contraction coupling (CaMKIIδ, Cav1.2, RyR2). In conclusion, LQTS-associated CaM mutation E140G severely impacts the structure-function relationship of CaM and its regulation of CaMKIIδ and Cav1.2. This provides a crucial insight into the molecular factors contributing to CaM-mediated arrhythmias with a central role for CaMKIIδ.

Tenuous transcriptional threshold of human sex determination. II. SRY exploits water-mediated clamp at the edge of ambiguity.

Y-encoded transcription factor SRY initiates male differentiation in therian mammals. This factor contains a high-mobility-group (HMG) box, which mediates sequence-specific DNA binding with sharp DNA bending. A companion article in this issue described sex-reversal mutations at box position 72 (residue 127 in human SRY), invariant as Tyr among mammalian orthologs. Although not contacting DNA, the aromatic ring seals the domain's minor wing at a solvent-exposed junction with a basic tail. A seeming paradox was posed by the native-like biochemical properties of inherited Swyer variant Y72F: its near-native gene-regulatory activity is consistent with the father's male development, but at odds with the daughter's XY female somatic phenotype. Surprisingly, aromatic rings (Y72, F72 or W72) confer higher transcriptional activity than do basic or polar side chains generally observed at solvated DNA interfaces (Arg, Lys, His or Gln). Whereas biophysical studies (time-resolved fluorescence resonance energy transfer and heteronuclear NMR spectroscopy) uncovered only subtle perturbations, dissociation of the Y72F complex was markedly accelerated relative to wild-type. Studies of protein-DNA solvation by molecular-dynamics (MD) simulations of an homologous high-resolution crystal structure (SOX18) suggest that Y72 para-OH anchors a network of water molecules at the tail-DNA interface, perturbed in the variant in association with nonlocal conformational fluctuations. Loss of the Y72 anchor among SRY variants presumably "unclamps" its basic tail, leading to (a) rapid DNA dissociation despite native affinity and (b) attenuated transcriptional activity at the edge of sexual ambiguity. Conservation of Y72 suggests that this water-mediated clamp operates generally among SRY and metazoan SOX domains.

Discovery, Structure, and Mechanism of a Class II Sesquiterpene Cyclase.

Terpene cyclases (TCs), extraordinary enzymes that create the structural diversity seen in terpene natural products, are traditionally divided into two classes, class I and class II. Although the structural and mechanistic features of class I TCs are well-known, the corresponding details in class II counterparts have not been fully characterized. Here, we report the genome mining discovery and structural characterization of two class II sesquiterpene cyclases (STCs) from Streptomyces. These drimenyl diphosphate synthases (DMSs) are the first STCs shown to possess β,γ-didomain architecture. High-resolution X-ray crystal structures of DMS from Streptomyces showdoensis (SsDMS) in complex with both a farnesyl diphosphate and Mg2+ unveiled an induced-fit mechanism, with an unprecedented Mg2+ binding mode, finally solving one of the lingering questions in class II TC enzymology. This study supports continued genome mining for novel bacterial TCs and provides new mechanistic insights into canonical class II TCs that will lead to advances in TC engineering and synthetic biology.

A backbone-dependent rotamer library with high (ϕ, ψ) coverage using metadynamics simulations.

Backbone-dependent rotamer libraries are commonly used to assign the side chain dihedral angles of amino acids when modeling protein structures. Most rotamer libraries are created by curating protein crystal structure data and using various methods to extrapolate the existing data to cover all possible backbone conformations. However, these rotamer libraries may not be suitable for modeling the structures of cyclic peptides and other constrained peptides because these molecules frequently sample backbone conformations rarely seen in the crystal structures of linear proteins. To provide backbone-dependent side chain information beyond the α-helix, β-sheet, and PPII regions, we used explicit-solvent metadynamics simulations of model dipeptides to create a new rotamer library that has high coverage in the (ϕ, ψ) space. Furthermore, this approach can be applied to build high-coverage rotamer libraries for noncanonical amino acids. The resulting Metadynamics of Dipeptides for Rotamer Distribution (MEDFORD) rotamer library predicts the side chain conformations of high-resolution protein crystal structures with similar accuracy (~80%) to a state-of-the-art rotamer library. Our ability to test the accuracy of MEDFORD at predicting the side chain dihedral angles of amino acids in noncanonical backbone conformation is restricted by the limited structural data available for cyclic peptides. For the cyclic peptide data that are currently available, MEDFORD and the state-of-the-art rotamer library perform comparably. However, the two rotamer libraries indeed make different rotamer predictions in noncanonical (ϕ, ψ) regions. For noncanonical amino acids, the MEDFORD rotamer library predicts the χ1 values with approximately 75% accuracy.

Structural basis for the carotenoid binding and transport function of a START domain

STARD3, a steroidogenic acute regulatory lipid transfer protein, was identified as a key xanthophyll-binding protein in the human retina. STARD3 and its homologs in invertebrates are known to bind and transport carotenoids, but this lacks structural elucidation. Here, we report high-resolution crystal structures of the apo- and zeaxanthin (ZEA)-bound carotenoid-binding protein from silkworm Bombyx mori (BmCBP). Having a STARD3-like fold, BmCBP features novel elements, including the Ω1-loop that, in the apoform, is uniquely fixed on the α4-helix by an R173-D279 salt bridge. We exploit absorbance, Raman and dichroism spectroscopy, and calorimetry to describe how ZEA and BmCBP mutually affect each other in the complex. We identify key carotenoid-binding residues, confirm their roles by ZEA-binding capacity and X-ray structures of BmCBP mutants, and also demonstrate that markedly different carotenoid-binding capacities of BmCBP and human STARD3 stem from differences in the structural organization of their carotenoid-binding cavity.

Interaction between DLC-1 and SAO-1 facilitates CED-4 translocation during apoptosis in the Caenorhabditis elegans germline

Apoptosis is one of the major forms of programmed cell death, and it serves vital biological functions in multicellular animal and plant cells. The core mechanism of apoptosis is highly conserved in metazoans, where the translocation of CED-4/Apaf-1 from mitochondria to the nuclear membrane is required to initiate and execute apoptosis. However, the underlying molecular mechanisms of this translocation are poorly understood. In this study, we showed that SAO-1 binds DLC-1 and prevents its degradation to promote apoptosis in C. elegans germ cells. We demonstrated that SAO-1 and DLC-1 regulate CED-4/Apaf-1 nuclear membrane accumulation during apoptosis. Isothermal titration calorimetry-based assay and high-resolution crystal structure analysis further revealed that SAO-1 interacted with DLC-1 to form a 2:4 complex: each of the two β-sheets in the SAO-1 peptide interacted with two DLC-1 dimers. Point mutations at the SAO-1-DLC-1 binding interface significantly inhibited apoptotic corpse formation and CED-4 nuclear membrane accumulation within C. elegans germ cells. In conclusion, our study provides a new perspective on the regulation of CED-4-mediated apoptosis.

High-Resolution Crystal Structure of Muscle Phosphoglycerate Mutase Provides Insight into Its Nuclear Import and Role.

Phosphoglycerate mutase (PGAM) is a glycolytic enzyme converting 3-phosphoglycerate to 2-phosphoglycerate, which in mammalian cells is expressed in two isoforms: brain (PGAM1) and muscle (PGAM2). Recently, it was shown that besides its enzymatic function, PGAM2 can be imported to the cell nucleus where it co-localizes with the nucleoli. It was suggested that it functions there to stabilize the nucleolar structure, maintain mRNA expression, and assist in the assembly of new pre-ribosomal subunits. However, the precise mechanism by which the protein translocates to the nucleus is unknown. In this study, we present the first crystal structure of PGAM2, identify the residues involved in the nuclear localization of the protein and propose that PGAM contains a "quaternary nuclear localization sequence (NLS)", i.e., one that consists of residues from different protein chains. Additionally, we identify potential interaction partners for PGAM2 in the nucleoli and demonstrate that 14-3-3ζ/δ is indeed an interaction partner of PGAM2 in the nucleus. We also present evidence that the insulin/IGF1-PI3K-Akt-mTOR signaling pathway is responsible for the nuclear localization of PGAM2.

Identification of New L-Heptanoylphosphatidyl Inositol Pentakisphosphate Derivatives Targeting the Interaction with HIV-1 Gag by Molecular Modelling Studies.

The HIV-1 Gag protein binds to the host cell membrane and assembles into immature particles. Then, in the course of immature virion budding, activated protease cleaves Gag into its main components: MA, CA, NC, and p6 proteins. The highly basic residues of MA predominantly interact with the acidic head of phosphatidyl-inositol-4,5-bisphosphate (PI(4,5)P2) inserted into the membrane. Our research group developed L-Heptanoylphosphatidyl Inositol Pentakisphosphate (L-HIPPO) and previously confirmed that this compound bound to the MA more strongly than PI(4,5)P2 and inositol hexakisphosphate (IP6) did. Therefore, herein we rationally designed eight new L-HIPPO derivatives based on the fact that the most changeable parts of L-HIPPO were two acyl chains. After that, we employed molecular docking for eight compounds via Maestro software using high-resolution crystal structures of MA in complex with IP6 (PDB IDs: 7E1I, 7E1J, and 7E1K), which were recently elucidated by our research group. The most promising docking scores were obtained with benzene-inserted compounds. Thus, we generated a library containing 213 new aromatic group-inserted L-HIPPO derivatives and performed the same molecular docking procedure. According to the results, we determined the nine new L-HIPPO derivatives most effectively binding to the MA with the most favorable scoring functions and pharmacokinetic properties for further exploration.

Discovery, evaluation and mechanism study of WDR5-targeted small molecular inhibitors for neuroblastoma.

Neuroblastoma is the most common and deadliest tumor in infancy. WDR5 (WD Repeat Domain 5), a critical factor supporting an N-myc transcriptional complex via its WBM site and interacting with chromosome via its WIN site, promotes the progression of neuroblastoma, thus making it a potential anti-neuroblastoma drug target. So far, a few WIN site inhibitors have been reported, and the WBM site disruptors are rare to see. In this study we conducted virtual screening to identify candidate hit compounds targeting the WBM site of WDR5. As a result, 60 compounds were selected as candidate WBM site inhibitors. Cell proliferation assay demonstrated 6 structurally distinct WBM site inhibitors, numbering as compounds 4, 7, 11, 13, 19 and 22, which potently suppressed 3 neuroblastoma cell lines (MYCN-amplified IMR32 and LAN5 cell lines, and MYCN-unamplified SK-N-AS cell line). Among them, compound 19 suppressed the proliferation of IMR32 and LAN5 cells with EC50 values of 12.34 and 14.89 μM, respectively, and exerted a moderate inhibition on SK-N-AS cells, without affecting HEK293T cells at 20 μM. Analysis of high-resolution crystal complex structure of compound 19 against WDR5 revealed that it competitively occupied the hydrophobic pocket where V264 was located, which might disrupt the interaction of MYC with WDR5 and further MYC-medicated gene transcription. By performing RNA-seq analysis we demonstrated the differences in molecular action mechanisms of the compound 19 and a WIN site inhibitor OICR-9429. Most interestingly, we established the particularly high synergy rate by combining WBM site inhibitor 19 and the WIN site inhibitor OICR-9429, providing a novel therapeutic avenue for neuroblastoma.

LY3041658/ interleukin-8 complex structure as targets for IL-8 small molecule inhibitors discovery using a combination of in silico methods

ABSTRACT Since interleukin-8 (IL-8/CXCL8) and its receptor, CXCR1 and CXCR2, were known in the early 1990s, biological pathways related to these proteins were proven to have high clinical value in cancer and inflammatory/autoimmune conditions treatment. Recently, IL-8 has been identified as biomarker for severe COVID-19 patients and COVID-19 prognosis. Boyles et al. (mAbs 12 (2020), pp. 1831880) have published a high-resolution X-ray crystal structure of the LY3041658 Fab in a complex human CXCL8. They described the ability to bind to IL-8 and the blocking of IL-8/its receptors interaction by the LY3041658 monoclonal antibody. Therefore, the study has been designed to identify potential small molecules inhibiting interleukin-8 by targeting LY3041658/IL-8 complex structure using an in silico approach. A structure‑based pharmacophore and molecular docking models of the protein active site cavity were generated to identify possible candidates, followed by virtual screening with the ZINC database. ADME analysis of hit compounds was also conducted. Molecular dynamics simulations were then performed to survey the behaviour and stability of the ligand-protein complexes. Furthermore, the MM/PBSA technique has been utilized to evaluate the free binding energy. The final data confirmed that one newly obtained compound, ZINC21882765, may serve as the best potential inhibitor for IL-8.

Viruses inhibit TIR gcADPR signalling to overcome bacterial defence.

The Toll/interleukin-1 receptor (TIR) domain is a key component of immune receptors that identify pathogen invasion in bacteria, plants and animals1-3. In the bacterial antiphage system Thoeris, as well as in plants, recognition of infection stimulates TIR domains to produce an immune signalling molecule whose molecular structure remains elusive. This molecule binds and activates the Thoeris immune effector, which then executes the immune function1. We identified a large family of phage-encoded proteins, denoted here as Thoeris anti-defence 1 (Tad1), that inhibit Thoeris immunity. We found that Tad1 proteins are 'sponges' that bind and sequester the immune signalling molecule produced by TIR-domain proteins, thus decoupling phage sensing from immune effector activation and rendering Thoeris inactive. Tad1 can also efficiently sequester molecules derived from a plant TIR-domain protein, and a high-resolution crystal structure of Tad1 bound to a plant-derived molecule showed a unique chemical structure of 1 ''-2' glycocyclic ADPR (gcADPR). Our data furthermore suggest that Thoeris TIR proteins produce a closely related molecule, 1''-3' gcADPR, which activates ThsA an order of magnitude more efficiently than the plant-derived 1''-2' gcADPR. Our results define the chemical structure of a central immune signalling molecule and show a new mode of action by which pathogens can suppress host immunity.

Structure of the SHOC2-MRAS-PP1C complex provides insights into RAF activation and Noonan syndrome.

SHOC2 acts as a strong synthetic lethal interactor with MEK inhibitors in multiple KRAS cancer cell lines. SHOC2 forms a heterotrimeric complex with MRAS and PP1C that is essential for regulating RAF and MAPK-pathway activation by dephosphorylating a specific phosphoserine on RAF kinases. Here we present the high-resolution crystal structure of the SHOC2-MRAS-PP1C (SMP) complex and apo-SHOC2. Our structures reveal that SHOC2, MRAS, and PP1C form a stable ternary complex in which all three proteins synergistically interact with each other. Our results show that dephosphorylation of RAF substrates by PP1C is enhanced upon interacting with SHOC2 and MRAS. The SMP complex forms only when MRAS is in an active state and is dependent on SHOC2 functioning as a scaffolding protein in the complex by bringing PP1C and MRAS together. Our results provide structural insights into the role of the SMP complex in RAF activation and how mutations found in Noonan syndrome enhance complex formation, and reveal new avenues for therapeutic interventions.

Comparative structural analyses of the NHL domains from the human E3 ligase TRIM-NHL family.

Tripartite motif (TRIM) proteins constitute one of the largest subfamilies of the RING-type E3 ubiquitin ligases that play a role in diverse processes from homeostasis and immune response to viral restriction. While TRIM proteins typically harbor an N-terminal RING finger, a B-box and a coiled-coil domain, a high degree of diversity lies in their C termini that contain diverse protein interaction modules, most of which, both structures and their roles in intermolecular interactions, remain unknown. Here, high-resolution crystal structures of the NHL domains of three of the four human TRIM-NHL proteins, namely TRIM2, TRIM3 and TRIM71, are presented. Comparative structural analyses revealed that, despite sharing an evolutionarily conserved six-bladed β-propeller architecture, the low sequence identities resulted in distinct properties of these interaction domains at their putative binding sites for macromolecules. Interestingly, residues lining the binding cavities represent a hotspot for genetic mutations linked to several diseases. Thus, high sequence diversity within the conserved NHL domains might be essential for differentiating binding partners among TRIM-NHL proteins.

High-resolution crystal structure and chemical screening reveal pantothenate kinase as a new target for antifungal development

Fungal infections are the leading cause of mortality by eukaryotic pathogens, with an estimated 150 million severe life-threatening cases and 1.7 million deaths reported annually. The rapid emergence of multidrug-resistant fungal isolates highlights the urgent need for new drugs with new mechanisms of action. In fungi, pantothenate phosphorylation, catalyzed by PanK enzyme, is the first step in the utilization of pantothenic acid and coenzyme A biosynthesis. In all fungi sequenced so far, this enzyme is encoded by a single PanK gene. Here, we report the crystal structure of a fungal PanK alone as well as with high-affinity inhibitors from a single chemotype identified through a high-throughput chemical screen. Structural, biochemical, and functional analyses revealed mechanisms governing substrate and ligand binding, dimerization, and catalysis and helped identify new compounds that inhibit the growth of several Candida species. The data validate PanK as a promising target for antifungal drug development.

Dioxygen Activation and Nδ,Nε-Dihydroxylation Mechanism Involved in the Formation of N-Nitrosourea Pharmacophore in Streptozotocin Catalyzed by Nonheme Diiron Enzyme SznF.

SznF is a nonheme diiron-dependent enzyme that catalyzes the critical N-nitrosation involved in the formation of the N-nitrosourea moiety in the pancreatic cancer drug streptozotocin. The N-nitrosation contains two successive N-hydroxylation and N-nitrosation steps, which are carried out by two separate active sites, namely, the central domain and cupin domain. Recently, the crystal structure of SznF was obtained, and the central domain was proved to contain a diiron cofactor to catalyze the N-hydroxylation. In this work, to gain insights into the O2 activation and the successive N-hydroxylation mechanism, on the basis of the high-resolution crystal structure, the enzyme-substrate complex models were constructed, and a series of combined QM/MM calculations were performed. Based on our calculations, the activation of O2 starts from the diiron(II,III)-superoxo (S) to generate the diiron(IV)-oxo species (Q) via a diiron(III,III)-peroxo (P)-like transition state or unstable intermediate (P'), and species P' can be described as a hybridization of diiron(IV)-oxo species and diiron(III,III)-peroxo (P) owing to the long distances of Fe1-Fe2 (4.22 Å) and O1-O2 (1.89 Å), which is different from those of other nonheme diiron enzymes. In the following hydroxylation of Nδ and Nε, the Nδ-hydroxylation was confirmed to occur first, agreeing with the experimental observations. Because the diiron(IV)-oxo species (Q) is responsible for hydroxylation, the reaction follows the H-abstraction/OH rebound mechanism, and the first abstraction occurs on the Nδ-H rather than Nε-H, which may be attributed to the different orientation of Fe(IV)-oxo relative to N-H as well as the bond dissociation enthalpies of two N-H bonds. The hydroxylation of N-methyl-L-arginine does not employ the diiron(III,III)-hydroperoxo (P″) to trigger the electrophilic attack of the guanidine to directly form the N-O bond, as previously suggested. In addition, our calculations also revealed that the direct attack of the Fe(IV)═O unit to the Nδ of the substrate corresponds to a higher barrier than that in the H-abstraction/OH rebound mechanism. These results may provide useful information for understanding the formation of the di-hydroxylation intermediate involved in the biosynthesis of N-nitrosation.

How a 10-epi-Cubebol Synthase Avoids Premature Reaction Quenching to Form a Tricyclic Product at High Purity.

Terpenes are the largest class of natural products and are attractive targets in the fuel, fragrance, pharmaceutical, and flavor industries. Harvesting terpenes from natural sources is environmentally intensive and often gives low yields and purities, requiring further downstream processing. Engineered terpene synthases (TSs) offer a solution to these problems, but the low sequence identity and high promiscuity among TSs are major challenges for targeted engineering. Rational design of TSs requires identification of key structural and chemical motifs that steer product outcomes. Producing the sesquiterpenoid 10-epi-cubebol from farnesyl pyrophosphate (FPP) requires many steps and some of Nature’s most difficult chemistry. 10-epi-Cubebol synthase from Sorangium cellulosum (ScCubS) guides a highly reactive carbocationic substrate through this pathway, preventing early quenching and ensuring correct stereochemistry at every stage. The cyclizations carried out by ScCubS potentially represent significant evolutionary expansions in the chemical space accessible by TSs. Here, we present the high-resolution crystal structure of ScCubS in complex with both a trinuclear magnesium cluster and pyrophosphate. Computational modeling, experiment, and bioinformatic analysis identified residues important in steering the reaction chemistry. We show that S206 is crucial in 10-epi-cubebol synthesis by enlisting the nearby F211 to shape the active site contour and prevent the formation of early escape cadalane products. We also show that N327 and F104 control the distribution between several early-stage cations and whether the final product is derived from the germacrane, cadalane, or cubebane hydrocarbon scaffold. Using these insights, we reengineered ScCubS so that its main product was germacradien-4-ol, which derives from the germacrane, rather than the cubebane, scaffold. Our work emphasizes that mechanistic understanding of cation stabilization in TSs can be used to guide catalytic outcomes.