Factors affecting the viscosity in high concentration solutions of different monoclonal antibodies

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.

Landscape of Subcutaneous Administration Strategies for Monoclonal Antibodies in Oncology

In recent decades, subcutaneous (SC) administration of monoclonal antibodies (mAbs) has emerged as a promising alternative to intravenous delivery in oncology, offering comparable therapeutic efficacy while addressing patient preferences. This perspective article provides an in-depth analysis of the technological landscape surrounding SC mAb administration in oncology. It outlines various technologies under evaluation across developmental stages, spanning from preclinical investigations to the integration of established methodologies in clinical practice. Additionally, this perspective article explores emerging trends and prospective trajectories, shedding light on the evolving landscape of SC mAb administration. Furthermore, it emphasizes key checkpoints related to quality attributes essential for optimizing mAb delivery via the SC route. This review serves as a valuable resource for researchers, clinicians, and healthcare policymakers, offering insights into the advancement of SC mAb administration in oncology and its implications for patient care.

Flow Activation Energy of High-Concentration Monoclonal Antibody Solutions and Protein-Protein Interactions Influenced by NaCl and Sucrose

The solution viscosity and protein-protein interactions (PPIs) as a function of temperature (4-40 °C) were measured at a series of protein concentrations for a monoclonal antibody (mAb) with different formulation conditions, which include NaCl and sucrose. The flow activation energy (Eη) was extracted from the temperature dependence of solution viscosity using the Arrhenius equation. PPIs were quantified via the protein diffusion interaction parameter (kD) measured by dynamic light scattering, together with the osmotic second virial coefficient and the structure factor obtained through small-angle X-ray scattering. Both viscosity and PPIs were found to vary with the formulation conditions. Adding NaCl introduces an attractive interaction but leads to a significant reduction in the viscosity. However, adding sucrose enhances an overall repulsive effect and leads to a slight decrease in viscosity. Thus, the averaged (attractive or repulsive) PPI information is not a good indicator of viscosity at high protein concentrations for the mAb studied here. Instead, a correlation based on the temperature dependence of viscosity (i.e., Eη) and the temperature sensitivity in PPIs was observed for this specific mAb. When kD is more sensitive to the temperature variation, it corresponds to a larger value of Eη and thus a higher viscosity in concentrated protein solutions. When kD is less sensitive to temperature change, it corresponds to a smaller value of Eη and thus a lower viscosity at high protein concentrations. Rather than the absolute value of PPIs at a given temperature, our results show that the temperature sensitivity of PPIs may be a more useful metric for predicting issues with high viscosity of concentrated solutions. In addition, we also demonstrate that caution is required in choosing a proper protein concentration range to extract kD. In some excipient conditions studied here, the appropriate protein concentration range needs to be less than 4 mg/mL, remarkably lower than the typical concentration range used in the literature.

Elucidation of the Reversible Self-Association Interface of a Diabody-Interleukin Fusion Protein Using Hydrogen-Exchange Mass Spectrometry and In Silico Modeling

Reversible self-association (RSA) of therapeutic proteins presents major challenges in the development of high-concentration formulations, especially those intended for subcutaneous administration. Understanding self-association mechanisms is therefore critical to the design and selection of candidates with acceptable developability to advance to clinical trials. The combination of experiments and in silico modeling presents a powerful tool to elucidate the interface of self-association. RSA of monoclonal antibodies has been studied extensively under different solution conditions and have been shown to involve interactions for both the antigen-binding fragment and the crystallizable fragment. Novel modalities such as bispecific antibodies, antigen-binding fragments, single-chain-variable fragments, and diabodies constitute a fast-growing class of antibody-based therapeutics that have unique physiochemical properties compared to monoclonal antibodies. In this study, the RSA interface of a diabody-interleukin 22 fusion protein (FP-1) was studied using hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS) in combination with in silico modeling. Taken together, the results show that a complex solution behavior underlies the self-association of FP-1 and that the interface thereof can be attributed to a specific segment in the variable light chain of the diabody. These findings also demonstrate that the combination of HDX-MS with in silico modeling is a powerful tool to guide the design and candidate selection of novel biotherapeutic modalities.

A First Insight into the Developability of an Immunoglobulin G3: A Combined Computational and Experimental Approach

Immunoglobulin G 3 (IgG3) monoclonal antibodies (mAbs) are high-value scaffolds for developing novel therapies. Despite their wide-ranging therapeutic potential, IgG3 physicochemical properties and developability characteristics remain largely under-characterized. Protein-protein interactions elevate solution viscosity in high-concentration formulations, impacting physicochemical stability, manufacturability, and the injectability of mAbs. Therefore, in this manuscript, the key molecular descriptors and biophysical properties of a model anti-IL-8 IgG1 and its IgG3 ortholog are characterized. A computational and experimental framework was applied to measure molecular descriptors impacting their downstream developability. Findings from this approach underpin a detailed understanding of the molecular characteristics of IgG3 mAbs as potential therapeutic entities. This work is the first report examining the manufacturability of IgG3 for high-concentration mAb formulations. While poorer conformational and colloidal stability and elevated solution viscosity were observed for IgG3, future efforts controlling surface potential through sequence-engineering of solvent-accessible patches can be used to improve biophysical parameters that dictate mAb developability.

Small-Angle X-ray Scattering as a Powerful Tool for Phase and Crystallinity Assessment of Monoclonal Antibody Crystallites in Support of Batch Crystallization

Crystalline suspensions of monoclonal antibodies (mAbs) have great potential to improve drug substance isolation and purification on a large scale and to be used for drug delivery via high-concentration formulations. Crystalline mAb suspensions are expected to have enhanced chemical and physical properties relative to mAb solutions delivered intravenously, making them attractive candidates for subcutaneous delivery. In contrast to small molecules, the development of protein crystalline suspensions is not a widely used approach in the pharmaceutical industry. This is mainly due to the challenges in finding crystalline hits and the suboptimal physical properties of the resulting crystallites when hits are found. Modern advances in instrumentation and increased knowledge of mAb crystallization have, however, resulted in higher probabilities of discovering crystal forms and improving their particle properties and characterization. In this regard, physical, analytical characterization plays a central role in the initial steps of understanding and later optimizing the crystallization of mAbs and requires careful selection of the appropriate tools. This contribution describes a novel crystal structure of the antibody pembrolizumab and demonstrates the usefulness of small-angle X-ray scattering (SAXS) for characterizing its crystalline suspensions. It illustrates the advantages of SAXS when used to (i) confirm crystallinity and crystal phase of crystallites produced in batch mode; (ii) confirm crystallinity under various conditions and detect variations in crystal phases, enabling fine-tuning of the crystallizations for phase control across multiple batches; (iii) monitor the physical response and stability of the crystallites in suspension with regard to filtration and washing; and (iv) monitor the physical stability of the crystallites upon drying. Overall, this work highlights how SAXS is an essential tool for mAb crystallization characterization.

Concentration-Dependent Diffusion of Monoclonal Antibodies: Underlying Mechanisms of Anomalous Diffusion

The development, storage, transport, and subcutaneous delivery of highly concentrated monoclonal antibody formulations pose significant challenges due to the high solution viscosity and low diffusion of the antibody molecules in crowded environments. These issues often stem from the self-associating behavior of the antibody molecules, potentially leading to aggregation. In this work, we used a dissipative particle dynamics-based coarse-grained model to investigate the diffusion behavior of IgG1 antibody molecules in aqueous solutions with 15 and 32 mM NaCl and antibody concentrations ranging from 10 to 400 mg/mL. We determined the coarse-grained interaction parameters by matching the calculated structure factor with the computational and experimental data from the literature. Our results indicate Fickian diffusion for antibody concentrations of 10 and 25 mg/mL and anomalous diffusion for concentrations exceeding 50 mg/mL. The anomalous diffusion was observed for ∼0.33 to 0.4 μs, followed by Fickian diffusion for all antibody concentrations. We observed a strong linear correlation between the diffusion behavior of the antibody molecules (diffusion coefficient D and anomalous diffusion exponent α) and the amount of aggregates present in the solution and between the amount of aggregates and the Coulomb interaction energy. The investigation of underlying mechanisms for anomalous diffusion revealed that in crowded environments at high antibody concentrations, the attractive interaction between electrostatically complementary regions of the antibody molecules could further bring the neighboring molecules closer to one another, ultimately resulting in aggregate formation. Further, the Coulomb attraction can continue to draw more molecules together, forming larger aggregates.

The search for novel proline analogs for viscosity reduction and stabilization of highly concentrated monoclonal antibody solutions

Administration of monoclonal antibodies (mAbs) is currently focused on subcutaneous injection associated with increased patient adherence and reduced treatment cost, leading to sustainable healthcare. The main bottleneck is low volume that can be injected, requiring highly concentrated mAb solutions. The latter results in increased solution viscosity with pronounced mAb aggregation propensity because of intensive protein-protein interactions. Small molecule excipients have been proposed to restrict the protein-protein interactions, contributing to reduced viscosity. The aim of the study was to discover novel compounds that reduce the viscosity of highly concentrated mAb solution. First, the chemical space of proline analogs was explored and 35 compounds were determined. Viscosity measurements revealed that 18 proline analogs reduced the mAb solution viscosity similar to or more than proline. The compounds forming both electrostatic and hydrophobic interactions with mAb reduced the viscosity of the formulation more efficiently without detrimentally effecting mAb physical stability. A correlation between the level of interaction and viscosity-reducing effect was confirmed with molecular dynamic simulations. Structure rigidity of the compounds and aromaticity contributed to their viscosity-reducing effect, dependent on molecule size. The study results highlight the novel proline analogs as an effective approach in viscosity reduction in development of biopharmaceuticals for subcutaneous administration.

Quantifying Protein Shape to Elucidate Its Influence on Solution Viscosity in High-Concentration Electrolyte Solutions

Therapeutic proteins with a high concentration and low viscosity are highly desirable for subcutaneous and certain local injections. The shape of a protein is known to influence solution viscosity; however, the precise quantification of protein shape and its relative impact compared to other factors like charge-charge interactions remains unclear. In this study, we utilized seven model proteins of varying shapes and experimentally determined their shape factors (v) based on Einstein's viscosity theory, which correlate strongly with the ratios of the proteins' surface area to the 2/3 power of their respective volumes, based on protein crystal structures resolved experimentally or predicted by AlphaFold. This finding confirms the feasibility of computationally estimating protein shape factors from amino acid sequences alone. Furthermore, our results demonstrated that, in high-concentration electrolyte solutions, a more spherical protein shape increases the protein's critical concentration (C*), the transition concentration beyond which protein viscosity increases exponentially relative to concentration increases. In summary, our work elucidates protein shape as a key determinant of solution viscosity through quantitative analysis and comparison with other contributing factors. This provides insights into molecular engineering strategies to optimize the molecular design of therapeutic proteins, thus optimizing their viscosity.

Antibody association in solution: cluster distributions and mechanisms

Understanding factors that affect the clustering and association of antibodies molecules in solution is critical to their development as therapeutics. For 19 different monoclonal antibody (mAb) solutions, we measured the viscosities, the second virial coefficients, the Kirkwood-Buff integrals, and the cluster distributions of the antibody molecules as functions of protein concentration. Solutions were modeled using the statistical-physics Wertheim liquid-solution theory, representing antibodies as Y-shaped molecular structures of seven beads each. We found that high-viscosity solutions result from more antibody molecules per cluster. Multi-body properties such as viscosity are well predicted experimentally by the 2-body Kirkwood-Buff quantity, G22, but not by the second virial coefficient, B22, and well-predicted theoretically from the Wertheim protein-protein sticking energy. Weakly interacting antibodies are rate-limited by nucleation; strongly interacting ones by propagation. This approach gives a way to relate micro to macro properties of solutions of associating proteins.

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.

Gel characteristics of low-acetyl spruce galactoglucomannans

Galactoglucomannans (GGM) recovered from abundant forest industry side-streams has been widely recognized as a renewable hydrocolloid. The low molar mass and presence of O-acetyl side-groups results in low viscous dispersions and weak intermolecular interactions that make GGM unsuitable for hydrogel formation, unless forcefully chemically derivatized and/or crosslinked with other polymers. Here we present the characterization of hydrogels prepared from GGM after tailoring the degree of acetylation by alkaline treatment during its recovery. Specifically, we investigated gel characteristics of low-acetyl GGM dispersions prepared at varied solid concentrations (5, 10 and 15 %) and pH (4, 7 and 10), and then subjected to ultrasonication. The results indicated that low-acetyl GGM dispersions formed gels (G' > G″) at all other studied solid concentration and pH level combinations except 5 % and pH 4. High pH levels, leading to further removal of acetyl groups, and high solid concentration facilitated the gel formation. GGM hydrogels were weak gels with strong shear-thinning behavior and thixotropic properties, and high hardness and water holding capacity; which were enhanced with increased pH and solid concentration, and prolonged storage time. Our study showed the possibility to utilize low-acetyl GGM as mildly processed gelling or thickening agents, and renewable materials for bio-based hydrogels.

Predicting Colloidal Stability of High-Concentration Monoclonal Antibody Formulations in Common Pharmaceutical Buffers Using Improved Polyethylene Glycol Induced Protein Precipitation Assay

Colloidal stability is an important consideration when developing high concentration mAb formulations. PEG-induced protein precipitation is a commonly used assay to assess the colloidal stability of protein solutions. However, the practical usefulness and the current theoretical model for this assay have yet to be verified over a large formulation space across multiple mAbs and mAb-based modalities. In the present study, we used PEG-induced protein precipitation assays to evaluate colloidal stability of 3 mAbs in 24 common formulation buffers at 20 and 5 °C. These prediction assays were conducted at low protein concentration (1 mg/mL). We also directly characterized high concentration (100 mg/mL) formulations for cold-induced phase separation, turbidity, and concentratibility by ultrafiltration. This systematic study allowed analysis of the correlation between the results of low concentration assays and the high concentration attributes. The key findings of this study include the following: (1) verification of the usefulness of three different parameters (Cmid, μB, and Tcloud) from PEG-induced protein precipitation assays for ranking colloidal stability of high concentration mAb formulations; (2) a new method to implement PEG-induced protein precipitation assay suitable for high throughput screening with low sample consumption; (3) improvement in the theoretical model for calculating robust thermodynamic parameters of colloidal stability (μB and εB) that are independent of specific experimental settings; (4) systematic evaluation of the effects of pH and buffer salts on colloidal stability of mAbs in common formulation buffers. These findings provide improved theoretical and practical tools for assessing the colloidal stability of mAbs and mAb-based modalities during formulation development.

Application of Formulation Principles to Stability Issues Encountered During Processing, Manufacturing, and Storage of Drug Substance and Drug Product Protein Therapeutics

The field of formulation and stabilization of protein therapeutics has become rather extensive. However, most of the focus has been on stabilization of the final drug product. Yet, proteins experience stress and degradation through the manufacturing process, starting with fermentaition. This review describes how formulation principles can be applied to stabilize biopharmaceutical proteins during bioprocessing and manufacturing, considering each unit operation involved in prepration of the drug substance. In addition, the impact of the container on stabilty is discussed as well.

The Evolution of Commercial Antibody Formulations

It has been nearly four decades since the first therapeutic monoclonal antibodies were approved and made available for widespread human use. Herein, US and EU approved antibody formulations are reviewed, and their nature and compositions are evaluated over time. From 1986 through Jan 2023, significant formulation trends have occurred and to represent this, 165 commercial antibody therapeutic formulations were binned into 5 different periods of time. Overall, we have observed the following: 1) The average formulation pH has decreased in recent years by over 0.5 units along with a decrease in variability that is largely driven by non-high concentration liquid in vial presentations for IV administration, 2) The use of certain excipients and buffers such as histidine, sucrose, metal chelators, arginine and methionine has become significantly more common, whereas formulations that contain phosphate, salt, no sugar or no surfactant have fallen out of favor, 3) Overall formulation space has increasingly become more homogenous and has converged in terms of formulation pH and excipient preferences regardless of formulation concentration, drug product presentation, and route of administration, 4) The average calculated isoelectric point (pI) has decreased 0.26 units, and 5) Overall, the average formulation pH and calculated pI for all commercial antibodies surveyed was 6.0 and 8.4, respectively. These trends and formulation convergence may be driven by multiple factors such as advancements in high-throughput computational and analytical technologies, the increased emphasis and understanding of certain developability attributes and formulation principles during lead selection and formulation development, and the adoption of low-risk development platform approaches.

Predictive Nature of High-Throughput Assays in ADC Formulation Screening

Utilization of high-throughput biophysical screening techniques during early screening studies is warranted due to the limited amount of material and large number of samples. But the predictability of the data to longer-term storage stability is critical as the high-throughput methods assist in defining the design space for the longer-term studies. In this study, the biophysical properties of two ADCs in 16 formulation conditions were evaluated using high-throughput techniques. Conformational stability and colloidal stability were evaluated by determining Tm values, kD, B22, and Tagg. In addition, the samples were placed on stability and the extent of aggregate formation over the 8-week interval was determined. The rank order of the 16 different formulations in the high-throughput assays was compared to the rank order observed during the stability studies to assess the predictive capabilities of the screening methods. It was demonstrated that similar rank orders can be expected between high-throughput physical stability indicating assays such as Tagg and B22 and traditional aggregation by SEC data, whereas conformational stability read-outs (Tm) are less predictive. In addition, the high-throughput assays appropriately identified the poor performing formulation conditions, which is ultimately what is desired of screening assays.

An Intercompany Perspective on Practical Experiences of Predicting, Optimizing and Analyzing High Concentration Biologic Therapeutic Formulations

Developing high-dose biologic drugs for subcutaneous injection often requires high-concentration formulations and optimizing viscosity, solubility, and stability while overcoming analytical, manufacturing, and administration challenges. To understand industry approaches for developing high-concentration formulations, the Formulation Workstream of the BioPhorum Development Group, an industry-wide consortium, conducted an inter-company collaborative exercise which included several surveys. This collaboration provided an industry perspective, experience, and insight into the practicalities for developing high-concentration biologics. To understand solubility and viscosity, companies desire predictive tools, but experience indicates that these are not reliable and experimental strategies are best. Similarly, most companies prefer accelerated and stress stability studies to in-silico or biophysical-based prediction methods to assess aggregation. In addition, optimization of primary container-closure and devices are pursued to mitigate challenges associated with high viscosity of the formulation. Formulation strategies including excipient selection and application of studies at low concentration to high-concentration formulations are reported. Finally, analytical approaches to high concentration formulations are presented. The survey suggests that although prediction of viscosity, solubility, and long-term stability is desirable, the outcome can be inconsistent and molecule dependent. Significant experimental studies are required to confirm robust product definition as modeling at low protein concentrations will not necessarily extrapolate to high concentration formulations.

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.

A systematic review of commercial high concentration antibody drug products approved in the US: formulation composition, dosage form design and primary packaging considerations

Three critical aspects that define high concentration antibody products (HCAPs) are as follows: 1) formulation composition, 2) dosage form, and 3) primary packaging configuration. HCAPs have become successful in the therapeutic sector due to their unique advantage of allowing subcutaneous self-administration. Technical challenges, such as physical and chemical instability, viscosity, delivery volume limitations, and product immunogenicity, can hinder successful development and commercialization of HCAPs. Such challenges can be overcome by robust formulation and process development strategies, as well as rational selection of excipients and packaging components. We compiled and analyzed data from US Food and Drug Administration-approved and marketed HCAPs that are ≥100 mg/mL to identify trends in formulation composition and quality target product profile. This review presents our findings and discusses novel formulation and processing technologies that enable the development of improved HCAPs at ≥200 mg/mL. The observed trends can be used as a guide for further advancements in the development of HCAPs as more complex antibody-based modalities enter biologics product development.

Blueprint for antibody biologics developability

Large-molecule antibody biologics have revolutionized medicine owing to their superior target specificity, pharmacokinetic and pharmacodynamic properties, safety and toxicity profiles, and amenability to versatile engineering. In this review, we focus on preclinical antibody developability, including its definition, scope, and key activities from hit to lead optimization and selection. This includes generation, computational and in silico approaches, molecular engineering, production, analytical and biophysical characterization, stability and forced degradation studies, and process and formulation assessments. More recently, it is apparent these activities not only affect lead selection and manufacturability, but ultimately correlate with clinical progression and success. Emerging developability workflows and strategies are explored as part of a blueprint for developability success that includes an overview of the four major molecular properties that affect all developability outcomes: 1) conformational, 2) chemical, 3) colloidal, and 4) other interactions. We also examine risk assessment and mitigation strategies that increase the likelihood of success for moving the right candidate into the clinic.

Liquid crystal phase formation and non-Newtonian behavior of oligonucleotide formulations

Viscosity behavior of liquid oligonucleotide therapeutics and its dependence on formulation properties has been poorly studied to date. We observed a high increase in viscosity and solidification of therapeutic oligonucleotide formulations with increasing oligonucleotide concentration creating challenges during drug product manufacturing. In this study, we characterized the viscosity behavior of three different single strand DNA oligonucleotides based on oligonucleotide concentration and formulation composition. We subsequently studied the underlying mechanism for increased viscosity at higher oligonucleotide concentrations by dynamic light scattering (DLS), 1H nuclear magnetic resonance (NMR), differential scanning calorimetry (DSC), and polarized light microscopy. Viscosity was highly dependent on formulation composition, oligonucleotide sequence, and concentration, and especially dependent on the presence and combination of different individual ions, such as the presence of sodium chloride in the formulation. In samples with elevated viscosity, the viscosity behavior was characterized by non-Newtonian, shear-thinning flow behavior. We further studied these samples by DLS and 1H NMR, which revealed the presence of supra-molecular assemblies, and further characterization by polarized light and DSC characterized these assemblies as liquid crystals in the formulation. The present study links the macroscopic viscosity behavior of oligonucleotide formulations to the formation of supra-molecular assemblies and to the presence of liquid crystals, and highlights the importance of formulation composition selection for these therapeutics.

Technology development to evaluate the effectiveness of viscosity reducing excipients

Addition of pharmaceutical excipients is a commonly used approach to decrease the viscosity of highly concentrated protein formulations, which otherwise could not be subcutaneously injected or processed. The variety of protein-protein interactions, which are responsible for increased viscosities, makes a portfolio approach necessary. Screening of several excipients to develop such a portfolio is time and money consuming in industrial settings. Responsible protein-protein interactions were investigated using the interaction parameter k_D obtained from dynamic light scattering measurements in the studies presented herein. Together with in-silico calculated excipient parameter, k_D could be used as a screening tool accelerating screening and formulation development as k_D is suitable to high-throughput formats using small quantities of protein and low concentrations. A qualitative correlation between k_D and high-concentration viscosity behavior could be shown in our case.

Predicting Experimental B22 Values and the Effects of Histidine Charge States for Monoclonal Antibodies Using Coarse-Grained Molecular Simulations

Static light scattering (SLS) was used to characterize five monoclonal antibodies (MAbs) as a function of total ionic strength (TIS) at pH values between 5.5 and 7.0. Second osmotic virial coefficient (B₂₂) values were determined experimentally for each MAb as a function of TIS using low protein concentration SLS data. Coarse-grained molecular simulations were performed to predict the B₂₂ values for each MAb at a given pH and TIS. To include the effect of charge fluctuations of titratable residues in the B₂₂ calculations, a statistical approach was introduced in the Monte Carlo algorithm based on the protonation probability based on a given pH value and the Henderson-Hasselbalch equation. The charged residues were allowed to fluctuate individually, based on the sampled microstates and the influence of electrostatic interactions on net protein-protein interactions during the simulations. Compared to static charge simulations, the new approach provided improved results compared to experimental B₂₂ values at pH conditions near the pK_a of titratable residues.

Application of a Simple Short-Range Attraction and Long-Range Repulsion Colloidal Model toward Predicting the Viscosity of Protein Solutions

Some hard-sphere colloidal models have been criticized for inaccurately predicting the solution viscosity of complex biological molecules like proteins. Competing short-range attractions and long-range repulsions, also known as short-range attraction and long-range repulsion (SALR) interactions, have been thought to affect the microstructure of a protein solution at low to moderate ionic strength. However, such interactions have been implicated primarily in causing phase transition, protein gelation, or reversible cluster formation, and their effect on protein solution viscosity change is not fully understood. In this work, we show the application of a hard-sphere colloidal model with SALR interactions toward predicting the viscosity of dilute to semi-dilute protein solutions. The comparison is performed for a globular-shaped albumin and Y-shaped therapeutic monoclonal antibody that are not explained by previous colloidal models. The model predictions show that it is the coupling between attractions and repulsions that gives rise to the observed experimental trends in solution viscosity as a function of pH, concentration, and ionic strength. The parameters of the model are obtained from measurements of the second virial coefficient and net surface charge/zeta-potential, without additional fitting of the viscosity.

Visual Monitoring of Disintegration of Solid Oral Dosage Forms in Simulated Gastric Fluids Using Low-Field NMR Imaging

Compared to traditional drug release monitoring with manual sampling and testing procedures, low-field nuclear magnetic resonance (LF-NMR) imaging is a one-step, visual, non-destructive, and non-invasive measurement method. Here, we reported the application of LF-NMR to image the morphology, component, sub-diffusion, and spatial distribution of a solid oral formulation, Biyankang tablets, during dissolution in vitro. The drug ingredients with characteristic relaxation times were distinguished and localized based on the signal of standards, such as patchouli oil, Xanthium strumarium extract, and starch. The hydration, swelling, disintegration, and sub-diffusion of tablets in simulated gastric fluids (SGF) were visualized statically. All tablets showed similar expansion (37.4-42.0%) along the direction of thickness at 25 min and reached a full disintegration at 145 min, at pH 1.80-6.15, indicating pH-independent swelling and disintegration. Compared to that static immersion within 20 mL SGF, the tablet disintegration time was shortened by ~ 11% in 30 mL SGF. The application of shear reduced the time by ~ 28%, suggesting a major role of hydrodynamic condition in tablet dissolution. The ability to simultaneously visualize, distinguish, and localize drug ingredients using LF-NMR is expected to provide valuable information to develop drug release monitoring systems in vitro and potentially in vivo using small animal studies.

Use of Debye-Hückel-Henry charge measurements in early antibody development elucidates effects of non-specific association

The diffusion interaction parameter (kD) has been demonstrated to be a high-throughput technique for characterizing interactions between proteins in solution. kD reflects both attractive and repulsive interactions, including long-ranged electrostatic repulsions. Here, we plot the mutual diffusion coefficient (Dm) as a function of the experimentally determined Debye-Hückel-Henry surface charge (ZDHH) for seven human monoclonal antibodies (mAbs) in 15 mM histidine, pH 6. We find that graphs of Dm versus ZDHH intersect at ZDHH, ~ 2.6, independent of protein concentration. The same data plotted as kD vs. ZDHH shows a transition from net attractive to net repulsive interactions in the same region of the ZDHH intersection point. These data suggest that there is a minimum surface charge necessary on these mAbs needed to overcome attractive interactions.

Particle-Based Microrheology As a Tool for Characterizing Protein-Based Materials

Microrheology based on video microscopy of embedded tracer particles has the potential to be used for high-throughput protein-based materials characterization. This potential is due to a number of characteristics of the techniques, including the suitability for measurement of low sample volumes, noninvasive and noncontact measurements, and the ability to set up a large number of samples for facile, sequential measurement. In addition to characterization of the bulk rheological properties of proteins in solution, for example, viscosity, microrheology can provide insight into the dynamics and self-assembly of protein-based materials as well as heterogeneities in the microenvironment being probed. Specifically, passive microrheology in the form of multiple particle tracking and differential dynamic microscopy holds promise for applications in high-throughput characterization because of the lack of user interaction required while making measurements. Herein, recent developments in the use of multiple particle tracking and differential dynamic microscopy are reviewed for protein characterization and their potential to be applied in a high-throughput, automatable setting.

DeepSCM: An efficient convolutional neural network surrogate model for the screening of therapeutic antibody viscosity

Predicting high concentration antibody viscosity is essential for developing subcutaneous administration. Computer simulations provide promising tools to reach this aim. One such model is the spatial charge map (SCM) proposed by Agrawal and coworkers (mAbs. 2015, 8(1):43-48). SCM applies molecular dynamics simulations to calculate a score for the screening of antibody viscosity at high concentrations. However, molecular dynamics simulations are computationally costly and require structural information, a significant application bottleneck. In this work, high throughput computing was performed to calculate the SCM scores for 6596 nonredundant antibody variable regions. A convolutional neural network surrogate model, DeepSCM, requiring only sequence information, was then developed based on this dataset. The linear correlation coefficient of the DeepSCM and SCM scores achieved 0.9 on the test set (N = 1320). The DeepSCM model was applied to screen the viscosity of 38 therapeutic antibodies that SCM correctly classified and resulted in only one misclassification. The DeepSCM model will facilitate high concentration antibody viscosity screening. The code and parameters are freely available at https://github.com/Lailabcode/DeepSCM.

Opalescence Arising from Network Assembly in Antibody Solution

Opalescence of therapeutic antibody solutions is one of the concerns in drug formulation. However, the mechanistic insights into the opalescence of antibody solutions remain unclear. Here, we investigated the assembly states of antibody molecules as a function of antibody concentration. The solutions of bovine gamma globulin and human immunoglobulin G at around 100 mg/mL showed the formation of submicron-scale network assemblies. The network assembly resulted in the appearance of opalescence with a transparent blue color without the precipitates of antibodies. Furthermore, the addition of trehalose and arginine, previously known to act as protein stabilizers and protein aggregation suppressors, was able to suppress the opalescence arising from the network assembly. These results will provide an important information for evaluating and improving protein formulations.

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.

Utility of High Resolution 2D NMR Fingerprinting in Assessing Viscosity of Therapeutic Monoclonal Antibodies

Purpose

The viscosity of highly concentrated therapeutic monoclonal antibody (mAb) formulations at concentrations ≥ 100 mg/mL can significantly affect the stability, processing, and drug product development for subcutaneous delivery. An early identification of a viscosity prone mAb during candidate selection stages are often beneficial for downstream processes. Higher order structure of mAbs may often dictate their viscosity behavior at high concentration. Thus it is beneficial to gauge or rank-order their viscosity behavior using noninvasive structural fingerprinting methods and to potentially screen for suitable viscosity lowering excipients.

Methods

In this study, Dynamic Light Scattering (DLS) and 2D NMR based methyl fingerprinting were used to correlate viscosity behavior of a set of Pfizer mAbs. The viscosities of mAbs were determined. Respective Fab and Fc domains were generated for studies.

Result

Methyl fingerprinting of intact mAbs allows for differentiation of viscosity prone mAbs from well behaved ones even at 30–40 mg/ml, where bulk viscosity of the solutions are near identical. For viscosity prone mAbs, peak broadening and or distinct chemical shift changes were noted in intact and fragment fingerprints, unlike the well-behaved mAbs, indicative of protein protein interactions (PPI).

Conclusion

Fab-Fab or Fab-Fc interactions may lead to formation of protein networks at high concentration. The early transients to these network formation may be manifested through peak broadening or peak shift in the 2D NMR spectrum of mAb/mAb fragments. Such insights go beyond rank ordering mAbs based on viscosity behavior, which can be obtained by other methods as well..

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-022-03200-6.

Investigating Crystalline Protein Suspension Formulations of Pembrolizumab from MAS NMR Spectroscopy

Developing biological formulations to maintain the chemical and structural integrity of therapeutic antibodies remains a significant challenge. Monoclonal antibody (mAb) crystalline suspension formulation is a promising alternative for high concentration subcutaneous drug delivery. It demonstrates many merits compared to the solution formulation to reach a high concentration at the reduced viscosity and enhanced stability. One main challenge in drug development is the lack of high-resolution characterization of the crystallinity and stability of mAb microcrystals in the native formulations. Conventional analytical techniques often cannot evaluate structural details of mAb microcrystals in the native suspension due to the presence of visible particles, relatively small crystal size, high protein concentration, and multicomponent nature of a liquid formulation. This study demonstrates the first high-resolution characterization of mAb microcrystalline suspension using magic angle spinning (MAS) NMR spectroscopy. Crystalline suspension formulation of pembrolizumab (Keytruda, Merck & Co., Inc., Kenilworth, NJ 07033, U.S.) is utilized as a model system. Remarkably narrow ¹³C spectral linewidth of approximately 29 Hz suggests a high order of crystallinity and conformational homogeneity of pembrolizumab crystals. The impact of thermal stress and dehydration on the structure, dynamics, and stability of these mAb crystals in the formulation environment is evaluated. Moreover, isotopic labeling and heteronuclear ¹³C and ¹⁵N spectroscopies have been utilized to identify the binding of caffeine in the pembrolizumab crystal lattice, providing molecular insights into the cocrystallization of the protein and ligand. Our study provides valuable structural details for facilitating the design of crystalline suspension formulation of Keytruda and demonstrates the high potential of MAS NMR as an advanced tool for biophysical characterization of biological therapeutics.

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.

Differences in human IgG1 and IgG4 S228P monoclonal antibodies viscosity and self-interactions: Experimental assessment and computational predictions of domain interactions

Human/humanized IgG4 antibodies have reduced effector function relative to IgG1 antibodies, which is desirable for certain therapeutic purposes. However, the developability and biophysical properties for IgG4 antibodies are not well understood. This work focuses on the head-to-head comparison of key biophysical properties, such as self-interaction and viscosity, for 14 human/humanized, and chimeric IgG1 and IgG4 S228P monoclonal antibody pairs that contain the identical variable regions. Experimental measurements showed that the IgG4 S228P antibodies have similar or higher self-interaction and viscosity than that of IgG1 antibodies in 20 mM sodium acetate, pH 5.5. We report sequence and structural drivers for the increased viscosity and self-interaction detected in IgG4 S228P antibodies through a combination of experimental data and computational models. Further, we applied and extended a previously established computational model for IgG1 antibodies to predict the self-interaction and viscosity behavior for each antibody pair, providing insight into the structural characteristics and differences of these two isotypes. Interestingly, we observed that the IgG4 S228P swapped variants, where the CH3 domain was swapped for that of an IgG1, showed reduced self-interaction behavior. These domain swapped IgG4 S228P molecules also showed reduced viscosity from experiment and coarse-grained simulations. We also observed that experimental diffusion interaction parameter (kD) values have a high correlation with computational diffusivity prediction for both IgG1 and IgG4 S228P isotypes.

Abbreviations: AHc, constant region Hamaker constant; AHv, variable region Hamaker constant; CDRs, Complementarity-determining regions; CG, Coarse-grained model; CH1, Constant heavy chain 1; CH2 Constant heavy chain 2; CH3 Constant heavy chain 3; chgCH3 Effective charge on the CH3 region; CL Constant light chain; cP, Centipoise; DLS, Dynamic light scattering; Fab, Fragment antigen-binding; Fc, Fragment crystallizable; Fv, Variable domaing; (r) Radial distribution function; H1 CDR1 of Heavy Chain; H2 CDR2 of Heavy Chain; H3 CDR3 of Heavy Chain; HVI, High viscosity index; IgG1 human immunoglobulin of IgG1 subclass; IgG4 human immunoglobulin of IgG4 subclass; kD, Diffusion interaction parameter; L1 CDR1 of Light Chain; L2 CDR2 of Light Chain; L3 CDR3 of Light Chain; mAb, Monoclonal antibody; MD, Molecular dynamics; PPI Protein–protein interactions; SCM, Spatial charge map; UP-SEC, Ultra-high-performance size-exclusion chromatography; VH, Variable domain of Heavy Chain; VL, Variable domain of Light Chain

Towards an improved prediction of concentrated antibody solution viscosity using the Huggins coefficient

The viscosity of a monoclonal antibody solution must be monitored and controlled as it can adversely affect product processing, packaging and administration. Engineering low viscosity mAb formulations is challenging as prohibitive amounts of material are required for concentrated solution analysis, and it is difficult to predict viscosity from parameters obtained through low-volume, high-throughput measurements such as the interaction parameter, k_D, and the second osmotic virial coefficient, B₂₂. As a measure encompassing the effect of intermolecular interactions on dilute solution viscosity, the Huggins coefficient, k_h, is a promising candidate as a parameter measureable at low concentrations, but indicative of concentrated solution viscosity. In this study, a differential viscometry technique is developed to measure the intrinsic viscosity, [η], and the Huggins coefficient, k_h, of protein solutions. To understand the effect of colloidal protein-protein interactions on the viscosity of concentrated protein formulations, the viscometric parameters are compared to k_D and B₂₂ of two mAbs, tuning the contributions of repulsive and attractive forces to the net protein-protein interaction by adjusting solution pH and ionic strength. We find a strong correlation between the concentrated protein solution viscosity and the k_h but this was not observed for the k_D or the b₂₂, which have been previously used as indicators of high concentration viscosity. Trends observed in [η] and k_h values as a function of pH and ionic strength are rationalised in terms of protein-protein interactions.

Insights into the Conformation and Self-Association of a Concentrated Monoclonal Antibody using Isothermal Chemical Denaturation and Nuclear Magnetic Resonance

The purpose of this investigation was to highlight the utility of nuclear magnetic resonance (NMR) as a multi-attribute method for the characterization of therapeutic antibodies. In this case study, we compared results from isothermal chemical denaturation (ICD) and NMR with standard methods to relate conformational states of a model monoclonal antibody (mAb1) with protein-protein interactions (PPI) that lead to self - association in concentrated solutions. The increase in aggregation rate and relative viscosity for mAb1 was found to be both concentration and pH dependent. The free energy of unfolding (∆G⁰) from ICD and thermal analysis in dilute solutions indicated that although the native state predominated between pH 4 - pH 7, it was disrupted at the CH₂ and unfolded noncooperatively under acidic conditions. One-dimensional (1D) ¹H NMR and two-dimensional (2D) ¹³C-¹H NMR performed, in concentrated solutions, confirmed that PPI between pH 4-7 occurred while mAb1 was in the native state. NMR corroborated that mAb1 maintained a dominant native state at formulation-relevant conditions at the tested pH range, had increased global molecular tumbling dynamics at lower pH and confirmed increased PPI at higher pH conditions. This report aligns and compares typical characterization of an IgG1 with assessment of structure by NMR and provided a more precise assessment and deeper insight into the conformation of an IgG1 in concentrated solutions.

An accelerated surface-mediated stress assay of antibody instability for developability studies

High physical stability is required for the development of monoclonal antibodies (mAbs) into successful therapeutic products. Developability assays are used to predict physical stability issues such as high viscosity and poor conformational stability, but protein aggregation remains a challenging property to predict. Among different types of stresses, air–water and solid–liquid interfaces are well known to potentially trigger protein instability and induce aggregation. Yet, in contrast to the increasing number of developability assays to evaluate bulk properties, there is still a lack of experimental methods to evaluate antibody stability against interfaces. Here, we investigate the potential of a hydrophobic nanoparticle surface-mediated stress assay to assess the stability of mAbs during the early stages of development. We evaluate this surface-mediated accelerated stability assay on a rationally designed library of 14 variants of a humanized IgG4, featuring a broad span of solubility values and other developability properties. The assay could identify variants characterized by high instability against agitation in the presence of air–water interfaces. Remarkably, for the set of investigated molecules, we observe strong correlations between the extent of aggregation induced by the surface-mediated stress assay and other developability properties of the molecules, such as aggregation upon storage at 45°C, self-association (evaluated by affinity-capture self-interaction nanoparticle spectroscopy) and nonspecific interactions (estimated by cross-interaction chromatography, stand-up monolayer chromatography (SMAC), SMAC*). This highly controlled surface-mediated stress assay has the potential to complement and increase the ability of the current set of screening techniques to assess protein aggregation and developability potential of mAbs during the early stages of drug development.

Abbreviations:AC-SINS: Affinity-Capture Self-Interaction Nanoparticle Spectroscopy; AMS: Ammonium sulfate precipitation; ANS: 1-anilinonaphtalene-8-sulfonate; CIC: Cross-interaction chromatography; DLS: Dynamic light scattering; HIC: Hydrophobic interaction chromatography; HNSSA: Hydrophobic nanoparticles surface-stress assay; mAb: Monoclonal antibody; NP: Nanoparticle; SEC: Size exclusion chromatography; SMAC: Stand-up monolayer chromatography; WT: Wild type

Caffeine as a Viscosity Reducer for Highly Concentrated Monoclonal Antibody Solutions

Many monoclonal antibody (mAb) solutions exhibit high viscosity at elevated concentrations, which prevents manufacturing and injecting of concentrated mAb drug products at the small volumes needed for subcutaneous (SC) administration. Addition of excipients that interrupt intermolecular interactions is a common approach to reduce viscosity of high concentration mAb formulations. However, in some cases widely used excipients can fail to lower viscosity. Here, using infliximab and ipilimumab as model proteins, we show that caffeine effectively lowers the viscosity of both mAb formulations, whereas other common viscosity-reducing excipients, sodium chloride and arginine, do not. Furthermore, stability studies under accelerated conditions show that caffeine has no impact on stability of lyophilized infliximab or liquid ipilimumab formulations. In addition, presence of caffeine in the formulations does not affect in vitro bioactivities of infliximab or ipilimumab. Results from this study suggest that caffeine could be a useful viscosity reducing agent that complements other traditional excipients and provides viscosity reduction to a wider range of mAb drug products.

Towards in silico Process Modeling for Vaccines

Chemical, manufacturing, and control development timelines occupy a significant part of vaccine end-to-end development. In the on-going race for accelerating timelines, in silico process development constitutes a viable strategy that can be achieved through an artificial intelligence (AI)-driven or a mechanistically oriented approach. In this opinion, we focus on the mechanistic option and report on the modeling competencies required to achieve it. By inspecting the most frequent vaccine process units, we identify fluid mechanics, thermodynamics and transport phenomena, intracellular modeling, hybrid modeling and data science, and model-based design of experiments as the pillars for vaccine development. In addition, we craft a generic pathway for accommodating the modeling competencies into an in silico process development strategy.

Comparison of Huggins Coefficients and Osmotic Second Virial Coefficients of Buffered Solutions of Monoclonal Antibodies

The Huggins coefficient kH is a well-known metric for quantifying the increase in solution viscosity arising from intermolecular interactions in relatively dilute macromolecular solutions, and there has been much interest in this solution property in connection with developing improved antibody therapeutics. While numerous kH measurements have been reported for select monoclonal antibodies (mAbs) solutions, there has been limited study of kH in terms of the fundamental molecular interactions that determine this property. In this paper, we compare measurements of the osmotic second virial coefficient B22, a common metric of intermolecular and interparticle interaction strength, to measurements of kH for model antibody solutions. This comparison is motivated by the seminal work of Russel for hard sphere particles having a short-range "sticky" interparticle interaction, and we also compare our data with known results for uncharged flexible polymers having variable excluded volume interactions because proteins are polypeptide chains. Our observations indicate that neither the adhesive hard sphere model, a common colloidal model of globular proteins, nor the familiar uncharged flexible polymer model, an excellent model of intrinsically disordered proteins, describes the dependence of kH of these antibodies on B22. Clearly, an improved understanding of protein and ion solvation by water as well as dipole-dipole and charge-dipole effects is required to understand the significance of kH from the standpoint of fundamental protein-protein interactions. Despite shortcomings in our theoretical understanding of kH for antibody solutions, this quantity provides a useful practical measure of the strength of interprotein interactions at elevated protein concentrations that is of direct significance for the development of antibody formulations that minimize the solution viscosity.

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.

Monoclonal Antibody Aggregation near Silicone Oil-Water Interfaces

In this work, we study the hydrodynamic behavior of monoclonal antibodies in the presence of silicone oil-water interfaces. We model the antibody molecules using a coarse-grained 24-bead model, where two beads are used to represent each antibody domain. We consider the spatial variation of the antibody polarity in our model as each bead represents a set of hydrophilic or hydrophobic amino acids. We use the dissipative particle dynamics scheme to represent the coarse-grained force field which governs the motion of the beads. In addition, interprotein interactions are modeled using an electrostatic force field. The model parameters are determined by comparing the structure factor against experimental structure factor data ranging from a low concentration regime (10 mg/mL) to a high concentration regime (150 mg/mL). Next, we conduct simulations for a suspension of antibody molecules in the presence of silicone oil-water interfaces. Protein loss from the bulk solution is noticed as the molecules adsorb at the interface. We observe dynamic cluster formation in the solution bulk and at the interface, as the antibody molecules self-associate along their trajectories. We quantify the aggregation using a density clustering algorithm and investigate the effect of the antibody concentration on the diffusivity of the antibody solution, aggregation propensity, and protein loss from the bulk. Our study shows that numerical simulations can be an important tool for understanding the molecular mechanisms driving protein aggregation near hydrophobic interfaces.

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.

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.

Generation of mAb Variants with Less Attractive Self-Interaction but Preserved Target Binding by Well-Directed Mutation

Strongly attractive self-interaction of therapeutic protein candidates can impose challenges for manufacturing, filling, stability, and administration due to elevated viscosity or aggregation propensity. Suitable formulations can mitigate these issues to a certain extent. Understanding the self-interaction mechanism on a molecular basis and rational protein engineering provides a more fundamental approach, and it can save costs and efforts as well as alleviate risks at later stages of development. In this study, we used computational methods for the identification of aggregation-prone regions in a mAb and generated mutants based on these findings. We applied hydrogen-deuterium exchange mass spectrometry to identify distinct self-interaction hot spots. Ultimately, we generated mAb variants based on a combination of both approaches and identified mutants with low attractive self-interaction propensity, minimal off-target binding, and even improved target binding. Our data show that the introduction of arginine in spatial proximity to hydrophobic patches is highly beneficial on all these levels. For our mAb, variants that contain more than one aspartate residue flanking to the hydrophobic HCDR3 show decreased attractive self-interaction at unaffected off-target and target binding. The combined engineering strategy described here underlines the high potential of understanding self-interaction in the early stages of development to predict and reduce the risk of failure in subsequent development.

Rheological investigation of collagen, fibrinogen, and thrombin solutions for drop-on-demand 3D bioprinting

Collagen, fibrinogen, and thrombin proteins in aqueous buffer solutions are widely used as precursors of natural biopolymers in three-dimensional (3D) bioprinting applications. The proteins are sourced from animals and their quality may vary from batch to batch, inducing differences in the rheological properties of such solutions. In this work, we investigate the rheological response of collagen, fibrinogen, and thrombin protein solutions in bulk and at the solution/air interface. Interfacial rheological measurements show that fibrous collagen, fibrinogen and globular thrombin proteins adsorb and aggregate at the solution/air interface, forming a viscoelastic solid film at the interface. The viscoelastic film corrupts the bulk rheological measurements in rotational rheometers by contributing to an apparent yield stress, which increases the apparent bulk viscosity up to shear rates as high as 1000 s^-1. The addition of a non-ionic surfactant, such as polysorbate 80 (PS80) in small amounts between 0.001 and 0.1 v/v%, prevents the formation of the interfacial layer, allowing the estimation of true bulk viscosity of the solutions. The estimation of viscosity not only helps in identifying those protein solutions that are potentially printable with drop-on-demand (DOD) inkjet printing but also detects inconsistencies in flow behavior among the batches.