A self-stabilized model of the chymotrypsin catalytic pocket. The energy profile of the overall catalytic cycle

A model of the catalytic triad of chymotrypsin is built assuring the arrangement and properties as they are within the complete enzyme. The model contains 18 amino acid residues of chymotrypsin and its substrate. A total of 135 atoms (including 70 heavy atoms) were subjected to full ab initio geometry optimizations through 127 individual steps along the reaction coordinate of the complete catalytic mechanism. It was shown that the described model of the catalytic apparatus forms a self-stabilized molecule ensemble without the rest of the enzyme and substrate. According to the calculations, the formations of the first and second tetrahedral intermediates in the model have 20.3 and 15.7 kcal/mol activation energy barriers, respectively. Removing elements of the catalytic apparatus such as the (1) catalytic aspartate or (2) the anion hole, as well as (3) inserting a water molecule "early" in the catalytic process, or (4) introducing conformational rigidity of the substrate, results in an increase of the above energy barrier of the first catalytic step in the model by 6.4, 13.7, 3.7, and 4.1 kcal/mol, respectively. Based on the calculated process one can conclude that the catalytic reaction in this model is much more similar to the reaction in the enzyme than to the reference reaction. To our knowledge, this is the first model system that mimics the complete catalytic mechanism.     

Structure and mechanism of human cytosolic phospholipase A(2)

cPLA(2) is an 85-kDa enzyme whose primary function, the release of arachidonic acid from phospholipid membranes, is a crucial reaction in the metabolism of lipid mediators of inflammation. cPLA(2) consists of two domains: an N-terminal, C2-type unit analogous to those present in other membrane-targeting molecules, and a catalytic domain harboring an active site dyad at the bottom of a deep, mostly hydrophobic catalytic funnel. The absence of a third active site residue in the cPLA(2) cleft, as observed in other lipases, suggests that the enzyme proceeds through a novel catalytic mechanism. Crystallographic and biochemical studies of cPLA(2) will provide essential information for the development of small molecule inhibitors which may be employed in the control of inflammatory and other highly regulated processes.     

Systematic comparison of catalytic mechanisms of hydrolysis and transfer reactions classified in the EzCatDB database

Catalytic mechanisms of 270 enzymes from 131 superfamilies, mainly hydrolases and transferases, were analyzed based on their enzyme structures. A method of systematic comparison and classification of the catalytic reactions was developed. Hydrolysis and transfer reactions closely resemble one another, displaying common mechanisms, single displacement, and double displacement. These displacement mechanisms might be further subclassified according to the type of catalytic factors and nucleophilic substitution involved. Several types of catalytic factors exist: nucleophile, acid, base, stabilizer, modulator, cofactors. Nucleophilic substitution might be categorized as S(N)1/S(N)2 (or dissociative/associative) reactions. The classification indicates that some mechanisms favor particular types of catalytic factors. In hydrolyses of amide bonds and phosphoric ester bonds, mechanisms with single displacement tend to use inorganic cofactors such as zinc and magnesium ions as important catalysts, whereas those with double displacement frequently do not use such cofactors. In contrast, hydrolyses of O-glycoside bond rarely use such cofactors, with one exception. The trypsin-like hydrolytic reaction, which is catalyzed by the classic catalytic triad comprising serine/histidine/aspartate, can be considered as a "super-reaction" because it is observed in at least three nonhomologous enzymes, whereas most reactions are singlets without any nonhomologous enzymes. By dividing complex reactions into several reactions, correlations between active site structures and catalytic functions can be suggested. This classification method is applicable to other reactions such as elimination and isomerization. Furthermore, it will facilitate annotation of enzyme functions from 3D patterns of enzyme active sites. The classification is available at http://mbs.cbrc.jp/EzCatDB/RLCP/index.html.     

Reaction mechanism of the membrane-bound ATPase of submitochondrial particles from beef heart

Submitochondrial particles from beef heart, washed with dilute solutions of KCl so as to activate the latent, membrane-bound ATPase, F1, may be used to study single site catalysis by the enzyme. [gamma-32P]ATP, incubated with a molar excess of catalytic sites, a condition which favors binding of substrate in only a single catalytic site on the enzyme, is hydrolyzed via a four-step reaction mechanism. The mechanism includes binding in a high affinity catalytic site, Ka = 10(12)M-1, a hydrolytic step for which the equilibrium constant is near unity, and two product release steps in which Pi dissociates from catalytic sites about 10 times more rapidly than ADP. Catalysis by the membrane-bound ATPase also is characterized by a 10(6)-fold acceleration in the rate of net hydrolysis of [gamma-32P]ATP, bound in the high affinity catalytic site, that occurs when substrate is made available to additional catalytic sites on the enzyme. These aspects of the reaction mechanism of the ATPase of submitochondrial particles closely parallel the reaction mechanism determined for solubilized, homogeneous F1 (Grubmeyer, C., Cross, R. L., and Penefsky, H. S. (1982) J. Biol. Chem. 257, 12092-12100). The finding that removal of the enzyme from the membrane does not significantly alter the properties of single site catalysis lends support to models of ATP synthesis in oxidative phosphorylation, catalyzed by membrane-bound F1, that have been based on the study of the soluble enzyme.     

The chemistry of protein catalysis

We report, for the first time, on the statistics of chemical mechanisms and amino acid residue functions that occur in enzyme reaction sequences using the MACiE database of 202 distinct enzyme reaction mechanisms as a knowledge base. MACiE currently holds representatives from each Enzyme Commission sub-subclass where there is an available crystal structure and sufficient evidence in the primary literature for a mechanism. Each catalytic step of every reaction sequence in MACiE is fully annotated, so that it includes the function of the catalytic residues involved in the reaction and the chemical mechanisms by which substrates are transformed into products. We show that the most catalytic amino acid residues are histidine, cysteine and aspartate, which are also the residues whose side-chains are more likely to serve as reactants, and that have the greatest versatility of function. We show that electrophilic reactions in enzymes are very rare, and the majority of enzyme reactions rely upon nucleophilic and general acid/base chemistry. However, although rare, radical (homolytic) reactions are much more common than electrophilic reactions. Thus, the majority of amino acid residues perform stabilisation roles (as spectators) or proton shuttling roles (as reactants). The analysis presented provides a better understanding of the mechanisms of enzyme catalysis and may act as an initial step in the validation and prediction of mechanism in an enzyme active site.     

Mechanistic Investigations of Unsaturated Glucuronyl Hydrolase from Clostridium perfringens

Experiments were carried out to probe the details of the hydration-initiated hydrolysis catalyzed by the Clostridium perfringens unsaturated glucuronyl hydrolase of glycoside hydrolase family 88 in the CAZy classification system. Direct ¹H NMR monitoring of the enzymatic reaction detected no accumulated reaction intermediates in solution, suggesting that rearrangement of the initial hydration product occurs on-enzyme. An attempt at mechanism-based trapping of on-enzyme intermediates using a 1,1-difluoro-substrate was unsuccessful because the probe was too deactivated to be turned over by the enzyme. Kinetic isotope effects arising from deuterium-for-hydrogen substitution at carbons 1 and 4 provide evidence for separate first-irreversible and overall rate-determining steps in the hydration reaction, with two potential mechanisms proposed to explain these results. Based on the positioning of catalytic residues in the enzyme active site, the lack of efficient turnover of a 2-deoxy-2-fluoro-substrate, and several unsuccessful attempts at confirmation of a simpler mechanism involving a covalent glycosyl-enzyme intermediate, the most plausible mechanism is one involving an intermediate bearing an epoxide on carbons 1 and 2.     

Identification of the catalytic residue of Thermococcus litoralis 4-alpha-glucanotransferase through mechanism-based labeling

Thermococcus litoralis 4-alpha-glucanotransferase (TLGT) belongs to family 57 of glycoside hydrolases and catalyzes the disproportionation and cycloamylose synthesis reactions. Family 57 glycoside hydrolases have not been well investigated, and even the catalytic mechanism involving the active site residues has not been studied. Using 3-ketobutylidene-beta-2-chloro-4-nitrophenyl maltopentaoside (3KBG5CNP) as a donor and glucose as an acceptor, we showed that the disproportionation reaction of TLGT involves a ping-pong bi-bi mechanism. On the basis of this reaction mechanism, the glycosyl-enzyme intermediate, in which a donor substrate was covalently bound to the catalytic nucleophile, was trapped by treating the enzyme with 3KBG5CNP in the absence of an acceptor and was detected by matrix-assisted laser desorption ionization time-of-flight mass spectrometry after peptic digestion. Postsource decay analysis suggested that either Glu-123 or Glu-129 was the catalytic nucleophile of TLGT. Glu-123 was completely conserved between family 57 enzymes, and the catalytic activity of the E123Q mutant enzyme was greatly decreased. On the other hand, Glu-129 was a variable residue, and the catalytic activity of the E129Q mutant enzyme was not decreased. These results indicate that Glu-123 is the catalytic nucleophile of TLGT. Sequence alignment of TLGT and family 38 enzymes (class II alpha-mannosidases) revealed that Glu-123 of TLGT corresponds to the nucleophilic aspartic acid residue of family 38 glycoside hydrolases, suggesting that family 57 and 38 glycoside hydrolases may have had a common ancestor.     

An Acyl-covalent Enzyme Intermediate of Lecithin:Retinol Acyltransferase

Synthesis of fatty acid retinyl esters determines systemic vitamin A levels and provides substrate for production of visual chromophore (11-cis-retinal) in vertebrates. Lecithin:retinol acyltransferase (LRAT), the main enzyme responsible for retinyl ester formation, catalyzes the transfer of an acyl group from the sn-1 position of phosphatidylcholine to retinol. To delineate the catalytic mechanism of this reaction, we expressed and purified a fully active, soluble form of this enzyme and used it to examine the possible formation of a transient acyl-enzyme intermediate. Detailed mass spectrometry analyses revealed that LRAT undergoes spontaneous, covalent modification upon incubation with a variety of phosphatidylcholine substrates. The addition of an acyl chain occurs at the Cys¹⁶¹ residue, indicating formation of a thioester intermediate. This observation provides the first direct experimental evidence of thioester intermediate formation that constitutes the initial step in the proposed LRAT catalytic reaction. Additionally, we examined the effect of increasing fatty acyl side chain length in phosphatidylcholine on substrate accessibility in this reaction, which provided insights into the function of the single membrane-spanning domain of LRAT. These observations are critical to understanding the catalytic mechanism of LRAT protein family members as well as other lecithin:acyltransferases wherein Cys residues are required for catalysis.     

Yeast cystathionine beta-synthase is a pyridoxal phosphate enzyme but, unlike the human enzyme, is not a heme protein

Our studies of cystathionine beta-synthase from Saccharomyces cerevisiae (yeast) are aimed at clarifying the cofactor dependence and catalytic mechanism and obtaining a system for future investigations of the effects of mutations that cause human disease (homocystinuria or coronary heart disease). We report methods that yielded high expression of the yeast gene in Escherichia coli and of purified yeast cystathionine beta-synthase. The absorption and circular dichroism spectra of the homogeneous enzyme were characteristic of a pyridoxal phosphate enzyme and showed the absence of heme, which is found in human and rat cystathionine beta-synthase. The absence of heme in the yeast enzyme facilitates spectroscopic studies to probe the catalytic mechanism. The reaction of the enzyme with L-serine in the absence of L-homocysteine produced the aldimine of aminoacrylate, which absorbed at 460 nm and had a strong negative circular dichroism band at 460 nm. The formation of this intermediate from the product, L-cystathionine, demonstrates the partial reversibility of the reaction. Our results establish the overall catalytic mechanism of yeast cystathionine beta-synthase and provide a useful system for future studies of structure and function. The absence of heme in the functional yeast enzyme suggests that heme does not play an essential catalytic role in the rat and human enzymes. The results are consistent with the absence of heme in the closely related enzymes O-acetylserine sulfhydrylase, threonine deaminase, and tryptophan synthase.     

Site-directed mutagenesis of glutamate 317 of bovine alpha-1,3Galactosyltransferase and its effect on enzyme activity: implications for reaction mechanism

Bovine alpha1,3galactosyltransferase (alpha1,3GalT) transfers galactose from UDP-alpha-galactose to terminal beta-linked galactosyl residues with retention of configuration of the incoming galactose residue. The epitope synthesized has been shown to be critical for xenotransplantation. According to a proposed double-displacement reaction mechanism, glutamate-317 (E317) is thought to be the catalytic nucleophile. The proposed catalytic role of E317 involves an initial nucleophilic attack with inversion of configuration and formation of a covalent sugar-enzyme intermediate between E317 and galactose from the donor substrate, followed by a second nucleophilic attack performed by the acceptor substrate with a second inversion of configuration. To determine whether E317 of alpha1,3GalT is critical for enzyme activity, site-directed mutagenesis was used to substitute alanine, aspartic acid, cysteine and histidine for E317. If the proposed reaction mechanism for the alpha1,3GalT enzyme is correct, E317D and E317H would produce active enzymes since they can act as nucleophiles. The non-conservative mutation E317A and conservative mutation E317C are predicted to produce inactive or very low activity enzymes since the E317A mutant cannot engage in a nucleophilic attack, and the E317C mutant would trap the galactose residue. The results obtained demonstrate that E317D and E317H mutants retained activity, albeit significantly less than the wild-type enzyme. Additionally, both E317A and E317C mutant also retained enzyme activity, suggesting that E317 is not the catalytic nucleophile proposed in the double-displacement mechanism. Therefore, either a different amino acid may act as the catalytic nucleophile or the reaction must proceed by a different mechanism.     

The crystal structure of the formiminotransferase domain of formiminotransferase-cyclodeaminase: implications for substrate channeling in a bifunctional enzyme

BACKGROUND: The bifunctional enzyme formiminotransferase-cyclodeaminase (FTCD) contains two active sites at different positions on the protein structure. The enzyme binds a gamma-linked polyglutamylated form of the tetrahydrofolate substrate and channels the product of the transferase reaction from the transferase active site to the cyclodeaminase active site. Structural studies of this bifunctional enzyme and its monofunctional domains will provide insight into the mechanism of substrate channeling and the two catalytic reactions. RESULTS: The crystal structure of the formiminotransferase (FT) domain of FTCD has been determined in the presence of a product analog, folinic acid. The overall structure shows that the FT domain comprises two subdomains that adopt a novel alpha/beta fold. Inspection of the folinic acid binding site reveals an electrostatic tunnel traversing the width of the molecule. The distribution of charged residues in the tunnel provides insight into the possible mode of substrate binding and channeling. The electron density reveals that the non-natural stereoisomer, (6R)-folinic acid, binds to the protein; this observation suggests a mechanism for product release. In addition, a single molecule of glycerol is bound to the enzyme and indicates a putative binding site for formiminoglutamate. CONCLUSIONS: The structure of the FT domain in the presence of folinic acid reveals a possible novel mechanism for substrate channeling. The position of the folinic acid and a bound glycerol molecule near to the sidechain of His82 suggests that this residue may act as the catalytic base required for the formiminotransferase mechanism.     

Site-directed mutagenesis of glutamate 317 of bovine α-1,3Galactosyltransferase and its effect on enzyme activity: Implications for reaction mechanism

Bovine α1,3galactosyltransferase (α1,3GalT) transfers galactose from UDP-α-galactose to terminal β-linked galactosyl residues with retention of configuration of the incoming galactose residue. The epitope synthesized has been shown to be critical for xenotransplantation. According to a proposed double-displacement reaction mechanism, glutamate-317 (E317) is thought to be the catalytic nucleophile. The proposed catalytic role of E317 involves an initial nucleophilic attack with inversion of configuration and formation of a covalent sugar–enzyme intermediate between E317 and galactose from the donor substrate, followed by a second nucleophilic attack performed by the acceptor substrate with a second inversion of configuration. To determine whether E317 of α1,3GalT is critical for enzyme activity, site-directed mutagenesis was used to substitute alanine, aspartic acid, cysteine and histidine for E317. If the proposed reaction mechanism for the α1,3GalT enzyme is correct, E317D and E317H would produce active enzymes since they can act as nucleophiles. The non-conservative mutation E317A and conservative mutation E317C are predicted to produce inactive or very low activity enzymes since the E317A mutant cannot engage in a nucleophilic attack, and the E317C mutant would trap the galactose residue. The results obtained demonstrate that E317D and E317H mutants retained activity, albeit significantly less than the wild-type enzyme. Additionally, both E317A and E317C mutant also retained enzyme activity, suggesting that E317 is not the catalytic nucleophile proposed in the double-displacement mechanism. Therefore, either a different amino acid may act as the catalytic nucleophile or the reaction must proceed by a different mechanism.     

Inactivation of peroxidase by hydrogen peroxide and its protection by a reductant agent

Hydrogen peroxide, the oxidant substrate of peroxidase, is also an inactivating agent of this enzyme. The reductant substrates protect the enzyme from the inactivating process. A reaction mechanism is proposed, in which two competitive routes exist for Compound I of peroxidase; one catalytic and one inactivating. The analytical solution produced at the end of the reaction supports the proposed mechanism and shows the dependence between the number of turnovers of the enzyme (r) and the ratio of both substrates.     

Evolutionary optimization of the catalytic efficiency of enzymes.

1. The rate equation for a generalized Michaelian type of enzymic reaction mechanism has been analyzed in order to establish how the mechanism should be kinetically designed in order to optimize the catalytic efficiency of the enzyme for a given average magnitude of true and apparent first-order rate constants in the mechanism at given concentrations of enzyme, substrate and product. 2. As long as on-velocity constants for substrate and product binding to the enzyme have not reached the limiting value for a diffusion-controlled association process, the optimal state of enzyme operation will be characterized by forward (true and apparent) first-order rate constants of equal magnitude and reverse rate constants of equal magnitude. The drop in free energy driving the catalysed reaction will occur to an equal extent for each reaction step in the mechanism. All internal equilibrium constants will be of equal magnitude and reflect only the closeness of the catalysed reaction to equilibrium conditions. 3. When magnitudes of on-velocity constants for substrate and product binding have reached their upper limits, the optimal kinetic design of the reaction mechanism becomes more complex and has to be established by numerical methods. Numerical solutions, calculated for triosephosphate isomerase, indicate that this particular enzyme may or may not be considered to exhibit close to maximal efficiency, depending on what value is assigned to the upper limit for a ligand association rate constant. 4. Arguments are presented to show that no useful information on the evolutionary optimization of the catalytic efficiency of enzymes can be obtained by previously taken approaches that are based on the application of linear free-energy relationships for rate and equilibrium constants in the reaction mechanism.     

On the mechanism of the metallo-beta-lactamase from Bacteroides fragilis.

The catalytic mechanism of metallo-beta-lactamase from Bacteroides fragilis, a dinuclear Zn(II)-containing enzyme responsible for multiple antibiotic resistance, has been investigated by using nitrocefin as a substrate. Rapid-scanning and single-wavelength stopped-flow studies revealed the accumulation during turnover of an enzyme-bound intermediate with intense absorbance at 665 nm (epsilon = 30 000 M(-1) cm(-1)). The proposed minimum kinetic mechanism for the B. fragilis metallo-beta-lactamase-catalyzed nitrocefin hydrolysis [Wang, Z., and Benkovic, S. J. (1998) J. Biol. Chem. 273, 22402-22408] was confirmed, and more accurate kinetic parameters were obtained from computer simulations and fitting. The intermediate was shown to be a novel anionic species bound to the enzyme through a Zn-acyl linkage and contains a negatively charged nitrogen leaving group. This is the first time such an intermediate was observed in the catalytic cycle of a Zn(II)-containing hydrolase and is evidence for a unique beta-lactam hydrolysis mechanism in which the amine can leave as an anion; prior protonation is not required. The electrostatic interaction between the negatively charged intermediate and the positively charged dinuclear Zn(II) center of the enzyme is important for stabilization of the intermediate. The catalytic reaction was accelerated in the presence of exogenous nucleophiles or anions, and neither the product nor the enzyme was modified during turnover, indicating that a Zn-bound hydroxide (rather than Asp-103) is the active site nucleophile. On the basis of all the information on hand, a catalytic mechanism of the B. fragilis metallo-beta-lactamase is proposed.     

Some properties of the catalytic subunit of adenosine 3':5'-monophosphate-dependent protein kinase from pigeon breast muscle]

A method for the preparation of a homogenous catalytic subunit of adenosine 3':5'-monophosphate-dependent protein kinase from pigeon breast muscle was developed. The molecular weight of the enzyme as determined by electrophoresis in the presence of sodium dodecyl sulfate was found to be 42000. The pH optimum of the catalytic subunit was around 8.0. The active site of the catalytic subunit was studied using some derivatives of ATP, containing different reactive groups in the triphosphate chain of the molecule. It may be assumed that the pH optimum of the enzyme inactivation by adenosine 5'-chloromethylpyrophosphonate and the protective effect of ATP suggest covalent binding of the imidazole ring in the enzyme active site. The kinetic mechanism of the protein kinase reaction was studied using the initial rate experiments and reaction product inhibition. The results obtained were consistent with a random Bi-Bi kinetic mechanism.     

Comparison of the Intrinsic Dynamics of Aminoacyl-tRNA Synthetases

Aminoacyl-tRNA synthetases (AARSs) are an important family of enzymes that catalyze tRNA aminoacylation reaction (Ibba and Soll in Annu Rev Biochem 2000, 69:617–650) [1]. AARSs are grouped into two broad classes (class I and II) based on sequence/structural homology and mode of their interactions with the tRNA molecule (Ibba and Soll in Annu Rev Biochem 2000, 69:617–650) [1]. As protein dynamics play an important role in enzyme function, we explored the intrinsic dynamics of these enzymes using normal mode analysis and investigated if the two classes and six subclasses (Ia–c and IIa–c) of AARSs exhibit any distinct patterns of motion. The present study found that the intrinsic dynamics-based classification of these enzymes is similar to that obtained based on sequence/structural homology for most enzymes. However, the classification of seryl-tRNA synthetase was not straightforward; the internal mobility patterns of this enzyme are comparable to both IIa and IIb AARSs. This study revealed only a few general mobility patterns in these enzymes—(1) the insertion domain is generally engaged in anticorrelated motion with respect to the catalytic domain for both classes of AARSs and (2) anticodon binding domain dynamics are partly correlated and partly anticorrelated with respect to other domains for class I enzymes. In most of the class II AARSs, the anticodon binding domain is predominately engaged in anticorrelated motion with respect to the catalytic domain and correlated to the insertion domain. This study supports the notion that dynamic-based classification could be useful for functional classification of proteins.     

Identification of the catalytic mechanism and estimation of kinetic parameters for fumarase

The enzyme fumarase catalyzes the reversible hydration of fumarate to malate. The reaction catalyzed by fumarase is critical for cellular energetics as a part of the tricarboxylic acid cycle, which produces reducing equivalents to drive oxidative ATP synthesis. A catalytic mechanism for the fumarase reaction that can account for the kinetic behavior of the enzyme observed in both isotope exchange studies and initial velocity studies has not yet been identified. In the present study, we develop an 11-state kinetic model of the enzyme based on the current consensus on its catalytic mechanism and design a series of experiments to estimate the model parameters and identify the major flux routes through the mechanism. The 11-state mechanism accounts for competitive binding of inhibitors and activation by different anions, including phosphate and fumarate. The model is identified from experimental time courses of the hydration of fumarate to malate obtained over a wide range of buffer and substrate concentrations. Further, the 11-state model is found to effectively reduce to a five-state model by lumping certain successive steps together to yield a mathematically less complex representation that is able to match the data. Analysis suggests the primary reaction route of the catalytic mechanism, with fumarate binding to the free unprotonated enzyme and a proton addition prior to malate release in the fumarate hydration reaction. In the reverse direction (malate dehydration), malate binds the protonated form of the enzyme, and a proton is generated before fumarate is released from the active site.     

Clinical Use of Aromatase Inhibitors: Current Data and Future Perspectives

Abstract

A systematic procedure for the kinetic study of irreversible inhibition when the enzyme is consumed in the reaction which it catalyses, has been developed and analysed. Whereas in most reactions the enzymes are regenerated after each catalytic event and serve as reusable transacting effectors, in the consumed enzymes each catalytic center participates only once and there is no enzyme turnover. A systematic kinetic analysis of irreversible inhibition of these enzyme reactions is presented. Based on the algebraic criteria proposed in this work, it should be possible to evaluate either the mechanism of inhibition (complexing or non-complexing), or the type of inhibition (competitive, non-competitive, uncompetitive, mixed non-competitive). In addition, all kinetic constants involved in each case could be calculated. An experimental application of this analysis is also presented, concerning peptide bond formation in vitro. Using the puromycin reaction, which is a model reaction for the study of peptide bond formation in vitro and which follows the same kinetic law as the enzymes under study, we have found that: (i) the antibiotic spiramycin inhibits the puromycin reaction as a competitive irreversible inhibitor in a one step mechanism with an association rate constant equal to 1.3 × 10⁴M^-1s^-1 and, (ii) hydroxylamine inhibits the same reaction as an irreversible non-competitive inhibitor also in a one step mechanism with a rate constant equal to 1.6 × 10^-3 M^-1s^-1.     

Potential of mean force calculation for the proton and hydride transfer reactions catalyzed by medium-chain acyl-CoA dehydrogenase: effect of mutations on enzyme catalysis

Potential of mean force calculations have been performed on the wild-type medium-chain acyl-CoA dehydrogenase (MCAD) and two of its mutant forms. Initial simulation and analysis of the active site of the enzyme reveal that an arginine residue (Arg256), conserved in the substrate-binding domain of this group of enzymes, exists in two alternate conformations, only one of which makes the enzyme active. This active conformation was used in subsequent computations of the enzymatic reactions. It is known that the catalytic alpha,beta-dehydrogenation of fatty acyl-CoAs consists of two C-H bond dissociation processes: a proton abstraction and a hydride transfer. Energy profiles of the two reaction steps in the wild-type MCAD demonstrate that the reaction proceeds by a stepwise mechanism with a transient species. The activation barriers of the two steps differ by only approximately 2 kcal/mol, indicating that both may contribute to the rate-limiting process. Thus this may be a stepwise dissociation mechanism whose relative barriers can be tuned by suitable alterations of the substrate and/or enzyme. Analysis of the structures along the reaction path reveals that Arg256 plays a key role in maintaining the reaction center hydrogen-bonding network involving the thioester carbonyl group, which stabilizes transition states as well as the intervening transient species. Mutation of this arginine residue to glutamine increases the activation barrier of the hydride transfer reaction by approximately 5 kcal/mol, and the present simulations predict a substantial loss of catalytic activity for this mutant. Structural analysis of this mutant reveals that the orientation of the thioester moiety of the substrate has been changed significantly as compared to that in the wild-type enzyme. In contrast, simulation of the active site of the Thr168Ala mutant shows no significant change in the relative orientation of the substrate and the cofactor in the active site; as a result, this mutation has very little effect on the overall reaction barrier, and this is consistent with the experimental data. This study demonstrates that significant insights into the catalytic mechanism can be obtained from simulation studies, and the results can be used to design novel mechanistic probes for the enzyme.