Structural and thermodynamic effects of post-translational modifications in mutant and wild type Cu, Zn superoxide dismutase

Direct observation of prion-like propagation of protein misfolding templated by pathogenic mutants

Many neurodegenerative diseases feature misfolded proteins that propagate via templated conversion of natively folded molecules. However, crucial questions about how such prion-like conversion occurs and what drives it remain unsolved, partly because technical challenges have prevented direct observation of conversion for any protein. We observed prion-like conversion in single molecules of superoxide dismutase-1 (SOD1), whose misfolding is linked to amyotrophic lateral sclerosis. Tethering pathogenic misfolded SOD1 mutants to wild-type molecules held in optical tweezers, we found that the mutants vastly increased misfolding of the wild-type molecule, inducing multiple misfolded isoforms. Crucially, the pattern of misfolding was the same in the mutant and converted wild-type domains and varied when the misfolded mutant was changed, reflecting the templating effect expected for prion-like conversion. Ensemble measurements showed decreased enzymatic activity in tethered heterodimers as conversion progressed, mirroring the single-molecule results. Antibodies sensitive to disease-specific epitopes bound to the converted protein, implying that conversion produced disease-relevant misfolded conformers.

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.

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.

Computational Insights into the Unfolding of a Destabilized Superoxide Dismutase 1 Mutant

In this work, we investigate the β-barrel of superoxide dismutase 1 (SOD1) in a mutated form, the isoleucine 35 to alanine (I35A) mutant, commonly used as a model system to decipher the role of the full-length apoSOD1 protein in amyotrophic lateral sclerosis (ALS). It is known from experiments that the mutation reduces the stability of the SOD1 barrel and makes it largely unfolded in the cell at 37 degrees Celsius. We deploy state-of-the-art computational machinery to examine the thermal destabilization of the I35A mutant by comparing two widely used force fields, Amber a99SB-disp and CHARMM36m. We find that only the latter force field, when combined with the Replica Exchange with Solute Scaling (REST2) approach, reproduces semi-quantitatively the experimentally observed shift in the melting between the original and the mutated SOD1 barrel. In addition, we analyze the unfolding process and the conformational landscape of the mutant, finding these largely similar to those of the wildtype. Nevertheless, we detect an increased presence of partially misfolded states at ambient temperatures. These states, featuring conformational changes in the region of the β-strands β4-β6, might provide a pathway for nonnative aggregation.

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.

β-Methylamino-L-alanine substitution of serine in SOD1 suggests a direct role in ALS etiology

Exposure to the environmental toxin β-methylamino-L-alanine (BMAA) is linked to amyotrophic lateral sclerosis (ALS), but its disease-promoting mechanism remains unknown. We propose that incorporation of BMAA into the ALS-linked protein Cu,Zn superoxide dismutase (SOD1) upon translation promotes protein misfolding and aggregation, which has been linked to ALS onset and progression. Using molecular simulation and predictive energetic computation, we demonstrate that substituting any serine with BMAA in SOD1 results in structural destabilization and aberrant dynamics, promoting neurotoxic SOD1 aggregation. We propose that translational incorporation of BMAA into SOD1 is directly responsible for its toxicity in neurodegeneration, and BMAA modification of SOD1 may serve as a biomarker of ALS.

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.

Molecular binding response of naringin and naringenin to H46R mutant SOD1 protein in combating protein aggregation using density functional theory and discrete molecular dynamics

Amyotrophic lateral sclerosis (ALS) is a calamitous neurodegenerative disorder characterized by denervation of upper and lower motor neurons. Numerous hypotheses suggest that toxic protein misfolding and aggregation cause ALS, similar to that of other neurodegenerative diseases, such as Alzheimers and Parkinsons. Protruding causes of familial ALS are mutations in the gene encoding Cu/Zn superoxide dismutase-1 (SOD1), which decrease protein stability and endorse protein aggregation. Thus, the interference concerning aggregate formation and destabilization in SOD1 is considered to be an impending therapeutic strategy. In this work, we utilized computational chemistry methods to initially study the effect of substitution mutation, His46Arg on SOD1 protein. Further, we described the interaction of two naturally occurring polyphenol compounds, naringin and naringenin with mutant SOD1 that is regarded to hinder the protein aggregation. Subsequently, the use of quantum chemical and molecular mechanics calculations speculated that naringin had a strong binding affinity with mutant SOD1 and impeded the formation of toxic aggregates than that of naringenin. Ultimately, we could conjecture that ingesting of polyphenol-rich foods in ALS patients may be regarded to improvise their living. Moreover, the findings from our study could pave a way in the field of structure-based drug design in developing potential anti-aggregation inhibitors against incurable ALS, affecting the human population.

RNA Fibers as Optimized Nanoscaffolds for siRNA Coordination and Reduced Immunological Recognition

RNA is a versatile biomaterial that can be used to engineer nanoassemblies for personalized treatment of various diseases. Despite promising advancements, the design of RNA nanoassemblies with minimal recognition by the immune system remains a major challenge. Here, an approach is reported to engineer RNA fibrous structures to operate as a customizable platform for efficient coordination of siRNAs and for maintaining low immunostimulation. Functional RNA fibers are studied in silico and their formation is confirmed by various experimental techniques and visualized by atomic force microscopy (AFM). It is demonstrated that the RNA fibers offer multiple advantages among which are: i) programmability and modular design that allow for simultaneous controlled delivery of multiple siRNAs and fluorophores, ii) reduced immunostimulation when compared to other programmable RNA nanoassemblies, and iii) simple production protocol for endotoxin-free fibers with the option of their cotranscriptional assembly. Furthermore, it is shown that functional RNA fibers can be efficiently delivered with various organic and inorganic carriers while retaining their structural integrity in cells. Specific gene silencing triggered by RNA fibers is assessed in human breast cancer and melanoma cell lines, with the confirmed ability of functional fibers to selectively target single nucleotide mutations.

Partially native intermediates mediate misfolding of SOD1 in single-molecule folding trajectories

Prion-like misfolding of superoxide dismutase 1 (SOD1) is associated with the disease ALS, but the mechanism of misfolding remains unclear, partly because misfolding is difficult to observe directly. Here we study the most misfolding-prone form of SOD1, reduced un-metallated monomers, using optical tweezers to measure unfolding and refolding of single molecules. We find that the folding is more complex than suspected, resolving numerous previously undetected intermediate states consistent with the formation of individual β-strands in the native structure. We identify a stable core of the protein that unfolds last and refolds first, and directly observe several distinct misfolded states that branch off from the native folding pathways at specific points after the formation of the stable core. Partially folded intermediates thus play a crucial role mediating between native and non-native folding. These results suggest an explanation for SOD1's propensity for prion-like misfolding and point to possible targets for therapeutic intervention.

A Structural Model for Vinculin Insertion into PIP2-Containing Membranes and the Effect of Insertion on Vinculin Activation and Localization

Vinculin, a scaffolding protein that localizes to focal adhesions (FAs) and adherens junctions, links the actin cytoskeleton to the adhesive super-structure. While vinculin binds to a number of cytoskeletal proteins, it can also associate with phosphatidylinositol 4,5-bisphosphate (PIP₂) to drive membrane association. To generate a structural model for PIP₂-dependent interaction of vinculin with the lipid bilayer, we conducted lipid-association, nuclear magnetic resonance, and computational modeling experiments. We find that two basic patches on the vinculin tail drive membrane association: the basic collar specifically recognizes PIP₂, while the basic ladder drives association with the lipid bilayer. Vinculin mutants with defects in PIP₂-dependent liposome association were then expressed in vinculin knockout murine embryonic fibroblasts. Results from these analyses indicate that PIP₂ binding is not required for localization of vinculin to FAs or FA strengthening, but is required for vinculin activation and turnover at FAs to promote its association with the force transduction FA nanodomain.

Discovery of Age-Related Protein Folding Stability Differences in the Mouse Brain Proteome

Described here is the application of thermodynamic stability measurements to study age-related differences in the folding and stability of proteins in a rodent model of aging. Thermodynamic stability profiles were generated for 809 proteins in brain cell lysates from mice, aged 6 (n = 7) and 18 months (n = 9) using the Stability of Proteins from Rates of Oxidation (SPROX) technique. The biological variability of the protein stability measurements was low and within the experimental error of SPROX. A total of 83 protein hits were detected with age-related stability differences in the brain samples. Remarkably, the large majority of the brain protein hits were destabilized in the old mice, and the hits were enriched in proteins that have slow turnover rates (p < 0.07). Furthermore, 70% of the hits have been previously linked to aging or age-related diseases. These results help validate the use of thermodynamic stability measurements to capture relevant age-related proteomic changes and establish a new biophysical link between these proteins and aging.

A Phosphomimetic Mutation Stabilizes SOD1 and Rescues Cell Viability in the Context of an ALS-Associated Mutation

The majority of amyotrophic lateral sclerosis (ALS)-related mutations in the enzyme Cu,Zn superoxide dismutase (SOD1), as well as a post-translational modification, glutathionylation, destabilize the protein and lead to a misfolded oligomer that is toxic to motor neurons. The biophysical role of another physiological SOD1 modification, T2-phosphorylation, has remained a mystery. Here, we find that a phosphomimetic mutation, T2D, thermodynamically stabilizes SOD1 even in the context of a strongly SOD1-destabilizing mutation, A4V, one of the most prevalent and aggressive ALS-associated mutations in North America. This stabilization protects against formation of toxic SOD oligomers and positively impacts motor neuron survival in cellular assays. We solve the crystal structure of T2D-SOD1 and explain its stabilization effect using discrete molecular dynamics (DMD) simulations. These findings imply that T2-phosphorylation may be a plausible innate cellular protection response against SOD1-induced cytotoxicity, and stabilizing the SOD1 native conformation might offer us viable pharmaceutical strategies against currently incurable ALS.

Computational investigation of curcumin, a natural polyphenol that inhibits the destabilization and the aggregation of human SOD1 mutant (Ala4Val)

Curcumin inhibits the aberrant aggregation in mutant SOD1 protein, thereby decreasing the propensity of β-sheets and the toxicity level.

In silico prediction and characterization of protein post-translational modifications

This review outlines the computational approaches and procedures for predicting post translational modification (PTM)-induced changes in protein conformation and their influence on protein function(s), the latter being assessed as differential affinity in interaction with either low (ligands for receptors or transporters, substrates for enzymes) or high molecular mass molecules (proteins or nucleic acids in supramolecular assemblies). The scope for an in silico approach is discussed against a summary of the in vitro evidence on the structural and functional outcome of protein PTM.

Global analysis of protein folding thermodynamics for disease state characterization

Current methods for the large-scale characterization of disease states generally rely on the analysis of gene and/or protein expression levels. These existing methods fail to detect proteins with disease-related functions and unaltered expression levels. Here we describe the large-scale use of thermodynamic measurements of protein folding and stability for the characterization of disease states. Using the Stable Isotope Labeling with Amino Acids in Cell Culture and Stability of Proteins from Rates of Oxidation (SILAC-SPROX) technique, we assayed ∼800 proteins for protein folding and stability changes in three different cell culture models of breast cancer including the MCF-10A, MCF-7, and MDA-MB-231 cell lines. The thermodynamic stability profiles generated here created distinct molecular markers to differentiate the three cell lines, and a significant fraction (∼45%) of the differentially stabilized proteins did not have altered expression levels. Thus, the differential thermodynamic profiling strategy reported here created novel molecular signatures of breast cancer and provided additional insight into the molecular basis of the disease. Our results establish the utility of protein folding and stability measurements for the study of disease processes, and they suggest that such measurements may be useful for biomarker discovery in disease.

Genotype-property patient-phenotype relations suggest that proteome exhaustion can cause amyotrophic lateral sclerosis

Late-onset neurodegenerative diseases remain poorly understood as search continues for the perceived pathogenic protein species. Previously, variants in Superoxide Dismutase 1 (SOD1) causing Amyotrophic Lateral Sclerosis (ALS) were found to destabilize and reduce net charge, suggesting a pathogenic aggregation mechanism. This paper reports analysis of compiled patient data and experimental and computed protein properties for variants of human SOD1, a major risk factor of ALS. Both stability and reduced net charge correlate significantly with disease, with larger significance than previously observed. Using two independent methods and two data sets, a probability < 3% (t-statistical test) is found that ALS-causing mutations share average stability with all possible 2907 SOD1 mutations. Most importantly, un-weighted patient survival times correlate strongly with the misfolded/unfolded protein copy number, expressed as an exponential function of the experimental stabilities (R² = 0.31, p = 0.002), and this phenotype is further aggravated by charge (R² = 0.51, p = 1.8 x 10−5). This finding suggests that disease relates to the copy number of misfolded proteins. Exhaustion of motor neurons due to expensive protein turnover of misfolded protein copies is consistent with the data but can further explain e.g. the expression-dependence of SOD1 pathogenicity, the lack of identification of a molecular toxic mode, elevated SOD1 mRNA levels in sporadic ALS, bioenergetic effects and increased resting energy expenditure in ALS patients, genetic risk factors affecting RNA metabolism, and recent findings that a SOD1 mutant becomes toxic when proteasome activity is recovered after washout of a proteasome inhibitor. Proteome exhaustion is also consistent with energy-producing mitochondria accumulating at the neuromuscular junctions where ALS often initiates. If true, this exhaustion mechanism implies a complete change of focus in treatment of ALS towards actively nursing the energy state and protein turnover of the motor neurons.

Phosphorylation as conformational switch from the native to amyloid state: Trp-cage as a protein aggregation model

The 20 residue long Trp-cage miniprotein is an excellent model for both computational and experimental studies of protein folding and stability. Recently, great attention emerged to study disease-related protein misfolding, aggregation, and amyloid formation, with the aim of revealing their structural and thermodynamic background. Trp-cage is sensitive to both environmental and structure-modifying effects. It aggregates with ease upon structure destabilization, and thus it is suitable for modeling aggregation and amyloid formation. Here, we characterize the amyloid formation of several sequence modified and side-chain phosphorylated Trp-cage variants. We applied NMR, circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopies, molecular dynamics (MD) simulations, and transmission electron microscopy (TEM) in conjunction with thioflavin-T (ThT) fluorescence measurements to reveal the structural consequences of side-chain phosphorylation. We demonstrate that the native fold is destabilized upon serine phosphorylation, and the resultant highly dynamic structures form amyloid-like ordered aggregates with high intermolecular β-structure content. The only exception is the D9S(P) variant, which follows an alternative aggregation process by forming thin fibrils, presenting a CD spectrum of PPII helix, and showing low ThT binding capability. We propose a complex aggregation model for these Trp-cage miniproteins. This model assumes an additional aggregated state, a collagen triple helical form that can precede amyloid formation. The phosphorylation of a single serine residue serves as a conformational switch, triggering aggregation, otherwise mediated by kinases in cell. We show that Trp-cage miniprotein is indeed a realistic model of larger globular systems of composite folding and aggregation landscapes and helps us to understand the fundamentals of deleterious protein aggregation and amyloid formation.

Oxidation of the tryptophan 32 residue of human superoxide dismutase 1 caused by its bicarbonate-dependent peroxidase activity triggers the non-amyloid aggregation of the enzyme

The role of oxidative post-translational modifications of human superoxide dismutase 1 (hSOD1) in the amyotrophic lateral sclerosis (ALS) pathology is an attractive hypothesis to explore based on several lines of evidence. Among them, the remarkable stability of hSOD1(WT) and several of its ALS-associated mutants suggests that hSOD1 oxidation may precede its conversion to the unfolded and aggregated forms found in ALS patients. The bicarbonate-dependent peroxidase activity of hSOD1 causes oxidation of its own solvent-exposed Trp(32) residue. The resulting products are apparently different from those produced in the absence of bicarbonate and are most likely specific for simian SOD1s, which contain the Trp(32) residue. The aims of this work were to examine whether the bicarbonate-dependent peroxidase activity of hSOD1 (hSOD1(WT) and hSOD1(G93A) mutant) triggers aggregation of the enzyme and to comprehend the role of the Trp(32) residue in the process. The results showed that Trp(32) residues of both enzymes are oxidized to a similar extent to hSOD1-derived tryptophanyl radicals. These radicals decayed to hSOD1-N-formylkynurenine and hSOD1-kynurenine or to a hSOD1 covalent dimer cross-linked by a ditryptophan bond, causing hSOD1 unfolding, oligomerization, and non-amyloid aggregation. The latter process was inhibited by tempol, which recombines with the hSOD1-derived tryptophanyl radical, and did not occur in the absence of bicarbonate or with enzymes that lack the Trp(32) residue (bovine SOD1 and hSOD1(W32F) mutant). The results support a role for the oxidation products of the hSOD1-Trp(32) residue, particularly the covalent dimer, in triggering the non-amyloid aggregation of hSOD1.

Non-native soluble oligomers of Cu/Zn superoxide dismutase (SOD1) contain a conformational epitope linked to cytotoxicity in amyotrophic lateral sclerosis (ALS)

Soluble misfolded Cu/Zn superoxide dismutase (SOD1) is implicated in motor neuron death in amyotrophic lateral sclerosis (ALS); however, the relative toxicities of the various non-native species formed by SOD1 as it misfolds and aggregates are unknown. Here, we demonstrate that early stages of SOD1 aggregation involve the formation of soluble oligomers that contain an epitope specific to disease-relevant misfolded SOD1; this epitope, recognized by the C4F6 antibody, has been proposed as a marker of toxic species. Formation of potentially toxic oligomers is likely to be exacerbated by an oxidizing cellular environment, as evidenced by increased oligomerization propensity and C4F6 reactivity when oxidative modification by glutathione is present at Cys-111. These findings suggest that soluble non-native SOD1 oligomers, rather than native-like dimers or monomers, share structural similarity to pathogenic misfolded species found in ALS patients and therefore represent potential cytotoxic agents and therapeutic targets in ALS.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

Computing stability effects of mutations in human superoxide dismutase 1

Protein stability is affected in several diseases and is of substantial interest in efforts to correlate genotypes to phenotypes. Superoxide dismutase 1 (SOD1) is a suitable test case for such correlations due to its abundance, stability, available crystal structures and thermochemical data, and physiological importance. In this work, stability changes of SOD1 mutations were computed with five methods, CUPSAT, I-Mutant2.0, I-Mutant3.0, PoPMuSiC, and SDM, with emphasis on structural sensitivity as a potential issue in structure-based protein calculation. The large correlation between experimental literature data of SOD1 dimers and monomers (r = 0.82) suggests that mutations in separate protein monomers are mostly additive. PoPMuSiC was most accurate (typical MAE ~ 1 kcal/mol, r ~ 0.5). The relative performance of the methods was not very structure-dependent, and the more accurate methods also displayed less structural sensitivity, with the standard deviation from different high-resolution structures down to ~0.2 kcal/mol. Structures of variable resolution and number of protein copies locally affected specific sites, emphasizing the use of state-relevant crystal structures when such sites are of interest, but had little impact on overall batch estimates. Protein-interaction effects (as a mimic of crystal packing) were small for the more accurate methods. Thus, batch computations, relevant to, e.g., comparisons of disease/nondisease mutant sets or different clades in phylogenetic trees, are much more significant than single mutant calculations and may be the only meaningful way to computationally bridge the genotype-phenotype gap of proteomics. Finally, mutations involving glycine were most difficult to model, of relevance to future method improvement. This could be due to structure changes (glycine has a low structural propensity) or water colocalization with glycine.

Comparative structural and conformational studies on H43R and W32F mutants of copper-zinc superoxide dismutase by molecular dynamics simulation

Recently, mutations in copper-zinc superoxide dismutase (SOD1) have been linked to familial amyotrophic lateral sclerosis (fALS), a progressive neurodegenerative disease involving motor neuron loss, paralysis and death. It is mainly due to protein misfolding and aggregation resulting from the enhanced peroxidase activity of SOD1 mutants. In this study, we have carried out a 20 ns molecular dynamics simulation for wild type (WT), H43R and W32F mutated SOD1's dimer and compared their structure and conformational properties by extracting several quantitative properties from the trajectory to understand the pathology of fALS disease. Our results show considerable differences in H43R compared to WT and W32F mutated SOD1, such as increasing distances between the critical residues results in open conformation at the active site, strong fluctuations in the important loops (Zinc and electrostatic loops) and weakening of important hydrogen bonds especially between N (His 43/Arg 43) and carbonyl oxygen (His 120) in agreement with the experimental report. The calculated buried surface area of dimer interface for WT, H43R and W32F are 682, 726 and 657 Å(2) respectively, representing the loss of dimerization in H43R. Essential dynamics reveal that overall motions of WT and W32F are mainly involved in three to four eigenvectors, but in H43R the overall motions are mainly in the first eigenvector. These data thus provide a unifying description for the structural destabilization, enhanced peroxidase activity, loss of dismutation activity and increase in aggregation propensity in the pathology of fALS diseases.

Deamidation of asparagine to aspartate destabilizes Cu, Zn superoxide dismutase, accelerates fibrillization, and mirrors ALS-linked mutations

The reactivity of asparagine residues in Cu, Zn superoxide dismutase (SOD1) to deamidate to aspartate remains uncharacterized; its occurrence in SOD1 has not been investigated, and the biophysical effects of deamidation on SOD1 are unknown. Deamidation is, nonetheless, chemically equivalent to Asn-to-Asp missense mutations in SOD1 that cause amyotrophic lateral sclerosis (ALS). This study utilized computational methods to identify three asparagine residues in wild-type (WT) SOD1 (i.e., N26, N131, and N139) that are predicted to undergo significant deamidation (i.e., to >20%) on time scales comparable to the long lifetime (>1 year) of SOD1 in large motor neurons. Site-directed mutagenesis was used to successively substitute these asparagines with aspartate (to mimic deamidation) according to their predicted deamidation rate, yielding: N26D, N26D/N131D, and N26D/N131D/N139D SOD1. Differential scanning calorimetry demonstrated that the thermostability of N26D/N131D/N139D SOD1 is lower than WT SOD1 by ~2-8 °C (depending upon the state of metalation) and <3 °C lower than the ALS mutant N139D SOD1. The triply deamidated analog also aggregated into amyloid fibrils faster than WT SOD1 by ~2-fold (p < 0.008**) and at a rate identical to ALS mutant N139D SOD1 (p > 0.2). A total of 534 separate amyloid assays were performed to generate statistically significant comparisons of aggregation rates among WT and N/D SOD1 proteins. Capillary electrophoresis and mass spectrometry demonstrated that ~23% of N26 is deamidated to aspartate (iso-aspartate was undetectable) in a preparation of WT human SOD1 (isolated from erythrocytes) that has been used for decades by researchers as an analytical standard. The deamidation of asparagine--an analytically elusive, sub-Dalton modification--represents a plausible and overlooked mechanism by which WT SOD1 is converted to a neurotoxic isoform that has a similar structure, instability, and aggregation propensity as ALS mutant N139D SOD1.

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.

Competitive binding of natural amphiphiles with graphene derivatives

Understanding the transformation of graphene derivatives by natural amphiphiles is essential for elucidating the biological and environmental implications of this emerging class of engineered nanomaterials. Using rapid discrete-molecular-dynamics simulations, we examined the binding of graphene and graphene oxide with peptides, fatty acids, and cellulose, and complemented our simulations by experimental studies of Raman spectroscopy, FTIR, and UV-Vis spectrophotometry. Specifically, we established a connection between the differential binding and the conformational flexibility, molecular geometry, and hydrocarbon content of the amphiphiles. Importantly, our dynamics simulations revealed a Vroman-like competitive binding of the amphiphiles for the graphene oxide substrate. This study provides a mechanistic basis for addressing the transformation, evolution, transport, biocompatibility, and toxicity of graphene derivatives in living systems and the natural environment.

Energetics of oligomeric protein folding and association

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

Abnormal SDS-PAGE migration of cytosolic proteins can identify domains and mechanisms that control surfactant binding

The amino acid substitution or post-translational modification of a cytosolic protein can cause unpredictable changes to its electrophoretic mobility during SDS-PAGE. This type of "gel shifting" has perplexed biochemists and biologists for decades. We identify a mechanism for "gel shifting" that predominates among a set of ALS (amyotrophic lateral sclerosis) mutant hSOD1 (superoxide dismutase) proteins, post-translationally modified hSOD1 proteins, and homologous SOD1 proteins from different organisms. By first comparing how 39 amino acid substitutions throughout hSOD1 affected SDS-PAGE migration, we found that substitutions that caused gel shifting occurred within a single polyacidic domain (residues ~80-101), and were nonisoelectric. Substitutions that decreased the net negative charge of domain 80-101 increased migration; only one substitution increased net negative charge and slowed migration. Capillary electrophoresis, circular dichroism, and size exclusion chromatography demonstrated that amino acid substitutions increase migration during SDS-PAGE by promoting the binding of three to four additional SDS molecules, without significantly altering the secondary structure or Stokes radius of hSOD1-SDS complexes. The high negative charge of domain 80-101 is required for SOD1 gel shifting: neutralizing the polyacidic domain (via chimeric mouse-human SOD1 fusion proteins) inhibited amino acid substitutions from causing gel shifting. These results demonstrate that the pattern of gel shifting for mutant cytosolic proteins can be used to: (i) identify domains in the primary structure that control interactions between denatured cytosolic proteins and SDS and (ii) identify a predominant chemical mechanism for the interaction (e.g., hydrophobic vs. electrostatic).

Amyotrophic lateral sclerosis: new insights into underlying molecular mechanisms and opportunities for therapeutic intervention

Recent years have witnessed a renewed interest in the pathogenic mechanisms of amyotrophic lateral sclerosis (ALS), a late-onset progressive degeneration of motor neurons. The discovery of new genes associated with the familial form of the disease, along with a deeper insight into pathways already described for this disease, has led scientists to reconsider previous postulates. While protein misfolding, mitochondrial dysfunction, oxidative damage, defective axonal transport, and excitotoxicity have not been dismissed, they need to be re-examined as contributors to the onset or progression of ALS in the light of the current knowledge that the mutations of proteins involved in RNA processing, apparently unrelated to the previous "old partners," are causative of the same phenotype. Thus, newly envisaged models and tools may offer unforeseen clues on the etiology of this disease and hopefully provide the key to treatment.

Structural and dynamic determinants of protein-peptide recognition

Protein-peptide interactions play important roles in many cellular processes, including signal transduction, trafficking, and immune recognition. Protein conformational changes upon binding, an ill-defined peptide binding surface, and the large number of peptide degrees of freedom make the prediction of protein-peptide interactions particularly challenging. To address these challenges, we perform rapid molecular dynamics simulations in order to examine the energetic and dynamic aspects of protein-peptide binding. We find that, in most cases, we recapitulate the native binding sites and native-like poses of protein-peptide complexes. Inclusion of electrostatic interactions in simulations significantly improves the prediction accuracy. Our results also highlight the importance of protein conformational flexibility, especially side-chain movement, which allows the peptide to optimize its conformation. Our findings not only demonstrate the importance of sufficient sampling of the protein and peptide conformations, but also reveal the possible effects of electrostatics and conformational flexibility on peptide recognition.

An over-oxidized form of superoxide dismutase found in sporadic amyotrophic lateral sclerosis with bulbar onset shares a toxic mechanism with mutant SOD1

Recent studies suggest that Cu/Zn superoxide dismutase (SOD1) could be pathogenic in both familial and sporadic amyotrophic lateral sclerosis (ALS) through either inheritable or nonheritable modifications. The presence of a misfolded WT SOD1 in patients with sporadic ALS, along with the recently reported evidence that reducing SOD1 levels in astrocytes derived from sporadic patients inhibits astrocyte-mediated toxicity on motor neurons, suggest that WT SOD1 may acquire toxic properties similar to familial ALS-linked mutant SOD1, perhaps through posttranslational modifications. Using patients' lymphoblasts, we show here that indeed WT SOD1 is modified posttranslationally in sporadic ALS and is iper-oxidized (i.e., above baseline oxidation levels) in a subset of patients with bulbar onset. Derivatization analysis of oxidized carbonyl compounds performed on immunoprecipitated SOD1 identified an iper-oxidized SOD1 that recapitulates mutant SOD1-like properties and damages mitochondria by forming a toxic complex with mitochondrial Bcl-2. This study conclusively demonstrates the existence of an iper-oxidized SOD1 with toxic properties in patient-derived cells and identifies a common SOD1-dependent toxicity between mutant SOD1-linked familial ALS and a subset of sporadic ALS, providing an opportunity to develop biomarkers to subclassify ALS and devise SOD1-based therapies that go beyond the small group of patients with mutant SOD1.

Amyotrophic lateral sclerosis: update and new developments

Amyotrophic lateral sclerosis (ALS) is the most common form of motor neuron disease. It is typically characterized by adult-onset degeneration of the upper and lower motor neurons, and is usually fatal within a few years of onset. A subset of ALS patients has an inherited form of the disease, and a few of the known mutant genes identified in familial cases have also been found in sporadic forms of ALS. Precisely how the diverse ALS-linked gene products dictate the course of the disease, resulting in compromised voluntary muscular ability, is not entirely known. This review addresses the major advances that are being made in our understanding of the molecular mechanisms giving rise to the disease, which may eventually translate into new treatment options.

Local unfolding of Cu, Zn superoxide dismutase monomer determines the morphology of fibrillar aggregates

Aggregation of Cu, Zn superoxide dismutase (SOD1) is often found in amyotrophic lateral sclerosis patients. The fibrillar aggregates formed by wild type and various disease-associated mutants have recently been found to have distinct cores and morphologies. Previous computational and experimental studies of wild-type SOD1 suggest that the apo-monomer, highly aggregation prone, displays substantial local unfolding dynamics. The residual folded structure of locally unfolded apoSOD1 corresponds to peptide segments forming the aggregation core as identified by a combination of proteolysis and mass spectroscopy. Therefore, we hypothesize that the destabilization of apoSOD1 caused by various mutations leads to distinct local unfolding dynamics. The partially unfolded structure, exposing the hydrophobic core and backbone hydrogen bond donors and acceptors, is prone to aggregate. The peptide segments in the residual folded structures form the "building block" for aggregation, which in turn determines the morphology of the aggregates. To test this hypothesis, we apply a multiscale simulation approach to study the aggregation of three typical SOD1 variants: wild type, G37R, and I149T. Each of these SOD1 variants has distinct peptide segments forming the core structure and features different aggregate morphologies. We perform atomistic molecular dynamics simulations to study the conformational dynamics of apoSOD1 monomer and coarse-grained molecular dynamics simulations to study the aggregation of partially unfolded SOD1 monomers. Our computational studies of monomer local unfolding and the aggregation of different SOD1 variants are consistent with experiments, supporting the hypothesis of the formation of aggregation "building blocks" via apo-monomer local unfolding as the mechanism of SOD1 fibrillar aggregation.

Glutathionylation at Cys-111 induces dissociation of wild type and FALS mutant SOD1 dimers

Mutation of the ubiquitous cytosolic enzyme Cu/Zn superoxide dismutase (SOD1) is hypothesized to cause familial amyotrophic lateral sclerosis (FALS) through structural destabilization leading to misfolding and aggregation. Considering the late onset of symptoms as well as the phenotypic variability among patients with identical SOD1 mutations, it is clear that nongenetic factor(s) impact ALS etiology and disease progression. Here we examine the effect of Cys-111 glutathionylation, a physiologically prevalent post-translational oxidative modification, on the stabilities of wild type SOD1 and two phenotypically diverse FALS mutants, A4V and I112T. Glutathionylation results in profound destabilization of SOD1(WT) dimers, increasing the equilibrium dissociation constant K(d) to ~10-20 μM, comparable to that of the aggressive A4V mutant. SOD1(A4V) is further destabilized by glutathionylation, experiencing an ~30-fold increase in K(d). Dissociation kinetics of glutathionylated SOD1(WT) and SOD1(A4V) are unchanged, as measured by surface plasmon resonance, indicating that glutathionylation destabilizes these variants by decreasing association rate. In contrast, SOD1(I112T) has a modestly increased dissociation rate but no change in K(d) when glutathionylated. Using computational structural modeling, we show that the distinct effects of glutathionylation on different SOD1 variants correspond to changes in composition of the dimer interface. Our experimental and computational results show that Cys-111 glutathionylation induces structural rearrangements that modulate stability of both wild type and FALS mutant SOD1. The distinct sensitivities of SOD1 variants to glutathionylation, a modification that acts in part as a coping mechanism for oxidative stress, suggest a novel mode by which redox regulation and aggregation propensity interact in ALS.