Publications

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.

How individual enzymes evolved is relatively well understood. However, individual enzymes rarely confer a physiological advantage on their own. Judging by its current state, the emergence of metabolism seemingly demanded the simultaneous emergence of many enzymes. Indeed, how multicomponent interlocked systems, like metabolic pathways, evolved is largely an open question. This complexity can be unlocked if we assume that survival of the fittest applies not only to genes and enzymes but also to the metabolites they produce. This review develops our current knowledge of enzyme evolution into a wider hypothesis of pathway and network evolution. We describe the current models for pathway evolution and offer an integrative metabolite-enzyme coevolution hypothesis. Our hypothesis addresses the origins of new metabolites and of new enzymes and the order of their recruitment. We aim to not only survey established knowledge but also present open questions and potential ways of addressing them.

The linkage between regulatory elements of transcription, such as promoters, and their protein products is central to gene function. Promoter-protein coevolution is therefore expected, but rarely observed, and the manner by which these two regulatory levels are linked remains largely unknown. We study glutamate dehydrogenase-a 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 promoter-enzyme coevolution. Characterization of the proteins shows that, compared to RocG, GudB's enzymatic activity is highly dependent on glutamate and pH. Promoter-enzyme 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.

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, C23's 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 sequence-activity relationships along the pathway of PTE's 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(7) M-1 min(-1) 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 sequence-activity relationships.

Modifications of the bacterial ribosome regulate the function of the ribosome and modulate its susceptibility to antibiotics. By modifying a highly conserved adenosine A2503 in 23S rRNA, methylating enzyme Cfr confers resistance to a range of ribosome-targeting antibiotics. The same adenosine is also methylated by RlmN, an enzyme widely distributed among bacteria. While RlmN modifies C2, Cfr modifies the C8 position of A2503. Shared nucleotide substrate and phylogenetic relationship between RlmN and Cfr prompted us to investigate evolutionary origin of antibiotic resistance in this enzyme family. Using directed evolution of RlmN under antibiotic selection, we obtained RlmN variants that mediate low-level resistance. Surprisingly, these variants confer resistance not through the Cfr-like C8 methylation, but via inhibition of the endogenous RlmN C2 methylation of A2503. Detection of RlmN inactivating mutations in clinical resistance isolates suggests that the mechanism used by the in vitro evolved variants is also relevant in a clinical setting. Additionally, as indicated by a phylogenetic analysis, it appears that Cfr did not diverge from the RlmN family but from another distinct family of predicted radical SAM methylating enzymes whose function remains unknown.

The recent attacks with the nerve agent sarin in Syria reveal the necessity of effective countermeasures against highly toxic organophosphorus compounds. Multiple studies provide evidence that a rapid onset of antidotal therapy might be life-saving but current standard antidotal protocols comprising reactivators and competitive muscarinic antagonists show a limited efficacy for several nerve agents. We here set out to test the newly developed phosphotriesterase (PTE) mutant C23AL by intravenous (i.v.), intramuscular (i.m.; model for autoinjector) and intraosseous (i.o.; model for intraosseous insertion device) application in an in vivo guinea pig model after VX challenge (similar to 2LD(50)). C23AL showed a C-max of 0.63 mu mol L (1) after i.o. and i.v. administration of 2 mg kg (1) providing a stable plasma profile up to 180 min experimental duration with 0.41 and 0.37 mu mol L (1) respectively. The i.m. application of C23AL did not result in detectable plasma levels. All animals challenged with VX and subsequent i.o. or i.v. C23AL therapy survived although an in part substantial inhibition of erythrocyte, brain and diaphragm AChE was detected. Theoretical calculation of the time required to hydrolyze in vivo 96.75% of the toxic VX enantiomer is consistent with previous studies wherein similar activity of plasma containing catalytic scavengers of OPs resulted in non-lethal protection although accompanied with a variable severity of cholinergic symptoms. The relatively low C23AL plasma level observed immediately after its i.v. or i.o load, point at a possible volume of distribution greater than the guinea pig plasma content, and thus underlines the necessity of in vivo experiments in antidote research. In conclusion the i.o. application of PTE is efficient and resulted in comparable plasma levels to the i.v. application at a given time. Thus, i.o. vascular access systems could improve the post-exposure PTE therapy of nerve agent poisoning. (C) 2016 Elsevier Irela

Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase ( hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at similar to 2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20 degrees C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il.

Molecular evolution has focused on the divergence of molecular functions, yet we know little about how structurally distinct protein folds emerge de novo. We characterized the evolutionary trajectories and selection forces underlying emergence of b-propeller proteins, a globular and symmetric fold group with diverse functions. The identification of short propeller- like motifs (<50 amino acids) in natural genomes indicated that they expanded via tandem duplications to form extant propellers. We phylogenetically reconstructed 47-residue ancestral motifs that form five-bladed lectin propellers via oligomeric assembly. We demonstrate a functional trajectory of tandemduplications of these motifs leading to monomeric lectins. Foldability, i.e., higher efficiency of folding, was the main parameter leading to improved functionality along the entire evolutionary trajectory. However, folding constraints changed along the trajectory: initially, conflicts between monomer folding and oligomer assembly dominated, whereas subsequently, upon tandem duplication, tradeoffs between monomer stability and foldability took precedence.

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) end 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.

Organismal adaptation to extreme temperatures yields enzymes with distinct configurational stabilities, including thermophilic and psychrophilic enzymes, which are adapted to high and low temperatures, respectively. These enzymes are widely assumed to also have unique rate temperature dependencies. Thermophilic enzymes, for example, are considered optimal at high temperatures and effectively inactive at low temperatures due to excess rigidity. Surveying published data, we find that thermophilic, mesophilic, and psychrophilic enzymes exhibit indistinguishable rate temperature dependencies. Furthermore, given the nonenzymatic rate temperature dependency, all enzymes, regardless of their operation temperatures, become >10-fold less powerful catalysts per 25 C temperature increase. Among other factors, this loss of rate acceleration may be ascribed to thermally induced vibrations compromising the activesite catalytic configuration, suggesting that many enzymes are in fact insufficiently rigid.

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.

Keywords: HYDROLASES

Protein evolvability includes two elements-robustness (or neutrality, mutations having no effect) and innovability (mutations readily inducing new functions). How are these two conflicting demands bridged? Does the ability to bridge them relate to the observation that certain folds, such as TIM barrels, accommodate numerous functions, whereas other folds support only one? Here, we hypothesize that the key to innovability is polarity-an active site composed of flexible, loosely packed loops alongside a well-separated, highly ordered scaffold. We show that highly stabilized variants of TEM-1 beta-lactamase exhibit selective rigidification of the enzyme's scaffold while the active-site loops maintained their conformational plasticity. Polarity therefore results in stabilizing, compensatory mutations not trading off, but instead promoting the acquisition of new activities. Indeed, computational analysis indicates that in folds that accommodate only one function throughout evolution, for example, dihydrofolate reductase, >= 60% of the active-site residues belong to the scaffold. In contrast, folds associated with multiple functions such as the TIM barrel show high scaffold-active-site polarity (similar to 20% of the active site comprises scaffold residues) and >2-fold higher rates of sequence divergence at active-site positions. Our work suggests structural measures of fold polarity that appear to be correlated with innovability, thereby providing new insights regarding protein evolution, design, and engineering. (C) 2013 Elsevier Ltd. All rights reserved.

Recent advances in enzyme engineering enable dramatic improvements in catalytic efficiency and/or selectivity, as well as de novo engineering of enzymes to catalyze reactions where natural enzymes are not available. Can these capabilities be utilized to transform biosynthesis pathways? Metabolic engineering is traditionally based on combining existing enzymes to give new, or modified, pathways, within a new context and/or organism. How efficient, however, are the individual enzyme components? Is there room to improve pathway performance by enzyme engineering? We discuss the differences between enzymes in central versus specialized, or secondary metabolism and highlight unique features of specialized metabolism enzymes participating in the synthesis of natural products. We argue that, for the purpose of metabolic engineering, the catalytic efficiency and selectivity of many enzymes can be improved with the aim of achieving higher rates, yields and product purities. We also note the relative abundance of spontaneous reactions in specialized metabolism, and the potential advantage of engineering enzymes that will catalyze these steps. Specialized metabolism therefore offers new opportunities to integrate enzyme and pathway engineering, thereby achieving higher metabolic efficiencies, enhanced production rates and improved product purities.

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, eta) varied together, defining a scaling curve: eta(2) = b/m + eta(2)(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.

Arsenate and phosphate are abundant on Earth and have striking similarities: nearly identical pK(a) values(1,2), similarly charged oxygen atoms, and thermochemical radii that differ by only 4% (ref. 3). Phosphate is indispensable and arsenate is toxic, but this extensive similarity raises the question whether arsenate may substitute for phosphate in certain niches(4,5). However, whether it is used or excluded, discriminating phosphate from arsenate is a paramount challenge. Enzymes that utilize phosphate, for example, have the same binding mode and kinetic parameters as arsenate, and the latter's presence therefore decouples metabolism(6,7). Can proteins discriminate between these two anions, and how would they do so? In particular, cellular phosphate uptake systems face a challenge in arsenate-rich environments. Here we describe a molecular mechanism for this process. We examined the periplasmic phosphate-binding proteins (PBPs) of the ABC-type transport system that mediates phosphate uptake into bacterial cells, including two PBPs from the arsenate-rich Mono Lake Halomonas strain GFAJ-1. All PBPs tested are capable of discriminating phosphate over arsenate at least 500-fold. The exception is one of the PBPs of GFAJ-1 that shows roughly 4,500-fold discrimination and its gene is highly expressed under phosphate-limiting conditions. Sub-angstrom-resolution structures of Pseudomonas fluorescens PBP with both arsenate and phosphate show a unique mode of binding that mediates discrimination. An extensive network of dipole-anion interactions(8,9), and of repulsive interactions, results in the 4% larger arsenate distorting a unique low-barrier hydrogen bond. These features enable the phosphate transport system to bind phosphate selectively over arsenate (at least 10(3) excess) even in highly arsenate-rich environments.

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 (3150 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) Biochemistry 45, 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 X 10(4)). 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.

Computational design is a test of our understanding of enzyme catalysis and a means of engineering novel, tailor-made enzymes. While the de novo computational design of catalytically efficient enzymes remains a challenge, designed enzymes may comprise unique starting points for further optimization by directed evolution. Directed evolution of two computationally designed Kemp eliminases, KE07 and KE70, led to low to moderately efficient enzymes (k(cat)/K-m values of 2,000-fold increase in catalytic efficiency, mainly via higher k(cat) values. The best KE59 variants exhibited k(cat)/K-m values up to 0.6 x 10(6) M(-1)s(-1), and k(cat)/k(uncat) values of

A preferred strategy for preventing nerve agents intoxication is catalytic scavenging by enzymes that hydrolyze them before they reach their targets. Using directed evolution, we simultaneously enhanced the activity of a previously described serum paraoxonase 1 (PON1) variant for hydrolysis of the toxic S-P isomers of the most threatening G-type nerve agents. The evolved variants show 2,500 for the S-P isomer in an evolved variant. Given their ability to hydrolyze G-agents, these evolved variants may serve as broad-range G-agent prophylactics.

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. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

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 similar to 90% of proteins, the distribution of positional rates exhibits three peaks: a peak of slow evolving residues, with average log(2)[normalized rate], log(2)mu, of ca. -2, corresponding primarily to core residues; a peak of fast evolving residues (log(2)mu similar to 0.5) largely corresponding to surface residues; and a very fast peak (log(2)mu similar to 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(2)mu similar to-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.

Whether evolution is erratic due to random historical details, or is repeatedly directed along similar paths by certain constraints, remains unclear. Epistasis (i.e. non-additive interaction between mutations that affect fitness) is a mechanism that can contribute to both scenarios. Epistasis can constrain the type and order of selected mutations, but it can also make adaptive trajectories contingent upon the first random substitution. This effect is particularly strong under sign epistasis, when the sign of the fitness effects of a mutation depends on its genetic background. In the current study, we examine how epistatic interactions between mutations determine alternative evolutionary pathways, using in vitro evolution of the antibiotic resistance enzyme TEM-1 beta-lactamase. First, we describe the diversity of adaptive pathways among replicate lines during evolution for resistance to a novel antibiotic (cefotaxime). Consistent with the prediction of epistatic constraints, most lines increased resistance by acquiring three mutations in a fixed order. However, a few lines deviated from this pattern. Next, to test whether negative interactions between alternative initial substitutions drive this divergence, alleles containing initial substitutions from the deviating lines were evolved under identical conditions. Indeed, these alternative initial substitutions consistently led to lower adaptive peaks, involving more and other substitutions than those observed in the common pathway. We found that a combination of decreased enzymatic activity and lower folding cooperativity underlies negative sign epistasis in the clash between key mutations in the common and deviating lines (Gly238Ser and Arg164Ser, respectively). Our results demonstrate that epistasis contributes to contingency in protein evolution by amplifying the selective consequences of random mutations.

Organophosphate nerve agents are extremely lethal compounds. Rapid in vivo organophosphate clearance requires bioscavenging enzymes with catalytic efficiencies of >10(7) (M(-1) min(-1)). Although serum paraoxonase (PON1) is a leading candidate for such a treatment, it hydrolyzes the toxic S(p) isomers of G-agents with very slow rates. We improved PON1's catalytic efficiency by combining random and targeted mutagenesis with high-throughput screening using fluorogenic analogs in emulsion compartments. We thereby enhanced PON1's activity toward the coumarin analog of S(p)-cyclosarin by similar to 10(5)-fold. We also developed a direct screen for protection of acetylcholinesterase from inactivation by nerve agents and used it to isolate variants that degrade the toxic isomer of the coumarin analog and cyclosarin itself with k(cat)/K(M) similar to 10(7) M(-1) min(-1). We then demonstrated the in vivo prophylactic activity of an evolved variant. These evolved variants and the newly developed screens provide the basis for engineering PON1 for prophylaxis against other G-type agents.

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. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

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 similar to 100 amino acid fragments derived from tachylectin-2-a monomeric, 236 amino acids, five-bladed beta-propeller. Their structures reveal a unique pentameric assembly and novel beta-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 beta-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.

Serum paraoxonase (PON1) is all anti-atherogenic interfacially activated lipo-lactonase that Was shown to selectively bind high-density lipoprotein (HDL) carrying apolipoprotein A-1 (apoA-1). ApoA-1 binding occurs with nanomolar affinity and Induces a dramatic increase in enzyme stability and lactonase activity. This study examined the association of PON1 With reconstituted HDL (rHDL) carrying apolipoprotein E, and its consequences on the stability and enzymatic activity of PON1, and on its anti-atherogenic potential. The results indicate that reconstituted HDL particles prepared With two most common isoforms of apoE (apoE3 and apoE4) associate with rePON1 in a manner and affinity similar to those of apoA-1. Binding to apoE-HDL stimulates the lactonase activity and stabilizes the enzyme, although the latter Occurs to a > 10-rold lesser extent compared to apoA-I-HDL particles. The anti-atherogenic potential of PON1, measured by inhibition of LDL oxidation and stimulation of macrophage cholesterol efflux, was also stimulated by apoE-HDL, at levels of 40-96% compared to apoA-1-HDL. Overall, reconstituted apoE-HDL exhibits properties similar to those of apoA-1-HDL, but with a lower capacity to stabilize PON1 and to induce its anti-atherogenic functions. ApoE, apoA-1, and to a lesser degree apoA-IV show distinct structural and functional similarities but little sequence homology. That these apolipoproteins. but not apoA-11, bind PON1 with high affinity and stimulate its activity Suggests that PON1-HDL recognition is based primarily oil surface properties of the apolipoproteins and that specific protein-protein interactions may play only a secondary role.

The past several years have seen novel insights at the interface of protein biophysics and evolution. The accepted paradigm that proteins can tolerate nearly any amino acid substitution has been replaced by the view that the deleterious effects of mutations, and especially their tendency to undermine the thermodynamic and kinetic stability of protein, is a major constraint on protein evolvability-the ability of proteins to acquire changes in sequence and function. We summarize recent findings regarding how mutations affect protein stability, and how stability affects protein evolution. We describe ways of predicting and analyzing stability effects of mutations, and mechanisms that buffer or compensate for these destabilizing effects and thereby promote protein evolvabilty, in nature and in the laboratory.

Keywords: Biochemistry & Molecular Biology

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

The divergence of new gene functions is described by various scenarios that involve gene duplication, albeit, at fundamentally different stages. We performed experimental measurements and developed a subsequent model, aimed at predicting the rate of divergence under different scenarios. We used gene libraries of TEM-1 beta-lactamase that were drifted under purifying selection toward the original penicillinase activity or under no selection at all. The frequency of genes conferring a new function (degradation of a cephalosporin antibiotic) was measured at various stages of the drift, and a model that accounts for the differences in the observed adaptation dynamics of the drifting TEM-1 populations was derived. The results indicate that rapid nonfunctionalization in the population relieved from selection (Ohno's model) afforded only a narrow window of adaptation to cefotaxime (neofunctionalization). The trade-off between TEM-1's original function and the new evolving function also disfavored the "gene sharing" model. The rate of adaptation was maximal when selection for the original function was partially relieved to enable the accumulation of potentially adaptive mutations, while still purging a large fraction of otherwise deleterious mutations. Altogether, scenarios of subfunctionalization seem more feasible, whereby sustaining the original function by two copies facilitates the accumulation potentially adaptive mutations while purging nonfunctionalization mutations.

The design of new enzymes for reactions not catalysed by naturally occurring biocatalysts is a challenge for protein engineering and is a critical test of our understanding of enzyme catalysis. Here we describe the computational design of eight enzymes that use two different catalytic motifs to catalyse the Kemp elimination - a model reaction for proton transfer from carbon - with measured rate enhancements of up to 10 5 and multiple turnovers. Mutational analysis confirms that catalysis depends on the computationally designed active sites, and a high- resolution crystal structure suggests that the designs have close to atomic accuracy. Application of in vitro evolution to enhance the computational designs produced a >200-fold increase in k(cat)/K-m (k(cat)/K-m of 2,600 M(-1)s(-1) and k(cat)/k(uncat) of >10(6)). These results demonstrate the power of combining computational protein design with directed evolution for creating new enzymes, and we anticipate the creation of a wide range of useful new catalysts in the future.

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 FAGS (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 HDLapoA-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 5 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 as well as the levels and quality of the HDL particles with which the enzyme is associated.

How the thermodynamic stability effects of protein mutations (Delta Delta 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 Delta Delta 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 Delta Delta G (similar to 0.6 kcal/mol), whereas the core residues exhibit a wider distribution with a stronger destabilizing mean (similar to 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. (C) 2007 Elsevier Ltd. All rights reserved.

The directed evolution of proteins has benefited greatly from site-specific methods of diversification such as saturation mutagenesis. These techniques target diversity to a number of chosen positions that are usually non-contiguous in the protein's primary structure. However, the number of targeted positions can be large, thus leading to impractically large library size, wherein almost all library variants are inactive and the likelihood of selecting desirable properties is extremely small. We describe a versatile combinatorial method for the partial diversification of large sets of residues. Our library oligonucleotides comprise randomized codons that are flanked by wild-type sequences. Adding these oligonucleotides to an assembly PCR of wildtype gene fragments incorporates the randomized cassettes, at their target sites, into the reassembled gene. Varying the oligonucleotides concentration resulted in library variants that carry a different average number of mutated positions that comprise a random subset of the entire set of. diversified codons. This method, dubbed Incorporating Synthetic Oligos via Gene Reassembly (ISOR), was used to create libraries of a cytosine-C5 methyltransferase wherein 45 individual positions were randomized. One library, containing an average of 5.6 mutated residues per gene, was selected, and mutants with wild-type-like activities isolated. We also created libraries of serum paraoxonase PON1 harboring insertions and deletions (indels) in various areas surrounding the active site. Screening these libraries yielded a range of mutants with altered substrate specificities and indicated that certain regions of this enzyme have a surprisingly high tolerance to indels.

Serum paraoxonase (PON1) is a HDL-associated enzyme exhibiting potentially antiatherogenic properties. Here, we examined the common PON1-192R/Q human polymorphism. Despite numerous studies, the effect of this polymorphism on the antiatherogenic potential of PON1 is yet unresolved. Our structural model suggests that amino acid 192 constitutes part of the HDL-anchoring surface and active site of PON1. Based on our findings that PON1 is an interfacially activated lipolactonase that selectively binds HDL carrying apolipoprotein A-I (apoA-I) and is thereby greatly stabilized and catalytically activated, we examined the interaction of the PON1-192 isozymes with reconstituted HDL-apoA-I particles. We found that PON1 position 192 is indeed involved in HDL binding. The PON1-192Q binds HDL with a 3-fold lower affinity than the R isozyme and consequently exhibits significantly reduced stability, lipolactonase activity, and macrophage cholesterol efflux. We also observed the lower affinity and stability of the 192Q versus the 192R isozyme in sera of individuals belonging to the corresponding genotypes. The observed differences in the properties of PON1-192R/Q isozymes provide a basis for further analysis of the contribution of the 192R/Q polymorphism to the susceptibility to atherosclerosis, although other factors, such as the overall levels of PON1, may play a more significant role.

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 goal of in vitro compartmentalization (IVC) is to divide a large reaction between many microscopic compartments(1). This technique was first developed to generate 'artificial cells' for the directed evolution of proteins (Fig. 1). Typically, an aqueous solution of genes and an in vitro transcription-translation system is stirred (or homogenized) into an oil-surfactant mixture to create a water-in-oil (w/o) emulsion with similar to 10(10) aqueous droplets per ml of emulsion. The majority of droplets contain no more than a single gene along with all of the molecular machinery needed to express that gene. The expressed proteins and the products of their catalytic activities cannot leave the droplets, and so genotype is coupled to phenotype in vitro, making it possible to select very large libraries of genes (10(8)-10(11) genes). We describe the advantages and applications of IVC in Box 1. Here we present a protocol for performing a directed evolution experiment by IVC that makes use of one or more w/o emulsions. This procedure involves the generation of a gene library, the performance of a selection, and the subsequent recovery of the selected genes by PCR. We also describe two procedures for converting w/o emulsions to water-in-oil-in- water (w/o/w) emulsions for high-throughput screening using a fluorescence-activated cell sorter (FACS; Box 2 and Fig. 2, and Box 3). Finally, we describe two methods for delivering substrates, regulators and other compounds to the preformed aqueous droplets of a w/o emulsion (Box 4).

We addressed the ability of various organophosphorus (OP) hydrolases to catalytically scavenge toxic OP nerve agents. Mammalian paraoxonase (PON1) was found to be more active than Pseudomonas diminuta OP hydrolase (OPH) and squid O,O-di-isopropyl fluorophosphatase (DFPase) in detoxifying cyclosarin (O-cyclohexyl methylphosphonofluoridate) and soman (O-pinacolyl methylphosphonofluoridate). Subsequently, nine directly evolved PON1 variants, selected for increased hydrolytic rates with a fluorogenic diethylphosphate ester, were tested for detoxification of cyclosarin, soman, O-isopropyl-O-(p-nitrophenyl) methyl phosphonate (IMP-pNP), DFP, and chlorpyrifos-oxon (ChPo). Detoxification rates were determined by temporal acetylcholinesterase inhibition by residual nonhydrolyzed OP. As stereoisomers of cyclosarin and soman differ significantly in their acetylcholinesterase-inhibiting potency, we actually measured the hydrolysis of the more toxic stereoisomers. Cyclosarin detoxification was similar to 10-fold faster with PONI mutants V346A and L69V. V346A also exhibited fourfold and sevenfold faster hydrolysis of DFP and ChPo, respectively, compared with wild-type, and ninefold higher activity towards soman. L69V exhibited 100-fold faster hydrolysis of DFP than the wildtype. The active-site mutant H115W exhibited 270-380-fold enhancement toward hydrolysis of the P-S bond in parathiol, a phosphorothiolate analog of parathion. This study identifies three key positions in PON1 that affect OP hydrolysis, Leu69, Val346 and His115, and several amino-acid replacements that significantly enhance the hydrolysis of toxic OPs. GC/pulsed flame photometer detector analysis, compared with assay of residual acetylcholinesterase inhibition, displayed stereoselective hydrolysis of cyclosarin, soman, and IMP-pNP, indicating that PONI is less active toward the more toxic optical isomers.

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.

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 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 x 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

induced fit is a predominant phenomenon in protein-ligand interactions, yet it is invariably attributed without establishing the existence, let alone the structure, of the initial, low-affinity encounter complex. We determined the crystal structure of the encounter complex on the pathway of ligand binding by IgE antibody SPE7. We show that this complex is formed by a wide range of ligands that initially bind with identical affinity. Nonspecific ligands rapidly dissociate, whereupon the antibody isomerizes to a nonbinding isomer. Specific ligand complexes, however, slowly isomerize to give a high-affinity complex. This isomerization involves backbone and side-chain rearrangements of up to 14 angstrom and the formation of specific hydrogen bonds. The postbinding conformational switch, combined with the prebinding isomerization to an energetically favorable nonbinding isomer, results in a "kinetic discrimination" mechanism that mediates selective binding, by a factor of > 10(3), between highly related ligands that initially bind with the same affinity. This model may apply to proteins that bind multiple ligands in a specific manner or other proteins that, although capable of binding many ligands, are activated by only a few.

The availability of vast gene repertoires from both natural sources (genomic and cDNA libraries) and artificial sources (gene libraries) demands the development and application of novel technologies that enable the screening or selection of large libraries for a variety of enzymatic activities. We describe recent developments in the selection of enzyme-coding genes for directed evolution and functional genomics. We focus on HTS approaches that enable selection from large libraries (> 10(6) gene variants) with relatively humble means (i.e. non-robotic systems), and on in vitro compartmentalization in particular.

Keywords: Chemistry, Multidisciplinary

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.

How proteins with new functions (e. g., drug or antibiotic resistance or degradation of man-made chemicals) evolve in a matter of months or years is still unclear. This ability is dependent on the induction of new phenotypic traits by a small number of mutations (plasticity). But mutations often have deleterious effects on functions that are essential for survival. How are these seemingly conflicting demands met at the single-protein level? Results from directed laboratory evolution experiments indicate that the evolution of a new function is driven by mutations that have little effect on the native function but large effects on the promiscuous functions that serve as starting point. Thus, an evolving protein can initially acquire increased fitness for a new function without losing its original function. Gene duplication and the divergence of a completely new protein may then follow.

We present a new and facile method to evaluate w/o/w emulsions containing fluorescent markers by flow cytometry. Flow cytometry allows simultaneous measurement of w/o/w emulsion droplets "marked" with a fluorescent marker or "blank" without the need for complicated sample preparation. The yield of preparation of the w/o/w emulsion and the release rate of the fluorescent marker FITC-BSA were investigated by this new method. The release fraction (after 24 h) of FITC-BSA from the w/o/w emulsion decreased with increasing concentration of FITC-BSA inside the internal phase, just like the release fraction of NaCl as marker from the w/o/w emulsion. Flow cytometry results show that the yield and release behavior in w/o/w emulsions are in agreement with results reported by more complicated methods.

Serum paraoxonases (PONs) are a group of enzymes that play a key role in organophosphate (OP) detoxification and in prevention of atherosclerosis. However, their structure and mechanism of action are poorly understood. PONs seem like jacks-of-all-trades, acting on a very wide range of substrates, most of which are of no physiological relevance. Family shuffling and screening lead to the first PON1 variants that express in a soluble and active form in Escherichia coli. We describe variants with kinetic parameters similar to those reported for PONs purified from sera and others that show dramatically increased activities. In particular, we have evolved POW variants with OP-hydrolyzing activities 40-fold higher than wild type and a specificity switch of >2,000-fold, producing PONs specialized for OP rather than ester hydrolysis. Analysis of the newly evolved variants provides insights into the evolutionary relationships between different family members.

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'). Crystallog raphy 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.

A number of monoclonal antibodies elicited against a nitrobenzyl (Nbzl)-phosphonate transition-state analogue (TSA), and which were selected for the hydrolysis of the corresponding Nbzl-ester, were also found to catalyze the hydrolysis of the analogous p-nitrophenyl(Np) ester with notable efficiency and specificity. The activity towards the Np-ester is higher in terms of rates (k(cat); as expected from the higher intrinsic reactivity of Np-esters); however, the rate acceleration (k(cat)/k(uncat)) is close to or lower than that observed with the Nbzl-ester. Unexpectedly, the affinity to the Np-ester substrate (1/K-M) and therefore k(cat)/K-M are significantly higher. The best example is antibody D2.4 having a k(cat)/K-M value of 64 s(-1). M(-1) with the Nbzl-ester and 9400 s-1 . M(-1) with the Np-ester. Moreover, due to a lower product inhibition by p-nitrophenol relative to p-nitrobenzyl alcohol, these antibodies exhibit more than 1000 turnovers with the Np-ester. The differential affinity of these antibodies to the Nbzl phosphonate TSA versus the Nbzl-ester substrate (K-S/K-TSA or K-M/K-i) correlates well with the observed rate enhancement (k(cat)/k(uncat)). For the Np-ester, however, stabilisation of the transition state (as reflected by K-S/K-TSA and by the catalytic proficiencies, k(cat)/K-M/k(uncat)) does not fully account for the catalytic power (k(cat)/k(uncat)), indicating a more complex catalytic mechanism than simply transition-state stabilization. A comparison of the kinetic parameters of D2.4 with other Np-ester-hydrolyzing antibodies raised against Np-phosphonate haptens emphasizes the marked advantage of this antibody which was elicited against an Nbzl-phosphonate hapten. These results appear to be general: anti-(Nbzl-phosphonate TSA) antibodies obtained from other mouse strains and using different immunization protocols are also efficient Np-esterases. They demonstrate the use of an expanded TSA-hapten, where a spacer (a methylene group) mimics bon

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 angstrom 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. Conclusion: 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.

The low abundance and activity of catalytic antibodies are major obstacles to their selection from the virtually unlimited repertoire of antibody binding sites. The requirement for new screening methodologies is further emphasized by the availability of combinatorial libraries, in which a functional polypeptide has to be selected out of millions of possibilities. We present a simple and sensitive screening approach (termed catELISA) based on immobilized substrates and immunodetection of the end product of the catalyzed reaction. The feasibility of catELISA is demonstrated here by the generation of potent ester-hydrolyzing antibodies by direct screening of hybridoma supernatants. We show that this approach is not only facile but general: it is not limited by type of reaction, substrate, or catalyst (enzymes, catalytic antibodies, chemical catalysts). catELISA opens a route to catalytic antibodies that replaces existing lengthy and arduous methods, thus allowing us to expand their number and improve their quality and to address questions that would otherwise be difficult to answer.