Effects of additives on surfactant phase behavior relevant to bacteriorhodopsin crystallization

The interactions leading to crystallization of the integral membrane protein bacteriorhodopsin solubilized in n-octyl-beta-D-glucoside were investigated. Osmotic second virial coefficients (B(22)) were measured by self-interaction chromatography using a wide range of additives and precipitants, including polyethylene glycol (PEG) and heptane-1,2,3-triol (HT). In all cases, attractive protein-detergent complex (PDC) interactions were observed near the surfactant cloud point temperature, and there is a correlation between the surfactant cloud point temperatures and PDC B(22) values. Light scattering, isothermal titration calorimetry, and tensiometry reveal that although the underlying reasons for the patterns of interaction may be different for various combinations of precipitants and additives, surfactant phase behavior plays an important role in promoting crystallization. In most cases, solution conditions that led to crystallization fell within a similar range of slightly negative B(22) values, suggesting that weakly attractive interactions are important as they are for soluble proteins. However, the sensitivity of the cloud point temperatures and resultant coexistence curves varied significantly as a function of precipitant type, which suggests that different types of forces are involved in driving phase separation depending on the precipitant used.     

Static light scattering studies of OmpF porin: implications for integral membrane protein crystallization

Integral membrane proteins carry out some of the most important functions of living cells, yet relatively few details are known about their structures. This is due, in large part, to the difficulties associated with preparing membrane protein crystals suitable for X-ray diffraction analysis. Mechanistic studies of membrane protein crystallization may provide insights that will aid in determining future membrane protein structures. Accordingly, the solution behavior of the bacterial outer membrane protein OmpF porin was studied by static light scattering under conditions favorable for crystal growth. The second osmotic virial coefficient (B22) was found to be a predictor of the crystallization behavior of porin, as has previously been found for soluble proteins. Both tetragonal and trigonal porin crystals were found to form only within a narrow window of B22 values located at approximately -0.5 to -2 X 10(-4) mol mL g(-2), which is similar to the "crystallization slot" observed for soluble proteins. The B22 behavior of protein-free detergent micelles proved very similar to that of porin-detergent complexes, suggesting that the detergent's contribution dominates the behavior of protein-detergent complexes under crystallizing conditions. This observation implies that, for any given detergent, it may be possible to construct membrane protein crystallization screens of general utility by manipulating the solution properties so as to drive detergent B22 values into the crystallization slot. Such screens would limit the screening effort to the detergent systems most likely to yield crystals, thereby minimizing protein requirements and improving productivity.     

Osmotic second virial cross‐coefficient measurements for binary combination of lysozyme,ovalbumin, and α‐amylase in salt solutions

Interactions measurement is a valuable tool to predict equilibrium phase separation of a desired protein in the presence of unwanted macromolecules. In this study, cross‐interactions were measured as the osmotic second virial cross‐coefficients (B₂₃) for the three binary protein systems involving lysozyme, ovalbumin, and α‐amylase in salt solutions (sodium chloride and ammonium sulfate). They were correlated with solubility for the binary protein mixtures. The cross‐interaction behavior at different salt concentrations was interpreted by either electrostatic or hydrophobic interaction forces. At low salt concentrations, the protein surface charge dominates cross‐interaction behavior as a function of pH. With added ovalbumin, the lysozyme solubility decreased linearly at low salt concentration in sodium chloride and increased at high salt concentration in ammonium sulfate. The B₂₃ value was found to be proportional to the slope of the lysozyme solubility against ovalbumin concentration and the correlation was explained by preferential interaction theory. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1203–1211, 2013     

Correlation of second virial coefficient with solubility for proteins in salt solutions

In this work, osmotic second virial coefficients (B(22)) were determined and correlated with the measured solubilities for the proteins, α-amylase, ovalbumin, and lysozyme. The B(22) values and solubilities were determined in similar solution conditions using two salts, sodium chloride and ammonium sulfate in an acidic pH range. An overall decrease in the solubility of the proteins (salting out) was observed at high concentrations of ammonium sulfate and sodium chloride solutions. However, for α-amylase, salting-in behavior was also observed in low concentration sodium chloride solutions. In ammonium sulfate solutions, the B(22) are small and close to zero below 2.4 M. As the ammonium sulfate concentrations were further increased, B(22) values decreased for all systems studied. The effect of sodium chloride on B(22) varies with concentration, solution pH, and the type of protein studied. Theoretical models show a reasonable fit to the experimental derived data of B(22) and solubility. B(22) is also directly proportional to the logarithm of the solubility values for individual proteins in salt solutions, so the log-linear empirical models developed in this work can also be used to rapidly predict solubility and B(22) values for given protein-salt systems.     

The influence of molecular variation on protein interactions

We investigated the influence of solvation forces on protein-protein interactions for two forms of lysozyme: hen egg white (HEWL) and turkey egg white (TEWL). Turkey egg white has more surface exposed hydrophobic residues than HEWL and the protein-protein interactions of TEWL are shown to be more attractive than those of HEWL, for the conditions studied. The importance of including a solvation term in the potential of mean force model, to account for molecular variation in protein surface characteristics, is highlighted. We also show that the magnitude of this solvation term can be estimated using readily available data.     

Protein oligomerization through domain swapping: role of inter-molecular interactions and protein concentration

Domain swapping has been shown to be an important mechanism controlling multiprotein assembly and has been suggested recently as a possible mechanism underlying protein aggregation. Understanding oligomerization via domain swapping is therefore of theoretical and practical importance. By using a symmetrized structure-based (Gō) model, we demonstrate that in the free-energy landscape of domain swapping, a large free-energy barrier separates monomeric and domain-swapped dimeric configurations. We investigate the effect of finite monomer concentration, by implementing a new semi-analytical method, which involves computing the second virial coefficient, a thermodynamic indicator of inter-molecular interactions. This method, together with the symmetrized structure-based (Gō) model, minimizes the need for expensive many-protein simulations, providing a convenient framework to investigate concentration effect. Finally, we perform direct simulations of domain-swapped trimer formation, showing that this modeling approach can be used for higher-order oligomers.     

Interactions of lysozyme in guanidinium chloride solutions from static and dynamic light-scattering measurements

The interactions of partially unfolded proteins provide insight into protein folding and protein aggregation. In this work, we studied partially unfolded hen egg lysozyme interactions in solutions containing up to 7 M guanidinium chloride (GdnHCl). The osmotic second virial coefficient (B(22)) of lysozyme was measured using static light scattering in GdnHCl aqueous solutions at 20 degrees C and pH 4.5. B(22) is positive in all solutions, indicating repulsive protein-protein interactions. At low GdnHCl concentrations, B(22) decreases with rising ionic strength: in the absence of GdnHCl, B(22) is 1.1 x 10(-3) mLmol/g(2), decreasing to 3.0 x 10(-5) mLmol/g(2) in the presence of 1 M GdnHCl. Lysozyme unfolds in solutions at GdnHCl concentrations higher than 3 M. Under such conditions, B(22) increases with ionic strength, reaching 8.0 x 10(-4) mLmol/g(2) at 6.5 M GdnHCl. Protein-protein hydrodynamic interactions were evaluated from concentration-dependent diffusivity measurements, obtained from dynamic light scattering. At moderate GdnHCl concentrations, lysozyme interparticle interactions are least repulsive and hydrodynamic interactions are least attractive. The lysozyme hydrodynamic radius was calculated from infinite-dilution diffusivity and did not change significantly during protein unfolding. Our results contribute toward better understanding of protein interactions of partially unfolded states in the presence of a denaturant; they may be helpful for the design of protein refolding processes that avoid protein aggregation.     

Lysozyme-lysozyme self-interactions as assessed by the osmotic second virial coefficient: Impact for physical protein stabilization

The purpose of the presented study is to understand the physicochemical properties of proteins in aqueous solutions in order to identify solution conditions with reduced attractive protein-protein interactions, to avoid the formation of protein aggregates and to increase protein solubility. This is assessed by measuring the osmotic second virial coefficient (B²²), a parameter of solution non-ideality, which is obtained using self-interaction chromatography. The model protein is lysozyme. The influence of various solution conditions on B²² was investigated: protonation degree, ionic strength, pharmaceutical relevant excipients and combinations thereof. Under acidic solution conditions B²² is positive, favoring protein repulsion. A similar trend is observed for the variation of the NaCl concentration, showing that with increasing the ionic strength protein attraction is more likely. B²² decreases and becomes negative. Thus, solution conditions are obtained favoring attractive protein-protein interactions. The B²² parameter also reflects, in general, the influence of the salts of the Hofmeister series with regard to their salting-in/salting-out effect. It is also shown that B²² correlates with protein solubility as well as physical protein stability.     

Roles of conformational stability and colloidal stability in the aggregation of recombinant human granulocyte colony-stimulating factor

We studied the non-native aggregation of recombinant human granulocyte stimulating factor (rhGCSF) in solution conditions where native rhGCSF is both conformationally stable compared to its unfolded state and at concentrations well below its solubility limit. Aggregation of rhGCSF first involves the perturbation of its native structure to form a structurally expanded transition state, followed by assembly process to form an irreversible aggregate. The energy barriers of the two steps are reflected in the experimentally measured values of free energy of unfolding (DeltaG(unf)) and osmotic second virial coefficient (B(22)), respectively. Under solution conditions where rhGCSF conformational stability dominates (i.e., large DeltaG(unf) and negative B(22)), the first step is rate-limiting, and increasing DeltaG(unf) (e.g., by the addition of sucrose) decreases aggregation. In solutions where colloidal stability is high (i.e., large and positive B(22) values) the second step is rate-limiting, and solution conditions (e.g., low pH and low ionic strength) that increase repulsive interactions between protein molecules are effective at reducing aggregation. rhGCSF aggregation is thus controlled by both conformational stability and colloidal stability, and depending on the solution conditions, either could be rate-limiting.     

Managing and mining protein crystallization data

The crystallization of macromolecules remains a major bottleneck in structural biology. The routine screening of more than one thousand crystallization conditions and subsequent optimization by fine screening presents a challenge to conventional laboratory notebook keeping. In addition, the development of high-throughput robotic crystallization and imaging systems presents a pressing need for low-cost laboratory information management system (LIMS). Here we describe CLIMS2, a crystallization LIMS that features a simple, user-friendly graphical interface, allowing the storage, management, retrieval and mining of crystallization data. The CLIMS2 executable and documentation is freely available at http://clims.med.monash.edu.au.     

Influence of precipitating agents on thermodynamic parameters of protein crystallization solutions

X‐ray crystallography is the most powerful method for determining three‐dimensional structures of proteins to (near‐)atomic resolution, but protein crystallization is a poorly explained and often intractable phenomenon. Differential Scanning Calorimetry was used to measure the thermodynamic parameters (ΔG, ΔH, ΔS) of temperature‐driven unfolding of two globular proteins, lysozyme, and ribonuclease A, in various salt solutions. The mixtures were categorized into those that were conducive to crystallization of the protein and those that were not. It was found that even fairly low salt concentrations had very large effects on thermodynamic parameters. High concentrations of salts conducive to crystallization stabilized the native folded forms of proteins, whereas high concentrations of salts that did not crystallize them tended to destabilize them. Considering the ΔH and TΔS contributions to the ΔG of unfolding separately, high concentrations of crystallizing salts were found to enthalpically stabilize and entropically destabilize the protein, and vice‐versa for the noncrystallizing salts. These observations suggest an explanation, in terms of protein stability and entropy of hydration, of why some salts are good crystallization agents for a given protein and others are not. This in turn provides theoretical insight into the process of protein crystallization, suggesting ways of predicting and controlling it. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 642–652, 2016.     

Predictive crystallization of ribonuclease A via rapid screening of osmotic second virial coefficients

Important progress has been made in recent years toward developing a molecular-level understanding of protein phase behavior in terms of the osmotic second virial coefficient, a thermodynamic parameter that characterizes pairwise protein interactions. Yet there has been little practical application of this knowledge to the field of protein crystallization, largely because of the difficult and time-consuming nature of traditional techniques for characterizing protein interactions. Self-interaction chromatography has recently been proposed as a highly efficient method for measuring the osmotic second virial coefficient. The utility of the technique is examined in this work by characterizing virial coefficients for ribonuclease A under 59 solution conditions using several crystallization additives, including PEG, sodium chloride, ammonium sulfate, and propanol. The virial coefficient measurements show some counterintuitive trends and shed light on the previous difficulties in crystallizing ribonuclease A. Crystallization experiments at the corresponding solution conditions were conducted by using ultracentrifugal crystallization. Using this methodology, ribonuclease A crystals were obtained under conditions for which the virial coefficients fell within the "crystallization slot." Crystallographic characterization showed that the crystals diffract to high resolution. Metastable crystals were also obtained for conditions outside, but near, the "crystallization slot," and they could also be frozen and used to collect structural information.     

Donnan membrane equilibrium,sedimentation equilibrium,and coil expansion of DNA in salt solutions

This is a review of applications of the McMillan-Mayer-Hill virial theory and the ionic double-layer theory to dilute colloidal solutions, in particular, solutions of DNA. Interactions of highly charged colloidal rods are developed in terms of the second virial coefficients between two rods, and between one rod and one small co-ion. The relevant cluster integrals are evaluated with interaction potentials based on the Poisson-Boltzmann equation. The treatment is extended to the intrachain repulsion responsible for the statistical swelling of coiled DNA (excluded volume effect). The theory is compared with three sets of experimental data: The salt distribution in Donnan membrane equilibria of DNA-salt solutions, sedimentation equilibria of short DNA fragments at different ionic strengths, and the intrinsic viscosity of T7 DNA in NaCl solutions. In all cases the theory agrees well with the experiments. The agreement is not convincing for the sedimentation equilibrium at low ionic strength, because here the experimental DNA concentration is too high for the truncated dilute solution expansion of the DNA-salt repulsion.     

Kinetics and equilibria of lysozyme precipitation and crystallization in concentrated ammonium sulfate solutions

The kinetics and thermodynamics of lysozyme precipitation in ammonium sulfate solutions at pH 4 and 8 and room temperature were studied. X-ray powder diffraction (XRD) was used to characterize the structure of lysozyme precipitates. It was found that, if sufficient time was allowed, microcrystals developed following an induction period after initial lysozyme precipitation, even up to ionic strengths of 8 m and at acidic pH, where lysozyme is refractory to crystallization in ammonium sulfate. The full set of precipitation and crystallization data allowed construction of a phase diagram of lysozyme, showing the ammonium sulfate dependence. It suggests that precipitation may reflect a frustrated metastable liquid-liquid phase separation, which would allow this process to be understood within the framework of the generic phase diagram for proteins. The results also demonstrate that XRD, more frequently used for characterizing inorganic and organic polycrystalline materials, is useful both in characterizing the presence of crystals in the dense phase and in verifying the crystal form of proteins.     

Protein-protein interactions in concentrated electrolyte solutions

Protein-protein interactions were measured for ovalbumin and for lysozyme in aqueous salt solutions. Protein-protein interactions are correlated with a proposed potential of mean force equal to the free energy to desolvate the protein surface that is made inaccessible to the solvent due to the protein-protein interaction. This energy is calculated from the surface free energy of the protein that is determined from protein-salt preferential-interaction parameter measurements. In classical salting-out behavior, the protein-salt preferential interaction is unfavorable. Because addition of salt raises the surface free energy of the protein according to the surface-tension increment of the salt, protein-protein attraction increases, leading to a reduction in solubility. When the surface chemistry of proteins is altered by binding of a specific ion, salting-in is observed when the interactions between (kosmotrope) ion-protein complexes are more repulsive than those between the uncomplexed proteins. However, salting-out is observed when interactions between (chaotrope) ion-protein complexes are more attractive than those of the uncomplexed proteins.     

Selection of membrane protein targets for crystallization using PFO-PAGE electrophoresis

A method to rapidly assess the oligomeric composition of multimeric proteins is notably absent from reported schemes for high throughput production and crystallization of membrane proteins. In this report we have investigated the suitability of PFO-PAGE electrophoresis for this purpose and present examples where it proves highly informative in selecting conditions favouring the functional oligomeric state of the target protein. Features such as the ability to analyze several samples in parallel, including crude membrane extracts, suggest it will be highly adaptable to high throughput analysis of membrane proteins.     

Identifying protein construct variants with increased crystallization propensity--a case study

This study describes an efficient multiparallel automated workflow of cloning, expression, purification, and crystallization of a large set of construct variants for isolated protein domains aimed at structure determination by X-ray crystallography. This methodology is applied to MAPKAP kinase 2, a key enzyme in the inflammation pathway and thus an attractive drug target. The study reveals a distinct subset of truncation variants with improved crystallization properties. These constructs distinguish themselves by increased solubility and stability during a parallel automated multistep purification process including removal of the recombinant tag. High-throughput protein melting point analysis characterizes this subset of constructs as particularly thermostable. Both parallel purification screening and melting point determination clearly identify residue 364 as the optimal C terminus for the kinase domain. Moreover, all three constructs that ultimately crystallized feature this C terminus. At the N terminus, only three amino acids differentiate a noncrystallizing from a crystallizing construct. This study addresses the very common issues associated with difficult to crystallize proteins, those of solubility and stability, and the crucial importance of particular residues in the formation of crystal contacts. A methodology is suggested that includes biophysical measurements to efficiently identify and produce construct variants of isolated protein domains which exhibit higher crystallization propensity.     

Comparison of protein-protein interactions in serine protease-inhibitor and antibody-antigen complexes: implications for the protein docking problem

The protein-protein interaction energy of 12 nonhomologous serine protease-inhibitor and 15 antibody-antigen complexes is calculated using a molecular mechanics formalism and dissected in terms of the main-chain vs. side-chain contribution, nonrotameric side-chain contributions, and amino acid residue type involvement in the interface interaction. There are major differences in the interactions of the two types of protein-protein complex. Protease-inhibitor complexes interact predominantly through a main-chain-main-chain mechanism while antibody-antigen complexes interact predominantly through a side-chain-side-chain or a side-chain-main-chain mechanism. However, there is no simple correlation between the main-chain-main-chain interaction energy and the percentage of main-chain surface area buried on binding. The interaction energy is equally effected by the presence of nonrotameric side-chain conformations, which constitute approximately 20% of the interaction energy. The ability to reproduce the interface interaction energy of the crystal structure if original side-chain conformations are removed from the calculation is much greater in the protease-inhibitor complexes than the antibody-antigen complexes. The success of a rotameric model for protein-protein docking appears dependent on the extent of the main-chain-main-chain contribution to binding. Analysis of (1) residue type and (2) residue pair interactions at the interface show that antibody-antigen interactions are very restricted with over 70% of the antibody energy attributable to just six residue types (Tyr > Asp > Asn > Ser > Glu > Trp) in agreement with previous studies on residue propensity. However, it is found here that 50% of the antigen energy is attributable to just four residue types (Arg = Lys > Asn > Asp). On average just 12 residue pair interactions (6%) contribute over 40% of the favorable interaction energy in the antibody-antigen complexes, with charge-charge and charge/polar-tyrosine interactions being prominent. In contrast protease inhibitors use a diverse set of residue types and residue pair interactions.