Mutagenesis Study of Active Site Residues in Chorismate Mutase from Bacillus subtilis

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.

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.

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.

Electrostatic transition state stabilization rather than reactant destabilization provides the chemical basis for efficient chorismate mutase catalysis

For more than half a century, transition state theory has provided a useful framework for understanding the origins of enzyme catalysis. As proposed by Pauling, enzymes accelerate chemical reactions by binding transition states tighter than substrates, thereby lowering the activation energy compared with that of the corresponding uncatalyzed process. This paradigm has been challenged for chorismate mutase (CM), a well-characterized metabolic enzyme that catalyzes the rearrangement of chorismate to prephenate. Calculations have predicted the decisive factor in CM catalysis to be ground state destabilization rather than transition state stabilization. Using X-ray crystallography, we show, in contrast, that a sluggish variant of Bacillus subtilis CM, in which a cationic active-site arginine was replaced by a neutral citrulline, is a poor catalyst even though it effectively preorganizes chorismate for the reaction. A series of high-resolution molecular snapshots of the reaction coordinate, including the apo enzyme, and complexes with substrate, transition state analog and product, demonstrate that an active site, which is only complementary in shape to a reactive substrate conformer, is insufficient for effective catalysis. Instead, as with other enzymes, electrostatic stabilization of the CM transition state appears to be crucial for achieving high reaction rates.

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.

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.

pH Dependence of catalysis by Pseudomonas aeruginosa isochorismate-pyruvate lyase: implications for transition state stabilization and the role of lysine 42

An isochorismate-pyruvate lyase with adventitious chorismate mutase activity from Pseudomonas aerugionsa (PchB) achieves catalysis of both pericyclic reactions in part by the stabilization of reactive conformations and in part by electrostatic transition-state stabilization. When the active site loop Lys42 is mutated to histidine, the enzyme develops a pH dependence corresponding to a loss of catalytic power upon deprotonation of the histidine. Structural data indicate that the change is not due to changes in active site architecture, but due to the difference in charge at this key site. With loss of the positive charge on the K42H side chain at high pH, the enzyme retains lyase activity at ∼100-fold lowered catalytic efficiency but loses detectable mutase activity. We propose that both substrate organization and electrostatic transition state stabilization contribute to catalysis. However, the dominant reaction path for catalysis is dependent on reaction conditions, which influence the electrostatic properties of the enzyme active site amino acid side chains.

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.

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Attractive noncovalent interactions in asymmetric catalysis: links between enzymes and small molecule catalysts

Catalysis by neutral, organic, small molecules capable of binding and activating substrates solely via noncovalent interactions--particularly H-bonding--has emerged as an important approach in organocatalysis. The mechanisms by which such small molecule catalysts induce high enantioselectivity may be quite different from those used by catalysts that rely on covalent interactions with substrates. Attractive noncovalent interactions are weaker, less distance dependent, less directional, and more affected by entropy than covalent interactions. However, the conformational constraint required for high stereoinduction may be achieved, in principle, if multiple noncovalent attractive interactions are operating in concert. This perspective will outline some recent efforts to elucidate the cooperative mechanisms responsible for stereoinduction in highly enantioselective reactions promoted by noncovalent catalysts.

Structure-function analyses of isochorismate-pyruvate lyase from Pseudomonas aeruginosa suggest differing catalytic mechanisms for the two pericyclic reactions of this bifunctional enzyme

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions in a single active site. PchB physiologically produces salicylate and pyruvate from isochorismate for ultimate incorporation of the salicylate into the siderophore pyochelin. PchB also produces prephenate from chorismate, most likely due to structural homology to the Escherchia coli chorismate mutase. The molecular basis of catalysis among enzymatic pericyclic reactions is a matter of debate, one view holding that catalysis may be derived from electrostatic transition state stabilization and the opposing view that catalysis is derived from the generation of a reactive substrate conformation. Mutant forms of PchB were generated by site-directed mutagenesis at the site (K42) hypothesized to be key for electrostatic transition state stabilization (K42A, K42Q, K42E, and K42H). The loop containing K42 is mobile, and a mutant to slow loop dynamics was also designed (A43P). Finally, a previously characterized mutation (I87T) was also produced. Circular dichroism was used to assess the overall effect on secondary structure as a result of the mutations, and X-ray crystallographic structures are reported for K42A with salicylate and pyruvate bound and for apo-I87T. The data illustrate that the active site architecture is maintained in K42A-PchB, which indicates that differences in activity are not caused by secondary structural changes or by differences in active site loop conformation but rather by the chemical nature of this key residue. In contrast, the I87T structure demonstrates considerable mobility, suggesting that loop dynamics and conformational plasticity may be important for efficient catalysis. Finally, the mutational effects on k(cat) provide evidence that the two activities of PchB are not covariant and that a single hypothesis may not provide a sufficient explanation for catalysis.

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.

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.

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.

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.

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.

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.

Directed evolution of protein enzymes using nonhomologous random recombination

We recently reported the development of nonhomologous random recombination (NRR) as a method for nucleic acid diversification and applied NRR to the evolution of DNA aptamers. Here, we describe a modified method, protein NRR, that enables proteins to access diversity previously difficult or impossible to generate. We investigated the structural plasticity of protein folds and the ability of helical motifs to function in different contexts by applying protein NRR and in vivo selection to the evolution of chorismate mutase (CM) enzymes. Functional CM mutants evolved using protein NRR contained many insertions, deletions, and rearrangements. The distribution of these changes was not random but clustered in certain regions of the protein. Topologically rearranged but functional enzymes also emerged from these studies, indicating that multiple connectivities can accommodate a functional CM active site and demonstrating the ability to generate new domain connectivities through protein NRR. Protein NRR was also used to randomly recombine CM and fumarase, an unrelated but also alpha-helical protein. Whereas the resulting library contained fumarase fragments in many contexts before functional selection, library members surviving selection for CM activity invariably contained a CM core with fumarase sequences found only at the termini or in one loop. These results imply that internal helical fragments cannot be swapped between these proteins without the loss of nearly all CM activity. Our findings suggest that protein NRR will be useful in probing the functional requirements of enzymes and in the creation of new protein topologies.

A Comparative Study of Claisen and Cope Rearrangements Catalyzed by Chorismate Mutase. An Insight into Enzymatic Efficiency: Transition State Stabilization or Substrate Preorganization?

In this work we present a detailed analysis of the activation free energies and averaged interactions for the Claisen and Cope rearrangements of chorismate and carbachorismate catalyzed by Bacillus subtilischorismate mutase (BsCM) using quantum mechanics/molecular mechanics (QM/MM) simulation methods. In gas phase, both reactions are described as concerted processes, with the activation free energy for carbachorismate being about 10-15 kcal mol(-)(1) larger than for chorismate, at the AM1 and B3LYP/6-31G levels. Aqueous solution and BsCM active site environments reduce the free energy barriers for both reactions, due to the fact that in these media the two carboxylate groups can be approached more easily than in the gas phase. The enzyme specifically reduces the activation free energy of the Claisen rearrangement about 3 kcal mol(-)(1) more than that for the Cope reaction. This result is due to a larger transition state stabilization associated to the formation of a hydrogen bond between Arg90 and the ether oxygen. When this oxygen atom is changed by a methylene group, the interaction is lost and Arg90 moves inside the active site establishing stronger interactions with one of the carboxylate groups. This fact yields a more intense rearrangement of the substrate structure. Comparing two reactions in the same enzyme, we have been able to obtain conclusions about the relative magnitude of the substrate preorganization and transition state stabilization effects. Transition state stabilization seems to be the dominant effect in this case.

Quantitative evaluation of noncovalent chorismate mutase-inhibitor binding by ESI-MS

Electrospray time-of-flight mass spectrometry was used to quantitatively determine the dissociation constant of chorismate mutase and a transition state analogue inhibitor. This system presents a fairly complex stoichiometry because the native protein is a homotrimer with three equal and independent substrate binding sites. We can detect the chorismate mutase trimer as well as chorismate mutase-inhibitor complexes by choosing appropriate conditions in the ESI source. To verify that the protein-inhibitor complexes are specific, titration experiments with different enzyme variants and different inhibitors were performed. A plot of the number of bound inhibitors versus added inhibitor concentration revealed saturation behavior with 3:1 (inhibitor:functional trimer) stoichiometry for the TSA. The soft ESI conditions, the relatively high protein mass of 43.5 kDa, and the low charge state (high m/z) result in broad peaks, a typical problem in analyzing noncovalent protein complexes. Due to the low molecular weight of the TSA (226 Da) the peaks of the free protein and the protein with one, two or three inhibitors bound cannot be clearly resolved. For data analysis, relative peak areas of the deconvoluted spectra of chorismate mutase-inhibitor complexes were obtained by fitting appropriate peak shapes to the signals corresponding to the free enzyme and its complexes with one, two, or three inhibitor molecules. From the relative peak areas we were able to calculate a dissociation constant that agreed well with known solution-phase data. This method may be generally useful for interpreting mass spectra of noncovalent complexes that exhibit broad peaks in the high m/z range.

Apparent NAC effect in chorismate mutase reflects electrostatic transition state stabilization

The catalytic reaction of chorismate mutase (CM) has been the subject of major current attention. Nevertheless, the origin of the catalytic power of CM remains an open question. In particular, it has not been clear whether the enzyme works by providing electrostatic transition state stabilization (TSS), by applying steric strain, or by populating near attack conformation (NAC). The present work explores this issue by a systematic quantitative analysis. The overall catalytic effect is reproduced by the empirical valence bond (EVB) method. In addition, the binding free energy of the ground state and the transition state is evaluated, demonstrating that the enzyme works by TSS. Furthermore, the evaluation of the electrostatic contribution to the reduction of the activation energy establishes that the TSS results from electrostatic effects. It is also found that the apparent NAC effect is not the reason for the catalytic effect but the result of the TSS. It is concluded that in CM as in other enzymes the key catalytic effect is electrostatic TSS. However, since the charge distribution of the transition state and the reactant state is similar, the stabilization of the transition state leads to reduction in the distance between the reacting atoms in the reactant state.

Selective stabilization of the chorismate mutase transition state by a positively charged hydrogen bond donor

Citrulline was incorporated via chemical semisynthesis at position 90 in the active site of the AroH chorismate mutase from Bacillus subtilis. The wild-type arginine at this position makes hydrogen-bonding interactions with the ether oxygen of chorismate. Replacement of the positively charged guanidinium group with the isosteric but neutral urea has a dramatic effect on the ability of the enzyme to convert chorismate into prephenate. The Arg90Cit variant exhibits a >104-fold decrease in the catalytic rate constant kcat with a 2.7-fold increase in the Michaelis constant Km. In contrast, its affinity for a conformationally constrained inhibitor molecule that effectively mimics the geometry but not the dissociative character of the transition state is only reduced by a factor of approximately 6. These results show that an active site merely complementary to the reactive conformation of chorismate is insufficient for catalysis of the mutase reaction. Instead, electrostatic stabilization of the polarized transition state by provision of a cationic hydrogen bond donor proximal to the oxygen in the breaking C-O bond is essential for high catalytic efficiency.

Investigation of the enzymatic and nonenzymatic cope rearrangement of carbaprephenate to carbachorismate

The dimethyl esters of carbaprephenate and 4-epi-carbaprephenate were prepared by modification of published procedures. In methanol these compounds are converted quantitatively to isomeric 6-hydroxytricyclo[3.3.1.0(2,7)]non-3-en-1,3-dimethyl esters via a two-step sequence involving an initial Cope rearrangement, followed by intramolecular Diels-Alder reaction of the dimethyl carbachorismate or 4-epi-carbachorismate intermediates. Carbaprephenate and its epimer were obtained by alkaline hydrolysis of the corresponding dimethyl esters. These compounds, in contrast to their ester precursors, undergo spontaneous acid-catalyzed decarboxylation in aqueous solution. Only at high pH does the Cope rearrangement compete with decarboxylation. At pH 12 and 90 degrees C, carbaprephenate slowly rearranges to carbachorismate, which rapidly loses water to give 3-(2-carboxyallyl)benzoic acid as the major product. A small amount of the intramolecular Diels-Alder adduct derived from carbachorismate is also observed by NMR as a minor product. Carbaprephenate is not a substrate for the enzyme chorismate mutase from Bacillus subtilis (BsCM), nor does carbaprephenate inhibit the normal chorismate mutase activity of this enzyme, even when present in 200-fold excess over chorismate. Its low affinity for the enzyme-active site is presumably a consequence of placing a methylene group rather than an oxygen atom proximal to the essential cationic residue Arg90. Nevertheless, BsCM variants that lack this cation (R90G and R90A) do not accelerate the Cope rearrangement of carbaprephenate either, and a catalytic antibody 1F7, which exhibits modest chorismate mutase activity, is similarly inactive. Poor substrate binding and the relatively high barrier for the Cope compared to the Claisen rearrangement presumably account for the lack of detectable catalysis. Acceleration of this sigmatropic rearrangement apparently requires more than an active site that is complementary in shape to the reactive substrate conformer.

How much do enzymes really gain by restraining their reacting fragments?

The steric effect, exerted by enzymes on their reacting substrates, has been considered as a major factor in enzyme catalysis. In particular, it has been proposed that enzymes catalyze their reactions by pushing their reacting fragments to a catalytic configuration which is sometimes called near attack configuration (NAC). This work uses computer simulation approaches to determine the relative importance of the steric contribution to enzyme catalysis. The steric proposal is expressed in terms of well defined thermodynamic cycles that compare the reaction in the enzyme to the corresponding reaction in water. The S(N)2 reaction of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10, which was used in previous studies to support the strain concept is chosen as a test case for this proposal. The empirical valence bond (EVB) method provides the reaction potential surfaces in our studies. The reliability and efficiency of this method make it possible to obtain stable results for the steric free energy. Two independent strategies are used to evaluate the actual magnitude of the steric effect. The first applies restraints on the substrate coordinates in water in a way that mimics the steric effect of the protein active site. These restraints are then released and the free energy associated with the release process provides the desired estimate of the steric effect. The second approach eliminates the electrostatic interactions between the substrate and the surrounding in the enzyme and in water, and compares the corresponding reaction profiles. The difference between the resulting profiles provides a direct estimate of the nonelectrostatic contribution to catalysis and the corresponding steric effect. It is found that the nonelectrostatic contribution is about -0.7 kcal/mol while the full "apparent steric contribution" is about -2.2 kcal/mol. The apparent steric effect includes about -1.5 kcal/mol electrostatic contribution. The total electrostatic contribution is found to account for almost all the observed catalytic effect ( approximately -6.1 kcal/mol of the -6.8 calculated total catalytic effect). Thus, it is concluded that the steric effect is not the major source of the catalytic power of haloalkane dehalogenase. Furthermore, it is found that the largest component of the apparent steric effect is associated with the solvent reorganization energy. This solvent-induced effect is quite different from the traditional picture of balance between the repulsive interaction of the reactive fragments and the steric force of the protein.

Investigating and Engineering Enzymes by Genetic Selection

Natural enzymes have arisen over millions of years by the gradual process of Darwinian evolution. The fundamental steps of evolution-mutation, selection, and amplification-can also be exploited in the laboratory to create and characterize protein catalysts on a human timescale. In vivo genetic selection strategies enable the exhaustive analysis of protein libraries with 10(10) different members, and even larger ensembles can be studied with in vitro methods. Evolutionary approaches can consequently yield statistically meaningful insight into the complex and often subtle interactions that influence protein folding, structure, and catalytic mechanism. Such methods are also being used increasingly as an adjunct to design, thus providing access to novel proteins with tailored catalytic activities and selectivities.

Substrate conformational transitions in the active site of chorismate mutase: their role in the catalytic mechanism

Chorismate mutase acts at the first branch-point of aromatic amino acid biosynthesis and catalyzes the conversion of chorismate to prephenate. The results of molecular dynamics simulations of the substrate in solution and in the active site of chorismate mutase are reported. Two nonreactive conformers of chorismate are found to be more stable than the reactive pseudodiaxial chair conformer in solution. It is shown by QM/MM molecular dynamics simulations, which take into account the motions of the enzyme, that when these inactive conformers are bound to the active site, they are rapidly converted to the reactive chair conformer. This result suggests that one contribution of the enzyme is to bind the more prevalent nonreactive conformers and transform them into the active form in a step before the chemical reaction. The motion of the reactive chair conformer in the active site calculated by using the QM/MM potential generates transient structures that are closer to the transition state than is the stable CHAIR conformer.

A strategically positioned cation is crucial for efficient catalysis by chorismate mutase

Combinatorial mutagenesis and in vivo selection experiments previously afforded functional variants of the AroH class Bacillus subtilis chorismate mutase lacking the otherwise highly conserved active site residue Arg(90). Here, we present a detailed kinetic and crystallographic study of several such variants. Removing the arginine side chain (R90G and R90A) reduced catalytic efficiency by more than 5 orders of magnitude. Reintroducing a positive charge to the active site through lysine substitutions restored more than a factor of a thousand in k(cat). Remarkably, the lysine could be placed at position 90 or at the more remote position 88 provided a sterically suitable residue was present at the partner site. Crystal structures of the double mutants C88S/R90K and C88K/R90S show that the lysine adopts an extended conformation that would place its epsilon-ammonium group within hydrogen-bonding distance of the ether oxygen of bound chorismate in the transition state. These results provide support for the hypothesis that developing negative charge in the highly polarized transition state is stabilized electrostatically by a strategically placed cation. The implications of this finding for the mechanism of all natural chorismate mutases and for the design of artificial catalysts are discussed.

Probing the role of the C-terminus of Bacillus subtilis chorismate mutase by a novel random protein-termination strategy

A novel strategy combining random protein truncation and genetic selection has been developed to identify dispensable C-terminal segments of an enzyme. This approach, which entails the random introduction of premature termination codons, was applied to the last 17 residues of chorismate mutase from Bacillus subtilis (BsCM). Although structurally ill-defined, the C-terminus of BsCM has been proposed to cap the active site upon substrate binding and affect catalysis. However, sequence patterns of 178 selected gene variants show that the final 11 residues of the protein can be mutated and even removed without significantly impairing activity in vivo. In fact, none of the randomized residues is absolutely required, but a preference for wild-type Lys111, Ala112, Leu115, and Arg116 is apparent. These residues are part of a C-terminal 3(10)-helix and provide contacts with the rest of the protein or its ligands. The kinetic parameters of selected enzyme variants show that truncations and mutations do not significantly impair catalytic turnover (k(cat)) but substantially decrease k(cat)/K(m). Thus, while the 17 C-terminal residues of BsCM do not participate directly in the chemical rearrangement, they appear to contribute to enzymatic efficiency via uniform binding of the substrate and transition state.

Effective fragment potentials and the enzyme active site

Optimization of the binding conformation of a substrate in an enzyme active site using ab initio quantum chemistry methods are intractable since the active site comprises several hundred atoms. However, the active site can be decomposed into an active and spectator region where the spectator residues are represented by effective fragment potentials and reducing the number of all-electron atoms involved in the chemistry to a reasonable level. The effective fragment potentials for electrostatics and polarization are implemented in GAMESS but the repulsive and charge transfer potentials are fit to interaction energies of water with models of the residues. These repulsive/charge transfer potentials are generated for the protein residues and the EFP are then used to optimize binding of a transition state analogue to chorismate mutase (B. subtilis) and small dianions to ribonuclease A. For chorismate mutase the calculated binding conformation compares well to the comparable X-ray structure. The binding of the inhibitor to the glutamate/glutamine mutant active site is then predicted with the optimization including the glutamine residue constrained only at the C alpha atom. The binding conformations suggest important roles for tyr108 and arg63, which have not been noted earlier. The electrostatic stabilization of the transition state by the active site charge distribution has to be augmented by a specific electronic activation by glu78. In ribonuclease A, the protons are found to move to provide a clustering of the charges to bind the small dianions, phosphate, thiophosphate, and sulfate.

Bacillus subtilis chorismate mutase is partially diffusion-controlled

The effect of viscosogens on the enzyme-catalyzed rearrangement of chorismate to prephenate has been studied. The steady-state parameters kcat and kcat/Km for the monofunctional chorismate mutase from Bacillus subtilis (BsCM) decreased significantly with increasing concentrations of glycerol, whereas the 'sluggish' BsCM mutants C75A and C75S were insensitive to changes in microviscosity. The latter results rule out extraneous interactions of the viscosogen as an explanation for the effects observed with the wild-type enzyme. Additional control experiments show that neither viscosogen-induced shifts in the pH-dependence of the enzyme-catalyzed reaction nor small perturbations of the conformational equilibrium of chorismate can account for the observed effects. Instead, BsCM appears to be limited by substrate binding and product release at low and high substrate concentrations, respectively. Analysis of the kinetic data indicates that diffusive transition states are between 30 and 40% rate-determining in these concentration regimes; the chemical step must contribute to the remaining kinetic barrier. The relatively low value of the 'on' rates for chorismate and prephenate (approximately 2 x 106 m-1.s-1) probably reflects the need for a rare conformation of the enzyme, the ligand, or both for successful binding. Interestingly, the chorismate mutase domain of the bifunctional chorismate mutase-prephenate dehydratase from Escherichia coli, which has steady-state kinetic parameters comparable to those of BsCM but has a much less accessible active site, is insensitive to changes in viscosity and the reaction it catalyses is not diffusion-controlled.

Use of site-directed mutagenesis to identify residues specific for each reaction catalyzed by chorismate mutase-prephenate dehydrogenase from Escherichia coli

Site-directed mutagenesis was performed on the bifunctional enzyme chorismate mutase-prephenate dehydrogenase in order to identify groups important for each of the two reactions. We selected two residues for mutagenesis, Lys37 and His131, identified previously by differential peptide mapping to be essential for activity [Christendat, D., and Turnbull, J. (1996) Biochemistry 35, 4468-4479]. Kinetic studies reveal that K37Q exhibits no mutase activity while retaining wild-type dehydrogenase activity, verifying that Lys37 plays a key role in the mutase. By contrast His131 is not critical for the dehydrogenase; H131A is a reasonably efficient catalyst exhibiting 10% dehydrogenase and 30% mutase activity compared to the wild-type enzyme. Chemical modification of H131A by diethyl pyrocarbonate further inactivated the dehydrogenase, suggesting that a different histidine is now accessible to modification. To identify this group, the protein's remaining eight histidines were changed to alanine or asparagine. A single substitution, H197N, decreased the dehydrogenase activity by 5 orders of magnitude while full mutase activity was retained. In H197N, the Michaelis constants for prephenate and NAD+ and the mutant's elution profile from Sepharose-AMP were similar to those of wild-type enzyme, indicating that catalysis rather than substrate binding is altered. Log V for the dehydrogenase reaction catalyzed by H197N is pH-independent and is in contrast to wild-type enzyme, which shows a decrease in activity at low pH and pK of about 6.5. We conclude that His197 is an essential catalytic residue in the dehydrogenase reaction.

Thermodynamics of a transition state analogue inhibitor binding to Escherichia coli chorismate mutase: probing the charge state of an active site residue and its role in inhibitor binding and catalysis

Electrostatic interactions play important roles in the catalysis of chorismate to prephenate by chorismate mutase. Mutation of Gln88 to glutamate in the monofunctional chorismate mutase from Escherichia coli results in an enzyme with a pH profile of activity significantly different from that of the wild type protein. To investigate whether the mutation alters the substrate binding process or the catalysis, we have directly determined the thermodynamic parameters of a transition state analogue inhibitor binding to the wild-type chorismate mutase and its Q88E mutant using isothermal titration calorimetry. The results demonstrate that solvent reorganization and hydrophobic interactions contribute the predominant free energy to inhibitor binding. The charge state of Glu88 in the Q88E mutant was experimentally determined and was shown to be protonated at pH 4.5 and ionized at pH 7.8, consistent with earlier hypotheses. Most surprisingly, inhibitor binding energetics do not exhibit significant pH dependency for both enzymes. Our findings indicate that the charge state of Glu88 has a small impact on inhibitor binding but plays an important role in the catalytic process.

Enzymatic activity and partial purification of solanapyrone synthase: first enzyme catalyzing Diels-Alder reaction

In cell-free extracts of Alternaria solani, an enzymatic activity converting prosolanapyrone II to solanapyrones A and D via oxidation and subsequent Diels-Alder reaction has been found. Chromatography with DEAE-Sepharose provided two active fractions, pools 1 and 2. The former fraction converted prosolanapyrone II to solanapyrones A and D in a ratio of 2.2:1 with optical purities of 99% and 45% ee, respectively. The latter fraction did so in a ratio of 7.6:1 with 99% and nearly 0% ee, respectively. The enzyme partially purified from pool 2 native molecular weight of 40-62 kD and a pl of 4.25. The high reactivity of prosolanapyrone III in aqueous solution and the chromatographic behavior of the enzyme in pool 2 suggest that a single enzyme catalyzes both the oxidation and Diels-Alder reaction.