The limitations of an exclusively colloidal view of protein solution hydrodynamics and rheology

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

Characterizing Protein-Protein Interactions and Viscosity of a Monoclonal Antibody from Low to High Concentration Using Small-Angle X-ray Scattering and Molecular Dynamics Simulations

Understanding protein-protein interactions and formation of reversible oligomers (clusters) in concentrated monoclonal antibody (mAb) solutions is necessary for designing stable, low viscosity (η) concentrated formulations for processing and subcutaneous injection. Here we characterize the strength (K) of short-range anisotropic attractions (SRA) for 75-200 mg/mL mAb2 solutions at different pH and cosolute conditions by analyzing structure factors (Seff(q)) from small-angle X-ray scattering (SAXS) using coarse-grained molecular dynamics simulations. Best fit simulations additionally provide cluster size distributions, fractal dimensions, cluster occluded volume, and mAb coordination numbers. These equilibrium properties are utilized in a model to account for increases in viscosity caused by occluded volume in the clusters (packing effects) and dissipation of stress across lubricated fractal clusters. Seff(q) is highly sensitive to K at 75 mg/mL where mAbs can mutually align to form SRA contacts but becomes less sensitive at 200 mg/mL as steric repulsion due to packing becomes dominant. In contrast, η at 200 mg/mL is highly sensitive to SRA and the average cluster size from SAXS/simulation, which is observed to track the cluster relaxation time from shear thinning. By analyzing the distribution of sub-bead hot spots on the 3D mAb surface, we identify a strongly attractive hydrophobic patch in the complementarity determining region (CDR) at pH 4.5 that contributes to the high K and consequently large cluster sizes and high η. Adding NaCl screens electrostatic interactions and increases the impact of hydrophobic attraction on cluster size and raises η, whereas nonspecific binding of Arg attenuates all SRA, reducing η. The hydrophobic patch is absent at higher pH values, leading to smaller K, smaller clusters, and lower η. This work constitutes a first attempt to use SAXS and CG modeling to link both structural and rheological properties of concentrated mAb solutions to the energetics of specific hydrophobic patches on mAb surfaces. As such, our work opens an avenue for future research, including the possibility of designing coarse-grained models with physically meaningful interacting hot spots.

Hydrocolloids of Egg White and Gelatin as a Platform for Hydrogel-Based Tissue Engineering

Innovative materials are needed to produce scaffolds for various tissue engineering and regenerative medicine (TERM) applications, including tissue models. Materials derived from natural sources that offer low production costs, easy availability, and high bioactivity are highly preferred. Chicken egg white (EW) is an overlooked protein-based material. Whilst its combination with the biopolymer gelatin has been investigated in the food technology industry, mixed hydrocolloids of EW and gelatin have not been reported in TERM. This paper investigates these hydrocolloids as a suitable platform for hydrogel-based tissue engineering, including 2D coating films, miniaturized 3D hydrogels in microfluidic devices, and 3D hydrogel scaffolds. Rheological assessment of the hydrocolloid solutions suggested that temperature and EW concentration can be used to fine-tune the viscosity of the ensuing gels. Fabricated thin 2D hydrocolloid films presented globular nano-topography and in vitro cell work showed that the mixed hydrocolloids had increased cell growth compared with EW films. Results showed that hydrocolloids of EW and gelatin can be used for creating a 3D hydrogel environment for cell studies inside microfluidic devices. Finally, 3D hydrogel scaffolds were fabricated by sequential temperature-dependent gelation followed by chemical cross-linking of the polymeric network of the hydrogel for added mechanical strength and stability. These 3D hydrogel scaffolds displayed pores, lamellae, globular nano-topography, tunable mechanical properties, high affinity for water, and cell proliferation and penetration properties. In conclusion, the large range of properties and characteristics of these materials provide a strong potential for a large variety of TERM applications, including cancer models, organoid growth, compatibility with bioprinting, or implantable devices.

Microfluidic techniques for mechanical measurements of biological samples

The use of microfluidics to make mechanical property measurements is increasingly common. Fabrication of microfluidic devices has enabled various types of flow control and sensor integration at micrometer length scales to interrogate biological materials. For rheological measurements of biofluids, the small length scales are well suited to reach high rates, and measurements can be made on droplet-sized samples. The control of flow fields, constrictions, and external fields can be used in microfluidics to make mechanical measurements of individual bioparticle properties, often at high sampling rates for high-throughput measurements. Microfluidics also enables the measurement of bio-surfaces, such as the elasticity and permeability properties of layers of cells cultured in microfluidic devices. Recent progress on these topics is reviewed, and future directions are discussed.

Simultaneous Sizing and Refractive Index Analysis of Heterogeneous Nanoparticle Suspensions

Rapid and reliable characterization of heterogeneous nanoparticle suspensions is a key technology across the nanosciences. Although approaches exist for homogeneous samples, they are often unsuitable for polydisperse suspensions, as particles of different sizes and compositions can lead to indistinguishable signals at the detector. Here, we introduce holographic nanoparticle tracking analysis, holoNTA, as a straightforward methodology that decouples size and material refractive index contributions. HoloNTA is applicable to any heterogeneous nanoparticle sample and has the sensitivity to measure the intrinsic heterogeneity of the sample. Specifically, we combined high dynamic range k-space imaging with holographic 3D single-particle tracking. This strategy enables long-term tracking by extending the imaging volume and delivers precise and accurate estimates of both scattering amplitude and diffusion coefficient of individual nanoparticles, from which particle refractive index and hydrodynamic size are determined. We specifically demonstrate, by simulations and experiments, that irrespective of localization uncertainty and size, the sizing sensitivity is improved as our extended detection volume yields considerably longer particle trajectories than previously reported by comparable technologies. As validation, we measured both homogeneous and heterogeneous suspensions of nanoparticles in the 40-250 nm size range and further monitored protein corona formation, where we identified subtle differences between the nanoparticle-protein complexes derived from avidin, bovine serum albumin, and streptavidin. We foresee that our approach will find many applications of both fundamental and applied nature where routine quantification and sizing of nanoparticles are required.

Molecular rotors in haemoglobin and bovine serum albumin proteins

Molecular rotors are fluorescent viscosity probes and their response in simple fluids is known to be a Förster-Hoffman power law, allowing the viscosity of the medium to be quantified by its fluorescence intensity. They are attractive probes in biological media, usually consisting of proteins, but how does a molecular rotor behave in a protein solution? The response of the DASPI molecular rotor is compared in two globular protein solutions of similar size, haemoglobin (Hb) and bovine serum albumin, one absorbent, the other not. In absorbent Hb, a model validated by experiments in triangular geometry allows one to correct the absorbing effect and to compare the rotor response in both proteins. With concomitant microrheology measurements, we investigate the relation between the DASPI fluorescence intensity and solution viscosity. In protein solutions, we show that viscosity is no longer the parameter determining the rotor response in contrast to simple fluids. Varying the viscosity by concentration or temperature is not equivalent, and the Förster-Hoffmann power laws do not apply when the solution concentration varies. We show that the concentration regime of the protein solution, semi-dilute or concentrated, determines the sensitivity of the rotor to its environment.

Application of a Simple Short-Range Attraction and Long-Range Repulsion Colloidal Model toward Predicting the Viscosity of Protein Solutions

Some hard-sphere colloidal models have been criticized for inaccurately predicting the solution viscosity of complex biological molecules like proteins. Competing short-range attractions and long-range repulsions, also known as short-range attraction and long-range repulsion (SALR) interactions, have been thought to affect the microstructure of a protein solution at low to moderate ionic strength. However, such interactions have been implicated primarily in causing phase transition, protein gelation, or reversible cluster formation, and their effect on protein solution viscosity change is not fully understood. In this work, we show the application of a hard-sphere colloidal model with SALR interactions toward predicting the viscosity of dilute to semi-dilute protein solutions. The comparison is performed for a globular-shaped albumin and Y-shaped therapeutic monoclonal antibody that are not explained by previous colloidal models. The model predictions show that it is the coupling between attractions and repulsions that gives rise to the observed experimental trends in solution viscosity as a function of pH, concentration, and ionic strength. The parameters of the model are obtained from measurements of the second virial coefficient and net surface charge/zeta-potential, without additional fitting of the viscosity.

Colloid-like solution behavior of computationally designed coiled coil bundlemers

The use of isotropic potential models of simple colloids for describing complex protein-protein interactions is a topic of ongoing debate in the biophysical community. This contention stems from the unavailability of synthetic protein-like model particles that are amenable to systematic experimental characterization. In this article, we test the utility of colloidal theory to capture the solution structure, interactions and dynamics of novel globular protein-mimicking, computationally designed peptide assemblies called bundlemers that are programmable model systems at the intersection of colloids and proteins. Small-angle neutron scattering (SANS) measurements of semi-dilute bundlemer solutions in low and high ionic strength solution indicate that bundlemers interact locally via repulsive interactions that can be described by a screened repulsive potential. We also present neutron spin echo (NSE) spectroscopy results that show high-Q freely-diffusive dynamics of bundlemers. Importantly, formation of clusters due to short-range attractive, inter-bundlemer interactions is observed in SANS even at dilute bundlemer concentrations, which is indicative of the complexity of the bundlemer charged surface. The similarities and differences between bundlemers and simple colloidal as well as complex protein-protein interactions is discussed in detail.

Deconstruction and Reassembly of Renewable Polymers and Biocolloids into Next Generation Structured Materials

This review considers the most recent developments in supramolecular and supraparticle structures obtained from natural, renewable biopolymers as well as their disassembly and reassembly into engineered materials. We introduce the main interactions that control bottom-up synthesis and top-down design at different length scales, highlighting the promise of natural biopolymers and associated building blocks. The latter have become main actors in the recent surge of the scientific and patent literature related to the subject. Such developments make prominent use of multicomponent and hierarchical polymeric assemblies and structures that contain polysaccharides (cellulose, chitin, and others), polyphenols (lignins, tannins), and proteins (soy, whey, silk, and other proteins). We offer a comprehensive discussion about the interactions that exist in their native architectures (including multicomponent and composite forms), the chemical modification of polysaccharides and their deconstruction into high axial aspect nanofibers and nanorods. We reflect on the availability and suitability of the latter types of building blocks to enable superstructures and colloidal associations. As far as processing, we describe the most relevant transitions, from the solution to the gel state and the routes that can be used to arrive to consolidated materials with prescribed properties. We highlight the implementation of supramolecular and superstructures in different technological fields that exploit the synergies exhibited by renewable polymers and biocolloids integrated in structured materials.

A microfluidic approach to studying the injection flow of concentrated albumin solutions

Abstract

Subcutaneous injection by means of prefilled syringes allows patients to self-administrate high-concentration (100 g/L or more) protein-based drugs. Although the shear flow of concentrated globulins or monoclonal antibodies has been intensively studied and related to the injection force proper of SC processes, very small attention has been paid to the extensional behavior of this category of complex fluids. This work focuses on the flow of concentrated bovine serum albumin (BSA) solutions through a microfluidic “syringe-on-chip” contraction device which shares some similarities with the geometry of syringes used in SC self-injection. By comparing the velocity and pressure measurements in complex flow with rheometric shear measurements obtained by means of the “Rheo-chip” device, it is shown that the extensional viscosity plays an important role in the injection process of protinaceous drugs.

Article Highlights

A microfluidic “syringe on chip” device mimicking the injection flow of protinaceous drugs has been developed.

The velocity field of concentrated BSA solutions through the “syringe on chip” is Newtonian-like.

The extensional viscosity of concentrated protein solutions should also be considered when computing injection forces through needles.

Droplet-Based Microfluidic Tool to Quantify Viscosity of Concentrating Protein Solutions

PURPOSE

Measurement of the viscosity of concentrated protein solutions is vital for the manufacture and delivery of protein therapeutics. Conventional methods for viscosity measurements require large solution volumes, creating a severe limitation during the early stage of protein development. The goal of this work is to develop a robust technique that requires minimal sample.

METHODS

In this work, a droplet-based microfluidic device is developed to quantify the viscosity of protein solutions while concentrating in micrometer-scale droplets. The technique requires only microliters of sample. The corresponding viscosity is characterized by multiple particle tracking microrheology (MPT).

RESULTS

We show that the viscosities quantified in the microfluidic device are consistent with macroscopic results measured by a conventional rheometer for poly(ethylene) glycol (PEG) solutions. The technique was further applied to quantify viscosities of well-studied lysozyme and bovine serum albumin (BSA) solutions. Comparison to both macroscopic measurements and models (Krieger-Dougherty model) demonstrate the validity of the approach.

CONCLUSION

The droplet-based microfluidic device provides accurate quantitative values of viscosity over a range of concentrations for protein solutions with small sample volumes (~ μL) and high compositional resolution. This device will be extended to study the effect of different excipients and other additives on the viscosity of protein solutions.

Towards an improved prediction of concentrated antibody solution viscosity using the Huggins coefficient

The viscosity of a monoclonal antibody solution must be monitored and controlled as it can adversely affect product processing, packaging and administration. Engineering low viscosity mAb formulations is challenging as prohibitive amounts of material are required for concentrated solution analysis, and it is difficult to predict viscosity from parameters obtained through low-volume, high-throughput measurements such as the interaction parameter, k_D, and the second osmotic virial coefficient, B₂₂. As a measure encompassing the effect of intermolecular interactions on dilute solution viscosity, the Huggins coefficient, k_h, is a promising candidate as a parameter measureable at low concentrations, but indicative of concentrated solution viscosity. In this study, a differential viscometry technique is developed to measure the intrinsic viscosity, [η], and the Huggins coefficient, k_h, of protein solutions. To understand the effect of colloidal protein-protein interactions on the viscosity of concentrated protein formulations, the viscometric parameters are compared to k_D and B₂₂ of two mAbs, tuning the contributions of repulsive and attractive forces to the net protein-protein interaction by adjusting solution pH and ionic strength. We find a strong correlation between the concentrated protein solution viscosity and the k_h but this was not observed for the k_D or the b₂₂, which have been previously used as indicators of high concentration viscosity. Trends observed in [η] and k_h values as a function of pH and ionic strength are rationalised in terms of protein-protein interactions.

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

The Huggins coefficient kH is a well-known metric for quantifying the increase in solution viscosity arising from intermolecular interactions in relatively dilute macromolecular solutions, and there has been much interest in this solution property in connection with developing improved antibody therapeutics. While numerous kH measurements have been reported for select monoclonal antibodies (mAbs) solutions, there has been limited study of kH in terms of the fundamental molecular interactions that determine this property. In this paper, we compare measurements of the osmotic second virial coefficient B22, a common metric of intermolecular and interparticle interaction strength, to measurements of kH for model antibody solutions. This comparison is motivated by the seminal work of Russel for hard sphere particles having a short-range "sticky" interparticle interaction, and we also compare our data with known results for uncharged flexible polymers having variable excluded volume interactions because proteins are polypeptide chains. Our observations indicate that neither the adhesive hard sphere model, a common colloidal model of globular proteins, nor the familiar uncharged flexible polymer model, an excellent model of intrinsically disordered proteins, describes the dependence of kH of these antibodies on B22. Clearly, an improved understanding of protein and ion solvation by water as well as dipole-dipole and charge-dipole effects is required to understand the significance of kH from the standpoint of fundamental protein-protein interactions. Despite shortcomings in our theoretical understanding of kH for antibody solutions, this quantity provides a useful practical measure of the strength of interprotein interactions at elevated protein concentrations that is of direct significance for the development of antibody formulations that minimize the solution viscosity.

Protein intrinsic viscosity determination with the Viscosizer TD instrument: reaching beyond the initially expected applications

Intrinsic viscosity is a key hydrodynamic parameter to understand molecular structure and hydration, as well as intramolecular interactions. Commercially available instruments measure intrinsic viscosity by recording the macromolecular mobility in a capillary. These instruments monitor Taylor dispersion using an absorbance or fluorescence detector. By design, these instruments behave like U-tube viscometers. To our knowledge, there are no studies to date showing that the Viscosizer TD instrument (Malvern-Panalytical) is able to measure the intrinsic viscosity of macromolecules. In this study, we then performed our assays on the Poly(ethylene oxide) polymer (PEO), used classically as a standard for viscometry measurements and on three model proteins: the bovine serum albumin (BSA), the bevacizumab monoclonal antibody, and the RTX Repeat Domain (RD) of the adenylate cyclase toxin of Bordetella pertussis (CyaA). The presence of P20 in the samples is critical to get reliable results. The data obtained with our in-house protocol show a strong correlation with intrinsic viscosity values obtained using conventional techniques. However, with respect to them, our measurements could be performed at relatively low concentrations, between 2 and 5 mg/ml, using only 7 µL per injection. Altogether, our results show that the Viscosizer TD instrument is able to measure intrinsic viscosities in a straightforward manner. This simple and innovative approach should give a new boost to intrinsic viscosity measurements and should reignite the interest of biophysicists, immunologists, structural biologists and other researchers for this key physicochemical parameter.

Highland games: A benchmarking exercise in predicting biophysical and drug properties of monoclonal antibodies from amino acid sequences

Biopharmaceutical product and process development do not yet take advantage of predictive computational modeling to nearly the degree seen in industries based on smaller molecules. To assess and advance progress in this area, spirited coopetition (mutually beneficial collaboration between competitors) was successfully used to motivate industrial scientists to develop, share, and compare data and methods which would normally have remained confidential. The first "Highland Games" competition was held in conjunction with the October 2018 Recovery of Biological Products Conference in Ashville, NC, with the goal of benchmarking and assessment of the ability to predict development-related properties of six antibodies from their amino acid sequences alone. Predictions included purification-influencing properties such as isoelectric point and protein A elution pH, and biophysical properties such as stability and viscosity at very high concentrations. Essential contributions were made by a large variety of individuals, including companies which consented to provide antibody amino acid sequences and test materials, volunteers who undertook the preparation and experimental characterization of these materials, and prediction teams who attempted to predict antibody properties from sequence alone. Best practices were identified and shared, and areas in which the community excels at making predictions were identified, as well as areas presenting opportunities for considerable improvement. Predictions of isoelectric point and protein A elution pH were especially good with all-prediction average errors of 0.2 and 1.6 pH unit, respectively, while predictions of some other properties were notably less good. This manuscript presents the events, methods, and results of the competition, and can serve as a tutorial and as a reference for in-house benchmarking by others. Organizations vary in their policies concerning disclosure of methods, but most managements were very cooperative with the Highland Games exercise, and considerable insight into common and best practices is available from the contributed methods. The accumulated data set will serve as a benchmarking tool for further development of in silico prediction tools.

Evolution of Interactions in the Protein Solution As Induced by Mono and Multivalent Ions

The evolution of interactions in the bovine serum albumin (BSA) protein solution on addition of mono and multivalent (di, tri and tetra) counterions has been studied using small-angle neutron scattering (SANS), dynamic light scattering (DLS) and ζ-potential measurements. It is found that in the presence of mono and divalent counterions, protein behavior can be well explained by DLVO theory, combining the contributions of screened Coulomb repulsion with the van der Waals attraction. The addition of mono or divalent salts in protein solution reduces the repulsive barrier and hence the overall interaction becomes attractive, but the system remains in one-phase for the entire concentration range of the salts, added in the system. However, contrary to DLVO theory, the protein solution undergoes a reentrant phase transition from one-phase to a two-phase system and then back to the one-phase system in the presence of tri and tetravalent counterions. The results show that tri and tetravalent (unlike mono and divalent) counterions induce short-range attraction between the protein molecules, leading to the transformation from one-phase to two-phase system. The two-phase is characterized by the fractal structure of protein aggregates. The excess condensation of these higher-valent counterions in the double layer around the BSA causes the reversal of charge of the protein molecules resulting into reentrant of the one-phase, at higher salt concentrations. The complete phase behavior with mono and multivalent ions has been explained in terms of the interplay of electrostatic repulsion and ion-induced short-range attraction between the protein molecules.

Calorimetric studies of interactions between low molecular weight salts and bovine serum albumin in water at pH values below and above the isoionic point

Isothermal titration calorimetry was used to determine the temperature and salt concentration dependence of the enthalpy of mixing, Δ_mix H, of bovine serum albumin (BSA) in aqueous buffer solutions with several low molecular weight salts. Three buffers were used: acetate (pH = 4.0), MOPS (7.2), and borate (9.2). Since the isoionic point of BSA is at pI ≈ 4.7, the net charge of BSA in acetate buffer was positive (≈ +20), while in the other two buffer solutions it was negative (≈ -15 in MOPS and ≈ -25 in borate). The majority of the recorded heat effects were exothermic, while only at pH = 9.2 a weak endothermic effect upon mixing BSA with LiCl, NaCl, and KCl was observed. For all buffer solutions the absolute values of Δ_mix H of sodium salts followed the order: NaCl < NaBr < NaNO₃ < NaI < NaSCN, which is the reverse Hofmeister series for anions. The magnitude of the effects was the largest in acetate buffer and decreased with an increasing pH value of the solution. While the effect of varying the anion of the added salts was strongly pronounced at all pH values, the effect of the cation (LiCl, NaCl, KCl, RbCl and CsCl salts) was weak. The most interesting feature of the results obtained for pH > pI was the fact that Δ_mix H were considerably more sensitive to the anion (co-ion to the net BSA charge) than to the cation species. This indicated that anions interacted quite strongly with the BSA even at pH values where the net charge of the protein was negative. We showed that Δ_mix H at high addition of salts correlated well with the enthalpy of hydration of the corresponding salt anion. This finding suggested, consistently with some previous studies, that a part of the exothermic contribution to Δ_mix H originated from the hydration changes upon the protein-salt interaction. Theoretical analysis, based on the primitive model of highly asymmetric electrolyte solutions solved within the mean spherical approximation, was used to estimate Coulomb effects upon mixing.

Microliter viscometry using a bright-field microscope: η-DDM

The rheological properties of a medium can be inferred from the Brownian motion of colloidal tracer particles using the microrheology procedure. The tracer motion can be characterized by the mean-squared displacement (MSD). It can be calculated from the intermediate scattering function determined by Differential Dynamic Microscopy (DDM). Here we show that DDM together with the empirical Cox-Merz rule is particularly suited to measure the steady-shear viscosity, i.e. the viscosity towards zero frequency, due to its ability to provide reliable information on long time and length scales and hence small frequencies. This method, η-DDM, is tested and illustrated using three different systems: Newtonian fluids (glycerol-water mixtures), colloidal suspensions (protein samples) and a viscoelastic polymer solution (aqueous poly(ethylene oxide) solution). These tests show that common lab equipment, namely a bright-field optical microscope, can be used as a convenient and reliable microliter viscometer. Because η-DDM requires much smaller sample volumes than classical rheometry, only a few microliters, it is particularly useful for biological and soft matter systems.

Membrane proteins significantly restrict exosome mobility

Exosomes are membrane nanovesicles implicated in cell-to-cell signaling in which they transfer their molecular cargo from the parent to the recipient cells. This role essentially depends on the exosomes' small size, which is the prerequisite for their rapid migration through the crowded extracellular matrix and into and out of circulation. Here we report much lower exosome mobility than expected from the size of their vesicles, implicate membrane proteins in a substantially impeded rate of migration, and suggest an approach to quantifying the impact. The broadly distributed excess hydrodynamic resistance provided by surface proteins produces a highly heterogeneous and microenvironment-dependent hindrance to exosome mobility. The implications of the findings on exosome-mediated signaling are discussed.

Cluster Size and Quinary Structure Determine the Rheological Effects of Antibody Self-Association at High Concentrations

The question of how nonspecific reversible intermolecular protein interactions affect solution rheology at high concentrations is fundamentally rooted in the translation of nanometer-scale interactions into macroscopic properties. Well-defined solutions of purified monoclonal antibodies (mAbs) provide a useful system with which to investigate the manifold intricacies of weak protein interactions at high concentrations. Recently, characterization of self-associating IgG1 antibody (mAb2) solutions has established the direct role of protein clusters on concentrated mAb rheology. Expanding on our earlier work with three additional mAbs (mAb1, mAb3, and mAb4), the observed concentration-dependent static light scattering and rheological data present a substantially more complex relationship between protein interactions and solution viscosity at high concentrations. The four mAb systems exhibited divergent correlations between cluster formation (size) and concentrated solution viscosities dependent on mAb primary sequence and solution conditions. To address this challenge, well-established features of colloidal cluster phenomena could be applied as a framework for interpreting our observations. The initial stages of mAb cluster formation were investigated with small-angle X-ray scattering (SAXS) and ensemble-optimized fit methods, to uncover shifts in the dimer structure populations which are produced by changes in mAb interaction modes and association valence under the different solution conditions. Analysis of mAb average cluster number and effective hydrodynamic radii at high concentrations revealed cluster architectures can have a wide range of fractal dimensions. Collectively, the static light scattering, SAXS, and rheological characterization demonstrate that nonspecific and anisotropic attractive intermolecular interactions produce antibody clusters with different quinary structures to regulate the rheological properties of concentrated mAb solutions.

Viscosity Control of Protein Solution by Small Solutes: A Review

Viscosity of protein solution is one of the most troublesome issues for the high-concentration formulation of protein drugs. In this review, we summarize the practical methods that suppress the viscosity of protein solution using small molecular additives. The small amount of salts decreases the viscosity that results from electrostatic repulsion and attraction. The chaotrope suppresses the hydrophobic attraction and cluster formation, which can lower the solution viscosity. Arginine hydrochloride (ArgHCl) also suppresses the solution viscosity due to the hydrophobic and aromatic interactions between protein molecules. The small molecular additives are the simplest resolution of the high viscosity of protein solution as well as understanding of the primary cause in complex phenomena of protein interactions.

Detecting protein folding by thermal fluctuations of microcantilevers

The accurate characterization of proteins in both their native and denatured states is essential to effectively understand protein function, folding and stability. As a proof of concept, a micro rheological method is applied, based on the characterization of thermal fluctuations of a micro cantilever immersed in a bovine serum albumin solution, to assess changes in the viscosity associated with modifications in the protein’s structure under the denaturant effect of urea. Through modeling the power spectrum density of the cantilever’s fluctuations over a broad frequency band, it is possible to implement a fitting procedure to accurately determine the viscosity of the fluid, even at low volumes. Increases in viscosity during the denaturant process are identified using the assumption that the protein is a hard sphere, with a hydrodynamic radius that increases during unfolding. This is modeled accordingly through the Einstein-Batchelor formula. The Einstein-Batchelor formula estimates are verified through dynamic light scattering, which measures the hydrodynamic radius of proteins. Thus, this methodology is proven to be suitable for the study of protein folding in samples of small size at vanishing shear stresses.

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

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

Rheology of clustering protein solutions

High viscosity is a major challenge with protein therapeutics at extremely high concentrations. To overcome this obstacle, it is essential to understand the relationship between the concentration of a protein solution and its viscosity as a function of shear rate and temperature. Here, lysozyme is a model charged globular protein having both short-ranged attraction (SA) and long-ranged repulsion (LR) that promote the formation of dynamic clusters at high concentrations. We report viscosity measurements from a micro-capillary rheometer (using only several microliters of solution) over a wide range of lysozyme solution concentrations, shear rates, and temperatures. Solution structural relaxation dynamics are also probed by dynamic light scattering (DLS). As a result of lysozyme's SALR interactions, the viscosity increased dramatically across all shear rates with increasing concentration and decreasing temperature. While most of the solutions exhibited Newtonian behavior, shear thinning was exhibited at the highest concentration (480 g/l) and lowest temperatures at shear rates above approximately 10(4 )s(-1). The onset shear rate for thinning and a structural relaxation rate estimated from a slow-mode measured by DLS are compared. These measurements provide insights into the properties of protein solutions and their microscopic structural origins.

Coupling of isotropic and directional interactions and its effect on phase separation and self-assembly

The interactions of molecules and particles in solution often involve an interplay between isotropic and highly directional interactions that lead to a mutual coupling of phase separation and self-assembly. This situation arises, for example, in proteins interacting through hydrophobic and charged patch regions on their surface and in nanoparticles with grafted polymer chains, such as DNA. As a minimal model of complex fluids exhibiting this interaction coupling, we investigate spherical particles having an isotropic interaction and a constellation of five attractive patches on the particle's surface. Monte Carlo simulations and mean-field calculations of the phase boundaries of this model depend strongly on the relative strength of the isotropic and patch potentials, where we surprisingly find that analytic mean-field predictions become increasingly accurate as the directional interactions become increasingly predominant. We quantitatively account for this effect by noting that the effective interaction range increases with increasing relative directional to isotropic interaction strength. We also identify thermodynamic transition lines associated with self-assembly, extract the entropy and energy of association, and characterize the resulting cluster properties obtained from simulations using percolation scaling theory and Flory-Stockmayer mean-field theory. We find that the fractal dimension and cluster size distribution are consistent with those of lattice animals, i.e., randomly branched polymers swollen by excluded volume interactions. We also identify a universal functional form for the average molecular weight and a nearly universal functional form for a scaling parameter characterizing the cluster size distribution. Since the formation of branched clusters at equilibrium is a common phenomenon in nature, we detail how our analysis can be used in experimental characterization of such associating fluids.

The effect of protein concentration on the viscosity of a recombinant albumin solution formulation

The effect of protein concentration on solution viscosity in a commercially available biopharmaceutical formulation of recombinant albumin (rAlbumin) was studied.

Critical examination of the colloidal particle model of globular proteins

Recent studies of globular protein solutions have uniformly adopted a colloidal view of proteins as particles, a perspective that neglects the polymeric primary structure of these biological macromolecules, their intrinsic flexibility, and their ability to sample a large configurational space. While the colloidal perspective often serves as a useful idealization in many cases, the macromolecular identity of proteins must reveal itself under thermodynamic conditions in which the native state is no longer stable, such as denaturing solvents and high protein concentrations where macromolecules tend to have screened excluded volume, charge, and hydrodynamic interactions. Under extreme pH conditions, charge repulsion interactions within the protein chain can overcome the attractive hydrogen-bonding interactions, holding it in its native globular state. Conformational changes can therefore be expected to have great significance on the shear viscosity and other rheological properties of protein solutions. These changes are not envisioned in conventional colloidal protein models and we have initiated an investigation of the scattering and rheological properties of model proteins. We initiate this effort by considering bovine serum albumin because it is a globular protein whose solution properties have also been extensively investigated as a function of pH, temperature, ionic strength, and concentration. As we anticipated, near-ultraviolet circular dichroism measurements and intrinsic viscosity measurements clearly indicate that the bovine serum albumin tertiary structure changes as protein concentration and pH are varied. Our findings point to limited validity of the colloidal protein model and to the need for further consideration and quantification of the effects of conformational changes on protein solution viscosity, protein association, and the phase behavior. Small-angle Neutron Scattering measurements have allowed us to assess how these conformational changes influence protein size, shape, and interprotein interaction strength.

Explaining the non-newtonian character of aggregating monoclonal antibody solutions using small-angle neutron scattering

A monoclonal antibody solution displays an increase in low shear rate viscosity upon aggregation after prolonged incubation at 40°C. The morphology and interactions leading to the formation of the aggregates responsible for this non-Newtonian character are resolved using small-angle neutron scattering. Our data show a weak repulsive barrier before proteins aggregate reversibly, unless a favorable contact with high binding energy occurs. Two types of aggregates were identified after incubation at 40°C: oligomers with radius of gyration ∼10 nm and fractal submicrometer particles formed by a slow reaction-limited aggregation process, consistent with monomers colliding many times before finding a favorable strong interaction site. Before incubation, these antibody solutions are Newtonian liquids with no increase in low shear rate viscosity and no upturn in scattering at low wavevector, whereas aggregated solutions under the same conditions have both of these features. These results demonstrate that fractal submicrometer particles are responsible for the increase in low shear rate viscosity and low wavevector upturn in scattered intensity of aggregated antibody solutions; both are removed from aggregated samples by filtering.

Aggregate structure, morphology and the effect of aggregation mechanisms on viscosity at elevated protein concentrations

Non-native aggregation is a common issue in a number of degenerative diseases and during manufacturing of protein-based therapeutics. There is a growing interest to monitor protein stability at intermediate to high protein concentrations, which are required for therapeutic dosing of subcutaneous injections. An understanding of the impact of protein structural changes and interactions on the protein aggregation mechanisms and resulting aggregate size and morphology may lead to improved strategies to reduce aggregation and solution viscosity. This report investigates non-native aggregation of a model protein, α-chymotrypsinogen, under accelerated conditions at elevated protein concentrations. Far-UV circular dichroism and Raman scattering show structural changes during aggregation. Size exclusion chromatography and laser light scattering are used to monitor the progression of aggregate growth and monomer loss. Monomer loss is concomitant with increased β-sheet structures as monomers are added to aggregates, which illustrate a transition from a native monomeric state to an aggregate state. Aggregates grow predominantly through monomer-addition, resulting in a semi-flexible polymer morphology. Analysis of aggregation growth kinetics shows that pH strongly affects the characteristic timescales for nucleation (τn) and growth (τg), while the initial protein concentration has only minor effects on τn or τg. Low-shear viscosity measurements follow a common scaling relationship between average aggregate molecular weight (Mw(agg)) and concentration (σ), which is consistent with semi-dilute polymer-solution theory. The results establish a link between aggregate growth mechanisms, which couple Mw(agg) and σ, to increases in solution viscosity even at these intermediate protein concentrations (less than 3w/v %).