Atomistic Modeling of Macromolecular Crowding Predicts Modest Increases in Protein Folding and Binding Stability

Theoretical models predict that macromolecular crowding can increase protein folding stability, but depending on details of the models (e.g., how the denatured state is represented), the level of stabilization predicted can be very different. In this study, we represented the native and denatured states atomistically, with conformations sampled from explicit-solvent molecular dynamics simulations at room temperature and high temperature, respectively. We then designed an efficient algorithm to calculate the allowed fraction, f, when the protein molecule is placed inside a box of crowders. That a fraction of placements of the protein molecule is disallowed because of volume exclusion by the crowders leads to an increase in chemical potential, given by Δμ = −k_BT lnf. The difference in Δμ between the native and denatured states predicts the effect of crowding on the folding free energy. Even when the crowders occupied 35% of the solution volume, the stabilization reached only 1.5 kcal/mol for cytochrome b₅₆₂. The modest stabilization predicted is consistent with experimental studies. Interestingly, a mixture of different sized crowders was found to exert a greater effect than the sum of the individual species of crowders. The stabilization of crowding on the binding stability of barnase and barstar, based on atomistic modeling of the proteins, was similarly modest. These findings have profound implications for macromolecular crowding inside cells.     

The Folding Mechanism of BBL: Plasticity of Transition-State Structure Observed within an Ultrafast Folding Protein Family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Φ-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Φ-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Φ values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Φ values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.     

Macromolecular Crowding as a Suppressor of Human IAPP Fibril Formation and Cytotoxicity

The biological cell is known to exhibit a highly crowded milieu, which significantly influences protein aggregation and association processes. As several cell degenerative diseases are related to the self-association and fibrillation of amyloidogenic peptides, understanding of the impact of macromolecular crowding on these processes is of high biomedical importance. It is further of particular relevance as most in vitro studies on amyloid aggregation have been performed in diluted solution which does not reflect the complexity of their cellular surrounding. The study presented here focuses on the self-association of the type-2 diabetes mellitus related human islet amyloid polypeptide (hIAPP) in various crowded environments including network-forming macromolecular crowding reagents and protein crowders. It was possible to identify two competing processes: a crowder concentration and type dependent stabilization of globular off-pathway species and a – consequently - retarded or even inhibited hIAPP fibrillation reaction. The cause of these crowding effects was revealed to be mainly excluded volume in the polymeric crowders, whereas non-specific interactions seem to be most dominant in protein crowded environments. Specific hIAPP cytotoxicity assays on pancreatic β-cells reveal non-toxicity for the stabilized globular species, in contrast to the high cytotoxicity imposed by the normal fibrillation pathway. From these findings it can be concluded that cellular crowding is able to effectively stabilize the monomeric conformation of hIAPP, hence enabling the conduction of its normal physiological function and prevent this highly amyloidogenic peptide from cytotoxic aggregation and fibrillation.     

Protein Folding: Adding a Nucleus to Guide Helix Docking Reduces Landscape Roughness

The elongated three-helix‐bundle spectrin domains R16 and R17 fold and unfold unusually slowly over a rough energy landscape, in contrast to the homologue R15, which folds fast over a much smoother, more typical landscape. R15 folds via a nucleation–condensation mechanism that guides the docking of the A and C-helices. However, in R16 and R17, the secondary structure forms first and the two helices must then dock in the correct register. Here, we use variants of R16 and R17 to demonstrate that substitution of just five key residues is sufficient to alter the folding mechanism and reduce the landscape roughness. We suggest that, by providing access to an alternative, faster, folding route over their landscape, R16 and R17 can circumvent their slow, frustrated wild-type folding mechanism.     

Effect of Macromolecular Crowding on Protein Binding Stability: Modest Stabilization and Significant Biological Consequences

Macromolecular crowding has long been known to significantly affect protein oligomerization, and yet no direct quantitative measurements appear to have been made of its effects on the binding free energy of the elemental step of adding a single subunit. Here, we report the effects of two crowding agents on the binding free energy of two subunits in the Escherichia coli polymerase III holoenzyme. The crowding agents are found, paradoxically, to have only a modest stabilizing effect, of the order of 1 kcal/mol, on the binding of the two subunits. Systematic variations in the level of stabilization with crowder size are nevertheless observed. The data are consistent with theoretical predictions based on atomistic modeling of excluded-volume interactions with crowders. We reconcile the apparent paradox presented by our data by noting that the modest effects of crowding on elemental binding steps are cumulative, and thus lead to substantial stabilization of higher oligomers. Correspondingly, the effects of small variations in the level of crowding during the lifetime of a cell may be magnified, suggesting that crowding may play a role in increased susceptibility to protein aggregation-related diseases with aging.     

Denatured State Structural Property Determines Protein Stabilization by Macromolecular Crowding: A Thermodynamic and Structural Approach

Understanding of protein structure and stability gained to date has been acquired through investigations made under dilute conditions where total macromolecular concentration never surpasses 10 g l⁻¹. However, biological macromolecules are known to evolve and function under crowded intracellular environments that comprises of proteins, nucleic acids, ribosomes and carbohydrates etc. Crowded environment is known to result in altered biological properties including thermodynamic, structural and functional aspect of macromolecules as compared to the macromolecules present in our commonly used experimental dilute buffers (for example, Tris HCl or phosphate buffer). In this study, we have investigated the thermodynamic and structural consequences of synthetic crowding agent (Ficoll 70) on three different proteins (Ribonuclease-A, lysozyme and holo α-lactalbumin) at different pH values. We report here that the effect of crowding is protein dependent in terms of protein thermal stability and structure. We also observed that the structural characteristics of the denatured state determines if crowding will have an effect or not on the protein stability.     

Effects of Macromolecular Crowding on the Structure of a Protein Complex: A Small-Angle Scattering Study of Superoxide Dismutase

Macromolecular crowding can alter the structure and function of biological macromolecules. We used small-angle scattering to measure the effects of macromolecular crowding on the size of a protein complex, SOD (superoxide dismutase). Crowding was induced using 400 MW PEG (polyethylene glycol),TEG (triethylene glycol), α-MG (methyl-α-glucoside), and TMAO (trimethylamine n-oxide). Parallel small-angle neutron scattering and small-angle x-ray scattering allowed us to unambiguously attribute apparent changes in radius of gyration to changes in the structure of SOD. For a 40% PEG solution, we find that the volume of SOD was reduced by 9%. Considering the osmotic pressure due to PEG, this deformation corresponds to a highly compressible structure. Small-angle x-ray scattering done in the presence of TEG suggests that for further deformation—beyond a 9% decrease in volume—the resistance to deformation may increase dramatically.     

Minimal Effects of Macromolecular Crowding on an Intrinsically Disordered Protein: A Small-Angle Neutron Scattering Study

Small-angle neutron scattering was used to study the effects of macromolecular crowding by two globular proteins, i.e., bovine pancreatic trypsin inhibitor and equine metmyoglobin, on the conformational ensemble of an intrinsically disordered protein, the N protein of bacteriophage λ. The λ N protein was uniformly labeled with ²H, and the concentrations of D₂O in the samples were adjusted to match the neutron scattering contrast of the unlabeled crowding proteins, thereby masking their contribution to the scattering profiles. Scattering from the deuterated λ N was recorded for samples containing up to 0.12 g/mL bovine pancreatic trypsin inhibitor or 0.2 g/mL metmyoglobin. The radius of gyration of the uncrowded protein was estimated to be 30 Å and was found to be remarkably insensitive to the presence of crowders, varying by <2 Å for the highest crowder concentrations. The scattering profiles were also used to estimate the fractal dimension of λ N, which was found to be ∼1.8 in the absence or presence of crowders, indicative of a well-solvated and expanded random coil under all of the conditions examined. These results are contrary to the predictions of theoretical treatments and previous experimental studies demonstrating compaction of unfolded proteins by crowding with polymers such as dextran and Ficoll. A computational simulation suggests that some previous treatments may have overestimated the effective volumes of disordered proteins and the variation of these volumes within an ensemble. The apparent insensitivity of λ N to crowding may also be due in part to weak attractive interactions with the crowding proteins, which may compensate for the effects of steric exclusion.     

Minimal Effects of Macromolecular Crowding on an Intrinsically Disordered Protein: A Small-Angle Neutron Scattering Study

Small-angle neutron scattering was used to study the effects of macromolecular crowding by two globular proteins, i.e., bovine pancreatic trypsin inhibitor and equine metmyoglobin, on the conformational ensemble of an intrinsically disordered protein, the N protein of bacteriophage λ. The λ N protein was uniformly labeled with ²H, and the concentrations of D₂O in the samples were adjusted to match the neutron scattering contrast of the unlabeled crowding proteins, thereby masking their contribution to the scattering profiles. Scattering from the deuterated λ N was recorded for samples containing up to 0.12 g/mL bovine pancreatic trypsin inhibitor or 0.2 g/mL metmyoglobin. The radius of gyration of the uncrowded protein was estimated to be 30 Å and was found to be remarkably insensitive to the presence of crowders, varying by <2 Å for the highest crowder concentrations. The scattering profiles were also used to estimate the fractal dimension of λ N, which was found to be ∼1.8 in the absence or presence of crowders, indicative of a well-solvated and expanded random coil under all of the conditions examined. These results are contrary to the predictions of theoretical treatments and previous experimental studies demonstrating compaction of unfolded proteins by crowding with polymers such as dextran and Ficoll. A computational simulation suggests that some previous treatments may have overestimated the effective volumes of disordered proteins and the variation of these volumes within an ensemble. The apparent insensitivity of λ N to crowding may also be due in part to weak attractive interactions with the crowding proteins, which may compensate for the effects of steric exclusion.     

Connecting the Dots: The Effects of Macromolecular Crowding on Cell Physiology

The physicochemical properties of cellular environments with a high macromolecular content have been systematically characterized to explain differences observed in the diffusion coefficients, kinetics parameters, and thermodynamic properties of proteins inside and outside of cells. However, much less attention has been given to the effects of macromolecular crowding on cell physiology. Here, we review recent findings that shed some light on the role of crowding in various cellular processes, such as reduction of biochemical activities, structural reorganization of the cytoplasm, cytoplasm fluidity, and cellular dormancy. We conclude by presenting some unresolved problems that require the attention of biophysicists, biochemists, and cell physiologists. Although it is still underappreciated, macromolecular crowding plays a critical role in life as we know it.     

The Effect of Attractive Interactions and Macromolecular Crowding on Crystallins Association

In living systems proteins are typically found in crowded environments where their effective interactions strongly depend on the surrounding medium. Yet, their association and dissociation needs to be robustly controlled in order to enable biological function. Uncontrolled protein aggregation often causes disease. For instance, cataract is caused by the clustering of lens proteins, i.e., crystallins, resulting in enhanced light scattering and impaired vision or blindness. To investigate the molecular origins of cataract formation and to design efficient treatments, a better understanding of crystallin association in macromolecular crowded environment is needed. Here we present a theoretical study of simple coarse grained colloidal models to characterize the general features of how the association equilibrium of proteins depends on the magnitude of intermolecular attraction. By comparing the analytic results to the available experimental data on the osmotic pressure in crystallin solutions, we identify the effective parameters regimes applicable to crystallins. Moreover, the combination of two models allows us to predict that the number of binding sites on crystallin is small, i.e. one to three per protein, which is different from previous estimates. We further observe that the crowding factor is sensitive to the size asymmetry between the reactants and crowding agents, the shape of the protein clusters, and to small variations of intermolecular attraction. Our work may provide general guidelines on how to steer the protein interactions in order to control their association.     

Macromolecular Crowding Effects on Two Homologs of Ribosomal Protein S16: Protein-Dependent Structural Changes and Local Interactions

Proteins function in cellular environments that are crowded with biomolecules, and in this reduced available space, their biophysical properties may differ from those observed in dilute solutions in vitro. Here, we investigated the effects of a synthetic macromolecular crowding agent, dextran 20, on the folded states of hyperthermophilic (S16_Thermo) and mesophilic (S16_Meso) homologs of the ribosomal protein S16. As expected for an excluded-volume effect, the resistance of the mesophilic protein to heat-induced unfolding increased in the presence of dextran 20, and chemical denaturation experiments at different fixed temperatures showed the macromolecular crowding effect to be temperature-independent. Förster resonance energy transfer experiments show that intramolecular distances between an intrinsic Trp residue and BODIPY-labeled S16_Meso depend on the level of the crowding agent. The BODIPY group was attached at three specific positions in S16_Meso, allowing measurements of three intraprotein distances. All S16_Meso variants exhibited a decrease in the average Trp-BODIPY distance at up to 100 mg/mL dextran 20, whereas the changes in distance became anisotropic (one distance increased, two distances decreased) at higher dextran concentrations. In contrast, the two S16_Thermo mutants did not show any changes in Trp-BODIPY distances upon increase of dextran 20 concentrations. It should be noted that the fluorescence quantum yields and lifetimes of BODIPY attached to the two S16 homologs decreased gradually in the presence of dextran 20. To investigate the origin of this decrease, we studied the BODIPY quantum yield in three protein variants in the presence of a tyrosine-labeled dextran. The experiments revealed distinct tyrosine quenching behaviors of BODIPY in the three variants, suggesting a dynamic local interaction between dextran and one particular S16 variant.     

Macromolecular Crowding Effects on Two Homologs of Ribosomal Protein S16: Protein-Dependent Structural Changes and Local Interactions

Proteins function in cellular environments that are crowded with biomolecules, and in this reduced available space, their biophysical properties may differ from those observed in dilute solutions in vitro. Here, we investigated the effects of a synthetic macromolecular crowding agent, dextran 20, on the folded states of hyperthermophilic (S16_Thermo) and mesophilic (S16_Meso) homologs of the ribosomal protein S16. As expected for an excluded-volume effect, the resistance of the mesophilic protein to heat-induced unfolding increased in the presence of dextran 20, and chemical denaturation experiments at different fixed temperatures showed the macromolecular crowding effect to be temperature-independent. Förster resonance energy transfer experiments show that intramolecular distances between an intrinsic Trp residue and BODIPY-labeled S16_Meso depend on the level of the crowding agent. The BODIPY group was attached at three specific positions in S16_Meso, allowing measurements of three intraprotein distances. All S16_Meso variants exhibited a decrease in the average Trp-BODIPY distance at up to 100 mg/mL dextran 20, whereas the changes in distance became anisotropic (one distance increased, two distances decreased) at higher dextran concentrations. In contrast, the two S16_Thermo mutants did not show any changes in Trp-BODIPY distances upon increase of dextran 20 concentrations. It should be noted that the fluorescence quantum yields and lifetimes of BODIPY attached to the two S16 homologs decreased gradually in the presence of dextran 20. To investigate the origin of this decrease, we studied the BODIPY quantum yield in three protein variants in the presence of a tyrosine-labeled dextran. The experiments revealed distinct tyrosine quenching behaviors of BODIPY in the three variants, suggesting a dynamic local interaction between dextran and one particular S16 variant.     

Dynameomics: A Consensus View of the Protein Unfolding/Folding Transition State Ensemble across a Diverse Set of Protein Folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-α, all-β, and mixed α/β and α + β) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an α-helix maintained more structure in the TS ensemble than did residues starting in β-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.     

Localization of Protein Aggregation in Escherichia coli Is Governed by Diffusion and Nucleoid Macromolecular Crowding Effect

Aggregates of misfolded proteins are a hallmark of many age-related diseases. Recently, they have been linked to aging of Escherichia coli (E. coli) where protein aggregates accumulate at the old pole region of the aging bacterium. Because of the potential of E. coli as a model organism, elucidating aging and protein aggregation in this bacterium may pave the way to significant advances in our global understanding of aging. A first obstacle along this path is to decipher the mechanisms by which protein aggregates are targeted to specific intercellular locations. Here, using an integrated approach based on individual-based modeling, time-lapse fluorescence microscopy and automated image analysis, we show that the movement of aging-related protein aggregates in E. coli is purely diffusive (Brownian). Using single-particle tracking of protein aggregates in live E. coli cells, we estimated the average size and diffusion constant of the aggregates. Our results provide evidence that the aggregates passively diffuse within the cell, with diffusion constants that depend on their size in agreement with the Stokes-Einstein law. However, the aggregate displacements along the cell long axis are confined to a region that roughly corresponds to the nucleoid-free space in the cell pole, thus confirming the importance of increased macromolecular crowding in the nucleoids. We thus used 3D individual-based modeling to show that these three ingredients (diffusion, aggregation and diffusion hindrance in the nucleoids) are sufficient and necessary to reproduce the available experimental data on aggregate localization in the cells. Taken together, our results strongly support the hypothesis that the localization of aging-related protein aggregates in the poles of E. coli results from the coupling of passive diffusion-aggregation with spatially non-homogeneous macromolecular crowding. They further support the importance of “soft” intracellular structuring (based on macromolecular crowding) in diffusion-based protein localization in E. coli.     

Dependence of Protein Folding Stability and Dynamics on the Density and Composition of Macromolecular Crowders

We investigate the effect of macromolecular crowding on protein folding, using purely repulsive crowding particles and a self-organizing polymer model of protein folding. We find that the variation in folding stability with crowder size for typical α-, β-, and α/β-proteins is well described by an adaptation of the scaled particle theory. The native state, the transition state, and the unfolded protein are treated as effective hard spheres, with the folded and transition state radii independent of the size and concentration of the crowders. Remarkably, we find that, as the effective unfolded state radius is very weakly dependent on the crowder concentration, it can also be approximated by a single size. The same model predicts the effect of crowding on the folding barrier and therefore refolding rates with no adjustable parameters. A simple extension of the scaled-particle theory model, assuming additivity, can also describe the behavior of mixtures of crowding particles.