Effect of a concentrated "inert" macromolecular cosolute on the stability of a globular protein with respect to denaturation by heat and by chaotropes: a statistical-thermodynamic model

An equilibrium statistical-thermodynamic model for the effect of volume exclusion arising from high concentrations of stable macromolecules upon the stability of a trace globular protein with respect to denaturation by heat and by chaotropes is presented. The stable cosolute and the native form of the trace protein are modeled by effective hard spherical particles. The denatured state of the trace protein is represented as an ensemble of substates modeled by random coils having the same contour length but different rms end-to-end distances (i.e., different degrees of compaction). The excess or nonideal chemical potential of the native state and of each denatured substate is calculated as a function of the concentration of stable cosolute, leading to an estimate of the relative abundance of each state and substate, and the ensemble average free energy of the transition between native and denatured protein. The effect of the addition of stable cosolute upon the temperature of half-denaturation and upon the concentration of chaotrope required to half-denature the tracer at constant temperature is then estimated. At high cosolute concentration (>100 g/l) these effects are predicted to be large and readily measurable experimentally, provided that an experimental system exhibiting a fully reversible unfolding equilibrium at high total macromolecular concentration can be developed.     

Alanine Scan of Core Positions in Ubiquitin Reveals Links between Dynamics,Stability, and Function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.     

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg²⁺ ion concentrations are low, K⁺ concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo–like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg²⁺ (0.5–2 mM) and K⁺ (140 mM) if the solution is supplemented with physiological amounts (∼20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.     

Design and structural analysis of an engineered thermostable chicken lysozyme.

A hyperstable (hs) variant of chicken egg-white lysozyme with enhanced thermal (delta Tm approximately +10.5 degrees C) and chemical (delta Cm for guanidine hydrochloride denaturation = +1.3 M) stabilities relative to wild-type (WT) was constructed by combining several individual stabilizing substitutions. The free energy difference between the native and denatured states of the hs variant is 3.1 (GdnHCl, 25 degrees C) to 4.0 (differential scanning calorimetry, 74 degrees C) kcal mol-1 greater than that of WT. The specific activity of the hs variant is 2.5-fold greater than that of WT. The choice of mutations came from diverse sources: (1) The I55L/S91T core construct with delta Tm = 3.3 degrees C from WT was available from the accompanying study (Shih P, Holland DR, Kirsch JF, 1995, Protein Sci 4:2050-2062). (2) The A31V mutation was suggested by the better atomic packing in the human lysozyme structure where the Ala 31 equivalent is Leu. (3) The H15L and R114H substitutions were selected on the basis of sequence comparisons with pheasant lysozymes that are more stable than the chicken enzyme. (4) The D101S variant was identified from a screen of mutants previously prepared in this laboratory. The effects of the individual mutations on stability are cumulative and nearly additive.     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

The interplay between protein stability and dynamics in conformational diseases: The case of hPGK1 deficiency

Conformational diseases often show defective protein folding efficiency in vivo upon mutation, affecting protein properties such as thermodynamic stability and folding/unfolding/misfolding kinetics as well as the interactions of the protein with the protein homeostasis network. Human phosphoglycerate kinase 1 (hPGK1) deficiency is a rare inherited disease caused by mutations in hPGK1 that lead to loss-of-function. This disease offers an excellent opportunity to explore the complex relationships between protein stability and dynamics because of the different unfolding mechanisms displayed towards chemical and thermal denaturation. This work explores these relationships using two thermostable mutants (p.E252A and p.T378P) causing hPGK1 deficiency and WT hPGK1 using proteolysis and chemical denaturation. p.T378P is degraded ~ 30-fold faster at low protease concentrations (here, the proteolysis step is rate-limiting) and ~ 3-fold faster at high protease concentrations (where unfolding kinetics is rate-limiting) than WT and p.E252A, indicating that p.T378P is thermodynamically and kinetically destabilized. Urea denaturation studies support the decrease in thermodynamic stability and folding cooperativity for p.T378P, as well as changes in folding/unfolding kinetics. The present study reveals changes in the folding landscape of hPGK1 upon mutation that may affect protein folding efficiency and stability in vivo, also suggesting that native state stabilizers and protein homeostasis modulators may help to correct folding defects in hPGK1 deficiency. Moreover, detailed kinetic proteolysis studies are shown to be powerful and simple tools to provide deep insight into mutational effects on protein folding and stability in conformational diseases.     

Contribution of Thr29 to the thermodynamic stability of goat alpha-lactalbumin as determined by experimental and theoretical approaches

The Thr29 residue in the hydrophobic core of goat alpha-lactalbumin (alpha-LA) was substituted with Val (Thr29Val) and Ile (Thr29Ile) to investigate the contribution of Thr29 to the thermodynamic stability of the protein. We carried out protein stability measurements, X-ray crystallographic analyses, and free energy calculations based on molecular dynamics simulation. The equilibrium unfolding transitions induced by guanidine hydrochloride demonstrated that the Thr29Val and Thr29Ile mutants were, respectively, 1.9 and 3.2 kcal/mol more stable than the wild-type protein (WT). The overall structures of the mutants were almost identical to that of WT, in spite of the disruption of the hydrogen bonding between the side-chain O-H group of Thr29 and the main-chain C=O group of Glu25. To analyze the stabilization mechanism of the mutants, we performed free energy calculations. The calculated free energy differences were in good agreement with the experimental values. The stabilization of the mutants was mainly caused by solvation loss in the denatured state. Furthermore, the O-H group of Thr29 favorably interacts with the C=O group of Glu25 to form hydrogen bonds and, simultaneously, unfavorably interacts electrostatically with the main-chain C=O group of Thr29. The difference in the free energy profile of the unfolding path between WT and the Thr29Ile mutant is discussed in light of our experimental and theoretical results.     

Effect of amino acid replacements on the structural stability of fish myoglobin

Structural stabilities of myoglobin (Mb) from several fish (scombridae) species differ significantly, although their amino acid sequence identity is very high (>95%), suggesting that only a few substitutions greatly affect the stability of Mb. Accordingly, recombinant Mbs with point mutation(s) derived from bigeye tuna Mb cDNA were expressed as GST-fusion proteins in the soluble fractions of Escherichia coli. After removal of the GST segment, the stability of five mutants, namely, P13A, I21M, V57I, A62G, and I21M/V57I, together with the wild type (WT) were investigated, taking temperature dependency of alpha-helical content and denaturant concentration dependency of Soret band absorbance as parameters. As a result, the stability of P13A against denaturants and its alpha-helical content at 10 degrees C was found to be the highest among the mutants, whereas those of A62G were the lowest. The stabilities of V57I and I21M/V57I were higher than that of WT, though that of I21M was nearly the same as WT. These findings suggest that the structural stability of fish Mb is tuned up only by the substitutions of a few amino acid residues located in the alpha-helical segments forming the hydrophobic heme pocket.     

Tumor necrosis factor receptor-associated factor 6 (TRAF6) associates with huntingtin protein and promotes its atypical ubiquitination to enhance aggregate formation

Huntington disease (HD) is a neurodegenerative disorder caused by an expansion of polyglutamines in the first exon of huntingtin (HTT), which confers aggregation-promoting properties to amino-terminal fragments of the protein (N-HTT). Mutant N-HTT aggregates are enriched for ubiquitin and contain ubiquitin E3 ligases, thus suggesting a role for ubiquitination in aggregate formation. Here, we report that tumor necrosis factor receptor-associated factor 6 (TRAF6) binds to WT and polyQ-expanded N-HTT in vitro as well as to endogenous full-length proteins in mouse and human brain in vivo. Endogenous TRAF6 is recruited to cellular inclusions formed by mutant N-HTT. Transient overexpression of TRAF6 promotes WT and mutant N-HTT atypical ubiquitination with Lys(6), Lys(27), and Lys(29) linkage formation. Both interaction and ubiquitination seem to be independent from polyQ length. In cultured cells, TRAF6 enhances mutant N-HTT aggregate formation, whereas it has no effect on WT N-HTT protein localization. Mutant N-HTT inclusions are enriched for ubiquitin staining only when TRAF6 and Lys(6), Lys(27), and Lys(29) ubiquitin mutants are expressed. Finally, we show that TRAF6 is up-regulated in post-mortem brains from HD patients where it is found in the insoluble fraction. These results suggest that TRAF6 atypical ubiquitination warrants investigation in HD pathogenesis.     

Population of on-pathway intermediates in the folding of ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.     

Quantitative protein stability measurement in vivo

The equilibrium between the native and denatured states of a protein can be key to its function and regulation. Traditionally, the folding equilibrium constant has been measured in vitro using purified protein and simple buffers. However, the biological environment of proteins can differ from these in vitro conditions in ways that could significantly perturb stability. Here, we present the first quantitative comparison between the stability of a protein in vitro and in the cytoplasm of Escherichia coli using amide hydrogen exchange detected by MALDI mass spectrometry (SUPREX). The results indicate that the thermodynamic stability of monomeric lambda repressor within the cell is the same as its stability measured in a simple buffer in vitro. However, when the E. coli are placed in a hyperosmotic environment, the in vivo stability is greatly enhanced. The in vivo SUPREX method provides a general and quantitative way to measure protein stabilities in the cell and will be useful for applications where intracellular stability information provides important biological insights.     

Ultracentrifugal studies of the effect of molecular crowding by trimethylamine N-oxide on the self-association of muscle glycogen phosphorylase b.

The suitability of sedimentation equilibrium for characterizing the self-association of muscle glycogen phosphorylase b has been reappraised. Whereas sedimentation equilibrium distributions for phosphorylase b in 40 mM Hepes buffer (pH 6.8) supplemented with 1 mM AMP signify a lack of chemical equilibrium attainment, those in buffer supplemented additionally with potassium sulfate conform with the requirements of a dimerizing system in chemical as well as sedimentation equilibrium. Because the rate of attainment of chemical equilibrium under the former conditions is sufficiently slow to allow resolution of the dimeric and tetrameric enzyme species by sedimentation velocity, this procedure has been used to examine the effects of thermodynamic nonideality arising from molecular crowding by trimethylamine N-oxide on the self-association behaviour of phosphorylase b. In those terms the marginally enhanced extent of phosphorylase b self-association observed in the presence of high concentrations of the cosolute is taken to imply that the effects of thermodynamic nonideality on the dimer-tetramer equilibrium are being countered by those displacing the T<==>R isomerization equilibrium for dimer towards the smaller, nonassociating T state. Because the R state is the enzymically active form, an inhibitory effect is the predicted consequence of molecular crowding by high concentrations of unrelated solutes. Thermodynamic nonideality thus provides an alternative explanation for the inhibitory effects of high concentrations of glycerol, sucrose and ethylene glycol on phosphorylase b activity, phenomena that have been attributed to extremely weak interaction of these cryoprotectants with the T state of the enzyme.     

Energetics of side chain packing in staphylococcal nuclease assessed by systematic double mutant cycles.

All 44 possible double mutant permutations of isoleucine, leucine, and valine were constructed in 11 pairings of six sites in the core of staphylococcal nuclease. The stabilities of these mutants were determined by guanidine hydrochloride denaturation. Comparison of the stabilities of all double mutants with those expected from addition of the corresponding single mutants showed that the effects of the two single mutations are energetically independent of each other in 30 of the double mutants. However, a substantial minority, 14, of the double mutants have stability effects that are not additive. In these cases, it appears that direct van der Waals contacts between the two side chains are present. The requirement of direct van der Waals contact for the interdependence of mutational stability effects is somewhat surprising in light of results previously reported by others. In addition, it was found that double mutants that did not alter or lower the overall number of atoms in the core and that showed nonadditive behavior were more stable than expected from addition of the effects of the corresponding single mutants. A net increase in the number of atoms in the core usually, but not always, resulted in a mutant that was less stable than expected. In contrast to previous staphylococcal nuclease double mutants, energetically significant changes to the denatured state do not appear to be occurring in these packing mutants. These conclusions imply that attempts to engineer protein stability based on single mutant data will be generally successful if overall core size is preserved and if residues are not in van der Waals contact.     

Exploring structural and thermodynamic stabilities of human prion protein pathogenic mutants D202N, E211Q and Q217R

The central event in the pathogenesis of prion protein (PrP) is a profound conformational change from its α-helical (PrP(C)) to its β-sheet-rich isoform (PrP(Sc)). Many single amino acid mutations of PrP are associated with familial prion diseases, such as D202N, E211Q, and Q217R mutations located at the third native α-helix of human PrP. In order to explore the underlying structural and dynamic effects of these mutations, we performed all-atom molecular dynamics (MD) simulations for the wild-type (WT) PrP and its mutants. The obtained results indicate that these amino acid substitutions have subtle effects on the protein structures, but show large changes of the overall electrostatic potential distributions. We can infer that the changes of PrP electrostatic surface due to the studied mutations may influence the intermolecular interactions during the aggregation process. In addition, the mutations also affect the thermodynamic stabilities of PrP.     

Contribution to catalysis and stability of the five cysteines in Escherichia coli aspartate aminotransferase. Preparation and properties of a cysteine-free enzyme.

The five cysteines, at positions 82, 191, 192, 270, and 401, of Escherichia coli aspartate aminotransferase (AATase) were, individually and in some combinations, converted to alanine by site-directed mutagenesis (C82A, C191A, C192A, C270A, C401A). Cys-191, which is conserved in all AATase isozymes, was mutated to serine as well (C191S). A quintuple mutant, with all cysteines converted to alanines (Quint), was also constructed. The effects of these single and multiple mutations were examined by steady-state kinetics and urea denaturation. The thermal stabilities of Quint and of the wild-type enzyme (WT) were determined by differential scanning calorimetry. The mutants had kcat values up to 50% greater than that of WT and KMAsp and KM alpha-KG values up to 1.5- and 3.3-fold higher than that of WT. The mutants C82A and C191A exhibit nearly the same CM in urea denaturation experiments as WT, while the other single mutants and Quint are less stable, with CM differences of up to 0.7 M urea. Quint is also less thermostable than WT, with a delta TM of 3.3-4.4 degrees C. Thus the five cysteine replacements yield small, but significant, changes in catalytic and denaturation parameters, but none of the cysteines was found to be essential. The changes manifested in the mutation of the conserved Cys-191 to alanine are no greater than those observed with the four nonconserved cysteines. We consider the evolutionary implications of these findings.     

Interdomain communication between weak structural elements within a disease-related human tRNA

The structure of the human mitochondrial (hs mt) tRNALeu(UUR) features several domains that are predicted to exhibit limited thermodynamic stability. An elevated frequency of disease-related mutations within these domains suggests a link between structural instability and the functional effects of pathogenic mutations. A series of tRNAs featuring mutations within the D and anticodon stems were prepared and investigated using nuclease probing. Structural mapping studies indicated that these domains were partially denatured for the wild type (WT) hs mt tRNALeu(UUR) and were significantly stabilized by mutations introducing additional or stronger base pairs into the stem regions. In addition, trends in the aminoacylation activities of the D stem mutants suggested that the loose structure is required for function, with mutants displaying the most ordered structures exhibiting the lowest levels of aminoacylation activity. A pronounced interdependence of the structures of the anticodon and D stems was observed, with mutations strengthening the D stem stabilizing the anticodon stem and vice versa. The existence of strong interdomain communication was further elucidated with a mutant of hs mt tRNALeu(UUR) containing a stabilized D stem and a pathogenic mutation that disrupted the anticodon stem. Strengthening the structure of the D stem completely restored the function of the disease-related mutant to WT levels, indicating that propagated structural weaknesses contribute to the functional deactivation of this tRNA by mutations.     

Trigger factor assisted folding of green fluorescent protein

Guanidine induced equilibrium and kinetic folding of a variant of green fluorescent protein (F99S/M153T/V163A, GFPuv) was studied. Using manual mixing and stopped-flow techniques, we combined different probes, including tryptophan fluorescence, chromophore fluorescence and reactivity with DTNB, to trace the spontaneous and TF-assisted folding of guanidine denatured GFPuv. We found that both unfolding and refolding of GFPuv occurred in a stepwise manner and a stable intermediate was populated under equilibrium conditions. The thermodynamic parameters obtained show that the intermediate state of GFPuv is quite compact compared to the denatured state and most of the green fluorescence is retained in this state. By studying GFPuv folding assisted by TF and a number of TF mutants, we found that wild-type TF catalyzes proline isomerization and accelerates the folding rate at low TF concentrations, but retards GFPuv folding and decelerates the folding rate at high TF concentrations. This reflects the two activities of TF, as an enzyme and as a chaperone. A general mechanism of TF assisted protein folding is discussed.     

Allogeneic HLA-A*02-restricted WT1-specific T cells from mismatched donors are highly reactive but show off-target promiscuity

T cells recognizing tumor-associated Ags such as Wilms tumor protein (WT1) are thought to exert potent antitumor reactivity. However, no consistent high-avidity T cell responses have been demonstrated in vaccination studies with WT1 as target in cancer immunotherapy. The aim of this study was to investigate the possible role of negative thymic selection on the avidity and specificity of T cells directed against self-antigens. T cell clones directed against the HLA-A*0201-binding WT1(126-134) peptide were generated from both HLA-A*02-positive (self-HLA-restricted) and HLA-A*02-negative [nonself (allogeneic) HLA [allo-HLA]-restricted] individuals by direct ex vivo isolation using tetramers or after in vitro priming and selection. The functional avidity and specificity of these T cell clones was analyzed in-depth. Self-HLA-restricted WT1-specific clones only recognized WT1(126-134) with low avidities. In contrast, allo-HLA-restricted WT1 clones exhibited profound functional reactivity against a multitude of HLA-A*02-positive targets, even in the absence of exogenously loaded WT1 peptide, indicative of Ag-binding promiscuity. To characterize this potential promiscuity, reactivity of the T cell clones against 400 randomly selected HLA-A*0201-binding peptides was investigated. The self-HLA-restricted WT1-specific T cell clones only recognized the WT1 peptide. In contrast, the allo-HLA-restricted WT1-reactive clones recognized besides WT1 various other HLA-A*0201-binding peptides. In conclusion, allogeneic HLA-A*02-restricted WT1-specific T cells isolated from mismatched donors may be more tumor-reactive than their autologous counterparts but can show specific off-target promiscuity of potential clinical importance. As a result of this, administration of WT1-specific T cells generated from HLA-mismatched donors should be performed with appropriate precautions against potential off-target effects.     

Deactivation processes in PsbA1-Photosystem II and PsbA3-Photosystem II under photoinhibitory conditions in the cyanobacterium Thermosynechococcus elongatus

The sensitivity to high light conditions of Photosystem II with either PsbA1 (WT*1) or PsbA3 (WT*3) as the D1 protein was studied in whole cells of the thermophilic cyanobacterium Thermosynechococcus elongatus. When the cells are cultivated under high light conditions the following results were found: (i) The O(2) evolution activity decreases faster in WT*1 cells than in WT*3 cells both in the absence and in the presence of lincomycin, a protein synthesis inhibitor; (ii) In WT*1 cells, the rate constant for the decrease of the O(2) evolution activity is comparable in the presence and in the absence of lincomycin; (iii) The D1 content revealed by western blot analysis decays similarly in both WT*1 and WT*3 cells and much slowly than O(2) evolution; (iv) The faster decrease in O(2) evolution in WT*1 than in WT*3 cells correlates with a much faster inhibition of the S(2)-state formation; (v) The shape of the WT*1 cells is altered. All these results are in agreement with a photo-inhibition process resulting in the loss of the O(2) activity much faster than the D1 turnover in PsbA1-PSII and likely to a greater production of reactive oxygen species under high light conditions in WT*1 than in WT*3. This latter result is discussed in view of the known effects of the PsbA1 to PsbA3 substitution on the redox properties of the Photosystem II cofactors. The observation that under low light conditions WT*3 cells are able to express the psbA(3) gene, whereas under similar conditions wild type cells are expressing mainly the psbA(1) gene is also discussed. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.