Factors Defining Effects of Macromolecular Crowding on Protein Stability: An in Vitro/in Silico Case Study Using Cytochrome c

Myoglobin unfolding in crowding and confinement

Crowding and confinement have often been used synonymously with regard to their effect on the structure and dynamics of proteins. In this work, we have investigated the unfolding of the protein myoglobin (Mb) entrapped in the confinement of the water pool of AOT reverse micelles and in the presence of some commonly used macromolecular crowding agents (Ficoll 70, Dextran 70, and Dextran 40). Our results reveal that confinement effects can be quite destabilizing in nature for Mb with the extent of distortion depending on a host of factors apart from the size of the confining cage. Effects of the crowding agents on myoglobin also show a deviation from the general notion that synthetic macromolecular crowding agents are always stabilizing in nature. Ficoll 70 was observed to be particularly destabilizing in its influence on Mb unfolding. Moreover, tryptophan lifetime studies point to the fact that the Trp-heme distance in Mb might not always be a reliable probe of the secondary structural dissolution of the protein.

Macromolecular crowding and protein stability

An understanding of cellular chemistry requires knowledge of how crowded environments affect proteins. The influence of crowding on protein stability arises from two phenomena, hard-core repulsions and soft (i.e., chemical) interactions. Most efforts to understand crowding effects on protein stability, however, focus on hard-core repulsions, which are inherently entropic and stabilizing. We assessed these phenomena by measuring the temperature dependence of NMR-detected amide proton exchange and used these data to extract the entropic and enthalpic contributions of crowding to the stability of ubiquitin. Contrary to expectations, the contribution of chemical interactions is large and in many cases dominates the contribution from hardcore repulsions. Our results show that both chemical interactions and hard-core repulsions must be considered when assessing the effects of crowding and help explain previous observations about protein stability and dynamics in cells.

Effects of macromolecular crowding on the structural stability of human α-lactalbumin

The folding of protein, an important process for protein to fulfill normal functions, takes place in crowded physiological environments. α-Lactalbumin, as a model system for protein-folding studies, has been used extensively because it can form stable molten globule states under a range of conditions. Here we report that the crowding agents Ficoll 70, dextran 70, and polyethylene glycol (PEG) 2000 have different effects on the structural stability of human α-lactalbumin (HLA) represented by the transition to a molten globule state: dextran 70 dramatically enhances the thermal stability of Ca(2+)-depleted HLA (apo-HLA) and Ficoll 70 enhances the thermal stability of apo-HLA to some extent, while PEG 2000 significantly decreases the thermal stability of apo-HLA. Ficoll 70 and dextran 70 have no obvious effects on trypsin degradation of apo-HLA but PEG 2000 accelerates apo-HLA degradation by trypsin and destabilizes the native conformation of apo-HLA. Furthermore, no interaction is observed between apo-HLA and Ficoll 70 or dextran 70, but a weak, non-specific interaction between the apo form of the protein and PEG 2000 is detected, and such a weak, non-specific interaction could overcome the excluded-volume effect of PEG 2000. Our data are consistent with the results of protein stability studies in cells and suggest that stabilizing excluded-volume effects of crowding agents can be ameliorated by non-specific interactions between proteins and crowders.

A physics-based approach of coarse-graining the cytoplasm of Escherichia coli (CGCYTO)

We have investigated protein stability in an environment of Escherichia coli cytoplasm using coarse-grained computer simulations. To coarse-grain a small slide of E. coli's cytoplasm consisting of over 16 million atoms, we have developed a self-assembled clustering algorithm (CGCYTO). CGCYTO uses the shape parameter and asphericity as well as a parameter λ (ranging from 0 to 1) that measures the covolume of a test protein and a macromolecule against the covolume of a test protein and a sphere of equal volume as that of a macromolecule for the criteria of coarse-graining a cytoplasmic model. A cutoff λ(c) = 0.8 was chosen based on the size of a test protein and computational resources and it determined the resolution of a coarse-grained cytoplasm. We compared the results from a polydisperse cytoplasmic model (PD model) produced by CGCYTO with two other coarse-grained hard-sphere cytoplasmic models: 1), F70 model, macromolecules in the cytoplasm were modeled by homogeneous hard spheres with a radius of 55 Å, the size of Ficoll70 and 2), HS model, each macromolecule in the cytoplasm was modeled by a hard sphere of equal volume. It was found that the folding temperature T(f) of a test protein (apoazurin) in a PD model is ~5° greater than that in a F70 model. In addition, the deviation of T(f) in a PD model is twice as much as that in a HS model when an apoazurin is randomly placed at different voids formed by particle fluctuations in PD models.

The contrasting effect of macromolecular crowding on amyloid fibril formation

Background

Amyloid fibrils associated with neurodegenerative diseases can be considered biologically relevant failures of cellular quality control mechanisms. It is known that in vivo human Tau protein, human prion protein, and human copper, zinc superoxide dismutase (SOD1) have the tendency to form fibril deposits in a variety of tissues and they are associated with different neurodegenerative diseases, while rabbit prion protein and hen egg white lysozyme do not readily form fibrils and are unlikely to cause neurodegenerative diseases. In this study, we have investigated the contrasting effect of macromolecular crowding on fibril formation of different proteins.

Methodology/Principal Findings

As revealed by assays based on thioflavin T binding and turbidity, human Tau fragments, when phosphorylated by glycogen synthase kinase-3β, do not form filaments in the absence of a crowding agent but do form fibrils in the presence of a crowding agent, and the presence of a strong crowding agent dramatically promotes amyloid fibril formation of human prion protein and its two pathogenic mutants E196K and D178N. Such an enhancing effect of macromolecular crowding on fibril formation is also observed for a pathological human SOD1 mutant A4V. On the other hand, rabbit prion protein and hen lysozyme do not form amyloid fibrils when a crowding agent at 300 g/l is used but do form fibrils in the absence of a crowding agent. Furthermore, aggregation of these two proteins is remarkably inhibited by Ficoll 70 and dextran 70 at 200 g/l.

Conclusions/Significance

We suggest that proteins associated with neurodegenerative diseases are more likely to form amyloid fibrils under crowded conditions than in dilute solutions. By contrast, some of the proteins that are not neurodegenerative disease-associated are unlikely to misfold in crowded physiological environments. A possible explanation for the contrasting effect of macromolecular crowding on these two sets of proteins (amyloidogenic proteins and non-amyloidogenic proteins) has been proposed.

Conformational sampling of peptides in the presence of protein crowders from AA/CG-multiscale simulations

Macromolecular crowding is recognized as an important factor influencing folding and conformational dynamics of proteins and nucleic acids. Previous views of crowding have focused on the mostly entropic volume exclusion effect of crowding, but recent studies are indicating the importance of enthalpic effects, in particular, changes in electrostatic interactions due to crowding. Here, temperature replica exchange molecular dynamics simulations of trp-cage and melittin in the presence of explicit protein crowders are presented to further examine the effect of protein crowders on peptide dynamics. The simulations involve a three-component multiscale modeling scheme where the peptides are represented at an atomistic level, the crowder proteins at a coarse-grained level, and the surrounding aqueous solvent as implicit solvent. This scheme optimally balances a physically realistic description for the peptide with computational efficiency. The multiscale simulations were compared with simulations of the same peptides in different dielectric environments with dielectric constants ranging from 5 to 80. It is found that the sampling in the presence of the crowders resembles sampling with reduced dielectric constants between 10 and 40. Furthermore, diverse conformational ensembles are generated in the presence of crowders including partially unfolded states for trp-cage. These findings emphasize the importance of enthalpic interactions over volume exclusion effects in describing the effects of cellular crowding.

Comparison of chemical and thermal protein denaturation by combination of computational and experimental approaches. II

Chemical and thermal denaturation methods have been widely used to investigate folding processes of proteins in vitro. However, a molecular understanding of the relationship between these two perturbation methods is lacking. Here, we combined computational and experimental approaches to investigate denaturing effects on three structurally different proteins. We derived a linear relationship between thermal denaturation at temperature T(b) and chemical denaturation at another temperature T(u) using the stability change of a protein (ΔG). For this, we related the dependence of ΔG on temperature, in the Gibbs-Helmholtz equation, to that of ΔG on urea concentration in the linear extrapolation method, assuming that there is a temperature pair from the urea (T(u)) and the aqueous (T(b)) ensembles that produces the same protein structures. We tested this relationship on apoazurin, cytochrome c, and apoflavodoxin using coarse-grained molecular simulations. We found a linear correlation between the temperature for a particular structural ensemble in the absence of urea, T(b), and the temperature of the same structural ensemble at a specific urea concentration, T(u). The in silico results agreed with in vitro far-UV circular dichroism data on apoazurin and cytochrome c. We conclude that chemical and thermal unfolding processes correlate in terms of thermodynamics and structural ensembles at most conditions; however, deviations were found at high concentrations of denaturant.

Folding dynamics of Trp-cage in the presence of chemical interference and macromolecular crowding. I

Proteins fold and function in the crowded environment of the cell's interior. In the recent years it has been well established that the so-called "macromolecular crowding" effect enhances the folding stability of proteins by destabilizing their unfolded states for selected proteins. On the other hand, chemical and thermal denaturation is often used in experiments as a tool to destabilize a protein by populating the unfolded states when probing its folding landscape and thermodynamic properties. However, little is known about the complicated effects of these synergistic perturbations acting on the kinetic properties of proteins, particularly when large structural fluctuations, such as protein folding, have been involved. In this study, we have first investigated the folding mechanism of Trp-cage dependent on urea concentration by coarse-grained molecular simulations where the impact of urea is implemented into an energy function of the side chain and/or backbone interactions derived from the all-atomistic molecular dynamics simulations with urea through a Boltzmann inversion method. In urea solution, the folding rates of a model miniprotein Trp-cage decrease and the folded state slightly swells due to a lack of contact formation between side chains at the terminal regions. In addition, the equilibrium m-values of Trp-cage from the computer simulations are in agreement with experimental measurements. We have further investigated the combined effects of urea denaturation and macromolecular crowding on Trp-cage's folding mechanism where crowding agents are modeled as hard-spheres. The enhancement of folding rates of Trp-cage is most pronounced by macromolecular crowding effect when the extended conformations of Trp-cast dominate at high urea concentration. Our study makes quantitatively testable predictions on protein folding dynamics in a complex environment involving both chemical denaturation and macromolecular crowding effects.

Role of metal in folding and stability of copper proteins in vitro

Metal coordination is required for function of many proteins. For biosynthesis of proteins coordinating a metal, the question arises if the metal binds before, during or after folding of the polypeptide. Moreover, when the metal is bound to the protein, how does its coordination affect biophysical properties such as stability and dynamics? Understanding how metals are utilized by proteins in cells on a molecular level requires accurate descriptions of the thermodynamic and kinetic parameters involved in protein-metal complexes. Copper is one of the essential transition metals found in the active sites of many key proteins. To avoid toxicity of free copper ions, living systems have developed elaborate copper-transport systems that involve dedicated proteins that facilitate efficient and specific delivery of copper to target proteins. This review describes in vitro and in silico biophysical work assessing the role of copper in folding and stability of copper-binding proteins. Examples of proteins discussed are: a blue-copper protein (Pseudomonas aeruginosa azurin), members of copper-transport systems (bacterial CopZ, human Atox1 and ATP7B domains) and multi-copper ferroxidases (yeast Fet3p and human ceruloplasmin). The consequences of interactions between copper proteins and platinum-complexes are also discussed. This article is part of a Special Issue entitled: Cell Biology of Metals.

Three-dimensional stochastic off-lattice model of binding chemistry in crowded environments

Molecular crowding is one of the characteristic features of the intracellular environment, defined by a dense mixture of varying kinds of proteins and other molecules. Interaction with these molecules significantly alters the rates and equilibria of chemical reactions in the crowded environment. Numerous fundamental activities of a living cell are strongly influenced by the crowding effect, such as protein folding, protein assembly and disassembly, enzyme activity, and signal transduction. Quantitatively predicting how crowding will affect any particular process is, however, a very challenging problem because many physical and chemical parameters act synergistically in ways that defy easy analysis. To build a more realistic model for this problem, we extend a prior stochastic off-lattice model from two-dimensional (2D) to three-dimensional (3D) space and examine how the 3D results compare to those found in 2D. We show that both models exhibit qualitatively similar crowding effects and similar parameter dependence, particularly with respect to a set of parameters previously shown to act linearly on total reaction equilibrium. There are quantitative differences between 2D and 3D models, although with a generally gradual nonlinear interpolation as a system is extended from 2D to 3D. However, the additional freedom of movement allowed to particles as thickness of the simulation box increases can produce significant quantitative change as a system moves from 2D to 3D. Simulation results over broader parameter ranges further show that the impact of molecular crowding is highly dependent on the specific reaction system examined.

Crowding effects on the small, fast-folding protein lambda6-85

The microsecond folder lambda6-85 is a small (9.2 kDa = 9200 amu) five helix bundle protein. We investigated the stability of lambda6-85 in two different low-fluorescence crowding matrices: the large 70 kDa carbohydrate Ficoll 70, and the small 14 kDa thermophilic protein SubL. The same thermal stability of secondary structure was measured by circular dichroism in aqueous buffer, and at a crowding fraction phi = 15 +/- 1% of Ficoll 70. Tryptophan fluorescence detection (probing a tertiary contact) yielded the same thermal stability in Ficoll, but 4 degrees C lower in aqueous buffer. Temperature-jump kinetics revealed that the relaxation rate, corrected for bulk viscosity, was very similar in Ficoll and in aqueous buffer. Thus viscosity, hydrodynamics and crowding seem to compensate one another. However, a new fast phase was observed in Ficoll, attributed to crowding-induced downhill folding. We also measured the stability of lambda6-85 in phi = 14 +/- 1% SubL, which acts as a smaller more rigid crowder. Significantly greater stabilization (7 to 13 degrees C depending on probe) was observed than in the Ficoll matrix. The results highlight the importance of crowding agent choice for studies of small, fast-folding proteins amenable to comparison with molecular dynamics simulations.

Protein stability and folding kinetics in the nucleus and endoplasmic reticulum of eucaryotic cells

We measure the stability and folding relaxation rate of phosphoglycerate kinase (PGK) Förster resonance energy transfer (FRET) constructs localized in the nucleus or in the endoplasmic reticulum (ER) of eukaryotic cells. PGK has a more compact native state in the cellular compartments than in aqueous solution. Its native FRET signature is similar to that previously observed in a carbohydrate-crowding matrix, consistent with crowding being responsible for the compact native state of PGK in the cell. PGK folds through multiple states in vitro, but its folding kinetics is more two-state-like in the ER, so the folding mechanism can be modified by intracellular compartments. The nucleus increases PGK stability and folding rate over the cytoplasm and ER, even though the density of crowders in the nucleus is no greater than in the ER or cytoplasm. Nuclear folding kinetics (and to a lesser extent, thermodynamics) vary less from cell to cell than in the cytoplasm or ER, indicating a more homogeneous crowding and chemical environment in the nucleus.

Unified regression model of binding equilibria in crowded environments

Molecular crowding is a critical feature distinguishing intracellular environments from idealized solution-based environments and is essential to understanding numerous biochemical reactions, from protein folding to signal transduction. Many biochemical reactions are dramatically altered by crowding, yet it is extremely difficult to predict how crowding will quantitatively affect any particular reaction systems. We previously developed a novel stochastic off-lattice model to efficiently simulate binding reactions across wide parameter ranges in various crowded conditions. We now show that a polynomial regression model can incorporate several interrelated parameters influencing chemistry under crowded conditions. The unified model of binding equilibria accurately reproduces the results of particle simulations over a broad range of variation of six physical parameters that collectively yield a complicated, non-linear crowding effect. The work represents an important step toward the long-term goal of computationally tractable predictive models of reaction chemistry in the cellular environment.

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.

Proteins are workhorses for driving biological functions inside cells. Calmodulin (CaM) is a protein that can carry cellular signals by triggered conformational changes due to calcium binding that alters target binding. Interestingly, CaM is able to bind over 300 targets. One of the challenges in characterizing CaM's ability to bind multiple targets lies in that CaM is a flexible protein and its structure is easily modulated by the physicochemical changes in its surroundings, particularly inside a complex cellular milieu. In order to determine structure-function relationships of CaM, we employed a combined approach of experiments, computer simulations and statistical physics in the investigation of the effect of calcium-binding, salt concentration, and macromolecular crowding on CaM. The results revealed unique folding energy landscapes of CaM in the absence and presence of calcium ions and the structural implications of CaM are interpreted under cell-like conditions. Further, a large conformational change in CaM in response to environmental impacts, dictates the packing of local helices that may be critical to its function of target binding and recognition among vast target selections.

Modulation of functionally significant conformational equilibria in adenylate kinase by high concentrations of trimethylamine oxide attributed to volume exclusion

The effect of an inert small molecule osmolyte, trimethyl amine N-oxide (TMAO), upon the conformational equilibria of Escherichia coli adenylate kinase was studied using time-resolved FRET. The relative populations of open and closed clefts between the LID and the CORE domains were measured as functions of the concentrations of the substrate ATP over the concentration range 0-18 mM and TMAO over the concentration range 0-4 M. A model was constructed according to which the enzyme exists in equilibrium among four conformational states, corresponding to combinations of open and closed conformations of the LID-CORE and AMP-CORE clefts. ATP is assumed to bind only to those conformations with the closed LID-CORE cleft, and TMAO is assumed to be differentially excluded as a hard spherical particle from each of the four conformations in accordance with calculations based upon x-ray crystallographic structures. This model was found to describe quantitatively the dependence of the fraction of the closed LID-CORE cleft upon the concentrations of both ATP and TMAO over the entire range of concentrations with just five undetermined parameters.

The parallel G-quadruplex structure of vertebrate telomeric repeat sequences is not the preferred folding topology under physiological conditions

G-quadruplex topologies of telomeric repeat sequences from vertebrates were investigated in the presence of molecular crowding (MC) mimetics, namely polyethylene glycol 200 (PEG), Ficoll 70 as well as Xenopus laevis egg extract by CD and NMR spectroscopy and native PAGE. Here, we show that the conformational behavior of the telomeric repeats in X. laevis egg extract or in Ficoll is notably different from that observed in the presence of PEG. While the behavior of the telomeric repeat in X. laevis egg extract or in Ficoll resembles results obtained under dilute conditions, PEG promotes the formation of high-order parallel topologies. Our data suggest that PEG should not be used as a MC mimetic.

Power-law dependence of the melting temperature of ubiquitin on the volume fraction of macromolecular crowders

The dependence of the melting temperature increase (ΔT(m)) of the protein ubiquitin on the volume fraction (ϕ) of several commonly used macromolecular crowding agents (dextran 6, 40, and 70 and ficoll 70) was quantitatively examined and compared to a recently developed theoretical crowding model, i.e., ΔT(m) ∼ (R(g)∕R(c))(α)φ(α∕3). We found that in the current case this model correctly predicts the power-law dependence of ΔT(m) on φ but significantly overestimates the role of the size (i.e., R(c)) of the crowding agent. In addition, we found that for ubiquitin the exponent α is in the range of 4.1-6.5, suggesting that the relation of α=3∕(3ν-1) is a better choice for estimating α based on the Flory coefficient (ν) of the polypeptide chain. Taken together these findings highlight the importance of improving our knowledge and theoretical treatment of the microcompartmentalization of the commonly used model crowding agents.

Protein folding in the cell: challenges and progress

It is hard to imagine a more extreme contrast than that between the dilute solutions used for in vitro studies of protein folding and the crowded, compartmentalized, sticky, spatially inhomogeneous interior of a cell. This review highlights recent research exploring protein folding in the cell with a focus on issues that are generally not relevant to in vitro studies of protein folding, such as macromolecular crowding, hindered diffusion, cotranslational folding, molecular chaperones, and evolutionary pressures. The technical obstacles that must be overcome to characterize protein folding in the cell are driving methodological advances, and we draw attention to several examples, such as fluorescence imaging of folding in cells and genetic screens for in-cell stability.