Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

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.     

Uniform binding and negative catalysis at the origin of enzymes

Enzymes are well known for their catalytic abilities, some even reaching “catalytic perfection” in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate‐handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a “substrate‐handle‐binding‐first” mechanism (“binding‐first” for brevity) instead of “chemistry‐first” and we are, therefore, left to wonder what the role of non‐catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back‐to‐back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as “uniform binding” (parallel stabilization of all enzyme‐bound states to the same degree). Indeed, we show that a uniform‐binding proto‐catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate‐handle‐binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction''s course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of “binding‐first” enzyme families like the aldolase superfamily.     

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome‐scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence‐level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.     

Substrate-induced Changes in Protease Active Site Conformation Impact on Subsequent Reactions with Substrates

Enzymatic catalysis of biochemical reactions is essential to all living systems. The “lock and key” and “induced fit” models were early contributions to our understanding of the mechanisms involved in the reaction between an enzyme and its substrate. However, whether a given substrate-induced conformation is rigid or remains flexible has not yet been determined. By measuring the enzyme activity and intrinsic fluorescence of a nonspecific Eisenia fetida protease-I with different chromogenic substrates, we show that in subsequent reactions of protease with substrates, both the “lock and key” and “induced fit” mechanisms are used depending on the degree of conformational change required. Chromozym-Th- or chromosym-Ch-induced protease conformations were unable to bind chromozym-U. The chromosym-U-induced protease conformation remained flexible and could be further induced by chromozym-Th and chromozym-Ch. When low concentrations of guanidine HCl were used to disturb the conformation of the enzyme, only small changes in intrinsic fluorescence of the chromozym-Th-induced protease were detected, in contrast to the native enzyme whose intrinsic fluorescence markedly increased. This indicates that the substrate-induced enzyme was relatively rigid compared with the native protease. Utilizing a lock and key mechanism for secondary substrate reactions may have adaptive value in that it facilitates high efficiency in enzymatic reactions.     

A Hybrid Mechanism for the Synechocystis Arsenate Reductase Revealed by Structural Snapshots during Arsenate Reduction

Evolution of enzymes plays a crucial role in obtaining new biological functions for all life forms. Arsenate reductases (ArsC) are several families of arsenic detoxification enzymes that reduce arsenate to arsenite, which can subsequently be extruded from cells by specific transporters. Among these, the Synechocystis ArsC (SynArsC) is structurally homologous to the well characterized thioredoxin (Trx)-coupled ArsC family but requires the glutaredoxin (Grx) system for its reactivation, therefore classified as a unique Trx/Grx-hybrid family. The detailed catalytic mechanism of SynArsC is unclear and how the “hybrid” mechanism evolved remains enigmatic. Herein, we report the molecular mechanism of SynArsC by biochemical and structural studies. Our work demonstrates that arsenate reduction is carried out via an intramolecular thiol-disulfide cascade similar to the Trx-coupled family, whereas the enzyme reactivation step is diverted to the coupling of the glutathione-Grx pathway due to the local structural difference. The current results support the hypothesis that SynArsC is likely a molecular fossil representing an intermediate stage during the evolution of the Trx-coupled ArsC family from the low molecular weight protein phosphotyrosine phosphatase (LMW-PTPase) family.     

A Fundamental Trade-off in Covalent Switching and Its Circumvention by Enzyme Bifunctionality in Glucose Homeostasis

Covalent modification provides a mechanism for modulating molecular state and regulating physiology. A cycle of competing enzymes that add and remove a single modification can act as a molecular switch between “on” and “off” and has been widely studied as a core motif in systems biology. Here, we exploit the recently developed “linear framework” for time scale separation to determine the general principles of such switches. These methods are not limited to Michaelis-Menten assumptions, and our conclusions hold for enzymes whose mechanisms may be arbitrarily complicated. We show that switching efficiency improves with increasing irreversibility of the enzymes and that the on/off transition occurs when the ratio of enzyme levels reaches a value that depends only on the rate constants. Fluctuations in enzyme levels, which habitually occur due to cellular heterogeneity, can cause flipping back and forth between on and off, leading to incoherent mosaic behavior in tissues, that worsens as switching becomes sharper. This trade-off can be circumvented if enzyme levels are correlated. In particular, if the competing catalytic domains are on the same protein but do not influence each other, the resulting bifunctional enzyme can switch sharply while remaining coherent. In the mammalian liver, the switch between glycolysis and gluconeogenesis is regulated by the bifunctional 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2). We suggest that bifunctionality of PFK-2/FBPase-2 complements the metabolic zonation of the liver by ensuring coherent switching in response to insulin and glucagon.     

Enzyme promiscuity – A light on the “darker” side of enzyme specificity

Abstract

Enzyme promiscuity can be defined as the capability of enzymes to catalyse side reaction in addition to its main reaction. The side reaction of an enzyme is termed as promiscuous or sometimes as the “darker” side of enzyme cross-reactivity/specificity. This unique property of enzyme allows organisms to adapt under varying environmental conditions. Promiscuous enzymes can modify their catalytic activities with altered substrates and can adjust their catalytic and kinetic mechanisms according to substrate properties. This group of enzymes evolved from ancestral proteins found in primitive organisms like archaea that survive under extreme environmental conditions. Such ancestral proteins possessed the potential to catalyse a wide range of reactions at low levels, hence create families or superfamilies of highly specialized enzymes. Further, some enzymes were identified which have non-catalytic functions in addition to their major catalytic activities. These enzymes are referred to as moonlighting enzymes. The study of these enzymes will provide important information regarding enzyme evolution and will help in optimizing protein engineering applications.     

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.     

Enzymatic Detoxication,Conformational Selection,and the Role of Molten Globule Active Sites

The role of conformational ensembles in enzymatic reactions remains unclear. Discussion concerning “induced fit” versus “conformational selection” has, however, ignored detoxication enzymes, which exhibit catalytic promiscuity. These enzymes dominate drug metabolism and determine drug-drug interactions. The detoxication enzyme glutathione transferase A1–1 (GSTA1–1), exploits a molten globule-like active site to achieve remarkable catalytic promiscuity wherein the substrate-free conformational ensemble is broad with barrierless transitions between states. A quantitative index of catalytic promiscuity is used to compare engineered variants of GSTA1–1 and the catalytic promiscuity correlates strongly with characteristics of the thermodynamic partition function, for the substrate-free enzymes. Access to chemically disparate transition states is encoded by the substrate-free conformational ensemble. Pre-steady state catalytic data confirm an extension of the conformational selection model, wherein different substrates select different starting conformations. The kinetic liability of the conformational breadth is minimized by a smooth landscape. We propose that “local” molten globule behavior optimizes detoxication enzymes.     

Structural Bases for Stability-Function Tradeoffs in Antibiotic Resistance

Preorganization of enzyme active sites for substrate recognition typically comes at a cost to the stability of the folded form of the protein; consequently, enzymes can be dramatically stabilized by substitutions that attenuate the size and preorganization “strain” of the active site. How this stability-activity tradeoff constrains enzyme evolution has remained less certain, and it is unclear whether one should expect major stability insults as enzymes mutate towards new activities or how these new activities manifest structurally. These questions are both germane and easy to study in β-lactamases, which are evolving on the timescale of years to confer resistance to an ever-broader spectrum of β-lactam antibiotics. To explore whether stability is a substantial constraint on this antibiotic resistance evolution, we investigated extended-spectrum mutants of class C β-lactamases, which had evolved new activity versus third-generation cephalosporins. Five mutant enzymes had between 100-fold and 200-fold increased activity against the antibiotic cefotaxime in enzyme assays, and the mutant enzymes all lost thermodynamic stability (from 1.7 kcal mol^− 1 to 4.1 kcal mol^− 1), consistent with the stability-function hypothesis. Intriguingly, several of the substitutions were 10-20 Å from the catalytic serine; the question of how they conferred extended-spectrum activity arose. Eight structures, including complexes with inhibitors and extended-spectrum antibiotics, were determined by X-ray crystallography. Distinct mechanisms of action, including changes in the flexibility and ground-state structures of the enzyme, are revealed for each mutant. These results explain the structural bases for the antibiotic resistance conferred by these substitutions and their corresponding decrease in protein stability, which will constrain the evolution of new antibiotic resistance.     

Effects of magnesium and related divalent metal ions in topoisomerase structure and function

The catalytic steps through which DNA topoisomerases produce their biological effects and the interference of drug molecules with the enzyme–DNA cleavage complex have been thoroughly investigated both from the biophysical and the biochemical point of view. This provides the basic structural insight on how this family of essential enzymes works in living systems and how their functions can be impaired by natural and synthetic compounds. Besides other factors, the physiological environment is known to affect substantially the biological properties of topoisomerases, a key role being played by metal ion cofactors, especially divalent ions (Mg²⁺), that are crucial to bestow and modulate catalytic activity by exploiting distinctive chemical features such as ionic size, hardness and characteristics of the coordination sphere including coordination number and geometry. Indeed, metal ions mediate fundamental aspects of the topoisomerase-driven transphosphorylation process by affecting the kinetics of the forward and the reverse steps and by modifying the enzyme conformation and flexibility. Of particular interest in type IA and type II enzymes are ionic interactions involving the Toprim fold, a protein domain conserved through evolution that contains a number of acidic residues essential for catalysis. A general two-metal ion mechanism is widely accepted to account for the biophysical and biochemical data thus far available.     

Rapid Identification of Sequences for Orphan Enzymes to Power Accurate Protein Annotation

The power of genome sequencing depends on the ability to understand what those genes and their proteins products actually do. The automated methods used to assign functions to putative proteins in newly sequenced organisms are limited by the size of our library of proteins with both known function and sequence. Unfortunately this library grows slowly, lagging well behind the rapid increase in novel protein sequences produced by modern genome sequencing methods. One potential source for rapidly expanding this functional library is the “back catalog” of enzymology – “orphan enzymes,” those enzymes that have been characterized and yet lack any associated sequence. There are hundreds of orphan enzymes in the Enzyme Commission (EC) database alone. In this study, we demonstrate how this orphan enzyme “back catalog” is a fertile source for rapidly advancing the state of protein annotation. Starting from three orphan enzyme samples, we applied mass-spectrometry based analysis and computational methods (including sequence similarity networks, sequence and structural alignments, and operon context analysis) to rapidly identify the specific sequence for each orphan while avoiding the most time- and labor-intensive aspects of typical sequence identifications. We then used these three new sequences to more accurately predict the catalytic function of 385 previously uncharacterized or misannotated proteins. We expect that this kind of rapid sequence identification could be efficiently applied on a larger scale to make enzymology’s “back catalog” another powerful tool to drive accurate genome annotation.     

New algorithms and an in silico benchmark for computational enzyme design

The creation of novel enzymes capable of catalyzing any desired chemical reaction is a grand challenge for computational protein design. Here we describe two new algorithms for enzyme design that employ hashing techniques to allow searching through large numbers of protein scaffolds for optimal catalytic site placement. We also describe an in silico benchmark, based on the recapitulation of the active sites of native enzymes, that allows rapid evaluation and testing of enzyme design methodologies. In the benchmark test, which consists of designing sites for each of 10 different chemical reactions in backbone scaffolds derived from 10 enzymes catalyzing the reactions, the new methods succeed in identifying the native site in the native scaffold and ranking it within the top five designs for six of the 10 reactions. The new methods can be directly applied to the design of new enzymes, and the benchmark provides a powerful in silico test for guiding improvements in computational enzyme design.     

The Natural History of Biocatalytic Mechanisms

Phylogenomic analysis of the occurrence and abundance of protein domains in proteomes has recently showed that the α/β architecture is probably the oldest fold design. This holds important implications for the origins of biochemistry. Here we explore structure-function relationships addressing the use of chemical mechanisms by ancestral enzymes. We test the hypothesis that the oldest folds used the most mechanisms. We start by tracing biocatalytic mechanisms operating in metabolic enzymes along a phylogenetic timeline of the first appearance of homologous superfamilies of protein domain structures from CATH. A total of 335 enzyme reactions were retrieved from MACiE and were mapped over fold age. We define a mechanistic step type as one of the 51 mechanistic annotations given in MACiE, and each step of each of the 335 mechanisms was described using one or more of these annotations. We find that the first two folds, the P-loop containing nucleotide triphosphate hydrolase and the NAD(P)-binding Rossmann-like homologous superfamilies, were α/β architectures responsible for introducing 35% (18/51) of the known mechanistic step types. We find that these two oldest structures in the phylogenomic analysis of protein domains introduced many mechanistic step types that were later combinatorially spread in catalytic history. The most common mechanistic step types included fundamental building blocks of enzyme chemistry: “Proton transfer,” “Bimolecular nucleophilic addition,” “Bimolecular nucleophilic substitution,” and “Unimolecular elimination by the conjugate base.” They were associated with the most ancestral fold structure typical of P-loop containing nucleotide triphosphate hydrolases. Over half of the mechanistic step types were introduced in the evolutionary timeline before the appearance of structures specific to diversified organisms, during a period of architectural diversification. The other half unfolded gradually after organismal diversification and during a period that spanned ∼2 billion years of evolutionary history.     

Resource Uptake and the Evolution of Moderately Efficient Enzymes

Enzymes speed up reactions that would otherwise be too slow to sustain the metabolism of selfreplicators. Yet, most enzymes seem only moderately efficient, exhibiting kinetic parameters orders of magnitude lower than their expected physically achievable maxima and spanning over surprisingly large ranges of values. Here, we question how these parameters evolve using a mechanistic model where enzyme efficiency is a key component of individual competition for resources. We show that kinetic parameters are under strong directional selection only up to a point, above which enzymes appear to evolve under near-neutrality, thereby confirming the qualitative observation of other modeling approaches. While the existence of a large fitness plateau could potentially explain the extensive variation in enzyme features reported, we show using a population genetics model that such a widespread distribution is an unlikely outcome of evolution on a common landscape, as mutation–selection–drift balance occupy a narrow area even when very moderate biases towards lower efficiency are considered. Instead, differences in the evolutionary context encountered by each enzyme should be involved, such that each evolves on an individual, unique landscape. Our results point to drift and effective population size playing an important role, along with the kinetics of nutrient transporters, the tolerance to high concentrations of intermediate metabolites, and the reversibility of reactions. Enzyme concentration also shapes selection on kinetic parameters, but we show that the joint evolution of concentration and efficiency does not yield extensive variance in evolutionary outcomes when documented costs to protein expression are applied.     

Origins of specificity and promiscuity in metabolic networks

How enzymes have evolved to their present form is linked to the question of how pathways emerged and evolved into extant metabolic networks. To investigate this mechanism, we have explored the chemical diversity present in a largely unbiased data set of catalytic reactions processed by modern enzymes across the tree of life. In order to get a quantitative estimate of enzyme chemical diversity, we measure enzyme multispecificity or promiscuity using the reaction molecular signatures. Our main finding is that reactions that are catalyzed by a highly specific enzyme are shared by poorly divergent species, suggesting a later emergence of this function during evolution. In contrast, reactions that are catalyzed by highly promiscuous enzymes are more likely to appear uniformly distributed across species in the tree of life. From a functional point of view, promiscuous enzymes are mainly involved in amino acid and lipid metabolisms, which might be associated with the earliest form of biochemical reactions. In this way, results presented in this paper might assist us with the identification of primeval promiscuous catalytic functions contributing to life's minimal metabolism.     

MACiE: a database of enzyme reaction mechanisms

SUMMARY: MACiE (mechanism, annotation and classification in enzymes) is a publicly available web-based database, held in CMLReact (an XML application), that aims to help our understanding of the evolution of enzyme catalytic mechanisms and also to create a classification system which reflects the actual chemical mechanism (catalytic steps) of an enzyme reaction, not only the overall reaction. AVAILABILITY: http://www-mitchell.ch.cam.ac.uk/macie/.     

Convergent evolution of enzyme active sites is not a rare phenomenon

Since convergent evolution of enzyme active sites was first identified in serine proteases, other individual instances of this phenomenon have been documented. However, a systematic analysis assessing the frequency of this phenomenon across enzyme space is still lacking. This work uses the Query3d structural comparison algorithm to integrate for the first time detailed knowledge about catalytic residues, available through the Catalytic Site Atlas (CSA), with the evolutionary information provided by the Structural Classification of Proteins (SCOP) database. This study considers two modes of convergent evolution: (i) mechanistic analogues which are enzymes that use the same mechanism to perform related, but possibly different, reactions (considered here as sharing the first three digits of the EC number); and (ii) transformational analogues which catalyse exactly the same reaction (identical EC numbers), but may use different mechanisms. Mechanistic analogues were identified in 15% (26 out of 169) of the three-digit EC groups considered, showing that this phenomenon is not rare. Furthermore 11 of these groups also contain transformational analogues. The catalytic triad is the most widespread active site; the results of the structural comparison show that this mechanism, or variations thereof, is present in 23 superfamilies. Transformational analogues were identified for 45 of the 951 four-digit EC numbers present within the CSA and about half of these were also mechanistic analogues exhibiting convergence of their active sites. This analysis has also been extended to the whole Protein Data Bank to provide a complete and manually curated list of the all the transformational analogues whose structure is classified in SCOP. The results of this work show that the phenomenon of convergent evolution is not rare, especially when considering large enzymatic families.