The evolution of catalytic efficiency and substrate promiscuity in human theta class 1-1 glutathione transferase

Theta class glutathione transferases (GST) from various species exhibit markedly different catalytic activities in conjugating the tripeptide glutathione (GSH) to a variety of electrophilic substrates. For example, the human theta 1-1 enzyme (hGSTT1-1) is 440-fold less efficient than the rat theta 2-2 enzyme (rGSTT2-2) with the fluorogenic substrate 7-amino-4-chloromethyl coumarin (CMAC). Large libraries of hGSTT1-1 constructed by error-prone PCR, DNA shuffling, or saturation mutagenesis were screened for improved catalytic activity towards CMAC in a quantitative fashion using flow cytometry. An iterative directed evolution approach employing random mutagenesis in conjunction with homologous recombination gave rise to enzymes exhibiting up to a 20,000-fold increase in k(cat)/K(M) compared to hGSTT1-1. All highly active clones encoded one or more mutations at residues 32, 176, or 234. Combinatorial saturation mutagenesis was used to evaluate the full complement of natural amino acids at these positions, and resulted in the isolation of enzymes with catalytic rates comparable to those exhibited by the fastest mutants obtained via directed evolution. The substrate selectivities of enzymes resulting from random mutagenesis, DNA shuffling, and combinatorial saturation mutagenesis were evaluated using a series of distinct electrophiles. The results revealed that promiscuous substrate activities arose in a stochastic manner, as they did not correlate with catalytic efficiency towards the CMAC selection substrate. In contrast, chimeric enzymes previously constructed by homology-independent recombination of hGSTT-1 and rGSTT2-2 exhibited very different substrate promiscuity profiles, and showed a more defined relationship between evolved and promiscuous activities.     

Using experimental evolution to probe molecular mechanisms of protein function

Directed evolution is a powerful tool for engineering protein function. The process of directed evolution involves iterative rounds of sequence diversification followed by assaying activity of variants and selection. The range of sequence variants and linked activities generated in the course of an evolution are a rich information source for investigating relationships between sequence and function. Key residue positions determining protein function, combinatorial contributors to activity and even potential functional mechanisms have been revealed in directed evolutions. The recent application of high throughput sequencing substantially increases the information that can be retrieved from directed evolution experiments. Combined with computational analysis this additional sequence information has allowed high‐resolution analysis of individual residue contributions to activity. These developments promise to significantly enhance the depth of insight that experimental evolution provides into mechanisms of protein function.     

Systems-Wide Prediction of Enzyme Promiscuity Reveals a New Underground Alternative Route for Pyridoxal 5’-Phosphate Production in E. coli

Recent insights suggest that non-specific and/or promiscuous enzymes are common and active across life. Understanding the role of such enzymes is an important open question in biology. Here we develop a genome-wide method, PROPER, that uses a permissive PSI-BLAST approach to predict promiscuous activities of metabolic genes. Enzyme promiscuity is typically studied experimentally using multicopy suppression, in which over-expression of a promiscuous ‘replacer’ gene rescues lethality caused by inactivation of a ‘target’ gene. We use PROPER to predict multicopy suppression in Escherichia coli, achieving highly significant overlap with published cases (hypergeometric p = 4.4e-13). We then validate three novel predicted target-replacer gene pairs in new multicopy suppression experiments. We next go beyond PROPER and develop a network-based approach, GEM-PROPER, that integrates PROPER with genome-scale metabolic modeling to predict promiscuous replacements via alternative metabolic pathways. GEM-PROPER predicts a new indirect replacer (thiG) for an essential enzyme (pdxB) in production of pyridoxal 5’-phosphate (the active form of Vitamin B₆), which we validate experimentally via multicopy suppression. We perform a structural analysis of thiG to determine its potential promiscuous active site, which we validate experimentally by inactivating the pertaining residues and showing a loss of replacer activity. Thus, this study is a successful example where a computational investigation leads to a network-based identification of an indirect promiscuous replacement of a key metabolic enzyme, which would have been extremely difficult to identify directly.     

Catalytic and binding poly-reactivities shared by two unrelated proteins: The potential role of promiscuity in enzyme evolution

It is generally accepted that enzymes evolved via gene duplication of existing proteins. But duplicated genes can serve as a starting point for the evolution of a new function only if the protein they encode happens to exhibit some activity towards this new function. Although the importance of such catalytic promiscuity in enzyme evolution has been proposed, little is actually known regarding how common promiscuous catalytic activities are in proteins or their origins, magnitudes, and potential contribution to the survival of an organism. Here we describe a pattern of promiscuous activities in two completely unrelated proteins-serum albumins and a catalytic antibody (aldolase antibody 38C2). Despite considerable structural dissimilarities-in the shape of the cavities and the position of catalytic lysine residues-both active sites are able to catalyze the Kemp elimination, a model reaction for proton transfer from carbon. We also show that these different active sites can bind promiscuously an array of hydrophobic negatively charged ligands. We suggest that the basic active-site features of an apolar pocket and a lysine residue can act as a primitive active site allowing these promiscuous activities to take place. We also describe, by modelling product formation at different substrate concentrations, how promiscuous activities of this kind- inefficient and rudimentary as they are-can provide a considerable selective advantage and a starting point for the evolution of new functions.     

Thematic minireview series on enzyme evolution in the post-genomic era

The emergence of genomics; ongoing computational advances; and the development of large-scale sequence, structural, and functional databases have created important new interdisciplinary linkages between molecular evolution, molecular biology, and enzymology. The five minireviews in this series survey advances and challenges in this burgeoning field from complementary perspectives. The series has three major themes. The first is the evolution of enzyme superfamilies, in which members exhibit increasing sequence, structural, and functional divergence with increasing time of divergence from a common ancestor. The second is the evolutionary role of promiscuous enzymes, which, in addition to their primary function, have adventitious secondary activities that frequently provide the starting point for the evolution of new enzymes. The third is the importance of in silico approaches to the daunting challenge of assigning and predicting the functions of the many uncharacterized proteins in the large-scale sequence and structural databases that are now available. A recent computational advance, the use of protein similarity networks that map functional data onto proteins clustered by similarity, is presented as an approach that can improve functional insight and inference. The three themes are illustrated with several examples of enzyme superfamilies, including the amidohydrolase, metallo-β-lactamase, and enolase superfamilies.     

Employing directed evolution for the functional analysis of multi-specific proteins

Multi-specific proteins located at the heart of complex protein–protein interaction (PPI) networks play essential roles in the survival and fitness of the cell. In addition, multi-specific or promiscuous enzymes exhibit activity toward a wide range of substrates so as to increase cell evolvability and robustness. However, despite their high importance, investigating the in vivo function of these proteins is difficult, due to their complex nature. Typically, deletion of these proteins leads to the abolishment of large PPI networks, highlighting the difficulty in examining the contributions of specific interactions/activities to complex biological processes and cell phenotypes. Protein engineering approaches, including directed evolution and computational protein design, allow for the generation of multi-specific proteins in which certain activities remain intact while others are abolished. The generation and examination of these mutants both in vitro and in vivo can provide high-resolution analysis of biological processes and cell phenotypes and provide new insight into the evolution and molecular function of this important protein family.     

Switch of substrate specificity of hyperthermophilic acylaminoacyl peptidase by combination of protein and solvent engineering

The inherent evolvability of promiscuous enzymes endows them with great potential to be artificially evolved for novel functions. Previously, we succeeded in transforming a promiscuous acylaminoacyl peptidase (apAAP) from the hyperthermophilic archaeon Aeropyrum pernix K1 into a specific carboxylesterase by making a single mutation. In order to fulfill the urgent requirement of thermostable lipolytic enzymes, in this paper we describe how the substrate preference of apAAP can be further changed from p-nitrophenyl caprylate (pNP-C8) to p-nitrophenyl laurate (pNP-C12) by protein and solvent engineering. After one round of directed evolution and subsequent saturation mutagenesis at selected residues in the active site, three variants with enhanced activity towards pNP-C12 were identified. Additionally, a combined mutant W474V/F488G/R526V/T560W was generated, which had the highest catalytic efficiency (k_cat/K_m) for pNP-C12, about 71-fold higher than the wild type. Its activity was further increased by solvent engineering, resulting in an activity enhancement of 280-fold compared with the wild type in the presence of 30% DMSO. The structural basis for the improved activity was studied by substrate docking and molecular dynamics simulation. It was revealed that W474V and F488G mutations caused a significant change in the geometry of the active center, which may facilitate binding and subsequent hydrolysis of bulky substrates. In conclusion, the combination of protein and solvent engineering may be an effective approach to improve the activities of promiscuous enzymes and could be used to create naturally rare hyperthermophilic enzymes.     

Setting the stage for evolution of a new enzyme

The evolution of novel enzymes has fueled the diversification of life on earth for billions of years. Insights into events that set the stage for the evolution of a new enzyme can be obtained from ancestral reconstruction and laboratory evolution. Ancestral reconstruction can reveal the emergence of a promiscuous activity in a pre-existing protein and the impact of subsequent mutations that enhance a new activity. Laboratory evolution provides a more holistic view by revealing mutations elsewhere in the genome that indirectly enhance the level of a newly important enzymatic activity. This review will highlight recent studies that probe the early stages of the evolution of a new enzyme from these complementary points of view.     

An automated flow for directed evolution based on detection of promiscuous scaffolds using spatial and electrostatic properties of catalytic residues

The aspiration to mimic and accelerate natural evolution has fueled interest in directed evolution experiments, which endow or enhance functionality in enzymes. Barring a few de novo approaches, most methods take a template protein having the desired activity, known active site residues and structure, and proceed to select a target protein which has a pre-existing scaffold congruent to the template motif. Previously, we have established a computational method (CLASP) based on spatial and electrostatic properties to detect active sites, and a method to quantify promiscuity in proteins. We exploit the prospect of promiscuous active sites to serve as the starting point for directed evolution and present a method to select a target protein which possesses a significant partial match with the template scaffold (DECAAF). A library of partial motifs, constructed from the active site residues of the template protein, is used to rank a set of target proteins based on maximal significant matches with the partial motifs, and cull out the best candidate from the reduced set as the target protein. Considering the scenario where this 'incubator' protein lacks activity, we identify mutations in the target protein that will mirror the template motif by superimposing the target and template protein based on the partial match. Using this superimposition technique, we analyzed the less than expected gain of activity achieved by an attempt to induce β-lactamase activity in a penicillin binding protein (PBP) (PBP-A from T. elongatus), and attributed this to steric hindrance from neighboring residues. We also propose mutations in PBP-5 from E. coli, which does not have similar steric constraints. The flow details have been worked out in an example which aims to select a substitute protein for human neutrophil elastase, preferably related to grapevines, in a chimeric anti-microbial enzyme which bolsters the innate immune defense system of grapevines.     

Shared promiscuous activities and evolutionary features in various members of the amidohydrolase superfamily

The amidohydrolase superfamily comprises hundreds of hydrolytic enzymes of the (beta/alpha)8 barrel fold with mono- or binuclear active-site metal centers, and a diverse spectrum of substrates and reactions. Promiscuous activities, or cross-reactivities, between different members of the same superfamily may provide important hints regarding evolutionary and mechanistic relationships. We examined three members: dihydroorotase (DHO), phosphotriesterase (PTE), and PTE-homology protein (PHP). Of particular interest are PTE, which is thought to have evolved within the last several decades, and PHP, an amidohydrolase superfamily member of unknown function, and the closest known homologue of PTE. We found a diverse and partially overlapping pattern of promiscuous activities in these enzymes, including a significant lactonase activity in PTE, esterase activities in both PTE and PHP, and a weak PTE activity in DHO. Directed evolution was applied to improve the promiscuous esterase activities of PTE and PHP. Remarkably, the most recurrent mutation increasing esterase activity in PTE, or PHP, maps to the same location in their superposed 3D structures. The evolved variants also exhibit newly acquired promiscuous activities that were not selected for, including very weak, yet measurable, paraoxonase activity in PHP. Our results illustrate the mechanistic, structural, and evolutionary links between these enzymes, and highlight the importance of studying laboratory evolution intermediates that might resemble node intermediates along the evolutionary pathways leading to the divergence of enzyme superfamilies.     

Evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase family of the enolase superfamily

Understanding how proteins evolve to provide both exquisite specificity and proficient activity is a fundamental problem in biology that has implications for protein function prediction and protein engineering. To study this problem, we analyzed the evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase (OSBS/NAAAR) family, part of the mechanistically diverse enolase superfamily. Although all characterized members of the family catalyze the OSBS reaction, this family is extraordinarily divergent, with some members sharing <15% identity. In addition, a member of this family, Amycolatopsis OSBS/NAAAR, is promiscuous, catalyzing both dehydration and racemization. Although the OSBS/NAAAR family appears to have a single evolutionary origin, no sequence or structural motifs unique to this family could be identified; all residues conserved in the family are also found in enolase superfamily members that have different functions. Based on their species distribution, several uncharacterized proteins similar to Amycolatopsis OSBS/NAAAR appear to have been transmitted by lateral gene transfer. Like Amycolatopsis OSBS/NAAAR, these might have additional or alternative functions to OSBS because many are from organisms lacking the pathway in which OSBS is an intermediate. In addition to functional differences, the OSBS/NAAAR family exhibits surprising structural variations, including large differences in orientation between the two domains. These results offer several insights into protein evolution. First, orthologous proteins can exhibit significant structural variation, and specificity can be maintained with little conservation of ligand-contacting residues. Second, the discovery of a set of proteins similar to Amycolatopsis OSBS/NAAAR supports the hypothesis that new protein functions evolve through promiscuous intermediates. Finally, a combination of evolutionary, structural, and sequence analyses identified characteristics that might prime proteins, such as Amycolatopsis OSBS/NAAAR, for the evolution of new activities.     

An optimized ATP/PP(i)-exchange assay in 96-well format for screening of adenylation domains for applications in combinatorial biosynthesis

We report a new format for measuring ATP/[(32)P]pyrophosphate exchange in a higher throughput assay of adenylation domains (A-domains) of non-ribosomal peptide synthetases. These enzymes are key specificity determinants in the assembly line biosynthesis of non-ribosomal peptides, an important class of natural products with an activity spectrum ranging from antibiotic to antitumor activities. Our assay in 96-well format allows the rapid measurement of approximately 1000 data points per week as a basis for precise assessment of the kinetics of A-domains. The assay also allows quantitative high-throughput screening of the substrate specificity of A-domains identifying alternative, promiscuous substrates. We show that our assay is able to give high quality data for the T278A mutant of the A-domain of the tyrocidine synthetase module TycA with a 330-fold lower k(cat)/K(M). The large dynamic range of this assay will be useful for the screening of libraries of mutant A-domains. Finally we describe and evaluate a procedure for the high-throughput purification of A-domains in 96-well format for the latter purpose. Our approach will be of utility for mechanistic analysis, substrate profiling and directed evolution of the A-domains, to ultimately enable the combinatorial biosynthesis of non-natural analogues of non-ribosomal peptides that may have potential as alternative drug candidates.     

Neutral genetic drift can alter promiscuous protein functions, potentially aiding functional evolution

Background  

Many of the mutations accumulated by naturally evolving proteins are neutral in the sense that they do not significantly alter a protein's ability to perform its primary biological function. However, new protein functions evolve when selection begins to favor other, "promiscuous" functions that are incidental to a protein's original biological role. If mutations that are neutral with respect to a protein's primary biological function cause substantial changes in promiscuous functions, these mutations could enable future functional evolution.     

Divergence and convergence in enzyme evolution: parallel evolution of paraoxonases from quorum-quenching lactonases

We discuss the basic features of divergent versus convergent evolution and of the common scenario of parallel evolution. The example of quorum-quenching lactonases is subsequently described. Three different quorum-quenching lactonase families are known, and they belong to three different superfamilies. Their key active-site architectures have converged and are strikingly similar. Curiously, a promiscuous organophosphate hydrolase activity is observed in all three families. We describe the structural and mechanistic features that underline this converged promiscuity and how this promiscuity drove the parallel divergence of organophosphate hydrolases within these lactonase families by either natural or laboratory evolution.     

Phylogenetic analysis of the corbiculate Apinae based on morphology of the sting apparatus (Hymenoptera: Apidae)

This study aimed to test the various competing hypotheses regarding the relationships among the four tribes of corbiculate apine bees (Euglossini “orchid bees”, Bombini “bumble bees”, Meliponini “stingless bees”, and Apini “honey bees”) with a completely new set of previously unstudied morphological characters derived from the sting apparatus. The result was one most parsimonious tree of 49 steps, CI = 89, RI = 93 that is perfectly congruent with most studies based on morphological and combined morphological/molecular data, i.e., Euglossini + (Bombini + (Meliponini + Apini)), supporting a well accepted scenario of social evolution for these bees. This data matrix was then combined with other published matrices for this group in order to perform simultaneous analyses. The problem of how to best combine the multiple matrices that did not use the same exemplars was investigated. © The Willi Hennig Society 2007.     

Directed evolution of the promiscuous esterase activity of carbonic anhydrase II

A promiscuous activity of an existing enzyme can confer an evolutionary advantage by providing an immediate response to a new selection pressure and a starting point for the divergence of a new enzyme. This work seeks to examine how this process might take place. Human carbonic anhydrase II (hCAII) is an enzyme that evolved to catalyze the reversible hydration of CO(2) and performs this task at a remarkable rate (k(cat) approximately 10(6) s(-)(1)). hCAII also exhibits promiscuous activity toward highly activated esters such as 4-nitrophenyl acetate. We describe a much weaker esterase activity of hCAII toward the bulkier and much less activated ester substrate 2-naphthyl acetate (2NA). Directed evolution of hCAII produced a variant with 40-fold higher rates toward 2NA, owing to two mutations: one within the active site (Ala65Val) and one at its mouth (Thr200Ala). Structure-activity studies suggest that these mutations led to adaptation of the active site for bulkier substrates and for the catalysis of nonactivated esters. The mutations did not, however, significantly alter the native activity of hCAII. Our results support the notion that the evolution of a new function can be driven by mutations that increase a promiscuous function (which serves as the starting point for the evolutionary process) but do not harm the native function.     

Enhancing the Promiscuous Phosphotriesterase Activity of a Thermostable Lactonase (GkaP) for the Efficient Degradation of Organophosphate Pesticides

The phosphotriesterase-like lactonase (PLL) enzymes in the amidohydrolase superfamily hydrolyze various lactones and exhibit latent phosphotriesterase activities. These enzymes serve as attractive templates for in vitro evolution of neurotoxic organophosphates (OPs) with hydrolytic capabilities that can be used as bioremediation tools. Here, a thermostable PLL from Geobacillus kaustophilus HTA426 (GkaP) was targeted for joint laboratory evolution with the aim of enhancing its catalytic efficiency against OP pesticides. By a combination of site saturation mutagenesis and whole-gene error-prone PCR approaches, several improved variants were isolated. The most active variant, 26A8C, accumulated eight amino acid substitutions and demonstrated a 232-fold improvement over the wild-type enzyme in reactivity (k(cat)/K(m)) for the OP pesticide ethyl-paraoxon. Concomitantly, this variant showed a 767-fold decrease in lactonase activity with δ-decanolactone, imparting a specificity switch of 1.8 × 10(5)-fold. 26A8C also exhibited high hydrolytic activities (19- to 497-fold) for several OP pesticides, including parathion, diazinon, and chlorpyrifos. Analysis of the mutagenesis sites on the GkaP structure revealed that most mutations are located in loop 8, which determines substrate specificity in the amidohydrolase superfamily. Molecular dynamics simulation shed light on why 26A8C lost its native lactonase activity and improved the promiscuous phosphotriesterase activity. These results permit us to obtain further insights into the divergent evolution of promiscuous enzymes and suggest that laboratory evolution of GkaP may lead to potential biological solutions for the efficient decontamination of neurotoxic OP compounds.     

Computational methods for sequence mapping of large combinatorial libraries and deduced sequence signatures

Here we describe a computational approach for the high-throughput sequence mapping of combinatorial libraries obtained by DNA shuffling. Original algorithms and their software implementation were developed for the automated and reliable analysis of hybridization data of differentially labeled oligonucleotide probes with PCR products spotted on DNA microarrays. This novel approach allows a context-dependent sequence attribution tolerant to fluctuations in experimental conditions and is well adapted to hybridization signals of variable qualities resulting from high-throughput PCR amplification from colonies. In addition, the analysis permits the calculation of sequence signatures that are characteristic of combinatorial library structure, defects, and diversity. The approach is of interest for the characterization and the equalization (library reduction to nonredundant structures) of combinatorial libraries involved in directed evolution and could be extrapolated to high-throughput polymorphism analysis.     

A quantitative index of substrate promiscuity

Catalytic promiscuity is a widespread, but poorly understood, phenomenon among enzymes with particular relevance to the evolution of new functions, drug metabolism, and in vitro biocatalyst engineering. However, there is at present no way to quantitatively measure or compare this important parameter of enzyme function. Here we define a quantitative index of promiscuity (I) that can be calculated from the catalytic efficiencies of an enzyme toward a defined set of substrates. A weighted promiscuity index (J) that accounts for patterns of similarity and dissimilarity among the substrates in the set is also defined. Promiscuity indices were calculated for three different enzyme classes: eight serine and cysteine proteases, two glutathione S-transferase (GST) isoforms, and three cytochrome P450 (CYP) isoforms. The proteases ranged from completely specific (granzyme B, J = 0.00) to highly promiscuous (cruzain, J = 0.83). The four drug-metabolizing enzymes studied (GST A1-1 and the CYP isoforms) were highly promiscuous, with J values between 0.72 and 0.92; GST A4-4, involved in the clearance of lipid peroxidation products, is moderately promiscuous (J = 0.37). Promiscuity indices also allowed for studies of correlation between substrate promiscuity and an enzyme's activity toward its most-favored substrate, for each of the three enzyme classes.