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

Frustration in physiology and molecular medicine

Molecules provide the ultimate language in terms of which physiology and pathology must be understood. Myriads of proteins participate in elaborate networks of interactions and perform chemical activities coordinating the life of cells. To perform these often amazing tasks, proteins must move and we must think of them as dynamic ensembles of three dimensional structures formed first by folding the polypeptide chains so as to minimize the conflicts between the interactions of their constituent amino acids. It is apparent however that, even when completely folded, not all conflicting interactions have been resolved so the structure remains 'locally frustrated'. Over the last decades it has become clearer that this local frustration is not just a random accident but plays an essential part of the inner workings of protein molecules. We will review here the physical origins of the frustration concept and review evidence that local frustration is important for protein physiology, protein-protein recognition, catalysis and allostery. Also, we highlight examples showing how alterations in the local frustration patterns can be linked to distinct pathologies. Finally we explore the extensions of the impact of frustration in higher order levels of organization of systems including gene regulatory networks and the neural networks of the brain.

Frustration In Physiology And Molecular Medicine

Molecules provide the ultimate language in terms of which physiology and pathology must be understood. Myriads of proteins participate in elaborate networks of interactions and perform chemical activities coordinating the life of cells. To perform these often amazing tasks, proteins must move and we must think of them as dynamic ensembles of three dimensional structures formed first by folding the polypeptide chains so as to minimize the conflicts between the interactions of their constituent amino acids. It is apparent however that, even when completely folded, not all conflicting interactions have been resolved so the structure remains 'locally frustrated'. Over the last decades it has become clearer that this local frustration is not just a random accident but plays an essential part of the inner workings of protein molecules. We will review here the physical origins of the frustration concept and review evidence that local frustration is important for protein physiology, protein-protein recognition, catalysis and allostery. Also, we highlight examples showing how alterations in the local frustration patterns can be linked to distinct pathologies. Finally we explore the extensions of the impact of frustration in higher order levels of organization of systems including gene regulatory networks and the neural networks of the brain.

Multi-scale in silico analysis of the phase separation behavior of FUS mutants

Fused in sarcoma (FUS) is an intrinsically disordered RNA-binding protein that helps to regulate transcription and RNA transport while reversibly assembling into membraneless organelles (MLOs). Some mutations of FUS can promote irreversible aggregation, contributing to neurodegenerative diseases. We previously reported a multi-scale computational framework combining a series of molecular dynamics simulations (MD) followed by lattice Monte Carlo (MC) simulations to describe the tendency and dynamics of the assembly and disassembly of intrinsically disordered proteins (IDPs) using wild-type (WT)-FUS as an illustrative example. In this study, we utilized our computational model to simulate three FUS mutants widely experimented with glycine point mutation G156E, arginine point mutation R244C, and deletion of the C-terminal nuclear localization signal (ΔNLS). MD simulation results conveyed that G156E has improved sticker contact probability compared to WT-FUS, while R244C has slightly lower contact probability, which is also complemented by change of net interactions according to the molecular mechanics Poisson Boltzmann surface area (MMPBSA) method. The MC simulation results revealed that G156E has a higher aggregation propensity than the WT-FUS, while ΔNLS has more liquid-like assemblies. R244C demonstrated higher dynamics at the beginning, while over the evolution of MC simulations, it tends to aggregate compared to WT-FUS. In addition, the G156E mutant has more stable protein aggregates, lacking the rapid dynamics shown in all other scenarios. From the peak height of radial distribution functions (RDFs) of the assemblies, the phase separation propensity in ascending order is ΔNLS < FUS-WT < R244C < G156E. Moreover, interpreting the dynamic assembly propensity (DAP) parameter over time, the fluidity of the assemblies in ascending order is G156E < FUS-WT < R244C < ΔNLS. The results obtained from this study support that the computational model is able to predict the effect of mutation down to single amino acid substitution on the phase separation behavior of FUS. This efficient in silico method can be generalized to investigate the phase separation propensity of other IDPs and their mutants.

Exploring the Structural and Functional Consequences of Deleterious Missense Nonsynonymous SNPs in the EPOR Gene: A Computational Approach

BACKGROUND

Mutations in the EPOR gene can disrupt its normal signaling pathways, leading to hematological disorders such as polycythemia vera and other myeloproliferative diseases.

METHODOLOGY

In this study, a range of bioinformatics tools, including SIFT, PolyPhen-2, SNAP2, SNPs & Go, PhD-SNP, I-Mutant2.0, MuPro, MutPred, ConSurf, HOPE, and Interpro were used to assess the deleterious effects of missense nonsynonymous single nucleotide polymorphisms (nsSNPs) on protein structure and function. Furthermore, molecular dynamics simulations (MDS) were conducted to assess the structural deviations of the identified mutant variants in comparison to the wild type.

RESULTS

The results identified two nsSNPs, R223P and G302S, as deleterious, significantly affecting protein structure and function. Both substitutions occur in functionally conserved regions and are predicted to be pathogenic, associated with altered molecular mechanisms. The MDSs indicated that while the wild-type EPOR maintained optimal stability, the G302S and R223P variants exhibited substantial deviations, adversely affecting overall protein stability and compactness.

CONCLUSIONS

The computational analysis of missense nsSNPs in the EPOR gene identified two missense SNPs, R223P and G302S, as deleterious, occurring at highly conserved regions, and having substantial effects on erythropoietin receptor (EPO-R) protein structure and function, suggesting their potential pathogenic consequences.

Allergic potential & molecular mechanism of skin sensitization of cinnamaldehyde under environmental UVB exposure

Fragrance, a key ingredient in cosmetics, often triggers skin allergy causes rashes, itching, dryness, and cracked or scaly skin. Cinnamaldehyde (CA), derived from the bark of the cinnamon tree, used as a fragrance and is a moderate skin sensitizer. CA exhibits strong UVB absorption, its allergic potential and the molecular mechanisms underlying skin sensitization under UVB exposure remain largely unexplored. To investigate the allergic potential and molecular mechanisms of CA-induced skin sensitization under ambient UVB radiation, we employed various alternative in-silico, in-chemico and in-vitro tools. CA under ambient UVB isomerizes from trans to cis CA after 1hr of exposure. Furthermore, DPRA assay and docking with simulation studies demonstrated the enhanced allergic potential of cis-CA. Additionally, our study evaluated intracellular ROS levels and the expression of Nrf2, Catalase, and MMP-2, and 9 in KeratinoSens cells, showing significant upregulation under UVB exposure in the presence of CA. Moreover, our findings indicate that CA activates THP-1 cells co-stimulatory surface marker (CD86) via the activation of intracellular ROS, phagocytosis, and genes of the TLR4 pathway. These insights into the mechanisms uncovered by our study are crucial for managing triggers of allergic skin diseases caused by fragrance use and concurrent exposure to environmental UVB/sunlight.

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.

Computer-aided identification of potential inhibitors against Necator americanus glutathione S-transferase 3

Hookworm infection is caused by the blood-feeding hookworm gastrointestinal nematodes. Its harmful effects include anemia and retarded growth and are common in the tropics. A current control method involves the mass drug administration of synthetic drugs, mainly albendazole and mebendazole. There are however concerns of low efficacy and drug resistance due to their repeated and excessive use. Although, Necator americanus glutathione S-transferase 3 (Na-GST-3) is a notable target, using natural product libraries for computational elucidation of promising leads is underexploited. This study sought to use pharmacoinformatics techniques to identify compounds of natural origins with the potential to be further optimized as promising inhibitors. A compendium of 3182 African natural products together with five known helminth GST inhibitors including Cibacron blue was screened against the active sites of the Na-GST-3 structure (PDB ID: 3W8S). The hit compounds were profiled to ascertain the mechanisms of binding, anthelmintic bioactivity, physicochemical and pharmacokinetic properties. The AutoDock Vina docking protocol was validated by obtaining 0.731 as the area under the curve calculated via the receiver operating characteristics curve. Four compounds comprising ZINC85999636, ZINC35418176, ZINC14825190, and Dammarane Triterpene13 were identified as potential lead compounds with binding energies less than -9.0 kcal/mol. Furthermore, the selected compounds formed key intermolecular interactions with critical residues Tyr95, Gly13 and Ala14. Notably, ZINC85999636, ZINC14825190, and dammarane triterpene13 were predicted as anthelmintics, whilst all the four molecules shared structural similarities with known inhibitors. Molecular modelling showed that the compounds had reasonably good binding free energies. More so, they had high binding affinities when screened against other variants of the Na-GST, namely Na-GST-1 and Na-GST-2. Ligand quality assessment using ligand efficiency dependent lipophilicity, ligand efficiency, ligand efficiency scale and fit quality scale showed the molecules are worthy candidates for further optimization. The inhibitory potentials of the molecules warrant in vitro studies to evaluate their effect on the heme regulation mechanisms.