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.

Accurate Computation of Thermodynamic Activation Parameters in the Chorismate Mutase Reaction from Empirical Valence Bond Simulations

Chorismate mutase (CM) enzymes have long served as model systems for benchmarking new methods and tools in computational chemistry. Despite the enzymes' prominence in the literature, the extent of the roles that activation enthalpy and entropy play in catalyzing the conversion of chorismate to prephenate is still subject to debate. Knowledge of these parameters is a key piece in fully understanding the mechanism of chorismate mutases. Within this study, we utilize EVB/MD free energy perturbation calculations at a range of temperatures, allowing us to extract activation enthalpies and entropies from an Arrhenius plot of activation free energies of the reaction catalyzed by a monofunctional Bacillus subtilis CM and the promiscuous enzyme isochorismate pyruvate lyase of Pseudomonas aeruginosa. In comparison to the uncatalyzed reaction, our results show that both enzyme-catalyzed reactions exhibit a substantial reduction in activation enthalpy, while the effect on activation entropy is relatively minor, demonstrating that enzyme-catalyzed CM reactions are enthalpically driven. Furthermore, we observe that the monofunctional CM from B. subtilis more efficiently catalyzes this reaction than its promiscuous counterpart. This is supported by a structural analysis of the reaction pathway at the transition state, from which we identified key residues explaining the enthalpically driven nature of the reactions and also the difference in efficiencies between the two enzymes.

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.

Crystal structure of chorismate mutase from Burkholderia thailandensis

Burkholderia thailandensis is often used as a model for more virulent members of this genus of proteobacteria that are highly antibiotic-resistant and are potential agents of biological warfare that are infective by inhalation. As part of ongoing efforts to identify potential targets for the development of rational therapeutics, the structures of enzymes that are absent in humans, including that of chorismate mutase from B. thailandensis, have been determined by the Seattle Structural Genomics Center for Infectious Disease. The high-resolution structure of chorismate mutase from B. thailandensis was determined in the monoclinic space group P21 with three homodimers per asymmetric unit. The overall structure of each protomer has the prototypical AroQγ topology and shares conserved binding-cavity residues with other chorismate mutases, including those with which it has no appreciable sequence identity.

Crystal structure of chorismate mutase from Burkholderia phymatum

The bacterium Burkholderia phymatum is a promiscuous symbiotic nitrogen-fixating bacterium that belongs to one of the largest groups of Betaproteobacteria. Other Burkholderia species are known to cause disease in plants and animals, and some are potential agents for biological warfare. Structural genomics efforts include characterizing the structures of enzymes from pathways that can be targeted for drug development. As part of these efforts, chorismate mutase from B. phymatum was produced and crystallized, and a 1.95 Å resolution structure is reported. This enzyme shares less than 33% sequence identity with other homologs of known structure. There are two classes of chorismate mutase: AroQ and AroH. The bacterial subclass AroQγ has reported roles in virulence. Chorismate mutase from B. phymatum has the prototypical AroQγ topology and retains the characteristic chorismate mutase active site. This suggests that substrate-based chorismate mutase inhibitors will not be specific and are likely to affect beneficial bacteria such as B. phymatum.

[Cloning and bioinformatics analysis of chorismate mutase gene from Salvia miltiorrhiza]

Chorismate mutase catalyzes the conversion of chorismate to prephenate that is the first committed step in the biosynthesis of the aromatic amino acids phenylalanine and tyrosine. A chorismate mutase gene, designated SmCM1, was isolated from Salvia miltiorrhiza by using RT-PCR. The full length of SmCM1 cDNA consists of 948 nucleotides and has an open reading frame of 765 bp. The deduced amino acid sequence of SmCM1 has 255 amino acid residues which forms a 36.0 kD polypeptide with calculated pI of 6.41 as expected. The putative polypeptide contains a CM_2 super family function domain. Blast W results showed that SmCM1 had 70% of the similarity with Petunia x hybrid CM, 72% of the similarity with Arabidopsis thaliana CM, and 64% of similarity with Populus trichocarpa CM. The transcription level of SmCM1 in root, stem and leaf was analysed by realtime quantitative PCR. The results showed the expression level of the SmCM1 in leaf was highest, and lowest in root. Yeast extract and silver ion joint induction could markedly stimulate the increase of mRNA expression of SmCM1 and its upstream 3-deoxy-7- phosphoheptulonate synthase (DAHPS) and chorismate synthase (CS). It was 7.9, 5.5 and 9.8 times of control on 8 h after induction, respectively.

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.

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.

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.

Consensus protein design without phylogenetic bias

Consensus design is an appealing strategy for the stabilization of proteins. It exploits amino acid conservation in sets of homologous proteins to identify likely beneficial mutations. Nevertheless, its success depends on the phylogenetic diversity of the sequence set available. Here, we show that randomization of a single protein represents a reliable alternative source of sequence diversity that is essentially free of phylogenetic bias. A small number of functional protein sequences selected from binary-patterned libraries suffice as input for the consensus design of active enzymes that are easier to produce and substantially more stable than individual members of the starting data set. Although catalytic activity correlates less consistently with sequence conservation in these extensively randomized proteins, less extreme mutagenesis strategies might be adopted in practice to augment stability while maintaining function.

A novel noncovalent complex of chorismate mutase and DAHP synthase from Mycobacterium tuberculosis: protein purification, crystallization and X-ray diffraction analysis

Chorismate mutase catalyzes a key step in the shikimate-biosynthetic pathway and hence is an essential enzyme in bacteria, plants and fungi. Mycobacterium tuberculosis contains two chorismate mutases, a secreted and an intracellular one, the latter of which (MtCM; Rv0948c; 90 amino-acid residues; 10 kDa) is the subject of this work. Here are reported the gene expression, purification and crystallization of MtCM alone and of its complex with another shikimate-pathway enzyme, DAHP synthase (MtDS; Rv2178c; 472 amino-acid residues; 52 kDa), which has been shown to enhance the catalytic efficiency of MtCM. The MtCM-MtDS complex represents the first noncovalent enzyme complex from the common shikimate pathway to be structurally characterized. Soaking experiments with a transition-state analogue are also reported. The crystals of MtCM and the MtCM-MtDS complex diffracted to 1.6 and 2.1 A resolution, respectively.

Structural and functional investigation of a secreted chorismate mutase from the plant-parasitic nematode Heterodera schachtii in the context of related enzymes from diverse origins

In this article, we present the cloning of Hscm1, a gene for chorismate mutase (CM) from the beet cyst nematode Heterodera schachtii. CM is a key branch-point enzyme of the shikimate pathway, and secondary metabolites that arise from this pathway control developmental programmes and defence responses of the plant. By manipulating the plant's endogenous shikimate pathway, the nematode can influence the plant physiology for its own benefit. Hscm1 is a member of the CM gene family and is expressed during the pre-parasitic and parasitic stages of the nematode's life cycle. In situ mRNA hybridization reveals an expression pattern specific to the subventral and dorsal pharyngeal glands. The predicted protein has a signal peptide for secretion in addition to two domains. The N-terminal domain of the mature protein, which is only found in cyst nematodes, contains six conserved cysteine residues, which may reflect the importance of disulphide bond formation for protein stabilization. The C-terminal domain holds a single catalytic site and has similarity to secreted CMs of pathogenic bacteria, classifying HsCM1 as an AroQ(gamma) enzyme. The presumed catalytic residues are discussed in detail, and genetic complementation experiments indicate that the C-terminal domain is essential for enzyme activity. Finally, we show how the modular design of the protein is mirrored in the genomic sequence by the intron/exon organization, suggesting exon shuffling as a mechanism for the evolutionary assembly of this protein.

On the relationship between folding and chemical landscapes in enzyme catalysis

Elucidating the relationship between the folding landscape of enzymes and their catalytic power has been one of the challenges of modern enzymology. The present work explores this issue by using a simplified folding model to generate the free-energy landscape of an enzyme and then to evaluate the activation barriers for the chemical step in different regions of the landscape. This approach is used to investigate the recent finding that an engineered monomeric chorismate mutase exhibits catalytic efficiency similar to the naturally occurring dimer even though it exhibits the properties of an intrinsically disordered molten globule. It is found that the monomer becomes more confined than its native-like counterpart upon ligand binding but still retains a wider catalytic region. Although the overall rate acceleration is still determined by reduction of the reorganization energy, the detailed contribution of different barriers yields a more complex picture for the chemical process than that of a single path. This work provides insight into the relationship between folding landscapes and catalysis. The computational approach used here may also provide a powerful strategy for modeling single-molecule experiments and designing enzymes.

Exploring SCC-DFTB paths for mapping QM/MM reaction mechanisms

A new first-order procedure for locating transition structures (TS) that employs hybrid quantum mechanical/molecular mechanical (QM/MM) potentials has been developed. This new technique (RPATh+RESD) combines the replica path method (RPATh) and standard reaction coordinate driving (RCD) techniques in an approach that both efficiently determines reaction barriers and successfully eliminates two key weaknesses of RCD calculations (i.e., hysteresis/discontinuities in the path and the sequential nature of the RCD procedure). In addition, we have extended CHARMM's QM/MM reaction pathway methods, the RPATh and nudged elastic band (NEB) methods, to incorporate SCC-DFTB wave functions. This newly added functionality has been applied to the chorismate mutase-catalyzed interconversion of chorismate to prephenate, which is a key step in the shikimate pathway of bacteria, fungi, and other higher plants. The RPATh+RESD barrier height (DeltaE=5.7 kcal/mol) is in good agreement with previous results from full-energy surface mapping studies (Zhang, X.; Zhang, X.; Bruice, T. C. Biochemistry 2005, 44, 10443-10448). Full reaction paths were independently mapped with RPATh and NEB methods and showed good agreement with the final transition state from the RPATh+RESD "gold standard" and previous high-level QM/MM transition states (Woodcock, H. L.; Hodoscek, M.; Gilbert, T. B.; Gill, P. M. W.; Schaefer, H. F.; Brooks, B. R. J. Comput. Chem. 2007, 28, 1485-1502). The SCC-DFTB TS geometry most closely approximates the MP2/6-31+G(d) QM/MM result. However, the barrier height is underestimated and possibly points to an area for improvement in SCC-DFTB parametrization. In addition, the steepest descents (SD) minimizer for the NEB method was modified to uncouple the in-path and off-path degrees of freedom during the minimization, which significantly improved performance. The convergence behavior of the RPATh and NEB was examined for SCC-DFTB wave functions, and it was determined that, in general, both methods converge at about the same rate, although the techniques used for convergence may be different. For instance, RPATh can effectively use the adopted basis Newton-Raphson (ABNR) minimizer, where NEB seems to require a combination of SD and ABNR.

Exhaustive Mutagenesis of Six Secondary Active-Site Residues in Escherichia coli Chorismate Mutase Shows the Importance of Hydrophobic Side Chains and a Helix N-Capping Position for Stability and Catalysis

Secondary active-site residues in enzymes, including hydrophobic amino acids, may contribute to catalysis through critical interactions that position the reacting molecule, organize hydrogen-bonding residues, and define the electrostatic environment of the active site. To ascertain the tolerance of an important model enzyme to mutation of active-site residues that do not directly hydrogen bond with the reacting molecule, all 19 possible amino acid substitutions were investigated in six positions of the engineered chorismate mutase domain of the Escherichia coli chorismate mutase-prephenate dehydratase. The six secondary active-site residues were selected to clarify results of a previous test of computational enzyme design procedures. Five of the positions encode hydrophobic side chains in the wild-type enzyme, and one forms a helix N-capping interaction as well as a salt bridge with a catalytically essential residue. Each mutant was evaluated for its ability to complement an auxotrophic chorismate mutase deletion strain. Kinetic parameters and thermal stabilities were measured for variants with in vivo activity. Altogether, we find that the enzyme tolerated 34% of the 114 possible substitutions, with a few mutations leading to increases in the catalytic efficiency of the enzyme. The results show the importance of secondary amino acid residues in determining enzymatic activity, and they point to strengths and weaknesses in current computational enzyme design procedures.

A Transition Path Sampling Study of the Reaction Catalyzed by the Enzyme Chorismate Mutase

The study of the chemical steps in enzyme-catalyzed reactions represents a challenge for molecular simulation techniques. One concern is how to calculate paths for the reaction. Common techniques include the definition of a reaction coordinate in terms of a small set of (normally) geometrical variables or the determination of minimum energy paths on the potential energy surface of the reacting system. Both have disadvantages, the former because it presupposes knowledge of which variables are likely to be important for reaction and the latter because it provides a static picture and dynamical effects are ignored. In this paper, we employ the transition path sampling method developed by Chandler and co-workers, which overcomes some of these limitations. The reaction that we have chosen is the chorismate-mutase-catalyzed conversion of chorismate into prephenate, which has become something of a test case for simulation studies of enzyme mechanisms. We generated an ensemble of approximately 1000 independent transition paths for the reaction in the enzyme and another approximately 500 for the corresponding reaction in solution. A large variety of analyses of these paths was performed, but we have concentrated on characterizing the transition state ensemble, particularly the flexibility of its structures with respect to other ligands of the enzyme and the time evolution of various geometrical and energetic properties as the reaction proceeds. We have also devised an approximate technique for locating transition state structures along the paths.

Biochemical and structural characterization of the secreted chorismate mutase (Rv1885c) from Mycobacterium tuberculosis H37Rv: an *AroQ enzyme not regulated by the aromatic amino acids

The gene Rv1885c from the genome of Mycobacterium tuberculosis H37Rv encodes a monofunctional and secreted chorismate mutase (*MtCM) with a 33-amino-acid cleavable signal sequence; hence, it belongs to the *AroQ class of chorismate mutases. Consistent with the heterologously expressed *MtCM having periplasmic destination in Escherichia coli and the absence of a discrete periplasmic compartment in M. tuberculosis, we show here that *MtCM secretes into the culture filtrate of M. tuberculosis. *MtCM functions as a homodimer and exhibits a dimeric state of the protein at a concentration as low as 5 nM. *MtCM exhibits simple Michaelis-Menten kinetics with a Km of 0.5 +/- 0.05 mM and a k(cat) of 60 s(-1) per active site (at 37 degrees C and pH 7.5). The crystal structure of *MtCM has been determined at 1.7 A resolution (Protein Data Bank identifier 2F6L). The protein has an all alpha-helical structure, and the active site is formed within a single chain without any contribution from the second chain in the dimer. Analysis of the structure shows a novel fold topology for the protein with a topologically rearranged helix containing Arg134. We provide evidence by site-directed mutagenesis that the residues Arg49, Lys60, Arg72, Thr105, Glu109, and Arg134 constitute the catalytic site; the numbering of the residues includes the signal sequence. Our investigation on the effect of phenylalanine, tyrosine, and tryptophan on *MtCM shows that *MtCM is not regulated by the aromatic amino acids. Consistent with this observation, the X-ray structure of *MtCM does not have an allosteric regulatory site.

Investigation of the first shot phenomenon in MALDI mass spectrometry of protein complexes

Matrix-assisted laser desorption/ionization-mass spectrometry (MALDI-MS) of noncovalent protein complexes is difficult, due to the disruptive nature of processes occurring during MALDI sample preparation and ion formation. Sometimes the observation of intact noncovalent protein complexes with MALDI is only possible if data are acquired from the first laser shot fired at a fresh sample; this is called the 'first shot phenomenon'. To study the origin of the first shot phenomenon, we used MALDI-MS and confocal laser scanning microscopy (CLSM) to examine typical MALDI sample preparations with embedded protein complexes, labeled with fluorophores. Fluorescence energy transfer techniques allowed the differentiation between intact and dissociated protein complexes with CLSM. In cases where a first shot behavior was observed by MALDI-MS, it was found to be accompanied by localization of protein complexes at the exterior of the sample crystals. Segregation of the large protein complexes to the exterior and dissociation of the complexes in the crystal interior during sample crystallization can rationalize this observation.

Relative tolerance of mesostable and thermostable protein homologs to extensive mutation

Evolvability, designability, and plasticity of a protein are properties that are important to protein engineers, but difficult to quantify. Here, we directly compare homologous AroQ chorismate mutases from the thermophile Methanococcus jannaschii and the mesophile Escherichia coli with respect to their capacity to accommodate extensive mutation. The N-terminal helix comprising about 40% of these proteins was randomized at the genetic level using a binary pattern of hydrophobic and hydrophilic residues based on the respective wild-type sequences. Catalytically active library members were identified by a survival-selection assay in a chorismate mutase-deficient E. coli strain. Functional variants were found approximately approximately 10-times more frequently with the thermostable protein compared to its mesostable counterpart. Moreover, detailed sequence analysis revealed that functional M. jannaschii enzyme variants contained a smaller number of conserved residues and tolerated greater variability at individual sequence positions. Our results thus highlight the greater robustness of the thermostable protein with respect to amino acid substitution, while identifying specific sites important for constructing active enzymes. Overall, they support the notion that redesign projects will benefit from using a thermostable starting structure, even at very high mutational loads.

Combinatorial methods for small-molecule placement in computational enzyme design

The incorporation of small-molecule transition state structures into protein design calculations poses special challenges because of the need to represent the added translational, rotational, and conformational freedoms within an already difficult optimization problem. Successful approaches to computational enzyme design have focused on catalytic side-chain contacts to guide placement of small molecules in active sites. We describe a process for modeling small molecules in enzyme design calculations that extends previously described methods, allowing favorable small-molecule positions and conformations to be explored simultaneously with sequence optimization. Because all current computational enzyme design methods rely heavily on sampling of possible active site geometries from discrete conformational states, we tested the effects of discretization parameters on calculation results. Rotational and translational step sizes as well as side-chain library types were varied in a series of computational tests designed to identify native-like binding contacts in three natural systems. We find that conformational parameters, especially the type of rotamer library used, significantly affect the ability of design calculations to recover native binding-site geometries. We describe the construction and use of a crystallographic conformer library and find that it more reliably captures active-site geometries than traditional rotamer libraries in the systems tested.

Hybrid Quantum Mechanics/Molecular Mechanics Simulations with Two-Dimensional Interpolated Corrections: Application to Enzymatic Processes

Hybrid quantum mechanics/molecular mechanics (QM/MM) techniques are widely used to study chemical reactions in large systems. Because of the computational cost associated with the high dimensionality of these systems, the quantum description is usually restricted to low-level methods, such as semiempirical Hamiltonians. In some cases, the description obtained at this computational level is quite poor and corrections must be considered. We here propose a simple but efficient way to include higher-level corrections to be used in potential energy surface explorations and in the calculation of potentials of mean force. We evaluate a correction energy term as the difference between a high-level and a low-level calculation on the QM subsystem, employing either the polarized or the gas-phase wave function, obtained as a function of two geometrical coordinates relevant in the process considered. Through the use of two-dimensional bicubic splines this correction energy is included in the simulations, ensuring the continuity and derivability of the energy function. We have tested the proposed scheme with two prototypical examples: the chorismate to prephenate rearrangement catalyzed by Bacillus subtilis chorismate mutase and the catechol methylation catalyzed by catechol O-methyltransferase. In both cases the use of interpolated corrections clearly improves the energetic and geometric descriptions of the reaction.

The 2.15 A crystal structure of Mycobacterium tuberculosis chorismate mutase reveals an unexpected gene duplication and suggests a role in host-pathogen interactions

Chorismate mutase catalyzes the first committed step toward the biosynthesis of the aromatic amino acids, phenylalanine and tyrosine. While this biosynthetic pathway exists exclusively in the cell cytoplasm, the Mycobacterium tuberculosis enzyme has been shown to be secreted into the extracellular medium. The secretory nature of the enzyme and its existence in M. tuberculosis as a duplicated gene are suggestive of its role in host-pathogen interactions. We report here the crystal structure of homodimeric chorismate mutase (Rv1885c) from M. tuberculosis determined at 2.15 A resolution. The structure suggests possible gene duplication within each subunit of the dimer (residues 35-119 and 130-199) and reveals an interesting proline-rich region on the protein surface (residues 119-130), which might act as a recognition site for protein-protein interactions. The structure also offers an explanation for its regulation by small ligands, such as tryptophan, a feature previously unknown in the prototypical Escherichia coli chorismate mutase. The tryptophan ligand is found to be sandwiched between the two monomers in a dimer contacting residues 66-68. The active site in the "gene-duplicated" monomer is occupied by a sulfate ion and is located in the first half of the polypeptide, unlike in the Saccharomyces cerevisiae (yeast) enzyme, where it is located in the later half. We hypothesize that the M. tuberculosis chorismate mutase might have a role to play in host-pathogen interactions, making it an important target for designing inhibitor molecules against the deadly pathogen.

Temperature Dependence of the Structure of the Substrate and Active Site of the Thermus thermophilus Chorismate Mutase E·S Complex

Molecular dynamics (MD) simulations of Thermus thermophilus chorismate mutase substrate complex (TtCM x S) have been carried out at 298 K, 333 K, and the temperature of optimum activity: 343 K. The enzyme exists as trimeric subunits with active sites shared between two neighboring subunits. Two features distinguish intersubunit linkages of the thermophilic and mesophilic enzyme Bacillus subtilis chorismate mutase substrate complex (BsCM x S): (i) electrostatic interactions by intersubunit ion pairs (Arg3-Glu40*/41, Arg76-Glu51* and Arg69*-Asp101, residues labeled with an asterisk are from the neighboring subunit) in the TtCM x S are not present in the structure of the BsCM x S; and (ii) replacement of polar residues with short and nonpolar residues in the interstices of the TtCM x S tighten the intersubunit hydrophobic interactions compared to BsCM x S. Concerning the active site, electrostatic interactions of the critically placed Arg6 and Arg63* with the two carboxylates of chorismate place the latter in a reactive conformation to spontaneously undergo a Claisen rearrangement. The optimum geometry at the active site has the CZ atoms of the two arginines 11 A apart. With a decrease in temperature, Arg63* moves toward Arg6 and the average conformation structure of chorismate moves further away from the reactive ground state conformation. This movement is due to the decrease in distance separating the electrostatic (in the main) and hydrophobic interacting pairs holding the two subunits together.

Preliminary X-ray crystallographic analysis of the secreted chorismate mutase from Mycobacterium tuberculosis: a tricky crystallization problem solved

Chorismate mutase catalyzes the conversion of chorismate to prephenate in the biosynthesis of the aromatic amino acids tyrosine and phenylalanine in bacteria, fungi and plants. Here, the crystallization of the unusual secreted chorismate mutase from Mycobacterium tuberculosis (encoded by Rv1885c), a 37.2 kDa dimeric protein belonging to the AroQ(gamma) subclass of mutases, is reported. Crystal optimization was non-trivial and is discussed in detail. To obtain crystals of sufficient quality, it was critical to initiate crystallization at higher precipitant concentration and then transfer the drops to lower precipitant concentrations within 5-15 min, in an adaptation of a previously described technique [Saridakis & Chayen (2000), Protein Sci. 9, 755-757]. As a result of the optimization, diffraction improved from 3.5 to 1.3 A resolution. The crystals belong to space group P2(1), with unit-cell parameters a = 42.6, b = 72.6, c = 62.0 angstroms, beta = 104.5 degrees. The asymmetric unit contains one biological dimer, with 167 amino acids per protomer. A soak with a transition-state analogue is also described.

The catalytic power of enzymes: Conformational selection or transition state stabilization?

The mechanism by which enzymes produce enormous rate enhancements in the reactions they catalyze remains unknown. Two viewpoints, selection of ground state conformations and stabilization of the transition state, are present in the literature in apparent opposition. To provide more insight into current discussion about enzyme efficiency, a two-state model of enzyme catalysis was developed. The model was designed to include both the pre-chemical (ground state conformations) and the chemical (transition state) components of the process for the substrate both in water and in the enzyme. Although the model is of general applicability, the chorismate to prephenate reaction catalyzed by chorismate mutase was chosen for illustrative purposes. The resulting kinetic equations show that the catalytic power of enzymes, quantified as the k(cat)/k(uncat) ratio, is the product of two terms: one including the equilibrium constants for the substrate conformational states and the other including the rate constants for the uncatalyzed and catalyzed chemical reactions. The model shows that these components are not mutually exclusive and can be simultaneously present in an enzymic system, being their relative contribution a property of the enzyme. The developed mathematical expressions reveal that the conformational and reaction components of the process perform differently for the translation of molecular efficiency (changes in energy levels) into observed enzymic efficiency (changes in k(cat)), being, in general, more productive the component involving the transition state.

Crystal structure of prephenate dehydrogenase from Aquifex aeolicus. Insights into the catalytic mechanism

The enzyme prephenate dehydrogenase catalyzes the oxidative decarboxylation of prephenate to 4-hydroxyphenylpyruvate for the biosynthesis of tyrosine. Prephenate dehydrogenases exist as either monofunctional or bifunctional enzymes. The bifunctional enzymes are diverse, since the prephenate dehydrogenase domain is associated with other enzymes, such as chorismate mutase and 3-phosphoskimate 1-carboxyvinyltransferase. We report the first crystal structure of a monofunctional prephenate dehydrogenase enzyme from the hyper-thermophile Aquifex aeolicus in complex with NAD+. This protein consists of two structural domains, a modified nucleotide-binding domain and a novel helical prephenate binding domain. The active site of prephenate dehydrogenase is formed at the domain interface and is shared between the subunits of the dimer. We infer from the structure that access to the active site is regulated via a gated mechanism, which is modulated by an ionic network involving a conserved arginine, Arg250. In addition, the crystal structure reveals for the first time the positions of a number of key catalytic residues and the identity of other active site residues that may participate in the reaction mechanism; these residues include Ser126 and Lys246 and the catalytic histidine, His147. Analysis of the structure further reveals that two secondary structure elements, beta3 and beta7, are missing in the prephenate dehydrogenase domain of the bifunctional chorismate mutase-prephenate dehydrogenase enzymes. This observation suggests that the two functional domains of chorismate mutase-prephenate dehydrogenase are interdependent and explains why these domains cannot be separated.

1.6 A crystal structure of the secreted chorismate mutase from Mycobacterium tuberculosis: novel fold topology revealed

The presence of exported chorismate mutases produced by certain organisms such as Mycobacterium tuberculosis has been shown to correlate with their pathogenicity. As such, these proteins comprise a new group of promising selective drug targets. Here, we report the high-resolution crystal structure of the secreted dimeric chorismate mutase from M. tuberculosis (*MtCM; encoded by Rv1885c), which represents the first 3D-structure of a member of this chorismate mutase family, termed the AroQ(gamma) subclass. Structures are presented both for the unliganded enzyme and for a complex with a transition state analog. The protomer fold resembles the structurally characterized (dimeric) Escherichia coli chorismate mutase domain, but exhibits a new topology, with helix H4 of *MtCM carrying the catalytic site residue missing in the shortened helix H1. Furthermore, the structure of each *MtCM protomer is significantly more compact and only harbors one active site pocket, which is formed entirely by one polypeptide chain. Apart from the structural model, we present evidence as to how the substrate may enter the active site.

Isolation and characterization of another cDNA encoding a chorismate mutase from the phytoparasitic nematode Meloidogyne arenaria

A new cDNA, named Ma-cm-2, encoding a chorismate mutase (CM), has been isolated from Meloidogyne arenaria. The full-length cDNA, carrying the trans-spliced SL1 leader sequence, was 753-bp long with an open reading frame of 576 bp. The deduced protein MA-CM-2 including amino-terminal signal peptide shows significant similarity to CMs of Meloidogyne incognita, Meloidogyne javanica, and also bacteria. Secondary structure prediction of MA-CM-2 indicates the presence of the three conserved alpha-helix domains present in the Escherichia coli CMs. Reverse transcription and polymerase chain reaction analysis showed that its transcript abundance is high in the early developmental stages and low in later ones. In situ mRNA hybridization revealed that the transcripts of Ma-cm-2 accumulated specifically in the two subventral oesophageal gland cells of M. arenaria. The widespread existence of CMs in the sedentary endoparasitic nematodes implicates that this enzyme plays an important role in the host-parasite interaction.

The Allosteric Mechanism of Yeast Chorismate Mutase: A Dynamic Analysis

The effector-regulated allosteric mechanism of yeast chorismate mutase (YCM) was studied by normal mode analysis and targeted molecular dynamics. The normal mode analysis shows that the conformational change between YCM in the R state and in the T state can be represented by a relatively small number of low-frequency modes. This suggests that the transition is coded in the structure and is likely to have a low energetic barrier. Quantitative comparisons (i.e. frequencies) between the low-frequency modes of YCM with and without effectors (modeled structures) reveal that the binding of Trp increases the global flexibility, whereas Tyr decreases global flexibility. The targeted molecular dynamics simulation of substrate analog release from the YCM active site suggests that a series of residues are critical for orienting and "recruiting" the substrate. The simulation led to the switching of a series of substrate-release-coupled salt-bridge partners in the ligand-binding domain; similar changes occur in the transition between YCM R-state and T-state crystal structures. Thus, the normal mode analysis and targeted molecular dynamics results provide evidence that the effectors regulate YCM activity by influencing the global flexibility. The change in flexibility is coupled to the binding of substrate to the T state and release of the product from the R state, respectively.

Experimental and Computational Investigation of the Uncatalyzed Rearrangement and Elimination Reactions of Isochorismate

The versatile biosynthetic intermediate isochorismate decomposes in aqueous buffer by two competitive pathways, one leading to isoprephenate by a facile Claisen rearrangement and the other to salicylate via elimination of the enolpyruvyl side chain. Computation suggests that both processes are concerted but asynchronous pericyclic reactions, with considerable C-O cleavage in the transition state but relatively little C-C bond formation (rearrangement) or hydrogen atom transfer to the enolpyruvyl side chain (elimination). Kinetic experiments show that rearrangement is roughly 8-times more favorable than elimination. Moreover, transfer of the C2 hydrogen atom to C9 was verified by monitoring the decomposition of [2-(2)H]isochorismate, which was prepared chemoenzymatically from labeled shikimate, by (2)H NMR spectroscopy and observing the appearance of [3-(2)H]pyruvate. Finally, the isotope effects obtained with the C2 deuterated substrate are in good agreement with calculations assuming pericyclic reaction mechanisms. These results provide a benchmark for mechanistic investigations of isochorismate mutase and isochorismate pyruvate lyase, the enzymes that respectively catalyze the rearrangement and elimination reactions in plants and bacteria.

All Electron Quantum Chemical Calculation of the Entire Enzyme System Confirms a Collective Catalytic Device in the Chorismate Mutase Reaction

To elucidate the catalytic power of enzymes, we analyzed the reaction profile of Claisen rearrangement of Bacillus subtilis chorismate mutase (BsCM) by all electron quantum chemical calculations using the fragment molecular orbital (FMO) method. To the best of our knowledge, this is the first report of ab initio-based quantum chemical calculations of the entire enzyme system, where we provide a detailed analysis of the catalytic factors that accomplish transition-state stabilization (TSS). FMO calculations deliver an ab initio-level estimate of the intermolecular interaction between the substrate and the amino acid residues of the enzyme. To clarify the catalytic role of Arg90, we calculated the reaction profile of the wild-type BsCM as well as Lys90 and Cit90 mutant BsCMs. Structural refinement and the reaction path determination were performed at the ab initio QM/MM level, and FMO calculations were applied to the QM/MM refined structures. Comparison between three types of reactions established two collective catalytic factors in the BsCM reaction: (1) the hydrogen bonds connecting the Glu78-Arg90-substrate cooperatively control the stability of TS relative to the ES complex and (2) the positive charge on Arg90 polarizes the substrate in the TS region to gain more electrostatic stabilization.

A definitive mechanism for chorismate mutase

In previous research presentations, we have described the important features of the chorismate --> prephenate reaction using molecular dynamics (MD) and thermodynamic integration studies. This investigation of the reaction in Escherichia coli and water involves QM/MM procedures (SCCDFTB/MM two-dimensional reaction coordinates to identify transition state structures in the water, enzyme, and gas phase followed by B3LYP/6-31+G* single-point computations which allow the determination of activation energies in water and in the E. coli enzyme). Computed activation energies of 11.3 kcal/mol in enzyme and 20.3 kcal/mol in water may be compared to the experimental values of 12.7 and 20.7 kcal/mol, respectively. The transition state structures in the gas phase, water, and enzyme are much the same. The transition states are characteristic of a concerted pericyclic rearrangement. The very small differences in the partial charges of O13 in NAC and TS support only a small preferential (10%) electrostatic stabilization of TS. The free energy of NAC formation in water exceeds that in enzyme by 8.5 kcal/mol, and it is this favored formation of NAC that provides the major kinetic advantage to the enzymatic reaction. These findings compare most favorably with those previous observations of this laboratory employing molecular dynamics and thermodynamic integrations. A definitive mechanism for the chorismate mutase enzymes is provided.

Feedback inhibition of chorismate mutase/prephenate dehydrogenase (TyrA) of Escherichia coli: generation and characterization of tyrosine-insensitive mutants

In order to get insights into the feedback regulation by tyrosine of the Escherichia coli chorismate mutase/prephenate dehydrogenase (CM/PDH), which is encoded by the tyrA gene, feedback-inhibition-resistant (fbr) mutants were generated by error-prone PCR. The tyrA(fbr) mutants were selected by virtue of their resistance toward m-fluoro-D,L-tyrosine, and seven representatives were characterized on the biochemical as well as on the molecular level. The PDH activities of the purified His6-tagged TyrA proteins exhibited up to 35% of the enzyme activity of TyrA(WT), but tyrosine did not inhibit the mutant PDH activities. On the other hand, CM activities of the TyrA(fbr) mutants were similar to those of the TyrA(WT) protein. Analyses of the DNA sequences of the tyrA genes revealed that tyrA(fbr) contained amino acid substitutions either at Tyr263 or at residues 354 to 357, indicating that these two sites are involved in the feedback inhibition by tyrosine.

The proficiency of a thermophilic chorismate mutase enzyme is solely through an entropic advantage in the enzyme reaction

A study of the Thermus thermophilus chorismate mutase (TtCM) is described by using quantum mechanics (self-consistent-charge density-functional tight binding)/molecular mechanics, umbrella sampling, and the weighted histogram analysis method. The computed free energies of activation for the reactions in water and TtCM are comparable to the experimental values. The free energies for formation of near attack conformer have been determined to be 8.06 and 0.05 kcal/mol in water and TtCM, respectively. The near attack conformer stabilization contributes approximately 90% to the proficiency of the enzymatic reaction compared with the reaction in water. The transition state (TS) structures and partial atom charges are much the same in the enzymatic and water reactions. The difference in the electrostatic interactions of Arg-89 with O13 in the enzyme-substrate complex and enzyme-TS complex provides the latter with but 0.55 kcal/mol of 1.92 kcal/mol total TS stabilization. Differences in electrostatic interactions between components at the active site in the enzyme-substrate complex and enzyme-TS complex are barely significant, such that TS stabilization is of minor importance and the enzymatic catalysis is through an entropic advantage.

An Active Enzyme Constructed from a 9-Amino Acid Alphabet

Nature employs a set of 20 amino acids to produce a repertoire of protein structures endowed with sophisticated functions. Here, we combined design and selection to create an enzyme composed entirely from a set of only 9 amino acids that can rescue auxotrophic cells lacking chorismate mutase. The simplified protein captures key structural features of its natural counterpart but appears to be somewhat less stable and more flexible. The potential of a dramatically reduced amino acid alphabet to produce an active catalyst supports the notion that primordial enzymes may have possessed low amino acid diversity and suggests that combinatorial engineering strategies, such as the one used here, may be generally applied to create enzymes with novel structures and functions.

Mechanistic insights into the isochorismate pyruvate lyase activity of the catalytically promiscuous PchB from combinatorial mutagenesis and selection

PchB from Pseudomonas aeruginosa possesses isochorismate pyruvate lyase (IPL) and weak chorismate mutase (CM) activity. Homology modeling based on a structurally characterized CM, coupled with randomization of presumed key active site residues (Arg54, Glu90, Gln91) and in vivo selection for CM activity, was used to derive mechanistic insights into the IPL activity of PchB. Mutation of Arg54 was incompatible with viability, and the CM and IPL activities of an engineered R54K variant were reduced 1,000-fold each. The observation that position 90 was tolerant to substitution but position 91 was essentially confined to Gln or Glu in functional variants rules out involvement of Glu90 in general base catalysis. Counter to the generally accepted mechanistic hypothesis for pyruvate lyases, we propose for PchB a rare [1,5]-sigmatropic reaction mechanism that invokes electrostatic catalysis in analogy to the [3,3]-pericyclic rearrangement of chorismate in CMs. A common catalytic principle for both PchB functions is also supported by the covariance of the catalytic parameters for the CM and IPL activities and the shared functional requirement for a protonated Glu91 in Q91E variants. The experiments demonstrate that focusing directed evolution strategies on the readily accessible surrogate activity of an enzyme can provide valuable insights into the mechanism of the primary reaction.

Simultaneous optimization of enzyme activity and quaternary structure by directed evolution

Natural evolution has produced efficient enzymes of enormous structural diversity. We imitated this natural process in the laboratory to augment the efficiency of an engineered chorismate mutase with low activity and an unusual hexameric topology. By applying two rounds of DNA shuffling and genetic selection, we obtained a 400-fold more efficient enzyme, containing three non-active-site mutations. Detailed biophysical characterization of the evolved variant suggests that it exists predominantly as a trimer in solution, but is otherwise similarly stable as the parent hexamer. The dramatic structural and functional effects achieved by a small number of seemingly innocuous substitutions highlights the utility of directed evolution for modifying protein-protein interactions to produce novel quaternary states with optimized activities.

Away from the edge II: in-house Se-SAS phasing with chromium radiation

Recently, the demands of high-throughput macromolecular crystallography have driven continuous improvements in phasing methods, data-collection protocols and many other technologies. Single-wavelength anomalous scattering (SAS) phasing with chromium X-ray radiation opens a new possibility for phasing a protein with data collected in-house and has led to several successful examples of de novo structure solution using only weak anomalous scatterers such as sulfur. To further reduce data-collection time and make SAS phasing more robust, it is natural to combine selenomethionine-derivatized protein (SeMet protein) with Cr Kalpha radiation to take advantage of the larger anomalous scattering signal from selenium (f'' = 2.28 e(-)) compared with sulfur (f'' = 1.14 e(-)). As reported herein, the crystal structure of a putative chorismate mutase from Clostridium thermocellum was determined using Se-SAS with Cr Kalpha radiation. Each protein molecule contains eight selenomethionine residues in 148 amino-acid residues, providing a calculated Bijvoet ratio of about 3.5% at the Cr Kalpha wavelength. A single data set to 2.2 A resolution with approximately ninefold redundancy was collected using an imaging-plate detector coupled with a Cr source. Structure solution, refinement and deposition to the Protein Data Bank were performed within 9 h of the availability of the scaled diffraction data. The procedure used here is applicable to many other proteins and promises to become a routine pathway for in-house high-throughput crystallography.

Crystallization and preliminary X-ray crystallographic studies of Mycobacterium tuberculosis chorismate mutase

Chorismate mutase catalyzes the first committed step in the biosynthesis of the aromatic amino acids phenylalanine and tyrosine in bacteria, fungi and higher plants. The recent re-annotation of the Mycobacterium tuberculosis genome has revealed the presence of a duplicate set of genes coding for chorismate mutase. The mycobacterial gene Rv1885c bears <20% sequence homology to other bacterial chorismate mutases, thus serving as a potential target for the development of inhibitors specific to the pathogen. The M. tuberculosis chorismate mutase was crystallized in space group C2 and the crystals diffracted to a resolution of 2.2 A. Matthews coefficient and self-rotation function calculations revealed the presence of two monomers in the asymmetric unit.

Multiple-Steering QM−MM Calculation of the Free Energy Profile in Chorismate Mutase

A novel technique for computing free energy profiles in enzymatic reactions using the multiple steering molecular dynamics approach in the context of an efficient QM-MM density functional scheme is presented. The conversion reaction of chorismate to prephenate catalyzed by the Bacillus subtilis enzyme chorismate mutase has been chosen as an illustrative example.

Investigation of Ligand Binding and Protein Dynamics in Bacillus subtilis Chorismate Mutase by Transverse Relaxation Optimized Spectroscopy−Nuclear Magnetic Resonance,

The structural and dynamical consequences of ligand binding to a monofunctional chorismate mutase from Bacillus subtilis have been investigated by solution NMR spectroscopy. TROSY methods were employed to assign 98% of the backbone (1)H(N), (1)H(alpha), (15)N, (13)C', and (13)C(alpha) resonances as well as 86% of the side chain (13)C resonances of the 44 kDa trimeric enzyme at 20 degrees C. This information was used to map chemical shift perturbations and changes in intramolecular mobility caused by binding of prephenate or a transition state analogue to the X-ray structure. Model-free interpretation of backbone dynamics for the free enzyme and its complexes based on (15)N relaxation data measured at 600 and 900 MHz showed significant structural consolidation of the protein in the presence of a bound ligand. In agreement with earlier structural and biochemical studies, substantial ordering of 10 otherwise highly flexible residues at the C-terminus is particularly notable. The observed changes suggest direct contact between this protein segment and the bound ligand, providing support for the proposal that the C-terminus can serve as a lid for the active site, limiting diffusion into and out of the pocket and possibly imposing conformational control over substrate once bound. Other regions of the protein that experience substantial ligand-induced changes also border the active site or lie along the subunit interfaces, indicating that the enzyme adapts dynamically to ligands by a sort of induced fit mechanism. It is believed that the mutase-catalyzed chorismate-to-prephenate rearrangement is partially encounter controlled, and backbone motions on the millisecond time scale, as seen here, may contribute to the reaction barrier.

Characterization of the secreted chorismate mutase from the pathogen Mycobacterium tuberculosis

The gene encompassing ORF Rv1885c with weak sequence similarity to AroQ chorismate mutases (CMs) was cloned from the genome of Mycobacterium tuberculosis and expressed in Escherichia coli. The gene product (*MtCM) complements a CM-deficient E. coli strain, but only if produced without the predicted N-terminal signal sequence typical of M. tuberculosis. The mature *MtCM, which was purified by exploiting its resistance to irreversible thermal denaturation, possesses high CM activity in vitro. The enzyme follows simple Michaelis-Menten kinetics, having a k(cat) of 50 s(-1) and a K(m) of 180 microM (at 30 degrees C and pH 7.5). *MtCM was shown to be a dimer by analytical ultracentrifugation and size-exclusion chromatography. Secondary-structure prediction and CD spectroscopy confirmed that *MtCM is a member of the all-alpha-helical AroQ class of CMs, but it seems to have a topologically rearranged AroQ fold. Because CMs are normally intracellular metabolic enzymes required for the biosynthesis of phenylalanine and tyrosine, the existence of an exported CM in Gram-positive M. tuberculosis is puzzling. The observation that homologs of *MtCM with a predicted export sequence are generally only present in parasitic or pathogenic organisms suggests that secreted CMs may have evolved to participate in some aspect of parasitism or pathogenesis yet to be unraveled.

Computing Kinetic Isotope Effects for Chorismate Mutase with High Accuracy. A New DFT/MM Strategy

A novel procedure has been applied to compute experimentally unobserved intrinsic kinetic isotope effects upon the rearrangement of chorismate to prephenate catalyzed by B. subtilis chorismate mutase. In this modified QM/MM approach, the "low-level" QM description of the quantum region is corrected during the optimization procedure by means of a "high-level" calculation in vacuo, keeping the QM-MM interaction contribution at a quantum "low-level". This allows computation of energies, gradients, and Hessians including the polarization of the QM subsystem and its interaction with the MM environment, both terms calculated using the low-level method at a reasonable computational cost. New information on an important enzymatic transformation is provided with greater reliability than has previously been possible. The predicted kinetic isotope effects on Vmax/Km are 1.33 and 0.86 (at 30 degrees C) for 5-3H and 9-3H2 substitutions, respectively, and 1.011 and 1.055 (at 22 degrees C) for 1-13C and 7-18O substitutions, respectively.

Allosteric mechanisms in ACT domain containing enzymes involved in amino acid metabolism

An important sequence motif identified by sequence analysis is shared by the ACT domain family, which has been found in a number of diverse proteins. Most of the proteins containing the ACT domain seem to be involved in amino acid and purine synthesis and are in many cases allosteric enzymes with complex regulation enforced by the binding of ligands. Here we explore the current understanding of the ACT domain function including its role as an allosteric module in a selected group of enzymes. We will further describe in more detail three of the proteins where some understanding is available on function and structure: i) the archetypical ACT domain protein E. coli 3PGDH, which catalyzes the first step in the biosynthesis of L-Ser, ii) the bifunctional chorismate mutase/prephenate dehydratase (P-protein) from E. coli, which catalyzes the first two steps in the biosynthesis of L-Phe, and iii) the mammalian aromatic amino acid hydroxylases, with special emphasis on phenylalanine hydroxylase, which catalyzes the first step in the catabolic degradation of L-phenylalanine (L-Phe). The ACT domain is commonly involved in the binding of a small regulatory molecule, such as the amino acids L-Ser and L-Phe in the case of 3PGDH and P-protein, respectively. On the other hand, for PAH, and probably for other enzymes, this domain appears to have been incorporated as a handy, flexible small module with the potential to provide allosteric regulation via transmission of finely tuned conformational changes, not necessarily initiated by regulatory ligand binding at the domain itself.

Differential Transition-State Stabilization in Enzyme Catalysis: Quantum Chemical Analysis of Interactions in the Chorismate Mutase Reaction and Prediction of the Optimal Catalytic Field

Chorismate mutase is a key model system in the development of theories of enzyme catalysis. To analyze the physical nature of catalytic interactions within the enzyme active site and to estimate the stabilization of the transition state (TS) relative to the substrate (differential transition state stabilization, DTSS), we have carried out nonempirical variation-perturbation analysis of the electrostatic, exchange, delocalization, and correlation interactions of the enzyme-bound substrate and transition-state structures derived from ab initio QM/MM modeling of Bacillus subtilis chorismate mutase. Significant TS stabilization by approximately -23 kcal/mol [MP2/6-31G(d)] relative to the bound substrate is in agreement with that of previous QM/MM modeling and contrasts with suggestions that catalysis by this enzyme arises purely from conformational selection effects. The most important contributions to DTSS come from the residues, Arg90, Arg7, Glu78, a crystallographic water molecule, Arg116, and Arg63, and are dominated by electrostatic effects. Analysis of the differential electrostatic potential of the TS and substrate allows calculation of the catalytic field, predicting the optimal location of charged groups to achieve maximal DTSS. Comparison with the active site of the enzyme from those of several species shows that the positions of charged active site residues correspond closely to the optimal catalytic field, showing that the enzyme has evolved specifically to stabilize the TS relative to the substrate.

An enzymatic molten globule: efficient coupling of folding and catalysis

A highly active, monomeric chorismate mutase, obtained by topological redesign of a dimeric helical bundle enzyme from Methanococcus jannaschii, was investigated by NMR and various other biochemical techniques, including H/D exchange. Although structural disorder is generally considered to be incompatible with efficient catalysis, the monomer, unlike its natural counterpart, unexpectedly possesses all of the characteristics of a molten globule. Global conformational ordering, observed upon binding of a transition state analog, indicates that folding can be coupled to catalysis with minimal energetic penalty. These results support the suggestion that many modern enzymes might have evolved from molten globule precursors. Insofar as their structural plasticity confers relaxed substrate specificity and/or catalytic promiscuity, molten globules may also be attractive starting points for the evolution of new catalysts in the laboratory.

Deciphering enzymes. Genetic selection as a probe of structure and mechanism

The efficient engineering of enzymes with novel activities remains an ongoing challenge. Towards this end, genetic selection techniques provide a method for finding rare solutions to catalytic problems that requires only a limited foreknowledge of structure-function relationships. We have used genetic selections to extensively probe the structure and mechanism of chorismate mutases. The insights gained from these investigations will aid future enzyme design efforts.

Structure of the OmpA-like domain of RmpM from Neisseria meningitidis

RmpM is a putative peptidoglycan binding protein from Neisseria meningitidis that has been shown to interact with integral outer membrane proteins such as porins and TonB-dependent transporters. Here we report the 1.9 A crystal structure of the C-terminal domain of RmpM. The 150-residue domain adopts a betaalphabetaalphabetabeta fold, as first identified in Bacillus subtilis chorismate mutase. The C-terminal RmpM domain is homologous to the periplasmic, C-terminal domain of Escherichia coli OmpA; these domains are thought to be responsible for non-covalent interactions with peptidoglycan. From the structure of the OmpA-like domain of RmpM, we suggest a putative peptidoglycan binding site and identify residues that may be essential for binding. Both the crystal structure and solution experiments indicate that RmpM may exist as a dimer. This would promote more efficient peptidoglycan binding, by allowing RmpM to interact simultaneously with two glycan chains through its C-terminal, OmpA-like binding domain, while its (structurally uncharacterized) N-terminal domain could stabilize oligomers of porins and TonB-dependent transporters in the outer membrane.

Chorismate mutase of Thermus thermophilus is a monofunctional AroH class enzyme inhibited by tyrosine

aroG, encoding the monofunctional chorismate mutase (TtCM) of the thermophilic gram-negative bacterium Thermus thermophilus, was cloned and its gene product characterized. TtCM was purified to homogeneity on an SDS polyacrylamide gel as a His-fusion protein with a deduced molecular mass of 15.8 kDa. The enzyme belongs to the rare group of AroH-type chorismate mutases which are mainly found in gram-positive bacteria of the Bacillus/ Clostridia group and have recently also been described for gram-negative organisms. The native molecular mass is consistent with a pseudo-alpha/beta barrel enzyme that is organized as a trimer. Comparison of the enzyme's structure with that of its mesophilic counterpart from Bacillus revealed an increase in hydrophilicity on the protein's surface, greater hydrophobicity in cavities within the protein, and greater restriction of conformational freedom, features that contribute to the thermal stability of this chorismate mutase. The kinetic data show Michaelis-Menten substrate saturation with a Km of 290 microM, and a kcat/ Km value of 180 s(-1) mM(-1). TtCM was inhibited by tyrosine with a Ki =34 microM, possibly in a competitive manner.