Protein unfolding in crowded milieu: what crowding can do to a protein undergoing unfolding?

The natural environment of a protein inside a cell is characterized by the almost complete lack of unoccupied space, limited amount of free water, and the tightly packed crowd of various biological macromolecules, such as proteins, nucleic acids, polysaccharides, and complexes thereof. This extremely crowded natural milieu is poorly mimicked by slightly salted aqueous solutions containing low concentrations of a protein of interest. The accepted practice is to model crowded environments by adding high concentrations of various polymers that serve as model “crowding agents” to the solution of a protein of interest. Although studies performed under these model conditions revealed that macromolecular crowding might have noticeable influence on various aspects related to the protein structure, function, folding, conformational stability, and aggregation propensity, the complete picture describing conformational behavior of a protein under these conditions is missing as of yet. Furthermore, there is an accepted belief that the conformational stability of globular proteins increases in the presence crowding agents due to the excluded volume effects. The goal of this study was to conduct a systematic analysis of the effect of high concentrations of PEG-8000 and Dextran-70 on the unfolding behavior of eleven globular proteins belonging to different structural classes.     

Brownian dynamics simulation of protein solutions: structural and dynamical properties

The study of solutions of biomacromolecules provides an important basis for understanding the behavior of many fundamental cellular processes, such as protein folding, self-assembly, biochemical reactions, and signal transduction. Here, we describe a Brownian dynamics simulation procedure and its validation for the study of the dynamic and structural properties of protein solutions. In the model used, the proteins are treated as atomically detailed rigid bodies moving in a continuum solvent. The protein-protein interaction forces are described by the sum of electrostatic interaction, electrostatic desolvation, nonpolar desolvation, and soft-core repulsion terms. The linearized Poisson-Boltzmann equation is solved to compute electrostatic terms. Simulations of homogeneous solutions of three different proteins with varying concentrations, pH, and ionic strength were performed. The results were compared to experimental data and theoretical values in terms of long-time self-diffusion coefficients, second virial coefficients, and structure factors. The results agree with the experimental trends and, in many cases, experimental values are reproduced quantitatively. There are no parameters specific to certain protein types in the interaction model, and hence the model should be applicable to the simulation of the behavior of mixtures of macromolecules in cell-like crowded environments.     

Translational and rotational motions of proteins in a protein crowded environment

Fluorescence correlation spectroscopy (FCS) was used to measure the translational diffusion of labeled apomyoglobin (tracer) in concentrated solutions of ribonuclease A and human serum albumin (crowders), as a quantitative model system of protein diffusive motions in crowded physiological environments. The ratio of the diffusion coefficient of the tracer protein in the protein crowded solutions and its diffusion coefficient in aqueous solution has been interpreted in terms of local apparent viscosities, a molecular parameter characteristic for each tracer-crowder system. In all protein solutions studied in this work, local translational viscosity values were larger than the solution bulk viscosity, and larger than rotational viscosities estimated for apomyoglobin in the same crowding solutions. Here we propose a method to estimate local apparent viscosities for the tracer translational and rotational diffusion directly from the bulk viscosity of the concentrated protein solutions. As a result of this study, the identification of protein species and the study of hydrodynamic changes and interactions in model crowded protein solutions by means of FCS and time-resolved fluorescence depolarization techniques may be expected to be greatly simplified.     

Function,Structure, and Stability of Enzymes Confined in Agarose Gels

Research over the past few decades has attempted to answer how proteins behave in molecularly confined or crowded environments when compared to dilute buffer solutions. This information is vital to understanding in vivo protein behavior, as the average spacing between macromolecules in the cell cytosol is much smaller than the size of the macromolecules themselves. In our study, we attempt to address this question using three structurally and functionally different model enzymes encapsulated in agarose gels of different porosities. Our studies reveal that under standard buffer conditions, the initial reaction rates of the agarose-encapsulated enzymes are lower than that of the solution phase enzymes. However, the encapsulated enzymes retain a higher percentage of their activity in the presence of denaturants. Moreover, the concentration of agarose used for encapsulation had a significant effect on the enzyme functional stability; enzymes encapsulated in higher percentages of agarose were more stable than the enzymes encapsulated in lower percentages of agarose. Similar results were observed through structural measurements of enzyme denaturation using an 8-anilinonaphthalene-1-sulfonic acid fluorescence assay. Our work demonstrates the utility of hydrogels to study protein behavior in highly confined environments similar to those present in vivo; furthermore, the enhanced stability of gel-encapsulated enzymes may find use in the delivery of therapeutic proteins, as well as the design of novel strategies for biohybrid medical devices.     

Simulation and modeling of crowding effects on the thermodynamic and kinetic properties of proteins with atomic details

Recent experimental studies of protein folding and binding under crowded solutions suggest that crowding agents exert subtle influences on the thermodynamic and kinetic properties of the proteins. While some of the crowding effects can be understood qualitatively from simple models of the proteins, quantitative rationalization of these effects requires an atomistic representation of the protein molecules in modeling their interactions with crowders. A computational approach, known as postprocessing, has opened the door for atomistic modeling of crowding effects. This review summarizes the applications of the postprocessing approach for studying crowding effects on the thermodynamics and kinetics of protein folding, conformational transition, and binding. The integration of atomistic modeling with experiments in crowded solutions promises new insight into biochemical processes in cellular environments.     

Effect of crowding on protein-protein association rates: fundamental differences between low and high mass crowding agents

Physiological media constitutes a crowded environment that serves as the field of action for protein-protein interaction in vivo. Measuring protein-protein interaction in crowded solutions can mimic this environment. In this work we follow the process of protein-protein association and its rate constants (k(on)) of the beta-lactamase (TEM)-beta-lactamase inhibitor protein (BLIP) complex in crowded solution using both low and high molecular mass crowding agents. In all crowded solutions (0-40% (w/w) of ethylene glycol (EG), poly(ethylene glycol) (PEG) 200, 1000, 3350, 8000 Da Ficoll-70 and Haemaccel the measured absolute k(on), but not k(off) values, were found to be slower as compared to buffer. However, there is a fundamental difference between low and high mass crowding agents. In the presence of low mass crowding agents and Haemaccel k(on) depends inversely on the solution viscosity. In high mass polymer solutions k(on) changes only slightly, even at viscosities 12-fold higher than water. The border between low and high molecular mass polymers is sharp and is dictated by the ratio between the polymer length (L) and its persistence length (Lp). Polymers that are long enough to form a flexible coil (L/Lp > 2) behave as high molecular mass polymers and those who are unable to do so (L/Lp < 2) behave as low molecular mass polymers. We concluded that although polymers solution are crowded, this property is not uniform; i.e. there are areas in the solution that contain bulk water, and in these areas proteins can diffuse and associate almost as if they were in diluted environment. This porous medium may be taken as mimicking some aspects of the cellular environment, where many of the macromolecules are organized along membranes and the cytoskeleton. To determine the contribution of electrostatic attraction between proteins in crowded milieu, we followed k(on) of wt-TEM and three BLIP analogs with up to 100-fold increased values of k(on) due to electrostatic steering. Faster associating BLIP variants keep their relative advantage in all crowded solutions, including Haemaccel. This result suggests that faster associating protein complexes keep their advantage also in complex environment.     

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein “breathing”, that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.     

Eutectic interactions in binary systems containing cholesterol, cholesteryl esters and triacylglycerols

Binary phase behavior of saturated cholesterol esters with trilaurin or cholesterol, and cholesterol with trilaurin, is studied. The existence of specific molecular interactions is detected by comparing the liquidus curve of the eutectic with ideal theory of freezing point depression and correcting the theoretical curve with the Bragg-Williams model when necessary. X-ray data indicate that all eutectic solids are nearly totally fractionated. The phase diagrams are sometimes well-explained by ideal solution theory indicating that polar interactions (e.g., the hydrogen bonding of cholesterol) are much less important than van der Waals interactions between neighboring molecules. However, the hydrogen bonding networks of cholesterol can lead to nonideal solution behavior with other lipids, a phenomenon consistent with previous observations on simpler molecular binaries. An observed nonideal solution behavior of triacylglycerol with cholesterol esters, on the other hand, is unexpected since significant polar interactions are expected to be 'buried' in the predominant nonpolar volume of the molecules involved.     

Effects of crowding by mono-, di-, and tetrasaccharides on cytochrome c-cytochrome c peroxidase binding: comparing experiment to theory

In cells, protein-protein interactions occur in an environment that is crowded with other molecules, but in vitro studies are almost exclusively performed in dilute solution. To gain information about the effects of crowding on protein complex formation, we used isothermal titration calorimetry to measure the stoichiometry, the free energy change, and the enthalpy change for the binding of yeast iso-1-ferricytochrome c to yeast ferricytochrome c peroxidase in dilute solution and in solutions crowded with the sugars glucose, sucrose, and stachyose. The stoichiometry is 1:1 under all conditions. The sugars stabilize the complex, but by only 0.1-0.5 kcal.mol(-)(1), and the increased stability is not correlated with the change in enthalpy of complex formation. We then compared the measured stability changes to values obtained from several analyses that are currently used to predict crowding-induced changes in biomolecular equilibria. None of the analyses are completely successful by themselves, and the results suggest that a complete analysis must account for both excluded-volume and chemical interactions.     

A cell is more than the sum of its (dilute) parts: A brief history of quinary structure

Most knowledge of protein structure and function is derived from experiments performed with purified protein resuspended in dilute, buffered solutions. However, proteins function in the crowded, complex cellular environment. Although the first four levels of protein structure provide important information, a complete understanding requires consideration of quinary structure. Quinary structure comprises the transient interactions between macromolecules that provides organization and compartmentalization inside cells. We review the history of quinary structure in the context of several metabolic pathways, and the technological advances that have yielded recent insight into protein behavior in living cells. The evidence demonstrates that protein behavior in isolated solutions deviates from behavior in the physiological environment.     

Osmotic pressures of solutions of ATP and catecholamines relating to storage in chromaffin granules

The chromaffin granule, which is the catecholamine storage organelle of the adrenal medulla, contains at least 0.73 M ions, yet it is isotonic with 0.3 osM solutions. One hypothesis which accounts for this disparity is formation of a complex between major constituents of the granule: the catecholamines, the proteins, and the ATP. In this paper we show by vapor pressure osmometry, which affords a direct measure of colligative properties, that ATP-catecholamine mixtures form highly nonideal solutions. At 37 degrees C, solutions containing 0.6 M epinephrine and 0.15 M ATP show an effective osmotic pressure of only 0.25 osM. The existence of polymeric complexes is implied by the fact that the increase of osmotic pressure with increasing concentrations of ATP and catecholamine falls off substantially at concentrations approaching those in the chromaffin granules. Neither inorganic ions nor calcium chelators cause regain of ideal colligative behavior. Osmotic measurements on model compounds suggest that the primary interaction is between the phosphate and amino groups. There is also evidence that the effects are not wholly due to the formation of discrete complexes; factors of nonideal solution behavior also play a role in lowering the osmotic pressure. These observations show that the stability of the chromaffin granule in situ can be accounted for, perhaps entirely, by spontaneous interactions among nucleotides and catecholamines.     

Nonideality and protein thermal denaturation

We studied the thermal denaturation of eglin c by using CD spectropolarimetry and differential scanning calorimetry (DSC). At low protein concentrations, denaturation is consistent with the classical two-state model. At concentrations greater than several hundred microM, however, the calorimetric enthalpy and the midpoint transition temperature increase with increasing protein concentration. These observations suggested the presence of intermediates and/or native state aggregation. However, the transitions are symmetric, suggesting that intermediates are absent, the DSC data do not fit models that include aggregation, and analytical ultracentrifugation (AUC) data show that native eglin c is monomeric. Instead, the AUC data show that eglin c solutions are nonideal. Analysis of the AUC data gives a second virial coefficient that is close to values calculated from theory and the DSC data are consistent with the behavior expected for nonideal solutions. We conclude that the concentration dependence is caused by differential nonideality of the native and denatured states. The nondeality arises from the high charge of the protein at acid pH and is exacerbated by low buffer concentrations. Our conclusion may explain differences between van't Hoff and calorimetric denaturation enthalpies observed for other proteins whose behavior is otherwise consistent with the classical two-state model.     

The effect of macromolecular crowding, ionic strength and calcium binding on calmodulin dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.     

Anomalous diffusion of proteins due to molecular crowding

We have studied the diffusion of tracer proteins in highly concentrated random-coil polymer and globular protein solutions imitating the crowded conditions encountered in cellular environments. Using fluorescence correlation spectroscopy, we measured the anomalous diffusion exponent alpha characterizing the dependence of the mean-square displacement of the tracer proteins on time, r(2)(t) approximately t(alpha). We observed that the diffusion of proteins in dextran solutions with concentrations up to 400 g/l is subdiffusive (alpha < 1) even at low obstacle concentration. The anomalous diffusion exponent alpha decreases continuously with increasing obstacle concentration and molecular weight, but does not depend on buffer ionic strength, and neither does it depend strongly on solution temperature. At very high random-coil polymer concentrations, alpha reaches a limit value of alpha(l) approximately 3/4, which we take to be the signature of a coupling between the motions of the tracer proteins and the segments of the dextran chains. A similar, although less pronounced, subdiffusive behavior is observed for the diffusion of streptavidin in concentrated globular protein solutions. These observations indicate that protein diffusion in the cell cytoplasm and nucleus should be anomalous as well, with consequences for measurements of solute diffusion coefficients in cells and for the modeling of cellular processes relying on diffusion.     

Fluorescence anisotropy as a probe to study tracer proteins in crowded solutions

Fluorescence anisotropy has been widely used to study the dynamics and interactions of biomolecules in diluted solutions. Comparable studies on single tracer macromolecules at the cellular level are now feasible because of the recent development of non-invasive fluorescence markers, like the growing family of the green fluorescence proteins (GFPs), and the advances in time-resolved fluorescence microscopy instrumentation. The interpretation of fluorescence polarization data in terms of dynamics and biological function of the macromolecular complexes in these physiological environments requires a deep understanding of the tracer rotational diffusion in such complex media. In this work we have studied the rotational diffusion of a tracer protein, apomyoglobin labeled with 1-anilino-8-naphthalene sulfonate, in crowded solutions of an unrelated protein, ribonuclease A. We have evaluated the deviation of the different tracer rotational motions from the Stokes-Einstein-Debye diffusion behavior, and its relation to the properties of the transient molecular cavities where the tracer is rotating in the fluorescence lifetime window. Finally, we have analyzed the application of fluorescence polarization methods to determine the apparent equilibrium constants of homo and hetero-associations of macromolecules in crowded conditions.     

Soft interactions and crowding

The intracellular milieu is complex, heterogeneous and crowded—an environment vastly different from dilute solutions in which most biophysical studies are performed. The crowded cytoplasm excludes about a third of the volume available to macromolecules in dilute solution. This excluded volume is the sum of two parts: steric repulsions and chemical interactions, also called soft interactions. Until recently, most efforts to understand crowding have focused on steric repulsions. Here, we summarize the results and conclusions from recent studies on macromolecular crowding, emphasizing the contribution of soft interactions to the equilibrium thermodynamics of protein stability. Despite their non-specific and weak nature, the large number of soft interactions present under many crowded conditions can sometimes overcome the stabilizing steric, excluded volume effect.     

Validation of fractal-like kinetic models by time-resolved binding kinetics of dansylamide and carbonic anhydrase in crowded media

Kinetic studies of biochemical reactions are typically carried out in a dilute solution that rarely contains anything more than reactants, products, and buffers. In such studies, mass-action-based kinetic models are used to analyze the progress curves. However, intracellular compartments are crowded by macromolecules. Therefore, we investigated the adequacy of the proposed generalizations of the mass-action model, which are meant to describe reactions in crowded media. To validate these models, we measured time-resolved kinetics for dansylamide binding to carbonic anhydrase in solutions crowded with polyethylene glycol and Ficoll. The measured progress curves clearly show the effects of crowding. The fractal-like model proposed by Savageau was used to fit these curves. In this model, the association rate coefficient k_a allometrically depends on concentrations of reactants. We also considered the fractal kinetic model proposed by Schnell and Turner, in which k_a depends on time according to a Zipf-Mandelbrot distribution, and some generalizations of these models. We found that the generalization of the mass-action model, in which association and dissociation rate coefficients are concentration-dependent, represents the preferred model. Other models based on time-dependent rate coefficients were inadequate or not preferred by model selection criteria.     

Osmotic pressure measurements of ovalbumin and lysozyme mixtures

Ovalbumin and lysozyme have been reported to undergo a mixed association in solutions of low ionic strength. Osmotic pressure experiments were performed on ovalbumin and on lysozyme solutions in 0.06 M cacodylate buffer (I = 0.02, pH = 5.8) at 30 and at 37 degrees C. The individual proteins did not undergo any self-associations at either temperature; these measurements indicated that each of the solutions was nonideal. Osmotic pressure experiments on three blends of lysozyme and ovalbumin at 30 and 37 degrees C could be interpreted in two ways. One interpretation was that a nonideal, nonassociating mixture of A and B was present; for the three solutions the mixed nonideal term BAB was negative. A negative nonideal term is usually interpreted as indicating an association. The other interpretation of the data was as a quasi-ideal mixed association of the type A + B in equilibrium AB.     

Determination of molecular parameters by fitting sedimentation data to finite-element solutions of the Lamm equation.

A method for fitting experimental sedimentation velocity data to finite-element solutions of various models based on the Lamm equation is presented. The method provides initial parameter estimates and guides the user in choosing an appropriate model for the analysis by preprocessing the data with the G(s) method by van Holde and Weischet. For a mixture of multiple solutes in a sample, the method returns the concentrations, the sedimentation (s) and diffusion coefficients (D), and thus the molecular weights (MW) for all solutes, provided the partial specific volumes (v) are known. For nonideal samples displaying concentration-dependent solution behavior, concentration dependency parameters for s(sigma) and D(delta) can be determined. The finite-element solution of the Lamm equation used for this study provides a numerical solution to the differential equation, and does not require empirically adjusted correction terms or any assumptions such as infinitely long cells. Consequently, experimental data from samples that neither clear the meniscus nor exhibit clearly defined plateau absorbances, as well as data from approach-to-equilibrium experiments, can be analyzed with this method with enhanced accuracy when compared to other available methods. The nonlinear least-squares fitting process was accomplished by the use of an adapted version of the "Doesn't Use Derivatives" nonlinear least-squares fitting routine. The effectiveness of the approach is illustrated with experimental data obtained from protein and DNA samples. Where applicable, results are compared to methods utilizing analytical solutions of approximated Lamm equations.     

Quantitative characterization of weak self-association in concentrated solutions of immunoglobulin G via the measurement of sedimentation equilibrium and osmotic pressure

The sedimentation equilibrium of solutions of immunoglobulin G in saline buffer, over a concentration range up to 125 g/L, was measured and analyzed in the context of a model that takes into account the possibility of attractive intermolecular interaction leading to the reversible formation of oligomeric species and repulsive intermolecular interaction leading to nonideal solution behavior. Additionally, previously published data on the concentration dependence of the osmotic pressure of immunoglobulin G under similar conditions, over a concentration range up to 400 g/L, were analyzed in the context of a newly developed thermodynamic formalism describing the osmotic pressure of a solution containing multiple nondiffusible solute species at an arbitrary concentration. Both sets of data are quantitatively accounted for by a model in which IgG self-associates at very high concentration to form (predominantly) trimers under the conditions of these experiments.