Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis

Bisdemethoxycurcumin, a novel potent polyphenolic compound, effectively inhibits the formation of amyloid aggregates in ALS-associated hSOD1 mutant (L38R)

Protein misfolding is a biological process that leads to protein aggregation. Anomalous misfolding and aggregation of human superoxide dismutase (hSOD1) into amyloid aggregates is a characteristic feature of amyotrophic lateral sclerosis (ALS), a neurodegenerative illness. Thus, focusing on the L38R mutant may be a wise decision to comprehend the SOD1 disease process in ALS. We suggest that Bisdemethoxycurcumin (BDMC) may be a strong anti-amyloidogenic polyphenol against L38R mutant aggregation. Protein stability, hydrophobicity, and flexibility were altered when BDMC was bound to the L38R mutant, as shown by molecular dynamic (MD) simulations and molecular docking. FTIR data shows α-Helix dominance in BDMC-containing samples, with reduced β-sheet and disordered peaks, indicating the decrease of aggregate species. ThT aggregation kinetics curves show BDMC reduces L38R mutant aggregation dose-dependently, with higher BDMC concentrations yielding greater reductions. TEM images showed various quantities of amorphous aggregates, but notably, 60 μM BDMC markedly reduced aggregate density, underscoring BDMC's inhibitory effect. Hemolysis tests revealed aggregate species in BDMC-treated samples were less toxic than in L38R mutant samples alone at the same concentrations and exposure times. Overall, BDMC has substantial potential to develop highly effective inhibitors that mitigate the risk of fatal ALS.

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.

Molecular docking and molecular dynamics simulation of wheat gluten-derived antioxidant peptides acting through the Keap1-Nrf2 pathway

BACKGROUND

In our previous study, we successfully identified five peptides from wheat gluten: Ala-Pro-Ser-Tyr (APSY), Leu-Tyr (LY), Pro-Tyr (PY), Arg-Gly-Gly-Tyr (RGGY) and Tyr-Gln (YQ). Molecular docking and molecular dynamics simulation methods were employed to investigate the interaction between these antioxidant peptides and the Kelch-like ECH-associated protein 1 (Keap1 protein), revealing the molecular mechanism of their non-competitive binding. In addition, the total antioxidant capacity of the five peptides was determined using the 2,2'-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) method.

RESULTS

The affinities of APSY, LY, PY, RGGY and YQ were -8.9, -8.3, -8.5, -9.1 and - 7.9 kcal mol-1, respectively. The five peptides effectively bound to Keap1 protein through hydrogen, π-σ, π-alkyl and alkyl interactions. Significant roles were observed for the P1 pocket residue ARG-415 and the P3 pocket residue ALA-556 in the interactions of the Keap1-peptide complexes. Molecular dynamics simulations further elucidated the dynamic process of peptide binding to the Keap1 protein. All five peptides formed stable complexes with Keap1 protein, with van der Waals forces playing crucial roles in these complex systems, indicative of the peptides' strong binding ability to Keap1 protein. The van der Waals forces were -178.74, -123.11, -134.36, -132.59, and -121.44 kJ mol-1 for the Keap1-APSY, Keap1-LY, Keap1-PY, Keap1-RGGY and Keap1-YQ complexes, respectively. These peptides exhibited excellent antioxidant effects. Among them, the YQ peptide exhibited the highest total antioxidant capacity, with an activity value of 1.18 ± 0.06 mmol Trolox equivalent (TE) L-1 at a concentration of 0.10 mg mL-1. The RGGY, PY, LY and APSY peptides followed in descending order, with activity values of 0.91 ± 0.05, 0.72 ± 0.06, 0.62 ± 0.04 and 0.60 ± 0.05 mmol TE L-1, respectively.

CONCLUSION

These results unveiled the molecular mechanism by which the five antioxidant peptides act on active pockets through the Keap1-Nrf2 signaling pathway, providing a theoretical basis for the development of antioxidants. © 2024 Society of Chemical Industry.

Computational screening of pathogenic missense nsSNPs in heme oxygenase 1 (HMOX1) gene and their structural and functional consequences

Heme Oxygenase 1 (HMOX1) is a cytoprotective enzyme, exhibiting the highest activity in the spleen, catalyzing the heme ring breakdown into products of biological significance- biliverdin, CO, and Fe2+. In vascular cells, HMOX1 possesses strong anti-apoptotic, antioxidant, anti-proliferative, anti-inflammatory, and immunomodulatory actions. The majority of these activities are crucial for the prevention of atherogenesis. Single amino acid substitutions in proteins generated by missense non-synonymous single nucleotide polymorphism (nsSNPs) in the protein-encoding regions of genes are potent enough to cause significant medical challenges due to the alteration of protein structure and function. The current study aimed at characterizing and analyzing high-risk nsSNPs associated with the human HMOX1 gene. Preliminary screening of the total available 288 missense SNPs was performed through the lens of deleteriousness and stability prediction tools. Finally, a total of seven nsSNPs (Y58D, A131T, Y134H, F166S, F167S, R183S and M186V) were found to be most deleterious by all tools that are present at highly conserved positions. Molecular dynamics simulations (MDS) analysis explained the mutational effects on the dynamic action of the wild-type and mutant proteins. In a nutshell, R183S (rs749644285) was identified as a highly detrimental mutation that could significantly render the enzymatic activity of HMOX1. The finding of this computational analysis might help subject the experimental confirmatory analysis to characterize the role of nsSNPs in HMOX1.Communicated by Ramaswamy H. Sarma.

Chronological and Biological Aging in Amyotrophic Lateral Sclerosis and the Potential of Senolytic Therapies

Amyotrophic Lateral Sclerosis (ALS) is a group of sporadic and genetic neurodegenerative disorders that result in losses of upper and lower motor neurons. Treatment of ALS is limited, and survival is 2-5 years after disease onset. While ALS can occur in younger individuals, the risk significantly increases with advancing age. Notably, both sporadic and genetic forms of ALS share pathophysiological features overlapping hallmarks of aging including genome instability/DNA damage, mitochondrial dysfunction, inflammation, proteostasis, and cellular senescence. This review explores chronological and biological aging in the context of ALS onset and progression. Age-related muscle weakness and motor unit loss mirror aspects of ALS pathology and coincide with peak ALS incidence, suggesting a potential link between aging and disease development. Hallmarks of biological aging, including DNA damage, mitochondrial dysfunction, and cellular senescence, are implicated in both aging and ALS, offering insights into shared mechanisms underlying disease pathogenesis. Furthermore, senescence-associated secretory phenotype and senolytic treatments emerge as promising avenues for ALS intervention, with the potential to mitigate neuroinflammation and modify disease progression.

An In Silico Analysis of Genetic Variants and Structural Modeling of the Human Frataxin Protein in Friedreich's Ataxia

Friedreich's Ataxia (FRDA) stands out as the most prevalent form of hereditary ataxias, marked by progressive movement ataxia, loss of vibratory sensitivity, and skeletal deformities, severely affecting daily functioning. To date, the only medication available for treating FRDA is Omaveloxolone (Skyclarys®), recently approved by the FDA. Missense mutations within the human frataxin (FXN) gene, responsible for intracellular iron homeostasis regulation, are linked to FRDA development. These mutations induce FXN dysfunction, fostering mitochondrial iron accumulation and heightened oxidative stress, ultimately triggering neuronal cell death pathways. This study amalgamated 226 FXN genetic variants from the literature and database searches, with only 18 previously characterized. Predictive analyses revealed a notable prevalence of detrimental and destabilizing predictions for FXN mutations, predominantly impacting conserved residues crucial for protein function. Additionally, an accurate, comprehensive three-dimensional model of human FXN was constructed, serving as the basis for generating genetic variants I154F and W155R. These variants, selected for their severe clinical implications, underwent molecular dynamics (MD) simulations, unveiling flexibility and essential dynamic alterations in their N-terminal segments, encompassing FXN42, FXN56, and FXN78 domains pivotal for protein maturation. Thus, our findings indicate potential interaction profile disturbances in the FXN42, FXN56, and FXN78 domains induced by I154F and W155R mutations, aligning with the existing literature.

Mechanistic and structural insight into R2R3-MYB transcription factor in plants: molecular dynamics based binding free energy analysis

The MYB transcription factor (TF) family is essential for various plant growth and development processes, including responses to biotic and abiotic stresses. This study investigated the R2R3-MYB protein structure from five plants, including cereal crops. The R2R3-MYB protein structure was docked with the DNA structure, and the best complexes were selected for two runs of molecular dynamics (MD) simulations to investigate the key interacting residues and conformational changes in the R2R3-MYB proteins caused by DNA binding. The MM/PBSA method calculated the binding free energy for each R2R3-MYB protein-DNA complex, showing strong interaction. Hydrophobic and hydrogen bonds significantly stabilized the R2R3-MYB protein-DNA complexes. The principal component analysis showed high restrictions on the movement of protein atoms in the phase space. A similar MD simulation analysis was performed using the crystal structure of the R2R3-MYB protein-DNA complex from Arabidopsis thaliana, and the generated complexes resembled the X-ray crystal structure. This is the first detailed study on the R2R3-MYB protein-DNA complex in cereal crops, providing a cost-effective solution to identify the key interacting residues and analyze the conformational changes in the MYB domain before and after DNA binding.Communicated by Ramaswamy H. Sarma.

Molecular dynamics and simulation analysis against superoxide dismutase (SOD) target of Micrococcus luteus with secondary metabolites from Bacillus licheniformis recognized by genome mining approach

Micrococcus luteus, also known as M. luteus, is a bacterium that inhabits mucous membranes, human skin, and various environmental sources. It is commonly linked to infections, especially among individuals who have compromised immune systems. M. luteus is capable of synthesizing the enzyme superoxide dismutase (SOD) as a component of its protective response to reactive oxygen species (ROS). This enzyme serves as a promising target for drug development in various diseases. The current study utilized a subtractive genomics approach to identify potential therapeutic targets from M. luteus. Additionally, genome mining was employed to identify and characterize the biosynthetic gene clusters (BGCs) responsible for the production of secondary metabolites in Bacillus licheniformis (B. licheniformis), a bacterium known for its production of therapeutically relevant secondary metabolites. Subtractive genomics resulted in identification of important extracellular protein SOD as a drug target that plays a crucial role in shielding cells from damage caused by ROS. Genome mining resulted in identification of five potential ligands (secondary metabolites) from B. licheniformis such as, Bacillibactin (BAC), Paenibactin (PAE), Fengycin (FEN), Surfactin (SUR) and Lichenysin (LIC). Molecular docking was used to predict and analyze the binding interactions between these five ligands and target protein SOD. The resulting protein-ligand complexes were further analyzed for their motions and interactions of atoms and molecules over 250 ns using molecular dynamics (MD) simulation analysis. The analysis of MD simulations suggests, Bacillibactin as the probable candidate to arrest the activities of SOD. All the five compounds reported in this study were found to act by directly/indirectly interacting with ROS molecules, such as superoxide radicals (O2-) and hydrogen peroxide (H2O2), and transforming them into less reactive species. This antioxidant activity contributes to its protective effects against oxidative stress-induced damage in cells making them likely candidate for various applications, including in the development of antioxidant-based therapies, nutraceuticals, and functional foods.

Mechanism of the interaction of toxic SOD1 fibrils with two potent polyphenols: curcumin and quercetin

Amyotrophic lateral sclerosis (ALS) is a debilitating neurodegenerative disease commonly caused due to the aggregation of superoxide dismutase 1 (SOD1) protein. Finding inhibitors of SOD1 aggregation is of prime concern, but understanding the mechanistic action of inhibitors is equally important. Recent experiments found that two polyphenols, curcumin, and quercetin, have the ability to inhibit SOD1 aggregation. Quercetin was experimentally proven to break pre-formed fibrils into shorter segments, while curcumin did not significantly affect the pre-formed species. Here, we delve deeper into understanding the mechanism of action of quercetin and curcumin on pre-formed octameric fibrils of SOD1 (28PVKVWGSIKGL38: chains A-H) with the help of molecular dynamics (MD) simulations of a fibril docked polyphenol complex. Our results suggest that quercetin shows π-π stacking interaction with one of the key residues for toxic amyloid formation, Trp 32 of chains D, E, and F, and breaks the peptide chains G, and H from the rest of the fibril. On the other hand, curcumin binds to the hydrophobic amino acids of almost all the chains B-H and stabilizes the fibril rather than destabilizing it. Binding free energy calculations using MM/PBSA showed that curcumin binds more strongly to the SOD1 fibril due to greater van der Waals interactions compared to quercetin. These findings provide insights for the development of potential ALS treatments.

Myricetin: A Potent Anti-Amyloidogenic Polyphenol against Superoxide Dismutase 1 Aggregation

Amyotrophic lateral sclerosis (ALS) is believed to be caused by the aggregation of misfolded or mutated superoxide dismutase 1 (SOD1). As there is currently no treatment, research into aggregation inhibitors continues. Based on docking, molecular dynamics (MD) simulations, and experimental observations, we propose that myricetin, a plant flavonoid, can act as a potent anti-amyloidogenic polyphenol against SOD1 aggregation. Our MD simulation results showed that myricetin stabilizes the protein interface, destabilizes the preformed fibril, and decreases the rate of fibril elongation. Myricetin inhibits the aggregation of SOD1 in a dose-dependent manner as shown by the ThT aggregation kinetics curves. Our transmission electron microscopy, dynamic light scattering, and circular dichroism experiments indicate that fewer shorter fibrils have formed. Fluorescence spectroscopy results predict the involvement of a static quenching mechanism characterized by a strong binding between protein and myricetin. Importantly, size exclusion chromatography revealed the potential of myricetin for fibril destabilization and depolymerization. These experimental observations complement the MD results. Thus, myricetin is a potent SOD1 aggregation inhibitor that can reduce the fibril load. Using the structure of myricetin as a reference, it is possible to design more effective therapeutic inhibitors against ALS that prevent the disease and reverse its effects.

Role of charged residues of the "electrostatic loop" of hSOD1 in promotion of aggregation: Implications for the mechanism of ALS-associated mutations under amyloidogenic conditions

Protein misfolding and amyloid formation are hallmarks of numerous diseases, including amyotrophic lateral sclerosis (ALS), in which hSOD1 aggregation is involved in pathogenesis. We used two point mutations in the electrostatic loop, G138E and T137R, to analyze charge distribution under destabilizing circumstances to gain more about how ALS-linked mutations affect SOD1 protein stability or net repulsive charge. We show that protein charge is important in the ALS disease process using bioinformatics and experiments. The MD simulation findings demonstrate that the mutant protein differs significantly from WT SOD1, which is consistent with the experimental evidence. The specific activity of the wild type was 1.61 and 1.48 times higher than that of the G138E and T137R mutants, respectively. Under amyloid induction conditions, the intensity of intrinsic and ANS fluorescence in both mutants reduced. Increasing the content of β-sheet structures in mutants can be attributed to aggregation propensity, which was confirmed using CD polarimetry and FTIR spectroscopy. Our findings show that two ALS-related mutations promote the formation of amyloid-like aggregates at near physiological pH under destabilizing conditions, which were detected using spectroscopic probes such as Congo red and ThT fluorescence, and also further confirmation of amyloid-like species by TEM. Overall, our results provide evidence supporting the notion that negative charge changes combined with other destabilizing factors play an important role in increasing protein aggregation by reducing repulsive negative charges.

A double point mutation of SOD1 targeting net charge promotes aggregation under destabilizing conditions: Correlation of charge distribution and ALS-provoking mutation

A progressive neurodegenerative illness such as amyotrophic lateral sclerosis (ALS) is characterized by the misfolding and aggregation of human CuZn superoxide dismutase (hSOD1) into amyloid aggregates. Thus, designing strategies for the choice of WT-SOD1 and double mutant (G12D/G138E) with an increased net negative charge can be a good idea to elucidate the pathological mechanism of SOD1 in ALS under some destabilizing conditions. Consequently, we show evidence that protein charge, together with other destabilizing conditions, plays an important role in ALS disease. To achieve this purpose, we use methods, such as spectroscopy and transmission electron microscopy (TEM) to monitor the formation of amyloid aggregation. The specific activity of WT-SOD1 was approximately 1.72 times higher than that of the double mutant. Under amyloidogenic circumstances, structural properties such as local, secondary, oligomeric, and fibrillar structures were explored. The double mutant's far-UV CD spectra displayed a broad minimum peak in the region 213 to 218 nm, suggesting the production of β-rich amyloid fibrils. FTIR spectra of the double mutant samples at different incubation times showed a low-frequency peak around 1630-1640 cm^-1, attributed to a parallel β-sheet. Moreover, CR-binding assay and TEM analysis revealed and confirmed that mutation with an increased repulsive charge promotes the formation of fibrous aggregates. Consequently, ALS mutations with a higher repulsive charge are the apparent exceptions that validate the rule. This findings revealed that the double mutant increases protein aggregation through a novel mechanism, likely involving destabilization of structure and a change in the net negative charge.

In Silico Analyses of a Promising Drug Candidate for the Treatment of Amyotrophic Lateral Sclerosis Targeting Superoxide Dismutase I Protein

Amyotrophic lateral sclerosis (ALS) is the most prevalent motor neuron disorder in adults, which is associated with a highly disabling condition. To date, ALS remains incurable, and the only drugs approved by the FDA for its treatment confer a limited survival benefit. Recently, SOD1 binding ligand 1 (SBL-1) was shown to inhibit in vitro the oxidation of a critical residue for SOD1 aggregation, which is a central event in ALS-related neurodegeneration. In this work, we investigated the interactions between SOD1 wild-type and its most frequent variants, i.e., A4V (NP_000445.1:p.Ala5Val) and D90A (NP_000445.1:p.Asp91Val), with SBL-1 using molecular dynamics (MD) simulations. The pharmacokinetics and toxicological profile of SBL-1 were also characterized in silico. The MD results suggest that the complex SOD1-SBL-1 remains relatively stable and interacts within a close distance during the simulations. This analysis also suggests that the mechanism of action proposed by SBL-1 and its binding affinity to SOD1 may be preserved upon mutations A4V and D90A. The pharmacokinetics and toxicological assessments suggest that SBL-1 has drug-likeness characteristics with low toxicity. Our findings, therefore, suggested that SBL-1 may be a promising strategy to treat ALS based on an unprecedented mechanism, including for patients with these frequent mutations.

Genetics of amyotrophic lateral sclerosis: seeking therapeutic targets in the era of gene therapy

Amyotrophic lateral sclerosis (ALS) is an intractable disease that causes respiratory failure leading to mortality. The main locus of ALS is motor neurons. The success of antisense oligonucleotide (ASO) therapy in spinal muscular atrophy (SMA), a motor neuron disease, has triggered a paradigm shift in developing ALS therapies. The causative genes of ALS and disease-modifying genes, including those of sporadic ALS, have been identified one after another. Thus, the freedom of target choice for gene therapy has expanded by ASO strategy, leading to new avenues for therapeutic development. Tofersen for superoxide dismutase 1 (SOD1) was a pioneer in developing ASO for ALS. Improving protocols and devising early interventions for the disease are vital. In this review, we updated the knowledge of causative genes in ALS. We summarized the genetic mutations identified in familial ALS and their clinical features, focusing on SOD1, fused in sarcoma (FUS), and transacting response DNA-binding protein. The frequency of the C9ORF72 mutation is low in Japan, unlike in Europe and the United States, while SOD1 and FUS are more common, indicating that the target mutations for gene therapy vary by ethnicity. A genome-wide association study has revealed disease-modifying genes, which could be the novel target of gene therapy. The current status and prospects of gene therapy development were discussed, including ethical issues. Furthermore, we discussed the potential of axonal pathology as new therapeutic targets of ALS from the perspective of early intervention, including intra-axonal transcription factors, neuromuscular junction disconnection, dysregulated local translation, abnormal protein degradation, mitochondrial pathology, impaired axonal transport, aberrant cytoskeleton, and axon branching. We simultaneously discuss important pathological states of cell bodies: persistent stress granules, disrupted nucleocytoplasmic transport, and cryptic splicing. The development of gene therapy based on the elucidation of disease-modifying genes and early intervention in molecular pathology is expected to become an important therapeutic strategy in ALS.

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

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

Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study

More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.

In silico design of a TLR4-mediating multiepitope chimeric vaccine against amyotrophic lateral sclerosis via advanced immunoinformatics

Amyotrophic lateral sclerosis (ALS) is the most prevalent motor neuron disorder worldwide. In ALS, progressing disease can result from misfolding and aggregation of superoxide dismutase-1 (SOD1) or TAR DNA-binding protein 43 kDa (TDP43). An efficient immunotherapy for ALS should spare intact SOD1 while eliminating its dysfunctional variant. We utilized advanced immunoinformatics to suggest a potential vaccine candidate against ALS by proposing a model of dynamic TLR4 mediation and induction of a specific Th2-biased shift against mutant SOD1, TDP43, and TRAF6, a protein that specifically interacts with dysfunctional SOD1. SOD1, TDP43, and TRAF6 were retrieved in FASTA. Immune Epitopes Database and CTLpred suggested T/B-cell epitopes from disease-specific regions of selected antigens. A TLR4-mediating adjuvant, RS01, was used. Sequences were assembled via suitable linkers. Tertiary structure of the protein was calculated. Refined protein structure and physicochemical features of the 3D structure were verified in silico. Differential immune induction was assessed via C-ImmSim. GROningen MAchine for Chemical Simulation was used to assess evolution of the docked vaccine-TLR4 complex in blood. Our protein showed high structural quality and was nonallergenic and immune inducing. Also, the vaccine-TLR4 complex stability was verified by RMSD, RMSF, gyration, and visual analyses of the molecular dynamic trajectory. Contact residues in the vaccine-TLR4 complex showed favorable binding energies. Immune stimulation analyses of the proposed candidate demonstrated a sustained memory cell response and a strong adaptive immune reaction. We proposed a potential vaccine candidate against ALS and verified its physicochemical and immune inducing features. Future studies should assess this vaccine in animal studies.

Probing the polyphenolic flavonoid, morin as a highly efficacious inhibitor against amyloid(A4V) mutant SOD1 in fatal amyotrophic lateral sclerosis

Deposition of misfolded protein aggregates in key areas of human brain is the quintessential trait of various pertinent neurodegenerative disorders including amyotrophic lateral sclerosis (ALS). Genetic point mutations in Cu/Zn superoxide dismutase (SOD1) are found to be the most important contributing factor behind familial ALS. Especially, single nucleotide polymorphism (SNP) A4V is the most nocuous since it substantially decreases life expectancy of patients. Besides, the use of naturally occurring polyphenolic flavonoids is profoundly being advocated for palliating amyloidogenic behavior of proteopathic proteins. In the present analysis, through proficient computational tools, we have attempted to ascertain a pharmacodynamically promising flavonoid compound that effectively curbs the pathogenic behavior of A4V SOD1 mutant. Initial screening of flavonoids that exhibit potency against amyloids identified morin, myricetin and epigallocatechin gallate as promising leads. Further, with the help of feasible and yet adept protein-ligand interaction studies and stalwart molecular simulation analyses, we were able to observe that aforementioned flavonoids were able to considerably divert mutant A4V SOD1 from its distinct pathogenic behavior. Among which, morin showed the most curative potential against A4V SOD1. Therefore, morin holds a great therapeutic potential in contriving highly efficacious inhibitors in mitigating fatal and insuperable ALS.

In Silico Three-Dimensional (3D) Modeling of the SecY Protein of ‘Candidatus Phytoplasma Solani’ Strains Associated with Grapevine “Bois Noir” and Its Possible Relationship with Strain Virulence

Grapevine “bois noir”, related to the presence of ‘Candidatus Phytoplasma solani’ (‘Ca. P. solani’), represents a serious threat in several vine-growing areas worldwide. In surveys conducted over two years, mild and/or moderate symptoms and lower pathogen titer were mainly associated with ‘Ca. P. solani’ strains harboring a secY gene sequence variant (secY52), whereas severe symptoms and higher titer were mainly observed in grapevines infected by phytoplasma strains carrying any one of another four variants. A comparison of amino acid sequences of the protein SecY of ‘Ca. P. solani’ strains revealed the presence of conservative and semi-conservative substitutions. The deduced three-dimensional (3D) protein analysis unveiled that one semi-conservative substitution identified in the sequence variant secY52 is responsible for a structural disordered region that probably confers a flexibility for binding to distinct molecular complexes. In fact, the other analyzed variants show an organized structure and the 3D in silico prediction allowed the identification of β-sheets. Thus, differences in symptom severity and pathogen concentration observed in grapevines infected by ‘Ca. P. solani’ strains carrying distinct secY gene sequence variants suggest a possible relationship between SecY protein structure and phytoplasma strain virulence.

Intermolecular Mechanism and Dynamic Investigation of Avian Influenza H7N9 Virus' Susceptibility to E119V-Substituted Peramivir-Neuraminidase Complex

The H7N9 virus attaches itself to the human cell receptor protein containing the polysaccharide that terminates with sialic acid. The mutation of neuraminidase at residue E119 has been explored experimentally. However, there is no adequate information on the substitution with E119V in peramivir at the intermolecular level. Therefore, a good knowledge of the interatomic interactions is a prerequisite in understanding its transmission mode and subsequent effective inhibitions of the sialic acid receptor cleavage by neuraminidase. Herein, we investigated the mechanism and dynamism on the susceptibility of the E119V mutation on the peramivir-neuraminidase complex relative to the wildtype complex at the intermolecular level. This study aims to investigate the impact of the 119V substitution on the neuraminidase-peramivir complex and unveil the residues responsible for the complex conformations. We employed molecular dynamic (MD) simulations and extensive post-MD analyses in the study. These extensive computational investigations were carried out on the wildtype and the E119V mutant complex of the protein for holistic insights in unveiling the effects of this mutation on the binding affinity and the conformational terrain of peramivir-neuraminidase E119V mutation. The calculated total binding energy (ΔGbind) for the peramivir wildtype is -49.09 ± 0.13 kcal/mol, while the E119V mutant is -58.55 ± 0.15 kcal/mol. The increase in binding energy (9.46 kcal/mol) is consistent with other post-MD analyses results, confirming that E119V substitution confers a higher degree of stability on the protein complex. This study promises to proffer contributory insight and additional knowledge that would enhance future drug designs and help in the fight targeted at controlling the avian influenza H7N9 virus. Therefore, we suggest that experimentalists collaborate with computational chemists for all investigations of this topic, as we have done in our previous studies.

Molecular dynamics and protein frustration analysis of human fused in Sarcoma protein variants in Amyotrophic Lateral Sclerosis type 6: An In Silico approach

Amyotrophic lateral sclerosis (ALS) is the most frequent adult-onset motor neuron disorder. The disease is characterized by degeneration of upper and lower motor neurons, leading to death usually within five years after the onset of symptoms. While most cases are sporadic, 5%-10% of cases can be associated with familial inheritance, including ALS type 6, which is associated with mutations in the Fused in Sarcoma (FUS) gene. This work aimed to evaluate how the most frequent ALS-related mutations in FUS, R521C, R521H, and P525L affect the protein structure and function. We used prediction algorithms to analyze the effects of the non-synonymous single nucleotide polymorphisms and performed evolutionary conservation analysis, protein frustration analysis, and molecular dynamics simulations. Most of the prediction algorithms classified the three mutations as deleterious. All three mutations were predicted to reduce protein stability, especially the mutation R521C, which was also predicted to increase chaperone binding tendency. The protein frustration analysis showed an increase in frustration in the interactions involving the mutated residue 521C. Evolutionary conservation analysis showed that residues 521 and 525 of human FUS are highly conserved sites. The molecular dynamics results indicate that protein stability could be compromised in all three mutations. They also affected the exposed surface area and protein compactness. The analyzed mutations also displayed high flexibility in most residues in all variants, most notably in the interaction site with the nuclear import protein of FUS.

Bridging the Bridging Imidazolate in the Bimetallic Center of the Cu/Zn SOD1 and ALS

Metallation status of human Cu/Zn superoxide dismutase 1 (SOD1) plays a pivotal role in the pathogenesis of amyotrophic lateral sclerosis (ALS). All of the amino acids found in the bimetallic center have been associated with ALS except for two positions. H63 which forms the bridging imidazolate ion in the bimetallic center and K136 which is not directly involved in coordination but located in the bimetallic center were not reported to be mutated in any of the identified ALS cases. In this study, we investigated the structure and flexibility of five SOD1 variants by using classical molecular dynamics simulations. These variants include three substitutions on the non-ALS-linked positions; H63A, H63R, K136A and ALS-linked positions; G37R, H46R/H48D. We have generated four systems for each variant differing in metallation and presence of the intramolecular disulfide bond. Overall, a total of 24 different dimers including the wild-type were generated and simulated at two temperatures, 298 and 400 K. We have monitored backbone mobility, fluctuations and compactness of the dimer structures to assess whether the hypothetical mutations would behave similar to the ALS-linked variants. Results showed that particularly two mutants, H63R and K136A, drastically affected the dimer dynamics by increasing the fluctuations of the metal binding loops compared with the control mutations. Further, these variants resulted in demetallation of the dimers, highlighting probable ALS toxicity that could be elicited by the SOD1 variants of H63R and K136A. Overall, this study bridges two putative SOD1 positions in the metallic center and ALS, underlining the potential use of atomistic simulations for studying disease variants.