Computational Investigation on Electrostatic Loop Mutants Instigating Destabilization and Aggregation on Human SOD1 Protein Causing Amyotrophic Lateral Sclerosis

Theoretical insights into the binding of mono/di-ethyl phthalates to superoxide dismutase and associated structural changes impairing antioxidant activity: A coupled molecular docking and dynamics simulation approach

Phthalic acid esters (PAEs) are plasticizers known to increase oxidative stress by impairing antioxidant defenses, including superoxide dismutase (SOD) activity. Since oxidative stress plays a critical role in disease development, the disruption of SOD function by PAEs presents a significant concern. However, the precise molecular mechanisms underlying the regulation of SOD activity remain unclear. This study investigated how diethyl phthalate (DEP) and its major metabolite, monoethyl phthalate (MEP), affected SOD activity using molecular docking and dynamics simulations. The results revealed that both DEP and MEP bound to SOD through weak hydrophilic interactions and hydrogen bonds with residues Lys9, Thr17, Asn51, Thr52, and Arg141 in the bottom of the enzyme's two subunit cavities. These interactions triggered structural changes, particularly in the electrostatic loop and catalytic channels, destabilizing SOD. DEP and MEP increased the enzyme's radius of gyration and solvent-accessible surface area while disrupting intra-protein interactions. MEP showed a stronger inhibitory effect, significantly altering SOD's conformation. This change correlated with reduced catalytic activity (R2 > 0.9). Consequently, the inhibition of the enzyme is primarily due to the disruption of Arg141's conformation and function, which weakens SOD's antioxidant defense and potentially contributes to diseases related to oxidative damage. These results underscore the health risks posed by PAEs, especially following metabolic transformation, and highlight the importance of addressing their oxidative impact.

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.

Molecular docking and molecular dynamics of hypoxia-inducible factor (HIF-1alpha): towards potential inhibitors

HIF-1α is a primary regulator in the adaptation of cancer cells to hypoxia. The aim was to find out new inhibitors of the HIF-1α. A molecular dynamic (MD) simulation performed on HIF-1α showed stable dynamic features. Virtual screening of 217 anticancer drugs was performed along with a positive control (2-Methoxyestradiolm, 2-ME2) on an optimized HIF-1α and dynamically simulated structure. Docking results produced two compounds namely pycnidione and nilotinib of high binding affinity -9.34 kcal/mol and -9.04 kcal/mol respectively, whereas 2-ME2 displayed a relatively lower affinity (-6.68 kcal/mol). For the three complexes, MD of 200 ns simulation was run. Data analysis showed that the three medications behaved similarly in the MD simulation. Nilotinib had a lower RMSD and higher SASA than the other complexes. In addition, the Nilotinib-HIF-1α combination had a lower RMSF value, a flatter Rg, and a number of hydrogen bonds similar to other complexes. MM-GBSA analysis revealed that nilotinib, pycnidione and 2-ME2 compounds had free binding energy of -23.77 ± 5.29, -21.85 ± 4.24 and -7.53 ± 6.62 kcal/mol respectively. Nilotinib and pycnidione bind competitively to HIF-1α, with nilotinib showing consistent molecular-dynamic properties. They relatively pass the blood-brain barrier, non-carcinogenic, and have IV-category acute oral toxicity. They have low CYP inhibitory characteristics. Further investigations are therefore warranted to elucidate their implications in hypoxia pathways, cell proliferation, apoptosis, survival, and metastatic potential.

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.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Characterization of the p.L145F and p.S135N Mutations in SOD1: Impact on the Metabolism of Fibroblasts Derived from Amyotrophic Lateral Sclerosis Patients

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by the loss of the upper and lower motor neurons (MNs). About 10% of patients have a family history (familial, fALS); however, most patients seem to develop the sporadic form of the disease (sALS). SOD1 (Cu/Zn superoxide dismutase-1) is the first studied gene among the ones related to ALS. Mutant SOD1 can adopt multiple misfolded conformation, lose the correct coordination of metal binding, decrease structural stability, and form aggregates. For all these reasons, it is complicated to characterize the conformational alterations of the ALS-associated mutant SOD1, and how they relate to toxicity. In this work, we performed a multilayered study on fibroblasts derived from two ALS patients, namely SOD1^L145F and SOD1^S135N, carrying the p.L145F and the p.S135N missense variants, respectively. The patients showed diverse symptoms and disease progression in accordance with our bioinformatic analysis, which predicted the different effects of the two mutations in terms of protein structure. Interestingly, both mutations had an effect on the fibroblast energy metabolisms. However, while the SOD1^L145F fibroblasts still relied more on oxidative phosphorylation, the SOD1^S135N fibroblasts showed a metabolic shift toward glycolysis. Our study suggests that SOD1 mutations might lead to alterations in the energy metabolism.

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

Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.

Novel SOD1 monoclonal antibodies against the electrostatic loop preferentially detect misfolded SOD1 aggregates

Amyotrophic lateral sclerosis (ALS) is a progressive neurological disease that leads to motor neuron degeneration and paralysis. Superoxide dismutase (SOD1) mutations are the second most common cause of familial ALS and are responsible for up to 20 % of familial ALS cases. In ALS patients, SOD1 can form toxic misfolded aggregates that deposit in the brain and spinal cord. To better detect SOD1 aggregates and expand the repertoire of conformational SOD1 antibodies, SOD1 monoclonal antibodies were generated by immunizing SOD1 knockout mice with an SOD1 fragment consisting of amino acids 129-146, which make up part of the electrostatic loop. A series of hybridomas secreting antibodies were screened and five different SOD1 monoclonal antibodies (2C10, 2F8, 4B11, 5H5, and 5A10) were found to preferentially detect denatured or aggregated SOD1 by enzyme-linked immunosorbent assay (ELISA), filter trap assay, and immunohistochemical analysis in SOD1 mouse models. The staining with these antibodies was compared to Campbell-Switzer argyrophilic reactivity of pathological inclusions. These new conformational selective SOD1 antibodies will be useful for clinical diagnosis of SOD1 ALS and potentially for passive immunotherapy.

Exploring the effect of nsSNPs in human YPEL3 gene in cellular senescence

YPEL3 that induces cellular senescence in both normal and tumour cells of humans may show altered expression under the influence of incidental mutations. In this study, we proposed the first structure of Native YPEL3 protein and its five possible deleterious mutants—V40M, C61Y, G98R, G108S, and A131T and predicted their deleterious effects to alter stability, flexibility and conformational changes in the protein. The MD simulation (RMSD, RMSF, Rg, h-bond and SASA) analysis revealed that the variants V40M, G98R and G108S increased the flexibility in protein, and variant V40M imparted more compactness to the protein.. In general, variants attributed changes in the native conformation and structure of the YPEL3 protein which might affect the native function of cellular senescence. The study provides opportunities for health professionals and practitioners in formulating précised medicines to effectively cure various cancers. We propose in-vitro or in-vivo studies should consider these reported nsSNPs while examining any malfunction in the YPEL3 protein.

Evolutionary Analyses of Sequence and Structure Space Unravel the Structural Facets of SOD1

Superoxide dismutase (SOD) is the primary enzyme of the cellular antioxidant defense cascade. Misfolding, concomitant oligomerization, and higher order aggregation of human cytosolic SOD are linked to amyotrophic lateral sclerosis (ALS). Although, with two metal ion cofactors SOD1 is extremely robust, the de-metallated apo form is intrinsically disordered. Since the rise of oxygen-based metabolism and antioxidant defense systems are evolutionary coupled, SOD is an interesting protein with a deep evolutionary history. We deployed statistical analysis of sequence space to decode evolutionarily co-varying residues in this protein. These were validated by applying graph theoretical modelling to understand the impact of the presence of metal ion co-factors in dictating the disordered (apo) to hidden disordered (wild-type SOD1) transition. Contact maps were generated for different variants, and the selected significant residues were mapped on separate structure networks. Sequence space analysis coupled with structure networks helped us to map the evolutionarily coupled co-varying patches in the SOD1 and its metal-depleted variants. In addition, using structure network analysis, the residues with a major impact on the internal dynamics of the protein structure were investigated. Our results reveal that the bulk of these evolutionarily co-varying residues are localized in the loop regions and positioned differentially depending upon the metal residence and concomitant steric restrictions of the loops.

Structure and Function of an Inflammatory Cytokine, Interleukin-2, Analyzed Using the Bioinformatic Approach

The inflammatory cytokine, interleukin-2 (IL-2), is an important regulator of cellular functions. This relatively less studied member of the interleukin protein family is responsible for multiple immuno-modulatory and immuno-stimulatory tasks, like T cell activation, triggering of natural killer cells, inflammation, as well as proliferation and progression of autoimmune diseases and cancers. In this communication we report the temporally variant structural aspects of the IL-2 ligand and its receptor interfaces, based on the available crystal structures. The intended goal of this effort is to generate simulated results that could potentially aid the designs of novel structure based therapeutics.