De Novo Evolutionary Emergence of a Symmetrical Protein Is Shaped by Folding Constraints.

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 ß-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 tandem duplications 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.

Automated Design of Efficient and Functionally Diverse Enzyme Repertoires.

Substantial improvements in enzyme activity demand multiple mutations at spatially proximal positions in the active site. Such mutations, however, often exhibit unpredictable epistatic (non-additive) effects on activity. Here we describe FuncLib, an automated method for designing multipoint mutations at enzyme active sites using phylogenetic analysis and Rosetta design calculations. We applied FuncLib to two unrelated enzymes, a phosphotriesterase and an acetyl-CoA synthetase. All designs were active, and most showed activity profiles that significantly differed from the wild-type and from one another. Several dozen designs with only 3-6 active-site mutations exhibited 10- to 4,000-fold higher efficiencies with a range of alternative substrates, including hydrolysis of the toxic organophosphate nerve agents soman and cyclosarin and synthesis of butyryl-CoA. FuncLib is implemented as a web server (http://FuncLib.weizmann.ac.il); it circumvents iterative, high-throughput experimental screens and opens the way to designing highly efficient and diverse catalytic repertoires.

Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase.

Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.

PGR5-LIKE PHOTOSYNTHETIC PHENOTYPE1A redox states alleviate photoinhibition during changes in light intensity.

Plants have evolved photosynthetic regulatory mechanisms to maintain homeostasis in response to light changes during diurnal transitions and those caused by passing clouds or by wind. One such adaptation directs photosynthetic electron flow to a cyclic pathway to alleviate excess energy surges. Here, we assign a function to regulatory cysteines of PGR5-like protein 1A (PGRL1A), a constituent of the PROTON GRADIENT REGULATION5 (PGR5)-dependent cyclic electron flow (CEF) pathway. During step increases from darkness to low light intensity in Arabidopsis (Arabidopsis thaliana), the intermolecular disulfide of the PGRL1A 59-kDa complex was reduced transiently within seconds to the 28-kDa form. In contrast, step increases from darkness to high light stimulated a stable, partially reduced redox state in PGRL1A. Mutations of 2 cysteines in PGRL1A, Cys82 and Cys183, resulted in a constitutively pseudo-reduced state. The mutant displayed higher proton motive force (PMF) and nonphotochemical quenching (NPQ) than the wild type (WT) and showed altered donor and acceptor dynamic flow around PSI. These changes were found to correspond with the redox state of PGRL1A. Continuous light regimes did not affect mutant growth compared to the WT. However, under fluctuating regimes of high light, the mutant showed better growth than the WT. In contrast, in fluctuating regimes of low light, the mutant displayed a growth penalty that can be attributed to constant stimulation of CEF under low light. Treatment with photosynthetic inhibitors indicated that PGRL1A redox state control depends on the penultimate Fd redox state. Our results showed that redox state changes in PGRL1A are crucial to optimize photosynthesis.

Automated Structure- and Sequence-Based Design of Proteins for High Bacterial Expression and Stability.

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 ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°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.

Identification and characterization of the key enzyme in the biosynthesis of the neurotoxin ß-ODAP in grass pea.

Grass pea (Lathyrus sativus L.) is a grain legume commonly grown in Asia and Africa for food and forage. It is a highly nutritious and robust crop, capable of surviving both droughts and floods. However, it produces a neurotoxic compound, ß-N-oxalyl-L-α,ß-diaminopropionic acid (ß-ODAP), which can cause a severe neurological disorder when consumed as a primary diet component. While the catalytic activity associated with ß-ODAP formation was demonstrated more than 50 years ago, the enzyme responsible for this activity has not been identified. Here, we report on the identity, activity, 3D structure, and phylogenesis of this enzyme-ß-ODAP synthase (BOS). We show that BOS belongs to the benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/benzoyltransferase, deacetylvindoline 4-O-acetyltransferase superfamily of acyltransferases and is structurally similar to hydroxycinnamoyl transferase. Using molecular docking, we propose a mechanism for its catalytic activity, and using heterologous expression in tobacco leaves (Nicotiana benthamiana), we demonstrate that expression of BOS in the presence of its substrates is sufficient for ß-ODAP production in vivo. The identification of BOS may pave the way toward engineering ß-ODAP-free grass pea cultivars, which are safe for human and animal consumption.

Pore-modulating toxins exploit inherent slow inactivation to block K+ channels.

Voltage-dependent potassium channels (Kvs) gate in response to changes in electrical membrane potential by coupling a voltage-sensing module with a K+-selective pore. Animal toxins targeting Kvs are classified as pore blockers, which physically plug the ion conduction pathway, or as gating modifiers, which disrupt voltage sensor movements. A third group of toxins blocks K+ conduction by an unknown mechanism via binding to the channel turrets. Here, we show that Conkunitzin-S1 (Cs1), a peptide toxin isolated from cone snail venom, binds at the turrets of Kv1.2 and targets a network of hydrogen bonds that govern water access to the peripheral cavities that surround the central pore. The resulting ectopic water flow triggers an asymmetric collapse of the pore by a process resembling that of inherent slow inactivation. Pore modulation by animal toxins exposes the peripheral cavity of K+ channels as a novel pharmacological target and provides a rational framework for drug design.

Enzyme Evolution: An Epistatic Ratchet versus a Smooth Reversible Transition.

Evolutionary trajectories are deemed largely irreversible. In a newly diverged protein, reversion of mutations that led to the functional switch typically results in loss of both the new and the ancestral functions. Nonetheless, evolutionary transitions where reversions are viable have also been described. The structural and mechanistic causes of reversion compatibility versus incompatibility therefore remain unclear. We examined two laboratory evolution trajectories of mammalian paraoxonase-1, a lactonase with promiscuous organophosphate hydrolase (OPH) activity. Both trajectories began with the same active-site mutant, His115Trp, which lost the native lactonase activity and acquired higher OPH activity. A neo-functionalization trajectory amplified the promiscuous OPH activity, whereas the re-functionalization trajectory restored the native activity, thus generating a new lactonase that lacks His115. The His115 revertants of these trajectories indicated opposite trends. Revertants of the neo-functionalization trajectory lost both the evolved OPH and the original lactonase activity. Revertants of the trajectory that restored the original lactonase function were, however, fully active. Crystal structures and molecular simulations show that in the newly diverged OPH, the reverted His115 and other catalytic residues are displaced, thus causing loss of both the original and the new activity. In contrast, in the re-functionalization trajectory, reversion compatibility of the original lactonase activity derives from mechanistic versatility whereby multiple residues can fulfill the same task. This versatility enables unique sequence-reversible compositions that are inaccessible when the active site was repurposed toward a new function.

Correction: Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces.

[This corrects the article DOI: 10.1371/journal.pcbi.1007207.].

Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces.

Antibodies developed for research and clinical applications may exhibit suboptimal stability, expressibility, or affinity. Existing optimization strategies focus on surface mutations, whereas natural affinity maturation also introduces mutations in the antibody core, simultaneously improving stability and affinity. To systematically map the mutational tolerance of an antibody variable fragment (Fv), we performed yeast display and applied deep mutational scanning to an anti-lysozyme antibody and found that many of the affinity-enhancing mutations clustered at the variable light-heavy chain interface, within the antibody core. Rosetta design combined enhancing mutations, yielding a variant with tenfold higher affinity and substantially improved stability. To make this approach broadly accessible, we developed AbLIFT, an automated web server that designs multipoint core mutations to improve contacts between specific Fv light and heavy chains (http://AbLIFT.weizmann.ac.il). We applied AbLIFT to two unrelated antibodies targeting the human antigens VEGF and QSOX1. Strikingly, the designs improved stability, affinity, and expression yields. The results provide proof-of-principle for bypassing laborious cycles of antibody engineering through automated computational affinity and stability design.

Principles for computational design of binding antibodies.

Natural proteins must both fold into a stable conformation and exert their molecular function. To date, computational design has successfully produced stable and atomically accurate proteins by using so-called "ideal" folds rich in regular secondary structures and almost devoid of loops and destabilizing elements, such as cavities. Molecular function, such as binding and catalysis, however, often demands nonideal features, including large and irregular loops and buried polar interaction networks, which have remained challenging for fold design. Through five design/experiment cycles, we learned principles for designing stable and functional antibody variable fragments (Fvs). Specifically, we (i) used sequence-design constraints derived from antibody multiple-sequence alignments, and (ii) during backbone design, maintained stabilizing interactions observed in natural antibodies between the framework and loops of complementarity-determining regions (CDRs) 1 and 2. Designed Fvs bound their ligands with midnanomolar affinities and were as stable as natural antibodies, despite having >30 mutations from mammalian antibody germlines. Furthermore, crystallographic analysis demonstrated atomic accuracy throughout the framework and in four of six CDRs in one design and atomic accuracy in the entire Fv in another. The principles we learned are general, and can be implemented to design other nonideal folds, generating stable, specific, and precise antibodies and enzymes.

A de novo protein binding pair by computational design and directed evolution.

The de novo design of protein-protein interfaces is a stringent test of our understanding of the principles underlying protein-protein interactions and would enable unique approaches to biological and medical challenges. Here we describe a motif-based method to computationally design protein-protein complexes with native-like interface composition and interaction density. Using this method we designed a pair of proteins, Prb and Pdar, that heterodimerize with a Kd of 130 nM, 1000-fold tighter than any previously designed de novo protein-protein complex. Directed evolution identified two point mutations that improve affinity to 180 pM. Crystal structures of an affinity-matured complex reveal binding is entirely through the designed interface residues. Surprisingly, in the in vitro evolved complex one of the partners is rotated 180° relative to the original design model, yet still maintains the central computationally designed hotspot interaction and preserves the character of many peripheral interactions. This work demonstrates that high-affinity protein interfaces can be created by designing complementary interaction surfaces on two noninteracting partners and underscores remaining challenges.

Triple-Stage Mass Spectrometry Unravels the Heterogeneity of an Endogenous Protein Complex.

Protein complexes often represent an ensemble of different assemblies with distinct functions and regulation. This increased complexity is enabled by the variety of protein diversification mechanisms that exist at every step of the protein biosynthesis pathway, such as alternative splicing and post transcriptional and translational modifications. The resulting variation in subunits can generate compositionally distinct protein assemblies. These different forms of a single protein complex may comprise functional variances that enable response and adaptation to varying cellular conditions. Despite the biological importance of this layer of complexity, relatively little is known about the compositional heterogeneity of protein complexes, mostly due to technical barriers of studying such closely related species. Here, we show that native mass spectrometry (MS) offers a way to unravel this inherent heterogeneity of protein assemblies. Our approach relies on the advanced Orbitrap mass spectrometer capable of multistage MS analysis across all levels of protein organization. Specifically, we have implemented a two-step fragmentation process in the inject flatapole device, which was converted to a linear ion trap, and can now probe the intact protein complex assembly, through its constituent subunits, to the primary sequence of each protein. We demonstrate our approach on the yeast homotetrameric FBP1 complex, the rate-limiting enzyme in gluconeogenesis. We show that the complex responds differently to changes in growth conditions by tuning phosphorylation dynamics. Our methodology deciphers, on a single instrument and in a single measurement, the stoichiometry, kinetics, and exact position of modifications, contributing to the exposure of the multilevel diversity of protein complexes.

Fine Tuning of Chlorophyll Spectra by Protein-Induced Ring Deformation.

The ability to tune the light-absorption properties of chlorophylls by their protein environment is the key to the robustness and high efficiency of photosynthetic light-harvesting proteins. Unfortunately, the intricacy of the natural complexes makes it very difficult to identify and isolate specific protein-pigment interactions that underlie the spectral-tuning mechanisms. Herein we identify and demonstrate the tuning mechanism of chlorophyll spectra in type II water-soluble chlorophyll binding proteins from Brassicaceae (WSCPs). By comparing the molecular structures of two natural WSCPs we correlate a shift in the chlorophyll red absorption band with deformation of its tetrapyrrole macrocycle that is induced by changing the position of a nearby tryptophan residue. We show by a set of reciprocal point mutations that this change accounts for up to 2/3 of the observed spectral shift between the two natural variants.

AcpM, the meromycolate extension acyl carrier protein of Mycobacterium tuberculosis, is activated by the 4'-phosphopantetheinyl transferase PptT, a potential target of the multistep mycolic acid biosynthesis.

Modification of acyl carrier proteins (ACP) or domains by the covalent binding of a 4'-phosphopantetheine (4'-PP) moiety is a fundamental condition for activation of fatty acid synthases (FASes) and polyketide synthases (PKSes). Binding of 4'-PP is mediated by 4' phosphopantetheinyl transfersases (PPTases). Mycobacterium tuberculosis (Mtb) possesses two essential PPTases: acyl carrier protein synthase (Mtb AcpS), which activates the multidomain fatty acid synthase I (FAS I), and Mtb PptT, an Sfp-type broad spectrum PPTase that activates PKSes. To date, it has not been determined which of the two Mtb PPTases, AcpS or PptT, activates the meromycolate extension ACP, Mtb AcpM, en route to the production of mycolic acids, the main components of the mycobacterial cell wall. In this study, we tested the enzymatic activation of a highly purified Mtb apo-AcpM to Mtb holo-AcpM by either Mtb PptT or Mtb AcpS. By using SDS-PAGE band shift assay and mass spectrometry analysis, we found that Mtb PptT is the PPTase that activates Mtb AcpM. We measured the catalytic activity of Mtb PptT toward CoA, using an activation assay of a blue pigment synthase, BpsA (a nonribosomal peptide synthase, NRPS). BpsA activation by Mtb PptT was inhibited by Mtb apo-AcpM through competition for CoA, in accord with Mtb AcpM activation. A structural model of the putative interaction between Mtb PptT and Mtb AcpM suggests that both hydrophobic and electrostatic interactions stabilize this complex. To conclude, activation of Mtb AcpM by Mtb PptT reveals a potential target of the multistep mycolic acid biosynthesis.

Bridging the gaps in design methodologies by evolutionary optimization of the stability and proficiency of designed Kemp eliminase KE59.

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 ≤ 5 10(4) M(-1)s(-1)). Here we describe the optimization of a third design, KE59. Although KE59 was the most catalytically efficient Kemp eliminase from this design series (by k(cat)/K(m), and by catalyzing the elimination of nonactivated benzisoxazoles), its impaired stability prevented its evolutionary optimization. To boost KE59's evolvability, stabilizing consensus mutations were included in the libraries throughout the directed evolution process. The libraries were also screened with less activated substrates. Sixteen rounds of mutation and selection led to > 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 10(6) M(-1)s(-1), and k(cat)/k(uncat) values of ≤ 10(7) almost regardless of substrate reactivity. Biochemical, structural, and molecular dynamics (MD) simulation studies provided insights regarding the optimization of KE59. Overall, the directed evolution of three different designed Kemp eliminases, KE07, KE70, and KE59, demonstrates that computational designs are highly evolvable and can be optimized to high catalytic efficiencies.

Kemp elimination catalysts by computational enzyme design.

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.

VWD domain stabilization by autocatalytic Asp-Pro cleavage.

Domains known as von Willebrand factor type D (VWD) are found in extracellular and cell-surface proteins including von Willebrand factor, mucins, and various signaling molecules and receptors. Many VWD domains have a glycine-aspartate-proline-histidine (GDPH) amino-acid sequence motif, which is hydrolytically cleaved post-translationally between the aspartate (Asp) and proline (Pro). The Fc IgG binding protein (FCGBP), found in intestinal mucus secretions and other extracellular environments, contains 13 VWD domains, 11 of which have a GDPH cleavage site. In this study, we investigated the structural and biophysical consequences of Asp-Pro peptide cleavage in a representative FCGBP VWD domain. We found that endogenous Asp-Pro cleavage increases the resistance of the domain to exogenous proteolytic degradation. Tertiary structural interactions made by the newly generated chain termini, as revealed by a crystal structure of an FCGBP segment containing the VWD domain, may explain this observation. Notably, the Gly-Asp peptide bond, upstream of the cleavage site, assumed the cis configuration in the structure. In addition to these local features of the cleavage site, a global organizational difference was seen when comparing the FCGBP segment structure with the numerous other structures containing the same set of domains. Together, these data illuminate the outcome of GDPH cleavage and demonstrate the plasticity of proteins with VWD domains, which may contribute to their evolution for function in a dynamic extracellular environment.

Structure of ALD1, a plant-specific homologue of the universal diaminopimelate aminotransferase enzyme of lysine biosynthesis.

Diaminopimelate aminotransferase (DAP-AT) is an enzyme in the lysine-biosynthesis pathway. Conversely, ALD1, a close homologue of DAP-AT in plants, uses lysine as a substrate in vitro. Both proteins require pyridoxal-5'-phosphate (PLP) for their activity. The structure of ALD1 from the flowering plant Arabidopsis thaliana (AtALD1) was solved at a resolution of 2.3 Å. Comparison of AtALD1 with the previously solved structure of A. thaliana DAP-AT (AtDAP-AT) revealed similar interactions with PLP despite sequence differences within the PLP-binding site. However, sequence differences between the binding site of AtDAP-AT for malate, a purported mimic of substrate binding, and the corresponding site in AtALD1 led to different interactions. This suggests that either the substrate itself, or the substrate-binding mode, differs in the two proteins, supporting the known in vitro findings.