Directed Evolution of Protein Catalysts.

Directed evolution is a powerful technique for generating tailor-made enzymes for a wide range of biocatalytic applications. Following the principles of natural evolution, iterative cycles of mutagenesis and screening or selection are applied to modify protein properties, enhance catalytic activities, or develop completely new protein catalysts for non-natural chemical transformations. This review briefly surveys the experimental methods used to generate genetic diversity and screen or select for improved enzyme variants. Emphasis is placed on a key challenge, namely how to generate novel catalytic activities that expand the scope of natural reactions. Two particularly effective strategies, exploiting catalytic promiscuity and rational design, are illustrated by representative examples of successfully evolved enzymes. Opportunities for extending these approaches to more complex biocatalytic systems are also considered.

Tight and specific lanthanide binding in a de novo TIM barrel with a large internal cavity designed by symmetric domain fusion.

De novo protein design has succeeded in generating a large variety of globular proteins, but the construction of protein scaffolds with cavities that could accommodate large signaling molecules, cofactors, and substrates remains an outstanding challenge. The long, often flexible loops that form such cavities in many natural proteins are difficult to precisely program and thus challenging for computational protein design. Here we describe an alternative approach to this problem. We fused two stable proteins with C2 symmetry-a de novo designed dimeric ferredoxin fold and a de novo designed TIM barrel-such that their symmetry axes are aligned to create scaffolds with large cavities that can serve as binding pockets or enzymatic reaction chambers. The crystal structures of two such designs confirm the presence of a 420 cubic Ångström chamber defined by the top of the designed TIM barrel and the bottom of the ferredoxin dimer. We functionalized the scaffold by installing a metal-binding site consisting of four glutamate residues close to the symmetry axis. The protein binds lanthanide ions with very high affinity as demonstrated by tryptophan-enhanced terbium luminescence. This approach can be extended to other metals and cofactors, making this scaffold a modular platform for the design of binding proteins and biocatalysts.

Minor Actinides Can Replace Essential Lanthanides in Bacterial Life.

Certain f-block elements-the lanthanides-have biological relevance in the context of methylotrophic bacteria. The respective strains incorporate these 4 f elements into the active site of one of their key metabolic enzymes, a lanthanide-dependent methanol dehydrogenase. In this study, we investigated whether actinides, the radioactive 5 f elements, can replace the essential 4 f elements in lanthanide-dependent bacterial metabolism. Growth studies with Methylacidiphilum fumariolicum SolV and the Methylobacterium extorquens AM1 ΔmxaF mutant demonstrate that americium and curium support growth in the absence of lanthanides. Moreover, strain SolV favors these actinides over late lanthanides when presented with a mixture of equal amounts of lanthanides together with americium and curium. Our combined in vivo and in vitro results establish that methylotrophic bacteria can utilize actinides instead of lanthanides to sustain their one-carbon metabolism if they possess the correct size and a +III oxidation state.

RNA specificity and regulation of catalysis in the eukaryotic polynucleotide kinase Clp1.

RNA-specific polynucleotide kinases of the Clp1 subfamily are key components of various RNA maturation pathways. However, the structural basis explaining their substrate specificity and the enzymatic mechanism is elusive. Here, we report crystal structures of Clp1 from Caenorhabditis elegans (ceClp1) in a number of nucleotide- and RNA-bound states along the reaction pathway. The combined structural and biochemical analysis of ceClp1 elucidates the RNA specificity and lets us derive a general model for enzyme catalysis of RNA-specific polynucleotide kinases. We identified an RNA binding motif referred to as "clasp" as well as a conformational switch that involves the essential Walker A lysine (Lys127) and regulates the enzymatic activity of ceClp1. Structural comparison with other P loop proteins, such as kinases, adenosine triphosphatases (ATPases), and guanosine triphosphatases (GTPases), suggests that the observed conformational switch of the Walker A lysine is a broadly relevant mechanistic feature.

Emergence of a Negative Activation Heat Capacity during Evolution of a Designed Enzyme.

Temperature influences the reaction kinetics and evolvability of all enzymes. To understand how evolution shapes the thermodynamic drivers of catalysis, we optimized the modest activity of a computationally designed enzyme for an elementary proton-transfer reaction by nearly 4 orders of magnitude over 9 rounds of mutagenesis and screening. As theorized for primordial enzymes, the catalytic effects of the original design were almost entirely enthalpic in origin, as were the rate enhancements achieved by laboratory evolution. However, the large reductions in ΔH⧧ were partially offset by a decrease in TΔS⧧ and unexpectedly accompanied by a negative activation heat capacity, signaling strong adaptation to the operating temperature. These findings echo reports of temperature-dependent activation parameters for highly evolved natural enzymes and are relevant to explanations of enzymatic catalysis and adaptation to changing thermal environments.

Optimization of Enzyme Mechanism along the Evolutionary Trajectory of a Computationally Designed (Retro-)Aldolase.

De novo biocatalysts have been successfully generated by computational design and subsequent experimental optimization. Here, we examined the evolutionary history of the computationally designed (retro-)aldolase RA95. The modest activity of the starting enzyme was previously improved 105-fold over many rounds of mutagenesis and screening to afford a proficient biocatalyst for enantioselective cleavage and synthesis of ß-hydroxyketones. Using a set of representative RA95 variants, we probed individual steps in the multistep reaction pathway to determine which processes limit steady-state turnover and how mutations that accumulated along the evolutionary trajectory influenced the kinetic mechanism. We found that the overall rate-limiting step for aldol cleavage shifted from C-C bond scission (or an earlier step in the pathway) for the computational design to product release for the evolved enzymes. Specifically, interconversion of Schiff base and enamine intermediates, formed covalently between acetone and the catalytic lysine residue, was found to be the slowest step for the most active variants. A complex hydrogen bond network of four active site residues, which was installed in the late stages of laboratory evolution, apparently enhances lysine reactivity and facilitates efficient proton shuffling. This catalytic tetrad accounts for the tremendous rate acceleration observed for all steps of the mechanism, most notably Schiff base formation and hydrolysis. Comparison of our results with kinetic and structural studies on natural aldolases provides valuable feedback for computational enzyme design and laboratory evolution approaches alike.

Characterizing Active Site Conformational Heterogeneity along the Trajectory of an Enzymatic Phosphoryl Transfer Reaction.

States along the phosphoryl transfer reaction catalyzed by the nucleoside monophosphate kinase UmpK were captured and changes in the conformational heterogeneity of conserved active site arginine side-chains were quantified by NMR spin-relaxation methods. In addition to apo and ligand-bound UmpK, a transition state analog (TSA) complex was utilized to evaluate the extent to which active site conformational entropy contributes to the transition state free energy. The catalytically essential arginine side-chain guanidino groups were found to be remarkably rigid in the TSA complex, indicating that the enzyme has evolved to restrict the conformational freedom along its reaction path over the energy landscape, which in turn allows the phosphoryl transfer to occur selectively by avoiding side reactions.

trans-Acting arginine residues in the AAA+ chaperone ClpB allosterically regulate the activity through inter- and intradomain communication.

The molecular chaperone ClpB/Hsp104, a member of the AAA+ superfamily (ATPases associated with various cellular activities), rescues proteins from the aggregated state in collaboration with the DnaK/Hsp70 chaperone system. ClpB/Hsp104 forms a hexameric, ring-shaped complex that functions as a tightly regulated, ATP-powered molecular disaggregation machine. Highly conserved and essential arginine residues, often called arginine fingers, are located at the subunit interfaces of the complex, which also harbor the catalytic sites. Several AAA+ proteins, including ClpB/Hsp104, possess a pair of such trans-acting arginines in the N-terminal nucleotide binding domain (NBD1), both of which were shown to be crucial for oligomerization and ATPase activity. Here, we present a mechanistic study elucidating the role of this conserved arginine pair. First, we found that the arginines couple nucleotide binding to oligomerization of NBD1, which is essential for the activity. Next, we designed a set of covalently linked, dimeric ClpB NBD1 variants, carrying single subunits deficient in either ATP binding or hydrolysis, to study allosteric regulation and intersubunit communication. Using this well defined environment of site-specifically modified, cross-linked AAA+ domains, we found that the conserved arginine pair mediates the cooperativity of ATP binding and hydrolysis in an allosteric fashion.

The molecular mechanism of Hsp100 chaperone inhibition by the prion curing agent guanidinium chloride.

The Hsp100 chaperones ClpB and Hsp104 utilize the energy from ATP hydrolysis to reactivate aggregated proteins in concert with the DnaK/Hsp70 chaperone system, thereby playing an important role in protein quality control. They belong to the family of AAA+ proteins (ATPases associated with various cellular activities), possess two nucleotide binding domains per monomer (NBD1 and NBD2), and oligomerize into hexameric ring complexes. Furthermore, Hsp104 is involved in yeast prion propagation and inheritance. It is well established that low concentrations of guanidinium chloride (GdmCl) inhibit the ATPase activity of Hsp104, leading to so called "prion curing," the loss of prion-related phenotypes. Here, we present mechanistic details about the Hsp100 chaperone inhibition by GdmCl using the Hsp104 homolog ClpB from Thermus thermophilus. Initially, we demonstrate that NBD1 of ClpB, which was previously considered inactive as a separately expressed construct, is a fully active ATPase on its own. Next, we show that only NBD1, but not NBD2, is affected by GdmCl. We present a crystal structure of ClpB NBD1 in complex with GdmCl and ADP, showing that the Gdm(+) ion binds specifically to the active site of NBD1. A conserved essential glutamate residue is involved in this interaction. Additionally, Gdm(+) interacts directly with the nucleotide, thereby increasing the nucleotide binding affinity of NBD1. We propose that both the interference with the essential glutamate and the modulation of nucleotide binding properties in NBD1 is responsible for the GdmCl-specific inhibition of Hsp100 chaperones.

Elements in nucleotide sensing and hydrolysis of the AAA+ disaggregation machine ClpB: a structure-based mechanistic dissection of a molecular motor.

ATPases of the AAA+ superfamily are large oligomeric molecular machines that remodel their substrates by converting the energy from ATP hydrolysis into mechanical force. This study focuses on the molecular chaperone ClpB, the bacterial homologue of Hsp104, which reactivates aggregated proteins under cellular stress conditions. Based on high-resolution crystal structures in different nucleotide states, mutational analysis and nucleotide-binding kinetics experiments, the ATPase cycle of the C-terminal nucleotide-binding domain (NBD2), one of the motor subunits of this AAA+ disaggregation machine, is dissected mechanistically. The results provide insights into nucleotide sensing, explaining how the conserved sensor 2 motif contributes to the discrimination between ADP and ATP binding. Furthermore, the role of a conserved active-site arginine (Arg621), which controls binding of the essential Mg2+ ion, is described. Finally, a hypothesis is presented as to how the ATPase activity is regulated by a conformational switch that involves the essential Walker A lysine. In the proposed model, an unusual side-chain conformation of this highly conserved residue stabilizes a catalytically inactive state, thereby avoiding unnecessary ATP hydrolysis.

A sterile alpha-motif domain in NafY targets apo-NifDK for iron-molybdenum cofactor delivery via a tethered domain.

NafY participates in the final steps of nitrogenase maturation, having a dual role as iron-molybdenum cofactor (FeMo-co) carrier and as chaperone to the FeMo-co-deficient apo-NifDK (apo-dinitrogenase). NafY contains an N-terminal domain of unknown function (n-NafY) and a C-terminal domain (core-NafY) necessary for FeMo-co binding. We show here that n-NafY and core-NafY have very weak interactions in intact NafY. The NMR structure of n-NafY reveals that it belongs to the sterile α-motif (SAM) family of domains, which are frequently involved in protein-protein interactions. The presence of a SAM domain in NafY was unexpected and could not be inferred from its amino acid sequence. Although SAM domains are very commonly found in eukaryotic proteins, they have rarely been identified in prokaryotes. The n-NafY SAM domain binds apo-NifDK. As opposed to full-length NafY, n-NafY impaired FeMo-co insertion when present in molar excess relative to FeMo-co and apo-NifDK. The implications of these observations are discussed to offer a plausible mechanism of FeMo-co insertion. NafY domain structure, molecular tumbling, and interdomain motion, as well as NafY interaction with apo-NifDK are consistent with the function of NafY in FeMo-co delivery to apo-NifDK.

Coupling of oligomerization and nucleotide binding in the AAA+ chaperone ClpB.

Members of the family of ATPases associated with various cellular activities (AAA+) typically form homohexameric ring complexes and are able to remodel their substrates, such as misfolded proteins or protein-protein complexes, in an ATP-driven process. The molecular mechanism by which ATP hydrolysis is coordinated within the multimeric complex and the energy is converted into molecular motions, however, is poorly understood. This is partly due to the fact that the oligomers formed by AAA+ proteins represent a highly complex system and analysis depends on simplification and prior knowledge. Here, we present nucleotide binding and oligomer assembly kinetics of the AAA+ protein ClpB, a molecular chaperone that is able to disaggregate protein aggregates in concert with the DnaK chaperone system. ClpB bears two AAA+ domains (NBD1 and NBD2) on one subunit and forms homohexameric ring complexes. In order to dissect individual mechanistic steps, we made use of a reconstituted system based on two individual constructs bearing either the N-terminal (NBD1) or the C-terminal AAA+ domain (NBD2). In contrast to the C-terminal construct, the N-terminal construct does not bind the fluorescent nucleotide MANT-dADP in isolation. However, sequential mixing experiments suggest that NBD1 obtains nucleotide binding competence when incorporated into an oligomeric complex. These findings support a model in which nucleotide binding to NBD1 is dependent on and regulated by trans-acting elements from neighboring subunits, either by direct interaction with the nucleotide or by stabilization of a nucleotide binding-competent state. In this way, they provide a basis for intersubunit communication within the functional ClpB complex.

Design and engineering of artificial metalloproteins: from de novo metal coordination to catalysis.

Metalloproteins are essential to sustain life. Natural evolution optimized them for intricate structural, regulatory and catalytic functions that cannot be fulfilled by either a protein or a metal ion alone. In order to understand this synergy and the complex design principles behind the natural systems, simpler mimics were engineered from the bottom up by installing de novo metal sites in either natural or fully designed, artificial protein scaffolds. This review focuses on key challenges associated with this approach. We discuss how proteins can be equipped with binding sites that provide an optimal coordination environment for a metal cofactor of choice, which can be a single metal ion or a complex multinuclear cluster. Furthermore, we highlight recent studies in which artificial metalloproteins were engineered towards new functions, including electron transfer and catalysis. In this context, the powerful combination of de novo protein design and directed evolution is emphasized for metalloenzyme development.

Americium preferred: lanmodulin, a natural lanthanide-binding protein favors an actinide over lanthanides.

The separation and recycling of lanthanides is an active area of research with a growing demand that calls for more environmentally friendly lanthanide sources. Likewise, the efficient and industrial separation of lanthanides from the minor actinides (Np, Am-Fm) is one of the key questions for closing the nuclear fuel cycle; reducing costs and increasing safety. With the advent of the field of lanthanide-dependent bacterial metabolism, bio-inspired applications are in reach. Here, we utilize the natural lanthanide chelator lanmodulin and the luminescent probes Eu3+ and Cm3+ to investigate the inter-metal competition behavior of all lanthanides (except Pm) and the major actinide plutonium as well as three minor actinides neptunium, americium and curium to lanmodulin. Using time-resolved laser-induced fluorescence spectroscopy we show that lanmodulin has the highest relative binding affinity to Nd3+ and Eu3+ among the lanthanide series. When equimolar mixtures of Cm3+ and Am3+ are added to lanmodulin, lanmodulin preferentially binds to Am3+ over Cm3+ whilst Nd3+ and Cm3+ bind with similar relative affinity. The results presented show that a natural lanthanide-binding protein can bind a major and various minor actinides with high relative affinity, paving the way to bio-inspired separation applications. In addition, an easy and versatile method was developed, using the fluorescence properties of only two elements, Eu and Cm, for inter-metal competition studies regarding lanthanides and selected actinides and their binding to biological molecules.

Analysis of electrostatic coupling throughout the laboratory evolution of a designed retroaldolase.

The roles of local interactions in the laboratory evolution of a highly active, computationally designed retroaldolase (RA) are examined. Partial Order Optimum Likelihood (POOL) is used to identify catalytically important amino acid interactions in several RA95 enzyme variants. The series RA95.5, RA95.5-5, RA95.5-8, and RA95.5-8F, representing progress along an evolutionary trajectory with increasing activity, is examined. Computed measures of coupling between charged states of residues show that, as evolution proceeds and higher activities are achieved, electrostatic coupling between the biochemically active amino acids and other residues is increased. In silico residue scanning suggests multiple coupling partners for the catalytic lysine K83. The effects of two predicted partners, Y51 and E85, are tested using site-directed mutagenesis and kinetic analysis of the variants Y51F and E85Q. The Y51F variants show decreases in kcat relative to wild type, with the greatest losses observed for the more evolved constructs; they also exhibit significant decreases in kcat /KM across the series. Only modest decreases in kcat /KM are observed for the E85Q variants with little effect on kcat . Computed metrics of the degree of coupling between protonation states rise significantly as evolution proceeds and catalytic turnover rate increases. Specifically, the charge state of the catalytic lysine K83 becomes more strongly coupled to those of other amino acids as the enzyme evolves to a better catalyst.

Lysine acylation using conjugating enzymes for site-specific modification and ubiquitination of recombinant proteins.

Enzymes are powerful tools for protein labelling due to their specificity and mild reaction conditions. Many protocols, however, are restricted to modifications at protein termini, rely on non-peptidic metabolites or require large recognition domains. Here we report a chemoenzymatic method, which we call lysine acylation using conjugating enzymes (LACE), to site-specifically modify folded proteins at internal lysine residues. LACE relies on a minimal genetically encoded tag (four residues) recognized by the E2 small ubiquitin-like modifier-conjugating enzyme Ubc9, and peptide or protein thioesters. Together, this approach obviates the need for E1 and E3 enzymes, enabling isopeptide formation with just Ubc9 in a programmable manner. We demonstrate the utility of LACE by the site-specific attachment of biochemical probes, one-pot dual-labelling in combination with sortase, and the conjugation of wild-type ubiquitin and ISG15 to recombinant proteins.

Bap (Sil1) regulates the molecular chaperone BiP by coupling release of nucleotide and substrate.

BiP is the endoplasmic member of the Hsp70 family. BiP is regulated by several co-chaperones including the nucleotide-exchange factor (NEF) Bap (Sil1 in yeast). Bap is a two-domain protein. The interaction of the Bap C-terminal domain with the BiP ATPase domain is sufficient for its weak NEF activity. However, stimulation of the BiP ATPase activity requires full-length Bap, suggesting a complex interplay of these two factors. Here, single-molecule FRET experiments with mammalian proteins reveal that Bap affects the conformation of both BiP domains, including the lid subdomain, which is important for substrate binding. The largely unstructured Bap N-terminal domain promotes the substrate release from BiP. Thus, Bap is a conformational regulator affecting both nucleotide and substrate interactions. The preferential interaction with BiP in its ADP state places Bap at a late stage of the chaperone cycle, in which it coordinates release of substrate and ADP, thereby resetting BiP for ATP and substrate binding.

Efficient laboratory evolution of computationally designed enzymes with low starting activities using fluorescence-activated droplet sorting.

De novo biocatalysts with non-natural functionality are accessible by computational enzyme design. The catalytic activities obtained for the initial designs are usually low, but can be optimized significantly by directed evolution. Nevertheless, rate accelerations approaching the level of natural enzymes can only be achieved over many rounds of tedious and time-consuming laboratory evolution. In this work, we show that microfluidic-based screening using fluorescence-activated droplet sorting (FADS) is ideally suited for efficient optimization of designed enzymes with low starting activity, essentially straight out of the computer. We chose the designed retro-aldolase RA95.0, which had been previously evolved by conventional microtiter plate screening, as an example and reoptimized it using the microfluidic-based assay. Our results show that FADS is sufficiently sensitive to detect enzyme activities as low as kcat/Km = 0.5 M(-1)s(-1) The ultra-high throughput of this system makes screening of large mutant libraries possible in which clusters of up to five residues are randomized simultaneously. Thus, combinations of beneficial mutations can be identified directly, leading to large jumps in catalytic activity of up to 80-fold within a single round of evolution. By exploring several evolutionary trajectories in parallel, we identify alternative active site arrangements that exhibit comparably enhanced efficiency but opposite enantioselectivity.