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.

In silico analyses of acetylcholinesterase (AChE) and its genetic variants in interaction with the anti-Alzheimer drug Rivastigmine.

Alzheimer's disease (AD) is the leading cause of dementia worldwide. Despite causing great social and economic impact, there is currently no cure for AD. The most effective therapy to manage AD symptoms is based on acetylcholinesterase inhibitors (AChEi), from which rivastigmine presented numerous benefits. However, mutations in AChE, which affect approximately 5% of the population, can modify protein structure and function, changing the individual response to Alzheimer's treatment. In this study, we performed computer simulations of AChE wild type and variants R34Q, P135A, V333E, and H353N, identified by one or more genome-wide association studies, to evaluate their effects on protein structure and interaction with rivastigmine. The functional effects of AChE variants were predicted using eight machine learning algorithms, while the evolutionary conservation of AChE residues was analyzed using the ConSurf server. Autodock4.2.6 was used to predict the binding modes for the hAChE-rivastigmine complex, which is still unknown. Molecular dynamics (MD) simulations were performed in triplicates for the AChE wild type and mutants using the GROMACS packages. Among the analyzed variants, P135A was classified as deleterious by all the functional prediction algorithms, in addition to occurring at highly conserved positions, which may have harmful consequences on protein function. The molecular docking results suggested that rivastigmine interacts with hAChE at the upper active-site gorge, which was further confirmed by MD simulations. Our MD findings also suggested that the complex hAChE-rivastigmine remains stable over time. The essential dynamics revealed flexibility alterations at the active-site gorge upon mutations P135A, V333E, and H353N, which may lead to strong and nonintuitive consequences to hAChE binding. Nonetheless, similar binding affinities were registered in the MMPBSA analysis for the hAChE wild type and variants when complexed to rivastigmine. Finally, our findings indicated that the rivastigmine binding to hAChE is an energetically favorable process mainly driven by negatively charged amino acids.

Trehalose synthesis inhibitor: A molecular in silico drug design.

Infectious diseases are serious public health problems, affecting a large portion of the world's population. A molecule that plays a key role in pathogenic organisms is trehalose and recently has been an interest in the metabolism of this molecule for drug development. The trehalose-6-phosphate synthase (TPS1) is an enzyme responsible for the biosynthesis of trehalose-6-phosphate (T6P) in the TPS1/TPS2 pathway, which results in the formation of trehalose. Studies carried out by our group demonstrated the inhibitory capacity of T6P in the TPS1 enzyme from Saccharomyces cerevisiae, preventing the synthesis of trehalose. By in silico techniques, we compiled sequences and experimentally determined structures of TPS1. Sequence alignments and molecular modeling were performed. The generated structures were submitted in validation of algorithms, aligned structurally and analyzed evolutionarily. Molecular docking methodology was applied to analyze the interaction between T6P and TPS1 and ADMET properties of T6P were analyzed. The results demonstrated the models created presented sequence and structural similarities with experimentally determined structures. With the molecular docking, a cavity in the protein surface was identified and the molecule T6P was interacting with the residues TYR-40, ALA-41, MET-42, and PHE-372, indicating the possible uncompetitive inhibition mechanism provided by this ligand, which can be useful in directing the molecular design of inhibitors. In ADMET analyses, T6P had acceptable risk values compared with other compounds from World Drug Index. Therefore, these results may present a promising strategy to explore to develop a broad-spectrum antibiotic of this specific target with selectivity, potency, and reduced side effects, leading to a new way to treat infectious diseases like tuberculosis and candidiasis.

SOD1 in amyotrophic lateral sclerosis development - in silico analysis and molecular dynamics of A4F and A4V variants.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease that is characterized by the selective loss of motor neurons. Approximately 5% to 10% of patients with ALS have a family history of the disease, and approximately 20% of familial amyotrophic lateral sclerosis (fALS) cases are associated with mutations in Cu/Zn superoxide dismutase (SOD1). In this study, we evaluated the structural and functional effects of human A4F and A4V SOD1 protein mutations. We performed an in silico analysis using prediction algorithms of nonsynonymous single-nucleotide polymorphisms (nsSNPs) associated with the fALS development. Our structural conservation results show that the mutations analyzed (A4V and A4F) were in a highly conserved region. Molecular dynamics simulations using the Linux GROMACS package revealed how these mutations affect protein structure, protein stability, and aggregation. These results suggest that there might be an effect on the SOD1 function. Understanding the molecular basis of disease provides new insights useful for rational drug design and advancing our understanding of the ALS development.

In Vivo Characterization of I91T Sod2 Polymorphism of Saccharomyces cerevisiae.

The mitochondrial antioxidant enzyme Mn-Superoxide Dismutase (Sod2) is essential for mammalian survival. I82T mutation in human Sod2 has been linked to a wide variety of diseases, including Alzheimer's and Parkinson's diseases as well as some types of cancers. Yeast wild-type (WT) Sod2 and the mutant Sod2 I91T, which corresponds to the human mutant Sod2 I82T, were cloned in sod2Δ strain. Residue I82 is conserved among a variety of species, showing that it has a biological importance. To assess the functionality of Sod2 I91T under oxidative stress, yeast cells were shifted from glucose (fermentative metabolism) to glycerol growth medium (respiratory metabolism). Overexpression of both Sod2 WT and Sod2 I91T increased Sod activity, but in long-term, the mutation brought impairment to Sod function. Aconitase, a sensor of superoxide radical production in vivo, had its activity preserved by overexpressions of both Sod2, in lesser extent in sod2ΔSod2I91T. In respiratory metabolism, sod2ΔSod2WT and sod2ΔSod2I91T showed high viability; although, sod2ΔSod2I91T showed high percentage of cells with mitochondrial function compromised. Moreover, the fitness analysis of mixed cultures showed that sod2ΔSod2I91T was less robust than WT cells. Although overexpression of Sod2 containing I91T mutation allows higher cell viability, longevity of cells is hampered, showing that in long-term this mutation is not neutral. J. Cell. Biochem. 118: 1078-1086, 2017. © 2016 Wiley Periodicals, Inc.

Revisiting yeast trehalose metabolism.

Establishing the function of trehalose in yeast cells has led us, over the years, through a long path-from simple energy storage carbohydrate, then a stabilizer and protector of membranes and proteins, through a safety valve against damage caused by oxygen radicals, up to regulator of the glycolytic path. In addition, trehalose biosynthesis has been proposed as a target for novel drugs against several pathogens. Since this pathway is entirely absent in mammalian cells and makes use of highly specific enzymes, trehalose metabolism might be an interesting target for the development of novel therapies. In this review, we want to address some recent points investigated about trehalose metabolism in Saccharomyces cerevisiae, focusing mainly on the mechanism by which this simple disaccharide protects against stress and on the enzymes involved in its synthesis and breakdown. We believe that these concepts are of great importance for medical and biotechnological applications.

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.

Molecular mimicry between Zika virus and central nervous system inflammatory demyelinating disorders: the role of NS5 Zika virus epitope and PLP autoantigens.

BACKGROUND: Evidence indicates a strong link between Zika virus (ZikV) and neurological complications. Acute myelitis, optic neuritis, polyneuropathy, and encephalomyelitis that mimic inflammatory idiopathic demyelination disorders (IIDD) after ZikV infection have been reported in Brazil. OBJECTIVE: The present study aims to investigate the possible occurrence of molecular mimicry between ZikV antigens and Multiple Sclerosis (MS) autoantigens, the most frequent IIDD of the central nervous system (CNS). METHODS: A retrospective cohort study with 305 patients admitted due to suspected arbovirus infection in Rio de Janeiro was performed, all subjects were submitted to neurological examination, and a biological sample was collected for serologic and molecular diagnostic. Bioinformatics tools were used to analyze the peptides shared between ZikV antigens and MS autoantigens. RESULTS: Of 305 patients, twenty-six were positive for ZikV and 4 presented IDD patterns found in MS cases. Sequence homology comparisons by bioinformatics approach between NS5 ZikV and PLP MS protein revealed a homology of 5/6 consecutive amino acids (CSSVPV/CSAVPV) with 83% identity, deducing a molecular mimicry. Analysis of the 3D structures revealed a similar conformation with alpha helix presentation. CONCLUSIONS: Molecular mimicry between NS5 Zika virus antigen and PLP MS autoantigens emerge as a possible mechanism for IDD spectrum in genetically susceptible individuals.

ANTECEDENTES: Evidências indicam uma forte ligação entre o vírus Zika (ZikV) e complicações neurológicas. Mielite aguda, neurite óptica, polineuropatia e encefalomielite que mimetizam distúrbios inflamatórios de desmielinização idiopáticos (DDII) após infecção por ZikV têm sido relatadas no Brasil. OBEJTIVO: O presente estudo tem como objetivo investigar a possível ocorrência de mimetismo molecular entre antígenos do ZikV e autoantígenos da Esclerose Múltipla (EM), a DDII mais frequente do sistema nervoso central (SNC). MéTODOS: Foi realizado um estudo de coorte retrospectivo com 305 pacientes internados por suspeita de infecção por arbovírus no Rio de Janeiro, todos os indivíduos foram submetidos a exame neurológico e coleta de amostra biológica para diagnóstico sorológico e molecular. Ferramentas de bioinformática foram usadas para analisar os peptídeos compartilhados entre antígenos do ZikV e autoantígenos da EM. RESULTADOS: Dos 305 pacientes, vinte e seis foram positivos para ZikV e 4 apresentaram padrão IDD encontrado em casos de EM. As comparações de homologia de sequência por abordagem de bioinformática entre a proteína NS5 ZikV e PLP EM revelaram uma homologia de 5/6 aminoácidos consecutivos (CSSVPV/CSAVPV) com 83% de identidade, deduzindo um mimetismo molecular. A análise das estruturas 3D revelou uma conformação semelhante com apresentação em alfa-hélice. CONCLUSõES: O mimetismo molecular entre o antígeno NS5 do vírus Zika e o autoantígeno PLP da EM surge como um possível mecanismo para o espectro IDD em indivíduos geneticamente suscetíveis.

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.

Quinolines-Based SARS-CoV-2 3CLpro and RdRp Inhibitors and Spike-RBD-ACE2 Inhibitor for Drug-Repurposing Against COVID-19: An in silico Analysis.

The novel coronavirus SARS-CoV-2 disease "COVID-19" emerged in China and rapidly spread to other countries; due to its rapid worldwide spread, the WHO has declared this as a global emergency. As there is no specific treatment prescribed to treat COVID-19, the seeking of suitable therapeutics among existing drugs seems valuable. The structure availability of coronavirus macromolecules has encouraged the finding of conceivable anti-SARS-CoV-2 therapeutics through in silico analysis. The results reveal that quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) and saquinavir strongly interact with the active site (Cys-His catalytic dyad), thereby are predicted to hinder the activity of SARS-CoV-2 3CLpro. Out of 113 quinoline-drugs, elvitegravir and oxolinic acid are able to interact with the NTP entry-channel and thus interfere with the RNA-directed 5'-3' polymerase activity of SARS-CoV-2 RdRp. The bioactivity-prediction results also validate the outcome of the docking study. Moreover, as SARS-CoV-2 Spike-glycoprotein uses human ACE2-receptor for viral entry, targeting the Spike-RBD-ACE2 has been viewed as a promising strategy to control the infection. The result shows rilapladib is the only quinoline that can interrupt the Spike-RBD-ACE2 complex. In conclusion, owing to their ability to target functional macromolecules of SARS-CoV-2, along with positive ADMET properties, quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI), saquinavir, elvitegravir, oxolinic acid, and rilapladib are suggested for the treatment of COVID-19.

In silico analysis of the tryptophan hydroxylase 2 (TPH2) protein variants related to psychiatric disorders.

The tryptophan hydroxylase 2 (TPH2) enzyme catalyzes the first step of serotonin biosynthesis. Serotonin is known for its role in several homeostatic systems related to sleep, mood, and food intake. As the reaction catalyzed by TPH2 is the rate-limiting step of serotonin biosynthesis, mutations in TPH2 have been associated with several psychiatric disorders (PD). This work undertakes an in silico analysis of the effects of genetic mutations in the human TPH2 protein. Ten algorithms were used to predict the functional and stability effects of the TPH2 mutations. ConSurf was used to estimate the evolutionary conservation of TPH2 amino acids. GROMACS was used to perform molecular dynamics (MD) simulations of TPH2 WT and P260S, R303W, and R441H, which had already been associated with the development of PD. Forty-six TPH2 variants were compiled from the literature. Among the analyzed variants, those occurring at the catalytic domain were shown to be more damaging to protein structure and function. The ConSurf analysis indicated that the mutations affecting the catalytic domain were also more conserved throughout evolution. The variants S364K and S383F were predicted to be deleterious by all the functional algorithms used and occurred at conserved positions, suggesting that they might be deleterious. The MD analyses indicate that the mutations P206S, R303W, and R441H affect TPH2 flexibility and essential mobility at the catalytic and oligomerization domains. The variants P206S, R303W, and R441H also exhibited alterations in dimer binding affinity and stability throughout the simulations. Thus, these mutations may impair TPH2 functional interactions and, consequently, its function, leading to the development of PD. Furthermore, we developed a database, SNPMOL (http://www.snpmol.org/), containing the results presented in this paper. Understanding the effects of TPH2 mutations on protein structure and function may lead to improvements in existing treatments for PD and facilitate the design of further experiments.

In silico analysis of PFN1 related to amyotrophic lateral sclerosis.

Profilin 1 (PFN1) protein plays key roles in neuronal growth and differentiation, membrane trafficking, and regulation of the actin cytoskeleton. Four natural variants of PFN1 were described as related to ALS, the most common adult-onset motor neuron disorder. However, the pathological mechanism of PFN1 in ALS is not yet completely understood. The goal of this work is to thoroughly analyze the effects of the ALS-related mutations on PFN1 structure and function using computational simulations. Here, PhD-SNP, PMUT, PolyPhen-2, SIFT, SNAP, SNPS&GO, SAAP, nsSNPAnalyzer, SNPeffect4.0 and I-Mutant2.0 were used to predict the functional and stability effects of PFN1 mutations. ConSurf was used for the evolutionary conservation analysis, and GROMACS was used to perform the MD simulations. The mutations C71G, M114T, and G118V, but not E117G, were predicted as deleterious by most of the functional prediction algorithms that were used. The stability prediction indicated that the ALS-related mutations could destabilize PFN1. The ConSurf analysis indicated that the mutation C71G, M114T, E117G, and G118V occur in highly conserved positions. The MD results indicated that the studied mutations could affect the PFN1 flexibility at the actin and PLP-binding domains, and consequently, their intermolecular interactions. It may be therefore related to the functional impairment of PFN1 upon C71G, M114T, E117G and G118V mutations, and their involvement in ALS development. We also developed a database, SNPMOL (http://www.snpmol.org/), containing the results presented on this paper for biologists and clinicians to exploit PFN1 and its natural variants.

In silico analysis of the V66M variant of human BDNF in psychiatric disorders: An approach to precision medicine.

Brain-derived neurotrophic factor (BDNF) plays an important role in neurogenesis and synapse formation. The V66M is the most prevalent BDNF mutation in humans and impairs the function and distribution of BDNF. This mutation is related to several psychiatric disorders. The pro-region of BDNF, particularly position 66 and its adjacent residues, are determinant for the intracellular sorting and activity-dependent secretion of BDNF. However, it has not yet been fully elucidated. The present study aims to analyze the effects of the V66M mutation on BDNF structure and function. Here, we applied nine algorithms, including SIFT and PolyPhen-2, for functional and stability prediction of the V66M mutation. The complete theoretical model of BNDF was generated by Rosetta and validated by PROCHECK, RAMPAGE, ProSa, QMEAN and Verify-3D algorithms. Structural alignment was performed using TM-align. Phylogenetic analysis was performed using the ConSurf server. Molecular dynamics (MD) simulations were performed and analyzed using the GROMACS 2018.2 package. The V66M mutation was predicted as deleterious by PolyPhen-2 and SIFT in addition to being predicted as destabilizing by I-Mutant. According to SNPeffect, the V66M mutation does not affect protein aggregation, amyloid propensity, and chaperone binding. The complete theoretical structure of BDNF proved to be a reliable model. Phylogenetic analysis indicated that the V66M mutation of BDNF occurs at a non-conserved position of the protein. MD analyses indicated that the V66M mutation does not affect the BDNF flexibility and surface-to-volume ratio, but affects the BDNF essential motions, hydrogen-bonding and secondary structure particularly at its pre and pro-domain, which are crucial for its activity and distribution. Thus, considering that these parameters are determinant for protein interactions and, consequently, protein function; the alterations observed throughout the MD analyses may be related to the functional impairment of BDNF upon V66M mutation, as well as its involvement in psychiatric disorders.

Protective effect of neglected plant Diplocyclos palmatus on quorum sensing mediated infection of Serratia marcescens and UV-A induced photoaging in model Caenorhabditis elegans.

Plants are considered to be a leading source for possible human therapeutic agents. This holistic study has investigated the anti-quorum sensing (anti-QS), anti-infection, antioxidant and anti-photoaging properties of neglected plant Diplocyclos palmatus. The results showed that D. palmatus methanolic leaf extract (DPME) effectively inhibited the quorum sensing (QS) regulated virulence factor production as well as biofilm formation in Serratia marcescens. The transcriptomic analysis revealed that DPME significantly downed the expression of QS-regulated genes such as fimA, fimC, flhC, bsmB, pigP and shlA in S. marcescens, which supports the outcome of in vitro bioassays. Further, the docking study revealed that the presence of active compounds, namely tocopherols and phytol, DPME exhibited its anti-QS activity against S. marcescens. In addition, DPME treatment extended the lifespan of S. marcescens infected C. elegans by the action of dropping the internal accumulation. Further, qPCR analysis clearly revealed that DPME treatment significantly up-regulated the expression of the lifespan-related gene (daf-16) and immune-related genes (clec-60, clec-87, lys-7 and bec-1) in S. marcescens infected C.elegans. On the other hand, DPME extensively reduced the UV-A induced ROS stress, thereby, extended the lifespan in UV-A photoaged C. elegans. Further, the qPCR analysis also confirmed the up-regulation of daf-16, clec-60, clec-87 and col-19 genes which advocated the improvement of the lifespan, healthspan and collagen production in UV-A photoaged C. elegans. Further bioassays evidenced that that the lifespan extension of photoaged C. elegans was accomplished by the actions of antioxidants such as tocopherols and phytol in DPME.

Amyotrophic Lateral Sclerosis Type 20 - In Silico Analysis and Molecular Dynamics Simulation of hnRNPA1.

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease that affects the upper and lower motor neurons. 5-10% of cases are genetically inherited, including ALS type 20, which is caused by mutations in the hnRNPA1 gene. The goals of this work are to analyze the effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on hnRNPA1 protein function, to model the complete tridimensional structure of the protein using computational methods and to assess structural and functional differences between the wild type and its variants through Molecular Dynamics simulations. nsSNP, PhD-SNP, Polyphen2, SIFT, SNAP, SNPs&GO, SNPeffect and PROVEAN were used to predict the functional effects of nsSNPs. Ab initio modeling of hnRNPA1 was made using Rosetta and refined using KoBaMIN. The structure was validated by PROCHECK, Rampage, ERRAT, Verify3D, ProSA and Qmean. TM-align was used for the structural alignment. FoldIndex, DICHOT, ELM, D2P2, Disopred and DisEMBL were used to predict disordered regions within the protein. Amino acid conservation analysis was assessed by Consurf, and the molecular dynamics simulations were performed using GROMACS. Mutations D314V and D314N were predicted to increase amyloid propensity, and predicted as deleterious by at least three algorithms, while mutation N73S was predicted as neutral by all the algorithms. D314N and D314V occur in a highly conserved amino acid. The Molecular Dynamics results indicate that all mutations increase protein stability when compared to the wild type. Mutants D314N and N319S showed higher overall dimensions and accessible surface when compared to the wild type. The flexibility level of the C-terminal residues of hnRNPA1 is affected by all mutations, which may affect protein function, especially regarding the protein ability to interact with other proteins.

Structural modeling and in silico analysis of human superoxide dismutase 2.

Aging in the world population has increased every year. Superoxide dismutase 2 (Mn-SOD or SOD2) protects against oxidative stress, a main factor influencing cellular longevity. Polymorphisms in SOD2 have been associated with the development of neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, as well as psychiatric disorders, such as schizophrenia, depression and bipolar disorder. In this study, all of the described natural variants (S10I, A16V, E66V, G76R, I82T and R156W) of SOD2 were subjected to in silico analysis using eight different algorithms: SNPeffect, PolyPhen-2, PhD-SNP, PMUT, SIFT, SNAP, SNPs&GO and nsSNPAnalyzer. This analysis revealed disparate results for a few of the algorithms. The results showed that, from at least one algorithm, each amino acid substitution appears to harmfully affect the protein. Structural theoretical models were created for variants through comparative modelling performed using the MHOLline server (which includes MODELLER and PROCHECK) and ab initio modelling, using the I-Tasser server. The predicted models were evaluated using TM-align, and the results show that the models were constructed with high accuracy. The RMSD values of the modelled mutants indicated likely pathogenicity for all missense mutations. Structural phylogenetic analysis using ConSurf revealed that human SOD2 is highly conserved. As a result, a human-curated database was generated that enables biologists and clinicians to explore SOD2 nsSNPs, including predictions of their effects and visualisation of the alignment of both the wild-type and mutant structures. The database is freely available at http://bioinfogroup.com/database and will be regularly updated.

Structural and functional analysis of human SOD1 in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with familial inheritance (fALS) in 5% to 10% of cases; 25% of those are caused by mutations in the superoxide dismutase 1 (SOD1) protein. More than 100 mutations in the SOD1 gene have been associated with fALS, altering the geometry of the active site, protein folding and the interaction between monomers. We performed a functional analysis of non-synonymous single nucleotide polymorphisms (nsSNPs) in 124 fALS SOD1 mutants. Eleven different algorithms were used to estimate the functional impact of the replacement of one amino acid on protein structure: SNPs&GO, PolyPhen-2, SNAP, PMUT, Sift, PhD-SNP, nsSNPAnalyzer, TANGO, WALTZ, LIMBO and FoldX. For the structural analysis, theoretical models of 124 SNPs of SOD1 were created by comparative modeling using the MHOLline workflow, which includes Modeller and Procheck. Models were aligned with the native protein by the TM-align algorithm. A human-curated database was developed using the server side include in Java, JMOL. The results of this functional analysis indicate that the majority of the 124 natural mutants are harmful to the protein structure and thus corroborate the correlation between the reported mutations and fALS. In the structural analysis, all models showed conformational changes when compared to wild-type SOD1, and the degree of structural alignment varied between them. The SOD1 database converge structural and functional analyses of SOD1; it is a vast resource for the molecular analysis of amyotrophic lateral sclerosis, which allows the user to expand his knowledge on the molecular basis of the disease. The SOD1 database is available at http://bioinfogroup.com/database.