Crystal structures of the monofunctional chorismate mutase from Bacillus subtilis and its complex with a transition state analog

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.

Exploration of biochemical reactivity with a QM/MM growing string method

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions

In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

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 k_cat/K_m 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.

Nonenzymatic Reactions in Natural Product Formation

Biosynthetic mechanisms of natural products primarily depend on systems of protein catalysts. However, within the field of biosynthesis, there are cases in which the inherent chemical reactivity of metabolic intermediates and substrates evades the involvement of enzymes. These reactions are difficult to characterize based on their reactivity and occlusion within the milieu of the cellular environment. As we continue to build a strong foundation for how microbes and higher organisms produce natural products, therein lies a need for understanding how protein independent or nonenzymatic biosynthetic steps can occur. We have classified such reactions into four categories: intramolecular, multicomponent, tailoring, and light-induced reactions. Intramolecular reactions is one of the most well studied in the context of biomimetic synthesis, consisting of cyclizations and cycloadditions due to the innate reactivity of the intermediates. There are two subclasses that make up multicomponent reactions, one being homologous multicomponent reactions which results in dimeric and pseudodimeric natural products, and the other being heterologous multicomponent reactions, where two or more precursors from independent biosynthetic pathways undergo a variety of reactions to produce the mature natural product. The third type of reaction discussed are tailoring reactions, where postmodifications occur on the natural products after the biosynthetic machinery is completed. The last category consists of light-induced reactions involving ecologically relevant UV light rather than high intensity UV irradiation that is traditionally used in synthetic chemistry. This review will cover recent nonenzymatic biosynthetic mechanisms and include sources for those reviewed previously.

Multistructural microiteration combined with QM/MM-ONIOM electrostatic embedding

Multistructural microiteration (MSM) is a method to take account of contributions of multiple surrounding structures in a geometrical optimization or reaction path calculation using the quantum mechanics/molecular mechanics (QM/MM) ONIOM method. In this study, we combined MSM with the electrostatic embedding (EE) scheme of the QM/MM-ONIOM method by extending its original formulation for mechanical embedding (ME). MSM-EE takes account of the polarization in the QM region induced by point charges assigned to atoms in the multiple surrounding structures, where the point charges are scaled by the weight factor of each surrounding structure determined through MSM. The performance of MSM-EE was compared with that of the other methods, i.e., ONIOM-ME, ONIOM-EE, and MSM-ME, by applying them to three chemical processes: (1) chorismate-to-prephenate transformation in aqueous solution, (2) the same transformation as (1) in an enzyme, and (3) hydroxylation in p-hydroxybenzoate hydroxylase. These numerical tests of MSM-EE yielded barriers and reaction energies close to experimental values with computational costs comparable to those of the other three methods.

Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian

Molecular dynamics (MD) simulations employing ab initio quantum mechanical and molecular mechanical (ai-QM/MM) potentials are considered to be the state of the art, but the high computational cost associated with the ai-QM calculations remains a theoretical challenge for their routine application. Here, we present a modified protocol of the multiple time step (MTS) method for accelerating ai-QM/MM MD simulations of condensed-phase reactions. Within a previous MTS protocol [Nam J. Chem. Theory Comput. 2014, 10, 4175], reference forces are evaluated using a low-level (semiempirical QM/MM) Hamiltonian and employed at inner time steps to propagate the nuclear motions. Correction forces, which arise from the force differences between high-level (ai-QM/MM) and low-level Hamiltonians, are applied at outer time steps, where the MTS algorithm allows the time-reversible integration of the correction forces. To increase the outer step size, which is bound by the highest-frequency component in the correction forces, the semiempirical QM Hamiltonian is recalibrated in this work to minimize the magnitude of the correction forces. The remaining high-frequency modes, which are mainly bond stretches involving hydrogen atoms, are then removed from the correction forces. When combined with a Langevin or SIN(R) thermostat, the modified MTS-QM/MM scheme remains robust with an up to 8 (with Langevin) or 10 fs (with SIN(R)) outer time step (with 1 fs inner time steps) for the chorismate mutase system. This leads to an over 5-fold speedup over standard ai-QM/MM simulations, without sacrificing the accuracy in the predicted free energy profile of the reaction.

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.

Machine-Learning-Assisted Free Energy Simulation of Solution-Phase and Enzyme Reactions

Despite recent advances in the development of machine learning potentials (MLPs) for biomolecular simulations, there has been limited effort on developing stable and accurate MLPs for enzymatic reactions. Here we report a protocol for performing machine-learning-assisted free energy simulation of solution-phase and enzyme reactions at the ab initio quantum-mechanical/molecular-mechanical (ai-QM/MM) level of accuracy. Within our protocol, the MLP is built to reproduce the ai-QM/MM energy and forces on both QM (reactive) and MM (solvent/enzyme) atoms. As an alternative strategy, a delta machine learning potential (ΔMLP) is trained to reproduce the differences between the ai-QM/MM and semiempirical (se) QM/MM energies and forces. To account for the effect of the condensed-phase environment in both MLP and ΔMLP, the DeePMD representation of a molecular system is extended to incorporate the external electrostatic potential and field on each QM atom. Using the Menshutkin and chorismate mutase reactions as examples, we show that the developed MLP and ΔMLP reproduce the ai-QM/MM energy and forces with errors that on average are less than 1.0 kcal/mol and 1.0 kcal mol-1 Å-1, respectively, for representative configurations along the reaction pathway. For both reactions, MLP/ΔMLP-based simulations yielded free energy profiles that differed by less than 1.0 kcal/mol from the reference ai-QM/MM results at only a fraction of the computational cost.

Discovery and investigation of natural Diels-Alderases

It has been proposed that biosyntheses of many natural products involve pericyclic reactions, including Diels-Alder (DA) reaction. However, only a small set of enzymes have been proposed to catalyze pericyclic reactions. Most surprisingly, there has been no formal identification of natural enzymes that can be defined to catalyze DA reactions (DAases), despite the wide application of the reaction in chemical syntheses of complex organic compounds. However, recent studies began to accumulate a growing body of evidence that supports the notion that enzymes that formally catalyze DA reactions, in fact exist. In this review, I will begin by describing a short history behind the discovery and characterization of macrophomate synthase, one of the earliest enzymes that was proposed to catalyze an intermolecular DA reaction during the biosynthesis of a substituted benzoic acid in a phytopathogenic fungus Macrophoma commelinae. Then, I will discuss representative enzymes that have been chemically authenticated to catalyze DA reactions, with emphasis on more recent discoveries of DAases involved mainly in fungal secondary metabolite biosynthesis except for one example from a marine streptomycete. The current success in identification of a series of DAases and enzymes that catalyze other pericyclic reactions owes to the combined efforts from both the experimental and theoretical approaches in discovering natural products. Such efforts typically involve identifying the chemical features derived from cycloaddition reactions, isolating the biosynthetic genes that encode enzymes that generate such chemical features and deciphering the reaction mechanisms for the enzyme-catalyzed pericyclic reactions.

Using a genetic algorithm to find molecules with good docking scores

A graph-based genetic algorithm (GA) is used to identify molecules (ligands) with high absolute docking scores as estimated by the Glide software package, starting from randomly chosen molecules from the ZINC database, for four different targets: Bacillus subtilis chorismate mutase (CM), human β2-adrenergic G protein-coupled receptor (β2AR), the DDR1 kinase domain (DDR1), and β-cyclodextrin (BCD). By the combined use of functional group filters and a score modifier based on a heuristic synthetic accessibility (SA) score our approach identifies between ca 500 and 6,000 structurally diverse molecules with scores better than known binders by screening a total of 400,000 molecules starting from 8,000 randomly selected molecules from the ZINC database. Screening 250,000 molecules from the ZINC database identifies significantly more molecules with better docking scores than known binders, with the exception of CM, where the conventional screening approach only identifies 60 compounds compared to 511 with GA+Filter+SA. In the case of β2AR and DDR1, the GA+Filter+SA approach finds significantly more molecules with docking scores lower than −9.0 and −10.0. The GA+Filters+SA docking methodology is thus effective in generating a large and diverse set of synthetically accessible molecules with very good docking scores for a particular target. An early incarnation of the GA+Filter+SA approach was used to identify potential binders to the COVID-19 main protease and submitted to the early stages of the COVID Moonshot project, a crowd-sourced initiative to accelerate the development of a COVID antiviral.

Hybrid Peptide-Thiourea Catalyst for Asymmetric Michael Additions of Aldehydes to Heterocyclic Nitroalkenes

Bifunctional organocatalysis combining covalent and noncovalent activation is presented. The hybrid peptide-thiourea catalyst features a N-terminal proline moiety for aldehyde activation and a thiourea unit for electrophile activation. This catalyst effectively promotes asymmetric Michael additions of aldehydes to challenging but biologically relevant heterocycle-containing nitroalkenes. The catalyst can be used under solvent-free conditions. Spectroscopic and density functional theory studies elucidate the catalyst structure and mode of action.

Chorismate mutase and isochorismatase, two potential effectors of the migratory nematode Hirschmanniella oryzae, increase host susceptibility by manipulating secondary metabolite content of rice

Hirschmanniella oryzae is one of the most devastating nematodes on rice, leading to substantial yield losses. Effector proteins aid the nematode during the infection process by subduing plant defence responses. In this research we characterized two potential H. oryzae effector proteins, chorismate mutase (HoCM) and isochorismatase (HoICM), and investigated their enzymatic activity and their role in plant immunity. Both HoCM and HoICM proved to be enzymatically active in complementation tests in mutant Escherichia coli strains. Infection success by the migratory nematode H. oryzae was significantly higher in transgenic rice lines constitutively expressing HoCM or HoICM. Expression of HoCM, but not HoICM, increased rice susceptibility against the sedentary nematode Meloidogyne graminicola also. Transcriptome and metabolome analyses indicated reductions in secondary metabolites in the transgenic rice plants expressing the potential nematode effectors. The results presented here demonstrate that both HoCM and HoICM suppress the host immune system and that this may be accomplished by lowering secondary metabolite levels in the plant.

Biophysical and In-Silico Studies of Phytochemicals Targeting Chorismate Synthase from Drug-Resistant Moraxella Catarrhalis

Chorismate serves as a crucial precursor for the synthesis of many aromatic compounds essential for the survival and virulence in various bacteria and protozoans. Chorismate synthase, a vital enzyme in the shikimate pathway, is responsible for the formation of chorismate from enolpyruvylshikimate-3-phosphate (EPSP). Moraxella catarrhalis is reported to be resistant to many beta-lactam antibiotics and causes chronic ailments such as otitis media, sinusitis, laryngitis, and bronchopulmonary infections. Here, we have cloned the aroC gene from Moraxella catarrhalis in pET28c and heterologously produced the chorismate synthase (~ 43 kDa) in Escherichia coli BL21(DE3) cells. We have predicted the three-dimensional structure of this enzyme and used the refined model for ligand-based virtual screening against Supernatural Database using PyRx tool that led to the identification of the top three molecules (caffeic acid, gallic acid, and o-coumaric acid). The resultant protein-ligand complex structures were subjected to 50 ns molecular dynamics (MD) simulation using GROMACS. Further, the binding energy was calculated by MM/PBSA approach using the trajectory obtained from MD simulation. The binding affinities of these compounds were validated with ITC experiments, which suggest that gallic acid has the highest binding affinity amongst these three phytochemicals. Together, these results pave the way for the use of these phytochemicals as potential anti-bacterial compounds.

Combining Random Mutagenesis and Metabolic Engineering for Enhanced Tryptophan Production in Synechocystis sp. Strain PCC 6803

Tryptophan (Trp) is an essential aromatic amino acid that has value as an animal feed supplement, as the amount found in plant-based sources is insufficient. An alternative to production by engineered microbial fermentation is to have tryptophan biosynthesized by a photosynthetic microorganism that could replace or supplement both the plant and industrially used microbes. We selected Synechocystis sp. strain PCC 6803, a model cyanobacterium, as the host and studied metabolic engineering and random mutagenesis approaches. Previous work on engineering heterotrophic microbes for improved Trp titers has targeted allosteric feedback regulation in enzymes 3-deoxy-d-arabinoheptulosonate 7-phosphate synthase (DAHPS) and anthranilate synthase (AS) as major bottlenecks in the shikimate pathway. In this work, the genes encoding feedback-resistant enzymes from Escherichia coli, aroG^fbr and trpE^fbr , were overexpressed in the host wild-type (WT) strain. Separately, the WT strain was subjected to random mutagenesis and selection using an amino acid analog to isolate tryptophan-overproducing strains. The randomly mutagenized strains were sequenced in order to identify the mutations that resulted in the desirable phenotypes. Interestingly, the tryptophan overproducers had mutations in the gene encoding chorismate mutase (CM), which catalyzes the conversion of chorismate to prephenate. The best tryptophan overproducer from random mutagenesis was selected as a host for metabolic engineering where aroG^fbr and trpE^fbr were overexpressed. The best strain developed produced 212 ± 23 mg/liter of tryptophan after 10 days of photoautotrophic growth under 3% (vol/vol) CO₂ We demonstrated that a combination of random mutagenesis and metabolic engineering was superior to either individual approach.IMPORTANCE Aromatic amino acids such as tryptophan are primarily used as additives in the animal feed industry and are typically produced using genetically engineered heterotrophic organisms such as Escherichia coli This involves a two-step process, where the substrate such as molasses is first obtained from plants followed by fermentation by heterotrophic organisms. We have engineered photoautotrophic cyanobacterial strains by a combination of random mutagenesis and metabolic engineering. These strains grow on CO₂ as the sole carbon source and utilize light as the sole energy source to produce tryptophan, thus converting the two-step process into a single step. Our results show that combining random mutagenesis and metabolic engineering was superior to either approach alone. This study also builds a foundation for further engineering of cyanobacteria for industrial tryptophan production.

Accelerated computation of free energy profile at ab initio quantum mechanical/molecular mechanical accuracy via a semi-empirical reference potential. II. Recalibrating semi-empirical parameters with force matching

An efficient and accurate reference potential simulation protocol is proposed for producing ab initio quantum mechanical/molecular mechanical (AI-QM/MM) quality free energy profiles for chemical reactions in a solvent or macromolecular environment. This protocol involves three stages: (a) using force matching to recalibrate a semi-empirical quantum mechanical (SE-QM) Hamiltonian for the specific reaction under study; (b) employing the recalibrated SE-QM Hamiltonian (in combination with molecular mechanical force fields) as the reference potential to drive umbrella samplings along the reaction pathway; and (c) computing AI-QM/MM energy values for collected configurations from the sampling and performing weighted thermodynamic perturbation to acquire an AI-QM/MM corrected reaction free energy profile. For three model reactions (identity SN2 reaction, Menshutkin reaction, and glycine proton transfer reaction) in aqueous solution and one enzyme reaction (Claisen arrangement in chorismate mutase), our simulations using recalibrated PM3 SE-QM Hamiltonians well reproduced QM/MM free energy profiles at the B3LYP/6-31G* level of theory all within 1 kcal mol-1 with a 20 to 45 fold reduction in the computer time.

Structural basis for stereoselective dehydration and hydrogen-bonding catalysis by the SAM-dependent pericyclase LepI

LepI is an S-adenosylmethionine (SAM)-dependent pericyclase that catalyzes the formation of 2-pyridone natural product leporin C. Biochemical characterization showed LepI can catalyze the stereoselective dehydration to yield a reactive (E)-quinone methide that can undergo bifurcating intramolecular Diels-Alder (IMDA) and hetero-Diels-Alder (HDA) cyclizations from an ambimodal transition state, as well as a [3,3]-retro-Claisen rearrangement to recycle the IMDA product into leporin C. Here we solved the X-ray crystal structures of SAM-bound LepI and in complex with a substrate analog, the product leporin C, and a retro-Claisen reaction transition-state analog to understand the structural basis for the multitude of reactions. Structural and mutational analysis revealed how Nature evolves a classic methyltransferase active site into one that can serve as a dehydratase and a multifunctional pericyclase. Catalysis of both sets of reactions employs H133 and R295, two active site residues that are not found in canonical methyltransferases. An alternative role of SAM, which is not found to be in direct contact with the substrate, is also proposed.

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.

Rational Discovery of (+) (S) Abscisic Acid as a Potential Antifungal Agent: a Repurposing Approach

Fungal infections are spreading widely worldwide, and the types of treatment are limited due to the lack of diverse therapeutic agents and their associated side effects and toxicity. The discovery of new antifungal classes is vital and critical. We discovered the antifungal activity of abscisic acid through a rational drug design methodology that included the building of homology models for fungal chorismate mutases and a pharmacophore model derived from a transition state inhibitor. Ligand-based virtual screening resulted in some hits that were filtered using molecular docking and molecular dynamic simulations studies. Both in silico methods and in vitro antifungal assays were used as tools to select and validate the abscisic acid repurposing. Abscisic acid inhibition assays confirmed the inhibitory effect of abscisic acid on chorismate mutase through the inhibition of phenylpyruvate production. The repositioning of abscisic acid, the well-known and naturally occurring plant growth regulator, as a potential antifungal agent because of its suggested action as an inhibitor to several fungal chorismate mutases was the main result of this work.

Tyrosine biosynthesis, metabolism, and catabolism in plants

L-Tyrosine (Tyr) is an aromatic amino acid (AAA) required for protein synthesis in all organisms, but synthesized de novo only in plants and microorganisms. In plants, Tyr also serves as a precursor of numerous specialized metabolites that have diverse physiological roles as electron carriers, antioxidants, attractants, and defense compounds. Some of these Tyr-derived plant natural products are also used in human medicine and nutrition (e.g. morphine and vitamin E). While the Tyr biosynthesis and catabolic pathways have been extensively studied in microbes and animals, respectively, those of plants have received much less attention until recently. Accumulating evidence suggest that the Tyr biosynthetic pathways differ between microbes and plants and even within the plant kingdom, likely to support the production of lineage-specific plant specialized metabolites derived from Tyr. The interspecies variations of plant Tyr pathway enzymes can now be used to enhance the production of Tyr and Tyr-derived compounds in plants and other synthetic biology platforms.

Structure of Chorismate Mutase-like Domain of DAHPS from Bacillus subtilis Complexed with Novel Inhibitor Reveals Conformational Plasticity of Active Site

3-deoxy-D-arabino-heptulosonate-7-phosphate-synthase (DAHPS) is the first enzyme of the shikimate pathway and is responsible for the synthesis of aromatic amino acids in microorganisms. This pathway is an attractive target for antimicrobial drugs. In Bacillus subtilis, the N-terminal domain of the bifunctional DAHPS enzyme belongs to an AroQ class of chorismate mutase and is functionally homologous to the downstream AroH class chorismate mutase. This is the first structure of chorismate mutase, AroQ (BsCM_2) enzyme from Bacillus subtilis in complex with citrate and chlorogenic acid at 1.9 Å and 1.8 Å resolution, respectively. This work provides the structural basis of ligand binding into the active site of AroQ class of chorismate mutase, while accompanied by the conformational flexibility of active site loop. Molecular dynamics results showed that helix H2′ undergoes uncoiling at the first turn and increases the mobility of loop L1′. The side chains of Arg45, Phe46, Arg52 and Lys76 undergo conformational changes, which may play an important role in DAHPS regulation by the formation of the domain-domain interface. Additionally, binding studies showed that the chlorogenic acid binds to BsCM_2 with a higher affinity than chorismate. These biochemical and structural findings could lead to the development of novel antimicrobial drugs.

X-Ray snapshots of a pyridoxal enzyme: a catalytic mechanism involving concerted [1,5]-hydrogen sigmatropy in methionine γ-lyase

Pyridoxal 5′-phosphate (PLP)-enzymes are essentially involved in amino acid and amine metabolism of a wide variety of organisms. Despite their extensive biochemical studies, there are little evidence and structural data to comprehensively elaborate the catalytic mechanism. We obtained X-ray snapshots of l-methionine γ-lyase from Entamoeba histolytica (EhMGL), a PLP-enzyme catalyzing the γ-elimination reaction of methionine. Here, we suggest a catalytic mechanism of EhMGL by using the X-ray snapshots covering all stages of this multistep catalysis reaction. Initial formation of a Michaelis complex is followed by the migration of double bond from the C4′=Nα–Cα moiety in an intermediate PLP-methionine imine to C4′–Nα=Cα in pyridoxamine 5′-phosphate (PMP)-α,β-dehydromethionine imine without intervention of a putative quinonoid intermediate. The enzyme can facilitate the subsequent γ-elimination of methanethiol by the possible general acid-base catalysis of Tyr108 for the E1cB mechanism, enabling to form the ene-imine C4′–Nα=Cα–Cβ=Cγ structure with the s-cis conformation, which is prerequisite for the non-enzymatic symmetry-allowed suprafacial [1,5]-hydrogen shift to complete the catalytic cycle by releasing α-ketobutyrate. The mechanism based on the X-ray snapshots is consistent with the reactivity of MGL toward methionine analogues. The generality of such a mechanism involving non-enzymatic concerted reaction in other PLP enzymes is discussed.

Quantum chemical modeling of the reaction path of chorismate mutase based on the experimental substrate/product complex

Chorismate mutase is a well-known model enzyme, catalyzing the Claisen rearrangement of chorismate to prephenate. Recent high-resolution crystal structures along the reaction coordinate of this enzyme enabled computational analyses at unprecedented detail. Using quantum chemical simulations, we investigated how the catalytic reaction mechanism is affected by electrostatic and hydrogen-bond interactions. Our calculations showed that the transition state (TS) was mainly stabilized electrostatically, with Arg90 playing the leading role. The effect was augmented by selective hydrogen-bond formation to the TS in the wild-type enzyme, facilitated by a small-scale local induced fit. We further identified a previously underappreciated water molecule, which separates the negative charges during the reaction. The analysis includes the wild-type enzyme and a non-natural enzyme variant, where the catalytic arginine was replaced with an isosteric citrulline residue.

Mycobacterium tuberculosis chorismate mutase: A potential target for TB

Mycobacterium tuberculosis chorismate mutase (MtbCM) catalyzes the rearrangement of chorismate to prephenate in the shikimate biosynthetic pathway to form the essential amino acids, phenylalanine and tyrosine. Two genes encoding chorismate mutase have been identified in Mtb. The secretory form,∗MtbCM (encoded by Rv1885c) is assumed to play a key role in pathogenesis of tuberculosis. Also, the inhibition of MtbCM may hinder the supply of nutrients to the organism. Indeed, the existence of chorismate mutase (CM) in bacteria, fungi and higher plants but not in human and low sequence homology among known CM makes it an interesting target for the discovery of anti-tubercular agents. The present article mainly focuses on the recent developments in the structure, function and inhibition of MtbCM. The understanding of various aspects of MtbCM as presented in the current article may facilitate the design and subsequent chemical synthesis of new inhibitors against ∗MtbCM, that could lead to the discovery and development of novel and potent anti-tubercular agents in future.

The Cation-π Interaction in Small-Molecule Catalysis

Catalysis by small molecules (≤1000 Da, 10(-9)  m) that are capable of binding and activating substrates through attractive, noncovalent interactions has emerged as an important approach in organic and organometallic chemistry. While the canonical noncovalent interactions, including hydrogen bonding, ion pairing, and π stacking, have become mainstays of catalyst design, the cation-π interaction has been comparatively underutilized in this context since its discovery in the 1980s. However, like a hydrogen bond, the cation-π interaction exhibits a typical binding affinity of several kcal mol(-1) with substantial directionality. These properties render it attractive as a design element for the development of small-molecule catalysts, and in recent years, the catalysis community has begun to take advantage of these features, drawing inspiration from pioneering research in molecular recognition and structural biology. This Review surveys the burgeoning application of the cation-π interaction in catalysis.

Horizontal gene transfer of acetyltransferases, invertases and chorismate mutases from different bacteria to diverse recipients

BACKGROUND

Hoplolaimina plant-parasitic nematodes (PPN) are a lineage of animals with many documented cases of horizontal gene transfer (HGT). In a recent study, we reported on three likely HGT candidate genes in the soybean cyst nematode Heterodera glycines, all of which encode secreted candidate effectors with putative functions in the host plant. Hg-GLAND1 is a putative GCN5-related N-acetyltransferase (GNAT), Hg-GLAND13 is a putative invertase (INV), and Hg-GLAND16 is a putative chorismate mutase (CM), and blastp searches of the non-redundant database resulted in highest similarity to bacterial sequences. Here, we searched nematode and non-nematode sequence databases to identify all the nematodes possible that contain these three genes, and to formulate hypotheses about when they most likely appeared in the phylum Nematoda. We then performed phylogenetic analyses combined with model selection tests of alternative models of sequence evolution to determine whether these genes were horizontally acquired from bacteria.

RESULTS

Mining of nematode sequence databases determined that GNATs appeared in Hoplolaimina PPN late in evolution, while both INVs and CMs appeared before the radiation of the Hoplolaimina suborder. Also, Hoplolaimina GNATs, INVs and CMs formed well-supported clusters with different rhizosphere bacteria in the phylogenetic trees, and the model selection tests greatly supported models of HGT over descent via common ancestry. Surprisingly, the phylogenetic trees also revealed additional, well-supported clusters of bacterial GNATs, INVs and CMs with diverse eukaryotes and archaea. There were at least eleven and eight well-supported clusters of GNATs and INVs, respectively, from different bacteria with diverse eukaryotes and archaea. Though less frequent, CMs from different bacteria formed supported clusters with multiple different eukaryotes. Moreover, almost all individual clusters containing bacteria and eukaryotes or archaea contained species that inhabit very similar niches.

CONCLUSIONS

GNATs were horizontally acquired late in Hoplolaimina PPN evolution from bacteria most similar to the saprophytic and plant-pathogenic actinomycetes. INVs and CMs were horizontally acquired from bacteria most similar to rhizobacteria and Burkholderia soil bacteria, respectively, before the radiation of Hoplolaimina. Also, these three gene groups appear to have been frequent subjects of HGT from different bacteria to numerous, diverse lineages of eukaryotes and archaea, which suggests that these genes may confer important evolutionary advantages to many taxa. In the case of Hoplolaimina PPN, this advantage likely was an improved ability to parasitize plants.

Large-Scale Density Functional Theory Transition State Searching in Enzymes

Linear-scaling quantum mechanical density functional theory calculations have been applied to study the rearrangement of chorismate to prephenate in large-scale models of the Bacillus subtilis chorismate mutase enzyme. By treating up to 2000 atoms at a consistent quantum mechanical level of theory, we obtain an unbiased, almost parameter-free description of the transition state geometry and energetics. The activation energy barrier is calculated to be lowered by 10.5 kcal mol(-1) in the enzyme, compared with the equivalent reaction in water, which is in good agreement with experiment. Natural bond orbital analysis identifies a number of active site residues that are important for transition state stabilization in chorismate mutase. This benchmark study demonstrates that linear-scaling density functional theory techniques are capable of simulating entire enzymes at the ab initio quantum mechanical level of accuracy.

Structural evolution of differential amino acid effector regulation in plant chorismate mutases

Chorismate mutase converts chorismate into prephenate for aromatic amino acid biosynthesis. To understand the molecular basis of allosteric regulation in the plant chorismate mutases, we analyzed the three Arabidopsis thaliana chorismate mutase isoforms (AtCM1-3) and determined the x-ray crystal structures of AtCM1 in complex with phenylalanine and tyrosine. Functional analyses show a wider range of effector control in the Arabidopsis chorismate mutases than previously reported. AtCM1 is activated by tryptophan with phenylalanine and tyrosine acting as negative effectors; however, tryptophan, cysteine, and histidine activate AtCM3. AtCM2 is a nonallosteric form. The crystal structure of AtCM1 in complex with tyrosine and phenylalanine identifies differences in the effector sites of the allosterically regulated yeast enzyme and the other two Arabidopsis isoforms. Site-directed mutagenesis of residues in the effector site reveals key features leading to differential effector regulation in these enzymes. In AtCM1, mutations of Gly-213 abolish allosteric regulation, as observed in AtCM2. A second effector site position, Gly-149 in AtCM1 and Asp-132 in AtCM3, controls amino acid effector specificity in AtCM1 and AtCM3. Comparisons of chorismate mutases from multiple plants suggest that subtle differences in the effector site are conserved in different lineages and may lead to specialized regulation of this branch point enzyme.

Substrate distortion contributes to the catalysis of orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate (UMP) by 17 orders of magnitude. Eight new crystal structures with ligand analogues combined with computational analyses of the enzyme's short-lived intermediates and the intrinsic electronic energies to distort the substrate and other ligands improve our understanding of the still controversially discussed reaction mechanism. In their respective complexes, 6-methyl-UMP displays significant distortion of its methyl substituent bond, 6-amino-UMP shows the competition between the K72 and C6 substituents for a position close to D70, and the methyl and ethyl esters of OMP both induce rotation of the carboxylate group substituent out of the plane of the pyrimidine ring. Molecular dynamics and quantum mechanics/molecular mechanics computations of the enzyme-substrate complex also show the bond between the carboxylate group and the pyrimidine ring to be distorted, with the distortion contributing a 10-15% decrease of the ΔΔG(⧧) value. These results are consistent with ODCase using both substrate distortion and transition-state stabilization, primarily exerted by K72, in its catalysis of the OMP decarboxylation reaction.

Mapping enzymatic catalysis using the effective fragment molecular orbital method: towards all ab initio biochemistry

We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol⁻¹ for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively.

Computational enzyme design

Recent developments in computational chemistry and biology have come together in the "inside-out" approach to enzyme engineering. Proteins have been designed to catalyze reactions not previously accelerated in nature. Some of these proteins fold and act as catalysts, but the success rate is still low. The achievements and limitations of the current technology are highlighted and contrasted to other protein engineering techniques. On its own, computational "inside-out" design can lead to the production of catalytically active and selective proteins, but their kinetic performances fall short of natural enzymes. When combined with directed evolution, molecular dynamics simulations, and crowd-sourced structure-prediction approaches, however, computational designs can be significantly improved in terms of binding, turnover, and thermal stability.

Molecular dynamics simulation of the last step of a catalytic cycle: product release from the active site of the enzyme chorismate mutase from Mycobacterium tuberculosis

The protein chorismate mutase MtCM from Mycobacterium tuberculosis catalyzes one of the few pericyclic reactions known in biology: the transformation of chorismate to prephenate. Chorismate mutases have been widely studied experimentally and computationally to elucidate the transition state of the enzyme catalyzed reaction and the origin of the high catalytic rate. However, studies about substrate entry and product exit to and from the highly occluded active site of the enzyme have to our knowledge not been performed on this enzyme. Crystallographic data suggest a possible substrate entry gate, that involves a slight opening of the enzyme for the substrate to access the active site. Using multiple molecular dynamics simulations, we investigate the natural dynamic process of the product exiting from the binding pocket of MtCM. We identify a dominant exit pathway, which is in agreement with the gate proposed from the available crystallographic data. Helices H2 and H4 move apart from each other which enables the product to exit from the active site. Interestingly, in almost all exit trajectories, two residues arginine 72 and arginine 134, which participate in the burying of the active site, are accompanying the product on its exit journey from the catalytic site.

Modification of residue 42 of the active site loop with a lysine-mimetic side chain rescues isochorismate-pyruvate lyase activity in Pseudomonas aeruginosa PchB

PchB is an isochorismate-pyruvate lyase from Pseudomonas aeruginosa. A positively charged lysine residue is located in a flexible loop that behaves as a lid to the active site, and the lysine residue is required for efficient production of salicylate. A variant of PchB that lacks the lysine at residue 42 has a reduced catalytic free energy of activation of up to 4.4 kcal/mol. Construction of a lysine isosteric residue bearing a positive charge at the appropriate position leads to the recovery of 2.5-2.7 kcal/mol (about 60%) of the 4.4 kcal/mol by chemical rescue. Exogenous addition of ethylamine to the K42A variant leads to a neglible recovery of activity (0.180 kcal/mol, roughly 7% rescue), whereas addition of propylamine caused an additional modest loss in catalytic power (0.056 kcal/mol, or 2% loss). This is consistent with the view that (a) the lysine-42 residue is required in a specific conformation to stabilize the transition state and (b) the correct conformation is achieved for a lysine-mimetic side chain at site 42 in the course of loop closure, as expected for transition-state stabilization by the side chain ammonio function. That the positive charge is the main effector of transition state stabilization is shown by the construction of a lysine-isosteric residue capable of exerting steric effects and hydrogen bonding but not electrostatic effects, leading to a modest increase of catalytic power (0.267-0.505 kcal/mol of catalytic free energy, or roughly 6-11% rescue).

Hybrid schemes based on quantum mechanics/molecular mechanics simulations goals to success, problems, and perspectives

The development of characterization techniques, advanced synthesis methods, as well as molecular modeling has transformed the study of systems in a well-established research field. The current research challenges in biocatalysis and biotransformation evolve around enzyme discovery, design, and optimization. How can we find or create enzymes that catalyze important synthetic reactions, even reactions that may not exist in nature? What is the source of enzyme catalytic power? To answer these and other related questions, the standard strategies have evolved from trial-and-error methodologies based on chemical knowledge, accumulated experience, and common sense into a clearly multidisciplinary science that allows one to reach the molecular design of tailor-made enzyme catalysts. This is even more so when one refers to enzyme catalysts, for which the detailed structure and composition are known and can be manipulated to introduce well-defined residues which can be implicated in the chemical rearrangements taking place in the active site. The methods and techniques of theoretical and computational chemistry are becoming more and more important in both understanding the fundamental biological roles of enzymes and facilitating their utilization in biotechnology. Improvement of the catalytic function of enzymes is important from scientific and industrial viewpoints, and to put this fact in the actual perspective as well as the potentialities, we recommend the very recent report of Sanderson [Sanderson, K. (2011). Chemistry: enzyme expertise. Nature 471, 397.]. Great fundamental advances have been made toward the ab initio design of enzyme catalysts based on molecular modeling. This has been based on the molecular mechanistic knowledge of the reactions to be catalyzed, together with the development of advanced synthesis and characterization techniques. The corresponding molecular mechanism can be studied by means of powerful quantum chemical calculations. The catalytic active site can be optimized to improve the transition state analogues (TSA) and to enhance the catalytic activity, even improve the active site to favor a desired direction of some promiscuous enzymes. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. Current computational tools to assist experimentalists for the design and engineering of proteins with desired catalytic properties are described. The interplay between enzyme design, molecular simulations, and experiments will be presented to emphasize the interdisciplinary nature of this research field. This text highlights the recent advances and examples selected from our laboratory are shown, of how the applications of these tools are a first attempt to de novo design of protein active sites. Identification of neutral/advantageous/deleterious mutation platforms can be exploited to penetrate some of Nature's closely guarded secrets of chemical reactivity. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. The first part describes briefly how the molecular modeling is carried out. Then, we discuss the requirements of hybrid quantum mechanical/molecular mechanics molecular dynamics (QM/MM MD) simulations, analyzing what are the basis of these theoretical methodologies, how we can use them with a view to its application in the study of enzyme catalysis, and what are the best methodologies for assessing its catalytic potential. In the second part, we focus on some selected examples, taking as a common guide the chorismate to prephenate rearrangement, studying the corresponding molecular mechanism in vacuo, in solution and in an enzyme environment. In addition, examples involving catalytic antibodies (CAs) and promiscuous enzymes will be presented. Finally, a special emphasis is made to provide some hints about the logical evolution that can be anticipated in this research field. Moreover, it helps in understanding the open directions in this area of knowledge and highlights the importance of computational approaches in discovering specific drugs and the impact on the rational design of tailor-made enzymes.

The shikimate pathway and aromatic amino Acid biosynthesis in plants

L-tryptophan, L-phenylalanine, and L-tyrosine are aromatic amino acids (AAAs) that are used for the synthesis of proteins and that in plants also serve as precursors of numerous natural products, such as pigments, alkaloids, hormones, and cell wall components. All three AAAs are derived from the shikimate pathway, to which ≥30% of photosynthetically fixed carbon is directed in vascular plants. Because their biosynthetic pathways have been lost in animal lineages, the AAAs are essential components of the diets of humans, and the enzymes required for their synthesis have been targeted for the development of herbicides. This review highlights recent molecular identification of enzymes of the pathway and summarizes the pathway organization and the transcriptional/posttranscriptional regulation of the AAA biosynthetic network. It also identifies the current limited knowledge of the subcellular compartmentalization and the metabolite transport involved in the plant AAA pathways and discusses metabolic engineering efforts aimed at improving production of the AAA-derived plant natural products.

Pericyclic reactions catalyzed by chorismate-utilizing enzymes

One of the fundamental questions of enzymology is how catalytic power is derived. This review focuses on recent developments in the structure--function relationships of chorismate-utilizing enzymes involved in siderophore biosynthesis to provide insight into the biocatalysis of pericyclic reactions. Specifically, salicylate synthesis by the two-enzyme pathway in Pseudomonas aeruginosa is examined. The isochorismate-pyruvate lyase is discussed in the context of its homologues, the chorismate mutases, and the isochorismate synthase is compared to its homologues in the MST family (menaquinone, siderophore, or tryptophan biosynthesis) of enzymes. The tentative conclusion is that the activities observed cannot be reconciled by inspection of the active site participants alone. Instead, individual activities must arise from unique dynamic properties of each enzyme that are tuned to promote specific chemistries.

Entropic and enthalpic components of catalysis in the mutase and lyase activities of Pseudomonas aeruginosa PchB

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions, demonstrating the eponymous activity and also chorismate mutase activity. The thermodynamic parameters for these enzyme-catalyzed activities, as well as the uncatalyzed isochorismate decomposition, are reported from temperature dependence of k(cat) and k(uncat) data. The entropic effects do not contribute to enzyme catalysis as expected from previously reported chorismate mutase data. Indeed, an entropic penalty for the enzyme-catalyzed mutase reaction (ΔS(++) = -12.1 ± 0.6 cal/(mol K)) is comparable to that of the previously reported uncatalyzed reaction, whereas that of the enzyme-catalyzed lyase reaction (ΔS(++) = -24.3 ± 0.2 cal/(mol K)) is larger than that of the uncatalyzed lyase reaction (-15.77 ± 0.02 cal/(mol K)) documented here. With the assumption that chemistry is rate-limiting, we propose that a reactive substrate conformation is formed upon loop closure of the active site and that ordering of the loop contributes to the entropic penalty for converting the enzyme substrate complex to the transition state.

Transition-state charge stabilization through multiple non-covalent interactions in the guanidinium-catalyzed enantioselective Claisen rearrangement

The mechanism by which chiral arylpyrrole-substituted guanidinium ions promote the Claisen rearrangement of O-allyl α-ketoesters and induce enantioselectivity was investigated by experimental and computational methods. In addition to stabilization of the developing negative charge on the oxallyl fragment of the rearrangement transition state by hydrogen-bond donation, evidence was obtained for a secondary attractive interaction between the π-system of a catalyst aromatic substituent and the cationic allyl fragment. Across a series of substituted arylpyrrole derivatives, enantioselectivity was observed to vary predictably according to this proposal. This mechanistic analysis led to the development of a new p-dimethylaminophenyl-substituted catalyst, which afforded improvements in enantioselectivity relative to the parent phenyl catalyst for a representative set of substrates.

Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions

Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.

Structures of the first representatives of Pfam family PF06684 (DUF1185) reveal a novel variant of the Bacillus chorismate mutase fold and suggest a role in amino-acid metabolism

The crystal structures of BB2672 and SPO0826 were determined to resolutions of 1.7 and 2.1 Å by single-wavelength anomalous dispersion and multiple-wavelength anomalous dispersion, respectively, using the semi-automated high-throughput pipeline of the Joint Center for Structural Genomics (JCSG) as part of the NIGMS Protein Structure Initiative (PSI). These proteins are the first structural representatives of the PF06684 (DUF1185) Pfam family. Structural analysis revealed that both structures adopt a variant of the Bacillus chorismate mutase fold (BCM). The biological unit of both proteins is a hexamer and analysis of homologs indicates that the oligomer interface residues are highly conserved. The conformation of the critical regions for oligomerization appears to be dependent on pH or salt concentration, suggesting that this protein might be subject to environmental regulation. Structural similarities to BCM and genome-context analysis suggest a function in amino-acid synthesis.