Folding of Cu, Zn Superoxide Dismutase and Familial Amyotrophic Lateral Sclerosis

Biophysical and spectroscopical insights into structural modulation of species in the aggregation pathway of superoxide dismutase 1

Superoxide dismutase 1 (SOD1) aggregation is implicated in the development of Amyotrophic Lateral Sclerosis (ALS). Despite knowledge of the role of SOD1 aggregation, the mechanistic understanding remains elusive. Our investigation aimed to unravel the complex steps involved in SOD1 aggregation associated with ALS. Therefore, we probed the aggregation using ThT fluorescence, size-exclusion chromatography, and surface-enhanced Raman spectroscopy (SERS). The removal of metal ions and disulfide bonds resulted in the dimers rapidly first converting to an extended monomers then coming together slowly to form non-native dimers. The rapid onset of oligomerization happens above critical non-native dimer concentration. Structural features of oligomer was obtained through SERS. The kinetic data supported a fragmentation-dominant mechanism for the fibril formation. Quercetin acts as inhibitor by delaying the formation of non-native dimer and soluble oligomers by decreasing the elongation rate. Thus, results provide significant insights into the critical steps in oligomer formation and their structure.

Neurological Diseases can be Regulated by Phase Separation

Several neurological diseases arise from abnormal protein aggregation within neurones and this is closely regulated by phase separation. One such is motor neurone disease and aberrant aggregation of superoxide dismutase. Again these events are regulated by electrical forces that are examined.

Unveiling the double-edged sword: SOD1 trimers possess tissue-selective toxicity and bind septin-7 in motor neuron-like cells

Misfolded species of superoxide dismutase 1 (SOD1) are associated with increased death in amyotrophic lateral sclerosis (ALS) models compared to insoluble protein aggregates. The mechanism by which structurally independent SOD1 trimers cause cellular toxicity is unknown but may drive disease pathology. Here, we uncovered the SOD1 trimer interactome-a map of potential tissue-selective protein-binding partners in the brain, spinal cord, and skeletal muscle. We identified binding partners and key pathways associated with SOD1 trimers and found that trimers may affect normal cellular functions such as dendritic spine morphogenesis and synaptic function in the central nervous system and cellular metabolism in skeletal muscle. We discovered SOD1 trimer-selective enrichment of genes. We performed detailed computational and biochemical characterization of SOD1 trimer protein binding for septin-7. Our investigation highlights key proteins and pathways within distinct tissues, revealing a plausible intersection of genetic and pathophysiological mechanisms in ALS through interactions involving SOD1 trimers.

Big versus small: The impact of aggregate size in disease

Protein aggregation results in an array of different size soluble oligomers and larger insoluble fibrils. Insoluble fibrils were originally thought to cause neuronal cell deaths in neurodegenerative diseases due to their prevalence in tissue samples and disease models. Despite recent studies demonstrating the toxicity associated with soluble oligomers, many therapeutic strategies still focus on fibrils or consider all types of aggregates as one group. Oligomers and fibrils require different modeling and therapeutic strategies, targeting the toxic species is crucial for successful study and therapeutic development. Here, we review the role of different-size aggregates in disease, and how factors contributing to aggregation (mutations, metals, post-translational modifications, and lipid interactions) may promote oligomers opposed to fibrils. We review two different computational modeling strategies (molecular dynamics and kinetic modeling) and how they are used to model both oligomers and fibrils. Finally, we outline the current therapeutic strategies targeting aggregating proteins and their strengths and weaknesses for targeting oligomers versus fibrils. Altogether, we aim to highlight the importance of distinguishing the difference between oligomers and fibrils and determining which species is toxic when modeling and creating therapeutics for protein aggregation in disease.

Intrinsic structural vulnerability in the hydrophobic core induces species-specific aggregation of canine SOD1 with degenerative myelopathy-linked E40K mutation

Canine degenerative myelopathy (DM), a fatal neurodegenerative disease in dogs, shares clinical and genetic features with amyotrophic lateral sclerosis (ALS), a human motor neuron disease. Mutations in the SOD1 gene encoding Cu/Zn superoxide dismutase (SOD1) cause canine DM and a subset of inherited human ALS. The most frequent DM causative mutation is homozygous E40K mutation which induces the aggregation of canine SOD1 but not of human SOD1. However, the mechanism through which canine E40K mutation induces species-specific aggregation of SOD1 remains unknown. By screening human/canine chimeric SOD1s, we identified that the humanized mutation of the 117th residue (M117L), encoded by exon 4, significantly reduced aggregation propensity of canine SOD1^E40K. Conversely, introducing a mutation of leucine 117 to methionine, a residue homologous to canine, promoted E40K-dependent aggregation in human SOD1. M117L mutation improved protein stability and reduced cytotoxicity of canine SOD1^E40K. Furthermore, crystal structural analysis of canine SOD1 proteins revealed that M117L increased the packing within the hydrophobic core of the β-barrel structure, contributing to the increased protein stability. Our findings indicate that the structural vulnerability derived intrinsically from Met 117 in the hydrophobic core of the β-barrel structure induces E40K-dependent species-specific aggregation in canine SOD1.

A computational strategy for therapeutic development against superoxide dismutase (SOD1) amyloid formation: effect of polyphenols on the various events in the aggregation pathway

Pathology of superoxide dismutase 1 (SOD1) aggregation is linked to a neurodegenerative disease known as amyotrophic lateral sclerosis (ALS). Without suitable post-translational modifications (PTMs), the protein structure tends to become aggregation-prone. Understanding the role of PTMs and targeting the aggregation-prone SOD1 with small molecules can be used to design a strategy to inhibit its aggregation. Microsecond long molecular dynamics (MD) simulations followed by free energy surface (FES) analyses show that the loss of structure in the apo monomer happens locally and stepwise. Removing the disulfide bond from apoprotein leads to further instability in the zinc-binding loop, giving rise to non-native protein conformations. Further, it was found that these non-native conformations have a higher propensity to form a non-native dimer. We chose three structurally similar polyphenols based on their binding energies and investigated their impact on SOD1 aggregation kinetics. MD simulations of apo-SOD1^SH/corkscrew fibril-polyphenol complexes were also carried out. The effect of polyphenols was seen on fibril elongation as well. Based on the experiments and MD simulation results, it can be inferred that the choice of inhibitors is influenced not only by the binding energy but also by dimer interface stabilization, the proclivity to form non-native dimers, the propensity to break fibrils, and the propensity to decrease the rate of elongation. The polyphenols with 3' and 4' hydroxyl groups are better inhibitors of SOD1 aggregation.

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

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

Quercetin and Baicalein Act as Potent Antiamyloidogenic and Fibril Destabilizing Agents for SOD1 Fibrils

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease that has been associated with the deposition of aggregates of superoxide dismutase 1 (SOD1). Effective therapeutics against SOD1 fibrillation is still an area of active research. Herein, we demonstrate the potential of two naturally occurring flavonoids (quercetin and baicalein) to inhibit fibrillation of wild-type SOD1 with the aid of a series of biophysical techniques. Our seeding experiments reveal that both of these flavonoids significantly affect the fibril elongation. Interestingly, our ThT binding assay, TEM, and SDS-PAGE experiments suggest that these flavonoids also disintegrate the fibrils into shorter fragments but do not completely depolymerize them into monomers. Binding parameters obtained from the analysis of UV-vis spectra suggest that these flavonoids bind moderately to native SOD1 dimer and have different binding sites. Docking of these flavonoids with a non-native monomer, non-native trimer, and oligomer derived from the 11-residue segment of SOD1 indicates that both quercetin and baicalein can bind to these species and thus can arrest the elongation of fibrils by blocking the fibrillar core regions on the intermediate species formed during aggregation of SOD1. MTT assay data revealed that both the flavonoids reduced the cytotoxicity of SOD1 fibrils. Experimental data also show the antiamyloidogenic potential of both flavonoids against A4V SOD1 mutant fibrillation. Thus, our findings may provide a direction for designing effective therapeutic agents against ALS which can act as promising antiamyloidogenic and fibril destabilizing agents.

The Role of Data in Model Building and Prediction: A Survey Through Examples

The goal of Science is to understand phenomena and systems in order to predict their development and gain control over them. In the scientific process of knowledge elaboration, a crucial role is played by models which, in the language of quantitative sciences, mean abstract mathematical or algorithmical representations. This short review discusses a few key examples from Physics, taken from dynamical systems theory, biophysics, and statistical mechanics, representing three paradigmatic procedures to build models and predictions from available data. In the case of dynamical systems we show how predictions can be obtained in a virtually model-free framework using the methods of analogues, and we briefly discuss other approaches based on machine learning methods. In cases where the complexity of systems is challenging, like in biophysics, we stress the necessity to include part of the empirical knowledge in the models to gain the minimal amount of realism. Finally, we consider many body systems where many (temporal or spatial) scales are at play—and show how to derive from data a dimensional reduction in terms of a Langevin dynamics for their slow components.

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.

Predictive methods for computational metalloenzyme redesign - a test case with carboxypeptidase A

Computational metalloenzyme design is a multi-scale problem. It requires treating the metal coordination quantum mechanically, extensive sampling of the protein backbone, and additionally accounting for the polarization of the active site by both the metal cation and the surrounding protein (a phenomenon called electrostatic preorganization). We bring together a combination of theoretical methods that jointly offer these desired qualities: QM/DMD for mixed quantum-classical dynamic sampling, quantum theory of atoms in molecules (QTAIM) for the assessment of electrostatic preorganization, and Density Functional Theory (DFT) for mechanistic studies. Within this suite of principally different methods, there are both complementarity of capabilities and cross-validation. Using these methods, predictions can be made regarding the relative activities of related enzymes, as we show on the native Zn²⁺-dependent carboxypeptidase A (CPA), and its mutant proteins, which are hypothesized to hydrolyze modified substrates. For the native CPA, we replicated the catalytic mechanism and the rate in close agreement with the experiment, giving validity to the QM/DMD predicted structure, the DFT mechanism, and the QTAIM assessment of catalytic activity. For most sequences of the modified substrate and tried CPA mutants, substantially worsened activity is predicted. However, for the substrate mutant that contains Asp instead of Phe at the C-terminus, one CPA mutant exhibits a reasonable activity, as predicted across the theoretical methods. CPA is a well-studied system, and here it serves as a testing ground for the offered methods.

Packing energetics determine the folding routes of the RNase-H proteins

The RNase-H proteins show a diverse range of folding routes with structurally distinct folding nuclei.

Stabilizing Off-pathway Oligomers by Polyphenol Nanoassemblies for IAPP Aggregation Inhibition

Experimental studies have shown that many naturally occurring polyphenols have inhibitory effect on the aggregation of several proteins. Here, we use discrete molecular dynamics (DMD) simulations and high-throughput dynamic light scattering (DLS) experiments to study the anti-aggregation effects of two polyphenols, curcumin and resveratrol, on the aggregation of islet amyloid polypeptide (IAPP or amylin). Our DMD simulations suggest that the aggregation inhibition is caused by stabilization of small molecular weight IAPP off-pathway oligomers by the polyphenols. Our analysis indicates that IAPP-polyphenol hydrogen bonds and π-π stacking combined with hydrophobic interactions are responsible for the stabilization of oligomers. The presence of small oligomers is confirmed with DLS measurements in which nanometer-sized oligomers are found to be stable for up to 7.5 hours, the time frame within which IAPP aggregates in the absence of polyphenols. Our study offers a general anti-aggregation mechanism for polyphenols, and further provides a computational framework for the future design of anti-amyloid aggregation therapeutics.

Functionally Relevant Specific Packing Can Determine Protein Folding Routes

Functional residues can modulate the folding mechanisms of proteins. In some proteins, mutations to such residues can radically change the primary folding route. Is it possible then to learn more about the functional regions of a protein by investigating just its choice of folding route? The folding and the function of the protein Escherichia coli ribonuclease H (ecoRNase-H) have been extensively studied and its folding route is known to near-residue resolution. Here, we computationally study the folding of ecoRNase-H using molecular dynamics simulations of structure-based models of increasing complexity. The differences between a model that correctly predicts the experimentally determined folding route and a simpler model that does not can be attributed to a set of six aromatic residues clustered together in a region of the protein called CORE. This clustering, which we term "specific" packing, drives CORE to fold early and determines the folding route. Both the residues involved in specific packing and their packing are largely conserved across E. coli-like RNase-Hs from diverse species. Residue conservation is usually implicated in function. Here, the identified residues either are known to bind substrate in ecoRNase-H or pack against the substrate in the homologous human RNase-H where a substrate-bound crystal structure exists. Thus, the folding mechanism of ecoRNase-H is a byproduct of functional demands upon its sequence. Using our observations on specific packing, we suggest mutations to an engineered HIV RNase-H to make its function better. Our results show that understanding folding route choice in proteins can provide unexpected insights into their function.

Competition-cooperation relationship networks characterize the competition and cooperation between proteins

By analyzing protein-protein interaction (PPI) networks, one can find that a protein may have multiple binding partners. However, it is difficult to determine whether the interactions with these partners occur simultaneously from binary PPIs alone. Here, we construct the yeast and human competition-cooperation relationship networks (CCRNs) based on protein structural interactomes to clearly exhibit the relationship (competition or cooperation) between two partners of the same protein. If two partners compete for the same interaction interface, they would be connected by a competitive edge; otherwise, they would be connected by a cooperative edge. The properties of three kinds of hubs (i.e., competitive, modest, and cooperative hubs) are analyzed in the CCRNs. Our results show that competitive hubs have higher clustering coefficients and form clusters in the human CCRN, but these tendencies are not observed in the yeast CCRN. We find that the human-specific proteins contribute significantly to these differences. Subsequently, we conduct a series of computational experiments to investigate the regulatory mechanisms that avoid competition between proteins. Our comprehensive analyses reveal that for most yeast and human protein competitors, transcriptional regulation plays an important role. Moreover, the human-specific proteins have a particular preference for other regulatory mechanisms, such as alternative splicing.

Calcium binding to gatekeeper residues flanking aggregation-prone segments underlies non-fibrillar amyloid traits in superoxide dismutase 1 (SOD1)

Calcium deregulation is a central feature among neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Calcium accumulates in the spinal and brain stem motor neurons of ALS patients triggering multiple pathophysiological processes which have been recently shown to include direct effects on the aggregation cascade of superoxide dismutase 1 (SOD1). SOD1 is a Cu/Zn enzyme whose demetallated form is implicated in ALS protein deposits, contributing to toxic gain of function phenotypes. Here we undertake a combined experimental and computational study aimed at establishing the molecular details underlying the regulatory effects of Ca(2+) over SOD1 aggregation potential. Isothermal titration calorimetry indicates entropy driven low affinity association of Ca(2+) ions to apo SOD1, at pH7.5 and 37°C. Molecular dynamics simulations denote a noticeable loss of native structure upon Ca(2+) association that is especially prominent at the zinc-binding and electrostatic loops, whose decoupling is known to expose the central SOD1 β-barrel triggering aggregation. Structural mapping of the preferential apo SOD1 Ca(2+) binding locations reveals that among the most frequent ligands for Ca(2+) are negatively-charged gatekeeper residues located in boundary positions with respect to segments highly prone to edge-to-edge aggregation. Calcium interactions thus diminish gatekeeping roles of these residues, by shielding repulsive interactions via stacking between aggregating β-sheets, partly blocking fibril formation and promoting amyloidogenic oligomers such as those found in ALS inclusions. Interestingly, many fALS mutations occur at these positions, disclosing how Ca(2+) interactions recreate effects similar to those of genetic defects, a finding with relevance to understand sporadic ALS pathomechanisms.

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.

Cu,Zn-superoxide dismutase without Zn is folded but catalytically inactive

Amyotrophic lateral sclerosis has been linked to the gain of aberrant function of superoxide dismutase, Cu,Zn-SOD1 upon protein misfolding. The mechanism of SOD1 misfolding is thought to involve mutations leading to the loss of Zn, followed by protein unfolding and aggregation. We show that the removal of Zn from SOD1 may not lead to an immediate unfolding but immediately deactivates the enzyme through a combination of subtle structural and electronic effects. Using quantum mechanics/discrete molecular dynamics, we showed that both Zn-less wild-type (WT)-SOD1 and its D124N mutant that does not bind Zn have at least metastable folded states. In those states, the reduction potential of Cu increases, leading to the presence of detectable amounts of Cu(I) instead of Cu(II) in the active site, as confirmed experimentally. The Cu(I) protein cannot participate in the catalytic Cu(I)-Cu(II) cycle. However, even without the full reduction to Cu(I), the Cu site in the Zn-less variants of SOD1 is shown to be catalytically incompetent: unable to bind superoxide in a way comparable to the WT-SOD1. The changes are more radical and different in the D124N Zn-less mutant than in the Zn-less WT-SOD1, suggesting D124N being perhaps not the most adequate model for Zn-less SOD1. Overall, Zn in SOD1 appears to be influencing the Cu site directly by adjusting its reduction potential and geometry. Thus, the role of Zn in SOD1 is not just structural, as was previously thought; it is a vital part of the catalytic machinery.

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.

Structural propensities of human ubiquitination sites: accessibility, centrality and local conformation

The existence and function of most proteins in the human proteome are regulated by the ubiquitination process. To date, tens of thousands human ubiquitination sites have been identified from high-throughput proteomic studies. However, the mechanism of ubiquitination site selection remains elusive because of the complicated sequence pattern flanking the ubiquitination sites. In this study, we perform a systematic analysis of 1,330 ubiquitination sites in 505 protein structures and quantify the significantly high accessibility and unexpectedly high centrality of human ubiquitination sites. Further analysis suggests that the higher centrality of ubiquitination sites is associated with the multi-functionality of ubiquitination sites, among which protein-protein interaction sites are common targets of ubiquitination. Moreover, we demonstrate that ubiquitination sites are flanked by residues with non-random local conformation. Finally, we provide quantitative and unambiguous evidence that most of the structural propensities contain specific information about ubiquitination site selection that is not represented by the sequence pattern. Therefore, the hypothesis about the structural level of the ubiquitination site selection mechanism has been substantially approved.

A quantitative measure of electrostatic perturbation in holo and apo enzymes induced by structural changes

Biological pathways are subject to subtle manipulations that achieve a wide range of functional variation in differing physiological niches. In many instances, changes in the structure of an enzyme on ligand binding germinate electrostatic perturbations that form the basis of its changed catalytic or transcriptional efficiency. Computational methods that seek to gain insights into the electrostatic changes in enzymes require expertise to setup and computing prowess. In the current work, we present a fast, easy and reliable methodology to compute electrostatic perturbations induced by ligand binding (MEPP). The theoretical foundation of MEPP is the conserved electrostatic potential difference (EPD) in cognate pairs of active site residues in proteins with the same functionality. Previously, this invariance has been used to unravel promiscuous serine protease and metallo-β-lactamase scaffolds in alkaline phosphatases. Given that a similarity in EPD is significant, we expect differences in the EPD to be significant too. MEPP identifies residues or domains that undergo significant electrostatic perturbations, and also enumerates residue pairs that undergo significant polarity change. The gain in a certain polarity of a residue with respect to neighboring residues, or the reversal of polarity between two residues might indicate a change in the preferred ligand. The methodology of MEPP has been demonstrated on several enzymes that employ varying mechanisms to perform their roles. For example, we have attributed the change in polarity in residue pairs to be responsible for the loss of metal ion binding in fructose 1,6-bisphosphatases, and corroborated the pre-organized state of the active site of the enzyme with respect to functionally relevant changes in electric fields in ketosteroid isomerases.

Redox properties of the disulfide bond of human Cu,Zn superoxide dismutase and the effects of human glutaredoxin 1

The intramolecular disulfide bond in hSOD1 [human SOD1 (Cu,Zn superoxide dismutase 1)] plays a key role in maintaining the protein's stability and quaternary structure. In mutant forms of SOD1 that cause familial ALS (amyotrophic lateral sclerosis), this disulfide bond is more susceptible to chemical reduction, which may lead to destabilization of the dimer and aggregation. During hSOD1 maturation, disulfide formation is catalysed by CCS1 (copper chaperone for SOD1). Previous studies in yeast demonstrate that the yeast GSH/Grx (glutaredoxin) redox system promotes reduction of the hSOD1 disulfide in the absence of CCS1. In the present study, we probe further the interaction between hSOD1, GSH and Grxs to provide mechanistic insight into the redox kinetics and thermodynamics of the hSOD1 disulfide. We demonstrate that hGrx1 (human Grx1) uses a monothiol mechanism to reduce the hSOD1 disulfide, and the GSH/hGrx1 system reduces ALS mutant SOD1 at a faster rate than WT (wild-type) hSOD1. However, redox potential measurements demonstrate that the thermodynamic stability of the disulfide is not consistently lower in ALS mutants compared with WT hSOD1. Furthermore, the presence of metal cofactors does not influence the disulfide redox potential. Overall, these studies suggest that differences in the GSH/hGrx1 reaction rate with WT compared with ALS mutant hSOD1 and not the inherent thermodynamic stability of the hSOD1 disulfide bond may contribute to the greater pathogenicity of ALS mutant hSOD1.

The complex molecular biology of amyotrophic lateral sclerosis (ALS)

Amyotrophic lateral sclerosis (ALS) is an adult-onset neurodegenerative disorder that causes selective death of motor neurons followed by paralysis and death. A subset of ALS cases is caused by mutations in the gene for Cu, Zn superoxide dismutase (SOD1), which impart a toxic gain of function to this antioxidant enzyme. This neurotoxic property is widely believed to stem from an increased propensity to misfold and aggregate caused by decreased stability of the native homodimer or a tendency to lose stabilizing posttranslational modifications. Study of the molecular mechanisms of SOD1-related ALS has revealed a complex array of interconnected pathological processes, including glutamate excitotoxicity, dysregulation of neurotrophic factors and axon guidance proteins, axonal transport defects, mitochondrial dysfunction, deficient protein quality control, and aberrant RNA processing. Many of these pathologies are directly exacerbated by misfolded and aggregated SOD1 and/or cytosolic calcium overload, suggesting the primacy of these events in disease etiology and their potential as targets for therapeutic intervention.

Discrete molecular dynamics: an efficient and versatile simulation method for fine protein characterization

Until now it has been impractical to observe protein folding in silico for proteins larger than 50 residues. Limitations of both force field accuracy and computational efficiency make the folding problem very challenging. Here we employ discrete molecular dynamics (DMD) simulations with an all-atom force field to fold fast-folding proteins. We extend the DMD force field by introducing long-range electrostatic interactions to model salt-bridges and a sequence-dependent semiempirical potential accounting for natural tendencies of certain amino acid sequences to form specific secondary structures. We enhance the computational performance by parallelizing the DMD algorithm. Using a small number of commodity computers, we achieve sampling quality and folding accuracy comparable to the explicit-solvent simulations performed on high-end hardware. We demonstrate that DMD can be used to observe equilibrium folding of villin headpiece and WW domain, study two-state folding kinetics, and sample near-native states in ab initio folding of proteins of ∼100 residues.

Engineering a thermo-stable superoxide dismutase functional at sub-zero to >50°C, which also tolerates autoclaving

Superoxide dismutase (SOD) is a critical enzyme associated with controlling oxygen toxicity arising out of oxidative stress in any living system. A hyper-thermostable SOD isolated from a polyextremophile higher plant Potentilla atrosanguinea Lodd. var. argyrophylla (Wall. ex Lehm.) was engineered by mutation of a single amino acid that enhanced the thermostability of the enzyme to twofold. The engineered enzyme was functional from sub-zero temperature to >50°C, tolerated autoclaving (heating at 121°C, at a pressure of 1.1 kg per square cm for 20 min) and was resistant to proteolysis. The present work is the first example to enhance the thermostability of a hyper-thermostable protein and has potential to application to other proteins for enhancing thermostability.

Coarse-grained simulations of protein aggregation

Protein aggregation is believed to be responsible for a number of human diseases and limits the yields of pharmaceutical proteins during production. Computer simulations can be used to develop novel experimentally testable hypotheses pertaining to aggregation. While all-atom simulations with explicit solvent are too computationally intensive to address the multitude of relevant time scales, coarse-grained models make it possible to observe the transition of monomers to an equilibrium containing aggregates. Here, we provide the reader with background information and a list of steps for setting up, performing, and analyzing computer simulations of aggregating coarse-grained (CG) proteins.

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.

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

Aggregation of Cu,Zn superoxide dismutase (SOD1) is implicated in amyotrophic lateral sclerosis. Glutathionylation and phosphorylation of SOD1 is omnipresent in the human body, even in healthy individuals, and has been shown to increase SOD1 dimer dissociation, which is the first step on the pathway toward SOD1 aggregation. We found that post-translational modification of SOD1, especially glutathionylation, promotes dimer dissociation. We discovered an intermediate state in the pathway to dissociation, a conformational change that involves a "loosening" of the β-barrels and a loss or shift of dimer interface interactions. In modified SOD1, this intermediate state is stabilized as compared to unmodified SOD1. The presence of post-translational modifications could explain the environmental factors involved in the speed of disease progression. Because post-translational modifications such as glutathionylation are often induced by oxidative stress, post-translational modification of SOD1 could be a factor in the occurrence of sporadic cases of amyotrophic lateral sclerosis, which represent 90% of all cases of the disease.

Functional aspects of redox control during neuroinflammation

Neuroinflammation is a CNS reaction to injury in which some severe pathologies, regardless of their origin, converge. The phenomenon emphasizes crosstalk between neurons and glia and reveals a complex interaction with oxidizing agents through redox sensors localized in enzymes, receptors, and transcription factors. When oxidizing pressures cause reversible molecular changes, such as minimal or transitory proinflammatory cytokine overproduction, redox couples provide a means of translating the presence of reactive oxygen or nitrogen species into useful signals in the cell. Additionally, thiol-based redox sensors convey information about localized changes in redox potential induced by physiologic or pathologic situations. They are susceptible to oxidative changes and become key events during neuroinflammation, altering the course of a signaling response or the behavior of specific transcription factors. When oxidative stress augments the pressure on the intracellular environment, the effective reduction potential of redox pairs diminishes, and cell signaling shifts toward proinflammatory and proapoptotic signals, creating a vicious cycle between oxidative stress and neuroinflammation. In addition, electrophilic compounds derived from the oxidative cascade react with key protein thiols and interfere with redox signaling. This article reviews the relevant functional aspects of redox control during the neuroinflammatory process.

Lymphotactin: how a protein can adopt two folds

Metamorphic proteins such as lymphotactin are a notable exception of the empirical principle that structured natural proteins possess a unique three-dimensional structure. In particular, the human chemokine lymphotactin protein exists in two distinct conformations (one monomeric and one dimeric) under physiological conditions. In this work, we use a C(alpha) Go model to show how this very peculiar behavior can be reproduced. From the study of the thermodynamics and of the kinetics, we characterize the interconversion mechanism. In particular, this takes place through the docking of the two chains living in a third monomeric, partially unfolded, state which shows a residual structure involving a set of local contacts common to the two native conformations. The main feature of two fold proteins appears to be the sharing of a common set of local contacts between the two distinct folds as confirmed by the study of two designed two fold proteins. Metamorphic proteins may be more common than expected.

Structural changes to monomeric CuZn superoxide dismutase caused by the familial amyotrophic lateral sclerosis-associated mutation A4V

Amyotrophic lateral sclerosis (ALS) is a progressive motor neuron degenerative disease, and the inherited form, familial ALS (fALS), has been linked to over 100 different point mutations scattered throughout the Cu-Zn superoxide dismutase protein (SOD1). The disease is likely due to a toxic gain of function caused by the misfolding, oligomerization, and eventual aggregation of mutant SOD1, but it is not yet understood how the structurally diverse mutations result in a common disease phenotype. The behavior of the apo-monomer fALS-associated mutant protein A4V was explored using molecular-dynamics simulations to elucidate characteristic structural changes to the protein that may allow the mutant form to improperly associate with other monomer subunits. Simulations showed that the mutant protein is less stable than the WT protein overall, with shifts in residue-residue contacts that lead to destabilization of the dimer and metal-binding sites, and stabilization of nonnative contacts that leads to a misfolded state. These findings provide a unifying explanation for disparate experimental observations, allow a better understanding of alterations of residue contacts that accompany loss of SOD1 structural integrity, and suggest sites where compensatory changes may stabilize the mutant structure.

Multiscale approaches for studying energy transduction in dynein

Cytoplasmic dynein is an important motor that drives all minus-end directed movement along microtubules. Dynein is a complex motor whose processive motion is driven by ATP-hydrolysis. Dynein's run length has been measured to be several millimetres with typical velocities in the order of a few nanometres per second. Therefore, the average time between steps is a fraction of a second. When this time scale is compared with typical time scales for protein side chain and backbone movements (approximately 10(-9) s and approximately 10(-5) s, respectively), it becomes clear that a multi-timescale modelling approach is required to understand energy transduction in this protein. Here, we review recent efforts to use computational and mathematical modelling to understand various aspects of dynein's chemomechanical cycle. First, we describe a structural model of dynein's motor unit showing a heptameric organization of the motor subunits. Second, we describe our molecular dynamics simulations of the motor unit that are used to investigate the dynamics of the various motor domains. Third, we present a kinetic model of the coordination between the two dynein heads. Lastly, we investigate the various potential geometries of the dimer during its hydrolytic and stepping cycle.

Modifications of superoxide dismutase (SOD1) in human erythrocytes: a possible role in amyotrophic lateral sclerosis

Over 100 mutations in Cu/Zn-superoxide dismutase (SOD1) result in familial amyotrophic lateral sclerosis. Dimer dissociation is the first step in SOD1 aggregation, and studies suggest nearly every amino acid residue in SOD1 is dynamically connected to the dimer interface. Post-translational modifications of SOD1 residues might be expected to have similar effects to mutations, but few modifications have been identified. Here we show, using SOD1 isolated from human erythrocytes, that human SOD1 is phosphorylated at threonine 2 and glutathionylated at cysteine 111. A second SOD1 phosphorylation was observed and mapped to either Thr-58 or Ser-59. Cysteine 111 glutathionylation promotes SOD1 monomer formation, a necessary initiating step in SOD1 aggregation, by causing a 2-fold increase in the K(d). This change in the dimer stability is expected to result in a 67% increase in monomer concentration, 315 nm rather than 212 nm at physiological SOD1 concentrations. Because protein glutathionylation is associated with redox regulation, our finding that glutathionylation promotes SOD1 monomer formation supports a model in which increased oxidative stress promotes SOD1 aggregation.

Insights from coarse-grained Gō models for protein folding and dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.

Dynamical roles of metal ions and the disulfide bond in Cu, Zn superoxide dismutase folding and aggregation

Misfolding and aggregation of Cu, Zn superoxide dismutase (SOD1) is implicated in neuronal death in amyotrophic lateral sclerosis. Each SOD1 monomer binds to 1 copper and 1 zinc ion and maintains its disulfide bond (Cys-57-Cys-146) in the reducing cytoplasm of cell. Mounting experimental evidence suggests that metal loss and/or disulfide reduction are important for initiating misfolding and aggregation of SOD1. To uncover the role of metals and the disulfide bond in the SOD1 folding, we systemically study the folding thermodynamics and structural dynamics of SOD1 monomer and dimer with and without metal binding and under disulfide-intact or disulfide-reduced environments in computational simulations. We use all-atom discrete molecular dynamics for sampling. Our simulation results provide dynamical evidence to the stabilizing role of metal ions in both dimer and monomer SOD1. The disulfide bond anchors a loop (Glu-49 to Asn-53) that contributes to the dimer interface. The reduction of the disulfide bond in SOD1 with metal ions depleted results in a flexible Glu-49-Asn-53 loop, which, in turn, disrupts dimer formation. Interestingly, the disulfide bond reduction does not affect the thermostability of monomer SOD1 as significantly as the metal ions do. We further study the structural dynamics of metal-free SOD1 monomers, the precursor for aggregation, in simulations and find inhomogeneous local unfolding of beta-strands. The strands protected by the metal-binding and electrostatic loops are found to unfold first after metal loss, leading to a partially unfolded beta-sheet prone to aggregation. Our simulation study sheds light on the critical role of metals and disulfide bond in SOD1 folding and aggregation.

Probing protein aggregation using discrete molecular dynamics

Understanding the role of biomolecular dynamics in cellular processes leading to human diseases and the ability to rationally manipulate these processes is of fundamental importance in scientific research. The last decade has witnessed significant progress in probing biophysical behavior of proteins. However, we are still limited in understanding how changes in protein dynamics and inter-protein interactions occurring in short length- and time-scales lead to aberrations in their biological function. Bridging this gap in biology probed using computer simulations marks a challenging frontier in computational biology. Here we examine hypothesis-driven simplified protein models in conjunction with discrete molecular dynamics in the study of protein aggregation, implicated in series of neurodegenerative diseases, such as Alzheimer's and Huntington's diseases. Discrete molecular dynamics simulations of simplified protein models have emerged as a powerful methodology with its ability to bridge the gap in time and length scales from protein dynamics to aggregation, and provide an indispensable tool for probing protein aggregation.

Diminished self-chaperoning activity of the DeltaF508 mutant of CFTR results in protein misfolding

The absence of a functional ATP Binding Cassette (ABC) protein called the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) from apical membranes of epithelial cells is responsible for cystic fibrosis (CF). Over 90% of CF patients carry at least one mutant allele with deletion of phenylalanine at position 508 located in the N-terminal nucleotide binding domain (NBD1). Biochemical and cell biological studies show that the ΔF508 mutant exhibits inefficient biosynthetic maturation and susceptibility to degradation probably due to misfolding of NBD1 and the resultant misassembly of other domains. However, little is known about the direct effect of the Phe508 deletion on the NBD1 folding, which is essential for rational design strategies of cystic fibrosis treatment. Here we show that the deletion of Phe508 alters the folding dynamics and kinetics of NBD1, thus possibly affecting the assembly of the complete CFTR. Using molecular dynamics simulations, we find that meta-stable intermediate states appearing on wild type and mutant folding pathways are populated differently and that their kinetic accessibilities are distinct. The structural basis of the increased misfolding propensity of the ΔF508 NBD1 mutant is the perturbation of interactions in residue pairs Q493/P574 and F575/F578 found in loop S7-H6. As a proof-of-principle that the S7-H6 loop conformation can modulate the folding kinetics of NBD1, we virtually design rescue mutations in the identified critical interactions to force the S7-H6 loop into the wild type conformation. Two redesigned NBD1-ΔF508 variants exhibited significantly higher folding probabilities than the original NBD1-ΔF508, thereby partially rescuing folding ability of the NBD1-ΔF508 mutant. We propose that these observed defects in folding kinetics of mutant NBD1 may also be modulated by structures separate from the 508 site. The identified structural determinants of increased misfolding propensity of NBD1-ΔF508 are essential information in correcting this pathogenic mutant.

Deletion of a single residue, phenylalanine at position 508, in the first nucleotide binding domain (NBD1) of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is present in approximately 90% of cystic fibrosis (CF) patients. Experiments show that this mutant protein exhibits inefficient biosynthetic maturation and susceptibility to degradation probably due to misfolding of NBD1 and the resultant incorrect interactions of other domains. However, little is known about the direct effect of the Phe508 deletion on NBD1 folding. Here, using molecular dynamics simulations of NBD1-WT, NBD1-F508A, and NBD1-ΔF508, we show that the deletion of Phe508 indeed alters the kinetics of NBD1 folding. We also find that the intermediate states appearing on wild type and mutant folding pathways are populated differently and that their kinetic accessibilities are distinct. Moreover, we identified critical interactions not necessarily localized near position 508, such as Q493/P574 and F575/F587, to be significant structural elements influencing the kinetic difference between wild type and mutant NBD1. We propose that these observed alterations in folding kinetics of mutant NBD1 result in misassembly of the whole multi-domain protein, thereby causing its premature degradation.

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.

Deciphering protein dynamics from NMR data using explicit structure sampling and selection

Perhaps one of the most prominent realizations of recent years is the critical role that protein dynamics plays in many facets of cellular function. While characterization of protein dynamics is fundamental to our understanding of protein function, the ability to explicitly detect an ensemble of protein conformations from dynamics data is a paramount challenge in structural biology. Here, we report a new computational method, Sample and Select, for determining the ensemble of protein conformations consistent with NMR dynamics data. This method can be generalized and extended to different sources of dynamics data, enabling broad applicability in deciphering protein dynamics at different timescales. The structural ensemble derived from Sample and Select will provide structural and dynamic information that should aid us in understanding and manipulating protein function.

Multiscale modeling of nucleosome dynamics

Nucleosomes form the fundamental building blocks of chromatin. Subtle modifications of the constituent histone tails mediate chromatin stability and regulate gene expression. For this reason, it is important to understand structural dynamics of nucleosomes at atomic levels. We report a novel multiscale model of the fundamental chromatin unit, a nucleosome, using a simplified model for rapid discrete molecular dynamics simulations and an all-atom model for detailed structural investigation. Using a simplified structural model, we perform equilibrium simulations of a single nucleosome at various temperatures. We further reconstruct all-atom nucleosome structures from simulation trajectories. We find that histone tails bind to nucleosomal DNA via strong salt-bridge interactions over a wide range of temperatures, suggesting a mechanism of chromatin structural organization whereby histone tails regulate inter- and intranucleosomal assemblies via binding with nucleosomal DNA. We identify specific regions of the histone core H2A/H2B-H4/H3-H3/H4-H2B/H2A, termed "cold sites", which retain a significant fraction of contacts with adjoining residues throughout the simulation, indicating their functional role in nucleosome organization. Cold sites are clustered around H3-H3, H2A-H4 and H4-H2A interhistone interfaces, indicating the necessity of these contacts for nucleosome stability. Essential dynamics analysis of simulation trajectories shows that bending across the H3-H3 is a prominent mode of intranucleosomal dynamics. We postulate that effects of salts on mononucleosomes can be modeled in discrete molecular dynamics by modulating histone-DNA interaction potentials. Local fluctuations in nucleosomal DNA vary significantly along the DNA sequence, suggesting that only a fraction of histone-DNA contacts make strong interactions dominating mononucleosomal dynamics. Our findings suggest that histone tails have a direct functional role in stabilizing higher-order chromatin structure, mediated by salt-bridge interactions with adjacent DNA.

Reversible pH-Driven Conformational Switching of Tethered Superoxide Dismutase with Gold Nanoparticle Enhanced Surface Plasmon Resonance Spectroscopy

A new class of surface-immobilized protein nanomachines can be reversibly actuated by cycling the solution pH between 2.5 and 12.3, which induces a conformational change, thereby modulating the thickness of superoxide dismutase (SOD1) tethered to the Au thin film. By placing Au nanoparticles (AuNP) atop the immobilized SOD1 by means of a gold-thiol assembly, the nanoscale motion of SOD1 at the interface produces mechanical work to lift and then lower the AuNP from the Au substrate by a distance of ca. 3 nm and transduces this motion into an easily measurable reflectivity change in the surface plasmon resonance (SPR) spectrum. As-made supported conjugate consisting of SOD1 and AuNP is quite robust and stable, and its operation in response to pH variations, which mirrors the conformational changes of responsive SOD1 at the interface, is found to be highly reversible and reproducible. This is the first demonstration of the development of novel solid-state sensors and/or switching devices based on substrate-bound protein conformational changes and AuNP enhanced SPR spectroscopy.

Common dynamical signatures of familial amyotrophic lateral sclerosis-associated structurally diverse Cu, Zn superoxide dismutase mutants

More than 100 structurally diverse point mutations leading to aggregation in the dimeric enzyme Cu, Zn superoxide dismutase (SOD1) are implicated in familial amyotrophic lateral sclerosis (FALS). Although SOD1 dimer dissociation is a known requirement for its aggregation, the common structural basis for diverse FALS mutations resulting in aggregation is not fully understood. In molecular dynamics simulations of wild-type SOD1 and three structurally diverse FALS mutants (A4V, G37R, and H46R), we find that a common effect of mutations on SOD1 dimer is the mutation-induced disruption of dynamic coupling between monomers. In the wild-type dimer, the principal coupled motion corresponds to a "breathing motion" of the monomers around an axis parallel to the dimer interface, and an opening-closing motion of the distal metal-binding loops. These coupled motions are disrupted in all three mutants independent of the mutation location. Loss of coupled motions in mutant dimers occurs with increased disruption of a key stabilizing structural element (the beta-plug) leading to the de-protection of edge strands. To rationalize disruption of coupling, which is independent of the effect of the mutation on global SOD1 stability, we analyze the residue-residue interaction network formed in SOD1. We find that the dimer interface and metal-binding loops, both involved in coupled motions, are regions of high connectivity in the network. Our results suggest that independent of the effect on protein stability, altered protein dynamics, due to long-range communication within its structure, may underlie the aggregation of mutant SOD1 in FALS.

Intracellular conformational alterations of mutant SOD1 and the implications for fALS-associated SOD1 mutant induced motor neuron cell death

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by the selective death of motor neurons. Approximately 10% of ALS cases are familial (fALS) and about 25% of fALS patients inherit autosomal dominant mutations in the gene encoding copper-zinc superoxide dismutase (SOD1). Over 90 different SOD1 mutations have been identified in fALS patients. It has been established that the ALS-linked SOD1 mutations provoke a new toxic function, the nature of which remains unclear. In vitro studies using various biophysical techniques have demonstrated that the SOD1 mutants share a reduced conformational stability. However, conformational alterations of the ALS mutants have not been directly demonstrated in vivo. We employed an SOD1-GFP fusion protein system in this study to monitor the intracellular protein conformation. We demonstrate that the ALS-linked SOD1 mutants adopt different conformations from the wild-type (WT) protein in living cells. Moreover, the conformational alterations of mutant SOD1 render the mutants susceptible to the formation of high-molecular-weight complexes prior to the appearance of detergent-resistant aggregates. Finally, we show that the motor neuron-like cells expressing mutant SOD1 are more susceptible to H2O2 induced cell death compared to the cells expressing WT SOD1. This study provides direct evidence of in vivo conformational differences between WT and mutant SOD1. In addition, the SOD1-GFP system can be exploited in future studies to investigate how conformational alterations of mutant SOD1 lead to protein aggregation and to study the potential toxicity of such aggregates in familial ALS.

Simple but predictive protein models

The traditional approach to computational biophysics studies of molecular systems is brute force molecular dynamics simulations under the conditions of interest. The disadvantages of this approach are that the time and length scales that are accessible to computer simulations often do not reach biologically relevant scales. An alternative approach, which we call intuitive modeling, is hypothesis-driven and based on tailoring simplified protein models to the systems of interest. Using intuitive modeling, the length and time scales that can be achieved using simplified protein models exceed those of traditional molecular-dynamic simulations. Here, we describe several recent studies that signify the predictive power of simplified protein models within the intuitive-modeling approach.

A Single Disulfide Bond Differentiates Aggregation Pathways of ß2-Microglobulin

Deposition of wild-type beta2-microglobulin (beta2m) into amyloid fibrils is a complication in patients undergoing long-term hemodialysis. The native beta-sandwich fold of beta2m has a highly conserved disulfide bond linking Cys25 and Cys80. Oxidized beta2m forms needle-like amyloid fibrils at pH 2.5 in vitro, whereas reduced beta2m, at acid pH, in which the intra-chain disulfide bond is disrupted, cannot form typical fibrils. Instead, reduced beta2m forms thinner and more flexible filaments. To uncover the difference in molecular mechanisms underlying the aggregation of the oxidized and reduced beta2m, we performed molecular dynamics simulations of beta2m oligomerization under oxidized and reduced conditions. We show that, consistent with experimental observations, the oxidized beta2m forms domain-swapped dimer, in which the two proteins exchange their N-terminal segments complementing each other. In contrast, both dimers and trimers, formed by reduced beta2m, are comprised of parallel beta-sheets between monomers and stabilized by the hydrogen bond network along the backbone. The oligomerized monomers are in extended conformations, capable of further aggregation. We find that both reduced and oxidized dimers are thermodynamically less stable than their corresponding monomers, indicating that beta2m oligomerization is not accompanied by the formation of a thermodynamically stable dimer. Our studies suggest that the different aggregation pathways of oxidized and reduced beta2m are dictated by the formation of distinct precursor oligomeric species that are modulated by Cys25-Cys80 disulfide-bonds. We propose that the propagation of domain swapping is the aggregation mechanism for the oxidized beta2m, while "parallel stacking" of partially unfolded beta2m is the aggregation mechanism for the reduced beta2m.