Small-angle neutron scattering characterization of monoclonal antibody conformations and interactions at high concentrations

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

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.

Combining Scattering Experiments and Colloid Theory to Characterize Charge Effects in Concentrated Antibody Solutions

Charges and their contribution to protein-protein interactions are essential for the key structural and dynamic properties of monoclonal antibody (mAb) solutions. In fact, they influence the apparent molecular weight, the static structure factor, the collective diffusion coefficient, or the relative viscosity, and their concentration dependence. Further, charges play an important role in the colloidal stability of mAbs. There exist standard experimental tools to characterize mAb net charges, such as the measurement of the electrophoretic mobility, the second virial coefficient, or the diffusion interaction parameter. However, the resulting values are difficult to directly relate to the actual overall net charge of the antibody and to theoretical predictions based on its known molecular structure. Here, we report the results of a systematic investigation of the solution properties of a charged IgG1 mAb as a function of concentration and ionic strength using a combination of electrophoretic measurements, static and dynamic light scattering, small-angle X-ray scattering, and tracer particle-based microrheology. We analyze and interpret the experimental results using established colloid theory and coarse-grained computer simulations. We discuss the potential and limits of colloidal models for the description of the interaction effects of charged mAbs, in particular pointing out the importance of incorporating shape and charge anisotropy when attempting to predict structural and dynamic solution properties at high concentrations.

Impact of Glycosylation on Protein-Protein Self-Interactions of Monoclonal Antibodies

Protein self-interactions measured via second osmotic virial coefficients (B22) and dynamic light scattering interaction parameter values (kD) are often used as metrics for assessing the favorability of protein candidates and different formulations during monoclonal antibody (MAb) product development. Model predictions of B22 or kD typically do not account for glycans, though glycosylation can potentially impact experimental MAb self-interactions. To the best of our knowledge, the impact of MAb glycosylation on the experimentally measured B22 and kD values has not yet been reported. B22 and kD values of two fully deglycosylated MAbs and their native (i.e., fully glycosylated) counterparts were measured by light scattering over a range of pH and ionic strength conditions. Significant differences between B22 and kD of the native and deglycosylated forms were observed at a range of low to high ionic strengths used to modulate the effect of electrostatic contributions. Differences were most pronounced at low ionic strength, indicating that electrostatic interactions are a contributing factor. Though B22 and kD values were statistically equivalent at high ionic strengths where electrostatics were fully screened, we observed protein-dependent qualitative differences, which indicate that steric interactions may also play a role in the observed B22 and kD differences. A domain-level coarse-grained molecular model accounting for charge differences was considered to potentially provide additional insight but was not fully predictive of the behavior across all of the solution conditions investigated. This highlights that both the level of modeling and lack of inclusion of glycans may limit existing models in making quantitatively accurate predictions of self-interactions.

Electrostatically Mediated Attractive Self-Interactions and Reversible Self-Association of Fc-Fusion Proteins

Attractive self-interactions and reversible self-association are implicated in many problematic solution behaviors for therapeutic proteins, such as irreversible aggregation, elevated viscosity, phase separation, and opalescence. Protein self-interactions and reversible oligomerization of two Fc-fusion proteins (monovalent and bivalent) and the corresponding fusion partner protein were characterized experimentally with static and dynamic light scattering as a function of pH (5 and 6.5) and ionic strength (10 mM to at least 300 mM). The fusion partner protein and monovalent Fc-fusion each displayed net attractive electrostatic self-interactions at pH 6.5 and net repulsive electrostatic self-interactions at pH 5. Solutions of the bivalent Fc-fusion contained higher molecular weight species that prevented quantification of typical interaction parameters (B22 and kD). All three of the proteins displayed reversible self-association at pH 6.5, where oligomers dissociated with increased ionic strength. Coarse-grained molecular simulations were used to model the self-interactions measured experimentally, assess net self-interactions for the bivalent Fc-fusion, and probe the specific electrostatic interactions between charged amino acids that were involved in attractive electrostatic self-interactions. Mayer-weighted pairwise electrostatic energies from the simulations suggested that attractive electrostatic self-interactions at pH 6.5 for the two Fc-fusion proteins were due to cross-domain interactions between the fusion partner domain(s) and the Fc domain.

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

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

Modulating shape transition in surfactant stabilized reverse microemulsions

The formation of reverse microemulsions (RMs) of spherical shape in the oil/water/surfactant ternary mixture at high molar ratio of water to surfactant (ω) is well established. Using dynamic light scattering, small-angle X-ray and neutron scattering, we elucidate the formation of non-spherical reverse microemulsions stabilised by sodium bis(2-ethylhexyl) sulfosuccinate (AOT) at ω = 10 and volume fractions of the dispersed phase, Φ, ranging from 0.005 to 0.20. In addition, we propose a strategy to tune the aspect ratio of non-spherical droplets and colloidal interactions by (i) varying the volume fraction of the dispersed phase (ii) changing the temperature, and (iii) by substituting the aliphatic oil with a mixture of aliphatic and aromatic hydrocarbons. This tunability of anisotropy along with a precise control of the interactions in the RMs, their ability to form spontaneously and their thermodynamic stability is crucial to provide a handle on reaction kinetics, synthesis of anisotropic nanoparticles as well as for their application as lubricants and viscosity modifiers.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

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

Raman Spectroscopic Analysis of Highly-Concentrated Antibodies under the Acid-Treated Conditions

PURPOSE

Antibody drugs are usually formulated as highly-concentrated solutions, which would easily generate aggregates, resulting in loss of efficacy. Although low pH increases the colloidal dispersion of antibodies, acid denaturation can be an issue. Therefore, knowing the physical properties at low pH under high concentration conditions is important.

METHODS

Raman spectroscopy was used to investigate pH-induced conformational changes of antibodies at 50 mg/ml. Experiments in pH 3 to 7 were performed for human serum IgG and recombinant rituximab.

RESULTS

We detected the evident changes at pH 3 in Tyr and Trp bands, which are the sensitive markers of intermolecular interactions. Thermal transition analysis over the pH range demonstrated that the thermal transition temperature (Tm) was highest at pH 3. Acid-treated and neutralized one showed higher Tm than that of pH 7, indicating that their extent of intermolecular interactions correlated with the Tm values. Onset temperature was clearly different between concentrated and diluted samples. Colloidal analyses confirmed the findings of the Raman analysis.

CONCLUSION

Our studies demonstrated the positive correlation between Raman analysis and colloidal information, validating as a method for evaluating antibody conformation associated with aggregation propensities.

Enabling Efficient Design of Biological Formulations Through Advanced Characterization

The present review summarizes the use of differential scanning calorimetry (DSC) and scattering techniques in the context of protein formulation design and characterization. The scattering techniques include wide angle X-ray diffractometry (XRD), small-angle neutron scattering (SANS) and small-angle X-ray scattering (SAXS). While DSC is valuable for understanding thermal behavior of the excipients, XRD provides critical information about physical state of solutes during freezing, annealing and in the final lyophile. However, as these techniques lack the sensitivity to detect biomolecule-related transitions, complementary characterization techniques such as small-angle scattering can provide valuable insights.

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

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

Modelling network formation in folded protein hydrogels by cluster aggregation kinetics

Globular folded protein-based hydrogels are becoming increasingly attractive due to their specific biological functionality, as well as their responsiveness to stimuli. By modelling folded proteins as colloids, there are rich opportunities to explore network formation mechanisms in protein hydrogels that negate the need for computationally expensive simulations which capture the full complexity of proteins. Here we present a kinetic lattice-based model which simulates the formation of irreversibly chemically crosslinked, folded protein-based hydrogels. We identify the critical point of gel percolation, explore the range of network regimes covering diffusion-limited to reaction-limited cluster aggregation (DLCA and RLCA, respectively) network formation mechanisms and predict the final network structure, fractal dimensions and final gel porosity. We reveal a crossover between DLCA and RLCA mechanisms as a function of protein volume fraction and show how the final network structure is governed by the structure at the percolation point, regardless of the broad variation of non-percolating cluster masses observed across all systems. An analysis of the pore size distribution in the final network structures reveals that, approaching RLCA, gels have larger maximal pores than the DLCA counterparts for both volume fractions studied. This general kinetic model and the analysis tools generate predictions of network structure and concurrent porosity over a broad range of experimentally controllable parameters that are consistent with current expectations and understanding of experimental results.

Analytical Techniques for the Characterization and Quantification of Monoclonal Antibodies

Monoclonal antibodies (mAbs) are a fast-growing class of biopharmaceuticals. They are widely used in the identification and detection of cell makers, serum analytes, and pathogenic agents, and are remarkably used for the cure of autoimmune diseases, infectious diseases, or malignancies. The successful application of therapeutic mAbs is based on their ability to precisely interact with their appropriate target sites. The precision of mAbs rely on the isolation techniques delivering pure, consistent, stable, and safe lots that can be used for analytical, diagnostic, or therapeutic applications. During the creation of a biologic, the key quality features of a particular mAb, such as structure, post-translational modifications, and activities at the biomolecular and cellular levels, must be characterized and profiled in great detail. This implies the requirement of powerful state of the art analytical techniques for quality control and characterization of mAbs. Until now, various analytical techniques have been developed to characterize and quantify the mAbs according to the regulatory guidelines. The present review summarizes the major techniques used for the analyses of mAbs which include chromatographic, electrophoretic, spectroscopic, and electrochemical methods in addition to the modifications in these methods for improving the quality of mAbs. This compilation of major analytical techniques will help students and researchers to have an overview of the methodologies employed by the biopharmaceutical industry for structural characterization of mAbs for eventual release of therapeutics in the drug market.

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

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

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.

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.

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.

How neutron scattering techniques benefit investigating structures and dynamics of monoclonal antibody

Over the past several decades, great progresses have been made for the pharmaceutical industry of monoclonal antibody (mAb). More and more mAb products were approved for human therapeutics. This review describes the state of art of utilizing neutron scattering to investigate mAbs, in the aspects of structures, dynamics, physicochemical stability, functionality, etc. Firstly, brief histories of mAbs and neutron scattering, as well as some basic knowledges and principles of neutron scattering were introduced. Then specific examples were demonstrated. For the structure and structural evolution investigation of in dilute and concentrated mAbs solution, in situ small angle neutron scattering (SANS) was frequently utilized. Neutron reflectometry (NR) is powerful to probe the absorption behaviors of mAbs on various surfaces and interfaces. While for dynamic investigation, quasi-elastic scattering techniques such as neutron spin echo (NSE) demonstrate the capabilities. With this review, how to utilize and take advantages of neutron scattering on investigating structures and dynamics of mAbs were demonstrated and discussed.

Spray-Dried Monoclonal Antibody Suspension for High-Concentration and Low-Viscosity Subcutaneous Injection

Administration of highly concentrated monoclonal antibodies (mAbs) through injection is often not possible as the viscosity can be readily above 50 mPa·s when the concentration exceeds 150 mg/mL. Besides, highly concentrated mAb solutions always exhibit increased aggregation propensity and lower stability, which raise the difficulty for the successful development of highly concentrated mAb formulations. We hereby explored the possibility of suspension as another formulation form for high-concentration proteins to reduce viscosity and maintain stability. Specifically, we demonstrated that spray drying can serve as a process to prepare particles for suspension. Particles prepared from formulations with different mAb/trehalose mass ratios displayed good physical stability and antibody binding affinity, as indicated by circular dichroism, fluorescence spectroscopy, and surface plasmon resonance (SPR)-based bioassay analyses. During spray drying, a surface tension-dominated enrichment of mAb on the particle surface was observed, but this did not show a significant negative impact on mAb stability. Spray-dried particles were subsequently suspended into benzyl benzoate, and the resulting suspension showed good stability and a lower viscosity when compared to its counterpart solution. Furthermore, mAbs recovered from the suspension maintained their conformational structure. Our study demonstrated that the suspension displayed low viscosity and good physical stability, so it may offer novel opportunities for the preparation of highly concentrated protein formulations.

Effects of Monovalent Salt on Protein-Protein Interactions of Dilute and Concentrated Monoclonal Antibody Formulations

In this study, we used sodium chloride (NaCl) to extensively modulate non-specific protein-protein interactions (PPI) of a humanized anti-streptavidin monoclonal antibody class 2 molecule (ASA-IgG2). The changes in PPI with varying NaCl (C_NaCl) and monoclonal antibody (mAb) concentration (C_mAb) were assessed using the diffusion interaction parameter k_D and second virial coefficient B₂₂ measured from solutions with low to moderate C_mAb. The effective structure factor S(q)_eff measured from concentrated mAb solutions using small-angle X-ray and neutron scattering (SAXS/SANS) was also used to characterize the PPI. Our results found that the nature of net PPI changed not only with C_NaCl, but also with increasing C_mAb. As a result, parameters measured from dilute and concentrated mAb samples could lead to different predictions on the stability of mAb formulations. We also compared experimentally determined viscosity results with those predicted from interaction parameters, including k_D and S(q)_eff. The lack of a clear correlation between interaction parameters and measured viscosity values indicates that the relationship between viscosity and PPI is concentration-dependent. Collectively, the behavior of flexible mAb molecules in concentrated solutions may not be correctly predicted using models where proteins are considered to be uniform colloid particles defined by parameters derived from low C_mAb.

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.

Peptide Design and Self-assembly into Targeted Nanostructure and Functional Materials

Peptides have been extensively utilized to construct nanomaterials that display targeted structure through hierarchical assembly. The self-assembly of both rationally designed peptides derived from naturally occurring domains in proteins as well as intuitively or computationally designed peptides that form β-sheets and helical secondary structures have been widely successful in constructing nanoscale morphologies with well-defined 1-d, 2-d, and 3-d architectures. In this review, we discuss these successes of peptide self-assembly, especially in the context of designing hierarchical materials. In particular, we emphasize the differences in the level of peptide design as an indicator of complexity within the targeted self-assembled materials and highlight future avenues for scientific and technological advances in this field.

Molecular Flexibility of Antibodies Preserved Even in the Dense Phase after Macroscopic Phase Separation

Antibody therapies are typically based on high-concentration formulations that need to be administered subcutaneously. These conditions induce several challenges, inter alia a viscosity suitable for injection, sufficient solution stability, and preservation of molecular function. To obtain systematic insights into the molecular factors, we study the dynamics on the molecular level under strongly varying solution conditions. In particular, we use solutions of antibodies with poly(ethylene glycol), in which simple cooling from room temperature to freezing temperatures induces a transition from a well-dispersed solution into a phase-separated and macroscopically arrested system. Using quasi-elastic neutron scattering during in situ cooling ramps and in prethermalized measurements, we observe a strong decrease in antibody diffusion, while internal flexibility persists to a significant degree, thus ensuring the movement necessary for the preservation of molecular function. These results are relevant for a more dynamic understanding of antibodies in high-concentration formulations, which affects the formation of transient clusters governing the solution viscosity.

Reversible Self-Association in Lactate Dehydrogenase during Freeze-Thaw in Buffered Solutions Using Neutron Scattering

The aims of this work were to evaluate the effect of freezing and thawing stresses on lactate dehydrogenase (LDH) stability under three conditions. (i) In a solution buffered with sodium phosphate (NaP; 10 and 100 mM). The selective crystallization of disodium hydrogen phosphate during freezing caused a pronounced pH shift. (ii) In a solution buffered with histidine, where there was no pH shift due to buffer salt crystallization. (iii) At different concentrations of LDH so as to determine the self-stabilizing ability of LDH. The change in LDH tetrameric conformation was measured by small-angle neutron scattering (SANS). The pH of the phosphate buffer solutions was monitored as a function of temperature to quantify the pH shift. The conditions of buffer component crystallization from solution were identified using low-temperature X-ray diffractometry. Dynamic light scattering (DLS) enabled us to determine the effect of freeze-thawing on the protein aggregation behavior. LDH, at a high concentration (1000 μg/mL; buffer concentration 10 mM), has a pronounced self-stabilizing effect and did not aggregate after five freeze-thaw cycles. At lower LDH concentrations (10 and 100 μg/mL), only with the selection of an appropriate buffer, irreversible aggregation could be avoided. While SANS provided qualitative information with respect to protein conformation, the insights from DLS were quantitative with respect to the particle size of the aggregates. SANS is the only technique which can characterize the protein both in the frozen and thawed states.

New Model to Predict the Concentration-Dependent Viscosity of Monoclonal Antibody Solutions

Concentrated monoclonal antibody solutions exhibit high solution viscosity, which is experimentally measured to be ∼1-2 orders of magnitude higher than the viscosity of water. However, physical processes responsible for the high antibody viscosity are not fully understood. We show that fluid occlusion due to the trapped solvent molecules within the boundaries formed by the aggregated antibodies is responsible for the elevated solution viscosity. We develop a theory to predict the viscosity of monoclonal antibodies based on the geometry of the antibody molecule and the aggregate morphology. We validate our theory with experiments and highlight useful insights obtained from the viscosity equation which could help in controlling the drug viscosity at the molecular design stage itself.

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.

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.

Toward a Suite of Coarse-Grained Models for Molecular Simulation of Monoclonal Antibodies and Therapeutic Proteins

A series of coarse-grained models for molecular simulation of proteins are considered, with emphasis on the application of predicting protein-protein self-interactions for monoclonal antibodies (MAbs). As an illustrative example and for quantitative comparison, the models are used to predict osmotic virial coefficients over a broad range of attractive and repulsive self-interactions and solution conditions for a series of MAbs where the second osmotic virial coefficient has been experimentally determined in prior work. The models are compared based on how well they can predict experimental behavior, their computational burdens, and scalability. An intermediate-resolution model is also introduced that can capture specific electrostatic interactions with improved efficiency and similar or improved accuracy when compared to the previously published models. Guidance is included for the selection of coarse-grained models more generally for capturing a balance of electrostatic, steric, and short-ranged nonelectrostatic interactions for proteins from low to high concentrations.

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.

Temperature Dependence of Protein Solution Viscosity and Protein-Protein Interactions: Insights into the Origins of High-Viscosity Protein Solutions

Protein solution viscosity (η) as a function of temperature was measured at a series of protein concentrations under a range of formulation conditions for two monoclonal antibodies (MAbs) and a globular protein (aCgn). Based on theoretical arguments, a strong temperature dependence for protein-protein interactions (PPI) indicates highly anisotropic, short-ranged attractions that could lead to higher solution viscosities. The semi-empirical Ross-Minton model was used to determine the apparent intrinsic viscosity, shape, and "crowding" factors for each protein as a function of temperature and formulation conditions. The apparent intrinsic viscosity was independent of temperature for aCgn, while a slight decrease with increasing temperature was observed for the MAbs. The temperature dependence of solution viscosity was analyzed using the Andrade-Eyring equation to determine the effective activation energy of viscous flow (E_a,η). While E_a,η values were different for each protein, they were independent of formulation conditions for a given protein. PPI were quantified via the osmotic second virial coefficient (B₂₂) and the protein diffusion interaction parameter (k_D) as a function of temperature under the same formulation conditions as the viscosity measurements. Net interactions ranged from strongly attractive to repulsive by changing formulation pH and ionic strength for each protein. Overall, larger activation energies for PPI corresponded to larger activation energies for η, and those were predictive of the highest η values at higher protein concentrations.

Coarse-Grained Molecular Dynamics Simulations for Understanding the Impact of Short-Range Anisotropic Attractions on Structure and Viscosity of Concentrated Monoclonal Antibody Solutions

Understanding protein-protein interactions in concentrated therapeutic monoclonal antibody (mAb) solutions is desirable for improved drug discovery, processing, and administration. Here, we deduce both the net protein charge and the magnitude and geometry of short-ranged, anisotropic attractions of a mAb across multiple concentrations and cosolute conditions by comparing structure factors S(q) obtained from small-angle X-ray scattering experiments with those from molecular dynamics (MD) simulations. The simulations, which utilize coarse-grained 12-bead models exhibiting a uniform van der Waals attraction, uniform electrostatic repulsion, and short-range attractions between specific beads, are versatile enough to fit S(q) of a wide range of protein concentrations and ionic strength with the same charge on each bead and a single anisotropic short-range attraction strength. Cluster size distributions (CSDs) obtained from best fit simulations reveal that the experimental structure is consistent with small reversible oligomers in even low viscosity systems and help quantify the impact of these clusters on viscosity. The ability to systematically use experimental S(q) data together with MD simulations to discriminate between different possible protein-protein interactions, as well as to predict viscosities from protein CSDs, is beneficial for designing mAbs and developing formulation strategies that avoid high viscosities and aggregation at high concentration.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.

Self-diffusion of a highly concentrated monoclonal antibody by fluorescence correlation spectroscopy: insight into protein-protein interactions and self-association

The dynamic behavior of monoclonal antibodies (mAbs) at high concentration provides insight into protein microstructure and protein-protein interactions (PPI) that influence solution viscosity and protein stability. At high concentration, interpretation of the collective-diffusion coefficient Dc, as determined by dynamic light scattering (DLS), is highly challenging given the complex hydrodynamics and PPI at close spacings. In contrast, self-diffusion of a tracer particle by Brownian motion is simpler to understand. Herein, we develop fluorescence correlation spectroscopy (FCS) for the measurement of the long-time self-diffusion of mAb2 over a wide range of concentrations and viscosities in multiple co-solute formulations with varying PPI. The normalized self-diffusion coefficient D0/Ds (equal to the microscopic relative viscosity ηeff/η0) was found to be smaller than η/η0. Smaller ratios of the microscopic to macroscopic viscosity (ηeff/η) are attributed to a combination of weaker PPI and less self-association. The interaction parameters extracted from fits of D0/Ds with a length scale dependent viscosity model agree with previous measurements of PPI by SLS and SAXS. Trends in the degree of self-association, estimated from ηeff/η with a microviscosity model, are consistent with oligomer sizes measured by SLS. Finally, measurements of collective diffusion and osmotic compressibility were combined with FCS data to demonstrate that the changes in self-diffusion between formulations are due primarily to changes in the protein-protein friction in these systems, and not to protein-solvent friction. Thus, FCS is a robust and accessible technique for measuring mAb self-diffusion, and, by extension, microviscosity, PPI and self-association that govern mAb solution dynamics.

Surface heterogeneity affects percolation and gelation of colloids: dynamic simulations with random patchy spheres

Surface heterogeneity of colloidal particles has a significant impact on their structure in solution and their rheological properties. During particle synthesis, heterogeneous chemical functionalization, processes of self-assembly, or phase separation, can all lead to heterogeneous colloidal surfaces which impart anisotropic interactions to suspended particles. Additionally, an important class of colloids, biological macromolecules, exhibit similar localized, short-ranged, anisotropic interactions, which have a significant impact on their solution properties. Therefore, understanding the assembly and rheology of such colloids can provide insight into a wide variety of relevant physical systems. In this computational study, we investigate dispersions of particles having surface patches with randomized functionality as a model for heterogeneous colloids. We use Brownian dynamics simulations with hydrodynamic interactions to explore the differences between these random patchy particles and homogeneous (or isotropic) particles. The common basis used for comparing dispersions of particles with different surface functionality is equality of the second virial coefficient, so that dispersions of particles with different patterns of surface heterogeneity are similar thermodynamically at low particle concentrations. We show that at modest particle concentrations, significant deviations from the isotropic model are evident in the dispersion micro-structure, giving drastically different percolation transition points depending on the degree of surface heterogeneity. However, these deviations can be rationalized and a universal percolation criteria derived in terms of the osmotic pressure of the dispersion. Heterogeneous interactions also impose extra constraints on the relative translation and rotation between neighboring particles, which increase the viscosity and elastic modulus of aggregated dispersions and gels built from heterogeneous colloids and shifts the gel point measurably.

A Colloid Approach to Self-Assembling Antibodies

Concentrated solutions of monoclonal antibodies have attracted considerable attention due to their importance in pharmaceutical formulations; yet, their tendency to aggregate and the resulting high viscosity pose considerable problems. Here we tackle this problem by a soft condensed matter physics approach, which combines a variety of experimental measurements with a patchy colloid model, amenable of analytical solution. We thus report results of structural antibodies and dynamic properties obtained through scattering methods and microrheological experiments. We model the data using a colloid-inspired approach, explicitly taking into account both the anisotropic shape of the molecule and its charge distribution. Our simple patchy model is able to disentangle self-assembly and intermolecular interactions and to quantitatively describe the concentration-dependence of the osmotic compressibility, collective diffusion coefficient, and zero shear viscosity. Our results offer new insights on the key problem of antibody formulations, providing a theoretical and experimental framework for a quantitative assessment of the effects of additional excipients or chemical modifications and a prediction of the resulting viscosity.

Web-based display of protein surface and pH-dependent properties for assessing the developability of biotherapeutics

Protein instability leads to reversible self-association and irreversible aggregation which is a major concern for developing new biopharmaceutical leads. Protein solution behaviour is dictated by the physicochemical properties of the protein and the solution. Optimising protein solutions through experimental screens and targeted protein engineering can be a difficult and time consuming process. Here, we describe development of the protein-sol web server, which was previously restricted to protein solubility prediction from amino acid sequence. Tools are presented for calculating and mapping patches of hydrophobicity and charge on the protein surface. In addition, predictions of folded state stability and net charge are displayed as a heatmap for a range of pH and ionic strength conditions. Tools are evaluated in the context of antibodies, their fragments and interactions. Surprisingly, antibody-antigen interfaces are, on average, at least as polar as Fab surfaces. This benchmarking process provides the user with thresholds with which to assess non-polar surface patches, and possible solubility implications, in proteins of interest. Stability heatmaps compare favourably with experimental data for CH2 and CH3 domains. Display and quantification of surface polarity and pH/ionic strength dependence will be useful generally for investigation of protein biophysics.

Colloidal systems with a short-range attraction and long-range repulsion: Phase diagrams, structures, and dynamics

Colloidal systems with both a short-range attraction and long-range repulsion (SALR) have rich phases compared with the traditional hard sphere systems or sticky hard sphere systems. The competition between the short-range attraction and long-range repulsion results in the frustrated phase separation, which leads to the formation of intermediate range order (IRO) structures and introduces new phases to both equilibrium and nonequilibrium phase diagrams, such as clustered fluid, cluster percolated fluid, Wigner glass, and cluster glass. One hallmark feature of many SALR systems is the appearance of the IRO peak in the interparticle structure factor, which is associated with different types of IRO structures. The relationship between the IRO peak and the clustered fluid state has been careful investigated. Not surprisingly, the morphology of clusters in solutions can be affected and controlled by the SALR potential. And the effect of the SALR potential on the dynamic properties is also reviewed here. Even though much progress has been made in understanding SALR systems, many future works are still needed to have quantitative comparisons between experiments and simulations/theories and understand the differences from different experimental systems. Owing to the large parameter space available for SALR systems, many exciting features of SALR systems are not fully explored yet. Because proteins in low-salinity solutions have SALR interactions, the understanding of SALR systems can greatly help understand protein behavior in concentrated solutions or crowded conditions.

Protein-Protein Interactions of Highly Concentrated Monoclonal Antibody Solutions via Static Light Scattering and Influence on the Viscosity

The ability to design and formulate mAbs to minimize attractive interactions at high concentrations is important for protein processing, stability, and administration, particularly in subcutaneous delivery, where high viscosities are often challenging. The strength of protein-protein interactions (PPIs) of an IgG1 and IgG4 monoclonal antibody (mAb) from low to high concentration was determined by static light scattering (SLS) and used to understand viscosity data. The PPI were tuned using NaCl and five organic ionic co-solutes. The PPI strength was quantified by the normalized structure factor S(0)/ S(0)_HS and Kirkwood-Buff integral G₂₂/ G_22,HS (HS = hard sphere) determined from the SLS data and also by fits with (1) a spherical Yukawa potential and (2) an interacting hard sphere (IHS) model, which describes attraction in terms of hypothetical oligomers. The IHS model was better able to capture the scattering behavior of the more strongly interacting systems (mAb and/or co-solute) than the spherical Yukawa potential. For each descriptor of PPI, linear correlations were obtained between the viscosity at high concentration (200 mg/mL) and the interaction strengths evaluated both at low (20 mg/mL) and high concentrations (200 mg/mL) for a given mAb. However, the only parameter that provided a correlation across both mAbs was the oligomer mass ratio ( m_oligomer/ m_{monomer+dimer}) from the IHS model, indicating the importance of self-association (in addition to the direct influence of the attractive PPI) on the viscosity.

Structure and Relaxation in Solutions of Monoclonal Antibodies

Reversible self-association of therapeutic antibodies is a key factor in high protein solution viscosities. In the present work, a coarse-grained computational model accounting for electrostatic, dispersion, and long-ranged hydrodynamic interactions of two model monoclonal antibodies is applied to understand the nature of self-association, predicting the solution microstructure and resulting transport properties of the solution. For the proteins investigated, the structure factor across a range of solution conditions shows quantitative agreement with neutron-scattering experiments. We observe a homogeneous, dynamical association of the antibodies with no evidence of phase separation. Calculations of self-diffusivity and viscosity from coarse-grained dynamic simulations show the appropriate trends with concentration but, respectively, over- and under-predict the experimentally measured values. By adding constraints to the self-associated clusters that rigidify them under flow, prediction of the transport properties is significantly improved with respect to experimental measurements. We hypothesize that these rigidity constraints are associated with missing degrees of freedom in the coarse-grained model resulting from patchy and heterogeneous interactions among coarse-grained domains. These results demonstrate how structural anisotropy and anisotropy of interactions generated by features at the 2-5 nm length scale in antibodies are sufficient to recover the dynamics and rheological properties of these important macromolecular solutions.