Probing the free energy landscapes of ALS disease mutants of SOD1 by NMR spectroscopy

MS785-MS27 Reactive Misfolded/Non-Native Zn-Deficient SOD1 Species Exhibit Cytotoxicity and Adopt Heterozygous Conformations in Motor Neurons

Misfolding of superoxide dismutase-1 (SOD1) is a pathological hallmark of amyotrophic lateral sclerosis (ALS) with SOD1 mutations. The development of antibodies specific for misfolded SOD1 deepens our understanding of how the protein participates in ALS pathogenesis. Since the term "misfolding" refers to various disordered conformers other than the natively folded one, which misfolded species are recognized by specific antibodies should be determined. Here, we molecularly characterized the recognition by MS785-MS27, an antibody cocktail experimentally confirmed to recognize over 100 ALS-linked SOD1 mutants. Indirect ELISA revealed that the antibody cocktail recognized Zn-deficient wild-type and mutated SOD1 species. It also recognized conformation-disordered wild-type and mutated SOD1 species, such as unfolded and oligomeric forms, but had less affinity for the aggregated form. Antibody-reactive SOD1 exhibited cytotoxicity to a motor neuron cell model, which was blocked by Zn treatment with Zn-deficient SOD1. Immunohistochemistry revealed antibody-reactive SOD1 mainly in spinal motor neurons of SOD1G93A mice throughout the disease course, and the distribution after symptomatic stages differed from that of other misfolded SOD1 species. This suggests that misfolded/non-native SOD1 species exist as heterogeneous populations. In conclusion, MS785-MS27 recognizes various conformation-disordered SOD1 species lacking the Zn ion.

Molecular dynamics analysis of superoxide dismutase 1 mutations suggests decoupling between mechanisms underlying ALS onset and progression

Mutations in the superoxide dismutase 1 (SOD1) gene are the second most common known cause of ALS. SOD1 variants express high phenotypic variability and over 200 have been reported in people with ALS. It was previously proposed that variants can be broadly classified in two groups, 'wild-type like' (WTL) and 'metal binding region' (MBR) variants, based on their structural location and biophysical properties. MBR variants, but not WTL variants, were associated with a reduction of SOD1 enzymatic activity. In this study we used molecular dynamics and large clinical datasets to characterise the differences in the structural and dynamic behaviour of WTL and MBR variants with respect to the wild-type SOD1, and how such differences influence the ALS clinical phenotype. Our study identified marked structural differences, some of which are observed in both variant groups, while others are group specific. Moreover, collecting clinical data of approximately 500 SOD1 ALS patients carrying variants, we showed that the survival time of patients carrying an MBR variant is generally longer (∼6 years median difference, p < 0.001) with respect to patients with a WTL variant. In conclusion, our study highlighted key differences in the dynamic behaviour between WTL and MBR SOD1 variants, and between variants and wild-type SOD1 at an atomic and molecular level, that could be further investigated to explain the associated phenotypic variability. Our results support the hypothesis of a decoupling between mechanisms of onset and progression of SOD1 ALS, and an involvement of loss-of-function of SOD1 with the disease progression.

Mutant SOD1 aggregates formed in vitro and in cultured cells are polymorphic and differ from those arising in the CNS

Mutations in the human Superoxide dismutase 1 (hSOD1) gene are well-established cause of the motor neuron disease ALS. Patients and transgenic (Tg) ALS model mice carrying mutant variants develop hSOD1 aggregates in the CNS. We have identified two hSOD1 aggregate strains, which both transmit spreading template-directed aggregation and premature fatal paralysis when inoculated into adult transgenic mice. This prion-like spread of aggregation could be a primary disease mechanism in SOD1-induced ALS. Human SOD1 aggregation has been studied extensively both in cultured cells and under various conditions in vitro. To determine how the structure of aggregates formed in these model systems related to disease-associated aggregates in the CNS, we used a binary epitope-mapping assay to examine aggregates of hSOD1 variants G93A, G85R, A4V, D90A, and G127X formed in vitro, in four different cell lines and in the CNS of Tg mice. We found considerable variability between replicate sets of in vitro-generated aggregates. In contrast, there was a high similarity between replicates of a given hSOD1 mutant in a given cell line, but pronounced variations between different hSOD1 mutants and different cell lines in both structures and amounts of aggregates formed. The aggregates formed in vitro or in cultured cells did not replicate the aggregate strains that arise in the CNS. Our findings suggest that the distinct aggregate morphologies in the CNS could result from a micro-environment with stringent quality control combined with second-order selection by spreading ability. Explorations of pathogenesis and development of therapeutics should be conducted in models that replicate aggregate structures forming in the CNS.

Quenched hydrogen-deuterium amide exchange optimization for high-resolution structural analysis of cellular protein aggregates

Inclusion bodies (IBs) are large, insoluble aggregates that often form during the overexpression of proteins in bacteria. These aggregates are of broad fundamental and practical significance, for recombinant protein preparation and due to their relevance to aggregation-related medical conditions and their recent emergence as promising functional nanomaterials. Despite their significance, high resolution knowledge of IB structure remains very limited. Such knowledge will advance understanding and control of IB formation and properties in myriad practical applications. Here, we report a detailed quenched hydrogen-deuterium amide exchange (qHDX) method with NMR readout to define the structure of IBs at the level of individual residues throughout the protein. Applying proper control of experimental conditions, such as sample pH, water content, temperature, and intrinsic rate of amide exchange, yields in depth results for these cellular protein aggregates. qHDX results illustrated for Cu, Zn superoxide dismutase 1 (SOD1) and Adnectins show their IBs include native-like structure and some but not all mutations alter IB structure.

High-Resolution NMR H/D Exchange of Human Superoxide Dismutase Inclusion Bodies Reveals Significant Native Features Despite Structural Heterogeneity

Protein aggregation is central to aging, disease and biotechnology. While there has been recent progress in defining structural features of cellular protein aggregates, many aspects remain unclear due to heterogeneity of aggregates presenting obstacles to characterization. Here we report high-resolution analysis of cellular inclusion bodies (IBs) of immature human superoxide dismutase (SOD1) mutants using NMR quenched amide hydrogen/deuterium exchange (qHDX), FTIR and Congo red binding. The extent of aggregation is correlated with mutant global stability and, notably, the free energy of native dimer dissociation, indicating contributions of native-like monomer associations to IB formation. This is further manifested by a common pattern of extensive protection against H/D exchange throughout nine mutant SOD1s despite their diverse characteristics. These results reveal multiple aggregation-prone regions in SOD1 and illuminate how aggregation may occur via an ensemble of pathways.

NMR methods for exploring 'dark' states in ligand binding and protein-protein interactions

A survey, primarily based on work in the authors' laboratory during the last 10 years, is provided of recent developments in NMR studies of exchange processes involving protein-ligand and protein-protein interactions. We start with a brief overview of the theoretical background of Dark state Exchange Saturation Transfer (DEST) and lifetime line-broadening (ΔR₂) NMR methodology. Some limitations of the DEST/ΔR₂ methodology in applications to molecular systems with intermediate molecular weights are discussed, along with the means of overcoming these limitations with the help of closely related exchange NMR techniques, such as the measurements of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion, exchange-induced chemical shifts or rapidly-relaxing components of relaxation decays. Some theoretical underpinnings of the quantitative description of global dynamics of proteins on the surface of very high molecular weight particles (nanoparticles) are discussed. Subsequently, several applications of DEST/ΔR₂ methodology are described from a methodological perspective with an emphasis on providing examples of how kinetic and relaxation parameters for exchanging systems can be reliably extracted from NMR data for each particular model of exchange. Among exchanging systems that are not associated with high molecular weight species, we describe several exchange NMR-based studies that focus on kinetic modelling of transient pre-nucleation oligomerization of huntingtin peptides that precedes aggregation and fibril formation.

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.

Measuring radiofrequency fields in NMR spectroscopy using offset-dependent nutation profiles

The application of NMR spectroscopy for studying molecular and reaction dynamics relies crucially on the measurement of the magnitude of radiofrequency (RF) fields that are used to nutate or lock the nuclear magnetization. Here, we report a method for measuring RF field amplitudes that leverages the intrinsic modulations observed in offset-dependent NMR nutation profiles of small molecules. Such nutation profiles are exquisitely sensitive to the magnitude of the RF field, and B₁ values ranging from 1 to 2000 Hz, as well the inhomogeneity in B₁ distributions, can be determined with high accuracy and precision using this approach. In order to measure B₁ fields associated with NMR experiments carried out on protein or nucleic acids, where these modulations are obscured by the large transverse relaxation rate constants of the analyte, our approach can be used in conjunction with a suitable external small molecule standard, expanding the scope of the method for large biomolecules.

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

With over 150 heritable mutations identified as disease-causative, superoxide dismutase 1 (SOD1) has been a main target of amyotrophic lateral sclerosis (ALS) research and therapeutic efforts. However, recent evidence has suggested that neither loss of function nor protein aggregation is responsible for promoting neurotoxicity. Furthermore, there is no clear pattern to the nature or the location of these mutations that could suggest a molecular mechanism behind SOD1-linked ALS. Here, we utilize reliable and accurate computational techniques to predict the perturbations of 10 such mutations to the free energy changes of SOD1 as it matures from apo monomer to metallated dimer. We find that the free energy perturbations caused by these mutations strongly depend on maturational progress, indicating the need for state-specific therapeutic targeting. We also find that many mutations exhibit similar patterns of perturbation to native and non-native maturation, indicating strong thermodynamic coupling between the dynamics at various sites of maturation within SOD1. These results suggest the presence of an allosteric network in SOD1 which is vulnerable to disruption by these mutations. Analysis of these perturbations may contribute to uncovering a unifying molecular mechanism which explains SOD1-linked ALS and help to guide future therapeutic efforts.

Superoxide Dismutase 1 in Health and Disease: How a Frontline Antioxidant Becomes Neurotoxic

Cu/Zn superoxide dismutase (SOD1) is a frontline antioxidant enzyme catalysing superoxide breakdown and is important for most forms of eukaryotic life. The evolution of aerobic respiration by mitochondria increased cellular production of superoxide, resulting in an increased reliance upon SOD1. Consistent with the importance of SOD1 for cellular health, many human diseases of the central nervous system involve perturbations in SOD1 biology. But far from providing a simple demonstration of how disease arises from SOD1 loss-of-function, attempts to elucidate pathways by which atypical SOD1 biology leads to neurodegeneration have revealed unexpectedly complex molecular characteristics delineating healthy, functional SOD1 protein from that which likely contributes to central nervous system disease. This review summarises current understanding of SOD1 biology from SOD1 genetics through to protein function and stability.

NMR spectroscopy captures the essential role of dynamics in regulating biomolecular function

Biomolecules are in constant motion. To understand how they function, and why malfunctions can cause disease, it is necessary to describe their three-dimensional structures in terms of dynamic conformational ensembles. Here, we demonstrate how nuclear magnetic resonance (NMR) spectroscopy provides an essential, dynamic view of structural biology that captures biomolecular motions at atomic resolution. We focus on examples that emphasize the diversity of biomolecules and biochemical applications that are amenable to NMR, such as elucidating functional dynamics in large molecular machines, characterizing transient conformations implicated in the onset of disease, and obtaining atomic-level descriptions of intrinsically disordered regions that make weak interactions involved in liquid-liquid phase separation. Finally, we discuss the pivotal role that NMR has played in driving forward our understanding of the biomolecular dynamics-function paradigm.

Exposure of β6/β7-Loop in Zn/Cu Superoxide Dismutase (SOD1) Is Coupled to Metal Loss and Is Transiently Reversible During Misfolding

Upon losing its structural integrity (misfolding), SOD1 acquires neurotoxic properties to become a pathogenic protein in ALS, a neurodegenerative disease targeting motor neurons; understanding the mechanism of misfolding may enable new treatment strategies for ALS. Here, we reported a monoclonal antibody, SE21, targeting the β6/β7-loop region of SOD1. The exposure of this region is coupled to metal loss and is entirely reversible during the early stages of misfolding. By using SE21 mAb, we demonstrated that, in apo-SOD1 incubated under the misfolding-promoting conditions, the reversible phase, during which SOD1 is capable of restoring its nativelike conformation in the presence of metals, is followed by an irreversible structural transition, autocatalytic in nature, which takes place prior to the onset of SOD1 aggregation and results in the formation of atypical apo-SOD1 that is unable to bind metals. The reversible phase defines a window of opportunity for pharmacological intervention using metal mimetics that stabilize SOD1 structure in its nativelike conformation to attenuate the spreading of the misfolding signal and disease progression by preventing the exposure of pathogenic SOD1 epitopes. Phenotypically similar apo-SOD1 species with impaired metal binding properties may also be produced via oxidation of Cys¹¹¹, underscoring the diversity of SOD1 misfolding pathways.

Visualizing and trapping transient oligomers in amyloid assembly pathways

Oligomers which form during amyloid fibril assembly are considered to be key contributors towards amyloid disease. However, understanding how such intermediates form, their structure, and mechanisms of toxicity presents significant challenges due to their transient and heterogeneous nature. Here, we discuss two different strategies for addressing these challenges: use of (1) methods capable of detecting lowly-populated species within complex mixtures, such as NMR, single particle methods (including fluorescence and force spectroscopy), and mass spectrometry; and (2) chemical and biological tools to bias the amyloid energy landscape towards specific oligomeric states. While the former methods are well suited to following the kinetics of amyloid assembly and obtaining low-resolution structural information, the latter are capable of producing oligomer samples for high-resolution structural studies and inferring structure-toxicity relationships. Together, these different approaches should enable a clearer picture to be gained of the nature and role of oligomeric intermediates in amyloid formation and disease.

•

Methods to study structure, toxicity, and kinetics of transient amyloid oligomers.

•

NMR and single particle methods can characterize lowly-populated oligomers.

•

Chemical tools/antibodies stabilize oligomers for structural and toxicity studies

•

A combination of methods is needed to fully characterize amyloid assembly pathways.

Amyotrophic Lateral Sclerosis: Proteins, Proteostasis, Prions, and Promises

Amyotrophic lateral sclerosis (ALS) is characterized by the progressive degeneration of the motor neurons that innervate muscle, resulting in gradual paralysis and culminating in the inability to breathe or swallow. This neuronal degeneration occurs in a spatiotemporal manner from a point of onset in the central nervous system (CNS), suggesting that there is a molecule that spreads from cell-to-cell. There is strong evidence that the onset and progression of ALS pathology is a consequence of protein misfolding and aggregation. In line with this, a hallmark pathology of ALS is protein deposition and inclusion formation within motor neurons and surrounding glia of the proteins TAR DNA-binding protein 43, superoxide dismutase-1, or fused in sarcoma. Collectively, the observed protein aggregation, in conjunction with the spatiotemporal spread of symptoms, strongly suggests a prion-like propagation of protein aggregation occurs in ALS. In this review, we discuss the role of protein aggregation in ALS concerning protein homeostasis (proteostasis) mechanisms and prion-like propagation. Furthermore, we examine the experimental models used to investigate these processes, including in vitro assays, cultured cells, invertebrate models, and murine models. Finally, we evaluate the therapeutics that may best prevent the onset or spread of pathology in ALS and discuss what lies on the horizon for treating this currently incurable disease.

Insight into the Autosomal-Dominant Inheritance Pattern of SOD1-Associated ALS from Native Mass Spectrometry

About 20% of all familial amyotrophic lateral sclerosis (ALS) cases are associated with mutations in superoxide dismutase (SOD1), a homodimeric protein. The disease has an autosomal-dominant inheritance pattern. It is, therefore, important to determine whether wild-type and mutant SOD1 subunits self-associate randomly or preferentially. A measure for the extent of bias in subunit association is the coupling constant determined in a double-mutant cycle type analysis. Here, cell lysates containing co-expressed wild-type and mutant SOD1 subunits were analyzed by native mass spectrometry to determine these coupling constants. Strikingly, we find a linear positive correlation between the coupling constant and the reported average duration of the disease. Our results indicate that inter-subunit communication and a preference for heterodimerization greatly increase the disease severity.

Friction-Limited Folding of Disulfide-Reduced Monomeric SOD1

The folding reaction of a stable monomeric variant of Cu/Zn superoxide dismutase (mSOD1), an enzyme responsible for the conversion of superoxide free radicals into hydrogen peroxide and oxygen, is known to be among the slowest folding processes that adhere to two-state behavior. The long lifetime, ∼10 s, of the unfolded state presents ample opportunities for the polypeptide chain to transiently sample nonnative structures before the formation of the productive folding transition state. We recently observed the formation of a nonnative structure in a peptide model of the C-terminus of SOD1, a sequence that might serve as a potential source of internal chain friction-limited folding. To test for friction-limited folding, we performed a comprehensive thermodynamic and kinetic analysis of the folding mechanism of mSOD1 in the presence of the viscogens glycerol and glucose. Using a, to our knowledge, novel analysis of the folding reactions, we found the disulfide-reduced form of the protein that exposes the C-terminal sequence, but not its disulfide-oxidized counterpart that protects it, experiences internal chain friction during folding. The sensitivity of the internal friction to the disulfide bond status suggests that one or both of the cross-linked regions play a critical role in driving the friction-limited folding. We speculate that the molecular mechanisms giving rise to the internal friction of disulfide-reduced mSOD1 might play a role in the amyotrophic lateral sclerosis-linked aggregation of SOD1.

The biophysics of superoxide dismutase-1 and amyotrophic lateral sclerosis

Few proteins have come under such intense scrutiny as superoxide dismutase-1 (SOD1). For almost a century, scientists have dissected its form, function and then later its malfunction in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). We now know SOD1 is a zinc and copper metalloenzyme that clears superoxide as part of our antioxidant defence and respiratory regulation systems. The possibility of reduced structural integrity was suggested by the first crystal structures of human SOD1 even before deleterious mutations in the sod1 gene were linked to the ALS. This concept evolved in the intervening years as an impressive array of biophysical studies examined the characteristics of mutant SOD1 in great detail. We now recognise how ALS-related mutations perturb the SOD1 maturation processes, reduce its ability to fold and reduce its thermal stability and half-life. Mutant SOD1 is therefore predisposed to monomerisation, non-canonical self-interactions, the formation of small misfolded oligomers and ultimately accumulation in the tell-tale insoluble inclusions found within the neurons of ALS patients. We have also seen that several post-translational modifications could push wild-type SOD1 down this toxic pathway. Recently we have come to view ALS as a prion-like disease where both the symptoms, and indeed SOD1 misfolding itself, are transmitted to neighbouring cells. This raises the possibility of intervention after the initial disease presentation. Several small-molecule and biologic-based strategies have been devised which directly target the SOD1 molecule to change the behaviour thought to be responsible for ALS. Here we provide a comprehensive review of the many biophysical advances that sculpted our view of SOD1 biology and the recent work that aims to apply this knowledge for therapeutic outcomes in ALS.

Nonnative structure in a peptide model of the unfolded state of superoxide dismutase 1 (SOD1): Implications for ALS-linked aggregation

Dozens of mutations throughout the sequence of the gene encoding superoxide dismutase 1 (SOD1) have been linked to toxic protein aggregation in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). A parsimonious explanation for numerous genotypes resulting in a common phenotype would be mutation-induced perturbation of the folding free-energy surface that increases the populations of high-energy states prone to aggregation. The absence of intermediates in the folding of monomeric SOD1 suggests that the unfolded ensemble is a potential source of aggregation. To test this hypothesis, here we dissected SOD1 into a set of peptides end-labeled with FRET probes to model the local behavior of the corresponding sequences in the unfolded ensemble. Using time-resolved FRET, we observed that the peptide corresponding to the Loop VII-β8 sequence at the SOD1 C terminus was uniquely sensitive to denaturant. Utilizing a two-dimensional form of maximum entropy modeling, we demonstrate that the sensitivity to denaturant is the surprising result of a two-state-like transition from a compact to an expanded state. Variations of the peptide sequence revealed that the compact state involves a nonnative interaction between the disordered N terminus and the hydrophobic C terminus of the peptide. This nonnative intramolecular structure could serve as a precursor for intermolecular association and result in aggregation associated with ALS. We propose that this precursor would provide a common molecular target for therapeutic intervention in the dozens of ALS-linked SOD1 mutations.

The molecular pathogenesis of superoxide dismutase 1-linked ALS is promoted by low oxygen tension

Mutations in superoxide dismutase 1 (SOD1) cause amyotrophic lateral sclerosis (ALS). Disease pathogenesis is linked to destabilization, disorder and aggregation of the SOD1 protein. However, the non-genetic factors that promote disorder and the subsequent aggregation of SOD1 have not been studied. Mainly located to the reducing cytosol, mature SOD1 contains an oxidized disulfide bond that is important for its stability. Since O₂ is required for formation of the bond, we reasoned that low O₂ tension might be a risk factor for the pathological changes associated with ALS development. By combining biochemical approaches in an extensive range of genetically distinct patient-derived cell lines, we show that the disulfide bond is an Achilles heel of the SOD1 protein. Culture of patient-derived fibroblasts, astrocytes, and induced pluripotent stem cell-derived mixed motor neuron and astrocyte cultures (MNACs) under low O₂ tensions caused reductive bond cleavage and increases in disordered SOD1. The effects were greatest in cells derived from patients carrying ALS-linked mutations in SOD1. However, significant increases also occurred in wild-type SOD1 in cultures derived from non-disease controls, and patients carrying mutations in other common ALS-linked genes. Compared to fibroblasts, MNACs showed far greater increases in SOD1 disorder and even aggregation of mutant SOD1s, in line with the vulnerability of the motor system to SOD1-mediated neurotoxicity. Our results show for the first time that O₂ tension is a principal determinant of SOD1 stability in human patient-derived cells. Furthermore, we provide a mechanism by which non-genetic risk factors for ALS, such as aging and other conditions causing reduced vascular perfusion, could promote disease initiation and progression.

An NMR View of Protein Dynamics in Health and Disease

Biological molecules are often highly dynamic, and this flexibility can be critical for function. The large range of sampled timescales and the fact that many of the conformers that are continually explored are only transiently formed and sparsely populated challenge current biophysical approaches. Solution nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for characterizing biomolecular dynamics in detail, even in cases where excursions involve short-lived states. Here, we briefly review a number of NMR experiments for studies of biomolecular dynamics on the microsecond-to-second timescale and focus on applications to protein and nucleic acid systems that clearly illustrate the functional relevance of motion in both health and disease.

Tryptophan 32 mediates SOD1 toxicity in a in vivo motor neuron model of ALS and is a promising target for small molecule therapeutics

SOD1 misfolding, toxic gain of function, and spread are proposed as a pathological basis of amyotrophic lateral sclerosis (ALS), but the nature of SOD1 toxicity has been difficult to elucidate. Uniquely in SOD1 proteins from humans and other primates, and rarely in other species, a tryptophan residue at position 32 (W32) is predicted to be solvent exposed and to participate in SOD1 misfolding. We hypothesized that W32 is influential in SOD1 acquiring toxicity, as it is known to be important in template-directed misfolding. We tested if W32 contributes to SOD1 cytotoxicity and if it is an appropriate drug target to ameliorate ALS-like neuromuscular deficits in a zebrafish model of motor neuron axon morphology and function (swimming). Embryos injected with human SOD1 variant with W32 substituted for a serine (SOD1^W32S) had reduced motor neuron axonopathy and motor deficits compared to those injected with wildtype or disease-associated SOD1. A library of FDA-approved small molecules was ranked with virtual screening based on predicted binding to W32, and subsequently filtered for analogues using a pharmacophore model based on molecular features of the uracil moiety of a small molecule previously predicted to interact with W32 (5'-fluorouridine or 5'-FUrd). Along with testing 5'-FUrd and uridine, a lead candidate from this list was selected based on its lower toxicity and improved blood brain barrier penetrance; telbivudine significantly rescued SOD1 toxicity in a dose-dependent manner. The mechanisms whereby the small molecules ameliorated motor neuron phenotypes were specifically mediated through human SOD1 and its residue W32, because these therapeutics had no measurable impact on the effects of UBQLN4^D90A, EtOH, or tryptophan-deficient human SOD1^W32S. By substituting W32 for a more evolutionarily conserved residue (serine), we confirmed the significant influence of W32 on human SOD1 toxicity to motor neuron morphology and function; further, we performed pharmaceutical targeting of the W32 residue for rescuing SOD1 toxicity. This unique residue offers future novel insights into SOD1 stability and toxic gain of function, and therefore poses an potential target for drug therapy.

A Hotspot for Disease-Associated Variants of Human PGM1 Is Associated with Impaired Ligand Binding and Loop Dynamics

Human phosphoglucomutase 1 (PGM1) plays a central role in cellular glucose homeostasis, catalyzing the conversion of glucose 1-phosphate and glucose 6-phosphate. Recently, missense variants of this enzyme were identified as causing an inborn error of metabolism, PGM1 deficiency, with features of a glycogen storage disease and a congenital disorder of glycosylation. Previous studies of selected PGM1 variants have revealed various mechanisms for enzyme dysfunction, including regions of structural disorder and side-chain rearrangements within the active site. Here, we examine variants within a substrate-binding loop in domain 4 (D4) of PGM1 that cause extreme impairment of activity. Biochemical, structural, and computational studies demonstrate multiple detrimental impacts resulting from these variants, including loss of conserved ligand-binding interactions and reduced mobility of the D4 loop, due to perturbation of its conformational ensemble. These potentially synergistic effects make this conserved ligand-binding loop a hotspot for disease-related variants in PGM1 and related enzymes.

General Expressions for Carr-Purcell-Meiboom-Gill Relaxation Dispersion for N-Site Chemical Exchange

The Carr-Purcell-Meiboom-Gill (CPMG) nuclear magnetic resonance experiment is widely used to characterize chemical exchange phenomena in biological macromolecules. Theoretical expressions for the nuclear spin relaxation rate constant for two-site chemical exchange during CPMG pulse trains valid for all time scales are well-known as are descriptions of N-site exchange in the fast limit. We have obtained theoretical expressions for N-site exchange outside of the fast limit by using approximations to an average Liouvillian describing the decay of magnetization during a CPMG pulse train. We obtain general expressions for CPMG experiments for any N-site scheme and all experimentally accessible time scales. For sufficiently slow chemical exchange, we obtain closed-form expressions for the relaxation rate constant and a general characteristic polynomial for arbitrary kinetic schemes. Furthermore, we highlight features that qualitatively characterize CPMG curves obtained for various N-site kinetic topologies, quantitatively characterize CPMG curves obtained from systems in various N-site exchange situations, and test distinguishability of kinetic models.

Large SOD1 aggregates, unlike trimeric SOD1, do not impact cell viability in a model of amyotrophic lateral sclerosis

Aberrant accumulation of misfolded Cu, Zn superoxide dismutase (SOD1) is a hallmark of SOD1-associated amyotrophic lateral sclerosis (ALS), an invariably fatal neurodegenerative disease. While recent discovery of nonnative trimeric SOD1-associated neurotoxicity has suggested a potential pathway for motor neuron impairment, it is yet unknown whether large, insoluble aggregates are cytotoxic. Here we designed SOD1 mutations that specifically stabilize either the fibrillar form or the trimeric state of SOD1. The designed mutants display elevated populations of fibrils or trimers correspondingly, as demonstrated by gel filtration chromatography and electron microscopy. The trimer-stabilizing mutant, G147P, promoted cell death, even more potently in comparison with the aggressive ALS-associated mutants A4V and G93A. In contrast, the fibril-stabilizing mutants, N53I and D101I, positively impacted the survival of motor neuron-like cells. Hence, we conclude the SOD1 oligomer and not the mature form of aggregated fibril is critical for the neurotoxic effects in the model of ALS. The formation of large aggregates is in competition with trimer formation, suggesting that aggregation may be a protective mechanism against formation of toxic oligomeric intermediates.

Phylogenetic spread of sequence data affects fitness of SOD1 consensus enzymes: Insights from sequence statistics and structural analyses

Non-natural protein sequences with native-like structures and functions can be constructed successfully using consensus design. This design strategy is relatively well understood in repeat proteins with simple binding function, however detailed studies are lacking in globular enzymes. The SOD1 family is a good model for such studies due to the availability of large amount of sequence and structure data motivated by involvement of human SOD1 in the fatal motor neuron disease amyotrophic lateral sclerosis (ALS). We constructed two consensus SOD1 enzymes from multiple sequence alignments from all organisms and eukaryotic organisms. A significant difference in their catalytic activities shows that the phylogenetic spread of the sequences used affects the fitness of the construct obtained. A mutation in an electrostatic loop and overall design incompatibilities between bacterial and eukaryotic sequences were implicated in this disparity. Based on this analysis, a bioinformatics approach was used to classify mutations thought to cause familial ALS providing a unique high level view of the physical basis of disease-causing aggregation of human SOD1.

Effects of maturation on the conformational free-energy landscape of SOD1

Amyotrophic lateral sclerosis (ALS) is a devastating fatal syndrome characterized by very rapid degeneration of motor neurons. A leading hypothesis is that ALS is caused by toxic protein misfolding and aggregation, as also occurs in many other neurodegenerative disorders, such as prion, Alzheimer's, Parkinson's, and Huntington's diseases. A prominent cause of familial ALS is mutations in the protein superoxide dismutase (SOD1), which promote the formation of misfolded SOD1 conformers that are prone to aberrant interactions both with each other and with other cellular components. We have shown previously that immature SOD1, lacking bound Cu and Zn metal ions and the intrasubunit disulfide bond (apoSOD1^2SH), has a rugged free-energy surface (FES) and exchanges with four other conformations (excited states) that have millisecond lifetimes and sparse populations on the order of a few percent. Here, we examine further states of SOD1 along its maturation pathway, as well as those off-pathway resulting from metal loss that have been observed in proteinaceous inclusions. Metallation and disulfide bond formation lead to structural transformations including local ordering of the electrostatic loop and native dimerization that are observed in rare conformers of apoSOD1^2SH; thus, SOD1 maturation may occur via a population-switch mechanism whereby posttranslational modifications select for preexisting structures on the FES. Metallation and oxidation of SOD1 stabilize the native, mature conformation and decrease the number of detected excited conformational states, suggesting that it is the immature forms of the protein that contribute to misfolded conformations in vivo rather than the highly stable enzymatically active dimer.

Disulfide driven folding for a conditionally disordered protein

Conditionally disordered proteins are either ordered or disordered depending on the environmental context. The substrates of the mitochondrial intermembrane space (IMS) oxidoreductase Mia40 are synthesized on cytosolic ribosomes and diffuse as intrinsically disordered proteins to the IMS, where they fold into their functional conformations; behaving thus as conditionally disordered proteins. It is not clear how the sequences of these polypeptides encode at the same time for their ability to adopt a folded structure and to remain unfolded. Here we characterize the disorder-to-order transition of a Mia40 substrate, the human small copper chaperone Cox17. Using an integrated real-time approach, including chromatography, fluorescence, CD, FTIR, SAXS, NMR, and MS analysis, we demonstrate that in this mitochondrial protein, the conformational switch between disordered and folded states is controlled by the formation of a single disulfide bond, both in the presence and in the absence of Mia40. We provide molecular details on how the folding of a conditionally disordered protein is tightly regulated in time and space, in such a way that the same sequence is competent for protein translocation and activity.

Disulfide driven folding for a conditionally disordered protein

Conditionally disordered proteins are either ordered or disordered depending on the environmental context. The substrates of the mitochondrial intermembrane space (IMS) oxidoreductase Mia40 are synthesized on cytosolic ribosomes and diffuse as intrinsically disordered proteins to the IMS, where they fold into their functional conformations; behaving thus as conditionally disordered proteins. It is not clear how the sequences of these polypeptides encode at the same time for their ability to adopt a folded structure and to remain unfolded. Here we characterize the disorder-to-order transition of a Mia40 substrate, the human small copper chaperone Cox17. Using an integrated real-time approach, including chromatography, fluorescence, CD, FTIR, SAXS, NMR, and MS analysis, we demonstrate that in this mitochondrial protein, the conformational switch between disordered and folded states is controlled by the formation of a single disulfide bond, both in the presence and in the absence of Mia40. We provide molecular details on how the folding of a conditionally disordered protein is tightly regulated in time and space, in such a way that the same sequence is competent for protein translocation and activity.

A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations

Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.

Utilizing tagged paramagnetic shift reagents to monitor protein dynamics by NMR

Calmodulin is a ubiquitous calcium sensor protein, known to serve as a critical interaction hub with a wide range of signaling partners. While the holo form of calmodulin (CaM-4Ca²⁺) has a well-defined ground state structure, it has been shown to undergo exchange, on a millisecond timescale, to a conformation resembling that of the peptide bound state. Tagged paramagnetic relaxation agents have been previously used to identify long-range dipolar interactions through relaxation effects on nuclear spins of interest. In the case of calmodulin, this lead to the determination of the relative orientation of the N- and C-terminal domains and the presence of a weakly populated peptide bound like state. Here, we make use of pseudocontact shifts from a tagged paramagnetic shift reagent which allows us to define minor states both in ¹³C and ¹⁵N NMR spectra and through ¹³C- and ¹⁵N-edited ¹H-CPMG relaxation dispersion measurements. This is validated by pulsed EPR (DEER) spectroscopy which reveals an ensemble consisting of a compact peptide-bound like conformer, an intermediate peptide-bound like conformer, and a (dumbbell-like) extended ground state conformer of CaM-4Ca²⁺, where addition of the MLCK peptide increases the population of the peptide-bound conformers. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Parsing disease-relevant protein modifications from epiphenomena: perspective on the structural basis of SOD1-mediated ALS

Conformational change and modification of proteins are involved in many cellular functions. However, they can also have adverse effects that are implicated in numerous diseases. How structural change promotes disease is generally not well-understood. This perspective illustrates how mass spectrometry (MS), followed by toxicological and epidemiological validation, can discover disease-relevant structural changes and therapeutic strategies. We (with our collaborators) set out to characterize the structural and toxic consequences of disease-associated mutations and post-translational modifications (PTMs) of the cytosolic antioxidant protein Cu/Zn-superoxide dismutase (SOD1). Previous genetic studies discovered >180 different mutations in the SOD1 gene that caused familial (inherited) amyotrophic lateral sclerosis (fALS). Using hydrogen-deuterium exchange with mass spectrometry, we determined that diverse disease-associated SOD1 mutations cause a common structural defect - perturbation of the SOD1 electrostatic loop. X-ray crystallographic studies had demonstrated that this leads to protein aggregation through a specific interaction between the electrostatic loop and an exposed beta-barrel edge strand. Using epidemiology methods, we then determined that decreased SOD1 stability and increased protein aggregation are powerful risk factors for fALS progression, with a combined hazard ratio > 300 (for comparison, a lifetime of smoking is associated with a hazard ratio of ~15 for lung cancer). The resulting structural model of fALS etiology supported the hypothesis that some sporadic ALS (sALS, ~80% of ALS is not associated with a gene defect) could be caused by post-translational protein modification of wild-type SOD1. We developed immunocapture antibodies and high sensitivity top-down MS methods and characterized PTMs of wild-type SOD1 using human tissue samples. Using global hydrogen-deuterium exchange, X-ray crystallography and neurotoxicology, we then characterized toxic and protective subsets of SOD1 PTMs. To cap this perspective, we present proof-of-concept that post-translational modification can cause disease. We show that numerous mutations (N➔D; Q➔E), which result in the same chemical structure as the PTM deamidation, cause multiple diseases. Copyright © 2017 John Wiley & Sons, Ltd.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

Free energy landscape remodeling of the cardiac pacemaker channel explains the molecular basis of familial sinus bradycardia

The hyperpolarization-activated and cyclic nucleotide-modulated ion channel (HCN) drives the pacemaker activity in the heart, and its malfunction can result in heart disorders. One such disorder, familial sinus bradycardia, is caused by the S672R mutation in HCN, whose electrophysiological phenotypes include a negative shift in the channel activation voltage and an accelerated HCN deactivation. The outcomes of these changes are abnormally low resting heart rates. However, the molecular mechanism underlying these electrophysiological changes is currently not fully understood. Crystallographic investigations indicate that the S672R mutation causes limited changes in the structure of the HCN intracellular gating tetramer, but its effects on protein dynamics are unknown. Here, we utilize comparative S672R versus WT NMR analyses to show that the S672R mutation results in extensive perturbations of the dynamics in both apo- and holo-forms of the HCN4 isoform, reflecting how S672R remodels the free energy landscape for the modulation of HCN4 by cAMP, i.e. the primary cyclic nucleotide modulator of HCN channels. We show that the S672R mutation results in a constitutive shift of the dynamic auto-inhibitory equilibrium toward inactive states of HCN4 and broadens the free-energy well of the apo-form, enhancing the millisecond to microsecond dynamics of the holo-form at sites critical for gating cAMP binding. These S672R-induced variations in dynamics provide a molecular basis for the electrophysiological phenotypes of this mutation and demonstrate that the pathogenic effects of the S672R mutation can be rationalized primarily in terms of modulations of protein dynamics.