Catalysing new reactions during evolution: economy of residues and mechanism

The diversity of function in some enzyme superfamilies shows that during evolution, enzymes have evolved to catalyse different reactions on the same structure scaffold. In this analysis, we examine in detail how enzymes can modify their chemistry, through a comparison of the catalytic residues and mechanisms in 27 pairs of homologous enzymes of totally different functions. We find that evolution is very economical. Enzymes retain structurally conserved residues to aid catalysis, including residues that bind catalytic metal ions and modulate cofactor chemistry. We examine the conservation of residue type and residue function in these structurally conserved residue pairs. Additionally, enzymes often retain common mechanistic steps catalyzed by structurally conserved residues. We have examined these steps in the context of their overall reactions.     

Substrate recognition and catalysis by UCH-L1

Deubiquitinating enzymes regulate essential cellular processes, and their dysregulation is implicated in multiple disease states. Ubiquitin carboxy-terminal hydrolase L1 (UCH-L1) has garnered attention for its links with Parkinson's disease and cancer; however, the mechanism of action of this enzyme in cells remains poorly understood. In order to advance our understanding of UCH-L1 function, we have been developing small molecule modulators of the enzyme for use as tools to probe its role in cells. In support of these efforts, an investigation of the mechanism of UCH-L1 catalysis was previously reported. Here, we extend this mechanistic evaluation and examine substrate recognition by UCH-L1. We developed a panel of ubiquitin fusions to test the contribution of specific residues of ubiquitin to binding and catalysis by the enzyme, and determined the activation parameters of selected variants to gain additional mechanistic insight. Ubiquitin side chains critical for establishing the Michaelis complex and enabling catalysis were identified, and features of this complex that differ between UCH-L1 and a homologue, UCH-L3, were revealed. These data provide dramatic examples of differences in substrate specificity between these enzymes. The implications of our experiments with UCH-L1 for selective inhibitor design and the relationship to disease are discussed.     

Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein''s designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Author''s Summary

Enzymes are nature''s molecular machines that catalyze biochemical reactions with remarkable efficiency. Recent evidence suggests that enzyme function may involve not only direct structural interactions between the enzyme and its substrate, but also internal motions of the enzyme itself. Here, we describe a computational investigation of three classes of enzymes that catalyze completely different biochemical reactions. Remarkably, the mobile enzyme regions and the nature of these motions are the same across species ranging from single-celled organisms to complex life-forms. Also surprisingly, non-homologous enzymes that catalyze the same chemical reaction but do not share sequence or structural similarity reveal a similar impact of enzyme motions on their reaction mechanisms. Flexible enzyme regions are found to be connected by conserved networks of coupled interactions that connect surface regions to active-site residues. These networks may provide a mechanism for the solvent on an enzyme''s surface to couple to the reaction catalyzed by the enzyme. These results have implications for understanding the mechanism of allostery (long-range effects), and for protein engineering and drug design.     

Structural and functional effects of mutations altering the subunit interface of mitochondrial malate dehydrogenase

Among highly conserved residues in eucaryotic mitochondrial malate dehydrogenases are those with roles in maintaining the interactions between identical monomeric subunits that form the dimeric enzymes. The contributions of two of these residues, Asp-43 and His-46, to structural stability and catalytic function were investigated by construction of mutant enzymes containing Asn-43 and Leu-46 substitutions using in vitro mutagenesis of the Saccharomyces cerevisiae gene (MDH1) encoding mitochondrial malate dehydrogenase. The mutant enzymes were expressed in and purified from a yeast strain containing a disruption of the chromosomal MDH1 locus. The enzyme containing the H46L substitution, as compared to the wild type enzyme, exhibits a dramatic shift in the pH profile for catalysis toward an optimum at low pH values. This shift corresponds with an increased stability of the dimeric form of the mutant enzyme, suggesting that His-46 may be the residue responsible for the previously described pH-dependent dissociation of mitochondrial malate dehydrogenase. The D43N substitution results in a mutant enzyme that is essentially inactive in in vitro assays and that tends to aggregate at pH 7.5, the optimal pH for catalysis for the dimeric wild type enzyme.     

Conservation of electrostatic properties within enzyme families and superfamilies

Electrostatic interactions play a key role in enzyme catalytic function. At long range, electrostatics steer the incoming ligand/substrate to the active site, and at short distances, electrostatics provide the specific local interactions for catalysis. In cases in which electrostatics determine enzyme function, orthologs should share the electrostatic properties to maintain function. Often, electrostatic potential maps are employed to depict how conserved surface electrostatics preserve function. We expand on previous efforts to explain conservation of function, using novel electrostatic sequence and structure analyses of four enzyme families and one enzyme superfamily. We show that the spatial charge distribution is conserved within each family and superfamily. Conversely, phylogenetic analysis of key electrostatic residues provide the evolutionary origins of functionality.     

Yeast epiarginase regulation,an enzyme-enzyme activity control: identification of residues of ornithine carbamoyltransferase and arginase responsible for enzyme catalytic and regulatory activities

In the presence of ornithine and arginine, ornithine carbamoyltransferase (OTCase) and arginase form a one-to-one enzyme complex in which the activity of OTCase is inhibited whereas arginase remains catalytically active. The mechanism by which these nonallosteric enzymes form a stable complex triggered by the binding of their respective substrates raises the question of how such a cooperative association is induced. Analyses of mutations in both enzymes identify residues that are required for their association, some of them being important for catalysis. In arginase, two cysteines at the C terminus of the protein are crucial for its epiarginase function but not for its catalytic activity and trimeric structure. In OTCase, mutations of putative ornithine binding residues, Asp-182, Asn-184, Asn-185, Cys-289, and Glu-256 greatly reduced the affinity for ornithine and impaired the interaction with arginase. The four lysine residues located in the SMG loop, Lys-260, Lys-263, Lys-265, and Lys-268, also play an important role in mediating the sensitivity of OTCase to ornithine and to arginase and appear to be involved in transducing and enhancing the signal given by ornithine for the closure of the catalytic domain.     

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.     

Heterologous expression of the vanadium-containing chloroperoxidase from Curvularia inaequalis in Saccharomyces cerevisiae and site-directed mutagenesis of the active site residues His(496), Lys(353), Arg(360), and Arg(490).

The vanadium-containing chloroperoxidase from the fungus Curvularia inaequalis is heterologously expressed to high levels in the yeast Saccharomyces cerevisiae. Characterization of the recombinant enzyme reveals that this behaves very similar to the native chloroperoxidase. Site-directed mutagenesis is performed on four highly conserved active site residues to examine their role in catalysis. When the vanadate-binding residue His(496) is changed into an alanine, the mutant enzyme loses the ability to bind vanadate covalently resulting in an inactive enzyme. The negative charges on the vanadate oxygens are compensated by hydrogen bonds with the residues Arg(360), Arg(490), and Lys(353). When these residues are changed into alanines the mutant enzymes lose the ability to effectively oxidize chloride but can still function as bromoperoxidases. A general mechanism for haloperoxidase catalysis is proposed that also correlates the kinetic properties of the mutants with the charge and the hydrogen-bonding network in the vanadate-binding site.     

Mechanisms of enzyme action and inhibition: Transition state analogues for acid-base catalysis

The theory of absolute reaction rates implies that the grip of a catalyst on a substrate tightens with substrate activation, relaxing later as products are formed and released. Analogs that mimic different kinds of substrate activation can, through the structural details of their complexes with enzymes, indicate how active site residues are involved in the enhancement of reactions rates. In several cases, bonds involved in general acid-base catalysis have been identified tentatively; and recent evidence points to a hydrogen bond of remarkable stability in the transition state in enzymatic deamination of adenosine. Similar approaches have been used to enzymes that act primarily by substrate distortion, nucleophilic catalysis, solvent removal and catalysis by approximation. Two recurring observations, that were not expected in theory, have been the binding of inhibitors in ionized forms that are rare in solution, and changes in enzyme configuration that accompany binding of transition state analogs. Origins and implications of these findings will be discussed with specific reference to the role of solvent water in catalytic phenomena.Supported by Grant GM-18325 from the National Institutes of Health.     

Monoterpene and sesquiterpene synthases and the origin of terpene skeletal diversity in plants

The multitude of terpene carbon skeletons in plants is formed by enzymes known as terpene synthases. This review covers the monoterpene and sesquiterpene synthases presenting an up-to-date list of enzymes reported and evidence for their ability to form multiple products. The reaction mechanisms of these enzyme classes are described, and information on how terpene synthase proteins mediate catalysis is summarized. Correlations between specific amino acid motifs and terpene synthase function are described, including an analysis of the relationships between active site sequence and cyclization type and a discussion of whether specific protein features might facilitate multiple product formation.     

Protein dynamics and enzymatic catalysis: investigating the peptidyl-prolyl cis-trans isomerization activity of cyclophilin A

A growing body of evidence suggests a connection between protein dynamics and enzymatic catalysis. In this paper, we present a variety of computational studies designed to investigate the role of protein dynamics in the detailed mechanism of peptidyl-prolyl cis-trans isomerization catalyzed by human cyclophilin A. The results identify a network of protein vibrations, extending from surface regions of the enzyme to the active site and coupled to substrate turnover. Indications are that this network may have a role in promoting catalysis. Crucial parts of this network are found to be conserved in 10 cyclophilin structures from six different species. Experimental evidence for the existence of this network comes from previous NMR relaxation studies, where motions in several residues, forming parts of this network, were detected only during substrate turnover. The high temperature factors (from X-ray crystal structures) associated with the network residues provide further evidence of these vibrations. Along with the knowledge of enzyme structure, this type of network could provide new insights into enzymatic catalysis and the effect of distant ligand binding on protein function. The procedure outlined in this paper is general and can be applied to other enzymatic systems as well. This presents an interesting opportunity; collaborative experimental and theoretical investigations designed to characterize in detail the nature and function of this type of network could enhance the understanding of protein dynamics in enzymatic catalysis.     

Amphipatic molecules affect the kinetic profile of Pseudomonas putida chlorocatechol 1,2-dioxygenase

Dioxygenases are nonheme iron enzymes that biodegrade recalcitrant compounds, such as catechol and derivatives, released into the environment by modern industry. Intradiol dioxygenases have attracted much attention due to the interest in their use for bioremediation, which has demanded efforts towards understanding their action mechanism and also how to control it. The role of unexpected amphipatic molecules, observed in crystal structures of intradiol dioxygenases, during catalysis has been poorly explored. We report results obtained with the intradiol enzyme chlorocatechol 1,2-dioxygenase (1,2-CCD) from Pseudomonas putida subjected to delipidation. The delipidated enzyme is more stable and shows more cooperative thermal denaturation. The kinetics changes from Michaelis–Menten to a cooperative scheme, indicating that conformational changes propagate between monomers in the absence of amphipatic molecules. Furthermore, these molecules inhibit catalysis, yielding lower v _max values. To the best of our knowledge, this is the first report concerning the effects of amphipatic molecules on 1,2-CCD function.     

The guanidino-group modifying enzymes: structural basis for their diversity and commonality

The guanidino-group modifying enzyme (GME) superfamily contains many drug targets, including metabolic enzymes from pathogenic microorganisms as well as key regulatory proteins from higher eukaryotes. These enzymes, despite their diverse sequences, adopt the common alpha/beta propeller fold and catalyze the modification of (methylated) guanidino groups. Our structural superposition and structure-based alignment for the GMEs have identified key residues that are involved in the catalysis and substrate binding. We have shown that conserved guanidino-carboxyl interactions are utilized in two different ways; the acidic residues in the catalytic site form hydrogen bonds to the substrate guanidino group, and the enzyme Arg residues at several key positions recognize the carboxyl group of the substrate and fix its orientation. Based on this observation, we have proposed rules for classifying the GME sequences and predicting their molecular function from the conservation of the key acidic and Arg residues. Other novel motifs have been identified, which involve residues that are not in direct contact with the substrate but are likely to stabilize the active-site conformation through hydrogen-bonding networks. In addition, we have examined the domain architecture of the GMEs. Although most members consist of a single catalytic domain, fold recognition analysis has identified a likely bifunctional enzyme from a cyanobacterium. It has also revealed common immunoglobulin-like beta-sandwich domains found in the enzymes that recognize protein substrates. These findings will be useful for predicting the precise mechanism of action for potential novel targets and designing therapeutic compounds against them.     

Understanding the Functional Roles of Amino Acid Residues in Enzyme Catalysis

The MACiE database contains 223 distinct step-wise enzyme reaction mechanisms and holds representatives from each EC sub-subclass where there is a crystal structure and sufficient evidence in the literature to support 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 mechanism by which substrates are transformed into products. Using MACiE as a knowledge base, we have seen that the top 10 most catalytic residues are histidine, aspartate, glutamate, lysine, cysteine, arginine, serine, threonine, tyrosine and tryptophan. Of these only seven (cysteine, histidine, aspartate, lysine, serine, threonine and tyrosine) dominate catalysis and provide essentially five functional roles that are essential. Stabilisation is the most common and essential role for all classes of enzyme, followed by general acid/base (proton acceptor and proton donor) functionality, with nucleophilic addition following closely behind (nucleophile and nucleofuge). We investigated the occurrence of these residues in MACiE and the Catalytic Site Atlas and found that, as expected, certain residue types are associated with each functional role, with some residue types able to perform diverse roles. In addition, it was seen that different EC classes of enzyme have a tendency to employ different residues for catalysis. Further, we show that whilst the differences between EC classes in catalytic residue composition are not immediately obvious from the general classes of Ingold mechanisms, there is some weak correlation between the mechanisms involved in a given EC class and the functions that the catalytic amino acid residues are performing. The analysis presented here provides a valuable insight into the functional roles of catalytic amino acid residues, which may have applications in many aspects of enzymology, from the design of novel enzymes to the prediction and validation of enzyme reaction mechanisms.     

Mechanism of proton transfer in [FeFe]-hydrogenase from Clostridium pasteurianum

[FeFe]-Hydrogenases are complex metalloproteins that catalyze the reversible reduction of protons to molecular hydrogen utilizing a unique diiron subcluster bridged to a [4Fe4S] subcluster. Extensive studies have concentrated on the nature and catalytic activity of the active site, yet relatively little information is available concerning the mechanism of proton transport that is required for this activity. Previously, structural characterization of [FeFe]-hydrogenase from Clostridium pasteurianum indicated a potential proton transport pathway involving four residues (Cys-299, Glu-279, Ser-319, and Glu-282) that connect the active site to the enzyme surface. Here, we demonstrate that substitution of any of these residues resulted in a drastic reduction in hydrogenase activity relative to the native enzyme, supporting the importance of these residues in catalysis. Inhibition studies of native and amino acid-substituted enzymes revealed that Zn(2+) specifically blocked proton transfer by binding to Glu-282, confirming the role of this residue in the identified pathway. In addition, all four of these residues are strictly conserved, suggesting that they may form a proton transport pathway that is common to all [FeFe]-hydrogenases.     

Residues Required for Activity in Escherichia coli o-Succinylbenzoate Synthase (OSBS) Are Not Conserved in All OSBS Enzymes

Understanding how enzyme specificity evolves will provide guiding principles for protein engineering and function prediction. The o-succinylbenzoate synthase (OSBS) family is an excellent model system for elucidating these principles because it has many highly divergent amino acid sequences that are <20% identical, and some members have evolved a second function. The OSBS family belongs to the enolase superfamily, members of which use a set of conserved residues to catalyze a wide variety of reactions. These residues are the only conserved residues in the OSBS family, so they are not sufficient to determine reaction specificity. Some enzymes in the OSBS family catalyze another reaction, N-succinylamino acid racemization (NSAR). NSARs cannot be segregated into a separate family because their sequences are highly similar to those of known OSBSs, and many of them have both OSBS and NSAR activities. To determine how such divergent enzymes can catalyze the same reaction and how NSAR activity evolved, we divided the OSBS family into subfamilies and compared the divergence of their active site residues. Correlating sequence conservation with the effects of mutations in Escherichia coli OSBS identified two nonconserved residues (R159 and G288) at which mutations decrease efficiency ≥200-fold. These residues are not conserved in the subfamily that includes NSAR enzymes. The OSBS/NSAR subfamily binds the substrate in a different orientation, eliminating selective pressure to retain arginine and glycine at these positions. This supports the hypothesis that specificity-determining residues have diverged in the OSBS family and provides insight into the sequence changes required for the evolution of NSAR activity.     

Kinetic properties of Streptococcus pneumoniae hyaluronate lyase

Streptococcus pneumoniae hyaluronate lyase is a surface antigen of this bacterial pathogen, which causes significant mortality and morbidity in human populations worldwide. The primary function of this enzyme is the degradation of hyaluronan, a major component of the extracellular matrix of the tissues of practically all vertebrates. The enzyme uses a processive mode of action to degrade hyaluronan to a final product, an unsaturated disaccharide hyaluronan unit. This catalysis proceeds via a five-step proton acceptance and donation mechanism that includes substrate binding, catalysis, release of the disaccharide product, translocation of the remaining hyaluronan substrate, and proton exchange with microenvironment. Based on the analysis of the three-dimensional structure of the native enzyme and its complexes with hexasaccharide substrate and disaccharide product, several residues have been chosen for mutation studies. These mutated residues included the catalytic residues Asn349, His399, Tyr408, and residues responsible for substrate binding and translocation, Arg243 and Asn580. The comparison of the kinetic properties of the wild-type with the mutant enzymes allowed for the characterization of every mutant and the correlation of the kinetic properties of the enzyme with its structure. The comparison of the wild-type hyaluronate lyase with other polysaccharide-degrading enzymes, the hydrolases endonuclease and glucoamylase, shows striking similarity of K(m)s for all of these different enzymes.     

Evidence for distinct roles in catalysis for residues of the serine-serine-lysine catalytic triad of fatty acid amide hydrolase

Fatty acid amide hydrolase (FAAH) is a mammalian amidase signature enzyme that inactivates neuromodulatory fatty acid amides, including the endogenous cannabinoid anandamide and the sleep-inducing substance oleamide. The recent determination of the three-dimensional structures of FAAH and two distantly related bacterial amidase signature enzymes indicates that these enzymes employ an unusual serine-serine-lysine triad for catalysis (Ser-241/Ser-217/Lys-142 in FAAH). Mutagenesis of each of the triad residues in FAAH has been shown to severely reduce amidase activity; however, how these residues contribute, both individually and in cooperation, to catalysis remains unclear. Here, through a combination of site-directed mutagenesis, enzyme kinetics, and chemical labeling experiments, we provide evidence that each FAAH triad residue plays a distinct role in catalysis. In particular, the mutation of Lys-142 to alanine indicates that this residue functions as both a base involved in the activation of the Ser-241 nucleophile and an acid that participates in the protonation of the substrate leaving group. This latter property appears to support the unusual ability of FAAH to hydrolyze amides and esters at equivalent rates. Interestingly, although structural evidence indicates that the impact of Lys-142 on catalysis probably occurs through the bridging Ser-217, the mutation of this latter residue to alanine impaired catalytic activity but left the amide/ester hydrolysis ratios of FAAH intact. Collectively, these findings suggest that FAAH possesses a specialized active site structure dedicated to a mechanism for competitive amide and ester hydrolysis where nucleophile attack and leaving group protonation occur in a coordinated manner dependent on Lys-142.     

Metal ions in biological catalysis: from enzyme databases to general principles

We analysed the roles and distribution of metal ions in enzymatic catalysis using available public databases and our new resource Metal-MACiE (). In Metal-MACiE, a database of metal-based reaction mechanisms, 116 entries covering 21% of the metal-dependent enzymes and 70% of the types of enzyme-catalysed chemical transformations are annotated according to metal function. We used Metal-MACiE to assess the functions performed by metals in biological catalysis and the relative frequencies of different metals in different roles, which can be related to their individual chemical properties and availability in the environment. The overall picture emerging from the overview of Metal-MACiE is that redox-inert metal ions are used in enzymes to stabilize negative charges and to activate substrates by virtue of their Lewis acid properties, whereas redox-active metal ions can be used both as Lewis acids and as redox centres. Magnesium and zinc are by far the most common ions of the first type, while calcium is relatively less used. Magnesium, however, is most often bound to phosphate groups of substrates and interacts with the enzyme only transiently, whereas the other metals are stably bound to the enzyme. The most common metal of the second type is iron, which is prevalent in the catalysis of redox reactions, followed by manganese, cobalt, molybdenum, copper and nickel. The control of the reactivity of redox-active metal ions may involve their association with organic cofactors to form stable units. This occurs sometimes for iron and nickel, and quite often for cobalt and molybdenum. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Effect of dimer dissociation on activity and thermostability of the alpha-glucuronidase from Geobacillus stearothermophilus: dissecting the different oligomeric forms of family 67 glycoside hydrolases

The oligomeric organization of enzymes plays an important role in many biological processes, such as allosteric regulation, conformational stability and thermal stability. alpha-Glucuronidases are family 67 glycosidases that cleave the alpha-1,2-glycosidic bond between 4-O-methyl-D-glucuronic acid and xylose units as part of an array of hemicellulose-hydrolyzing enzymes. Currently, two crystal structures of alpha-glucuronidases are available, those from Geobacillus stearothermophilus (AguA) and from Cellvibrio japonicus (GlcA67A). Both enzymes are homodimeric, but surprisingly their dimeric organization is different, raising questions regarding the significance of dimerization for the enzymes' activity and stability. Structural comparison of the two enzymes suggests several elements that are responsible for the different dimerization organization. Phylogenetic analysis shows that the alpha-glucuronidases AguA and GlcA67A can be classified into two distinct subfamilies of bacterial alpha-glucuronidases, where the dimer-forming residues of each enzyme are conserved only within its own subfamily. It seems that the different dimeric forms of AguA and GlcA67A represent the two alternative dimeric organizations of these subfamilies. To study the biological significance of the dimerization in alpha-glucuronidases, we have constructed a monomeric form of AguA by mutating three of its interface residues (W328E, R329T, and R665N). The activity of the monomer was significantly lower than the activity of the wild-type dimeric AguA, and the optimal temperature for activity of the monomer was around 35 degrees C, compared to 65 degrees C of the wild-type enzyme. Nevertheless, the melting temperature of the monomeric protein, 72.9 degrees C, was almost identical to that of the wild-type, 73.4 degrees C. It appears that the dimerization of AguA is essential for efficient catalysis and that the dissociation into monomers results in subtle conformational changes in the structure which indirectly influence the active site region and reduce the activity. Structural and mechanistic explanations for these effects are discussed.