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.

A TRAV26-1-encoded recognition motif focuses the biased T cell response in celiac disease.

The semi-public T-cell response towards the gluten epitope DQ2.5-glia-α2 uses a prototypic TCR encoded by the germline segments TRAV26-1 and TRBV7-2. Through mutagenesis experiments, we show that a TRAV26-1encoded recognition motif contacts the MHC ß-chain and the TCR CDR3ß loop underpinning this conserved T-cell response restricted to the prototypic TCRs.

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.

Consolidation of proteomics data in the Cancer Proteomics database.

Cancer is a class of diseases characterized by abnormal cell growth and one of the major reasons for human deaths. Proteins are involved in the molecular mechanisms leading to cancer, furthermore they are affected by anti-cancer drugs, and protein biomarkers can be used to diagnose certain cancer types. Therefore, it is important to explore the proteomics background of cancer. In this report, we developed the Cancer Proteomics database to re-interrogate published proteome studies investigating cancer. The database is divided in three sections related to cancer processes, cancer types, and anti-cancer drugs. Currently, the Cancer Proteomics database contains 9778 entries of 4118 proteins extracted from 143 scientific articles covering all three sections: cell death (cancer process), prostate cancer (cancer type) and platinum-based anti-cancer drugs including carboplatin, cisplatin, and oxaliplatin (anti-cancer drugs). The detailed information extracted from the literature includes basic information about the articles (e.g., PubMed ID, authors, journal name, publication year), information about the samples (type, study/reference, prognosis factor), and the proteomics workflow (Subcellular fractionation, protein, and peptide separation, mass spectrometry, quantification). Useful annotations such as hyperlinks to UniProt and PubMed were included. In addition, many filtering options were established as well as export functions. The database is freely available at http://cancerproteomics.uio.no.

Computational design of nanomolar-binding antibodies specific to multiple SARS-CoV-2 variants by engineering a specificity switch of antibody 80R using RosettaAntibodyDesign (RAbD) results in potential generalizable therapeutic antibodies for novel SARS-CoV-2 virus.

The human infectious disease COVID-19 caused by the SARS-CoV-2 virus has become a major threat to global public health. Developing a vaccine is the preferred prophylactic response to epidemics and pandemics. However, for individuals who have contracted the disease, the rapid design of antibodies that can target the SARS-CoV-2 virus fulfils a critical need. Further, discovering antibodies that bind multiple variants of SARS-CoV-2 can aid in the development of rapid antigen tests (RATs) which are critical for the identification and isolation of individuals currently carrying COVID-19. Here we provide a proof-of-concept study for the computational design of high-affinity antibodies that bind to multiple variants of the SARS-CoV-2 spike protein using RosettaAntibodyDesign (RAbD). Well characterized antibodies that bind with high affinity to the SARS-CoV-1 (but not SARS-CoV-2) spike protein were used as templates and re-designed to bind the SARS-CoV-2 spike protein with high affinity, resulting in a specificity switch. A panel of designed antibodies were experimentally validated. One design bound to a broad range of variants of concern including the Omicron, Delta, Wuhan, and South African spike protein variants.

Corrigendum to "Computational design of nanomolar-binding antibodies specific to multiple SARS-CoV-2 variants by engineering a specificity switch of antibody 80R using RosettaAntibodyDesign (RAbD) results in potential generalizable therapeutic antibodies for novel SARS-CoV-2 virus" [Heliyon 9(4) (April 2023) e15032].

[This corrects the article DOI: 10.1016/j.heliyon.2023.e15032.].

Hallucinating structure-conditioned antibody libraries for target-specific binders.

Antibodies are widely developed and used as therapeutics to treat cancer, infectious disease, and inflammation. During development, initial leads routinely undergo additional engineering to increase their target affinity. Experimental methods for affinity maturation are expensive, laborious, and time-consuming and rarely allow the efficient exploration of the relevant design space. Deep learning (DL) models are transforming the field of protein engineering and design. While several DL-based protein design methods have shown promise, the antibody design problem is distinct, and specialized models for antibody design are desirable. Inspired by hallucination frameworks that leverage accurate structure prediction DL models, we propose the FvHallucinator for designing antibody sequences, especially the CDR loops, conditioned on an antibody structure. Such a strategy generates targeted CDR libraries that retain the conformation of the binder and thereby the mode of binding to the epitope on the antigen. On a benchmark set of 60 antibodies, FvHallucinator generates sequences resembling natural CDRs and recapitulates perplexity of canonical CDR clusters. Furthermore, the FvHallucinator designs amino acid substitutions at the VH-VL interface that are enriched in human antibody repertoires and therapeutic antibodies. We propose a pipeline that screens FvHallucinator designs to obtain a library enriched in binders for an antigen of interest. We apply this pipeline to the CDR H3 of the Trastuzumab-HER2 complex to generate in silico designs predicted to improve upon the binding affinity and interfacial properties of the original antibody. Thus, the FvHallucinator pipeline enables generation of inexpensive, diverse, and targeted antibody libraries enriched in binders for antibody affinity maturation.

Affinity maturation of TCR-like antibodies using phage display guided by structural modeling.

TCR-like antibodies represent a unique type of engineered antibodies with specificity toward pHLA, a ligand normally restricted to the sensitive recognition by T cells. Here, we report a phage display-based sequential development path of such antibodies. The strategy goes from initial lead identification through in silico informed CDR engineering in combination with framework engineering for affinity and thermostability optimization, respectively. The strategy allowed the identification of HLA-DQ2.5 gluten peptide-specific TCR-like antibodies with low picomolar affinity. Our method outlines an efficient and general method for development of this promising class of antibodies, which should facilitate their utility including translation to human therapy.

Biophysical differences in IgG1 Fc-based therapeutics relate to their cellular handling, interaction with FcRn and plasma half-life.

Antibody-based therapeutics (ABTs) are used to treat a range of diseases. Most ABTs are either full-length IgG1 antibodies or fusions between for instance antigen (Ag)-binding receptor domains and the IgG1 Fc fragment. Interestingly, their plasma half-life varies considerably, which may relate to how they engage the neonatal Fc receptor (FcRn). As such, there is a need for an in-depth understanding of how different features of ABTs affect FcRn-binding and transport behavior. Here, we report on how FcRn-engagement of the IgG1 Fc fragment compare to clinically relevant IgGs and receptor domain Fc fusions, binding to VEGF or TNF-α. The results reveal FcRn-dependent intracellular accumulation of the Fc, which is in line with shorter plasma half-life than that of full-length IgG1 in human FcRn-expressing mice. Receptor domain fusion to the Fc increases its half-life, but not to the extent of IgG1. This is mirrored by a reduced cellular recycling capacity of the Fc-fusions. In addition, binding of cognate Ag to ABTs show that complexes of similar size undergo cellular transport at different rates, which could be explained by the biophysical properties of each ABT. Thus, the study provides knowledge that should guide tailoring of ABTs regarding optimal cellular sorting and plasma half-life.

Antibody variable sequences have a pronounced effect on cellular transport and plasma half-life.

Monoclonal IgG antibodies are the fastest growing class of biologics, but large differences exist in their plasma half-life in humans. Thus, to design IgG antibodies with favorable pharmacokinetics, it is crucial to identify the determinants of such differences. Here, we demonstrate that the variable region sequences of IgG antibodies greatly affect cellular uptake and subsequent recycling and rescue from intracellular degradation by endothelial cells. When the variable sequences are masked by the cognate antigen, it influences both their transport behavior and binding to the neonatal Fc receptor (FcRn), a key regulator of IgG plasma half-life. Furthermore, we show how charge patch differences in the variable domains modulate both binding and transport properties and that a short plasma half-life, due to unfavorable charge patches, may partly be overcome by Fc-engineering for improved FcRn binding.

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.

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.

Targeting the MHC Ligandome by Use of TCR-Like Antibodies.

Monoclonal antibodies (mAbs) are valuable as research reagents, in diagnosis and in therapy. Their high specificity, the ease in production, favorable biophysical properties and the opportunity to engineer different properties make mAbs a versatile class of biologics. mAbs targeting peptide-major histocompatibility molecule (pMHC) complexes are often referred to as "TCR-like" mAbs, as pMHC complexes are generally recognized by T-cell receptors (TCRs). Presentation of self- and non-self-derived peptide fragments on MHC molecules and subsequent activation of T cells dictate immune responses in health and disease. This includes responses to infectious agents or cancer but also aberrant responses against harmless self-peptides in autoimmune diseases. The ability of TCR-like mAbs to target specific peptides presented on MHC allows for their use to study peptide presentation or for diagnosis and therapy. This extends the scope of conventional mAbs, which are generally limited to cell-surface or soluble antigens. Herein, we review the strategies used to generate TCR-like mAbs and provide a structural comparison with the analogous TCR in pMHC binding. We further discuss their applications as research tools and therapeutic reagents in preclinical models as well as challenges and limitations associated with their use.

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 VL-VH 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 TCRα framework-centered codon shapes a biased T cell repertoire through direct MHC and CDR3ß interactions.

Selection of biased T cell receptor (TCR) repertoires across individuals is seen in both infectious diseases and autoimmunity, but the underlying molecular basis leading to these shared repertoires remains unclear. Celiac disease (CD) occurs primarily in HLA-DQ2.5+ individuals and is characterized by a CD4+ T cell response against gluten epitopes dominated by DQ2.5-glia-α1a and DQ2.5-glia-α2. The DQ2.5-glia-α2 response recruits a highly biased TCR repertoire composed of TRAV26-1 paired with TRBV7-2 harboring a semipublic CDR3ß loop. We aimed to unravel the molecular basis for this signature. By variable gene segment exchange, directed mutagenesis, and cellular T cell activation studies, we found that TRBV7-3 can substitute for TRBV7-2, as both can contain the canonical CDR3ß loop. Furthermore, we identified a pivotal germline-encoded MHC recognition motif centered on framework residue Y40 in TRAV26-1 engaging both DQB1*02 and the canonical CDR3ß. This allowed prediction of expanded DQ2.5-glia-α2-reactive TCR repertoires, which were confirmed by single-cell sorting and TCR sequencing from CD patient samples. Our data refine our understanding of how HLA-dependent biased TCR repertoires are selected in the periphery due to germline-encoded residues.