Original Research ARTICLE
Susceptibility of Mutant SOD1 to Form a Destabilized Monomer Predicts Cellular Aggregation and Toxicity but Not In vitro Aggregation Propensity
- 1Lab 210, Illawarra Health and Medical Research Institute, Wollongong, NSW, Australia
- 2Science Medicine and Health Faculty, School of Biological Sciences, University of Wollongong, Wollongong, NSW, Australia
Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by the rapid and progressive degeneration of upper and lower motor neurons in the spinal cord, brain stem and motor cortex. The first gene linked to ALS was the gene encoding the free radical scavenging enzyme superoxide dismutase-1 (SOD1) that currently has over 180, mostly missense, ALS-associated mutations identified. SOD1-associated fALS patients show remarkably broad mean survival times (<1 year to ~17 years death post-diagnosis) that are mutation dependent. A hallmark of SOD1-associated ALS is the deposition of SOD1 into large insoluble aggregates in motor neurons. This is thought to be a consequence of mutation induced structural destabilization and/or oxidative damage leading to the misfolding and aggregation of SOD1 into a neurotoxic species. Here we aim to understand the relationship between SOD1 variant toxicity, structural stability, and aggregation propensity using a combination of cell culture and purified protein assays. Cell based assays indicated that aggregation of SOD1 variants correlate closely to cellular toxicity. However, the relationship between cellular toxicity and disease severity was less clear. We next utilized mass spectrometry to interrogate the structural consequences of metal loss and disulfide reduction on fALS-associated SOD1 variant structure. All variants showed evidence of unfolded, intermediate, and compact conformations, with SOD1G37R, SOD1G93A and SOD1V148G having the greatest abundance of intermediate and unfolded SOD1. SOD1G37R was an informative outlier as it had a high propensity to unfold and form oligomeric aggregates, but it did not aggregate to the same extent as SOD1G93A and SOD1V148G in in vitro aggregation assays. Furthermore, seeding the aggregation of DTT/EDTA-treated SOD1G37R with preformed SOD1G93A fibrils elicited minimal aggregation response, suggesting that the arginine substitution at position-37 blocks the templating of SOD1 onto preformed fibrils. We propose that this difference may be explained by multiple strains of SOD1 aggregate and this may also help explain the slow disease progression observed in patients with SOD1G37R.
Amyotrophic lateral sclerosis (ALS) is a disease characterized by the progressive degeneration of both upper and lower motor neurons, ultimately resulting in death. The symptoms of ALS typically manifest in people aged 55–65 years, with younger patients rarely being identified (Gordon, 2011). Post-diagnosis, the average disease duration is approximately 3 years, although patients suffering different familial variants can have shorter and longer durations. Genetic studies of ALS patients have revealed that 90% of cases are sporadic (sALS) and the remaining 10% of cases are familial (fALS) (Renton et al., 2014). Several of the fALS associated genes (TDP-43, C9orf72, FUS, ALS2) have roles in the biogenesis and trafficking of RNA species, and another clear grouping are associated with protein degradation (CCNF, VCP, SQSTM1, UBQLN2, OPTN, TBK1); most are also linked to frontotemporal dementia (Cirulli et al., 2015; Peters et al., 2015; Williams et al., 2016). The first fALS linked gene discovered, and most studied, is the gene encoding superoxide dismutase-1 (SOD1) which is only associated with a pure motor neuron disease phenotype.
Superoxide dismutase-1 (SOD1) is a cytosolic 32 kDa homodimeric enzyme that catalyzes the dismutation of oxygen radicals to either molecular oxygen or hydrogen peroxide (McCord and Fridovich, 1969). Each subunit contains a catalytically active copper ion, a zinc ion, and an intramolecular disulfide bond. These post-translational modifications confer substantial stability to the native fold of SOD1, endowing it with resistance to heat denaturation and maintenance of catalytic activity in denaturing conditions (Senoo et al., 1988; Rodriguez et al., 2002). Zinc binding occurs spontaneously following synthesis and folding, whereas copper binding and disulfide bond formation are mediated by association with the copper-chaperone for SOD1 (Rae et al., 2000; Furukawa et al., 2004). Currently, there are known to be over 180 natural mutations that can occur within the SOD1 sequence that are associated with fALS [http://alsod.iop.kcl.ac.uk - (Abel et al., 2013)]. These mutations occur throughout the SOD1 sequence, and it is thought that all have the ability to destabilize the native fold, leading to an increased propensity to misfold and/or aggregate (Khare et al., 2006; Nordlund and Oliveberg, 2008; Byrstrom et al., 2010; Munch and Bertolotti, 2010). Owing to the growing evidence for prion-like propagation associated with SOD1 fALS (Ayers et al., 2014, 2016; Grad et al., 2014; Zeineddine et al., 2015), and the fact that a misfolded/aggregated wild-type SOD1 species has been found associated with sALS (Bosco et al., 2010), it has been suggested that there is a common pathogenic conformer or aggregate species that exists for all SOD1 variants, although it is yet to be elucidated.
Metal coordination and disulfide bond formation are important for stability of SOD1 tertiary and quaternary structures (Hayward et al., 2002; Tiwari and Hayward, 2003; Lynch et al., 2004; Svensson et al., 2010), and their absence has been shown to increase the aggregation propensity of SOD1 (Arnesano et al., 2004; Ding and Dokholyan, 2008). From this, it has been inferred that an immature/unfolded conformation(s) of SOD1 is responsible for aggregation. Investigation of these possible unfolded SOD1 conformations has been difficult due to the heterogeneity of non-native protein structures. High resolution analyses, such as X-ray crystallography, are impeded by the existence of disordered protein regions that can assume multiple conformations (Banci et al., 2009), and typical separation techniques, such as size exclusion chromatography (SEC) do not provide sufficient resolution to separate different conformations for downstream analysis (Hong et al., 2012), although use of small-angle X-ray scattering coupled with SEC can provide valuable insight into the structural plasticity of proteins in solution by separating different conformers and oligomeric states (Wright et al., 2011, 2016; Furukawa et al., 2015). Nuclear magnetic resonance (NMR) spectroscopy has been used to great extent on disordered SOD1, monitoring unfolding, refolding, and even identifying a zinc-deficient form of SOD1 proposed as a precursor to the toxic species (Luchinat et al., 2014; Szpryngiel et al., 2015). Although NMR analysis is an extremely useful tool in structural biology, it requires high concentrations of protein (0.1−3 mM), which can lead to undesirable artifacts, especially in the case of protein-protein interactions.
In comparison, native mass spectrometry (MS) is a powerful tool for the analysis of structurally heterogeneous protein samples as it allows for the separation of different protein species based upon mass, charge, and drift time when coupled to ion mobility (IM-MS) (Mehmood et al., 2015), allowing for the resolution of species in solution at variable concentrations. This study details the use of native MS to investigate the possible misfolded/unfolded conformations SOD1 can occupy when metal binding and disulfide bonding are disrupted, and the relationship of these species to SOD1 aggregation. We first determined if there was a correlation between isolated recombinant protein aggregation and cellular aggregation of different SOD1 variants, focusing on SOD1-EGFP fusion proteins. Next we attempted to assess the oligomeric state of our SOD1 variants under similar aggregation conditions, finding that SOD1 variants had differential mobility in both SEC and native PAGE. This was followed by native MS analysis to separate and analyse the variably destabilized SOD1 species. Our observations indicate that the oligomeric state of the SOD1 variants was primarily monomeric and that the monomer assumed 3 conformations based upon mass-to-charge distributions. The relative abundances of the conformations were dependent upon the SOD1 variant being analyzed and the formation of destabilized SOD1 monomer correlated with cellular aggregation. Lastly, the slow progressing SOD1G37R mutant had a high proportion of unfolded monomer, high propensity to form oligomeric aggregates and insoluble aggregates, but could not be seeded by G93A aggregates in a plate-based aggregation assay, suggesting that it may exist in a different strain to other aggregation prone mutants.
Expression and Purification of Recombinant SOD1
Plasmids encoding SOD1WT and SOD1G93A were a gift from the Oliveberg group (Stockholm, Sweden). Plasmids encoding SOD1G37R, SOD1H46R, SOD1D90A, and SOD1V148G were designed in-house and generated by Genscript (New Jersey, USA). Protein expression and purification were performed according to (Lindberg et al., 2002) with slight modifications. Briefly, plasmids containing genes for the expression of SOD1 and yeast-CCS were transformed into chemically competent BL21 (DE3) E. coli using heat-shock. Following transformation expression cultures were induced using IPTG in the presence of copper and zinc. Following lysis and ammonium-sulfate precipitation, the expressed protein was purified using size exclusion chromatography (Hiload 16/60 Superdex 75 PG, GE USA) and anion exchange chromatography (HiTrap DEAE, GE USA). Purity was assessed by SDS-PAGE and mass spectrometry, with pure samples snap frozen and stored in 1 × PBS at −20°C. Misfolded/unfolded SOD1 was generated by incubating purified SOD1 variants at a concentration of 30 μM in 1 × PBS with 5 mM EDTA and 20 mM DTT at 37°C for 2 h. All protein concentrations were determined using a bicinchoninic acid assay.
Isolated Recombinant Protein Aggregation Assays
Superoxide dismutase-1 (SOD1) variants were aggregated at a concentration of 30 μM (dimer) in 1 × PBS (pH 7.4) containing 20 mM DTT, 5 mM EDTA, with 10 μM thioflavin T (ThT). Plate-reader assays were performed on a POLARstar Omega (BMG labtech) in clear bottomed 384-well plates (Greiner), with a final well volume of 50 μl. Following addition of DTT/EDTA to wells containing SOD1 protein, the plate was incubated at 37°C for 30 min before being covered with an adhesive slip. Aggregation was induced with double-orbital shaking at 300 rpm for 330 s at the start of a 900 s cycle. ThT was excited at 440 nm and its fluorescence was measured at 490 nm. Seeded aggregation assays were carried out similarly, with the exception that the reaction mixture was 30 μM SOD1, 10 mM DTT, 1 mM EDTA, 10 μM ThT, in 1 × PBS (pH 7.4). Seeded assays also contained 0.3% (w/w) of total protein as seed from a previous aggregation assay. Analysis of aggregation kinetics was carried out as described by Cox et al. (2016).
Native Page Analysis
Superoxide dismutase-1 (SOD1) variants at a concentration of 30 μM (10 μg total protein per well) were loaded into stain-free any kD gradient gels (Biorad, USA) in a tris-glycine buffer. Gels were electrophoresed at 100 V for 3 h at 4°C with stirring. Following electrophoresis, gels were imaged on a stain-free imager (Biorad, USA).
Analytical Gel-Filtration and Buffer Exchange
Analytical gel-filtration chromatography was carried out in conjunction with buffer exchange of SOD1 variants into 200 mM NH4OAc (pH 6.8). SOD1 variants were concentrated using microfuge concentrators (Vivaspin 25 10 kDa MWCO, GE USA) to a concentration of ~450 μM prior to analytical gel filtration. Concentrated native and DTT/EDTA-treated SOD1 variants were loaded onto either a Superdex-75 10/300 column or Superdex-200 10/300 column (GE, USA) at a flow rate of 0.5 ml/min with the elution profile being measured at 280 nm, and 0.5 ml fractions being collected for immediate MS analysis. Fractions were placed on ice as they eluted.
Mass spectrometry analysis was performed using a SYNAPT G1 HDMS (Waters, UK) with parameters set according to previous work (McAlary et al., 2013). Briefly, SOD1 samples at 10 μM in 200 mM NH4 OAc were loaded into gold-coated borosilicate capillaries (made in-house) and subjected to nano-electrospray ionization. All spectra were externally calibrated using 10 mg/ml caesium-iodide in 50% n-propanol, and were processed using Masslynx 4.1. For determination of the abundances of the observed conformations MATLAB R2014b (Version 8.4) was used to fit Gaussians to the data.
Secondary Tissue Culture
Mouse neuroblastoma/motor neuron hybrid cells (NSC-34) (Cashman et al., 1992) were cultured in Dulbecco's modified Eagles medium-F12 (DMEM-F12) (Invitrogen, Australia), supplemented with 10% (v/v) heat inactivated fetal bovine serum (FBS) (Bovogen Biologicals, Australia). In order to passage and plate cells, they were washed once with DMEM-F12 and treated with 0.25% trypsin, 0.02% EDTA (Sigma-Aldrich, Australia) to lift off the adherent cells. The cells were pelleted via centrifugation (500 × g, 5 min) and resuspended in pre-warmed DMEM-F12, supplemented with 10% FBS. Following washing, plates and chamber-slides were seeded at a confluency of 40% and cultured at 37°C in a humidified incubator with 5% atmospheric CO2 for 24 h prior to transfection.
Transfection of NSC-34 Cells
NSC-34 cells were transfected, 24 h post-plating into either multi-well plates or chamber-slides (Ibidi Germany), using Lipofectamine 3000 (Invitrogen, USA) according to the manufacturer's instructions. Lipofectamine 3000 reagent was added to pre-warmed DMEM-F12 to a concentration of 6% (v/v). P3000 reagent and plasmid DNA were added to a separate DMEM-F12 volume to concentrations of 4% (v/v) and 20 ng/μL, respectively. The DNA/lipofectamine complex was generated by mixing one volume of 6% (v/v) lipofectamine in DMEM-F12, with one volume of DNA/P3000 mix in DMEM-F12. The DNA/lipofectamine mixture was incubated for 30 min at room temperature, after which it was dispensed onto cultured cells. For cells that were co-transfected, the total amount of DNA was halved between plasmid constructs.
All cell imaging was performed using a Leica TCS SP5 Confocal microscope (Leica, Germany). Cells were imaged live in 5% atmospheric CO2, with the incubation unit set to a temperature of 37°C. For all confocal experiments, the argon laser power was at 20%. All green-fluorescent protein (EGFP) tagged constructs were excited with 10% laser power at 488 nm and emission was analyzed from 507 to 550 nm. Images were acquired using a HCX PL APO 63.0 × 1.20 water immersion objective, using PMT detectors, with a pin-hole of 100 μm.
To quantify inclusion formation of SOD1-EGFP constructs, z-stacks of transfected cells were acquired. Cells were cultured into μSlide 8-well chamber slides (Ibidi, Germany) and transfected with c-terminal SOD1-EGFP constructs. Imaging was performed 48 h post-transfection. Z-stacks series were acquired at 512 × 512 pixels, laser scanning speed 400 Hz, line and frame averaging = 1, where each image was a 400-500 nm slice of the z-axis. Each Z-stack was systematically acquired in a new area of the well to avoid overlap. At least five z-stack series were acquired per well of a chamber slide, accounting for at least 100 transfected cells imaged per well. Z-stacks were processed into a single image using LAS-AF Lite software (Leica, Germany), where cells containing inclusions were determined by the presence of bright fluorescent puncta in the cytoplasmic region greater than 2 μm in diameter.
Time Resolved Fluorescent Imaging
To assess the toxicity of SOD1-EGFP constructs in NSC-34 cells over a time course, an IncuCyte® automated fluorescent microscope (Essen BioScience, USA) was used. NSC-34 cells were plated into 12-wells plates at a confluency of 60% and were transfected 24 h post-plating. Cells were dissociated 24 h post-transfection and plated into 96-well plates at a confluency of 20% in phenol-red-free DMEM-F12 supplemented with 10% FBS. At least 3 images were acquired per-well at 2 h time points for 68 h in both phase and green channels with the green channel acquisition time at 400 ms. The processing definition generated to analyse the images utilized top-hat back ground subtraction [radius = 100 μm, threshold = 0.5 general calibration units (GCU)], edge-splitting (edge sensitivity = 0), filters (minimum area = 175 μm2, minimum mean intensity = 1.2). The equations below describe data analysis and presentation.
First, the number of GFP positive cells at each time point (GFPtx) for each transfection was normalized to the intial value (GFPt0) determined in the first scan after plating.
Then, the normalized values of the SOD1 mutants at each time point (Normalised GFPtx) were divided by the normalized SOD1WT data at the same time points to determine the proportion of GFP positive cells relative to SOD1WT.
FloIT Assay by Flow Cytometry
The FloIT assay was performed as described in (Whiten et al., 2016). Briefly, NSC-34 cells were transiently transfected in a 24-well plate with SOD1-EGFP constructs as described above. Cell were harvested, washed twice and resuspended in 1 × PBS (0.5 mL/tube). A 150 μL aliquot was taken and used to analyse transfection efficiency. The remaining cells were pelleted (300 × g, 5 min) and placed on ice. Supernatant was removed from each sample and the cells were lysed with the addition of 0.5% Triton X-100 in 1 × PBS containing a 1:1000 dilution of RedDot2 immediately prior to flow cytometric analysis. Events were analyzed using a LSRFortessa X-20 Cell Analyzer (BD Biosciences) with the excitation set to 488 nm and the emission collected with a 525/50 band-pass filter for EGFP, and an excitation of 640 nm with the emission collected using a 670/30 nm band-pass filter for RedDot2. All parameters were set to log10 during acquisition from cell lysates, which were gated according to forward and side scatter. The forward scatter threshold was set to the minimum value (200 AU) to minimize the exclusion of small protein inclusions. Nuclei were identified based on RedDot2 fluorescence and forward scatter for removal from further analysis. The remaining particles were analyzed for presence of inclusions based on EGFP fluorescence, forward scatter, and comparison to lysates of NSC-34 cells expressing EGFP only. The number of inclusions in the population of particles was normalized to the number of nuclei, and reported as inclusions/100 transfected cells (iFloIT) according to the equation
where ni represents the number of inclusions acquired, nnuc is the number of nuclei acquired, and γ is the transfection efficiency (expressed as a fraction).
Cellular Aggregation Correlates with Cell Loss
It has been previously shown that SOD1 variants expressed in cultured cells have different propensities to aggregate (Turner et al., 2005; Prudencio et al., 2009). This variation in aggregation may offer an explanation to the large disparity in average length of disease duration between SOD1 mutations. We first sought to determine the aggregation propensities of our chosen SOD1 variants as EGFP fusion proteins in cultured cells using microscopy. To this end, mouse neuroblastoma/motor neuron hybrid cells (NSC-34) (Cashman et al., 1992) were transfected with SOD1-EGFP constructs and imaged 48 h post-transfection and the number of transfected cells that formed inclusions counted. We chose to work with variants SOD1WT, SOD1G37R, SOD1H46R, SOD1D90A, SOD1G93A, and SOD1V148G owing to their variable patient survival time (SOD1G37R = ~17, SOD1H46R = ~17.6, SOD1D90A ~8, SOD1G93A ~3.1, SOD1V148G = ~2.1 years) (Wang et al., 2008). We observed that variants, SOD1G37R, SOD1G93A, and SOD1V148G had significantly increased populations of transfected cells containing inclusions when compared to SOD1WT (Figures 1A,B). In contrast, there was no significant difference between SOD1WT and the variants SOD1H46R and SOD1D90A (Figure 1B). The flow cytometry based FloIT assay showed results similar to manual counting (Figure 1C). While the FloIT technique finds 2–3 times more inclusions per 100 cells, this is due to the fact that multiple inclusions can occur in a single cell (Whiten et al., 2016). Regardless, the outcome from manual counting and FloIT analysis demonstrate that SOD1G37R, SOD1G93A, and SOD1V148G aggregate significantly more in cells that the other variants tested. We next examined the effects of overexpression of SOD1-EGFP constructs on cell survival through time using real-time imaging. The number of surviving cells was followed and plotted relative to SOD1WT survival (Figure 1D). As previously observed (Turner et al., 2005), the survival of cells expressing SOD1 variants correlated with the number of cells with aggregates (Figure 1E). That is the variants SOD1G37R, SOD1G93A, and SOD1V148G had significant levels of cell death, while SOD1H46R had an intermediate level of cell death and SOD1D90A was no different to SOD1WT. This suggests that the cellular aggregation propensity of the variant predicts cell death in NSC-34 cells.
Figure 1. SOD1 variants have differential aggregation propensities in NSC-34 cells. NSC-34 cells were transiently transfected with SOD1-EGFP constructs and imaged 48 h following transfection to examine the formation of inclusions. (A) Representative images of NSC-34 cells transfected with different SOD1-EGFP variants with arrows marking cells with inclusions. Scale bar represents 50 μM. (B) Manual counting of inclusions from imaged cells showed that SOD1variants had differential aggregation propensities when expressed in NSC-34 cells. (C) FloIT analysis of transfected NSC-34 cells. Cells were lysed and the inclusions counted using flow cytometry 48 h post-transfection. (D) The toxicity of SOD1-EGFP constructs compared to SOD1WT in NSC-34 cells. (E) Cell death is plotted against cellular aggregation propensity, there is a plausible correlation between aggregate load and cell death at 48 h post-transfection. Statistical significance was determined via a one-way ANOVA with a Tukey's multiple comparison post-test (***P <0.001, **P <0.01, ns = non-significant).
Aggregation Propensity of Recombinant SOD1 Does Not Always Predict Its Cellular Aggregation
The variation in cellular aggregation propensity between SOD1 variants could be due to their intrinsic propensity to aggregate, or given that cells actively generate inclusions in a manner thought to represent a quality control compartment (Weisberg et al., 2012; Farrawell et al., 2015), the level of aggregation may represent an ability of a particular variant to evade cellular protein quality control machinery (Yerbury et al., 2013, 2016). To distinguish between these two possibilities, we sought to assess if there was a correlation between aggregation of purified recombinant SOD1 and SOD1-GFP aggregation in cells. Aggregation of recombinant SOD1 was performed as previously described by us and others (Furukawa et al., 2008; Roberts et al., 2013; Yerbury et al., 2013), making use of DTT (to reduce the intramolecular disulfide) and EDTA (to chelate Cu and Zn), as well as shaking, to promote aggregation of the SOD1 variants. Aggregation was measured via Thioflavin T fluorescence, which has been used as a measure of amyloid formation and is thought to fluoresce when it binds to the β-sheet structure characteristic of amyloid. The aggregation propensities differed between SOD1 variants (Figures 2A,B), with SOD1V148G, SOD1G93A, SOD1G37R, and SOD1D90A obeying Boltzmann-sigmoidal kinetics typical of fibrillar aggregation. SOD1 aggregation is noted to be a highly stochastic process, with large variance between technical and biological replicates (Abdolvahabi et al., 2016); however, a comparison of the rates of fibril elongation (measured as change in ThT fluorescence; Figure 2C) between two assays here showed that inter-assay variation was minimal with the exception of SOD1G93A. SOD1WT, and SOD1H46R did not fit standard fibrillar aggregation kinetics, and as such analysis of the fibril elongation phase was not performed on these variants. Rather, SOD1WT and SOD1H46R did not seem to aggregate at all under our assay conditions (20 mM DTT, 5 mM EDTA 1 × PBS 37°C, shaking agitation) or was not aggregating into fibrils. For the most part, our assays suggested that propensity to aggregate was dependent upon mutation, and when we plotted the relative final fluorescence of each SOD1 variant in the first isolated aggregation assay (Figure 2A) against the inclusion formation data from transfected mammalian cells (Figure 2D) we observed a relatively weak correlation (r2 = 0.6). However, if we omit SOD1G37R as an outlier, due to its interestingly high propensity to form inclusions in cells comparative to its patient survival time, the correlation is much stronger (r2 = 0.86). These data taken together suggest that in our case, with the exception of SOD1G37R, aggregation propensity of isolated SOD1 can predict its cellular aggregation.
Figure 2. Isolated recombinant SOD1 aggregation. Aggregation of purified recombinant SOD1 as promoted with 20 mM DTT 5 mM EDTA in 1 × PBS, with shaking at 37°C, and was monitored by ThT fluorescence. (A,B) Plate-based aggregation assay of recombinant SOD1 variants showing the differential aggregation profiles of each variant where each panel represents a separate assay under the same conditions. Dotted lines represent SEM of 4 technical replicates. (C) Analysis of the fibril elongation rate (measured as ΔThT fluorescence) of assays (A,B). SOD1WT and SOD1H46R are not included due to the sigmoidal fit being below r2 = 0.75. (D) Plotting the percentage of cells with inclusions against the relative maximum ThT signal from assay A gave an r2 value of 0.60 (solid black line) which was strengthened to 0.86 (dotted line) when SOD1G37R was omitted from analysis.
DTT and EDTA Cause Changes in SOD1 Oligomeric State
A proposed common characteristic of SOD1 variants is dimer destabilization and dissociation into monomers, which can subsequently lead to the formation of higher order aggregates in vitro (Nordlund and Oliveberg, 2008). Indeed, antibodies detecting epitopes exposed only in the misfolded SOD1 monomer show immunoreactivity in motor neurons from human tissue as well as transgenic rat and mouse (SOD1A4V, SOD1G93A, and SOD1G37R) models of disease (Rakhit et al., 2007). To determine if propensity to dissociate to monomer was related to cellular aggregation we used the same destabilizing conditions used above (30 μM SOD1 20 mM DTT, 5 mM EDTA in 1 × PBS) but examined oligomeric state after 2 h (i.e., prior to insoluble aggregates being formed). Native PAGE of our variants at a concentration of 15 μM resolved SOD1WT, SOD1H46R, and SOD1D90A into clear monomer and dimer populations (Figure 3). Interestingly, we observed no separate monomer or dimer populations for SOD1G37R, SOD1G93A, or SOD1V148G. Rather, the bands seem to indicate that a large proportion of the protein remained dimeric under aggregation inducing conditons. In particular, SOD1G37R showed a substantial difference compared to SOD1WT as there was a large proportion of protein that had a lower mobility when compared to the dimer band, suggesting the presence of larger oligomeric species. Native PAGE has been used previously to investigate the electrophoretic mobility of SOD1 variants following disulfide reduction (Tiwari and Hayward, 2003), where it was proposed that increased mobility was most likely a result of metal loss, and retarded mobility was most likely a result of a conformational change in the protein structure. We next used analytical SEC after incubation in DTT and EDTA, in order to further investigate oligomeric distribution. We observed that the elution profiles of the SOD1 variants corresponded closely to the mobility observed for each variant in the native PAGE assay (Figure 3), suggesting similar changes were detected by each technique. The dimer and monomer populations were, again, apparent for SOD1WT, SOD1H46R, and SOD1D90A, with the dimer eluting at ~12.5 mL, and the monomer at ~14.3 mL. However, the elution profile of the DTT/EDTA-treated SOD1 dimer was distinct to the primarily dimeric non-treated controls (~13.6 mL), indicating a change in the hydrodynamic volume of the SOD1 dimer. This shift to a lower elution volume was observed for all SOD1 variants, with SOD1G37R, SOD1G93A, and SOD1V148G having major elution peaks at ~12.5 ml, but no peak at ~14.3 mL in contrast to SOD1WT, indicating a lack of monomer. Interestingly, SOD1G37R showed substantial protein at the void volume indicating the presence of oligomeric species with masses greater than or equal to ~100 kDa (S-75 10/300 GL column). SOD1G93A, SOD1V148G, and to a lesser extent SOD1D90A showed the presence of oligomeric species with peak shoulders ranging from ~9 to 11 mL. The analytical SEC, in conjunction with the native PAGE analysis, suggest that after a 2 h treatment with DTT and EDTA, SOD1 undergoes a conformational change that affects the hydrodynamic volume and the packing of the protein. It is well known that changes to the structure of molecules can make SEC data ambiguous (Hong et al., 2012), potentially making the interpretation of these data complicated. Due to this, native mass spectrometry analysis was performed on fractions eluted from analytical SEC assays.
Figure 3. Analytical size-exclusion chromatography of SOD1 variants. SOD1 variants were incubated at 37°C, in 1 × PBS, with 20 mM DTT, and 5 mM EDTA, for 2 h. After incubation, samples were concentrated at 4°C to 450 μM in a volume of 100 μl and loaded onto an analytical S-75 column equilibrated in 200 mM NH4OAc, and eluted at 0.5 ml/min into 0.5 ml fractions. Black is control native SOD1 which was not incubated in the presence of DTT or EDTA. Blue is SOD1 + DTT + EDTA, which showed different elution profiles for all SOD1 variants, with peaks at ~12 mL being labeled “a” and the peak at ~14 mL being labeled “b.” Inset above each chromatogram is native PAGE of the 2 h incubation time point.
Putative SOD1 Dimer from SEC Is Unfolded Monomer
The native PAGE and analytical SEC data indicated the variable presence of SOD1 dimer and monomer populations depending on the variant. In addition, a potential change in dimer conformation was observed as the dimer peak had a lower elution volume after treatment. Native mass spectrometry (MS) can be used to simultaneously investigate the suite of potential species in this dynamic system, giving high-resolution mass data that can also be used to determine modification states (disulfide status, metal binding, oxidative modification). MS was used to analyse the main peaks eluting from the analytical SEC from Figure 3 (MS analysis in Figure 4), finding that in the case of SOD1WT, (Figures 4A,B) the oligomeric distribution in the mass spectra matched the oligomeric state determined from the SEC (i.e., the peak at 12.5 mL was predominantly dimer and the peak at 14.3 mL predominantly monomer). The charge state distributions of both monomer and dimer in this assay were similar to that observed previously for native SOD1 variants (dimer 9 – 12+, monomer 7 – 8+) (McAlary et al., 2013). Analytical SEC of SOD1G93A resulted in only one major peak from DTT/EDTA treatment at 12.5 mL (peak a), consistent with the elution of the DTT/EDTA-treated SOD1WT dimer peak (Figure 3). However, when this was analyzed via MS, it was observed that a substantially higher abundance of monomer ranging from charge states 7+ to 21+ was present (Figure 4C). This broad monomeric charge state distribution was also observed in the ~12 mL peak for SOD1G37R (Figure 4D, peak a), the ~14.2 mL peak for SOD1D90A (Supplementary Figure 1A), and the ~12.1 mL peak for SOD1V148G (Supplementary Figure 1B), although each variant comprised different relative abundances of the highly charged monomer ions (9 – 21+). The charge state distribution of proteins in mass spectrometry is dependent upon conformation (Prakash and Mazumdar, 2005), and the presence of these highly charged monomeric ions is indicative of a variably unfolded protein in the solution phase, due to greater access of solvent to chargeable sites during ionization. SOD1 contains 23 basic amino acids (His = 8, Lys = 11, Arg = 4), and when recombinantly expressed and purified from E. coli, is methionine cleaved but not acetylated at its N-terminus (Leinweber et al., 2004), meaning it can theoretically carry 23 positive charges if fully unfolded. Our data shows a maximum of 21 charges applied to SOD1G93A (Figure 4), and even in the case of SOD1G37R (Figure 4), which has 24 chargeable sites, we only observed a maximum of 21 charges, suggesting that 2 (3 for SOD1G37R) of the basic amino acids are masked from protonation. Taken together we conclude that the monomer ions ranging from 7 to 11+ have maintained some structure after DTT/EDTA treatment, and that 12 – 21+ ions are representative of a highly unfolded monomer.
Figure 4. Native mass spectrometry of peaks from analytical SEC. Fractions from under peak apexes of DTT/EDTA-treated SOD1 were subject to mass spectrometry under gentle conditions optimized to maintain a “solution-like” conformation. Oligomeric and conformational state is represented by ellipsoids and lines above peaks with charge states marked. (A) Spectra of peak “a” from SOD1WT showing a primarily dimeric conformation similar to what is observed in native SOD1 (McAlary et al., 2013). (B) Spectra of peak ‘b’ from SOD1WT showing a primarily monomeric conformation at charge states similar to native monomeric SOD1WT (McAlary et al., 2013). (C) Spectra of peak “a” from SOD1G93A showing a primarily monomeric distribution with peaks corresponding to monomeric SOD1 with substantially higher charge states than what was observed in SOD1WT. (D) Spectra of peak′a′ from SOD1G37R showing a similar oligomeric and conformational distribution to SOD1G93A with greater abundance of highly charged (9–21+) monomer.
SOD1G37R Has a High Propensity to Form Oligomeric Aggregates
The elution profile of SOD1G37R, from analytical SEC, showed a large abundance of protein in the void volume, indicating high molecular weight oligomeric species. In order to identify oligomeric species mass spectrometry was performed on the void volume fraction as well as a second analytical size exclusion analysis with freshly prepared DTT/EDTA treated SOD1G37R using an S-200 column (1.3 MDa exclusion limit). The S-200 SEC of SOD1G37R (Figure 5A) showed a primary elution peak for DTT/EDTA treated SOD1G37R at ~15.6 mL compared to the untreated control SOD1G37R at ~16.5 mL. This lower elution volume was consistent with the peaks observed in the S-75 SEC (Figure 3). The absorbance at 280 nm was broad ranging from 8 to 21 mL, with 8–14 mL corresponding to high molecular weight species, and 18–21 mL likely corresponding to residual DTT from the treatment (MS of this peak showed no protein signal–data not shown). To identify the main oligomeric species of the major peak (Figure 5A, peak a), mass spectrometry was performed on this fraction (Figure 5B), which showed the charge state distribution, of primarily monomeric SOD1G37R, similar to that observed from the S-75 SEC (Figure 4). The fractions at lower elution volumes, containing oligomeric SOD1G37R, were also subject to MS, however low signal was acquired, most likely due to the very low abundance of the oligomers present. To overcome this, the void volume fraction of the SOD1G37R S-75 SEC analysis was subject to MS. We observed a broad distribution of oligomers ranging from 1 to 12 subunits (Figure 5C), with the monomer being the most abundant species, indicating that the oligomers eluted from the column were highly labile. Generation of a theoretical mass spectrum containing the potentially identified species allowed for the correct assignment of each ion species (Figure 5D). We plotted the abundance of each oligomer to determine if there was a pattern of oligomerization, and found that the large oligomers were assembling in patterns of 4 subunits, as shown by the abundance of oligomers composed of 4, 8, and 12 subunits, and the absence of oligomers made from 7, 9, 10, or 11 subunits (Figure 5E). Since the monomer and dimer ions detected from this fraction were showing a charge state distribution typical for SOD1 (monomer = 7−8+, dimer = 9–12+), we plotted the average charge of each oligomeric state against the number of subunits in that oligomer, and found that the charge added per subunit addition was consistent at ~2 charges per subunit added to the complex, indicating that the oligomers were forming in an ordered fashion (Figure 5F).
Figure 5. Analysis of the void volume peak from DTT/EDTA-treated SOD1G37R. Owing to the observation that a large portion of SOD1G37R eluted in the void volume, this fraction was analyzed with native MS. (A) S-200 column chromatogram of native SOD1G37R (black) and DTT/EDTA-treated SOD1G37R (blue) showing differences in elution profile where DTT/EDTA-treated is eluting at lower volumes showing a broad oligomeric distribution. The main peak is labeled “a” and markers above the chromatogram represent molecular weight standards. (B) Spectra of peak “a” from the S-200 chromatogram showing a similar distribution of SOD1G37R to that seen previously in the S-75 column analysis. (C) Spectra of the void volume peak from the SOD1G37R S-75 column analysis showing a broad oligomeric distribution of SOD1G37R ranging from 1 to 12 subunits. (D) Theoretical spectra and deconvolution of the mass spectra observed in (C), showing where peak overlap may occur between oligomers, as well as the major charge state of each oligomer. (E) Graphing the relative abundance of each oligomer shows a pattern whereby oligomers are formed as groups of 4 subunits due to the relatively low abundance of 2–3 and 5–7 mers and no 9–11 mers. (F) Plotting the average charge of each oligomer against the number of subunits in that oligomer gives are strong r2 value of 0.99, indicating that the association of these oligomers proceeds in a ordered fashion.
Destabilised Monomer Abundance Predicts Cellular Aggregation
Given the identification of various states of unfolded monomer from SEC peaks we next wanted to perform a direct comparison of the relative abundances of charge states, and oligomeric distributions for each variant. This was achieved by analyzing the total protein content of each experiment to allow for comparison. Here we observed, that for each SOD1 variant after DTT and EDTA treatment the major species was monomeric, ranging from charges of 7+ to 21+ (Figure 6), as previously observed (Figures 4, 5). The dimer abundance was substantially decreased in all cases, compared to previous spectra that were acquired from peak a in analytical SEC. This indicates that the pooling and dilution of the samples from analytical SEC resulted in dissociation of the dimer. In the case of SOD1G37R (Figure 6B), there was no evidence of any oligomeric species, even though fractions thought to contain them were included into the pooled sample. This is most likely due to a combination of the biases involved in mass spectrometry based analyses, such as this, where system parameters determine sensitivity at specific mass ranges, and the previously observed labile nature of the oligomers. In contrast to previously acquired spectra from fractions corresponding to peaks in the analytical SEC (Figures 4, 5), it was observed that highly charged monomer was present for all SOD1 variants, although at different abundances for each SOD1 variant. SODWT, SOD1H46R, and SOD1D90A (Figures 6A,C,D, respectively) showed relatively low abundances of highly charged SOD1 monomer, instead the major ion species was the monomeric 8+ charge state. SOD1G37R, SOD1G93A, and SOD1V148G (Figures 6B,E,F, respectively) had considerably increased relative abundances of highly charged monomeric species, in comparison to the typical 7+ and 8+ monomeric charge states.
Figure 6. Mass spectra of pooled DTT/EDTA-treated SOD1 variants from SEC. Fractions across the elution profile containing SOD1 were pooled, diluted to a concentration of 10 μM, and subject to mass spectrometry. (A) SOD1WT showing a high abundance of monomer (gray) compared to dimer (black). Inset is an expansion of the high charged 10 – 21+ monomer ions, normalized to the 10+ ion intensity. (B) SOD1G37R spectrum showing that highly charged monomer is the major species present. (C) SOD1H46R showing a primarily low charged monomeric distribution. Inset is an expansion of the high charge region normalized to the 9+ ion intensity. (D) SOD1D90A showing that the most abundant ion species is the 8+ charge state with low abundance of highly charged monomer. Inset is an expansion of the high charged region, normalized to the intensity of the 16+ ion. (E) SOD1G93A showing a substantial abundance of high charged monomer. (F) SOD1V148G showing a substantial abundance of high charged monomer.
Following the analysis and identification of protein species evident from a close inspection of different monomer ions, it was determined that there were three main monomeric conformations that were being detected by MS. These conformations were termed compact, intermediate, and unfolded due to the major species across the monomeric ions being reduced apo-SOD1 for all variants i.e., no holo-SOD1 was observed. It was also based on the Gaussian distribution of peaks (Figure 7A) where we determined the relative abundance of each peak and plotted the abundance of each charge state. From these Gaussian fits it was possible to determine the relative abundance of each conformation for all SOD1 variants (Figure 7B). We found that SOD1WT, SOD1H46R, and SOD1D90A monomers primarily inhabited the compact conformation, with much lower abundances of intermediate and unfolded conformations. SOD1G37R, SOD1G93A, and SOD1V148G had far higher abundances of intermediate and unfolded monomer, compared to SOD1WT, and out of the three variants, only SOD1G37R showed a higher abundance of unfolded monomer compared to compact monomer. Since SOD1 aggregation has been previously proposed to be dependent upon the amount of unfolded/destabilized protein (Byrstrom et al., 2010), the relative ThT intensity from previous aggregation data (Figure 2A) was plotted against the percentage of monomer that was destabilzed (sum of intermediate and unfolded abundances) (Figure 7C). The correlation was relatively weak (r2 = 0.31), however when SOD1G37R was omitted from the linear fit, the correlation became much stronger (r2 = 0.70). SOD1G37R is a curious outlier as it is defined as “wild-type like” due to its structural similarities with SOD1WT (Hart et al., 1998) and has a long variable disease duration of ~17 years (Wang et al., 2008), yet in our assays presents severe structural destabilization on similar levels as the aggressive SOD1G93A and SOD1V148G variants (~2.2 and ~2.3 years patient survival, respectively). Given that aggregation propensity did not very well predict cellular aggregation we next asked whether the relative abundance of destabilized monomer correlated with cellular aggregation (Figure 7D). Here the correlation was very strong (r2 = 0.87) and as such it seems that the main predictor of cellular aggregation is the propensity of the SOD1 variants to unfold as previously suggested (Nordlund and Oliveberg, 2008).
Figure 7. Determination of the abundance of monomeric conformations detected by mass spectrometry. (A) Plotting the abundance of each charge state allowed for the determination of 3 Gaussian distributions that correspond to increasing charge states and decreasing peak width, indicative of 3 conformations (compact = blue, intermediate = green, unfolded = red). (B) The abundance of each specific conformation for each SOD1 variant shows the monomeric distribution where SOD1WT, SOD1H46R, and SOD1D90A present a primarily compact monomeric distribution, compared to the other variants which showed high abundances of intermediate and unfolded monomer conformations. Error bars represent SD of the mean from 3 separate DTT/EDTA treatments of 3 separate nano-electrosprays. (C) Plotting the percentage abundance of destabilized monomer (unfolded and intermediate conformation) against the relative ThT intensity from aggregation assay A (Figure 2) showed a weak correlation (r2 = 0.31), however omission of SOD1G37R from the linear fit gave a much stronger correlation (r2 = 0.70). (D) Plotting the percentage abundance of destabilized monomer against the percentage of cells containing inclusions after transfection showed a strong correlation (r2 = 0.87). Error bars represent SD of the mean from at least 3 separate experiments.
Unfolded SOD1G37R Does Not Template on to SOD1G93A Seeds
While we observe a strong correlation between propensity to unfold, cellular aggregation and resultant cell loss, translating this to human disease severity is limited by the fact that one of our variants SOD1G37R has a high propensity to aggregate and is toxic in cells but has a long mean disease duration of 17 years. The concepts of aggregate strains and prion-like seeding have recently become prominent within the study of neurodegenerative disease, especially relating to ALS (Grad et al., 2014; Bergh et al., 2015; Lang et al., 2015; Bidhendi et al., 2016). Owing to the fact that spinal cord material from SOD1G37R mice were unable to propagate disease in SOD1G85R mice while SOD1G93A material was (Ayers et al., 2014), we next performed an experiment to elucidate if this unfolded species could template onto other preformed SOD1 aggregates. SOD1G93A was chosen to be the ‘seed’ variant, as we had previously shown it could seed SOD1WT aggregation (Grad et al., 2014), in this experiment a minimal amount of preformed aggregate seed (seed was 0.3% total protein by mass) was added to the assay. Remarkably, it was observed that SOD1G37R did not respond to the addition of seed, whereas SOD1WT showed an appreciable response (Figures 8B,D) and SOD1G93A showed significantly altered aggregation kinetics (Figures 8A,D). SOD1G37R showed no difference in either fibril elongation rate (measured as change in ThT fluorescence) (Figures 8C,D), or lag phase (Figure 8D), indicating that the addition of the seed was not affecting either the nucleation time, or the generation of fibrillar structures. The kinetics of SOD1WT were difficult to determine due to a weak sigmoidal fit, however it was clear from the assay (Figure 8A) that the addition of SOD1G93A seed increased aggregation. Addition of seed to SOD1G93A aggregation significantly increased the elongation rate (Figure 8D) compared to no seed, and also significantly decreased the lag phase (Figure 8E), indicating that the addition of seed was pushing the aggregation reaction toward fibril elongation, overcoming the need to form aggregate nuclei. Next we compared the fold increase in the mean maximum ThT fluorescence between variants (Figure 8F) after seeding and found that SOD1G37R had only a slight increase in max fluorescence, compared to SOD1WT which was higher than SOD1G37R and SOD1G93A which was substantially higher than both other variants. Taken together, this work suggests that SOD1G37R may be aggregating as a different strain that is dependent upon the substitution of arginine at position-37.
Figure 8. SOD1G37R is resistant to seeding by preformed SOD1G93A aggregates. Aggregation of SOD1 was promoted with DTT, EDTA, and shaking, with the formation of fibrillar structures being monitored by ThT fluorescence. Preformed SOD1G93A aggregates were added (or not) at a final amount of 0.3% total protein. (A) Aggregation assay of SOD1G93A normalized to the maximum signal showing the ThT fluorescence over time for each sample. Dotted lines represent the maximum fluorescence of SOD1WT (red) and SOD1G37R (blue) assays shown in panels B and C, respectively. (B) Aggregation assay of SOD1WT with and without seed, normalized to the maximum fluorescent signal. (C) Aggregation assay of SOD1G37R with and without seed, normalized to the maximum fluorescence signal (D) Kinetic analysis gave the fibril elongation rate (measured as change in ThT fluorescence), showing that seeded SOD1G93A had a significantly increased fibril elongation rate compared to other treatments. SOD1WT is not shown here due to it not aggregating without seed. (E) Kinetic analysis gave the lag phase of aggregation, showing that addition of seed significantly decreased the lag phase of SOD1G93A but not SOD1G37R. (F) Comparing the fold increase in the max fluorescence (Fmax) values when aggregation was seeded for each variant, showing that SOD1G93A and SOD1WT had greater increases compared to SOD1G37R. Each assay presented here is representative of 2 separate assays and were performed in the same microplate where error bars represent SEM of 4 technical replicates. Statistical significance was determined using a one-way ANOVA with a Tukey's post-test (***P <0.001).
In this work we have shown that the reduction and chelation of SOD1 variants induces unfolding into multiple conformational states, and that the major predictor for cellular aggregation of SOD1 variants is the abundance of destabilized SOD1 states. It is well known that disruption of the metal ions and the intramolecular disulfide has severe effects upon SOD1 structural stability (Hayward et al., 2002; Tiwari and Hayward, 2003; Lynch et al., 2004; Banci et al., 2009; Svensson et al., 2010). Indeed, in the native conformation (metal bound, disulfide bonded form) most ALS-associated SOD1 variants have high melting temperatures, however, when reduced and metal-free, some variants, such as SOD1G93A, melt below physiological temperature (Furukawa and O'Halloran, 2005; Rodriguez et al., 2005). Taking this into account, it is clear that the conformational distribution we observed in MS assays is a result of the lower thermal stabilities of the SOD1 variants when in their reduced apo-form, which was induced via incubation in DTT and EDTA at physiological temperature. Take SOD1G93A for example, which in its reduced apo-form, has a reported melting temperature of 31.2 ± 0.3°C (Furukawa and O'Halloran, 2005). We observed that SOD1G93A was significantly destabilized following DTT/EDTA treatment at physiological temperature, in agreement with this previous work. The most interesting aspect of the MS-based assay however, was that it was able to further separate the products of destabilization, determining the existence and relative abundance, of several conformational states. Previously published thermal stability data has not detected these states (Furukawa and O'Halloran, 2005; Rodriguez et al., 2005), most likely due to the low resolution provided by the technique, or that these conformations are in a constant state of flux. Higher resolution analyses of disordered SOD1 have been conducted using NMR, where reduced apo-SOD1G93A was found to be highly unstructured, with little propensity to form secondary structure (Luchinat et al., 2014), and reduced apo-SOD1WT was found to assume a monomeric native-like β-barrel structure with disordered loops (Sekhar et al., 2015). This model of reduced apo-SOD1WT matches the analysis presented here where the relatively low charge state distribution is indicative of a more compact monomeric species, close to what we observed for disulfide oxidized holo-SOD1WT.
Likewise, other work has considered reduced apo-SOD1G93A to be representative of several other SOD1 variants, including SOD1G37R, which is in partial agreement with our results (Luchinat et al., 2014). SOD1G37R was the most destabilized variant observed, with SOD1G93A a close second. In-cell NMR spectroscopy revealed that the levels of unstructured reduced apo-SOD1G37R and reduced apo-SOD1G93A were similar when overexpressed in mammalian cells, failing to detect any tertiary structure (Luchinat et al., 2014). In contrast, the MS-based assay performed here detected reduced apo-SOD1 species inhabiting quaternary, tertiary, and unstructured states, as determined by IM-MS, across all of our SOD1 variants. The discrepancy is potentially due to the experimental design, i.e., the work present here has investigated SOD1 that has been unfolded from its native conformation, whereas Luchinat et al. (2014) and colleagues have investigated SOD1 structures that may not yet be folded. Regardless of the order, the removal of post-translational modifications from SOD1 is widely thought to result in a product that is similar to the immature polypeptide that can take part in aggregation.
The isolated protein and cell-based aggregation assays showed a clear difference between specific SOD1 variant aggregation propensities. For the most part, aggregation in the cellular context correlated with aggregation of isolated protein, with the exception of SOD1G37R. Ivanova et al. (2014) determined that the SOD1 sequence segment 30KVWGSIKGL38 was a fibril forming segment but the same peptide containing the SOD1G37R mutation did not form fibrillar structures, suggesting that the SOD1G37R mutation may lead to aggregation through other fibril forming segments or to amorphous aggregation. SOD1G37R was a consistent outlier throughout this work, and even more so when taking into account the interesting correlation between patient survival time and SOD1 mutation, where in general SOD1 variants more prone to aggregation are correlated with a more rapid disease progression (Wang et al., 2008). Given the association with a long disease course, one might expect SOD1G37R to have a low aggregation propensity in cells, as well as a low isolated protein aggregation propensity. However, its in-cell aggregation was significantly increased compared to SOD1WT, similar to that observed for SOD1G93A and SOD1V148G. A potential explanation for this is the difference between the intracellular environment compared to that of an isolated protein assay. The isolated protein assay contains only pure SOD1 protein at a relatively low concentration (uncrowded) of 30 μM dimer, and is optimized to produce fibrillar aggregates at a rapid rate with vigorous shaking. The cytoplasm contains an abundance of protein and other macromolecules in high concentration (~400 mg/mL), resulting in a crowded environment (Fulton, 1982; Zeskind et al., 2007). This crowding effect, along with other cellular processes, are likely changing the dynamics of SOD1 aggregation, when compared to isolated protein. Indeed, transfection of cells typically results in variable expression levels of the protein of interest, across the population of cells that have taken up the vector. In the case presented here, the highest expressing cells are the ones most likely to develop inclusions, consistent with the idea that protein concentration is an important factor contributing to aggregation. Inducing high expression levels of SOD1 in the absence of similarly high expression of the copper-chaperone for SOD1 (CCS) is likely to result in a large pool of reduced apo-SOD1 that is more susceptible to aggregation. Recent work has highlighted a faulty interaction between SOD1A4V and human CCS (Wright et al., 2016). This suggests that SOD1 aggregation may result from an overloading of its specific quality control systems in motor neurons. Another explanation is the energy driven aggregation that occurs in cells. Previous work has shown that large inclusions (similar to the inclusions we examine here) are dependent upon microtubule transport (Farrawell et al., 2015) and are not a simply a product of random aggregation. It is possible that aggregation occurs after misfolded SOD1 has been transported and concentrated at the site of an inclusion.
The analytical SEC and MS analysis lends credence to the concept of crowding, where destabilized SOD1 was concentrated prior to SEC (~450 μM), however, the presence of a substantial amount of high molecular weight species was observed only in the case of SOD1G37R, indicating a rapid formation of oligomeric species for this variant. These oligomeric species were sensitive to dilution, as evidenced by the most abundant species in our MS analysis of the SOD1G37R void volume being monomeric. This potentially translates into a greater cellular aggregation propensity due to the high concentration of SOD1 with certain cells, but a lower propensity of isolated protein to aggregate due to the comparatively decreased concentration of SOD1 in our assay. Furthermore, there is a possibility that SOD1G37R forms aggregate nuclei rapidly, leading to a decreased pool of available monomer to template onto extending fibrils, although we cannot rule out the possibility of a competing amorphous aggregation pathway for SOD1G37R. Interestingly, while in many neurodegenerative diseases soluble hydrophobic oligomeric aggregates are implicated as the most toxic aggregate species (Bolognesi et al., 2010), here we find the opposite; SOD1G37R has the highest propensity to form oligomers, but is the amongst the slowest progressing mutants in humans.
Recent work has emphasized the role that different aggregation pathways, leading to the formation of different aggregate strains, may play in the role of amyotrophic lateral sclerosis (Bergh et al., 2015; Lang et al., 2015; Bidhendi et al., 2016). Using binary epitope mapping, Bergh et al. (2015) determined that mice overexpressing human SOD1D90A developed SOD1 aggregates that were found to be competing strains termed strain A and B. Each strain demonstrated differential kinetics and mortality in mice. A follow on from this saw the injection of mice expressing human SOD1WT with either of these strains developing ALS-like symptom at different rates, and also presented different pathophysiology (Bidhendi et al., 2016). The strain hypothesis may explain, at least in part, why different SOD1 mutants present variable disease progression in patients. Indeed, the work presented here indicates that SOD1G37R may be aggregating into a strain distinct from SOD1G93A similar to previous work that suggested SOD1G37R mice spinal cord homogenates could not transmit disease to a SOD1G85R mouse, but that SOD1G93A mice homogenates could transmit disease (Ayers et al., 2014). Given the MS data showed that SOD1G37R had the highest abundance of unfolded monomer, it might be predicted to have substantially altered aggregation dynamics in the presence of preformed aggregate seed (decreased lag phase, increased elongation rate), however this was not the case. Previous work has shown that the SOD1 sequence segment 30KVWGSIKGL38 is fibril forming (Ivanova et al., 2014), and in conjunction with our observations presented here suggests it is possible that SOD1G93A fibrils are polymerizing based on this segment, and therefore the seeded aggregation of SOD1G37R did not take place due to the arginine substitution at position 37. This combined with the observation of high levels of unfolded SOD1G37R in the cell-free assays would infer SOD1G37R fibril formation is dependent upon other sequence segments, however, we still observe a low in vitro aggregation propensity for this SOD1 variant, indicating that the sequence segment 30KVWGSIKGL38 may indeed be the primary fibril forming segment within SOD1. This SOD1G37R aggregate strain may be, for the most part, benign, explaining the comparatively long disease duration of this SOD1 variant, although further characterization of in vitro and in vivo SOD1G37R aggregates would need to be performed to confirm this interpretation. Interestingly, of the variables measured here, including in vitro aggregation, cellular aggregation, structural destabilization, and cell death, the strongest correlation to human SOD1-associated fALS patient survival was in vitro aggregation (Figure 9A). Cellular inclusion formation, cell death, and structural destabilization all had very weak correlations with patient survival (Figures 9B,C,D) again consistent with the idea that a strain of aggregate reactive with ThT might best predict disease severity.
Figure 9. Correlation of measured variables to human disease duration. Since sufferers of SOD1-associated fALS cases have variable survival times post-diagnosis, the measured variables in this work were plotted against mean disease duration. (A) The relative maximum ThT signal from in vitro aggregation assay plotted against disease duration gave a plausible correlation (r2 = 0.77). (B) The percentage of cells containing inclusions 48 h post-transfection plotted against disease duration gave a very weak correlation (r2 = 0.16). (C) Plotting the cell death comparative to SOD1WT 48 h post-transfection gave a very weak correlation (r2 = 0.01). (D) Plotting the relative abundance of unfolded SOD1 measured from MS assays against disease duration gave a very weak correlation (r2 = 0.03).
In conclusion, the work presented here underscores the necessity of combining mass analysis techniques to size determination techniques, as the relationship between size and mass may be drastically altered by conformational changes to molecules. Native mass spectrometry is an excellent technique as it provides a means by which to not only separate different conformers, but identify modification states, as we have shown above. The conformers we have identified related to different monomeric conformational states of several ALS-associated SOD1 variants, showing that disulfide reduction and metal chelation resulted in differential conformational states, which correlates with cellular aggregation and toxicity. SOD1G37R presented itself as an outlier when taking into account its high abundance of unfolded protein and oligomeric species, yet a comparatively low in vitro aggregation propensity. Seeding SOD1G37R aggregation with preformed SOD1G93A aggregates had no effect on its aggregation kinetics, suggesting that SOD1G37R may aggregate as a different strain and thus potentially explain its slow disease progression. Understanding the interplay between unfolding and subsequent aggregation, as well as the potential for the formation of variable aggregate strains of more SOD1 variants may shed light on the phenomenon of SOD1 variant related differential patient survival times, and may also provide information on the disease causing species in SOD1-associated ALS.
LM performed experiments, interpreted and analyzed data, as well as wrote the initial manuscript. JA designed experiments, interpreted and analyzed data. JJY designed experiments, interpreted and analyzed data, and wrote and edited the manuscript.
JJY is supported by an NHMRC Career Development Fellowship (1084144) and a Dementia Teams Grant (1095215). LM was supported by a University of Wollongong matching scholarship.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank the MSURRF group for access to mass spectrometry facilities. We thank Prof. Neil Cashman and Prof. Jenny Beck for insightful comments and review of material presented here. We would also like to thank Prof. Mikael Oliveberg for kindly providing SOD1 in vitro expression plasmids.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnins.2016.00499/full#supplementary-material
Abdolvahabi, A., Shi, Y., Chuprin, A., and Shaw, B. F. (2016). Stochastic formation of fibrillar and amorphous superoxide dismutase oligomers linked to amyotrophic lateral sclerosis. ACS Chem. Neurosci. 7, 799–810. doi: 10.1021/acschemneuro.6b00048
Abel, O., Shatunov, A., Jones, R. A., Andersen, M. P., Powell, F. J., and Al-Chalabi, A. (2013). Development of a smartphone app for a genetics website: the amyotrophic lateral sclerosis online genetics database (ALSoD). JMIR Mhealth Uhealth 1, e18. doi: 10.2196/mhealth.2706
Arnesano, F., Banci, L., Bertini, I., Martinelli, M., Furukawa, Y., and O'Halloran, T. V. (2004). The unusually stable quaternary structure of Human Cu,Zn-superoxide Dismutase 1 is controlled by both metal occupancy and disulfide status. J. Biol. Chem. 279, 47998–48003. doi: 10.1074/jbc.M406021200
Ayers, J. I., Fromholt, S. E., O'Neal, V. M., Diamond, J. H., and Borchelt, D. R. (2016). Prion-like propagation of mutant SOD1 misfolding and motor neuron disease spread along neuroanatomical pathways. Acta Neuropathol. 131, 103–114. doi: 10.1007/s00401-015-1514-0
Ayers, J. I., Fromholt, S., Koch, M., DeBosier, A., McMahon, B., Xu, G., et al. (2014). Experimental transmissibility of mutant SOD1 motor neuron disease. Acta Neuropathol. 128, 791–803. doi: 10.1007/s00401-014-1342-7
Banci, L., Bertini, I., Boca, M., Calderone, V., Cantini, F., Girotto, S., et al. (2009). Structural and dynamic aspects related to oligomerization of apo SOD1 and its mutants. Proc. Natl. Acad. Sci. 106, 6980–6985. doi: 10.1073/pnas.0809845106
Bergh, J., Zetterström, P., Andersen, P. M., Brännström, T., Graffmo, K. S., Jonsson, P, A., et al. (2015). Structural and kinetic analysis of protein-aggregate strains in vivo using binary epitope mapping. Proc. Natl. Acad. Sci. 112, 4489–4494. doi: 10.1073/pnas.1419228112
Bidhendi, E. E., Bergh, J., Zetterstrom, P., Andersen, P. M., Marklund, S. L., and Brännström, T. (2016). Two superoxide dismutase prion strains transmit amyotrophic lateral sclerosis-like disease. J. Clin. Invest. 126, 2249–2253. doi: 10.1172/JCI84360
Bolognesi, B., Kumita, J. R., Barros, T. P., Esbjorner, E. K., Luheshi, L. M., Crowther, D. C., et al. (2010). ANS binding reveals common features of cytotoxic amyloid species. ACS Chem. Biol. 5, 735–740. doi: 10.1021/cb1001203
Bosco, D. A., Morfini, G., Karabacak, N. M., Song, Y., Gros-Louis, F., Pasinelli, P., et al. (2010). Wild-type and mutant SOD1 share an aberrant conformation and a common pathogenic pathway in ALS. Nat. Neurosci. 13, 1396–1403. doi: 10.1038/nn.2660
Byrström, R., Andersen, P., Gröbner, G., and Oliveberg, M. (2010). SOD1 mutations targeting surface hydrogen bonds promote amyotrophic lateral sclerosis without reducing apo-state stability. J. Biol. Chem. 285, 19544–19552. doi: 10.1074/jbc.M109.086074
Cashman, N. R., Durham, H. D., Blusztajn, J. K., Oda, K., Tabira, T., Shaw, I. T., et al. (1992). Neuroblastoma x spinal cord (NSC) hybrid cell lines resemble developing motor neurons. Dev. Dynamics 194, 209–221. doi: 10.1002/aja.1001940306
Cirulli, E. T., Lasseigne, B. N., Petrovski, S., Sapp, P. C., Dion, P. A., Leblond, C. S., et al. (2015). Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways. Science 347, 1436–1441. doi: 10.1126/science.aaa3650
Cox, D., Selig, E., Griffin, M. D. W., Carver, J. A., and Ecroyd, H. (2016). Small heat shock proteins prevent alpha-synuclein aggregation via transient interactions and their efficacy is affected by the rate of aggregation. J. Biol. Chem. 291, 22618–22629. doi: 10.1074/jbc.M116.739250
Ding, F., and Dokholyan, N. V. (2008). Dynamical roles of metal ions and the disulfide bond in Cu, Zn Superoxide dismutase folding and aggregation. Proc. Natl. Acad. Sci. U.S.A. 105, 19696–19701. doi: 10.1073/pnas.0803266105
Farrawell, N. E., Lambert-Smith, I. A., Warraich, S. T., Blair, I. P., Saunders, D. N., Hatters, D. M., et al. (2015). Distinct partitioning of ALS associated TDP-43, FUS and SOD1 mutants into cellular inclusions. Sci. Rep. 5:13416. doi: 10.1038/srep13416
Furukawa, Y., Anzai, I., Akiyama, S., Imai, M., Cruz, F., Tomohide, S., et al. (2015). Conformational disorder of the most Immature Cu,Zn-Superoxide dismutase leading to amyotrophic lateral sclerosis. J. Biol. Chem. 291, 4144–4155. doi: 10.1074/jbc.M115.683763
Furukawa, Y., Kaneko, K., Yamanaka, K., O'Halloran, T. V., and Nukina, N. (2008). Complete loss of post-translational modifications triggers fibrillar aggregation of SOD1 in the familial form of amyotrophic lateral sclerosis. J. Biol. Chem. 283, 24167–24176. doi: 10.1074/jbc.M802083200
Furukawa, Y., and O'Halloran, T. (2005). Amyotrophic lateral sclerosis mutations have the greatest destabilizing effect on the apo- and reduced form of SOD1, leading to unfolding and oxidative aggregation. J. Biol. Chem. 280, 17266–17274. doi: 10.1074/jbc.M500482200
Furukawa, Y., Torres, A. S., and O'Halloran, T. V. (2004). Oxygen-induced maturation of SOD1: a key role for disulfide formation by the Copper Chaperone (CCS). EMBO J. 23, 2872–2881. doi: 10.1038/sj.emboj.7600276
Grad, L. I., Yerbury, J. J., Turner, B. J., Guest, W. C., Pokrishevsky, E., O'Neill, M. A., et al. (2014). Intercellular propagated misfolding of wild-type Cu/Zn superoxide dismutase occurs via exosome-dependent and -independent mechanisms. Proc. Natl. Acad. Sci. U.S.A. 111, 3620–3625. doi: 10.1073/pnas.1312245111
Hart, P. J., Liu, H., Pellegrini, M., Nersissian, A. M., Gralla, E. B., Valentine, J. S., et al. (1998). Subunit asymmetry in the three-dimensional structure of a human CuZnSOD mutant found in familial amyotrophic lateral sclerosis. Protein Sci. 7, 545–555.
Hayward, L. J., Rodriguez, J. A., Kim, J. W., Tiwari, A., Goto, J. J., Cabelli, D. E., et al. (2002). Decreased metallation and activity in subsets of mutant superoxide dismutases associated with familial amyotrophic lateral sclerosis. J. Biol. Chem. 277, 15923–15931. doi: 10.1074/jbc.M112087200
Hong, P., Koza, S., and Bouvier, E. S. P. (2012). Size-exclusion chromatography for the analysis of protein biotherapeutics and their aggregates. J. Liq. Chromatogr. Relat. Technol. 35, 2923–2950. doi: 10.1080/10826076.2012.743724
Ivanova, M. I., Sievers, S. A., Guenther, E. L., Johnson, L. M., Winkler, D. D., Galaleldeen, A., et al. (2014). Aggregation-triggering segments of SOD1 fibril formation support a common pathway for familial and sporadic ALS. Proc. Natl. Acad. Sci. U.S.A. 111, 197–201. doi: 10.1073/pnas.1320786110
Khare, S. D., Caplow, M., and Dokholyan, N. V. (2006). FALS Mutations in Cu, Zn superoxide dismutase destabilize the dimer and increase dimer dissociation propensity: a large-scale thermodynamic analysis. Amyloid 13, 226–235. doi: 10.1080/13506120600960486
Lang, L., Zetterstrom, P., Brannstrom, T., Marklund, S., Danielsson, J., and Oliveberg, M. (2015). SOD1 aggregation in ALS mice shows simplistic test tube behaviour. Proc. Natl. Acad. Sci. 112, 9878–9883. doi: 10.1073/pnas.1503328112
Leinweber, B., Barofsky, E., Barofsky, D. F., Ermilov, V., Nylin, K., and Beckman, J. S. (2004). Aggregation of ALS mutant superoxide dismutase expressed in Escherichia coli. Free Rad. Biol. Med. 36, 911–918. doi: 10.1016/j.freeradbiomed.2003.12.021
Lindberg, M. J., Tibell, L., and Oliveberg, M. (2002). Common denominator of Cu/Zn superoxide dismutase mutants associated with amyotrophic lateral sclerosis: decreased stability of the apo state. Proc. Natl. Acad. Sci. U.S.A. 99, 16607–16612. doi: 10.1073/pnas.262527099
Luchinat, E., Barbieri, L., Rubino, J. T., Kozyreva, T., Cantini, F., and Banci, L. (2014). In-cell NMR reveals potential precursor of toxic species from SOD1 fALS mutants. Nat. Commun. 5, 5502. doi: 10.1038/ncomms6502
Lynch, S. M., Boswell, S. A., and Colón, W. (2004). Kinetic stability of Cu/Zn superoxide dismutase is dependent on its metal ligands: implications for ALS. Biochemistry 43, 16525–16531. doi: 10.1021/bi048831v
McAlary, L., Yerbury, J. J., and Aquilina, J. A. (2013). Glutathionylation potentiates benign superoxide dismutase 1 variants to the toxic forms associated with amyotrophic lateral sclerosis. Sci. Rep. 3:3275. doi: 10.1038/srep03275
Mehmood, S., Allison, T. M., and Robinson, C. V. (2015). Mass spectrometry of protein complexes: from origins to applications. Annu. Rev. Phys. Chem. 66, 453–474. doi: 10.1146/annurev-physchem-040214-121732
Münch, C., and Bertolotti, A. (2010). Exposure of hydrophobic surfaces initiates aggregation of diverse ALS-Causing Superoxide Dismutase-1 mutants. J. Mol. Biol. 399, 512–525. doi: 10.1016/j.jmb.2010.04.019
Prakash, H., and Mazumdar, S. (2005). Direct correlation of the crystal structure of proteins with the maximum positive and negative charge states of gaseous protein ions produced by electrospray ionization. J. Am. Soc. Mass Spectrom. 16, 1409–1421. doi: 10.1016/j.jasms.2005.04.002
Prudencio, M., Hart, P. J., Borchelt, D. R., and Andersen, P. M. (2009). Variation in aggregation propensities among ALS-associated variants of SOD1: correlation to human disease. Hum. Mol. Genet. 18, 3217–3226. doi: 10.1093/hmg/ddp260
Rae, T. D., Torres, A. S., Pufahl, R., and O'Halloran, T. (2000). Mechanism of Cu,Zn-superoxide dismutase activation by the human metallchaperone hCCS. J. Biol. Chem. 276, 5166–5176. doi: 10.1074/jbc.M008005200
Rakhit, R., Robertson, J., Vande Velde, C., Horne, P., Ruth, D. M., Griffin, J., et al. (2007). An immunological epitope selective for pathological monomer-misfolded SOD1 in ALS. Nat. Med. 13, 754–759. doi: 10.1038/nm1559
Roberts, K., Zeineddine, R., Corcoran, L., Li, W., Campbell, I. L., and Yerbury, J. J. (2013). Extracellular aggregated Cu/Zn superoxide dismutase activates microglia to give a cytotoxic phenotype. Glia 61, 409–419. doi: 10.1002/glia.22444
Rodriguez, J. A., Shaw, B. F., Durazo, A., Sohn, S. H., Doucette, P. A., Nersissian, A. M., et al. (2005). Destabilization of apoprotein is insufficient to explain Cu,Zn-Superoxide dismutase-linked ALS pathogenesis. Proc. Natl. Acad. Sci. U.S.A. 102, 10516–10521. doi: 10.1073/pnas.0502515102
Rodriguez, J. A., Valentine, J. S., Eggers, D. K., Roe, J. A., Tiwari, A., Brown, R. H., et al. (2002). Familial amyotrophic lateral sclerosis-associated mutations decrease the thermal stability of distinctly metallated species of human Copper/Zinc Superoxide dismutase. J. Biol. Chem. 277, 15932–15937. doi: 10.1074/jbc.M112088200
Sekhar, A., Rumfeldt, J. A., Broom, H. R., Doyle, C. M., Bouvignies, G., Meiering, E. M., et al. (2015). Thermal fluctuations of immature SOD1 lead to separate folding and misfolding pathways. Elife 4, e07296. doi: 10.1038/srep31138
Senoo, Y., Katoh, K., Nakai, Y., Hashimoto, Y., Bando, K., and Teramoto, S. (1988). Activity and stability of recombinant human superoxide dismutase in buffer solutions and hypothermic perfusates. Acta Med. Okayama 42, 169–174.
Svensson, A. E., Bilsel, O., Kayatekin, C., Adefusika, J., Zitzewitz, J. A., and Matthews, C. R. (2010). Metal-free ALS variants of dimeric human Cu,Zn-Superoxide dismutase have enhanced populations of monomeric species. PLoS ONE 5:e10064. doi: 10.1371/journal.pone.0010064
Tiwari, A., and Hayward, L. J. (2003). Familial amyotrophic lateral sclerosis mutants of Copper/Zinc superoxide dismutase are susceptible to disulfide reduction. J. Biol. Chem. 278, 5984–5992. doi: 10.1074/jbc.M210419200
Turner, B. J., Atkin, J. D., Farg, M. A., Zang, D. W., Rembach, A., Lopes, E. C., et al. (2005). Impaired extracellular secretion of mutant superoxide dismutase 1 associates with neurotoxicity in familial amyotrophic lateral sclerosis. J. Neurosci. 25, 108–117. doi: 10.1523/JNEUROSCI.4253-04.2005
Wang, Q., Johnson, J. L., Agar, N. Y., and Agar, J. N. (2008). Protein aggregation and protein instability govern familial amyotrophic lateral sclerosis patient survival. PLoS Biol. 6:e170. doi: 10.1371/journal.pbio.0060170
Weisberg, S. J., Lyakhovetsky, R., Werdiger, A. C., Gitler, A. D., Soen, Y., and Kaganovich, D. (2012). Compartmentalization of superoxide dismutase 1 (SOD1G93A) aggregates determines their toxicity. Proc. Natl. Acad. Sci. U.S.A. 109, 15811–15816. doi: 10.1073/pnas.1205829109
Whiten, D. R., San Gil, R., McAlary, L., Yerbury, J. J., Ecroyd, H., and Wilson, M. R. (2016). Rapid flow cytometric measurement of protein inclusions and nuclear trafficking. Sci. Rep. 6:31138. doi: 10.1038/srep31138
Williams, K. L., Topp, S., Yang, S., Smith, B., Fifita, J. A., Warraich, S. T., et al. (2016). CCNF mutations in amyotrophic lateral sclerosis and frontotemporal dementia. Nat. Commun. 7:11253. doi: 10.1038/ncomms11253
Wright, G. S., Hasnain, S. S., and Grossmann, J. G. (2011). The structural plasticity of the human Copper Chaperone for SOD1: insights from combined size-exclusion chromatographic and SOlution X-Ray scattering studies. Biochem. J. 439, 39–44. doi: 10.1042/BJ20110948
Yerbury, J. J., Gower, D., Vanags, L., Roberts, K., Lee, J. A., and Ecroyd, H. (2013). The small heat shock proteins alphaB-crystallin and Hsp27 suppress SOD1 aggregation in vitro. Cell Stress Chaperones 18, 251–257. doi: 10.1007/s12192-012-0371-1
Yerbury, J. J., Ooi, L., Dillin, A., Saunders, D. N., Hatters, D. M., Beart, P. M., et al. (2016). Walking the tightrope: proteostasis and neurodegenerative disease. J. Neurochem. 137, 489–505. doi: 10.1111/jnc.13575
Zeineddine, R., Pundavela, J. F., Corcoran, L., Stewart, E. M., Do-Ha, D., Bax, M., et al. (2015). SOD1 protein aggregates stimulate macropinocytosis in neurons to facilitate their propagation. Mol. Neurodegener. 10, 57. doi: 10.1186/s13024-015-0053-4
Keywords: ALS, SOD1, protein misfolding, protein aggregation, mass spectrometry
Citation: McAlary L, Aquilina JA and Yerbury JJ (2016) Susceptibility of Mutant SOD1 to Form a Destabilized Monomer Predicts Cellular Aggregation and Toxicity but Not In vitro Aggregation Propensity. Front. Neurosci. 10:499. doi: 10.3389/fnins.2016.00499
Received: 28 September 2016; Accepted: 20 October 2016;
Published: 04 November 2016.
Edited by:Wendy Noble, King's College London, UK
Reviewed by:Patrizia Longone, Fondazione Santa Lucia (IRCCS), Italy
Ian Paul Johnson, University of Adelaide, Australia
Copyright © 2016 McAlary, Aquilina and Yerbury. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Justin J. Yerbury, firstname.lastname@example.org