Modeling and mitigation of high-concentration antibody viscosity through structure-based computer-aided protein design

A framework for the biophysical screening of antibody mutations targeting solvent-accessible hydrophobic and electrostatic patches for enhanced viscosity profiles

The formulation of high-concentration monoclonal antibody (mAb) solutions in low dose volumes for autoinjector devices poses challenges in manufacturability and patient administration due to elevated solution viscosity. Often many therapeutically potent mAbs are discovered, but their commercial development is stalled by unfavourable developability challenges. In this work, we present a systematic experimental framework for the computational screening of molecular descriptors to guide the design of 24 mutants with modified viscosity profiles accompanied by experimental evaluation. Our experimental observations using a model anti-IL8 mAb and eight engineered mutant variants reveal that viscosity reduction is influenced by the location of hydrophobic interactions, while targeting positively charged patches significantly increases viscosity in comparison to wild-type anti-IL-8 mAb. We conclude that most predicted in silico physicochemical properties exhibit poor correlation with measured experimental parameters for antibodies with suboptimal developability characteristics, emphasizing the need for comprehensive case-by-case evaluation of mAbs. This framework combining molecular design and triage via computational predictions with experimental evaluation aids the agile and rational design of mAbs with tailored solution viscosities, ensuring improved manufacturability and patient convenience in self-administration scenarios.

Intramodule pressure profiles and protein accumulation during tangential flow filtration

Tangential flow filtration (TFF) through a 30 kDa nominal molecular weight cut-off (MWCO) ultrafiltration membrane is widely employed to concentrate purified monoclonal antibodies (mAbs) to levels required for their formulation into injectable biologics. While TFF has been used to remove casein from milk for cheese production for over 35 years, and in pharmaceutical manufacture of biotherapeutic proteins for 20 years, the rapid decline in filtration rate (i.e., flux) at high protein concentrations is a limitation that still needs to be addressed. This is particularly important for mAbs, many of which are 140-160 kDa immunoglobulin G (IgG) type proteins recovered at concentrations of 200 mg/mL or higher. This work reports the direct measurement of local transmembrane pressure drops and off-line confocal imaging of protein accumulation in stagnant regions on the surface of a 30 kDa regenerated cellulose membrane in a flat-sheet configuration widely used in manufacture of biotherapeutic proteins. These first-of-a-kind measurements using 150 kDa bovine IgG show that while axial pressure decreases by 58 psi across a process membrane cassette, the decrease in transmembrane pressure drop is constant at about 1.2 psi/cm along the 20.7 cm length of the membrane. Confocal laser scanning microscopy of the membrane surface at the completion of runs where retentate protein concentration exceeds 200 mg/mL, shows a 50 μm thick protein layer is uniformly deposited. The localized measurements made possible by the modified membrane system confirm the role of protein deposition on limiting ultrafiltration rate and indicate possible targets for improving membrane performance.

Molecular surface descriptors to predict antibody developability: sensitivity to parameters, structure models, and conformational sampling

In silico assessment of antibody developability during early lead candidate selection and optimization is of paramount importance, offering a rapid and material-free screening approach. However, the predictive power and reproducibility of such methods depend heavily on the selection of molecular descriptors, model parameters, accuracy of predicted structure models, and conformational sampling techniques. Here, we present a set of molecular surface descriptors specifically designed for predicting antibody developability. We assess the performance of these descriptors by benchmarking their correlations with an extensive array of experimentally determined biophysical properties, including viscosity, aggregation, hydrophobic interaction chromatography, human pharmacokinetic clearance, heparin retention time, and polyspecificity. Further, we investigate the sensitivity of these surface descriptors to methodological nuances, such as the choice of interior dielectric constant, hydrophobicity scales, structure prediction methods, and the impact of conformational sampling. Notably, we observe systematic shifts in the distribution of surface descriptors depending on the structure prediction method used, driving weak correlations of surface descriptors across structure models. Averaging the descriptor values over conformational distributions from molecular dynamics mitigates the systematic shifts and improves the consistency across different structure prediction methods, albeit with inconsistent improvements in correlations with biophysical data. Based on our benchmarking analysis, we propose six in silico developability risk flags and assess their effectiveness in predicting potential developability issues for a set of case study molecules.

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.

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.

Characterizing Protein-Protein Interactions and Viscosity of a Monoclonal Antibody from Low to High Concentration Using Small-Angle X-ray Scattering and Molecular Dynamics Simulations

Understanding protein-protein interactions and formation of reversible oligomers (clusters) in concentrated monoclonal antibody (mAb) solutions is necessary for designing stable, low viscosity (η) concentrated formulations for processing and subcutaneous injection. Here we characterize the strength (K) of short-range anisotropic attractions (SRA) for 75-200 mg/mL mAb2 solutions at different pH and cosolute conditions by analyzing structure factors (Seff(q)) from small-angle X-ray scattering (SAXS) using coarse-grained molecular dynamics simulations. Best fit simulations additionally provide cluster size distributions, fractal dimensions, cluster occluded volume, and mAb coordination numbers. These equilibrium properties are utilized in a model to account for increases in viscosity caused by occluded volume in the clusters (packing effects) and dissipation of stress across lubricated fractal clusters. Seff(q) is highly sensitive to K at 75 mg/mL where mAbs can mutually align to form SRA contacts but becomes less sensitive at 200 mg/mL as steric repulsion due to packing becomes dominant. In contrast, η at 200 mg/mL is highly sensitive to SRA and the average cluster size from SAXS/simulation, which is observed to track the cluster relaxation time from shear thinning. By analyzing the distribution of sub-bead hot spots on the 3D mAb surface, we identify a strongly attractive hydrophobic patch in the complementarity determining region (CDR) at pH 4.5 that contributes to the high K and consequently large cluster sizes and high η. Adding NaCl screens electrostatic interactions and increases the impact of hydrophobic attraction on cluster size and raises η, whereas nonspecific binding of Arg attenuates all SRA, reducing η. The hydrophobic patch is absent at higher pH values, leading to smaller K, smaller clusters, and lower η. This work constitutes a first attempt to use SAXS and CG modeling to link both structural and rheological properties of concentrated mAb solutions to the energetics of specific hydrophobic patches on mAb surfaces. As such, our work opens an avenue for future research, including the possibility of designing coarse-grained models with physically meaningful interacting hot spots.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.

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.

Subclass Effects on Self-Association and Viscosity of Monoclonal Antibodies at High Concentrations

The effects of a subclass of monoclonal antibodies (mAbs) on protein-protein interactions, formation of reversible oligomers (clusters), and viscosity (η) are not well understood at high concentrations. Herein, we quantify a short-range anisotropic attraction between the complementarity-determining region (CDR) and CH3 domains (K_CDR-CH3) for vedolizumab IgG1, IgG2, or IgG4 subclasses by fitting small-angle X-ray scattering (SAXS) structure factor S_eff(q) data with an extensive library of 12-bead coarse-grained (CG) molecular dynamics simulations. The K_CDR-CH3 bead attraction strength was isolated from the strength of long-range electrostatic repulsion for the full mAb, which was determined from the theoretical net charge and a scaling parameter ψ to account for solvent accessibility and ion pairing. At low ionic strength (IS), the strongest short-range attraction (K_CDR-CH3) and consequently the largest clusters and highest η were observed with IgG1, the subclass with the most positively charged CH3 domain. Furthermore, the trend in K_CDR-CH3 with the subclass followed the electrostatic interaction energy between the CDR and CH3 regions calculated with the BioLuminate software using the 3D mAb structure and molecular interaction potentials. Whereas the equilibrium cluster size distributions and fractal dimensions were determined from fits of SAXS with the MD simulations, the degree of cluster rigidity under flow was estimated from the experimental η with a phenomenological model. For the systems with the largest clusters, especially IgG1, the inefficient packing of mAbs in the clusters played the largest role in increasing η, whereas for other systems, the relative contribution from stress produced by the clusters was more significant. The ability to relate η to short-range attraction from SAXS measurements at high concentrations and to theoretical characterization of electrostatic patches on the 3D surface is not only of fundamental interest but also of practical value for mAb discovery, processing, formulation, and subcutaneous delivery.

Accelerating therapeutic protein design with computational approaches toward the clinical stage

Therapeutic protein, represented by antibodies, is of increasing interest in human medicine. However, clinical translation of therapeutic protein is still largely hindered by different aspects of developability, including affinity and selectivity, stability and aggregation prevention, solubility and viscosity reduction, and deimmunization. Conventional optimization of the developability with widely used methods, like display technologies and library screening approaches, is a time and cost-intensive endeavor, and the efficiency in finding suitable solutions is still not enough to meet clinical needs. In recent years, the accelerated advancement of computational methodologies has ushered in a transformative era in the field of therapeutic protein design. Owing to their remarkable capabilities in feature extraction and modeling, the integration of cutting-edge computational strategies with conventional techniques presents a promising avenue to accelerate the progression of therapeutic protein design and optimization toward clinical implementation. Here, we compared the differences between therapeutic protein and small molecules in developability and provided an overview of the computational approaches applicable to the design or optimization of therapeutic protein in several developability issues.

Low-data interpretable deep learning prediction of antibody viscosity using a biophysically meaningful representation

Deep learning, aided by the availability of big data sets, has led to substantial advances across many disciplines. However, many scientific problems of practical interest lack sufficiently large datasets amenable to deep learning. Prediction of antibody viscosity is one such problem where deep learning methods have not yet been explored due to the relative scarcity of relevant training data. In this work, we overcome this limitation using a biophysically meaningful representation that enables us to develop generalizable models even under limited training data. We present, PfAbNet-viscosity, a 3D convolutional neural network architecture, to predict high-concentration viscosity of therapeutic antibodies. We show that with the electrostatic potential surface of the antibody variable region as the only input to the network, the models trained on as few as couple dozen datapoints can generalize with high accuracy. Our feature attribution analysis shows that PfAbNet-viscosity has learned key biophysical drivers of viscosity. The applicability of our approach to other biological systems is discussed.

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.

Predictive modeling of concentration-dependent viscosity behavior of monoclonal antibody solutions using artificial neural networks

Solutions of monoclonal antibodies (mAbs) can show increased viscosity at high concentration, which can be a disadvantage during protein purification, filling, and administration. The viscosity is determined by protein-protein-interactions, which are influenced by the antibody's sequence as well as solution conditions, like pH, buffer type, or the presence of salts and other excipients. To predict viscosity, experimental parameters, like the diffusion interaction parameter (kD), or computational tools harnessing information derived from primary sequence, are often used, but a reliable predictive tool is still missing. We present a modeling approach employing artificial neural networks (ANNs) using experimental factors combined with simulation-derived parameters plus viscosity data from 27 highly concentrated (180 mg/mL) mAbs. These ANNs can be used to predict if mAbs exhibit problematic viscosity at distinct concentrations or to model viscosity-concentration-curves.

Structure-Based Optimization of Antibody-Based Biotherapeutics for Improved Developability: A Practical Guide for Molecular Modelers

A great effort to avoid known developability risks is now more often being made earlier during the lead candidate discovery and optimization phase of biotherapeutic drug development. Predictive computational strategies, used in the early stages of antibody discovery and development, to mitigate the risk of late-stage failure of antibody candidates, are highly valuable. Various structure-based methods exist for accurately predicting properties critical to developability, and, in this chapter, we discuss the history of their development and demonstrate how they can be used to filter large sets of candidates arising from target affinity screening and to optimize lead candidates for developability. Methods for modeling antibody structures from sequence and detecting post-translational modifications and chemical degradation liabilities are also discussed.

Non-specificity as the sticky problem in therapeutic antibody development

Antibodies are highly potent therapeutic scaffolds with more than a hundred different products approved on the market. Successful development of antibody-based drugs requires a trade-off between high target specificity and target binding affinity. In order to better understand this problem, we here review non-specific interactions and explore their fundamental physicochemical origins. We discuss the role of surface patches - clusters of surface-exposed amino acid residues with similar physicochemical properties - as inducers of non-specific interactions. These patches collectively drive interactions including dipole-dipole, π-stacking and hydrophobic interactions to complementary moieties. We elucidate links between these supramolecular assembly processes and macroscopic development issues, such as decreased physical stability and poor in vivo half-life. Finally, we highlight challenges and opportunities for optimizing protein binding specificity and minimizing non-specificity for future generations of therapeutics.

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

The widespread interest in antibody therapeutics has led to much focus on identifying antibody candidates with favorable developability properties. In particular, there is broad interest in identifying antibody candidates with highly repulsive self-interactions in standard formulations (e.g., low ionic strength buffers at pH 5-6) for high solubility and low viscosity. Likewise, there is also broad interest in identifying antibody candidates with low levels of non-specific interactions in physiological solution conditions (PBS, pH 7.4) to promote favorable pharmacokinetic properties. To what extent antibodies that possess both highly repulsive self-interactions in standard formulations and weak non-specific interactions in physiological solution conditions can be systematically identified remains unclear and is a potential impediment to successful therapeutic drug development. Here, we evaluate these two properties for 42 IgG1 variants based on the variable fragments (Fvs) from four clinical-stage antibodies and complementarity-determining regions from 10 clinical-stage antibodies. Interestingly, we find that antibodies with the strongest repulsive self-interactions in a standard formulation (pH 6 and 10 mM histidine) display the strongest non-specific interactions in physiological solution conditions. Conversely, antibodies with the weakest non-specific interactions under physiological conditions display the least repulsive self-interactions in standard formulations. This behavior can be largely explained by the antibody isoelectric point, as highly basic antibodies that are highly positively charged under standard formulation conditions (pH 5-6) promote repulsive self-interactions that mediate high colloidal stability but also mediate strong non-specific interactions with negatively charged biomolecules at physiological pH and vice versa for antibodies with negatively charged Fv regions. Therefore, IgG1s with weakly basic isoelectric points between 8 and 8.5 and Fv isoelectric points between 7.5 and 9 typically display the best combinations of strong repulsive self-interactions and weak non-specific interactions. We expect that these findings will improve the identification and engineering of antibody candidates with drug-like biophysical properties.

Separating clinical antibodies from repertoire antibodies, a path to in silico developability assessment

Approaches for antibody discovery have seen substantial improvement and success in recent years. Yet, advancing antibodies into the clinic remains difficult because therapeutic developability concerns are challenging to predict. We developed a computational model to simplify antibody developability assessment and enable accelerated early-stage screening. To this end, we quantified the ability of hundreds of sequence- and structure-based descriptors to differentiate clinical antibodies that have undergone rigorous screening and characterization for drug-like properties from antibodies in the human repertoire that are not natively paired. This analysis identified 144 descriptors capable of distinguishing clinical from repertoire antibodies. Five descriptors were selected and combined based on performance and orthogonality into a single model referred to as the Therapeutic Antibody Developability Analysis (TA-DA). On a hold-out test set, this tool separated clinical antibodies from repertoire antibodies with an AUC = 0.8, demonstrating the ability to identify developability attributes unique to clinical antibodies. Based on our results, the TA-DA score may serve as an approach for selecting lead antibodies for further development.

Abbreviations: Affinity-Capture Self-Interaction Nanoparticle Spectroscopy (AC-SINS), Area Under the Curve (AUC), Complementary-Determining Region (CDR), Clinical-Stage Therapeutics (CST), Framework (FR), Monoclonal Antibodies (mAbs), Observed Antibody Space (OAS), Receiver Operating Characteristic (ROC), Size-Exclusion Chromatography (SEC), Structural Aggregation Propensity (SAP), Therapeutic Antibody Developability Analysis (TA-DA), Therapeutic Antibody Profiler (TAP), Therapeutic Structural Antibody Database (Thera-SAbDab), Variable Heavy (VH), Variable Light (VL).

Ongoing Challenges to Develop High Concentration Monoclonal Antibody-based Formulations for Subcutaneous Administration: Quo Vadis?

Although many subcutaneously (s.c.) delivered, high-concentration antibody formulations (HCAF) have received regulatory approval and are widely used commercially, formulation scientists are still presented with many ongoing challenges during HCAF development with new mAb and mAb-based candidates. Depending on the specific physicochemical and biological properties of a particular mAb-based molecule, such challenges vary from pharmaceutical attributes e.g., stability, viscosity, manufacturability, to clinical performance e.g., bioavailability, immunogenicity, and finally to patient experience e.g., preference for s.c. vs. intravenous delivery and/or preferred interactions with health-care professionals. This commentary focuses on one key formulation obstacle encountered during HCAF development: how to maximize the dose of the drug? We examine methodologies for increasing the protein concentration, increasing the volume delivered, or combining both approaches together. We discuss commonly encountered hurdles, i.e., physical protein instability and solution volume limitations, and we provide recommendations to formulation scientists to facilitate their development of s.c. administered HCAF with new mAb-based product candidates.

Structure-based charge calculations for predicting isoelectric point, viscosity, clearance, and profiling antibody therapeutics

The effect of hydrophobicity on antibody aggregation is well understood, and it has been shown that charge calculations can be useful for high-concentration viscosity and pharmacokinetic (PK) clearance predictions. In this work, structure-based charge descriptors are evaluated for their predictive performance on recently published antibody pI, viscosity, and clearance data. From this, we devised four rules for therapeutic antibody profiling which address developability issues arising from hydrophobicity and charged-based solution behavior, PK, and the ability to enrich for those that are approved by the U.S. Food and Drug Administration. Differences in strategy for optimizing the solution behavior of human IgG1 antibodies versus the IgG2 and IgG4 isotypes and the impact of pH alterations in formulation are discussed.

Identifying biophysical assays and in silico properties that enrich for slow clearance in clinical-stage therapeutic antibodies

Understanding the pharmacokinetic (PK) properties of a drug, such as clearance, is a crucial step for evaluating efficacy. The PK of therapeutic antibodies can be complex and is influenced by interactions with the target, Fc-receptors, anti-drug antibodies, and antibody intrinsic factors. A growing body of literature has linked biophysical properties of antibodies, particularly nonspecific-binding propensity, hydrophobicity and charged regions to rapid clearance in preclinical species and selected human PK studies. A clear understanding of the connection between biophysical properties and their impact on PK would allow for early selection and optimization of antibodies and reduce costly attrition during clinical trials due to sub-optimal human clearance. Due to the difficulty in obtaining large and unbiased human PK data, previous studies have focused mostly on preclinical PK. For this study, we obtained and curated the most comprehensive clinical PK dataset to date and calculated accurate estimates of linear clearance for 64 monoclonal antibodies ranging from investigational candidates in Phase 2 trials to marketed products. This allows for the first time a deep analysis of the influence of biophysical and sequence-based in silico properties directly on human clearance. We use statistical analysis and a Random Forest classifier to identify properties that have the greatest influence in our dataset. Our findings indicate that in vitro poly-specificity assay and in silico estimated isoelectric point can discriminate fast and slow clearing antibodies, extending previous observations on preclinical clearance. This provides a simple yet powerful approach to select antibodies with desirable PK during early-stage screening.

Characterization and Modeling of Reversible Antibody Self-Association Provide Insights into Behavior, Prediction, and Correction

Reversible antibody self-association, while having major developability and therapeutic implications, is not fully understood or readily predictable and correctable. For a strongly self-associating humanized mAb variant, resulting in unacceptable viscosity, the monovalent affinity of self-interaction was measured in the low μM range, typical of many specific and biologically relevant protein-protein interactions. A face-to-face interaction model extending across both the heavy-chain (HC) and light-chain (LC) Complementary Determining Regions (CDRs) was apparent from biochemical and mutagenesis approaches as well as computational modeling. Light scattering experiments involving individual mAb, Fc, Fab, and Fab'2 domains revealed that Fabs self-interact to form dimers, while bivalent mAb/Fab'2 forms lead to significant oligomerization. Site-directed mutagenesis of aromatic residues identified by homology model patch analysis and self-docking dramatically affected self-association, demonstrating the utility of these predictive approaches, while revealing a highly specific and tunable nature of self-binding modulated by single point mutations. Mutagenesis at these same key HC/LC CDR positions that affect self-interaction also typically abolished target binding with notable exceptions, clearly demonstrating the difficulties yet possibility of correcting self-association through engineering. Clear correlations were also observed between different methods used to assess self-interaction, such as Dynamic Light Scattering (DLS) and Affinity-Capture Self-Interaction Nanoparticle Spectroscopy (AC-SINS). Our findings advance our understanding of therapeutic protein and antibody self-association and offer insights into its prediction, evaluation and corrective mitigation to aid therapeutic development.

Machine Learning Applied to Determine the Molecular Descriptors Responsible for the Viscosity Behavior of Concentrated Therapeutic Antibodies

Predicting the solution viscosity of monoclonal antibody (mAb) drug products remains as one of the main challenges in antibody drug design, manufacturing, and delivery. In this work, the concentration-dependent solution viscosity of 27 FDA-approved mAbs was measured at pH 6.0 in 10 mM histidine-HCl. Six mAbs exhibited high viscosity (>30 cP) in solutions at 150 mg/mL mAb concentration. Combining molecular modeling and machine learning feature selection, we found that the net charge in the mAbs and the amino acid composition in the Fv region are key features which govern the viscosity behavior. For mAbs whose behavior was not dominated by charge effects, we observed that high viscosity is correlated with more hydrophilic and fewer hydrophobic residues in the Fv region. A predictive model based on the net charges of mAbs and a high viscosity index is presented as a fast screening tool for classifying low- and high-viscosity mAbs.

Discovery-stage identification of drug-like antibodies using emerging experimental and computational methods

There is intense and widespread interest in developing monoclonal antibodies as therapeutic agents to treat diverse human disorders. During early-stage antibody discovery, hundreds to thousands of lead candidates are identified, and those that lack optimal physical and chemical properties must be deselected as early as possible to avoid problems later in drug development. It is particularly challenging to characterize such properties for large numbers of candidates with the low antibody quantities, concentrations, and purities that are available at the discovery stage, and to predict concentrated antibody properties (e.g., solubility, viscosity) required for efficient formulation, delivery, and efficacy. Here we review key recent advances in developing and implementing high-throughput methods for identifying antibodies with desirable in vitro and in vivo properties, including favorable antibody stability, specificity, solubility, pharmacokinetics, and immunogenicity profiles, that together encompass overall drug developability. In particular, we highlight impressive recent progress in developing computational methods for improving rational antibody design and prediction of drug-like behaviors that hold great promise for reducing the amount of required experimentation. We also discuss outstanding challenges that will need to be addressed in the future to fully realize the great potential of using such analysis for minimizing development times and improving the success rate of antibody candidates in the clinic.

Toward Drug-Like Multispecific Antibodies by Design

The success of antibody therapeutics is strongly influenced by their multifunctional nature that couples antigen recognition mediated by their variable regions with effector functions and half-life extension mediated by a subset of their constant regions. Nevertheless, the monospecific IgG format is not optimal for many therapeutic applications, and this has led to the design of a vast number of unique multispecific antibody formats that enable targeting of multiple antigens or multiple epitopes on the same antigen. Despite the diversity of these formats, a common challenge in generating multispecific antibodies is that they display suboptimal physical and chemical properties relative to conventional IgGs and are more difficult to develop into therapeutics. Here we review advances in the design and engineering of multispecific antibodies with drug-like properties, including favorable stability, solubility, viscosity, specificity and pharmacokinetic properties. We also highlight emerging experimental and computational methods for improving the next generation of multispecific antibodies, as well as their constituent antibody fragments, with natural IgG-like properties. Finally, we identify several outstanding challenges that need to be addressed to increase the success of multispecific antibodies in the clinic.