A kiwellin disarms the metabolic activity of a secreted fungal virulence factor.

Fungi-induced plant diseases affect global food security and plant ecology. The biotrophic fungus Ustilago maydis causes smut disease in maize (Zea mays) plants by secreting numerous virulence effectors that reprogram plant metabolism and immune responses1,2. The secreted fungal chorismate mutase Cmu1 presumably affects biosynthesis of the plant immune signal salicylic acid by channelling chorismate into the phenylpropanoid pathway3. Here we show that one of the 20 maize-encoded kiwellins (ZmKWL1) specifically blocks the catalytic activity of Cmu1. ZmKWL1 hinders substrate access to the active site of Cmu1 through intimate interactions involving structural features that are specific to fungal Cmu1 orthologues. Phylogenetic analysis suggests that plant kiwellins have a versatile scaffold that can specifically counteract pathogen effectors such as Cmu1. We reveal the biological activity of a member of the kiwellin family, a widely conserved group of proteins that have previously been recognized only as important human allergens.

Metabolic engineering of Synechocystis sp. PCC 6803 for the improved production of phenylpropanoids.

BACKGROUND: Phenylpropanoids are a large group of plant secondary metabolites with various biological functions, derived from aromatic amino acids. Cyanobacteria are promising host organisms for sustainable production of plant phenylpropanoids. We have previously engineered Synechocystis sp. PCC 6803 to produce trans-cinnamic acid (tCA) and p-coumaric acid (pCou), the first intermediates of phenylpropanoid pathway, by overexpression of phenylalanine- and tyrosine ammonia lyases. In this study, we aimed to enhance the production of the target compounds tCA and pCou in Synechocystis. RESULTS: We eliminated the 4-hydroxyphenylpyruvate dioxygenase (HPPD) activity, which is a competing pathway consuming tyrosine and, possibly, phenylalanine for tocopherol synthesis. Moreover, several genes of the terminal steps of the shikimate pathway were overexpressed alone or in operons, such as aromatic transaminases, feedback insensitive cyclohexadienyl dehydrogenase (TyrC) from Zymomonas mobilis and the chorismate mutase (CM) domain of the fused chorismate mutase/prephenate dehydratase enzyme from Escherichia coli. The obtained engineered strains demonstrated nearly 1.5 times enhanced tCA and pCou production when HPPD was knocked out compared to the parental production strains, accumulating 138 ± 3.5 mg L-1 of tCA and 72.3 ± 10.3 mg L-1 of pCou after seven days of photoautotrophic growth. However, there was no further improvement when any of the pathway genes were overexpressed. Finally, we used previously obtained AtPRM8 and TsPRM8 Synechocystis strains with deregulated shikimate pathway as a background for the overexpression of synthetic constructs with ppd knockout. CONCLUSIONS: HPPD elimination enhances the tCA and pCou productivity to a similar extent. The use of PRM8 based strains as a background for overexpression of synthetic constructs, however, did not promote tCA and pCou titers, which indicates a tight regulation of the terminal steps of phenylalanine and tyrosine synthesis. This work contributes to establishing cyanobacteria as hosts for phenylpropanoid production.

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study.

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

What Drives Chorismate Mutase to Top Performance? Insights from a Combined In Silico and In Vitro Study.

Unlike typical chorismate mutases, the enzyme from Mycobacterium tuberculosis (MtCM) has only low activity on its own. Remarkably, its catalytic efficiency kcat/Km can be boosted more than 100-fold by complex formation with a partner enzyme. Recently, an autonomously fully active MtCM variant was generated using directed evolution, and its structure was solved by X-ray crystallography. However, key residues were involved in crystal contacts, challenging the functional interpretation of the structural changes. Here, we address these challenges by microsecond molecular dynamics simulations, followed up by additional kinetic and structural analyses of selected sets of specifically engineered enzyme variants. A comparison of wild-type MtCM with naturally and artificially activated MtCMs revealed the overall dynamic profiles of these enzymes as well as key interactions between the C-terminus and the active site loop. In the artificially evolved variant of this model enzyme, this loop is preorganized and stabilized by Pro52 and Asp55, two highly conserved residues in typical, highly active chorismate mutases. Asp55 stretches across the active site and helps to appropriately position active site residues Arg18 and Arg46 for catalysis. The role of Asp55 can be taken over by another acidic residue, if introduced at position 88 close to the C-terminus of MtCM, as suggested by molecular dynamics simulations and confirmed by kinetic investigations of engineered variants.

Mechanism of chorismate dehydratase MqnA, the first enzyme of the futalosine pathway, proceeds via substrate-assisted catalysis.

MqnA, the only chorismate dehydratase known so far, catalyzes the initial step in the biosynthesis of menaquinone via the futalosine pathway. Details of the MqnA reaction mechanism remain unclear. Here, we present crystal structures of Streptomyces coelicolor MqnA and its active site mutants in complex with chorismate and the product 3-enolpyruvyl-benzoate, produced during heterologous expression in Escherichia coli. Together with activity studies, our data are in line with dehydration proceeding via substrate assisted catalysis, with the enol pyruvyl group of chorismate acting as catalytic base. Surprisingly, structures of the mutant Asn17Asp with copurified ligand suggest that the enzyme converts to a hydrolase by serendipitous positioning of the carboxyl group. All complex structures presented here exhibit a closed Venus flytrap fold, with the enzyme exploiting the characteristic ligand binding properties of the fold for specific substrate binding and catalysis. The conformational rearrangements that facilitate complete burial of substrate/product, with accompanying topological changes to the enzyme surface, could foster substrate channeling within the biosynthetic pathway.

A secreted form of chorismate mutase (Rv1885c) in Mycobacterium bovis BCG contributes to pathogenesis by inhibiting mitochondria-mediated apoptotic cell death of macrophages.

BACKGROUND: Mycobacterium tuberculosis is the causative agent of tuberculosis (TB), and its pathogenicity is associated with its ability to evade the host defense system. The secretory form of the chorismate mutase of M. tuberculosis (TBCM, encoded by Rv1885c) is assumed to play a key role in the pathogenesis of TB; however, the mechanism remains unknown. METHODS: A tbcm deletion mutant (B∆tbcm) was generated by targeted gene knockout in BCG to investigate the pathogenic role of TBCM in mice or macrophages. We compared the pathogenesis of B∆tbcm and wild-type BCG in vivo by measuring the bacterial clearance rate and the degree of apoptosis. Promotion of the intrinsic apoptotic pathway was evaluated in infected bone marrow-derived macrophages (BMDMs) by measuring apoptotic cell death, loss of mitochondrial membrane potential and translocation of pore-forming proteins. Immunocytochemistry, western blotting and real-time PCR were also performed to assess the related protein expression levels after infection. Furthermore, these findings were validated by complementation of tbcm in BCG. RESULTS: Deletion of the tbcm gene in BCG leads to reduced pathogenesis in a mouse model, compared to wild type BCG, by promoting apoptotic cell death and bacterial clearance. Based on these findings, we found that intrinsic apoptosis and mitochondrial impairment were promoted in B∆tbcm-infected BMDMs. B∆tbcm down-regulates the expression of Bcl-2, which leads to mitochondrial outer membrane permeabilization (MOMP), culminating in cytochrome c release from mitochondria. Consistent with this, transcriptome profiling also indicated that B∆tbcm infection is more closely related to altered mitochondrial-related gene expression than wild-type BCG infection, suggesting an inhibitory role of TBCM in mitochondrial dysfunction. Moreover, genetic complementation of B∆tbcm (C∆tbcm) restored its capacity to inhibit mitochondria-mediated apoptotic cell death. CONCLUSIONS: Our findings demonstrate the contribution of TBCM to bacterial survival, inhibiting intrinsic apoptotic cell death of macrophages as a virulence factor of M. tuberculosis complex (MTBC) strains, which could be a potential target for the development of TB therapy.

Best Practices on QM/MM Simulations of Biological Systems.

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

Wang resin catalysed sonochemical synthesis of pyrazolo[4,3-d]pyrimidinones and 2,3-dihydroquinazolin-4(1H)-ones: Identification of chorismate mutase inhibitors having effects on Mycobacterium tuberculosis cell viability.

The enzyme chorismate mutase (or CM that is vital for the survival of bacteria) is an interesting pharmacological target for the identification of new anti-tubercular agents. The 5,5-disibstituted pyrazolo[4,3-d]pyrimidinone derivatives containing the fragment based on 4-amino-1-methyl-3-propyl-1H-pyrazole-5-carboxamide were designed and explored as the potential inhibitors of chorismate mutase. Based on encouraging docking results of two representative molecules evaluated in silico against MtbCM (PDB: 2FP2) the Wang resin catalysed sonochemical synthesis of target N-heteroarenes were undertaken. The methodology involved the reaction of 4-amino-1-methyl-3-propyl-1H-pyrazole-5-carboxamide with the appropriate cyclic/acyclic ketones to afford the desired products in acceptable (51-94%) yields. The methodology was also extended successfully towards the synthesis of 2,2-disubstituted 2,3-dihydroquinazolin-4(1H)-ones in excellent (85-90%) yields. In vitro MTT assay against the RAW 264.7 cell line followed by enzymatic assay against MtbCM identified 3b and 3c as active compounds that showed two H-bonding via their NH (at position 6) and CO group with MtbCM in silico and encouraging (54-57%) inhibition at 30 µM in vitro. Notably, none of the 2,2-disubstituted 2,3-dihydroquinazolin-4(1H)-ones showed any significant inhibition of MtbCM suggesting the favourable role of the pyrazole moiety in case of pyrazolo[4,3-d]pyrimidinones. The favourable role of cyclopentyl ring attached to the pyrazolo[4,3-d]pyrimidinone moiety and that of two methyl groups in place of cyclopentyl ring was also indicated by the SAR study. Besides showing effects against MtbCM in the concentration response study, 3b and 3c showed little or no effects on mammalian cell viability up to 100 µM in an MTT assay but decreased the % Mtb cell viability at 10-30 µM with > 20% decrease at 30 µM in an Alamar Blue Assay. Moreover, no adverse effects were noted for these compounds when tested for teratogenicity and hepatotoxicity in zebrafish at various concentrations. Overall, being the only example of MtbCM inhibitors that showed effects on Mtb cell viability the compound 3b and 3c are of further interest form the view point of discovery and development of new anti-tubercular agents.

Can the local electric field be a descriptor of catalytic activity? A case study on chorismate mutase.

The current theoretical perception of enzymatic activity is highly reliant on the determination of the activation energy of the reactions, which is often calculated using computationally demanding quantum mechanical calculations. With the ever-increasing use of bioengineering techniques that produce too many variants of the same enzyme, a fast and accurate way to study the relative efficiency of enzymes is currently in high demand. Here, we propose the local electric field (LEF) of the enzyme along the reaction axis as a descriptor for the enzymatic activity using the example of chorismate mutase in its native form and several variants (R90A, R90G, and R90K/C88S). The study shows a direct correlation between the calculated enzymatic EF and the enzymatic activity for all the complexes. MD simulations of the Michaelis complex and the transition state analog (TSA) show a stabilizing force on the TSA due to the enzymatic EF. QM/MM and QM-only DFT calculations in the presence of an external electric field (EEF) oriented along the reaction axis show that the electric field can interact with the dipole moment of the TS, thereby stabilizing it and thus lowering the activation energy.

Evolving the naturally compromised chorismate mutase from Mycobacterium tuberculosis to top performance.

Chorismate mutase (CM), an essential enzyme at the branch-point of the shikimate pathway, is required for the biosynthesis of phenylalanine and tyrosine in bacteria, archaea, plants, and fungi. MtCM, the CM from Mycobacterium tuberculosis, has less than 1% of the catalytic efficiency of a typical natural CM and requires complex formation with 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase for high activity. To explore the full potential of MtCM for catalyzing its native reaction, we applied diverse iterative cycles of mutagenesis and selection, thereby raising kcat/Km 270-fold to 5 × 105m-1s-1, which is even higher than for the complex. Moreover, the evolutionarily optimized autonomous MtCM, which had 11 of its 90 amino acids exchanged, was stabilized compared with its progenitor, as indicated by a 9 °C increase in melting temperature. The 1.5 Å crystal structure of the top-evolved MtCM variant reveals the molecular underpinnings of this activity boost. Some acquired residues (e.g. Pro52 and Asp55) are conserved in naturally efficient CMs, but most of them lie beyond the active site. Our evolutionary trajectories reached a plateau at the level of the best natural enzymes, suggesting that we have exhausted the potential of MtCM. Taken together, these findings show that the scaffold of MtCM, which naturally evolved for mediocrity to enable inter-enzyme allosteric regulation of the shikimate pathway, is inherently capable of high activity.

A Secreted Chorismate Mutase from Xanthomonas arboricola pv. juglandis Attenuates Virulence and Walnut Blight Symptoms.

Walnut blight is a significant above-ground disease of walnuts caused by Xanthomonas arboricola pv. juglandis (Xaj). The secreted form of chorismate mutase (CM), a key enzyme of the shikimate pathway regulating plant immunity, is highly conserved between plant-associated beta and gamma proteobacteria including phytopathogens belonging to the Xanthomonadaceae family. To define its role in walnut blight disease, a dysfunctional mutant of chorismate mutase was created in a copper resistant strain Xaj417 (XajCM). Infections of immature walnut Juglans regia (Jr) fruit with XajCM were hypervirulent compared with infections with the wildtype Xaj417 strain. The in vitro growth rate, size and cellular morphology were similar between the wild-type and XajCM mutant strains, however the quantification of bacterial cells by dPCR within walnut hull tissues showed a 27% increase in XajCM seven days post-infection. To define the mechanism of hypervirulence, proteome analysis was conducted to compare walnut hull tissues inoculated with the wild type to those inoculated with the XajCM mutant strain. Proteome analysis revealed 3296 Jr proteins (five decreased and ten increased with FDR ≤ 0.05) and 676 Xaj417 proteins (235 increased in XajCM with FDR ≤ 0.05). Interestingly, the most abundant protein in Xaj was a polygalacturonase, while in Jr it was a polygalacturonase inhibitor. These results suggest that this secreted chorismate mutase may be an important virulence suppressor gene that regulates Xaj417 virulence response, allowing for improved bacterial survival in the plant tissues.

Unravelling the Allosteric Targeting of PHGDH at the ACT-Binding Domain with a Photoactivatable Diazirine Probe and Mass Spectrometry Experiments.

The serine biosynthetic pathway is a key element contributing to tumor proliferation. In recent years, targeting of phosphoglycerate dehydrogenase (PHGDH), the first enzyme of this pathway, intensified and revealed to be a promising strategy to develop new anticancer drugs. Among attractive PHGDH inhibitors are the α-ketothioamides. In previous work, we have demonstrated their efficacy in the inhibition of PHGDH in vitro and in cellulo. However, the precise site of action of this series, which would help the rational design of new inhibitors, remained undefined. In the present study, the detailed mechanism-of-action of a representative α-ketothioamide inhibitor is reported using several complementary experimental techniques. Strikingly, our work led to the identification of an allosteric site on PHGDH that can be targeted for drug development. Using mass spectrometry experiments and an original α-ketothioamide diazirine-based photoaffinity probe, we identified the 523Q-533F sequence on the ACT regulatory domain of PHGDH as the binding site of α-ketothioamides. Mutagenesis experiments further documented the specificity of our compound at this allosteric site. Our results thus pave the way for the development of new anticancer drugs using a completely novel mechanism-of-action.

Different Solvent and Conformational Entropy Contributions to the Allosteric Activation and Inhibition Mechanisms of Yeast Chorismate Mutase.

Allosteric regulation is important in many biological processes, including cell signaling, gene regulation, and metabolism. Saccharomyces cerevisiae chorismate mutase (ScCM) is a key homodimeric enzyme in the shikimate pathway responsible for the generation of aromatic amino acids, where it is allosterically inhibited and activated by Tyr and Trp, respectively. Our previous studies indicated that binding of both allosteric effectors is negatively cooperative, that is binding at one allosteric binding site discourages binding at the other, due to the entropic penalty of binding the second allosteric effector. We utilized variable temperature isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR) experiments to better understand the entropic contributions to allosteric effector binding, including changes to solvent entropy and protein conformational entropy. Upon binding either Tyr or Trp, ScCM experiences a quenching of motions on the picosecond-to-nanosecond time scale, which we could relate to a loss of protein conformational entropy. Further ITC and NMR studies were consistent with the Tyr-bound form of ScCM being associated with more water molecules compared to the Trp-bound form and Tyr binding being associated with a less positive solvent entropy change. These studies provide insight into the role of structural dynamics in ScCM function and highlight the importance of solvent entropy changes in allosteric regulation, a historically underappreciated concept.

The lineage and diversity of putative amino acid sensor ACR proteins in plants.

Amino acid metabolic enzymes often contain a regulatory ACT domain, named for aspartate kinase, chorismate mutase, and TyrA (prephenate dehydrogenase). Arabidopsis encodes 12 putative amino acid sensor ACT repeat (ACR) proteins, all containing ACT repeats but no identifiable catalytic domain. Arabidopsis ACRs comprise three groups based on domain composition and sequence: group I and II ACRs contain four ACTs each, and group III ACRs contain two ACTs. Previously, all three groups had been documented only in Arabidopsis. Here, we extended this to algae and land plants, showing that all three groups of ACRs are present in most, if not all, land plants, whereas among algal ACRs, although quite diverse, only group III is conserved. The appearance of canonical group I and II ACRs thus accompanied the evolution of plants from living in water to living on land. Alignment of ACTs from plant ACRs revealed a conserved motif, DRPGLL, at the putative ligand-binding site. Notably, the unique features of the DRPGLL motifs in each ACT domain are conserved in ACRs from algae to land plants. The conservation of plant ACRs is reminiscent of that of human cellular arginine sensor for mTORC1 (CASTOR1), a member of a small protein family highly conserved in animals. CASTOR proteins also have four ACT domains, although the sequence identities between ACRs and CASTORs are very low. Thus, plant ACRs and animal CASTORs may have adapted the regulatory ACT domains from a more ancient metabolic enzyme, and then evolved independently.

De novo resveratrol production through modular engineering of an Escherichia coli-Saccharomyces cerevisiae co-culture.

BACKGROUND: Resveratrol is a plant secondary metabolite with diverse, potential health-promoting benefits. Due to its nutraceutical merit, bioproduction of resveratrol via microbial engineering has gained increasing attention and provides an alternative to unsustainable chemical synthesis and straight extraction from plants. However, many studies on microbial resveratrol production were implemented with the addition of water-insoluble phenylalanine or tyrosine-based precursors to the medium, limiting in the sustainable development of bioproduction. RESULTS: Here we present a novel coculture platform where two distinct metabolic background species were modularly engineered for the combined total and de novo biosynthesis of resveratrol. In this scenario, the upstream Escherichia coli module is capable of excreting p-coumaric acid into the surrounding culture media through constitutive overexpression of codon-optimized tyrosine ammonia lyase from Trichosporon cutaneum (TAL), feedback-inhibition-resistant 3-deoxy-d-arabinoheptulosonate-7-phosphate synthase (aroGfbr) and chorismate mutase/prephenate dehydrogenase (tyrAfbr) in a transcriptional regulator tyrR knockout strain. Next, to enhance the precursor malonyl-CoA supply, an inactivation-resistant version of acetyl-CoA carboxylase (ACC1S659A,S1157A) was introduced into the downstream Saccharomyces cerevisiae module constitutively expressing codon-optimized 4-coumarate-CoA ligase from Arabidopsis thaliana (4CL) and resveratrol synthase from Vitis vinifera (STS), and thus further improve the conversion of p-coumaric acid-to-resveratrol. Upon optimization of the initial inoculation ratio of two populations, fermentation temperature, and culture time, this co-culture system yielded 28.5 mg/L resveratrol from glucose in flasks. In further optimization by increasing initial net cells density at a test tube scale, a final resveratrol titer of 36 mg/L was achieved. CONCLUSIONS: This is first study that demonstrates the use of a synthetic E. coli-S. cerevisiae consortium for de novo resveratrol biosynthesis, which highlights its potential for production of other p-coumaric-acid or resveratrol derived biochemicals.

Energy and Enzyme Activity Landscapes of Yeast Chorismate Mutase at Cellular Concentrations of Allosteric Effectors.

In solution, proteins fluctuate among many conformational substates, with their relative free energies determining substate populations and energy barriers determining conformational exchange kinetics. It has been suggested that members of the conformational ensemble may be responsible for different protein functions, although it is generally difficult to test such a proposal in most systems. A model protein for deciphering individual substate contributions is the homodimeric Saccharomyces cerevisiae chorismate mutase (ScCM) enzyme, which is negatively and positively regulated by tyrosine and tryptophan, respectively. Previous X-ray crystallography structures revealed two equivalent allosteric binding pockets that can be occupied by either tryptophan or tyrosine. We proposed that under cellular conditions there are six potential states of ScCM: no allosteric effector bound, a single tyrosine bound, a single tryptophan bound, two tyrosines bound, two tryptophans bound, and a mixed bound state in which tyrosine and tryptophan occupy different allosteric sites. We used isothermal titration calorimetry and solution-state nuclear magnetic resonance spectroscopy to confirm the existence of all six states and construct the complete six-state equilibrium binding profile. We were also able to assign enzyme activities to each state, which allowed us to derive the enzyme activity landscape across the range of cellular concentrations of tyrosine and tryptophan. Surprisingly, the mixed bound state had the highest enzyme activity, which suggested that the shikimate pathway is shunted toward tyrosine production under most conditions.

Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution.

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.

Synthesis of 3-indolylmethyl substituted (pyrazolo/benzo)triazinone derivatives under Pd/Cu-catalysis: Identification of potent inhibitors of chorismate mutase (CM).

The chorismate mutase (CM) is considered as an attractive target for the identification of potential antitubercular agents due to its absence in animals but not in bacteria. A series of 3-indolylmethyl substituted pyrazolotriazinone derivatives were designed and docked into CM in silico as potential inhibitors. These compounds were efficiently synthesized using the Pd/Cu-catalyzed coupling-cyclization in a single pot involving the construction of indole ring. The methodology was later extended to the preparation of corresponding benzo analogs of pyrazolotriazinones i.e. 3-indolylmethyl substituted benzotriazinone derivatives. Several of these novel compounds showed significant inhibition of CM when tested in vitro at 30 µM. The SAR (Structure-Activity-Relationship) studies suggested that benzotriazinone moiety was more favorable over the pyrazolotriazinone ring. The two best active compounds showed IC50 ∼ 0.4-0.9 µM (better than the reference/known compounds used) and no toxicity till 30 µM in vitro.

Inter-Enzyme Allosteric Regulation of Chorismate Mutase in Corynebacterium glutamicum: Structural Basis of Feedback Activation by Trp.

Corynebacterium glutamicum is widely used for the industrial production of amino acids, nucleotides, and vitamins. The shikimate pathway enzymes DAHP synthase (CgDS, Cg2391) and chorismate mutase (CgCM, Cgl0853) play a key role in the biosynthesis of aromatic compounds. Here we show that CgCM requires the formation of a complex with CgDS to achieve full activity, and that both CgCM and CgDS are feedback regulated by aromatic amino acids binding to CgDS. Kinetic analysis showed that Phe and Tyr inhibit CgCM activity by inter-enzyme allostery, whereas binding of Trp to CgDS strongly activates CgCM. Mechanistic insights were gained from crystal structures of the CgCM homodimer, tetrameric CgDS, and the heterooctameric CgCM-CgDS complex, refined to 1.1, 2.5, and 2.2 Å resolution, respectively. Structural details from the allosteric binding sites reveal that DAHP synthase is recruited as the dominant regulatory platform to control the shikimate pathway, similar to the corresponding enzyme complex from Mycobacterium tuberculosis.

Evolution of allosteric regulation in chorismate mutases from early plants.

Plants, fungi, and bacteria synthesize the aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan. Chorismate mutase catalyzes the branch point reaction of phenylalanine and tyrosine biosynthesis to generate prephenate. In Arabidopsis thaliana, there are two plastid-localized chorismate mutases that are allosterically regulated (AtCM1 and AtCM3) and one cytosolic isoform (AtCM2) that is unregulated. Previous analysis of plant chorismate mutases suggested that the enzymes from early plants (i.e. bryophytes/moss, lycophytes, and basal angiosperms) formed a clade distinct from the isoforms found in flowering plants; however, no biochemical information on these enzymes is available. To understand the evolution of allosteric regulation in plant chorismate mutases, we analyzed a basal lineage of plant enzymes homologous to AtCM1 based on sequence similarity. The chorismate mutases from the moss/bryophyte Physcomitrella patens (PpCM1 and PpCM2), the lycophyte Selaginella moellendorffii (SmCM), and the basal angiosperm Amborella trichopoda (AmtCM1 and AmtCM2) were characterized biochemically. Tryptophan was a positive effector for each of the five enzymes examined. Histidine was a weak positive effector for PpCM1 and AmtCM1. Neither tyrosine nor phenylalanine altered the activity of SmCM; however, tyrosine was a negative regulator of the other four enzymes. Phenylalanine down-regulates both moss enzymes and AmtCM2. The 2.0 ŠX-ray crystal structure of PpCM1 in complex with the tryptophan identified the allosteric effector site and reveals structural differences between the R- (more active) and T-state (less active) forms of plant chorismate mutases. Molecular insight into the basal plant chorismate mutases guides our understanding of the evolution of allosteric regulation in these enzymes.