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

Selenium- and zinc-biofortified bean sprouts improve cognitive dysfunction in aging mice by reducing oxidative stress

Currently, age-related diseases are a serious global health problem. In this study, selenium- and zinc-biofortified bean sprouts (ZSB) were cultivated, and the intervention effect on cognitive dysfunction was explored in aging mice model established by combining AlCl3 and d-galactose treatment. The results showed that compared with the model group, ZSB could decrease the accumulation of Al and increase the accumulation of Zn and Se in mice. The levels of serum superoxide dismutase, glutathione, and glutathione peroxidase in liver and glutathione and acetylcholine in brain were significantly increased (p < 0.05) after the treatments of biofortified bean sprouts. Additionally, biofortified bean sprouts treatment reduced the level of lipid peroxidation and inhibited the deposition of Aβ1-40 in brain and the activity of acetylcholinesterase. It could also increase the abundance of lactic acid bacteria in cecum and inhibit the proliferation of Escherichia coli, and then the content of short-chain fatty acids was increased. More importantly, hippocampal neuron damage was reversed, and the correct rate of spontaneous behavior alternations was increased by 35% after the intervention of biofortified bean sprouts. Overall, ZSB can effectively slow down aging process through a variety of ways, which provide countermeasures for delaying aging.

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

Nicotinamide Adenine Dinucleotide Precursor Supplementation Modulates Neurite Complexity and Survival in Motor Neurons from Amyotrophic Lateral Sclerosis Models

Aims: Increasing nicotinamide adenine dinucleotide (NAD+) availability has been proposed as a therapeutic approach to prevent neurodegeneration in amyotrophic lateral sclerosis (ALS). Accordingly, NAD+ precursor supplementation appears to exert neuroprotective effects in ALS patients and mouse models. The mechanisms mediating neuroprotection remain uncertain but could involve changes in multiple cell types. We investigated a potential direct effect of the NAD+ precursor nicotinamide mononucleotide (NMN) on the health of cultured induced pluripotent stem cell (iPSC)-derived human motor neurons and in motor neurons isolated from two ALS mouse models, that is, mice overexpressing wild-type transactive response DNA binding protein-43 (TDP-43) or the ALS-linked human superoxide dismutase 1 with the G93A mutation (hSOD1G93A). Results: NMN treatment increased the complexity of neuronal processes in motor neurons isolated from both mouse models and in iPSC-derived human motor neurons. In addition, NMN prevented neuronal death induced by trophic factor deprivation. In mouse and human motor neurons expressing ALS-linked mutant superoxide dismutase 1, NMN induced an increase in glutathione levels, but this effect was not observed in nontransgenic or TDP-43 overexpressing motor neurons. In contrast, NMN treatment normalized the TDP-43 cytoplasmic mislocalization induced by its overexpression. Innovation: NMN can directly act on motor neurons to increase the growth and complexity of neuronal processes and prevent the death induced by trophic factor deprivation. Conclusion: Our results support a direct beneficial effect of NAD+ precursor supplementation on the maintenance of the neuritic arbor in motor neurons. Importantly, this was observed in motor neurons isolated from two different ALS models, with and without involvement of TDP-43 pathology, supporting its therapeutic potential in sporadic and familial ALS. Antioxid. Redox Signal. 41, 573-589.

Interactions of binary mixtures of metals on rainbow trout (Oncorhynchus mykiss) heart mitochondrial H2O2 homeodynamics

For continuous pumping of blood, the heart needs a constant supply of energy (ATP) that is primarily met via oxidative phosphorylation in the mitochondria of cardiomyocytes. However, sustained high rates of electron transport for energy conversion redox reactions predisposes the heart to the production of reactive oxygen species (ROS) and oxidative stress. Mitochondrial ROS are fundamental drivers of responses to environmental stressors including metals but knowledge of how combinations of metals alter mitochondrial ROS homeodynamics remains sparse. We explored the effects and interactions of binary mixtures of copper (Cu), cadmium (Cd), and zinc (Zn), metals that are common contaminants of aquatic systems, on ROS (hydrogen peroxide, H2O2) homeodynamics in rainbow trout (Oncorhynchus mykiss) heart mitochondria. Isolated mitochondria were energized with glutamate-malate or succinate and exposed to a range of concentrations of the metals singly and in equimolar binary concentrations. Speciation analysis revealed that Cu was highly complexed by glutamate or Tris resulting in Cu2+ concentrations in the picomolar to nanomolar range. The concentration of Cd2+ was 7.2-7.5 % of the total while Zn2+ was 15 % and 21 % of the total during glutamate-malate and succinate oxidation, respectively. The concentration-effect relationships for Cu and Cd on mitochondrial H2O2 emission depended on the substrate while those for Zn were similar during glutamate-malate and succinate oxidation. Cu + Zn and Cu + Cd mixtures exhibited antagonistic interactions wherein Cu reduced the effects of both Cd and Zn, suggesting that Cu can mitigate oxidative distress caused by Cd or Zn. Binary combinations of the metals acted additively to reduce the rate constant and increase the half-life of H2O2 consumption while concomitantly suppressing thioredoxin reductase and stimulating glutathione peroxidase activities. Collectively, our study indicates that binary mixtures of Cu, Zn, and Cd act additively or antagonistically to modulate H2O2 homeodynamics in heart mitochondria.

The physiological and pathophysiological roles of copper in the nervous system

Copper is a critical trace element in biological systems due the vast number of essential enzymes that require the metal as a cofactor, including cytochrome c oxidase, superoxide dismutase and dopamine-β-hydroxylase. Due its key role in oxidative metabolism, antioxidant defence and neurotransmitter synthesis, copper is particularly important for neuronal development and proper neuronal function. Moreover, increasing evidence suggests that copper also serves important functions in synaptic and network activity, the regulation of circadian rhythms, and arousal. However, it is important to note that because of copper's ability to redox cycle and generate reactive species, cellular levels of the metal must be tightly regulated to meet cellular needs while avoiding copper-induced oxidative stress. Therefore, it is essential that the intricate system of copper transporters, exporters, copper chaperones and copper trafficking proteins function properly and in coordinate fashion. Indeed, disorders of copper metabolism such as Menkes disease and Wilson disease, as well as diseases linked to dysfunction of copper-requiring enzymes, such as SOD1-linked amyotrophic lateral sclerosis, demonstrate the dramatic neurological consequences of altered copper homeostasis. In this review, we explore the physiological importance of copper in the nervous system as well as pathologies related to improper copper handling.

Discovery of a novel homozygous SOD1 truncating variant bolsters infantile SOD1 deficiency syndrome

OBJECTIVE

Superoxide dismutase 1 (SOD1) is an important antioxidant enzyme whose main function is to neutralise superoxide free radicals in the cytoplasm. Heterozygous variants in SOD1 are responsible for a substantial percentage of familial amyotrophic lateral sclerosis (ALS) cases. Recently, several reports have shown that biallelic loss of SOD1 function results in a novel phenotype called infantile SOD1 deficiency syndrome, which is consistent with a recessive pattern of inheritance and can be distinguished from typical (adult-onset) ALS.

METHODS

We documented detailed family histories and clinical data, followed by whole-exome sequencing and family co-segregation analysis through Sanger sequencing. To facilitate comparisons, relevant data from fifteen previously reported patients with SOD1-related neurodevelopmental disorders were included.

RESULTS

This study presents a new Turkish family with two affected children exhibiting severe delayed motor development, infancy-onset loss of motor skills, axial hypotonia, tetraspasticity, and impaired cognitive functions. Genetic analysis revealed a novel homozygous frameshift variant in SOD1 (c.248dupG [p.Asp84Argfs*8]), with computational biochemical studies shedding light on the mechanistic aspects of SOD1 dysfunction.

CONCLUSIONS

Our findings contribute an affirmative report of a fourth biallelic variant resulting in a severe clinical phenotype, reminiscent of those induced by previously identified homozygous loss-of-function SOD1 variants. This research not only advances our understanding of the pathogenesis of this debilitating neurological syndrome but also aligns with ongoing intensive efforts to comprehend and address SOD1-linked ALS.

Extreme diversity of 12 cations in folding ALS-linked hSOD1 unveils novel hSOD1-dependent mechanisms for Fe2+/Cu2+-induced cytotoxicity

153-Residue copper-zinc superoxide dismutase 1 (hSOD1) is the first gene whose mutation was linked to FALS. To date, > 180 ALS-causing mutations have been identified within hSOD1, yet the underlying mechanism still remains mysterious. Mature hSOD1 is exceptionally stable constrained by a disulfide bridge to adopt a Greek-key β-barrel fold that accommodates copper/zinc cofactors. Conversely, nascent hSOD1 is unfolded and susceptible to aggregation and amyloid formation, requiring Zn2+ to initiate folding to a coexistence of folded and unfolded states. Recent studies demonstrate mutations that disrupt Zn2+-binding correlate with their ability to form toxic aggregates. Therefore, to decode the role of cations in hSOD1 folding provides not only mechanistic insights, but may bear therapeutic implications for hSOD1-linked ALS. Here by NMR, we visualized the effect of 12 cations: 8 essential for humans (Na+, K+, Ca2+, Zn2+, Mg2+, Mn2+, Cu2+, Fe2+), 3 mimicking zinc (Ni2+, Cd2+, Co2+), and environmentally abundant Al3+. Surprisingly, most cations, including Zn2+-mimics, showed negligible binding or induction for folding of nascent hSOD1. Cu2+ exhibited extensive binding to the unfolded state but led to severe aggregation. Unexpectedly, for the first time Fe2+ was deciphered to have Zn2+-like folding-inducing capacity. Zn2+ was unable to induce folding of H80S/D83S-hSOD1, while Fe2+ could. In contrast, Zn2+ could trigger folding of G93A-hSOD1, but Fe2+ failed. Notably, pre-existing Fe2+ disrupted the Zn2+-induced folding of G93A-hSOD1. Comparing with the ATP-induced folded state, our findings delineate that hSOD1 maturation requires: (1) intrinsic folding capacity encoded by the sequence; (2) specific Zn2+-coordination; (3) disulfide formation and Cu-load catalyzed by hCCS. This study unveils a previously-unknown interplay of cations in governing the initial folding of hSOD1, emphasizing the pivotal role of Zn2+ in hSOD1-related ALS and implying new hSOD1-dependent mechanisms for Cu2+/Fe2+-induced cytotoxicity, likely relevant to aging and other diseases.

Misfolding-Associated Exposure of Natively Buried Residues in Mutant SOD1 Facilitates Binding to TRAF6

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease primarily impacting motor neurons. Mutations in superoxide dismutase 1 (SOD1) are the second most common cause of familial ALS. Several of these mutations lead to misfolding or toxic gain of function in the SOD1 protein. Recently, we reported that misfolded SOD1 interacts with TNF receptor-associated factor 6 (TRAF6) in the SOD1^G93A rat model of ALS. Further, we showed in cultured cells that several mutant SOD1 proteins, but not wildtype SOD1 protein, interact with TRAF6 via the MATH domain. Here, we sought to uncover the structural details of this interaction through molecular dynamics (MD) simulations of a dimeric model system, coarse grained using the AWSEM force field. We used direct MD simulations to identify buried residues, and predict binding poses by clustering frames from the trajectories. Metadynamics simulations were also used to deduce preferred binding regions on the protein surfaces from the potential of the mean force in orientation space. Well-folded SOD1 was found to bind TRAF6 via co-option of its native homodimer interface. However, if loops IV and VII of SOD1 were disordered, as typically occurs in the absence of stabilizing Zn²⁺ ion binding, these disordered loops now participated in novel interactions with TRAF6. On TRAF6, multiple interaction hot-spots were distributed around the equatorial region of the MATH domain beta barrel. Expression of TRAF6 variants with mutations in this region in cultured cells demonstrated that TRAF6^T475 facilitates interaction with different SOD1 mutants. These findings contribute to our understanding of the disease mechanism and uncover potential targets for the development of therapeutics.

Neural relational inference to learn long-range allosteric interactions in proteins from molecular dynamics simulations

Protein allostery is a biological process facilitated by spatially long-range intra-protein communication, whereby ligand binding or amino acid change at a distant site affects the active site remotely. Molecular dynamics (MD) simulation provides a powerful computational approach to probe the allosteric effect. However, current MD simulations cannot reach the time scales of whole allosteric processes. The advent of deep learning made it possible to evaluate both spatially short and long-range communications for understanding allostery. For this purpose, we applied a neural relational inference model based on a graph neural network, which adopts an encoder-decoder architecture to simultaneously infer latent interactions for probing protein allosteric processes as dynamic networks of interacting residues. From the MD trajectories, this model successfully learned the long-range interactions and pathways that can mediate the allosteric communications between distant sites in the Pin1, SOD1, and MEK1 systems. Furthermore, the model can discover allostery-related interactions earlier in the MD simulation trajectories and predict relative free energy changes upon mutations more accurately than other methods.

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.

Determining Copper and Zinc Content in Superoxide Dismutase Using Electron Capture Dissociation Under Native Spray Conditions

Localizing metal binding to specific sites in proteins remains a challenging analytical problem in vitro and in vivo. Although metal binding can be maintained by "native" electrospray ionization with intact proteins for quantitation by mass spectrometry, subsequent fragmentation of proteins with slow-heating methods like collision-induced dissociation (CID) can scramble and detach metals. In contrast, electron capture dissociation (ECD) fragmentation produces highly localized bond cleavage that is well known to preserve posttranslational modifications. We show how a newly available ECD tool that can be retrofitted on standard QTOF mass spectrometers allows the sites of copper and zinc binding to be localized in the antioxidant enzyme Cu, Zn superoxide dismutase (SOD1). The loss of zinc from Cu, Zn SOD1 has been shown to induce motor neuron death and could have a causal role in the fatal neurodegenerative disease, amyotrophic lateral sclerosis (ALS). The methods described enable copper loss to be distinguished from zinc using distinct ECD fragments of SOD1 and are broadly applicable to other metalloproteins.

Hydrogen Peroxide and Amyotrophic Lateral Sclerosis: From Biochemistry to Pathophysiology

Free radicals are unstable chemical reactive species produced during Redox dyshomeostasis (RDH) inside living cells and are implicated in the pathogenesis of various neurodegenerative diseases. One of the most complicated and life-threatening motor neurodegenerative diseases (MND) is amyotrophic lateral sclerosis (ALS) because of the poor understanding of its pathophysiology and absence of an effective treatment for its cure. During the last 25 years, researchers around the globe have focused their interest on copper/zinc superoxide dismutase (Cu/Zn SOD, SOD1) protein after the landmark discovery of mutant SOD1 (mSOD1) gene as a risk factor for ALS. Substantial evidence suggests that toxic gain of function due to redox disturbance caused by reactive oxygen species (ROS) changes the biophysical properties of native SOD1 protein thus, instigating its fibrillization and misfolding. These abnormal misfolding aggregates or inclusions of SOD1 play a role in the pathogenesis of both forms of ALS, i.e., Sporadic ALS (sALS) and familial ALS (fALS). However, what leads to a decrease in the stability and misfolding of SOD1 is still in question and our scientific knowledge is scarce. A large number of studies have been conducted in this area to explore the biochemical mechanistic pathway of SOD1 aggregation. Several studies, over the past two decades, have shown that the SOD1-catalyzed biochemical reaction product hydrogen peroxide (H₂O₂) at a pathological concentration act as a substrate to trigger the misfolding trajectories and toxicity of SOD1 in the pathogenesis of ALS. These toxic aggregates of SOD1 also cause aberrant localization of TAR-DNA binding protein 43 (TDP-43), which is characteristic of neuronal cytoplasmic inclusions (NCI) found in ALS. Here in this review, we present the evidence implicating the pivotal role of H₂O₂ in modulating the toxicity of SOD1 in the pathophysiology of the incurable and highly complex disease ALS. Also, highlighting the role of H₂O₂ in ALS, we believe will encourage scientists to target pathological concentrations of H₂O₂ thereby halting the misfolding of SOD1.

Free energy calculations of ALS-causing SOD1 mutants reveal common perturbations to stability and dynamics along the maturation pathway

With over 150 heritable mutations identified as disease-causative, superoxide dismutase 1 (SOD1) has been a main target of amyotrophic lateral sclerosis (ALS) research and therapeutic efforts. However, recent evidence has suggested that neither loss of function nor protein aggregation is responsible for promoting neurotoxicity. Furthermore, there is no clear pattern to the nature or the location of these mutations that could suggest a molecular mechanism behind SOD1-linked ALS. Here, we utilize reliable and accurate computational techniques to predict the perturbations of 10 such mutations to the free energy changes of SOD1 as it matures from apo monomer to metallated dimer. We find that the free energy perturbations caused by these mutations strongly depend on maturational progress, indicating the need for state-specific therapeutic targeting. We also find that many mutations exhibit similar patterns of perturbation to native and non-native maturation, indicating strong thermodynamic coupling between the dynamics at various sites of maturation within SOD1. These results suggest the presence of an allosteric network in SOD1 which is vulnerable to disruption by these mutations. Analysis of these perturbations may contribute to uncovering a unifying molecular mechanism which explains SOD1-linked ALS and help to guide future therapeutic efforts.

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.

Utilizing Ion Mobility-Mass Spectrometry to Investigate the Unfolding Pathway of Cu/Zn Superoxide Dismutase

Native mass spectrometry has emerged as a powerful tool for structural biology as it enables the evaluation of molecules as they occur in their physiological conditions. Ion mobility spectrometry-mass spectrometry (IMS-MS) has shown essential in these analyses as it allows the measurement of the shape of a molecule, denoted as its collision cross section (CCS), and mass. The structural information garnered from native IMS-MS provides insight into the tertiary and quaternary structure of proteins and can be used to validate NMR or crystallographic X-ray structures. Additionally, due to the rapid nature (millisecond measurements) and ability of IMS-MS to analyze heterogeneous solutions, it can be used to address structural questions not possible with traditional structural approaches. Herein, we applied multiple solution conditions to systematically denature bovine Cu/Zn-superoxide dismutase (SOD1) and assess its unfolding pathway from the holo-dimer to the holo-monomer, single-metal monomer, and apo-monomer. Additionally, we compared and noted 1–2% agreement between CCS values from both drift tube IMS and trapped IMS for the SOD1 holo-monomer and holo-dimer. The observed CCS values were in excellent agreement with computational CCS values predicted from the homo-dimer crystal structure, showcasing the ability to use both IMS-MS platforms to provide valuable structural information for molecular modeling of protein interactions and structural assessments.

An Allosteric Pathway in Copper, Zinc Superoxide Dismutase Unravels the Molecular Mechanism of the G93A Amyotrophic Lateral Sclerosis-Linked Mutation

Several different mutations of the protein copper, zinc superoxide dismutase (SOD1) produce the neurodegenerative disorder amyotrophic lateral sclerosis (ALS). The molecular mechanism by which the diverse mutations converge to a similar pathology is currently unknown. The electrostatic loop (EL) of SOD1 is known to be affected in all of the studied ALS-linked mutations of SOD1. In this work, we employ a multiscale simulation approach to show that this perturbation corresponds to an increased probability of the EL detaching from its native position, exposing the metal site of the protein to water. From extensive atomistic and coarse-grained molecular dynamics (MD) simulations, we identify an allosteric pathway that explains the action of the distant G93A mutation on the EL. Finally, we employ quantum mechanics/molecular mechanics MD simulations to show that the opening of the EL decreases the Zn(II) affinity of the protein. As the loss of Zn(II) is at the center of several proposed pathogenic mechanisms in SOD1-linked ALS, the structural effect identified here not only is in agreement with the experimental data but also places the opening of the electrostatic loop as the possible main pathogenic effect for a significant number of ALS-linked SOD1 mutations.

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.

Folding of copper proteins: role of the metal?

Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.

Copper ion / H2O2 oxidation of Cu/Zn-Superoxide dismutase: Implications for enzymatic activity and antioxidant action

Copper ion-catalyzed oxidation of yeast SOD1 (ySOD1) was examined to determine early oxidative modifications, including oxidation of a crucial disulfide bond, and the structural and functional repercussions of these events. The study used distinct oxidative conditions: Cu²⁺/H₂O₂, Cu²⁺/H₂O₂/AscH⁻ and Cu²⁺/H₂O₂/glucose. Capillary electrophoresis experiments and quantification of protein carbonyls indicate that ySOD1 is highly susceptible to oxidative modification and that changes can be detected within 0.1 min of the initiation of the reaction. Oxidation-induced structural perturbations, characterized by circular dichroism, revealed the formation of partially-unfolded ySOD1 species in a dose-dependent manner. Consistent with these structural changes, pyrogallol assay indicates a partial loss of enzymatic activity. ESI-MS analyses showed seven distinct oxidized ySOD1 species under mild oxidation within 0.1 min. LC/MS analysis after proteolytic digestion demonstrated that the copper-coordinating active site histidine residues, His47 and His49, were converted into 2-oxo-histidine. Furthermore, the Cu and Zn bridging residue, His64 is converted into aspartate/asparagine. Importantly, the disulfide-bond Cys58-Cys147 which is critical for the structural and functional integrity of ySOD1 was detected as being oxidized at Cys147. We propose, based on LC/MS analyses, that disulfide-bond oxidation occurs without disulfide bond cleavage. Modifications were also detected at Met85 and five surface-exposed Lys residues. Based on these data we propose that the Cys58-Cys147 bond may act as a sacrificial target for oxidants and protect ySOD1 from oxidative inactivation arising from exposure to Cu²⁺/H₂O₂ and auto-inactivation during extended enzymatic turnover.

•

Oxidation of yeast superoxide dismutase (ySOD1) by Cu²⁺/H₂O₂ is examined.

•

Rapid modification of His, Met, Cys and Lys residues detected by LC-MS methods.

•

Oxidation of active site His residues and partial protein unfolding are early events.

•

The Cys58-Cys147 disulfide bond is oxidized and may act as a sacrificial target.

•

Excess exogenous Cu²⁺ decreases protein damage and can reverse loss of activity.

Neuroprotective effect of CuATSM on neurotoxin-induced motor neuron loss in an ALS mouse model

CuATSM is a PET-imaging agent that has recently received attention for its success in extending the lifespan in animals in several neurodegenerative disease models. In the SOD1^G93A model of ALS, CuATSM prolonged mouse longevity far longer than any previously tested therapeutic agents. The mechanism underlying this outcome has not been fully understood, but studies suggest that this copper complex contributes to maintaining copper homeostasis in mitochondria. More specifically for the SOD1 model, the molecule supplies copper back to the SOD1 protein. Additionally, CuATSM demonstrated similar protective effects in various in vivo Parkinson's disease mouse models. In the current pilot study, we utilized a neurodegenerative mouse model of motor neuron degeneration induced by the neurotoxin β-sitosterol β-D-glucoside. In this model, slow but distinct and progressive features of sporadic ALS occur. Treatment with CuATSM kept animal behavioural performance on par with the controls and prevented the extensive motor neuron degeneration and microglia activation seen in the untreated animals. These outcomes support a broader neuroprotective role for CuATSM beyond mutant SOD models of ALS.

A benchmark for the size of the QM system required for accurate hybrid QM/MM calculations on the metal site of the protein copper, zinc superoxide dismutase

The protein superoxide dismutase 1 (SOD1) is a copper and zinc-binding protein that has been implicated in the pathogenesis of amyotrophic lateral sclerosis (ALS). The Zn(II) binding to SOD1 is critical for the stability of the protein, and has been by itself implicated in ALS pathogenesis. Hence, the quantum mechanical (QM) study of the Zn(II)-site of SOD1 is relevant for understanding ALS. The hybrid QM-molecular mechanics (QM/MM) approach commonly employed for the QM study of proteins is highly dependent on the size of the sub-system treated quantum-mechanically. The size of the QM system also determines the computational feasibility of a given method. In the present work, we compare optimized geometries for the metal site and Zn(II) dissociation energies obtained with a QM/MM methodology employing different sizes for the QM sub-system. We find that geometries converge rapidly to RMSDs of around 0.3 Å, and fails to converge further, while a QM system of 480 atoms was required for converging the Zn(II) interaction energy of SOD1 to within 5 kcal*mol^-1, and a 611-atoms QM system for a 1 kcal*mol^-1 convergence with respect to our reference, 1280 QM-atoms system. Graphical Abstract The size of the QM system is critical for both the accuracy and the computational cost of a QM/MM calculation. We have identified a optimum balance for the study of the active site of the coppper, zinc superoxide dismutase.

TFE-induced local unfolding and fibrillation of SOD1: bridging the experiment and simulation studies

Misfolding and aggregation of Cu, Zn Superoxide dismutase (SOD1) is involved in the neurodegenerative disease, amyotrophic lateral sclerosis. Many studies have shown that metal-depleted, monomeric form of SOD1 displays substantial local unfolding dynamics and is the precursor for aggregation. Here, we have studied the structure and dynamics of different apo monomeric SOD1 variants associated with unfolding and aggregation in aqueous trifluoroethanol (TFE) through experiments and simulation. TFE induces partially unfolded β-sheet-rich extended conformations in these SOD1 variants, which subsequently develops aggregates with fibril-like characteristics. Fibrillation was achieved more easily in disulfide-reduced monomeric SOD1 when compared with wild-type and mutant monomeric SOD1. At higher concentrations of TFE, a native-like structure with the increase in α-helical content was observed. The molecular dynamics simulation results illustrate distinct structural dynamics for different regions of SOD1 variants and show uniform local unfolding of β-strands. The strands protected by the zinc-binding and electrostatic loops were found to unfold first in 20% (v/v) TFE, leading to a partial unfolding of β-strands 4, 5, and 6 which are prone to aggregation. Our results thus shed light on the role of local unfolding and conformational dynamics in SOD1 misfolding and aggregation.

Effects of maturation on the conformational free-energy landscape of SOD1

Amyotrophic lateral sclerosis (ALS) is a devastating fatal syndrome characterized by very rapid degeneration of motor neurons. A leading hypothesis is that ALS is caused by toxic protein misfolding and aggregation, as also occurs in many other neurodegenerative disorders, such as prion, Alzheimer's, Parkinson's, and Huntington's diseases. A prominent cause of familial ALS is mutations in the protein superoxide dismutase (SOD1), which promote the formation of misfolded SOD1 conformers that are prone to aberrant interactions both with each other and with other cellular components. We have shown previously that immature SOD1, lacking bound Cu and Zn metal ions and the intrasubunit disulfide bond (apoSOD1^2SH), has a rugged free-energy surface (FES) and exchanges with four other conformations (excited states) that have millisecond lifetimes and sparse populations on the order of a few percent. Here, we examine further states of SOD1 along its maturation pathway, as well as those off-pathway resulting from metal loss that have been observed in proteinaceous inclusions. Metallation and disulfide bond formation lead to structural transformations including local ordering of the electrostatic loop and native dimerization that are observed in rare conformers of apoSOD1^2SH; thus, SOD1 maturation may occur via a population-switch mechanism whereby posttranslational modifications select for preexisting structures on the FES. Metallation and oxidation of SOD1 stabilize the native, mature conformation and decrease the number of detected excited conformational states, suggesting that it is the immature forms of the protein that contribute to misfolded conformations in vivo rather than the highly stable enzymatically active dimer.

The Role of Metal Binding in the Amyotrophic Lateral Sclerosis-Related Aggregation of Copper-Zinc Superoxide Dismutase

Protein misfolding and conformational changes are common hallmarks in many neurodegenerative diseases involving formation and deposition of toxic protein aggregates. Although many players are involved in the in vivo protein aggregation, physiological factors such as labile metal ions within the cellular environment are likely to play a key role. In this review, we elucidate the role of metal binding in the aggregation process of copper-zinc superoxide dismutase (SOD1) associated to amyotrophic lateral sclerosis (ALS). SOD1 is an extremely stable Cu-Zn metalloprotein in which metal binding is crucial for folding, enzymatic activity and maintenance of the native conformation. Indeed, demetalation in SOD1 is known to induce misfolding and aggregation in physiological conditions in vitro suggesting that metal binding could play a key role in the pathological aggregation of SOD1. In addition, this study includes recent advances on the role of aberrant metal coordination in promoting SOD1 aggregation, highlighting the influence of metal ion homeostasis in pathologic aggregation processes.

Solvent sensitivity of protein aggregation in Cu, Zn superoxide dismutase: a molecular dynamics simulation study

Misfolding and aggregation of Cu, Zn Superoxide Dismutase (SOD1) is often found in amyotrophic lateral sclerosis (ALS) patients. The central apo SOD1 barrel was involved in protein maturation and pathological aggregation in ALS. In this work, we employed atomistic molecular dynamics (MD) simulations to study the conformational dynamics of SOD1^barrel monomer in different concentrations of trifluoroethanol (TFE). We find concentration dependence unusual structural and dynamical features, characterized by the local unfolding of SOD1^barrel. This partially unfolded structure is characterized by the exposure of hydrophobic core, is highly dynamic in nature, and is the precursor of aggregation seen in SOD1^barrel. Our computational studies supports the hypothesis of the formation of aggregation 'building blocks' by means of local unfolding of apo monomer as the mechanism of SOD1 fibrillar aggregation. The non-monotonic TFE concentration dependence of protein conformational changes was explored through simulation studies. Our results suggest that altered protein conformation and dynamics within its structure may underlie the aggregation of SOD1 in ALS.

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.

Copper delivery to the CNS by CuATSM effectively treats motor neuron disease in SOD(G93A) mice co-expressing the Copper-Chaperone-for-SOD

Over-expression of mutant copper, zinc superoxide dismutase (SOD) in mice induces ALS and has become the most widely used model of neurodegeneration. However, no pharmaceutical agent in 20 years has extended lifespan by more than a few weeks. The Copper-Chaperone-for-SOD (CCS) protein completes the maturation of SOD by inserting copper, but paradoxically human CCS causes mice co-expressing mutant SOD to die within two weeks of birth. Hypothesizing that co-expression of CCS created copper deficiency in spinal cord, we treated these pups with the PET-imaging agent CuATSM, which is known to deliver copper into the CNS within minutes. CuATSM prevented the early mortality of CCSxSOD mice, while markedly increasing Cu, Zn SOD protein in their ventral spinal cord. Remarkably, continued treatment with CuATSM extended the survival of these mice by an average of 18 months. When CuATSM treatment was stopped, these mice developed ALS-related symptoms and died within 3 months. Restoring CuATSM treatment could rescue these mice after they became symptomatic, providing a means to start and stop disease progression. All ALS patients also express human CCS, raising the hope that familial SOD ALS patients could respond to CuATSM treatment similarly to the CCSxSOD mice.

The Disulfide Bond, but Not Zinc or Dimerization, Controls Initiation and Seeded Growth in Amyotrophic Lateral Sclerosis-linked Cu,Zn Superoxide Dismutase (SOD1) Fibrillation

Aggregation of copper-zinc superoxide dismutase (SOD1) is a defining feature of familial ALS caused by inherited mutations in the sod1 gene, and misfolded and aggregated forms of wild-type SOD1 are found in both sporadic and familial ALS cases. Mature SOD1 owes its exceptional stability to a number of post-translational modifications as follows: formation of the intramolecular disulfide bond, binding of copper and zinc, and dimerization. Loss of stability due to the failure to acquire one or more of these modifications is proposed to lead to aggregation in vivo. Previously, we showed that the presence of apo-, disulfide-reduced SOD1, the most immature form of SOD1, results in initiation of fibrillation of more mature forms that have an intact Cys-57-Cys-146 disulfide bond and are partially metallated. In this study, we examine the ability of each of the above post-translational modifications to modulate fibril initiation and seeded growth. Cobalt or zinc binding, despite conferring great structural stability, neither inhibits the initiation propensity of disulfide-reduced SOD1 nor consistently protects disulfide-oxidized SOD1 from being recruited into growing fibrils across wild-type and a number of ALS mutants. In contrast, reduction of the disulfide bond, known to be necessary for fibril initiation, also allows for faster recruitment during seeded amyloid growth. These results identify separate factors that differently influence seeded growth and initiation and indicate a lack of correlation between the overall thermodynamic stability of partially mature SOD1 states and their ability to initiate fibrillation or be recruited by a growing fibril.

Aberrant zinc binding to immature conformers of metal-free copper-zinc superoxide dismutase triggers amorphous aggregation

Superoxide dismutase 1 (SOD1) is a Cu/Zn metalloenzyme that aggregates in amyotrophic lateral sclerosis (ALS), a fatal neurodegenerative disorder. Correct metal insertion during SOD1 biosynthesis is critical to prevent misfolding; however Zn(2+) can bind to the copper-site leading to an aberrantly metallated protein. These effects of Zn(2+) misligation on SOD1 aggregation remain to be explored, even though Zn(2+) levels are upregulated in ALS motor neurons. Here we use complementary biophysical methods to investigate Zn(2+) binding and its effects on the aggregation of three immature metal-free SOD1 conformers that represent biogenesis intermediates: dimeric, monomeric and reduced monomeric SOD1. Using isothermal titration calorimetry we determined that Zn(2+) binds to all conformers both at the zinc- as well as to the copper-site; however Zn(2+) binding mechanisms to the zinc-site have distinct characteristics across immature conformers. We show that this 'zinc overload' of immature SOD1 promotes intermolecular interactions, as evidenced by dynamic light scattering and ThT fluorescence kinetic studies. Analysis of aged zinc-induced aggregates by energy-dispersive X-ray and electron energy-loss spectroscopy shows that aggregates integrate some Zn(2+). In addition, electron diffraction analysis identifies nano-scaled crystalline materials and amyloid fibril-like reflections. Transmission electron microscopy reveals that Zn(2+) diverts the SOD1 aggregation pathway from fibrils to amorphous aggregate, and electrophoretic analysis evidences an increase in insoluble materials. Overall, we provide evidence that aberrant zinc coordination to immature conformers broadens the population of SOD1 misfolded species at early aggregation stages and provide evidence for a high structural polymorphism and heterogeneity of SOD1 aggregates.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

Diffuse binding of Zn(2+) to the denatured ensemble of Cu/Zn superoxide dismutase 1

•

Nearly complete backbone assignments for monomeric SOD1 were obtained.

•

Monomeric SOD1 is unstructured in 9 M urea.

•

Zn²⁺-binding to monomeric SOD1 occurs through diffuse coordination to most His residues.

•

The binding geometry of Zn²⁺ is different than in the native, folded SOD1.

The stability and structural properties of the metalloprotein superoxide dismutase 1 (SOD1) are found to depend critically on metal ions. Native SOD1 monomers coordinate one structural Zn²⁺ and one redox-active Cu^2+/1+ to the active site. To do this, the Zn²⁺ ions need to interact with the SOD1 protein on the denatured side of the folding barrier, prior to the formation of the folding nucleus. In this study, we have examined at residue level the nature of this early Zn²⁺ binding by NMR studies on the urea denatured-state of SOD1. Nearly complete backbone chemical shift assignments were obtained in 9 M urea at physiological pH, conditions at which NMR studies are scarce. Our results demonstrate that SOD1 is predominantly unstructured under these conditions. Chemical-shift changes upon Zn²⁺ titration show that denatured SOD1 retains a significant affinity to Zn²⁺ ions, even in 9 M urea. However, the Zn²⁺ interactions are not limited to the native metal-binding ligands in the two binding sites, but are seen for all His residues. Moreover, the native Cu^2+/1+ ligand H46 seems not to bind as well as the other His residues, while the nearby non-native H43 does bind, indicating that the binding geometry is relaxed. The result suggests that the Zn²⁺-binding observed to catalyze folding of SOD1 in physiological buffer is initiated by diffuse, non-specific coordination to the coil, which subsequently funnels by ligand exchange into the native coordination geometry of the folded monomer. Altogether, this diffuse binding is a result with fundamental implications for folding of metalloproteins in general.

In-cell NMR reveals potential precursor of toxic species from SOD1 fALS mutants

Mutations in the superoxide dismutase 1 (SOD1) gene are related to familial cases of amyotrophic lateral sclerosis (fALS). Here we exploit in-cell NMR to characterize the protein folding and maturation of a series of fALS-linked SOD1 mutants in human cells and to obtain insight into their behaviour in the cellular context, at the molecular level. The effect of various mutations on SOD1 maturation are investigated by changing the availability of metal ions in the cells, and by coexpressing the copper chaperone for SOD1, hCCS. We observe for most of the mutants the occurrence of an unstructured SOD1 species, unable to bind zinc. This species may be a common precursor of potentially toxic oligomeric species, that are associated with fALS. Coexpression of hCCS in the presence of copper restores the correct maturation of the SOD1 mutants and prevents the formation of the unstructured species, confirming that hCCS also acts as a molecular chaperone.

Insights into SOD1-linked amyotrophic lateral sclerosis from NMR studies of Ni(2+)- and other metal-ion-substituted wild-type copper-zinc superoxide dismutases

The dimeric Cu-Zn superoxide dismutase (SOD1) is a particularly interesting system for biological inorganic chemical studies because substitutions of the native Cu and/or Zn ions by a nonnative metal ion cause minimal structural changes and result in high enzymatic activity for those derivatives with Cu remaining in the Cu site. The pioneering NMR studies of the magnetically coupled derivative Cu2Co2SOD1 by Ivano Bertini and coworkers are of particular importance in this regard. In addition to Co(2+), Ni(2+) is a versatile metal ion for substitution into SOD1, showing very little disturbance of the structure in Cu2Ni2SOD1 and acting as a very good mimic of the native Cu ion in Ni2Zn2SOD1. The NMR studies presented here were inspired by and are indebted to Ivano Bertini's paramagnetic NMR pursuits of metalloproteins. We report Ni(2+) binding to apo wild-type SOD1 and a time-dependent Ni(2+) migration from the Zn site to the Cu site, and the preparation and characterization of Ni2Ni2SOD1, which shows coordination properties similar to those of Cu2Cu2SOD1, namely, an anion-binding property different from that of the wild type and a possibly broken bridging His. Mutations in the human SOD1 gene can cause familial amyotrophic lateral sclerosis (ALS), and mutant SOD1 proteins with significantly altered metal-binding behaviors are implicated in causing the disease. We conclude by discussing the effects of the ALS mutations on the remarkable stabilities and metal-binding properties of wild-type SOD1 proteins and the implications concerning the causes of SOD1-linked ALS.

An emerging role for misfolded wild-type SOD1 in sporadic ALS pathogenesis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder that targets motor neurons, leading to paralysis and death within a few years of disease onset. While several genes have been linked to the inheritable, or familial, form of ALS, much less is known about the cause(s) of sporadic ALS, which accounts for ~90% of ALS cases. Due to the clinical similarities between familial and sporadic ALS, it is plausible that both forms of the disease converge on a common pathway and, therefore, involve common factors. Recent evidence suggests the Cu,Zn-superoxide dismutase (SOD1) protein to be one such factor that is common to both sporadic and familial ALS. In 1993, mutations were uncovered in SOD1 that represent the first known genetic cause of familial ALS. While the exact mechanism of mutant-SOD1 toxicity is still not known today, most evidence points to a gain of toxic function that stems, at least in part, from the propensity of this protein to misfold. In the wild-type SOD1 protein, non-genetic perturbations such as metal depletion, disruption of the quaternary structure, and oxidation, can also induce SOD1 to misfold. In fact, these aforementioned post-translational modifications cause wild-type SOD1 to adopt a “toxic conformation” that is similar to familial ALS-linked SOD1 variants. These observations, together with the detection of misfolded wild-type SOD1 within human post-mortem sporadic ALS samples, have been used to support the controversial hypothesis that misfolded forms of wild-type SOD1 contribute to sporadic ALS pathogenesis. In this review, we present data from the literature that both support and contradict this hypothesis. We also discuss SOD1 as a potential therapeutic target for both familial and sporadic ALS.

The role of solvent exclusion in the interaction between D124 and the metal site in SOD1: implications for ALS

Structural changes in the metal site of the copper-zinc superoxide dismutase (SOD1) are involved in the various mechanisms proposed for the pathogenesis of the SOD1-linked familial form of amyotrophic lateral sclerosis (ALS). Elucidating how the metal site of SOD1 can be disrupted by ALS-linked mutations is important for a better understanding of the pathogenesis of the disease and for developing more efficient treatments. Residue D124, a second-sphere ligand of the copper and zinc ions, is known from experimental studies to be essential for the integrity of the metal-site structure. In this work, we used density functional theory calculations and molecular dynamics simulations to elucidate which factors keep D124 attached to the metal site and how structural changes may disrupt the binding between D124 and the metal first-sphere ligands. The calculations show that D124 is kept attached to the metal site in a kinetic trap. The exclusion of solvent molecules by the electrostatic loop of the protein is found to create the binding of D124 to the metal site. The calculations also indicate that changes in the structure of the electrostatic loop of the protein can weaken the D124-metal site interaction, lowering the affinity of the zinc site for the metal. Destabilization of the electrostatic loop of SOD1 has been previously shown to be a common property of ALS-linked variants of the protein, but its role in the pathogenesis of SOD1-linked ALS has not been elucidated.

Early steps in thermal unfolding of superoxide dismutase 1 are similar to the conformational changes associated with the ALS-associated A4V mutation

There are over 100 mutations in Cu/Zn superoxide dismutase (SOD1) that result in a subset of familial amyotrophic lateral sclerosis (fALS) cases. The hypothesis that dissociation of the dimer, misfolding of the monomer and subsequent aggregation of mutant SOD1 leads to fALS has been gaining support as an explanation for how these disparate missense mutations cause the same disease. These forms are only responsible for a fraction of the ALS cases; however, the rest are sporadic. Starting with a folded apo monomer, the species considered most likely to be involved in misfolding, we used high-temperature all-atom molecular dynamics simulations to explore the events of the wild-type protein unfolding through the denatured state. All simulations showed early loss of structure along the β5-β6 edge of the β-sandwich, supporting earlier findings of instability in this region. Transition state structures identified from the simulations are in good agreement with experiment, providing detailed, validated molecular models for this elusive state. Furthermore, we compare the process of thermal unfolding investigated here to that of the lethal A4V mutant-induced unfolding at physiological temperature and find that the pathways are very similar.

Identification of human monoclonal antibodies specific for human SOD1 recognizing distinct epitopes and forms of SOD1

Mutations in the gene encoding human SOD1 (hSOD1) can cause amyotrophic lateral sclerosis (ALS) yet the mechanism by which mutant SOD1 can induce ALS is not fully understood. There is currently no cure for ALS or treatment that significantly reduces symptoms or progression. To develop tools to understand the protein conformations present in mutant SOD1-induced ALS and as possible immunotherapy, we isolated and characterized eleven unique human monoclonal antibodies specific for hSOD1. Among these, five recognized distinct linear epitopes on hSOD1 that were not available in the properly-folded protein but were available on forms of protein with some degree of misfolding. The other six antibodies recognized conformation-dependent epitopes that were present in the properly-folded protein with two different recognition profiles: three could bind hSOD1 dimer or monomer and the other three were specific for hSOD1 dimer only. Antibodies with the capacity to bind hSOD1 monomer were able to prevent increased hydrophobicity when mutant hSOD1 was exposed to increased temperature and EDTA, suggesting that the antibodies stabilized the native structure of hSOD1. Two antibodies were tested in a G93A mutant hSOD1 transgenic mouse model of ALS but did not yield a statistically significant increase in overall survival. It may be that the two antibodies selected for testing in the mouse model were not effective for therapy or that the model and/or route of administration were not optimal to produce a therapeutic effect. Therefore, additional testing will be required to determine therapeutic potential for SOD1 mutant ALS and potentially some subset of sporadic ALS.

Polymer uncrossing and knotting in protein folding, and their role in minimal folding pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold [Formula: see text], [Formula: see text], [Formula: see text], and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation "alignment". The consensus minimal pathway is constructed and shown schematically for representative cases of an [Formula: see text], [Formula: see text], and knotted protein. An overlap parameter is defined between pathways; we find that [Formula: see text] proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and [Formula: see text] proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.

Ligand binding and aggregation of pathogenic SOD1

Mutations in the gene encoding Cu/Zn superoxide dismutase-1 cause amyotrophic lateral sclerosis. Superoxide dismutase-1 mutations decrease protein stability and promote aggregation. The mutant monomer is thought to be an intermediate in the pathway from the superoxide dismutase-1 dimer to aggregate. Here we find that the monomeric copper-apo, zinc-holo protein is structurally perturbed and the apo-protein aggregates without reattainment of the monomer-dimer equilibrium. Intervention to stabilize the superoxide dismutase-1 dimer and inhibit aggregation is regarded as a potential therapeutic strategy. We describe protein-ligand interactions for two compounds, Isoproterenol and 5-fluorouridine, highlighted as superoxide dismutase-1 stabilizers. We find both compounds interact with superoxide dismutase-1 at a key region identified at the core of the superoxide dismutase-1 fibrillar aggregates, β-barrel loop II-strand 3, rather than the proposed dimer interface site. This illustrates the need for direct structural observations when developing compounds for protein-targeted therapeutics.

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

Composition of soluble misfolded superoxide dismutase-1 in murine models of amyotrophic lateral sclerosis

A common cause of amyotrophic lateral sclerosis is mutations in superoxide dismutase-1, which provoke the disease by an unknown mechanism. We have previously found that soluble hydrophobic misfolded mutant human superoxide dismutase-1 species are enriched in the vulnerable spinal cords of transgenic model mice. The levels were broadly inversely correlated with life spans, suggesting involvement in the pathogenesis. Here, we used methods based on antihuman superoxide dismutase-1 peptide antibodies specific for misfolded species to explore the composition and amounts of soluble misfolded human superoxide dismutase-1 in tissue extracts. Mice expressing 5 different human superoxide dismutase-1 variants with widely variable structural characteristics were examined. The levels were generally higher in spinal cords than in other tissues. The major portion of misfolded superoxide dismutase-1 was shown to be monomers lacking the C57-C146 disulfide bond with large hydrodynamic volume, indicating a severely disordered structure. The remainder of the misfolded protein appeared to be non-covalently associated in 130- and 250-kDa complexes. The malleable monomers should be prone to aggregate and associate with other cellular components, and should be easily translocated between compartments. They may be the primary cause of toxicity in superoxide dismutase-1-induced amyotrophic lateral sclerosis.

Enthalpic barriers dominate the folding and unfolding of the human Cu, Zn superoxide dismutase monomer

The rate-limiting step in the formation of the native dimeric state of human Cu, Zn superoxide dismutase (SOD1) is a very slow monomer folding reaction that governs the lifetime of its unfolded state. Mutations at dozens of sites in SOD1 are known to cause a fatal motor neuron disease, amyotrophic lateral sclerosis, and recent experiments implicate the unfolded state as a source of soluble oligomers and histologically observable aggregates thought to be responsible for toxicity. To determine the thermodynamic properties of the transition state ensemble (TSE) limiting the folding of this high-contact-order β-sandwich motif, we performed a combined thermal/urea denaturation thermodynamic/kinetic analysis. The barriers to folding and unfolding are dominated by the activation enthalpy at 298 K and neutral pH; the activation entropy is favorable and reduces the barrier height for both reactions. The absence of secondary structure formation or large-scale chain collapse prior to crossing the barrier for folding led to the conclusion that dehydration of nonpolar surfaces in the TSE is responsible for the large and positive activation enthalpy. Although the activation entropy favors the folding reaction, the transition from the unfolded state to the native state is entropically disfavored at 298 K. The opposing entropic contributions to the free energies of the TSE and the native state during folding provide insights into structural properties of the TSE. The results also imply a crucial role for water in governing the productive folding reaction and enhancing the propensity for the aggregation of SOD1.