Highly active rubiscos discovered by systematic interrogation of natural sequence diversity.

CO2 is converted into biomass almost solely by the enzyme rubisco. The poor carboxylation properties of plant rubiscos have led to efforts that made it the most kinetically characterized enzyme, yet these studies focused on < 5% of its natural diversity. Here, we searched for fast-carboxylating variants by systematically mining genomic and metagenomic data. Approximately 33,000 unique rubisco sequences were identified and clustered into ≈ 1,000 similarity groups. We then synthesized, purified, and biochemically tested the carboxylation rates of 143 representatives, spanning all clusters of form-II and form-II/III rubiscos. Most variants (> 100) were active in vitro, with the fastest having a turnover number of 22 ± 1 s-1 -sixfold faster than the median plant rubisco and nearly twofold faster than the fastest measured rubisco to date. Unlike rubiscos from plants and cyanobacteria, the fastest variants discovered here are homodimers and exhibit a much simpler folding and activation kinetics. Our pipeline can be utilized to explore the kinetic space of other enzymes of interest, allowing us to get a better view of the biosynthetic potential of the biosphere.

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.

Communication between viruses guides lysis-lysogeny decisions.

Temperate viruses can become dormant in their host cells, a process called lysogeny. In every infection, such viruses decide between the lytic and the lysogenic cycles, that is, whether to replicate and lyse their host or to lysogenize and keep the host viable. Here we show that viruses (phages) of the SPbeta group use a small-molecule communication system to coordinate lysis-lysogeny decisions. During infection of its Bacillus host cell, the phage produces a six amino-acids-long communication peptide that is released into the medium. In subsequent infections, progeny phages measure the concentration of this peptide and lysogenize if the concentration is sufficiently high. We found that different phages encode different versions of the communication peptide, demonstrating a phage-specific peptide communication code for lysogeny decisions. We term this communication system the 'arbitrium' system, and further show that it is encoded by three phage genes: aimP, which produces the peptide; aimR, the intracellular peptide receptor; and aimX, a negative regulator of lysogeny. The arbitrium system enables a descendant phage to 'communicate' with its predecessors, that is, to estimate the amount of recent previous infections and hence decide whether to employ the lytic or lysogenic cycle.

Two closed ATP- and ADP-dependent conformations in yeast Hsp90 chaperone detected by Mn(II) EPR spectroscopic techniques.

Hsp90 plays a central role in cell homeostasis by assisting folding and maturation of a large variety of clients. It is a homo-dimer, which functions via hydrolysis of ATP-coupled to conformational changes. Hsp90's conformational cycle in the absence of cochaperones is currently postulated as apo-Hsp90 being an ensemble of "open"/"closed" conformations. Upon ATP binding, Hsp90 adopts an active ATP-bound closed conformation where the N-terminal domains, which comprise the ATP binding site, are in close contact. However, there is no consensus regarding the conformation of the ADP-bound Hsp90, which is considered important for client release. In this work, we tracked the conformational states of yeast Hsp90 at various stages of ATP hydrolysis in frozen solutions employing electron paramagnetic resonance (EPR) techniques, particularly double electron-electron resonance (DEER) distance measurements. Using rigid Gd(III) spin labels, we found the C domains to be dimerized with same distance distribution at all hydrolysis states. Then, we substituted the ATPase Mg(II) cofactor with paramagnetic Mn(II) and followed the hydrolysis state using hyperfine spectroscopy and measured the inter-N-domain distance distributions via Mn(II)-Mn(II) DEER. The point character of the Mn(II) spin label allowed us resolve 2 different closed states: The ATP-bound (prehydrolysis) characterized by a distance distribution having a maximum of 4.3 nm, which broadened and shortened, shifting the mean to 3.8 nm at the ADP-bound state (posthydrolysis). This provides experimental evidence to a second closed conformational state of Hsp90 in solution, referred to as "compact." Finally, the so-called high-energy state, trapped by addition of vanadate, was found structurally similar to the posthydrolysis state.

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

Efficient Targeted Degradation via Reversible and Irreversible Covalent PROTACs.

Proteolysis targeting chimeras (PROTACs) represent an exciting inhibitory modality with many advantages, including substoichiometric degradation of targets. Their scope, though, is still limited to date by the requirement for a sufficiently potent target binder. A solution that proved useful in tackling challenging targets is the use of electrophiles to allow irreversible binding to the target. However, such binding will negate the catalytic nature of PROTACs. Reversible covalent PROTACs potentially offer the best of both worlds. They possess the potency and selectivity associated with the formation of the covalent bond, while being able to dissociate and regenerate once the protein target is degraded. Using Bruton's tyrosine kinase (BTK) as a clinically relevant model system, we show efficient degradation by noncovalent, irreversible covalent, and reversible covalent PROTACs, with <10 nM DC50's and >85% degradation. Our data suggest that part of the degradation by our irreversible covalent PROTACs is driven by reversible binding prior to covalent bond formation, while the reversible covalent PROTACs drive degradation primarily by covalent engagement. The PROTACs showed enhanced inhibition of B cell activation compared to ibrutinib and exhibit potent degradation of BTK in patient-derived primary chronic lymphocytic leukemia cells. The most potent reversible covalent PROTAC, RC-3, exhibited enhanced selectivity toward BTK compared to noncovalent and irreversible covalent PROTACs. These compounds may pave the way for the design of covalent PROTACs for a wide variety of challenging targets.

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.

Tracking Conformational Changes in Calmodulin in vitro, in Cell Extract, and in Cells by Electron Paramagnetic Resonance Distance Measurements.

It is an open question whether the conformations of proteins sampled in dilute solutions are the same as in the cellular environment. Here we address this question by double electron-electron resonance (DEER) distance measurements with Gd(III) spin labels to probe the conformations of calmodulin (CaM) in vitro, in cell extract, and in human HeLa cells. Using the CaM mutants N53C/T110C and T34C/T117C labeled with maleimide-DOTA-Gd(III) in the N- and C-terminal domains, we observed broad and varied interdomain distance distributions. The in vitro distance distributions of apo-CaM and holo-CaM in the presence and absence of the IQ target peptide can be described by combinations of closed, open, and collapsed conformations. In cell extract, apo- and holo-CaM bind to target proteins in a similar way as apo- and holo-CaM bind to IQ peptide in vitro. In HeLa cells, however, in the presence or absence of elevated in-cell Ca2+ levels CaM unexpectedly produced more open conformations and very broad distance distributions indicative of many different interactions with in-cell components. These results show-case the importance of in-cell analyses of protein structures.

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.

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.

A Recombinant Collagen-mRNA Platform for Controllable Protein Synthesis.

We have developed a collagen-mRNA platform for controllable protein production that is intended to be less prone to the problems associated with commonly used mRNA therapy as well as with collagen skin-healing procedures. A collagen mimic was constructed according to a recombinant method and was used as scaffold for translating mRNA chains into proteins. Cysteines were genetically inserted into the collagen chain at positions allowing efficient ribosome translation activity while minimizing mRNA misfolding and degradation. Enhanced green fluorescence protein (eGFP) mRNA bound to collagen was successfully translated by cell-free Escherichia coli ribosomes. This system enabled an accurate control of specific protein synthesis by monitoring expression time and level. Luciferase-mRNA was also translated on collagen scaffold by eukaryotic cell extracts. Thus we have demonstrated the feasibility of controllable protein synthesis on collagen scaffolds by ribosomal machinery.

Sensing Protein Surfaces with Targeted Fluorescent Receptors.

Invited for the cover of this issue is the group of David Margulies at the Weizmann Institute of Science (Israel). The image highlights the analogy between fluorescent molecular sensors and a miniaturized camera that can capture changes that occur at the nanoscale and shed light on the structural state of proteins. Read the full text of the article at 10.1002/chem.201502069.

Sensing Protein Surfaces with Targeted Fluorescent Receptors.

A methodology for creating fluorescent molecular sensors that respond to changes that occur on the surfaces of specific proteins is presented. This approach, which relies on binding cooperatively between a specific His-tag binder and a nonspecific protein-surface receptor, enabled the development of a sensor that can track changes on the surface of a His-tag-labeled calmodulin (His-CaM) upon interacting with metal ions, small molecules, and protein binding partners. The way this approach was used to detect dephosphorylation of an unlabeled calmodulin-dependent protein kinase II (CaMKII), and the binding of Bax BH3 to His-tagged B-cell lymphoma 2 (Bcl-2) protein is also presented.

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.

GTP binding to the ROC domain of DAP-kinase regulates its function through intramolecular signalling.

Death-associated protein kinase (DAPk) was recently suggested by sequence homology to be a member of the ROCO family of proteins. Here, we show that DAPk has a functional ROC (Ras of complex proteins) domain that mediates homo-oligomerization and GTP binding through a defined P-loop motif. Upon binding to GTP, the ROC domain negatively regulates the catalytic activity of DAPk and its cellular effects. Mechanistically, GTP binding enhances an inhibitory autophosphorylation at a distal site that suppresses kinase activity. This study presents a new mechanism of intramolecular signal transduction, by which GTP binding operates in cis to affect the catalytic activity of a distal domain in the protein.