Milestones in directed enzyme evolution

Directed evolution has now been used for over two decades as an alternative to rational design for protein engineering. Protein function, however, is complex, and modifying enzyme activity is a tall order. We can now improve existing enzyme activity, change enzyme selectivity and evolve function de novo using directed evolution. Although directed evolution is now used routinely to improve existing enzyme activity, there are still only a handful of examples where substrate selectivity has been modified sufficiently for practical application, and the de novo evolution of function largely eludes us.     

How well is enzyme function conserved as a function of pairwise sequence identity?

Enzyme function conservation has been used to derive the threshold of sequence identity necessary to transfer function from a protein of known function to an unknown protein. Using pairwise sequence comparison, several studies suggested that when the sequence identity is above 40%, enzyme function is well conserved. In contrast, Rost argued that because of database bias, the results from such simple pairwise comparisons might be misleading. Thus, by grouping enzyme sequences into families based on sequence similarity and selecting representative sequences for comparison, he showed that enzyme function starts to diverge quickly when the sequence identity is below 70%. Here, we employ a strategy similar to Rost's to reduce the database bias; however, we classify enzyme families based not only on sequence similarity, but also on functional similarity, i.e. sequences in each family must have the same four digits or the same first three digits of the enzyme commission (EC) number. Furthermore, instead of selecting representative sequences for comparison, we calculate the function conservation of each enzyme family and then average the degree of enzyme function conservation across all enzyme families. Our analysis suggests that for functional transferability, 40% sequence identity can still be used as a confident threshold to transfer the first three digits of an EC number; however, to transfer all four digits of an EC number, above 60% sequence identity is needed to have at least 90% accuracy. Moreover, when PSI-BLAST is used, the magnitude of the E-value is found to be weakly correlated with the extent of enzyme function conservation in the third iteration of PSI-BLAST. As a result, functional annotation based on the E-values from PSI-BLAST should be used with caution. We also show that by employing an enzyme family-specific sequence identity threshold above which 100% functional conservation is required, functional inference of unknown sequences can be accurately accomplished. However, this comes at a cost: those true positive sequences below this threshold cannot be uniquely identified.     

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.     

Optimization of batch processes involving simultaneous enzymatic and microbial reactions

Two general models for batch simultaneous enzymatic and microbial reaction (SEMR) processes are presented, the second derived from and simpler than the first and accounting for enzyme denaturation. Using the second model and parameter values from the literature, simulation was used to examine a range of enzyme addition rate strategies (in which the rate was a linear function of time) for a relatively fast ethanol fermentation and for a longer duration citric acid fermentation, both using cellulose as the substrate. For the ethanol process it is optimal (for a specific objective function which accounts for product value and enzyme cost) to add all the enzyme at the beginning of the process. But for the citric acid process a linearly decreasing enzyme addition rate, coupled with the addition of a small fraction of the enzyme at time zero, is better than pure batch operation or operation with the best constant enzyme feed rate.     

Understanding enzyme structure and function in terms of the shifting specificity model

The purpose of this paper is to suggest that the prominence of Haldane's explanation for enzyme catalysis significantly hinders investigations in understanding enzyme structure and function. This occurs despite the existence of much evidence that the Haldane model cannot embrace. Some of the evidence, in fact, disproves the model. A brief history of the explanation of enzyme catalysis is presented. The currently accepted view of enzyme catalysis--the Haldane model--is examined in terms of its strengths and weaknesses. An alternate model for general enzyme catalysis (the Shifting Specificity model) is reintroduced and an assessment of why it may be superior to the Haldane model is presented. Finally, it is proposed that a re-examination of many current aspects in enzyme structure and function (specifically, protein folding, x-ray and NMR structure analyses, enzyme stability curves, enzyme mimics, catalytic antibodies, and the loose packing of enzyme folded forms) in terms of the new model may offer crucial insights.     

Evolutionarily conserved substrate substructures for automated annotation of enzyme superfamilies

The evolution of enzymes affects how well a species can adapt to new environmental conditions. During enzyme evolution, certain aspects of molecular function are conserved while other aspects can vary. Aspects of function that are more difficult to change or that need to be reused in multiple contexts are often conserved, while those that vary may indicate functions that are more easily changed or that are no longer required. In analogy to the study of conservation patterns in enzyme sequences and structures, we have examined the patterns of conservation and variation in enzyme function by analyzing graph isomorphisms among enzyme substrates of a large number of enzyme superfamilies. This systematic analysis of substrate substructures establishes the conservation patterns that typify individual superfamilies. Specifically, we determined the chemical substructures that are conserved among all known substrates of a superfamily and the substructures that are reacting in these substrates and then examined the relationship between the two. Across the 42 superfamilies that were analyzed, substantial variation was found in how much of the conserved substructure is reacting, suggesting that superfamilies may not be easily grouped into discrete and separable categories. Instead, our results suggest that many superfamilies may need to be treated individually for analyses of evolution, function prediction, and guiding enzyme engineering strategies. Annotating superfamilies with these conserved and reacting substructure patterns provides information that is orthogonal to information provided by studies of conservation in superfamily sequences and structures, thereby improving the precision with which we can predict the functions of enzymes of unknown function and direct studies in enzyme engineering. Because the method is automated, it is suitable for large-scale characterization and comparison of fundamental functional capabilities of both characterized and uncharacterized enzyme superfamilies.     

Kinetics of enzyme action of Cartazyme NS-10 prior to bleaching of kraft pulp

The kinetics of enzyme treatment of hardwood kraft pulp with commercial xylanase Cartazyme NS-10 was investigated. The enzyme treatment process was found to follow closely the topochemical modified equation of Prout–Tompkins. The influence of the initial enzyme concentration was studied and the applicability of the power kinetic equation was established for the initial rate of the process. An equation of practical use was obtained, which provides the temperature function of the amount of the reducing substances, depending on the enzyme initial concentration. This function made it possible to control the process of prior enzyme treatment of unbleached kraft pulp.     

Enzyme function less conserved than anticipated

The level of sequence similarity that implies similarity in protein structure is well established. Recently, many groups proposed thresholds for similarity in sequence implying similarity in enzymatic function. All previous results suggest the strong conservation of enzymatic function above levels of 50% pairwise sequence identity. Here, I argue that all groups substantially overestimated the conservation of enzyme function because their data sets were either too biased, or too small. An unbiased analysis suggested that less than 30% of the pair fragments above 50% sequence identity have entirely identical EC numbers. Another surprising finding was that even BLAST E-values below 10(-50) did not suffice to automatically transfer enzyme function without errors. As expected, most misclassifications originated from similarities in relatively short regions and/or from transferring annotations for different domains. Both problems cannot be corrected easily by adjusting the thresholds for automatic transfer of genome annotations. A score relating sequence identity to alignment length (distance from HSSP-threshold) outperformed statistical BLAST scores for high sequence similarity. In particular, the distance score allowed error-free transfer of enzyme function for the 10% most similar enzyme pairs. The results illustrated how difficult it is to assess the conservation of protein function and to guarantee error-free genome annotations, in general: sets with millions of pair comparisons might not suffice to arrive at statistically significant conclusions. In practice, the revised detailed estimates for the sequence conservation of enzyme function may provide important benchmarks for everyday sequence analysis and for more cautious automatic genome annotations.     

Properties and subcellular distribution of two partially purified ornithine transcarbamoylases in cell suspensions of sugarcane

The spatially separated forms of ornithine transcarbamoylase (EC 2.1.3.3) of different molecular weights coexist in sugarcane (Saccharum sp.). The smaller form of the enzyme (mol wt 79,000) appears to be cytoplasmic, while a larger form (mol wt 224,000) sedimented with mitochondria. The Km of the cytoplasmic enzyme for ornithine was 3.11 mm, while the enzyme in the mitochondrial fraction had a Km of 0.50 mm for this substrate; both enzymes had similar affinity for carbamoyl phosphate (0.12 mm). Characteristics of the smaller ornithine transcarbamoylase are in keeping with a predominantly catabolic function, those of the enzyme which sediments with mitochondria, with an anabolic function. Only the mitochondrial enzyme was regulated in vivo by exogenous arginine.     

Substrate specificity of nicotinamide methyltransferase isolated from porcine liver

Nicotinamide methyltransferase (EC 2.1.1.1) has been purified over 1300-fold from porcine liver. The enzyme is electrophoretically homogeneous, exhibiting a relative molecular mass of 27,000. In addition to acting on nicotinamide and close structural analogs such as thionicotinamide and 3-acetylpyridine, the enzyme actively accommodates poor analogs such as quinoline, isoquinoline, and 1,2,3,4-tetrahydroisoquinoline as methyl group acceptors. The enzyme may thus have the function of detoxicating numerous alkaloids in vivo. In some cases, the action of the enzyme might paradoxically increase the toxicities of substrates, but the hepatotoxic antibiotic pyrazinamide, which we considered as potentially such an enzyme-activated electrophile, did not function detectably as a substrate for the isolated enzyme.     

Influence of pH on enzyme stabilization: an analysis using a series-type mechanism

A series-type model is utilized to show the influence of pH on enzyme inactivation kinetics and stability. Examples of enzyme inactivations involving both single-step and series-type mechanisms are presented. Empirical relations for the inactivation rate constant for the first step and the residual activity as a function of pH are presented. This provides physical insights into the enzyme inactivation processes. The analysis forms the beginning of a framework within which one could quantitatively manipulate the inactivation rate constants and the residual activity for enzymes in desired directions as a function of pH.     

The carbamoyl-phosphate synthetase of Pyrococcus furiosus is enzymologically and structurally a carbamate kinase.

The hyperthermophiles Pyrococcus furiosus and Pyrococcus abyssi make pyrimidines and arginine from carbamoyl phosphate (CP) synthesized by an enzyme that differs from other carbamoyl-phosphate synthetases and that resembles carbamate kinase (CK) in polypeptide mass, amino acid sequence, and oligomeric organization. This enzyme was reported to use ammonia, bicarbonate, and two ATP molecules as carbamoyl-phosphate synthetases to make CP and to exhibit bicarbonatedependent ATPase activity. We have reexamined these findings using the enzyme of P. furiosus expressed in Escherichia coli from the corresponding gene cloned in a plasmid. We show that the enzyme uses chemically made carbamate rather than ammonia and bicarbonate and catalyzes a reaction with the stoichiometry and equilibrium that are typical for CK. Furthermore, the enzyme catalyzes actively full reversion of the CK reaction and exhibits little bicarbonate-dependent ATPase. In addition, it cross-reacts with antibodies raised against CK from Enterococcus faecium, and its three-dimensional structure, judged by x-ray crystallography of enzyme crystals, is very similar to that of CK. Thus, the enzyme is, in all respects other than its function in vivo, a CK. Because in other organisms the function of CK is to make ATP from ADP and CP derived from arginine catabolism, this is the first example of using CK for making rather than using CP. The reasons for this use and the adaptation of the enzyme to this new function 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.     

酶的分子设计、改造与工程应用

酶工程的研究已经发展到分子水平 ,在体外通过基因工程、化学、物理等手段改造酶分子结构与功能 ,大幅提高了酶分子的进化效率和催化效率 ,生产有价值的非天然酶。对酶工程学若干“热点”和前沿课题的研究、应用进行了概述 ,分析了国际上酶工程研究及应用技术、手段、方法 ,包括体外分子进化、核酶和抗体酶的设计、酶分子的定向固定化技术、酶蛋白分子的化学修饰、融合酶、人工合成及模拟酶等技术 ,并展望了酶工程的技术进步和应用的新进展。     

Experimental Evolution of a New Enzymatic Function. II. Evolution of Multiple Functions for EBG Enzyme in E. COLI

The evolution of ebgo enzyme of Escherichia coli, an enzyme which is unable to hydrolyze lactose, lactulose, lactobionate, or galactose-arabinoside effectively, has been directed in successive steps so that the evolved enzyme is able to hydrolyze these galactosides effectively. I show that in order for a strain of E. coli with a lacZ deletion to evolve the ability to use lactobionate as a carbon source, a series of mutations must occur in the ebg genes, and that these mutations must be selected in a particular order. The ordered series of mutations constitutes an obligatory evolutionary pathway for the acquisition of a new function for ebgo enzyme. A comparison of newly evolved strains with parental strains shows that when ebg enzyme acquires a new function, its old functions often suffer; but that in several cases old functions are either unaffected or are improved. I conclude that divergence of functions catalyzed by an enzyme need not require gene duplication.     

Increasing the thermal stability of an oligomeric protein, beta-glucuronidase.

The reporter enzyme beta-glucuronidase was mutagenized and evolved for thermostability. After four cycles of screening the best variant was more active than the wild-type enzyme, and retained function at 70 degrees C, whereas the wild-type enzyme lost function at 65 degrees C. Variants derived from sequential mutagenesis were shuffled together, and re-screened for thermostability. The best variants retained activities at even higher temperatures (80 degrees C), but had specific activities that were now less than that of the wild-type enzyme. The mutations clustered near the tetramer interface of the enzyme, and many of the evolved variants showed much greater resistance to quaternary structure disruption at high temperatures, which is also a characteristic of naturally thermostable enzymes. Together, these results suggest a pathway for the evolution of thermostability in which enzymes initially become stable at high temperatures without loss of activity at low temperatures, while further evolution leads to enzymes that have kinetic parameters that are optimized for high temperatures.     

谷氨酰-脯氨酰-tRNA合成酶的非经典功能——炎症相关基因特异性的翻译沉默

氨基酰-tRNA合成酶（aaRS）催化氨基酰化反应,为生物体内的蛋白质合成提供原料。哺乳动物细胞质中一种双功能aaRS谷氨酰-脯氨酰-tRNA合成酶（EPRS）通常负责将谷氨酸和脯氨酸分别接载到对应的tRNA的3’末端参与蛋白质翻译;此外,它还具有与氨基酰化经典功能无关的单核/巨噬细胞特异性炎症相关基因翻译沉默的非经典功能。在过去的十五年间,对于EPRS参与炎症反应相关基因表达调控的功能研究取得了一系列重要进展,揭示了看家基因EPRS与人类疾病发生和发展的潜在联系。     

N-terminal domains of the human telomerase catalytic subunit required for enzyme activity in vivo

Most tumor cells depend upon activation of the ribonucleoprotein enzyme telomerase for telomere maintenance and continual proliferation. The catalytic activity of this enzyme can be reconstituted in vitro with the RNA (hTR) and catalytic (hTERT) subunits. However, catalytic activity alone is insufficient for the full in vivo function of the enzyme. In addition, the enzyme must localize to the nucleus, recognize chromosome ends, and orchestrate telomere elongation in a highly regulated fashion. To identify domains of hTERT involved in these biological functions, we introduced a panel of 90 N-terminal hTERT substitution mutants into telomerase-negative cells and assayed the resulting cells for catalytic activity and, as a marker of in vivo function, for cellular proliferation. We found four domains to be essential for in vitro and in vivo enzyme activity, two of which were required for hTR binding. These domains map to regions defined by sequence alignments and mutational analysis in yeast, indicating that the N terminus has also been functionally conserved throughout evolution. Additionally, we discovered a novel domain, DAT, that "dissociates activities of telomerase," where mutations left the enzyme catalytically active, but was unable to function in vivo. Since mutations in this domain had no measurable effect on hTERT homomultimerization, hTR binding, or nuclear targeting, we propose that this domain is involved in other aspects of in vivo telomere elongation. The discovery of these domains provides the first step in dissecting the biological functions of human telomerase, with the ultimate goal of targeting this enzyme for the treatment of human cancers.     

Testing for diffusion limitations in salt-activated enzyme catalysts operating in organic solvents

The dramatic activation of serine proteases in nonaqueous media resulting from lyophilization in the presence of KCl is shown to be unrelated to relaxation of potential substrate diffusional limitations. Specifically, lyophilizing subtilisin Carlsberg in the presence of KCl and phosphate buffer in different proportions, ranging from 99% (w/w) enzyme to 1% (w/w) enzyme in the final lyophilized solids, resulted in biocatalyst preparations that were not influenced by substrate diffusion. This result was made evident through use of a classical analysis whereby initial catalytic rates, normalized per weight of total enzyme in the catalyst material, were measured as a function of active enzyme for biocatalyst preparations containing different ratios of active to inactive enzyme. The active enzyme content of a given biocatalyst preparation was controlled by mixing native subtilisin with subtilisin preinactivated with PMSF, a serine protease inhibitor, and lyophilizing the enzyme mixture in the presence of different fractions of KCl and phosphate buffer. Plots of initial reaction rates as a function of percent active subtilisin in the biocatalyst were linear for all biocatalyst preparations. Thus, enzyme activation (reported elsewhere to be as high as 3750-fold in hexane for the transesterification of N-Ac-L-Phe-OEt with n-PrOH) is a manifestation of intrinsic enzyme activation and not relaxation of diffusional limitations resulting from diluted enzyme preparations. Similar activation is reported herein for thermolysin, a nonserine protease, thereby demonstrating that enzyme activation due to lyophilization in the presence of KCl may be a general phenomenon for proteolytic enzymes.     

Categorization of enzyme deactivations using a series-type mechanism

A series-type enzyme deactivation model involving an active enzyme precursor and a final enzyme state with possible non-zero activity is proposed to categorize enzyme deactivation curves. The enzyme activity is a weighted function of the active enzyme states. The deactivation curves may be broadly classified into two major categories wherein the activity is either always less than or it may be more than the initial activity for some time period. Data taken from the literature may be classified into 14 cases. Complex enzyme deactivation curves exhibiting enzyme stabilization and a flex are some of the features that are classified.