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

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.

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.'

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.

Exploring HCl-HCl interactions: QZVPP calculations, improved Lennard-Jones potential, and second virial coefficient analysis for thermodynamics and industrial applications

In this paper, we present a comprehensive analysis of HCl-HCl interactions, including QZVPP calculations, energy fitting, conformation validation, and the determination of the second virial coefficient B using improved Lennard-Jones (ILJ) potential parameters. To acquire accurate interaction energies, initial QZVPP calculations are performed on approximately 1851 randomly generated HCl-HCl conformations. Then, these energies are used to fit an improved Lennard-Jones potential energy surface, allowing for a robust description of HCl-HCl interactions. The ILJ potential parameters are then used to validate particular HCl dimer conformations, ensuring their stability and consistency with experimental observations. The correlation between calculated and experimental conformations strengthens the validity of the ILJ potential parameters. In addition, the second viral coefficient B is calculated at various temperatures using the ILJ potential. The obtained B values are compared to experimental data, demonstrating close agreement, and validating the ILJ potential's ability to accurately capture the intermolecular interactions and gas-phase behavior of the HCl-HCl system. The results of this study demonstrate the effective implementation of QZVPP calculations, energy fitting, and ILJ potential parameters in validating HCl-HCl conformations and accurately determining the second virial coefficient B. The high degree of concordance between calculated B values and experimental data demonstrates the validity of the ILJ potential and its suitability for modeling HCl-HCl interactions. This research contributes to a greater comprehension of HCl-HCl interactions and their implications for numerous chemical and atmospheric processes. The validated conformations, energy fitting method, and calculated second virial coefficients provide valuable instruments for future research and pave the way for more accurate modeling and simulations of HCl-HCl systems.

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.

New Facet in Viscometry of Charged Associating Polymer Systems in Dilute Solutions

The peculiarities of viscosity data treatment for two series of polymer systems exhibiting associative properties: brush-like amphiphilic copolymers-charged alkylated N-methyl-N-vinylacetamide and N-methyl-N-vinylamine copolymer (MVAA-co-MVAC_nH_2n+1) and charged chains of sodium polystyrene-4-sulfonate (PSSNa) in large-scale molecular masses (MM) and in extreme-scale of the ionic strength of solutions were considered in this study. The interest in amphiphilic macromolecular systems is explained by the fact that they are considered as micellar-forming structures in aqueous solutions, and these structures are able to carry hydrophobic biologically active compounds. In the case of appearing the hydrophobic interactions, attention was paid to discussing convenient ways to extract the correct value of intrinsic viscosity η from the combined analysis of Kraemer and Huggins plots, which were considered as twin plots. Systems and situations were demonstrated where intrachain hydrophobic interactions occurred. The obtained data were discussed in terms of lnηr vs. cη plots as well as in terms of normalized scaling relationships where ηr was the relative viscosity of the polymer solution. The first plot allowed for the detection and calibration of hydrophobic interactions in polymer chains, while the second plot allowed for the monitoring of the change in the size of charged chains depending on the ionic strength of solutions.

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.

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.

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

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

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.