In silico analysis of PFN1 related to amyotrophic lateral sclerosis

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 analysis of alpha-synuclein protein variants and posttranslational modifications related to Parkinson's disease

Parkinson's disease (PD) is among the most prevalent neurodegenerative disorders, affecting over 10 million people worldwide. The protein encoded by the SNCA gene, alpha-synuclein (ASYN), is the major component of Lewy body (LB) aggregates, a histopathological hallmark of PD. Mutations and posttranslational modifications (PTMs) in ASYN are known to influence protein aggregation and LB formation, possibly playing a crucial role in PD pathogenesis. In this work, we applied computational methods to characterize the effects of missense mutations and PTMs on the structure and function of ASYN. Missense mutations in ASYN were compiled from the literature/databases and underwent a comprehensive predictive analysis. Phosphorylation and SUMOylation sites of ASYN were retrieved from databases and predicted by algorithms. ConSurf was used to estimate the evolutionary conservation of ASYN amino acids. Molecular dynamics (MD) simulations of ASYN wild-type and variants A30G, A30P, A53T, and G51D were performed using the GROMACS package. Seventy-seven missense mutations in ASYN were compiled. Although most mutations were not predicted to affect ASYN stability, aggregation propensity, amyloid formation, and chaperone binding, the analyzed mutations received relatively high rates of deleterious predictions and predominantly occurred at evolutionarily conserved sites within the protein. Moreover, our predictive analyses suggested that the following mutations may be possibly harmful to ASYN and, consequently, potential targets for future investigation: K6N, T22I, K34E, G36R, G36S, V37F, L38P, G41D, and K102E. The MD analyses pointed to remarkable flexibility and essential dynamics alterations at nearly all domains of the studied variants, which could lead to impaired contact between NAC and the C-terminal domain triggering protein aggregation. These alterations may have functional implications for ASYN and provide important insight into the molecular mechanism of PD, supporting the design of future biomedical research and improvements in existing therapies for the disease.

Comparison between PFN1 and SOD1 mutations in amyotrophic lateral sclerosis

BACKGROUND

The objective of this study was to characterize the prototypical phenotype of patients with amyotrophic lateral sclerosis (ALS) associated with PFN1 mutations in profilin 1 (PFN1) and to determine clinical indications to test for mutations in this gene.

MATERIAL AND METHODS

The phenotype of three relatives carrying the M114V PFN1 mutation are detailed here and are compared with those of patients with ALS linked to PFN1 previously reported in the literature.

RESULTS

In this pedigree and in the literature, the main clinical findings which best describe familial ALS linked to PFN1 might be the following characteristics: pedigrees over five cases, age of onset around 50 years, site of onset systematically lower limbs and the absence of cognitive impairment.

CONCLUSION

First, the infrequent incidence of patients with ALS linked to PFN1 mutation supports the pursuit of a precise characterization of the phenotype linked to PFN1 mutations. Then, the numerous similarities between the phenotype amongst patients linked to SOD1 and PFN1 mutations and between histological features amongst both mice models prompts a review of the current ALS classifications, taking into consideration both phenotype and genotype.

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.

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.

PSnpBind: a database of mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow

A key concept in drug design is how natural variants, especially the ones occurring in the binding site of drug targets, affect the inter-individual drug response and efficacy by altering binding affinity. These effects have been studied on very limited and small datasets while, ideally, a large dataset of binding affinity changes due to binding site single-nucleotide polymorphisms (SNPs) is needed for evaluation. However, to the best of our knowledge, such a dataset does not exist. Thus, a reference dataset of ligands binding affinities to proteins with all their reported binding sites' variants was constructed using a molecular docking approach. Having a large database of protein-ligand complexes covering a wide range of binding pocket mutations and a large small molecules' landscape is of great importance for several types of studies. For example, developing machine learning algorithms to predict protein-ligand affinity or a SNP effect on it requires an extensive amount of data. In this work, we present PSnpBind: A large database of 0.6 million mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow. It provides a web interface to explore and visualize the protein-ligand complexes and a REST API to programmatically access the different aspects of the database contents. PSnpBind is open source and freely available at https://psnpbind.org .

Biochemical characterization of actin assembly mechanisms with ALS-associated profilin variants

Eight separate mutations in the actin-binding protein profilin-1 have been identified as a rare cause of amyotrophic lateral sclerosis (ALS). Profilin is essential for many neuronal cell processes through its regulation of lipids, nuclear signals, and cytoskeletal dynamics, including actin filament assembly. Direct interactions between profilin and actin monomers inhibit actin filament polymerization. In contrast, profilin can also stimulate polymerization by simultaneously binding actin monomers and proline-rich tracts found in other proteins. Whether the ALS-associated mutations in profilin compromise these actin assembly functions is unclear. We performed a quantitative biochemical comparison of the direct and formin mediated impact for the eight ALS-associated profilin variants on actin assembly using classic protein-binding and single-filament microscopy assays. We determined that the binding constant of each profilin for actin monomers generally correlates with the actin nucleation strength associated with each ALS-related profilin. In the presence of formin, the A20T, R136W, Q139L, and C71G variants failed to activate the elongation phase of actin assembly. This diverse range of formin-activities is not fully explained through profilin-poly-L-proline (PLP) interactions, as all ALS-associated variants bind a formin-derived PLP peptide with similar affinities. However, chemical denaturation experiments suggest that the folding stability of these profilins impact some of these effects on actin assembly. Thus, changes in profilin protein stability and alterations in actin filament polymerization may both contribute to the profilin-mediated actin disruptions in ALS.

In silico studies reveal structural deviations of mutant profilin-1 and interaction with riluzole and edaravone in amyotrophic lateral sclerosis

This study aimed to investigate four of the eight PFN-1 mutations that are located near the actin-binding domain and determine the structural changes due to each mutant and unravel how these mutations alter protein structural behavior. Swapaa's command in UCSF chimera for generating mutations, FTMAP were employed and the data was analyzed by RMSD, RMSF graphs, Rg, hydrogen bonding analysis, and RRdisMaps utilizing Autodock4 and GROMACS. The functional changes and virtual screening, structural dynamics, and chemical bonding behavior changes, molecular docking simulation with two current FDA-approved drugs for ALS were investigated. The highest reduction and increase in Rg were found to exist in the G117V and M113T mutants, respectively. The RMSF data consistently shows changes nearby to this site. The in silico data described indicate that each of the mutations is capable of altering the structure of PFN-1 in vivo. The potential effect of riluzole and edaravone two FDA approved drugs for ALS, impacting the structural deviations and stabilization of the mutant PFN-1 is evaluated using in silico tools. Overall, the analysis of data collected reveals structural changes of mutant PFN-1 protein that may explain the neurotoxicity and the reason(s) for possible loss and gain of function of PFN-1 in the neurotoxic model of ALS.

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.