Transition state variation in enzymatic reactions

Experimental analysis of enzymatic transition states by kinetic isotope effect methods has established geometric variation in related transition state structures. Differences are apparent in development of the reaction coordinate, in solvolytic transition states relative to those in enzymatic catalytic sites, in the stereochemistry of related substrates at the transition state, and in reactions catalyzed by related enzymes.     

Enzymatic transition states: thermodynamics, dynamics and analogue design

Kinetic isotope effects and computational chemistry have defined the transition state structures for several members of the N-ribosyltransferase family. Transition state analogues designed to mimic their cognate transition state structures are among the most powerful enzyme inhibitors. In complexes of N-ribosyltransferases with their transition state analogues, the dynamic nature of the transition state is converted to an ordered, thermodynamic structure closely related to the transition state. This phenomenon is documented by peptide bond H/D exchange, crystallography and computational chemistry. Complexes with substrate, transition state and product analogues reveal reaction coordinate motion and catalytic interactions. Isotope-edited spectroscopic analysis and binding specificity of these complexes provides information about specific enzyme-transition state contacts. In combination with protein dynamic QM/MM models, it is proposed that the transition state is reached by stochastic dynamic excursions of the protein groups near the substrates in the closed conformation. Examples from fully dissociated (D(N) *A(N)), hybrid (D(N)A(N)) and symmetric nucleophilic displacement (A(N)D(N)) transition states are found in the N-ribosyltransferases. The success of transition state analogue inhibitor design based on kinetic isotope effects validates this approach to understanding enzymatic transition states.     

An aspartate transcarbamylase lacking catalytic subunit interactions. Study of conformational changes by ultraviolet absorbance and circular dichroism spectroscopy.

A modified form of aspartate transcarbamylase is synthesized by Escherichia coli in the presence of 2-thiouracil which does not exhibit homotropic cooperative interactions between active sites yet retains heterotropic cooperative interactions due to nucleotide binding. The conformational changes induced in the modified enzyme by the binding of different ligands (substrates, substrate analogs, a transition state analog, and nucleotide effectors) were studied using ultraviolet absorbance and circular dichroism difference spectroscopy. Comparison of the results for the modified enzyme and its isolated subunits to those for the native enzyme and its isolated subunits showed that the conformational changes detected by these methods are qualitatively similar in the two enzymes. Comparison of the absorbance difference spectra due to the binding of a transition substrate analog to the intact native or modified enzymes to the corresponding results for the isolated subunits suggested that ligand binding causes an increased exposure to solvent of certain tyrosyl and phenylalanyl residues in the intact enzymes but not in the isolated subunits. This result is consistent with a diminution of subunit contacts due to substrate binding in the course of homotropic interactions in the native enzyme. Such conformational changes, though perhaps necessary for homotropic cooperativity, are not sufficient to cause homotropic cooperativity since the modified enzyme gave identical perturbations. Interactions of the transition state analog, N-(phosphonacetyl)-L-aspartate, with the modified enzyme were studied. Enzyme kinetic data obtained at low aspartate concentrations showed that this transition state analog does not stimulate activity, but rather exhibits the inhibition predicted for the total absence of homotropic cooperative interactions in the modified enzyme. Spectrophotometric titrations of the number of catalytic sites with the transition state analog showed that the modified enzyme and its isolated subunits possess, respectively, four and two high affinity sites for the inhibitor instead of six and three observed in the case of the normal enzyme and its isolated catalytic subunits. These results are correlated with the lower specific enzymatic activities of the modified enzyme and its catalytic subunits compared to the normal corresponding enzymatic species.     

Selection and characterization of lipase abzyme from phage displayed antibody libraries

Antibodies with enzymatic activity were named abzymes or catalytic antibodies. In the present study, the lipolytic abzymes were selected from the phage displayed antibody libraries against a transition state analog (TSA) of lipases/esterases. After three rounds of selection, four monoclonal phage particles capable of binding significantly with the TSA were obtained. The soluble scFv antibody fragments were further expressed and obtained using Escherichia coli strain HB2151. The binding capabilities and the apparent enzymatic activities of the purified antibody proteins were measured. The 3D structures of the expressed antibodies were also predicted through homology modeling and binding-site prediction algorithm. The present method demonstrates that selection from phage displayed antibody libraries is an efficient and convenient means to find new abzymes.     

Binding of the transition state analog MgADP-fluoroaluminate to F1-ATPase

Escherichia coli F1-ATPase from mutant betaY331W was potently inhibited by fluoroaluminate plus MgADP but not by MgADP alone. beta-Trp-331 fluorescence was used to measure MgADP binding to catalytic sites. Fluoroaluminate induced a very large increase in MgADP binding affinity at catalytic site one, a smaller increase at site two, and no effect at site three. Mutation of either of the critical catalytic site residues beta-Lys-155 or beta-Glu-181 to Gln abolished the effects of fluoroaluminate on MgADP binding. The results indicate that the MgADP-fluoroaluminate complex is a transition state analog and independently demonstrate that residues beta-Lys-155 and (particularly) beta-Glu-181 are important for generation and stabilization of the catalytic transition state. Dicyclohexylcarbodiimide-inhibited enzyme, with 1% residual steady-state ATPase, showed normal transition state formation as judged by fluoroaluminate-induced MgADP binding affinity changes, consistent with a proposed mechanism by which dicyclohexylcarbodiimide prevents a conformational interaction between catalytic sites but does not affect the catalytic step per se. The fluorescence technique should prove valuable for future transition state studies of F1-ATPase.     

Correlated conformational fluctuations during enzymatic catalysis: Implications for catalytic rate enhancement

Correlated enzymatic conformational fluctuations are shown to contribute to the rate of enhancement achieved during catalysis. Cytidine deaminase serves as a model system. Crystallographic temperature factor data for this enzyme complexed with substrate analog, transition-state analog, and product are available, thereby establishing a measure of atomic scale conformational fluctuations along the (approximate) reaction coordinate. First, a neural network-based algorithm is used to visualize the decreased conformational fluctuations at the transition state. Second, a dynamic diffusion equation along the reaction coordinate is solved and shows that the flux velocity through the associated enzymatic conformation space is greatest at the transition state. These results suggest (1) that there are both dynamic and energetic restrictions to conformational fluctuations at the transition state, (2) that enzymatic catalysis occurs on a fluctuating potential energy surface, and (3) a form for the potential energy. The Michaelis-Menten equations are modified to describe catalysis on this fluctuating potential energy profile, leading to enhanced catalytic rates when fluctuations along the reaction coordinate are appropriately correlated. This represents a dynamic tuning of the enzyme for maximally effective transformation of the ES complex into EP.     

How similar are enzyme active site geometries derived from quantum mechanical theozymes to crystal structures of enzyme‐inhibitor complexes? Implications for enzyme design

Quantum mechanical optimizations of theoretical enzymes (theozymes), which are predicted catalytic arrays of biological functionalities stabilizing a transition state, have been carried out for a set of nine diverse enzyme active sites. For each enzyme, the theozyme for the rate-determining transition state plus the catalytic groups modeled by side-chain mimics was optimized using B3LYP/6-31G(d) or, in one case, HF/3-21G(d) quantum mechanical calculations. To determine if the theozyme can reproduce the natural evolutionary catalytic geometry, the positions of optimized catalytic atoms, i.e., covalent, partial covalent, or stabilizing interactions with transition state atoms, are compared to the positions of the atoms in the X-ray crystal structure with a bound inhibitor. These structure comparisons are contrasted to computed substrate-active site structures surrounded by the same theozyme residues. The theozyme/transition structure is shown to predict geometries of active sites with an average RMSD of 0.64 A from the crystal structure, while the RMSD for the bound intermediate complexes are significantly higher at 1.42 A. The implications for computational enzyme design are discussed.     

X-ray, neutron and NMR studies of the catalytic mechanism of aspartic proteinases

Current proposals for the catalytic mechanism of aspartic proteinases are largely based on X-ray structures of bound oligopeptide inhibitors possessing non-hydrolysable analogues of the scissile peptide bond. Until recent years, the positions of protons on the catalytic aspartates and the ligand in these complexes had not been determined with certainty due to the inadequate resolution of these analyses. There has been much interest in locating the catalytic protons at the active site of aspartic proteinases since this has major implications for detailed understanding of the mechanism of action and the design of improved transition state mimics for therapeutic applications. In this review we discuss the results of studies which have shed light on the locations of protons at the catalytic centre. The first direct determination of the proton positions stemmed from neutron diffraction data collected from crystals of the fungal aspartic proteinase endothiapepsin bound to a transition state analogue (H261). The neutron structure of the complex at a resolution of 2.1 A provided evidence that Asp 215 is protonated and that Asp 32 is the negatively charged residue in the transition state complex. Atomic resolution X-ray studies of inhibitor complexes have corroborated this finding. A similar study of the native enzyme established that it, unexpectedly, has a dipeptide bound at the catalytic site which is consistent with classical reports of inhibition by short peptides and the ability of pepsins to catalyse transpeptidation reactions. Studies by NMR have confirmed the findings of low-barrier and single-well hydrogen bonds in the complexes with transition state analogues.     

Atomic resolution analysis of the catalytic site of an aspartic proteinase and an unexpected mode of binding by short peptides

The X-ray structures of native endothiapepsin and a complex with a hydroxyethylene transition state analog inhibitor (H261) have been determined at atomic resolution. Unrestrained refinement of the carboxyl groups of the enzyme by using the atomic resolution data indicates that both catalytic aspartates in the native enzyme share a single negative charge equally; that is, in the crystal, one half of the active sites have Asp 32 ionized and the other half have Asp 215 ionized. The electron density map of the native enzyme refined at 0.9 A resolution demonstrates that there is a short peptide (probably Ser-Thr) bound noncovalently in the active site cleft. The N-terminal nitrogen of the dipeptide interacts with the aspartate diad of the enzyme by hydrogen bonds involving the carboxyl of Asp 215 and the catalytic water molecule. This is consistent with classical findings that the aspartic proteinases can be inhibited weakly by short peptides and that these enzymes can catalyze transpeptidation reactions. The dipeptide may originate from autolysis of the N-terminal Ser-Thr sequence of the enzyme during crystallization.     

The 2.0 A structure of human hypoxanthine-guanine phosphoribosyltransferase in complex with a transition-state analog inhibitor.

The structure of human HGPRT bound to the transition-state analog immucillinGP and Mg2+-pyrophosphate has been determined to 2.0 A resolution. ImmucillinGP was designed as a stable analog with the stereoelectronic features of the transition state. Bound inhibitor at the catalytic site indicates that the oxocarbenium ion of the transition state is stabilized by neighboring-group participation from MgPPi and O5'. A short hydrogen bond forms between Asp 137 and the purine ring analog. Two Mg2+ ions sandwich the pyrophosphate and contact both hydroxyls of the ribosyl analog. The transition-state analog is shielded from bulk solvent by a catalytic loop that moves approximately 25 A to cover the active site and becomes an ordered antiparallel beta-sheet.     

Protein kinase structure and function analysis with chemical tools

Protein kinases are the largest enzyme superfamily involved in cell signal transduction and represent therapeutic targets for a range of diseases. There have been intensive efforts from many labs to understand their catalytic mechanisms, discover inhibitors and discern their cellular functions. In this review, we will describe two approaches developed to analyze protein kinases: bisubstrate analog inhibition and phosphonate analog utilization. Both of these methods have been used in combination with the protein semisynthesis method expressed protein ligation to advance our understanding of kinase-substrate interactions and functional elucidation of phosphorylation. Previous work on the nature of the protein kinase mechanism suggests it follows a dissociative transition state. A bisubstrate analog was designed against the insulin receptor kinase to mimic the geometry of a dissociative transition state reaction coordinate distance. This bisubstrate compound proved to be a potent inhibitor against the insulin receptor kinase and occupied both peptide and nucleotide binding sites. Bisubstrate compounds with altered hydrogen bonding potential as well as varying spacers between the adenine and the peptide demonstrate the importance of the original design features. We have also shown that related bisubstrate analogs can be used to potently block serine/threonine kinases including protein kinase A. Since many protein kinases recognize folded protein substrates for efficient phosphorylation, it was advantageous to incorporate the peptide-ATP conjugates into protein structures. Using expressed protein ligation, a Src-ATP conjugate was produced and shown to be a high affinity ligand for the Csk tyrosine kinase. Nonhydrolyzable mimics of phosphoSer/phosphoTyr can be useful in examining the functionality of phosphorylation events. Using expressed protein ligation, we have employed phosphonomethylene phenylalanine and phosphonomethylene alanine to probe the phosphorylation of Tyr and Ser, respectively. These tools have permitted an analysis of the SH2-phosphatases (SHP1 and SHP2), revealing a novel intramolecular stimulation of catalytic activity mediated by the corresponding phosphorylation events. They have also been used to characterize the cellular regulation of the melatonin rhythm enzyme by phosphorylation.     

New probes of the F1-ATPase catalytic transition state reveal that two of the three catalytic sites can assume a transition state conformation simultaneously

MgADP in combination with fluoroscandium (ScFx) is shown to form a potently inhibitory, tightly bound, noncovalent complex at the catalytic sites of F(1)-ATPase. The F(1).MgADP.ScFx complex mimics a catalytic transition state. Notably, ScFx caused large enhancement of MgADP binding affinity at both catalytic sites 1 and 2, with little effect at site 3. These results indicate that sites 1 and 2 may form a transition state conformation. A new direct optical probe of F(1)-ATPase catalytic transition state conformation is also reported, namely, substantial enhancement of fluorescence emission of residue beta-Trp-148 observed upon binding of MgADP.ScFx or MgIDP. ScFx. Using this fluorescence signal, titrations were performed with MgIDP.ScFx which demonstrated that catalytic sites 1 and 2 can both form a transition state conformation but site 3 cannot. Supporting data were obtained using MgIDP-fluoroaluminate. Current models of the MgATP hydrolysis mechanism uniformly make the assumption that only one catalytic site hydrolyzes MgATP at any one time. The fluorometal analogues demonstrate that two sites have the capability to form the transition state simultaneously.     

High resolution crystal structures of human Rab5a and five mutants with substitutions in the catalytically important phosphate-binding loop

GTPase domain crystal structures of Rab5a wild type and five variants with mutations in the phosphate-binding loop are reported here at resolutions up to 1.5 A. Of particular interest, the A30P mutant was crystallized in complexes with GDP, GDP+AlF(3), and authentic GTP, respectively. The other variant crystals were obtained in complexes with a non-hydrolyzable GTP analog, GppNHp. All structures were solved in the same crystal form, providing an unusual opportunity to compare structures of small GTPases with different catalytic rates. The A30P mutant exhibits dramatically reduced GTPase activity and forms a GTP-bound complex stable enough for crystallographic analysis. Importantly, the A30P structure with bound GDP plus AlF(3) has been solved in the absence of a GTPase-activating protein, and it may resemble that of a transition state intermediate. Conformational changes are observed between the GTP-bound form and the transition state intermediate, mainly in the switch II region containing the catalytic Gln(79) residue and independent of A30P mutation-induced local alterations in the P-loop. The structures suggest an important catalytic role for a P-loop backbone amide group, which is eliminated in the A30P mutant, and support the notion that the transition state of GTPase-mediated GTP hydrolysis is of considerable dissociative character.     

The magnitude of enzyme transition state analog binding constants

Transition state binding theory utilizes non-enzymic and enzymic rate ratios to predict the ratio of transition state analog dissociation constants to substrate dissociation constants. In this paper we show that enzyme rate accelerations due solely to lessened entropy requirements, arising from the juxtaposition of a catalytic group and a substrate binding site at an enzyme active site, will result in a ratio of transition state and substrate dissociation constants which is different, in general, from the ratio of non-enzymic and enzymic rate constants. The arguments presented in this paper provide a possible explanation for the frequently observed large discrepancy between the measured and predicted values for transition state analog dissociation constants.     

The dynamic role of distal side residues in heme hydroperoxidase catalysis. Interplay between X-ray crystallography and ab initio MD simulations

The enzymatic cycle of hydroperoxidases involves the resting Fe(III) state of the enzyme and the high-valent iron intermediates Compound I and Compound II. These states might be characterized by X-ray crystallography and the transition pathways between each state can be investigated using atomistic simulations. Here we review our recent work in the modeling of two key steps of the enzymatic reaction of hydroperoxidases: the formation of Cpd I in peroxidase and the reduction of Cpd I in catalase. It will be shown that small conformational motions of distal side residues (His in peroxidases and His/Asn in catalases), not,or only partially, revealed by the available X-ray structures, play an important role in the catalytic processes examined.     

A computational method for design of connected catalytic networks in proteins

Computational design of new active sites has generally proceeded by geometrically defining interactions between the reaction transition state(s) and surrounding side‐chain functional groups which maximize transition‐state stabilization, and then searching for sites in protein scaffolds where the specified side‐chain–transition‐state interactions can be realized. A limitation of this approach is that the interactions between the side chains themselves are not constrained. An extensive connected hydrogen bond network involving the catalytic residues was observed in a designed retroaldolase following directed evolution. Such connected networks could increase catalytic activity by preorganizing active site residues in catalytically competent orientations, and enabling concerted interactions between side chains during catalysis, for example, proton shuffling. We developed a method for designing active sites in which the catalytic side chains, in addition to making interactions with the transition state, are also involved in extensive hydrogen bond networks. Because of the added constraint of hydrogen‐bond connectivity between the catalytic side chains, to find solutions, a wider range of interactions between these side chains and the transition state must be considered. Our new method starts from a ChemDraw‐like two‐dimensional representation of the transition state with hydrogen‐bond donors, acceptors, and covalent interaction sites indicated, and all placements of side‐chain functional groups that make the indicated interactions with the transition state, and are fully connected in a single hydrogen‐bond network are systematically enumerated. The RosettaMatch method can then be used to identify realizations of these fully‐connected active sites in protein scaffolds. The method generates many fully‐connected active site solutions for a set of model reactions that are promising starting points for the design of fully‐preorganized enzyme catalysts.     

Role of protein conformational dynamics in the catalysis by 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

Enzymatic catalysis has conflicting structural requirements of the enzyme. In order for the enzyme to form a Michaelis complex, the enzyme must be in an open conformation so that the substrate can get into its active center. On the other hand, in order to maximize the stabilization of the transition state of the enzymatic reaction, the enzyme must be in a closed conformation to maximize its interactions with the transition state. The conflicting structural requirements can be resolved by a flexible active center that can sample both open and closed conformational states. For a bisubstrate enzyme, the Michaelis complex consists of two substrates in addition to the enzyme. The enzyme must remain flexible upon the binding of the first substrate so that the second substrate can get into the active center. The active center is fully assembled and stabilized only when both substrates bind to the enzyme. However, the side-chain positions of the catalytic residues in the Michaelis complex are still not optimally aligned for the stabilization of the transition state, which lasts only approximately 10(-13) s. The instantaneous and optimal alignment of catalytic groups for the transition state stabilization requires a dynamic enzyme, not an enzyme which undergoes a large scale of movements but an enzyme which permits at least a small scale of adjustment of catalytic group positions. This review will summarize the structure, catalytic mechanism, and dynamic properties of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase and examine the role of protein conformational dynamics in the catalysis of a bisubstrate enzymatic reaction.     

Catalytic mechanism of the inverting N-acetylglucosaminyltransferase I: DFT quantum mechanical model of the reaction pathway and determination of the transition state structure

The complex N-glycan structures on glycoproteins play important roles in cell adhesion and recognition events in metazoan organisms. A critical step in the biosynthetic pathway leading from high mannose to these complex structures includes the transfer of N-acetylglucosamine (GlcNAc) to a mannose residue by the inverting N-acetylglucosaminyltransferase I (GnT-I). The catalytic mechanism of this enzymatic reaction is explored herein using DFT quantum chemical methods. The computational model used to follow the reaction is based on the X-ray crystallographic structure of GnT-I and contains 127 atoms that represent fragments of residues critical for the substrate binding and catalysis. The mechanism of the catalytic reaction was monitored by means of a 2D potential energy map calculated as a function of predefined reaction coordinates at the B3LYP/6-31G** level. This potential energy surface revealed one transition state associated with a reaction pathway following a concerted mechanism. The reaction barrier was estimated, and the structure of the transition state was characterized at the B3LYP/6-311++G**// B3LYP/6-31G** level.