Using mechanism similarity to understand enzyme evolution

Enzyme reactions take place in the active site through a series of catalytic steps, which are collectively termed the enzyme mechanism. The catalytic step is thereby the individual unit to consider for the purposes of building new enzyme mechanisms — i.e. through the mix and match of individual catalytic steps, new enzyme mechanisms and reactions can be conceived. In the case of natural evolution, it has been shown that new enzyme functions have emerged through the tweaking of existing mechanisms by the addition, removal, or modification of some catalytic steps, while maintaining other steps of the mechanism intact. Recently, we have extracted and codified the information on the catalytic steps of hundreds of enzymes in a machine-readable way, with the aim of automating this kind of evolutionary analysis. In this paper, we illustrate how these data, which we called the “rules of enzyme catalysis”, can be used to identify similar catalytic steps across enzymes that differ in their overall function and/or structural folds. A discussion on a set of three enzymes that share part of their mechanism is used as an exemplar to illustrate how this approach can reveal divergent and convergent evolution of enzymes at the mechanistic level.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12551-022-01022-9.     

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.     

Using reaction mechanism to measure enzyme similarity

The concept of reaction similarity has been well studied in terms of the overall transformation associated with a reaction, but not in terms of mechanism. We present the first method to give a quantitative measure of the similarity of reactions based upon their explicit mechanisms. Two approaches are presented to measure the similarity between individual steps of mechanisms: a fingerprint-based approach that incorporates relevant information on each mechanistic step; and an approach based only on bond formation, cleavage and changes in order. The overall similarity for two reaction mechanisms is then calculated using the Needleman-Wunsch alignment algorithm. An analysis of MACiE, a database of enzyme mechanisms, using our measure of similarity identifies some examples of convergent evolution of chemical mechanisms. In many cases, mechanism similarity is not reflected by similarity according to the EC system of enzyme classification. In particular, little mechanistic information is conveyed by the class level of the EC system.     

The enzymatic nature of an anonymous protein sequence cannot reliably be inferred from superfamily level structural information alone

As the largest fraction of any proteome does not carry out enzymatic functions, and in order to leverage 3D structural data for the annotation of increasingly higher volumes of sequence data, we wanted to assess the strength of the link between coarse grained structural data (i.e., homologous superfamily level) and the enzymatic versus non-enzymatic nature of protein sequences. To probe this relationship, we took advantage of 41 phylogenetically diverse (encompassing 11 distinct phyla) genomes recently sequenced within the GEBA initiative, for which we integrated structural information, as defined by CATH, with enzyme level information, as defined by Enzyme Commission (EC) numbers. This analysis revealed that only a very small fraction (about 1%) of domain sequences occurring in the analyzed genomes was found to be associated with homologous superfamilies strongly indicative of enzymatic function. Resorting to less stringent criteria to define enzyme versus non-enzyme biased structural classes or excluding highly prevalent folds from the analysis had only modest effect on this proportion. Thus, the low genomic coverage by structurally anchored protein domains strongly associated to catalytic activities indicates that, on its own, the power of coarse grained structural information to infer the general property of being an enzyme is rather limited.     

Hierarchical classification of hydrolases catalytic sites

Universal ontology of catalytic sites is required to systematize enzyme catalytic sites, their evolution as well as relations between catalytic sites and protein families, organisms and chemical reactions. Here we present a classification of hydrolases catalytic sites based on hierarchical organization. The web-accessible database provides information on the catalytic sites, protein folds, EC numbers and source organisms of the enzymes and includes software allowing for analysis and visualization of the relations between them. AVAILABILITY: http://www.enzyme.chem.msu.ru/hcs/     

Quantitative Comparison of Catalytic Mechanisms and Overall Reactions in Convergently Evolved Enzymes: Implications for Classification of Enzyme Function

Functionally analogous enzymes are those that catalyze similar reactions on similar substrates but do not share common ancestry, providing a window on the different structural strategies nature has used to evolve required catalysts. Identification and use of this information to improve reaction classification and computational annotation of enzymes newly discovered in the genome projects would benefit from systematic determination of reaction similarities. Here, we quantified similarity in bond changes for overall reactions and catalytic mechanisms for 95 pairs of functionally analogous enzymes (non-homologous enzymes with identical first three numbers of their EC codes) from the MACiE database. Similarity of overall reactions was computed by comparing the sets of bond changes in the transformations from substrates to products. For similarity of mechanisms, sets of bond changes occurring in each mechanistic step were compared; these similarities were then used to guide global and local alignments of mechanistic steps. Using this metric, only 44% of pairs of functionally analogous enzymes in the dataset had significantly similar overall reactions. For these enzymes, convergence to the same mechanism occurred in 33% of cases, with most pairs having at least one identical mechanistic step. Using our metric, overall reaction similarity serves as an upper bound for mechanistic similarity in functional analogs. For example, the four carbon-oxygen lyases acting on phosphates (EC 4.2.3) show neither significant overall reaction similarity nor significant mechanistic similarity. By contrast, the three carboxylic-ester hydrolases (EC 3.1.1) catalyze overall reactions with identical bond changes and have converged to almost identical mechanisms. The large proportion of enzyme pairs that do not show significant overall reaction similarity (56%) suggests that at least for the functionally analogous enzymes studied here, more stringent criteria could be used to refine definitions of EC sub-subclasses for improved discrimination in their classification of enzyme reactions. The results also indicate that mechanistic convergence of reaction steps is widespread, suggesting that quantitative measurement of mechanistic similarity can inform approaches for functional annotation.     

Evolutionary trace analysis of the alpha-D-phosphohexomutase superfamily

The alpha-D-phosphohexomutase superfamily is composed of four related enzymes that catalyze a reversible, intramolecular phosphoryl transfer on their sugar substrates. The enzymes in this superfamily play important and diverse roles in carbohydrate metabolism in organisms from bacteria to humans. Recent structural and mechanistic studies of one member of this superfamily, phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa, have provided new insights into enzyme mechanism and substrate recognition. Here we use sequence-sequence and sequence-structure comparisons via evolutionary trace analysis to examine 71 members of the alpha-D-phosphohexomutase superfamily. These analyses show that key residues in the active site, including many of those involved in substrate contacts in the P. aeruginosa PMM/PGM complexes, are conserved throughout the enzyme family. Several important regions show class-specific differences in sequence that appear to be correlated with differences in substrate specificity exhibited by subgroups of the family. In addition, we describe the translocation of a 20-residue segment containing the catalytic phosphoserine of phosphoacetylglucosamine mutase, which uniquely identifies members of this subgroup.     

Multisite promiscuity in the processing of endogenous substrates by human carboxylesterase 1

Human carboxylesterase 1 (hCE1) is a drug and endobiotic-processing serine hydrolase that exhibits relatively broad substrate specificity. It has been implicated in a variety of endogenous cholesterol metabolism pathways including the following apparently disparate reactions: cholesterol ester hydrolysis (CEH), fatty acyl Coenzyme A hydrolysis (FACoAH), acyl-Coenzyme A:cholesterol acyltransfer (ACAT), and fatty acyl ethyl ester synthesis (FAEES). The structural basis for the ability of hCE1 to perform these catalytic actions involving large substrates and products has remained unclear. Here we present four crystal structures of the hCE1 glycoprotein in complexes with the following endogenous substrates or substrate analogues: Coenzyme A, the fatty acid palmitate, and the bile acids cholate and taurocholate. While the active site of hCE1 was known to be promiscuous and capable of interacting with a variety of chemically distinct ligands, these structures reveal that the enzyme contains two additional ligand-binding sites and that each site also exhibits relatively non-specific ligand-binding properties. Using this multisite promiscuity, hCE1 appears structurally capable of assembling several catalytic events depending, apparently, on the physiological state of the cellular environment. These results expand our understanding of enzyme promiscuity and indicate that, in the case of hCE1, multiple non-specific sites are employed to perform distinct catalytic actions.     

Molecular Architecture of Photoreceptor Phosphodiesterase Elucidated by Chemical Cross-Linking and Integrative Modeling

Photoreceptor phosphodiesterase (PDE6) is the central effector enzyme in visual excitation pathway in rod and cone photoreceptors. Its tight regulation is essential for the speed, sensitivity, recovery and adaptation of visual detection. Although major steps in the PDE6 activation/deactivation pathway have been identified, mechanistic understanding of PDE6 regulation is limited by the lack of knowledge about the molecular organization of the PDE6 holoenzyme (αβγγ). Here, we characterize the PDE6 holoenzyme by integrative structural determination of the PDE6 catalytic dimer (αβ), based primarily on chemical cross-linking and mass spectrometric analysis. Our models built from high-density cross-linking data elucidate a parallel organization of the two catalytic subunits, with juxtaposed α-helical segments within the tandem regulatory GAF domains to provide multiple sites for dimerization. The two catalytic domains exist in an open configuration when compared to the structure of PDE2 in the apo state. Detailed structural elements for differential binding of the γ-subunit to the GAFa domains of the α- and β-subunits are revealed, providing insight into the regulation of the PDE6 activation/deactivation cycle.     

Plasticity of enzyme active sites

The expectation is that any similarity in reaction chemistry shared by enzyme homologues is mediated by common functional groups conserved through evolution. However, detailed enzyme studies have revealed the flexibility of many active sites, in that different functional groups, unconserved with respect to position in the primary sequence, mediate the same mechanistic role. Nevertheless, the catalytic atoms might be spatially equivalent. More rarely, the active sites have completely different locations in the protein scaffold. This variability could result from: (1) the hopping of functional groups from one position to another to optimize catalysis; (2) the independent specialization of a low-activity primordial enzyme in different phylogenetic lineages; (3) functional convergence after evolutionary divergence; or (4) circular permutation events.     

Sequence and Structural Features of Enzymes and their Active Sites by EC Class

We have analysed a non-redundant set of 294 enzymes for differences in sequence and structural features between the six main Enzyme Commission (EC) classification groups. This systematic study of enzymes, and their active sites in particular, aims to increase understanding of how the structure of an enzyme relates to its functional role. Many features showed significant differences between the EC classes, including active-site polarity, enzyme size and active-site amino acid propensities. Many attributes correlate with each other to form clusters of related features from which we chose representative features for further analysis. Oxidoreductases have more non-polar active sites, which can be attributed to cofactor binding and a preference for Glu over Asp in active sites in comparison to the other classes. Lyases form a significantly higher proportion of oligomers than any other class, whilst the hydrolases form the largest proportion of monomers. These features were then used in a prediction model that classified each enzyme into its top EC class with an accuracy of 33.1%, which is an increase of 16.4% over random classification. Understanding the link between structure and function is critical to improving enzyme design and the prediction of protein function from structure without transfer of annotation from alignments.     

Using a library of structural templates to recognise catalytic sites and explore their evolution in homologous families

Catalytic site structure is normally highly conserved between distantly related enzymes. As a consequence, templates representing catalytic sites have the potential to succeed at function prediction in cases where methods based on sequence or overall structure fail. There are many methods for searching protein structures for matches to structural templates, but few validated template libraries to use with these methods. We present a library of structural templates representing catalytic sites, based on information from the scientific literature. Furthermore, we analyse homologous template families to discover the diversity within families and the utility of templates for active site recognition. Templates representing the catalytic sites of homologous proteins mostly differ by less than 1A root mean square deviation, even when the sequence similarity between the two proteins is low. Within these sets of homologues there is usually no discernible relationship between catalytic site structure similarity and sequence similarity. Because of this structural conservation of catalytic sites, the templates can discriminate between matches to related proteins and random matches with over 85% sensitivity and predictive accuracy. Templates based on protein backbone positions are more discriminating than those based on side-chain atoms. These analyses show encouraging prospects for prediction of functional sites in structural genomics structures of unknown function, and will be of use in analyses of convergent evolution and exploring relationships between active site geometry and chemistry. The template library can be queried via a web server at and is available for download.     

Modeling active RNA structures using the intersection of conformational space: application to the lead-activated ribozyme.

The Pb2+ cleavage of a specific phosphodiester bond in yeast tRNA(Phe) is the classical model of metal-assisted RNA catalysis. In vitro selection experiments have identified a tRNA(Phe) variant, the leadzyme, that is very active in cleavage by Pb2+. We present here a three-dimensional modeling protocol that was used to propose a structure for this ribozyme, and is based on the computation of the intersection of conformational space of sequence variants and the use of chemical modification data. Sequence and secondary structure data were used in a first round of computer modeling that allowed identification of conformations compatible with all known leadzyme variants. Common conformations were then tested experimentally by evaluating the activity of analogues containing modified nucleotides in the catalytic core. These experiments led to a new structural hypothesis that was tested in a second round of computer modeling. The resulting proposal for the active conformation of the leadzyme is consistent with all known structural data. The final model suggests an in-line SN2 attack mechanism and predicts two Pb2+ binding sites. The protocol presented here is generally applicable in modeling RNAs whenever the catalytic or binding activity of structural analogues is known.     

Exploring Protein Space: From Hydrolase to Ligase by Substitution

The understanding of how proteins evolve to perform novel functions has long been sought by biologists. In this regard, two homologous bacterial enzymes, PafA and Dop, pose an insightful case study, as both rely on similar mechanistic properties, yet catalyze different reactions. PafA conjugates a small protein tag to target proteins, whereas Dop removes the tag by hydrolysis. Given that both enzymes present a similar fold and high sequence similarity, we sought to identify the differences in the amino acid sequence and folding responsible for each distinct activity. We tackled this question using analysis of sequence–function relationships, and identified a set of uniquely conserved residues in each enzyme. Reciprocal mutagenesis of the hydrolase, Dop, completely abolished the native activity, at the same time yielding a catalytically active ligase. Based on the available Dop and PafA crystal structures, this change of activity required a conformational change of a critical loop at the vicinity of the active site. We identified the conserved positions essential for stabilization of the alternative loop conformation, and tracked alternative mutational pathways that lead to a change in activity. Remarkably, all these pathways were combined in the evolution of PafA and Dop, despite their redundant effect on activity. Overall, we identified the residues and structural elements in PafA and Dop responsible for their activity differences. This analysis delineated, in molecular terms, the changes required for the emergence of a new catalytic function from a preexisting one.     

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.     

Interaction of aminoacyl-tRNA-synthetases with amino acids

The review deals with interactions of the key enzymes of the protein biosynthesis-aminoacyl-tRNA synthetases (EC 6.1.1.) with amino acids and their analogues, considering the contribution of different groups in the process of specific complex formation and catalysis. The important role of alpha-amino group of amino acid in the enzyme recognition has been revealed. Modification of the carboxylic group does not change significantly the analogues complex formation with aminoacyl-tRNA synthetases. However this group is essential for amino acid rearrangement in the specific complex with the enzyme. The structural organization of the enzyme binding sites specific for amino acids and the enzyme interaction with the analogues of aminoacyladenylates are discussed.     

Genome-wide analysis of enzyme structure-function combination across three domains of life

To investigate diverse enzyme structure-function combination (SFC) types in different species, 34 different genome sequences were annotated using the protein catalytic domain database SCOPEC (http://www.enzome.com/enzome/), in which both the structure and function for each entry are known. Annotated enzymes with catalytic domains from the same SCOP superfamily are considered to have an identical structure. Annotated enzymes sharing the identical three-digit EC number are considered to have the same enzymatic function. Results reveal that the different SFC types for enzymes identified in archaea, bacteria and eukaryota are 137, 300 and 313, respectively. About 80% of the SFCs identified in archaea can be consistently found in bacteria and eukaryota species, whereas 28% and 35% combination types in bacteria and eukaryota respectively are unique to their corresponding groups. The number of functions per structure and the number of structures per function for the annotated sequences were measured in different species. Furthermore, a new concept was proposed to represent enzymatic structures as a functional similarity network. Thus, the current study will be helpful to enhance the global view on the evolution of enzymatic structure and function.     

Substrate Shuttling between Active Sites of Uroporphyrinogen Decarboxylase Is Not Required to Generate Coproporphyrinogen

Uroporphyrinogen decarboxylase (URO-D; EC 4.1.1.37), the fifth enzyme of the heme biosynthetic pathway, is required for the production of heme, vitamin B12, siroheme, and chlorophyll precursors. URO-D catalyzes the sequential decarboxylation of four acetate side chains in the pyrrole groups of uroporphyrinogen to produce coproporphyrinogen. URO-D is a stable homodimer, with the active-site clefts of the two subunits adjacent to each other. It has been hypothesized that the two catalytic centers interact functionally, perhaps by shuttling of reaction intermediates between subunits. We tested this hypothesis by construction of a single-chain protein (single-chain URO-D) in which the two subunits were connected by a flexible linker. The crystal structure of this protein was shown to be superimposable with wild-type activity and to have comparable catalytic activity. Mutations that impaired one or the other of the two active sites of single-chain URO-D resulted in approximately half of wild-type activity. The distributions of reaction intermediates were the same for mutant and wild-type sequences and were unaltered in a competition experiment using I and III isomer substrates. These observations indicate that communication between active sites is not required for enzyme function and suggest that the dimeric structure of URO-D is required to achieve conformational stability and to create a large active-site cleft.     

Catalytic Turnover-Based Phage Selection for Engineering the Substrate Specificity of Sfp Phosphopantetheinyl Transferase

We report a high-throughput phage selection method to identify mutants of Sfp phosphopantetheinyl transferase with altered substrate specificities from a large library of the Sfp enzyme. In this method, Sfp and its peptide substrates are co-displayed on the M13 phage surface as fusions to the phage capsid protein pIII. Phage-displayed Sfp mutants that are active with biotin-conjugated coenzyme A (CoA) analogues would covalently transfer biotin to the peptide substrates anchored on the same phage particle. Affinity selection for biotin-labeled phages would enrich Sfp mutants that recognize CoA analogues for carrier protein modification. We used this method to successfully change the substrate specificity of Sfp and identified mutant enzymes with more than 300-fold increase in catalytic efficiency with 3′-dephospho CoA as the substrate. The method we developed in this study provides a useful platform to display enzymes and their peptide substrates on the phage surface and directly couples phage selection with enzyme catalysis. We envision this method to be applied to engineering the catalytic activities of other protein posttranslational modification enzymes.     

Glycerol dehydrogenase. structure, specificity, and mechanism of a family III polyol dehydrogenase

BACKGROUND: Bacillus stearothermophilus glycerol dehydrogenase (GlyDH) (glycerol:NAD(+) 2-oxidoreductase, EC 1.1.1.6) catalyzes the oxidation of glycerol to dihydroxyacetone (1,3-dihydroxypropanone) with concomitant reduction of NAD(+) to NADH. Analysis of the sequence of this enzyme indicates that it is a member of the so-called iron-containing alcohol dehydrogenase family. Despite this sequence similarity, GlyDH shows a strict dependence on zinc for activity. On the basis of this, we propose to rename this group the family III metal-dependent polyol dehydrogenases. To date, no structural data have been reported for any enzyme in this group. RESULTS: The crystal structure of B. stearothermophilus glycerol dehydrogenase has been determined at 1.7 A resolution to provide structural insights into the mechanistic features of this family. The enzyme has 370 amino acid residues, has a molecular mass of 39.5 kDa, and is a homooctamer in solution. CONCLUSIONS: Analysis of the crystal structures of the free enzyme and of the binary complexes with NAD(+) and glycerol show that the active site of GlyDH lies in the cleft between the enzyme's two domains, with the catalytic zinc ion playing a role in stabilizing an alkoxide intermediate. In addition, the specificity of this enzyme for a range of diols can be understood, as both hydroxyls of the glycerol form ligands to the enzyme-bound Zn(2+) ion at the active site. The structure further reveals a previously unsuspected similarity to dehydroquinate synthase, an enzyme whose more complex chemistry shares a common chemical step with that catalyzed by glycerol dehydrogenase, providing a striking example of divergent evolution. Finally, the structure suggests that the NAD(+) binding domain of GlyDH may be related to that of the classical Rossmann fold by switching the sequence order of the two mononucleotide binding folds that make up this domain.