Changes in hydrophobicity mainly promotes the aggregation tendency of ALS associated SOD1 mutants

Delineating the Mechanistic Insight of Inhibition of α-Synuclein Fibrillation by Neuro Metabolite, Myo-inositol: Implications in Synucleopathies-Related Disorders

The fibrillation of α-synuclein (α-syn) is a major factor contributing to neuronal damage and is critical in developing synucleopathies-related disorders. Considering this, the discovery of new compounds that can inhibit or modulate α-syn aggregation is a significant area of research. While polyol osmolytes have been shown to reduce α-syn fibrillation, the impact of brain metabolites such as myo-inositol (MI) on α-syn aggregation has not yet been explored. This study is the first to examine the effects of MI on α-syn aggregation, utilizing spectroscopic, microscopic, and cell cytotoxicity assay. Various aggregation assays revealed that MI inhibits the α-syn fibrillation in a dose-dependent manner. Fluorescence microscopy observations suggest that MI inhibits the α-syn fibrillation by forming amorphous aggregates. MTT assay revealed that α-syn aggregates in the presence of different concentrations of MI were not toxic as compared to α-syn fibrils. Thus, the mechanistic insight of inhibition of α-syn fibrillation by MI was explored by employing interaction studies using spectroscopic, calorimetric, and in silico approaches. Surface plasmon resonance and isothermal titration calorimetry suggest that MI-α-syn interacted with significant binding affinity, and the reaction was spontaneous. Molecular docking results depict that MI interacted with the aggregation-prone residues (36-42) at the N-terminal of α-syn, thereby stabilizing the α-syn and preventing the fibril formation. Molecular dynamics simulation results demonstrate the stability of the complex formation of MI with α-syn. This study highlighted the mechanistic insight of MI on preventing the α-syn from forming amyloid fibril, which could be further explored for therapeutic management of synucleopathies-related disorders.

A random forest-based predictive model for classifying BRCA1 missense variants: a novel approach for evaluating the missense mutations effect

The right classification of variants is the key to pre-symptomatic detection of disease and conducting preventive actions. Since BRCA1 has a high incidence and penetrance in breast and ovarian cancers, a high-performance predictive tool can be employed to classify the clinical significance of its variants. Several tools have previously been developed for this purpose which poorly classify the significance in specific cases. The proposed tools commonly assign a score without providing any interpretation behind it. To reach an accurate predictive tool with interpretation abilities, in this study, we propose BRCA1-Forest which works based on random forest as a well-known machine learning technique for making interpretable decisions with high specificity and sensitivity in variants classification. The method involves narrowing down available options until reaching the final decision. To this end, a set of BRCA1 benign and pathogenic missense variants was collected first, and then, the dataset was prepared based on the effect of each variant on the protein sequence. The dataset was enriched by adding physicochemical changes and the conservation score of the amino acid position as pathogenicity criteria. The proposed model was trained based on the dataset to classify the clinical significance of variants. The performance of BRCA1-Forest was compared to four state-of-the-art methods, SIFT, PolyPhen2, CADD, and DANN, in terms of different evaluation metrics including precision, recall, false positive rate (FPR), the area under the receiver operator curve (AUC ROC), the area under the precision-recall curve (AUC-PR), and Mathew correlation coefficient (MCC). The results reveal that the proposed model outperforms the abovementioned tools in all metrics except for recall. The software of BRCA1-Forest is available at https://github.com/HamedKAAC/BRCA1Forest .

Structural insights into SOD1: from in silico and molecular dynamics to experimental analyses of ALS-associated E49K and R115G mutants

Protein stability is a crucial characteristic that influences both protein activity and structure and plays a significant role in several diseases. Cu/Zn superoxide dismutase 1 (SOD1) mutations serve as a model for elucidating the destabilizing effects on protein folding and misfolding linked to the lethal neurological disease, amyotrophic lateral sclerosis (ALS). In the present study, we have examined the structure and dynamics of the SOD1 protein upon two ALS-associated point mutations at the surface (namely, E49K and R115G), which are located in metal-binding loop IV and Greek key loop VI, respectively. Our analysis was performed through multiple algorithms on the structural characterization of the hSOD1 protein using computational predictions, molecular dynamics (MD) simulations, and experimental studies to understand the effects of amino acid substitutions. Predictive results of computational analysis predicted the deleterious and destabilizing effect of mutants on hSOD1 function and stability. MD outcomes also indicate that the mutations result in structural destabilization by affecting the increased content of β-sheet structures and loss of hydrogen bonds. Moreover, comparative intrinsic and extrinsic fluorescence results of WT-hSOD1 and mutants indicated structural alterations and increased hydrophobic surface pockets, respectively. As well, the existence of β-sheet-dominated structures was observed under amyloidogenic conditions using FTIR spectroscopy. Overall, our findings suggest that mutations in the metal-binding loop IV and Greek key loop VI lead to significant structural and conformational changes that could affect the structure and stability of the hSOD1 molecule, resulting in the formation of toxic intermediate species that cause ALS.

A Practical Guide to Computational Tools for Engineering Biocatalytic Properties

The growing demand for efficient, selective, and stable enzymes has fueled advancements in computational enzyme engineering, a field that complements experimental methods to accelerate enzyme discovery. With a plethora of software and tools available, researchers from different disciplines often face challenges in selecting the most suitable method that meets their requirements and available starting data. This review categorizes the computational tools available for enzyme engineering based on their capacity to enhance the following specific biocatalytic properties of biotechnological interest: (i) protein-ligand affinity/selectivity, (ii) catalytic efficiency, (iii) thermostability, and (iv) solubility for recombinant enzyme production. By aligning tools with their respective scoring functions, we aim to guide researchers, particularly those new to computational methods, in selecting the appropriate software for the design of protein engineering campaigns. De novo enzyme design, involving the creation of novel proteins, is beyond this review's scope. Instead, we focus on practical strategies for fine-tuning enzymatic performance within an established reference framework of natural proteins.

Fighting against amyotrophic lateral sclerosis (ALS) with flavonoids: a computational approach to inhibit superoxide dismutase (SOD1) mutant aggregation

Protein aggregation is a biological process that occurs when proteins misfold. Misfolding and aggregation of human superoxide dismutase (hSOD1) cause a neurodegenerative disease called amyotrophic lateral sclerosis (ALS). Among the mutations occurring, targeting the E21K mutation could be a good choice to understand the pathological mechanism of SOD1 in ALS, whereof it significantly reduces life hopefulness in patients. Naturally occurring polyphenolic flavonoids have been suggested as a way to alleviate the amyloidogenic behavior of proteins. In this study, computational tools were used to identify promising flavonoid compounds that effectively inhibit the pathogenic behavior of the E21K mutant. Initial screening identified Pelargonidin, Curcumin, and Silybin as promising leads. Molecular dynamics (MD) simulations showed that the binding of flavonoids to the mutated SOD1 caused changes in the protein stability, hydrophobicity, flexibility, and restoration of lost hydrogen bonds. Secondary structure analysis indicated that the protein destabilization and the increased propensity of β-sheet caused by the mutation were restored to the wild-type state upon binding of flavonoids. Free energy landscape (FEL) analysis was also used to differentiate aggregation, and results showed that Silybin followed by Pelargonidin had the most therapeutic efficacy against the E21K mutant SOD1. Therefore, these flavonoids hold great potential as highly effective inhibitors in mitigating ALS's fatal and insuperable effects.Communicated by Ramaswamy H. Sarma.

Implications of ALS-Associated Mutations on Biochemical and Biophysical Features of hSOD1 and Aggregation Formation

One of the recognized motor neuron degenerative disorders is amyotrophic lateral sclerosis (ALS). By now, several mutations have been reported and linked to ALS patients, some of which are induced by mutations in the human superoxide dismutase (hSOD1) gene. The ALS-provoking mutations are located throughout the structure of hSOD1 and promote the propensity to aggregate. Despite numerous investigations, the underlying mechanism related to the toxicity of mutant hSOD1 through the gain of a toxic function is still vague. We surveyed two mutant forms of hSOD1 by removing and adding cysteine at positions 146 and 72, respectively, to investigate the biochemical characterization and amyloid formation. Our findings predicted the harmful and destabilizing impact of two SOD1 mutants using multiple programs. The specific activity of the wild-type form was about 1.42- and 1.92-fold higher than that of C146R and G72C mutants, respectively. Comparative structural studies using CD spectropolarimetry, and intrinsic and ANS fluorescence showed alterations in secondary structure content, exposure of hydrophobic patches, and structural compactness of WT-hSOD1 vs. mutants. We demonstrated that two mutants were able to promote amyloid-like aggregates under amyloid induction circumstances (50-mM Tris-HCl pH 7.4, 0.2-M KSCN, 50-mM DTT, 37 °C, 190 rpm). Monitoring aggregates were done using an enhancement in thioflavin T fluorescence and alterations in Congo red absorption. The mutants accelerated fibrillation with subsequently greater fluorescence amplitude and a shorter lag time compared to WT-SOD1. These findings support the aggregation of ALS-associated SOD1 mutants as an integral part of ALS pathology.

Mutation/metal deficiency in the "electrostatic loop" enhanced aggregation process in apo/holo SOD1 variants: implications for ALS diseases

Despite the many mechanisms it has created to prevent unfolding and aggregation of proteins, many diseases are caused by abnormal folding of proteins, which are called misfolding diseases. During this process, proteins undergo structural changes and become stable, insoluble beta-sheet aggregates called amyloid fibrils. Mutations/disruptions in metal ion homeostasis in the ALS-associated metalloenzyme superoxide dismutase (SOD1) reduce conformational stability, consistent with the protein aggregation hypothesis for neurodegenerative diseases. However, the exact mechanism of involvement is not well understood. Hence, to understand the role of mutation/ metal deficiency in SOD1 misfolding and aggregation, we investigated the effects of apo/holo SOD1 variants on structural properties using biophysical/experimental techniques. The MD results support the idea that the mutation/metal deficiency can lead to a change in conformation. The increased content of β-sheet structures in apo/holo SOD1 variants can be attributed to the aggregation tendency, which was confirmed by FTIR spectroscopy and dictionary of secondary structure in proteins (DSSP) results. Thermodynamic studies of GdnHCl showed that metal deficiency/mutation/intramolecular S-S reduction together are required to initiate misfolding/aggregation of SOD1. The results showed that apo/holo SOD1 variants under destabilizing conditions induced amyloid aggregates at physiological pH, which were detected by ThT/ANS fluorescence, as well as further confirmation of amyloid/amorphous species by TEM. This study confirms that mutations in the electrostatic loop of SOD1 lead to structural abnormalities, including changes in hydrophobicity, reduced disulfide bonds, and an increased propensity for protein denaturation. This process facilitates the formation of amyloid/amorphous aggregates ALS-associated.

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Polyphenolic flavonoid compounds act as the inhibitory potential of aggregation process: Implications for the prevention and therapeutics against FALS-associated D101G SOD1 mutant

Aggregation of proteins is a biological phenomenon caused by misfolded proteins. Human superoxide dismutase (hSOD1) misfolding and aggregation underlie the neurological illness amyotrophic lateral sclerosis (ALS). The most significant contributing factor to ALS is genetic point mutations in SOD1. particularly, D101G mutant is the most harmful because it significantly reduces the life expectancy of patients. Subsequently, the use of natural polyphenolic flavonoids is strongly recommended to reduce the amyloidogenic behavior of protopathic proteins. In this study, using computational parameters such as protein-ligand interaction and molecular dynamics (MD) simulation analyses, we are trying to identify a pharmacodynamically promising flavonoid compound that can effectively inhibit the pathogenic behavior of the D101G mutant. Epigallocatechin-gallate (EGCG), Hesperidin, Isorhamnetin, and Diosmetin were identified as potential leads in a preliminary screening of flavonoids to anti-amyloid action. The results of MD showed that the binding of flavonoids to D101G mutant caused changes in stability, hydrophobicity of protein, and flexibility, as well as significantly led to the restoration of lost hydrogen bonds. Secondary structure analysis showed that protein destabilization and the increased propensity of β-sheet caused by the mutation were restored to the wild-type state upon binding of flavonoids. Besides, to differentiate aggregation, we elucidated alterations in the free energy landscape (FEL) and dynamic cross-correlation matrix (DCCM) of WT-SOD1 and mutant (unbound /bound) states. Among flavonoids, Epigallocatechin-gallate and Hesperidin had the most therapeutic efficacy against the D101G mutant. Therefore, Epigallocatechin-gallate and Hesperidin promise considerable therapeutic potential to develop highly effective inhibitors in reducing fatal and irreversible ALS.

Characterizing proteomic and transcriptomic features of missense variants in amyotrophic lateral sclerosis genes

Within recent years, there has been a growing number of genes associated with amyotrophic lateral sclerosis (ALS), resulting in an increasing number of novel variants, particularly missense variants, many of which are of unknown clinical significance. Here, we leverage the sequencing efforts of the ALS Knowledge Portal (3864 individuals with ALS and 7839 controls) and Project MinE ALS Sequencing Consortium (4366 individuals with ALS and 1832 controls) to perform proteomic and transcriptomic characterization of missense variants in 24 ALS-associated genes. The two sequencing datasets were interrogated for missense variants in the 24 genes, and variants were annotated with gnomAD minor allele frequencies, ClinVar pathogenicity classifications, protein sequence features including Uniprot functional site annotations, and PhosphoSitePlus post-translational modification site annotations, structural features from AlphaFold predicted monomeric 3D structures, and transcriptomic expression levels from Genotype-Tissue Expression. We then applied missense variant enrichment and gene-burden testing following binning of variation based on the selected proteomic and transcriptomic features to identify those most relevant to pathogenicity in ALS-associated genes. Using predicted human protein structures from AlphaFold, we determined that missense variants carried by individuals with ALS were significantly enriched in β-sheets and α-helices, as well as in core, buried or moderately buried regions. At the same time, we identified that hydrophobic amino acid residues, compositionally biased protein regions and regions of interest are predominantly enriched in missense variants carried by individuals with ALS. Assessment of expression level based on transcriptomics also revealed enrichment of variants of high and medium expression across all tissues and within the brain. We further explored enriched features of interest using burden analyses and identified individual genes were indeed driving certain enrichment signals. A case study is presented for SOD1 to demonstrate proof-of-concept of how enriched features may aid in defining variant pathogenicity. Our results present proteomic and transcriptomic features that are important indicators of missense variant pathogenicity in ALS and are distinct from features associated with neurodevelopmental disorders.

Computational identification and experimental validation of anti-filarial lead molecules targeting metal binding/substrate channel residues of Cu/Zn SOD1 from Wuchereria bancrofti

Lymphatic filariasis (LF) is a neglected mosquito-borne parasitic disease, widely caused by Wuchereria bancrofti (Wb) in tropical and sub-tropical countries. During a blood meal, the filarial nematodes are transmitted to humans by the infected mosquito. To counter attack the invaded nematodes, the human immune system produces reactive oxygen species. However, the anti-oxidant enzymes of nematodes counteract the host oxidative cytotoxicity. Cu/Zn Superoxide dismutase (SOD1), a member of antioxidant enzymes and are widely used by the nematodes to sustain the host oxidative stress across its lifecycle, hence targeting SOD1 to develop suitable drug molecules would help to overcome the problems related to efficacy and activity of drugs upon different stages of nematodes. In order to find the potent inhibitors, a three-dimensional structure of Cu/Zn WbSOD1 was modelled and the structural stability was analysed through simulation studies. The structure-guided virtual screening approach has been used to identify lead molecules from the ChemBridge based on the docking score, ADMET properties and protein-ligand complex stability analysis. The identified compounds were observed to interact with the copper, metal binding residues (His48, His63, His80 and His120) and catalytically important residue Arg146, which play a crucial role in the disproportionation of incoming superoxide radicals of Cu/Zn WbSOD1. Further, in vitro validation of the selected leads in the filarial worm Setaria digitata exhibited higher inhibition and better IC50 compared to the standard drug ivermectin. Thus, the identified leads could potentially inhibit enzyme activity, which could subsequently act as drug candidates to control LF.Communicated by Ramaswamy H. Sarma.

Self-Assembly of Platelet Lysates Proteins into Microparticles by Unnatural Disulfide Bonds for Bottom-Up Tissue Engineering

There is a demand to design microparticles holding surface topography while presenting inherent bioactive cues for applications in the biomedical and biotechnological fields. Using the pool of proteins present in human-derived platelet lysates (PLs), the production of protein-based microparticles via a simple and cost-effective method is reported, exploring the prone redox behavior of cysteine (Cy-SH) amino acid residues. The forced formation of new intermolecular disulfide bonds results in the precipitation of the proteins as spherical, pompom-like microparticles with adjustable sizes (15-50 µm in diameter) and surface topography consisting of grooves and ridges. These PL microparticles exhibit extraordinary cytocompatibility, allowing cell-guided microaggregates to form, while also working as injectable systems for cell support. Early studies also suggest that the surface topography provided by these PL microparticles can support osteogenic behavior. Consequently, these PL microparticles may find use to create live tissues via bottom-up procedures or injectable tissue-defect fillers, particularly for bone regeneration, with the prospect of working under xeno-free conditions.

SOD1-ALS-Browser: a web-utility for investigating the clinical phenotype in SOD1 amyotrophic lateral sclerosis

Objective: Variants in the superoxide dismutase (SOD1) gene are among the most common genetic causes of amyotrophic lateral sclerosis. Reflecting the wide spectrum of putatively deleterious variants that have been reported to date, it has become clear that SOD1-linked ALS presents a highly variable age at symptom onset and disease duration.Methods: Here we describe an open access web tool for comparative phenotype analysis in ALS: https://sod1-als-browser.rosalind.kcl.ac.uk/. The tool contains a built-in dataset of clinical information from 1383 people with ALS harboring a SOD1 variant resulting in one of 162 unique amino acid sequence alterations and from a non-SOD1 comparator ALS cohort of 13,469 individuals. We present two examples of analyses possible with this tool, testing how the ALS phenotype relates to SOD1 variants that alter amino acid residue hydrophobicity and to distinct variants at the 94th residue of SOD1, where six are sampled.Results and conclusions: The tool provides immediate access to the datasets and enables bespoke analysis of phenotypic trends associated with different protein variants, including the option for users to upload their own datasets for integration with the server data. The tool can be used to study SOD1-ALS and provides an analytical framework to study the differences between other user-uploaded ALS groups and our large reference database of SOD1 and non-SOD1 ALS. The tool is designed to be useful for clinicians and researchers, including those without programming expertise, and is highly flexible in the analyses that can be conducted.

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

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

Computational methods to predict protein aggregation

In most cases, protein aggregation stems from the establishment of non-native intermolecular contacts. The formation of insoluble protein aggregates is associated with many human diseases and is a major bottleneck for the industrial production of protein-based therapeutics. Strikingly, fibrillar aggregates are naturally exploited for structural scaffolding or to generate molecular switches and can be artificially engineered to build up multi-functional nanomaterials. Thus, there is a high interest in rationalizing and forecasting protein aggregation. Here, we review the available computational toolbox to predict protein aggregation propensities, identify sequential or structural aggregation-prone regions, evaluate the impact of mutations on aggregation or recognize prion-like domains. We discuss the strengths and limitations of these algorithms and how they can evolve in the next future.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Aberrant Stress Granule Dynamics and Aggrephagy in ALS Pathogenesis

Stress granules are conserved cytosolic ribonucleoprotein (RNP) compartments that undergo dynamic assembly and disassembly by phase separation in response to stressful conditions. Gene mutations may lead to aberrant phase separation of stress granules eliciting irreversible protein aggregations. A selective autophagy pathway called aggrephagy may partially alleviate the cytotoxicity mediated by these protein aggregates. Cells must perceive when and where the stress granules are transformed into toxic protein aggregates to initiate autophagosomal engulfment for subsequent autolysosomal degradation, therefore, maintaining cellular homeostasis. Indeed, defective aggrephagy has been causally linked to various neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). In this review, we discuss stress granules at the intersection of autophagy and ALS pathogenesis.