Mechanical Probes of SOD1 Predict Systematic Trends in Metal and Dimer Affinity of ALS-Associated Mutants

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.

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.

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.

A pH Switch Controls Zinc Binding in Tomato Copper-Zinc Superoxide Dismutase

Copper-zinc superoxide dismutase (SOD1) is a major antioxidant metalloenzyme that protects cells from oxidative damage by superoxide anions (O2-). Structural, biophysical, and other characteristics have in the past been compiled for mammalian SOD1s and for the highly homologous fungal and bovine SOD1s. Here, we characterize the biophysical properties of a plant SOD1 from tomato chloroplasts and present several of its crystal structures. The most unusual of these structures is a structure at low pH in which tSOD1 harbors zinc in the copper-binding site but contains no metal in the zinc-binding site. The side chain of D83, normally a zinc ligand, adopts an alternate rotameric conformation to form an unusual bidentate hydrogen bond with the side chain of D124, precluding metal binding in the zinc-binding site. This alternate conformation of D83 appears to be responsible for the previously observed pH-dependent loss of zinc from the zinc-binding site of SOD1. Titrations of cobalt into apo tSOD1 at a similar pH support the lack of an intact zinc-binding site. Further characterization of tSOD1 reveals that it is a weaker dimer relative to human SOD1 and that it can be activated in vivo through a copper chaperone for the SOD1-independent mechanism.

Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.

When a foreign gene meets its native counterpart: computational biophysics analysis of two PgiC loci in the grass Festuca ovina

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the “native” PgiC1 and “foreign” PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

Structural fluctuations and mechanical stabilities of the metamorphic protein RfaH

RfaH is a compact two-domain bacterial transcription factor that functions both as a regulator of transcription and an enhancer of translation. Underpinning the dual functional roles of RfaH is a partial but dramatic fold switch, which completely transforms the ~50-amino acid C-terminal domain (CTD) from an all-α state to an all-β state. The fold switch of the CTD occurs when RfaH binds to RNA polymerase (RNAP), however, the details of how this structural transformation is triggered is not well understood. Here we use all-atom Monte Carlo simulations to characterize structural fluctuations and mechanical stability properties of the full-length RfaH and the CTD as an isolated fragment. In agreement with experiments, we find that interdomain contacts are crucial for maintaining a stable, all-α CTD in free RfaH. To probe mechanical properties, we use pulling simulations to measure the work required to inflict local deformations at different positions along the chain. The resulting mechanical stability profile reveals that free RfaH can be divided into a "rigid" part and a "soft" part, with a boundary that nearly coincides with the boundary between the two domains. We discuss the potential role of this feature for how fold switching may be triggered by interaction with RNAP.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

CuATSM Protects Against the In Vitro Cytotoxicity of Wild-Type-Like Copper-Zinc Superoxide Dismutase Mutants but not Mutants That Disrupt Metal Binding

Mutations in the SOD1 gene are associated with some forms of familial amyotrophic lateral sclerosis (fALS). There are more than 150 different mutations in the SOD1 gene that have various effects on the copper-zinc superoxide dismutase (SOD1) enzyme structure, including the loss of metal binding and a decrease in dimer affinity. The copper-based therapeutic CuATSM has been proven to be effective at rescuing neuronal cells from SOD1 mutant toxicity and has also increased the life expectancy of mice expressing the human transgenes SOD1^G93A and SOD1^G37R. Furthermore, CuATSM is currently the subject of a phase I/II clinical trial in Australia as a treatment for ALS. To determine if CuATSM protects against a broad variety of SOD1 mutations, we used a well-established cell culture model of SOD1-fALS. NSC-34 cells expressing SOD1-EGFP constructs were treated with CuATSM and examined by time-lapse microscopy. Our results show a concentration-dependent protection of cells expressing mutant SOD1^A4V over the experimental time period. We tested the efficacy of CuATSM on 10 SOD1-fALS mutants and found that while protection was observed in cells expressing pathogenic wild-type-like mutants, cells expressing a truncation mutant or metal binding region mutants were not. We also show that CuATSM rescue is associated with an increase in human SOD1 activity and a decrease in the level of SOD1 aggregation in vitro. In conclusion, CuATSM has shown to be a promising therapeutic for SOD1-associated ALS; however, our in vitro results suggest that the protection afforded varies depending on the SOD1 variant, including negligible protection to mutants with deficient copper binding.

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.

Prediction of Misfolding-Specific Epitopes in SOD1 Using Collective Coordinates

We introduce a global, collective coordinate bias into molecular dynamics simulations that partially unfolds a protein, in order to predict misfolding-specific epitopes based on the regions that locally unfold. Several metrics are used to measure local disorder, including solvent exposed surface area (SASA), native contacts ( Q), and root mean squared fluctuations (RMSF). The method is applied to Cu, Zn superoxide dismutase (SOD1). For this protein, the processes of monomerization, metal loss, and conformational unfolding due to microenvironmental stresses are all separately taken into account. Several misfolding-specific epitopes are predicted, and consensus epitopes are calculated. These predicted epitopes are consistent with the "lower-resolution" peptide sequences used to raise disease-specific antibodies, but the epitopes derived from collective coordinates contain shorter, more refined sequences for the key residues constituting the epitope.

Quantum chemical and molecular mechanics studies on the assessment of interactions between resveratrol and mutant SOD1 (G93A) protein

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease that has been associated with mutations in metalloenzyme superoxide dismutase (SOD1) causing protein structural destabilization and aggregation. However, the mechanistic action and the cure for the disease still remain obscure. Herein, we initially studied the conformational preferences of SOD1 protein structures upon substitution of Ala at Gly93 in comparison with that of wild type. Our results corroborated with the previous experimental studies on the aggregation and the destabilizing activity of mutant SOD1 protein G93A. On the therapeutic point of view, we computationally analyzed the influence of resveratrol, a natural polyphenol widely found in red wine on mutant SOD1 relative to wild type, using molecular docking studies. Further, FMO calculations were performed, using GAMESS to study the pair residual interaction on the wild type and mutant complex systems. Consequently, the resveratrol showed greater interaction with mutant than the wild type. Subsequently, we evaluated the conformational preferences of wild type and mutant complex systems, where the protein conformational structures of mutant that were earlier found to lose their conformational stability was regained, upon binding with resveratrol. Similar trend of results were found on the 2-D free energy landscapes of both the wild type and mutant systems. Hence, the combined biophysical and quantum chemical studies in our study supported the results of previous experimental studies, thereby stipulating an action of resveratrol on mutant SOD1 and paving a way for the design of highly potent effective inhibitors against fALS affecting the mankind.

Molecular mechanisms underlying the impact of mutations in SOD1 on its conformational properties associated with amyotrophic lateral sclerosis as revealed with molecular modelling

Background

So far, little is known about the molecular mechanisms of amyotrophic lateral sclerosis onset and progression caused by SOD1 mutations. One of the hypotheses is based on SOD1 misfolding resulting from mutations and subsequent deposition of its cytotoxic aggregates. This hypothesis is complicated by the fact that known SOD1 mutations of similar clinical effect could be distributed over the whole protein structure.

Results

In this work, a measure of hydrogen bond stability in conformational states was studied with elastic network analysis of 35 SOD1 mutants. Twenty-eight hydrogen bonds were detected in nine of 35 mutants with their stability being significantly different from that with the wild-type. These hydrogen bonds were formed by the amino acid residues known from the literature to be located in contact between SOD1 aggregates. Additionally, residues disposed between copper binding sites of both protein subunits were found from the models to form a stiff core, which can be involved in mechanical impulse transduction between these active centres.

Conclusions

The modelling highlights that both stability of the copper binding site and stability of the dimer can play an important role in ALS progression.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0080-9) contains supplementary material, which is available to authorized users.

A faulty interaction between SOD1 and hCCS in neurodegenerative disease

A proportion of Amyotrophic lateral sclerosis (ALS) cases result from impaired mutant superoxide dismutase-1 (SOD1) maturation. The copper chaperone for SOD1 (hCCS) forms a transient complex with SOD1 and catalyses the final stages of its maturation. We find that a neurodegenerative disease-associated hCCS mutation abrogates the interaction with SOD1 by inhibiting hCCS zinc binding. Analogously, SOD1 zinc loss has a detrimental effect on the formation, structure and disassociation of the hCCS-SOD1 heterodimer. This suggests that hCCS functionality is impaired by ALS mutations that reduce SOD1 zinc affinity. Furthermore, stabilization of wild-type SOD1 by chemical modification including cisplatination, inhibits complex formation. We hypothesize that drug molecules designed to stabilize ALS SOD1 mutants that also target the wild-type form will lead to characteristics common in SOD1 knock-outs. Our work demonstrates the applicability of chromatographic SAXS when studying biomolecules predisposed to aggregation or dissociation; attributes frequently reported for complexes involved in neurodegenerative disease.

Effect of Oxidative Damage on the Stability and Dimerization of Superoxide Dismutase 1

During their life cycle, proteins are subject to different modifications involving reactive oxygen species. Such oxidative damage to proteins may lead to the formation of insoluble aggregates and cytotoxicity and is associated with age-related disorders including neurodegenerative diseases, cancer, and diabetes. Superoxide dismutase 1 (SOD1), a key antioxidant enzyme in human cells, is particularly susceptible to such modifications. Moreover, this homodimeric metalloenzyme has been directly linked to both familial and sporadic amyotrophic lateral sclerosis (ALS), a devastating, late-onset motor neuronal disease, with more than 150 ALS-related mutations in the SOD1 gene. Importantly, oxidatively damaged SOD1 aggregates have been observed in both familial and sporadic forms of the disease. However, the molecular mechanisms as well as potential implications of oxidative stress in SOD1-induced cytotoxicity remain elusive. In this study, we examine the effects of oxidative modification on SOD1 monomer and homodimer stability, the key molecular properties related to SOD1 aggregation. We use molecular dynamics simulations in combination with thermodynamic integration to study microscopic-level site-specific effects of oxidative "mutations" at the dimer interface, including lysine, arginine, proline and threonine carbonylation, and cysteine oxidation. Our results show that oxidative damage of even single residues at the interface may drastically destabilize the SOD1 homodimer, with several modifications exhibiting a comparable effect to that of the most drastic ALS-causing mutations known. Additionally, we show that the SOD1 monomer stability decreases upon oxidative stress, which may lead to partial local unfolding and consequently to increased aggregation propensity. Importantly, these results suggest that oxidative stress may play a key role in development of ALS, with the mutations in the SOD1 gene being an additional factor.

Destabilization of the dimer interface is a common consequence of diverse ALS-associated mutations in metal free SOD1

Neurotoxic misfolding of Cu, Zn-superoxide dismutase (SOD1) is implicated in causing amyotrophic lateral sclerosis, a devastating and incurable neurodegenerative disease. Disease-linked mutations in SOD1 have been proposed to promote misfolding and aggregation by decreasing protein stability and increasing the proportion of less folded forms of the protein. Here we report direct measurement of the thermodynamic effects of chemically and structurally diverse mutations on the stability of the dimer interface for metal free (apo) SOD1 using isothermal titration calorimetry and size exclusion chromatography. Remarkably, all mutations studied, even ones distant from the dimer interface, decrease interface stability, and increase the population of monomeric SOD1. We interpret the thermodynamic data to mean that substantial structural perturbations accompany dimer dissociation, resulting in the formation of poorly packed and malleable dissociated monomers. These findings provide key information for understanding the mechanisms and energetics underlying normal maturation of SOD1, as well as toxic SOD1 misfolding pathways associated with disease. Furthermore, accurate prediction of protein-protein association remains very difficult, especially when large structural changes are involved in the process, and our findings provide a quantitative set of data for such cases, to improve modelling of protein association.

Metal-deficient SOD1 in amyotrophic lateral sclerosis

Mutations to the ubiquitous antioxidant enzyme Cu/Zn superoxide dismutase (SOD1) were the first established genetic cause of the fatal, adult-onset neurodegenerative disease amyotrophic lateral sclerosis (ALS). It is widely accepted that these mutations do not cause ALS via a loss of antioxidant function, but elucidating the alternate toxic gain of function has proven to be elusive. Under physiological conditions, SOD1 binds one copper ion and one zinc ion per monomer to form a highly stable and functional homodimer, but there is now ample evidence to indicate aberrant persistence of SOD1 in an intermediate metal-deficient state may contribute to the protein’s involvement in ALS. This review briefly discusses some of the data to support a role for metal-deficient SOD1 in the development of ALS and some of the outcomes from drug development studies that have aimed to modify the symptoms of ALS by targeting the metal state of SOD1. The implications for the metal state of SOD1 in cases of sporadic ALS that do not involve mutant SOD1 are also discussed.

Interaction between dimer interface residues of native and mutated SOD1 protein: a theoretical study

Cu-Zn superoxide dismutase 1 (SOD1) is a highly conserved bimetallic protein enzyme, used for the scavenging the superoxide radicals (O2 (-)) produced due to aerobic metabolism in the mitochondrial respiratory chain. Over 100 mutations have been identified and found to be in the homodimeric structure of SOD1. The enzyme has to be maintained in its dimeric state for the structural stability and enzymatic activity. From our investigation, we found that the mutations apart from the dimer interface residues are found to affect the dimer stability of protein and hence enhancing the aggregation and misfolding tendency of mutated protein. The homodimeric state of SOD1 is found to be held together by the non-covalent interactions. The molecular dynamics simulation has been used to study the hydrogen bond interactions between the dimer interface residues of the monomers in native and mutated forms of SOD1 in apo- and holo-states. The results obtained by this analysis reveal the fact that the loss of hydrogen bond interactions between the monomers of the dimer is responsible for the reduced stability of the apo- and holo-mutant forms of SOD1. The conformers with dimer interface residues in native and mutated protein obtained by the molecular dynamics simulation is subjected to quantum mechanical study using M052X/6-31G(d) level of theory. The charge transfer between N-H···O interactions in the dimer interface residues were studied. The weak interaction between the monomers of the dimer accounts for the reduced dimerization and enhanced deformation energy in the mutated SOD1 protein.

Structural investigation on the electrostatic loop of native and mutated SOD1 and their interaction with therapeutic compounds

The electrostatic loop of the native and mutated SOD1 protein with single point mutation in the loop is subjected to MD simulation. The structure and electrostatic properties of the native and mutated loops before/after interacting with small compounds are compared.

Altered thiol chemistry in human amyotrophic lateral sclerosis-linked mutants of superoxide dismutase 1

Neurodegenerative diseases share a common characteristic, the presence of intracellular or extracellular deposits of protein aggregates in nervous tissues. Amyotrophic Lateral Sclerosis (ALS) is a severe and fatal neurodegenerative disorder, which affects preferentially motoneurons. Changes in the redox state of superoxide dismutase 1 (SOD1) are associated with the onset and development of familial forms of ALS. In human SOD1 (hSOD1), a conserved disulfide bond and two free cysteine residues can engage in anomalous thiol/disulfide exchange resulting in non-native disulfides, a hallmark of ALS that is related to protein misfolding and aggregation. Because of the many competing reaction pathways, traditional bulk techniques fall short at quantifying individual thiol/disulfide exchange reactions. Here, we adapt recently developed single-bond chemistry techniques to study individual disulfide isomerization reactions in hSOD1. Mechanical unfolding of hSOD1 leads to the formation of a polypeptide loop held by the disulfide. This loop behaves as a molecular jump rope that brings reactive Cys-111 close to the disulfide. Using force-clamp spectroscopy, we monitor nucleophilic attack of Cys-111 at either sulfur of the disulfide and determine the selectivity of the reaction. Disease-causing mutations G93A and A4V show greatly altered reactivity patterns, which may contribute to the progression of familial ALS.

Local unfolding and aggregation mechanisms of SOD1: a Monte Carlo exploration

Copper, zinc superoxide dismutase 1 (SOD1) is a ubiquitous homodimeric enzyme, whose misfolding and aggregation play a potentially key role in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). SOD1 aggregation is thought to be preceded by dimer dissociation and metal loss, but the mechanisms by which the metal-free monomer aggregates remain incompletely understood. Here we use implicit solvent all-atom Monte Carlo (MC) methods to investigate the local unfolding dynamics of the β-barrel-forming SOD1 monomer. Although event-to-event variations are large, on average, we find clear differences in dynamics among the eight strands forming the β-barrel. Most dynamic is the eighth strand, β8, which is located in the dimer interface of native SOD1. For the four strands in or near the dimer interface (β1, β2, β7, and β8), we perform aggregation simulations to assess the propensity of these chain segments to self-associate. We find that β1 and β2 readily self-associate to form intermolecular parallel β-sheets, whereas β8 shows a very low aggregation propensity.

Hypothesis: The unfolding power of protein dielectricity

A hypothesis is proposed on a potential role of protein dielectricity as an unfolding factor in protein-protein interactions. It is suggested that large protein complexes and aggregation seeds can unfold target proteins by virtue of their effect on the dielectric properties of water at the protein-solvent interface. Here, similar to the effect of membrane surfaces, protein surface can cause decrease in the local dielectric constant of solvent and thereby induce structural changes in a target protein approaching this surface. Some potential implementations of this hypothetical mechanism are also discussed.

SOD1 exhibits allosteric frustration to facilitate metal binding affinity

Superoxide dismutase-1 (SOD1) is a ubiquitous, Cu and Zn binding, free-radical defense enzyme whose misfolding and aggregation play a potential key role in amyotrophic lateral sclerosis, an invariably fatal neurodegenerative disease. Over 150 mutations in SOD1 have been identified with a familial form of the disease, but it is presently not clear what unifying features, if any, these mutants share to make them pathogenic. Here, we develop several unique computational assays for probing the thermo-mechanical properties of both ALS-associated and rationally designed SOD1 variants. Allosteric interaction-free energies between residues and metals are calculated, and a series of atomic force microscopy experiments are simulated with variable tether positions to quantify mechanical rigidity "fingerprints" for SOD1 variants. Mechanical fingerprinting studies of a series of C-terminally truncated mutants, along with an analysis of equilibrium dynamic fluctuations while varying native constraints, potential energy change upon mutation, frustratometer analysis, and analysis of the coupling between local frustration and metal binding interactions for a glycine scan of 90 residues together, reveal that the apo protein is internally frustrated, that these internal stresses are partially relieved by mutation but at the expense of metal-binding affinity, and that the frustration of a residue is directly related to its role in binding metals. This evidence points to apo SOD1 as a strained intermediate with "self-allostery" for high metal-binding affinity. Thus, the prerequisites for the function of SOD1 as an antioxidant compete with apo state thermo-mechanical stability, increasing the susceptibility of the protein to misfold in the apo state.