MoFvAb: Modeling the Fv region of 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.

Structures and immune recognition of Env trimers from two Asia prevalent HIV-1 CRFs

Structure-guided immunofocusing HIV-1 vaccine design entails a comprehensive understanding of Envs from diverse HIV-1 subtypes, including circulating recombinant forms (CRFs). Here, we present the cryo-EM structures of Envs from two Asia prevalent CRFs (CRF01_AE and CRF07_BC) at 3.0 and 3.5 Å. We compare the structures and glycosylation patterns of Envs from different subtypes and perform cross-clade statistical analyses to reveal the unique features of CRF01_AE V1 region, which are associated with the resistance to certain bNAbs. We also solve a 4.1 Å cryo-EM structure of CRF01_AE Env in complex with F6, the first bNAb from CRF01_AE-infected individuals. F6 recognizes a gp120-gp41 spanning epitope to allosterically destabilize the Env trimer apex and weaken inter-protomer packing, which in turn hinders the receptor binding and induces Env trimer disassembly, demonstrating a dual mechanism of neutralization. These findings broaden our understanding of CRF Envs and shed lights on immunofocusing HIV-1 vaccine design.

Broadly neutralizing humanized SARS-CoV-2 antibody binds to a conserved epitope on Spike and provides antiviral protection through inhalation-based delivery in non-human primates

The COVID-19 pandemic represents a global challenge that has impacted and is expected to continue to impact the lives and health of people across the world for the foreseeable future. The rollout of vaccines has provided highly anticipated relief, but effective therapeutics are required to further reduce the risk and severity of infections. Monoclonal antibodies have been shown to be effective as therapeutics for SARS-CoV-2, but as new variants of concern (VoC) continue to emerge, their utility and use have waned due to limited or no efficacy against these variants. Furthermore, cumbersome systemic administration limits easy and broad access to such drugs. As well, concentrations of systemically administered antibodies in the mucosal epithelium, a primary site of initial infection, are dependent on neonatal Fc receptor mediated transport and require high drug concentrations. To reduce the viral load more effectively in the lung, we developed an inhalable formulation of a SARS-CoV-2 neutralizing antibody binding to a conserved epitope on the Spike protein, ensuring pan-neutralizing properties. Administration of this antibody via a vibrating mesh nebulization device retained antibody integrity and resulted in effective distribution of the antibody in the upper and lower respiratory tract of non-human primates (NHP). In comparison with intravenous administration, significantly higher antibody concentrations can be obtained in the lung, resulting in highly effective reduction in viral load post SARS-CoV-2 challenge. This approach may reduce the barriers of access and uptake of antibody therapeutics in real-world clinical settings and provide a more effective blueprint for targeting existing and potentially emerging respiratory tract viruses.

Identifying developability risks for clinical progression of antibodies using high-throughput in vitro and in silico approaches

With the growing significance of antibodies as a therapeutic class, identifying developability risks early during development is of paramount importance. Several high-throughput in vitro assays and in silico approaches have been proposed to de-risk antibodies during early stages of the discovery process. In this review, we have compiled and collectively analyzed published experimental assessments and computational metrics for clinical antibodies. We show that flags assigned based on in vitro measurements of polyspecificity and hydrophobicity are more predictive of clinical progression than their in silico counterparts. Additionally, we assessed the performance of published models for developability predictions on molecules not used during model training. We find that generalization to data outside of those used for training remains a challenge for models. Finally, we highlight the challenges of reproducibility in computed metrics arising from differences in homology modeling, in vitro assessments relying on complex reagents, as well as curation of experimental data often used to assess the utility of high-throughput approaches. We end with a recommendation to enable assay reproducibility by inclusion of controls with disclosed sequences, as well as sharing of structural models to enable the critical assessment and improvement of in silico predictions.

Comparison of hydrophobicity scales for predicting biophysical properties of antibodies

While antibody-based therapeutics have grown to be one of the major classes of novel medicines, some antibody development candidates face significant challenges regarding expression levels, solubility, as well as stability and aggregation, under physiological and storage conditions. A major determinant of those properties is surface hydrophobicity, which promotes unspecific interactions and has repeatedly proven problematic in the development of novel antibody-based drugs. Multiple computational methods have been devised for in-silico prediction of antibody hydrophobicity, often using hydrophobicity scales to assign values to each amino acid. Those approaches are usually validated by their ability to rank potential therapeutic antibodies in terms of their experimental hydrophobicity. However, there is significant diversity both in the hydrophobicity scales and in the experimental methods, and consequently in the performance of in-silico methods to predict experimental results. In this work, we investigate hydrophobicity of monoclonal antibodies using hydrophobicity scales. We implement several scoring schemes based on the solvent-accessibility and the assigned hydrophobicity values, and compare the different scores and scales based on their ability to predict retention times from hydrophobic interaction chromatography. We provide an overview of the strengths and weaknesses of several commonly employed hydrophobicity scales, thereby improving the understanding of hydrophobicity in antibody development. Furthermore, we test several datasets, both publicly available and proprietary, and find that the diversity of the dataset affects the performance of hydrophobicity scores. We expect that this work will provide valuable guidelines for the optimization of biophysical properties in future drug discovery campaigns.

Investigating brain uptake of a non-targeting monoclonal antibody after intravenous and intracerebroventricular administration

Monoclonal antibodies play an important role in the treatment of various diseases. However, the development of these drugs against neurological disorders where the drug target is located in the brain is challenging and requires a good understanding of the local drug concentration in the brain. In this original research, we investigated the systemic and local pharmacokinetics in the brain of healthy rats after either intravenous (IV) or intracerebroventricular (ICV) administration of EGFRvIII-T-Cell bispecific (TCB), a bispecific monoclonal antibody. We established an experimental protocol that allows serial sampling in serum, cerebrospinal fluid (CSF) and interstitial fluid (ISF) of the prefrontal cortex in freely moving rats. For detection of drug concentration in ISF, a push-pull microdialysis technique with large pore membranes was applied. Brain uptake into CSF and ISF was characterized and quantified with a reduced brain physiologically-based pharmacokinetic model. The model allowed us to interpret the pharmacokinetic processes of brain uptake after different routes of administration. The proposed model capturing the pharmacokinetics in serum, CSF and ISF of the prefrontal cortex suggests a barrier function between the CSF and ISF that impedes free antibody transfer. This finding suggests that ICV administration may not be better suited to reach higher local drug exposure as compared to IV administration. The model enabled us to quantify the relative contribution of the blood-brain barrier (BBB) and Blood-CSF-Barrier to the uptake into the interstitial fluid of the brain. In addition, we compared the brain uptake of three monoclonal antibodies after IV dosing. In summary, the presented approach can be applied to profile compounds based on their relative uptake in the brain and provides quantitative insights into which pathways are contributing to the net exposure in the brain.

VHH Structural Modelling Approaches: A Critical Review

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

Inhibitory Monoclonal Antibodies and Their Recombinant Derivatives Targeting Surface-Exposed Carbonic Anhydrase XII on Cancer Cells

Monoclonal and recombinant antibodies are widely used for the diagnostics and therapy of cancer. They are generated to interact with cell surface proteins which are usually involved in the development and progression of cancer. Carbonic anhydrase XII (CA XII) contributes to the survival of tumors under hypoxic conditions thus is considered a candidate target for antibody-based therapy. In this study, we have generated a novel collection of monoclonal antibodies (MAbs) against the recombinant extracellular domain of CA XII produced in HEK-293 cells. Eighteen out of 24 MAbs were reactive with cellular CA XII on the surface of live kidney and lung cancer cells as determined by flow cytometry. One MAb 14D6 also inhibited the enzymatic activity of recombinant CA XII as measured by the stopped-flow assay. MAb 14D6 showed the migrastatic effect on human lung carcinoma A549 and renal carcinoma A498 cell lines in a ‘wound healing’ assay. It did not reduce the growth of multicellular lung and renal cancer spheroids but reduced the cell viability by the ATP Bioluminescence assay. Epitope mapping revealed the surface-exposed amino acid sequence (35-FGPDGENS-42) close to the catalytic center of CA XII recognized by the MAb 14D6. The variable regions of the heavy and light chains of MAb 14D6 were sequenced and their complementarity-determining regions were defined. The obtained variable sequences were used to generate recombinant antibodies in two formats: single-chain fragment variable (scFv) expressed in E. coli and scFv fused to human IgG1 Fc fragment (scFv-Fc) expressed in Chinese Hamster Ovary (CHO) cells. Both recombinant antibodies maintained the same specificity for CA XII as the parental MAb 14D6. The novel antibodies may represent promising tools for CA XII-related cancer research and immunotherapy.

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.

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.

Back-to-Germline (B2G) Procedure for Antibody Devolution

Bispecific antibodies (bsAbs) with avidity-enhanced specificity can be used to address target cells with increased specificity, ideally binding efficiently to cells that express two cognate antigens, yet not to cells that express only one of those. Building blocks required to generate such bsAbs are binders that recognize the two antigens with high specificity yet with various (including very low monovalent) affinities. The herein described ‘back-to-germline’ (B2G) procedure defines such derivatives. It converts parent antibodies with high specificity to derivatives that retain specificity but modulate affinity. The approach defines mutations to be introduced into antibody complementarity-determining regions (CDRs) regions without requiring structures of antibody-antigen complexes. Instead, it reverses the B-cell maturation process that increases affinities, with preference on CDR residues with high antigen contact probability. Placing germline residues at those positions generates VH and VL domains and Fv-combinations thereof that retain specificities but are ‘de-matured’ to different degrees. De-maturation influences on-rates and off-rates, and can produce entities with extremely low affinity for which binding can only be detected in bivalent formats. A comparison with alanine replacement in CDRs (so far, the most frequently applied technology) indicates that B2G may be more reliable/predictable without introduction of stickiness or poly-reactivity. The applicability for generating sets of affinity-modulated monospecific variants is exemplarily shown for antibodies that bind CD138, Her2/neu, and EGFR.

High-Throughput Antibody Structure Modeling and Design Using ABodyBuilder

Antibodies are proteins of the adaptive immune system; they can be designed to bind almost any molecule, and are increasingly being used as biotherapeutics. Experimental antibody design is an expensive and time-consuming process, and computational antibody design methods can now be used to help develop new therapeutics and diagnostics. Within the design pipeline, accurate antibody structure modeling is essential, as it provides the basis for antibody-antigen docking, binding affinity prediction, and estimating thermal stability. Ideally, models should be rapidly generated, allowing the exploration of the breadth of antibody space. This allows methods to replicate the natural processes of antibody diversification (e.g., V(D)J recombination and somatic hypermutation), and cope with large volumes of data that are typical of next-generation sequencing datasets. Here we describe ABodyBuilder and PEARS, algorithms that build and mutate antibody model structures. These methods take ~30 s to generate a model antibody structure.

Bioluminescent Antibodies for Point-of-Care Diagnostics

We introduce a general method to transform antibodies into ratiometric, bioluminescent sensor proteins for the no‐wash quantification of analytes. Our approach is based on the genetic fusion of antibody fragments to NanoLuc luciferase and SNAP‐tag, the latter being labeled with a synthetic fluorescent competitor of the antigen. Binding of the antigen, here synthetic drugs, by the sensor displaces the tethered fluorescent competitor from the antibody and disrupts bioluminescent resonance energy transfer (BRET) between the luciferase and fluorophore. The semisynthetic sensors display a tunable response range (submicromolar to submillimolar) and large dynamic range (ΔR_max>500 %), and they permit the quantification of analytes through spotting of the samples onto paper followed by analysis with a digital camera.

DynaDom: structure-based prediction of T cell receptor inter-domain and T cell receptor-peptide-MHC (class I) association angles

Background

T cell receptor (TCR) molecules are involved in the adaptive immune response as they distinguish between self- and foreign-peptides, presented in major histocompatibility complex molecules (pMHC). Former studies showed that the association angles of the TCR variable domains (Vα/Vβ) can differ significantly and change upon binding to the pMHC complex. These changes can be described as a rotation of the domains around a general Center of Rotation, characterized by the interaction of two highly conserved glutamine residues.

Methods

We developed a computational method, DynaDom, for the prediction of TCR Vα/Vβ inter-domain and TCR/pMHC orientations in TCRpMHC complexes, which allows predicting the orientation of multiple protein-domains. In addition, we implemented a new approach to predict the correct orientation of the carboxamide endgroups in glutamine and asparagine residues, which can also be used as an external, independent tool.

Results

The approach was evaluated for the remodeling of 75 and 53 experimental structures of TCR and TCRpMHC (class I) complexes, respectively. We show that the DynaDom method predicts the correct orientation of the TCR Vα/Vβ angles in 96 and 89% of the cases, for the poses with the best RMSD and best interaction energy, respectively. For the concurrent prediction of the TCR Vα/Vβ and pMHC orientations, the respective rates reached 74 and 72%. Through an exhaustive analysis, we could show that the pMHC placement can be further improved by a straightforward, yet very time intensive extension of the current approach.

Conclusions

The results obtained in the present remodeling study prove the suitability of our approach for interdomain-angle optimization. In addition, the high prediction rate obtained specifically for the energetically highest ranked poses further demonstrates that our method is a powerful candidate for blind prediction. Therefore it should be well suited as part of any accurate atomistic modeling pipeline for TCRpMHC complexes and potentially other large molecular assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-016-0071-7) contains supplementary material, which is available to authorized users.

Characterization of a re-engineered, mesothelin-targeted Pseudomonas exotoxin fusion protein for lung cancer therapy

Mesothelin overexpression in lung adenocarcinomas correlates with the presence of activating KRAS mutations and poor prognosis. Hence SS1P, a mesothelin-targeted immunotoxin, could offer valuable treatment options for these patients, but its use in solid tumor therapy is hampered by high immunogenicity and non-specific toxicity. To overcome both obstacles we developed RG7787, a de-immunized cytotoxic fusion protein comprising a humanized SS1 Fab fragment and a truncated, B-cell epitope silenced, 24 kD fragment of Pseudomonas exotoxin A (PE24). Reactivity of RG7787 with sera from immunotoxin-treated patients was >1000 fold reduced. In vitro RG7787 inhibited cell viability of lung cancer cell lines with picomolar potency. The pharmacokinetic properties of RG7787 in rodents were comparable to SS1P, yet it was tolerated up to 10 fold better without causing severe vascular leak syndrome or hepatotoxicity. A pharmacokinetic/pharmacodynamic model developed based on NCI-H596 xenograft studies showed that for RG7787 and SS1P, their in vitro and in vivo potencies closely correlate. At optimal doses of 2-3 mg/kg RG7787 is more efficacious than SS1P. Even large, well established tumors (600 mm(3)) underwent remission during three treatment cycles with RG7787. Also in two patient-derived lung cancer xenograft models, Lu7336 and Lu7187, RG7787 showed anti-tumor efficacy. In monotherapy two treatment cycles were moderately efficacious in the Lu7336 model but showed good anti-tumor activity in the KRAS mutant Lu7187 model (26% and 80% tumor growth inhibition, respectively). Combination of RG7787 with standard chemotherapies further enhanced efficacy in both models achieving near complete eradication of Lu7187 tumors.

ABodyBuilder: Automated antibody structure prediction with data-driven accuracy estimation

Computational modeling of antibody structures plays a critical role in therapeutic antibody design. Several antibody modeling pipelines exist, but no freely available methods currently model nanobodies, provide estimates of expected model accuracy, or highlight potential issues with the antibody's experimental development. Here, we describe our automated antibody modeling pipeline, ABodyBuilder, designed to overcome these issues. The algorithm itself follows the standard 4 steps of template selection, orientation prediction, complementarity-determining region (CDR) loop modeling, and side chain prediction. ABodyBuilder then annotates the 'confidence' of the model as a probability that a component of the antibody (e.g., CDRL3 loop) will be modeled within a root-mean square deviation threshold. It also flags structural motifs on the model that are known to cause issues during in vitro development. ABodyBuilder was tested on 4 separate datasets, including the 11 antibodies from the Antibody Modeling Assessment-II competition. ABodyBuilder builds models that are of similar quality to other methodologies, with sub-Angstrom predictions for the 'canonical' CDR loops. Its ability to model nanobodies, and rapidly generate models (∼30 seconds per model) widens its potential usage. ABodyBuilder can also help users in decision-making for the development of novel antibodies because it provides model confidence and potential sequence liabilities. ABodyBuilder is freely available at http://opig.stats.ox.ac.uk/webapps/abodybuilder .