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

Deep mutational scanning for therapeutic antibody engineering

The biophysical and functional properties of monoclonal antibody (mAb) drug candidates are often improved by protein engineering methods to increase the probability of clinical efficacy. One emerging method is deep mutational scanning (DMS) which combines the power of exhaustive protein mutagenesis and functional screening with deep sequencing and bioinformatics. The application of DMS has yielded significant improvements to the affinity, specificity, and stability of several preclinical antibodies alongside novel applications such as introducing multi-specific binding properties. DMS has also been applied directly on target antigens to precisely map antibody-binding epitopes and notably to profile the mutational escape potential of viral targets (e.g., SARS-CoV-2 variants). Finally, DMS combined with machine learning is enabling advances in the computational screening and engineering of therapeutic antibodies.

Direct-MS analysis of antibody-antigen complexes

In recent decades, antibodies (Abs) have attracted the attention of academia and the biopharmaceutical industry due to their therapeutic properties and versatility in binding a vast spectrum of antigens. Different engineering strategies have been developed for optimizing Ab specificity, efficacy, affinity, stability and production, enabling systematic screening and analysis procedures for selecting lead candidates. This quality assessment is critical but usually demands time-consuming and labor-intensive purification procedures. Here, we harnessed the direct-mass spectrometry (direct-MS) approach, in which the analysis is carried out directly from the crude growth media, for the rapid, structural characterization of designed Abs. We demonstrate that properties such as stability, specificity and interactions with antigens can be defined, without the need for prior purification.

Delicate balance among thermal stability, binding affinity, and conformational space explored by single-domain VHH antibodies

The high binding affinities and specificities of antibodies have led to their use as drugs and biosensors. Single-domain V_HH antibodies exhibit high specificity and affinity but have higher stability and solubility than conventional antibodies as they are single-domain proteins. In this work, based on physicochemical measurements and molecular dynamics (MD) simulations, we have gained insight that will facilitate rational design of single-chain V_HH antibodies. We first assessed two homologous V_HH antibodies by differential scanning calorimetry (DSC); one had a high (64.8 °C) and the other a low (58.6 °C) melting temperature. We then generated a series of the variants of the low stability antibody and analyzed their thermal stabilities by DSC and characterized their structures through MD simulations. We found that a single mutation that resulted in 8.2 °C improvement in melting temperature resulted in binding affinity an order of magnitude lower than the parent antibody, likely due to a shift of conformational space explored by the single-chain V_HH antibody. These results suggest that the delicate balance among conformational stability, binding capability, and conformational space explored by antibodies must be considered in design of fully functional single-chain V_HH antibodies.

Binding Characteristic of Various Antibody Formats Against Aflatoxins

The application of recombinant antibodies for the analysis of foods and food contaminants is now a major focus, given their capacity to be engineered to tailor their specificity, enhance their stability, and modify their structural formats to fit the desired analytical platform. In this study, human scFv antibody fragments generated against aflatoxin B1 (AFB1) were selected as the model antibody to explore the effect of antibody formats on their binding activity and to evaluate their potential use as immunoreagents for food contaminant analysis. Four human scFv antibody fragments against aflatoxin B1 (AFB1), previously isolated and engineered by chain shuffling, were converted into various formats, that is, scFv-AP fusions, scFv-Fc, and whole IgG molecules. The result indicated that the effects of the antibody format on the binding property varied, depending on the sequence of scFv. For all of the scFv clones, the scFv-AP fusion format showed the highest sensitivity by competitive ELISA, while the effects on the binding activity after conversion to scFv-Fc or IgG format varied, depending on the amino acid sequence of the antibodies. The sAFH-3e3 antibodies that showed the best performance by competitive ELISA were selected for further investigation. The sAFH-3e3 was converted to the scFv-GFP format and tested by fluorescence-linked immunosorbent assay (FLISA), which showed that its binding property was equivalent to those of scFv-Fc and IgG formats. The potential applications of the sAFH-3e3 in a rapid test kit format based on ELISA (scFv-AP) and in a lateral flow immunochromatography assay (LFIA) (IgG) were demonstrated. A comparison of methods for the extraction of AFB1 from matrices for use with these assay formats indicated that PBS and TBST are better than 70% methanol.

CDR1 Composition Can Affect Nanobody Recombinant Expression Yields

The isolation of nanobodies from pre-immune libraries by means of biopanning is a straightforward process. Nevertheless, the recovered candidates often require optimization to improve some of their biophysical characteristics. In principle, CDRs are not mutated because they are likely to be part of the antibody paratope, but in this work, we describe a mutagenesis strategy that specifically addresses CDR1. Its sequence was identified as an instability hot spot by the PROSS program, and the available structural information indicated that four CDR1 residues bound directly to the antigen. We therefore modified the loop flexibility with the addition of an extra glycine rather than by mutating single amino acids. This approach significantly increased the nanobody yields but traded-off with moderate affinity loss. Accurate modeling coupled with atomistic molecular dynamics simulations enabled the modifications induced by the glycine insertion and the rationale behind the engineering design to be described in detail.

Intrinsic physicochemical profile of marketed antibody-based biotherapeutics

Successful biologic drug discovery and development involves finding functional as well as developable candidates. Once a candidate has been demonstrated to be functional, the next step is to determine whether it can be translated into a drug product. This requires that the candidate can withstand stresses encountered during manufacturing, shipping, and storage. Additionally, it must be safe, efficacious, and possess good pharmacology. In silico analyses of the variable regions of 77 marketed antibody-based biotherapeutics have revealed five nonredundant physicochemical descriptors. Distributions of these descriptors, observed for marketed biotherapeutics, can help prioritize a drug candidate for experimental testing at early discovery stages, guide engineering efforts to further optimize it, and help increase the productivity of biologic drug discovery and development.

Feeding biopharma pipelines with biotherapeutic candidates that possess desirable developability profiles can help improve the productivity of biologic drug discovery and development. Here, we have derived an in silico profile by analyzing computed physicochemical descriptors for the variable regions (Fv) found in 77 marketed antibody-based biotherapeutics. Fv regions of these biotherapeutics demonstrate significant diversities in their germlines, complementarity determining region loop lengths, hydrophobicity, and charge distributions. Furthermore, an analysis of 24 physicochemical descriptors, calculated using homology-based molecular models, has yielded five nonredundant descriptors whose distributions represent stability, isoelectric point, and molecular surface characteristics of their Fv regions. Fv regions of candidates from our internal discovery campaigns, human next-generation sequencing repertoires, and those in clinical-stages (CST) were assessed for similarity with the physicochemical profile derived here. The Fv regions in 33% of CST antibodies show physicochemical properties that are dissimilar to currently marketed biotherapeutics. In comparison, physicochemical characteristics of ∼29% of the Fv regions in human antibodies and ∼27% of our internal hits deviated significantly from those of marketed biotherapeutics. The early availability of this information can help guide hit selection, lead identification, and optimization of biotherapeutic candidates. Insights from this work can also help support portfolio risk assessment, in-licensing, and biopharma collaborations.

Computationally designed pyocyanin demethylase acts synergistically with tobramycin to kill recalcitrant Pseudomonas aeruginosa biofilms

Pseudomonas aeruginosa is an opportunistic human pathogen that develops difficult-to-treat biofilms in immunocompromised individuals, cystic fibrosis patients, and in chronic wounds. P. aeruginosa has an arsenal of physiological attributes that enable it to evade standard antibiotic treatments, particularly in the context of biofilms where it grows slowly and becomes tolerant to many drugs. One of its survival strategies involves the production of the redox-active phenazine, pyocyanin, which promotes biofilm development. We previously identified an enzyme, PodA, that demethylated pyocyanin and disrupted P. aeruginosa biofilm development in vitro. Here, we asked if this protein could be used as a potential therapeutic for P. aeruginosa infections together with tobramycin, an antibiotic typically used in the clinic. A major roadblock to answering this question was the poor yield and stability of wild-type PodA purified from standard Escherichia coli overexpression systems. We hypothesized that the insufficient yields were due to poor packing within PodA's obligatory homotrimeric interfaces. We therefore applied the protein design algorithm, AffiLib, to optimize the symmetric core of this interface, resulting in a design that incorporated five mutations leading to a 20-fold increase in protein yield from heterologous expression and purification and a substantial increase in stability to environmental conditions. The addition of the designed PodA with tobramycin led to increased killing of P. aeruginosa cultures under oxic and hypoxic conditions in both the planktonic and biofilm states. This study highlights the potential for targeting extracellular metabolites to assist the control of P. aeruginosa biofilms that tolerate conventional antibiotic treatment.

Deep geometric representations for modeling effects of mutations on protein-protein binding affinity

Modeling the impact of amino acid mutations on protein-protein interaction plays a crucial role in protein engineering and drug design. In this study, we develop GeoPPI, a novel structure-based deep-learning framework to predict the change of binding affinity upon mutations. Based on the three-dimensional structure of a protein, GeoPPI first learns a geometric representation that encodes topology features of the protein structure via a self-supervised learning scheme. These representations are then used as features for training gradient-boosting trees to predict the changes of protein-protein binding affinity upon mutations. We find that GeoPPI is able to learn meaningful features that characterize interactions between atoms in protein structures. In addition, through extensive experiments, we show that GeoPPI achieves new state-of-the-art performance in predicting the binding affinity changes upon both single- and multi-point mutations on six benchmark datasets. Moreover, we show that GeoPPI can accurately estimate the difference of binding affinities between a few recently identified SARS-CoV-2 antibodies and the receptor-binding domain (RBD) of the S protein. These results demonstrate the potential of GeoPPI as a powerful and useful computational tool in protein design and engineering. Our code and datasets are available at: https://github.com/Liuxg16/GeoPPI.

Estimating the binding affinities of protein-protein interactions (PPIs) is crucial to understand protein function and design new functional proteins. Since the experimental measurement in wet-labs is labor-intensive and time-consuming, fast and accurate in silico approaches have received much attention. Although considerable efforts have been made in this direction, predicting the effects of mutations on the protein-protein binding affinity is still a challenging research problem. In this work, we introduce GeoPPI, a novel computational approach that uses deep geometric representations of protein complexes to predict the effects of mutations on the binding affinity. The geometric representations are first learned via a self-supervised learning scheme and then integrated with gradient-boosting trees to accomplish the prediction. We find that the learned representations encode meaningful patterns underlying the interactions between atoms in protein structures. Also, extensive tests on major benchmark datasets show that GeoPPI has made an important improvement over the existing methods in predicting the effects of mutations on the binding affinity.

Leishmania tarentolae cell-free based approach for rapid anitbody-antigen interaction analysis

Rapid techniques for producing high-quality recombinant proteins are essential for fast protein functional analysis, as well as various screening applications. Cell-free protein expression is an enabling tool in protein research capable of producing high-quality proteins within a few hours. In this chapter, we describe the use of a Leishmania tarentolae-based cell-free expression system to produce antibody fragments coupled to the analysis of their interaction with their ligands. Interaction analysis is performed using the scalable and sensitive AlphaLISA bead proximity assay. The method presented in this chapter offers a rapid and inexpensive approach for production of putative interacting protein pairs, as well as a multiplexable approach for their rapid interaction analysis.

Biomolecular Recognition of the Glycan Neoantigen CA19-9 by Distinct Antibodies

Glycans decorate the cell surface, secreted glycoproteins and glycolipids, and altered glycans are often found in cancers. Despite their high diagnostic and therapeutic potential, however, glycans are polar and flexible molecules that are quite challenging for the development and design of high-affinity binding antibodies. To understand the mechanisms by which glycan neoantigens are specifically recognized by antibodies, we analyze the biomolecular recognition of the tumor-associated carbohydrate antigen CA19-9 by two distinct antibodies using X-ray crystallography. Despite the potential plasticity of glycans and the very different antigen-binding surfaces presented by the antibodies, both structures reveal an essentially identical extended CA19-9 conformer, suggesting that this conformer's stability selects the antibodies. Starting from the bound structure of one of the antibodies, we use the AbLIFT computational algorithm to design a variant with seven core mutations in the variable domain's light-heavy chain interface that exhibits tenfold improved affinity for CA19-9. The results reveal strategies used by antibodies to specifically recognize glycan antigens and show how automated antibody-optimization methods may be used to enhance the clinical potential of existing antibodies.

An antibody engineering platform using amino acid networks: A case study in development of antiviral therapeutics

We present here a case study of an antibody-engineering platform that selects, modifies, and assembles antibody parts to construct novel antibodies. A salient feature of this platform includes the role of amino acid networks in optimizing framework regions (FRs) and complementarity determining regions (CDRs) to engineer new antibodies with desired structure-function relationships. The details of this approach are described in the context of its utility in engineering ZAb_FLEP, a potent anti-Zika virus antibody. ZAb_FLEP comprises of distinct parts, including heavy chain and light chain FRs and CDRs, with engineered features such as loop lengths and optimal epitope-paratope contacts. We demonstrate, with different test antibodies derived from different FR-CDR combinations, that despite these test antibodies sharing high overall sequence similarity, they yield diverse functional readouts. Furthermore, we show that strategies relying on one dimensional sequence similarity-based analyses of antibodies miss important structural nuances of the FR-CDR relationship, which is effectively addressed by the amino acid networks approach of this platform.

Design of a Trispecific Checkpoint Inhibitor and Natural Killer Cell Engager Based on a 2 + 1 Common Light Chain Antibody Architecture

Natural killer cell engagers gained enormous interest in recent years due to their potent anti-tumor activity and favorable safety profile. Simultaneously, chicken-derived antibodies entered clinical studies paving the way for avian-derived therapeutics. In this study, we describe the affinity maturation of a common light chain (cLC)-based, chicken-derived antibody targeting EGFR, followed by utilization of the same light chain for the isolation of CD16a- and PD-L1-specific monoclonal antibodies. The resulting binders target their respective antigen with single-digit nanomolar affinity while blocking the ligand binding of all three respective receptors. Following library-based humanization, bispecific and trispecific variants in a standard 1 + 1 or a 2 + 1 common light chain format were generated, simultaneously targeting EGFR, CD16a, and PD-L1. The trispecific antibody mediated an elevated antibody-dependent cellular cytotoxicity (ADCC) in comparison to the EGFR×CD16a bispecific variant by effectively bridging EGFR/PD-L1 double-positive cancer cells with CD16a-positive effector cells. These findings represent, to our knowledge, the first detailed report on the generation of a trispecific 2 + 1 antibodies exhibiting a common light chain and illustrate synergistic effects of trispecific antigen binding. Overall, this generic procedure paves the way for the engineering of tri- and oligospecific therapeutic antibodies derived from avian immunizations.

UEP: an open-source and fast classifier for predicting the impact of mutations in protein-protein complexes

MOTIVATION

Single protein residue mutations may reshape the binding affinity of protein-protein interactions. Therefore, predicting its effects is of great interest in biotechnology and biomedicine. Unfortunately, the availability of experimental data on binding affinity changes upon mutation is limited, which hampers the development of new and more precise algorithms. Here, we propose UEP, a classifier for predicting beneficial and detrimental mutations in protein-protein complexes trained on interactome data.

RESULTS

Regardless of the simplicity of the UEP algorithm, which is based on a simple three-body contact potential derived from interactome data, we report competitive results with the gold standard methods in this field with the advantage of being faster in terms of computational time. Moreover, we propose a consensus selection procedure by involving the combination of three predictors that showed higher classification accuracy in our benchmark: UEP, pyDock and EvoEF1/FoldX. Overall, we demonstrate that the analysis of interactome data allows predicting the impact of protein-protein mutations using UEP, a fast and reliable open-source code.

AVAILABILITY AND IMPLEMENTATION

UEP algorithm can be found at: https://github.com/pepamengual/UEP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Community-Wide Experimental Evaluation of the PROSS Stability-Design Method

Recent years have seen a dramatic improvement in protein-design methodology. Nevertheless, most methods demand expert intervention, limiting their widespread adoption. By contrast, the PROSS algorithm for improving protein stability and heterologous expression levels has been successfully applied to a range of challenging enzymes and binding proteins. Here, we benchmark the application of PROSS as a stand-alone tool for protein scientists with no or limited experience in modeling. Twelve laboratories from the Protein Production and Purification Partnership in Europe (P4EU) challenged the PROSS algorithm with 14 unrelated protein targets without support from the PROSS developers. For each target, up to six designs were evaluated for expression levels and in some cases, for thermal stability and activity. In nine targets, designs exhibited increased heterologous expression levels either in prokaryotic and/or eukaryotic expression systems under experimental conditions that were tailored for each target protein. Furthermore, we observed increased thermal stability in nine of ten tested targets. In two prime examples, the human Stem Cell Factor (hSCF) and human Cadherin-Like Domain (CLD12) from the RET receptor, the wild type proteins were not expressible as soluble proteins in E. coli, yet the PROSS designs exhibited high expression levels in E. coli and HEK293 cells, respectively, and improved thermal stability. We conclude that PROSS may improve stability and expressibility in diverse cases, and that improvement typically requires target-specific expression conditions. This study demonstrates the strengths of community-wide efforts to probe the generality of new methods and recommends areas for future research to advance practically useful algorithms for protein science.

Mutational fitness landscapes reveal genetic and structural improvement pathways for a vaccine-elicited HIV-1 broadly neutralizing antibody

Vaccine-based elicitation of broadly neutralizing antibodies holds great promise for preventing HIV-1 transmission. However, the key biophysical markers of improved antibody recognition remain uncertain in the diverse landscape of potential antibody mutation pathways, and a more complete understanding of anti-HIV-1 fusion peptide (FP) antibody development will accelerate rational vaccine designs. Here we survey the mutational landscape of the vaccine-elicited anti-FP antibody, vFP16.02, to determine the genetic, structural, and functional features associated with antibody improvement or fitness. Using site-saturation mutagenesis and yeast display functional screening, we found that 1% of possible single mutations improved HIV-1 envelope trimer (Env) affinity, but generally comprised rare somatic hypermutations that may not arise frequently in vivo. We observed that many single mutations in the vFP16.02 Fab could enhance affinity >1,000-fold against soluble FP, although affinity improvements against the HIV-1 trimer were more measured and rare. The most potent variants enhanced affinity to both soluble FP and Env, had mutations concentrated in antibody framework regions, and achieved up to 37% neutralization breadth compared to 28% neutralization of the template antibody. Altered heavy- and light-chain interface angles and conformational dynamics, as well as reduced Fab thermal stability, were associated with improved HIV-1 neutralization breadth and potency. We also observed parallel sets of mutations that enhanced viral neutralization through similar structural mechanisms. These data provide a quantitative understanding of the mutational landscape for vaccine-elicited FP-directed broadly neutralizing antibody and demonstrate that numerous antigen-distal framework mutations can improve antibody function by enhancing affinity simultaneously toward HIV-1 Env and FP.

A matrix of structure-based designs yields improved VRC01-class antibodies for HIV-1 therapy and prevention

Passive transfer of broadly neutralizing antibodies is showing promise in the treatment and prevention of HIV-1. One class of antibodies, the VRC01 class, appears especially promising. To improve VRC01-class antibodies, we combined structure-based design with a matrix-based approach to generate VRC01-class variants that filled an interfacial cavity, used diverse third-complementarity-determining regions, reduced potential steric clashes, or exploited extended contacts to a neighboring protomer within the envelope trimer. On a 208-strain panel, variant VRC01.23LS neutralized 90% of the panel at a geometric mean IC80 less than 1 μg/ml, and in transgenic mice with human neonatal-Fc receptor, the serum half-life of VRC01.23LS was indistinguishable from that of the parent VRC01LS, which has a half-life of 71 d in humans. A cryo-electron microscopy structure of VRC01.23 Fab in complex with BG505 DS-SOSIP.664 Env trimer determined at 3.4-Å resolution confirmed the structural basis for its ~10-fold improved potency relative to VRC01. Another variant, VRC07-523-F54-LS.v3, neutralized 95% of the 208-isolated panel at a geometric mean IC80 of less than 1 μg/ml, with a half-life comparable to that of the parental VRC07-523LS. Our matrix-based structural approach thus enables the engineering of VRC01 variants for HIV-1 therapy and prevention with improved potency, breadth, and pharmacokinetics.

Discovery-stage identification of drug-like antibodies using emerging experimental and computational methods

There is intense and widespread interest in developing monoclonal antibodies as therapeutic agents to treat diverse human disorders. During early-stage antibody discovery, hundreds to thousands of lead candidates are identified, and those that lack optimal physical and chemical properties must be deselected as early as possible to avoid problems later in drug development. It is particularly challenging to characterize such properties for large numbers of candidates with the low antibody quantities, concentrations, and purities that are available at the discovery stage, and to predict concentrated antibody properties (e.g., solubility, viscosity) required for efficient formulation, delivery, and efficacy. Here we review key recent advances in developing and implementing high-throughput methods for identifying antibodies with desirable in vitro and in vivo properties, including favorable antibody stability, specificity, solubility, pharmacokinetics, and immunogenicity profiles, that together encompass overall drug developability. In particular, we highlight impressive recent progress in developing computational methods for improving rational antibody design and prediction of drug-like behaviors that hold great promise for reducing the amount of required experimentation. We also discuss outstanding challenges that will need to be addressed in the future to fully realize the great potential of using such analysis for minimizing development times and improving the success rate of antibody candidates in the clinic.

The AbDesign computational pipeline for modular backbone assembly and design of binders and enzymes

The functional sites of many protein families are dominated by diverse backbone regions that lack secondary structure (loops) but fold stably into their functionally competent state. Nevertheless, the design of structured loop regions from scratch, especially in functional sites, has met with great difficulty. We therefore developed an approach, called AbDesign, to exploit the natural modularity of many protein families and computationally assemble a large number of new backbones by combining naturally occurring modular fragments. This strategy yielded large, atomically accurate, and highly efficient proteins, including antibodies and enzymes exhibiting dozens of mutations from any natural protein. The combinatorial backbone‐conformation space that can be accessed by AbDesign even for a modestly sized family of homologs may exceed the diversity in the entire PDB, providing the sub‐Ångstrom level of control over the positioning of active‐site groups that is necessary for obtaining highly active proteins. This manuscript describes how to implement the pipeline using code that is freely available at https://github.com/Fleishman‐Lab/AbDesign_for_enzymes.

Toward Drug-Like Multispecific Antibodies by Design

The success of antibody therapeutics is strongly influenced by their multifunctional nature that couples antigen recognition mediated by their variable regions with effector functions and half-life extension mediated by a subset of their constant regions. Nevertheless, the monospecific IgG format is not optimal for many therapeutic applications, and this has led to the design of a vast number of unique multispecific antibody formats that enable targeting of multiple antigens or multiple epitopes on the same antigen. Despite the diversity of these formats, a common challenge in generating multispecific antibodies is that they display suboptimal physical and chemical properties relative to conventional IgGs and are more difficult to develop into therapeutics. Here we review advances in the design and engineering of multispecific antibodies with drug-like properties, including favorable stability, solubility, viscosity, specificity and pharmacokinetic properties. We also highlight emerging experimental and computational methods for improving the next generation of multispecific antibodies, as well as their constituent antibody fragments, with natural IgG-like properties. Finally, we identify several outstanding challenges that need to be addressed to increase the success of multispecific antibodies in the clinic.

Mass Spectrometry Analysis of Intact Proteins from Crude Samples

Analysis of intact proteins by native mass spectrometry has emerged as a powerful tool for obtaining insight into subunit diversity, post-translational modifications, stoichiometry, structural arrangement, stability, and overall architecture. Typically, such an analysis is performed following protein purification procedures, which are time consuming, costly, and labor intensive. As this technology continues to move forward, advances in sample handling and instrumentation have enabled the investigation of intact proteins in situ and in crude samples, offering rapid analysis and improved conservation of the biological context. This emerging field, which involves various ion source platforms such as matrix-assisted laser desorption ionization (MALDI) and electrospray ionization (ESI) for both spatial imaging and solution-based analysis, is expected to impact many scientific fields, including biotechnology, pharmaceuticals, and clinical sciences. In this Perspective, we discuss the information that can be retrieved by such experiments as well as the current advantages and technical challenges associated with the different sampling strategies. Furthermore, we present future directions of these MS-based methods, including current limitations and efforts that should be made to make these approaches more accessible. Considering the vast progress we have witnessed in recent years, we anticipate that the advent of further innovations enabling minimal handling of MS samples will make this field more robust, user friendly, and widespread.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Practically useful protein-design methods combining phylogenetic and atomistic calculations

Our ability to design new or improved biomolecular activities depends on understanding the sequence-function relationships in proteins. The large size and fold complexity of most proteins, however, obscure these relationships, and protein-optimization methods continue to rely on laborious experimental iterations. Recently, a deeper understanding of the roles of stability-threshold effects and biomolecular epistasis in proteins has led to the development of hybrid methods that combine phylogenetic analysis with atomistic design calculations. These methods enable reliable and even single-step optimization of protein stability, expressibility, and activity in proteins that were considered outside the scope of computational design. Furthermore, ancestral-sequence reconstruction produces insights on missing links in the evolution of enzymes and binders that may be used in protein design. Through the combination of phylogenetic and atomistic calculations, the long-standing goal of general computational methods that can be universally applied to study and optimize proteins finally seems within reach.

Inhibition of fibroblast secreted QSOX1 perturbs extracellular matrix in the tumor microenvironment and decreases tumor growth and metastasis in murine cancer models

Extracellular matrix (ECM) plays an important role in tumor development and dissemination, but few points of therapeutic intervention targeting ECM of the tumor microenvironment have been exploited to date. Recent observations suggest that the enzymatic introduction of disulfide bond cross-links into the ECM may be modulated to affect cancer progression. Specifically, the disulfide bond-forming activity of the enzyme Quiescin sulfhydryl oxidase 1 (QSOX1) is required by fibroblasts to assemble ECM components for adhesion and migration of cancer cells. Based on this finding and the increased QSOX1 expression in the stroma of aggressive breast carcinomas, we developed monoclonal antibody inhibitors with the aim of preventing QSOX1 from participating in pro-metastatic ECM remodeling. Here we show that QSOX1 inhibitory antibodies decreased tumor growth and metastasis in murine cancer models and had added benefits when provided together with chemotherapy. Mechanistically, the inhibitors dampened stromal participation in tumor development, as the tumors of treated animals showed fewer myofibroblasts and poorer ECM organization. Thus, our findings demonstrate that specifically targeting excess stromal QSOX1 secreted in response to tumor-cell signaling provides a means to modulate the tumor microenvironment and may complement other therapeutic approaches in cancer.

A lipophilicity-based energy function for membrane-protein modelling and design

Membrane-protein design is an exciting and increasingly successful research area which has led to landmarks including the design of stable and accurate membrane-integral proteins based on coiled-coil motifs. Design of topologically more complex proteins, such as most receptors, channels, and transporters, however, demands an energy function that balances contributions from intra-protein contacts and protein-membrane interactions. Recent advances in water-soluble all-atom energy functions have increased the accuracy in structure-prediction benchmarks. The plasma membrane, however, imposes different physical constraints on protein solvation. To understand these constraints, we recently developed a high-throughput experimental screen, called dsTβL, and inferred apparent insertion energies for each amino acid at dozens of positions across the bacterial plasma membrane. Here, we express these profiles as lipophilicity energy terms in Rosetta and demonstrate that the new energy function outperforms previous ones in modelling and design benchmarks. Rosetta ab initio simulations starting from an extended chain recapitulate two-thirds of the experimentally determined structures of membrane-spanning homo-oligomers with <2.5Å root-mean-square deviation within the top-predicted five models (available online: http://tmhop.weizmann.ac.il). Furthermore, in two sequence-design benchmarks, the energy function improves discrimination of stabilizing point mutations and recapitulates natural membrane-protein sequences of known structure, thereby recommending this new energy function for membrane-protein modelling and design.