Kaplan-Meier Meets Chemical Kinetics: Intrinsic Rate of SOD1 Amyloidogenesis Decreased by Subset of ALS Mutations and Cannot Fully Explain Age of Disease Onset

Toxic SOD1 trimers are off-pathway in the formation of amyloid-like fibrils in ALS

Accumulation of insoluble amyloid fibrils is widely studied as a critical factor in the pathology of multiple neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS), a fatal neurodegenerative disease. Misfolded Cu, Zn superoxide dismutase (SOD1) was the first protein linked to ALS, and non-native SOD1 trimeric oligomers were recently linked to cytotoxicity, while larger oligomers were protective to cells. The balance between trimers and larger aggregates in the process of SOD1 aggregation is, thus, a critical determinant of potential therapeutic approaches to treat ALS. Yet, it is unknown whether these trimeric oligomers are a necessary intermediate for larger aggregate formation or a distinct off-pathway species competing with fibril formation. Depending on the on- or off-pathway scenario of trimer formation, we expect drastically different therapeutic approaches. Here, we show that the toxic SOD1 trimer is an off-pathway intermediate competing with protective fibril formation. We design mutant SOD1 constructs that remain in a trimeric state (super stable trimers) and show that stabilizing the trimeric SOD1 prevents formation of fibrils in vitro and in a motor neuron like cell model (NSC-34). Using size exclusion chromatography, we track the aggregation kinetics of purified SOD1 and show direct competition of trimeric SOD1 with larger oligomer and fibril formation. Finally, we show the trimer is structurally independent of both larger soluble oligomers and insoluble fibrils using circular dichroism spectroscopy and limited proteolysis.

First Principles Calculation of Protein-Protein Dimer Affinities of ALS-Associated SOD1 Mutants

Cu,Zn superoxide dismutase (SOD1) is a 32 kDa homodimer that converts toxic oxygen radicals in neurons to less harmful species. The dimerization of SOD1 is essential to the stability of the protein. Monomerization increases the likelihood of SOD1 misfolding into conformations associated with aggregation, cellular toxicity, and neuronal death in familial amyotrophic lateral sclerosis (fALS). The ubiquity of disease-associated mutations throughout the primary sequence of SOD1 suggests an important role of physicochemical processes, including monomerization of SOD1, in the pathology of the disease. Herein, we use a first-principles statistical mechanics method to systematically calculate the free energy of dimer binding for SOD1 using molecular dynamics, which involves sequentially computing conformational, orientational, and separation distance contributions to the binding free energy. We consider the effects of two ALS-associated mutations in SOD1 protein on dimer stability, A4V and D101N, as well as the role of metal binding and disulfide bond formation. We find that the penalty for dimer formation arising from the conformational entropy of disordered loops in SOD1 is significantly larger than that for other protein–protein interactions previously considered. In the case of the disulfide-reduced protein, this leads to a bound complex whose formation is energetically disfavored. Somewhat surprisingly, the loop free energy penalty upon dimerization is still significant for the holoprotein, despite the increased structural order induced by the bound metal cations. This resulted in a surprisingly modest increase in dimer binding free energy of only about 1.5 kcal/mol upon metalation of the protein, suggesting that the most significant stabilizing effects of metalation are on folding stability rather than dimer binding stability. The mutant A4V has an unstable dimer due to weakened monomer-monomer interactions, which are manifested in the calculation by a separation free energy surface with a lower barrier. The mutant D101N has a stable dimer partially due to an unusually rigid β-barrel in the free monomer. D101N also exhibits anticooperativity in loop folding upon dimerization. These computational calculations are, to our knowledge, the most quantitatively accurate calculations of dimer binding stability in SOD1 to date.

Supercharging Prions via Amyloid-Selective Lysine Acetylation

Repulsive electrostatic forces between prion-like proteins are a barrier against aggregation. In neuropharmacology, however, a prion's net charge (Z) is not a targeted parameter. Compounds that selectively boost prion Z remain unreported. Here, we synthesized compounds that amplified the negative charge of misfolded superoxide dismutase-1 (SOD1) by acetylating lysine-NH₃ ⁺ in amyloid-SOD1, without acetylating native-SOD1. Compounds resembled a "ball and chain" mace: a rigid amyloid-binding "handle" (benzothiazole, stilbene, or styrylpyridine); an aryl ester "ball"; and a triethylene glycol chain connecting ball to handle. At stoichiometric excess, compounds acetylated up to 9 of 11 lysine per misfolded subunit (ΔZ_fibril =-8100 per 10³ subunits). Acetylated amyloid-SOD1 seeded aggregation more slowly than unacetylated amyloid-SOD1 in vitro and organotypic spinal cord (these effects were partially due to compound binding). Compounds exhibited reactivity with other amyloid and non-amyloid proteins (e.g., fibrillar α-synuclein was peracetylated; serum albumin was partially acetylated; carbonic anhydrase was largely unacetylated).

Variation in the vulnerability of mice expressing human superoxide dismutase 1 to prion-like seeding: a study of the influence of primary amino acid sequence

Misfolded forms of superoxide dismutase 1 (SOD1) with mutations associated with familial amyotrophic lateral sclerosis (fALS) exhibit prion characteristics, including the ability to act as seeds to accelerate motor neuron disease in mouse models. A key feature of infectious prion seeding is that the efficiency of transmission is governed by the primary sequence of prion protein (PrP). Isologous seeding, where the sequence of the PrP in the seed matches that of the host, is generally much more efficient than when there is a sequence mis-match. Here, we used paradigms in which mutant SOD1 seeding homogenates were injected intraspinally in newborn mice or into the sciatic nerve of adult mice, to assess the influence of SOD1 primary sequence on seeding efficiency. We observed a spectrum of seeding efficiencies depending upon both the SOD1 expressed by mice injected with seeds and the origin of the seed preparations. Mice expressing WT human SOD1 or the disease variant G37R were resistant to isologous seeding. Mice expressing G93A SOD1 were also largely resistant to isologous seeding, with limited success in one line of mice that express at low levels. By contrast, mice expressing human G85R-SOD1 were highly susceptible to isologous seeding but resistant to heterologous seeding by homogenates from paralyzed mice over-expressing mouse SOD1-G86R. In other seeding experiments with G85R SOD1:YFP mice, we observed that homogenates from paralyzed animals expressing the H46R or G37R variants of human SOD1 were less effective than seeds prepared from mice expressing the human G93A variant. These sequence mis-match effects were less pronounced when we used purified recombinant SOD1 that had been fibrilized in vitro as the seeding preparation. Collectively, our findings demonstrate diversity in the abilities of ALS variants of SOD1 to initiate or sustain prion-like propagation of misfolded conformations that produce motor neuron disease.

Amyotrophic Lateral Sclerosis: Proteins, Proteostasis, Prions, and Promises

Amyotrophic lateral sclerosis (ALS) is characterized by the progressive degeneration of the motor neurons that innervate muscle, resulting in gradual paralysis and culminating in the inability to breathe or swallow. This neuronal degeneration occurs in a spatiotemporal manner from a point of onset in the central nervous system (CNS), suggesting that there is a molecule that spreads from cell-to-cell. There is strong evidence that the onset and progression of ALS pathology is a consequence of protein misfolding and aggregation. In line with this, a hallmark pathology of ALS is protein deposition and inclusion formation within motor neurons and surrounding glia of the proteins TAR DNA-binding protein 43, superoxide dismutase-1, or fused in sarcoma. Collectively, the observed protein aggregation, in conjunction with the spatiotemporal spread of symptoms, strongly suggests a prion-like propagation of protein aggregation occurs in ALS. In this review, we discuss the role of protein aggregation in ALS concerning protein homeostasis (proteostasis) mechanisms and prion-like propagation. Furthermore, we examine the experimental models used to investigate these processes, including in vitro assays, cultured cells, invertebrate models, and murine models. Finally, we evaluate the therapeutics that may best prevent the onset or spread of pathology in ALS and discuss what lies on the horizon for treating this currently incurable disease.

Examination of SOD1 aggregation modulators and their effect on SOD1 enzymatic activity as a proxy for potential toxicity

Small-molecule inhibitors of abnormal protein self-assembly are promising candidates for developing therapy against proteinopathies. Such compounds have been examined primarily as inhibitors of amyloid β-protein (Aβ), whereas testing of inhibitors of other amyloidogenic proteins has lagged behind. An important issue with screening compound libraries is that although an inhibitor suitable for therapy must be both effective and nontoxic, typical screening focuses on efficacy, whereas safety typically is tested at a later stage using cells and/or animals. In addition, typical thioflavin T (ThT)-fluorescence-based screens use the final fluorescence value as a readout, potentially missing important kinetic information. Here, we examined potential inhibitors of superoxide dismutase 1 (SOD1) using ThT-fluorescence including the different phases of fluorescence change and added a parallel screen of SOD1 activity as a potential proxy for compound toxicity. Some compounds previously reported to inhibit other amyloidogenic proteins also inhibited SOD1 aggregation at low micromolar concentrations, whereas others were ineffective. Analysis of the lag phase and exponential slope added important information that could help exclude false-positive or false-negative results. SOD1 was highly resistant to inhibition of its activity, and therefore, did not have the necessary sensitivity to serve as a proxy for examining potential toxicity.

Loss of angiogenin function is related to earlier ALS onset and a paradoxical increase in ALS duration

0.5-1% of ALS (Amyotrophic Lateral Sclerosis) and Parkinson's disease (PD) are associated with mutations in the angiogenin (ANG). These mutations are thought to cause disease through a loss of ANG function, but this hypothesis has not been evaluated statistically. In addition, the potential for ANG to promote disease has not been considered. With the goal of better defining the etiology of ANG-ALS, we assembled all clinical onset and disease duration data and determined if these were correlated with biochemical properties of ANG variants. Loss of ANG stability and ribonuclease activity were found to correlate with early ALS onset, confirming an aspect of the prevailing model of ANG-ALS. Conversely, loss of ANG stability and ribonuclease activity correlated with longer survival following diagnosis, which is inconsistent with the prevailing model. These results indicate that functional ANG appears to decrease the risk of developing ALS but exacerbate ALS once in progress. These findings are rationalized in terms of studies demonstrating that distinct mechanisms contribute to ALS onset and progression and propose that ANG replacement or stabilization would benefit pre-symptomatic ANG-ALS patients. However, this study challenges the prevailing hypothesis that augmenting ANG will benefit symptomatic ANG-ALS patients. Instead, our results suggest that silencing of ANG activity may be beneficial for symptomatic ALS patients. This study will serve as a call-to-arms for neurologists to consistently publish ALS and PD patient's clinical data-if all ANG-ALS patients' data were available our findings could be tested with considerable statistical power.

Kinetic Variability in Seeded Formation of ALS-Linked SOD1 Fibrils Across Multiple Generations

The unseeded aggregation of superoxide dismutase-1 (SOD1) into amyloid-like fibrils occurs stochastically in vitro and in vivo, that is, isolated populations of SOD1 proteins (within microplate wells or living cells) self-assemble into amyloid at rates that span a probability distribution. This stochasticity has been attributed to variable degrees of monomer depletion by competing pathways of amorphous and fibrillar aggregation (inter alia). Here, microplate-based thioflavin-T (ThT) fluorescence assays were performed at high iteration (∼300) to establish whether this observed stochasticity persists when progenitor ("parent") SOD1 fibrils are used to seed the formation of multiple generations of progeny fibrils (daughter, granddaughter, and great-granddaughter fibrils). Populations of progenitor fibrils formed stochastically at different rates and fluorescence intensity, however, progeny fibrils formed at more similar rates regardless of the formation rate of the progenitor fibril. For example, populations of progenitor fibrils that formed with a lag time of ∼30 h or ∼15 h both produced progeny fibrils with lag times of ∼8 h. Likewise, populations of progenitor fibrils with high or low maximum fluorescence (e.g., ∼450 or ∼75 A.U.) both produced progeny fibrils with more similar maximum fluorescence (∼125 A.U.). The rate of propagation was found to be more dependent on monomer concentration than seed concentration. These results can be rationalized by classical rate laws for primary nucleation and monomer-dependent secondary nucleation. We also find that the seeding propensity of some "families" of in vitro grown fibrils exhibit a finite lifetime (similar to that observed in the seeding of small molecule crystals and colloids). The single biological takeaway of this study is that the concentration of native SOD1 in a cell can have a stronger effect on rates of seeded aggregation than the concentration of prion-like seed that infected the cell.

The biophysics of superoxide dismutase-1 and amyotrophic lateral sclerosis

Few proteins have come under such intense scrutiny as superoxide dismutase-1 (SOD1). For almost a century, scientists have dissected its form, function and then later its malfunction in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). We now know SOD1 is a zinc and copper metalloenzyme that clears superoxide as part of our antioxidant defence and respiratory regulation systems. The possibility of reduced structural integrity was suggested by the first crystal structures of human SOD1 even before deleterious mutations in the sod1 gene were linked to the ALS. This concept evolved in the intervening years as an impressive array of biophysical studies examined the characteristics of mutant SOD1 in great detail. We now recognise how ALS-related mutations perturb the SOD1 maturation processes, reduce its ability to fold and reduce its thermal stability and half-life. Mutant SOD1 is therefore predisposed to monomerisation, non-canonical self-interactions, the formation of small misfolded oligomers and ultimately accumulation in the tell-tale insoluble inclusions found within the neurons of ALS patients. We have also seen that several post-translational modifications could push wild-type SOD1 down this toxic pathway. Recently we have come to view ALS as a prion-like disease where both the symptoms, and indeed SOD1 misfolding itself, are transmitted to neighbouring cells. This raises the possibility of intervention after the initial disease presentation. Several small-molecule and biologic-based strategies have been devised which directly target the SOD1 molecule to change the behaviour thought to be responsible for ALS. Here we provide a comprehensive review of the many biophysical advances that sculpted our view of SOD1 biology and the recent work that aims to apply this knowledge for therapeutic outcomes in ALS.

High-Throughput Microplate-Based Fluorescence Assays for Studying Stochastic Aggregation of Superoxide Dismutase-1

Investigating in vitro kinetics of superoxide dismutase-1 (SOD1) aggregation with high-throughput microplate-based assays provides valuable information regarding SOD1 pathogenesis in amyotrophic lateral sclerosis (ALS) and opens venues for the development of effective therapies. In this chapter, we first explain the step-by-step purification and demetallation of wild-type (WT) and ALS-variant SOD1 proteins from Saccharomyces cerevisiae (baker's yeast). We then describe the methodology for a microplate-based fluorescence assay that is used to study real-time kinetics of metal-free (apo)-SOD1 aggregation. This technique is highly sensitive, semiautomated, requires minimum modifications to protein, and produces a plethora of data in a short period of time. We also describe a new approach for extracting clinically relevant information from SOD1 aggregation data using Kaplan-Meier estimators.

Investigation of Anti-SOD1 Antibodies Yields New Structural Insight into SOD1 Misfolding and Surprising Behavior of the Antibodies Themselves

Mutations in Cu/Zn-superoxide dismutase (SOD1) gene are linked to 10-20% of familial amyotrophic lateral sclerosis (fALS) cases. The mutations cause misfolding and self-assembly of SOD1 into toxic oligomers and aggregates, resulting in motor neuron degeneration. The molecular mechanisms underlying SOD1 aggregation and toxicity are unclear. Characterization of misfolded SOD1 is particularly challenging because of its metastable nature. Antibodies against misfolded SOD1 are useful tools for this purpose, provided their specificity and selectivity are well-characterized. Here, we characterized three recently introduced antimisfolded SOD1 antibodies and compared them with two commercial, antimisfolded SOD1 antibodies raised against the fALS-linked variant G93A-SOD1. As controls, we compared the reactivity of these antibodies to two polyclonal anti-SOD1 antibodies expected to be insensitive to misfolding. We asked to what extent the antibodies could distinguish between WT and variant SOD1 and between native and misfolded conformations. WT, G93A-SOD1, or E100K-SOD1 were incubated under aggregation-promoting conditions and monitored using thioflavin-T fluorescence, electron microscopy, and dot blots. WT and G93A-SOD1 also were analyzed using native-PAGE/Western blot. The new antimisfolded SOD1 and the commercial antibody B8H10 showed variable reactivity using dot blots but generally showed maximum reactivity at the time misfolded SOD1 oligomers were expected to be most abundant. In contrast, only B8H10 and the control antibodies were reactive in Western blots. Unexpectedly, the polyclonal antibodies showed strong preference for the misfolded form of G93A-SOD1 in dot blots. Surprisingly, antimisfolded SOD1 antibody C4F6 was specific for the apo form of G93A-SOD1 but insensitive to misfolding. Antibody 10C12 showed preference for early misfolded structures, whereas 3H1 bound preferentially to late structures. These new antibodies allow distinction between putative early- and late-forming prefibrillar SOD1 oligomers.

Glycerolipid Headgroups Control Rate and Mechanism of Superoxide Dismutase-1 Aggregation and Accelerate Fibrillization of Slowly Aggregating Amyotrophic Lateral Sclerosis Mutants

Interactions between superoxide dismutase-1 (SOD1) and lipid membranes might be directly involved in the toxicity and intercellular propagation of aggregated SOD1 in amyotrophic lateral sclerosis (ALS), but the chemical details of lipid-SOD1 interactions and their effects on SOD1 aggregation remain unclear. This paper determined the rate and mechanism of nucleation of fibrillar apo-SOD1 catalyzed by liposomal surfaces with identical hydrophobic chains (RCH₂(O₂C₁₈H₃₃)₂), but headgroups of different net charge and hydrophobicity (i.e., R(CH₂)N⁺(CH₃)₃, RPO₄^-(CH₂)₂N⁺(CH₃)₃, and RPO₄^-). Under semiquiescent conditions (within a 96 well microplate, without a gyrating bead), the aggregation of apo-SOD1 into thioflavin-T-positive (ThT(+)) amyloid fibrils did not occur over 120 h in the absence of liposomal surfaces. Anionic liposomes triggered aggregation of apo-SOD1 into ThT(+) amyloid fibrils; cationic liposomes catalyzed fibrillization but at slower rates and across a narrower lipid concentration; zwitterionic liposomes produced nonfibrillar (amorphous) aggregates. The inability of zwitterionic liposomes to catalyze fibrillization and the dependence of fibrillization rate on anionic lipid concentration suggests that membranes catalyze SOD1 fibrillization by a primary nucleation mechanism. Membrane-catalyzed fibrillization was also examined for eight ALS variants of apo-SOD1, including G37R, G93R, D90A, and E100G apo-SOD1 that nucleate slower than or equal to WT SOD1 in lipid-free, nonquiescent amyloid assays. All ALS variants (with one exception) nucleated faster than WT SOD1 in the presence of anionic liposomes, wherein the greatest acceleratory effects were observed among variants with lower net negative surface charge (G37R, G93R, D90A, E100G). The exception was H46R apo-SOD1, which did not form ThT(+) species.

Large SOD1 aggregates, unlike trimeric SOD1, do not impact cell viability in a model of amyotrophic lateral sclerosis

Aberrant accumulation of misfolded Cu, Zn superoxide dismutase (SOD1) is a hallmark of SOD1-associated amyotrophic lateral sclerosis (ALS), an invariably fatal neurodegenerative disease. While recent discovery of nonnative trimeric SOD1-associated neurotoxicity has suggested a potential pathway for motor neuron impairment, it is yet unknown whether large, insoluble aggregates are cytotoxic. Here we designed SOD1 mutations that specifically stabilize either the fibrillar form or the trimeric state of SOD1. The designed mutants display elevated populations of fibrils or trimers correspondingly, as demonstrated by gel filtration chromatography and electron microscopy. The trimer-stabilizing mutant, G147P, promoted cell death, even more potently in comparison with the aggressive ALS-associated mutants A4V and G93A. In contrast, the fibril-stabilizing mutants, N53I and D101I, positively impacted the survival of motor neuron-like cells. Hence, we conclude the SOD1 oligomer and not the mature form of aggregated fibril is critical for the neurotoxic effects in the model of ALS. The formation of large aggregates is in competition with trimer formation, suggesting that aggregation may be a protective mechanism against formation of toxic oligomeric intermediates.

Misfolded SOD1 Accumulation and Mitochondrial Association Contribute to the Selective Vulnerability of Motor Neurons in Familial ALS: Correlation to Human Disease

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder, with a 10% genetic linkage, of which 20% of these cases may be attributed to mutations in superoxide dismutase (SOD1). Specific mutations in SOD1 have been associated with disease duration, which can be highly variable ranging from a life expectancy of 3 to beyond 10 years. SOD1 neurotoxicity has been attributed to aberrant accumulation of misfolded SOD1, which in its soluble form binds to intracellular organelles disrupting their function or forms insoluble toxic aggregates. To understand whether these biophysical properties of the mutant protein may influence disease onset and duration, we generated 19 point mutations in the SOD1 gene, based on available clinical data of disease onset and progression from patients. By overexpressing these mutants in motor-neuron-like NSC-34 cells, we demonstrate a variability in misfolding capacity between the different mutants with a correlation between the degree of protein misfolding and mutation severity. We also show a clear variation of the different SOD1 mutants to associate with mitochondrial-enriched fractions with a correlation between mutation severity and this association. In summary, these findings reveal a correlation between the accumulation of misfolded SOD1 species and their mitochondrial association with disease duration but not with disease onset, and they have implications for the potential therapeutic role of suppressing the accumulation of misfolded SOD1.

Lysine acylation in superoxide dismutase-1 electrostatically inhibits formation of fibrils with prion-like seeding

The acylation of lysine residues in superoxide dismutase-1 (SOD1) has been previously shown to decrease its rate of nucleation and elongation into amyloid-like fibrils linked to amyotrophic lateral sclerosis. The chemical mechanism underlying this effect is unclear, i.e. hydrophobic/steric effects versus electrostatic effects. Moreover, the degree to which the acylation might alter the prion-like seeding of SOD1 in vivo has not been addressed. Here, we acylated a fraction of lysine residues in SOD1 with groups of variable hydrophobicity, charge, and conformational entropy. The effect of each acyl group on the rate of SOD1 fibril nucleation and elongation were quantified in vitro with thioflavin-T (ThT) fluorescence, and we performed 594 iterate aggregation assays to obtain statistically significant rates. The effect of the lysine acylation on the prion-like seeding of SOD1 was assayed in spinal cord extracts of transgenic mice expressing a G85R SOD1-yellow fluorescent protein construct. Acyl groups with >2 carboxylic acids diminished self-assembly into ThT-positive fibrils and instead promoted the self-assembly of ThT-negative fibrils and amorphous complexes. The addition of ThT-negative, acylated SOD1 fibrils to organotypic spinal cord failed to produce the SOD1 inclusion pathology that typically results from the addition of ThT-positive SOD1 fibrils. These results suggest that chemically increasing the net negative surface charge of SOD1 via acylation can block the prion-like propagation of oligomeric SOD1 in spinal cord.

The unfolding mechanism of monomeric mutant SOD1 by simulated force spectroscopy

Mechanical unfolding of mutated apo, disulfide-reduced, monomeric superoxide dismutase 1 protein (SOD1) has been simulated via force spectroscopy techniques, using both an all-atom (AA), explicit solvent model and a coarse-grained heavy-atom Gō (HA-Gō) model. The HA-Gō model was implemented at two different pulling speeds for comparison. The most-common sequence of unfolding in the AA model agrees well with the most-common unfolding sequence of the HA-Gō model, when the same normalized pulling rate was used. Clustering of partially-native structures as the protein unfolds shows that the AA and HA-Gō models both exhibit a dominant pathway for early unfolding, which eventually bifurcates repeatedly to multiple branches after the protein is about half-unfolded. The force-extension curve exhibits multiple force drops, which are concomitant with jumps in the local interaction potential energy between specific β-strands in the protein. These sudden jumps in the potential energy coincide with the dissociation of specific pairs of β-strands, and thus intermediate unfolding events. The most common sequence of β-strand dissociation in the unfolding pathway of the AA model is β-strands 5, 4, 8, 7, 1, 2, then finally β-strands 3 and 6. The observation that β-strand 5 is among the first to unfold here, but the last to unfold in simulations of loop-truncated SOD1, could imply the existence of an evolutionary compensation mechanism, which would stabilize β-strands flanking long loops against their entropic penalty by strengthening intramolecular interactions. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.