Disulfide-Reduced ALS Variants of Cu, Zn Superoxide Dismutase Exhibit Increased Populations of Unfolded Species

Cu,Zn superoxide dismutase (SOD1) is a dimeric metal-binding enzyme responsible for the dismutation of toxic superoxide to hydrogen peroxide and oxygen in cells. Mutations at dozens of sites in SOD1 induce amyotrophic lateral sclerosis (ALS), a fatal gain-of-function neurodegenerative disease whose molecular basis is unknown. To obtain insights into effects of the mutations on the folded and unfolded populations of immature monomeric forms whose aggregation or self-association may be responsible for ALS, the thermodynamic and kinetic folding properties of a set of disulfide-reduced and disulfide-oxidized Zn-free and Zn-bound stable monomeric SOD1 variants were compared to properties of the wild-type (WT) protein. The most striking effect of the mutations on the monomer stability was observed for the disulfide-reduced metal-free variants. Whereas the WT and S134N monomers are > 95% folded at neutral pH and 37 °C, A4V, L38V, G93A, and L106V ranged from 50% to ∼ 90% unfolded. The reduction of the disulfide bond was also found to reduce the apparent Zn affinity of the WT monomer by 750-fold, into the nanomolar range, where it may be unable to compete for free Zn in the cell. With the exception of the S134N metal-binding variant, the Zn affinity of disulfide-oxidized SOD1 monomers showed little sensitivity to amino acid replacements. These results suggest a model for SOD1 aggregation where the constant synthesis of ALS variants of SOD1 on ribosomes provides a pool of species in which the increased population of the unfolded state may favor aggregation over productive folding to the native dimeric state.     

Zinc binding modulates the entire folding free energy surface of human Cu,Zn superoxide dismutase

Over 100 amino acid replacements in human Cu,Zn superoxide dismutase (SOD) are known to cause amyotrophic lateral sclerosis, a gain-of-function neurodegenerative disease that destroys motor neurons. Supposing that aggregates of partially folded states are primarily responsible for toxicity, we determined the role of the structurally important zinc ion in defining the folding free energy surface of dimeric SOD by comparing the thermodynamic and kinetic folding properties of the zinc-free and zinc-bound forms of the protein. The presence of zinc was found to decrease the free energies of a peptide model of the unfolded monomer, a stable variant of the folded monomeric intermediate, and the folded dimeric species. The unfolded state binds zinc weakly with a micromolar dissociation constant, and the folded monomeric intermediate and the native dimeric form both bind zinc tightly, with subnanomolar dissociation constants. Coupled with the strong driving force for the subunit association reaction, the shift in the populations toward more well-folded states in the presence of zinc decreases the steady-state populations of higher-energy states in SOD under expected in vivo zinc concentrations (approximately nanomolar). The significant decrease in the population of partially folded states is expected to diminish their potential for aggregation and account for the known protective effect of zinc. The ∼ 100-fold increase in the rate of folding of SOD in the presence of micromolar concentrations of zinc demonstrates a significant role for a preorganized zinc-binding loop in the transition-state ensemble for the rate-limiting monomer folding reaction in this β-barrel protein.     

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.     

The Folding Free-Energy Surface of HIV-1 Protease: Insights into the Thermodynamic Basis for Resistance to Inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric β-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR^?, was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.     

Reversible denaturation of the gene V protein of bacteriophage f1

The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.     

Backbone resonance assignments of monomeric SOD1 in dilute and crowded environments

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease that leads to movement disorders. In motor neurons of ALS patients, intracellular aggregates of superoxide dismutase 1 (SOD1) have often been observed. To elucidate the aggregation mechanism, it is important to analyze the folding equilibrium of SOD1 between folded and aggregation-prone unfolded states. However, in most cases, this folding equilibrium has been studied in dilute solution even though the aggregate formation occurs in a highly crowded intracellular environment. Indeed, a recent study reported that the folding stability of SOD1 decreased in an environment containing protein crowder molecules. To understand such a destabilization effect due to protein crowders, it is necessary to obtain more precise structural information on SOD1 in the presence of protein crowders. Here, we report the ¹H, ¹³C, and ¹⁵N backbone resonance assignments of monomeric SOD1 in the absence and presence of the protein crowder lysozyme. The chemical shift differences caused by addition of lysozyme suggest that SOD1 associated with lysozyme via negatively charged surfaces. Based on the assigned chemical shifts, the presence of lysozyme has a limited influence on the secondary structure of SOD1. We anticipate that our assignments will provide an important basis for elucidation of the crowding-induced folding destabilization of SOD1.     

Increased stability and DNA site discrimination of "single chain" variants of the dimeric beta-barrel DNA binding domain of the human papillomavirus E2 transcriptional regulator

Human papillomavirus infects millions of people worldwide and is a causal agent of cervical cancer in women. The HPV E2 protein controls the expression of all viral genes through binding of its dimeric C-terminal domain (E2C) to its target DNA site. We engineered monomeric versions of the HPV16 E2C, in order to probe the link of the dimeric beta-barrel fold to stability, dimerization, and DNA binding. Two single-chain variants, with 6 and 12 residue linkers (scE2C-6 and scE2C-12), were purified and characterized. Spectroscopy and crystallography show that the native structure is unperturbed in scE2C-12. The single chain variants are stabilized with respect to E2C, with effective concentrations of 0.6 to 6 mM. The early folding events of the E2C dimer and scE2C-12 are very similar and include formation of a compact species in the submillisecond time scale and a non-native monomeric intermediate with a half-life of 25 ms. However, monomerization changes the unfolding mechanism of the linked species from two-state to three-state, with a high-energy intermediate. Binding to the specific target site is up to 5-fold tighter in the single chain variants. Nonspecific DNA binding is up to 7-fold weaker in the single chain variants, leading to an overall 10-fold increased site discrimination capacity, the largest described so far for linked DNA binding domains. Titration calorimetric binding analysis, however, shows almost identical behavior for dimer and single-chain species, suggesting very subtle changes behind the increased specificity. Global analysis of the mechanisms probed suggests that the dynamics of the E2C domain, rather than the structure, are responsible for the differential properties. Thus, the plastic and dimeric nature of the domain did not evolve for a maximum affinity, specificity, and stability of the quaternary structure, likely because of regulatory reasons and for roles other than DNA binding played by partly folded dimeric or monomeric conformers.     

The Ribosome Can Prevent Aggregation of Partially Folded Protein Intermediates: Studies Using the Escherichia coli Ribosome

Background

Molecular chaperones that support de novo folding of proteins under non stress condition are classified as chaperone ‘foldases’ that are distinct from chaperone’ holdases’ that provide high affinity binding platform for unfolded proteins and prevent their aggregation specifically under stress conditions. Ribosome, the cellular protein synthesis machine can act as a foldase chaperone that can bind unfolded proteins and release them in folding competent state. The peptidyl transferase center (PTC) located in the domain V of the 23S rRNA of Escherichia coli ribosome (bDV RNA) is the chaperoning center of the ribosome. It has been proposed that via specific interactions between the RNA and refolding proteins, the chaperone provides information for the correct folding of unfolded polypeptide chains.

Results

We demonstrate using Escherichia coli ribosome and variants of its domain V RNA that the ribosome can bind to partially folded intermediates of bovine carbonic anhydrase II (BCAII) and lysozyme and suppress aggregation during their refolding. Using mutants of domain V RNA we demonstrate that the time for which the chaperone retains the bound protein is an important factor in determining its ability to suppress aggregation and/or support reactivation of protein.

Conclusion

The ribosome can behave like a ‘holdase’ chaperone and has the ability to bind and hold back partially folded intermediate states of proteins from participating in the aggregation process. Since the ribosome is an essential organelle that is present in large numbers in all living cells, this ability of the ribosome provides an energetically inexpensive way to suppress cellular aggregation. Further, this ability of the ribosome might also be crucial in the context that the ribosome is one of the first chaperones to be encountered by a large nascent polypeptide chains that have a tendency to form partially folded intermediates immediately following their synthesis.     

Local unfolding of Cu, Zn superoxide dismutase monomer determines the morphology of fibrillar aggregates

Aggregation of Cu, Zn superoxide dismutase (SOD1) is often found in amyotrophic lateral sclerosis patients. The fibrillar aggregates formed by wild type and various disease-associated mutants have recently been found to have distinct cores and morphologies. Previous computational and experimental studies of wild-type SOD1 suggest that the apo-monomer, highly aggregation prone, displays substantial local unfolding dynamics. The residual folded structure of locally unfolded apoSOD1 corresponds to peptide segments forming the aggregation core as identified by a combination of proteolysis and mass spectroscopy. Therefore, we hypothesize that the destabilization of apoSOD1 caused by various mutations leads to distinct local unfolding dynamics. The partially unfolded structure, exposing the hydrophobic core and backbone hydrogen bond donors and acceptors, is prone to aggregate. The peptide segments in the residual folded structures form the "building block" for aggregation, which in turn determines the morphology of the aggregates. To test this hypothesis, we apply a multiscale simulation approach to study the aggregation of three typical SOD1 variants: wild type, G37R, and I149T. Each of these SOD1 variants has distinct peptide segments forming the core structure and features different aggregate morphologies. We perform atomistic molecular dynamics simulations to study the conformational dynamics of apoSOD1 monomer and coarse-grained molecular dynamics simulations to study the aggregation of partially unfolded SOD1 monomers. Our computational studies of monomer local unfolding and the aggregation of different SOD1 variants are consistent with experiments, supporting the hypothesis of the formation of aggregation "building blocks" via apo-monomer local unfolding as the mechanism of SOD1 fibrillar aggregation.     

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the β-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the β-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the β-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the β-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Separating instability from aggregation propensity in γS-crystallin variants

Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.     

Effect of backbone cyclization on protein folding stability: chain entropies of both the unfolded and the folded states are restricted

Circular versions of a large number of proteins have been designed by connecting the N and C termini via peptide linkers. A motivation for these designs is the assumed enhancement in folding stability, because backbone cyclization reduces the chain entropy of the unfolded state. Here, it is recognized that backbone cyclization also reduces the chain entropy of a flexible peptide linker in the folded state. Specifically, the end-to-end distance of the linker is restricted to fluctuations around the average displacement between the N and C termini of the folded protein. The balance of the chain-entropy reductions in the folded and unfolded states is used to predict the change in the unfolding free energy, deltadeltaG(cycl), by backbone cyclization. Predicted values of deltadeltaG(cycl) are in quantitative agreement with results of a careful study on cyclizing the 34 residue PIN1 WW domain by linkers with two to seen residues. The experimental results of an optimal linker length l=4 and a maximum stabilization of 1.7 kcal/mol are reproduced. Calculations of deltadeltaG(cycl) for a broad selection of circular proteins suggest that the stabilizing effect of backbone cyclization is modest, reflecting entropy reductions in both the unfolded and the folded states.     

Abp1p and Fyn SH3 domains fold through similar low-populated intermediate states

Src homology 3 (SH3) domains are small modules that are thought to fold via a two-state mechanism, without the accumulation of significant populations of intermediate states. Relaxation dispersion NMR studies of the folding of G48V and G48M mutants of the Fyn SH3 domain have established that, at least for these modules, folding proceeds through the formation of a transient on-pathway intermediate with an equilibrium population of 1-2% that can be readily detected [Korzhnev, D. M., et al. (2004) Nature 430, 586-590]. To investigate the generality of this result, we present an (15)N relaxation dispersion NMR study of a pair of additional SH3 domains, including a G48V mutant of a stabilized Abp1p SH3 domain that shares 36% sequence identity with the Fyn SH3 module, and a A39V/N53P/V55L mutant Fyn SH3 domain. A transient folding intermediate is detected for both of the proteins studied here, and the dispersion data are well fit to a folding model of the form F <--> I <--> U, where F, I, and U correspond to folded, intermediate, and unfolded states, respectively. The temperature dependencies of the folding/unfolding rate constants were obtained so that the thermodynamic properties of each of F, I, and U could be established. The detection of I states in folding pathways of all SH3 domains examined to date via relaxation dispersion NMR spectroscopy indicates that such intermediates may well be a conserved feature in the folding of such domains in general but that their transient nature along with their low population makes detection difficult using more well-established approaches to the study of folding.     

Long-range electrostatic effects on peptide folding.

The reversible folding/unfolding of a short peptide in solution is studied by molecular dynamics simulations. The effects of long-range electrostatic interactions are examined and found to be important both for the equilibrium between folded and unfolded states and the dynamics of the folding process. The neglect of long-range electrostatics leads to an increased population of unfolded states and increased structural fluctuations. When such interactions are taken into account, the peptide unfolds and folds to the experimentally determined structure several times during a 25 ns trajectory, with approximately equal populations of folded and unfolded states in the neighborhood of its proposed melting temperature. The effect of using spherical boundary conditions rather than periodic ones does not appear to have any major effect on the folding dynamics.     

Flipping the switch from monomeric to dimeric CV-N has little effect on antiviral activity

Cyanovirin-N can exist in solution in monomeric and domain-swapped dimeric forms, with HIV-antiviral activity being reported for both. Here we present results for CV-N variants that form stable solution dimers: the obligate dimer [DeltaQ50]CV-N and the preferential dimer [S52P]CV-N. These variants exhibit comparable DeltaG values (10.6 +/- 0.5 and 9.4 +/- 0.5 kcal.mol(-1), respectively), similar to that of stabilized, monomeric [P51G]CV-N (9.8 +/- 0.5 kcal.mol(-1)), but significantly higher than wild-type CV-N (4.1 +/- 0.2 kcal.mol(-1)). During folding/unfolding, no stably folded monomer was observed under any condition for the obligate dimer [DeltaQ50]CV-N, whereas two monomeric, metastable species were detected for [S52P]CV-N at low concentrations. This is in contrast to our previous results for [P51G]CV-N and wild-type CV-N, for which the dimeric forms were found to be the metastable species. The dimeric mutants exhibit comparable antiviral activity against HIV and Ebola, similar to that of wild-type CV-N and the stabilized [P51G]CV-N variant.     

Unfolding and folding kinetics of amyotrophic lateral sclerosis-associated mutant Cu,Zn superoxide dismutases

More than 110 mutations in dimeric, Cu,Zn superoxide dismutase (SOD) have been linked to the fatal neurodegenerative disease, amyotrophic lateral sclerosis (ALS). In both human patients and mouse model studies, protein misfolding has been implicated in disease pathogenesis. A central step in understanding the misfolding/aggregation mechanism of this protein is the elucidation of the folding pathway of SOD. Here we report a systematic analyses of unfolding and folding kinetics using single- and double-jump experiments as well as measurements as a function of guanidium chloride, protein, and metal concentration for fully metallated (holo) pseudo wild-type and ALS-associated mutant (E100G, G93R, G93A, and metal binding mutants G85R and H46R) SODs. The kinetic mechanism for holo SODs involves native dimer, monomer intermediate, and unfolded monomer, with variable metal dissociation from the monomeric states depending on solution conditions. The effects of the ALS mutations on the kinetics of the holoproteins in guanidium chloride are markedly different from those observed previously for acid-induced unfolding and for the unmetallated (apo) forms of the proteins. The mutations decrease the stability of holo SOD mainly by increasing unfolding rates, which is particularly pronounced for the metal-binding mutants, and have relatively smaller effects on the observed folding kinetics. Mutations also seem to favour increased formation of a Zn-free monomer intermediate, which has been implicated in the formation of toxic aggregates. The results reveal the kinetic basis for the extremely high stability of wild-type holo SOD and the possible consequences of kinetic changes for disease.