Yeast chorismate mutase in the R state: simulations of the active site

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.

Mechanism of histone methylation catalyzed by protein lysine methyltransferase SET7/9 and origin of product specificity

Methylation of certain lysine residues in the N-terminal tails of core histone proteins in nucleosome is of fundamental importance in the regulation of chromatin structure and gene expression. Such histone modification is catalyzed by protein lysine methyltransferases (PKMTs). PKMTs contain a conserved SET domain in almost all of the cases and may transfer one to three methyl groups from S-adenosyl-L-methionine (AdoMet) to the epsilon-amino group of the target lysine residue. Here, quantum mechanical/molecular mechanical molecular dynamics and free-energy simulations are performed on human PKMT SET7/9 and its mutants to understand two outstanding questions for the reaction catalyzed by PKMTs: the mechanism for deprotonation of positively charged methyl lysine (lysine) and origin of product specificity. The results of the simulations suggest that Tyr-335 (an absolute conserved residue in PKMTs) may play the role as the general base for the deprotonation after dissociation of AdoHcy (S-adenosyl-L-homocysteine) and before binding of AdoMet. It is shown that conformational changes could bring Y335 to the target methyl lysine (lysine) for proton abstraction. This mechanism provides an explanation why methyl transfers could be catalyzed by PKMTs processively. The free-energy profiles for methyl transfers are reported and analyzed for wild type and certain mutants (Y305F and Y335F) and the active-site interactions that are of importance for the enzyme's function are discussed. The results of the simulations provide important insights into the catalytic process and lead to a better understanding of experimental observations concerning the origin of product specificity for PKMTs.

Using reaction mechanism to measure enzyme similarity

The concept of reaction similarity has been well studied in terms of the overall transformation associated with a reaction, but not in terms of mechanism. We present the first method to give a quantitative measure of the similarity of reactions based upon their explicit mechanisms. Two approaches are presented to measure the similarity between individual steps of mechanisms: a fingerprint-based approach that incorporates relevant information on each mechanistic step; and an approach based only on bond formation, cleavage and changes in order. The overall similarity for two reaction mechanisms is then calculated using the Needleman-Wunsch alignment algorithm. An analysis of MACiE, a database of enzyme mechanisms, using our measure of similarity identifies some examples of convergent evolution of chemical mechanisms. In many cases, mechanism similarity is not reflected by similarity according to the EC system of enzyme classification. In particular, little mechanistic information is conveyed by the class level of the EC system.

A theoretical analysis of rate constants and kinetic isotope effects corresponding to different reactant valleys in lactate dehydrogenase

In some enzymatic systems large conformational changes are coupled to the chemical step, in such a way that dispersion of rate constants can be observed in single-molecule experiments, each corresponding to the reaction from a different reactant valley. Under this perspective here we present a computational study of pyruvate to lactate transformation catalyzed by lactate dehydrogenase. The reaction consists of a hydride transfer and a proton transfer that seem to take place concertedly. The degree of asynchronicity and the energy barrier depend on the particular starting reactant valley. In order to estimate rate constants we used a free energy perturbation technique adapted to follow the intrinsic reaction coordinate for several possible reaction paths. Tunneling effects are also obtained with a slightly modified version of the ensemble-averaged variational transition state theory with multidimensional tunneling contributions. According to our results the closure of the active site by means of a flexible loop can lead to the formation of different reactant complexes displaying different features in the disposition of some key residues (such as Arg109), interactions with the substrate and number of water molecules in the active site. The chemical step of the reaction takes place with a different reaction rate from each of these complexes. Finally, primary kinetic isotope effects for replacement of the transferring hydrogen of the cofactor with a deuteride are in good agreement with experimental observations, thus validating our methodology.

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 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.

Molecular dynamics analysis of a buckyball-antibody complex

This is a multinanosecond molecular dynamics study of a bio-nano complex formed by a carbon nanoparticle, a buckyball C(60), and a biological molecule, an antibody, with high binding affinity and specificity. In the simulation, the ball is completely desolvated by the binding site of the antibody by means of a nearly perfect shape complementarity and extensive side-chain interactions, with the exception that about 17% of the surface is persistently exposed to solvent and could be used for functional derivatization. The interactions are predominantly hydrophobic, but significant polar interactions occur as well. There exists a rich body of various pi-stacking interactions. Aromatic side chains are involved in both double and triple stackings with the ball. Some ionic side chains, such as the guanidinium group of arginine, also form pi-stackings with the ball. The results suggest that pi-stackings are very efficient and common modes of biological recognition of pi-electron-rich carbon nanoparticles. Most importantly, the results demonstrate that, in general, an ordinary protein binding site, such as that of an antibody, can readily bind to a carbon nanoparticle with high affinity and specificity through recognition modes that are common in protein-ligand recognition.

Interactions of Streptomyces griseus aminopeptidase with amino acid reaction products and their implications toward a catalytic mechanism

Streptomyces griseus aminopeptidase (SGAP) is a double-zinc exopeptidase with a high preference toward large hydrophobic amino-terminus residues. It is a monomer of a relatively low molecular weight (30 kDa), it is heat stable, it displays a high and efficient catalytic turnover, and its activity is modulated by calcium ions. The small size, high activity, and heat stability make SGAP a very attractive enzyme for various biotechnological applications, among which is the processing of recombinant DNA proteins and fusion protein products. Several free amino acids, such as phenylalanine, leucine, and methionine, were found to act as weak inhibitors of SGAP and hence were chosen for structural studies. These inhibitors can potentially be regarded as product analogs because one of the products obtained in a normal enzymatic reaction is the cleaved amino terminal amino acid of the substrate. The current study includes the X-ray crystallographic analysis of the SGAP complexes with methionine (1.53 A resolution), leucine (1.70 A resolution), and phenylalanine (1.80 A resolution). These three high-resolution structures have been used to fully characterize the SGAP active site and to identify some of the functional groups of the enzyme that are involved in enzyme-substrate and enzyme-product interactions. A unique binding site for the terminal amine group of the substrate (including the side chains of Glu131 and Asp160, as well as the carbonyl group of Arg202) is indicated to play an important role in the binding and orientation of both the substrate and the product of the catalytic reaction. These studies also suggest that Glu131 and Tyr246 are directly involved in the catalytic mechanism of the enzyme. Both of these residues seem to be important for substrate binding and orientation, as well as the stabilization of the tetrahedral transition state of the enzyme-substrate complex. Glu131 is specifically suggested to function as a general base during catalysis by promoting the nucleophilic attack of the zinc-bound water/hydroxide on the substrate carbonyl carbon. The structures of the three SGAP complexes are compared with recent structures of three related aminopeptidases: Aeromonas proteolytica aminopeptidase (AAP), leucine aminopeptidase (LAP), and methionine aminopeptidase (MAP) and their complexes with corresponding inhibitors and analogs. These structural results have been used for the simulation of several species along the reaction coordinate and for the suggestion of a general scheme for the proteolytic reaction catalyzed by SGAP.

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.

HARO7 encodes chorismate mutase of the methylotrophic yeast Hansenula polymorpha and is derepressed upon methanol utilization

The HARO7 gene of the methylotrophic, thermotolerant yeast Hansenula polymorpha was cloned by functional complementation. HARO7 encodes a monofunctional 280-amino-acid protein with chorismate mutase (EC 5.4. 99.5) activity that catalyzes the conversion of chorismate to prephenate, a key step in the biosynthesis of aromatic amino acids. The HARO7 gene product shows strong similarities to primary sequences of known eukaryotic chorismate mutase enzymes. After homologous overexpression and purification of the 32-kDa protein, its kinetic parameters (k(cat) = 319.1 s(-1), n(H) = 1.56, [S](0.5) = 16.7 mM) as well as its allosteric regulatory properties were determined. Tryptophan acts as heterotropic positive effector; tyrosine is a negative-acting, heterotropic feedback inhibitor of enzyme activity. The influence of temperature on catalytic turnover and the thermal stability of the enzyme were determined and compared to features of the chorismate mutase enzyme of Saccharomyces cerevisiae. Using the Cre-loxP recombination system, we constructed mutant strains carrying a disrupted HARO7 gene that showed tyrosine auxotrophy and severe growth defects. The amount of the 0.9-kb HARO7 mRNA is independent of amino acid starvation conditions but increases twofold in the presence of methanol as the sole carbon source, implying a catabolite repression system acting on HARO7 expression.

The aroC gene of Aspergillus nidulans codes for a monofunctional, allosterically regulated chorismate mutase

The cDNA and the chromosomal locus of the aroC gene of Aspergillus nidulans were cloned and is the first representative of a filamentous fungal gene encoding chorismate mutase (EC 5.4.99.5), the enzyme at the first branch point of aromatic amino acid biosynthesis. The aroC gene complements the Saccharomyces cerevisiae aro7Delta as well as the A. nidulans aroC mutation. The gene consists of three exons interrupted by two short intron sequences. The expressed mRNA is 0.96 kilobases in length and aroC expression is not regulated on the transcriptional level under amino acid starvation conditions. aroC encodes a monofunctional polypeptide of 268 amino acids. Purification of this 30-kDa enzyme allowed determination of its kinetic parameters (k(cat) = 82 s(-1), n(H) = 1. 56, [S](0.5) = 2.3 mM), varying pH dependence of catalytic activity in different regulatory states, and an acidic pI value of 4.7. Tryptophan acts as heterotropic activator and tyrosine as negative acting, heterotropic feedback-inhibitor with a K(i) of 2.8 microM. Immunological data, homology modeling, as well as electron microscopy studies, indicate that this chorismate mutase has a dimeric structure like the S. cerevisiae enzyme. Site-directed mutagenesis of a crucial residue in loop220s (Asp(233)) revealed differences concerning the intramolecular signal transduction for allosteric regulation of enzymatic activity.