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.     

1H R1ρ relaxation dispersion experiments in aromatic side chains

Journal of Biomolecular NMR - Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on...     

Conformational Variation in Enzyme Catalysis: A Structural Study on Catalytic Residues

Conformational variation in catalytic residues can be captured as alternative snapshots in enzyme crystal structures. Addressing the question of whether active site flexibility is an intrinsic and essential property of enzymes for catalysis, we present a comprehensive study on the 3D variation of active sites of 925 enzyme families, using explicit catalytic residue annotations from the Mechanism and Catalytic Site Atlas and structural data from the Protein Data Bank. Through weighted pairwise superposition of the functional atoms of active sites, we captured structural variability at single-residue level and examined the geometrical changes as ligands bind or as mutations occur. We demonstrate that catalytic centres of enzymes can be inherently rigid or flexible to various degrees according to the function they perform, and structural variability most often involves a subset of the catalytic residues, usually those not directly involved in the formation or cleavage of bonds. Moreover, data suggest that 2/3 of active sites are flexible, and in half of those, flexibility is only observed in the side chain. The goal of this work is to characterise our current knowledge of the extent of flexibility at the heart of catalysis and ultimately place our findings in the context of the evolution of catalysis as enzymes evolve new functions and bind different substrates.     

Understanding nature's catalytic toolkit

Enzymes catalyse numerous reactions in nature, often causing spectacular accelerations in the catalysis rate. One aspect of understanding how enzymes achieve these feats is to explore how they use the limited set of residue side chains that form their 'catalytic toolkit'. Combinations of different residues form 'catalytic units' that are found repeatedly in different unrelated enzymes. Most catalytic units facilitate rapid catalysis in the enzyme active site either by providing charged groups to polarize substrates and to stabilize transition states, or by modifying the pKa values of other residues to provide more effective acids and bases. Given recent efforts to design novel enzymes, the rise of structural genomics and subsequent efforts to predict the function of enzymes from their structure, these units provide a simple framework to describe how nature uses the tools at her disposal, and might help to improve techniques for designing and predicting enzyme function.     

Conformational Dynamics of Protein Side Chains and Enzyme-Substrate Interaction

Abstract

The electric interaction of charged and polar protein side chains with dipole moments of the bonds of substrate molecules bound in enzyme active sites is considered. The conformational motion of the side chains leads to the fact that the electric interaction, besides a constant (electrostatic) part, contains a fluctuating one, which is a random force (noise) exerted on the substrate molecule. On the time scale of enzyme turnovers this noise can be considered as the white one with good approximation. The noise is external, and the explicit expression for its intensity has been obtained. The possible functional role of the noise as an activating factor in enzyme catalysis is discussed.     

Study of trypsin-substrate and trypsin-inhibitor complexes. 1. Conformation of Asp-102, His-57 and Ser-195 residues in the trypsin active center

Using the semi-empirical method of conformational analysis low energetic conformations were found for trypsin's catalytic triad Asp102, His57, Ser195 in the field of the active centre that consisted of about 800 atoms. It was shown that within the enzyme structure conformations which are preferential for isolated residues are formed. Lability obtained for the side chains is in accord with their functions in enzymatic catalysis. Results of the calculations carried out correlate with the X-ray data for trypsin complexes.     

Mutational and structural studies of the diisopropylfluorophosphatase from Loligo vulgaris shed new light on the catalytic mechanism of the enzyme

The active site, the substrate binding site, and the metal binding sites of the diisopropylfluorophosphatase (DFPase) from Loligo vulgaris have been modified by means of site-directed mutagenesis to improve our understanding of the reaction mechanism. Enzymatic characterization of mutants located in the major groove of the substrate binding pocket indicates that large hydrophobic side chains at these positions are favorable for substrate turnover. Moreover, the active site residue His287 proved to be beneficial, but not essential, for DFP hydrolysis. In most cases, hydrophobic side chains at position 287 led to significant catalytic activities although reduced relative to the wild-type enzyme. With respect to the Ca-1 binding site, where catalysis occurs, various mutants indicated that the net charge at this calcium-binding site as well as the relative positions of the charged calcium ligands is crucial for catalytic activity. The importance of the electrostatic potential at the active site was furthermore revealed by various mutations of residues lining the interior of the central water-filled tunnel, which traverses the entire protein structure. In this respect, the structural features of residue His181, which is located at the opposite end of the DFPase tunnel relative to the active site, were characterized extensively. It was concluded that a tunnel-spanning hydrogen bond network, which includes a large number of apparently slow exchanging water molecules, relays any modifications in the electrostatics of the system to the active site, thus affecting the catalytic reactivity of the enzyme.     

Altered orientation of active site residues in variants of human ferrochelatase. Evidence for a hydrogen bond network involved in catalysis

Ferrochelatase catalyzes the terminal step in heme biosynthesis, the insertion of ferrous iron into protoporphyrin to form protoheme IX. The crystal structures of human ferrochelatase both with and without the protoporphyrin substrate bound have been determined previously. The substrate-free enzyme has an open active site pocket, while in the substrate-bound enzyme, the active site pocket is closed around the porphyrin macrocycle and a number of active site residues have reoriented side chains. To understand how and why these structural changes occur, we have substituted three amino acid residues (H263, H341, and F337) whose side chains occupy different spatial positions in the substrate-free versus substrate-bound ferrochelatases. The catalytic and structural properties of ferrochelatases containing the amino acid substitutions H263C, H341C, and F337A were examined. It was found that in the H263C and H341C variants, but not the F337A variant enzymes, the side chains of N75, M76, R164, H263, F337, H341, and E343 are oriented in a fashion similar to what is found in ferrochelatase with the bound porphyrin substrate. However, all of the variant forms possess open active site pockets which are found in the structure of porphyrin-free ferrochelatase. Thus, while the interior walls of the active site pocket are remodeled in these variants, the exterior lips remain unaltered in position. One possible explanation for this collective reorganization of active site side chains is the presence of a hydrogen bond network among H263, H341, and E343. This network is disrupted in the variants by alteration of H263C or H341C. In the substrate-bound enzyme, the formation of a hydrogen bond between H263 and a pyrrole nitrogen results in disruption of the network. The possible role of this network in catalysis is discussed.     

Crystal structure of human uroporphyrinogen III synthase.

Uroporphyrinogen III synthase, U3S, the fourth enzyme in the porphyrin biosynthetic pathway, catalyzes cyclization of the linear tetrapyrrole, hydroxymethylbilane, to the macrocyclic uroporphyrino gen III, which is used in several different pathways to form heme, siroheme, chlorophyll, F(430) and vitamin B(12). U3S activity is essential in all organisms, and decreased activity in humans leads to the autosomal recessive disorder congenital erythropoetic porphyria. We have determined the crystal structure of recombinant human U3S at 1.85 A resolution. The protein folds into two alpha/beta domains connected by a beta-ladder. The active site appears to be located between the domains, and variations in relative domain positions observed between crystallographically independent molecules indicates the presence of flexibility that may be important in the catalytic cycle. Possible mechanisms of catalysis were probed by mutating each of the four invariant residues in the protein that have titratable side chains. Additionally, six other highly conserved and titratable side chains were also mutated. In no case, however, did one of these mutations abolish enzyme activity, suggesting that the mechanism does not require acid/base catalysis.     

Crystal structure of the reduced Schiff-base intermediate complex of transaldolase B from Escherichia coli: mechanistic implications for class I aldolases.

Transaldolase catalyzes transfer of a dihydroxyacetone moiety from a ketose donor to an aldose acceptor. During catalysis, a Schiff-base intermediate between dihydroxyacetone and the epsilon-amino group of a lysine residue at the active site of the enzyme is formed. This Schiff-base intermediate has been trapped by reduction with potassium borohydride, and the crystal structure of this complex has been determined at 2.2 A resolution. The overall structures of the complex and the native enzyme are very similar; formation of the intermediate induces no large conformational changes. The dihydroxyacetone moiety is covalently linked to the side chain of Lys 132 at the active site of the enzyme. The Cl hydroxyl group of the dihydroxyacetone moiety forms hydrogen bonds to the side chains of residues Asn 154 and Ser 176. The C3 hydroxyl group interacts with the side chain of Asp 17 and Asn 35. Based on the crystal structure of this complex a reaction mechanism for transaldolase is proposed.     

Altering of the metal specificity of Escherichia coli alkaline phosphatase

Analysis of sequence alignments of alkaline phosphatases revealed a correlation between metal specificity and certain amino acid side chains in the active site that are metal-binding ligands. The Zn(2+)-requiring Escherichia coli alkaline phosphatase has an Asp at position 153 and a Lys at position 328. Co(2+)-requiring alkaline phosphatases from Thermotoga maritima and Bacillus subtilis have a His and a Trp at these positions, respectively. The mutations D153H, K328W, and D153H/K328W were induced in E. coli alkaline phosphatase to determine whether these residues dictate the metal dependence of the enzyme. The wild-type and D153H enzymes showed very little activity in the presence of Co(2+), but the K328W and especially the D153H/K328W enzymes effectively use Co(2+) for catalysis. Isothermal titration calorimetry experiments showed that in all cases except for the D153H/K328W enzyme, a possible conformation change occurs upon binding Co(2+). These data together indicate that the active site of the D153H/K328W enzyme has been altered significantly enough to allow the enzyme to utilize Co(2+) for catalysis. These studies suggest that the active site residues His and Trp at the E. coli enzyme positions 153 and 328, respectively, at least partially dictate the metal specificity of alkaline phosphatase.     

Crystal structures of human mitochondrial branched chain aminotransferase reaction intermediates: ketimine and pyridoxamine phosphate forms

The three-dimensional structures of the isoleucine ketimine and the pyridoxamine phosphate forms of human mitochondrial branched chain aminotransferase (hBCATm) have been determined crystallographically at 1.9 A resolution. The hBCATm-catalyzed transamination can be described in molecular terms together with the earlier solved pyridoxal phosphate forms of the enzyme. The active site lysine, Lys202, undergoes large conformational changes, and the pyridine ring of the cofactor tilts by about 18 degrees during catalysis. A major determinant of the enzyme's substrate and stereospecificity for L-branched chain amino acids is a group of hydrophobic residues that form three hydrophobic surfaces and lock the side chain in place. Short-chain aliphatic amino acid side chains are unable to interact through van der Waals contacts with any of the surfaces whereas bulky aromatic side chains would result in significant steric hindrance. As shown by modeling, and in agreement with previous biochemical data, glutamate but not aspartate can form hydrogen bond interactions. The carboxylate group of the bound isoleucine is on the same side as the phosphate group of the cofactor. These active site interactions are largely retained in a model of the human cytosolic branched chain aminotransferase (hBCATc), suggesting that residues in the second tier of interactions are likely to determine the specificity of hBCATc for the drug gabapentin. Finally, the structures reveal a unique role for cysteine residues in the mammalian BCAT. Cys315 and Cys318, which immediately follow a beta-turn (residues 311-314) and are located just outside the active site, form an unusual thiol-thiolate hydrogen bond. This beta-turn positions Thr313 for its interaction with the pyridoxal phosphate oxygens and substrate alpha-carboxylate group.     

Analysis of catalytic residues in enzyme active sites

We present an analysis of the residues directly involved in catalysis in 178 enzyme active sites. Specific criteria were derived to define a catalytic residue, and used to create a catalytic residue dataset, which was then analysed in terms of properties including secondary structure, solvent accessibility, flexibility, conservation, quaternary structure and function. The results indicate the dominance of a small set of amino acid residues in catalysis and give a picture of a general active site environment. It is hoped that this information will provide a better understanding of the molecular mechanisms involved in catalysis and a heuristic basis for predicting catalytic residues in enzymes of unknown function.     

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.     

Structural requirements for catalysis and dimerization of human methionine adenosyltransferase I/III

We have used site-directed mutagenesis to probe the structural requirements for catalysis and dimerization of human hepatic methionine adenosyltransferase (hMAT). We built a homology model of the dimeric hMAT III inferred by the crystal structure of the highly homologous Escherichia coli MAT dimer. The active sites of both enzymes comprise the same amino acids and are located in the inter-subunit interface. All of the amino acids predicted to be in the hMAT III active site were mutated, as well as residues in a conserved ATP binding region. All of the mutations except one severely affected catalytic activity. On the other hand, dimerization was affected only by single mutations of three different residues, all on one monomer. The homology model suggested that the side chains of these residues stabilized the monomer and participated in a bridge between subunits consisting of a network of metal and phosphate ions. In agreement with this observation, we demonstrated that dimerization cannot occur in the absence of phosphate.     

Arginine coordination in enzymatic phosphoryl transfer: evaluation of the effect of Arg166 mutations in Escherichia coli alkaline phosphatase

Arginine residues are commonly found in the active sites of enzymes catalyzing phosphoryl transfer reactions. Numerous site-directed mutagenesis experiments establish the importance of these residues for efficient catalysis, but their role in catalysis is not clear. To examine the role of arginine residues in the phosphoryl transfer reaction, we have measured the consequences of mutations to arginine 166 in Escherichia coli alkaline phosphatase on hydrolysis of ethyl phosphate, on individual reaction steps in the hydrolysis of the covalent enzyme-phosphoryl intermediate, and on thio substitution effects. The results show that the role of the arginine side chain extends beyond its positive charge, as the Arg166Lys mutant is as compromised in activity as Arg166Ser. Through measurement of individual reaction steps, we construct a free energy profile for the hydrolysis of the enzyme-phosphate intermediate. This analysis indicates that the arginine side chain strengthens binding by approximately 3 kcal/mol and provides an additional 1-2 kcal/mol stabilization of the chemical transition state. A 2.1 A X-ray diffraction structure of Arg166Ser AP is presented, which shows little difference in enzyme structure compared to the wild-type enzyme but shows a significant reorientation of the bound phosphate. Altogether, these results support a model in which the arginine contributes to catalysis through binding interactions and through additional transition state stabilization that may arise from complementarity of the guanidinum group to the geometry of the trigonal bipyramidal transition state.     

Crystallographic analysis of a complex between human immunodeficiency virus type 1 protease and acetyl-pepstatin at 2.0-A resolution

The mode of binding of acetyl-pepstatin to the protease from the human immunodeficiency virus type 1 (HIV-1) has been determined by x-ray diffraction analysis. Crystals of an acetyl-pepstatin-HIV-1 protease complex were obtained in space group P2(1)2(1)2 (unit cell dimensions a = 58.39 A, b = 86.70 A, c = 46.27 A) by precipitation with sodium chloride. The structure was phased by molecular replacement methods, and a model for the structure was refined using diffraction data to 2.0 A resolution (R = 0.176 for 12901 reflections with I greater than sigma (I); deviation of bond distances from ideal values = 0.018 A; 172 solvent molecules included). The structure of the protein in the complex has been compared with the structure of the enzyme without the ligand. A core of 44 amino acids in each monomer, including residues in the active site and residues at the dimer interface, remains unchanged on binding of the inhibitor (root mean square deviation of alpha carbon positions = 0.39 A). The remaining 55 residues in each monomer undergo substantial rearrangement, with the most dramatic changes occurring at residues 44-57 (these residues comprise the so-called flaps of the enzyme). The flaps interact with one another and with the inhibitor so as to largely preserve the 2-fold symmetry of the protein. The inhibitor is bound in two approximately symmetric orientations. In both orientations the peptidyl backbone of the inhibitor is extended; a network of hydrogen bonds is formed between the inhibitor and the main body of the protein as well as between the inhibitor and the flaps. Hydrophobic side chains of residues in the body of the protein form partial binding sites for the side chains of the inhibitor; hydrophobic side chains of residues in the flaps complete these binding sites.     

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.     

Reconstruction by site-directed mutagenesis of the transition state for the activation of tyrosine by the tyrosyl-tRNA synthetase: a mobile loop envelopes the transition state in an induced-fit mechanism

Site-directed mutagenesis of the tyrosyl-tRNA synthetase followed by kinetic studies has shown that residues which are distant from the active site of the free enzyme are brought into play as the structure of the enzyme changes during catalysis. Positively charged side chains which are in mobile loops of the enzyme envelope the negatively charged pyrophosphate moiety during the transition state for the formation of tyrosyl adenylate in an induced-fit mechanism. Residues Lys-82 and Arg-86, which are on one side of the rim of the binding site pocket, and Lys-230 and Lys-233, which are on the other side, have been mutated to alanine residues and also to asparagine or glutamine. The resultant mutants still form 1 mol of tyrosyl adenylate/mol of dimer but with rate constants up to 8000 times lower. Construction of difference energy diagrams reveals that all the residues specifically interact with the transition state for the reaction and with pyrophosphate in the E.Tyr-AMP.PPi complex. Yet, the epsilon-NH3+ groups of Lys-230 and Lys-233 in the crystalline enzyme are at least 8 A too far away to interact with the pyrophosphate moiety in the transition state at the same time as do Lys-82 and Arg-86. Binding of substrates must, therefore, induce a conformational change in the enzyme that brings these residues into range. Consistent with this proposal is the observation that all four residues are in flexible regions of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Investigation of transition-state stabilization by residues histidine-45 and threonine-40 in the tyrosyl-tRNA synthetase

We have analyzed various mutations involving residues Thr-40 and His-45 in the tyrosyl-tRNA synthetase of Bacillus stearothermophilus. The utilization of binding energy in catalysis of tyrosyl adenylate formation from tyrosine and ATP was determined from the free energy profiles for the mutant enzymes. Our results confirm that the side chains of Thr-40 and His-45 provide a binding site for the pyrophosphoryl portion of the transition state of this reaction and for pyrophosphate in the reverse reaction. Deletion of these side chains destabilizes the transition-state by 4.9 and 4.1 kcal mol-1, respectively, consistent with a charged hydrogen-bonding interaction. To examine the role of His-45 further, we constructed the potentially conservative mutations His----Gln-45 and His----Asn-45. Both mutant enzymes are debilitated compared with the native enzyme. The His----Gln-45 enzyme is more active than enzyme in which the complete side chain is deleted (His----Ala-45), and so in this location glutamine is a semiconservative replacement. In contrast, the His----Asn-45 mutation is significantly worse than simple deletion of the side chain, indicating that asparagine at this position causes active destabilization of the transition state compared to His----Ala-45. The amide-NH2 of glutamine may be considered stereochemically equivalent to the epsilon-NH of histidine whereas the amide-NH2 of asparagine is comparable to the delta-NH of histidine. The results suggest that the epsilon-NH rather than the delta-NH group of His-45 is involved in the transition-state stabilization.(ABSTRACT TRUNCATED AT 250 WORDS)