Publications

Controlled degradation of proteins is necessary for ensuring their abundance and sustaining a healthy and accurately functioning proteome. One of the degradation routes involves the uncapped 20S proteasome, which cleaves proteins with a partially unfolded region, including those that are damaged or contain intrinsically disordered regions. This degradation route is tightly controlled by a recently discovered family of proteins named Catalytic Core Regulators (CCRs). Here, we show that CCRs function through an allosteric mechanism, coupling the physical binding of the PSMB4 β-subunit with attenuation of the complex's three proteolytic activities. In addition, by dissecting the structural properties that are required for CCR-like function, we could recapitulate this activity using a designed protein that is half the size of natural CCRs. These data uncover an allosteric path that does not involve the proteasome's enzymatic subunits but rather propagates through the non-catalytic subunit PSMB4. This way of 20S proteasome-specific attenuation opens avenues for decoupling the 20S and 26S proteasome degradation pathways as well as for developing selective 20S proteasome inhibitors.

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

Anthropogenic organophosphorus compounds (AOPCs), such as phosphotriesters, are used extensively as plasticizers, flame retardants, nerve agents, and pesticides. To date, only a handful of soil bacteria bearing a phosphotriesterase (PTE), the key enzyme in the AOPC degradation pathway, have been identified. Therefore, the extent to which bacteria are capable of utilizing AOPCs as a phosphorus source, and how widespread this adaptation may be, remains unclear. Marine environments with phosphorus limitation and increasing levels of pollution by AOPCs may drive the emergence of PTE activity. Here, we report the utilization of diverse AOPCs by four model marine bacteria and 17 bacterial isolates from the Mediterranean Sea and the Red Sea. To unravel the details of AOPC utilization, two PTEs from marine bacteria were isolated and characterized, with one of the enzymes belonging to a protein family that, to our knowledge, has never before been associated with PTE activity. When expressed in Escherichia coli with a phosphodiesterase, a PTE isolated from a marine bacterium enabled growth on a pesticide analog as the sole phosphorus source. Utilization of AOPCs may provide bacteria a source of phosphorus in depleted environments and offers a prospect for the bioremediation of a pervasive class of anthropogenic pollutants.

Across the Tree of Life (ToL), the complexity of proteomes varies widely. Our systematic analysis depicts that from the simplest archaea to mammals, the total number of proteins per proteome expanded ∼200-fold. Individual proteins also became larger, and multidomain proteins expanded ∼50-fold. Apart from duplication and divergence of existing proteins, completely new proteins were born. Along the ToL, the number of different folds expanded ∼5-fold and fold combinations ∼20-fold. Proteins prone to misfolding and aggregation, such as repeat and beta-rich proteins, proliferated ∼600-fold and, accordingly, proteins predicted as aggregation-prone became 6-fold more frequent in mammalian compared with bacterial proteomes. To control the quality of these expanding proteomes, core chaperones, ranging from heat shock proteins 20 (HSP20s) that prevent aggregation to HSP60, HSP70, HSP90, and HSP100 acting as adenosine triphosphate (ATP)-fueled unfolding and refolding machines, also evolved. However, these core chaperones were already available in prokaryotes, and they comprise ∼0.3% of all genes from archaea to mammals. This challenge-roughly the same number of core chaperones supporting a massive expansion of proteomes-was met by 1) elevation of messenger RNA (mRNA) and protein abundances of the ancient generalist core chaperones in the cell, and 2) continuous emergence of new substrate-binding and nucleotide-exchange factor cochaperones that function cooperatively with core chaperones as a network.

Production of molecular oxygen was a turning point in the Earths history. The geological record indicates the Great Oxidation Event, which marked a permanent transition to an oxidizing atmosphere around 2.4Ga. However, the degree to which oxygen was available to life before oxygenation of the atmosphere remains unknown. Here, phylogenetic analysis of all known oxygen-utilizing and -producing enzymes (O2-enzymes) indicates that oxygen became widely available to living organisms well before the Great Oxidation Event. About 60% of the O2-enzyme families whose birth can be dated appear to have emerged at the separation of terrestrial and marine bacteria (22families, compared to two families assigned to the last universal common ancestor). This node, dubbed the last universal oxygen ancestor, coincides with a burst of emergence of both oxygenases and other oxidoreductases, thus suggesting a wider availability of oxygen around 3.1Ga.

Natural product biosynthesis (NPB) is the Panda's Thumb of evolutionary biochemistry. Arm races between organisms, and ever-changing environments, result in relentless innovation. This review focusses on enzyme evolution in NPB. First, we review cases of de novo emergence, whereby a completely new enzymatic activity arose in a ligand-binding protein, or a new enzyme emerged including a completely new scaffold. Second, we briefly review the current models for enzyme evolution, and how they explain the inherent promiscuity of NPB enzymes and their tendency to produce multiple related products. We thus suggest that NPB enzymes a priori evolved to generate a specific product; they are, however, trapped in a multifunctional, generalist evolutionary state and thereby produce a diversity of products.

Oligomeric proteins are central to life. Duplication and divergence of their genes is a key evolutionary driver, also because duplications can yield very different outcomes. Given a homomeric ancestor, duplication can yield two paralogs that form two distinct homomeric complexes, or a heteromeric complex comprising both paralogs. Alternatively, one paralog remains a homomer while the other acquires a new partner. However, so far, conflicting trends have been noted with respect to which fate dominates, primarily because different methods and criteria are being used to assign the interaction status of paralogs. Here, we systematically analyzed all Saccharomyces cerevisiae and Escherichia coli oligomeric complexes that include paralogous proteins. We found that the proportions of homo-hetero duplication fates strongly depend on a variety of factors, yet that nonetheless, rigorous filtering gives a consistent picture. In E. coli about 50%, of the paralogous pairs appear to have retained the ancestral homomeric interaction, whereas in S. cerevisiae only ~10% retained a homomeric state. This difference was also observed when unique complexes were counted instead of paralogous gene pairs. We further show that this difference is accounted for by multiple cases of heteromeric yeast complexes that share common ancestry with homomeric bacterial complexes. Our analysis settles contradicting trends and conflicting previous analyses, and provides a systematic and rigorous pipeline for delineating the fate of duplicated oligomers in any organism for which protein-protein interaction data are available.

Traditionally, enzymes are referred to as remarkably fast and specific catalysts. The fact that many enzymes are capable of catalyzing other reactions, besides the one they physiologically specialize in, or evolved for, is definitely not new. Since a long time, breaches of specificity, or promiscuity, of enzymes have been recognized. Early examples of enzyme promiscuity include pyruvate decarboxylase,¹ carbonic anhydrase,² pepsin,³ chymotrypsin,⁴ and L-asparaginase.⁵ However, these early discussions of enzymatic versatility were scarce, and until recently, the wider implications of this darker side of enzymes were largely ignored. During the last decade, protein, and especially enzyme promiscuity, received considerable attention, and their importance in various contexts was systematically studied. Reviews by OBrien and Herschlag,⁶ and later Copley,⁷ were the first to highlight the mechanistic and evolutionary implications of promiscuity. Other, more recent reviews have focused on the practical implications of promiscuity in organic synthesis,^1,8,9 on promiscuity and divergence in certain enzyme families,¹⁰¹³ on mechanistic aspects of promiscuity,^1,14 and on promiscuity in the context of protein evolution¹⁴ and design.¹⁵ The primary focus of this chapter is the role of promiscuity in the evolution of new enzyme functions. New enzymes have constantly emerged throughout the natural history of this planet. Over the past decades, enzymes that degrade synthetic chemicals were introduced to the biosystem,¹⁶²⁰ and enzymes associated with drug resistance,²¹²⁴ provide vivid examples of how rapid the evolution of new enzymatic functions can be. The first direct connection between protein evolution and promiscuity was made in 1976 by Jensen.

Aminoacyl-tRNA synthetases (AARSs) charge tRNA with their cognate amino acids. Many other enzymes use amino acids as substrates, yet discrimination against noncognate amino acids that threaten the accuracy of protein translation is a hallmark of AARSs. Comparing AARSs to these other enzymes allowed us to recognize patterns in molecular recognition and strategies used by evolution for exercising selectivity. Overall, AARSs are 23 orders of magnitude more selective than most other amino acid utilizing enzymes. AARSs also reveal the physicochemical limits of molecular discrimination. For example, amino acids smaller by a single methyl moiety present a discrimination ceiling of ~200, while larger ones can be discriminated by up to 10
⁵-fold. In contrast, substrates larger by a hydroxyl group challenge AARS selectivity, due to promiscuous H-bonding with polar active site groups. This hydroxyl paradox is resolved by editing. Indeed, when the physicochemical discrimination limits are reached, post-transfer editing hydrolysis of tRNAs charged with noncognate amino acids, evolved. The editing site often selectively recognizes the edited noncognate substrate using the very same feature that the synthetic site could not efficiently discriminate against. Finally, the comparison to other enzymes also reveals that the selectivity of AARSs is an explicitly evolved trait, showing some clear examples of how selection acted not only to optimize catalytic efficiency with the target substrate, but also to abolish activity with noncognate threat substrates (negative selection).

The ubiquity of phospho-ligands suggests that phosphate binding emerged at the earliest stage of protein evolution. To evaluate this hypothesis and unravel its details, we identified all phosphate-binding protein lineages in the Evolutionary Classification of Protein Domains database. We found at least 250 independent evolutionary lineages that bind small molecule cofactors and metabolites with phosphate moieties. For many lineages, phosphate binding emerged later as a niche functionality, but for the oldest protein lineages, phosphate binding was the founding function. Across some 4 billion y of protein evolution, side-chain binding, in which the phosphate moiety does not interact with the backbone at all, emerged most frequently. However, in the oldest lineages, and most characteristically in αβα sandwich enzyme domains, N-helix binding sites dominate, where the phosphate moiety sits atop the N terminus of an α-helix. This discrepancy is explained by the observation that N-helix binding is uniquely realized by short, contiguous sequences with reduced amino acid diversity, foremost Gly, Ser, and Thr. The latter two amino acids preferentially interact with both the backbone amide and the side-chain hydroxyl (bidentate interaction) to promote binding by short sequences. We conclude that the first αβα sandwich domains emerged from shorter and simpler polypeptides that bound phospho-ligands via N-helix sites.

Chalcone isomerases are plant enzymes that perform enantioselective oxa-Michael cyclizations of 2'-hydroxychalcones into flavanones. An X-ray crystal structure of an enzyme-product complex combined with molecular dynamics simulations reveal an enzyme mechanism wherein the guanidinium ion of a conserved arginine positions the nucleophilic phenoxide and activates the electrophilic enone for cyclization through Bronsted and Lewis acid interactions. The reaction terminates by asymmetric protonation of the carbanion intermediate syn to the guanidinium. Interestingly, bifunctional guanidine- and urea-based chemical reagents, increasingly used for asymmetric organocatalytic applications, share mechanistic similarities with this natural system. Comparative protein crystal structures and molecular dynamics simulations further demonstrate how two active site water molecules coordinate a hydrogen bond network that enables expanded substrate reactivity for 6'-deoxychalcones in more recently evolved type-2 chalcone isomerases.

How does environmental complexity affect the evolution of single genes? Here, we measured the effects of a set of Bacillus subtilis glutamate dehydrogenase mutants across 19 different environments-from phenotypically homogeneous single-cell populations in liquid media to heterogeneous biofilms, plant roots and soil populations. The effects of individual gene mutations on organismal fitness were highly reproducible in liquid cultures. However, 84% of the tested alleles showed opposing fitness effects under different growth conditions (sign environmental pleiotropy). In colony biofilms and soil samples, different alleles dominated in parallel replica experiments. Accordingly, we found that in these heterogeneous cell populations the fate of mutations was dictated by a combination of selection and drift. The latter relates to programmed prophage excisions that occurred during biofilm development. Overall, for each condition, a wide range of glutamate dehydrogenase mutations persisted and sometimes fixated as a result of the combined action of selection, pleiotropy and chance. However, over longer periods and in multiple environments, nearly all of this diversity would be lost-across all the environments and conditions that we tested, the wild type was the fittest allele.

Protein Science, 27: 12191226. doi:10.1002/pro.2928 Equation 1 has been erroneously printed, and should read as follows: (Formula presented.) whereby f is the fraction of an allele competing against wild-type, f
₀ is the fraction of this allele at generation zero, s is its selection coefficient, and n is the number of generations. Note that if s > 0, as n increases, f/f
₀ increases exponentially; however, given a finite population size, this value will plateau at fixation (i.e., complete takeover of the competing allele). Similarly, if s

Chemists, including the most prominent ones, have always had a profound interest in living systems and biomolecules. Here, I describe few historical endeavors at the interface of Chemistry and Biology that inspired me and shaped my research interests. These examples highlight the uniqueness of implementing a chemist's mindset to fundamental biological questions, by obtaining biological insight via structure-function analyses and by implementing the fundamentals of physical chemistry. I therefore see Chemical Biology as the study of biological systems and molecules with the mindset of a chemist, and in my case, of a physical organic chemist. To me, addressing biological questions with the toolset and mindset of a chemist seems so obvious and natural that there may be no need to define a specialized sub-discipline.

Photorespiration recycles ribulose-1,5-bisphosphate carboxylase/oxygenase ( Rubisco) oxygenation product, 2-phosphoglycolate, back into the Calvin Cycle. Natural photorespiration, however, limits agricultural productivity by dissipating energy and releasing CO2. Several photorespiration bypasses have been previously suggested but were limited to existing enzymes and pathways that release CO2. Here, we harness the power of enzyme and metabolic engineering to establish synthetic routes that bypass photorespiration without CO2 release. By defining specific reaction rules, we systematically identified promising routes that assimilate 2-phosphoglycolate into the Calvin Cycle without carbon loss. We further developed a kinetic-stoichiometric model that indicates that the identified synthetic shunts could potentially enhance carbon fixation rate across the physiological range of irradiation and CO2, even if most of their enzymes operate at a tenth of Rubisco's maximal carboxylation activity. Glycolate reduction to glycolaldehyde is essential for several of the synthetic shunts but is not known to occur naturally. We, therefore, used computational design and directed evolution to establish this activity in two sequential reactions. An acetyl-CoA synthetase was engineered for higher stability and glycolyl-CoA synthesis. A propionyl-CoA reductase was engineered for higher selectivity for glycolyl-CoA and for use of NADPH over NAD(+), thereby favoring reduction over oxidation. The engineered glycolate reduction module was then combined with downstream condensation and assimilation of glycolaldehyde to ribulose 1,5-bisphosphate, thus providing proof of principle for a carbonconserving photorespiration pathway.

Enzymes catalyze a vast range of reactions. Their catalytic performances, mechanisms, global folds, and active-site architectures are also highly diverse, suggesting that enzymes are shaped by an entire range of physiological demands and evolutionary constraints, as well as by chemical and physicochemical constraints. We have attempted to identify signatures of these shaping demands and constraints. To this end, we describe a bird's-eye view of the enzyme space from two angles: evolution and chemistry. We examine various chemical reaction parameters that may have shaped the catalytic performances and active-site architectures of enzymes. We test and weigh these considerations against physiological and evolutionary factors. Although the catalytic properties of the "average" enzyme correlate with cellular metabolic demands and enzyme expression levels, at the level of individual enzymes, a multitude of physiological demands and constraints, combined with the coincidental nature of evolutionary processes, result in a complex picture. Indeed, neither reaction type (a chemical constraint) nor evolutionary origin alone can explain enzyme rates. Nevertheless, chemical constraints are apparent in the convergence of active site architectures in independently evolved enzymes, although significant variations within an architecture are common.

Marine organisms release dimethylsulfide (DMS) via cleavage of dimethylsulfoniopropionate (DMSP). Different genes encoding proteins with DMSP lyase activity are known, yet these exhibit highly variable levels of activity. Most assigned bacterial DMSP lyases, including DddK, DddL, DddQ, DddW, and DddY, appear to belong to one, cupin-like superfamily. Here, we attempted to define and map this superfamily dubbed cupin-DLL (DMSP lyases and lyase-like). To this end, we have pursued the characterization of various recombinant DMSP lyases belonging to this superfamily of metalloenzymes, and especially of DddY and DddL that seem to be the most active DMSP lyases in this superfamily. We identified two conserved sequence motifs that characterize this superfamily. These motifs include the metal-ligating residues that are absolutely essential and other residues including an active site tyrosine that seems to play a relatively minor role in DMSP lysis. We also identified a transition metal chelator, N,N,N',N'-tetrakis(2-pyridylmethyl)ethane-1,2-diamine (TPEN), that selectively inhibits all known members of the cupin-DLL superfamily that exhibit DMSP lyase activity. A phylogenetic analysis indicated that the known DMSP lyase families are sporadically distributed suggesting that DMSP lyases evolved within this superfamily multiple times. However, unusually low specific DMSP lyase activity and genome context analysis suggest that DMSP lyase is not the native function of most cupin-DLL families. Indeed, a systematic profiling of substrate selectivity with a series of DMSP analogues indicated that some members, most distinctly DddY and DddL, are bona fide DMSP lyases, while others, foremost DddQ, may only exhibit promiscuous DMSP lyase activity.

Rational design and directed evolution have proved to be successful approaches to increase catalytic efficiencies of both natural and artificial enzymes. Protein dynamics is recognized as important, but due to the inherent flexibility of biological macromolecules it is often difficult to distinguish which conformational changes are directly related to function. Here, we use directed evolution on an impaired mutant of the proline isomerase CypA and identify two second-shell mutations that partially restore its catalytic activity. We show both kinetically, using NMR spectroscopy, and structurally, by room-temperature X-ray crystallography, how local perturbations propagate through a large allosteric network to facilitate conformational dynamics. The increased catalysis selected for in the evolutionary screen is correlated with an accelerated interconversion between the two catalytically essential conformational substates, which are both captured in the high-resolution X-ray ensembles. Our data provide a glimpse of an evolutionary trajectory and show how subtle changes can fine-tune enzyme function.

The linkage between regulatory elements of transcription, such as promoters, and their protein products is central to gene function. Promoterprotein coevolution is therefore expected, but rarely observed, and the manner by which these two regulatory levels are linked remains largely unknown. We study glutamate dehydrogenasea hub of carbon and nitrogen metabolism. In Bacillus subtilis, two paralogues exist: GudB is constitutively transcribed whereas RocG is tightly regulated. In their active, oligomeric states, both enzymes show similar enzymatic rates. However, swaps of enzymes and promoters cause severe fitness losses, thus indicating promoterenzyme coevolution. Characterization of the proteins shows that, compared to RocG, GudB's enzymatic activity is highly dependent on glutamate and pH. Promoterenzyme swaps therefore result in excessive glutamate degradation when expressing a constitutive enzyme under a constitutive promoter, or insufficient activity when both the enzyme and its promoter are tightly regulated. Coevolution of transcriptional and enzymatic regulation therefore underlies paralogue-specific spatio-temporal control, especially under diverse growth conditions.

Determining the properties of proteins prior to purification saves time and labor. Here, we demonstrate a native mass spectrometry approach for rapid characterization of overexpressed proteins directly in crude cell lysates. The method provides immediate information on the identity, solubility, oligomeric state, overall structure, and stability, as well as ligand binding, without the need for purification.

S-Adenosylmethionine (SAM) is an essential methylation cofactor. The origins of SAM methylation are complex, seemingly demanding the simultaneous emergence of an enzyme that makes SAM and enzyme(s) that utilize it. We report that both ATP and adenosine spontaneously react with methionine to yield SAM, thus suggesting that SAM could have emerged by chance. SAM methylation thus exemplifies how metabolites and pathways can co-emerge through the gradual recruitment of individual enzymes in reverse order.

DNA methyltransferases (MTases) constitute an attractive target for protein engineering, thus opening the road to new ways of manipulating DNA in a unique and selective manner. Here, we review various aspects of MTase engineering, both methodological and conceptual, and also discuss future directions and challenges. Bacterial MTases that are part of restriction/modification (R/M) systems offer a convenient way for the selection of large gene libraries, both in vivo and in vitro. We review these selection methods, their strengths and weaknesses, and also the prospects for new selection approaches that will enable the directed evolution of mammalian DNA methyltransferases (Dnmts). We explore various properties of MTases that may be subject to engineering. These include engineering for higher stability and soluble expression (MTases, including bacterial ones, are prone to misfolding), engineering of the DNA target specificity, and engineering for the usage of S-adenosyl-L-methionine (AdoMet) analogs. Directed evolution of bacterial MTases also offers insights into how these enzymes readily evolve in nature, thus yielding MTases with a huge spectrum of DNA target specificities. Engineering for alternative cofactors, on the other hand, enables modification of DNA with various groups other than methyl and thus can be employed to map and redirect DNA epigenetic modifications.

The nearly 200,000 fatalities following exposure to organophosphorus (OP) pesticides each year and the omnipresent danger of a terroristic attack with OP nerve agents emphasize the demand for the development of effective OP antidotes. Standard treatments for intoxicated patients with a combination of atropine and an oxime are limited in their efficacy. Thus, research focuses on developing catalytic bioscavengers as an alternative approach using OP-hydrolyzing enzymes such as Brevundimonas diminuta phosphotriesterase (PTE). Recently, a PTE mutant dubbed C23 was engineered, exhibiting reversed stereoselectivity and high catalytic efficiency (k_cat/K_M) for the hydrolysis of the toxic enantiomers of VX, CVX, and VR. Additionally, C23s ability to prevent systemic toxicity of VX using a low protein dose has been shown in vivo. In this study, the catalytic efficiencies of V-agent hydrolysis by two newly selected PTE variants were determined. Moreover, in order to establish trends in sequenceactivity relationships along the pathway of PTEs laboratory evolution, we examined k_cat/K_M values of several variants with a number of V-type and G-type nerve agents as well as with different OP pesticides. Although none of the new PTE variants exhibited k_cat/K_M values >10⁷ M⁻¹ min⁻¹ with V-type nerve agents, which is required for effective prophylaxis, they were improved with VR relative to previously evolved variants. The new variants detoxify a broad spectrum of OPs and provide insight into OP hydrolysis and sequenceactivity relationships.

Ancestral reconstruction provides instrumental insights regarding the biochemical and biophysical characteristics of past proteins. A striking observation relates to the remarkably high thermostability of reconstructed ancestors. The latter has been linked to high environmental temperatures in the Precambrian era, the era relating to most reconstructed proteins. We found that inferred ancestors of the serum paraoxonase (PON) enzyme family, including the mammalian ancestor, exhibit dramatically increased thermostabilities compared with the extant, human enzyme (up to 30 °C higher melting temperature). However, the environmental temperature at the time of emergence of mammals is presumed to be similar to the present one. Additionally, the mammalian PON ancestor has superior folding properties (kinetic stability) - unlike the extant mammalian PONs, it expresses in E. coli in a soluble and functional form, and at a high yield. We discuss two potential origins of this unexpectedly high stability. First, ancestral stability may be overestimated by a "consensus effect," whereby replacing amino acids that are rare in contemporary sequences with the amino acid most common in the family increases protein stability. Comparison to other reconstructed ancestors indicates that the consensus effect may bias some but not all reconstructions. Second, we note that high stability may relate to factors other than high environmental temperature such as oxidative stress or high radiation levels. Foremost, intrinsic factors such as high rates of genetic mutations and/or of transcriptional and translational errors, and less efficient protein quality control systems, may underlie the high kinetic and thermodynamic stability of past proteins.

The last decade has seen a growing number of experiments aimed at systematically mapping the effects of mutations in different proteins, and of attempting to correlate their biophysical and biochemical effects with organismal fitness. While insightful, systematic laboratory measurements of fitness effects present challenges and difficulties. Here, we discuss the limitations associated with such measurements, and in particular the challenge of correlating the effects of mutations at the single protein level (\u201cprotein fitness\u201d) with their effects on organismal fitness. A variety of experimental setups are used, with some measuring the direct effects on protein function and others monitoring the growth rate of a model organism carrying the protein mutants. The manners by which fitness effects are calculated and presented also vary, and the conclusions, including the derived distributions of fitness effects of mutations, vary accordingly. The comparison of the effects of mutations in the laboratory to the natural protein diversity, namely to amino acid changes that have fixed in the course of millions of years of evolution, is also debatable. The results of laboratory experiments may, therefore, be less relevant to understanding long-term inter-species variations yet insightful with regard to short-term polymorphism, for example, in the study of the effects of human SNPs.

When asked to edit an issue of Current Opinion in Chemical Biology focused on biocatalysis, we felt that it was high time to cover the esoteric namely, the enzymatic reactions leading to biosynthesis of natural products in ecologically diverse niches and/or in somewhat unusual organisms. High time not because we have exhausted the study of standard enzymes, but because specialized enzymology is revealing new and beautiful aspects of Nature's chemistry, and because several of these niches are actually quite central with respect to life on this planet.

Systematic mappings of the effects of protein mutations are becoming increasingly popular. Unexpectedly, these experiments often find that proteins are tolerant to most amino acid substitutions, including substitutions in positions that are highly conserved in nature. To obtain a more realistic distribution of the effects of protein mutations, we applied a laboratory drift comprising 17 rounds of random mutagenesis and selection of M.HaeIII, a DNA methyltransferase. During this drift, multiple mutations gradually accumulated. Deep sequencing of the drifted gene ensembles allowed determination of the relative effects of all possible single nucleotide mutations. Despite being averaged across many different genetic backgrounds, about 67% of all nonsynonymous, missense mutations were evidently deleterious, and an additional 16% were likely to be deleterious. In the early generations, the frequency of most deleterious mutations remained high. However, by the 17th generation, their frequency was consistently reduced, and those remaining were accepted alongside compensatory mutations. The tolerance to mutations measured in this laboratory drift correlated with sequence exchanges seen in M.HaeIIIs natural orthologs. The biophysical constraints dictating purging in nature and in this laboratory drift also seemed to overlap. Our experiment therefore provides an improved method for measuring the effects of protein mutations that more closely replicates the natural evolutionary forces, and thereby a more realistic view of the mutational space of proteins.

The pioneering model of Henri, Michaelis, and Menten was based on the fast equilibrium assumption: the substrate binds its enzyme reversibly, and substrate dissociation is much faster than product formation. Here, we examine this assumption from a somewhat different point of view, asking what fraction of enzymesubstrate complexes are futile, i.e., result in dissociation rather than product formation. In Knowles notion of a \u201cperfect\u201d enzyme, all encounters of the enzyme with its substrate result in conversion to product. Thus, the perfect enzymes catalytic efficiency, kcat/KM, is constrained by only the diffusion on-rate, and the fraction of futile encounters (defined as φ) approaches zero. The available data on >1000 different enzymes suggest that for ≥90% of enzymes φ > 0.99 and for the \u201caverage enzyme\u201d φ ≥ 0.9999; namely,

Algal blooms produce large amounts of dimethyl sulfide (DMS), a volatile with a diverse signaling role in marine food webs that is emitted to the atmosphere, where it can affect cloud formation. The algal enzymes responsible for forming DMS from dimethylsulfoniopropionate (DMSP) remain unidentified despite their critical role in the global sulfur cycle. We identified and characterized Alma1, a DMSP lyase from the bloom-forming algae Emiliania huxleyi. Alma1 is a tetrameric, redox-sensitive enzyme of the aspartate racemase superfamily. Recombinant Alma1 exhibits biochemical features identical to the DMSP lyase in E. huxleyi, and DMS released by various E. huxleyi isolates correlates with their Alma1 levels. Sequence homology searches suggest that Alma1 represents a gene family present in major, globally distributed phytoplankton taxa and in other marine organisms.

The highly toxic organophosphorus (OP) nerve agent VX is characterized by a remarkable biological persistence which limits the effectiveness of standard treatment with atropine and oximes. Existing OP hydrolyzing enzymes show low activity against VX and hydrolyze preferentially the less toxic P(+)-VX enantiomer. Recently, a phosphotriesterase (PTE) mutant, C23, was engineered towards the hydrolysis of the toxic P(-) isomers of VX and other V-type agents with relatively high in vitro catalytic efficiency (k(cat)/K-M = 5 x 10(6) M-1 min(-1)). To investigate the suitability of the PTE mutant C23 as a catalytic scavenger, an in vivo guinea pig model was established to determine the efficacy of post-exposure treatment with C23 alone against VX intoxication. Injection of C23 (5 mg kg(-1) i.v.) 5 min after s.c. challenge with VX (similar to 2LD(50)) prevented systemic toxicity. A lower C23 dose (2 mg kg(-1)) reduced systemic toxicity and prevented mortality. Delayed treatment (i.e., 15 min post VX) with 5 mg kg(-1) C23 resulted in survival of all animals and only in moderate systemic toxicity. Although C23 did not prevent inhibition of erythrocyte acetylcholinesterase (AChE) activity, it partially preserved brain AChE activity. C23 therapy resulted in a rapid decrease of racemic VX blood concentration which was mainly due to the rate of degradation of the toxic P(-)-VX enantiomer that correlates with the C23 blood levels and its k(cat)/K-M value. Although performed under anesthesia, this proof-of-concept study demonstrated for the first time the ability of a catalytic bioscavenger to prevent systemic VX toxicity when given alone as a single postexposure treatment, and enables an initial assessment of a time window for this approach. In conclusion, the PTE mutant C23 may be considered as a promising starting point for the development of highly effective catalytic bioscavengers for post-exposure treatment of V-agents intoxication. (C) 2014 Elsevier Ireland Ltd. All rights

I discuss some physico-chemical and evolutionary aspects of enzyme accuracy (selectivity, specificity) and speed (turnover rate, processivity). Accuracy can be a beneficial side-product of active-sites being refined to proficiently convert a given substrate into one product. However, exclusion of undesirable, non-cognate substrates is also an explicitly evolved trait that may come with a cost. I define two schematic mechanisms. Ground-state discrimination applies to enzymes where selectivity is achieved primarily at the level of substrate binding. Exemplified by DNA methyltransferases and the ribosome, ground-state discrimination imposes strong accuracy-rate tradeoffs. Alternatively, transition-state discrimination, applies to relatively small substrates where substrate binding and chemistry are efficiently coupled, and evokes weaker tradeoffs. Overall, the mechanistic, structural and evolutionary basis of enzymatic accuracy-rate tradeoffs merits deeper understanding.

Assignment of protein folds to functions indicates that >60% of folds carry out one or two enzymatic functions, while few folds, for example, the TIM-barrel and Rossmann folds, exhibit hundreds. Are there structural features that make a fold amenable to functional innovation (innovability)? Do these features relate to robustness-the ability to readily accumulate sequence changes? We discuss several hypotheses regarding the relationship between the architecture of a protein and its evolutionary potential. We describe how, in a seemingly paradoxical manner, opposite properties, such as high stability and rigidity versus conformational plasticity and structural order versus disorder, promote robustness and/or innovability. We hypothesize that polarity. - differentiation and low connectivity between a protein's scaffold and its active-site. -is a key prerequisite for innovability.

Protein engineering by directed evolution relies on the use of libraries enriched with beneficial variants. Such libraries should explore large mutational diversities while avoiding high loads of deleterious mutations. Here we describe a simple protocol for incorporating synthetic oligonucleotides that encode designed, site-specific mutations by assembly PCR. This protocol enables a researcher to "hedge the bets," namely, to explore a large number of potentially beneficial mutations in a combinatorial manner such that individual library variants carry a limited number of mutations.

Short insertions and deletions (InDels) comprise an important part of the natural mutational repertoire. InDels are, however, highly deleterious, primarily because two-thirds result in frame-shifts. Bypass through slippage over homonucleotide repeats by transcriptional and/or translational infidelity is known to occur sporadically. However, the overall frequency of bypass and its relation to sequence composition remain unclear. Intriguingly, the occurrence of InDels and the bypass of frame-shifts are mechanistically related - occurring through slippage over repeats by DNA or RNA polymerases, or by the ribosome, respectively. Here, we show that the frequency of frame-shifting InDels, and the frequency by which they are bypassed to give full-length, functional proteins, are indeed highly correlated. Using a laboratory genetic drift, we have exhaustively mapped all InDels that occurred within a single gene. We thus compared the naive InDel repertoire that results from DNA polymerase slippage to the frame-shifting InDels tolerated following selection to maintain protein function. We found that InDels repeatedly occurred, and were bypassed, within homonucleotide repeats of 3-8 bases. The longer the repeat, the higher was the frequency of InDels formation, and the more frequent was their bypass. Besides an expected 8A repeat, other types of repeats, including short ones, and G and C repeats, were bypassed. Although obtained in vitro, our results indicate a direct link between the genetic occurrence of InDels and their phenotypic rescue, thus suggesting a potential role for frame-shifting InDels as bridging evolutionary intermediates.

Alignments of orthologous protein sequences convey a complex picture. Some positions are utterly conserved whilst others have diverged to variable degrees. Amongst the latter, many are non-exchangeable between extant sequences. How do functionally critical and highly conserved residues diverge? Why and how did these exchanges become incompatible within contemporary sequences? Our model is phosphoglycerate kinase (PGK), where lysine 219 is an essential active-site residue completely conserved throughout Eukaryota and Bacteria, and serine is found only in archaeal PGKs. Contemporary sequences tested exhibited complete loss of function upon exchanges at 219. However, a directed evolution experiment revealed that two mutations were sufficient for human PGK to become functional with serine at position 219. These two mutations made position 219 permissive not only for serine and lysine, but also to a range of other amino acids seen in archaeal PGKs. The identified trajectories that enabled exchanges at 219 show marked sign epistasis - a relatively small loss of function with respect to one amino acid (lysine) versus a large gain with another (serine, and other amino acids). Our findings support the view that, as theoretically described, the trajectories underlining the divergence of critical positions are dominated by sign epistatic interactions. Such trajectories are an outcome of rare mutational combinations. Nonetheless, as suggested by the laboratory enabled K219S exchange, given enough time and variability in selection levels, even utterly conserved and functionally essential residues may change.

Backbone modifications via insertions and deletions (InDels) may exert dramatic effects, for better (mediating new functions) and for worse (causing loss of structure and/or function). However, contrary to point mutations (substitutions), our knowledge of the evolution and structural-functional effects of InDels is limited and so is our capability to engineer them. We sought to assess how deleterious InDels are relative to point mutations and understand the mechanisms that mediate their acceptance. Analysis of the evolution of InDels in orthologous protein phylogenies indicated that their rate of purging is 9- to 100-fold higher than for point mutations. In yeast, for example, the substitutions-to-InDels ratio is approximately 14-fold higher in protein coding than in noncoding regions. The incorporation of InDels relative to substitutions is not only slow but also nonlinear. On average, ≥50 substitutions accumulate before the appearance of the first InDel. We also found enriched substitutions in sequential and spatial proximity to InDels, suggesting that certain substitutions are correlated with InDels. As indicated by the lag in InDels accumulation, some of these correlated substitutions may have occurred first, as apparently neutral mutations, and later enabled the accumulation of InDels that would be otherwise purged. Thus, compensatory substitutions may follow InDels in an "adaptive walk" as traditionally assumed, but might also accumulate first, by "neutral roaming." The dynamics of InDels accumulation also depends on their genomic frequencies - InDels in flies are 4-fold more frequent than in yeast and tend to be compensated rather than enabled.

Although largely deemed as structurally conserved, catalytic metal ion sites can rearrange, thereby contributing to enzyme evolvability. Here, we show that in paraoxonase-1, a lipo-lactonase, catalytic promiscuity and divergence into an organophosphate hydrolase are correlated with an alternative mode of the catalytic Ca^{2 +}. We describe the crystal structures of active-site mutants bearing mutations at position 115. The histidine at this position acts as a base to activate the lactone-hydrolyzing water molecule. Mutations to Trp or Gln indeed diminish paraoxonase-1's lactonase activity; however, the promiscuous organophosphate hydrolase activity is enhanced. The structures reveal a 1.8-Å upward displacement towards the enzyme's surface of the catalytic Ca^{2 +} in the His115 mutants and configurational changes in the ligating side chains and water molecules, relative to the wild-type enzyme. Biochemical analysis and molecular dynamics simulations suggest that this alternative, upward metal mode mediates the promiscuous hydrolysis of organophosphates. The upward Ca^{2 +} mode observed in the His115 mutants also appears to mediate the wild type's paraoxonase activity. However, whereas the upward mode dominates in the Trp115 mutant, it is scarcely populated in wild type. Thus, the plasticity of active-site metal ions may permit alternative, latent, promiscuous activities and also provide the basis for the divergence of new enzymatic functions.

Gene expression depends on the frequency of transcription events (burst frequency) and on the number of mRNA molecules made per event (burst size). Both processes are encoded in promoter sequence, yet their dependence on mutations is poorly understood. Theory suggests that burst size and frequency can be distinguished by monitoring the stochastic variation (noise) in gene expression: Increasing burst size will increase mean expression without changing noise, while increasing burst frequency will increase mean expression and decrease noise. To reveal principles by which promoter sequence regulates burst size and frequency, we randomly mutated 22 yeast promoters chosen to span a range of expression and noise levels, generating libraries of hundreds of sequence variants. In each library, mean expression (m) and noise (coefficient of variation, η) varied together, defining a scaling curve: η² = b/m + η_ext². This relation is expected if sequence mutations modulate burst frequency primarily. The estimated burst size (b) differed between promoters, being higher in promoter containing a TATA box and lacking a nucleosome-free region. The rare variants that significantly decreased b were explained by mutations in TATA, or by an insertion of an out-of-frame translation start site. The decrease in burst size due to mutations in TATA was promoter-dependent, but independent of other mutations. These TATA box mutations also modulated the responsiveness of gene expression to changing conditions. Our results suggest that burst size is a promoter-specific property that is relatively robust to sequence mutations but is strongly dependent on the interaction between the TATA box and promoter nucleosomes.

In nature, the evolution of new protein functions is driven not only by side-chain substitutions (point mutations), but also by backbone modifications (insertions and deletions). The current laboratory diversification methods, however, are largely limited to point mutations. Of particular interest are short insertions-by-duplication that are frequent in nature but cannot be introduced in vitro in a library format (i.e. in random locations and lengths). Here, we describe a new procedure that allows the generation of tandem repeats of random fragments of the target gene via rolling-circle amplification, and the concurrent incorporation of these repeats into the target gene. This procedure, dubbed tandem repeat insertion, or TRINS, results in a library of genes carrying insertions-by-duplication of variable lengths (3-150 bp) at random positions. This diversification pattern allows sampling of sequence space regions that are not readily accessible by other protocols. We demonstrate this method by constructing three different gene libraries, and by selecting insertion variants of TEM-1 β-lactamase.

Only decades after the introduction of organophosphate pesticides, bacterial phosphotriesterases (PTEs) have evolved to catalyze their degradation with remarkable efficiency. Their closest known relatives, lactonases, with promiscuous phosphotriasterase activity, dubbed PTE-like lactonases (PLLs), share only 30% sequence identity and also differ in the configuration of their active-site loops. PTE was therefore presumed to have evolved from a yet unknown PLL whose primary activity was the hydrolysis of quorum sensing homoserine lactones (HSLs) (Afriat et al. (2006) Biochemistry45, 13677-13686). However, how PTEs diverged from this presumed PLL remains a mystery. In this study we investigated loop remodeling as a means of reconstructing a homoserine lactonase ancestor that relates to PTE by few mutational steps. Although, in nature, loop remodeling is a common mechanism of divergence of enzymatic functions, reproducing this process in the laboratory is a challenge. Structural and phylogenetic analyses enabled us to remodel one of PTE's active-site loops into a PLL-like configuration. A deletion in loop 7, combined with an adjacent, highly epistatic, point mutation led to the emergence of an HSLase activity that is undetectable in PTE (k_cat/K_M values of up to 2 × 10⁴). The appearance of the HSLase activity was accompanied by only a minor decrease in PTE's paraoxonase activity. This specificity change demonstrates the potential role of bifunctional intermediates in the divergence of new enzymatic functions and highlights the critical contribution of loop remodeling to the rapid divergence of new enzyme functions.

The origins of enzyme specificity are well established. However, the molecular details underlying the ability of a single active site to promiscuously bind different substrates and catalyze different reactions remain largely unknown. To better understand the molecular basis of enzyme promiscuity, we studied the mammalian serum paraoxonase 1 (PON1) whose native substrates are lipophilic lactones. We describe the crystal structures of PON1 at a catalytically relevant pH and of its complex with a lactone analogue. The various PON1 structures and the analysis of active-site mutants guided the generation of docking models of the various substrates and their reaction intermediates. The models suggest that promiscuity is driven by coincidental overlaps between the reactive intermediate for the native lactonase reaction and the ground and/or intermediate states of the promiscuous reactions. This overlap is also enabled by different active-site conformations: the lactonase activity utilizes one active-site conformation whereas the promiscuous phosphotriesterase activity utilizes another. The hydrolysis of phosphotriesters, and of the aromatic lactone dihydrocoumarin, is also driven by an alternative catalytic mode that uses only a subset of the active-site residues utilized for lactone hydrolysis. Indeed, PON1's active site shows a remarkable level of networking and versatility whereby multiple residues share the same task and individual active-site residues perform multiple tasks (e.g., binding the catalytic calcium and activating the hydrolytic water). Overall, the coexistence of multiple conformations and alternative catalytic modes within the same active site underlines PON1's promiscuity and evolutionary potential.

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⁴ 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.

An ex vivo protocol was developed to assay the antidotal capacity of rePON1 variants to protect endogenous acetylcholinesterase and butyrylcholinesterase in human whole blood against OP nerve agents. This protocol permitted us to address the relationship between blood rePON1 concentrations, their kinetic parameters, and the level of protection conferred by rePON1 on the cholinesterases in human blood, following a challenge with cyclosarin (GF). The experimental data thus obtained were in good agreement with the predicted percent residual activities of blood cholinesterases calculated on the basis of the rate constants for inhibition of human acetylcholinesterase and butyrylcholinesterase by GF, the concentration of the particular rePON1 variant, and its k_cat/K_m value for GF. This protocol thus provides a rapid and reliable ex vivo screening tool for identification of rePON1 bioscavenger candidates suitable for protection of humans against organophosphorus-based toxicants. The results also permitted the refinement of a mathematical model for estimating the efficacious dose of rePON1s variants required for prophylaxis in humans.

Why do certain proteins evolve much slower than others? We compared not only rates per protein, but also rates per position within individual proteins. For ∼90% of proteins, the distribution of positional rates exhibits three peaks: a peak of slow evolving residues, with average log ₂[normalized rate], log ₂μ, of ca. -2, corresponding primarily to core residues; a peak of fast evolving residues (log ₂μ ∼ 0.5) largely corresponding to surface residues; and a very fast peak (log ₂μ ∼ 2) associated with disordered segments. However, a unique fraction of proteins that evolve very slowly exhibit not only a negligible fast peak, but also a peak with a log ₂μ ∼ -4, rather than the standard core peak of -2. Thus, a "freeze" of a protein's surface seems to stop core evolution as well. We also observed a much higher fraction of substitutions in potentially interacting residues than expected by chance, including substitutions in pairs of contacting surface-core residues. Overall, the data suggest that accumulation of surface substitutions enables the acceptance of substitutions in core positions. The underlying reason for slow evolution might therefore be a highly constrained surface due to protein - protein interactions or the need to prevent misfolding or aggregation. If the surface is inaccessible to substitutions, so becomes the core, thus resulting in very slow overall rates.

Two different scenarios for the recruitment of evolutionary starting points and their subsequent divergence to give new enzymes have been described. The coincidental, promiscuous starting activity may regard the same reaction chemistry on a new substrate (substrate ambiguity). Alternatively, substrate binding guides the recruitment of an enzyme whose reaction chemistry differs from that of the newly evolving one (catalytic promiscuity). While substrate ambiguity seems to underlie the divergence of most enzyme families, the relative levels of occurrence of these scenarios remain unknown. Screening the Escherichia coli proteome with a comparative series of xenobiotic substrates, we found that substrate ambiguity was, as anticipated, more frequent than reaction promiscuity. However, for at least one unnatural reaction (phosphonoesterase), a promiscuous enzyme was identified only when the substrate was decorated with the naturally abundant phosphate group. These findings support the prevailing hypothesis of chemistry-driven divergence but also suggest that recognition of familiar substrate motifs plays a role. In the absence of enzymes catalyzing the same chemistry, having a familiar, naturally occurring substrate motif (chemophore) such as phosphate may increase the likelihood of catalytic promiscuity. Chemophore anchoring may also find practical applications in identifying catalysts for unnatural reactions.

A newly identified bacterial strain that can grow in the presence of arsenate and possibly in the absence of phosphate, has raised much interest, but also fueled an active debate. Can arsenate substitute for phosphate in some or possibly in most of the absolutely essential phosphate-based biomolecules, including DNA? If so, then the possibility of alternative, arsenic-based life forms must be considered. The physicochemical similarity of these two oxyanions speaks in favor of this idea. However, arsenate-esters and arsenate-diesters in particular are extremely unstable in aqueous media. Here, we explore the potential of arsenate to be used as substrate by phosphate-utilizing enzymes. We review the existing literature on arsenate enzymology, that intriguingly, dates back to the 1930s. We address the issue of how and to what degree proteins can distinguish between arsenate and phosphate and what is known in general about oxyanion specificity. We also discuss how phosphate-arsenate promiscuity may affect evolutionary transitions between phosphate- and arsenate-based biochemistry. Finally, we highlight potential applications of arsenate as a structural and mechanistic probe of enzymes whose catalyzed reactions involve the making or breaking of phosphoester bonds.

Fluorogenic organophosphate inhibitors of acetylcholinesterase (AChE) homologous in structure to nerve agents provide useful probes for high throughput screening of mammalian paraoxonase (PON1) libraries generated by directed evolution of an engineered PON1 variant with wild-type like specificity (rePON1). Wt PON1 and rePON1 hydrolyze preferentially the less-toxic R_P enantiomers of nerve agents and of their fluorogenic surrogates containing the fluorescent leaving group, 3-cyano-7-hydroxy-4-methylcoumarin (CHMC). To increase the sensitivity and reliability of the screening protocol so as to directly select rePON1 clones displaying stereo-preference towards the toxic S_P enantiomer, and to determine accurately K_m and k_cat values for the individual isomers, two approaches were used to obtain the corresponding S_P and R_P isomers: (a) stereo-specific synthesis of the O-ethyl, O-n-propyl, and O-i-propyl analogs and (b) enzymic resolution of a racemic mixture of O-cyclohexyl methylphosphonylated CHMC. The configurational assignments of the S_P and R_P isomers, as well as their optical purity, were established by X-ray diffraction, reaction with sodium fluoride, hydrolysis by selected rePON1 variants, and inhibition of AChE. The S_P configuration of the tested surrogates was established for the enantiomer with the more potent anti-AChE activity, with S_P/R_P inhibition ratios of 10-100, whereas the R_P isomers of the O-ethyl and O-n-propyl were hydrolyzed by wt rePON1 about 600- and 70-fold faster, respectively, than the S_P counterpart. Wt rePON1-induced R_P/S_P hydrolysis ratios for the O-cyclohexyl and O-i-propyl analogs are estimated to be 1000. The various S_P enantiomers of O-alkyl-methylphosphonyl esters of CHMC provide suitable ligands for screening rePON1 libraries, and can expedite identification of variants with enhanced catalytic proficiency towards the toxic nerve agents.

The primary sequence of proteins usually dictates a single tertiary and quaternary structure. However, certain proteins undergo reversible backbone rearrangements. Such metamorphic proteins provide a means of facilitating the evolution of new folds and architectures. However, because natural folds emerged at the early stages of evolution, the potential role of metamorphic intermediates in mediating evolutionary transitions of structure remains largely unexplored. We evolved a set of new proteins based on ∼100 amino acid fragments derived from tachylectin-2 - a monomeric, 236 amino acids, five-bladed β-propeller. Their structures reveal a unique pentameric assembly and novel β-propeller structures. Although identical in sequence, the oligomeric subunits adopt two, or even three, different structures that together enable the pentameric assembly of two propellers connected via a small linker. Most of the subunits adopt a wild-type-like structure within individual five-bladed propellers. However, the bridging subunits exhibit domain swaps and asymmetric strand exchanges that allow them to complete the two propellers and connect them. Thus, the modular and metamorphic nature of these subunits enabled dramatic changes in tertiary and quaternary structure, while maintaining the lectin function. These oligomers therefore comprise putative intermediates via which β-propellers can evolve from smaller elements. Our data also suggest that the ability of one sequence to equilibrate between different structures can be evolutionary optimized, thus facilitating the emergence of new structures.

Understanding enzyme catalysis through the analysis of natural enzymes is a daunting challenge-their active sites are complex and combine numerous interactions and catalytic forces that are finely coordinated. Study of more rudimentary (wo)man-made enzymes provides a unique opportunity for better understanding of enzymatic catalysis. KE07, a computationally designed Kemp eliminase that employs a glutamate side chain as the catalytic base for the critical proton abstraction step and an apolar binding site to guide substrate binding, was optimized by seven rounds of random mutagenesis and selection, resulting in a >200-fold increase in catalytic efficiency. Here, we describe the directed evolution process in detail and the biophysical and crystallographic studies of the designed KE07 and its evolved variants. The optimization of KE07's activity to give a k_cat/K_M value of ~2600 s^-1 M^-1 and an ~10⁶-fold rate acceleration (k_cat/k_uncat) involved the incorporation of up to eight mutations. These mutations led to a marked decrease in the overall thermodynamic stability of the evolved KE07s and in the configurational stability of their active sites. We identified two primary contributions of the mutations to KE07's improved activity: (i) the introduction of new salt bridges to correct a mistake in the original design that placed a lysine for leaving-group protonation without consideration of its "quenching" interactions with the catalytic glutamate, and (ii) the tuning of the environment, the pK_a of the catalytic base, and its interactions with the substrate through the evolution of a network of hydrogen bonds consisting of several charged residues surrounding the active site.

How do intricate multi-residue features such as protein-protein interfaces evolve? To address this question, we evolved a new colicin-immunity binding interaction. We started with Im9, which inhibits its cognate DNase ColE9 at 10(-14) M affinity, and evolved it toward ColE7, which it inhibits promiscuously (K(d) > 10(-8) M). Iterative rounds of random mutagenesis and selection toward higher affinity for ColE7, and selectivity ( against ColE9 inhibition), led to an similar to 10(5)-fold increase in affinity and a 10(8)-fold increase in selectivity. Analysis of intermediates along the evolved variants revealed that changes in the binding configuration of the Im protein uncovered a latent set of interactions, thus providing the key to the rapid divergence of new Im7 variants. Overall, protein-protein interfaces seem to share the evolvability features of enzymes, that is, the exploitation of promiscuous interactions and alternative binding configurations via 'generalist' intermediates, and the key role of compensatory stabilizing mutations in facilitating the divergence of new functions.

Most protein mutations, and mutations that alter protein functions in particular, undermine stability and are therefore deleterious. Chaperones, or heat-shock proteins, are often implicated in buffering mutations, and could thus facilitate the acquisition of neutral genetic diversity and the rate of adaptation. We examined the ability of the Escherichia coli GroEL/GroES chaperonins to buffer destabilizing and adaptive mutations. Here we show that mutational drifts performed in vitro with four different enzymes indicated that GroEL/GroES overexpression doubled the number of accumulating mutations, and promoted the folding of enzyme variants carrying mutations in the protein core and/or mutations with higher destabilizing effects (destabilization energies of >3.5kcalmol^-1, on average, versus 1kcalmol ^-1 in the absence of GroEL/GroES). The divergence of modified enzymatic specificity occurred much faster under GroEL/GroES overexpression, in terms of the number of adapted variants (2-fold) and their improved specificity and activity (10-fold). These results indicate that protein stability is a major constraint in protein evolution, and buffering mechanisms such as chaperonins are key in alleviating this constraint.

The traditional view that proteins possess absolute functional specificity and a single, fixed structure conflicts with their marked ability to adapt and evolve new functions and structures. We consider an alternative,"avant- garde view" in which proteins are conformational dynamic and exhibit functional promiscuity. We surmise that these properties are the foundation stones of protein evolvability; they facilitate the divergence of new functions within existing folds and the evolution of entirely new folds. Packing modes of proteins also affect their evolvability, and poorly packed, disordered, and conformationally diverse proteins may exhibit high evolvability. This dynamic view of protein structure, function, and evolvability is extrapolated to describe hypothetical scenarios for the evolution of the early proteins and future research directions in the area of protein dynamism and evolution.

What changes occur when a natural protein that had been under low mutation rates is subjected to a neutral drift at high mutational loads, thus generating genetically diverse (polymorphic) gene ensembles that all maintain the protein's original function and structure? To address this question we subjected large populations of TEM-1 β-lactamase to a prolonged neutral drift, applying high mutation rates and purifying selection to maintain TEM-1's existing penicillinase activity. Purging of deleterious mutations and enrichment of beneficial ones maintained the sequence of these ensembles closer to TEM-1's family consensus and inferred ancestor. In particular, back-to-consensus/ancestor mutations that increase TEM-1's kinetic and thermodynamic stability were enriched. These acted as global suppressors and enabled the tolerance of a broad range of deleterious mutations, thus further increasing the genetic diversity of the drifting populations. The probability of a new function emerging (cefotaxime degradation) was also substantially increased in these ensembles owing to the presence of many gene variants carrying the global suppressors. Our findings indicate the unique features of large, polymorphic neutral ensembles generated under high mutational loads and prompt the speculation that the progenitors of today's proteins may have evolved under high mutational loads. The results also suggest that predictable back-to-consensus/ancestor changes can be used in the laboratory to generate highly diverse and evolvable gene libraries.

We address recent developments in the area of laboratory, or directed evolution, with a focus on enzymes and on new methodologies of generic potential. We survey three main areas: (i) library making techniques, including the application of computational and rational methods for library design; (ii) screening and selection techniques, including recent applications of enzyme screening by FACS (fluorescence activated cell sorter); (iii) new approaches for performing directed evolution, and in particular, the application of 'neutral drifts' (libraries generated by rounds of mutation and selection for the enzyme's original function) and of consensus mutations to generate highly evolvable starting points for directed evolution.

Serum paraoxonase (PON1) is a lipolactonase that associates with HDL-apolipoprotein A-I (HDL-apoA-I) and thereby plays a role in the prevention of atherosclerosis. Current sera tests make use of promiscuous substrates and provide no indications regarding HDL-PON1 complex formation. We developed new enzymatic tests that detect total PON1 levels, irrespective of HDL status and R/Q polymorphism, as well as the degree of catalytic stimulation and increased stability that follow PON1's tight binding to HDL-apoA-I. The tests are based on measuring total PON1 levels with a fluorogenic phosphotriester, measuring the lipolactonase activity with a chromogenic lactone, and assaying the enzyme's chelator-mediated inactivation rate. The latter two are affected by tight HDL binding and thereby derive the levels of the serum PON1-HDL complex. We demonstrate these new tests with a group of healthy individuals (n = 54) and show that the levels of PON1-HDL vary by a factor of 12. Whereas the traditionally applied paraoxonase and arylesterase tests weakly reflect PON1-HDL levels (R = 0.64), the lipolactonase test provides better correlation (R = 0.80). These new tests indicate the levels and activity of PON1 in a physiologically relevant context aswell as the levels and quality of the HDL particles with which the enzyme is associated.

How the thermodynamic stability effects of protein mutations (ΔΔG) are distributed is a fundamental property related to the architecture, tolerance to mutations (mutational robustness), and evolutionary history of proteins. The stability effects of mutations also dictate the rate and dynamics of protein evolution, with deleterious mutations being the main inhibitory factor. Using the FoldX algorithm that attempts to computationally predict ΔΔG effects of mutations, we deduced the overall distributions of stability effects for all possible mutations in 21 different globular, single domain proteins. We found that these distributions are strikingly similar despite a range of sizes and folds, and largely follow a bi-Gaussian function: The surface residues exhibit a narrow distribution with a mildly destabilizing mean ΔΔG (∼ 0.6 kcal/mol), whereas the core residues exhibit a wider distribution with a stronger destabilizing mean (∼ 1.4 kcal/mol). Since smaller proteins have a higher fraction of surface residues, the relative weight of these single distributions correlates with size. We also found that proteins evolved in the laboratory follow an essentially identical distribution, whereas de novo designed folds show markedly less destabilizing distributions (i.e. they seem more robust to the effects of mutations). This bi-Gaussian model provides an analytical description of the predicted distributions of mutational stability effects. It comprises a novel tool for analyzing proteins and protein models, for simulating the effect of mutations under evolutionary processes, and a quantitative description of mutational robustness.

Biological systems exhibit mutational robustness, or neutrality, whereby the impact of mutations is minimized. Does neutrality hamper their ability to adapt in the face of changing environments? We monitored changes in genotype and phenotype that occur within a neutral mutational network of an enzyme, experimentally and computationally, see accompanying article.. Using the enzyme PON1 as a model, we performed random mutagenesis and purifying selection to purge deleterious mutations. We characterized similar to 300 variants that are apparently neutral, or close to neutral, with respect to PON1's levels of expression and native lactonase activity. Their activities with promiscuous substrates and ligands indicated significant changes in adaptive potentials. Almost half of the variants exhibited changes in promiscuous activities, specificities, or inhibition, and several of these were found to be one or two mutations, closer to potentially new phenotypes: aryl esterase, thiolactonase, phosphotriesterase, or drug resistance. This empirical measure of phenotypic changes under neutrality supports the notion that sequence changes that are neutral, i.e., non-adaptive, in a current context can facilitate adaptation under changing circumstances, by both expanding the activity range of existing enzymes and thus providing an immediate advantage, and by reducing the number of mutations required for divergence of new functions.

In essence, evolutionary processes occur gradually, while maintaining fitness throughout. Along this line, it has been proposed that the ability of a progenitor to promiscuously catalyze a low level of the evolving activity could facilitate the divergence of a new function by providing an immediate selective advantage. To directly establish a role for promiscuity in the divergence of natural enzymes, we attempted to trace the origins of a bacterial phosphotriesterase (PTE), an enzyme thought to have evolved for the purpose of degradation of a synthetic insecticide introduced in the 20th century. We surmised that PTE's promiscuous lactonase activity may be a vestige of its progenitor and tested homologues annotated as "putative PTEs" for lactonase and phosphotriesterase activity. We identified three genes that define a new group of microbial lactonases dubbed PTE-like lactonases (PLLs). These enzymes proficiently hydrolyze various lactones, and in particular quorum-sensing N-acyl homoserine lactones (AHLs), and exhibit much lower promiscuous phosphotriesterase activities. PLLs share key sequence and active site features with PTE and differ primarily by an insertion in one surface loop. Given their biochemical and biological function, PLLs are likely to have existed for many millions of years. PTE could have therefore evolved from a member of the PLL family while utilizing its latent promiscuous paraoxonase activity as an essential starting point.

The past few years have seen significant advances in research related to the 'latent skills' of enzymes - namely, their capacity to promiscuously catalyze reactions other than the ones they evolved for. These advances regard (i) the mechanism of catalytic promiscuity - how enzymes, that generally exert exquisite specificity, promiscuously catalyze other, and sometimes barely related, reactions; (ii) the evolvability of promiscuous functions - namely, how latent activities evolve further, and in particular, how promiscuous activities can firstly evolve without severely compromising the original activity. These findings have interesting implications on our understanding of how new enzymes evolve. They support the key role of catalytic promiscuity in the natural history of enzymes, and suggest that today's enzymes diverged from ancestral proteins catalyzing a whole range of activities at low levels, to create families and superfamilies of potent and highly specialized enzymes.

Keywords: Toxicology

The efficient amplification of genomic libraries, cDNA libraries and other complex mixtures of genes by PCR is impeded by two phenomena: firstly, short fragments tend to be amplified in preference to larger ones; and, secondly, artifactual fragments are generated by recombination between homologous regions of DNA. Recombination in this case occurs when a primer is partially extended on one template during one cycle of PCR and further extended on another template during a later cycle. Thus, chimeric molecules are generated, the short ones of which are then preferentially amplified as described in Figure 1. A variety of PCR protocols have been proposed to minimize these problems, most of which rely on high template concentrations and low numbers of PCR cycles. Clearly, however, such an approach is not viable if little template DNA is available. Here we describe a protocol for amplifying complex DNA mixtures, based on the compartmentalization of genes in a water-in-oil (w/o) emulsion. Template fragments are segregated in the minute aqueous droplets of the emulsion and amplified by PCR in isolation (Fig. 1). This approach alleviates the problems described above while enabling the use of small amounts of template DNA and high numbers of PCR cycles. Box 1 describes an alternative method for generating very stable emulsions for emulsion PCR using the surfactant ABIL EM 90 (Fig. 2).

Homocysteine thiolactone (tHcy) is deemed a risk factor for cardiovascular diseases and strokes, presumably because it acylates the side chain of protein lysine residues ("N-homocysteinylation"), thereby causing protein damage and autoimmune responses. We analysed the kinetics of hydrolysis and aminolysis of tHcy and two related thiolactones (γ-thiobutyrolactone and N-trimethyl-tHcy), and we have thereby described the first detailed mechanism of thiolactone aminolysis. As opposed to the previously studied (thio and oxo)esters and (oxo)lactones, aminolysis of thiolactones was found to be first order with respect to amine concentration. Anchimeric assistance by the α-amino group of tHcy (through general acid/base catalysis) could not be detected, and the Brønsted plot (nucleophilicity versus pK_a) for aminolysis yielded a slope (β^nuc) value of 0.66. These data support a mechanism of aminolysis where the rate-determining step is the formation of a zwitterionic tetrahedral intermediate. The β^nuc value and steric factors dictate a regime whereby, at physiological pH values (pH7.4), maximal reactivity of tHcy is exhibited with primary amine groups with a pK_α value of 7.7; this allows the reactivity of various protein amino groups towards N-homocysteinylation to be predicted.

Serum paraoxonases (PONs) are calcium-dependent lactonases that catalyze the hydrolysis and formation of a variety of lactones, with a clear preference for lipophilic lactones. However, the lactonase mechanism of mammalian PON1, a high density lipoprotein-associated enzyme that is the most studied family member, remains unclear, and other family members have not been examined at all. We present a kinetic and site-directed mutagenesis study aimed at deciphering the lactonase mechanism of PON1 and PON3. The pH-rate profile determined for the lactonase activity of PON1 indicated an apparent pK(a) of similar to 7.4. We thus explored the role of all amino acids in the PON1 active site that are not directly ligated to the catalytic calcium and that possess an imidazolyl or carboxyl side chain (His(115), His(134), His(184), His(285), Asp(183), and Asp(269)). Extensive site-directed mutagenesis studies in which each amino acid candidate was replaced with all other 19 amino acids were conducted to identify the residue(s) that mediate the lactonase activity of PONs. The results indicate that the lactonase activity of PON1 and PON3 and the esterase activity of PON1 are mediated by the His(115)-His(134) dyad. Notably, the phosphotriesterase activity of PON1, which is a promiscuous activity of this enzyme, is mediated by other residues. To our knowledge, this is one of few examples of a histidine dyad in enzyme active sites and the first example of a hydrolytic enzyme with such a dyad.

An enzyme with completely new function can be created in the lab by mimicking natural evolutionary processes that both alter and preserve protein architecture.

The amidohydrolase superfamily comprises hundreds of hydrolytic enzymes of the (β/α)₈ 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.

Serum paraoxonase (PON1) is a high-density lipoprotein (HDL)-associated enzyme exhibiting antiatherogenic properties. This study examined the interaction of recombinant PON1 with reconstituted HDL comprised of PC, cholesterol, and various apolipoproteins (apoA-I, -II, and -IV). The affinity, stability, and lactonase activity were strongly correlated, with apoA-I exhibiting the strongest effects, apoA-IV exhibiting weaker yet significant effects, and apoA-II having a negative effect relative to protein-free particles. We found that PON1 binds apoA-I HDL with sub-nanomolar affinities (K_d ≪ 10^-9 M) and slow dissociation rates (t _1/2 > 80 min), while binding affinity for other particles was dramatically lower. A truncated form of PON1 lacking the N-terminal helix maintains considerable binding to apoA-I HDL (K_d = 1.2 × 10^-7 M), validating the structural model which indicates additional parts of the enzyme involved in HDL binding. Kinetic inactivation assays revealed the existence of an equilibrium between two forms of PON1 differing in their stability by a factor of 100. Various lipoproteins and detergent preparations shift this equilibrium toward the more stable conformation. Consistent with its highest affinity, only apoA-I HDL is capable of totally shifting the equilibrium toward the stable form. The paraoxonase and arylesterase activities were stimulated by HDL by 2-5-fold as previously reported, almost independently of the apoliporotein content. In contrast, only apoA-I is capable of stimulating the lactonase activity by ≤20-fold to k _cat/K_M values of 10⁶-10⁷ M ^-1 s^-1, while apoA-IV and apoA-II have almost no effect. Overall, the results indicate the high stability, selectivity, and catalytic proficiency of PON1 when anchored onto apoA-I HDL, toward lactone substrates, and lipophilic lactones in particular.

PON1 is the best-studied member of a family of enzymes called serum paraoxonases, or PONs, identified in mammals (including humans) and other vertebrates as well as in invertebrates. PONs exhibit a range of important activities, including drug metabolism and detoxification of organophosphates such as nerve agents. PON1 resides on HDL (the "good cholesterol") and is also involved in the prevention of atherosclerosis. Despite this wealth of activities, the identity of PON1's native substrate, namely, the substrate for which this enzyme and other enzymes from the PON family evolved, remains unknown. To elucidate the substrate preference and other details of PON1 mechanism of catalysis, structure - activity studies were performed with three groups of substrates that are known to be hydrolyzed by PON1: phosphotriesters, esters, and lactones. We found that the hydrolysis of aryl esters is governed primarily by steric factors and not the pK_a of the leaving group. The rates of hydrolysis of aliphatic esters are much slower and show a similar dependence on the pK_a of the leaving group to that of the nonenzymatic reactions in solution, while the aryl phosphotriesters show much higher dependence than the respective nonenzymatic reaction. PON1-catalyzed lactone hydrolysis shows almost no dependence on the pK_a of the leaving group, and unlike all other substrates, lactones seem to differ in their K_M rather than k_cat values. These, and the relatively high rates measured with several lactone substrates (k_cat/K _M ≈ 10⁶ M^-1 s^-1) imply that PON1 is in fact a lactonase.

Recent developments have been made in the application of directed evolution to achieve the efficient heterologous expression of proteins in Escherichia coli and yeast by increasing the stability and solubility of the protein in the host environment. One interesting conclusion that emerges is that the evolutionary process often improves the stability and solubility of an intermediate (apoprotein, proprotein or folding intermediate) that otherwise constitutes a bottleneck to functional expression, rather than altering the protein's final state.

Phosphotriesterase from Pseudomonas diminuta (PTE) is an extremely efficient metalloenzyme that hydrolyses a variety of compounds including organophosphorus nerve agents. Study of PTE has been hampered by difficulties with efficient expression of the recombinant form of this highly interesting and potentially useful enzyme. We identified a low-level esterolytic activity of PTE and then screened PTE gene libraries for improvements in 2-naphthyl acetate hydrolysis. However, the attempt to evolve this promiscuous esterase activity led to a variant (S5) containing three point mutations that resulted in a 20-fold increase in functional expression. Interestingly, the zinc holoenzyme form of S5 appears to be more sensitive than wild-type PTE to both thermal denaturation and addition of metal chelators. Higher functional expression of the S5 variant seems to lie in a higher stability of the metal-free apoenzyme. The results obtained in this work point out another-and often overlooked-possible determinant of protein expression and purification yields, i.e. the stability of intermediates during protein folding and processing.

Members of the serum paraoxonase (PON) family have been identified in mammals and other vertebrates, and in invertebrates. PONs exhibit a wide range of physiologically important hydrolytic activities, including drug metabolism and detoxification of nerve agents. PON1 and PON3 reside on high-density lipoprotein (HDL, 'good cholesterol') and are involved in the prevention of atherosclerosis. We describe the first crystal structure of a PON family member, a variant of PON1 obtained by directed evolution, at a resolution of 2.2 Å. PON1 is a six-bladed β-propeller with a unique active site lid that is also involved in HDL binding. The three-dimensional structure and directed evolution studies permit a detailed description of PON1's active site and catalytic mechanism, which are reminiscent of secreted phospholipase A2, and of the routes by which PON family members diverged toward different substrate and reaction selectivities.

Water-in-oil (w/o) emulsions can be used to compartmentalize and select large gene libraries for a predetermined function. The aqueous droplets of the w/o emulsion function as cell-like compartments in each of which a single gene is transcribed and translated to give multiple copies of the protein (e.g., an enzyme) it encodes. While compartmentalization ensures that the gene, the protein it encodes, and the products of the activity of this protein remain linked, it does not directly afford a way of selecting for the desired activity. Here we show that re-emulsification of w/o emulsions gives water-in-oil-in- water (w/o/w) emulsions with an external (continuous) water phase through which droplets containing fluorescent markers can be isolated by fluorescence- activated cell sorting (FACS). These w/o/w emulsions can be sorted by FACS, while the content of the aqueous droplets of the primary w/o emulsion remains intact. Consequently, genes embedded in these water droplets together with a fluorescent marker can be isolated and enriched from an excess of genes embedded in water droplets without a fluorescent marker. The ability of FACS instruments to sort up to 40,000 events per second may endow this technology a wide potential in the area of high-throughput screening and the directed evolution of enzymes.

Complex organisms have evolved from a limited number of primordial genes and proteins. However, the mechanisms by which the earliest proteins evolved and then served as the origin for the present diversity of protein function are unknown. Here, we outline a hypothesis based on the 'new view' of proteins whereby one sequence can adopt multiple structures and functions. We suggest that such conformational diversity could increase the functional diversity of a limited repertoire of sequences and, thereby, facilitate the evolution of new proteins and functions from old ones.

In vitro compartmentalisation (IVC), a technique for selecting genes encoding enzymes based on compartmentalising gene translation and enzymatic reactions in emulsions, was used to investigate the interaction of the DNA cytosine-5 methyltransferase M.HhaI with its target DNA (5-GCGC-3). Crystallography shows that the active site loop from the large domain of M.HhaI interacts with a flipped-out cytosine (the target for methylation) and two target recognition loops (loops I and II) from the small domain make almost all the other base-specific interactions. A library of M.HhaI genes was created by randomising all the loop II residues thought to make base-specific interactions and directly determine target specificity. The library was selected for 5-GCGC-3 methylation. Interestingly, in 11 selected active clones, 10 different sequences were found and none were wild-type. At two of the positions mutated (Ser252 and Tyr254) a number of different amino acids could be tolerated. At the third position, however, all active mutants had a glycine, as in wild-type M.HhaI, suggesting that Gly257 is crucial for DNA recognition and enzyme activity. Our results suggest that recognition of base pairs 3 and 4 of the target site either relies entirely on main chain interactions or that different residues from those identified in the crystal structure contribute to DNA recognition.

The x-ray structures of three esterase-like catalytic antibodies identified by screening for catalytic activity the entire hybridoma repertoire, elicited in response to a phosphonate transition state analog (TSA) hapten, were analyzed. The high resolution structures account for catalysis by transition state stabilization, and in all three antibodies a tyrosine residue participates in the oxyanion hole. Despite significant conformational differences in their combining sites, the three antibodies, which are the most efficient among those elicited, achieve catalysis in essentially the same mode, suggesting that evolution for binding to a single TSA followed by screening for catalysis lead to antibodies with structural convergence.

Background Antibodies with catalytic properties can be prepared by eliciting an antibody response against 'transition state analog' haptens. The specificity, rate and number of reaction cycles observed with these antibodies more closely resemble the properties of enzymes than any of the many other known enzyme- mimicking systems. Results We have determined to 3 å resolution the first X- ray structure of a catalytic antibody Fab. This antibody catalyzes the hydrolysis of a p-nitrophenyl ester. In conjunction with binding studies in solution, this structure of the uncomplexed site suggests a model for transition state fixation where two tyrosines mimic the oxyanion binding hole of serine proteases. A comparison with the structures of known Fabs specific for low molecular weight haptens reveals that this catalytic antibody has an unusually long groove at its combining site. Conclusions Since transition state analogs contain elements of the desired product, product inhibition is a severe problem in antibody catalysis. The observation of a long groove at the combining site may relate to the ability of this catalytic antibody to achieve multiple cycles of reaction.