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.

A double point mutation of SOD1 targeting net charge promotes aggregation under destabilizing conditions: Correlation of charge distribution and ALS-provoking mutation

A progressive neurodegenerative illness such as amyotrophic lateral sclerosis (ALS) is characterized by the misfolding and aggregation of human CuZn superoxide dismutase (hSOD1) into amyloid aggregates. Thus, designing strategies for the choice of WT-SOD1 and double mutant (G12D/G138E) with an increased net negative charge can be a good idea to elucidate the pathological mechanism of SOD1 in ALS under some destabilizing conditions. Consequently, we show evidence that protein charge, together with other destabilizing conditions, plays an important role in ALS disease. To achieve this purpose, we use methods, such as spectroscopy and transmission electron microscopy (TEM) to monitor the formation of amyloid aggregation. The specific activity of WT-SOD1 was approximately 1.72 times higher than that of the double mutant. Under amyloidogenic circumstances, structural properties such as local, secondary, oligomeric, and fibrillar structures were explored. The double mutant's far-UV CD spectra displayed a broad minimum peak in the region 213 to 218 nm, suggesting the production of β-rich amyloid fibrils. FTIR spectra of the double mutant samples at different incubation times showed a low-frequency peak around 1630-1640 cm^-1, attributed to a parallel β-sheet. Moreover, CR-binding assay and TEM analysis revealed and confirmed that mutation with an increased repulsive charge promotes the formation of fibrous aggregates. Consequently, ALS mutations with a higher repulsive charge are the apparent exceptions that validate the rule. This findings revealed that the double mutant increases protein aggregation through a novel mechanism, likely involving destabilization of structure and a change in the net negative charge.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

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.

Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway

TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.

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.

Delineating the Aggregation-Prone Hotspot Regions (Peptides) in the Human Cu/Zn Superoxide Dismutase 1

Amyotrophic lateral sclerosis (ALS) is a fatal, incurable neurodegenerative disease described by progressive degeneration of motor neurons. The most common familial form of ALS (fALS) has been associated with mutations in the Cu/Zn superoxide dismutase (SOD1) gene. Mutation-induced misfolding and aggregation of SOD1 is often found in ALS patients. In this work, we probe the aggregation properties of peptides derived from the SOD1. To examine the source of SOD1 aggregation, we have employed a computational algorithm to identify four peptides from the SOD1 protein sequence that aggregates into a fibril. Aided by computational algorithms, we identified four peptides likely involved in SOD1 fibrillization. These four aggregation-prone peptides were 14VQGIINFE21, 30KVWGSIKGL38, 101DSVISLS107, and 147GVIGIAQ153. In addition, the formation of fibril propensities from the identified peptides was investigated through different biophysical techniques. The atomic structures of two fibril-forming peptides from the C-terminal SOD1 showed that the steric zippers formed by 101DSVISLS107 and 147GVIGIAQ153 vary in their arrangement. We also discovered that fALS mutations in the peptide 147GVIGIAQ153 increased the fibril-forming propensity and altered the steric zipper's packing. Thus, our results suggested that the C-terminal peptides of SOD1 have a central role in amyloid formation and might be involved in forming the structural core of SOD1 aggregation observed in vivo.

Genomics-guided targeting of stress granule proteins G3BP1/2 to inhibit SARS-CoV-2 propagation

SARS-CoV-2 nucleocapsid (N) protein undergoes RNA-induced phase separation (LLPS) and sequesters the host key stress granule (SG) proteins, Ras-GTPase-activating protein SH3-domain-binding protein 1 and 2 (G3BP1 and G3BP2) to inhibit SG formation. This will allow viral packaging and propagation in host cells. Based on a genomic-guided meta-analysis, here we identify upstream regulatory elements modulating the expression of G3BP1 and G3BP2 (collectively called G3BP1/2). Using this strategy, we have identified FOXA1, YY1, SYK, E2F-1, and TGFBR2 as activators and SIN3A, SRF, and AKT-1 as repressors of G3BP1/2 genes. Panels of the activators and repressors were then used to identify drugs that change their gene expression signatures. Two drugs, imatinib, and decitabine have been identified as putative modulators of G3BP1/2 genes and their regulators, suggesting their role as COVID-19 mitigation agents. Molecular docking analysis suggests that both drugs bind to G3BP1/2 with a much higher affinity than the SARS-CoV-2 N protein. This study reports imatinib and decitabine as candidate drugs against N protein and G3BP1/2 protein.

Computer analysis of the relation between hydrogen bond stability in SOD1 mutants and the survival time of amyotrophic lateral sclerosis patients

BACKGROUND AND OBJECTIVE

Mutations in the SOD1 protein can lead to the death of motor neurons, which, in turn, causes an incurable disease called amyotrophic lateral sclerosis (ALS). At the same time, the mechanism of the onset and development of this disease is not fully understood and is often contradictory.

METHODS

Accelerated Molecular Dynamics as implemented in the OpenMM library, principal component analysis, regression analysis, random forest method.

RESULTS

The stability of hydrogen bonds in 72 mutants of the SOD1 protein was calculated. Principal component analysis was carried out. Based on ten principal components acting as predictors, a multiple linear regression model was constructed. An analysis of the correlation of these ten principal components with the initial values of the stability of hydrogen bonds made it possible to characterize the contribution of known structurally and functionally important sites in the SOD1 to the scatter of ALS patients' survival time.

CONCLUSION

Such an analysis made it possible to put forward hypotheses about the relationship between the stabilizing and destabilizing effects of mutations in different structurally and functionally important regions of SOD1 with the patients's survival time.

Targeting cathepsins: A potential link between COVID-19 and associated neurological manifestations

Many studies have shown that the lysosomal cathepsins, especially cathepsins B/L (CTSB/L) are required for SARS-CoV-2 entry into host cells. Lysosomal proteases, cathepsins are indispensable for normal health and are involved in several brain disorders occurring at different development age periods. On the other hand, it has been well known that COVID-19 infection is largely associated with several neurological disorders. Taken together these findings and given the high levels of expression of CTSB/L in the brain, we here proposed a reasonable hypothesis about the involvement of CTSB/L in the neurological manifestations linked to COVID-19. Pharmacological inhibitions of the CTSB/L could be a potential therapeutic target to block the virus entry as well as to mitigate the brain disorders. To this end, we utilized the network-based drug repurposing analyses to identify the possible drugs that can target CTSB/L. This study identifies the molecules like cyclosporine, phenytoin, and paclitaxel as potential drugs with binding ability to the CTSB/L. Further, we have performed molecular docking and all-atom molecular dynamics (MD) simulations to investigate the stability of CTSL-drug complexes. The results showed strong and stable binding of drugs with CTSL.

The structural, functional, and dynamic effect of Tau tubulin kinase1 upon a mutation: A neuro-degenerative hotspot

Alzheimer's disease (AD) is a progressive disorder that causes brain cells to degenerate and die. AD is one of the common causes of dementia that leads to a decline in thinking, behavioral and social skills that disrupts a person's ability to function independently. Tau-tubulin kinase1 (TTBK1) is a crucial disease regulating AD protein, which is majorly responsible for the phosphorylation and accumulation of tau protein at specific Serine/Threonine residues found in paired helical filaments, suggesting its role in tauopathy. TTBK1 involvement in many diseases and the restricted expression of TTBK1 to the central nervous system (CNS) makes TTBK1 an attractive therapeutic target for tauopathies. The genetic variations in TTBK1 are primarily involved in the TTBK1 pathogenesis. This study highlighted the destabilizing, damaging and deleterious effect of the mutation R142Q on TTBK1 structure through computational predictions and molecular dynamics simulations. The protein deviation, fluctuations, conformational dynamics, solvent accessibility, hydrogen bonding, and the residue-residue mapping confirmed the mutant effect to cause structural aberrations, suggesting overall destabilization due to the protein mutation. The presence of well-defined free energy minima was observed in TTBK1-wild type, as opposed to that in the R142Q mutant, reflecting structural deterioration. The overall findings from the study reveal that the presence of R142Q mutation on TTBK1 is responsible for the structural instability, leading to disruption of its biological functions. The mutation could be used as future diagnostic markers in treating AD.

Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine

The COVID-19 pandemic caused by SARS-CoV-2 represents a global public health emergency. The entry of SARS-CoV-2 into host cells requires the activation of its spike protein by host cell proteases. The serine protease, TMPRSS2, and cysteine proteases, Cathepsins B/L, activate spike protein and enable SARS-CoV-2 entry to the host cell through two completely different and independent pathways. Therefore, inhibiting either TMPRSS2 or cathepsin B/L may not sufficiently block the virus entry. We here hypothesized that simultaneous targeting of both the entry pathways would be more efficient to block the virus entry rather than targeting the entry pathways individually. To this end, we utilized the network-based drug repurposing analyses to identify the possible common drugs that can target both the entry pathways. This study, for the first time, reports the molecules like cyclosporine, calcitriol, and estradiol as candidate drugs with the binding ability to the host proteases, TMPRSS2, and cathepsin B/L. Next, we analyzed drug-gene and gene-gene interaction networks using 332 human targets of SARS-CoV-2 proteins. The network results indicate that, out of 332 human proteins, cyclosporine interacts with 216 (65%) proteins. Furthermore, we performed molecular docking and all-atom molecular dynamics (MD) simulations to explore the binding of drug with TMPRSS2 and cathepsin L. The molecular docking and MD simulation results showed strong and stable binding of cyclosporine A (CsA) with TMPRSS2 and CTSL genes. The above results indicate cyclosporine as a potential drug molecule, as apart from interacting with SARS-CoV-2 entry receptors, it also interacts with most of SARS-CoV-2 target host genes; thus it could potentially interfere with functions of SARS-CoV-2 proteins in human cells. We here also suggest that these antiviral drugs alone or in combination can simultaneously target both the entry pathways and thus can be considered as a potential treatment option for COVID-19.

"Identification of Nafamostat and VR23 as COVID-19 drug candidates by targeting 3CLpro and PLpro."

The sudden increase in the COVID-19 epidemic affected by novel coronavirus 2019 has jeopardized public health worldwide. Hence the necessities of a drug or therapeutic agent that heal SARS-CoV-2 infections are essential requirements. The viral genome encodes a large Polyprotein, further processed by the main protease/ 3C-like protease (3CL^pro) and papain-like proteases (PL^pro) into 16 nonstructural proteins to form a viral replication complex. These essential functions of 3CL^pro and PL^pro in virus duplication make these proteases a promising target for discovering potential therapeutic candidates and possible treatment for SARS-CoV-2 infection.

This study aimed to screen a unique set of protease inhibitors library against 3CL^pro and PL^pro of the SARS-CoV-2. A molecular docking study was performed using PyRx to reveal the binding affinity of the selected ligands and molecular dynamic simulations were executed to assess the three-dimensional stability of protein-ligand complexes. The pharmacodynamics parameters of the inhibitors were predicted using admetSAR. The top two ligands (Nafamostat and VR23) based on docking scores were selected for further studies. Selected ligands showed excellent pharmacokinetic properties with proper absorption, bioavailability and minimal toxicity. Due to the emerging and efficiency of remdesivir and dexamethasone in healing COVID-19 patients, ADMET properties of the selected ligands were thus compared with it. MD Simulation studies up to 100 ns revealed the ligands' stability at the target proteins' binding site residues. Therefore, Nafamostat and VR23 may provide potential treatment options against SARS-CoV-2 infections by potentially inhibiting virus duplication though more research is warranted.

A computational study of Tat-CDK9-Cyclin binding dynamics and its implication in transcription-dependent HIV latency

HIV is a virus that attacks the T cells. HIV may either actively replicate or become latent within host cells for years. Since HIV uses its own protein Tat to hijack the host CDK9-Cyclin complex for transcription, Tat is implicated in transcription-dependent HIV latency. To quantify the impact of Tat binding, we propose a computational framework to probe the dynamics of the CDK9-Cyclin interface and the ATP pocket reorganization upon binding by different Tat mutants. Specifically, we focus on mutations at three Tat residues P10, W11, and N12 that are known to interact directly with CDK9 based on the crystal structure of the Tat-CDK9-Cyclin complex. Our molecular dynamics simulations show that the CDK9-Cyclin interface becomes slightly weaker for P10S and W11R mutants but tighter for the K12N mutant. Furthermore, the side chain orientation of residue K48 in the ATP pocket of CDK9 is similar to the inactive state in P10S and W11R simulations, but similar to the active state in K12N simulations. These are consistent with some existing but puzzling observations of latency for these mutants. This framework may hence help gain a better understanding of the role of Tat in the transcription-dependent HIV latency establishment.

Trifluoroethanol Partially Unfolds G93A SOD1 Leading to Protein Aggregation: A Study by Native Mass Spectrometry and FPOP Protein Footprinting

Misfolding of Cu, Zn superoxide dismutase (SOD1) variants may lead to protein aggregation and ultimately amyotrophic lateral sclerosis (ALS). The mechanism and protein conformational changes during this process are complex and remain unclear. To study SOD1 variant aggregation at the molecular level and in solution, we chemically induced aggregation of a mutant variant (G93A SOD1) with trifluoroethanol (TFE) and used both native mass spectrometry (MS) to analyze the intact protein and fast photochemical oxidation of proteins (FPOP) to characterize the structural changes induced by TFE. We found partially unfolded G93A SOD1 monomers prior to oligomerization and identified regions of the N-terminus, C-terminus, and strands β5, β6 accountable for the partial unfolding. We propose that exposure of hydrophobic interfaces of these unstructured regions serves as a precursor to aggregation. Our results provide a possible mechanism and molecular basis for ALS-linked SOD1 misfolding and aggregation.

Haplotype Analysis of the First A4V-SOD1 Spanish Family: Two Separate Founders or a Single Common Founder?

Despite the genetic heterogeneity reported in familial amyotrophic lateral sclerosis (ALS) (fALS), Cu/Zn superoxide-dismutase (SOD1) gene mutations are the second most common cause of the disease, accounting for around 20% of all families (ALS1) and isolated sporadic cases (sALS). At least 186 different mutations in the SOD1 gene have been reported to date. The possibility of a single founder and separate founders have been investigated for D90A (p.D91A) and A4V (p.A5V), the most common mutations worldwide. High-throughput single nucleotide polymorphism genotyping studies have suggested two founders for A4V (one for the Amerindian population and another for the European population) although the possibility that the two populations are descended from a single ancient founder cannot be ruled out. We used 15 genetic variants spanning the human chromosome 21 from the SOD1 gene to the SCAF4 gene, comparing them with the population reference panels, to demonstrate that the first A4V Spanish pedigree shared the genetic background reported in the European population.

Rational design of linear tripeptides against the aggregation of human mutant SOD1 protein causing amyotrophic lateral sclerosis

Formation of protein aggregation is considered a hallmark feature of various neurological diseases. Amyotrophic lateral sclerosis is one such devastating neurodegenerative disorder characterized by mutation in Cu/Zn superoxide dismutase protein (SOD1). In our study, we contemplated the most aggregated and pathogenic mutant A4V in a viewpoint of finding a therapeutic regime by inhibiting the formation of the aggregates with the aid of tripeptides since new perspectives in the field of drug design in the current era are being focused on peptide-based drugs. Reports from the experimental study have stipulated that the SOD1 derived peptide, "LSGDHCIIGRTLVVHEKADD" was found to have the inhibitory activity against aggregated SOD1 protein. Moreover, it was determined that the hexapeptide, "LSGDHC" was the key factor in inhibiting the aggregates of SOD1. Accordingly, we utilized the computerized algorithms and programs on determining the binding efficiency and inhibitory activity of hexapeptide on mutant SOD1. Following that, we incorporated a cutting-edge methodology with the use of molecular docking, affinity predictions, alanine scanning, steered molecular dynamics (SMD) and discrete molecular dynamics (DMD) in designing the de novo tripeptides, which could act against the aggregated mutant SOD1 protein. Upon examining the results from the various conformational studies, we identified that CGH had an enhanced binding affinity and inhibitory activity against the aggregated mutant SOD1 protein than other tripeptides and hexapeptide. Thus, our study could be a lead for state-of-the-art design in peptide-based drugs for doctoring the cureless ALS disorder.

Unfolding Transitions of Peripheral Subunit Binding Domains Show Cooperative Behavior

Characterization of native, intermediate, and denatured states is crucial for understanding the factors influencing the stability of proteins. We have carried out molecular dynamics simulations to study the unfolding of three peripheral subunit binding domains (PSBDs): E. coli BBL, Bacillus stearothermophilus E3BD, and human hbSBD, at three different temperatures: 300, 330, and 400 K, and in the presence of two solvents: water and 5 M guanidinium hydrochloride (GndCl) solution. These proteins share similar folds, with two parallel helices, maintained via a hydrophobic core comprising residues from their interconnecting loop. BBL is more sensitive to thermal and chemical denaturation in comparison to hbSBD, and E3BD is the most stable of all of the three proteins. The effect of temperature on the stability of these proteins is more pronounced in "water-only" simulations compared to that in the presence of guanidium hydrochloride in high concentrations. Our results show cooperative unfolding transitions of these proteins, which are triggered by an initial melting of the C-terminal helix H2. The consequent loss of interhelical interactions or native contacts, as observed, leads to the subsequent melting of the N-terminal helix H1.

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.

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.