Nature and consequences of protein–protein interactions in high protein concentration solutions

Protein-protein interactions in dilute to concentrated solutions: α-chymotrypsinogen in acidic conditions

Protein-protein interactions were investigated for α-chymotrypsinogen by static and dynamic light scattering (SLS and DLS, respectively), as well as small-angle neutron scattering (SANS), as a function of protein and salt concentration at acidic conditions. Net protein-protein interactions were probed via the Kirkwood-Buff integral G22 and the static structure factor S(q) from SLS and SANS data. G22 was obtained by regressing the Rayleigh ratio versus protein concentration with a local Taylor series approach, which does not require one to assume the underlying form or nature of intermolecular interactions. In addition, G22 and S(q) were further analyzed by traditional methods involving fits to effective interaction potentials. Although the fitted model parameters were not always physically realistic, the numerical values for G22 and S(q → 0) were in good agreement from SLS and SANS as a function of protein concentration. In the dilute regime, fitted G22 values agreed with those obtained via the osmotic second virial coefficient B22 and showed that electrostatic interactions are the dominant contribution for colloidal interactions in α-chymotrypsinogen solutions. However, as protein concentration increases, the strength of protein-protein interactions decreases, with a more pronounced decrease at low salt concentrations. The results are consistent with an effective "crowding" or excluded volume contribution to G22 due to the long-ranged electrostatic repulsions that are prominent even at the moderate range of protein concentrations used here (<40 g/L). These apparent crowding effects were confirmed and quantified by assessing the hydrodynamic factor H(q → 0), which is obtained by combining measurements of the collective diffusion coefficient from DLS data with measurements of S(q → 0). H(q → 0) was significantly less than that for a corresponding hard-sphere system and showed that hydrodynamic nonidealities can lead to qualitatively incorrect conclusions regarding B22, G22, and static protein-protein interactions if one uses only DLS to assess protein interactions.

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

Small-angle neutron scattering (SANS) is used to probe the solution structure of two protein therapeutics (monoclonal antibodies 1 and 2 (MAb1 and MAb2)) and their protein-protein interaction (PPI) at high concentrations. These MAbs differ by small sequence alterations in the complementarity-determining region but show very large differences in solution viscosity. The analyses of SANS patterns as a function of different solution conditions suggest that the average intramolecular structure of both MAbs in solution is not significantly altered over the studied protein concentrations and experimental conditions. Even though a strong repulsive interaction is expected for both MAbs due to their net charges and low solvent ionic strength, analysis of the SANS data shows that the effective PPI for MAb1 is dominated by a very strong attraction at small volume fraction that becomes negligible at large concentrations. The MAb1 PPI cannot be modeled simply by a spherically symmetric central forces model. It is proposed that an anisotropic attraction strongly affects the local interprotein structure and leads to an anomalously large viscosity of concentrated MAb1 solutions. Conversely, MAb2 displays a repulsive interaction potential throughout the concentration series probed and a comparatively small solution viscosity.

Coarse-grained model for colloidal protein interactions, B(22), and protein cluster formation

Reversible protein cluster formation is an important initial step in the processes of native and non-native protein aggregation, but involves relatively long time and length scales for detailed atomistic simulations and extensive mapping of free energy landscapes. A coarse-grained (CG) model is presented to semiquantitatively characterize the thermodynamics and key configurations involved in the landscape for protein oligomerization, as well as experimental measures of interactions such as the osmotic second virial coefficient (B22). Based on earlier work (Grüenberger et al., J. Phys. Chem. B 2013, 117, 763), this CG model treats proteins as rigid bodies composed of one bead per amino acid, with each amino acid having specific parameters for its size, hydrophobicity, and charge. The net interactions are a combination of steric repulsions, short-range attractions, and screened long-range charge-charge interactions. Model parametrization was done by fitting simulation results against experimental value of B22 as a function of solution ionic strength for α-chymotrypsinogen A and γD-Crystallin (gD-Crys). The CG model is applied to characterize the pairwise interactions and dimerization of gD-Crys and the dependence on temperature, protein concentration, and ionic strength. The results illustrate that at experimentally relevant conditions where stable dimers do not form, the entropic contributions are predominant in the free-energy of protein cluster formation and colloidal protein interactions, arguing against interpretations that treat B22 primarily from energetic considerations alone. Additionally, the results suggest that electrostatic interactions help to modulate the population of the different stable configurations for protein nearest-neighbor pairs, while short-range attractions determine the relative orientations of proteins within these configurations. Finally, simulation results are combined with Principal Component Analysis to identify those amino-acids/surface patches that form interprotein contacts at conditions that favor dimerization of gD-Crys. The resulting regions agree with previously found aggregation-prone sites, as well as suggesting new ones that may be important.

High-throughput biophysical analysis of protein therapeutics to examine interrelationships between aggregate formation and conformational stability

Stabilization and formulation of therapeutic proteins against physical instability, both structural alterations and aggregation, is particularly challenging not only due to each protein's unique physicochemical characteristics but also their susceptibility to the surrounding milieu (pH, ionic strength, excipients, etc.) as well as various environmental stresses (temperature, agitation, lyophilization, etc.). The use of high-throughput techniques can significantly aid in the evaluation of stabilizing solution conditions by permitting a more rapid evaluation of a large matrix of possible combinations. In this mini-review, we discuss both key physical degradation pathways observed for protein-based drugs and the utility of various high-throughput biophysical techniques to aid in protein formulation development to minimize their occurrence. We then focus on four illustrative case studies with therapeutic protein candidates of varying sizes, shapes and physicochemical properties to explore different analytical challenges in monitoring protein physical instability. These include an IgG2 monoclonal antibody, an albumin-fusion protein, a recombinant pentameric plasma glycoprotein, and an antibody fragment (Fab). Future challenges and opportunities to improve and apply high-throughput approaches to protein formulation development are also discussed.

High throughput detection of antibody self-interaction by bio-layer interferometry

Self-interaction of an antibody may lead to aggregation, low solubility or high viscosity. Rapid identification of highly developable leads remains challenging, even though progress has been made with the introduction of techniques such as self-interaction chromatography (SIC) and cross-interaction chromatography (CIC). Here, we report a high throughput method to detect antibody clone self-interaction (CSI) using bio-layer interferometry (BLI) technology. Antibodies with strong self-interaction responses in the CSI-BLI assay also show delayed retention times in SIC and CIC. This method allows hundreds of candidates to be screened in a matter of hours with minimal material consumption.

Application of a PEG precipitation method for solubility screening: a tool for developing high protein concentration formulations

Previous publications demonstrated that the extrapolated solubility by polyethylene glycol (PEG) precipitation method (Middaugh et al., J Biol Chem 1979; 254:367-370; Juckes, Biochim Biophys Acta 1971; 229:535-546; Foster et al., Biochim Biophys Acta 1973; 317:505; Mahadevan and Hall, AIChE J 1990; 36:1517-1528; Stevenson and Hageman, Pharm Res 1995; 12:1671-1676) has a strong correlation to experimentally measured solubility of proteins. Here, we explored the utility of extrapolated solubility as a method to compare multiple protein drug candidates when nonideality of a highly soluble protein prohibits accurate quantitative solubility prediction. To achieve high efficiency and reduce the amount of protein required, the method is miniaturized to microwell plate format for high-throughput screening application. In this simplified version of the method, comparative solubility of proteins can be obtained without the need of concentration measurement of the supernatant following the precipitation step in the conventional method. The monoclonal antibodies with the lowest apparent solubilities determined by this method are the most difficult to be concentrated, indicating a good correlation between the prediction and empirical observations. This study also shows that the PEG precipitation method gives results for opalescence prediction that favorably compares to experimentally determined opalescence levels at high concentration. This approach may be useful in detecting proteins with potential solubility and opalescence problems prior to the time-consuming and expensive development process of high concentration formulation.

Aggregation of biopharmaceuticals in human plasma and human serum: implications for drug research and development

Analytical methods based on light microscopy, 90° light-scattering and surface plasmon resonance (SPR) allowed the characterization of aggregation that can occur when antibodies are mixed with human plasma. Light microscopy showed that aggregates formed when human plasma was mixed with 5% dextrose solutions of Herceptin^® (trastuzumab) or Avastin^® (bevacizumab) but not Remicade^® (infliximab). The aggregates in the plasma-Herceptin^®-5% dextrose solution were globular, size range 0.5–9 μm, with a mean diameter of 4 μm. The aggregates in the plasma-Avastin^®-5% dextrose samples had a mean size of 2 μm. No aggregation was observed when 0.9% NaCl solutions of Herceptin^®, Avastin^® and Remicade^® were mixed with human plasma. 90° light-scattering measurements showed that aggregates were still present 2.5 h after mixing Herceptin^® or Avastin^® with 5% dextrose-plasma solution. A SPR method was utilized to qualitatively describe the extent of interactions of surface-bound antibodies with undiluted human serum. Increased binding was observed in the case of Erbitux^® (cetuximab), whereas no binding was measured for Humira^® (adalimumab). The binding of sera components to 13 monoclonal antibodies was measured and correlated with known serum binding properties of the antibodies. The data presented in this paper provide analytical methods to study the intrinsic and buffer-dependent aggregation tendencies of therapeutic proteins when mixed with human plasma and serum.

Effects of ionic strength and sugars on the aggregation propensity of monoclonal antibodies: influence of colloidal and conformational stabilities

PURPOSE

To develop a general strategy for optimizing monoclonal antibody (MAb) formulations.

METHODS

Colloidal stabilities of four representative MAbs solutions were assessed based on the second virial coefficient (B 2) at 20°C and 40°C, and net charges at different NaCl concentrations, and/or in the presence of sugars. Conformational stabilities were evaluated from the unfolding temperatures. The aggregation propensities were determined at 40°C and after freeze-thawing. The electrostatic potential of antibody surfaces was simulated for the development of rational formulations.

RESULTS

Similar B 2 values were obtained at 20°C and 40°C, implying little dependence on temperature. B 2 correlated quantitatively with aggregation propensities at 40°C. The net charge partly correlated with colloidal stability. Salts stabilized or destabilized MAbs, depending on repulsive or attractive interactions. Sugars improved the aggregation propensity under freeze-thaw stress through improved conformational stability. Uneven and even distributions of potential surfaces were attributed to attractive and strong repulsive electrostatic interactions.

CONCLUSIONS

Assessment of colloidal stability at the lowest ionic strength is particularly effective for the development of formulations. If necessary, salts are added to enhance the colloidal stability. Sugars further improved aggregation propensities by enhancing conformational stability. These behaviors are rationally predictable according to the surface potentials of MAbs.

Capillary viscometer for fully automated measurement of the concentration and shear dependence of the viscosity of macromolecular solutions

The construction and operation of a novel viscometer/rheometer are described. The instrument is designed to measure the viscosity of a macromolecular solution while automatically varying both solute concentration and shear rate. Viscosity is calculated directly from Poiseuille's law, given the measured difference in pressure between two ends of a capillary tube through which the solution is flowing at a known rate. The instrument requires as little as 0.75 mL of a solution to provide a full profile of viscosity as a function of concentration and shear rate, and it can measure viscosities as high as 500 cP and as low as 1 cP, at shear rates between 10 and 2 × 10(3) s(-1). The results of control experiments are presented to document the accuracy and precision of measurement at both low and high concentration of synthetic polymers and proteins.

Kinetically controlled nanostructure formation in self-assembled globular protein-polymer diblock copolymers

Aqueous processing of globular protein-polymer diblock copolymers into solid-state materials and subsequent solvent annealing enables kinetic and thermodynamic control of nanostructure formation to produce block copolymer morphologies that maintain a high degree of protein fold and function. When model diblock copolymers composed of mCherry-b-poly(N-isopropylacrylamide) are used, orthogonal control over solubility of the protein block through changes in pH and the polymer block through changes in temperature is demonstrated during casting and solvent annealing. Hexagonal cylinders, perforated lamellae, lamellae, or hexagonal and disordered micellar phases are observed, depending on the coil fraction of the block copolymer and the kinetic pathway used for self-assembly. Good solvents for the polymer block produce ordered structures reminiscent of coil-coil diblock copolymers, while an unfavorable solvent results in kinetically trapped micellar structures. Decreasing solvent quality for the protein improves long-range ordering, suggesting that the strength of protein interactions influences nanostructure formation. Subsequent solvent annealing results in evolution of the nanostructures, with the best ordering and the highest protein function observed when annealing in a good solvent for both blocks. While protein secondary structure was found to be almost entirely preserved for all processing pathways, UV-vis spectroscopy of solid-state films indicates that using a good solvent for the protein block enables up to 70% of the protein to be retained in its functional form.

Coarse-grained modeling of the self-association of therapeutic monoclonal antibodies

Coarse-grained computational models of two therapeutic monoclonal antibodies are constructed to understand the effect of domain-level charge-charge electrostatics on the self-association phenomena at high protein concentrations. The coarse-grained representations of the individual antibodies are constructed using an elastic network normal-mode analysis. Two different models are constructed for each antibody for a compact Y-shaped and an extended Y-shaped configuration. The resulting simulations of these coarse-grained antibodies that interact through screened electrostatics are done at six different concentrations. It is observed that a particular monoclonal antibody (hereafter referred to as MAb1) forms three-dimensional heterogeneous structures with dense regions or clusters compared to a different monoclonal antibody (hereafter referred to as MAb2) that forms more homogeneous structures (no clusters). These structures, together with the potential mean force (PMF) and radial distribution functions (RDF) between pairs of coarse-grained regions on the MAbs, are qualitatively consistent with the experimental observation that MAb1 has a significantly higher viscosity compared to MAb2, especially at concentrations >50 mg/mL, even though the only difference between the MAbs lies with a few amino acids at the antigen-binding loops (CDRs). It is also observed that the structures in MAb1 are formed due to stronger Fab-Fab interactions in corroboration with experimental observations. Evidence is also shown that Fab-Fc interactions can be equally important in addition to Fab-Fab interactions. The coarse-grained representations are effective in picking up differences based on local charge distributions of domains and make predictions on the self-association characteristics of these protein solutions. This is the first computational study of its kind to show that there are differences in structures formed by two different monoclonal antibodies at high concentrations.

Antibody nanoparticle dispersions formed with mixtures of crowding molecules retain activity and in vivo bioavailability

Monoclonal antibodies continue to command a large market for treatment of a variety of diseases. In many cases, the doses required for therapeutic efficacy are large, limiting options for antibody delivery and administration. We report a novel formulation strategy based on dispersions of antibody nanoclusters that allows for subcutaneous injection of highly concentrated antibody (≈ 190 mg/mL). A solution of monoclonal antibody 1B7 was rapidly frozen and lyophilized using a novel spiral-wound in-situ freezing technology to generate amorphous particles. Upon gentle stirring, a translucent dispersion of approximately 430 nm protein clusters with low apparent viscosity (≈ 24 cp) formed rapidly in buffer containing the pharmaceutically acceptable crowding agents such as trehalose, polyethylene glycol, and n-methyl-2-pyrrolidone. Upon in vitro dilution of the dispersion, the nanoclusters rapidly reverted to monomeric protein with full activity, as monitored by dynamic light scattering and antigen binding. When administered to mice as an intravenous solution, subcutaneous solution, or subcutaneous dispersion at similar (4.6-7.3 mg/kg) or ultra-high dosages (51.6 mg/kg), the distribution and elimination kinetics were within error and the protein retained full activity. Overall, this method of generating high-concentration, low-viscosity dispersions of antibody nanoclusters could lead to improved administration and patient compliance, providing new opportunities for the biotechnology industry.

High-dose monoclonal antibodies via the subcutaneous route: challenges and technical solutions, an industry perspective

This review summarizes the various challenges in product development involved in subcutaneous administration of high-dose monoclonal antibodies and attempts to provide an industry perspective of some of the available technologies and potential avenues to overcome these challenges.

Advancements in high throughput biophysical technologies: applications for characterization and screening during early formulation development of monoclonal antibodies

The formulation development of monoclonal antibodies is extremely challenging, due to the diversity and complexity contained within this class of molecules. The physical and chemical properties of a monoclonal antibody dictate the behavior of the protein drug during manufacturing, storage and clinical administration. In the past few years, the use of high throughput technologies has been widely adapted to delineate unique properties of individual immunoglobulin G's (IgG's) important for their development. Numerous screening techniques have been designed to reveal physical and chemical characteristics of a protein relevant to stability under production, formulation and delivery conditions. In addition, protein stability under accelerated stresses has been utilized to predict long-term storage behavior for monoclonal antibodies in the formulation. In this review, we summarize the recent advancements in the field of biophysical technology, with a specific focus on the techniques that can be directly applied to the formulation development of monoclonal antibodies. Several case studies are also presented here to provide examples of combining existing biophysical methods with high throughput screening technology in the formulation development of monoclonal antibody drugs.

Determination of the dipole moments of RNAse SA wild type and a basic mutant

In this study, we report the effects of acidic to basic residue point mutations (5K) on the dipole moment of RNAse SA at different pHs. Dipole moments were determined by measuring solution capacitance of the wild type (WT) and the 5K mutant with an impedance analyzer. The dipole moments were then (1) compared with theoretically calculated dipole moments, (2) analyzed to determine the effect of the point mutations, and (3) analyzed for their contribution to overall protein-protein interactions (PPI) in solution as quantitated by experimentally derived second virial coefficients. We determined that experimental and calculated dipoles were in reasonable agreement. Differences are likely due to local motions of residue side chains, which are not accounted for by the calculated dipole. We observed that the proteins' dipole moments increase as the pH is shifted further from their isoelectric points and that the wild-type dipole moments were greater than those of the 5K. This is likely due to an increase in the proportion of one charge (either negative or positive) relative to the other. A greater charge disparity corresponded to a larger dipole moment. Finally, the larger dipole moments of the WT resulted in greater attractive overall PPI for that protein as compared to the 5K.

Viscosity of concentrated therapeutic protein compositions

The use of monoclonal antibodies as therapeutic agents has been increasing steadily over the last decade for the treatment of various conditions. There is often a need to deliver a large dose of the protein, so there is a trend toward developing commercially viable liquid formulations of highly concentrated antibodies. Such concentrated solutions are associated with a number of challenges, including optimization of production processes, plus chemical and physical stability of the final product where solution viscosity becomes a critical quality attribute. Assessment of the rheological characteristics of concentrated compositions is essential as are development strategies to reduce the viscosity. This review covers the state-of-the-art rheology measurement techniques, focusing particularly on concentrated protein solutions. Current understanding of the mechanisms leading to high viscosity and control by formulation parameters is discussed.

Rheological and syringeability properties of highly concentrated human polyclonal immunoglobulin solutions

This study of highly concentrated polyvalent immunoglobulin solutions, IgG, aimed at analyzing the relationships between protein concentration and aggregation on the one hand and viscosity on the other hand. Viscosity variations as a function of IgG concentration showed two well-defined behaviours: a Newtonian behaviour for low-concentrated solutions and a shear-thinning behaviour for highly concentrated ones. The viscosity data fitted very well with the Mooney model, suggesting the absence of intermolecular interactions in the IgG solutions that behaved like a non-interacting suspension of hard particles. The polyclonal nature of IgG seems to prevent intermolecular interaction. The shape factor, determined from Mooney fitting, revealed a non-spherical shape of the polyclonal IgG molecules. The rheological properties were also correlated with the injection force (F) through hypodermic needles by syringeability tests. Here, F was mainly affected by three parameters: the solution viscosity, the injection flow rate, and the needle characteristics. In fact, syringeability tests showed that F increased with IgG concentration and flow rate and decreased with the internal diameter of the needle. A zone for optimal injection conditions was then identified taking into account the different affecting parameters and mainly a maximum force for manual injection, which was fixed at 30N.

Low viscosity highly concentrated injectable nonaqueous suspensions of lysozyme microparticles

Subcutaneous injection of concentrated protein and peptide solutions, in the range of 100-400 mg/mL, is often not possible with a 25- to 27-gauge needle, as the viscosity can be well above 50 cP. Apparent viscosities below this limit are reported for suspensions of milled lysozyme microparticles up to nearly 400 mg/mL in benzyl benzoate or benzyl benzoate mixtures with safflower oils through a syringe with a 25- to 27-gauge needle at room temperature. These apparent viscosities were confirmed using a cone-and-plate rheometer. The intrinsic viscosity regressed from the Kreiger-Dougherty model was only slightly above the Einstein value of 2.5, indicating the increase in viscosity relative to that of the solvent was caused primarily by the excluded volume. Thus, the increases in viscosity from electrical double layer interactions (electroviscous effects), solvation of the particles, or deviations of the particle shape from a spherical geometry were minimal, and much smaller than typically observed for proteins dissolved in aqueous solutions. The small electroviscous effects are expected given the negligible zeta potential and thin double layers in the low dielectric constant organic solvent. The suspensions were resuspendable after a year, with essentially constant particle size after two months as measured by static light scattering. The lower apparent viscosities for highly concentrated protein suspensions relative to protein solutions, coupled with these favorable characteristics upon resuspension, may offer novel opportunities for subcutaneous injection of therapeutic proteins.

Long- and short-range electrostatic interactions affect the rheology of highly concentrated antibody solutions

PURPOSE

To explain the differences in protein-protein interactions (PPI) of concentrated versus dilute formulations of a model antibody.

METHODS

High frequency rheological measurements from pH 3.0 to 12.0 quantitated viscoelasticity and PPI at high concentrations. Dynamic light scattering (DLS) characterized PPI in dilute solutions.

RESULTS

For concentrated solutions at low ionic strength, the storage modulus, a viscosity component and a measure of PPI, is highest at the isoelectric point (pH 9.0) and lowest at pH 5.4. This profile flattens at higher ionic strength but not completely, indicating PPI consist of long-range electrostatics and other short-range attractions. At low concentrations, PPI are near zero at pI but become repulsive as the pH is shifted. Higher salt concentrations completely flatten this profile to zero, indicating that these PPI are mainly electrostatic.

CONCLUSIONS

This discrepancy occurs because long-range interactions are significant at low concentrations, whereas both long- and short-range interactions are significant at higher concentrations. Computer modeling was used to calculate antibody properties responsible for long- and short-range interactions, i.e. net charge and dipole moment. Charge-charge interactions are repulsive while dipole-dipole interactions are attractive. Their net effect correlated with the storage modulus profile. However, only charge-charge repulsions correlated with PPI determined by DLS.