Computational optimization of antibody humanness and stability by systematic energy-based ranking.

Conventional methods for humanizing animal-derived antibodies involve grafting their complementarity-determining regions onto homologous human framework regions. However, this process can substantially lower antibody stability and antigen-binding affinity, and requires iterative mutational fine-tuning to recover the original antibody properties. Here we report a computational method for the systematic grafting of animal complementarity-determining regions onto thousands of human frameworks. The method, which we named CUMAb (for computational human antibody design; available at http://CUMAb.weizmann.ac.il ), starts from an experimental or model antibody structure and uses Rosetta atomistic simulations to select designs by energy and structural integrity. CUMAb-designed humanized versions of five antibodies exhibited similar affinities to those of the parental animal antibodies, with some designs showing marked improvement in stability. We also show that (1) non-homologous frameworks are often preferred to highest-homology frameworks, and (2) several CUMAb designs that differ by dozens of mutations and that use different human frameworks are functionally equivalent.

Stable Mammalian Serum Albumins Designed for Bacterial Expression.

Albumin is the most abundant protein in the blood serum of mammals and has essential carrier and physiological roles. Albumins are also used in a wide variety of molecular and cellular experiments and in the cultivated meat industry. Despite their importance, however, albumins are challenging for heterologous expression in microbial hosts, likely due to 17 conserved intramolecular disulfide bonds. Therefore, albumins used in research and biotechnological applications either derive from animal serum, despite severe ethical and reproducibility concerns, or from recombinant expression in yeast or rice. We use the PROSS algorithm to stabilize human and bovine serum albumins, finding that all are highly expressed in E. coli. Design accuracy is verified by crystallographic analysis of a human albumin variant with 16 mutations. This albumin variant exhibits ligand binding properties similar to those of the wild type. Remarkably, a design with 73 mutations relative to human albumin exhibits over 40 °C improved stability and is stable beyond the boiling point of water. Our results suggest that proteins with many disulfide bridges have the potential to exhibit extreme stability when subjected to design. The designed albumins may be used to make economical, reproducible, and animal-free reagents for molecular and cell biology. They also open the way to high-throughput screening to study and enhance albumin carrier properties.

Evolutionary paths that link orthogonal pairs of binding proteins.

Some protein binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if their affinity exceeds the threshold required for function1-4. Thus, homologous and high-specificity binding pairs bring to light an evolutionary conundrum: how does a new specificity evolve while maintaining the required affinity in each intermediate5,6? Until now, a fully functional single-mutation path that connects two orthogonal pairs has only been described where the pairs were mutationally close thus enabling experimental enumeration of all intermediates2. We present an atomistic and graph-theoretical framework for discovering low molecular strain single-mutation paths that connect two extant pairs, enabling enumeration beyond experimental capability. We apply it to two orthogonal bacterial colicin endonuclease-immunity pairs separated by 17 interface mutations7. We were not able to find a strain-free and functional path in the sequence space defined by the two extant pairs. But including mutations that bridge amino acids that cannot be exchanged through single-nucleotide mutations led us to a strain-free 19-mutation trajectory that is completely viable in vivo. Our experiments show that the specificity switch is remarkably abrupt, resulting from only one radical mutation on each partner. Furthermore, each of the critical specificity-switch mutations increases fitness, demonstrating that functional divergence could be driven by positive Darwinian selection. These results reveal how even radical functional changes in an epistatic fitness landscape may evolve.

Allosteric regulation of the 20S proteasome by the Catalytic Core Regulators (CCRs) family.

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

Designed active-site library reveals thousands of functional GFP variants.

Mutations in a protein active site can lead to dramatic and useful changes in protein activity. The active site, however, is sensitive to mutations due to a high density of molecular interactions, substantially reducing the likelihood of obtaining functional multipoint mutants. We introduce an atomistic and machine-learning-based approach, called high-throughput Functional Libraries (htFuncLib), that designs a sequence space in which mutations form low-energy combinations that mitigate the risk of incompatible interactions. We apply htFuncLib to the GFP chromophore-binding pocket, and, using fluorescence readout, recover >16,000 unique designs encoding as many as eight active-site mutations. Many designs exhibit substantial and useful diversity in functional thermostability (up to 96 °C), fluorescence lifetime, and quantum yield. By eliminating incompatible active-site mutations, htFuncLib generates a large diversity of functional sequences. We envision that htFuncLib will be used in one-shot optimization of activity in enzymes, binders, and other proteins.

Computational design and molecular dynamics simulations suggest the mode of substrate binding in ceramide synthases.

Until now, membrane-protein stabilization has relied on iterations of mutations and screening. We now validate a one-step algorithm, mPROSS, for stabilizing membrane proteins directly from an AlphaFold2 model structure. Applied to the lipid-generating enzyme, ceramide synthase, 37 designed mutations lead to a more stable form of human CerS2. Together with molecular dynamics simulations, we propose a pathway by which substrates might be delivered to the ceramide synthases.

Computational design of BclxL inhibitors that target transmembrane domain interactions.

Several methods have been developed to explore interactions among water-soluble proteins or regions of proteins. However, techniques to target transmembrane domains (TMDs) have not been examined thoroughly despite their importance. Here, we developed a computational approach to design sequences that specifically modulate protein-protein interactions in the membrane. To illustrate this method, we demonstrated that BclxL can interact with other members of the B cell lymphoma 2 (Bcl2) family through the TMD and that these interactions are required for BclxL control of cell death. Next, we designed sequences that specifically recognize and sequester the TMD of BclxL. Hence, we were able to prevent BclxL intramembrane interactions and cancel its antiapoptotic effect. These results advance our understanding of protein-protein interactions in membranes and provide a means to modulate them. Moreover, the success of our approach may trigger the development of a generation of inhibitors targeting interactions between TMDs.

Design of a stable human acid-ß-glucosidase: towards improved Gaucher disease therapy and mutation classification.

Acid-ß-glucosidase (GCase, EC3.2.1.45), the lysosomal enzyme which hydrolyzes the simple glycosphingolipid, glucosylceramide (GlcCer), is encoded by the GBA1 gene. Biallelic mutations in GBA1 cause the human inherited metabolic disorder, Gaucher disease (GD), in which GlcCer accumulates, while heterozygous GBA1 mutations are the highest genetic risk factor for Parkinson's disease (PD). Recombinant GCase (e.g., Cerezyme® ) is produced for use in enzyme replacement therapy for GD and is largely successful in relieving disease symptoms, except for the neurological symptoms observed in a subset of patients. As a first step toward developing an alternative to the recombinant human enzymes used to treat GD, we applied the PROSS stability-design algorithm to generate GCase variants with enhanced stability. One of the designs, containing 55 mutations compared to wild-type human GCase, exhibits improved secretion and thermal stability. Furthermore, the design has higher enzymatic activity than the clinically used human enzyme when incorporated into an AAV vector, resulting in a larger decrease in the accumulation of lipid substrates in cultured cells. Based on stability-design calculations, we also developed a machine learning-based approach to distinguish benign from deleterious (i.e., disease-causing) GBA1 mutations. This approach gave remarkably accurate predictions of the enzymatic activity of single-nucleotide polymorphisms in the GBA1 gene that are not currently associated with GD or PD. This latter approach could be applied to other diseases to determine risk factors in patients carrying rare mutations.

Repertoire of Computationally Designed Peroxygenases for Enantiodivergent C-H Oxyfunctionalization Reactions.

The generation of enantiodivergent biocatalysts for C-H oxyfunctionalizations is ever more important in modern synthetic chemistry. Here, we have applied the FuncLib algorithm based on phylogenetic and Rosetta calculations to design a diverse repertoire of active, stable, and enantiodivergent fungal peroxygenases. 24 designs, each carrying 4-5 mutations in the catalytic core, were expressed functionally in yeast and benchmarked against characteristic model compounds. Several designs were active and stable in a range of temperature and pH, displaying unprecedented enantiodivergence, changing regioselectivity from alkyl to aromatic hydroxylation, and increasing catalytic efficiencies up to 10-fold, with 15-fold improvements in total turnover numbers over the parental enzyme. We find that this dramatic functional divergence stems from beneficial epistasis among the mutations and an extensive reorganization of the heme channel. Our work demonstrates that FuncLib can rapidly design highly functional libraries enriched in enantioselective peroxygenases not seen in nature for a range of biotechnological applications.

Protein quaternary structures in solution are a mixture of multiple forms.

Over half the proteins in the E. coli cytoplasm form homo or hetero-oligomeric structures. Experimentally determined structures are often considered in determining a protein's oligomeric state, but static structures miss the dynamic equilibrium between different quaternary forms. The problem is exacerbated in homo-oligomers, where the oligomeric states are challenging to characterize. Here, we re-evaluated the oligomeric state of 17 different bacterial proteins across a broad range of protein concentrations and solutions by native mass spectrometry (MS), mass photometry (MP), size exclusion chromatography (SEC), and small-angle X-ray scattering (SAXS), finding that most exhibit several oligomeric states. Surprisingly, some proteins did not show mass-action driven equilibrium between the oligomeric states. For approximately half the proteins, the predicted oligomeric forms described in publicly available databases underestimated the complexity of protein quaternary structures in solution. Conversely, AlphaFold multimer provided an accurate description of the potential multimeric states for most proteins, suggesting that it could help resolve uncertainties on the solution state of many proteins.

Designed High-Redox Potential Laccases Exhibit High Functional Diversity.

White-rot fungi secrete an impressive repertoire of high-redox potential laccases (HRPLs) and peroxidases for efficient oxidation and utilization of lignin. Laccases are attractive enzymes for the chemical industry due to their broad substrate range and low environmental impact. Since expression of functional recombinant HRPLs is challenging, however, iterative-directed evolution protocols have been applied to improve their expression, activity, and stability. We implement a rational, stabilize-and-diversify strategy to two HRPLs that we could not functionally express. First, we use the PROSS stability-design algorithm to allow functional expression in yeast. Second, we use the stabilized enzymes as starting points for FuncLib active-site design to improve their activity and substrate diversity. Four of the FuncLib-designed HRPLs and their PROSS progenitor exhibit substantial diversity in reactivity profiles against high-redox potential substrates, including lignin monomers. Combinations of 3-4 subtle mutations that change the polarity, solvation, and sterics of the substrate-oxidation site result in orders of magnitude changes in reactivity profiles. These stable and versatile HRPLs are a step toward generating an effective lignin-degrading consortium of enzymes that can be secreted from yeast. The stabilize-and-diversify strategy can be applied to other challenging enzyme families to study and expand the utility of natural enzymes.

Computationally designed hyperactive Cas9 enzymes.

The ability to alter the genomes of living cells is key to understanding how genes influence the functions of organisms and will be critical to modify living systems for useful purposes. However, this promise has long been limited by the technical challenges involved in genetic engineering. Recent advances in gene editing have bypassed some of these challenges but they are still far from ideal. Here we use FuncLib to computationally design Cas9 enzymes with substantially higher donor-independent editing activities. We use genetic circuits linked to cell survival in yeast to quantify Cas9 activity and discover synergistic interactions between engineered regions. These hyperactive Cas9 variants function efficiently in mammalian cells and introduce larger and more diverse pools of insertions and deletions into targeted genomic regions, providing tools to enhance and expand the possible applications of CRISPR-based gene editing.

De novo-designed transmembrane domains tune engineered receptor functions.

De novo-designed receptor transmembrane domains (TMDs) present opportunities for precise control of cellular receptor functions. We developed a de novo design strategy for generating programmed membrane proteins (proMPs): single-pass α-helical TMDs that self-assemble through computationally defined and crystallographically validated interfaces. We used these proMPs to program specific oligomeric interactions into a chimeric antigen receptor (CAR) that we expressed in mouse primary T cells and found that both in vitro CAR T cell cytokine release and in vivo antitumor activity scaled linearly with the oligomeric state encoded by the receptor TMD, from monomers up to tetramers. All programmed CARs stimulated substantially lower T cell cytokine release relative to the commonly used CD28 TMD, which we show elevated cytokine release through lateral recruitment of the endogenous T cell costimulatory receptor CD28. Precise design using orthogonal and modular TMDs thus provides a new way to program receptor structure and predictably tune activity for basic or applied synthetic biology.

Computer-aided engineering of staphylokinase toward enhanced affinity and selectivity for plasmin.

Cardio- and cerebrovascular diseases are leading causes of death and disability, resulting in one of the highest socio-economic burdens of any disease type. The discovery of bacterial and human plasminogen activators and their use as thrombolytic drugs have revolutionized treatment of these pathologies. Fibrin-specific agents have an advantage over non-specific factors because of lower rates of deleterious side effects. Specifically, staphylokinase (SAK) is a pharmacologically attractive indirect plasminogen activator protein of bacterial origin that forms stoichiometric noncovalent complexes with plasmin, promoting the conversion of plasminogen into plasmin. Here we report a computer-assisted re-design of the molecular surface of SAK to increase its affinity for plasmin. A set of computationally designed SAK mutants was produced recombinantly and biochemically characterized. Screening revealed a pharmacologically interesting SAK mutant with ∼7-fold enhanced affinity toward plasmin, ∼10-fold improved plasmin selectivity and moderately higher plasmin-generating efficiency in vitro. Collectively, the results obtained provide a framework for SAK engineering using computational affinity-design that could pave the way to next-generation of effective, highly selective, and less toxic thrombolytics.

Stabilization of the SARS-CoV-2 receptor binding domain by protein core redesign and deep mutational scanning.

Stabilizing antigenic proteins as vaccine immunogens or diagnostic reagents is a stringent case of protein engineering and design as the exterior surface must maintain recognition by receptor(s) and antigen-specific antibodies at multiple distinct epitopes. This is a challenge, as stability enhancing mutations must be focused on the protein core, whereas successful computational stabilization algorithms typically select mutations at solvent-facing positions. In this study, we report the stabilization of SARS-CoV-2 Wuhan Hu-1 Spike receptor binding domain using a combination of deep mutational scanning and computational design, including the FuncLib algorithm. Our most successful design encodes I358F, Y365W, T430I, and I513L receptor binding domain mutations, maintains recognition by the receptor ACE2 and a panel of different anti-receptor binding domain monoclonal antibodies, is between 1 and 2°C more thermally stable than the original receptor binding domain using a thermal shift assay, and is less proteolytically sensitive to chymotrypsin and thermolysin than the original receptor binding domain. Our approach could be applied to the computational stabilization of a wide range of proteins without requiring detailed knowledge of active sites or binding epitopes. We envision that this strategy may be particularly powerful for cases when there are multiple or unknown binding sites.

Highly Specific Monoclonal Antibody Targeting the Botulinum Neurotoxin Type E Exposed SNAP-25 Neoepitope.

Botulinum neurotoxin type E (BoNT/E), the fastest acting toxin of all BoNTs, cleaves the 25 kDa synaptosomal-associated protein (SNAP-25) in motor neurons, leading to flaccid paralysis. The specific detection and quantification of the BoNT/E-cleaved SNAP-25 neoepitope can facilitate the development of cell-based assays for the characterization of anti-BoNT/E antibody preparations. In order to isolate highly specific monoclonal antibodies suitable for the in vitro immuno-detection of the exposed neoepitope, mice and rabbits were immunized with an eight amino acid peptide composed of the C-terminus of the cleaved SNAP-25. The immunized rabbits developed a specific and robust polyclonal antibody response, whereas the immunized mice mostly demonstrated a weak antibody response that could not discriminate between the two forms of SNAP-25. An immune scFv phage-display library was constructed from the immunized rabbits and a panel of antibodies was isolated. The sequence alignment of the isolated clones revealed high similarity between both heavy and light chains with exceptionally short HCDR3 sequences. A chimeric scFv-Fc antibody was further expressed and characterized, exhibiting a selective, ultra-high affinity (pM) towards the SNAP-25 neoepitope. Moreover, this antibody enabled the sensitive detection of cleaved SNAP-25 in BoNT/E treated SiMa cells with no cross reactivity with the intact SNAP-25. Thus, by applying an immunization and selection procedure, we have isolated a novel, specific and high-affinity antibody against the BoNT/E-derived SNAP-25 neoepitope. This novel antibody can be applied in in vitro assays that determine the potency of antitoxin preparations and reduce the use of laboratory animals for these purposes.

Computationally designed dual-color MRI reporters for noninvasive imaging of transgene expression.

Imaging of gene-expression patterns in live animals is difficult to achieve with fluorescent proteins because tissues are opaque to visible light. Imaging of transgene expression with magnetic resonance imaging (MRI), which penetrates to deep tissues, has been limited by single reporter visualization capabilities. Moreover, the low-throughput capacity of MRI limits large-scale mutagenesis strategies to improve existing reporters. Here we develop an MRI system, called GeneREFORM, comprising orthogonal reporters for two-color imaging of transgene expression in deep tissues. Starting from two promiscuous deoxyribonucleoside kinases, we computationally designed highly active, orthogonal enzymes ('reporter genes') that specifically phosphorylate two MRI-detectable synthetic deoxyribonucleosides ('reporter probes'). Systemically administered reporter probes exclusively accumulate in cells expressing the designed reporter genes, and their distribution is displayed as pseudo-colored MRI maps based on dynamic proton exchange for noninvasive visualization of transgene expression. We envision that future extensions of GeneREFORM will pave the way to multiplexed deep-tissue mapping of gene expression in live animals.

Structure and receptor recognition by the Lassa virus spike complex.

Lassa virus (LASV) is a human pathogen, causing substantial morbidity and mortality1,2. Similar to other Arenaviridae, it presents a class-I spike complex on its surface that facilitates cell entry. The virus's cellular receptor is matriglycan, a linear carbohydrate that is present on α-dystroglycan3,4, but the molecular mechanism that LASV uses to recognize this glycan is unknown. In addition, LASV and other arenaviruses have a unique signal peptide that forms an integral and functionally important part of the mature spike5-8; yet the structure, function and topology of the signal peptide in the membrane remain uncertain9-11. Here we solve the structure of a complete native LASV spike complex, finding that the signal peptide crosses the membrane once and that its amino terminus is located in the extracellular region. Together with a double-sided domain-switching mechanism, the signal peptide helps to stabilize the spike complex in its native conformation. This structure reveals that the LASV spike complex is preloaded with matriglycan, suggesting the mechanism of binding and rationalizing receptor recognition by α-dystroglycan-tropic arenaviruses. This discovery further informs us about the mechanism of viral egress and may facilitate the rational design of novel therapeutics that exploit this binding site.

Stable and Functionally Diverse Versatile Peroxidases Designed Directly from Sequences.

White-rot fungi secrete a repertoire of high-redox potential oxidoreductases to efficiently decompose lignin. Of these enzymes, versatile peroxidases (VPs) are the most promiscuous biocatalysts. VPs are attractive enzymes for research and industrial use but their recombinant production is extremely challenging. To date, only a single VP has been structurally characterized and optimized for recombinant functional expression, stability, and activity. Computational enzyme optimization methods can be applied to many enzymes in parallel but they require accurate structures. Here, we demonstrate that model structures computed by deep-learning-based ab initio structure prediction methods are reliable starting points for one-shot PROSS stability-design calculations. Four designed VPs encoding as many as 43 mutations relative to the wildtype enzymes are functionally expressed in yeast, whereas their wildtype parents are not. Three of these designs exhibit substantial and useful diversity in their reactivity profiles and tolerance to environmental conditions. The reliability of the new generation of structure predictors and design methods increases the scale and scope of computational enzyme optimization, enabling efficient discovery and exploitation of the functional diversity in natural enzyme families directly from genomic databases.