Developability Assessment of Engineered Monoclonal Antibody Variants with a Complex Self-Association Behavior Using Complementary Analytical and in Silico Tools

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.

The role of intermolecular interactions on monoclonal antibody filtration through virus removal membranes

The removal of viruses by filtration is a critical unit operation to ensure the overall safety of monoclonal antibody (mAb) products. Many mAbs show very low filtrate flux during virus removal filtration, although there are still significant uncertainties regarding both the mechanisms and antibody properties that determine the filtration behavior. Experiments were performed with three highly purified mAbs through three different commercial virus filters (Viresolve Pro, Viresolve NFP, and Pegasus SV4) with different pore structures and chemistries. The flux decline observed during mAb filtration was largely reversible, even under conditions where the filtrate flux with the mAb was more than 100-fold smaller than the corresponding buffer flux. The extent of flux decline was highly correlated with the hydrodynamic diameter of the mAb as determined by dynamic light scattering (DLS). The mAb with the lowest filtrate flux for all three membranes showed the largest attractive intermolecular interactions and the greatest hydrophobicity, with the latter determined by binding to a butyl resin in an analytical hydrophobic interaction chromatography (HIC) column. These results strongly suggest that the flux behavior is dominated by reversible self-association of the mAbs, providing important insights into the design of more effective virus filtration processes and in the early identification of problematic mAbs/solution conditions.

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.

Automated optimisation of solubility and conformational stability of antibodies and proteins

Biologics, such as antibodies and enzymes, are crucial in research, biotechnology, diagnostics, and therapeutics. Often, biologics with suitable functionality are discovered, but their development is impeded by developability issues. Stability and solubility are key biophysical traits underpinning developability potential, as they determine aggregation, correlate with production yield and poly-specificity, and are essential to access parenteral and oral delivery. While advances for the optimisation of individual traits have been made, the co-optimization of multiple traits remains highly problematic and time-consuming, as mutations that improve one property often negatively impact others. In this work, we introduce a fully automated computational strategy for the simultaneous optimisation of conformational stability and solubility, which we experimentally validate on six antibodies, including two approved therapeutics. Our results on 42 designs demonstrate that the computational procedure is highly effective at improving developability potential, while not affecting antigen-binding. We make the method available as a webserver at www-cohsoftware.ch.cam.ac.uk.

Multidimensional Protein Solubility Optimization with an Ultrahigh-Throughput Microfluidic Platform

Protein-based biologics are highly suitable for drug development as they exhibit low toxicity and high specificity for their targets. However, for therapeutic applications, biologics must often be formulated to elevated concentrations, making insufficient solubility a critical bottleneck in the drug development pipeline. Here, we report an ultrahigh-throughput microfluidic platform for protein solubility screening. In comparison with previous methods, this microfluidic platform can make, incubate, and measure samples in a few minutes, uses just 20 μg of protein (>10-fold improvement), and yields 10,000 data points (1000-fold improvement). This allows quantitative comparison of formulation excipients, such as sodium chloride, polysorbate, histidine, arginine, and sucrose. Additionally, we can measure how solubility is affected by the combinatorial effect of multiple additives, find a suitable pH for the formulation, and measure the impact of mutations on solubility, thus enabling the screening of large libraries. By reducing material and time costs, this approach makes detailed multidimensional solubility optimization experiments possible, streamlining drug development and increasing our understanding of biotherapeutic solubility and the effects of excipients.

Molecular engineering of a cryptic epitope in Spike RBD improves manufacturability and neutralizing breadth against SARS-CoV-2 variants

There is a continued need for sarbecovirus vaccines that can be manufactured and distributed in low- and middle-income countries (LMICs). Subunit protein vaccines are manufactured at large scales at low costs, have less stringent temperature requirements for distribution in LMICs, and several candidates have shown protection against SARS-CoV-2. We previously reported an engineered variant of the SARS-CoV-2 Spike protein receptor binding domain antigen (RBD-L452K-F490W; RBD-J) with enhanced manufacturability and immunogenicity compared to the ancestral RBD. Here, we report a second-generation engineered RBD antigen (RBD-J6) with two additional mutations to a hydrophobic cryptic epitope in the RBD core, S383D and L518D, that further improved expression titers and biophysical stability. RBD-J6 retained binding affinity to human convalescent sera and to all tested neutralizing antibodies except antibodies that target the class IV epitope on the RBD core. K18-hACE2 transgenic mice immunized with three doses of a Beta variant of RBD-J6 displayed on a virus-like particle (VLP) generated neutralizing antibodies (nAb) to nine SARS-CoV-2 variants of concern at similar levels as two doses of Comirnaty. The vaccinated mice were also protected from challenge with Alpha or Beta SARS-CoV-2. This engineered antigen could be useful for modular RBD-based subunit vaccines to enhance manufacturability and global access, or for further development of variant-specific or broadly acting booster vaccines.

Robust production of monovalent bispecific IgG antibodies through novel electrostatic steering mutations at the CH1-Cλ interface

Bispecific antibodies represent an increasingly large fraction of biologics in therapeutic development due to their expanded scope in functional capabilities. Asymmetric monovalent bispecific IgGs (bsIgGs) have the additional advantage of maintaining a native antibody-like structure, which can provide favorable pharmacology and pharmacokinetic profiles. The production of correctly assembled asymmetric monovalent bsIgGs, however, is a complex engineering endeavor due to the propensity for non-cognate heavy and light chains to mis-pair. Previously, we introduced the DuetMab platform as a general solution for the production of bsIgGs, which utilizes an engineered interchain disulfide bond in one of the CH1-CL domains to promote orthogonal chain pairing between heavy and light chains. While highly effective in promoting cognate heavy and light chain pairing, residual chain mispairing could be detected for specific combinations of Fv pairs. Here, we present enhancements to the DuetMab design that improve chain pairing and production through the introduction of novel electrostatic steering mutations at the CH1-CL interface with lambda light chains (CH1-Cλ). These mutations work together with previously established charge-pair mutations at the CH1-CL interface with kappa light chains (CH1-Cκ) and Fab disulfide engineering to promote cognate heavy and light chain pairing and enable the reliable production of bsIgGs. Importantly, these enhanced DuetMabs do not require engineering of the variable domains and are robust when applied to a panel of bsIgGs with diverse Fv sequences. We present a comprehensive biochemical, biophysical, and functional characterization of the resulting DuetMabs to demonstrate compatibility with industrial production benchmarks. Overall, this enhanced DuetMab platform substantially streamlines process development of these disruptive biotherapeutics.

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.

Advanced Formulations/Drug Delivery Systems for Subcutaneous Delivery of Protein-Based Biotherapeutics

Multiple advanced formulations and drug delivery systems (DDSs) have been developed to deliver protein-based biotherapeutics via the subcutaneous (SC) route. These formulations/DDSs include high-concentration solution, co-formulation of two or more proteins, large volume injection, protein cluster/complex, suspension, nanoparticle, microparticle, and hydrogel. These advanced systems provide clinical benefits related to efficacy and safety, but meanwhile, have more complicated formulations and manufacturing processes compared to conventional solution formulations. To develop a fit-for-purpose formulation/DDS for SC delivery, scientists need to consider multiple factors, such as the primary indication, targeted site, immunogenicity, compatibility, biopharmaceutics, patient compliance, etc. Next, they need to develop appropriate formulation (s) and manufacturing processes using the QbD principle and have a control strategy. This paper aims to provide a comprehensive review of advanced formulations/DDSs recently developed for SC delivery of proteins, as well as some knowledge gaps and potential strategies to narrow them through future research.

Development of an In Vitro System To Emulate an In Vivo Subcutaneous Environment: Small Molecule Drug Assessment

A reliable in vitro system can support and guide the development of subcutaneous (SC) drug products. Although several in vitro systems have been developed, they have some limitations, which may hinder them from getting more engaged in SC drug product development. This study sought to develop a novel in vitro system, namely, Emulator of SubCutaneous Absorption and Release (ESCAR), to better emulate the in vivo SC environment and predict the fate of drugs in SC delivery. ESCAR was designed using computer-aided design (CAD) software and fabricated using the three-dimensional (3D) printing technique. ESCAR has a design of two acceptor chambers representing the blood uptake pathway and the lymphatic uptake pathway, respectively, although only the blood uptake pathway was investigated for small molecules in this study. Via conducting a DoE factor screening study using acetaminophen solution, the relationship of the output (drug release from the "SC" chamber to the "blood circulation" chamber) and the input parameters could be modeled using a variety of methods, including polynomial equations, machine learning methods, and Monte Carlo simulation-based methods. The results suggested that the hyaluronic acid (HA) concentration was a critical parameter, whereas the influence of the injection volume and injection position was not substantial. An in vitro-in vivo correlation (IVIVC) study was developed using griseofulvin suspension to explore the feasibility of applying ESCAR in formulation development and bioequivalence studies. The developed LEVEL A IVIVC model demonstrated that the in vivo PK profile could be correlated with the in vitro release profile. Therefore, using this model, for new formulations, only in vitro studies need to be conducted in ESCAR, and in vivo studies might be waived. In conclusion, ESCAR had important implications for research and development and quality control of SC drug products. Future work would be focused on further optimizing ESCAR and expanding its applications via assessing more types of molecules and formulations.

Two physics‐based models for pH ‐dependent calculations of protein solubility

When engineering a protein for its biological function, many physicochemical properties are also optimized throughout the engineering process, and the protein's solubility is among the most important properties to consider. Here, we report two novel computational methods to calculate the pH-dependent protein solubility, and to rank the solubility of mutants. The first is an empirical method developed for fast ranking of the solubility of a large number of mutants of a protein. It takes into account electrostatic solvation energy term calculated using Generalized Born approximation, hydrophobic patches, protein charge, and charge asymmetry, as well as the changes of protein stability upon mutation. This method has been tested on over 100 mutations for 17 globular proteins, as well as on 44 variants of five different antibodies. The prediction rate is over 80%. The antibody tests showed a Pearson correlation coefficient, R, with experimental data from .83 to .91. The second method is based on a novel, completely force-field-based approach using CHARMm program modules to calculate the binding energy of the protein to a part of the crystal lattice, generated from X-ray structure. The method predicted with very high accuracy the solubility of Ribonuclease SA and its 3K and 5K mutants as a function of pH without any parameter adjustments of the existing BIOVIA Discovery Studio binding affinity model. Our methods can be used for rapid screening of large numbers of design candidates based on solubility, and to guide the design of solution conditions for antibody formulation.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

Assessment of Therapeutic Antibody Developability by Combinations of In Vitro and In Silico Methods

Although antibodies have become the fastest-growing class of therapeutics on the market, it is still challenging to develop them for therapeutic applications, which often require these molecules to withstand stresses that are not present in vivo. We define developability as the likelihood of an antibody candidate with suitable functionality to be developed into a manufacturable, stable, safe, and effective drug that can be formulated to high concentrations while retaining a long shelf life. The implementation of reliable developability assessments from the early stages of antibody discovery enables flagging and deselection of potentially problematic candidates, while focussing available resources on the development of the most promising ones. Currently, however, thorough developability assessment requires multiple in vitro assays, which makes it labor intensive and time consuming to implement at early stages. Furthermore, accurate in vitro analysis at the early stage is compromised by the high number of potential candidates that are often prepared at low quantities and purity. Recent improvements in the performance of computational predictors of developability potential are beginning to change this scenario. Many computational methods only require the knowledge of the amino acid sequences and can be used to identify possible developability issues or to rank available candidates according to a range of biophysical properties. Here, we describe how the implementation of in silico tools into antibody discovery pipelines is increasingly offering time- and cost-effective alternatives to in vitro experimental screening, thus streamlining the drug development process. We discuss in particular the biophysical and biochemical properties that underpin developability potential and their trade-offs, review various in vitro assays to measure such properties or parameters that are predictive of developability, and give an overview of the growing number of in silico tools available to predict properties important for antibody development, including the CamSol method developed in our laboratory.

Deciphering the high viscosity of a therapeutic monoclonal antibody in high concentration formulations by microdialysis-hydrogen/deuterium exchange mass spectrometry

High concentration formulations of therapeutic monoclonal antibodies (mAbs) are highly desired for subcutaneous injection. However, high concentration formulations can exhibit unusual molecular behaviors, such as high viscosity or aggregation, that present challenges for manufacturing and administration. To understand the molecular mechanism of the high viscosity exhibited by high concentration protein formulations, we analyzed a human IgG4 (mAb1) at high protein concentrations using sedimentation velocity analytical ultracentrifugation (SV-AUC), X-ray crystallography, hydrogen/deuterium exchange mass spectrometry (HDX-MS), and protein surface patches analysis. Particularly, we developed a microdialysis HDX-MS method to determine intermolecular interactions at different protein concentrations. SV-AUC revealed that mAb1 displayed a propensity for self-association of Fab-Fab, Fab-Fc, and Fc-Fc. mAb1 crystal structure and HDX-MS results demonstrated self-association between complementarity-determining regions (CDRs) and Fc through electrostatic interactions. HDX-MS also indicated Fab-Fab interactions through hydrophobic surface patches constructed by mAb1 CDRs. Our multi-method approach, including fast screening of SV-AUC as well as interface analysis by X-ray crystallography and HDX-MS, helped to elucidate the high viscosity of mAb1 at high concentrations as induced by self-associations of Fab-Fc and Fab-Fab.

Measuring Self-Association of Antibody Lead Candidates with Dynamic Light Scattering

In this method chapter, we provide a brief overview of the key methods available to measure self-association of monoclonal antibodies (mAbs) and explain for which experimental throughputs they are usually applied. We then focus on dynamic light scattering (DLS) and describe experimental details on how to measure the diffusion interaction parameter (k_D) which is occasionally referred to as the gold standard for measuring self-association of proteins. The k_D is a well-established parameter to predict solution viscosity, which is one of the most critical developability parameters of mAbs. Finally, we present a pH and excipient screen that is designed to measure self-association with DLS under conditions that are relevant for bioprocessing and formulation of mAbs. The presented light scattering methods are well suited for lead candidate selections where it is essential to select mAbs with high developability potential for progression toward first human dose.

Improving antibody drug development using bionanotechnology

Monoclonal antibodies are being used to treat a remarkable breadth of human disorders. Nevertheless, there are several key challenges at the earliest stages of antibody drug development that need to be addressed using simple and widely accessible methods, especially related to generating antibodies against membrane proteins and identifying antibody candidates with drug-like biophysical properties (high solubility and low viscosity). Here we highlight key bionanotechnologies for preparing functional and stable membrane proteins in diverse types of lipoparticles that are being used to improve antibody discovery and engineering efforts. We also highlight key bionanotechnologies for high-throughput and ultra-dilute screening of antibody biophysical properties during antibody discovery and optimization that are being used for identifying antibodies with superior combinations of in vitro (formulation) and in vivo (half-life) properties.

An open-source automated PEG precipitation assay to measure the relative solubility of proteins with low material requirement

The solubility of proteins correlates with a variety of their properties, including function, production yield, pharmacokinetics, and formulation at high concentrations. High solubility is therefore a key requirement for the development of protein-based reagents for applications in life sciences, biotechnology, diagnostics, and therapeutics. Accurate solubility measurements, however, remain challenging and resource intensive, which limits their throughput and hence their applicability at the early stages of development pipelines, when long-lists of candidates are typically available in minute amounts. Here, we present an automated method based on the titration of a crowding agent (polyethylene glycol, PEG) to quantitatively assess relative solubility of proteins using about 200 µg of purified material. Our results demonstrate that this method is accurate and economical in material requirement and costs of reagents, which makes it suitable for high-throughput screening. This approach is freely-shared and based on a low cost, open-source liquid-handling robot. We anticipate that this method will facilitate the assessment of the developability of proteins and make it substantially more accessible.

Effect of charged mutation on aggregation of a pentapeptide: Insights from molecular dynamics simulations

Aggregation of therapeutic monoclonal antibodies (mAbs) can negatively affect their chemistry, manufacturing, and control attributes and lead to undesirable immune responses in patients. Therefore, optimization of lead mAb drug candidates during discovery stages to mitigate aggregation is increasingly becoming an integral part of their developability assessments. The disruption of short sequence motifs called aggregation prone regions (APRs) found in amino acid sequences of mAb candidates can potentially mitigate their aggregation. In this work, we have performed molecular dynamics simulations to study the aggregation of an APR (VLVIY) found in λ light chains of human antibodies and its single point mutant KLVIY. Eighteen different multicopy peptide simulation systems of "VLVIY" and "KLVIY" were constructed by varying their concentrations, temperatures, termini capping, and flanking gate-keeper regions. Within 20 ns of the simulation, peptide "VLVIY" formed an aggregate of 100 peptides at ~0.1 M concentration with a 60% reduction in solvent accessible surface area (SASA). Furthermore, analysis of the SASA change, peptide cluster distribution, and water residence time demonstrated how Val ➔ Lys mutation resists aggregation and improves solubility. Presence of Lys slows down aggregation kinetics via charge-charge repulsions and by raising the kinetic barrier to formation of large oligomers. However, the effect of the Val ➔ Lys mutation is dependent on sequence and structural contexts around the APR. This mutation also alters the solvation shell around the peptide by favoring solute-solvent interactions, thereby increasing its solubility. This work has provided a detailed mechanistic explanation of how APR disruption can mitigate aggregation in biotherapeutics and improve their developability.

Rational optimization of a monoclonal antibody improves the aggregation propensity and enhances the CMC properties along the entire pharmaceutical process chain

The discovery of therapeutic monoclonal antibodies (mAbs) primarily focuses on their biological activity favoring the selection of highly potent drug candidates. These candidates, however, may have physical or chemical attributes that lead to unfavorable chemistry, manufacturing, and control (CMC) properties, such as low product titers, conformational and colloidal instabilities, or poor solubility, which can hamper or even prevent development and manufacturing. Hence, there is an urgent need to consider the developability of mAb candidates during lead identification and optimization. This work provides a comprehensive proof of concept study for the significantly improved developability of a mAb variant that was optimized with the help of sophisticated in silico tools relative to its difficult-to-develop parental counterpart. Interestingly, a single amino acid substitution in the variable domain of the light chain resulted in a three-fold increased product titer after stable expression in Chinese hamster ovary cells. Microscopic investigations revealed that wild type mAb-producing cells displayed potential antibody inclusions, while the in silico optimized variant-producing cells showed a rescued phenotype. Notably, the drug substance of the in silico optimized variant contained substantially reduced levels of aggregates and fragments after downstream process purification. Finally, formulation studies unraveled a significantly enhanced colloidal stability of the in silico optimized variant while its folding stability and potency were maintained. This study emphasizes that implementation of bioinformatics early in lead generation and optimization of biotherapeutics reduces failures during subsequent development activities and supports the reduction of project timelines and resources.

Evaluation of in silico tools for the prediction of protein and peptide aggregation on diverse datasets

Several prediction algorithms and tools have been developed in the last two decades to predict protein and peptide aggregation. These in silico tools aid to predict the aggregation propensity and amyloidogenicity as well as the identification of aggregation-prone regions. Despite the immense interest in the field, it is of prime importance to systematically compare these algorithms for their performance. In this review, we have provided a rigorous performance analysis of nine prediction tools using a variety of assessments. The assessments were carried out on several non-redundant datasets ranging from hexapeptides to protein sequences as well as amyloidogenic antibody light chains to soluble protein sequences. Our analysis reveals the robustness of the current prediction tools and the scope for improvement in their predictive performances. Insights gained from this work provide critical guidance to the scientific community on advantages and limitations of different aggregation prediction methods and make informed decisions about their research needs.

Combining Unfolding Reversibility Studies and Molecular Dynamics Simulations to Select Aggregation-Resistant Antibodies

The efficient development of new therapeutic antibodies relies on developability assessment with biophysical and computational methods to find molecules with drug-like properties such as resistance to aggregation. Despite the many novel approaches to select well-behaved proteins, antibody aggregation during storage is still challenging to predict. For this reason, there is a high demand for methods that can identify aggregation-resistant antibodies. Here, we show that three straightforward techniques can select the aggregation-resistant antibodies from a dataset with 13 molecules. The ReFOLD assay provided information about the ability of the antibodies to refold to monomers after unfolding with chemical denaturants. Modulated scanning fluorimetry (MSF) yielded the temperatures that start causing irreversible unfolding of the proteins. Aggregation was the main reason for poor unfolding reversibility in both ReFOLD and MSF experiments. We therefore performed temperature ramps in molecular dynamics (MD) simulations to obtain partially unfolded antibody domains in silico and used CamSol to assess their aggregation potential. We compared the information from ReFOLD, MSF, and MD to size-exclusion chromatography (SEC) data that shows whether the antibodies aggregated during storage at 4, 25, and 40 °C. Contrary to the aggregation-prone molecules, the antibodies that were resistant to aggregation during storage at 40 °C shared three common features: (i) higher tendency to refold to monomers after unfolding with chemical denaturants, (ii) higher onset temperature of nonreversible unfolding, and (iii) unfolding of regions containing aggregation-prone sequences at higher temperatures in MD simulations.

Pragmatic mAb lead molecule engineering from a developability perspective

As the number of antibody drugs being approved and marketed increases, our knowledge of what makes potential drug candidates a successful product has increased tremendously. One of the critical parameters that have become clear in the field is the importance of mAb "developability." Efforts are being increasingly focused on simultaneously selecting molecules that exhibit both desirable biological potencies and manufacturability attributes. In the current study mutations to improve the developability profile of a problematic antibody that inconsistently precipitates in a batch scale-dependent fashion using a standard platform purification process are described. Initial bioinformatic analysis showed the molecule has no obvious sequence or structural liabilities that might lead it to precipitate. Subsequent analysis of the molecule revealed the presence of two unusual positively charged mutations on the light chain at the interface of VH and VL domains, which were hypothesized to be the primary contributor to molecule precipitation during process development. To investigate this hypothesis, straightforward reversion to the germline of these residues was carried out. The resulting mutants have improved expression titers and recovered stability within a forced precipitation assay, without any change to biological activity. Given the time pressures of drug development in industry, process optimization of the lead molecule was carried out in parallel to the "retrospective" mutagenesis approach. Bespoke process optimization for large-scale manufacturing was successful. However, we propose that such context-dependent sequence liabilities should be included in the arsenal of in silico developability screening early in development; particularly since this specific issue can be efficiently mitigated without the requirement for extensive screening of lead molecule variants.

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.

Predicting Antibody Developability Profiles Through Early Stage Discovery Screening

Monoclonal antibodies play an increasingly important role for the development of new drugs across multiple therapy areas. The term 'developability' encompasses the feasibility of molecules to successfully progress from discovery to development via evaluation of their physicochemical properties. These properties include the tendency for self-interaction and aggregation, thermal stability, colloidal stability, and optimization of their properties through sequence engineering. Selection of the best antibody molecule based on biological function, efficacy, safety, and developability allows for a streamlined and successful CMC phase. An efficient and practical high-throughput developability workflow (100 s-1,000 s of molecules) implemented during early antibody generation and screening is crucial to select the best lead candidates. This involves careful assessment of critical developability parameters, combined with binding affinity and biological properties evaluation using small amounts of purified material (<1 mg), as well as an efficient data management and database system. Herein, a panel of 152 various human or humanized monoclonal antibodies was analyzed in biophysical property assays. Correlations between assays for different sets of properties were established. We demonstrated in two case studies that physicochemical properties and key assay endpoints correlate with key downstream process parameters. The workflow allows the elimination of antibodies with suboptimal properties and a rank ordering of molecules for further evaluation early in the candidate selection process. This enables any further engineering for problematic sequence attributes without affecting program timelines.

Dissecting the molecular basis of high viscosity of monospecific and bispecific IgG antibodies

Some antibodies exhibit elevated viscosity at high concentrations, making them poorly suited for therapeutic applications requiring administration by injection such as subcutaneous or ocular delivery. Here we studied an anti-IL-13/IL-17 bispecific IgG₄ antibody, which has anomalously high viscosity compared to its parent monospecific antibodies. The viscosity of the bispecific IgG₄ in solution was decreased by only ~30% in the presence of NaCl, suggesting electrostatic interactions are insufficient to fully explain the drivers of viscosity. Intriguingly, addition of arginine-HCl reduced the viscosity of the bispecific IgG₄ by ~50% to its parent IgG level. These data suggest that beyond electrostatics, additional types of interactions such as cation-π and/or π-π may contribute to high viscosity more significantly than previously understood. Molecular dynamics simulations of antibody fragments in the mixed solution of free arginine and explicit water were conducted to identify hotspots involved in self-interactions. Exposed surface aromatic amino acids displayed an increased number of contacts with arginine. Mutagenesis of the majority of aromatic residues pinpointed by molecular dynamics simulations effectively decreased the solution's viscosity when tested experimentally. This mutational method to reduce the viscosity of a bispecific antibody was extended to a monospecific anti-GCGR IgG₁ antibody with elevated viscosity. In all cases, point mutants were readily identified that both reduced viscosity and retained antigen-binding affinity. These studies demonstrate a new approach to mitigate high viscosity of some antibodies by mutagenesis of surface-exposed aromatic residues on complementarity-determining regions that may facilitate some clinical applications.

Highly sensitive detection of antibody nonspecific interactions using flow cytometry

The rapidly evolving nature of antibody drug development has resulted in technologies that generate vast numbers (hundreds to thousands) of lead antibody candidates during early discovery. These candidates must be rapidly pared down to identify the most drug-like candidates for in-depth analysis of their safety and efficacy, which can only be performed on a limited number of antibodies due to time and resource requirements. One key biophysical property of successful antibody therapeutics is high specificity, defined as low levels of nonspecific binding or polyspecificity. Although there has been some progress in developing assays for detecting antibody polyspecificity, most of these assays are limited by poor sensitivity or assay formats that require proprietary antibody surface display methods, and some of these assays use complex and poorly defined polyspecificity reagents. Here we report the PolySpecificity Particle (PSP) assay, a sensitive flow cytometry assay for evaluating antibody nonspecific interactions that overcomes previous limitations and can be used for evaluating diverse types of IgGs, multispecific antibodies and Fc-fusion proteins. Our approach uses micron-sized magnetic beads coated with Protein A to capture antibodies at extremely dilute concentrations (<0.02 mg/mL). Flow cytometry analysis of polyspecificity reagent binding to these conjugates results in sensitive detection of differences in nonspecific interactions for clinical-stage antibodies. Our PSP assay strongly discriminates between antibodies with different levels of polyspecificity using previously reported polyspecificity reagents that are either well-defined proteins or highly complex protein mixtures. Moreover, we also find that a unique reagent, namely ovalbumin, results in the best assay sensitivity and specificity. Importantly, our assay is much more sensitive than standard assays such as ELISAs. We expect that our simple, sensitive, and high-throughput PSP assay will accelerate the development of safe and effective antibody therapeutics.

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.

Stability liabilities of biotherapeutic proteins: Early assessment as mitigation strategy

Identification of molecular liabilities and implementation of mitigation strategies are key aspects that need to be considered by pharmaceutical companies developing therapeutic proteins. In the field of monoclonal antibodies, an efficient and streamlined process known as developability assessment is used for the selection of the "fittest" candidate. Other protein modalities, have in most cases only a limited number of possible candidates, requiring a paradigm change to a concept of candidate enabling. The assessment of liabilities at early project phases with the possibility to re-engineer candidates becomes essential for the success of these projects. Each protein possesses a unique stability profile resulting from the interplay of conformational, colloidal, chemical and physical stability attributes. All of these attributes strongly depend on external factors. Conformational and colloidal stability profiles of three non-immunoglobulin domain based proteins, namely Carbonic anhydrase, Ovalbumin and Thyroglobulin, and of two monoclonal antibodies were assessed in dependence of solution pH, ionic strength and varying buffering agents. The impact of screened external factors on proteins' stability attributes varied significantly, indicating presence of molecule specific liabilities. Screening of such a broad space of conditions at early project phases is only feasible using low-material consuming, high-throughput analytical methods as exemplified in this study.

Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Charge and hydrophobicity are key features in sequence-trained machine learning models for predicting the biophysical properties of clinical-stage antibodies

Improved understanding of properties that mediate protein solubility and resistance to aggregation are important for developing biopharmaceuticals, and more generally in biotechnology and synthetic biology. Recent acquisition of large datasets for antibody biophysical properties enables the search for predictive models. In this report, machine learning methods are used to derive models for 12 biophysical properties. A physicochemical perspective is maintained in analysing the models, leading to the observation that models cluster largely according to charge (cross-interaction measurements) and hydrophobicity (self-interaction methods). These two properties also overlap in some cases, for example in a new interpretation of variation in hydrophobic interaction chromatography. Since the models are developed from differences of antibody variable loops, the next stage is to extend models to more diverse protein sets.

AVAILABILITY

The web application for the sequence-based algorithms are available on the protein-sol webserver, at https://protein-sol.manchester.ac.uk/abpred, with models and virtualisation software available at https://protein-sol.manchester.ac.uk/software.

Identifying Key Residues That Drive Strong Electrostatic Attractions between Therapeutic Antibodies

Attractive electrostatic protein-protein interactions (PPI) necessarily involve identifying oppositely charged regions of the protein surface that interact favorably. This cannot be done reliably if one only considers a single protein in isolation unless there are obvious charge "patches" that result in extreme molecular dipoles. Prior work [ J. Pharm. Sci. 2019 , 108 , 120 - 132 ] identified three monoclonal antibodies (MAbs) that displayed experimental behavior ranging from net repulsive to strongly attractive electrostatic interactions. The present work provides a systematic computational approach for identifying the origin of diverse PPI, in terms of which sets of amino acids or individual amino acids are most influential, and determining if there are different patterns of pairwise amino acid interaction "maps" that result in different behaviors. The charge was eliminated computationally, one by one, for each charged residue in the wild-type sequences, which resulted in predicted changes in the second osmotic virial coefficient. The results highlight interaction "maps" that correspond to cases with qualitatively different net electrostatic PPI for the different MAbs and solution conditions, as well as key sets of residues that contribute to strongly attractive PPI. A more computationally efficient method is also proposed to identify key amino acids based on Mayer-weighted interaction energies.

Protein Solubility Predictions Using the CamSol Method in the Study of Protein Homeostasis

One of the major functions of the protein homeostasis system is to maintain proteins in their soluble states, and indeed several human disorders are associated with the aberrant aggregation of proteins. An active involvement of the protein homeostasis system is necessary to avoid aggregation because proteins are expressed at levels close to their solubility limits, hence being poorly soluble. The mechanisms by which the protein homeostasis system acts to control protein aggregation are, however, still not known in much detail. To facilitate systematic investigations of these mechanisms, we describe here the CamSol method of predicting protein solubility, and illustrate its initial applications. We anticipate that with the advent of powerful proteomics and transcriptomic methods, in combination with the use of the CamSol method and related approaches to predict the solubility and other biophysical properties of proteins, it will become possible to increase our understanding of the principles of protein homeostasis related to the maintenance of the proteome in its soluble form.

Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies

The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.

Computer aided analysis of disease linked protein networks

Proteins can interact in various ways, ranging from direct physical relationships to indirect interactions in a formation of protein-protein interaction network. Diagnosis of the protein connections is critical to identify various cellular pathways. Today constructing and analyzing the protein interaction network is being developed as a powerful approach to create network pharmacology toward detecting unknown genes and proteins associated with diseases. Discovery drug targets regarding therapeutic decisions are exciting outcomes of studying disease networks. Protein connections may be identified by experimental and recent new computational approaches. Due to difficulties in analyzing in-vivo proteins interactions, many researchers have encouraged improving computational methods to design protein interaction network. In this review, the experimental and computational approaches and also advantages and disadvantages of these methods regarding the identification of new interactions in a molecular mechanism have been reviewed. Systematic analysis of complex biological systems including network pharmacology and disease network has also been discussed in this review.

Binding symmetry and surface flexibility mediate antibody self-association

Solution stability is an important factor in the optimization of engineered biotherapeutic candidates such as monoclonal antibodies because of its possible effects on manufacturability, pharmacology, efficacy and safety. A detailed atomic understanding of the mechanisms governing self-association of natively folded protein monomers is required to devise predictive tools to guide screening and re-engineering along the drug development pipeline. We investigated pairs of affinity-matured full-size antibodies and observed drastically different propensities to aggregate from variants differing by a single amino-acid. Biophysical testing showed that antigen-binding fragments (Fabs) from the aggregating antibodies also reversibly associated with equilibrium dissociation constants in the low-micromolar range. Crystal structures (PDB accession codes 6MXR, 6MXS, 6MY4, 6MY5) and bottom-up hydrogen-exchange mass spectrometry revealed that Fab self-association occurs in a symmetric mode that involves the antigen complementarity-determining regions. Subtle local conformational changes incurred upon point mutation of monomeric variants foster formation of complementary polar interactions and hydrophobic contacts to generate a dimeric Fab interface. Testing of popular in silico tools generally indicated low reliabilities for predicting the aggregation propensities observed. A structure-aggregation data set is provided here in order to stimulate further improvements of in silico tools for prediction of native aggregation. Incorporation of intermolecular docking, conformational flexibility, and short-range packing interactions may all be necessary features of the ideal algorithm.

Science and art of protein formulation development

Protein pharmaceuticals have become a significant class of marketed drug products and are expected to grow steadily over the next decade. Development of a commercial protein product is, however, a rather complex process. A critical step in this process is formulation development, enabling the final product configuration. A number of challenges still exist in the formulation development process. This review is intended to discuss these challenges, to illustrate the basic formulation development processes, and to compare the options and strategies in practical formulation development.