High-accuracy modeling of antibody structures by a search for minimum-energy recombination of backbone fragments

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.

Computational Methods in Immunology and Vaccinology: Design and Development of Antibodies and Immunogens

The design of new biomolecules able to harness immune mechanisms for the treatment of diseases is a prime challenge for computational and simulative approaches. For instance, in recent years, antibodies have emerged as an important class of therapeutics against a spectrum of pathologies. In cancer, immune-inspired approaches are witnessing a surge thanks to a better understanding of tumor-associated antigens and the mechanisms of their engagement or evasion from the human immune system. Here, we provide a summary of the main state-of-the-art computational approaches that are used to design antibodies and antigens, and in parallel, we review key methodologies for epitope identification for both B- and T-cell mediated responses. A special focus is devoted to the description of structure- and physics-based models, privileged over purely sequence-based approaches. We discuss the implications of novel methods in engineering biomolecules with tailored immunological properties for possible therapeutic uses. Finally, we highlight the extraordinary challenges and opportunities presented by the possible integration of structure- and physics-based methods with emerging Artificial Intelligence technologies for the prediction and design of novel antigens, epitopes, and antibodies.

Binding affinity prediction for antibody-protein antigen complexes: A machine learning analysis based on interface and surface areas

Specific antibodies can bind to protein antigens with high affinity and specificity, and this property makes them one of the best protein-based therapeutics. Accurate prediction of antibody‒protein antigen binding affinity is crucial for designing effective antibodies. The current predictive methods for protein‒protein binding affinity usually fail to predict the binding affinity of an antibody‒protein antigen complex with a comparable level of accuracy. Here, new models specific for antibody‒antigen binding affinity prediction are developed according to the different types of interface and surface areas present in antibody‒antigen complex. The contacts-based descriptors are also employed to construct or train different models specific for antibody‒protein antigen binding affinity prediction. The results of this study show that (i) the area-based descriptors are slightly better than the contacts-based descriptors in terms of the predictive power; (ii) the new models specific for antibody‒protein antigen binding affinity prediction are superior to the previously-used general models for predicting the protein‒protein binding affinities; (iii) the performances of the best area-based and contacts-based models developed in this work are better than the performances of a recently-developed graph-based model (i.e., CSM-AB) specific for antibody‒protein antigen binding affinity prediction. The new models developed in this work would not only help understand the mechanisms underlying antibody‒protein antigen interactions, but would also be of some applicable utility in the design and virtual screening of antibody-based therapeutics.

Multi-state modeling of antibody-antigen complexes with SAXS profiles and deep-learning models

Antibodies are an established class of human therapeutics. Epitope characterization is an important part of therapeutic antibody discovery. However, structural characterization of antibody-antigen complexes remains challenging. On the one hand, X-ray crystallography or cryo-electron microscopy provide atomic resolution characterization of the epitope, but the data collection process is typically long and the success rate is low. On the other hand, computational methods for modeling antibody-antigen structures from the individual components frequently suffer from a high false positive rate, rarely resulting in a unique solution. Recent deep learning models for structure prediction are also successful in predicting protein-protein complexes. However, they do not perform well for antibody-antigen complexes. Small Angle X-ray Scattering (SAXS) is a reliable technique for rapid structural characterization of protein samples in solution albeit at low resolution. Here, we present an integrative approach for modeling antigen-antibody complexes using the antibody sequence, antigen structure, and experimentally determined SAXS profiles of the antibody, antigen, and the complex. The method models antibody structures using a novel deep-learning approach, NanoNet. The structures of the antibodies and antigens are represented using multiple 3D conformations to account for compositional and conformational heterogeneity of the protein samples that are used to collect the SAXS data. The complexes are predicted by integrating the SAXS profiles with scoring functions for protein-protein interfaces that are based on statistical potentials and antibody-specific deep-learning models. We validated the method via application to four Fab:EGFR and one Fab:PCSK9 antibody:antigen complexes with experimentally available SAXS datasets. The integrative approach returns accurate predictions (interface RMSD<4Å) in the top five predictions for four out of five complexes (respective interface RMSD values of 1.95, 2.18, 2.66 and 3.87Å), providing support for the utility of such a computational pipeline for epitope characterization during therapeutic antibody discovery.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

VHH Structural Modelling Approaches: A Critical Review

VHH, i.e., VH domains of camelid single-chain antibodies, are very promising therapeutic agents due to their significant physicochemical advantages compared to classical mammalian antibodies. The number of experimentally solved VHH structures has significantly improved recently, which is of great help, because it offers the ability to directly work on 3D structures to humanise or improve them. Unfortunately, most VHHs do not have 3D structures. Thus, it is essential to find alternative ways to get structural information. The methods of structure prediction from the primary amino acid sequence appear essential to bypass this limitation. This review presents the most extensive overview of structure prediction methods applied for the 3D modelling of a given VHH sequence (a total of 21). Besides the historical overview, it aims at showing how model software programs have been shaping the structural predictions of VHHs. A brief explanation of each methodology is supplied, and pertinent examples of their usage are provided. Finally, we present a structure prediction case study of a recently solved VHH structure. According to some recent studies and the present analysis, AlphaFold 2 and NanoNet appear to be the best tools to predict a structural model of VHH from its sequence.

The neutralization potency of anti-SARS-CoV-2 therapeutic human monoclonal antibodies is retained against viral variants

A wide range of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) neutralizing monoclonal antibodies (mAbs) have been reported, most of which target the spike glycoprotein. Therapeutic implementation of these antibodies has been challenged by emerging SARS-CoV-2 variants harboring mutated spike versions. Consequently, re-assessment of previously identified mAbs is of high priority. Four previously selected mAbs targeting non-overlapping epitopes are now evaluated for binding potency to mutated RBD versions, reported to mediate escape from antibody neutralization. In vitro neutralization potencies of these mAbs, and two NTD-specific mAbs, are evaluated against two frequent SARS-CoV-2 variants of concern, the B.1.1.7 Alpha and the B.1.351 Beta. Furthermore, we demonstrate therapeutic potential of three selected mAbs by treatment of K18-human angiotensin-converting enzyme 2 (hACE2) transgenic mice 2 days post-infection with each virus variant. Thus, despite the accumulation of spike mutations, the highly potent MD65 and BL6 mAbs retain their ability to bind the prevalent viral mutants, effectively protecting against B.1.1.7 and B.1.351 variants.

RosettaCM for antibodies with very long HCDR3s and low template availability

Antibody-antigen co-crystal structures are a valuable resource for the fundamental understanding of antibody-mediated immunity. Determination of structures with antibodies in complex with their antigens, however, is a laborious task without guarantee of success. Therefore, homology modeling of antibodies and docking to their respective antigens has become a very important technique to drive antibody and vaccine design. The quality of the antibody modeling process is critical for the success of these endeavors. Here, we compare different computational protocols for predicting antibody structure from sequence in the biomolecular modeling software Rosetta-all of which use multiple existing antibody structures to guide modeling. Specifically, we compare protocols developed solely to predict antibody structure (RosettaAntibody, AbPredict) with a universal homology modeling protocol (RosettaCM). Following recent advances in homology modeling with multiple templates simultaneously, we propose that the use of multiple templates over the same antibody regions may improve modeling performance. To evaluate whether multi-template comparative modeling with RosettaCM can improve the modeling accuracy of antibodies over existing methods, this study compares the performance of the three modeling algorithms when modeling human antibodies taken from antibody-antigen co-crystal structures. In these benchmarking experiments, RosettaCM outperformed other methods when modeling antibodies with long HCDR3s and few available templates.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

Toward Biotherapeutics Formulation Composition Engineering using Site-Identification by Ligand Competitive Saturation (SILCS)

Formulation of protein-based therapeutics employ advanced formulation and analytical technologies for screening various parameters such as buffer, pH, and excipients. At a molecular level, physico-chemical properties of a protein formulation depend on self-interaction between protein molecules, protein-solvent and protein-excipient interactions. This work describes a novel in silico approach, SILCS-Biologics, for structure-based modeling of protein formulations. SILCS Biologics is based on the Site-Identification by Ligand Competitive Saturation (SILCS) technology and enables modeling of interactions among different components of a formulation at an atomistic level while accounting for protein flexibility. It predicts potential hotspot regions on the protein surface for protein-protein and protein-excipient interactions. Here we apply SILCS-Biologics on a Fab domain of a monoclonal antibody (mAbN) to model Fab-Fab interactions and interactions with three amino acid excipients, namely, arginine HCl, proline and lysine HCl. Experiments on 100 mg/ml formulations of mAbN showed that arginine increased, lysine reduced, and proline did not impact viscosity. We use SILCS-Biologics modeling to explore a structure-based hypothesis for the viscosity modulating effect of these excipients. Current efforts are aimed at further validation of this novel computational framework and expanding the scope to model full mAb and other protein therapeutics.

Macromolecular modeling and design in Rosetta: recent methods and frameworks

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

AbPredict 2: a server for accurate and unstrained structure prediction of antibody variable domains

SUMMARY

Methods for antibody structure prediction rely on sequence homology to experimentally determined structures. Resulting models may be accurate but are often stereochemically strained, limiting their usefulness in modeling and design workflows. We present the AbPredict 2 web-server, which instead of using sequence homology, conducts a Monte Carlo-based search for low-energy combinations of backbone conformations to yield accurate and unstrained antibody structures.

AVAILABILITY AND IMPLEMENTATION

We introduce several important improvements over the previous AbPredict implementation: (i) backbones and sidechains are now modeled using ideal bond lengths and angles, substantially reducing stereochemical strain, (ii) sampling of the rigid-body orientation at the light-heavy chain interface is improved, increasing model accuracy and (iii) runtime is reduced 20-fold without compromising accuracy, enabling the implementation of AbPredict 2 as a fully automated web-server (http://abpredict.weizmann.ac.il). Accurate and unstrained antibody model structures may in some cases obviate the need for experimental structures in antibody optimization workflows.

Polyclonal and convergent antibody response to Ebola virus vaccine rVSV-ZEBOV

Recombinant vesicular stomatitis virus-Zaire Ebola virus (rVSV-ZEBOV) is the most advanced Ebola virus vaccine candidate and is currently being used to combat the outbreak of Ebola virus disease (EVD) in the Democratic Republic of the Congo (DRC). Here we examine the humoral immune response in a subset of human volunteers enrolled in a phase 1 rVSV-ZEBOV vaccination trial by performing comprehensive single B cell and electron microscopy structure analyses. Four studied vaccinees show polyclonal, yet reproducible and convergent B cell responses with shared sequence characteristics. EBOV-targeting antibodies cross-react with other Ebolavirus species, and detailed epitope mapping revealed overlapping target epitopes with antibodies isolated from EVD survivors. Moreover, in all vaccinees, we detected highly potent EBOV-neutralizing antibodies with activities comparable or superior to the monoclonal antibodies currently used in clinical trials. These include antibodies combining the IGHV3-15/IGLV1-40 immunoglobulin gene segments that were identified in all investigated individuals. Our findings will help to evaluate and direct current and future vaccination strategies and offer opportunities for novel EVD therapies.

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.

A combined computational-experimental approach to define the structural origin of antibody recognition of sialyl-Tn, a tumor-associated carbohydrate antigen

Anti-carbohydrate monoclonal antibodies (mAbs) hold great promise as cancer therapeutics and diagnostics. However, their specificity can be mixed, and detailed characterization is problematic, because antibody-glycan complexes are challenging to crystallize. Here, we developed a generalizable approach employing high-throughput techniques for characterizing the structure and specificity of such mAbs, and applied it to the mAb TKH2 developed against the tumor-associated carbohydrate antigen sialyl-Tn (STn). The mAb specificity was defined by apparent K_D values determined by quantitative glycan microarray screening. Key residues in the antibody combining site were identified by site-directed mutagenesis, and the glycan-antigen contact surface was defined using saturation transfer difference NMR (STD-NMR). These features were then employed as metrics for selecting the optimal 3D-model of the antibody-glycan complex, out of thousands plausible options generated by automated docking and molecular dynamics simulation. STn-specificity was further validated by computationally screening of the selected antibody 3D-model against the human sialyl-Tn-glycome. This computational-experimental approach would allow rational design of potent antibodies targeting carbohydrates.

Understanding and overcoming trade-offs between antibody affinity, specificity, stability and solubility

The widespread use of monoclonal antibodies for therapeutic applications has led to intense interest in optimizing several of their natural properties (affinity, specificity, stability, solubility and effector functions) as well as introducing non-natural activities (bispecificity and cytotoxicity mediated by conjugated drugs). A common challenge during antibody optimization is that improvements in one property (e.g., affinity) can lead to deficits in other properties (e.g., stability). Here we review recent advances in understanding trade-offs between different antibody properties, including affinity, specificity, stability and solubility. We also review new approaches for co-optimizing multiple antibody properties and discuss how these methods can be used to rapidly and systematically generate antibodies for a wide range of applications.

The path to VICTORy - a beginner's guide to success using commercial research antibodies

Commercial research antibodies are crucial tools in modern cell biology and biochemistry. In the USA some $2 billion a year are spent on them, but many are apparently not fit-for-purpose, and this may contribute to the 'reproducibility crisis' in biological sciences. Inadequate antibody validation and characterization, lack of user awareness, and occasional incompetence amongst suppliers have had immense scientific and personal costs. In this Opinion, I suggest some paths to make the use of these vital tools more successful. I have attempted to summarize and extend expert views from the literature to suggest that sustained routine efforts should made in: (1) the validation of antibodies, (2) their identification, (3) communication and controls, (4) the training of potential users, (5) the transparency of original equipment manufacturer (OEM) marketing agreements, and (5) in a more widespread use of recombinant antibodies (together denoted the 'VICTOR' approach).

Improved antibody-based ricin neutralization by affinity maturation is correlated with slower off-rate values

While potent monoclonal antibodies against ricin were introduced over the years, the question whether increasing antibody affinity enables better toxin neutralization was not fully addressed yet. The aim of this study was to characterize the contribution of antibody affinity to the ricin neutralization potential of the antibody. cHD23 monoclonal antibody that targets the toxin B-subunit and interferes with its binding to membranal receptors, was isolated. In order to create antibody clones with improved affinity toward ricin, a scFv-phage display library containing mutated versions of the variable regions of cHD23 was constructed and clones with improved binding of ricin were isolated. Structural modeling of these mutants suggests that the inserted mutations may increase the antibody conformational flexibility thus improving its ability to bind ricin. While it was found that the selected clones exhibited improved neutralization of ricin, the correlation between the KD values and potency was only minor (r = 0.55). However, a positive correlation (r = 0.84) exist between the off-rate values (koff) of the affinity matured clones and their ability to neutralize ricin. As cell membranes display inordinately large amounts of potential surface binding sites for ricin, it is suggested that antibodies with improved off-rate values block the ability of the toxin to bind to target receptors, in a highly efficient manner. Currently, antibody-based therapy is the most effective treatment for ricin intoxication and it is anticipated that the findings of this study will provide useful information and a possible strategy to design an improved antibody-based therapy for the toxin.

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.