Antibodies with Weakly Basic Isoelectric Points Minimize Trade-offs between Formulation and Physiological Colloidal Properties

Reduction of monoclonal antibody viscosity using interpretable machine learning

Early identification of antibody candidates with drug-like properties is essential for simplifying the development of safe and effective antibody therapeutics. For subcutaneous administration, it is important to identify candidates with low self-association to enable their formulation at high concentration while maintaining low viscosity, opalescence, and aggregation. Here, we report an interpretable machine learning model for predicting antibody (IgG1) variants with low viscosity using only the sequences of their variable (Fv) regions. Our model was trained on antibody viscosity data (>100 mg/mL mAb concentration) obtained at a common formulation pH (pH 5.2), and it identifies three key Fv features of antibodies linked to viscosity, namely their isoelectric points, hydrophobic patch sizes, and numbers of negatively charged patches. Of the three features, most predicted antibodies at risk for high viscosity, including antibodies with diverse antibody germlines in our study (79 mAbs) as well as clinical-stage IgG1s (94 mAbs), are those with low Fv isoelectric points (Fv pIs < 6.3). Our model identifies viscous antibodies with relatively high accuracy not only in our training and test sets, but also for previously reported data. Importantly, we show that the interpretable nature of the model enables the design of mutations that significantly reduce antibody viscosity, which we confirmed experimentally. We expect that this approach can be readily integrated into the drug development process to reduce the need for experimental viscosity screening and improve the identification of antibody candidates with drug-like properties.

Human Dose and Pharmacokinetic Predictions for Biologics at Boehringer Ingelheim: A Retrospective Analysis

INTRODUCTION

Accurate predictions of pharmacokinetics and efficacious doses for biologics in humans are critical for selecting appropriate first-in-human starting doses and dose ranges and for estimating clinical material needs and cost of goods. This also impacts clinical feasibility, particularly for subcutaneously administered biologics.

METHODS

We performed a comprehensive comparison between predicted and observed clearances and doses in humans for a set of 22 biologic drugs developed at Boehringer Ingelheim (BI) over the last 2 decades. The analysis included biologics across three therapeutic areas comprising a wide variety of modalities: mono- and bispecific monoclonal antibodies (mAbs) and nanobodies and a Fab fragment.

RESULTS

Our analysis showed that observed clearances in humans were within twofold of predicted clearances for 17 out of 20 biologics (85%). Six biologics had uncharacteristically high observed human clearances (range 32-280 mL/h) for their respective molecular classes, impacting their clinical developability. For three molecules, molecular characteristics contributed to the high clearance. Clinically selected doses were within twofold of predicted for 58% of projects. With 42% and 25% of projects selecting clinical doses higher than two- or threefold the predicted value, respectively, the importance of better understanding not only the pharmacokinetic (PK) but also the predictivity of pharmacodynamic models is highlighted.

CONCLUSIONS

We provide a clinical pharmacology perspective on the commonly accepted twofold range of human clearance predictions as well as the implications of higher than predicted targeted efficacious plasma concentration on clinical development. Finally, an analysis of key success factors for biologics at BI was conducted, which may be relevant for the entire pharmaceutical industry. This is one of the largest retrospective analyses for biologics and provides further evidence that successful predictions of human PK and efficacious dose will be further facilitated by gathering key translational data early in research.

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤ 2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~ 200 anti-GCGR IgG1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2-10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

Optimization of therapeutic antibodies for reduced self-association and non-specific binding via interpretable machine learning

Antibody development, delivery, and efficacy are influenced by antibody-antigen affinity interactions, off-target interactions that reduce antibody bioavailability and pharmacokinetics, and repulsive self-interactions that increase the stability of concentrated antibody formulations and reduce their corresponding viscosity. Yet identifying antibody variants with optimal combinations of these three types of interactions is challenging. Here we show that interpretable machine-learning classifiers, leveraging antibody structural features descriptive of their variable regions and trained on experimental data for a panel of 80 clinical-stage monoclonal antibodies, can identify antibodies with optimal combinations of low off-target binding in a common physiological-solution condition and low self-association in a common antibody-formulation condition. For three clinical-stage antibodies with suboptimal combinations of off-target binding and self-association, the classifiers predicted variable-region mutations that optimized non-affinity interactions while maintaining high-affinity antibody-antigen interactions. Interpretable machine-learning models may facilitate the optimization of antibody candidates for therapeutic applications.

The origins of nonideality exhibited by monoclonal antibodies and Fab fragments in human serum

The development of therapeutic antibodies remains challenging, time-consuming, and expensive. A key contributing factor is a lack of understanding of how proteins are affected by complex biological environments such as serum and plasma. Nonideality due to attractive or repulsive interactions with cosolutes can alter the stability, aggregation propensity, and binding interactions of proteins in solution. Fluorescence correlation spectroscopy (FCS) can be used to measure apparent second virial coefficient (B2,app ) values for therapeutic and model monoclonal antibodies (mAbs) that capture the nature and strength of interactions with cosolutes directly in undiluted serum and similar complex biological media. Here, we use FCS-derived B2,app measurements to identify the components of human serum responsible for nonideal interactions with mAbs and Fab fragments. Most mAbs exhibit neutral or slightly attractive interactions with intact serum. Generally, mAbs display repulsive interactions with albumin and mildly attractive interactions with IgGs in the context of whole serum. Crucially, however, these attractive interactions are much stronger with pooled IgGs isolated from other serum components, indicating that the effects of serum nonideality can only be understood by studying the intact medium (rather than isolated components). Moreover, Fab fragments universally exhibited more attractive interactions than their parental mAbs, potentially rendering them more susceptible to nonideality-driven perturbations. FCS-based B2,app measurements have the potential to advance our understanding of how physiological environments impact protein-based therapeutics in general. Furthermore, incorporating such assays into preclinical biologics development may help de-risk molecules and make for a faster and more efficient development process.

The role of acidification of the medium on the erythrocyte agglutinating and adsorbing properties of blood group-specific monoclonal antibodies

Background

Studies on the specificity of ABH antigen-antibody interactions at different pH values are rare. Therefore, the aim of this study was to estimate the effect of acidification of the reacting medium on the agglutinating ability of anti-A monoclonal antibodies (mABs) and their inhibition by glycoconjugates of red blood cell membranes.

Methods

Anti-A mABs were obtained from the Fourth International Workshop on Monoclonal Antibodies Against Human Red Blood Cell and Related Antigens (on July 19-20, 2002, in Paris, France). The glycoconjugates were isolated from erythrocytes' membranes. The inhibition of the lipid and protein isotypes of the blood group A antigen was assessed.

Results

The acidic medium decreased the agglutinating ability of acid immunoglobulin M (IgM) anti-A mABs, in contrast to alkaline immunoglobulin G (IgG) mABs. Meanwhile, at pH 6.5, IgM antibodies exhibited a high adsorption capacity, while IgG antibodies demonstrated a strong adsorption capacity at an alkaline pH. mABs 2-19, 2-28, 2-22, and 2-8 were inhibited by the acidic lipid and protein glycotopes of erythrocyte membranes. However, mAB 2-8 was inhibited by both acidic and alkaline glycotope variants.

Conclusions

The agglutinating and adsorbing abilities of mABs, which are revealed at opposing pH values, should be taken into account during studies of antigen-antibody interactions.

AI/ML combined with next-generation sequencing of VHH immune repertoires enables the rapid identification of de novo humanized and sequence-optimized single domain antibodies: a prospective case study

Introduction: In this study, we demonstrate the feasibility of yeast surface display (YSD) and nextgeneration sequencing (NGS) in combination with artificial intelligence and machine learning methods (AI/ML) for the identification of de novo humanized single domain antibodies (sdAbs) with favorable early developability profiles. Methods: The display library was derived from a novel approach, in which VHH-based CDR3 regions obtained from a llama (Lama glama), immunized against NKp46, were grafted onto a humanized VHH backbone library that was diversified in CDR1 and CDR2. Following NGS analysis of sequence pools from two rounds of fluorescence-activated cell sorting we focused on four sequence clusters based on NGS frequency and enrichment analysis as well as in silico developability assessment. For each cluster, long short-term memory (LSTM) based deep generative models were trained and used for the in silico sampling of new sequences. Sequences were subjected to sequence- and structure-based in silico developability assessment to select a set of less than 10 sequences per cluster for production. Results: As demonstrated by binding kinetics and early developability assessment, this procedure represents a general strategy for the rapid and efficient design of potent and automatically humanized sdAb hits from screening selections with favorable early developability profiles.

EPR Viscometric Measurements Using a 13C-Labeled Triarylmethyl Radical in Protein-based Biotherapeutics and Human Synovial Fluids

The viscosity measurements are of clinical significance for evaluation of the potential pathological conditions of biological lubricants such as synovial fluids of joints, and for formulation and characterization of peptide- and protein-based biotherapeutics. Due to inherent potential therapeutic activity, protein drugs have proven to be one of the most efficient therapeutic agents in treatment of several life-threatening disorders, such as diabetes and autoimmune diseases. However, home-use applications for treating chronic inflammatory diseases, such as diabetes and rheumatoid arthritis, necessitate the development of high-concentration insulin and monoclonal antibodies formulations for patient self-administration. High protein concentrations can affect viscosity of the corresponding drug solutions complicating their manufacture and administration. The measurements of the viscosity of new insulin analogs and monoclonal antibodies solutions under development is of practical importance to avoid unwanted highly viscous, and therefore, painful for injection drug formulations. Recently, we have demonstrated capability of the electron paramagnetic resonance (EPR) viscometry using viscosity-sensitive 13C-labeled trityl spin probe (13C1-dFT) to report the viscosity of human blood, and interstitial fluids measured in various organs in mice ex-vivo and in anesthetized mice, in vivo. In the present work, we demonstrate utility of the EPR viscometry using 13C1-dFT to measure microviscosity of commercial insulin samples, antibodies solution, and human synovial fluids using small microliter volume samples (5-50 μL). This viscometry analysis approach provides useful tool to control formulations and administration of new biopharmaceuticals, and for evaluation of the state of synovial fluids of importance for clinical applications.

Physiologically Based Pharmacokinetic Modeling to Characterize the Effect of Molecular Charge on Whole-Body Disposition of Monoclonal Antibodies

Motivated by a series of work demonstrating the effect of molecular charge on antibody pharmacokinetics (PK), physiological-based pharmacokinetic (PBPK) models are emerging that relate in silico calculated charge or in vitro measures of polyspecificity to antibody PK parameters. However, only plasma data has been used for model development in these studies, leading to unvalidated assumptions. Here, we present an extended platform PBPK model for antibodies that incorporate charge-dependent endothelial cell pinocytosis rate and nonspecific off-target binding in the interstitial space and on circulating blood cells, to simultaneously characterize whole-body disposition of three antibody charge variants. Predictive potential of various charge metrics was also explored, and the difference between positive charge patches and negative charge patches (i.e., PPC-PNC) was used as the charge parameter to establish quantitative relationships with nonspecific binding affinities and endothelial cell uptake rate. Whole-body disposition of these charge variants was captured well by the model, with less than 2-fold predictive error in area under the curve of most plasma and tissue PK data. The model also predicted that with greater positive charge, nonspecific binding was more substantial, and pinocytosis rate increased especially in brain, heart, kidney, liver, lung, and spleen, but remained unchanged in adipose, bone, muscle, and skin. The presented PBPK model contributes to our understanding of the mechanisms governing the disposition of charged antibodies and can be used as a platform to guide charge engineering based on desired plasma and tissue exposures.

Design of a highly sensitive and versatile membrane-based immunosensor using a Cu-free click reaction

A highly sensitive immunosensor is developed using membrane pores as the recognition interface. In this sensor, a Cu-free click reaction is used to efficiently immobilize antibodies, and the sensor inhibits the adsorption of nonspecific proteins that degrade sensitivity. Furthermore, the sensor demonstrates rapid interleukin-6 detection in the picogram per milliliter range.

Characterizing Experimental Monoclonal Antibody Interactions and Clustering Using a Coarse-Grained Simulation Library and a Viscosity Model

Attractive protein-protein interactions in concentrated monoclonal antibody (mAb) solutions may lead to the formation of clusters that increase viscosity. Here, we propose an analytical model that relates mAb solution viscosity to clustering by accounting for the contributions of suboptimal mAb packing within a cluster and cluster fractal dimension. The influence of short-range, anisotropic attractions and long-range Coulombic repulsion on cluster properties is investigated by analyzing the cluster-size distributions, cluster fractal dimensions, radial distribution functions, and static structure factors from a library of coarse-grained molecular dynamics simulations. The library spans a vast range of mAb charges and attractive interactions in solutions of varying ionic strength. We present a framework for combining the viscosity model and simulation library to successfully characterize the attraction, repulsion, and clustering of an experimental mAb in three different pH and cosolute conditions by fitting the measured viscosity or structure factor from small-angle X-ray scattering. At low ionic strength, the cluster-size distribution is impacted by strong charges, and both the viscosity and net charge or structure factor and net charge must be considered to deconvolute the effects of short-range attraction and long-range repulsion.

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.

SUMO: In Silico Sequence Assessment Using Multiple Optimization Parameters

To select the most promising screening hits from antibody and VHH display campaigns for subsequent in-depth profiling and optimization, it is highly desirable to assess and select sequences on properties beyond only their binding signals from the sorting process. In addition, developability risk criteria, sequence diversity, and the anticipated complexity for sequence optimization are relevant attributes for hit selection and optimization. Here, we describe an approach for the in silico developability assessment of antibody and VHH sequences. This method not only allows for ranking and filtering multiple sequences with regard to their predicted developability properties and diversity, but also visualizes relevant sequence and structural features of potentially problematic regions and thereby provides rationales and starting points for multi-parameter sequence optimization.

Embedding Dynamics in Intrinsic Physicochemical Profiles of Market-Stage Antibody-Based Biotherapeutics

Adequate stability, manufacturability, and safety are crucial to bringing an antibody-based biotherapeutic to the market. Following the concept of holistic in silico developability, we introduce a physicochemical description of 91 market-stage antibody-based biotherapeutics based on orthogonal molecular properties of variable regions (Fvs) embedded in different simulation environments, mimicking conditions experienced by antibodies during manufacturing, formulation, and in vivo. In this work, the evaluation of molecular properties includes conformational flexibility of the Fvs using molecular dynamics (MD) simulations. The comparison between static homology models and simulations shows that MD significantly affects certain molecular descriptors like surface molecular patches. Moreover, the structural stability of a subset of Fv regions is linked to changes in their specific molecular interactions with ions in different experimental conditions. This is supported by the observation of differences in protein melting temperatures upon addition of NaCl. A DEvelopability Navigator In Silico (DENIS) is proposed to compare mAb candidates for their similarity with market-stage biotherapeutics in terms of physicochemical properties and conformational stability. Expanding on our previous developability guidelines (Ahmed et al. Proc. Natl. Acad. Sci. 2021, 118 (37), e2020577118), the hydrodynamic radius and the protein strand ratio are introduced as two additional descriptors that enable a more comprehensive in silico characterization of biotherapeutic drug candidates. Test cases show how this approach can facilitate identification and optimization of intrinsically developable lead candidates. DENIS represents an advanced computational tool to progress biotherapeutic drug candidates from discovery into early development by predicting drug properties in different aqueous environments.

Reduction of therapeutic antibody self-association using yeast-display selections and machine learning

Self-association governs the viscosity and solubility of therapeutic antibodies in high-concentration formulations used for subcutaneous delivery, yet it is difficult to reliably identify candidates with low self-association during antibody discovery and early-stage optimization. Here, we report a high-throughput protein engineering method for rapidly identifying antibody candidates with both low self-association and high affinity. We find that conjugating quantum dots to IgGs that strongly self-associate (pH 7.4, PBS), such as lenzilumab and bococizumab, results in immunoconjugates that are highly sensitive for detecting other high self-association antibodies. Moreover, these conjugates can be used to rapidly enrich yeast-displayed bococizumab sub-libraries for variants with low levels of immunoconjugate binding. Deep sequencing and machine learning analysis of the enriched bococizumab libraries, along with similar library analysis for antibody affinity, enabled identification of extremely rare variants with co-optimized levels of low self-association and high affinity. This analysis revealed that co-optimizing bococizumab is difficult because most high-affinity variants possess positively charged variable domains and most low self-association variants possess negatively charged variable domains. Moreover, negatively charged mutations in the heavy chain CDR2 of bococizumab, adjacent to its paratope, were effective at reducing self-association without reducing affinity. Interestingly, most of the bococizumab variants with reduced self-association also displayed improved folding stability and reduced nonspecific binding, revealing that this approach may be particularly useful for identifying antibody candidates with attractive combinations of drug-like properties.Abbreviations: AC-SINS: affinity-capture self-interaction nanoparticle spectroscopy; CDR: complementarity-determining region; CS-SINS: charge-stabilized self-interaction nanoparticle spectroscopy; FACS: fluorescence-activated cell sorting; Fab: fragment antigen binding; Fv: fragment variable; IgG: immunoglobulin; QD: quantum dot; PBS: phosphate-buffered saline; V_H: variable heavy; V_L: variable light.