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

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.

mmCSM-AB: guiding rational antibody engineering through multiple point mutations

While antibodies are becoming an increasingly important therapeutic class, especially in personalized medicine, their development and optimization has been largely through experimental exploration. While there have been many efforts to develop computational tools to guide rational antibody engineering, most approaches are of limited accuracy when applied to antibody design, and have largely been limited to analysing a single point mutation at a time. To overcome this gap, we have curated a dataset of 242 experimentally determined changes in binding affinity upon multiple point mutations in antibody-target complexes (89 increasing and 153 decreasing binding affinity). Here, we have shown that by using our graph-based signatures and atomic interaction information, we can accurately analyse the consequence of multi-point mutations on antigen binding affinity. Our approach outperformed other available tools across cross-validation and two independent blind tests, achieving Pearson's correlations of up to 0.95. We have implemented our new approach, mmCSM-AB, as a web-server that can help guide the process of affinity maturation in antibody design. mmCSM-AB is freely available at http://biosig.unimelb.edu.au/mmcsm_ab/.

Learning context-aware structural representations to predict antigen and antibody binding interfaces

MOTIVATION

Understanding how antibodies specifically interact with their antigens can enable better drug and vaccine design, as well as provide insights into natural immunity. Experimental structural characterization can detail the 'ground truth' of antibody-antigen interactions, but computational methods are required to efficiently scale to large-scale studies. To increase prediction accuracy as well as to provide a means to gain new biological insights into these interactions, we have developed a unified deep learning-based framework to predict binding interfaces on both antibodies and antigens.

RESULTS

Our framework leverages three key aspects of antibody-antigen interactions to learn predictive structural representations: (i) since interfaces are formed from multiple residues in spatial proximity, we employ graph convolutions to aggregate properties across local regions in a protein; (ii) since interactions are specific between antibody-antigen pairs, we employ an attention layer to explicitly encode the context of the partner; (iii) since more data are available for general protein-protein interactions, we employ transfer learning to leverage this data as a prior for the specific case of antibody-antigen interactions. We show that this single framework achieves state-of-the-art performance at predicting binding interfaces on both antibodies and antigens, and that each of its three aspects drives additional improvement in the performance. We further show that the attention layer not only improves performance, but also provides a biologically interpretable perspective into the mode of interaction.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available on github at https://github.com/vamships/PECAN.git.

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.

In silico analysis of SARS-CoV-2 spike glycoprotein and insights into antibody binding

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in Wuhan, China in December 2019. Since then, COVID-19, the disease caused by SARS-CoV-2, has become a rapidly spreading pandemic that has reached most countries in the world. So far, there are no vaccines or therapeutics to fight this virus. Here, I present an in silico analysis of the virus spike glycoprotein (recently determined at atomic resolution) and provide insights into how antibodies against the 2002 virus SARS-CoV might be modified to neutralize SARS-CoV-2. I ran docking experiments with Rosetta Dock to determine which substitutions in the 80R and m396 antibodies might improve the binding of these to SARS-CoV-2 and used molecular visualization and analysis software, including UCSF Chimera and Rosetta Dock, as well as other bioinformatics tools, including SWISS-MODEL. Supercomputers, including Bridges Large, Stampede and Frontera, were used for macromolecular assemblies and large scale analysis and visualization.

Molecular interactions between monoclonal oligomer-specific antibody 5E3 and its amyloid beta cognates

Oligomeric amyloid β (Aβ) is currently considered the most neurotoxic form of the Aβ peptide implicated in Alzheimer’s disease (AD). The molecular structures of the oligomers have remained mostly unknown due to their transient nature. As a result, the molecular mechanisms of interactions between conformation-specific antibodies and their Aβ oligomer (AβO) cognates are not well understood. A monoclonal conformation-specific antibody, m5E3, was raised against a structural epitope of Aβ oligomers. m5E3 binds to AβOs with high affinity, but not to Aβ monomers or fibrils. In this study, a computational model of the variable fragment (Fv) of the m5E3 antibody (Fv5E3) is introduced. We further employ docking and molecular dynamics simulations to determine the molecular details of the antibody-oligomer interactions, and to classify the AβOs as Fv5E3-positives and negatives, and to provide a rationale for the low affinity of Fv5E3 for fibrils. This information will help us to perform site-directed mutagenesis on the m5E3 antibody to improve its specificity and affinity toward oligomeric Aβ species. We also provide evidence for the possible capability of the m5E3 antibody to disaggregate AβOs and to fragment protofilaments.

How repertoire data are changing antibody science

Antibodies are vital proteins of the immune system that recognize potentially harmful molecules and initiate their removal. Mammals can efficiently create vast numbers of antibodies with different sequences capable of binding to any antigen with high affinity and specificity. Because they can be developed to bind to many disease agents, antibodies can be used as therapeutics. In an organism, after antigen exposure, antibodies specific to that antigen are enriched through clonal selection, expansion, and somatic hypermutation. The antibodies present in an organism therefore report on its immune status, describe its innate ability to deal with harmful substances, and reveal how it has previously responded. Next-generation sequencing technologies are being increasingly used to query the antibody, or B-cell receptor (BCR), sequence repertoire, and the amount of BCR data in public repositories is growing. The Observed Antibody Space database, for example, currently contains over a billion sequences from 68 different studies. Repertoires are available that represent both the naive state (i.e. antigen-inexperienced) and that after immunization. This wealth of data has created opportunities to learn more about our immune system. In this review, we discuss the many ways in which BCR repertoire data have been or could be exploited. We highlight its utility for providing insights into how the naive immune repertoire is generated and how it responds to antigens. We also consider how structural information can be used to enhance these data and may lead to more accurate depictions of the sequence space and to applications in the discovery of new therapeutics.

Characterization of a human monoclonal antibody generated from a B-cell specific for a prefusion-stabilized spike protein of Middle East respiratory syndrome coronavirus

Middle East respiratory syndrome coronavirus (MERS-CoV) causes severe respiratory infection and continues to infect humans, thereby contributing to a high mortality rate (34.3% in 2019). In the absence of an available licensed vaccine and antiviral agent, therapeutic human antibodies have been suggested as candidates for treatment. In this study, human monoclonal antibodies were isolated by sorting B cells from patient's PBMC cells with prefusion stabilized spike (S) probes and a direct immunoglobulin cloning strategy. We identified six receptor-binding domain (RBD)-specific and five S1 (non-RBD)-specific antibodies, among which, only the RBD-specific antibodies showed high neutralizing potency (IC50 0.006-1.787 μg/ml) as well as high affinity to RBD. Notably, passive immunization using a highly potent antibody (KNIH90-F1) at a relatively low dose (2 mg/kg) completely protected transgenic mice expressing human DPP4 against MERS-CoV lethal challenge. These results suggested that human monoclonal antibodies isolated by using the rationally designed prefusion MERS-CoV S probe could be considered potential candidates for the development of therapeutic and/or prophylactic antiviral agents for MERS-CoV human infection.

Modeling Antibody-Antigen Complexes by Information-Driven Docking

Antibodies are Y-shaped proteins essential for immune response. Their capability to recognize antigens with high specificity makes them excellent therapeutic targets. Understanding the structural basis of antibody-antigen interactions is therefore crucial for improving our ability to design efficient biological drugs. Computational approaches such as molecular docking are providing a valuable and fast alternative to experimental structural characterization for these complexes. We investigate here how information about complementarity-determining regions and binding epitopes can be used to drive the modeling process, and present a comparative study of four different docking software suites (ClusPro, LightDock, ZDOCK, and HADDOCK) providing specific options for antibody-antigen modeling. Their performance on a dataset of 16 complexes is reported. HADDOCK, which includes information to drive the docking, is shown to perform best in terms of both success rate and quality of the generated models in both the presence and absence of information about the epitope on the antigen.

BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes

Antibodies have become an indispensable tool for many biotechnological and clinical applications. They bind their molecular target (antigen) by recognizing a portion of its structure (epitope) in a highly specific manner. The ability to predict epitopes from antigen sequences alone is a complex task. Despite substantial effort, limited advancement has been achieved over the last decade in the accuracy of epitope prediction methods, especially for those that rely on the sequence of the antigen only. Here, we present BepiPred-2.0 (http://www.cbs.dtu.dk/services/BepiPred/), a web server for predicting B-cell epitopes from antigen sequences. BepiPred-2.0 is based on a random forest algorithm trained on epitopes annotated from antibody-antigen protein structures. This new method was found to outperform other available tools for sequence-based epitope prediction both on epitope data derived from solved 3D structures, and on a large collection of linear epitopes downloaded from the IEDB database. The method displays results in a user-friendly and informative way, both for computer-savvy and non-expert users. We believe that BepiPred-2.0 will be a valuable tool for the bioinformatics and immunology community.

Plasma Cells Are the Most Abundant Gluten Peptide MHC-expressing Cells in Inflamed Intestinal Tissues From Patients With Celiac Disease

BACKGROUND & AIMS

Development of celiac disease is believed to involve the transglutaminase-dependent response of CD4⁺ T cells toward deamidated gluten peptides in the intestinal mucosa of individuals with specific HLA-DQ haplotypes. We investigated the antigen presentation process during this mucosal immune response.

METHODS

We generated monoclonal antibodies (mAbs) specific for the peptide-MHC (pMHC) complex of HLA-DQ2.5 and the immunodominant gluten epitope DQ2.5-glia-α1a using phage display. We used these mAbs to assess gluten peptide presentation and phenotypes of presenting cells by flow cytometry and enzyme-linked immune absorbent spot (ELISPOT) in freshly prepared single-cell suspensions from intestinal biopsies from 40 patients with celiac disease (35 untreated and 5 on a gluten-free diet) as well as 18 subjects with confirmed noninflamed gut mucosa (controls, 12 presumed healthy, 5 undergoing pancreatoduodenectomy, and 1 with potential celiac disease).

RESULTS

Using the mAbs, we detected MHC complexes on cells from intestinal biopsies from patients with celiac disease who consume gluten, but not from patients on gluten-free diets. We found B cells and plasma cells to be the most abundant cells that present DQ2.5-glia-α1a in the inflamed mucosa. We identified a subset of plasma cells that expresses B-cell receptors (BCR) specific for gluten peptides or the autoantigen transglutaminase 2 (TG2). Expression of MHC class II (MHCII) was not restricted to these specific plasma cells in patients with celiac disease but was observed in an average 30% of gut plasma cells from patients and controls.

CONCLUSIONS

A population of plasma cells from intestinal biopsies of patients with celiac disease express MHCII; this is the most abundant cell type presenting the immunodominant gluten peptide DQ2.5-glia-α1a in the tissues from these patients. These results indicate that plasma cells in the gut can function as antigen-presenting cells and might promote and maintain intestinal inflammation in patients with celiac disease or other inflammatory disorders.

Can We Improve Vaccine Efficacy by Targeting T and B Cell Repertoire Convergence?

Traditional vaccine development builds on the assumption that healthy individuals have virtually unlimited antigen recognition repertoires of receptors in B cells and T cells [the B cell receptor (BCR) and TCR respectively]. However, there are indications that there are "holes" in the breadth of repertoire diversity, where no or few B or T cell are able to bind to a given antigen. Repertoire diversity may in these cases be a limiting factor for vaccine efficacy. Assuming that it is possible to predict which B and T cell receptors will respond to a given immunogen, vaccine strategies could be optimized and personalized. In addition, vaccine testing could be simplified if we could predict responses through sequencing BCR and TCRs. Bulk sequencing has shown putatively specific converging sequences after infection or vaccination. However, only single cell technologies have made it possible to capture the sequence of both heavy and light chains of a BCR or the alpha and beta chains the TCR. This has enabled the cloning of receptors and the functional validation of a predicted specificity. This review summarizes recent evidence of converging sequences in infectious diseases. Current and potential future applications of single cell technology in immune repertoire analysis are then discussed. Finally, possible short- and long- term implications for vaccine research are highlighted.

Mapping discontinuous epitopes for MRK-16, UIC2 and 4E3 antibodies to extracellular loops 1 and 4 of human P-glycoprotein

P-glycoprotein (P-gp), an ATP-dependent efflux pump, is associated with the development of multidrug resistance in cancer cells. Antibody-mediated blockade of human P-gp activity has been shown to overcome drug resistance by re-sensitizing resistant cancer cells to anticancer drugs. Despite the potential clinical application of this finding, the epitopes of the three human P-gp-specific monoclonal antibodies MRK-16, UIC2 and 4E3, which bind to the extracellular loops (ECLs) have not yet been mapped. By generating human-mouse P-gp chimeras, we mapped the epitopes of these antibodies to ECLs 1 and 4. We then identified key amino acids in these regions by replacing mouse residues with homologous human P-gp residues to recover binding of antibodies to the mouse P-gp. We found that changing a total of ten residues, five each in ECL1 and ECL4, was sufficient to recover binding of both MRK-16 and 4E3 antibodies, suggesting a common epitope. However, recovery of the conformation-sensitive UIC2 epitope required replacement of thirteen residues in ECL1 and the same five residues replaced in the ECL4 for MRK-16 and 4E3 binding. These results demonstrate that discontinuous epitopes for MRK-16, UIC2 and 4E3 are located in the same regions of ECL1 and 4 of the multidrug transporter.

In Silico Methods in Antibody Design

Antibody therapies with high efficiency and low toxicity are becoming one of the major approaches in antibody therapeutics. Based on high-throughput sequencing and increasing experimental structures of antibodies/antibody-antigen complexes, computational approaches can predict antibody/antigen structures, engineering the function of antibodies and design antibody-antigen complexes with improved properties. This review summarizes recent progress in the field of in silico design of antibodies, including antibody structure modeling, antibody-antigen complex prediction, antibody stability evaluation, and allosteric effects in antibodies and functions. We listed the cases in which these methods have helped experimental studies to improve the affinities and physicochemical properties of antibodies. We emphasized how the molecular dynamics unveiled the allosteric effects during antibody-antigen recognition and antibody-effector recognition.

Construction and Characterization of a Humanized Anti-Epstein-Barr Virus gp350 Antibody with Neutralizing Activity in Cell Culture

Acute Epstein-Barr virus (EBV) infection in immunosuppressed transplant patients can give rise to a malignant B-cell proliferation known as post-transplant lymphoproliferative disease (PTLD). The EBV major virion surface glycoprotein (gp)350 is a principal target of naturally occurring neutralizing antibodies and is viewed as the best target to prevent acute infection and PTLD in at-risk transplant recipients. We have constructed a humanized (hu) version of the murine anti-gp350 neutralizing monoclonal antibody 72a1. The hu72a1 IgG1 antibody displayed no significant anti-mouse activity, recognized both gp350 and its splice variant gp220 as well as a gp350 peptide that was shown to constitute the principal EBV gp350 neutralizing epitope when tested in immunoassays. Hu72a1 antibody blocked in vitro EBV infection of B cells at a level which equaled that of a mouse-human chimeric 72a1 antibody construct. This work provides a further structural and immunological understanding of the 72a1 antibody interaction with EBV gp350, and constitutes a launch point for future anti-EBV therapeutic antibodies designed to block EBV infection and prevent PTLD while eliminating the deleterious antigenic murine features of the original 72a1 antibody.

Mass spectrometric epitope mapping

Mass spectrometric epitope mapping has become a versatile method to precisely determine a soluble antigen's partial structure that directly interacts with an antibody in solution. Typical lengths of investigated antigens have increased up to several 100 amino acids while experimentally determined epitope peptides have decreased in length to on average 10-15 amino acids. Since the early 1990s more and more sophisticated methods have been developed and have forwarded a bouquet of suitable approaches for epitope mapping with immobilized, temporarily immobilized, and free-floating antibodies. While up to now monoclonal antibodies have been mostly used in epitope mapping experiments, the applicability of polyclonal antibodies has been proven. The antibody's resistance towards enzymatic proteolysis has been of key importance for the two mostly applied methods: epitope excision and epitope extraction. Sample consumption has dropped to low pmol amounts on both, the antigen and the antibody. While adequate in-solution sample handling has been most important for successful epitope mapping, mass spectrometric analysis has been found the most suitable read-out method from early on. The rapidity by which mass spectrometric epitope mapping nowadays is executed outperforms all alternative methods. Thus, it can be asserted that mass spectrometric epitope mapping has reached a state of maturity, which allows it to be used in any mass spectrometry laboratory. After 25 years of constant and steady improvements, its application to clinical samples, for example, for patient characterization and stratification, is anticipated in the near future. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 37:229-241, 2018.

How B-Cell Receptor Repertoire Sequencing Can Be Enriched with Structural Antibody Data

Next-generation sequencing of immunoglobulin gene repertoires (Ig-seq) allows the investigation of large-scale antibody dynamics at a sequence level. However, structural information, a crucial descriptor of antibody binding capability, is not collected in Ig-seq protocols. Developing systematic relationships between the antibody sequence information gathered from Ig-seq and low-throughput techniques such as X-ray crystallography could radically improve our understanding of antibodies. The mapping of Ig-seq datasets to known antibody structures can indicate structurally, and perhaps functionally, uncharted areas. Furthermore, contrasting naïve and antigenically challenged datasets using structural antibody descriptors should provide insights into antibody maturation. As the number of antibody structures steadily increases and more and more Ig-seq datasets become available, the opportunities that arise from combining the two types of information increase as well. Here, we review how these data types enrich one another and show potential for advancing our knowledge of the immune system and improving antibody engineering.

Computationally-driven identification of antibody epitopes

Understanding where antibodies recognize antigens can help define mechanisms of action and provide insights into progression of immune responses. We investigate the extent to which information about binding specificity implicitly encoded in amino acid sequence can be leveraged to identify antibody epitopes. In computationally-driven epitope localization, possible antibody–antigen binding modes are modeled, and targeted panels of antigen variants are designed to experimentally test these hypotheses. Prospective application of this approach to two antibodies enabled epitope localization using five or fewer variants per antibody, or alternatively, a six-variant panel for both simultaneously. Retrospective analysis of a variety of antibodies and antigens demonstrated an almost 90% success rate with an average of three antigen variants, further supporting the observation that the combination of computational modeling and protein design can reveal key determinants of antibody–antigen binding and enable efficient studies of collections of antibodies identified from polyclonal samples or engineered libraries.

Web-accessible molecular modeling with Rosetta: The Rosetta Online Server that Includes Everyone (ROSIE)

The Rosetta molecular modeling software package provides a large number of experimentally validated tools for modeling and designing proteins, nucleic acids, and other biopolymers, with new protocols being added continually. While freely available to academic users, external usage is limited by the need for expertise in the Unix command line environment. To make Rosetta protocols available to a wider audience, we previously created a web server called Rosetta Online Server that Includes Everyone (ROSIE), which provides a common environment for hosting web-accessible Rosetta protocols. Here we describe a simplification of the ROSIE protocol specification format, one that permits easier implementation of Rosetta protocols. Whereas the previous format required creating multiple separate files in different locations, the new format allows specification of the protocol in a single file. This new, simplified protocol specification has more than doubled the number of Rosetta protocols available under ROSIE. These new applications include pK_a determination, lipid accessibility calculation, ribonucleic acid redesign, protein-protein docking, protein-small molecule docking, symmetric docking, antibody docking, cyclic toxin docking, critical binding peptide determination, and mapping small molecule binding sites. ROSIE is freely available to academic users at http://rosie.rosettacommons.org.

In silico methods for design of biological therapeutics

It has been twenty years since the first rationally designed small molecule drug was introduced into the market. Since then, we have progressed from designing small molecules to designing biotherapeutics. This class of therapeutics includes designed proteins, peptides and nucleic acids that could more effectively combat drug resistance and even act in cases where the disease is caused because of a molecular deficiency. Computational methods are crucial in this design exercise and this review discusses the various elements of designing biotherapeutic proteins and peptides. Many of the techniques discussed here, such as the deterministic and stochastic design methods, are generally used in protein design. We have devoted special attention to the design of antibodies and vaccines. In addition to the methods for designing these molecules, we have included a comprehensive list of all biotherapeutics approved for clinical use. Also included is an overview of methods that predict the binding affinity, cell penetration ability, half-life, solubility, immunogenicity and toxicity of the designed therapeutics. Biotherapeutics are only going to grow in clinical importance and are set to herald a new generation of disease management and cure.

Structure-based cross-docking analysis of antibody-antigen interactions

Antibody–antigen interactions are critical to our immune response, and understanding the structure-based biophysical determinants for their binding specificity and affinity is of fundamental importance. We present a computational structure-based cross-docking study to test the identification of native antibody–antigen interaction pairs among cognate and non-cognate complexes. We picked a dataset of 17 antibody–antigen complexes of which 11 have both bound and unbound structures available, and we generated a representative ensemble of cognate and non-cognate complexes. Using the Rosetta interface score as a classifier, the cognate pair was the top-ranked model in 80% (14/17) of the antigen targets using bound monomer structures in docking, 35% (6/17) when using unbound, and 12% (2/17) when using the homology-modeled backbones to generate the complexes. Increasing rigid-body diversity of the models using RosettaDock’s local dock routine lowers the discrimination accuracy with the cognate antibody–antigen pair ranking in bound and unbound models but recovers additional top-ranked cognate complexes when using homology models. The study is the first structure-based cross-docking attempt aimed at distinguishing antibody–antigen binders from non-binders and demonstrates the challenges to address for the methods to be widely applicable to supplement high-throughput experimental antibody sequencing workflows.

T-dependent B cell responses to Plasmodium induce antibodies that form a high-avidity multivalent complex with the circumsporozoite protein

The repeat region of the Plasmodium falciparum circumsporozoite protein (CSP) is a major vaccine antigen because it can be targeted by parasite neutralizing antibodies; however, little is known about this interaction. We used isothermal titration calorimetry, X-ray crystallography and mutagenesis-validated modeling to analyze the binding of a murine neutralizing antibody to Plasmodium falciparum CSP. Strikingly, we found that the repeat region of CSP is bound by multiple antibodies. This repeating pattern allows multiple weak interactions of single F_AB domains to accumulate and yield a complex with a dissociation constant in the low nM range. Because the CSP protein can potentially cross-link multiple B cell receptors (BCRs) we hypothesized that the B cell response might be T cell independent. However, while there was a modest response in mice deficient in T cell help, the bulk of the response was T cell dependent. By sequencing the BCRs of CSP-repeat specific B cells in inbred mice we found that these cells underwent somatic hypermutation and affinity maturation indicative of a T-dependent response. Last, we found that the BCR repertoire of responding B cells was limited suggesting that the structural simplicity of the repeat may limit the breadth of the immune response.

Vaccines aim to protect by inducing the immune system to make molecules called antibodies that can recognize molecules on the surface of invading pathogens. In the case of malaria, our most advanced vaccine candidates aim to promote the production of antibodies that recognize the circumsporozoite protein (CSP) molecule on the surface of the invasive parasite stage called the sporozoite. In this report we use X-ray crystallography to determine the structure of CSP-binding antibodies at the atomic level. We use other techniques such as isothermal titration calorimetry and structural modeling to examine how this antibody interacts with the CSP molecule. Strikingly, we found that each CSP molecule could bind 6 antibodies. This finding has implications for the immune response and may explain why high titers of antibody are needed for protection. Moreover, because the structure of the CSP repeat is quite simple we determined that the number of different kinds of antibodies that could bind this molecule are quite small. However a high avidity interaction between those antibodies and CSP can result from a process called affinity maturation that allows the body to learn how to make improved antibodies specific for pathogen molecules. These data show that while it is challenging for the immune system to recognize and neutralize CSP, it should be possible to generate viable vaccines targeting this molecule.

Modeling oblong proteins and water-mediated interfaces with RosettaDock in CAPRI rounds 28-35

The 28th-35th rounds of the Critical Assessment of PRotein Interactions (CAPRI) served as a practical benchmark for our RosettaDock protein-protein docking protocols, highlighting strengths and weaknesses of the approach. We achieved acceptable or better quality models in three out of 11 targets. For the two α-repeat protein-green fluorescent protein (αrep-GFP) complexes, we used a novel ellipsoidal partial-global docking method (Ellipsoidal Dock) to generate models with 2.2 Å/1.5 Å interface RMSD, capturing 49%/42% of the native contacts, for the 7-/5-repeat αrep complexes. For the DNase-immunity protein complex, we used a new predictor of hydrogen-bonding networks, HBNet with Bridging Waters, to place individual water models at the complex interface; models were generated with 1.8 Å interface RMSD and 12% native water contacts recovered. The targets for which RosettaDock failed to create an acceptable model were typically difficult in general, as six had no acceptable models submitted by any CAPRI predictor. The UCH-L5-RPN13 and UCH-L5-INO80G de-ubiquitinating enzyme-inhibitor complexes comprised inhibitors undergoing significant structural changes upon binding, with the partners being highly interwoven in the docked complexes. Our failure to predict the nucleosome-enzyme complex in Target 95 was largely due to tight constraints we placed on our model based on sparse biochemical data suggesting two specific cross-interface interactions, preventing the correct structure from being sampled. While RosettaDock's three successes show that it is a state-of-the-art docking method, the difficulties with highly flexible and multi-domain complexes highlight the need for better flexible docking and domain-assembly methods. Proteins 2017; 85:479-486. © 2016 Wiley Periodicals, Inc.

Computational design of an epitope-specific Keap1 binding antibody using hotspot residues grafting and CDR loop swapping

Therapeutic and diagnostic applications of monoclonal antibodies often require careful selection of binders that recognize specific epitopes on the target molecule to exert a desired modulation of biological function. Here we present a proof-of-concept application for the rational design of an epitope-specific antibody binding with the target protein Keap1, by grafting pre-defined structural interaction patterns from the native binding partner protein, Nrf2, onto geometrically matched positions of a set of antibody scaffolds. The designed antibodies bind to Keap1 and block the Keap1-Nrf2 interaction in an epitope-specific way. One resulting antibody is further optimised to achieve low-nanomolar binding affinity by in silico redesign of the CDRH3 sequences. An X-ray co-crystal structure of one resulting design reveals that the actual binding orientation and interface with Keap1 is very close to the design model, despite an unexpected CDRH3 tilt and V_H/V_L interface deviation, which indicates that the modelling precision may be improved by taking into account simultaneous CDR loops conformation and V_H/V_L orientation optimisation upon antibody sequence change. Our study confirms that, given a pre-existing crystal structure of the target protein-protein interaction, hotspots grafting with CDR loop swapping is an attractive route to the rational design of an antibody targeting a pre-selected epitope.

Modeling and docking of antibody structures with Rosetta

We describe Rosetta-based computational protocols for predicting the 3D structure of an antibody from sequence (RosettaAntibody) and then docking the antibody to protein antigens (SnugDock). Antibody modeling leverages canonical loop conformations to graft large segments from experimentally determined structures, as well as offering (i) energetic calculations to minimize loops, (ii) docking methodology to refine the V_L-V_H relative orientation and (iii) de novo prediction of the elusive complementarity determining region (CDR) H3 loop. To alleviate model uncertainty, antibody-antigen docking resamples CDR loop conformations and can use multiple models to represent an ensemble of conformations for the antibody, the antigen or both. These protocols can be run fully automated via the ROSIE web server (http://rosie.rosettacommons.org/) or manually on a computer with user control of individual steps. For best results, the protocol requires roughly 1,000 CPU-hours for antibody modeling and 250 CPU-hours for antibody-antigen docking. Tasks can be completed in under a day by using public supercomputers.

A benchmark testing ground for integrating homology modeling and protein docking

Protein docking procedures carry out the task of predicting the structure of a protein-protein complex starting from the known structures of the individual protein components. More often than not, however, the structure of one or both components is not known, but can be derived by homology modeling on the basis of known structures of related proteins deposited in the Protein Data Bank (PDB). Thus, the problem is to develop methods that optimally integrate homology modeling and docking with the goal of predicting the structure of a complex directly from the amino acid sequences of its component proteins. One possibility is to use the best available homology modeling and docking methods. However, the models built for the individual subunits often differ to a significant degree from the bound conformation in the complex, often much more so than the differences observed between free and bound structures of the same protein, and therefore additional conformational adjustments, both at the backbone and side chain levels need to be modeled to achieve an accurate docking prediction. In particular, even homology models of overall good accuracy frequently include localized errors that unfavorably impact docking results. The predicted reliability of the different regions in the model can also serve as a useful input for the docking calculations. Here we present a benchmark dataset that should help to explore and solve combined modeling and docking problems. This dataset comprises a subset of the experimentally solved 'target' complexes from the widely used Docking Benchmark from the Weng Lab (excluding antibody-antigen complexes). This subset is extended to include the structures from the PDB related to those of the individual components of each complex, and hence represent potential templates for investigating and benchmarking integrated homology modeling and docking approaches. Template sets can be dynamically customized by specifying ranges in sequence similarity and in PDB release dates, or using other filtering options, such as excluding sets of specific structures from the template list. Multiple sequence alignments, as well as structural alignments of the templates to their corresponding subunits in the target are also provided. The resource is accessible online or can be downloaded at http://cluspro.org/benchmark, and is updated on a weekly basis in synchrony with new PDB releases. Proteins 2016; 85:10-16. © 2016 Wiley Periodicals, Inc.

Computational Tools for Aiding Rational Antibody Design

Antibodies are a group of proteins responsible for mediating immune reactions in vertebrates. They are able to bind a variety of structural motifs on noxious molecules tagging them for elimination from the organism. As a result of their versatile binding properties, antibodies are currently one of the most important classes of biopharmaceuticals. In this chapter, we discuss how knowledge-based computational methods can aid experimentalists in the development of potent antibodies. When using common experimental methods for antibody development, we often know the sequence of an antibody that binds to our molecule, antigen, of interest. We may also have a structure or model of the antigen. In these cases, computational methods can help by both modeling the antibody and identifying the antibody-antigen contact residues. This information can then play a key role in the rational design of more potent antibodies.

Accurate Structure Prediction of CDR H3 Loops Enabled by a Novel Structure-Based C-Terminal Constraint

Ab structure prediction has made great strides, but accurately modeling CDR H3 loops remains elusive. Unlike the other five CDR loops, CDR H3 does not adopt canonical conformations and must be modeled de novo. During Antibody Modeling Assessment II, we found that biasing simulations toward kinked conformations enables generating low-root mean square deviation models (Weitzner et al. 2014. Proteins 82: 1611-1623), and since then, we have presented new geometric parameters defining the kink conformation (Weitzner et al. 2015. Structure 23: 302-311). In this study, we use these parameters to develop a new biasing constraint. When applied to a benchmark set of high-quality CDR H3 loops, the average minimum root mean square deviation sampled is 0.93 Å, compared with 1.34 Å without the constraint. We then test the performance of the constrained de novo method for homology modeling and rigid-body docking and present the results for 1) the Antibody Modeling Assessment II targets, 2) the 2009 RosettaAntibody benchmark set, and 3) the high-quality set.

A Framework to Simplify Combined Sampling Strategies in Rosetta

A core task in computational structural biology is the search of conformational space for low energy configurations of a biological macromolecule. Because conformational space has a very high dimensionality, the most successful search methods integrate some form of prior knowledge into a general sampling algorithm to reduce the effective dimensionality. However, integrating multiple types of constraints can be challenging. To streamline the incorporation of diverse constraints, we developed the Broker: an extension of the Rosetta macromolecular modeling suite that can express a wide range of protocols using constraints by combining small, independent modules, each of which implements a different set of constraints. We demonstrate expressiveness of the Broker through several code vignettes. The framework enables rapid protocol development in both biomolecular design and structural modeling tasks and thus is an important step towards exposing the rich functionality of Rosetta's core libraries to a growing community of users addressing a diverse set of tasks in computational biology.

High-resolution protein-protein docking by global optimization: recent advances and future challenges

A computational protein-protein docking method that predicts atomic details of protein-protein interactions from protein monomer structures is an invaluable tool for understanding the molecular mechanisms of protein interactions and for designing molecules that control such interactions. Compared to low-resolution docking, high-resolution docking explores the conformational space in atomic resolution to provide predictions with atomic details. This allows for applications to more challenging docking problems that involve conformational changes induced by binding. Recently, high-resolution methods have become more promising as additional information such as global shapes or residue contacts are now available from experiments or sequence/structure data. In this review article, we highlight developments in high-resolution docking made during the last decade, specifically regarding global optimization methods employed by the docking methods. We also discuss two major challenges in high-resolution docking: prediction of backbone flexibility and water-mediated interactions.

Antibody specific epitope prediction-emergence of a new paradigm

The development of accurate tools for predicting B-cell epitopes is important but difficult. Traditional methods have examined which regions in an antigen are likely binding sites of an antibody. However, it is becoming increasingly clear that most antigen surface residues will be able to bind one or more of the myriad of possible antibodies. In recent years, new approaches have emerged for predicting an epitope for a specific antibody, utilizing information encoded in antibody sequence or structure. Applying such antibody-specific predictions to groups of antibodies in combination with easily obtainable experimental data improves the performance of epitope predictions. We expect that further advances of such tools will be possible with the integration of immunoglobulin repertoire sequencing data.

Peptides designed to spatially depict the Epstein-Barr virus major virion glycoprotein gp350 neutralization epitope elicit antibodies that block virus-neutralizing antibody 72A1 interaction with the native gp350 molecule

Epstein-Barr virus (EBV) is the etiologic agent of infectious mononucleosis and the root cause of B-cell lymphoproliferative disease in individuals with a weakened immune system, as well as a principal cofactor in nasopharyngeal carcinoma, various lymphomas, and other cancers. The EBV major virion surface glycoprotein gp350 is viewed as the best vaccine candidate to prevent infectious mononucleosis in healthy EBV-naive persons and EBV-related cancers in at-risk individuals. Previous epitope mapping of gp350 revealed only one dominant neutralizing epitope, which has been shown to be the target of the monoclonal antibody 72A1. Computer modeling of the 72A1 antibody interaction with the gp350 amino terminus was used to identify gp350 amino acids that could form strong ionic, electrostatic, or hydrogen bonds with the 72A1 antibody. Peptide DDRTTLQLAQNPVYIPETYPYIKWDN (designated peptide 2) and peptide GSAKPGNGSYFASVKTEMLGNEID (designated peptide 3) were designed to spatially represent the gp350 amino acids predicted to interact with the 72A1 antibody paratope. Peptide 2 bound to the 72A1 antibody and blocked 72A1 antibody recognition of the native gp350 molecule. Peptide 2 and peptide 3 were recognized by human IgG and shown to elicit murine antibodies that could target gp350 and block its recognition by the 72A1 antibody. This work provides a structural mapping of the interaction between the EBV-neutralizing antibody 72A1 and the major virion surface protein gp350. gp350 mimetic peptides that spatially depict the EBV-neutralizing epitope would be useful as a vaccine to focus the immune system exclusively to this important virus epitope.

IMPORTANCE

The production of virus-neutralizing antibodies targeting the Epstein-Barr virus (EBV) major surface glycoprotein gp350 is important for the prevention of infectious mononucleosis and EBV-related cancers. The data presented here provide the first in silico map of the gp350 interaction with a virus-blocking monoclonal antibody. Immunization with gp350 peptides identified by in silico mapping generated antibodies that cross-react with the EBV gp350 molecule and block recognition of the gp350 molecule by a virus-neutralizing antibody. Through its ability to focus the immune system exclusively on the gp350 sequence important for viral entry, these peptides may form the basis of an EBV vaccine candidate. This strategy would sidestep the production of other irrelevant gp350 antibodies that divert the immune system from generating a protective antiviral response or that impede access to the virus-blocking epitope by protective antibodies.

Structural and electrostatic analysis of HLA B-cell epitopes: inference on immunogenicity and prediction of humoral alloresponses

PURPOSE OF REVIEW

The immunogenic capacity of donor human leukocyte antigen (HLA) to induce humoral immune responses is not an intrinsic property of the mismatched alloantigen but depends on the HLA phenotype of the recipient. In recent years, advances in molecular sequence technology and information from X-ray crystallography have enabled structural comparison of donor and recipient HLA type providing an opportunity for a more rational approach for determining HLA compatibility. In this article, we review studies investigating the molecular basis of antibody-antigen interactions and present computational approaches to determine the complex physiochemical and structural properties of B-cell epitopes.

RECENT FINDINGS

The relative immunogenicity of individual HLA mismatches may be predicted from analysis of polymorphic amino acids at continuous and discontinuous HLA sequence positions. The use of alloantigen sequence information alone, however, provides limited insight into key determinants of B-cell epitope immunogenicity, such as the orientation, accessibility and physiochemical properties of amino acid side chains. Advances in computational molecular modelling techniques now enable assessment of HLA-alloantibody interactions at the atomic level. Recent evidence supports a strong link between HLA B-cell epitope surface electrostatic potential and their immunogenicity.

SUMMARY

Assessment of the surface electrostatic properties of HLA alloantigens and computational analyses of HLA-alloantibody interactions represent a promising area for future research into the molecular basis of HLA immunogenicity and antigenicity.

Prediction of Antibody Epitopes

Antibodies recognize their cognate antigens in a precise and effective way. In order to do so, they target regions of the antigenic molecules that have specific features such as large exposed areas, presence of charged or polar atoms, specific secondary structure elements, and lack of similarity to self-proteins. Given the sequence or the structure of a protein of interest, several methods exploit such features to predict the residues that are more likely to be recognized by an immunoglobulin. Here, we present two methods (BepiPred and DiscoTope) to predict linear and discontinuous antibody epitopes from the sequence and/or the three-dimensional structure of a target protein.

Antibodies: Computer-Aided Prediction of Structure and Design of Function

With the advent of high-throughput sequencing, and the increased availability of experimental structures of antibodies and antibody-antigen complexes, comes the improvement of computational approaches to predict the structure and design the function of antibodies and antibody-antigen complexes. While antibodies pose formidable challenges for protein structure prediction and design due to their large size and highly flexible loops in the complementarity-determining regions, they also offer exciting opportunities: the central importance of antibodies for human health results in a wealth of structural and sequence information that—as a knowledge base—can drive the modeling algorithms by limiting the conformational and sequence search space to likely regions of success. Further, efficient experimental platforms exist to test predicted antibody structure or designed antibody function, thereby leading to an iterative feedback loop between computation and experiment. We briefly review the history of computer-aided prediction of structure and design of function in the antibody field before we focus on recent methodological developments and the most exciting application examples.

A flexible docking approach for prediction of T cell receptor-peptide-MHC complexes

T cell receptors (TCRs) are immune proteins that specifically bind to antigenic molecules, which are often foreign peptides presented by major histocompatibility complex proteins (pMHCs), playing a key role in the cellular immune response. To advance our understanding and modeling of this dynamic immunological event, we assembled a protein-protein docking benchmark consisting of 20 structures of crystallized TCR/pMHC complexes for which unbound structures exist for both TCR and pMHC. We used our benchmark to compare predictive performance using several flexible and rigid backbone TCR/pMHC docking protocols. Our flexible TCR docking algorithm, TCRFlexDock, improved predictive success over the fixed backbone protocol, leading to near-native predictions for 80% of the TCR/pMHC cases among the top 10 models, and 100% of the cases in the top 30 models. We then applied TCRFlexDock to predict the two distinct docking modes recently described for a single TCR bound to two different antigens, and tested several protein modeling scoring functions for prediction of TCR/pMHC binding affinities. This algorithm and benchmark should enable future efforts to predict, and design of uncharacterized TCR/pMHC complexes.

Blind prediction performance of RosettaAntibody 3.0: grafting, relaxation, kinematic loop modeling, and full CDR optimization

Antibody Modeling Assessment II (AMA-II) provided an opportunity to benchmark RosettaAntibody on a set of 11 unpublished antibody structures. RosettaAntibody produced accurate, physically realistic models, with all framework regions and 42 of the 55 non-H3 CDR loops predicted to under an Ångström. The performance is notable when modeling H3 on a homology framework, where RosettaAntibody produced the best model among all participants for four of the 11 targets, two of which were predicted with sub-Ångström accuracy. To improve RosettaAntibody, we pursued the causes of model errors. The most common limitation was template unavailability, underscoring the need for more antibody structures and/or better de novo loop methods. In some cases, better templates could have been found by considering residues outside of the CDRs. De novo CDR H3 modeling remains challenging at long loop lengths, but constraining the C-terminal end of H3 to a kinked conformation allows near-native conformations to be sampled more frequently. We also found that incorrect VL -VH orientations caused models with low H3 RMSDs to score poorly, suggesting that correct VL -VH orientations will improve discrimination between near-native and incorrect conformations. These observations will guide the future development of RosettaAntibody.

Increasing stability of antibody via antibody engineering: stability engineering on an anti-hVEGF

Antibody stability is very important for expression, activity, specificity, and storage. This knowledge of antibody structure has made it possible for a computer-aided molecule design to be used to optimize and increase antibody stability. Many computational methods have been built based on knowledge or structure, however, a good integrated engineering system has yet to be developed that combines these methods. In the current study, we designed an integrated computer-aided engineering protocol, which included several successful methods. Mutants were designed considering factors that affected stability and multiwall filter screening was used to improve the design accuracy. Using this protocol, the thermo-stability of an anti-hVEGF antibody was significantly improved. Nearly 40% of the single-point mutants proved to be more stable than the parent antibody and most of the mutations could be stacked effectively. The T₅₀ also improved about 7°C by combinational mutation of seven sites in the light chain and three sites in the heavy chain. Data indicate that the protocol is an effective method for optimization of antibody structure, especially for improving thermo-stability. This protocol could also be used to enhance the stability of other antibodies.

Automated antibody structure prediction using Accelrys tools: results and best practices

We describe the methodology and results from our participation in the second Antibody Modeling Assessment experiment. During the experiment we predicted the structure of eleven unpublished antibody Fv fragments. Our prediction methods centered on template-based modeling; potential templates were selected from an antibody database based on their sequence similarity to the target in the framework regions. Depending on the quality of the templates, we constructed models of the antibody framework regions either using a single, chimeric or multiple template approach. The hypervariable loop regions in the initial models were rebuilt by grafting the corresponding regions from suitable templates onto the model. For the H3 loop region, we further refined models using ab initio methods. The final models were subjected to constrained energy minimization to resolve severe local structural problems. The analysis of the models submitted show that Accelrys tools allow for the construction of quite accurate models for the framework and the canonical CDR regions, with RMSDs to the X-ray structure on average below 1 Å for most of these regions. The results show that accurate prediction of the H3 hypervariable loops remains a challenge. Furthermore, model quality assessment of the submitted models show that the models are of quite high quality, with local geometry assessment scores similar to that of the target X-ray structures. Proteins 2014; 82:1583–1598. © 2014 The Authors. Proteins published by Wiley Periodicals, Inc.

Improving B-cell epitope prediction and its application to global antibody-antigen docking

Motivation: Antibodies are currently the most important class of biopharmaceuticals. Development of such antibody-based drugs depends on costly and time-consuming screening campaigns. Computational techniques such as antibody–antigen docking hold the potential to facilitate the screening process by rapidly providing a list of initial poses that approximate the native complex.

Results: We have developed a new method to identify the epitope region on the antigen, given the structures of the antibody and the antigen—EpiPred. The method combines conformational matching of the antibody–antigen structures and a specific antibody–antigen score. We have tested the method on both a large non-redundant set of antibody–antigen complexes and on homology models of the antibodies and/or the unbound antigen structure. On a non-redundant test set, our epitope prediction method achieves 44% recall at 14% precision against 23% recall at 14% precision for a background random distribution. We use our epitope predictions to rescore the global docking results of two rigid-body docking algorithms: ZDOCK and ClusPro. In both cases including our epitope, prediction increases the number of near-native poses found among the top decoys.

Availability and implementation: Our software is available from http://www.stats.ox.ac.uk/research/proteins/resources.

Contact:deane@stats.ox.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Antibody engineering and therapeutics, The Annual Meeting of the Antibody Society: December 8-12, 2013, Huntington Beach, CA

The 24^th Antibody Engineering & Therapeutics meeting brought together a broad range of participants who were updated on the latest advances in antibody research and development. Organized by IBC Life Sciences, the gathering is the annual meeting of The Antibody Society, which serves as the scientific sponsor. Preconference workshops on 3D modeling and delineation of clonal lineages were featured, and the conference included sessions on a wide variety of topics relevant to researchers, including systems biology; antibody deep sequencing and repertoires; the effects of antibody gene variation and usage on antibody response; directed evolution; knowledge-based design; antibodies in a complex environment; polyreactive antibodies and polyspecificity; the interface between antibody therapy and cellular immunity in cancer; antibodies in cardiometabolic medicine; antibody pharmacokinetics, distribution and off-target toxicity; optimizing antibody formats for immunotherapy; polyclonals, oligoclonals and bispecifics; antibody discovery platforms; and antibody-drug conjugates.

Extending RosettaDock with water, sugar, and pH for prediction of complex structures and affinities for CAPRI rounds 20-27

Rounds 20-27 of the Critical Assessment of PRotein Interactions (CAPRI) provided a testing platform for computational methods designed to address a wide range of challenges. The diverse targets drove the creation of and new combinations of computational tools. In this study, RosettaDock and other novel Rosetta protocols were used to successfully predict four of the 10 blind targets. For example, for DNase domain of Colicin E2-Im2 immunity protein, RosettaDock and RosettaLigand were used to predict the positions of water molecules at the interface, recovering 46% of the native water-mediated contacts. For α-repeat Rep4-Rep2 and g-type lysozyme-PliG inhibitor complexes, homology models were built and standard and pH-sensitive docking algorithms were used to generate structures with interface RMSD values of 3.3 Å and 2.0 Å, respectively. A novel flexible sugar-protein docking protocol was also developed and used for structure prediction of the BT4661-heparin-like saccharide complex, recovering 71% of the native contacts. Challenges remain in the generation of accurate homology models for protein mutants and sampling during global docking. On proteins designed to bind influenza hemagglutinin, only about half of the mutations were identified that affect binding (T55: 54%; T56: 48%). The prediction of the structure of the xylanase complex involving homology modeling and multidomain docking pushed the limits of global conformational sampling and did not result in any successful prediction. The diversity of problems at hand requires computational algorithms to be versatile; the recent additions to the Rosetta suite expand the capabilities to encompass more biologically realistic docking problems.

Antibody i-Patch prediction of the antibody binding site improves rigid local antibody-antigen docking

Antibodies are a class of proteins indispensable for the vertebrate immune system. The general architecture of all antibodies is very similar, but they contain a hypervariable region which allows millions of antibody variants to exist, each of which can bind to different molecules. This binding malleability means that antibodies are an increasingly important category of biopharmaceuticals and biomarkers. We present Antibody i-Patch, a method that annotates the most likely antibody residues to be in contact with the antigen. We show that our predictions correlate with energetic importance and thus we argue that they may be useful in guiding mutations in the artificial affinity maturation process. Using our predictions as constraints for a rigid-body docking algorithm, we are able to obtain high-quality results in minutes. Our annotation method and re-scoring system for docking achieve their predictive power by using antibody-specific statistics. Antibody i-Patch is available from http://www.stats.ox.ac.uk/research/proteins/resources.

Application of asymmetric statistical potentials to antibody-protein docking

MOTIVATION

An effective docking algorithm for antibody-protein antigen complex prediction is an important first step toward design of biologics and vaccines. We have recently developed a new class of knowledge-based interaction potentials called Decoys as the Reference State (DARS) and incorporated DARS into the docking program PIPER based on the fast Fourier transform correlation approach. Although PIPER was the best performer in the latest rounds of the CAPRI protein docking experiment, it is much less accurate for docking antibody-protein antigen pairs than other types of complexes, in spite of incorporating sequence-based information on the location of the paratope. Analysis of antibody-protein antigen complexes has revealed an inherent asymmetry within these interfaces. Specifically, phenylalanine, tryptophan and tyrosine residues highly populate the paratope of the antibody but not the epitope of the antigen.

RESULTS

Since this asymmetry cannot be adequately modeled using a symmetric pairwise potential, we have removed the usual assumption of symmetry. Interaction statistics were extracted from antibody-protein complexes under the assumption that a particular atom on the antibody is different from the same atom on the antigen protein. The use of the new potential significantly improves the performance of docking for antibody-protein antigen complexes, even without any sequence information on the location of the paratope. We note that the asymmetric potential captures the effects of the multi-body interactions inherent to the complex environment in the antibody-protein antigen interface.

AVAILABILITY

The method is implemented in the ClusPro protein docking server, available at http://cluspro.bu.edu.

Serverification of molecular modeling applications: the Rosetta Online Server that Includes Everyone (ROSIE)

The Rosetta molecular modeling software package provides experimentally tested and rapidly evolving tools for the 3D structure prediction and high-resolution design of proteins, nucleic acids, and a growing number of non-natural polymers. Despite its free availability to academic users and improving documentation, use of Rosetta has largely remained confined to developers and their immediate collaborators due to the code's difficulty of use, the requirement for large computational resources, and the unavailability of servers for most of the Rosetta applications. Here, we present a unified web framework for Rosetta applications called ROSIE (Rosetta Online Server that Includes Everyone). ROSIE provides (a) a common user interface for Rosetta protocols, (b) a stable application programming interface for developers to add additional protocols, (c) a flexible back-end to allow leveraging of computer cluster resources shared by RosettaCommons member institutions, and (d) centralized administration by the RosettaCommons to ensure continuous maintenance. This paper describes the ROSIE server infrastructure, a step-by-step 'serverification' protocol for use by Rosetta developers, and the deployment of the first nine ROSIE applications by six separate developer teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance, Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated by the number and diversity of these applications, ROSIE offers a general and speedy paradigm for serverification of Rosetta applications that incurs negligible cost to developers and lowers barriers to Rosetta use for the broader biological community. ROSIE is available at http://rosie.rosettacommons.org.

Structure-based non-canonical amino acid design to covalently crosslink an antibody-antigen complex

Engineering antibodies to utilize non-canonical amino acids (NCAA) should greatly expand the utility of an already important biological reagent. In particular, introducing crosslinking reagents into antibody complementarity determining regions (CDRs) should provide a means to covalently crosslink residues at the antibody-antigen interface. Unfortunately, finding the optimum position for crosslinking two proteins is often a matter of iterative guessing, even when the interface is known in atomic detail. Computer-aided antibody design can potentially greatly restrict the number of variants that must be explored in order to identify successful crosslinking sites. We have therefore used Rosetta to guide the introduction of an oxidizable crosslinking NCAA, l-3,4-dihydroxyphenylalanine (l-DOPA), into the CDRs of the anti-protective antigen scFv antibody M18, and have measured crosslinking to its cognate antigen, domain 4 of the anthrax protective antigen. Computed crosslinking distance, solvent accessibility, and interface energetics were three factors considered that could impact the efficiency of l-DOPA-mediated crosslinking. In the end, 10 variants were synthesized, and crosslinking efficiencies were generally 10% or higher, with the best variant crosslinking to 52% of the available antigen. The results suggest that computational analysis can be used in a pipeline for engineering crosslinking antibodies. The rules learned from l-DOPA crosslinking of antibodies may also be generalizable to the formation of other crosslinked interfaces and complexes.