Stability for function trade-offs in the enolase superfamily "catalytic module"

Enzyme catalysis reflects a dynamic interplay between charged and polar active site residues that facilitate function, stabilize transition states, and maintain overall protein stability. Previous studies show that substituting neutral for charged residues in the active site often significantly stabilizes a protein, suggesting a stability trade-off for functionality. In the enolase superfamily, a set of conserved active site residues (the "catalytic module") has repeatedly been used in nature in the evolution of many different enzymes for the performance of unique overall reactions involving a chemically diverse set of substrates. This catalytic module provides a robust solution for catalysis that delivers the common underlying partial reaction that supports all of the different overall chemical reactions of the superfamily. As this module has been so broadly conserved in the evolution of new functions, we sought to investigate the extent to which it follows the stability-function trade-off. Alanine substitutions were made for individual residues, groups of residues, and the entire catalytic module of o-succinylbenzoate synthase (OSBS), a member of the enolase superfamily from Escherichia coli. Of six individual residue substitutions, four (K131A, D161A, E190A, and D213A) substantially increased protein stability (by 0.46-4.23 kcal/mol), broadly consistent with prediction of a stability-activity trade-off. The residue most conserved across the superfamily, E190, is by far the most destabilizing. When the individual substitutions were combined into groups (as they are structurally and functionally organized), nonadditive stability effects emerged, supporting previous observations that residues within the module interact as two functional groups within a larger catalytic system. Thus, whereas the multiple-mutant enzymes D161A/E190A/D213A and K131A/K133A/D161A/E190A/D213A/K235A (termed 3KDED) are stabilized relative to the wild-type enzyme (by 1.77 and 3.68 kcal/mol, respectively), the net stabilization achieved in both cases is much weaker than what would be predicted if their stability contributions were additive. Organization of the catalytic module into systems that mitigate the expected stability cost due to the presence of highly charged active site residues may help to explain its repeated use for the evolution of many different functions.     

The evolutionary origins and catalytic importance of conserved electrostatic networks within TIM-barrel proteins

Conservation of function is the basic tenet of protein evolution. Conservation of key electrostatic properties is a frequently employed mechanism that leads to conserved function. In a previous report, we identified several conserved electrostatic properties in four protein families and one functionally diverse enzyme superfamily. In this report, we demonstrate the evolutionary and catalytic importance of electrostatic networks in three ubiquitous metabolic enzymes: triosephosphate isomerase, enolase, and transaldolase. Evolutionary importance is demonstrated using phylogenetic motifs (sequence fragments that parallel the overall familial phylogeny). Phylogenetic motifs frequently correspond to both catalytic residues and conserved interactions that fine-tune catalytic residue pKa values. Further, in the case of triosephosphate isomerase, quantitative differences in the catalytic Glu169 pKa values parallel subfamily differentiation. Finally, phylogenetic motifs are shown to structurally cluster around the active sites of eight different TIM-barrel families. Depending upon the mechanistic requisites of each reaction catalyzed, interruptions to the canonical fold may or may not be identified as phylogenetic motifs.     

Persistently conserved positions in structurally similar, sequence dissimilar proteins: Roles in preserving protein fold and function

Many protein pairs that share the same fold do not have any detectable sequence similarity, providing a valuable source of information for studying sequence-structure relationship. In this study, we use a stringent data set of structurally similar, sequence-dissimilar protein pairs to characterize residues that may play a role in the determination of protein structure and/or function. For each protein in the database, we identify amino-acid positions that show residue conservation within both close and distant family members. These positions are termed "persistently conserved". We then proceed to determine the "mutually" persistently conserved (MPC) positions: those structurally aligned positions in a protein pair that are persistently conserved in both pair mates. Because of their intra- and interfamily conservation, these positions are good candidates for determining protein fold and function. We find that 45% of the persistently conserved positions are mutually conserved. A significant fraction of them are located in critical positions for secondary structure determination, they are mostly buried, and many of them form spatial clusters within their protein structures. A substitution matrix based on the subset of MPC positions shows two distinct characteristics: (i) it is different from other available matrices, even those that are derived from structural alignments; (ii) its relative entropy is high, emphasizing the special residue restrictions imposed on these positions. Such a substitution matrix should be valuable for protein design experiments.     

Force-clamp spectroscopy detects residue co-evolution in enzyme catalysis

Understanding how the catalytic mechanisms of enzymes are optimized through evolution remains a major challenge in molecular biology. The concept of co-evolution implicates that compensatory mutations occur to preserve the structure and function of proteins. We have combined statistical analysis of protein sequences with the sensitivity of single molecule force-clamp spectroscopy to probe how catalysis is affected by structurally distant correlated mutations in Escherichia coli thioredoxin. Our findings show that evolutionary anti-correlated mutations have an inhibitory effect on enzyme catalysis, whereas positively correlated mutations rescue the catalytic activity. We interpret these results in terms of an evolutionary tuning of both the enzyme-substrate binding process and the chemistry of the active site. Our results constitute a direct observation of distant residue co-evolution in enzyme catalysis.     

Contribution of active site residues to substrate hydrolysis by USP2: insights into catalysis by ubiquitin specific proteases

The ubiquitin-specific protease (USP) structural class represents the largest and most diverse family of deubiquitinating enzymes (DUBs). Many USPs assume important biological roles and emerge as potential targets for therapeutic intervention. A clear understanding of USP catalytic mechanism requires a functional evaluation of the proposed key active site residues. Crystallographic data of ubiquitin aldehyde adducts of USP catalytic cores provided structural details on the catalytic triad residues, namely the conserved Cys and His, and a variable putative third residue, and inferred indirect structural roles for two other conserved residues (Asn and Asp), in stabilizing via a bridging water molecule the oxyanion of the tetrahedral intermediate (TI). We have expressed the catalytic domain of USP2 and probed by site-directed mutagenesis the role of these active site residues in the hydrolysis of peptide and isopeptide substrates, including a synthetic K48-linked diubiquitin substrate for which a label-free, mass spectrometry based assay has been developed to monitor cleavage. Hydrolysis of ubiquitin-AMC, a model substrate, was not affected by the mutations. Molecular dynamics simulations of USP2, free and complexed with the TI of a bona fide isopeptide substrate, were carried out. We found that Asn271 is structurally poised to directly stabilize the oxyanion developed in the acylation step, while being structurally supported by the adjacent absolutely conserved Asp575. Mutagenesis data functionally confirmed this structural role independent of the nature (isopeptide vs peptide) of the bond being cleaved. We also found that Asn574, structurally located as the third member of the catalytic triad, does not fulfill this role functionally. A dual supporting role is inferred from double-point mutation and structural data for the absolutely conserved residue Asp575, in oxyanion hole formation, and in maintaining the correct alignment and protonation of His557 for catalytic competency.     

Enzyme-mononucleotide interactions: three different folds share common structural elements for ATP recognition

Three ATP-dependent enzymes with different folds, cAMP-dependent protein kinase, D-Ala:D-Ala ligase and the alpha-subunit of the alpha2beta2 ribonucleotide reductase, have a similar organization of their ATP-binding sites. The most meaningful similarity was found over 23 structurally equivalent residues in each protein and includes three strands each from their beta-sheets, in addition to a connecting loop. The equivalent secondary structure elements in each of these enzymes donate four amino acids forming key hydrogen bonds responsible for the common orientation of the "AMP" moieties of their ATP-ligands. One lysine residue conserved throughout the three families binds the alpha-phosphate in each protein. The common fragments of structure also position some, but not all, of the equivalent residues involved in hydrophobic contacts with the adenine ring. These examples of convergent evolution reinforce the view that different proteins can fold in different ways to produce similar structures locally, and nature can take advantage of these features when structure and function demand it, as shown here for the common mode of ATP-binding by three unrelated proteins.     

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.     

Distinctive structural motifs co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole in the alpha/beta‐hydrolase fold enzymes

The alpha/beta‐hydrolases (ABH) are among the largest structural families of proteins that are found in nature. Although they vary in their sequence and function, the ABH enzymes use a similar acid–base‐nucleophile catalytic mechanism to catalyze reactions on different substrates. Because ABH enzymes are biocatalysts with a wide range of potential applications, protein engineering has taken advantage of their catalytic versatility to develop enzymes with industrial applications. This study is a comprehensive analysis of 40 ABH enzyme families focusing on two identified substructures: the nucleophile zone and the oxyanion zone, which co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole, and independently reported as critical for the enzymatic activity. We also frequently observed an aromatic cluster near the nucleophile and oxyanion zones, and opposite the ligand‐binding site. The nucleophile zone, the oxyanion zone and the residue cluster enriched in aromatic side chains comprise a three‐dimensional structural organization that shapes the active site of ABH enzymes and plays an important role in the enzymatic function by structurally stabilizing the catalytic nucleophile and the residues of the oxyanion hole. The structural data support the notion that the aromatic cluster can participate in co‐ordination of the catalytic histidine loop, and properly place the catalytic histidine next to the catalytic nucleophile.     

Composite active site of chondroitin lyase ABC accepting both epimers of uronic acid

Enzymes have evolved as catalysts with high degrees of stereospecificity. When both enantiomers are biologically important, enzymes with two different folds usually catalyze reactions with the individual enantiomers. In rare cases a single enzyme can process both enantiomers efficiently, but no molecular basis for such catalysis has been established. The family of bacterial chondroitin lyases ABC comprises such enzymes. They can degrade both chondroitin sulfate (CS) and dermatan sulfate (DS) glycosaminoglycans at the nonreducing end of either glucuronic acid (CS) or its epimer iduronic acid (DS) by a beta-elimination mechanism, which commences with the removal of the C-5 proton from the uronic acid. Two other structural folds evolved to perform these reactions in an epimer-specific fashion: (alpha/alpha)(5) for CS (chondroitin lyases AC) and beta-helix for DS (chondroitin lyases B); their catalytic mechanisms have been established at the molecular level. The structure of chondroitinase ABC from Proteus vulgaris showed surprising similarity to chondroitinase AC, including the presence of a Tyr-His-Glu-Arg catalytic tetrad, which provided a possible mechanism for CS degradation but not for DS degradation. We determined the structure of a distantly related Bacteroides thetaiotaomicron chondroitinase ABC to identify additional structurally conserved residues potentially involved in catalysis. We found a conserved cluster located approximately 12 A from the catalytic tetrad. We demonstrate that a histidine in this cluster is essential for catalysis of DS but not CS. The enzyme utilizes a single substrate-binding site while having two partially overlapping active sites catalyzing the respective reactions. The spatial separation of the two sets of residues suggests a substrate-induced conformational change that brings all catalytically essential residues close together.     

Protein Topology Determines Cysteine Oxidation Fate: The Case of Sulfenyl Amide Formation among Protein Families

Cysteine residues have a rich chemistry and play a critical role in the catalytic activity of a plethora of enzymes. However, cysteines are susceptible to oxidation by Reactive Oxygen and Nitrogen Species, leading to a loss of their catalytic function. Therefore, cysteine oxidation is emerging as a relevant physiological regulatory mechanism. Formation of a cyclic sulfenyl amide residue at the active site of redox-regulated proteins has been proposed as a protection mechanism against irreversible oxidation as the sulfenyl amide intermediate has been identified in several proteins. However, how and why only some specific cysteine residues in particular proteins react to form this intermediate is still unknown. In the present work using in-silico based tools, we have identified a constrained conformation that accelerates sulfenyl amide formation. By means of combined MD and QM/MM calculation we show that this conformation positions the NH backbone towards the sulfenic acid and promotes the reaction to yield the sulfenyl amide intermediate, in one step with the concomitant release of a water molecule. Moreover, in a large subset of the proteins we found a conserved beta sheet-loop-helix motif, which is present across different protein folds, that is key for sulfenyl amide production as it promotes the previous formation of sulfenic acid. For catalytic activity, in several cases, proteins need the Cysteine to be in the cysteinate form, i.e. a low pK_a Cys. We found that the conserved motif stabilizes the cysteinate by hydrogen bonding to several NH backbone moieties. As cysteinate is also more reactive toward ROS we propose that the sheet-loop-helix motif and the constraint conformation have been selected by evolution for proteins that need a reactive Cys protected from irreversible oxidation. Our results also highlight how fold conservation can be correlated to redox chemistry regulation of protein function.     

The evolution of function in strictosidine synthase-like proteins

The exponential growth of sequence data provides abundant information for the discovery of new enzyme reactions. Correctly annotating the functions of highly diverse proteins can be difficult, however, hindering use of this information. Global analysis of large superfamilies of related proteins is a powerful strategy for understanding the evolution of reactions by identifying catalytic commonalities and differences in reaction and substrate specificity, even when only a few members have been biochemically or structurally characterized. A comparison of >2500 sequences sharing the six-bladed β-propeller fold establishes sequence, structural, and functional links among the three subgroups of the functionally diverse N6P superfamily: the arylesterase-like and senescence marker protein-30/gluconolactonase/luciferin-regenerating enzyme-like (SGL) subgroups, representing enzymes that catalyze lactonase and related hydrolytic reactions, and the so-called strictosidine synthase-like (SSL) subgroup. Metal-coordinating residues were identified as broadly conserved in the active sites of all three subgroups except for a few proteins from the SSL subgroup, which have been experimentally determined to catalyze the quite different strictosidine synthase (SS) reaction, a metal-independent condensation reaction. Despite these differences, comparison of conserved catalytic features of the arylesterase-like and SGL enzymes with the SSs identified similar structural and mechanistic attributes between the hydrolytic reactions catalyzed by the former and the condensation reaction catalyzed by SS. The results also suggest that despite their annotations, the great majority of these >500 SSL sequences do not catalyze the SS reaction; rather, they likely catalyze hydrolytic reactions typical of the other two subgroups instead. This prediction was confirmed experimentally for one of these proteins.     

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.     

Anatomy of the E2 ligase fold: implications for enzymology and evolution of ubiquitin/Ub-like protein conjugation

The configuration of the active site of E2 ligases, central enzymes in the ubiquitin/ubiquitin-like protein (Ub/Ubl) conjugation systems, has long puzzled researchers. Taking advantage of the wealth of newly available structures and sequences of E2s from diverse organisms, we performed a large-scale comparative analysis of these proteins. As a result we identified a previously under-appreciated diversity in the active site of these enzymes, in particular, the spatial location of the catalytic cysteine and a constellation of associated conserved residues that potentially contributes to catalysis. We observed structural innovations of differing magnitudes occurring in various families across the E2 fold that might correlate in part with differences in target interaction. A key finding was the independent emergence on multiple occasions of a polar residue, often a histidine, in the vicinity of the catalytic cysteine in different E2 families. We propose that these convergently emerging polar residues have a common function, such as in the stabilization of oxyanion holes during Ub/Ubl transfer and spatial localization of the Ub/Ubl tails in the active site. Thus, the E2 ligases represent a rare example in enzyme evolution of high structural diversity of the active site and position of the catalytic residue despite all characterized members catalyzing a similar reaction. Our studies also indicated certain evolutionarily conserved features in all active members of the E2 superfamily that stabilize the unusual flap-like structure in the fold. These features are likely to form a critical mechanical element of the fold required for catalysis. The results presented here could aid in new experiments to understand E2 catalysis.     

Mutational and crystallographic analyses of the active site residues of the Bacillus circulans xylanase.

Using site-directed mutagenesis we have investigated the catalytic residues in a xylanase from Bacillus circulans. Analysis of the mutants E78D and E172D indicated that mutations in these conserved residues do not grossly alter the structure of the enzyme and that these residues participate in the catalytic mechanism. We have now determined the crystal structure of an enzyme-substrate complex to 108 A resolution using a catalytically incompetent mutant (E172C). In addition to the catalytic residues, Glu 78 and Glu 172, we have identified 2 tyrosine residues, Tyr 69 and Tyr 80, which likely function in substrate binding, and an arginine residue, Arg 112, which plays an important role in the active site of this enzyme. On the basis of our work we would propose that Glu 78 is the nucleophile and that Glu 172 is the acid-base catalyst in the reaction.     

Predicting and annotating catalytic residues: an information theoretic approach.

We introduce a computational method to predict and annotate the catalytic residues of a protein using only its sequence information, so that we describe both the residues' sequence locations (prediction) and their specific biochemical roles in the catalyzed reaction (annotation). While knowing the chemistry of an enzyme's catalytic residues is essential to understanding its function, the challenges of prediction and annotation have remained difficult, especially when only the enzyme's sequence and no homologous structures are available. Our sequence-based approach follows the guiding principle that catalytic residues performing the same biochemical function should have similar chemical environments; it detects specific conservation patterns near in sequence to known catalytic residues and accordingly constrains what combination of amino acids can be present near a predicted catalytic residue. We associate with each catalytic residue a short sequence profile and define a Kullback-Leibler (KL) distance measure between these profiles, which, as we show, effectively captures even subtle biochemical variations. We apply the method to the class of glycohydrolase enzymes. This class includes proteins from 96 families with very different sequences and folds, many of which perform important functions. In a cross-validation test, our approach correctly predicts the location of the enzymes' catalytic residues with a sensitivity of 80% at a specificity of 99.4%, and in a separate cross-validation we also correctly annotate the biochemical role of 80% of the catalytic residues. Our results compare favorably to existing methods. Moreover, our method is more broadly applicable because it relies on sequence and not structure information; it may, furthermore, be used in conjunction with structure-based methods.     

A matching algorithm for catalytic residue site selection in computational enzyme design

A loop closure-based sequential algorithm, PRODA_MATCH, was developed to match catalytic residues onto a scaffold for enzyme design in silico. The computational complexity of this algorithm is polynomial with respect to the number of active sites, the number of catalytic residues, and the maximal iteration number of cyclic coordinate descent steps. This matching algorithm is independent of a rotamer library that enables the catalytic residue to take any required conformation during the reaction coordinate. The catalytic geometric parameters defined between functional groups of transition state (TS) and the catalytic residues are continuously optimized to identify the accurate position of the TS. Pseudo-spheres are introduced for surrounding residues, which make the algorithm take binding into account as early as during the matching process. Recapitulation of native catalytic residue sites was used as a benchmark to evaluate the novel algorithm. The calculation results for the test set show that the native catalytic residue sites were successfully identified and ranked within the top 10 designs for 7 of the 10 chemical reactions. This indicates that the matching algorithm has the potential to be used for designing industrial enzymes for desired reactions.     

The geometry of interactions between catalytic residues and their substrates

The process of deducing the catalytic mechanism of an enzyme from its structure is highly complex and requires extensive experimental work to validate a proposed mechanism. As one step towards improving the reliability of this process, we have gathered statistics describing the typical geometry of catalytic residues with regard to the substrate and one another. In order to analyse residue-substrate interactions, we have assembled a dataset of structures of enzymes of known mechanism bound to substrate, product, or a substrate analogue. Despite the challenges presented in obtaining such experimental data, we were able to include 42 enzyme structures. We have also assembled a separate dataset of catalytic residues which act upon other catalytic residues, using a set of 60 enzyme structures. For both datasets, we have extracted the distances between residues with a given catalytic function and their target moieties. The geometry of residues whose function involves the transfer or sharing of hydrogens (either with substrate or another residue) was analysed more closely. The results showed that the geometry for such productive interactions (prior to the transition state) closely resembles that seen in non-catalytic hydrogen bonds, with distances and angles in the normal expected range. Such statistics provide limits on "expected geometries" for catalytic residues, which will help to identify these residues and elucidate enzyme mechanisms.     

Identification of active site residues in mevalonate diphosphate decarboxylase: implications for a family of phosphotransferases

A combination of sequence homology analyses of mevalonate diphosphate decarboxylase (MDD) proteins and structural information for MDD leads to the hypothesis that Asp 302 and Lys 18 are active site residues in MDD. These residues were mutated to replace acidic/basic side chains and the mutant proteins were isolated and characterized. Binding and competitive displacement studies using trinitrophenyl-ATP, a fluorescent analog of substrate ATP, indicate that these mutant enzymes (D302A, D302N, K18M) retain the ability to stoichiometrically bind nucleotide triphosphates at the active site. These observations suggest the structural integrity of the mutant MDD proteins. The functional importance of mutated residues was evaluated by kinetic analysis. The 10(3) and 10(5)-fold decreases in k(cat) observed for the Asp 302 mutants (D302N and D302A, respectively) support assignment of a crucial catalytic role to Asp 302. A 30-fold decrease in activity and a 16-fold inflation of the K(m) for ATP is documented for the K18M mutant, indicating that Lys 18 influences the active site but is not crucial for reaction chemistry. Demonstration of the influence of conserved aspartate 302 appears to represent the first documentation of the functional importance of a residue in the MDD catalytic site and affords insight into phosphotransferase reactions catalyzed by a variety of enzymes in the galactokinase, homoserine kinase, mevalonate kinase, phosphom-evalonate kinase (GHMP kinase) family.