Static light scattering from concentrated protein solutions II: experimental test of theory for protein mixtures and weakly self-associating proteins

Engineered polymeric excipients for enhancing the stability of protein biologics: Poly(N-isopropylacrylamide)-poly(ethylene glycol) (PNIPAM-PEG) block copolymers

Protein therapeutics, particularly antibodies, depend on maintaining their native structures for optimal function. Hydrophobic interfaces, such as the air-water interface, can trigger protein aggregation and denaturation. While completely avoiding such interfacial exposures during manufacturing and storage is impractical, minimizing them is crucial for enhancing protein drug stability and extending shelf life. In the biologics industry, surfactants like polysorbates are commonly used as additives (excipients) to mitigate these undesirable interfacial exposures. However, polysorbates, the most prevalent choice, have recognized limitations in terms of polydispersity, purity, and stability, prompting the exploration of alternative excipients. The present study identifies poly(N-isopropylacrylamide)-poly(ethylene glycol) (PNIPAM-PEG) block copolymers as a promising alternative to polysorbates. Due to its stronger affinity for the air-water interface, PNIPAM-PEG significantly outperforms polysorbates in enhancing protein stability. This claim is supported by results from multiple tests. Accelerated dynamic light scattering (DLS) experiments demonstrate PNIPAM-PEG's exceptional efficacy in preserving IgG stability against surface-induced aggregation, surpassing conventional polysorbate excipients (Tween 80 and Tween 20) under high-temperature conditions. Additionally, circular dichroism (CD) spectroscopy results reveal conformational alterations associated with aggregation, with PNIPAM-PEG consistently demonstrates a greater protective effect by mitigating negative shifts at λ ≅ 220 nm, indicative of changes in secondary structure. Overall, this study positions PNIPAM-PEG as a promising excipient for antibody therapeutics, facilitating the development of more stable and effective biopharmaceuticals.

In Silico View of Crowding: Biomolecular Processes to Nanomaterial Design

It is widely accepted that deciphering biomolecular structure and function requires going beyond the single-molecule or single-complex paradigm. The densely packed macromolecules, cosolutes, and metabolites in the living cell impose crowding effects on the biomolecular structure and dynamics that need to be accounted for. Molecular simulations have proven to be a powerful tool to advance the current molecular-level understanding of such a highly concentrated, complex milieu. This Mini-Review focuses on summarizing the understanding achieved so far for the effects of crowding on biomolecular processes using computational methods, along with highlighting a new direction in employing crowding as a tool for tunable nanomaterial design. The two schools of thought that form the pillars of the current understanding of crowding effects are discussed. The investigation of crowded solutions using physics-based models that encompass different time and length scales to mimic the intracellular environment are described. The limitations and challenges faced by the current models and simulation methods are addressed, highlighting the gaps to be filled for better agreement with experiments. Crowding can also act as an effective tool to modulate the structure-property-function relationships of nanomaterials, leading to the development of novel functional materials. A few recent studies, mostly experimental, have been summarized in this direction. The Mini-Review concludes with an outlook for future developments in this field in order to enable accurate mimicking of the intracellular environment using simulations and to bridge the gap between biological processes and nanomaterial design.

Simplified Equilibrium Model for Exploring the Combined Influences of Concentration, Aggregate Shape, Excluded Volume, and Surface Adsorption upon Aggregation Propensity and Distribution of Globular Macromolecules

A mesoscopic model for the equilibrium self-association of a globular macromolecule that may form oligomers of various shapes and unlimited sizes is presented. Allowance is made within this model for the effects of variation in the free energy of subunit contact within an oligomer of specified size and different shapes, the free energy of adsorption of an oligomer of specified size and shape to a planar surface, and the free energy of nonspecific excluded volume interaction between an oligomer of specified size and shape and an inert species occupying a specified fraction of total volume. The model is analytically soluble and permits rapid calculation and analysis of the effects of variation in each of the three free energy parameters upon the concentration dependence of the weight-average stoichiometry of the oligomer, the fraction of total macromolecule that is adsorbed, and the fraction of differently shaped oligomers that are adsorbed and in free solution.

Predictive Nature of High-Throughput Assays in ADC Formulation Screening

Utilization of high-throughput biophysical screening techniques during early screening studies is warranted due to the limited amount of material and large number of samples. But the predictability of the data to longer-term storage stability is critical as the high-throughput methods assist in defining the design space for the longer-term studies. In this study, the biophysical properties of two ADCs in 16 formulation conditions were evaluated using high-throughput techniques. Conformational stability and colloidal stability were evaluated by determining Tm values, kD, B22, and Tagg. In addition, the samples were placed on stability and the extent of aggregate formation over the 8-week interval was determined. The rank order of the 16 different formulations in the high-throughput assays was compared to the rank order observed during the stability studies to assess the predictive capabilities of the screening methods. It was demonstrated that similar rank orders can be expected between high-throughput physical stability indicating assays such as Tagg and B22 and traditional aggregation by SEC data, whereas conformational stability read-outs (Tm) are less predictive. In addition, the high-throughput assays appropriately identified the poor performing formulation conditions, which is ultimately what is desired of screening assays.

Effects of Cascading Optical Processes: Part I: Impacts on Quantification of Sample Scattering Extinction, Intensity, and Depolarization

Light scattering is a universal matter property that is especially prominent in nanoscale or larger materials. However, the effects of scattering-based cascading optical processes on experimental quantification of sample absorption, scattering, and emission intensities, as well as scattering and emission depolarization, have not been adequately addressed. Using a series of polystyrene nanoparticles (PSNPs) of different sizes as model analytes, we present a computational and experimental study on the effects of cascading light scattering on experimental quantification of NP scattering activities (scattering cross-section or molar coefficient), intensity, and depolarization. Part II and Part III of this series of companion articles explore the effects of cascading optical processes on sample absorption and fluorescence measurements, respectively. A general theoretical model is developed on how forward scattered light complicates the general applicability of Beer's law to the experimental UV-vis spectrum of scattering samples. The correlation between the scattering intensity and PSNP concentration is highly complicated with no robust linearity even when the scatterers' concentration is very low. Such complexity arises from the combination of concentration-dependence of light scattering depolarization and the scattering inner filter effects (IFEs). Scattering depolarization increases with the PSNP scattering extinction (thereby, its concentration) but can never reach unity (isotropic) due to the polarization dependence of the scattering IFE. The insights from this study are important for understanding the strengths and limitations of various scattering-based techniques for material characterization including nanoparticle quantification. They are also foundational for quantitative mechanistic understanding on the effects of light scattering on sample absorption and fluorescence measurements.

Near-Infrared Dye-Aptamer Assay for Small Molecule Detection in Complex Specimens

Aptamers are single-stranded oligonucleotides isolated in vitro that bind specific targets with high affinity and are commonly used as receptors in biosensors. Aptamer-based dye-displacement assays are a promising sensing platform because they are label-free, sensitive, simple, and rapid. However, these assays can exhibit impaired sensitivity in biospecimens, which contain numerous interferents that cause unwanted absorbance, scattering, and fluorescence in the UV-vis region. Here, this problem is overcome by utilizing near-infrared (NIR) signatures of the dye 3,3'-diethylthiadicarbocyanine iodide (Cy5). Cy5 initially complexes with aptamers as monomers and dimers; aptamer-target binding displaces the dye into solution, resulting in the formation of J-aggregates that provide a detectable NIR signal. The generality of our assay is demonstrated by detecting three different small-molecule analytes with their respective DNA aptamers at clinically relevant concentrations in serum and urine. These successful demonstrations show the utility of dye-aptamer NIR biosensors for high-throughput detection of analytes in clinical specimens.

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

Influence of Nonspecific Interactions on Protein Associations: Implications for Biochemistry In Vivo

The cellular interior is composed of a variety of microenvironments defined by distinct local compositions and composition-dependent intermolecular interactions. We review the various types of nonspecific interactions between proteins and between proteins and other macromolecules and supramolecular structures that influence the state of association and functional properties of a given protein existing within a particular microenvironment at a particular point in time. The present state of knowledge is summarized, and suggestions for fruitful directions of research are offered. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Bayesian analysis of static light scattering data for globular proteins

Static light scattering is a popular physical chemistry technique that enables calculation of physical attributes such as the radius of gyration and the second virial coefficient for a macromolecule (e.g., a polymer or a protein) in solution. The second virial coefficient is a physical quantity that characterizes the magnitude and sign of pairwise interactions between particles, and hence is related to aggregation propensity, a property of considerable scientific and practical interest. Estimating the second virial coefficient from experimental data is challenging due both to the degree of precision required and the complexity of the error structure involved. In contrast to conventional approaches based on heuristic ordinary least squares estimates, Bayesian inference for the second virial coefficient allows explicit modeling of error processes, incorporation of prior information, and the ability to directly test competing physical models. Here, we introduce a fully Bayesian model for static light scattering experiments on small-particle systems, with joint inference for concentration, index of refraction, oligomer size, and the second virial coefficient. We apply our proposed model to study the aggregation behavior of hen egg-white lysozyme and human γS-crystallin using in-house experimental data. Based on these observations, we also perform a simulation study on the primary drivers of uncertainty in this family of experiments, showing in particular the potential for improved monitoring and control of concentration to aid inference.

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.

The Cumulative Effect of Surface Adsorption and Excluded Volume in 2D and 3D on Protein Fibrillation

A refined mesoscopic model for the cumulative effect of repulsive excluded volume protein-protein interaction and attractive protein-surface interaction upon the properties of a trace protein capable of fiber formation is presented. The model predicts that very small changes in the magnitude of bulk volume occupancy or the strength of protein-surface attraction may result in very large changes in the extent of trace protein fibrillation and the distribution of trace protein between bulk and adsorbed phases.

Dynamic Light Scattering (DLS)

Focus a laser on dissolved particles and analyze the scattered light to reveal their size. This well established principle is used in dynamic light scattering (DLS), or also called photon-correlation spectroscopy, which is a widely popular and highly adaptable analytical method applied in different fields of life and material sciences, as well as in industrial quality control processes.

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

Non-specific Interactions Between Macromolecular Solutes in Concentrated Solution: Physico-Chemical Manifestations and Biochemical Consequences

A general thermodynamic formulation of the effect of hard and soft non-specific intermolecular interactions upon reaction equilibria is summarized. A highly simplified quantitative model for non-specific intermolecular interaction is introduced. This model is used to illustrate how the magnitudes of attractive and repulsive components of the overall intermolecular interaction, and the balance between them, influence the concentration-dependent properties of a highly concentrated solution of a single macromolecular solute. The properties calculated using the results of computer simulation and an approximate analytical model are found to agree qualitatively with the results of experimental measurements on protein solutions over a broad range of concentration.

Toward an understanding of biochemical equilibria within living cells

Four types of environmental effects that can affect macromolecular reactions in a living cell are defined: nonspecific intermolecular interactions, side reactions, partitioning between microenvironments, and surface interactions. Methods for investigating these interactions and their influence on target reactions in vitro are reviewed. Methods employed to characterize conformational and association equilibria in vivo are reviewed and difficulties in their interpretation cataloged. It is concluded that, in order to be amenable to unambiguous interpretation, in vivo studies must be complemented by in vitro studies carried out in well-characterized and controllable media designed to contain key elements of selected intracellular microenvironments.

Explicit Incorporation of Hard and Soft Protein-Protein Interactions into Models for Crowding Effects in Protein Mixtures. 2. Effects of Varying Hard and Soft Interactions upon Prototypical Chemical Equilibria

Previously derived approximate analytical relations for the activity coefficient of each solute in a mixture of up to three spherical solutes in a highly nonideal solution interacting via square well potentials of mean force (Hoppe, T.; Minton, A. P. J Phys Chem B. 2016, 120, 11866-11872) were used to explore the effect of heterogeneity in volume occupancy and intermolecular interactions upon prototypical schemes representing solubility, partitioning, conformational isomerization, and self-association in crowded solutions. Results generally indicate that all of the equilibria explored are exquisitely sensitive to variations in both volume occupancy and intermolecular interaction and have important implications for the design and execution of more detailed simulations of complex media.

Relating Protein-Protein Interactions and Aggregation Rates From Low to High Concentrations

At low protein concentrations (c2), non-native protein aggregation rates are known to be sensitive to changes in conformational stability and "weak" or "colloidal" protein-protein interactions. Protein-protein interactions are also known to be strong functions of c2. In the present work, protein-protein interactions and rates of aggregation were quantified systematically for a monoclonal antibody (MAb) across a broad range of c2 at pH 5.1 and 6.5, with or without 5 wt/wt % sucrose or 100 mM NaCl present. Aggregation rates were determined from initial-rate analysis with size-exclusion chromatography, and interactions were quantified with static and dynamic laser light scattering. A number of hypotheses were tested regarding whether changes in protein-protein interactions can be predictive of changes in aggregation rates versus c2. Hypotheses were based on (i) changes in thermodynamic activity; (ii) statistical mechanical fluctuation theory; and (iii) surface-contact probabilities. Arguments based on (i) and (ii) were qualitatively inconsistent with experimental rates and scattering. Hypothesis (iii) was reasonably successful and resulted in a semiquantitative correlation between rates and protein-protein interactions across almost 2 orders of magnitude in c2. However, (iii) requires one to assume that the concentration-dependent protein-protein Kirkwood-Buff integral is a reasonable surrogate for contact probabilities.

Macromolecular Crowding In Vitro, In Vivo, and In Between

Biochemical processes take place in heterogeneous and highly volume-occupied or crowded environments that can considerably influence the reactivity and distribution of participating macromolecules. We summarize here the thermodynamic consequences of excluded-volume and long-range nonspecific intermolecular interactions for macromolecular reactions in volume-occupied media. In addition, we summarize and compare the information content of studies of crowding in vitro and in vivo. We emphasize the importance of characterizing the behavior not only of labeled tracer macromolecules but also the composition and behavior of unlabeled macromolecules in the immediate vicinity of the tracer. Finally, we propose strategies for extending quantitative analyses of crowding in simple model systems to increasingly complex media up to and including intact cells.

Characterizing Reversible Protein Association at Moderately High Concentration Via Composition-Gradient Static Light Scattering

Analysis of weakly self-associating macromolecules at concentrations beyond a few g/L is challenging on account of the confounding effect of thermodynamic nonideality on the association signal. When the reversible association comprises only 1 or 2 oligomeric species in equilibrium with the monomer, the nonideality may be accounted for in a relatively rigorous manner, but if more association states are involved, the analysis becomes quite complex. We show that under reasonable assumptions, the nonideality in a composition-gradient static light scattering measurement may be accounted for in a simple fashion. The correction is applied to determining the stoichiometry and binding affinity of a protein previously characterized via sedimentation equilibrium and dynamic light scattering. The results of the new analysis are remarkably self-consistent and in line with the expectations for the form of self-association predicted previously from analysis of the surface residuals, establishing composition-gradient multi-angle static light scattering with nonideality corrections as a critical technology for characterizing associative interactions in concentrated solutions.

Charged Molecules Modulate the Volume Exclusion Effects Exerted by Crowders on FtsZ Polymerization

We have studied the influence of protein crowders, either combined or individually, on the GTP-induced FtsZ cooperative assembly, crucial for the formation of the dynamic septal ring and, hence, for bacterial division. It was earlier demonstrated that high concentrations of inert polymers like Ficoll 70, used to mimic the crowded cellular interior, favor the assembly of FtsZ into bundles with slow depolymerization. We have found, by fluorescence anisotropy together with light scattering measurements, that the presence of protein crowders increases the tendency of FtsZ to polymerize at micromolar magnesium concentration, being the effect larger with ovomucoid, a negatively charged protein. Neutral polymers and a positively charged protein also diminished the critical concentration of assembly, the extent of the effect being compatible with that expected according to pure volume exclusion models. FtsZ polymerization was also observed to be strongly promoted by a negatively charged polymer, DNA, and by some unrelated polymers like PEGs at concentrations below the crowding regime. The influence of mixed crowders mimicking the heterogeneity of the intracellular environment on the tendency of FtsZ to assemble was also studied and nonadditive effects were found to prevail. Far from exactly reproducing the bacterial cytoplasm environment, this approach serves as a simplified model illustrating how its intrinsically crowded and heterogeneous nature may modulate FtsZ assembly into a functional Z-ring.

Quantitative characterization of nonspecific self- and hetero-interactions of proteins in nonideal solutions via static light scattering

The dependence of static light scattering upon the compositions of solutions including hen egg white ovalbumin, hen egg white ovomucoid, ribonuclease A, and binary mixtures of these proteins at total concentrations of up to about 40 g/L were measured at different values of the pH and ionic strength. At the pH values of measurement, ovalbumin and ovomucoid have a net negative charge and ribonuclease A has a net positive charge. The observed dependence of scattering intensity upon solution composition may be accounted for by an extension of previously formulated equivalent hard particle models that allows for the presence of both repulsive interactions between like species and attractive interactions between unlike species in mixtures of positively and negatively charged proteins.

Combined dynamic light scattering and Raman spectroscopy approach for characterizing the aggregation of therapeutic proteins

Determination of the physicochemical properties of protein therapeutics and their aggregates is critical for developing formulations that enhance product efficacy, stability, safety and manufacturability. Analytical challenges are compounded for materials: (1) that are formulated at high concentration, (2) that are formulated with a variety of excipients, and (3) that are available only in small volumes. In this article, a new instrument is described that measures protein secondary and tertiary structure, as well as molecular size, over a range of concentrations and formulation conditions of low volume samples. Specifically, characterization of colloidal and conformational stability is obtained through a combination of two well-established analytical techniques: dynamic light scattering (DLS) and Raman spectroscopy, respectively. As the data for these two analytical modalities are collected on the same sample at the same time, the technique enables direct correlation between them, in addition to the more straightforward benefit of minimizing sample usage by providing multiple analytical measurements on the same aliquot non-destructively. The ability to differentiate between unfolding and aggregation that the combination of these techniques provides enables insights into underlying protein aggregation mechanisms. The article will report on mechanistic insights for aggregation that have been obtained from the application of this technique to the characterization of lysozyme, which was evaluated as a function of concentration and pH.

Determination of the second virial coefficient of bovine serum albumin under varying pH and ionic strength by composition-gradient multi-angle static light scattering

Composition-gradient multi-angle static light scattering (CG-MALS) is an emerging technique for the determination of intermolecular interactions via the second virial coefficient B22. With CG-MALS, detailed studies of the second virial coefficient can be carried out more accurately and effectively than with traditional methods. In addition, automated mixing, delivery and measurement enable high speed, continuous, fluctuation-free sample delivery and accurate results. Using CG-MALS we measure the second virial coefficient of bovine serum albumin (BSA) in aqueous solutions at various values of pH and ionic strength of a univalent salt (NaCl). The systematic variation of the second virial coefficient as a function of pH and NaCl strength reveals the net charge change and the isoelectric point of BSA under different solution conditions. The magnitude of the second virial coefficient decreases to 1.13 x 10(-5) ml*mol/g(2) near the isoelectric point of pH 4.6 and 25 mM NaCl. These results illuminate the role of fundamental long-range electrostatic and van der Waals forces in protein-protein interactions, specifically their dependence on pH and ionic strength.

Optimized UV detection of high-concentration antibody formulations using high-throughput SE-HPLC

High-concentration antibody solutions (>100 mg/mL) present significant challenges for formulation and process development, including formulation attributes such as increased solution viscosity, and the propensity for self-association. An additional challenge comes from the adaptation of analytical methods designed for low-concentration formulations to the high-concentration regime. The oligomeric state is a good example: it is a quality attribute monitored during pharmaceutical development and is one that can be affected by dilution; a typical first step in the analysis of high-concentration solutions. The objective of this work was to develop a size-exclusion HPLC (SE-HPLC) method that would allow the injection of high-concentration antibody formulations without the need for dilution prior to injection and their analysis in a high-throughput manner that does not create a bottleneck for the execution of complex formulation development studies. It was found that changing the UV detection wavelength from 215 to 235 nm simplified sample preparation by allowing for an approximately fivefold increase in injection load while maintaining the signal within the linear range of detection. In addition, the chromatographic peak properties (i.e., peak symmetry, resolution, and sensitivity) were determined to be consistent when compared with analytical methods developed for formulations with lower antibody concentrations.

Role of anisotropic interactions for proteins and patchy nanoparticles

Protein–protein interactions are inherently anisotropic to some degree, with orientation-dependent interactions between repulsive and attractive or complementary regions or “patches” on adjacent proteins. In some cases it has been suggested that such patch–patch interactions dominate the thermodynamics of dilute protein solutions, as captured by the osmotic second virial coefficient (B₂₂), but delineating when this will or will not be the case remains an open question. A series of simplified but exactly solvable models are first used to illustrate that a delicate balance exists between the strength of attractive patch–patch interactions and the patch size, and that repulsive patch–patch interactions contribute significantly to B₂₂ for only those conditions where the repulsions are long-ranged. Finally, B₂₂ is reformulated, without approximations, in terms of the density of states for a given interaction energy and particle–particle distance. Doing so illustrates the inherent balance of entropic and energetic contributions to B₂₂. It highlights that simply having strong patch–patch interactions will only cause anisotropic interactions to dominate B₂₂ solution properties if the unavoidable entropic penalties are overcome, which cannot occur if patches are too small. The results also indicate that the temperature dependence of B₂₂ may be a simple experimental means to assess whether a small number of strongly attractive configurations dominate the dilute solution behavior.

Protein aggregation and its impact on product quality

Protein pharmaceutical products are typically active as folded monomers that are composed of one or more protein chains, such as the heavy and light chains in monoclonal antibodies that are a mainstay of current drug pipelines. There are numerous possible aggregated states for a given protein, some of which are potentially useful, while most of which are considered deleterious from the perspective of pharmaceutical product quality and performance. This review provides an overview of how and why different aggregated states of proteins occur, how this potentially impacts product quality and performance, fundamental approaches to control aggregate formation, and the practical approaches that are currently used in the pharmaceutical industry.

Light-scattering-based analysis of biomolecular interactions

While light scattering has long been applied to the analysis of biomolecular interactions, recent advances have extended the practical use of light scattering techniques to cover a rather broad range of phenomena. In this paper I review essential light scattering theory as applied to specific interactions under thermodynamically ideal conditions and present examples showing how light scattering elucidates the dynamic equilibrium and kinetic behavior of proteins and other biomacromolecules.

Quantitative assessment of the relative contributions of steric repulsion and chemical interactions to macromolecular crowding

The term "macromolecular crowding" denotes the combined effects of high volume fractions of nominally unrelated macromolecules upon the equilibrium and transport properties of all macrosolutes, dilute as well as concentrated, in the crowded medium. We present a formal partitioning of the total crowding effect into contributions from steric exclusion (excluded volume) and weak, nonspecific attractive interactions between a concentrated "crowding agent" and reactant and product species present at trace concentration. A numerical example of the combined effect of both steric and chemical interactions between crowder and tracer upon the reversible dimerization of tracer is presented, based upon reasonable estimates of the magnitude of both repulsive and attractive interactions between tracer and crowder species.

Quantitative characterization of the interaction between sucrose and native proteins via static light scattering

The composition-dependent static light scattering of binary mixtures of each of four dilute globular proteins and sucrose were measured over a broad range of sucrose concentrations. A conventional analysis of the dependence of excess scattering of a single macrosolute in a continuum solvent yields unphysical results. The data are reanalyzed in the context of multicomponent light scattering theory to yield the dependence of the free energy of solvation of each protein upon the concentration of sucrose. The results could be satisfactorily accounted for by an effective hard particle model that indicates the nature of the underlying interactions between sucrose and each protein.

Reexamining protein-protein and protein-solvent interactions from Kirkwood-Buff analysis of light scattering in multi-component solutions

The classic analysis of Rayleigh light scattering (LS) is re-examined for multi-component protein solutions, within the context of Kirkwood-Buff (KB) theory as well as a more generalized canonical treatment. Significant differences arise when traditional treatments that approximate constant pressure and neglect concentration fluctuations in one or more (co)solvent/co-solute species are compared with more rigorous treatments at constant volume and with all species free to fluctuate. For dilute solutions, it is shown that LS can be used to rigorously and unambiguously obtain values for the osmotic second virial coefficient (B(22)), in contrast with recent arguments regarding protein interactions deduced from LS experiments. For more concentrated solutions, it is shown that conventional analysis over(under)-estimates the magnitude of B(22) for significantly repulsive(attractive) conditions, and that protein-protein KB integrals (G(22)) are the more relevant quantity obtainable from LS. Published data for α-chymotrypsinogen A and a series of monoclonal antibodies at different pH and salt concentrations are re-analyzed using traditional and new treatments. The results illustrate that while traditional analysis may be sufficient if one is interested in only the sign of B(22) or G(22), the quantitative values can be significantly in error. A simple approach is illustrated for determining whether protein concentration (c(2)) is sufficiently dilute for B(22) to apply, and for correcting B(22) values from traditional LS regression at higher c(2) values. The apparent molecular weight M(2, app) obtained from LS is shown to generally not be equal to the true molecular weight, with the differences arising from a combination of protein-solute and protein-cosolute interactions that may, in principle, also be determined from LS.

Quantitative characterization of temperature-independent and temperature-dependent protein-protein interactions in highly nonideal solutions

The interaction among each of three dilute "tracer" proteins (bovine serum albumin, superoxide dismutase, and ovomucoid) at a concentration of 2 mg/mL and each of two "crowder" proteins (ovomucoid and BSA) at concentrations up to 100 mg/mL was characterized by analysis of dependence of the equilibrium gradients of both tracer and crowder upon the concentration of crowder. The equilibrium gradients of both crowder proteins were found to be independent of temperature over the range 5-37 °C. The equilibrium gradients of tracer BSA and ovomucoid in the complementary crowder species were likewise found to be independent of temperature over this range, indicating that interaction among these tracers and crowders is predominantly repulsive and essentially entirely entropic in nature. The equilibrium gradient of tracer SOD in BSA was also found to be independent of temperature over this range, but the gradient of tracer SOD in ovomucoid depended significantly upon temperature in a manner indicating a significant enthalpic (attractive) component of the overall interaction between SOD and ovomucoid. The experimental data are analyzed using model-free expansions of the thermodynamic activity coefficients of tracer and crowder in powers of the concentration of crowder and using approximate statistical thermodynamic models based upon highly simplified descriptions of molecular structure and interactions. Detailed analysis of the results indicates a relatively small contribution of nonspecific attraction to the total protein-protein interaction, which is dominated by steric repulsion.

On the distribution of protein refractive index increments

The protein refractive index increment, dn/dc, is an important parameter underlying the concentration determination and the biophysical characterization of proteins and protein complexes in many techniques. In this study, we examine the widely used assumption that most proteins have dn/dc values in a very narrow range, and reappraise the prediction of dn/dc of unmodified proteins based on their amino acid composition. Applying this approach in large scale to the entire set of known and predicted human proteins, we obtain, for the first time, to our knowledge, an estimate of the full distribution of protein dn/dc values. The distribution is close to Gaussian with a mean of 0.190 ml/g (for unmodified proteins at 589 nm) and a standard deviation of 0.003 ml/g. However, small proteins <10 kDa exhibit a larger spread, and almost 3000 proteins have values deviating by more than two standard deviations from the mean. Due to the widespread availability of protein sequences and the potential for outliers, the compositional prediction should be convenient and provide greater accuracy than an average consensus value for all proteins. We discuss how this approach should be particularly valuable for certain protein classes where a high dn/dc is coincidental to structural features, or may be functionally relevant such as in proteins of the eye.

Stability of seasonal influenza vaccines investigated by spectroscopy and microscopy methods

The stability of different seasonal influenza vaccines was investigated by spectroscopy and microscopy methods before and after the following stress-conditions: (i) 2 and 4 weeks storage at 25°C, (ii) 1 day storage at 37°C and (iii) one freeze-thaw cycle. The subunit vaccine Influvac (Solvay Pharma) and the split vaccine Mutagrip (Sanofi Pasteur) were affected by all stresses. The split vaccine Fluarix (GlaxoSmithKline) was affected only by storage at 25°C. The virosomal vaccine Inflexal V (Berna Biotech) was stable after the temperature stresses but aggregated after one freeze-thaw cycle. This study provides new insights into commercial vaccines of low antigen concentration and highlights the importance of using multiple techniques to assess vaccine stability.

Beyond the second virial coefficient: Sedimentation equilibrium in highly non-ideal solutions

The general theory of sedimentation equilibrium (SE), applicable to mixtures of interacting sedimentable solutes at arbitrary concentration, is summarized. Practical techniques for the acquisition of SE data suitable for analysis are described. Experimental measurements and analyses of SE in concentrated protein solutions are reviewed. The method of non-ideal tracer sedimentation equilibrium (NITSE) is described. Experimental studies using NITSE to detect and quantitatively characterize intermolecular interactions in mixtures of dilute tracer species and concentrated proteins or polymers are reviewed.

Predicting solution aggregation rates for therapeutic proteins: approaches and challenges

Non-native aggregation is a common concern during therapeutic protein product development and manufacturing, particularly for liquid dosage forms. Because aggregates are often net irreversible under the conditions that they form, controlling aggregate levels requires control of aggregation rates across a range of solution conditions. Rational design of product formulation(s) would therefore benefit greatly from methods to accurately predict aggregation rates. This article focuses on the principles underlying current rate-prediction approaches for non-native aggregation, the limitations and strengths of different approaches, and illustrative examples from the authors' laboratories. The analysis highlights a number of reasons why accurate prediction of aggregation rates remains an outstanding challenge, and suggests some of the important areas for research to ultimately enable improved predictive capabilities in the future.

Indefinite noncooperative self-association of chicken deoxy hemoglobin D

The minor tetrameric hemoglobin (Hb), Hb D, of chicken red blood cells self-associates upon deoxygenation. This self-association enhances the cooperativity of oxygen binding. The maximal Hill coefficient is greater than 4 at high Hb concentrations. Previous measurements at low Hb concentrations were consistent with a monomer-to-dimer equilibrium and an association constant of ∼1.3-1.6 × 10(4) M(-1). Here, the Hb tetramer is considered as the monomer. However, new results indicate that the association extends beyond the dimer. We show by combination of Hb oligomer modeling and sedimentation velocity analyses that the data can be well described by an indefinite noncooperative or isodesmic association model. In this model, the deoxy Hb D associates noncooperatively to give a linear oligomeric chain with an equilibrium association constant of 1.42 × 10(4) M(-1) at 20°C for each step. The data are also well described by a monomer-dimer-tetramer equilibrium model with monomer-to-dimer and dimer-to-tetramer association constants of 1.87 and 1.03 × 10(4) M(-1) at 20°C, respectively. A hybrid recombinant Hb D was prepared with recombinant α(D)-globin and native β-globin to give a Hb D tetramer (α(2)(D)β(2)). This rHb D undergoes decreased deoxygenation-dependent self-association compared with the native Hb D. Residue glutamate 138 has previously been proposed to influence intertetramer interactions. Our results with recombinant Hb D show that Glu138 plays no role in deoxy Hb D intertetramer interactions.

Intermolecular interactions of IgG1 monoclonal antibodies at high concentrations characterized by light scattering

Light scattering intensity measurements of solutions of two purified monoclonal antibodies were performed over a wide range of concentrations (0.5-275 mg/mL) and ionic strengths (0.02 to 0.6 M). Despite extensive sequence homology between these mAbs, alteration of ∼20 amino acids in the complementarity determining regions resulted in different net intermolecular interactions and responses to solution ionic strength. The concentration dependence of scattering was analyzed by comparison with the predictions of three models, allowing for intermolecular interaction of various types. In order of increasing complexity, the three models account for: (1) steric repulsions (simple hard-sphere model), (2) steric repulsion with short-ranged attractive interactions of varying magnitude (adhesive hard-sphere model), and (3) steric and nonsteric repulsive interactions between several species whose relative concentrations may change as a function of total protein concentration as dictated by equilibrium self-association (effective hard-sphere mixture model). Simple scattering models of noninteracting and adhesive hard-sphere species permitted qualitative interpretation of contributions from excluded volume, electrostatic, and van der Waals interactions on net mAb interactions at high concentration as a function of ionic strength. mAb2 electrostatic interactions were repulsive, whereas mAb1 interactions were net attractive at low ionic strengths, attributed to an anisotropic distribution of molecular charge. The effective hard-sphere mixture model can account quantitatively for the dependence of scattering for both antibodies over the entire concentration range and at salt concentrations exceeding 40 mM. This analysis showed that at high ionic strength both mAbs self-associate weakly to form dimer with an affinity that varies little with salt concentration at concentrations exceeding 75 mM. In addition, mAb1 appears to self-associate further to form oligomers with stoichiometry of 4-6 and an affinity that declines substantially with increasing ionic strength. All three models lead to the conclusion that at high concentrations repulsive interactions are predominantly due to excluded volume, whereas additional features are salt-dependent and reflect a substantial electrostatic contribution to intermolecular interactions of both mAbs.

Effect of nonadditive repulsive intermolecular interactions on the light scattering of concentrated protein-osmolyte mixtures

The static light scattering of three globular proteins, bovine serum albumin, ovalbumin, and ovomucoid, and binary mixtures of each protein and trimethylamine oxide (TMAO) containing between 10 and 70% protein, were measured as a function of total weight per volume concentration up to 100 g/L. The observed dependence of scattering upon concentration may be accounted for quantitatively by an effective hard sphere model incorporating an extension that takes into account the nonadditive nature of the repulsive intermolecular interaction between protein and TMAO.

Quantitative characterization of polymer-polymer, protein-protein, and polymer-protein interaction via tracer sedimentation equilibrium

Quantitative analysis of the composition dependence of the concentration gradient of each species of macromolecule within a solution mixture at sedimentation equilibrium permits the quantitative characterization of self- and heterointeractions between sedimenting solutes. Sedimentation equilibrium experiments were conducted on solutions containing a trace concentration of FITC-labeled BSA in varying concentrations of Ficoll 70 and on solutions containing a trace concentration of FITC-labeled Ficoll 70 in varying concentrations of BSA. The equilibrium gradient of each solute component in each mixture was measured independently. Analysis of the resulting gradients resulted in evaluation of the dependence of the activity coefficient of Ficoll upon the concentrations of Ficoll and BSA at concentrations of up to 100 g/L and the dependence of the activity coefficient of BSA upon the concentrations of Ficoll and BSA at concentrations of up to 100 g/L. The activity coefficients of both species increase significantly with increasing Ficoll and BSA concentration and do not vary with temperature, to within the precision of measurement, over the temperature range of 5-37 degrees C, indicating that the dominant interaction between Ficoll molecules and between BSA and Ficoll molecules is repulsive and probably due to steric volume exclusion. The measured dependences may be accounted for quantitatively by a simple model in which BSA and Ficoll 70 are represented by equivalent rigid particles.

Allowance for the effect of protein charge in the characterization of nonideal solute self-association by sedimentation equilibrium

This theoretical investigation explores the use of statistical-mechanical approaches to characterize the reversible tetramerization of a protein monomer with the size and charge characteristics of serum albumin under conditions where consideration of nearest-neighbor interactions suffices to describe effects of thermodynamic non-ideality. Such analysis of simulated sedimentation equilibrium distributions points to the adequacy of both the scaled particle theory and potential-of-mean-force methods for determining the self-association constant. Although the latter method usually entails the assignment of a magnitude to monomer net charge, this requirement can be obviated to some extent by repeating the analysis for a range of monomer charges and identifying the most appropriate value as that associated with a minimum in the sum-of-squares-of-residuals (SSR) of the best-fit descriptions of the sedimentation equilibrium distribution. Reasonable estimates of the association constant are usually obtained from corresponding analyses of the same sedimentation equilibrium distributions with activity coefficients obtained by scaled particle theory, an approach which also involves the identification of parameters on the basis of a minimum in SSR. However, the value of monomer charge determined must be regarded as a curve-fitting parameter rather than a true measure of monomer charge. Similar qualifications are shown to prevail in the scaled particle theory approach, which also involves the identification of parameters (the effective monomer volume and the polymer/monomer volume ratio) on the basis of a minimum in SSR. We therefore recommend discontinuation of the practice whereby quite precise distinction between modes of self-association has been attempted on the grounds of the physical credibility of the magnitudes of these additional curve-fitting parameters.

Models of macromolecular crowding effects and the need for quantitative comparisons with experiment

In recent years significant effort has been devoted to exploring the potential effects of macromolecular crowding on protein folding and association phenomena. Theoretical calculations and molecular simulations have, in particular, been exploited to describe aspects of protein behavior in crowded and confined conditions and many aspects of the simulated behavior have reflected, at least at a qualitative level, the behavior observed in experiments. One major and immediate challenge for the theorists is to now produce models capable of making quantitatively accurate predictions of in vitro behavior. A second challenge is to derive models that explain results obtained from experiments performed in vivo, the results of which appear to call into question the assumed dominance of excluded-volume effects in vivo.

Direct allowance for the effects of thermodynamic nonideality in the quantitative characterization of protein self-association by osmometry

A procedure is described for the direct analysis of osmotic pressure data for reversibly dimerizing proteins that makes allowance for effects of thermodynamic nonideality on the statistical-mechanical basis of the potential-of-mean-force between molecules. Detailed consideration is also given to calculation of the magnitudes of the required virial coefficients. After illustration of the approach with analysis of simulated osmotic pressure data, the method is used to obtain dimerization constants from published osmotic pressure data for soybean proteinase inhibitor, hemoglobin and alpha-chymotrypsin.