Nonmonotonic variation with salt concentration of the second virial coefficient in protein solutions

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

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

Toward a Suite of Coarse-Grained Models for Molecular Simulation of Monoclonal Antibodies and Therapeutic Proteins

A series of coarse-grained models for molecular simulation of proteins are considered, with emphasis on the application of predicting protein-protein self-interactions for monoclonal antibodies (MAbs). As an illustrative example and for quantitative comparison, the models are used to predict osmotic virial coefficients over a broad range of attractive and repulsive self-interactions and solution conditions for a series of MAbs where the second osmotic virial coefficient has been experimentally determined in prior work. The models are compared based on how well they can predict experimental behavior, their computational burdens, and scalability. An intermediate-resolution model is also introduced that can capture specific electrostatic interactions with improved efficiency and similar or improved accuracy when compared to the previously published models. Guidance is included for the selection of coarse-grained models more generally for capturing a balance of electrostatic, steric, and short-ranged nonelectrostatic interactions for proteins from low to high concentrations.

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.

Anisotropy of the Coulomb interaction between folded proteins: consequences for mesoscopic aggregation of lysozyme

Toward quantitative description of protein aggregation, we develop a computationally efficient method to evaluate the potential of mean force between two folded protein molecules that allows for complete sampling of their mutual orientation. Our model is valid at moderate ionic strengths and accounts for the actual charge distribution on the surface of the molecules, the dielectric discontinuity at the protein-solvent interface, and the possibility of protonation or deprotonation of surface residues induced by the electric field due to the other protein molecule. We apply the model to the protein lysozyme, whose solutions exhibit both mesoscopic clusters of protein-rich liquid and liquid-liquid separation; the former requires that protein form complexes with typical lifetimes of approximately milliseconds. We find the electrostatic repulsion is typically lower than the prediction of the Derjaguin-Landau-Verwey-Overbeek theory. The Coulomb interaction in the lowest-energy docking configuration is nonrepulsive, despite the high positive charge on the molecules. Typical docking configurations barely involve protonation or deprotonation of surface residues. The obtained potential of mean force between folded lysozyme molecules is consistent with the location of the liquid-liquid coexistence, but produces dimers that are too short-lived for clusters to exist, suggesting lysozyme undergoes conformational changes during cluster formation.

Computational methods to predict therapeutic protein aggregation

Protein based biotherapeutics have emerged as a successful class of pharmaceuticals. However, these macromolecules endure a variety of physicochemical degradations during manufacturing, shipping, and storage, which may adversely impact the drug product quality. Of these degradations, the irreversible self-association of therapeutic proteins to form aggregates is a major challenge in the formulation of these molecules. Tools to predict and mitigate protein aggregation are, therefore, of great interest to biopharmaceutical research and development. In this chapter, a number of such computational tools developed to understand and predict the various steps involved in protein aggregation are described. These tools can be grouped into three general classes: unfolding kinetics and native state thermal stability, colloidal stability, and sequence/structure based aggregation liabilities. Chapter sections introduce each class by discussing how these predictive tools provide insight into the molecular events leading to protein aggregation. The computational methods are then explained in detail along with their advantages and limitations.

DDFT calibration and investigation of an anisotropic phase-field crystal model

The anisotropic phase-field crystal model recently proposed and used by Prieler et al (2009 J. Phys.: Condens. Matter 21 464110) is derived from microscopic density functional theory for anisotropic particles with fixed orientation. Its morphology diagram is also explored. In particular we have investigated the influence of anisotropy and undercooling on the process of nucleation and microstructure formation from the atomic to the microscale. To that end numerical simulations were performed varying those dimensionless parameters which represent anisotropy and undercooling in our anisotropic phase-field crystal model. The results from these numerical simulations are summarized in terms of a morphology diagram of the stable state phases. These stable phases are also investigated with respect to their kinetics and characteristic morphological features.

Nonadditivity in the effective interactions of binary charged colloidal suspensions

Based on primitive model computer simulations with explicit microions, we calculate the effective interactions in a binary mixture of charged colloids with species A and B for different size and charge ratios. An optimal pairwise interaction is obtained by fitting the many-body effective forces. This interaction is close to a Yukawa (or Derjaguin-Landau-Verwey-Overbeek (DLVO)) pair potential but the AB cross-interaction is different from the geometric mean of the two direct AA and BB interactions. As a function of charge asymmetry, the corresponding nonadditivity parameter is first positive, then significantly negative and is then positive again. We finally show that an inclusion of nonadditivity within an optimal effective Yukawa model gives better predictions for the fluid pair structure than DLVO theory.

Electrostatics in soft matter

Recent progress in understanding the effect of electrostatics in soft matter is presented. A vast number of materials contain ions, ranging from the molecular scale (e.g. electrolyte) to the meso/macroscopic one (e.g. charged colloidal particles or polyelectrolytes). Their (micro)structure and physico-chemical properties are especially dictated by the famous and redoubtable long-ranged Coulomb interaction. In particular, theoretical and simulational aspects, including the experimental motivations, will be discussed.

Molecular dynamics characterization of protein crystal contacts in aqueous solutions

We employ nonequilibrium molecular dynamics simulation to characterize the effective interactions between lysozyme molecules involved in the formation of two hydrophobic crystal contacts. We show that the effective interactions between crystal contacts do not exceed a few kT, the range of the attractive part of the potential is less than 4 angstroms, and, within this range, there is a significant depletion of water density between two protein contacts. Our findings highlight the different natures of protein crystallization and protein recognition processes.

Toward the description of electrostatic interactions between globular proteins: Potential of mean force in the primitive model

Monte Carlo simulations are used to calculate the exact potential of mean force between charged globular proteins in aqueous solution. The aim of the present paper is to study the influence of the ions of the added salt on the effective interaction between these nanoparticles. The charges of the model proteins, either identical or opposite, are either central or distributed on a discrete pattern. Contrarily to Poisson-Boltzmann predictions, attractive, and repulsive direct forces between proteins are not screened similarly. Moreover, it has been shown that the relative orientations of the charge patterns strongly influence salt-mediated interactions. More precisely, for short distances between the proteins, ions enhance the difference of the effective forces between (i) like-charged and oppositely charged proteins, (ii) attractive and repulsive relative orientations of the proteins, which may affect the selectivity of protein/protein recognition. Finally, such results observed with the simplest models are applied to a more elaborate one to demonstrate their generality.

Tail-induced attraction between nucleosome core particles

We study a possible electrostatic mechanism underlying the compaction of DNA inside the nuclei of eucaryotes: the tail-bridging effect between nucleosomes, the fundamental DNA packaging units of the chromatin complex. As a simple model of the nucleosome we introduce the eight-tail colloid, a charged sphere with eight oppositely charged, flexible, grafted chains that represent the terminal histone tails. We show that our complexes attract each other via the formation of chain bridges and contrast this to the effect of attraction via charge patches. We demonstrate that the attraction between eight-tail colloids can be tuned by changing the fraction of charged monomers on the tails. This suggests a physical mechanism of chromatin compaction where the degree of DNA condensation is controlled via biochemical means, namely the acetylation and deacetylation of lysines in the histone tails.

The nucleosome: a transparent, slippery, sticky and yet stable DNA-protein complex

Roughly three quarters of eucaryotic DNA are tightly wrapped onto protein cylinders organized in so-called nucleosomes. Despite this fact, the wrapped DNA cannot be inert since DNA is at the heart of many crucial life processes. We focus here on physical mechanisms that might allow nucleosomes to perform a great deal of such processes, specifically 1) on unwrapping fluctuations that give DNA-binding proteins access to the wrapped DNA portions without disrupting the nucleosome as a whole, 2) on corkscrew sliding along DNA and some implications and on 3) tail-bridging-induced attraction between nucleosomes as a means of controlling higher-order folding.

Second virial coefficient determination of a therapeutic peptide by self-interaction chromatography

Self-interaction of macromolecules has been shown to play an important role in a number of physical processes, including crystallization, solubility, viscosity, and aggregation. Peptide self-interaction is not as well studied as for larger proteins, but should play an equally important role. The osmotic second virial coefficient, B, can be used to quantify peptide and protein self-interaction. B values are typically measured using static light scattering (SLS). Peptides, however, do not scatter enough light to allow such measurements. This study describes the first use of self-interaction chromatography (SIC) for the measurement of peptide B values because SIC does not have the molecular size limitations of SLS. In the present work, SIC was used to measure B for enfuvirtide, a 36-amino acid therapeutic peptide, as a function of salt concentration, salt type, and pH. B was found to correlate strongly with solubility and apparent molecular weight. In general, the solubility of enfuvirtide increases with pH from 6 to 10 and decreases as the salt concentration increases from 0 to 0.5M for three different salts. The effect of peptide concentration on B was also investigated and shown to have a significant effect, but only at high concentrations (>80 mg/mL).

Diffusion and microstructural properties of solutions of charged nanosized proteins: Experiment versus theory

We have reanalyzed our former static small-angle x-ray scattering and photon correlation spectroscopy results on dense solutions of charged spherical apoferritin proteins using theories recently developed for studies of colloids. The static structure factors S(q), and the small-wave-number collective diffusion coefficient D(c) determined from those experiments are interpreted now in terms of a theoretical scheme based on a Derjaguin-Landau-Verwey-Overbeek-type continuum model of charged colloidal spheres. This scheme accounts, in an approximate way, for many-body hydrodynamic interactions. Stokesian dynamics computer simulations of the hydrodynamic function have been performed for the first time for dense charge-stabilized dispersions to assess the accuracy of the theoretical scheme. We show that the continuum model allows for a consistent description of all experimental results, and that the effective particle charge is dependent upon the protein concentration relative to the added salt concentration. In addition, we discuss the consequences of small ions dynamics for the collective protein diffusion within the framework of the coupled-mode theory.

Monte Carlo simulation of an inhomogeneous dielectric continuum model for B-DNA

Thermodynamic and structural properties of the counterion atmosphere surrounding B-DNA are calculated by Monte Carlo simulation in a spatially inhomogeneous, but piecewise uniform, dielectric continuum cell model - the "barbarous" model. A boundary element formulation is implemented to study the sensitivity of these properties with respect to perturbations in the location of discontinuous dielectric boundaries relative to fixed and mobile charges. High concentrations are considered corresponding to the liquid crystalline hexagonally ordered phase of DNA. Primitive model results are verified against other simulation reports and a comparison of barbarous model predictions with experimental data is discussed. The internal energy, osmotic coefficient, radial distributions and the population ratio of counterions in the geometrically resolved major and minor grooves are all found to strongly depend on the dielectric boundary position. This suggests that a self-consistent development of the model should consider a free surface problem where the boundary is not specified a priori.

Thermodynamic behaviors of polyampholytes at low temperatures

The thermodynamic properties of 6-mer polyampholyte chains were simulated with the multicanonical Monte Carlo method in a three-dimensional off-lattice model. The results demonstrate that the number of the charged monomers on the chain and the sequence significantly affect the transition behaviors of polyampholyte chains. In the case of the random distribution of the charged monomers, two oppositely charged monomers act as two stickers during a collapse process, due to which the movement of monomers is suppressed and the chains can more directly reach their solid state. Interestingly, in a proper arrangement of the charged monomers, the chains may have a proteinlike nondegenerate state, i.e., the chain is designable. The transition of such polyampholyte chains appears to be of a two-state, all-or-none type.

Phase behavior of aqueous solutions containing dipolar proteins from second-order perturbation theory

Due to the interplay of Coulombic repulsion and attractive dipolar and van der Waals interactions, solutions of globular proteins display a rich variety of phase behavior featuring fluid-fluid and fluid-solid transitions that strongly depend on solution pH and salt concentration. Using a simple model for charge, dispersion and dipole-related contributions to the interprotein potential, we calculate phase diagrams for protein solutions within the framework of second-order perturbation theory. For each phase, we determine the Helmholtz energy as the sum of a hard-sphere reference term and a perturbation term that reflects both the electrostatic and dispersion interactions. Dipolar effects can induce fluid-fluid phase separation or crystallization even in the absence of any significant dispersion attraction. Because dissolved electrolytes screen the charge-charge repulsion more strongly than the dipolar attraction, the ionic strength dependence of the potential of mean force can feature a minimum at intermediate ionic strengths offering an explanation for the observed nonmonotonic dependence of the phase behavior on salt concentration. Inclusion of correlations between charge-dipole and dipole-dipole interactions is essential for a reliable calculation of phase diagrams for systems containing charged dipolar proteins and colloids.

Attraction between DNA molecules mediated by multivalent ions

The effective force between two parallel DNA molecules is calculated as a function of their mutual separation for different valencies of counterion and salt ions and different salt concentrations. Computer simulations of the primitive model are used and the shape of the DNA molecules is accurately modeled using different geometrical shapes. We find that multivalent ions induce a significant attraction between the DNA molecules whose strength can be tuned by the averaged valency of the ions. The physical origin of the attraction is traced back either to electrostatics or to entropic contributions. For multivalent counterions and monovalent salt ions, we find a salt-enhanced repulsion effect: the force is first attractive but gets repulsive with increasing salt concentration. Furthermore, we show that the multivalent-ion-induced attraction does not necessarily correlate with DNA overcharging.

Thermodynamic instability in supersaturated lysozyme solutions: effect of salt and role of concentration fluctuations

Experimental and theoretical work has suggested that protein crystal nucleation can be affected by the separation of two metastable liquid phases with different local concentrations, or more specifically by critical density fluctuations. We measure the amplitude and correlation length of local concentration fluctuations by light scattering for supersaturated solutions of hen egg-white lysozyme (at pH 4.5 and at different NaCl concentrations, up to 7% w/v). By extrapolating the critical divergent behavior of concentration fluctuation amplitude versus temperature, we determine the spinodal line, that is the limit of stability. Cloud-point measurements are used to determine liquid-liquid coexistence, consistent with previous work. In the present work, which is an extensive study of off-critical fluctuations in supersaturated protein solution, we observe a nonclassical scaling divergent behavior of the correlation length of concentration fluctuations, thus suggesting that off-critical fluctuations may have a role in crystallization kinetics. To appropriately fit the spinodal data, an entropic term must be added to the van der Waals or to the adhesive hard-sphere model. We interpret this contribution as due to the salt-induced modulation of protein hydration.