Challenges in Predicting Protein-Protein Interactions from Measurements of Molecular Diffusivity

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

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

Understanding the Phase Behavior of a Multistimuli-Responsive Elastin-like Polymer: Insights from Dynamic Light Scattering Analysis

Elastin-like polymers are a class of stimuli-responsive protein polymers that hold immense promise in applications such as drug delivery, hydrogels, and biosensors. Yet, understanding the intricate interplay of factors influencing their stimuli-responsive behavior remains a challenging frontier. Using temperature-controlled dynamic light scattering and zeta potential measurements, we investigate the interactions between buffer, pH, salt, water, and protein using an elastin-like polymer containing ionizable lysine residues. We observed the elevation of transition temperature in the presence of the common buffering agent HEPES at low concentrations, suggesting a "salting-in" effect of HEPES as a cosolute through weak association with the protein. Our findings motivate a more comprehensive investigation of the influence of buffer and other cosolute molecules on elastin-like polymer behavior.

Investigation on the Combined Effect of Hydroxypropyl Beta-Cyclodextrin (HPβCD) and Polysorbate in Monoclonal Antibody Formulation

Monoclonal antibodies require careful formulation due to their inherent stability limitations. Polysorbates are commonly used to stabilize mAbs, but they are prone to degradation, which results in unwanted impurities. KLEPTOSE® HPβCD (hydroxypropyl beta-cyclodextrin) has functioned as a stable stabilizer for protein formulations in our previous research. The current study investigates the collaborative impact of combining polysorbates and HPβCD as excipients in protein formulations. The introduction of HPβCD in formulations showed it considerably reduced aggregation in two model proteins, bevacizumab and ipilimumab, following exposure to various stress conditions. The diffusion interaction parameter revealed a reduction in protein-protein interactions by HPβCD. In bevacizumab formulations, the subvisible particle counts per 0.4 mL of samples in commercial formulations vs. formulations containing both HPβCD and polysorbates subjected to distinct stressors were as follows: agitation, 87,308 particles vs. 15,350 particles; light, 25,492 particles vs. 6765 particles; and heat, 1775 particles vs. 460 particles. Isothermal titration calorimetry (ITC) measurement indicated a weak interaction between PS 80 and HPβCD, with a KD value of 74.7 ± 7.5 µM and binding sites of 5 × 10-3. Surface tension measurements illustrated that HPβCD enhanced the surface activity of polysorbates. The study suggests that combining these excipients can improve mAb stability in formulations, offering an alternative for the biopharmaceutical industry.

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

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

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

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

Comprehensive physicochemical and functional analysis of pembrolizumab based on controlled degradation studies: Impact on antigen-antibody binding

Monoclonal antibodies-based medicines are widely used in the treatment of different diseases. These medicines are very sensitive to exposure to different environmental conditions and their handling in hospitals may affect their safety and efficacy. This is the case for pembrolizumab (Keytruda®, 25 mg/mL), for which there is not yet much information on its risk behaviour associated with routine handling or unintentional mishandling. Here we performed a wider physicochemical and functional analysis of pembrolizumab medicine including controlled degradation studies: heat, freeze/thaw, agitation, accelerated light exposure and high hypertonic solution. After that, the samples were analysed by a set of analytical techniques to evaluated critical quality attributes: Far-UV CD, IT-FS, DLS, RP/UHPLC(UV)-MS, SE/UHPLC(UV), RP/UHPLC(UV)-MS/MS and ELISA. The results provide an in-depth understanding of the biochemical and biophysical properties of pembrolizumab, showing that the medicine is affected by accelerated light exposure and temperature of 60 °C, demonstrated by the detection of non-natural dimers and HMWS. Light exposure also revealed different isoform profile and increase in oxidations. Regarding functionality by means of the interaction antigen-antibody binding, all the stressors promoted a decrease in pembrolizumab capacity to bind to PD-1 receptor, although the biological activity remained still high for all of them, being 60 °C and accelerated light exposure the most affected.

A comprehensive evaluation of arginine and its derivatives as protein formulation stabilizers

Arginine and its derivatives (such as arginine ethyl ester and acetyl arginine) have varying degrees of protein aggregation suppressor effect across different protein solutions. To understand this performance ambiguity, we evaluated the activity of arginine, acetyl arginine, and arginine ethyl ester for aggregation suppressor effect against human intravenous immunoglobulin G (IgG) solution at pH 4.8. Both arginine and its cationic derivative arginine ethyl ester in their hydrochloride salt forms significantly reduced the colloidal and conformational stability (reduced kd and Tm) of IgG. Consequently, the monomer content was decreased with an increase in subvisible particulates after agitation or thermal stress. Furthermore, compared to arginine, arginine ethyl ester with one more cationic charge and hydrochloride salt form readily precipitated IgG at temperatures higher than 25 °C. On the contrary, acetyl arginine, which mostly exists in a neutral state at pH 4.8, efficiently suppressed the formation of subvisible particles retaining a high amount of monomer owing to its higher colloidal and conformational stability. Concisely, the charged state of additives significantly impacts protein stability. This study demonstrated that contrary to popular belief, arginine and its derivatives may either enhance or suppress protein aggregation depending on their net charge and concentration.

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

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

The effect of mAb and excipient cryoconcentration on long-term frozen storage stability – part 2: Aggregate formation and oxidation

We examined the impact of monoclonal antibody (mAb) and buffer concentration, mimicking the cryoconcentration found upon freezing in a 2 L bottle, on mAb stability during frozen storage. Upon cryoconcentration, larger protein molecules and small excipient molecules freeze-concentrate differently, resulting in different protein to stabiliser ratios within a container. Understanding the impact of these shifted ratios on protein stability is essential. For two mAbs a set of samples with constant mAb (5 mg/mL) or buffer concentration (medium histidine/adipic acid) was prepared and stored for 6 months at −10 °C. Stability was evaluated via size-exclusion chromatography, flow imaging microscopy, UV/Vis spectroscopy at 350 nm, and protein A chromatography. Dynamic light scattering was used to determine k_D values. Soluble aggregate levels were unaffected by mAb concentration, but increased with histidine concentration. No trend in optical density could be identified. In contrast, increasing mAb or buffer concentration facilitated the formation of subvisible particles. A trend towards attractive protein-protein interactions was seen with higher ionic strength. MAb oxidation levels were negatively affected by increasing histidine concentration, but became less with higher mAb concentration. Small changes in mAb and buffer composition had a significant impact on stability during six-month frozen storage. Thus, preventing cryoconcentration effects in larger freezing containers may improve long-term stability.

Proteins on the rack: Mechanistic studies on protein particle formation during peristaltic pumping

Peristaltic pumping can cause protein particle formation. The expected causes were unfolding by heat in the pump head, oxidative stress by cavitation generated during roller movement, interfacial adsorption to the tubing wall and mechanical stress by stretching of the tubing itself. The pump head reached 28°C during experiments stayed well below the onset of the melting point of the proteins. Thus, heat may only be a relevant root cause for proteins containing domains with very low unfolding temperature. Analysis by terephthalic acid dosimetry and protein oxidation via RP-HPLC ruled out major induction of reactive hydroxyl radicals by pumping, indicating that cavitation does not play a significant role in particle generation. Addition of surfactants suppresses protein adsorption to the tubing wall and drastically reduced protein particle formation. This indicates that interfacial protein adsorption is a key element. Repeated stretching of tubing filled with protein solution led to the formation of protein particles, demonstrating that expansion and compression of the protein film on the tubing surface is the second key component for particle formation. Thus, protein particle generation during peristaltic pumping originates from the formation of a protein film on the tubing surface which gets stretched and compressed, leading to film fragments entering the bulk solution. This interplay of protein film formation and its rupture has been also observed at liquid/liquid or liquid/air interfaces.

Dipolar Interactions and Protein Hydration in Highly Concentrated Antibody Formulations

Molecular interaction mechanisms in high-concentrated protein systems are of fundamental importance for the rational development of biopharmaceuticals such as monoclonal antibody (mAb) formulations. In such high-concentrated protein systems, the intermolecular distances between mAb molecules are reduced to the size of the protein diameter (approx. 10 nm). Thus, protein-protein interactions are more pronounced at high concentrations; so a direct extrapolation of physicochemical properties obtained from measurements at a low protein concentration of the corresponding properties at a high protein concentration is highly questionable. Besides the charge-charge interaction, the effects of molecular crowding, dipolar interaction, changes in protein hydration, and self-assembling tendency become more relevant. Here, protein hydration, protein dipole moment, and protein-protein interactions were studied in protein concentrations up to 200 mg/mL (= 1.3 mM) in different formulations for selected mAbs using dielectric relaxation spectroscopy (DRS). These data are correlated with the second virial coefficient, A₂, the diffusion interaction parameter, k_D, the elastic shear modulus, G', and the dynamic viscosity, η. When large contributions of dipolar protein-protein interactions were observed, the tendency of self-assembling and an increase in solution viscosity were detected. These effects were examined using specific buffer conditions. Furthermore, different types of protein-water interactions were identified via DRS, whereby the effect of high protein concentration on protein hydration was investigated for different high-concentrated liquid formulations (HCLFs).

Revisit PEG-Induced Precipitation Assay for Protein Solubility Assessment of Monoclonal Antibody Formulations

PURPOSE

Protein solubility is an important attribute of pharmaceutical monoclonal antibody (MAb) formulations, particularly at high MAb concentrations. PEG-induced protein precipitation has been routinely used to assess protein solubility. To provide insights for better understanding and implementation of PEG-induced protein precipitation assay, this work compares different solubility measures and examines their relevance to loss of protein solubility in concentrated formulations.

METHODS

Solubility of a MAb in 15 formulations was evaluated using PEG-induced precipitation assay. Three apparent protein solubility measures, the middle-point and onset PEG concentrations (c_mid and c_onset) as well as the binding free energy (μ_B), were obtained from the PEG-induced protein precipitation assay and compared to the DLS protein interaction parameter (k_D). Visual inspection of loss of protein solubility in concentrated formulations during storage was used to further examine the discrepancy of protein solubility ranking by these measures.

RESULTS

PEG-induced precipitation assay predicted overall protein solubility ranking similar to that by DLS k_D. However, for three formulations with ionic excipients NaCl, Arg·Cl, and Arg·Glu·Cl, PEG-induced precipitation assay yielded more accurate predictions compared to DLS k_D measurements. Furthermore, μ_B showed superior ability in distinguishing protein solubility for these formulations.

CONCLUSIONS

This study demonstrated good correlations between the protein solubility measures obtained from PEG-induced precipitation experiments and DLS k_D measurement. It also provides one example in which protein solubility ranking by binding free energy is more accurate than the other measures. The results support the theoretical proposition that μ_B has a potential to serve as standard protein solubility measure.

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.

High-Throughput Screening for Colloidal Stability of Peptide Formulations Using Dynamic and Static Light Scattering

Selection of an appropriate formulation to stabilize therapeutic proteins against aggregation is one of the most challenging tasks in early-stage drug product development. The amount of aggregates is more difficult to quantify in the case of peptides due to their small molecular size. Here, we investigated the suitability of diffusion self-interaction parameters (k_D) and osmotic second virial coefficients (B₂₂) for high-throughput (HT) screening of peptide formulations regarding their aggregation risk. These parameters were compared to the effect of thermal stress on colloidal stability. The formulation matrix comprised six buffering systems at two selected pH values, four tonicity agents, and a common preservative. The results revealed that electrostatic interactions are the main driver to control colloidal stability. Preferred formulations consisted of acetate and succinate buffer at pH 4.5 combined with glycerol or mannitol and optional m-cresol. k_D proved to be a suitable surrogate for B₂₂ as an indicator of high colloidal stability in the case of peptides as was previously described for globular proteins and antibodies. Formulation assessment solely based on k_D obtained by HT methods offers important insights into the optimization of colloidal stability during the early development of peptide-based liquid formulations and can be performed with a limited amount of peptide (∼360 mg).

Energetic Dissection of Mab-Specific Reversible Self-Association Reveals Unique Thermodynamic Signatures

PURPOSE

Reversible self-association (RSA) remains a challenge in the development of therapeutic monoclonal antibodies (mAbs). We recently analyzed the energetics of RSA for five IgG mAbs (designated as A-E) under matched conditions and using orthogonal methods. Here we examine the thermodynamics of RSA for two of the mAbs that showed the strongest evidence of RSA (mAbs C and E) to identify underlying mechanisms.

METHODS

Concentration-dependent dynamic light scattering and sedimentation velocity (SV) studies were carried out for each mAb over a range of temperatures. Because self-association was weak, the SV data were globally analyzed via direct boundary fitting to identify best-fit models, accurately determine interaction energetics, and account for the confounding effects of thermodynamic and hydrodynamic nonideality.

RESULTS

mAb C undergoes isodesmic self-association at all temperatures examined, with the energetics indicative of an enthalpically-driven reaction offset by a significant entropic penalty. By contrast, mAb E undergoes monomer-dimer self-association, with the reaction being entropically-driven and comprised of only a small enthalpic contribution.

CONCLUSIONS

Classical interpretations implicate van der Waals interactions and H-bond formation for mAb C RSA, and electrostatic interactions for mAb E. However, noting that RSA is likely coupled to additional equilibria, we also discuss the limitations of such interpretations.

Evaluation of biological activities, structural and conformational properties of bovine beta- and alpha-trypsin isoforms in aqueous-organic media

The study of the biological activity of trypsin isoforms in aqueous-organic media is of great interest to various fields of knowledge and biochemistry applications. Thus enzymatic, structural, and energetic properties of bovine β- and α-trypsin isoforms were compared in aqueous-organic media using 30 mg of each isoform. The results showed that the changes induced on the structure and activity of the same trypsin isoform occur at different concentrations. Better results for activity (ionic strength of 0.11 mol·L^-1, at 37 °C and pH 8.0) were found in 0-40% of ethanolic media in which the activity for β-trypsin was about 60% higher than ɑ-trypsin. The ethanolic system does not cause significant changes in the level of secondary structure but the β-trypsin isoform undergoes a major rearrangement. The use of until 60% (v/v) ethanol showed that β-trypsin presents a denaturation process 17% more cooperative. The organic solvent causes redistribution in the supramolecular arrangement of both isoforms: all concentrations used induced the β-trypsin molecules to rearrange into agglomerates. The ɑ-trypsin rearranges into agglomerates up to 60% (v/v) of ethanol and aggregates at 80% (v/v) of ethanol. Both isoforms keep the enzymatic activity up to 60% (v/v) of ethanol.

Monoclonal Antibody Aggregation near Silicone Oil-Water Interfaces

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

Microwell Plate-Based Dynamic Light Scattering as a High-Throughput Characterization Tool in Biopharmaceutical Development

High-throughput light scattering instruments are widely used in screening of biopharmaceutical formulations and can be easily incorporated into processes by utilizing multi-well plate formats. High-throughput plate readers are helpful tools to assess the aggregation tendency and colloidal stability of biological drug candidates based on the diffusion self-interaction parameter (k_D). However, plate readers evoke issues about the precision and variability of determined data. In this article, we report about the statistical evaluation of intra- and inter-plate variability (384-well plates) for the k_D analysis of protein and peptide solutions. ANOVA revealed no significant differences between the runs. In conclusion, the reliability and precision of k_D was dependent on the plate position of the sample replicates and k_D value. Positive k_D values (57.0 mL/g, coefficients of variation (CV) 8.9%) showed a lower variability compared to negative k_D values (−14.8 mL/g, CV 13.4%). The variability of k_D was not reduced using more data points (120 vs. 30). A k_D analysis exclusively based on center wells showed a lower CV (<2%) compared to edge wells (5–12%) or a combination of edge and center wells (2–5%). We present plate designs for k_D analysis within the early formulation development, screening up to 20 formulations consuming less than 50 mg of active pharmaceutical ingredient (API).

Computational Characterization of Antibody-Excipient Interactions for Rational Excipient Selection Using the Site Identification by Ligand Competitive Saturation-Biologics Approach

Protein therapeutics typically require a concentrated protein formulation, which can lead to self-association and/or high viscosity due to protein-protein interaction (PPI). Excipients are often added to improve stability, bioavailability, and manufacturability of the protein therapeutics, but the selection of excipients often relies on trial and error. Therefore, understanding the excipient-protein interaction and its effect on non-specific PPI is important for rational selection of formulation development. In this study, we validate a general workflow based on the site identification by ligand competitive saturation (SILCS) technology, termed SILCS-Biologics, that can be applied to protein therapeutics for rational excipient selection. The National Institute of Standards and Technology monoclonal antibody (NISTmAb) reference along with the CNTO607 mAb is used as model antibody proteins to examine PPIs, and NISTmAb was used to further examine excipient-protein interactions, in silico. Metrics from SILCS include the distribution and predicted affinity of excipients, buffer interactions with the NISTmAb Fab, and the relation of the interactions to predicted PPI. Comparison with a range of experimental data showed multiple SILCS metrics to be predictive. Specifically, the number of favorable sites to which an excipient binds and the number of sites to which an excipient binds that are involved in predicted PPIs correlate with the experimentally determined viscosity. In addition, a combination of the number of binding sites and the predicted binding affinity is indicated to be predictive of relative protein stability. Comparison of arginine, trehalose, and sucrose, all of which give the highest viscosity in combination with analysis of B₂₂ and k_D and the SILCS metrics, indicates that higher viscosities are associated with a low number of predicted binding sites, with lower binding affinity of arginine leading to its anomalously high impact on viscosity. The present study indicates the potential for the SILCS-Biologics approach to be of utility in the rational design of excipients during biologics formulation.

N-Acetylated-L-arginine (NALA) is an enhanced protein aggregation suppressor under interfacial stresses and elevated temperature for protein liquid formulations

Even though arginine hydrochloride has been recognized as a protein aggregation suppressor in the biopharmaceutical industry, its use has been questioned due to decreasing transition unfolding temperatures (T_m). Four compounds were designed to enhance the role of arginine by changing the length of the carbon chain with removal or N-acetylation of α-amino group. Biophysical properties were observed by differential scanning calorimetry (DSC), dynamic light scattering (DLS), size-exclusion chromatography (SEC), and flow imaging (FI). N-Acetyl-L-arginine (NALA) performed the best at minimizing decrease in T_m with arginine at different pH. NALA also demonstrated relatively higher colloidal stability than arginine hydrochloride, especially in the acidic pH, thereby reducing agitation stress of IgG. Moreover, NALA exhibited a cooperative effect with commercially used glycine buffer for IVIG to maintain the monomer contents with almost no change and suppressed larger particle formation after agitation with heat. The study concludes that the decreasing T_m of proteins by arginine hydrochloride is due to amide group in the α-carbon chain. Moreover, chemical modification on the group compared to removing it will be a breakthrough of arginine's limitations and optimize storage stability of protein therapeutics.

Quantifying Protein-Protein Interactions by Molecular Counting with Mass Photometry

Interactions between biomolecules control the processes of life in health and their malfunction in disease, making their characterization and quantification essential. Immobilization- and label-free analytical techniques are desirable because of their simplicity and minimal invasiveness, but they struggle with quantifying tight interactions. Here, we show that mass photometry can accurately count, distinguish by molecular mass, and thereby reveal the relative abundances of different unlabelled biomolecules and their complexes in mixtures at the single-molecule level. These measurements determine binding affinities over four orders of magnitude at equilibrium for both simple and complex stoichiometries within minutes, as well as the associated kinetics. These results introduce mass photometry as a rapid, simple and label-free method for studying sub-micromolar binding affinities, with potential for extension towards a universal approach for characterizing complex biomolecular interactions.

On the Aggregation and Nucleation Mechanism of the Monoclonal Antibody Anti-CD20 Near Liquid-Liquid Phase Separation (LLPS)

The crystallization of Anti-CD20, a full-length monoclonal antibody, has been studied in the PEG400/Na2SO4/Water system near Liquid-Liquid Phase Separation (LLPS) conditions by both sitting-drop vapour diffusion and batch methods. In order to understand the Anti-CD20 crystallization propensity in the solvent system of different compositions, we investigated some measurable parameters, normally used to assess protein conformational and colloidal stability in solution, with the aim to understand the aggregation mechanism of this complex biomacromolecule. We propose that under crystallization conditions a minor population of specifically aggregated protein molecules are present. While this minor species hardly contributes to the measured average solution behaviour, it induces and promotes crystal formation. The existence of this minor species is the result of the LLPS occurring concomitantly under crystallization conditions.

Amplification of Salt-Induced Protein Diffusiophoresis by Varying Salt from Potassium to Sodium to Magnesium Chloride in Water

Salt-induced diffusiophoresis is the migration of a macromolecule or a colloidal particle induced by a concentration gradient of salt in water. Here, the effect of salt type on salt-induced diffusiophoresis of the protein lysozyme at pH 4.5 and 25 °C was examined as a function of salt concentration for three chloride salts: NaCl, KCl, and MgCl₂. Diffusiophoresis coefficients were calculated from experimental ternary diffusion coefficients on lysozyme-salt-water mixtures. In all cases, diffusiophoresis of this positively charged protein occurs from high to low salt concentration. An appropriate mass transfer process was theoretically examined to show that concentration gradients of MgCl₂ produce significant lysozyme diffusiophoresis. This is attributed to the relatively low mobility of Mg²⁺ ions compared to that of Cl^- ions at low salt concentration and a strong thermodynamic nonideality of this salt at high salt concentration. These findings indicate that MgCl₂ concentration gradients could be exploited for protein manipulation in solution (e.g., using microfluidic technologies) with applications to protein adsorption and purification. The dependence of lysozyme diffusiophoresis on salt type was theoretically examined and linked to protein charge. The effect of salts on hydrogen-ion titration curves was experimentally characterized to understand the role of salt type on protein charge. Our results indicate that binding of Mg²⁺ ions to lysozyme further enhances protein diffusiophoresis.

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.

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.

The Investigation of Protein Diffusion via H-Cell Microfluidics

In this study, we developed a microfluidics method, using a so-called H-cell microfluidics device, for the determination of protein diffusion coefficients at different concentrations, pHs, ionic strengths, and solvent viscosities. Protein transfer takes place in the H-cell channels between two laminarly flowing streams with each containing a different initial protein concentration. The protein diffusion coefficients are calculated based on the measured protein mass transfer, the channel dimensions, and the contact time between the two streams. The diffusion rates of lysozyme, cytochrome c, myoglobin, ovalbumin, bovine serum albumin, and etanercept were investigated. The accuracy of the presented methodology was demonstrated by comparing the measured diffusion coefficients with literature values measured under similar solvent conditions using other techniques. At low pH and ionic strength, the measured lysozyme diffusion coefficient increased with the protein concentration gradient, suggesting stronger and more frequent intermolecular interactions. At comparable concentration gradients, the measured lysozyme diffusion coefficient decreased drastically as a function of increasing ionic strength (from zero onwards) and increasing medium viscosity. Additionally, a particle tracing numerical simulation was performed to achieve a better understanding of the macromolecular displacement in the H-cell microchannels. It was found that particle transfer between the two channels tends to speed up at low ionic strength and high concentration gradient. This confirms the corresponding experimental observation of protein diffusion measured via the H-cell microfluidics.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.

Hydrodynamic radius coincides with the slip plane position in the electrokinetic behavior of lysozyme

The zeta potential (ζ) is the effective charge energy of a solvated protein, describing the magnitude of electrostatic interactions in solution. It is commonly used in the assessment of adsorption processes and dispersion stability. Predicting ζ from molecular structure would be useful to the structure-based molecular design of drugs, proteins, and other molecules that hold charge-dependent function while remaining suspended in solution. One challenge in predicting ζ is identifying the location of the slip plane (X_SP ), a distance from the protein surface where ζ is theoretically defined. This study tests the hypothesis that the X_SP can be estimated by the Stokes-Einstein hydrodynamic radius (R_h ), using globular hen egg white lysozyme as a model system. Although the X_SP and R_h differ in their theoretical definitions, with the X_SP being the position of the ζ during electrokinetic phenomena (e.g., electrophoresis) and the R_h being a radius pertaining to the edge of solvation during diffusion, they both represent the point where water and ions no longer adhere to a molecule. This work identifies the limited range of ionic strengths in which the X_SP can be determined using diffusivity measurements and the Stokes-Einstein equation. In addition, a computational protocol is developed for determining the ζ from a protein crystal structure. At low ionic strengths, a hyperdiffusivity regime exists, requiring direct measurement of electrophoretic mobility to determine ζ. This work, therefore, supports a basic tenant of EDL theory that the electric double layer during diffusion and electrophoresis are equivalent in the Stokes-Einstein regime.