Relating Protein-Protein Interactions and Aggregation Rates From Low to 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.'

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.

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.

Injectable hydrogel particles for amorphous solid formulation of biologics

The fast pace of breakthroughs in cancer immunotherapy, combined with the new paradigm of moving toward high-concentration dosages and combinatorial treatments, is generating new challenges in the formulation of biologics. To address these challenges, we describe a method of formulation that enables high-concentration injectable and stable formulation of biologics as amorphous solids in aqueous suspension. This technology combines the benefits of liquid formulation with the stability of solid formulation and eliminates the need for drying and reconstitution. This widely applicable formulation integrates the amorphous solid forms of antibodies with the injectability, lubricity, and tunability of soft alginate hydrogel particles using a minimal process. The platform was evaluated for anti-PD-1 antibody pembrolizumab and human immunoglobulin G at concentrations up to 300 mg/mL with confirmed quality after release. The soft nature of the hydrogel matrix allowed packing the particles to high volume fractions.

Conformational Changes and Drivers of Monoclonal Antibody Liquid-Liquid Phase Separation

Liquid-liquid phase separation is a phenomenon within biology whereby proteins can separate into dense and more dilute phases with distinct properties. Three antibodies that undergo liquid-liquid phase separation were characterized in the protein-rich and protein-poor phases. In comparison to the protein-poor phase, the protein-rich phase demonstrates more blue-shift tryptophan emissions and red-shifted amide I absorbances. Large changes involving conformational isomerization around disulfide bonds were observed using Raman spectroscopy. Amide I and protein fluorescence differences between the phases persisted to temperatures above the critical temperature but ceased at the temperature at which aggregation occurred. In addition, large changes occurred in the structural organization of water molecules within the protein-rich phase for all three antibodies. It is hypothesized that as the proteins have the same chemical potential in both phases, the protein viscosity is higher in the protein-rich phase resulting in slowed diffusion dependent protein aggregation in this phase. For all three antibodies we performed accelerated stability studies and found that the protein-rich phase aggregated at the same rate or slower than the protein-poor phase.

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.

Roughness of Production Conditions: Does It Really Affect Stability of IgG-Based Antivenoms?

Antivenoms contain either pure animal IgGs or their fragments as an active substance, and are the only specific therapeutics against envenomation arising from snakebites. Although they are highly needed, the low sustainability of such preparations’ manufacture causes constant global shortages. One reason for this is the stability of the product, which contributes not only to the manufacture sustainability, but the product safety as well. It has been hypothesized that the roughness of conditions to which IgGs are exposed during downstream purification disturbs their conformation, making them prone to aggregation, particularly after exposure to secondary stress. The aim of this research was to investigate how the roughness of the downstream purification conditions influences the stability properties of purified IgGs. For this purpose, equine IgGs were extracted from unique hyperimmune plasma by two mild condition-based operational procedures (anion-exchange chromatography and caprylic acid precipitation) and three rougher ones (ammonium sulphate precipitation, cation-exchange chromatography and protein A affinity chromatography). The stability of the refined preparations was studied under non-optimal storage conditions (37 °C, 42 °C, and a transiently lower pH) by monitoring changes in the aggregate content and thermal stability of the pure IgGs. Mild purification protocols generated IgG samples with a lower aggregate share in comparison to the rougher ones. Their tendency for further aggregation was significantly associated with the initial aggregate share. The thermal stability of IgG molecules and the aggregate content in refined samples were inversely correlated. Since the initial proportion of aggregates in the samples was influenced by the operating conditions, we have shown a strong indication that each of them also indirectly affected the stability of the final preparations. This suggests that mild condition-based refinement protocols indeed generate more stable IgGs.

Off-label use of plastic syringes with silicone oil for intravenous infusion bags of antibodies

The formation of particulates in post-manufacture biopharmaceuticals continues to be a major concern in medical treatment. This study was designed to evaluate the content of micro-sized particles using flow imaging of antibodies in intravenous infusion bags. Intravenous immunoglobulin (IVIG) and Avastin® were selected as model drugs and plastic syringes with and without silicone oil (SO) were used to transfer the drugs into the bags (0.9% saline or 5% dextrose). Antibodies exposed to SO had significantly increased levels of microparticles in both diluents, suggesting SO accelerates particle formation, especially at a higher antibody concentration. Even before the drop stress, their count exceeded the USP guideline. Dropping the bags in the presence of SO produced larger microparticles. Meanwhile, air bubbles were retained longer in saline suggesting more protein film formation on its air-water interface. Overall, both drugs were conformationally stable and produced less particles in dextrose than in saline.

Size-based Degradation of Therapeutic Proteins - Mechanisms, Modelling and Control

Protein therapeutics are in great demand due to their effectiveness towards hard-to-treat diseases. Despite their high demand, these bio-therapeutics are very susceptible to degradation via aggregation, fragmentation, oxidation, and reduction, all of which are very likely to affect the quality and efficacy of the product. Mechanisms and modelling of these degradation (aggregation and fragmentation) pathways is critical for gaining a deeper understanding of stability of these products. This review aims to provide a summary of major developments that have occurred towards unravelling the mechanisms of size-based protein degradation (particularly aggregation and fragmentation), modelling of these size-based degradation pathways, and their control. Major caveats that remain in our understanding and control of size-based protein degradation have also been presented and discussed.

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.

Protein-protein interactions at high concentrations: Effects of ArgHCl and NaCl on the stability, viscosity and aggregation mechanisms of protein solution

The aim of this work was to use the diffusion coefficient ration (D_m/D_line) as a parameter to characterize the stability of protein at high concentration, to compare the effects of ArgHCl and NaCl on the interaction of highly concentrated proteins under different pH conditions, and to explore the correlation with protein stability. For this purpose, a high-concentration bovine serum albumin solution (BSA) was selected as the model system, and the diffusion coefficient, aggregation degree, conformational stability, and solution viscosity of the protein were studied by dynamic light scattering (DLS) and spectral detection techniques. The result showed that there was a significant correlation between the D_m/D_line and the protein aggregation. The D_m/D_line of the protein was minimum at pH 7.4, which corresponded to the maximum degree of aggregation and the highest solution viscosity. At pH 7.4, the hydrophobic interactions and the increased conformational stability of ArgHCl maximized the stability of the protein and reduced the viscosity of the solution by 69.3%. At pH 3.0, the strong charge shielding effect of ArgHCl and NaCl and the decreased conformational stability induced protein aggregation and the gel formation. These findings provided valuable insights into the mechanism of protein aggregation and the diffusion coefficient ration (D_m/D_line) could be a potential tool for the pre-formulation studies.

Global multi-method analysis of interaction parameters for reversibly self-associating macromolecules at high concentrations

Weak macromolecular interactions assume a dominant role in the behavior of highly concentrated solutions, and are at the center of a variety of fields ranging from colloidal chemistry to cell biology, neurodegenerative diseases, and manufacturing of protein drugs. They are frequently measured in different biophysical techniques in the form of second virial coefficients, and nonideality coefficients of sedimentation and diffusion, which may be related mechanistically to macromolecular distance distributions in solution and interparticle potentials. A problem arises for proteins where reversible self-association often complicates the concentration-dependent behavior, such that grossly inconsistent coefficients are measured in experiments based on different techniques, confounding quantitative conclusions. Here we present a global multi-method analysis that synergistically bridges gaps in resolution and sensitivity of orthogonal techniques. We demonstrate the method with a panel of monoclonal antibodies exhibiting different degrees of self-association. We show how their concentration-dependent behavior, examined by static and dynamic light scattering and sedimentation velocity, can be jointly described in a self-consistent framework that separates nonideality coefficients from self-association properties, and thereby extends the quantitative interpretation of nonideality coefficients to probe dynamics in highly concentrated protein solutions.

Meat analogs: Protein restructuring during thermomechanical processing

Increasing awareness of inefficient meat production and its future impact on global food security has led the food industry to look for a sustainable approach. Meat products have superior sensorial perception, because of their molecular composition and fibrous structure. Current understanding in the science of food structuring has enabled the utilization of alternative or nonmeat protein ingredients to create novel structured matrices that could resemble the textural functionality of real meat. The physicochemical and structural changes that occur in concentrated protein systems during thermomechanical processing lead to the creation of a fibrous or layered meat-like texture. Phase transitions in concentrated protein systems during protein-protein, protein-polysaccharide, protein-lipid, and protein-water interactions significantly influence the texture and the overall sensory quality of meat analogs. This review summarizes the roles of raw materials (moisture, protein type and concentration, lipids, polysaccharides, and air) and processing parameters (temperature, pH, and shear) in modulating the behavior of the protein phase during the restructuring process (structure-function-process relationship). The big challenge for the food industry is to manufacture concept-based (such as beef-like, chicken-like, etc.) meat analogs with controlled structural attributes. This information will be useful in developing superior meat analogs that fulfill consumer expectations when replacing meat in their diet.

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.

Relationship of PEG-induced precipitation with protein-protein interactions and aggregation rates of high concentration mAb formulations at 5 °C

Native protein-protein interactions can play an important role in determining the tendency of monoclonal antibodies (mAbs) to aggregate under storage conditions. In this context, phase separation of mAb solutions induced by the addition of neutral polymers such as poly(ethylene glycol) (PEG) represents a simple method to assess the tendency of proteins to self-associate in the native state. Here, we investigated their relationships between PEG-induced phase separation, protein-protein interactions and long-term aggregation rate of several formulations of four mAbs at 100 mg/mL and 5 °C over 12 weeks of storage. We observed that the location of the phase boundary correlated well with the osmotic second virial coefficient B₂₂ determined in absence of the polymer, indicating that for our solutions PEG primarily leads to depletion forces between protein molecules, which are additive to protein-protein interactions. However, limited correlation between aggregation rate at 5 °C and phase behavior was observed across different mAbs, pH values and ionic strengths, indicating that colloidal stability is not the only determinant of aggregation even at such low temperature and high protein concentration. Our results contribute to the growing realization that aggregation propensity in the context of antibody developability is a complex feature, which depends on a variety of biophysical properties rather than one single parameter.

Identifying Key Residues That Drive Strong Electrostatic Attractions between Therapeutic Antibodies

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

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.

Light Scattering to Quantify Protein-Protein Interactions at High Protein Concentrations

Static and dynamic (laser) light scattering (SLS and DLS, respectively) can be used to measure the so-called weak or colloidal protein-protein interactions in solution from low to high protein concentrations (c₂). This chapter describes a methodology to measure protein-protein self-interactions using SLS and DLS, with illustrative examples for monoclonal antibody solutions from low to high protein concentrations (c₂ ~ 1-10² g/L).

Isotonic concentrations of excipients control the dimerization rate of a therapeutic immunoglobulin G1 antibody during refrigerated storage based on their rank order of native-state interaction

Inert co-solutes, or excipients, are often included in protein biologic formulations to adjust the tonicity of liquid dosage forms intended for subcutaneous delivery. Despite the low concentration of their use, many of these excipients alter protein-protein interactions such as dimerization and aggregation rates of high concentration monoclonal antibody (mAb) therapeutics to varying extents during long-term refrigerated clinical storage, challenging the formulation scientist to make informed excipient selections at the earliest stages of development when protein supply and time are often limited. The objectives of this study were to better understand how isotonic concentrations of excipients influence the dimerization rates of a model mAb stored at refrigerated and room temperatures and explore protein sparing biophysical methods capable of predicting this dependence. Despite their prevalence of use in the biopharmaceutical industry, methods for assessing conformational stability such differential scanning calorimetry and isothermal equilibrium unfolding showed little predictive power and we highlight some of the assumptions and technical challenges of their use with mAbs. Conversely, measures of colloidal stability of the native-state such as preferential interaction coefficients measured by vapor pressure osmometry and solubility assessed by polyethylene-glycol induced precipitation correlated reasonably well with the mAb dimerization data and are most consistent with the excipients tested minimizing dimerization by interacting favorably with the residues comprising the protein-protein association interface.

Science and art of protein formulation development

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

Predicting High-Concentration Interactions of Monoclonal Antibody Solutions: Comparison of Theoretical Approaches for Strongly Attractive Versus Repulsive Conditions

Nonspecific protein-protein interactions of a monoclonal antibody were quantified experimentally using light scattering from low to high protein concentrations (c₂) and compared with prior work for a different antibody that yielded qualitatively different behavior. The c₂ dependence of the excess Rayleigh ratio (R^ex) provided the osmotic second virial coefficient (B₂₂) at low c₂ and the static structure factor (S_q=0) at high c₂, as a function of solution pH, total ionic strength (TIS), and sucrose concentration. Net repulsive interactions were observed at pH 5, with weaker repulsions at higher TIS. Conversely, attractive electrostatic interactions were observed at pH 6.5, with weaker attractions at higher TIS. Refined coarse-grained models were used to fit model parameters using experimental B₂₂ versus TIS data. The parameters were used to predict high-c₂ R^ex values via Monte Carlo simulations and separately with Mayer-sampling calculations of higher-order virial coefficients. For both methods, predictions for repulsive to mildly attractive conditions were quantitatively accurate. However, only qualitatively accurate predictions were practical for strongly attractive conditions. An alternative, higher resolution model was used to show semiquantitatively and quantitatively accurate predictions of strong electrostatic attractions at low c₂ and low ionic strength.

Chemical and physical instabilities in manufacturing and storage of therapeutic proteins

Development of a robust biologic drug product is accomplished by extensive formulation and process development screening studies; however, even in the most optimal formulation, a protein can undergo spontaneous degradation during manufacture, storage, and clinical use. Chemical changes to amino acid residues, such as oxidation of methionine or tryptophan, or changes in charge such as deamidation or carbonylation, can induce conformational changes in the overall protein structure, potentially leading to changes in physical - in addition to chemical - stability. Oxidation is often caused by light exposure or the presence of metal ions or peroxides. Asparagine deamidation is more likely to occur at higher pH and/or elevated temperature. Mechanical and interfacial stresses during manufacturing can lead to physical instabilities (i.e. various forms of aggregation). A well-defined manufacturing process and effective in-process controls are essential in minimizing chemical and physical instabilities, enabling robust production and distribution of a safe and efficacious drug product. In this work, the authors provide a review of developments in these areas over the past two years, with emphasis on manufacturability of therapeutically relevant proteins and protein-based drug products.

A Reappraisal of Sedimentation Nonideality Coefficients for the Analysis of Weak Interactions of Therapeutic Proteins

The study of weak or colloidal interactions of therapeutic proteins in different formulations allows prediction and optimization of protein stability. Various biophysical techniques have been applied to determine the second osmotic virial coefficient B₂ as it reflects on the macromolecular distance distribution that governs solution behavior at high concentration. In the present work, we exploit a direct link predicted by hydrodynamic theory between B₂ and the nonideality of sedimentation, commonly measured in sedimentation velocity analytical ultracentrifugation through the nonideality coefficient of sedimentation, k_S. Using sedimentation equilibrium analytical ultracentrifugation for independent measurement of B₂, we have examined the dependence of k_S on B₂ for model proteins in different buffers. The data exhibit the expected linear relationship and highlight the impact of protein shape on the magnitude of the nonideality coefficient k_S. Recently, measurements of k_S have been considerably simplified allowing higher throughput and simultaneous polydispersity assessment at higher protein concentrations. Thus, sedimentation velocity may offer a useful approach to compare the impact of formulation conditions on weak interactions and simultaneously on higher-order structure of therapeutic proteins.

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.

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.

Investigation of molecular aggregation mechanism of glipizide/cyclodextrin complexation by combined experimental and molecular modeling approaches

Cyclodextrin complexation is a wise strategy to enhance aqueous solubility of water-insoluble drugs. However, the aggregation mechanism of drug-cyclodextrin complexes is still unclear. This research aimed to investigate the molecular aggregation mechanism of glipizide/cyclodextrin complexation by the combination of experimental and modeling methods. Binding free energies between glipizide and cyclodextrins from modeling calculations were higher than those by the phase solubility diagram method. Both experimental and modeling results showed that methylated-β-cyclodextrin exhibited the best solubilizing capability to glipizide. Size-measurement results confirmed the aggregation between glipizide and all four cyclodextrins in high concentrations. Glipizide/γ-cyclodextrin and glipizide/β-cyclodextrin complexes showed stronger aggregation trend than HP-β-cyclodextrin and methylated-β-cyclodextrin. The substituted groups in the rim of HP-β-cyclodextrin and methylated-β-cyclodextrin lead to weak aggregation. This research provided us a clear molecular mechanism of glipizide/cyclodextrin complexation and aggregation. This research will also benefit the formulation development of cyclodextrin solubilization.

Quantitative NMR study of heat-induced aggregation of eye-lens crystallin proteins under crowding conditions

The eye lens contains a highly concentrated, polydisperse mixture of crystallins, and a loss in transparency during cataract formation is attributed to the aggregation of these proteins. Most biochemical and biophysical studies of crystallins have been performed in diluted samples because of various physical limitations of the respective method at physiological concentrations of up to 200-400 mg/mL. We introduce a straightforward proton NMR transverse relaxometry method to quantify simultaneously proteins in the dissolved and aggregated states at these elevated concentrations, because these states significantly differ in their transverse relaxation properties. The key feature of this method is a direct observation of the protein signal in a wide range of relaxation delays, from few microseconds up to few hundred milliseconds. We applied this method to follow heat-induced aggregation of bovine α- and γB-crystallin between 60 and 200 mg/mL. We find that at 60 °C, a temperature where both crystallins still comprise a native tertiary structure, γB-crystallin aggregated at these high protein concentrations with a time constant of about 30-40 h. α-crystallin remained soluble at 60 mg/mL but formed a transparent gel at 200 mg/mL. This quantitative NMR method can be applied to investigations of other proteins and their mixtures under various aggregation conditions.

Predicting Protein-Protein Interactions of Concentrated Antibody Solutions Using Dilute Solution Data and Coarse-Grained Molecular Models

Protein-protein interactions for solutions of an IgG1 molecule were quantified using static light scattering (SLS) measurements from low to high protein concentrations (c₂). SLS was used to determine second osmotic virial coefficients (B₂₂) at low c₂, and excess Rayleigh profiles (R^ex/K vs. c₂) and zero-q structure factors (S_q=0) as a function of c₂ at higher c₂ for a series of conditions (pH, sucrose concentration, and total ionic strength [TIS]). Repulsive (attractive) interactions were observed at low TIS (high TIS) for pH 5 and 6.5, with increasing repulsions when 5% w/w sucrose was also present. Previously developed and refined coarse-grained antibody models were used to fit model parameters from B₂₂ versus TIS data. The resulting parameters from low-c₂ conditions were used as the sole input to multiprotein Monte Carlo simulations to predict high-c₂R^ex/K and S_q=0 behavior up to 150 g/L. Experimental results at high-c₂ conditions were quantitatively predicted by the simulations for the coarse-grained models that treated antibody molecules as either 6 or 12 (sub) domains, which preserved the basic shape of a monoclonal antibody. Finally, preferential accumulation of sucrose around the protein surface was identified via high-precision density measurements, which self-consistently explained the simulation and experimental SLS results.

Charge Shielding Prevents Aggregation of Supercharged GFP Variants at High Protein Concentration

Understanding protein stability is central to combatting protein aggregation diseases and developing new protein therapeutics. At the high concentrations often present in biological systems, purified proteins can exhibit undesirable high solution viscosities and poor solubilities mediated by short-range electrostatic and hydrophobic protein-protein interactions. The interplay between protein amino acid sequence, protein structure, and solvent conditions to minimize protein-protein interactions is key to designing well-behaved pharmaceutical proteins. However, theoretical approaches have yet to yield a general framework to address these problems. Here, we analyzed the high concentration behavior of superfolder GFP (sfGFP) and two supercharged sfGFP variants engineered to have formal charges of -18 or +15. Under low cosolute conditions, sfGFP and the -18 variant formed a gel or phase separated at ∼10 mg/mL. Under conditions that screen surface charges, including formulations with high histidine or high NaCl concentrations, all three variants attained concentrations up to 250 mg/mL with moderate viscosities. Moreover, all three variants exhibited very similar viscosity-concentration profiles over this range. This effect was not mimicked by high sugar concentrations that exert excluded-volume effects without shielding charge. Collectively, these data demonstrate that charge shielding neutralizes not only long-range electrostatic interactions but also, surprisingly, short-range electrostatic effects due to surface charge anisotropy. This work shows that supercharged sfGFP behavior under high ionic strength is largely determined by particle geometry, a conclusion that is supported by colloid models and may be applicable to pharmaceutically relevant proteins.

Optimization of high-concentration endostatin formulation: Harmonization of excipients' contributions on colloidal and conformational stabilities

Recently, increasing research efforts have been devoted into developing high-concentration protein drugs for subcutaneous injection, especially for those with short half-lives and high-dose requirement. Proteins at high concentrations normally present increased colloidal and structural instability, such as aggregation, fibrillation and gelation, which significantly challenges the high-concentration formulation development of protein drugs. Here we used endostatin, a 20kD recombinant protein, as a model drug for high-concentration formulation optimization. The colloidal and conformational stability of endostatin at high concentration of 30mg/mL were investigated in formulations containing various excipients, including saccharides (mannitol, sorbitol and sucrose), salts (ArgHCl and NaCl), and surfactants (tween 20 and 80). Protein fibrillation was characterized and semi-quantified by optical polarized light microscopy and transmission electron microscopy, and the amount of fiber formation at elevated temperature of 40°C was determined. The soluble protein aggregates were characterized by dynamic and static light scattering before and after dilution. The conformational stability were characterized by polyacrylamide gel electrophoresis, fluorescence, circular dichroism, and differential scanning calorimetry. We observed that the soluble aggregation, fibrillation and gelation, induced by conformational and colloidal instabilities of the protein solution, could be substantially optimized by using suitable stabilizers such as combinations of saccharides and surfactants; while formation of gel and soluble aggregates at high protein concentration (e.g., 30mg/mL) and elevated temperature (40°C) could be prevented by avoiding the usage of salts. It's worth emphasizing that some stabilizers, such as salts and surfactants, could show opposite contributions in conformational and colloidal stabilities of endostatin. Therefore, cautions are needed when one attempts to correlate the colloidal stability of high-concentration proteins with their conformational stability, and the colloidal and conformational protein stabilities must be harmonized by a balanced selection of various types of excipients.

Optimizing the Bioavailability of Subcutaneously Administered Biotherapeutics Through Mechanochemical Drivers

The subcutaneous route offers myriad benefits for the administration of biotherapeutics in both acute and chronic diseases, including convenience, cost effectiveness and the potential for automation through closed-loop systems. Recent advances in parenteral administration devices and the use of additives which enhance drug dispersion have generated substantial additional interest in IV to SQ switching studies. Designing pre-clinical and clinical studies using SQ mediated delivery however requires deep understanding of complex inter-related physiologies and transport pathways governing the interstitial matrix, vascular system and lymphatic channels. This expert review will highlight key structural features which contribute to transport and biodistribution in the subcutaneous space and also assess the impact of drug formulations. Based on the rapidly growing interest in the SQ delivery route, a number of potential areas for future development are highlighted, which are likely to allow continued evolution and innovation in this important area.

Protein Partial Molar Volumes in Multicomponent Solutions from the Perspective of Inverse Kirkwood-Buff Theory

Inverse Kirkwood-Buff (KB) solution theory can be used to relate macroscopic quantities with molecular scale interactions and correlation functions, in the form of KB integrals. Protein partial specific volumes ([Formula: see text]) from high-precision density measurements can be used to quantify solvent-solute and solute-solute KB integrals. Currently, general expressions for [Formula: see text] as a function of cosolute concentration (c₃) have been provided for only binary and ternary solutions. We derive a general multicomponent expression for [Formula: see text] in terms of the relevant KB integrals for the case of low (infinite dilution) protein concentration but arbitrary cosolute concentrations. To test the utility of treating a quaternary system with a pseudoternary approximation, α-Chymotrypsinogen (aCgn) solutions with a series of solutes (NaCl, sucrose, and trehalose) were compared as a function of solute concentration with and without buffer present. Comparison between those ternary and quaternary solutions shows equivalent results within experimental uncertainty and suggests the pseudoternary approximation may be reasonable. In the case of aCgn, doing so also revealed that the preferential interactions can depend on pH. Analysis of steric contributions also provides an example that illustrates how KB integrals allow one to interpret [Formula: see text] in terms of contributions from molecular volume, excluded volume, and hydration/solvation effects.

Novel chemical degradation pathways of proteins mediated by tryptophan oxidation: tryptophan side chain fragmentation

Objectives

This minireview focuses on novel degradation pathways of proteins in solution via intermediary tryptophan (Trp) radical cations, which are generated via photo-induced electron transfer to suitable acceptors such as disulfide bonds.

Methods

Gas-phase mass spectrometry studies had indicated the potential for Trp radical cations to fragment via release of 3-methylene-3H-indol-1-ium from the side chain. HPLC-MS/MS analysis demonstrates that analogous fragmentation reactions occur during the exposure of peptides and proteins to light or accelerated stability testing.

Key findings

The light exposure of selected peptides and monoclonal antibodies leads to the conversion of Trp to glycine (Gly) or glycine hydroperoxide (GlyOOH), where GlyOOH could be reduced to hydroxyglycine, which undergoes subsequent cleavage. Product formation is consistent with Cα–Cβ fragmentation of intermediary Trp radical cations. For the peptide octreotide and specific glycoforms of IgG1 Fc domains, Trp side chain cleavage in aqueous solution is indicated by the formation of 3-methyleneindolenine (3-MEI), which adds to nucleophilic side chains, for example to Lys residues adjacent to the original Trp residues.

Conclusions

Trp side chain cleavage leads to novel reaction products on specific peptide and protein sequences, which may have consequences for potency and immunogenicity.