Evolutionary potential of (beta/alpha)8-barrels: functional promiscuity produced by single substitutions in the enolase superfamily

Computational redesign of a mononuclear zinc metalloenzyme for organophosphate hydrolysis

The ability to redesign enzymes to catalyze noncognate chemical transformations would have wide-ranging applications. We developed a computational method for repurposing the reactivity of metalloenzyme active site functional groups to catalyze new reactions. Using this method, we engineered a zinc-containing mouse adenosine deaminase to catalyze the hydrolysis of a model organophosphate with a catalytic efficiency (k(cat)/K(m)) of ~10(4) M(-1) s(-1) after directed evolution. In the high-resolution crystal structure of the enzyme, all but one of the designed residues adopt the designed conformation. The designed enzyme efficiently catalyzes the hydrolysis of the R(P) isomer of a coumarinyl analog of the nerve agent cyclosarin, and it shows marked substrate selectivity for coumarinyl leaving groups. Computational redesign of native enzyme active sites complements directed evolution methods and offers a general approach for exploring their untapped catalytic potential for new reactivities.

QM/MM analysis suggests that Alkaline Phosphatase (AP) and nucleotide pyrophosphatase/phosphodiesterase slightly tighten the transition state for phosphate diester hydrolysis relative to solution: implication for catalytic promiscuity in the AP superfamily

Several members of the Alkaline Phosphatase (AP) superfamily exhibit a high level of catalytic proffciency and promiscuity in structurally similar active sites. A thorough characterization of the nature of transition state for different substrates in these enzymes is crucial for understanding the molecular mechanisms that govern those remarkable catalytic properties. In this work, we study the hydrolysis of a phosphate diester, MpNPP(-), in solution, two experimentally well-characterized variants of AP (R166S AP, R166S/E322Y AP) and wild type Nucleotide pyrophosphatase/phosphodiesterase (NPP) by QM/MM calculations in which the QM method is an approximate density functional theory previously parametrized for phosphate hydrolysis (SCC-DFTBPR). The general agreements found between these calculations and available experimental data for both solution and enzymes support the use of SCC-DFTBPR/MM for a semiquantitative analysis of the catalytic mechanism and nature of transition state in AP and NPP. Although phosphate diesters are cognate substrates for NPP but promiscuous substrates for AP, the calculations suggest that their hydrolysis reactions catalyzed by AP and NPP feature similar synchronous transition states that are slightly tighter in nature compared to that in solution, due in part to the geometry of the bimetallic zinc motif. Therefore, this study provides the first direct computational support to the hypothesis that enzymes in the AP superfamily catalyze cognate and promiscuous substrates via similar transition states to those in solution. Our calculations do not support the finding of recent QM/MM studies by López-Canut and co-workers, who suggested that the same diester substrate goes through a much looser transition state in NPP/AP than in solution, a result likely biased by the large structural distortion of the bimetallic zinc site in their simulations. Finally, our calculations for different phosphate diester orientations and phosphorothioate diesters highlight that the interpretation of thio-substitution experiments is not always straightforward.

Miniaturized, microarray-based assays for chemical proteomic studies of protein function

Systematic analysis of protein and enzyme function typically requires scale-up of protein expression and purification prior to assay development; this can often be limiting. Miniaturization of assays provides an alternative approach, but simple, generic methods are in short supply. Here we show how custom microarrays can be adapted to this purpose. We discuss the different routes to array fabrication and describe in detail one facile approach in which the purification and immobilization procedures are combined into a single step, significantly simplifying the array fabrication process. We illustrate this approach by reference to the creation of arrays of human protein kinases and of human cytochrome P450s. We discuss methods for both ligand-binding and turnover-based assays, as well as data analysis on such arrays.

Divergent evolution in enolase superfamily: strategies for assigning functions

Nature's strategies for evolving catalytic functions can be deciphered from the information contained in the rapidly expanding protein sequence databases. However, the functions of many proteins in the protein sequence and structure databases are either uncertain (too divergent to assign function based on homology) or unknown (no homologs), thereby limiting the utility of the databases. The mechanistically diverse enolase superfamily is a paradigm for understanding the structural bases for evolution of enzymatic function. We describe strategies for assigning functions to members of the enolase superfamily that should be applicable to other superfamilies.

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.

Assessing directed evolution methods for the generation of biosynthetic enzymes with potential in drug biosynthesis

To address the synthesis of increasingly structurally diverse small-molecule drugs, methods for the generation of efficient and selective biological catalysts are becoming increasingly important. 'Directed evolution' is an umbrella term referring to a variety of methods for improving or altering the function of enzymes using a nature-inspired twofold strategy of mutagenesis followed by selection. This article provides an objective assessment of the effectiveness of directed evolution campaigns in generating enzymes with improved catalytic parameters for new substrates from the last decade, excluding studies that aimed to select for only improved physical properties and those that lack kinetic characterization. An analysis of the trends of methodologies and their success rates from 81 qualifying examples in the literature reveals the average fold improvement for k (cat) (or V (max)), K (m) and k (cat)/K (m) to be 366-, 12- and 2548-fold, respectively, whereas the median fold improvements are 5.4, 3 and 15.6. Further analysis by enzyme class, library-generation methodology and screening methodology explores relationships between successful campaigns and the methodologies employed.

Constructing de novo biosynthetic pathways for chemical synthesis inside living cells

Living organisms have evolved a vast array of catalytic functions that make them ideally suited for the production of medicinally and industrially relevant small-molecule targets. Indeed, native metabolic pathways in microbial hosts have long been exploited and optimized for the scalable production of both fine and commodity chemicals. Our increasing capacity for DNA sequencing and synthesis has revealed the molecular basis for the biosynthesis of a variety of complex and useful metabolites and allows the de novo construction of novel metabolic pathways for the production of new and exotic molecular targets in genetically tractable microbes. However, the development of commercially viable processes for these engineered pathways is currently limited by our ability to quickly identify or engineer enzymes with the correct reaction and substrate selectivity as well as the speed by which metabolic bottlenecks can be determined and corrected. Efforts to understand the relationship among sequence, structure, and function in the basic biochemical sciences can advance these goals for synthetic biology applications while also serving as an experimental platform for elucidating the in vivo specificity and function of enzymes and reconstituting complex biochemical traits for study in a living model organism. Furthermore, the continuing discovery of natural mechanisms for the regulation of metabolic pathways has revealed new principles for the design of high-flux pathways with minimized metabolic burden and has inspired the development of new tools and approaches to engineering synthetic pathways in microbial hosts for chemical production.

Protein function prediction: towards integration of similarity metrics

Genomic centers discover increasingly many protein sequences and structures, but not necessarily their full biological functions. Thus, currently, less than one percent of proteins have experimentally verified biochemical activities. To fill this gap, function prediction algorithms apply metrics of similarity between proteins on the premise that those sufficiently alike in sequence, or structure, will perform identical functions. Although high sensitivity is elusive, network analyses that integrate these metrics together hold the promise of rapid gains in function prediction specificity.

Systematic screening for catalytic promiscuity in 4-oxalocrotonate tautomerase: enamine formation and aldolase activity

The enzyme 4-oxalocrotonate tautomerase (4-OT) is part of a catabolic pathway for aromatic hydrocarbons in Pseudomonas putida mt-2, where it catalyzes the conversion of 2-hydroxy-2,4-hexadienedioate(1) to 2-oxo-3-hexenedioate(2). 4-OT is a member of the tautomerase superfamily, a group of homologous proteins that are characterized by a β-α-β structural fold and a catalytic amino-terminal proline. In the mechanism of 4-OT, Pro1 is a general base that abstracts the 2-hydroxyl proton of 1 for delivery to the C-5 position to yield 2. Here, 4-OT was explored for nucleophilic catalysis based on the mechanistic reasoning that its Pro1 residue has the correct protonation state (pK(a) ∼6.4) to be able to act as a nucleophile at pH 7.3. By using inhibition studies and mass spectrometry experiments it was first demonstrated that 4-OT can use Pro1 as a nucleophile to form an imine/enamine with various aldehyde and ketone compounds. The chemical potential of the smallest enamine (generated from acetaldehyde) was then explored for further reactions by using a small set of selected electrophiles. This systematic screening approach led to the discovery of a new promiscuous activity in wild-type 4-OT: the enzyme catalyzes the aldol condensation of acetaldehyde with benzaldehyde to form cinnamaldehyde. This low-level aldolase activity can be improved 16-fold with a single point mutation (L8R) in 4-OT's active site. The proposed mechanism of the reaction mimicks that used by natural class-I aldolases and designed catalytic aldolase antibodies. An important difference, however, is that these natural and designed aldolases use the primary amine of a lysine residue to form enamines with carbonyl substrates, whereas 4-OT uses the secondary amine of an active-site proline as the nucleophile catalyst. Further systematic screening of 4-OT and related proline-based biocatalysts might prove to be a useful approach to discover new promiscuous carbonyl transformation activities that could be exploited to develop new biocatalysts for carbon-carbon bond formation.

Protein function microarrays for customised systems-oriented proteome analysis

Protein microarrays have many potential applications in the systematic, quantitative analysis of protein function. However, simple, reproducible, and robust methods for array fabrication that are compatible with the study of large, custom collections of potentially unrelated proteins are required. Here, we discuss different routes to array fabrication and describe in detail one approach in which the purification and immobilisation procedures are combined into a single step, significantly simplifying the array fabrication process. We illustrate this approach by reference to the creation of an array of human protein kinases and discuss methods for assay and data analysis on such arrays.

Enhancement of the latent 3-isopropylmalate dehydrogenase activity of promiscuous homoisocitrate dehydrogenase by directed evolution

HICDH (homoisocitrate dehydrogenase), which is involved in lysine biosynthesis through α-aminoadipate, is a paralogue of IPMDH [3-IPM (3-isopropylmalate) dehydrogenase], which is involved in leucine biosynthesis. TtHICDH (Thermus thermophilus HICDH) can recognize isocitrate, as well as homoisocitrate, as the substrate, and also shows IPMDH activity, although at a considerably decreased rate. In the present study, the promiscuous TtHICDH was evolved into an enzyme showing distinct IPMDH activity by directed evolution using a DNA-shuffling technique. Through five repeats of DNA shuffling/screening, variants that allowed Escherichia coli C600 (leuB−) to grow on a minimal medium in 2 days were obtained. One of the variants LR5–1, with eight amino acid replacements, was found to possess a 65-fold increased kcat/Km value for 3-IPM, compared with TtHICDH. Introduction of a single back-replacement H15Y change caused a further increase in the kcat/Km value and a partial recovery of the decreased thermotolerance of LR5–1. Site-directed mutagenesis revealed that most of the amino acid replacements found in LR5–1 effectively increased IPMDH activity; replacements around the substrate-binding site contributed to the improved recognition for 3-IPM, and other replacements at sites away from the substrate-binding site enhanced the turnover number for the IPMDH reaction. The crystal structure of LR5–1 was determined at 2.4 Å resolution and revealed that helix α4 was displaced in a manner suitable for recognition of the hydrophobic γ-moiety of 3-IPM. On the basis of the crystal structure, possible reasons for enhancement of the turnover number are discussed.

Kinetic and structural characterization of a heterohexamer 4-oxalocrotonate tautomerase from Chloroflexus aurantiacus J-10-fl: implications for functional and structural diversity in the tautomerase superfamily

4-Oxalocrotonate tautomerase (4-OT) isozymes play prominent roles in the bacterial utilization of aromatic hydrocarbons as sole carbon sources. These enzymes catalyze the conversion of 2-hydroxy-2,4-hexadienedioate (or 2-hydroxymuconate) to 2-oxo-3-hexenedioate, where Pro-1 functions as a general base and shuttles a proton from the 2-hydroxyl group of the substrate to the C-5 position of the product. 4-OT, a homohexamer from Pseudomonas putida mt-2, is the most extensively studied 4-OT isozyme and the founding member of the tautomerase superfamily. A search of five thermophilic bacterial genomes identified a coded amino acid sequence in each that had been annotated as a tautomerase-like protein but lacked Pro-1. However, a nearby sequence has Pro-1, but the sequence is not annotated as a tautomerase-like protein. To characterize this group of proteins, two genes from Chloroflexus aurantiacus J-10-fl were cloned, and the corresponding proteins were expressed. Kinetic, biochemical, and X-ray structural analyses show that the two expressed proteins form a functional heterohexamer 4-OT (hh4-OT), composed of three alphabeta dimers. Like the P. putida enzyme, hh4-OT requires the amino-terminal proline and two arginines for the conversion of 2-hydroxymuconate to the product, implicating an analogous mechanism. In contrast to 4-OT, hh4-OT does not exhibit the low-level activity of another tautomerase superfamily member, the heterohexamer trans-3-chloroacrylic acid dehalogenase (CaaD). Characterization of hh4-OT enables functional assignment of the related enzymes, highlights the diverse ways the beta-alpha-beta building block can be assembled into an active enzyme, and provides further insight into the molecular basis of the low-level CaaD activity in 4-OT.

Enzyme promiscuity: a mechanistic and evolutionary perspective

Many, if not most, enzymes can promiscuously catalyze reactions, or act on substrates, other than those for which they evolved. Here, we discuss the structural, mechanistic, and evolutionary implications of this manifestation of infidelity of molecular recognition. We define promiscuity and related phenomena and also address their generality and physiological implications. We discuss the mechanistic enzymology of promiscuity--how enzymes, which generally exert exquisite specificity, catalyze other, and sometimes barely related, reactions. Finally, we address the hypothesis that promiscuous enzymatic activities serve as evolutionary starting points and highlight the unique evolutionary features of promiscuous enzyme functions.

Mechanisms of protein evolution and their application to protein engineering

Protein engineering holds great promise for the development of new biosensors, diagnostics, therapeutics, and agents for bioremediation. Despite some remarkable successes in experimental and computational protein design, engineered proteins rarely achieve the efficiency or specificity of natural enzymes. Current protein design methods utilize evolutionary concepts, including mutation, recombination, and selection, but the inability to fully recapitulate the success of natural evolution suggests that some evolutionary principles have not been fully exploited. One aspect of protein engineering that has received little attention is how to select the most promising proteins to serve as templates, or scaffolds, for engineering. Two evolutionary concepts that could provide a rational basis for template selection are the conservation of catalytic mechanisms and functional promiscuity. Knowledge of the catalytic motifs responsible for conserved aspects of catalysis in mechanistically diverse superfamilies could be used to identify promising templates for protein engineering. Second, protein evolution often proceeds through promiscuous intermediates, suggesting that templates which are naturally promiscuous for a target reaction could enhance protein engineering strategies. This review explores these ideas and alternative hypotheses concerning protein evolution and engineering. Future research will determine if application of these principles will lead to a protein engineering methodology governed by predictable rules for designing efficient, novel catalysts.

Engineering monolignol 4-O-methyltransferases to modulate lignin biosynthesis

Lignin is a complex polymer derived from the oxidative coupling of three classical monolignols. Lignin precursors are methylated exclusively at the meta-positions (i.e. 3/5-OH) of their phenyl rings by native O-methyltransferases, and are precluded from substitution of the para-hydroxyl (4-OH) position. Ostensibly, the para-hydroxyls of phenolics are critically important for oxidative coupling of phenoxy radicals to form polymers. Therefore, creating a 4-O-methyltransferase to substitute the para-hydroxyl of monolignols might well interfere with the synthesis of lignin. The phylogeny of plant phenolic O-methyltransferases points to the existence of a batch of evolutionarily "plastic" amino acid residues. Following one amino acid at a time path of directed evolution, and using the strategy of structure-based iterative site-saturation mutagenesis, we created a novel monolignol 4-O-methyltransferase from the enzyme responsible for methylating phenylpropenes. We show that two plastic residues in the active site of the parental enzyme are vital in dominating substrate discrimination. Mutations at either one of these separate the evolutionarily tightly linked properties of substrate specificity and regioselective methylation of native O-methyltransferase, thereby conferring the ability for para-methylation of the lignin monomeric precursors, primarily monolignols. Beneficial mutations at both sites have an additive effect. By further optimizing enzyme activity, we generated a triple mutant variant that may structurally constitute a novel phenolic substrate binding pocket, leading to its high binding affinity and catalytic efficiency on monolignols. The 4-O-methoxylation of monolignol efficiently impairs oxidative radical coupling in vitro, highlighting the potential for applying this novel enzyme in managing lignin polymerization in planta.

Directed evolution of a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 in the amidohydrolase superfamily

The PLL(PTE-like lactonase)-group of enzymes within the amidohydrolase superfamily hydrolyze N-acyl-homoserine lactones (AHLs) that are involved in bacterial quorum-sensing pathways. These enzymes possess the (beta/alpha)(8)-barrel fold and serve as attractive templates for in vitro evolution and engineering of quorum-quenching biological molecules that can serve as antivirulence therapeutic agents. Using a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 (GI: 41409766) as the initial template for in vitro evolution experiments, we enhanced the catalytic efficiency and increased the substrate range of the wild-type enzyme through a single point mutation on the loop at the C-terminal end of the eighth beta-strand. This N266Y mutant had an increased value of k(cat)/K(M) of 30- and 32-fold toward 3-oxo-N-octanoyl-l-homoserine lactone and N-hexanoyl-l-homoserine lactone, respectively; the evolved mutant also exhibited lactonase activity toward 3-oxo-N-hexanoyl-l-homoserine lactone and N-butyryl-l-homoserine lactone, AHLs that were previously not hydrolyzed by the wild-type enzyme. This article reinforces the evolutionary potential of the (beta/alpha)(8)-barrel fold and highlights the possibility of using quorum-quenching lactonases in the amidohydrolase superfamily as templates for engineering biomolecules of therapeutic use.

Evolution of enzymatic activities in the enolase superfamily: stereochemically distinct mechanisms in two families of cis,cis-muconate lactonizing enzymes

The mechanistically diverse enolase superfamily is a paradigm for elucidating Nature's strategies for divergent evolution of enzyme function. Each of the different reactions catalyzed by members of the superfamily is initiated by abstraction of the alpha-proton of a carboxylate substrate that is coordinated to an essential Mg(2+). The muconate lactonizing enzyme (MLE) from Pseudomonas putida, a member of a family that catalyzes the syn-cycloisomerization of cis,cis-muconate to (4S)-muconolactone in the beta-ketoadipate pathway, has provided critical insights into the structural bases for evolution of function within the superfamily. A second, divergent family of homologous MLEs that catalyzes anti-cycloisomerization has been identified. Structures of members of both families liganded with the common (4S)-muconolactone product (syn, Pseudomonas fluorescens, gi 70731221 ; anti, Mycobacterium smegmatis, gi 118470554 ) document that the conserved Lys at the end of the second beta-strand in the (beta/alpha)(7)beta-barrel domain serves as the acid catalyst in both reactions. The different stereochemical courses (syn and anti) result from different structural strategies for determining substrate specificity: although the distal carboxylate group of the cis,cis-muconate substrate attacks the same face of the proximal double bond, opposite faces of the resulting enolate anion intermediate are presented to the conserved Lys acid catalyst. The discovery of two families of homologous, but stereochemically distinct, MLEs likely provides an example of "pseudoconvergent" evolution of the same function from different homologous progenitors within the enolase superfamily, in which different spatial arrangements of active site functional groups and substrate specificity determinants support catalysis of the same reaction.

Enzyme (re)design: lessons from natural evolution and computation

The (re)design of enzymes to catalyze 'new' reactions is a topic of considerable practical and intellectual interest. Directed evolution (random mutagenesis followed by screening/selection) has been used widely to identify novel biocatalysts. However, 'rational' approaches using either natural divergent evolution or computational predictions based on chemical principles have been less successful. This review summarizes recent progress in evolution-based and computation-based (re)design.

Establishing wild-type levels of catalytic activity on natural and artificial (beta alpha)8-barrel protein scaffolds

The generation of high levels of new catalytic activities on natural and artificial protein scaffolds is a major goal of enzyme engineering. Here, we used random mutagenesis and selection in vivo to establish a sugar isomerisation reaction on both a natural (beta alpha)(8)-barrel enzyme and a catalytically inert chimeric (beta alpha)(8)-barrel scaffold, which was generated by the recombination of 2 (beta alpha)(4)-half barrels. The best evolved variants show turnover numbers and substrate affinities that are similar to those of wild-type enzymes catalyzing the same reaction. The determination of the crystal structure of the most proficient variant allowed us to model the substrate sugar in the novel active site and to elucidate the mechanistic basis of the newly established activity. The results demonstrate that natural and inert artificial protein scaffolds can be converted into highly proficient enzymes in the laboratory, and provide insights into the mechanisms of enzyme evolution.

The C-terminal domain of Escherichia coli dihydrodipicolinate synthase (DHDPS) is essential for maintenance of quaternary structure and efficient catalysis

Dihydrodipicolinate synthase (DHDPS) catalyses the first committed step in the biosynthesis of (S)-lysine, an essential constituent of bacterial cell walls. Escherichia coli DHDPS is homotetrameric, and each monomer contains an N-terminal (alpha/beta)(8)-barrel, responsible for catalysis and regulation, and three C-terminal alpha-helices, the function of which is unknown. This study investigated the C-terminal domain of E. coli DHDPS by characterising a C-terminal truncated DHDPS (DHDPS-H225*). DHDPS-H225* was unable to complement an (S)-lysine auxotroph, and showed significantly reduced solubility, stability, and maximum catalytic activity (k(cat)=1.20+/-0.01 s(-1)), which was only 1.6% of wild type E. coli DHDPS (DHDPS-WT). The affinity of DHDPS-H225* for substrates and the feedback inhibitor, (S)-lysine, remained comparable to DHDPS-WT. These changes were accompanied by disruption in the quaternary structure, which has previously been shown to be essential for efficient catalysis in this enzyme.

Mapping catalytic promiscuity in the alkaline phosphatase superfamily

"Promiscuous" enzymes possess activities in addition to their native ones. Promiscuous activities could be remnants from an evolutionary ancestor that has been adapted to fulfil a new function following gene duplication. Alternatively, the observation of promiscuity could indicate that an enzyme has the potential to evolve into a new catalyst. Thus, the observation of promiscuity defines functional relationships in enzyme superfamilies. Crosswise promiscuity can provide an additional layer of connectivity between members of a - usually structurally defined - superfamily to establish a system for tracking the emergence and interconversion of enzymatic function. The systematic analysis of measured promiscuous rates may serve as a basis for drawing up phylogenetic relationships based on the potential for catalysis and may be useful for active use in directed evolution, suggesting evolutionary "short cuts". We review recent observations of catalytic promiscuity in members of the alkaline phosphatase (AP) superfamily that exhibit reciprocal relationships of crosswise promiscuity with rate accelerations (kcat/KM)/k2 between 106 and 1018. Specifically, we focus on the mechanistic features that appear to form the basis of catalytic promiscuity in this superfamily.

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.

Evolution of function in the "two dinucleotide binding domains" flavoproteins

Structural and biochemical constraints force some segments of proteins to evolve more slowly than others, often allowing identification of conserved structural or sequence motifs that can be associated with substrate binding properties, chemical mechanisms, and molecular functions. We have assessed the functional and structural constraints imposed by cofactors on the evolution of new functions in a superfamily of flavoproteins characterized by two-dinucleotide binding domains, the “two dinucleotide binding domains” flavoproteins (tDBDF) superfamily. Although these enzymes catalyze many different types of oxidation/reduction reactions, each is initiated by a stereospecific hydride transfer reaction between two cofactors, a pyridine nucleotide and flavin adenine dinucleotide (FAD). Sequence and structural analysis of more than 1,600 members of the superfamily reveals new members and identifies details of the evolutionary connections among them. Our analysis shows that in all of the highly divergent families within the superfamily, these cofactors adopt a conserved configuration optimal for stereospecific hydride transfer that is stabilized by specific interactions with amino acids from several motifs distributed among both dinucleotide binding domains. The conservation of cofactor configuration in the active site restricts the pyridine nucleotide to interact with FAD from the re-side, limiting the flow of electrons from the re-side to the si-side. This directionality of electron flow constrains interactions with the different partner proteins of different families to occur on the same face of the cofactor binding domains. As a result, superimposing the structures of tDBDFs aligns not only these interacting proteins, but also their constituent electron acceptors, including heme and iron-sulfur clusters. Thus, not only are specific aspects of the cofactor-directed chemical mechanism conserved across the superfamily, the constraints they impose are manifested in the mode of protein–protein interactions. Overlaid on this foundation of conserved interactions, nature has conscripted different protein partners to serve as electron acceptors, thereby generating diversification of function across the superfamily.

The sequencing of genomes from different species has provided a unique opportunity for comparative analysis and opened the door to a higher level of understanding of living organisms. However, identifying the biochemical functions of the protein products coded by these genes has proved to be a major challenge. Computational methods that have been used to assign functions to such sequences often result in high levels of misannotations. Nature's strategy of evolving new function provides clues for formulating an accurate predictive scheme for functional annotation. Constraints associated with substrate binding properties and chemistry have been shown to be major determinants of guiding the evolution of new function. In this study, the authors have explored the functional and structural constraints imposed by complex cofactors on the evolution of new functions. Analysis of the large “two dinucleotide binding domains” flavoproteins (tDBDF) superfamily using structural comparisons and other bioinformatics approaches shows how structural requirements associated with cofactor reactivity constrain the mode of protein–protein interactions while providing the major route for evolution of functional diversification. The evolutionary framework established in this work may be generally useful for the analysis of functional divergence in other enzyme superfamilies that use complex cofactors.

Directed evolution of novel protein functions

Directed evolution has been successfully used to engineer proteins for basic and applied biological research. However, engineering of novel protein functions by directed evolution remains an overwhelming challenge. This challenge may come from the fact that multiple simultaneously or synergistic mutations are required for the creation of a novel protein function. Here we review the key developments in engineering of novel protein functions by using either directed evolution or a combined directed evolution and rational or computational design approach. Specific attention will be paid to a molecular evolution model for generation of novel proteins. The engineered novel proteins should not only broaden the range of applications of proteins but also provide new insights into protein structure-function relationship and protein evolution.

Determination of the stereochemistry of 2-succinyl-5-enolpyruvyl-6-hydroxy-3- cyclohexene-1-carboxylate, a key intermediate in menaquinone biosynthesis

The turnover product of the committed step of menaquinone biosynthesis was isolated and determined to be (1R,2S,5S,6S)-2-succinyl-5-enolpyruvyl-6-hydroxy-3-cyclohexene-1-carboxylate. Structural determination of this key intermediate represents a critical step to complete elucidation of the biosynthetic pathway.

Multicopy suppression underpins metabolic evolvability

Our understanding of the origins of new metabolic functions is based upon anecdotal genetic and biochemical evidence. Some auxotrophies can be suppressed by overexpressing substrate-ambiguous enzymes (i.e., those that catalyze the same chemical transformation on different substrates). Other enzymes exhibit weak but detectable catalytic promiscuity in vitro (i.e., they catalyze different transformations on similar substrates). Cells adapt to novel environments through the evolution of these secondary activities, but neither their chemical natures nor their frequencies of occurrence have been characterized en bloc. Here, we systematically identified multifunctional genes within the Escherichia coli genome. We screened 104 single-gene knockout strains and discovered that many (20%) of these auxotrophs were rescued by the overexpression of at least one noncognate E. coli gene. The deleted gene and its suppressor were generally unrelated, suggesting that promiscuity is a product of contingency. This genome-wide survey demonstrates that multifunctional genes are common and illustrates the mechanistic diversity by which their products enhance metabolic robustness and evolvability.

Functional Promiscuity Correlates with Conformational Heterogeneity in A-class Glutathione S-Transferases

The structurally related glutathione S-transferase isoforms GSTA1-1 and GSTA4-4 differ greatly in their relative catalytic promiscuity. GSTA1-1 is a highly promiscuous detoxification enzyme. In contrast, GSTA4-4 exhibits selectivity for congeners of the lipid peroxidation product 4-hydroxynonenal. The contribution of protein dynamics to promiscuity has not been studied. Therefore, hydrogen/deuterium exchange mass spectrometry (H/DX) and fluorescence lifetime distribution analysis were performed with glutathione S-transferases A1-1 and A4-4. Differences in local dynamics of the C-terminal helix were evident as expected on the basis of previous studies. However, H/DX demonstrated significantly greater solvent accessibility throughout most of the GSTA1-1 sequence compared with GSTA4-4. A Phe-111/Tyr-217 aromatic-aromatic interaction in A4-4, which is not present in A1-1, was hypothesized to increase core packing. "Swap" mutants that eliminate this interaction from A4-4 or incorporate it into A1-1 yield H/DX behavior that is intermediate between the wild type templates. In addition, the single Trp-21 residue of each isoform was exploited to probe the conformational heterogeneity at the intrasubunit domain-domain interface. Excited state fluorescence lifetime distribution analysis indicates that this core residue is more conformationally heterogeneous in GSTA1-1 than in GSTA4-4, and this correlates with greater stability toward urea denaturation for GSTA4-4. The fluorescence distribution and urea sensitivity of the mutant proteins were intermediate between the wild type templates. The results suggest that the differences in protein dynamics of these homologs are global. The results suggest also the possible importance of extensive conformational plasticity to achieve high levels of functional promiscuity, possibly at the cost of stability.

Minimalist active-site redesign: teaching old enzymes new tricks

Although nature evolves its catalysts over millions of years, enzyme engineers try to do it a bit faster. Enzyme active sites provide highly optimized microenvironments for the catalysis of biologically useful chemical transformations. Consequently, changes at these centers can have large effects on enzyme activity. The prediction and control of these effects provides a promising way to access new functions. The development of methods and strategies to explore the untapped catalytic potential of natural enzyme scaffolds has been pushed by the increasing demand for industrial biocatalysts. This Review describes the use of minimal modifications at enzyme active sites to expand their catalytic repertoires, including targeted mutagenesis and the addition of new reactive functionalities. Often, a novel activity can be obtained with only a single point mutation. The many successful examples of active-site engineering through minimal mutations give useful insights into enzyme evolution and open new avenues in biocatalyst research.

The subtle benefits of being promiscuous: adaptive evolution potentiated by enzyme promiscuity

In this commentary we discuss recent progress in our understanding of adaptive protein evolution. We begin with a general introduction to proteins and their evolution, quickly focusing on the question of how natural selection produces proteins with novel functions. We then summarize the theory of latent protein adaptation advanced by the joint articles by Amitai et al. (2007), HFSP J. 1, 67-78 and Wroe et al. (2007) HFSP J. 1, 79-87, published in the first issue edition of the HFSP Journal. This theory provides a biophysical framework linking the effects of individual mutations on promiscuous protein function, neutral genetic drift, and gene duplication to the process of adaptive protein evolution.

In vivo selection for the directed evolution of L-rhamnulose aldolase from L-rhamnulose-1-phosphate aldolase (RhaD)

Dihydroxyacetone phosphate (DHAP)-dependent aldolases have been widely used for organic synthesis. The major drawback of DHAP-dependent aldolases is their strict donor substrate specificity toward DHAP, which is expensive and unstable. Here we report the development of an in vivo selection system for the directed evolution of the DHAP-dependent aldolase, L-rhamnulose-1-phosphate aldolase (RhaD), to alter its donor substrate specificity from DHAP to dihydroxyacetone (DHA). We also report preliminary results on mutants that were discovered with this screen. A strain deficient in the L-rhamnose metabolic pathway in Escherichia coli (DeltarhaDAB, DE3) was constructed and used as a selection host strain. Co-expression of L-rhamnose isomerase (rhaA) and rhaD in the selection host did not restore its growth on minimal plate supplemented with L-rhamnose as a sole carbon source, because of the lack of L-rhamnulose kinase (RhaB) activity and the inability of WT RhaD aldolase to use unphosphorylated L-rhamnulose as a substrate. Use of this selection host and co-expression vector system gives us an in vivo selection for the desired mutant RhaD which can cleave unphosphorylated L-rhamnulose and allow the mutant to grow in the minimal media. An error-prone PCR (ep-PCR) library of rhaD gene on the co-expression vector was constructed and introduced into the rha-mutant, and survivors were selected in minimal media with l-rhamnose (MMRha media). An initial round of screening gave mutants allowing the selection strain to grow on MMRha plates. This in vivo selection system allows rapid screening of mutated aldolases that can utilize dihydroxyacetone as a donor substrate.

Novel enzymes through design and evolution

The generation of enzymes with new catalytic activities remains a major challenge. So far, several different strategies have been developed to tackle this problem, including site-directed mutagenesis, random mutagenesis (directed evolution), antibody catalysis, computational redesign, and de novo methods. Using these techniques, a broad array of novel enzymes has been created (aldolases, decarboxylases, dehydratases, isomerases, oxidases, reductases, and others), although their low efficiencies (10 to 100 M(-1) s(-l)) compared to those of the best natural enzymes (10(6) to 10(8) M(-1) s(-1)) remains a significant concern. Whereas rational design might be the most promising and versatile approach to generating new activities, directed evolution seems to be the best way to optimize the catalytic properties of novel enzymes. Indeed, impressive successes in enzyme engineering have resulted from a combination of rational and random design.

Enzyme promiscuity: mechanism and applications

Introductory courses in biochemistry teach that enzymes are specific for their substrates and the reactions they catalyze. Enzymes diverging from this statement are sometimes called promiscuous. It has been suggested that relaxed substrate and reaction specificities can have an important role in enzyme evolution; however, enzyme promiscuity also has an applied aspect. Enzyme condition promiscuity has, for a long time, been used to run reactions under conditions of low water activity that favor ester synthesis instead of hydrolysis. Together with enzyme substrate promiscuity, it is exploited in numerous synthetic applications, from the laboratory to industrial scale. Furthermore, enzyme catalytic promiscuity, where enzymes catalyze accidental or induced new reactions, has begun to be recognized as a valuable research and synthesis tool. Exploiting enzyme catalytic promiscuity might lead to improvements in existing catalysts and provide novel synthesis pathways that are currently not available.

Pathway engineering by designed divergent evolution

Designed divergent evolution is a proposed protein engineering methodology to redesign enzyme function. The methodology was developed on the basis of the theories of divergent molecular evolution: (i) enzymes with more active and specialized functions have evolved from ones with promiscuous functions; (ii) this process is driven by small numbers of amino acid substitutions (plasticity); and (iii) the effects of double or multiple mutations are often additive (quasi-additive assumption). Thus, in many cases the impact of multiple mutations can be calculated by first determining the effects of a mutation at a single position and subsequently summing these effects using the quasi-additive assumption. In this way, the shape of the fitness landscape of a particular enzyme function can be estimated. The combinations of mutations predicted to yield global optima for desired functions can then be selected and introduced into the enzymes. The methodology has been demonstrated to be very powerful to redesign enzyme function. The use of multiple redesigned enzymes in novel or reconstructed metabolic pathways will enable the production of natural and unnatural products that will find use as pharmaceuticals, agrochemicals and many other applications.

Evolution of enzymatic activity in the tautomerase superfamily: mechanistic and structural consequences of the L8R mutation in 4-oxalocrotonate tautomerase

4-Oxalocrotonate tautomerase (4-OT) and trans-3-chloroacrylic acid dehalogenase (CaaD) are members of the tautomerase superfamily, a group of structurally homologous proteins that share a beta-alpha-beta fold and a catalytic amino-terminal proline. 4-OT, from Pseudomonas putida mt-2, catalyzes the conversion of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate through the dienol intermediate 2-hydroxymuconate in a catabolic pathway for aromatic hydrocarbons. CaaD, from Pseudomonas pavonaceae 170, catalyzes the hydrolytic dehalogenation of trans-3-chloroacrylate in the trans-1,3-dichloropropene degradation pathway. Both reactions may involve an arginine-stabilized enediolate intermediate, a capability that may partially account for the low-level CaaD activity of 4-OT. Two active-site residues in 4-OT, Leu-8 and Ile-52, have now been mutated to the positionally conserved and catalytic ones in CaaD, alphaArg-8, and alphaGlu-52. The L8R and L8R/I52E mutants show improved CaaD activity (50- and 32-fold increases in k(cat)/K(m), respectively) and diminished 4-OT activity (5- and 1700-fold decreases in k(cat)/K(m), respectively). The increased efficiency of L8R-4-OT for the CaaD reaction stems primarily from an 8.8-fold increase in k(cat), whereas that of the L8R/I52E mutant is due largely to a 23-fold decrease in K(m). The presence of the additional arginine residue in the active site of L8R-4-OT does not alter the pK(a) of the Pro-1 amino group from that measured for the wild type (6.5 +/- 0.1 versus 6.4 +/- 0.2). Moreover, the crystal structure of L8R-4-OT is comparable to that of the wild type. Hence, the enhanced CaaD activity of L8R-4-OT is likely due to the additional arginine residue that can participate in substrate binding and/or stabilization of the putative enediolate intermediate. The results also suggest that the evolution of new functions within the tautomerase superfamily could be quite facile, requiring only a few strategically placed active-site mutations.

Mutation of Lys242 allows Delta3-Delta2-enoyl-CoA isomerase to acquire enoyl-CoA hydratase activity

We report here a novel example of generating hydratase activity through site-directed mutagenesis of a single residue Lys242 of rat liver mitochondrial Delta3-Delta2-enoyl-CoA isomerase, which is one of the key enzymes involved in fatty acid oxidation and a member of the crotonase superfamily. Lys242 is at the C-terminal of the enzyme, which is far from the active site in the crotonase superfamily and forms a salt bridge with Asp149. A variety of mutant expression plasmids were constructed, and it was observed that mutation of Lys242 to nonbasic residues allowed the mutants to have enoyl-CoA hydratase activity confirmed by HPLC analysis of the incubation mixture. Kinetic studies of these mutants were carried out for both isomerase and hydratase activities. Mutant K242C showed a k(cat) value of 1.0 s(-1) for hydration reaction. This activity constitutes about 10% of the total enzyme activity, and the remaining 90% is its natural isomerase activity. To the best of our knowledge, this is the first report on the generation of functional promiscuity through single amino acid mutation far from the active site. This may be a simple and efficient approach to designing a new enzyme based on an existing template. It could perhaps become a general methodology for facilitating an enzyme to acquire a type enzymatic activity that belongs to another member of the same superfamily, by interrupting a key structural element in order to introduce ambiguity, using site-directed mutagenesis.

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.