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

Mutation/metal deficiency in the "electrostatic loop" enhanced aggregation process in apo/holo SOD1 variants: implications for ALS diseases

Despite the many mechanisms it has created to prevent unfolding and aggregation of proteins, many diseases are caused by abnormal folding of proteins, which are called misfolding diseases. During this process, proteins undergo structural changes and become stable, insoluble beta-sheet aggregates called amyloid fibrils. Mutations/disruptions in metal ion homeostasis in the ALS-associated metalloenzyme superoxide dismutase (SOD1) reduce conformational stability, consistent with the protein aggregation hypothesis for neurodegenerative diseases. However, the exact mechanism of involvement is not well understood. Hence, to understand the role of mutation/ metal deficiency in SOD1 misfolding and aggregation, we investigated the effects of apo/holo SOD1 variants on structural properties using biophysical/experimental techniques. The MD results support the idea that the mutation/metal deficiency can lead to a change in conformation. The increased content of β-sheet structures in apo/holo SOD1 variants can be attributed to the aggregation tendency, which was confirmed by FTIR spectroscopy and dictionary of secondary structure in proteins (DSSP) results. Thermodynamic studies of GdnHCl showed that metal deficiency/mutation/intramolecular S-S reduction together are required to initiate misfolding/aggregation of SOD1. The results showed that apo/holo SOD1 variants under destabilizing conditions induced amyloid aggregates at physiological pH, which were detected by ThT/ANS fluorescence, as well as further confirmation of amyloid/amorphous species by TEM. This study confirms that mutations in the electrostatic loop of SOD1 lead to structural abnormalities, including changes in hydrophobicity, reduced disulfide bonds, and an increased propensity for protein denaturation. This process facilitates the formation of amyloid/amorphous aggregates ALS-associated.

Predicting binding affinity changes from long-distance mutations using molecular dynamics simulations and Rosetta

Computationally modeling how mutations affect protein-protein binding not only helps uncover the biophysics of protein interfaces, but also enables the redesign and optimization of protein interactions. Traditional high-throughput methods for estimating binding free energy changes are currently limited to mutations directly at the interface due to difficulties in accurately modeling how long-distance mutations propagate their effects through the protein structure. However, the modeling and design of such mutations is of substantial interest as it allows for greater control and flexibility in protein design applications. We have developed a method that combines high-throughput Rosetta-based side-chain optimization with conformational sampling using classical molecular dynamics simulations, finding significant improvements in our ability to accurately predict long-distance mutational perturbations to protein binding. Our approach uses an analytical framework grounded in alchemical free energy calculations while enabling exploration of a vastly larger sequence space. When comparing to experimental data, we find that our method can predict internal long-distance mutational perturbations with a level of accuracy similar to that of traditional methods in predicting the effects of mutations at the protein-protein interface. This work represents a new and generalizable approach to optimize protein free energy landscapes for desired biological functions.

Metal migration and subunit swapping in ALS-linked SOD1: Zn2+ transfer between mutant and wild-type occurs faster than the rate of heterodimerization

The heterodimerization of WT Cu, Zn superoxide dismutase-1 (SOD1), and mutant SOD1 might be a critical step in the pathogenesis of SOD1-linked amyotrophic lateral sclerosis (ALS). Rates and free energies of heterodimerization (ΔG_Het) between WT and ALS-mutant SOD1 in mismatched metalation states-where one subunit is metalated and the other is not-have been difficult to obtain. Consequently, the hypothesis that under-metalated SOD1 might trigger misfolding of metalated SOD1 by "stealing" metal ions remains untested. This study used capillary zone electrophoresis and mass spectrometry to track heterodimerization and metal transfer between WT SOD1, ALS-variant SOD1 (E100K, E100G, D90A), and triply deamidated SOD1 (modeled with N26D/N131D/N139D substitutions). We determined that rates of subunit exchange between apo dimers and metalated dimers-expressed as time to reach 30% heterodimer-ranged from t_30% = 67.75 ± 9.08 to 338.53 ± 26.95 min; free energies of heterodimerization ranged from ΔG_Het = -1.21 ± 0.31 to -3.06 ± 0.12 kJ/mol. Rates and ΔG_Het values of partially metalated heterodimers were more similar to those of fully metalated heterodimers than apo heterodimers, and largely independent of which subunit (mutant or WT) was metal-replete or metal-free. Mass spectrometry and capillary electrophoresis demonstrated that mutant or WT 4Zn-SOD1 could transfer up to two equivalents of Zn²⁺ to mutant or WT apo-SOD1 (at rates faster than the rate of heterodimerization). This result suggests that zinc-replete SOD1 can function as a chaperone to deliver Zn²⁺ to apo-SOD1, and that WT apo-SOD1 might increase the toxicity of mutant SOD1 by stealing its Zn²⁺.

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.