Critical examination of the colloidal particle model of globular proteins

Do Ionic Liquids Exhibit the Required Characteristics to Dissolve, Extract, Stabilize, and Purify Proteins? Past-Present-Future Assessment

Proteins are highly labile molecules, thus requiring the presence of appropriate solvents and excipients in their liquid milieu to keep their stability and biological activity. In this field, ionic liquids (ILs) have gained momentum in the past years, with a relevant number of works reporting their successful use to dissolve, stabilize, extract, and purify proteins. Different approaches in protein-IL systems have been reported, namely, proteins dissolved in (i) neat ILs, (ii) ILs as co-solvents, (iii) ILs as adjuvants, (iv) ILs as surfactants, (v) ILs as phase-forming components of aqueous biphasic systems, and (vi) IL-polymer-protein/peptide conjugates. Herein, we critically analyze the works published to date and provide a comprehensive understanding of the IL-protein interactions affecting the stability, conformational alteration, unfolding, misfolding, and refolding of proteins while providing directions for future studies in view of imminent applications. Overall, it has been found that the stability or purification of proteins by ILs is bispecific and depends on the structure of both the IL and the protein. The most promising IL-protein systems are identified, which is valuable when foreseeing market applications of ILs, e.g., in "protein packaging" and "detergent applications". Future directions and other possibilities of IL-protein systems in light-harvesting and biotechnology/biomedical applications are discussed.

Changes in Secondary Structure and Properties of Bovine Serum Albumin as a Result of Interactions with Gold Surface

Proteins can alter their shape when interacting with a surface. This study explores how bovine serum albumin (BSA) modifies structurally when it adheres to a gold surface, depending on the protein concentration and pH. We verified that the gold surface induces significant structural modifications to the BSA molecule using circular dichroism, infrared spectroscopy, and atomic force microscopy. Specifically, adsorbed molecules displayed increased levels of disordered structures and β-turns, with fewer α-helices than the native structure. MP-SPR spectroscopy demonstrated that the protein molecules preferred a planar orientation during adsorption. Molecular dynamics simulations revealed that the interaction between cysteines exposed to the outside of the molecule and the gold surface was vital, especially at pH=3.5. The macroscopic properties of the protein film observed by AFM and contact angles confirm the flexible nature of the protein itself. Notably, structural transformation is joined with the degree of hydration of protein layers.

Protein Association in Solution: Statistical Mechanical Modeling

Protein molecules associate in solution, often in clusters beyond pairwise, leading to liquid phase separations and high viscosities. It is often impractical to study these multi-protein systems by atomistic computer simulations, particularly in multi-component solvents. Instead, their forces and states can be studied by liquid state statistical mechanics. However, past such approaches, such as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, were limited to modeling proteins as spheres, and contained no microscopic structure-property relations. Recently, this limitation has been partly overcome by bringing the powerful Wertheim theory of associating molecules to bear on protein association equilibria. Here, we review these developments.

Ultrathin Micellar Foam Films of Sodium Caseinate Protein Solutions

Sodium caseinates (NaCas), derived from milk proteins called caseins, are often added to food formulations as emulsifiers, foaming agents, and ingredients for producing dairy products. In this contribution, we contrast the drainage behavior of single foam films made with micellar NaCas solutions with well-established features of stratification observed for the micellar sodium dodecyl sulfate (SDS) foam films. In reflected light microscopy, the stratified SDS foam films display regions with distinct gray colors due to differences in interference intensity from coexisting thick-thin regions. Using IDIOM (interferometry digital imaging optical microscopy) protocols we pioneered for mapping nanotopography of foam films, we showed that drainage via stratification in SDS films proceeds by the expansion of flat domains that are thinner than surrounding by a concentration-dependent step-size, and nonflat features (nanoridges and mesas) form at the moving front. Furthermore, stratifying SDS foam films show stepwise thinning, such that the step-size and terminal film thickness decrease with concentration. Here we visualize the nanotopography in protein films with high spatiotemporal resolution using IDIOM protocols to address two long-standing questions. Do protein foam films formulated with NaCas undergo drainage via stratification? Are thickness transitions and variations in protein foam films determined by intermicellar interactions and supramolecular oscillatory disjoining pressure? In contrast with foam films containing micellar SDS, we find that micellar NaCas foam films display just one step, nonflat and noncircular domains that expand without forming nanoridges and a terminal thickness that increases with NaCas concentration. We infer that the differences in adsorbing and self-assembling unimers triumph over any similarities in the structure and interactions of their micelles.

Modelling network formation in folded protein hydrogels by cluster aggregation kinetics

Globular folded protein-based hydrogels are becoming increasingly attractive due to their specific biological functionality, as well as their responsiveness to stimuli. By modelling folded proteins as colloids, there are rich opportunities to explore network formation mechanisms in protein hydrogels that negate the need for computationally expensive simulations which capture the full complexity of proteins. Here we present a kinetic lattice-based model which simulates the formation of irreversibly chemically crosslinked, folded protein-based hydrogels. We identify the critical point of gel percolation, explore the range of network regimes covering diffusion-limited to reaction-limited cluster aggregation (DLCA and RLCA, respectively) network formation mechanisms and predict the final network structure, fractal dimensions and final gel porosity. We reveal a crossover between DLCA and RLCA mechanisms as a function of protein volume fraction and show how the final network structure is governed by the structure at the percolation point, regardless of the broad variation of non-percolating cluster masses observed across all systems. An analysis of the pore size distribution in the final network structures reveals that, approaching RLCA, gels have larger maximal pores than the DLCA counterparts for both volume fractions studied. This general kinetic model and the analysis tools generate predictions of network structure and concurrent porosity over a broad range of experimentally controllable parameters that are consistent with current expectations and understanding of experimental results.

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.

Diffusive Dynamics of Bacterial Proteome as a Proxy of Cell Death

Temperature variations have a big impact on bacterial metabolism and death, yet an exhaustive molecular picture of these processes is still missing. For instance, whether thermal death is determined by the deterioration of the whole or a specific part of the proteome is hotly debated. Here, by monitoring the proteome dynamics of E. coli, we clearly show that only a minor fraction of the proteome unfolds at the cell death. First, we prove that the dynamical state of the E. coli proteome is an excellent proxy for temperature-dependent bacterial metabolism and death. The proteome diffusive dynamics peaks at about the bacterial optimal growth temperature, then a dramatic dynamical slowdown is observed that starts just below the cell's death temperature. Next, we show that this slowdown is caused by the unfolding of just a small fraction of proteins that establish an entangling interprotein network, dominated by hydrophobic interactions, across the cytoplasm. Finally, the deduced progress of the proteome unfolding and its diffusive dynamics are both key to correctly reproduce the E. coli growth rate.

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.

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

Unusual stability of protein molecules in the presence of multivalent counterions

Proteins are known to undergo denaturation and form different phases with varying physicochemical parameters. We report unusual stability of bovine serum albumin protein against commonly used denaturants (temperature and surfactant) in the charged reversal reentrant phase, caused by the multivalent counterions. Unlike monovalent counterions, which promote the denaturants' induced protein unfolding, the unfolding is restricted in the presence of multivalent ions. The observations are beyond the scope of general understanding of protein unfolding and are believed to be governed by ion-ion correlations driven strong condensation of the multivalent ions.

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

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

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

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

Biophysics of Phase Separation of Disordered Proteins Is Governed by Balance between Short- And Long-Range Interactions

Intrinsically disordered proteins play a crucial role in cellular phase separation, yet the diverse molecular forces driving phase separation are not fully understood. It is of utmost importance to understand how peptide sequence, and particularly the balance between the peptides’ short- and long-range interactions with other peptides, may affect the stability, structure, and dynamics of liquid–liquid phase separation in protein condensates. Here, using coarse-grained molecular dynamics simulations, we studied the liquid properties of the condensate in a series of polymers in which the ratio of short-range dispersion interactions to long-range electrostatic interactions varied. As the fraction of mutations that participate in short-range interactions increases at the expense of long-range electrostatic interactions, a significant decrease in the critical temperature of phase separation is observed. Nevertheless, sequences with a high fraction of short-range interactions exhibit stabilization, which suggests compensation for the loss of long-range electrostatic interactions. Decreased condensate stability is coupled with decreased translational diffusion of the polymers in the condensate, which may result in the loss of liquid characteristics in the presence of a high fraction of uncharged residues. The effect of exchanging long-range electrostatic interactions for short-range interactions can be explained by the kinetics of breaking intermolecular contacts with neighboring polymers and the kinetics of intramolecular fluctuations. While both time scales are coupled and increase as electrostatic interactions are lost, for sequences that are dominated by short-range interactions, the kinetics of intermolecular contact breakage significantly slows down. Our study supports the contention that different types of interactions can maintain protein condensates, however, long-range electrostatic interactions enhance its liquid-like behavior.

Amyloid Self-Assembly of Lysozyme in Self-Crowded Conditions: The Formation of a Protein Oligomer Hydrogel

A method is designed to quickly form protein hydrogels, based on the self-assembly of highly concentrated lysozyme solutions in acidic conditions. Their properties can be easily modulated by selecting the curing temperature. Molecular insights on the gelation pathway, derived by in situ FTIR spectroscopy, are related to calorimetric and rheological results, providing a consistent picture on structure–property correlations. In these self-crowded samples, the thermal unfolding induces the rapid formation of amyloid aggregates, leading to temperature-dependent quasi-stationary levels of antiparallel cross β-sheet links, attributed to kinetically trapped oligomers. Upon subsequent cooling, thermoreversible hydrogels develop by the formation of interoligomer contacts. Through heating/cooling cycles, the starting solutions can be largely recovered back, due to oligomer-to-monomer dissociation and refolding. Overall, transparent protein hydrogels can be easily formed in self-crowding conditions and their properties explained, considering the formation of interconnected amyloid oligomers. This type of biomaterial might be relevant in different fields, along with analogous systems of a fibrillar nature more commonly considered.

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.

Model for disordered proteins with strongly sequence-dependent liquid phase behavior

Phase separation of intrinsically disordered proteins is important for the formation of membraneless organelles or biomolecular condensates, which play key roles in the regulation of biochemical processes within cells. In this work, we investigated the phase separation of different sequences of a coarse-grained model for intrinsically disordered proteins and discovered a surprisingly rich phase behavior. We studied both the fraction of total hydrophobic parts and the distribution of hydrophobic parts. Not surprisingly, sequences with larger hydrophobic fractions showed conventional liquid-liquid phase separation. The location of the critical point was systematically influenced by the terminal beads of the sequence due to changes in interfacial composition and tension. For sequences with lower hydrophobicity, we observed not only conventional liquid-liquid phase separation but also re-entrant phase behavior in which the liquid phase density decreases at lower temperatures. For some sequences, we observed the formation of open phases consisting of aggregates, rather than a normal liquid. These aggregates had overall lower densities than the conventional liquid phases and exhibited complex geometries with large interconnected string-like or membrane-like clusters. Our findings suggest that minor alterations in the ordering of residues may lead to large changes in the phase behavior of the protein, a fact of significant potential relevance for biology.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.

Downstream Processing Technologies/Capturing and Final Purification : Opportunities for Innovation, Change, and Improvement. A Review of Downstream Processing Developments in Protein Purification

Increased pressure on upstream processes to maximize productivity has been crowned with great success, although at the cost of shifting the bottleneck to purification. As drivers were economical, focus is on now on debottlenecking downstream processes as the main drivers of high manufacturing cost. Devising a holistically efficient and economical process remains a key challenge. Traditional and emerging protein purification strategies with particular emphasis on methodologies implemented for the production of recombinant proteins of biopharmaceutical importance are reviewed. The breadth of innovation is addressed, as well as the challenges the industry faces today, with an eye to remaining impartial, fair, and balanced. In addition, the scope encompasses both chromatographic and non-chromatographic separations directed at the purification of proteins, with a strong emphasis on antibodies. Complete solutions such as integrated USP/DSP strategies (i.e., continuous processing) are discussed as well as gains in data quantity and quality arising from automation and high-throughput screening (HTS). Best practices and advantages through design of experiments (DOE) to access a complex design space such as multi-modal chromatography are reviewed with an outlook on potential future trends. A discussion of single-use technology, its impact and opportunities for further growth, and the exciting developments in modeling and simulation of DSP rounds out the overview. Lastly, emerging trends such as 3D printing and nanotechnology are covered. Graphical Abstract Workflow of high-throughput screening, design of experiments, and high-throughput analytics to understand design space and design space boundaries quickly. (Reproduced with permission from Gregory Barker, Process Development, Bristol-Myers Squibb).

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.

Soft-Matter Approaches for Controlling Food Protein Interactions and Assembly

Animal- and plant-based proteins are present in a wide variety of raw and processed foods. They play an important role in determining the final structure of food matrices. Food proteins are diverse in terms of their biological origin, molecular structure, and supramolecular assembly. This diversity has led to segmented experimental studies that typically focus on one or two proteins but hinder a more general understanding of food protein structuring as a whole. In this review, we propose a unified view of how soft-matter physics can be used to control food protein assembly. We discuss physical models from polymer and colloidal science that best describe and predict the phase behavior of proteins. We explore the occurrence of phase transitions along two axes: increasing protein concentration and increasing molecular attraction. This review provides new perspectives on the link between the interactions, phase transitions, and assembly of proteins that can help in designing new food products and innovative food processing operations.

Modeling the depletion effect caused by an addition of polymer to monoclonal antibody solutions

We present a theoretical study of colloidal stability of the model mixtures of monoclonal antibody molecules and non-adsorbing (no polymer-protein attraction) polymers. The antibodies are pictured as an assembly of seven hard spheres assuming a Y-like shape. Polymers present in the mixture are modeled as chain-like molecules having from 32 up to 128 monomers represented as hard spheres. We use Wertheim's thermodynamic perturbation theory to construct the two molecular species and to calculate measurable properties. The calculations are performed in the osmotic ensemble. In view that no direct attractive interaction is present in the model Hamiltonian, we only account for the entropic contribution to the phase equilibrium. We calculate chemical potentials and the equation of state for the model mixture to determine the liquid-liquid part of the phase diagram. We investigate how the critical antibody number density depends on the degree of polymerization and the bead size ratio of the polymer and protein components. The model mixture qualitatively correctly predicts some basic features of real systems. The effects of the model 'protein' geometry, that is the difference in results for the flexible Y-shaped protein versus the rigid spherical one, are also examined.

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.

Controlling the viscosities of antibody solutions through control of their binding sites

For biotechnological drugs, it is desirable to formulate antibody solutions with low viscosities. We go beyond previous colloid theories in treating protein-protein self-association of molecules that are antibody-shaped and flexible and have spatially specific binding sites. We consider interactions either through fragment antigen (Fab-Fab) or fragment crystalizable (Fab-Fc) binding. Wertheim's theory is adapted to compute the cluster-size distributions, viscosities, second virial coefficients, and Huggins coefficients, as functions of antibody concentration. We find that the aggregation properties of concentrated solutions can be anticipated from simpler-to-measure dilute solutions. A principal finding is that aggregation is controllable, in principle, through modifying the antibody itself, and not just the solution it is dissolved in. In particular: (i) monospecific antibodies having two identical Fab arms can form linear chains with intermediate viscosities. (ii) Bispecific antibodies having different Fab arms can, in some cases, only dimerize, having low viscosities. (iii) Arm-to-Fc binding allows for three binding partners, leading to networks and high viscosities.

A methodology to calculate small-angle scattering profiles of macromolecular solutions from molecular simulations in the grand-canonical ensemble

The theoretical framework to evaluate small-angle scattering (SAS) profiles for multi-component macromolecular solutions is re-examined from the standpoint of molecular simulations in the grand-canonical ensemble, where the chemical potentials of all species in solution are fixed. This statistical mechanical ensemble resembles more closely scattering experiments, capturing concentration fluctuations that arise from the exchange of molecules between the scattering volume and the bulk solution. The resulting grand-canonical expression relates scattering intensities to the different intra- and intermolecular pair distribution functions, as well as to the distribution of molecular concentrations on the scattering volume. This formulation represents a generalized expression that encompasses most of the existing methods to evaluate SAS profiles from molecular simulations. The grand-canonical SAS methodology is probed for a series of different implicit-solvent, homogeneous systems at conditions ranging from dilute to concentrated. These systems consist of spherical colloids, dumbbell particles, and highly flexible polymer chains. Comparison of the resulting SAS curves against classical methodologies based on either theoretical approaches or canonical simulations (i.e., at a fixed number of molecules) shows equivalence between the different scattering intensities so long as interactions between molecules are net repulsive or weakly attractive. On the other hand, for strongly attractive interactions, grand-canonical SAS profiles deviate in the low- and intermediate-q range from those calculated in a canonical ensemble. Such differences are due to the distribution of molecules becoming asymmetric, which yields a higher contribution from configurations with molecular concentrations larger than the nominal value. Additionally, for flexible systems, explicit discrimination between intra- and inter-molecular SAS contributions permits the implementation of model-free, structural analysis such as Guinier's plots at high molecular concentrations, beyond what the traditional limits are for such analysis.

Theory for the Liquid-Liquid Phase Separation in Aqueous Antibody Solutions

This study presents the theory for liquid-liquid phase separation for systems of molecules modeling monoclonal antibodies. Individual molecule is depicted as an assembly of seven hard spheres, organized to mimic the Y-shaped antibody. We consider the antibody-antibody interactions either through Fab, Fab' (two Fab fragments may be different), or Fc domain. Interaction between these three domains of the molecule (hereafter denoted as A, B, and C, respectively) is modeled by a short-range square-well attraction. To obtain numerical results for the model under study, we adapt Wertheim's thermodynamic perturbation theory. We use this model to calculate the liquid-liquid phase separation curve and the second virial coefficient B2. Various interaction scenarios are examined to see how the strength of the site-site interactions and their range shape the coexistence curve. In the asymmetric case, where an attraction between two sites is favored and the interaction energies for the other sites kept constant, critical temperature first increases and than strongly decreases. Some more microscopic information, for example, the probability for the particular two sites to be connected, has been calculated. Analysis of the experimental liquid-liquid phase diagrams, obtained from literature, is presented. In addition, we calculate the second virial coefficient under conditions leading to the liquid-liquid phase separation and present this quantity on the graph B2 versus protein concentration.

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.

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.

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

Dynamic light scattering can be used to measure the diffusivity of a protein within a formulation. The dependence of molecular diffusivity on protein concentration (traditionally expressed in terms of the interaction parameter k_D) is often used to infer whether protein-protein interactions are repulsive or attractive, resulting in solutions that are colloidally stable or unstable, respectively. However, a number of factors unrelated to intermolecular forces can also impact protein diffusion, complicating this interpretation. Here, we investigate the influence of multicomponent diffusion in a ternary protein-salt-water system on protein diffusion and k_D in the context of Nernst-Planck theory. This analysis demonstrates that large changes in protein diffusivity with protein concentration can result even for hard-sphere systems in the absence of protein-protein interactions. In addition, we show that dynamic light scattering measurements of diffusivity made at low ionic strength cannot be reliably used to detect protein conformational changes. We recommend comparing experimentally determined k_D values to theoretically predicted excluded-volume contributions, which will allow a more accurate assessment of protein-protein interactions.

Explicit-water theory for the salt-specific effects and Hofmeister series in protein solutions

Effects of addition of salts on stability of aqueous protein solutions are studied theoretically and the results are compared with experimental data. In our approach, all the interacting species, proteins, ions, and water molecules, are accounted for explicitly. Water molecules are modeled as hard spheres with four off-center attractive square-well sites. These sites serve to bind either another water or to solvate the ions or protein charges. The ions are represented as charged hard spheres, and decorated by attractive sites to allow solvation. Spherical proteins simultaneously possess positive and negative groups, represented by charged hard spheres, attached to the surface of the protein. The attractive square-well sites, mimicking the protein-protein van der Waals interaction, are located on the surface of the protein. To obtain numerical results, we utilized the energy route of Wertheim's associative mean spherical approximation. From measurable properties, we choose to calculate the second virial coefficient B2, which is closely related to the tendency of proteins to aggregate and eventually crystalize. Calculations are in agreement with experimental trends: (i) For low concentration of added salt, the alkali halide salts follow the inverse Hofmeister series. (ii) At higher concentration of added salt, the trend is reversed. (iii) When cations are varied, the salts follow the direct Hofmeister series. (iv) In contrast to the colloidal theories, our approach correctly predicts the non-monotonic behavior of B2 upon addition of salts. (v) With respect to anions, the theory predicts for the B2 values to follow different sequences below and above the iso-ionic point, as also confirmed experimentally. (vi) A semi-quantitative agreement between measured and calculated values for the second virial coefficient, as functions of pH of solution and added salt type and concentration, is obtained.

In-silico prediction of concentration-dependent viscosity curves for monoclonal antibody solutions

Early stage developability assessments of monoclonal antibody (mAb) candidates can help reduce risks and costs associated with their product development. Forecasting viscosity of highly concentrated mAb solutions is an important aspect of such developability assessments. Reliable predictions of concentration-dependent viscosity behaviors for mAb solutions in platform formulations can help screen or optimize drug candidates for flexible manufacturing and drug delivery options. Here, we present a computational method to predict concentration-dependent viscosity curves for mAbs solely from their sequence-structural attributes. This method was developed using experimental data on 16 different mAbs whose concentration-dependent viscosity curves were experimentally obtained under standardized conditions. Each concentration-dependent viscosity curve was fitted with a straight line, via logarithmic manipulations, and the values for intercept and slope were obtained. Intercept, which relates to antibody diffusivity, was found to be nearly constant. In contrast, slope, the rate of increase in solution viscosity with solute concentration, varied significantly across different mAbs, demonstrating the importance of intermolecular interactions toward viscosity. Next, several molecular descriptors for electrostatic and hydrophobic properties of the 16 mAbs derived using their full-length homology models were examined for potential correlations with the slope. An equation consisting of hydrophobic surface area of full-length antibody and charges on V_H, V_L, and hinge regions was found to be capable of predicting the concentration-dependent viscosity curves of the antibody solutions. Availability of this computational tool may facilitate material-free high-throughput screening of antibody candidates during early stages of drug discovery and development.

Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules

Artistic representation of limited valance units consisting of a soft core (in blue) and a small number of flexible bonding patches (in orange).

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.