Cu/Zn Superoxide Dismutase Forms Amyloid Fibrils under Near-Physiological Quiescent Conditions: The Roles of Disulfide Bonds and Effects of Denaturant

Unveiling local and global conformational changes and allosteric communications in SOD1 systems using molecular dynamics simulation and network analyses

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a serious neurodegenerative disorder affecting nerve cells in the brain and spinal cord that is caused by mutations in the superoxide dismutase 1 (SOD1) enzyme. ALS-related mutations cause misfolding, dimerisation instability, and increased formation of aggregates. The underlying allosteric mechanisms, however, remain obscure as far as details of their fundamental atomistic structure are concerned. Hence, this gap in knowledge limits the development of novel SOD1 inhibitors and the understanding of how disease-associated mutations in distal sites affect enzyme activity.

METHODS

We combined microsecond-scale based unbiased molecular dynamics (MD) simulation with network analysis to elucidate the local and global conformational changes and allosteric communications in SOD1 Apo (unmetallated form), Holo, Apo_CallA (mutant and unmetallated form), and Holo_CallA (mutant form) systems. To identify hotspot residues involved in SOD1 signalling and allosteric communications, we performed network centrality, community network, and path analyses.

RESULTS

Structural analyses showed that unmetallated SOD1 systems and cysteine mutations displayed large structural variations in the catalytic sites, affecting structural stability. Inter- and intra H-bond analyses identified several important residues crucial for maintaining interfacial stability, structural stability, and enzyme catalysis. Dynamic motion analysis demonstrated more balanced atomic displacement and highly correlated motions in the Holo system. The rationale for structural disparity observed in the disulfide bond formation and R143 configuration in Apo and Holo systems were elucidated using distance and dihedral probability distribution analyses.

CONCLUSION

Our study highlights the efficiency of combining extensive MD simulations with network analyses to unravel the features of protein allostery.

Role of charged residues of the "electrostatic loop" of hSOD1 in promotion of aggregation: Implications for the mechanism of ALS-associated mutations under amyloidogenic conditions

Protein misfolding and amyloid formation are hallmarks of numerous diseases, including amyotrophic lateral sclerosis (ALS), in which hSOD1 aggregation is involved in pathogenesis. We used two point mutations in the electrostatic loop, G138E and T137R, to analyze charge distribution under destabilizing circumstances to gain more about how ALS-linked mutations affect SOD1 protein stability or net repulsive charge. We show that protein charge is important in the ALS disease process using bioinformatics and experiments. The MD simulation findings demonstrate that the mutant protein differs significantly from WT SOD1, which is consistent with the experimental evidence. The specific activity of the wild type was 1.61 and 1.48 times higher than that of the G138E and T137R mutants, respectively. Under amyloid induction conditions, the intensity of intrinsic and ANS fluorescence in both mutants reduced. Increasing the content of β-sheet structures in mutants can be attributed to aggregation propensity, which was confirmed using CD polarimetry and FTIR spectroscopy. Our findings show that two ALS-related mutations promote the formation of amyloid-like aggregates at near physiological pH under destabilizing conditions, which were detected using spectroscopic probes such as Congo red and ThT fluorescence, and also further confirmation of amyloid-like species by TEM. Overall, our results provide evidence supporting the notion that negative charge changes combined with other destabilizing factors play an important role in increasing protein aggregation by reducing repulsive negative charges.

The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins

Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.

Tailoring the positive and negative solvatochromism for chalcone analogues to detect heterozygous protein co-aggregation

It is rare for one fluorophore scaffold to harbor both positive and negative solvatochromism. Herein, we tailor chalcone analogues to achieve both positive- and negative-polarity sensitivity of fluorescence intensity. We explore two chalcones of opposite solvatochromism to simultaneously detect the co-aggregation of wild-type and mutant superoxide dismutase that cause amyotrophic lateral sclerosis disease.

Identification of Distinct Soluble States During Fibril Formation Using Multilinear Analysis of NMR Diffusion Data

Protein misfolding and self-assembling into amyloid structures are associated with a number of diseases. Characterization of protein amyloid formation reactions is a challenging task as transient populations of multiple species are involved. Here we outline a method for identification and characterization of the individual soluble states during protein amyloid formation. The method combines NMR translational diffusion measurements with multilinear data analysis.

Influence of denaturants on amyloid β42 aggregation kinetics

Amyloid formation is linked to devastating neurodegenerative diseases, motivating detailed studies of the mechanisms of amyloid formation. For Aβ, the peptide associated with Alzheimer's disease, the mechanism and rate of aggregation have been established for a range of variants and conditions in vitro and in bodily fluids. A key outstanding question is how the relative stabilities of monomers, fibrils and intermediates affect each step in the fibril formation process. By monitoring the kinetics of aggregation of Aβ42, in the presence of urea or guanidinium hydrochloride (GuHCl), we here determine the rates of the underlying microscopic steps and establish the importance of changes in relative stability induced by the presence of denaturant for each individual step. Denaturants shift the equilibrium towards the unfolded state of each species. We find that a non-ionic denaturant, urea, reduces the overall aggregation rate, and that the effect on nucleation is stronger than the effect on elongation. Urea reduces the rate of secondary nucleation by decreasing the coverage of fibril surfaces and the rate of nucleus formation. It also reduces the rate of primary nucleation, increasing its reaction order. The ionic denaturant, GuHCl, accelerates the aggregation at low denaturant concentrations and decelerates the aggregation at high denaturant concentrations. Below approximately 0.25 M GuHCl, the screening of repulsive electrostatic interactions between peptides by the charged denaturant dominates, leading to an increased aggregation rate. At higher GuHCl concentrations, the electrostatic repulsion is completely screened, and the denaturing effect dominates. The results illustrate how the differential effects of denaturants on stability of monomer, oligomer and fibril translate to differential effects on microscopic steps, with the rate of nucleation being most strongly reduced.

Effects of Protein Crowders and Charge on the Folding of Superoxide Dismutase 1 Variants: A Computational Study

The neurodegenerative disease amyotrophic lateral sclerosis (ALS) is associated with the misfolding and aggregation of the metalloenzyme protein superoxide dismutase 1 (SOD1) via mutations that destabilize the monomer-dimer interface. In a cellular environment, crowding and electrostatic screening play essential roles in the folding and aggregation of the SOD1 monomers. Despite numerous studies on the effects of mutations on SOD1 folding, a clear understanding of the interplay between crowding, folding, and aggregation in vivo remains lacking. Using a structure-based minimal model for molecular dynamics simulations, we investigate the role of self-crowding and charge on the folding stability of SOD1 and the G41D mutant where experimentalists were intrigued by an alteration of the folding mechanism by a single point mutation from glycine to charged aspartic acid. We show that unfolded SOD1 configurations are significantly affected by charge and crowding, a finding that would be extremely costly to achieve with all-atom simulations, while the native state is not significantly altered. The mutation at residue 41 alters the interactions between proteins in the unfolded states instead of those within a protein. This paper suggests electrostatics may play an important role in the folding pathway of SOD1 and modifying the charge via mutation and ion concentration may change the dominant interactions between proteins, with potential impacts for aggregation of the mutants. This work provides a plausible reason for the alteration of the unfolded states to address why the mutant G41D causes the changes to the folding mechanism of SOD1 that have intrigued experimentalists.

Cryo-EM structure of an amyloid fibril formed by full-length human SOD1 reveals its conformational conversion

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease. Misfolded Cu, Zn-superoxide dismutase (SOD1) has been linked to both familial and sporadic ALS. SOD1 fibrils formed in vitro share toxic properties with ALS inclusions. Here we produced cytotoxic amyloid fibrils from full-length apo human SOD1 under reducing conditions and determined the atomic structure using cryo-EM. The SOD1 fibril consists of a single protofilament with a left-handed helix. The fibril core exhibits a serpentine fold comprising N-terminal segment (residues 3–55) and C-terminal segment (residues 86–153) with an intrinsic disordered segment. The two segments are zipped up by three salt bridge pairs. By comparison with the structure of apo SOD1 dimer, we propose that eight β-strands (to form a β-barrel) and one α-helix in the subunit of apo SOD1 convert into thirteen β-strands stabilized by five hydrophobic cavities in the SOD1 fibril. Our data provide insights into how SOD1 converts between structurally and functionally distinct states.

Misfolded SOD1 has been linked to both familial and sporadic ALS. Here the authors have determined the cryo-EM structure of SOD1 fibrils, providing insights into the conversion of SOD1 from its immature form into an aggregated form during pathogenesis of ALS.

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.

SOD1 in ALS: Taking Stock in Pathogenic Mechanisms and the Role of Glial and Muscle Cells

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the loss of motor neurons in the brain and spinal cord. While the exact causes of ALS are still unclear, the discovery that familial cases of ALS are related to mutations in the Cu/Zn superoxide dismutase (SOD1), a key antioxidant enzyme protecting cells from the deleterious effects of superoxide radicals, suggested that alterations in SOD1 functionality and/or aberrant SOD1 aggregation strongly contribute to ALS pathogenesis. A new scenario was opened in which, thanks to the generation of SOD1 related models, different mechanisms crucial for ALS progression were identified. These include excitotoxicity, oxidative stress, mitochondrial dysfunctions, and non-cell autonomous toxicity, also implicating altered Ca2+ metabolism. While most of the literature considers motor neurons as primary target of SOD1-mediated effects, here we mainly discuss the effects of SOD1 mutations in non-neuronal cells, such as glial and skeletal muscle cells, in ALS. Attention is given to the altered redox balance and Ca2+ homeostasis, two processes that are strictly related with each other. We also provide original data obtained in primary myocytes derived from hSOD1(G93A) transgenic mice, showing perturbed expression of Ca2+ transporters that may be responsible for altered mitochondrial Ca2+ fluxes. ALS-related SOD1 mutants are also responsible for early alterations of fundamental biological processes in skeletal myocytes that may impinge on skeletal muscle functions and the cross-talk between muscle cells and motor neurons during disease progression.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

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.

The Unfolding Journey of Superoxide Dismutase 1 Barrels under Crowding: Atomistic Simulations Shed Light on Intermediate States and Their Interactions with Crowders

The thermal stability of the superoxide dismutase 1 protein in a crowded solution is investigated by performing enhanced sampling molecular simulations. By complementing thermal unfolding experiments done close to physiological conditions (200 mg/mL), we provide evidence that the presence of the protein crowder bovine serum albumin in different packing states has only a minor, and essentially destabilizing, effect. The finding that quinary interactions counteract the pure stabilization contribution stemming from excluded volume is rationalized here by exploring the SOD1 unfolding mechanism in microscopic detail. In agreement with recent experiments, we unveil the importance of intermediate unfolded states as well as the correlation between protein conformations and local packing with the crowders. This link helps us to elucidate why certain SOD1 mutations involved in the ALS disease reverse the stability effect of the intracellular environment.

A Systematic and Comprehensive Review on Disease-Causing Genes in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder and is characterized by degeneration and axon loss from the upper motor neuron, that descends from the lower motor neuron in the brain. Over the period, assorted outcomes from medical findings, molecular pathogenesis, and structural and biophysical studies have abetted in providing thoughtful insights underlying the importance of disease-causing genes in ALS. Consequently, numerous mechanisms were proposed for the pathogenesis of ALS, considering protein mutations, aggregation, and misfolding. Besides, the answers to the majority of ALS cases that happen to be sporadic still remain obscure. The application in discovering susceptibility factors in ALS contemplating the genetic factors is to be further dissevered in the future years with innovation in research studies. Hence, this review targets in revisiting the breakthroughs on the disease-causing genes related with ALS.

Revealing Well-Defined Soluble States during Amyloid Fibril Formation by Multilinear Analysis of NMR Diffusion Data

Amyloid fibril formation is a hallmark of neurodegenerative disease caused by protein aggregation. Oligomeric protein states that arise during the process of fibril formation often coexist with mature fibrils and are known to cause cell death in disease model systems. Progress in this field depends critically on development of analytical methods that can provide information about the mechanisms and species involved in oligomerization and fibril formation. Here, we demonstrate how the powerful combination of diffusion NMR and multilinear data analysis can efficiently disentangle the number of involved species, their kinetic rates of formation or disappearance, spectral contributions, and diffusion coefficients, even without prior knowledge of the time evolution of the process or chemical shift assignments of the various species. Using this method we identify oligomeric species that form transiently during aggregation of human superoxide dismutase 1 (SOD1), which is known to form misfolded aggregates in patients with amyotrophic lateral sclerosis. Specifically, over a time course of 42 days, during which SOD1 fibrils form, we detect the disappearance of the native monomeric species, formation of a partially unfolded intermediate in the dimer to tetramer size range, subsequent formation of a distinct similarly sized species that dominates the final spectrum detected by solution NMR, and concomitant appearance of small peptide fragments.

Network mapping of the conformational heterogeneity of SOD1 by deploying statistical cluster analysis of FTIR spectra

A crucial contribution to the heterogeneity of the conformational landscape of a protein comes from the way an intermediate relates to another intermediate state in its journey from the unfolded to folded or misfolded form. Unfortunately, it is extremely hard to decode this relatedness in a quantifiable manner. Here, we developed an application of statistical cluster analyses to explore the conformational heterogeneity of a metalloenzyme, human cytosolic copper-zinc superoxide dismutase (SOD1), using the inputs from infrared spectroscopy. This study provides a quantifiable picture of how conformational information at one particular site (for example, the copper-binding pocket) is related to the information at the second site (for example, the zinc-binding pocket), and how this relatedness is transferred to the global conformational information of the protein. The distance outputs were used to quantitatively generate a network capturing the folding sub-stages of SOD1.

Stability and Local Unfolding of SOD1 in the Presence of Protein Crowders

Using NMR and Monte Carlo (MC) methods, we investigate the stability and dynamics of superoxide dismutase 1 (SOD1) in homogeneous crowding environments, where either bovine pancreatic trypsin inhibitor (BPTI) or the B1 domain of streptococcal protein G (PGB1) serves as a crowding agent. By NMR, we show that both crowders, and especially BPTI, cause a drastic loss in the overall stability of SOD1 in its apo monomeric form. Additionally, we determine chemical shift perturbations indicating that SOD1 interacts with the crowder proteins in a residue-specific manner that further depends on the identity of the crowding protein. Furthermore, the specificity of SOD1-crowder interactions is reciprocal: chemical shift perturbations on BPTI and PGB1 identify regions that interact preferentially with SOD1. By MC simulations, we investigate the local unfolding of SOD1 in the absence and presence of the crowders. We find that the crowders primarily interact with the long flexible loops of the folded SOD1 monomer. The basic mechanisms by which the SOD1 β-barrel core unfolds remain unchanged when adding the crowders. In particular, both with and without the crowders, the second β-sheet of the barrel is more dynamic and unfolding-prone than the first. Notably, the MC simulations (exploring the early stages of SOD1 unfolding) and the NMR experiments (under equilibrium conditions) identify largely the same set of PGB1 and BPTI residues as prone to form SOD1 contacts. Thus, contacts stabilizing the unfolded state of SOD1 in many cases appear to form early in the unfolding reaction.

Tryptophan 32 mediates SOD1 toxicity in a in vivo motor neuron model of ALS and is a promising target for small molecule therapeutics

SOD1 misfolding, toxic gain of function, and spread are proposed as a pathological basis of amyotrophic lateral sclerosis (ALS), but the nature of SOD1 toxicity has been difficult to elucidate. Uniquely in SOD1 proteins from humans and other primates, and rarely in other species, a tryptophan residue at position 32 (W32) is predicted to be solvent exposed and to participate in SOD1 misfolding. We hypothesized that W32 is influential in SOD1 acquiring toxicity, as it is known to be important in template-directed misfolding. We tested if W32 contributes to SOD1 cytotoxicity and if it is an appropriate drug target to ameliorate ALS-like neuromuscular deficits in a zebrafish model of motor neuron axon morphology and function (swimming). Embryos injected with human SOD1 variant with W32 substituted for a serine (SOD1^W32S) had reduced motor neuron axonopathy and motor deficits compared to those injected with wildtype or disease-associated SOD1. A library of FDA-approved small molecules was ranked with virtual screening based on predicted binding to W32, and subsequently filtered for analogues using a pharmacophore model based on molecular features of the uracil moiety of a small molecule previously predicted to interact with W32 (5'-fluorouridine or 5'-FUrd). Along with testing 5'-FUrd and uridine, a lead candidate from this list was selected based on its lower toxicity and improved blood brain barrier penetrance; telbivudine significantly rescued SOD1 toxicity in a dose-dependent manner. The mechanisms whereby the small molecules ameliorated motor neuron phenotypes were specifically mediated through human SOD1 and its residue W32, because these therapeutics had no measurable impact on the effects of UBQLN4^D90A, EtOH, or tryptophan-deficient human SOD1^W32S. By substituting W32 for a more evolutionarily conserved residue (serine), we confirmed the significant influence of W32 on human SOD1 toxicity to motor neuron morphology and function; further, we performed pharmaceutical targeting of the W32 residue for rescuing SOD1 toxicity. This unique residue offers future novel insights into SOD1 stability and toxic gain of function, and therefore poses an potential target for drug therapy.

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.