SnugDock: paratope structural optimization during antibody-antigen docking compensates for errors in antibody homology models

WUREN: Whole-modal union representation for epitope prediction

B-cell epitope identification plays a vital role in the development of vaccines, therapies, and diagnostic tools. Currently, molecular docking tools in B-cell epitope prediction are heavily influenced by empirical parameters and require significant computational resources, rendering a great challenge to meet large-scale prediction demands. When predicting epitopes from antigen-antibody complex, current artificial intelligence algorithms cannot accurately implement the prediction due to insufficient protein feature representations, indicating novel algorithm is desperately needed for efficient protein information extraction. In this paper, we introduce a multimodal model called WUREN (Whole-modal Union Representation for Epitope predictioN), which effectively combines sequence, graph, and structural features. It achieved AUC-PR scores of 0.213 and 0.193 on the solved structures and AlphaFold-generated structures, respectively, for the independent test proteins selected from DiscoTope3 benchmark. Our findings indicate that WUREN is an efficient feature extraction model for protein complexes, with the generalizable application potential in the development of protein-based drugs. Moreover, the streamlined framework of WUREN could be readily extended to model similar biomolecules, such as nucleic acids, carbohydrates, and lipids.

A comparison of antibody-antigen complex sequence-to-structure prediction methods and their systematic biases

The ability to accurately predict antibody-antigen complex structures from their sequences could greatly advance our understanding of the immune system and would aid in the development of novel antibody therapeutics. There have been considerable recent advancements in predicting protein-protein interactions (PPIs) fueled by progress in machine learning (ML). To understand the current state of the field, we compare six representative methods for predicting antibody-antigen complexes from sequence, including two deep learning approaches trained to predict PPIs in general (AlphaFold-Multimer and RoseTTAFold), two composite methods that initially predict antibody and antigen structures separately and dock them (using antibody-mode ClusPro), local refinement in Rosetta (SnugDock) of globally docked poses from ClusPro, and a pipeline combining homology modeling with rigid-body docking informed by ML-based epitope and paratope prediction (AbAdapt). We find that AlphaFold-Multimer outperformed other methods, although the absolute performance leaves considerable room for improvement. AlphaFold-Multimer models of lower quality display significant structural biases at the level of tertiary motifs (TERMs) toward having fewer structural matches in non-antibody-containing structures from the Protein Data Bank (PDB). Specifically, better models exhibit more common PDB-like TERMs at the antibody-antigen interface than worse ones. Importantly, the clear relationship between performance and the commonness of interfacial TERMs suggests that the scarcity of interfacial geometry data in the structural database may currently limit the application of ML to the prediction of antibody-antigen interactions.

De novo generation of SARS-CoV-2 antibody CDRH3 with a pre-trained generative large language model

Artificial Intelligence (AI) techniques have made great advances in assisting antibody design. However, antibody design still heavily relies on isolating antigen-specific antibodies from serum, which is a resource-intensive and time-consuming process. To address this issue, we propose a Pre-trained Antibody generative large Language Model (PALM-H3) for the de novo generation of artificial antibodies heavy chain complementarity-determining region 3 (CDRH3) with desired antigen-binding specificity, reducing the reliance on natural antibodies. We also build a high-precision model antigen-antibody binder (A2binder) that pairs antigen epitope sequences with antibody sequences to predict binding specificity and affinity. PALM-H3-generated antibodies exhibit binding ability to SARS-CoV-2 antigens, including the emerging XBB variant, as confirmed through in-silico analysis and in-vitro assays. The in-vitro assays validate that PALM-H3-generated antibodies achieve high binding affinity and potent neutralization capability against spike proteins of SARS-CoV-2 wild-type, Alpha, Delta, and the emerging XBB variant. Meanwhile, A2binder demonstrates exceptional predictive performance on binding specificity for various epitopes and variants. Furthermore, by incorporating the attention mechanism inherent in the Roformer architecture into the PALM-H3 model, we improve its interpretability, providing crucial insights into the fundamental principles of antibody design.

Vaccine design and development: Exploring the interface with computational biology and AI

Computational biology involves applying computer science and informatics techniques in biology to understand complex biological data. It allows us to collect, connect, and analyze biological data at a large scale and build predictive models. In the twenty first century, computational resources along with Artificial Intelligence (AI) have been widely used in various fields of biological sciences such as biochemistry, structural biology, immunology, microbiology, and genomics to handle massive data for decision-making, including in applications such as drug design and vaccine development, one of the major areas of focus for human and animal welfare. The knowledge of available computational resources and AI-enabled tools in vaccine design and development can improve our ability to conduct cutting-edge research. Therefore, this review article aims to summarize important computational resources and AI-based tools. Further, the article discusses the various applications and limitations of AI tools in vaccine development.

A Comparison of Antibody-Antigen Complex Sequence-to-Structure Prediction Methods and their Systematic Biases

The ability to accurately predict antibody-antigen complex structures from their sequences could greatly advance our understanding of the immune system and would aid in the development of novel antibody therapeutics. There have been considerable recent advancements in predicting protein-protein interactions (PPIs) fueled by progress in machine learning (ML). To understand the current state of the field, we compare six representative methods for predicting antibody-antigen complexes from sequence, including two deep learning approaches trained to predict PPIs in general (AlphaFold-Multimer, RoseTTAFold), two composite methods that initially predict antibody and antigen structures separately and dock them (using antibody-mode ClusPro), local refinement in Rosetta (SnugDock) of globally docked poses from ClusPro, and a pipeline combining homology modeling with rigid-body docking informed by ML-based epitope and paratope prediction (AbAdapt). We find that AlphaFold-Multimer outperformed other methods, although the absolute performance leaves considerable room for improvement. AlphaFold-Multimer models of lower-quality display significant structural biases at the level of tertiary motifs (TERMs) towards having fewer structural matches in non-antibody containing structures from the Protein Data Bank (PDB). Specifically, better models exhibit more common PDB-like TERMs at the antibody-antigen interface than worse ones. Importantly, the clear relationship between performance and the commonness of interfacial TERMs suggests that scarcity of interfacial geometry data in the structural database may currently limit application of machine learning to the prediction of antibody-antigen interactions.

A new era of antibody discovery: an in-depth review of AI-driven approaches

Given their high affinity and specificity for a range of macromolecules, antibodies are widely used in the treatment of autoimmune diseases, cancers, inflammatory diseases, and Alzheimer's disease (AD). Traditional experimental methods are time-consuming, expensive, and labor-intensive. Recent advances in artificial intelligence (AI) technologies provide complementary methods that can reduce the time and costs required for antibody design by minimizing failures and increasing the success rate of experimental tests. In this review, we scrutinize the plethora of AI-driven methodologies that have been deployed over the past 4 years for modeling antibody structures, predicting antibody-antigen interactions, optimizing antibody affinity, and generating novel antibody candidates. We also briefly address the challenges faced in integrating AI-based models with traditional antibody discovery pipelines and highlight the potential future directions in this burgeoning field.

A Dual Strategy-In Vitro and In Silico-To Evaluate Human Antitetanus mAbs Addressing Their Potential Protective Action on TeNT Endocytosis in Primary Rat Neuronal Cells

Tetanus disease, caused by C. tetani, starts with wounds or mucous layer contact. Prevented by vaccination, the lack of booster shots throughout life requires prophylactic treatment in case of accidents. The incidence of tetanus is high in underdeveloped countries, requiring the administration of antitetanus antibodies, usually derived from immunized horses or humans. Heterologous sera represent risks such as serum sickness. Human sera can carry unknown viruses. In the search for human monoclonal antibodies (mAbs) against TeNT (Tetanus Neurotoxin), we previously identified a panel of mAbs derived from B-cell sorting, selecting two nonrelated ones that binded to the C-terminal domain of TeNT (HCR/T), inhibiting its interaction with the cellular receptor ganglioside GT1b. Here, we present the results of cellular assays and molecular docking tools. TeNT internalization in neurons is prevented by more than 50% in neonatal rat spinal cord cells, determined by quantitative analysis of immunofluorescence punctate staining of Alexa Fluor 647 conjugated to TeNT. We also confirmed the mediator role of the Synaptic Vesicle Glycoprotein II (SV2) in TeNT endocytosis. The molecular docking assays to predict potential TeNT epitopes showed the binding of both antibodies to the HCR/T domain. A higher incidence was found between N1153 and W1297 when evaluating candidate residues for conformational epitope.

Prediction of Paratope-Epitope Pairs Using Convolutional Neural Networks

Antibodies play a central role in the adaptive immune response of vertebrates through the specific recognition of exogenous or endogenous antigens. The rational design of antibodies has a wide range of biotechnological and medical applications, such as in disease diagnosis and treatment. However, there are currently no reliable methods for predicting the antibodies that recognize a specific antigen region (or epitope) and, conversely, epitopes that recognize the binding region of a given antibody (or paratope). To fill this gap, we developed ImaPEp, a machine learning-based tool for predicting the binding probability of paratope-epitope pairs, where the epitope and paratope patches were simplified into interacting two-dimensional patches, which were colored according to the values of selected features, and pixelated. The specific recognition of an epitope image by a paratope image was achieved by using a convolutional neural network-based model, which was trained on a set of two-dimensional paratope-epitope images derived from experimental structures of antibody-antigen complexes. Our method achieves good performances in terms of cross-validation with a balanced accuracy of 0.8. Finally, we showcase examples of application of ImaPep, including extensive screening of large libraries to identify paratope candidates that bind to a selected epitope, and rescoring and refining antibody-antigen docking poses.

Unveiling the affinity-stability relationship in anti-measles virus antibodies: a computational approach for hotspots prediction

Recent years have seen an uptick in the use of computational applications in antibody engineering. These tools have enhanced our ability to predict interactions with antigens and immunogenicity, facilitate humanization, and serve other critical functions. However, several studies highlight the concern of potential trade-offs between antibody affinity and stability in antibody engineering. In this study, we analyzed anti-measles virus antibodies as a case study, to examine the relationship between binding affinity and stability, upon identifying the binding hotspots. We leverage in silico tools like Rosetta and FoldX, along with molecular dynamics (MD) simulations, offering a cost-effective alternative to traditional in vitro mutagenesis. We introduced a pattern in identifying key residues in pairs, shedding light on hotspots identification. Experimental physicochemical analysis validated the predicted key residues by confirming significant decrease in binding affinity for the high-affinity antibodies to measles virus hemagglutinin. Through the nature of the identified pairs, which represented the relative hydropathy of amino acid side chain, a connection was proposed between affinity and stability. The findings of the study enhance our understanding of the interactions between antibody and measles virus hemagglutinin. Moreover, the implications of the observed correlation between binding affinity and stability extend beyond the field of anti-measles virus antibodies, thereby opening doors for advancements in antibody research.

ABAG-docking benchmark: a non-redundant structure benchmark dataset for antibody-antigen computational docking

Accurate prediction of antibody-antigen complex structures is pivotal in drug discovery, vaccine design and disease treatment and can facilitate the development of more effective therapies and diagnostics. In this work, we first review the antibody-antigen docking (ABAG-docking) datasets. Then, we present the creation and characterization of a comprehensive benchmark dataset of antibody-antigen complexes. We categorize the dataset based on docking difficulty, interface properties and structural characteristics, to provide a diverse set of cases for rigorous evaluation. Compared with Docking Benchmark 5.5, we have added 112 cases, including 14 single-domain antibody (sdAb) cases and 98 monoclonal antibody (mAb) cases, and also increased the proportion of Difficult cases. Our dataset contains diverse cases, including human/humanized antibodies, sdAbs, rodent antibodies and other types, opening the door to better algorithm development. Furthermore, we provide details on the process of building the benchmark dataset and introduce a pipeline for periodic updates to keep it up to date. We also utilize multiple complex prediction methods including ZDOCK, ClusPro, HDOCK and AlphaFold-Multimer for testing and analyzing this dataset. This benchmark serves as a valuable resource for evaluating and advancing docking computational methods in the analysis of antibody-antigen interaction, enabling researchers to develop more accurate and effective tools for predicting and designing antibody-antigen complexes. The non-redundant ABAG-docking structure benchmark dataset is available at https://github.com/Zhaonan99/Antibody-antigen-complex-structure-benchmark-dataset.

Evaluation of AlphaFold antibody-antigen modeling with implications for improving predictive accuracy

High resolution antibody-antigen structures provide critical insights into immune recognition and can inform therapeutic design. The challenges of experimental structural determination and the diversity of the immune repertoire underscore the necessity of accurate computational tools for modeling antibody-antigen complexes. Initial benchmarking showed that despite overall success in modeling protein-protein complexes, AlphaFold and AlphaFold-Multimer have limited success in modeling antibody-antigen interactions. In this study, we performed a thorough analysis of AlphaFold's antibody-antigen modeling performance on 427 nonredundant antibody-antigen complex structures, identifying useful confidence metrics for predicting model quality, and features of complexes associated with improved modeling success. Notably, we found that the latest version of AlphaFold improves near-native modeling success to over 30%, versus approximately 20% for a previous version, while increased AlphaFold sampling gives approximately 50% success. With this improved success, AlphaFold can generate accurate antibody-antigen models in many cases, while additional training or other optimization may further improve performance.

Next generation of multispecific antibody engineering

Multispecific antibodies recognize two or more epitopes located on the same or distinct targets. This added capability through protein design allows these man-made molecules to address unmet medical needs that are no longer possible with single targeting such as with monoclonal antibodies or cytokines alone. However, the approach to the development of these multispecific molecules has been met with numerous road bumps, which suggests that a new workflow for multispecific molecules is required. The investigation of the molecular basis that mediates the successful assembly of the building blocks into non-native quaternary structures will lead to the writing of a playbook for multispecifics. This is a must do if we are to design workflows that we can control and in turn predict success. Here, we reflect on the current state-of-the-art of therapeutic biologics and look at the building blocks, in terms of proteins, and tools that can be used to build the foundations of such a next-generation workflow.

Simultaneous selection of nanobodies for accessible epitopes on immune cells in the tumor microenvironment

In the rapidly advancing field of synthetic biology, there exists a critical need for technology to discover targeting moieties for therapeutic biologics. Here we present INSPIRE-seq, an approach that utilizes a nanobody library and next-generation sequencing to identify nanobodies selected for complex environments. INSPIRE-seq enables the parallel enrichment of immune cell-binding nanobodies that penetrate the tumor microenvironment. Clone enrichment and specificity vary across immune cell subtypes in the tumor, lymph node, and spleen. INSPIRE-seq identifies a dendritic cell binding clone that binds PHB2. Single-cell RNA sequencing reveals a connection with cDC1s, and immunofluorescence confirms nanobody-PHB2 colocalization along cell membranes. Structural modeling and docking studies assist binding predictions and will guide nanobody selection. In this work, we demonstrate that INSPIRE-seq offers an unbiased approach to examine complex microenvironments and assist in the development of nanobodies, which could serve as active drugs, modified to become drugs, or used as targeting moieties.

Enhanced antibody-antigen structure prediction from molecular docking using AlphaFold2

Predicting the structure of antibody-antigen complexes has tremendous value in biomedical research but unfortunately suffers from a poor performance in real-life applications. AlphaFold2 (AF2) has provided renewed hope for improvements in the field of protein-protein docking but has shown limited success against antibody-antigen complexes due to the lack of co-evolutionary constraints. In this study, we used physics-based protein docking methods for building decoy sets consisting of low-energy docking solutions that were either geometrically close to the native structure (positives) or not (negatives). The docking models were then fed into AF2 to assess their confidence with a novel composite score based on normalized pLDDT and pTMscore metrics after AF2 structural refinement. We show benefits of the AF2 composite score for rescoring docking poses both in terms of (1) classification of positives/negatives and of (2) success rates with particular emphasis on early enrichment. Docking models of at least medium quality present in the decoy set, but not necessarily highly ranked by docking methods, benefitted most from AF2 rescoring by experiencing large advances towards the top of the reranked list of models. These improvements, obtained without any calibration or novel methodologies, led to a notable level of performance in antibody-antigen unbound docking that was never achieved previously.

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.

How can we discover developable antibody-based biotherapeutics?

Antibody-based biotherapeutics have emerged as a successful class of pharmaceuticals despite significant challenges and risks to their discovery and development. This review discusses the most frequently encountered hurdles in the research and development (R&D) of antibody-based biotherapeutics and proposes a conceptual framework called biopharmaceutical informatics. Our vision advocates for the syncretic use of computation and experimentation at every stage of biologic drug discovery, considering developability (manufacturability, safety, efficacy, and pharmacology) of potential drug candidates from the earliest stages of the drug discovery phase. The computational advances in recent years allow for more precise formulation of disease concepts, rapid identification, and validation of targets suitable for therapeutic intervention and discovery of potential biotherapeutics that can agonize or antagonize them. Furthermore, computational methods for de novo and epitope-specific antibody design are increasingly being developed, opening novel computationally driven opportunities for biologic drug discovery. Here, we review the opportunities and limitations of emerging computational approaches for optimizing antigens to generate robust immune responses, in silico generation of antibody sequences, discovery of potential antibody binders through virtual screening, assessment of hits, identification of lead drug candidates and their affinity maturation, and optimization for developability. The adoption of biopharmaceutical informatics across all aspects of drug discovery and development cycles should help bring affordable and effective biotherapeutics to patients more quickly.

Evaluation of AlphaFold Antibody-Antigen Modeling with Implications for Improving Predictive Accuracy

High resolution antibody-antigen structures provide critical insights into immune recognition and can inform therapeutic design. The challenges of experimental structural determination and the diversity of the immune repertoire underscore the necessity of accurate computational tools for modeling antibody-antigen complexes. Initial benchmarking showed that despite overall success in modeling protein-protein complexes, AlphaFold and AlphaFold-Multimer have limited success in modeling antibody-antigen interactions. In this study, we performed a thorough analysis of AlphaFold's antibody-antigen modeling performance on 429 nonredundant antibody-antigen complex structures, identifying useful confidence metrics for predicting model quality, and features of complexes associated with improved modeling success. We show the importance of bound-like component modeling in complex assembly accuracy, and that the current version of AlphaFold improves near-native modeling success to over 30%, versus approximately 20% for a previous version. With this improved success, AlphaFold can generate accurate antibody-antigen models in many cases, while additional training may further improve its performance.

Human-murine chimeric autoantibodies with high affinity and specificity for systemic sclerosis

Scleroderma 70 (Scl-70) is commonly used in the clinic for aiding systemic sclerosis (SSc) diagnosis due to its recognition as autoantibodies in the serum of SSc patients. However, obtaining sera positive for anti-Scl-70 antibody can be challenging; therefore, there is an urgent need to develop a specific, sensitive, and easily available reference for SSc diagnosis. In this study, murine-sourced scFv library were screened by phage display technology against human Scl-70, and the scFvs with high affinity were constructed into humanized antibodies for clinical application. Finally, ten high-affinity scFv fragments were obtained. Three fragments (2A, 2AB, and 2HD) were select for humanization. The physicochemical properties of the amino acid sequence, three-dimensional structural basis, and electrostatic potential distribution of the protein surface of different scFv fragments revealed differences in the electrostatic potential of their CDR regions determined their affinity for Scl-70 and expression. Notably, the specificity test showed the half-maximal effective concentration values of the three humanized antibodies were lower than that of positive patient serum. Moreover, these humanized antibodies showed high specificity for Scl-70 in diagnostic immunoassays for ANA. Among these three antibodies, 2A exhibited most positive electrostatic potential on the surface of the CDRs and highest affinity and specificity for Scl-70 but with least expression level; thus, it may provide new foundations for developing enhanced diagnostic strategies for SSc.

Fast, accurate antibody structure prediction from deep learning on massive set of natural antibodies

Antibodies have the capacity to bind a diverse set of antigens, and they have become critical therapeutics and diagnostic molecules. The binding of antibodies is facilitated by a set of six hypervariable loops that are diversified through genetic recombination and mutation. Even with recent advances, accurate structural prediction of these loops remains a challenge. Here, we present IgFold, a fast deep learning method for antibody structure prediction. IgFold consists of a pre-trained language model trained on 558 million natural antibody sequences followed by graph networks that directly predict backbone atom coordinates. IgFold predicts structures of similar or better quality than alternative methods (including AlphaFold) in significantly less time (under 25 s). Accurate structure prediction on this timescale makes possible avenues of investigation that were previously infeasible. As a demonstration of IgFold's capabilities, we predicted structures for 1.4 million paired antibody sequences, providing structural insights to 500-fold more antibodies than have experimentally determined structures.

The validation of artificial anti-monkeypox antibodies by in silico and experimental approaches

As a result of smallpox immunization programs that ended more than 40 years ago, a significant portion of the world's population is not immune. Moreover, due to the lack of anti-monkeypox drugs and vaccines against monkeypox, the spread of this virus may be the beginning of another challenge. In this study, novel antibodies against monkeypox virus were modeled based on a heavy chain of human antibody and a small peptide fragment. Docking of modeled antibodies with C19L protein showed the range of docking energy, and root-mean-square deviation (RMSD) was from -124 to -154 kcal/mL and 4-6 angstrom, respectively. Also, docking of modeled antibodies-C19L complex with gamma Fc receptor type I illustrated the range of docking energy, and RMSD was from -132 to -155 kcal/ml and 5-7 angstrom, respectively. Moreover, molecular dynamics simulation showed that antibody 62 had the highest stability with the lowest energy level and RMSD. Interestingly, no modeled antibodies had immunogenicity, allergenicity, and toxicity. Although all of them had good stability, only antibodies 25, 28, 54, and 62 had a half-life of >10 h. Moreover, the interaction between C19L protein and anti-C19L antibodies (wild-type and synthetic) was evaluated by the SPR method. We found that K_D in synthetic antibodies was lower than wild antibody. In terms of δH°, TδS°, and δG°, the results were consistent with binding parameters. Here, the lowest value of thermodynamic parameters was obtained for antibody 62. These data show that the synthetic antibodies, especially antibody 62, had a higher affinity than the wild-type antibody.

De novo design and Rosetta-based assessment of high-affinity antibody variable regions (Fv) against the SARS-CoV-2 spike receptor binding domain (RBD)

The continued emergence of new SARS-CoV-2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high-affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent-exposed hACE2-binding residues of the SARS-CoV-2 spike receptor binding domain (RBD) protein using the software tool OptMAVEn-2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9-mer library of "human Abs" based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade-off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced by OptMAVEn-2.0 using a Rosetta-based approach. We used Rosetta SnugDock for local docking of the designs to evaluate their potential to bind the spike RBD and performed "forward folding" with DeepAb to assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS-CoV-2 variants or other antigenic targets.

Toward real-world automated antibody design with combinatorial Bayesian optimization

Antibodies are multimeric proteins capable of highly specific molecular recognition. The complementarity determining region 3 of the antibody variable heavy chain (CDRH3) often dominates antigen-binding specificity. Hence, it is a priority to design optimal antigen-specific CDRH3 to develop therapeutic antibodies. The combinatorial structure of CDRH3 sequences makes it impossible to query binding-affinity oracles exhaustively. Moreover, antibodies are expected to have high target specificity and developability. Here, we present AntBO, a combinatorial Bayesian optimization framework utilizing a CDRH3 trust region for an in silico design of antibodies with favorable developability scores. The in silico experiments on 159 antigens demonstrate that AntBO is a step toward practically viable in vitro antibody design. In under 200 calls to the oracle, AntBO suggests antibodies outperforming the best binding sequence from 6.9 million experimentally obtained CDRH3s. Additionally, AntBO finds very-high-affinity CDRH3 in only 38 protein designs while requiring no domain knowledge.

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.

Prediction of Antibody-Antigen Binding via Machine Learning: Development of Data Sets and Evaluation of Methods

BACKGROUND

The mammalian immune system is able to generate antibodies against a huge variety of antigens, including bacteria, viruses, and toxins. The ultradeep DNA sequencing of rearranged immunoglobulin genes has considerable potential in furthering our understanding of the immune response, but it is limited by the lack of a high-throughput, sequence-based method for predicting the antigen(s) that a given immunoglobulin recognizes.

OBJECTIVE

As a step toward the prediction of antibody-antigen binding from sequence data alone, we aimed to compare a range of machine learning approaches that were applied to a collated data set of antibody-antigen pairs in order to predict antibody-antigen binding from sequence data.

METHODS

Data for training and testing were extracted from the Protein Data Bank and the Coronavirus Antibody Database, and additional antibody-antigen pair data were generated by using a molecular docking protocol. Several machine learning methods, including the weighted nearest neighbor method, the nearest neighbor method with the BLOSUM62 matrix, and the random forest method, were applied to the problem.

RESULTS

The final data set contained 1157 antibodies and 57 antigens that were combined in 5041 antibody-antigen pairs. The best performance for the prediction of interactions was obtained by using the nearest neighbor method with the BLOSUM62 matrix, which resulted in around 82% accuracy on the full data set. These results provide a useful frame of reference, as well as protocols and considerations, for machine learning and data set creation in the prediction of antibody-antigen binding.

CONCLUSIONS

Several machine learning approaches were compared to predict antibody-antigen interaction from protein sequences. Both the data set (in CSV format) and the machine learning program (coded in Python) are freely available for download on GitHub.

Characterization of Neutralizing Human Anti-Tetanus Monoclonal Antibodies Produced by Stable Cell Lines

Tetanus toxin (TeNT) is produced by C. tetani, a spore-forming bacillus broadly spread in the environment. Although an inexpensive and safe vaccine is available, tetanus persists because of a lack of booster shots and variable responses to vaccines due to immunocompromised status or age-decreased immune surveillance. Tetanus is most prevalent in low- and medium-income countries, where it remains a health problem. Neutralizing monoclonal antibodies (mAbs) can prevent the severity of illness and death caused by C. tetani infection. We identified a panel of mAbs that bind to TeNT, some of which were investigated in a preclinical assay, showing that a trio of mAbs that bind to different sites of TeNT can neutralize the toxin and prevent symptoms and death in mice. We also identified two mAbs that can impair the binding of TeNT to the GT1b ganglioside receptor in neurons. In this work, to generate a series of cell lines, we constructed vectors containing sequences encoding heavy and light constant regions that can receive the paired variable regions resulting from PCRs of human B cells. In this way, we generated stable cell lines for five mAbs and compared and characterized the antibody produced in large quantities, enabling the characterization experiments. We present the results regarding the cell growth and viability in a fed-batch culture, titer measurement, and specific productivity estimation. The affinity of purified mAbs was analyzed by kinetics and under steady-state conditions, as three mAbs could not dissociate from TeNT within 36,000 s. The binding of mAbs to TeNT was confirmed by ELISA and inhibition of toxin binding to GT1b. The use of the mAbs mixture confirmed the individual mAb contribution to inhibition. We also analyzed the binding of mAbs to FcγR by surface plasmon resonance (SPR) and the glycan composition. Molecular docking analyses showed the binding site of an anti-tetanus mAb.

Recent Progress in the Discovery and Development of Monoclonal Antibodies against Viral Infections

Monoclonal antibodies (mAbs), the new revolutionary class of medications, are fast becoming tools against various diseases thanks to a unique structure and function that allow them to bind highly specific targets or receptors. These specialized proteins can be produced in large quantities via the hybridoma technique introduced in 1975 or by means of modern technologies. Additional methods have been developed to generate mAbs with new biological properties such as humanized, chimeric, or murine. The inclusion of mAbs in therapeutic regimens is a major medical advance and will hopefully lead to significant improvements in infectious disease management. Since the first therapeutic mAb, muromonab-CD3, was approved by the U.S. Food and Drug Administration (FDA) in 1986, the list of approved mAbs and their clinical indications and applications have been proliferating. New technologies have been developed to modify the structure of mAbs, thereby increasing efficacy and improving delivery routes. Gene delivery technologies, such as non-viral synthetic plasmid DNA and messenger RNA vectors (DMabs or mRNA-encoded mAbs), built to express tailored mAb genes, might help overcome some of the challenges of mAb therapy, including production restrictions, cold-chain storage, transportation requirements, and expensive manufacturing and distribution processes. This paper reviews some of the recent developments in mAb discovery against viral infections and illustrates how mAbs can help to combat viral diseases and outbreaks.

Machine-designed biotherapeutics: opportunities, feasibility and advantages of deep learning in computational antibody discovery

Antibodies are versatile molecular binders with an established and growing role as therapeutics. Computational approaches to developing and designing these molecules are being increasingly used to complement traditional lab-based processes. Nowadays, in silico methods fill multiple elements of the discovery stage, such as characterizing antibody–antigen interactions and identifying developability liabilities. Recently, computational methods tackling such problems have begun to follow machine learning paradigms, in many cases deep learning specifically. This paradigm shift offers improvements in established areas such as structure or binding prediction and opens up new possibilities such as language-based modeling of antibody repertoires or machine-learning-based generation of novel sequences. In this review, we critically examine the recent developments in (deep) machine learning approaches to therapeutic antibody design with implications for fully computational antibody design.

Structural Features of Antibody-Peptide Recognition

Antibody recognition of antigens is a critical element of adaptive immunity. One key class of antibody-antigen complexes is comprised of antibodies targeting linear epitopes of proteins, which in some cases are conserved elements of viruses and pathogens of relevance for vaccine design and immunotherapy. Here we report a detailed analysis of the structural and interface features of this class of complexes, based on a set of nearly 200 nonredundant high resolution antibody-peptide complex structures that were assembled from the Protein Data Bank. We found that antibody-bound peptides adopt a broad range of conformations, often displaying limited secondary structure, and that the same peptide sequence bound by different antibodies can in many cases exhibit varying conformations. Propensities of contacts with antibody loops and extent of antibody binding conformational changes were found to be broadly similar to those for antibodies in complex with larger protein antigens. However, antibody-peptide interfaces showed lower buried surface areas and fewer hydrogen bonds than antibody-protein antigen complexes, while calculated binding energy per buried interface area was found to be higher on average for antibody-peptide interfaces, likely due in part to a greater proportion of buried hydrophobic residues and higher shape complementarity. This dataset and these observations can be of use for future studies focused on this class of interactions, including predictive computational modeling efforts and the design of antibodies or epitope-based vaccine immunogens.

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.

Optimal target saturation of ligand-blocking anti-GITR antibody IBI37G5 dictates FcγR-independent GITR agonism and antitumor activity

Glucocorticoid-induced tumor necrosis factor receptor (GITR) is a co-stimulatory receptor and an important target for cancer immunotherapy. We herein present a potent FcγR-independent GITR agonist IBI37G5 that can effectively activate effector T cells and synergize with anti-programmed death 1 (PD1) antibody to eradicate established tumors. IBI37G5 depends on both antibody bivalency and GITR homo-dimerization for efficient receptor cross-linking. Functional analyses reveal bell-shaped dose responses due to the unique 2:2 antibody-receptor stoichiometry required for GITR activation. Antibody self-competition is observed after concentration exceeded that of 100% receptor occupancy (RO), which leads to antibody monovalent binding and loss of activity. Retrospective pharmacokinetics/pharmacodynamics analysis demonstrates that the maximal efficacy is achieved at medium doses with drug exposure near saturating GITR occupancy during the dosing cycle. Finally, we propose an alternative dose-finding strategy that does not rely on the traditional maximal tolerated dose (MTD)-based paradigm but instead on utilizing the RO-function relations as biomarker to guide the clinical translation of GITR and similar co-stimulatory agonists.

Computational Structure Prediction for Antibody-Antigen Complexes From Hydrogen-Deuterium Exchange Mass Spectrometry: Challenges and Outlook

Although computational structure prediction has had great successes in recent years, it regularly fails to predict the interactions of large protein complexes with residue-level accuracy, or even the correct orientation of the protein partners. The performance of computational docking can be notably enhanced by incorporating experimental data from structural biology techniques. A rapid method to probe protein-protein interactions is hydrogen-deuterium exchange mass spectrometry (HDX-MS). HDX-MS has been increasingly used for epitope-mapping of antibodies (Abs) to their respective antigens (Ags) in the past few years. In this paper, we review the current state of HDX-MS in studying protein interactions, specifically Ab-Ag interactions, and how it has been used to inform computational structure prediction calculations. Particularly, we address the limitations of HDX-MS in epitope mapping and techniques and protocols applied to overcome these barriers. Furthermore, we explore computational methods that leverage HDX-MS to aid structure prediction, including the computational simulation of HDX-MS data and the combination of HDX-MS and protein docking. We point out challenges in interpreting and incorporating HDX-MS data into Ab-Ag complex docking and highlight the opportunities they provide to build towards a more optimized hybrid method, allowing for more reliable, high throughput epitope identification.

AbAdapt: an adaptive approach to predicting antibody-antigen complex structures from sequence

Motivation

The scoring of antibody-antigen docked poses starting from unbound homology models has not been systematically optimized for a large and diverse set of input sequences.

Results

To address this need, we have developed AbAdapt, a webserver that accepts antibody and antigen sequences, models their 3D structures, predicts epitope and paratope, and then docks the modeled structures using two established docking engines (Piper and Hex). Each of the key steps has been optimized by developing and training new machine-learning models. The sequences from a diverse set of 622 antibody-antigen pairs with known structure were used as inputs for leave-one-out cross-validation. The final set of cluster representatives included at least one 'Adequate' pose for 550/622 (88.4%) of the queries. The median (interquartile range) ranks of these 'Adequate' poses were 22 (5-77). Similar results were obtained on a holdout set of 100 unrelated antibody-antigen pairs. When epitopes were repredicted using docking-derived features for specific antibodies, the median ROC AUC increased from 0.679 to 0.720 in cross-validation and from 0.694 to 0.730 in the holdout set.

Availability and implementation

AbAdapt and related data are available at https://sysimm.org/abadapt/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Progress and challenges for the machine learning-based design of fit-for-purpose monoclonal antibodies

Although the therapeutic efficacy and commercial success of monoclonal antibodies (mAbs) are tremendous, the design and discovery of new candidates remain a time and cost-intensive endeavor. In this regard, progress in the generation of data describing antigen binding and developability, computational methodology, and artificial intelligence may pave the way for a new era of in silico on-demand immunotherapeutics design and discovery. Here, we argue that the main necessary machine learning (ML) components for an in silico mAb sequence generator are: understanding of the rules of mAb-antigen binding, capacity to modularly combine mAb design parameters, and algorithms for unconstrained parameter-driven in silico mAb sequence synthesis. We review the current progress toward the realization of these necessary components and discuss the challenges that must be overcome to allow the on-demand ML-based discovery and design of fit-for-purpose mAb therapeutic candidates.

Ensuring scientific reproducibility in bio-macromolecular modeling via extensive, automated benchmarks

Each year vast international resources are wasted on irreproducible research. The scientific community has been slow to adopt standard software engineering practices, despite the increases in high-dimensional data, complexities of workflows, and computational environments. Here we show how scientific software applications can be created in a reproducible manner when simple design goals for reproducibility are met. We describe the implementation of a test server framework and 40 scientific benchmarks, covering numerous applications in Rosetta bio-macromolecular modeling. High performance computing cluster integration allows these benchmarks to run continuously and automatically. Detailed protocol captures are useful for developers and users of Rosetta and other macromolecular modeling tools. The framework and design concepts presented here are valuable for developers and users of any type of scientific software and for the scientific community to create reproducible methods. Specific examples highlight the utility of this framework, and the comprehensive documentation illustrates the ease of adding new tests in a matter of hours.

Paratope states in solution improve structure prediction and docking

Structure-based antibody design and accurate predictions of antibody-antigen interactions remain major challenges in computational biology. By using molecular dynamics simulations, we show that a single static X-ray structure is not sufficient to identify determinants of antibody-antigen recognition. Here, we investigate antibodies that undergo substantial conformational changes upon antigen binding and have been classified as difficult cases in an extensive benchmark for antibody-antigen docking. We present thermodynamics and transition kinetics of these conformational rearrangements and show that paratope states can be used to improve antibody-antigen docking. By using the unbound antibody X-ray structure as starting structure for molecular dynamics simulations, we retain a binding competent conformation substantially different to the unbound antibody X-ray structure. We also observe that the kinetically dominant antibody paratope conformations are chosen by the bound antigen conformation with the highest probability. Thus, we show that paratope states in solution can improve antibody-antigen docking and structure prediction.

Design of peptides with high affinity binding to a monoclonal antibody as a basis for immunotherapy

About half of the US population is sensitized to one or more allergens, as found by a National Health and Nutrition Examination Survey (NHANES). The most common treatment for seasonal allergic responses is the daily use of oral antihistamines, which can control some of the symptoms, but are not effective for nasal congestion, and can be debilitating in many patients. Peptide immunotherapy is a promising new approach to treat allergic airway diseases. The small size of the immunogens cannot lead to an unwanted allergic reaction in sensitized patients, and the production of peptides with sufficient amounts for immunotherapy is time- and cost-effective. However, it is not known what peptides are the most effective for an immunotherapy of allergens. We previously produced a unique monoclonal antibody (mAb) E58, which can inhibit the binding of multiple groups of mAbs and human IgEs from patients affected by the major group 1 allergens of ragweed (Amb a 1) and conifer pollens (Jun a 1, Cup s 1, and Cry j 1). Here, we demonstrated that a combined approach, starting from two linear E58 epitopes of the tree pollen allergen Jun a 1 and the ragweed pollen allergen Amb a 1, and residue modifications suggested by molecular docking calculations and peptide design could identify a large number of high affinity binding peptides. We propose that this combined experimental and computational approach by structural analysis of linear IgE epitopes and peptide design, can lead to potential new candidates for peptide immunotherapy.

NbX: Machine Learning-Guided Re-Ranking of Nanobody-Antigen Binding Poses

Modeling the binding pose of an antibody is a prerequisite to structure-based affinity maturation and design. Without knowing a reliable binding pose, the subsequent structural simulation is largely futile. In this study, we have developed a method of machine learning-guided re-ranking of antigen binding poses of nanobodies, the single-domain antibody which has drawn much interest recently in antibody drug development. We performed a large-scale self-docking experiment of nanobody–antigen complexes. By training a decision tree classifier through mapping a feature set consisting of energy, contact and interface property descriptors to a measure of their docking quality of the refined poses, significant improvement in the median ranking of native-like nanobody poses by was achieved eightfold compared with ClusPro and an established deep 3D CNN classifier of native protein–protein interaction. We further interpreted our model by identifying features that showed relatively important contributions to the prediction performance. This study demonstrated a useful method in improving our current ability in pose prediction of nanobodies.

RAPID: A Rep-Seq Dataset Analysis Platform With an Integrated Antibody Database

The antibody repertoire is a critical component of the adaptive immune system and is believed to reflect an individual’s immune history and current immune status. Delineating the antibody repertoire has advanced our understanding of humoral immunity, facilitated antibody discovery, and showed great potential for improving the diagnosis and treatment of disease. However, no tool to date has effectively integrated big Rep-seq data and prior knowledge of functional antibodies to elucidate the remarkably diverse antibody repertoire. We developed a Rep-seq dataset Analysis Platform with an Integrated antibody Database (RAPID; https://rapid.zzhlab.org/), a free and web-based tool that allows researchers to process and analyse Rep-seq datasets. RAPID consolidates 521 WHO-recognized therapeutic antibodies, 88,059 antigen- or disease-specific antibodies, and 306 million clones extracted from 2,449 human IGH Rep-seq datasets generated from individuals with 29 different health conditions. RAPID also integrates a standardized Rep-seq dataset analysis pipeline to enable users to upload and analyse their datasets. In the process, users can also select set of existing repertoires for comparison. RAPID automatically annotates clones based on integrated therapeutic and known antibodies, and users can easily query antibodies or repertoires based on sequence or optional keywords. With its powerful analysis functions and rich set of antibody and antibody repertoire information, RAPID will benefit researchers in adaptive immune studies.

A high-affinity human TCR-like antibody detects celiac disease gluten peptide-MHC complexes and inhibits T cell activation

Antibodies specific for peptides bound to human leukocyte antigen (HLA) molecules are valuable tools for studies of antigen presentation and may have therapeutic potential. Here, we generated human T cell receptor (TCR)-like antibodies toward the immunodominant signature gluten epitope DQ2.5-glia-α2 in celiac disease (CeD). Phage display selection combined with secondary targeted engineering was used to obtain highly specific antibodies with picomolar affinity. The crystal structure of a Fab fragment of the lead antibody 3.C11 in complex with HLA-DQ2.5:DQ2.5-glia-α2 revealed a binding geometry and interaction mode highly similar to prototypic TCRs specific for the same complex. Assessment of CeD biopsy material confirmed disease specificity and reinforced the notion that abundant plasma cells present antigen in the inflamed CeD gut. Furthermore, 3.C11 specifically inhibited activation and proliferation of gluten-specific CD4⁺ T cells in vitro and in HLA-DQ2.5 humanized mice, suggesting a potential for targeted intervention without compromising systemic immunity.

IsAb: a computational protocol for antibody design

The design of therapeutic antibodies has attracted a large amount of attention over the years. Antibodies are widely used to treat many diseases due to their high efficiency and low risk of adverse events. However, the experimental methods of antibody design are time-consuming and expensive. Although computational antibody design techniques have had significant advances in the past years, there are still some challenges that need to be solved, such as the flexibility of antigen structure, the lack of antibody structural data and the absence of standard antibody design protocol. In the present work, we elaborated on an in silico antibody design protocol for users to easily perform computer-aided antibody design. First, the Rosetta web server will be applied to generate the 3D structure of query antibodies if there is no structural information available. Then, two-step docking will be used to identify the binding pose of an antibody-antigen complex when the binding information is unknown. ClusPro is the first method to be used to conduct the global docking, and SnugDock is applied for the local docking. Sequentially, based on the predicted binding poses, in silico alanine scanning will be used to predict the potential hotspots (or key residues). Finally, computational affinity maturation protocol will be used to modify the structure of antibodies to theoretically increase their affinity and stability, which will be further validated by the bioassays in the future. As a proof of concept, we redesigned antibody D44.1 and compared it with previously reported data in order to validate IsAb protocol. To further illustrate our proposed protocol, we used cemiplimab antibody, a PD-1 checkpoint inhibitor, as an example to showcase a step-by-step tutorial.

Robustification of RosettaAntibody and Rosetta SnugDock

In recent years, the observed antibody sequence space has grown exponentially due to advances in high-throughput sequencing of immune receptors. The rise in sequences has not been mirrored by a rise in structures, as experimental structure determination techniques have remained low-throughput. Computational modeling, however, has the potential to close the sequence–structure gap. To achieve this goal, computational methods must be robust, fast, easy to use, and accurate. Here we report on the latest advances made in RosettaAntibody and Rosetta SnugDock—methods for antibody structure prediction and antibody–antigen docking. We simplified the user interface, expanded and automated the template database, generalized the kinematics of antibody–antigen docking (which enabled modeling of single-domain antibodies) and incorporated new loop modeling techniques. To evaluate the effects of our updates on modeling accuracy, we developed rigorous tests under a new scientific benchmarking framework within Rosetta. Benchmarking revealed that more structurally similar templates could be identified in the updated database and that SnugDock broadened its applicability without losing accuracy. However, there are further advances to be made, including increasing the accuracy and speed of CDR-H3 loop modeling, before computational approaches can accurately model any antibody.

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.

An expanded benchmark for antibody-antigen docking and affinity prediction reveals insights into antibody recognition determinants

Accurate predictive modeling of antibody-antigen complex structures and structure-based antibody design remain major challenges in computational biology, with implications for biotherapeutics, immunity, and vaccines. Through a systematic search for high-resolution structures of antibody-antigen complexes and unbound antibody and antigen structures, in conjunction with identification of experimentally determined binding affinities, we have assembled a non-redundant set of test cases for antibody-antigen docking and affinity prediction. This benchmark more than doubles the number of antibody-antigen complexes and corresponding affinities available in our previous benchmarks, providing an unprecedented view of the determinants of antibody recognition and insights into molecular flexibility. Initial assessments of docking and affinity prediction tools highlight the challenges posed by this diverse set of cases, which includes camelid nanobodies, therapeutic monoclonal antibodies, and broadly neutralizing antibodies targeting viral glycoproteins. This dataset will enable development of advanced predictive modeling and design methods for this therapeutically relevant class of protein-protein interactions.

Deimmunization of protein therapeutics - Recent advances in experimental and computational epitope prediction and deletion

Biotherapeutics, and antimicrobial proteins in particular, are of increasing interest for human medicine. An important challenge in the development of such therapeutics is their potential immunogenicity, which can induce production of anti-drug-antibodies, resulting in altered pharmacokinetics, reduced efficacy, and potentially severe anaphylactic or hypersensitivity reactions. For this reason, the development and application of effective deimmunization methods for protein drugs is of utmost importance. Deimmunization may be achieved by unspecific shielding approaches, which include PEGylation, fusion to polypeptides (e.g., XTEN or PAS), reductive methylation, glycosylation, and polysialylation. Alternatively, the identification of epitopes for T cells or B cells and their subsequent deletion through site-directed mutagenesis represent promising deimmunization strategies and can be accomplished through either experimental or computational approaches. This review highlights the most recent advances and current challenges in the deimmunization of protein therapeutics, with a special focus on computational epitope prediction and deletion tools.

Discovery of Marburg virus neutralizing antibodies from virus-naïve human antibody repertoires using large-scale structural predictions

Marburg virus (MARV) disease is lethal, with fatality rates up to 90%. Neutralizing antibodies (Abs) are promising drug candidates to prevent or treat the disease. Current efforts are focused in part on vaccine development to induce such MARV-neutralizing Abs. We analyzed the antibody repertoire from healthy unexposed and previously MARV-infected individuals to assess if naïve repertoires contain suitable precursor antibodies that could become neutralizing with a limited set of somatic mutations. We computationally searched the human Ab variable gene repertoire for predicted structural homologs of the neutralizing Ab MR78 that is specific to the receptor binding site (RBS) of MARV glycoprotein (GP). Eight Ab heavy-chain complementarity determining region 3 (HCDR3) loops from MARV-naïve individuals and one from a previously MARV-infected individual were selected for testing as HCDR3 loop chimeras on the MR78 Ab framework. Three of these chimerized antibodies bound to MARV GP. We then tested a full-length native Ab heavy chain encoding the same 17-residue-long HCDR3 loop that bound to the MARV GP the best among the chimeric Abs tested. Despite only 57% amino acid sequence identity, the Ab from a MARV-naïve donor recognized MARV GP and possessed neutralizing activity against the virus. Crystallization of both chimeric and full-length native heavy chain-containing Abs provided structural insights into the mechanism of binding for these types of Abs. Our work suggests that the MARV GP RBS is a promising candidate for epitope-focused vaccine design to induce neutralizing Abs against MARV.

Novel murine mAbs define specific and cross-reactive epitopes on the latex profilin panallergen Hev b 8

The production of specific antibodies able to recognize allergens from different sources or block interactions between allergens and antibodies mediating allergic reactions is crucial for developing successful tools for diagnostics and therapeutics. Panallergens are highly conserved proteins present in widely different species, implicated in relevant cross-reactions. The panallergen latex profilin (Hev b 8) has been associated with the latex-food-pollen syndrome. We generated five monoclonal IgGs and one IgE from murine hybridomas against recombinant Hev b 8 and evaluated their interaction with this allergen using ELISA and biolayer interferometry (BLI). Affinity purified mAbs exhibited high binding affinities towards rHev b 8, with K_D1 values ranging from 10^-10 M to 10^-11 M. Some of these antibodies also recognized the recombinant profilins from maize and tomato (Zea m 12 and Sola l 1), and the ash tree pollen (Fra e 2). Competition ELISA demonstrated that some mAb pairs could bind simultaneously to rHev b 8. Using BLI, we detected competitive, non-competitive, and partial-competition interactions between pairs of mAbs with rHev b 8, suggesting the existence of at least two non-overlapping epitopes on the surface of this allergen. Three-dimensional models of the Fv of 1B4 and 2D10 IgGs and docking simulations of these Fvs with rHev b 8 revealed these epitopes. Furthermore, these two mAbs inhibited the interaction of polyclonal IgE and IgG4 antibodies from profilin-allergic patients with rHev b 8, indicating that the mAbs and the antibodies present in sera from allergic patients bind to overlapping epitopes on the allergen. These mAbs can be useful tools for immune-localization studies, immunoassay development, or standardization of allergenic products.

Computational approaches to therapeutic antibody design: established methods and emerging trends

Antibodies are proteins that recognize the molecular surfaces of potentially noxious molecules to mount an adaptive immune response or, in the case of autoimmune diseases, molecules that are part of healthy cells and tissues. Due to their binding versatility, antibodies are currently the largest class of biotherapeutics, with five monoclonal antibodies ranked in the top 10 blockbuster drugs. Computational advances in protein modelling and design can have a tangible impact on antibody-based therapeutic development. Antibody-specific computational protocols currently benefit from an increasing volume of data provided by next generation sequencing and application to related drug modalities based on traditional antibodies, such as nanobodies. Here we present a structured overview of available databases, methods and emerging trends in computational antibody analysis and contextualize them towards the engineering of candidate antibody therapeutics.

Characterizing Epitope Binding Regions of Entire Antibody Panels by Combining Experimental and Computational Analysis of Antibody: Antigen Binding Competition

Vaccines and immunotherapies depend on the ability of antibodies to sensitively and specifically recognize particular antigens and specific epitopes on those antigens. As such, detailed characterization of antibody-antigen binding provides important information to guide development. Due to the time and expense required, high-resolution structural characterization techniques are typically used sparingly and late in a development process. Here, we show that antibody-antigen binding can be characterized early in a process for whole panels of antibodies by combining experimental and computational analyses of competition between monoclonal antibodies for binding to an antigen. Experimental "epitope binning" of monoclonal antibodies uses high-throughput surface plasmon resonance to reveal which antibodies compete, while a new complementary computational analysis that we call "dock binning" evaluates antibody-antigen docking models to identify why and where they might compete, in terms of possible binding sites on the antigen. Experimental and computational characterization of the identified antigenic hotspots then enables the refinement of the competitors and their associated epitope binding regions on the antigen. While not performed at atomic resolution, this approach allows for the group-level identification of functionally related monoclonal antibodies (i.e., communities) and identification of their general binding regions on the antigen. By leveraging extensive epitope characterization data that can be readily generated both experimentally and computationally, researchers can gain broad insights into the basis for antibody-antigen recognition in wide-ranging vaccine and immunotherapy discovery and development programs.

Novel sampling strategies and a coarse-grained score function for docking homomers, flexible heteromers, and oligosaccharides using Rosetta in CAPRI rounds 37-45

Critical Assessment of PRediction of Interactions (CAPRI) rounds 37 through 45 introduced larger complexes, new macromolecules, and multistage assemblies. For these rounds, we used and expanded docking methods in Rosetta to model 23 target complexes. We successfully predicted 14 target complexes and recognized and refined near-native models generated by other groups for two further targets. Notably, for targets T110 and T136, we achieved the closest prediction of any CAPRI participant. We created several innovative approaches during these rounds. Since round 39 (target 122), we have used the new RosettaDock 4.0, which has a revamped coarse-grained energy function and the ability to perform conformer selection during docking with hundreds of pregenerated protein backbones. Ten of the complexes had some degree of symmetry in their interactions, so we tested Rosetta SymDock, realized its shortcomings, and developed the next-generation symmetric docking protocol, SymDock2, which includes docking of multiple backbones and induced-fit refinement. Since the last CAPRI assessment, we also developed methods for modeling and designing carbohydrates in Rosetta, and we used them to successfully model oligosaccharide-protein complexes in round 41. Although the results were broadly encouraging, they also highlighted the pressing need to invest in (a) flexible docking algorithms with the ability to model loop and linker motions and in (b) new sampling and scoring methods for oligosaccharide-protein interactions.

Methods for sequence and structural analysis of B and T cell receptor repertoires

B cell receptors (BCRs) and T cell receptors (TCRs) make up an essential network of defense molecules that, collectively, can distinguish self from non-self and facilitate destruction of antigen-bearing cells such as pathogens or tumors. The analysis of BCR and TCR repertoires plays an important role in both basic immunology as well as in biotechnology. Because the repertoires are highly diverse, specialized software methods are needed to extract meaningful information from BCR and TCR sequence data. Here, we review recent developments in bioinformatics tools for analysis of BCR and TCR repertoires, with an emphasis on those that incorporate structural features. After describing the recent sequencing technologies for immune receptor repertoires, we survey structural modeling methods for BCR and TCRs, along with methods for clustering such models. We review downstream analyses, including BCR and TCR epitope prediction, antibody-antigen docking and TCR-peptide-MHC Modeling. We also briefly discuss molecular dynamics in this context.