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

Discovery of a new Bcl-2 inhibitor through synthesis, anticancer activity, docking and MD simulations

A database of 300 compounds was virtually screened and docked against Bcl-2 protein; the stability of the best-formed complex was evaluated through Molecular dynamics, the top ten compounds with the best in-silico complexation affinities were synthesized, and their In-vitro cytotoxic activity was examined. Thiazolidinone (4e) and isoxazoline (4a-d) were evaluated in-silico. For further evaluation and examination, we designed and synthesized from naturally occurring (R)-carvone and characterized it via spectroscopic analysis, as well as tested for their anticancer activities towards human cancer cell lines such as HT-1080 (fibrosarcome cancer), MCF-7 and MDA-MB-231 (breast cancer) and A-549 (lung cancer) by using MTT method with Doxorubicin as standard drug. Among them, compound 4d showed the most promising anticancer activity against HT-1080, A-549, MCF-7, and MDA-MB-231 cell lines with IC50 values of 15.59 ± 3.21 µM; 18.32 ± 2.73 µM; 17.28 ± 0.33 µM and 19.27 ± 2.73 µM respectively.Communicated by Ramaswamy H. Sarma.

Unravelling the Relacatib activity against the CTSK proteins causing pycnodysostosis: a molecular docking and dynamics approach

Pycnodysostosis is an atypical autosomal recessive condition of Lysosomal storage disorder that originated due to the deficit of the enzyme Cathepsin K which is vital for normal osteoclast action in bone resorption. Abnormal degradation of type 1 collagen and accumulation of toxic undigested collagen fibers in lysosomes of the osteoclast cells resulting in high bone density, brittle bones, and a short stature is caused in CTSK protein-carrying individuals. The broad aim of this study is to identify the most significant variant through various computational pipelines. This study was initiated by retrieving a total number of thirty-six variants from NCBI, HGMD, and UniProt databases, and the Y283C variant was found to be more significant by various standard computational tools. A structural investigation was performed to understand and gain a better knowledge about the interaction profile for the native (1BY8) and variant (Y283C) with Relacatib (a small-molecule drug that blocks the function of Cathepsin K, an enzyme that has been linked to osteoporosis, osteoarthritis, and other bone-degrading diseases). The interaction profile was analyzed using molecular docking. Relacatib (ligand) had an average binding affinity for both native (-7.16 kcal/mol) and Y283C (-6.76 kcal/mol). Finally, Molecular dynamics simulations were done in duplicates to recognize the variant (Y283C) activity of the protein structure against Relacatib for 100 ns. This study assists in comprehending the most pathogenic amino-acid variant, the ligand interaction with the protein structure, and paves the way for understanding the steadiness of the ligand with the native and selected significant amino-acid variant.Communicated by Ramaswamy H. Sarma.

Neuropharmacological assessment and identification of possible lead compound (apomorphine) from Hygrophila spinosa through in-vivo and in-silico approaches

The aim of this research is to examine possible neurological activity of methanol, ethyl acetate, and aqueous extracts of Hygrophila spinosa and identify possible lead compounds through in silico analysis. In vivo, neuropharmacological activity was evaluated by using four distinct neuropharmacological assessment assays. Previously reported GC-MS data and earlier literature were utilized to identify the phytochemicals present in Hygrophila spinosa. Computational studies notably molecular docking and molecular dynamic simulations were conducted with responsible receptors to assess the stability of the best interacting compound. Pharmacokinetics properties like absorption, distribution, metabolism, excretion, and toxicity were considered to evaluate the drug likeliness properties of the identified compounds. All the in vivo results support the notion that different extracts (methanol, ethyl acetate, and aqueous) of Hygrophila spinosa have significant (*p = 0.05) sedative-hypnotic, anxiolytic, and anti-depressant activity. Among all the extracts, specifically methanol extracts of Hygrophila spinosa (MHS 400 mg/kg.b.w.) showed better sedative, anxiolytic and antidepressant activity than aqueous and ethyl acetate extracts. In silico molecular docking analysis revealed that among 53 compounds 7 compounds showed good binding affinities and one compound, namely apomorphine (CID: 6005), surprisingly showed promising binding affinity to all the receptors . An analysis of molecular dynamics simulations confirmed that apomorphine (CID: 6005) had a high level of stability at the protein binding site. Evidence suggests that Hygrophila spinosa has significant sedative, anxiolytic, and antidepressant activity. In silico analysis revealed that a particular compound (apomorphine) is responsible for this action. Further research is required in order to establish apomorphine as a drug for anxiety, depression, and sleep disorders.Communicated by Ramaswamy H. Sarma.

Computational study of the motor neuron protein KIF5A to identify nsSNPs, bioactive compounds, and its key regulators

Introduction: Kinesin family member 5A (KIF5A) is a motor neuron protein expressed in neurons and involved in anterograde transportation of organelles, proteins, and RNA. Variations in the KIF5A gene that interfere with axonal transport have emerged as a distinguishing feature in several neurodegenerative disorders, including hereditary spastic paraplegia (HSP10), Charcot-Marie-Tooth disease type 2 (CMT2), and Amyotrophic Lateral Sclerosis (ALS). Methods: In this study, we implemented a computational structural and systems biology approach to uncover the role of KIF5A in ALS. Using the computational structural biology method, we explored the role of non-synonymous Single Nucleotide Polymorphism (nsSNPs) in KIF5A. Further, to identify the potential inhibitory molecule against the highly destabilizing structure variant, we docked 24 plant-derived phytochemicals involved in ALS. Results: We found KIF5AS291F variant showed the most structure destabilizing behavior and the phytocompound "epigallocatechin gallate" showed the highest binding affinity (-9.0 Kcal/mol) as compared to wild KIF5A (-8.4 Kcal/mol). Further, with the systems biology approach, we constructed the KIF5A protein-protein interaction (PPI) network to identify the associated Kinesin Families (KIFs) proteins, modules, and their function. We also constructed a transcriptional and post-transcriptional regulatory network of KIF5A. With the network topological parameters of PPIN (Degree, Bottleneck, Closeness, and MNC) using CytoHubba and computational knock-out experiment using Network Analyzer, we found KIF1A, 5B, and 5C were the significant proteins. The functional modules were highly enriched with microtubule motor activity, chemical synaptic transmission in neurons, GTP binding, and GABA receptor activity. In regulatory network analysis, we found KIF5A post-transcriptionally down-regulated by miR-107 which is further transcriptionally up-regulated by four TFs (HIF1A, PPARA, SREBF1, and TP53) and down-regulated by three TFs (ZEB1, ZEB2, and LIN28A). Discussion: We concluded our study by finding a crucial variant of KIF5A and its potential therapeutic target (epigallocatechin gallate) and KIF5A associated significant genes with important regulators which could decrypt the novel therapeutics in ALS and other neurodegenerative diseases.

Mechanism Exploration of Amyloid-β-42 Disaggregation by Single-Chain Variable Fragments of Alzheimer's Disease Therapeutic Antibodies

Alzheimer's disease (AD) is a specific neurodegenerative disease. This study adopts single-chain variable fragments (scFvs) as a potential immunotherapeutic precursor for AD. According to the remarkable effects of monoclonal antibodies, such as the depolymerization or promotion of Aβ42 efflux by Crenezumab, Solanezumab, and 12B4, it is attractive to prepare corresponding scFvs targeting amyloid-β-42 protein (Aβ42) and investigate their biological activities. Crenezumab-like scFv (scFv-C), Solanezumab-like scFv (scFv-S), and 12B4-like scFv (scFv-12B4) were designed and constructed. The thermal stabilities and binding ability to Aβ42 of scFv-C, scFv-S, and scFv-12B4 were evaluated using unfolding profile and enzyme-linked immunosorbent assay. As the results indicated that scFv-C could recognize Aβ42 monomer/oligomer and promote the disaggregation of Aβ42 fiber as determined by the Thioflavin-T assay, the potential mechanism of its interaction with Aβ42 was investigated using molecular dynamics analysis. Interactions involving hydrogen bonds and salt bonds were predicted between scFv-C and Aβ42 pentamer, suggesting the possibility of inhibiting further aggregation of Aβ42. The successfully prepared scFvs, especially scFv-C, with favorable biological activity targeting Aβ42, might be developed for a potentially efficacious clinical application for AD.

New molecular insights into the A218V variant impact on the steroidogenic acute regulatory protein (STAR) associated with 46, XY disorders of sexual development

Disorders of sexual development (DSD) are an abnormal congenital conditions associated with atypical development of the urogenital tract and external genital structures. The steroidogenic acute regulatory (STAR) gene, associated with congenital lipoid adrenal hyperplasia (CLAH), is included in the targeted gene panel for the DSD diagnosis. Therefore, the genetic alterations of the STAR gene and their molecular effect were examined in the CLAH patients affected with DSD. Ten different Iranian families including twelve male pseudo-hermaphroditism patients with CLAH phenotype were studied using genetic linkage screening and STAR gene sequencing in the linked families to the STAR locus. Furthermore, the structural, dynamical, and functional impacts of the variants on the STAR in silico were analyzed. Sanger sequencing showed the pathogenic variant p.A218V in STAR gene, as the first report in Iranian population. Moreover, modeling and simulation analysis were performed using tools such as radius of gyration, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and molecular docking showed that p.A218V variant affects the residues interaction in cholesterol-binding site and the proper folding of STAR through increasing H-bound and the amount of α-Helix, deceasing total flexibility and changing fluctuations in some residues, resulting in reduced steroidogenic activity of the STAR protein. The study characterized the structural and functional changes of STAR caused by pathogenic variant p.A218V. It leads to limited cholesterol-binding activity of STAR, ultimately leading to the CLAH disease. Molecular dynamics simulation of STAR variants could help explain different clinical manifestations of CLAH disease.

Structure restoration and aggregate inhibition of V30M mutant transthyretin protein by potential quinoline molecules

Transthyretin (TTR) is a tetrameric protein found in human plasma and cerebrospinal fluid that functions as a transporter of thyroxine (T4) and retinol. A mutation resulting in the substitution of valine to methionine at position 30 (V30M) is the most common mutation that destabilizes the tetramer structure of TTR protein resulting in a fatal neuropathy known as TTR amyloidosis. The V30M TTR-induced neuropathy can be inhibited through stabilization of the TTR tetramer by the binding of small molecules. We accessed the potential of in-house synthesized quinoline molecules to stabilize the V30M TTR structure and analyzed the impact of protein-ligand interactions through molecular docking, molecular dynamics (MD) simulations, steered MD, and umbrella sampling simulations. This study revealed that the binding of quinoline molecules reverted back the structural changes including the residual flexibility, changes in secondary structural elements, and also restored the alterations in the electrostatic surface potential induced by the V30M mutation. Further, the top-most 4G and 4R molecules were compared with an FDA-approved drug (Tafamidis) and a reference quinoline molecule 14C. Here, we intend to suggest that the quinoline molecules could revert the structural changes, cease tetramer dissociation, prevent abnormal oligomerization and therefore could be developed as an effective therapeutics against TTR amyloidosis.

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, β-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of β-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (N^α-H---N^δ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

Virtual Screening-Based Drug Development for the Treatment of Nervous System Diseases

The incidence rate of nervous system diseases has increased in recent years. Nerve injury or neurodegenerative diseases usually cause neuronal loss and neuronal circuit damage, which seriously affect motor nerve and autonomic nervous function. Therefore, safe and effective treatment is needed. As traditional drug research becomes slower and more expensive, it is vital to enlist the help of cutting- edge technology. Virtual screening (VS) is an attractive option for the identification and development of promising new compounds with high efficiency and low cost. With the assistance of computer- aided drug design (CADD), VS is becoming more and more popular in new drug development and research. In recent years, it has become a reality to transform non-neuronal cells into functional neurons through small molecular compounds, which provides a broader application prospect than transcription factor-mediated neuronal reprogramming. This review mainly summarizes related theory and technology of VS and the drug research and development using VS technology in nervous system diseases in recent years, and focuses more on the potential application of VS technology in neuronal reprogramming, thus facilitating new drug design for both prevention and treatment of nervous system diseases.

Potential roles of the endoplasmic reticulum stress pathway in amyotrophic lateral sclerosis

The endoplasmic reticulum (ER) is a major organelle involved in protein quality control and cellular homeostasis. ER stress results from structural and functional dysfunction of the organelle, along with the accumulation of misfolded proteins and changes in calcium homeostasis, it leads to ER stress response pathway such as unfolded protein response (UPR). Neurons are particularly sensitive to the accumulation of misfolded proteins. Thus, the ER stress is involved in neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, prion disease and motor neuron disease (MND). Recently, the complex involvement of ER stress pathways has been demonstrated in experimental models of amyotrophic lateral sclerosis (ALS)/MND using pharmacological and genetic manipulation of the unfolded protein response (UPR), an adaptive response to ER stress. Here, we aim to provide recent evidence demonstrating that the ER stress pathway is an essential pathological mechanism of ALS. In addition, we also provide therapeutic strategies that can help treat diseases by targeting the ER stress pathway.

A novel causative functional mutation in GATA6 gene is responsible for familial dilated cardiomyopathy as supported by in silico functional analysis

Dilated cardiomyopathy (DCM), one of the most common types of cardiomyopathies has a heterogeneous nature and can be seen in Mendelian forms. Next Generation Sequencing is a powerful tool for identifying novel variants in monogenic disorders. We used whole-exome sequencing (WES) and Sanger sequencing techniques to identify the causative mutation of DCM in an Iranian pedigree. We found a novel variant in the GATA6 gene, leading to substituting Histidine by Tyrosine at position 329, observed in all affected family members in the pedigree, whereas it was not established in any of the unaffected ones. We hypothesized that the H329Y mutation may be causative for the familial pattern of DCM in this family. The predicted models of GATA6 and H329Y showed the high quality according to PROCHECK and ERRAT. Nonetheless, simulation results revealed that the protein stability decreased after mutation, while the flexibility may have been increased. Hence, the mutation led to the increased compactness of GATA6. Overall, these data indicated that the mutation could affect the protein structure, which may be related to the functional impairment of GATA6 upon H329Y mutation, likewise their involvement in pathologies. Further functional investigations would help elucidating the exact mechanism.

Comparative molecular dynamics analyses on PIK3CA hotspot mutations with PI3Kα specific inhibitors and ATP

PI3K pathway is heavily emphasized in cancer where PIK3CA, which encodes for the p110α subunit of PI3Kα, presents itself as the second most common mutated gene. A lot of effort has been put in developing PI3K inhibitors, opening promising avenues for the treatment of cancer. Among these, PI3Kα specific inhibitor alpelisib was approved by FDA for breast cancer and other α-isoform specific inhibitors such as inavolisib and serabelisib reached clinical trials. However, the mode of action of these inhibitors on mutated PI3Kα and how they interact with mutant structures has not been fully elucidated yet. In this study, we are revealing the calculated interactions and binding affinities of these inhibitors within the context of PIK3CA hotspot mutations (E542K, E545K and H1047R) by employing molecular dynamics (MD) simulations. We performed principal component analysis to understand the motions of the protein complex during our simulations and also checked the correlated motions of all amino acids. Binding affinity calculations with MM-PBSA confirmed the consistent binding of alpelisib across mutations and revealed relatively higher affinities for inavolisib towards wild-type and H1047R mutant structures in comparison to other inhibitors. On the other hand, E542K mutation significantly impaired the interaction of inavolisib and serabelisib with PI3Kα. We also investigated the structural relationship of the natural ligand ATP with PI3Kα, and interestingly realized a significant reduction in binding affinity for the mutants, with potentially unexpected implications on the mechanisms that render these mutations oncogenic. Moreover, correlated motions of all residues were generally higher for ATP except the H1047R mutation which exhibited a distinguishable reduction. The results presented here could be guiding for pre-clinical and clinical studies of personalized medicine where individual mutations are a strong consideration point.

Histidine tautomerism dependent conformational transitions driven aggregation of profilin-1: Implications in amyotrophic lateral sclerosis

Aggregation of profilin-1 (PFN1) causes a fatal neurodegenerative disease, familial amyotrophic lateral sclerosis (fALS). Histidine (His) tautomerism has been linked to the formation of fibril aggregation causing neurodegenerative disease. Characterization of intermediate species that form during aggregation is crucial, however, this has proven very challenging for experimentalists due to their transient nature. Hence, molecular dynamics (MD) simulations have been performed on the His tautomeric isomers εε, εδ, δε, and δδ of PFN1 to explain the structural changes and to correlate them with its aggregation propensity. MD simulations show that His133 presumably plays a major role in the aggregation of PFN1 upon His tautomerism compared to His119. Further, the formation of a new 3₁₀-helix is observed in εε and δε but 3₁₀-helix is not observed in δδ and εδ isomers. In addition, our findings unveil that β-sheet dominating conformations are observed in His119(δ)-His133(δ) δδ isomer of PFN1 with significant antiparallel β-sheets between residues T15-G23, S29-A33, L63-L65, Q68-S76, F83-T89, T97-T105, and K107-K115, suggesting a novel aggregation mechanism possibly occur for the formation of PFN1 aggregates. Overall, these results propose that MD simulations of PFN1 His tautomers can provide a detailed microscopic understanding of the aggregation mechanisms which are hard to probe through experiments.

Exploration of the Interactions between Maltase–Glucoamylase and Its Potential Peptide Inhibitors by Molecular Dynamics Simulation

Diabetes mellitus, a chronic metabolic disorder, represents a serious threat to human health. The gut enzyme maltase–glucoamylase (MGAM) has attracted considerable attention as a potential therapeutic target for the treatment of type 2 diabetes. Thus, developing novel inhibitors of MGAM holds the promise of improving clinical management. The dipeptides, Thr-Trp (TW) and Trp-Ala (WA), are known inhibitors of MGAM; however, studies on how they interact with MGAM are lacking. The work presented here explored these interactions by utilizing molecular docking and molecular dynamics simulations. Results indicate that the active center of the MGAM could easily accommodate the flexible peptides. Interactions involving hydrogen bonds, cation-π, and hydrophobic interactions are predicted between TW/WA and residues including Tyr1251, Trp1355, Asp1420, Met1421, Glu1423, and Arg1510 within MGAM. The electrostatic energy was recognized as playing a dominant role in both TW-MGAM and WA-MGAM systems. The binding locations of TW/WA are close to the possible acid-base catalytic residue Asp1526 and might be the reason for MGAM inhibition. These findings provide a theoretical structural model for the development of future inhibitors.

Synthesis of benzimidazole/triphenylamine-based compounds, evaluation of their bioactivities and an in silico study with receptor tyrosine kinases

The antiproliferative activity of AL-1 against various cancer cells indicated the applicability of the BI-TPA-based compound as a potential multi-cancer inhibitor.

Delineating binding potential, stability of Sulforaphane-N-acetyl-cysteine in the active site of histone deacetylase 2 and testing its cytotoxicity against distinct cancer lines through stringent molecular dynamics, DFT and cell-based assays

Histone deacetylase 2 (HDAC2), an isozyme of Class I HDACs has potent imputations in actuating neurodegenerative signaling. Currently, there are sizeable therapeutic disquiets with the use of synthetic histone deacetylase inhibitors in disease management. This strongly suggests the unfulfilled medical necessity of plant substitutes for therapeutic intervention. Sulforaphane-N-acetyl-cysteine (SFN-N-acetylcysteine or SFN-NAC), a sulforaphane metabolite has shown significantly worthier activity against HDACs under in vitro conditions. However, the atomistic studies of SFN-NAC against HDAC2 are currently lacking. Thus, the present study employed a hybrid strategy including extra-precision (XP) grid-based flexible molecular docking, molecular mechanics generalized born surface area (MM-GBSA), e-Pharmacophores method, and molecular dynamics simulation for exploring the binding strengh, mode of interaction, e-Pharmacophoric features, and stability of SFN-NAC towards HDAC2. Further, the globally acknowledged density functional theory (DFT) study was performed on SFN-NAC and entinostat individually in complex state with HDAC2. Apart from this, these inhibitors were tested against three distinct cancer cell models and one transformed cell line for cytotoxic activity. Moreover, double mutant of HDAC2 was generated and the binding orientation and interaction of SFN-NAC was scrutinized in this state. On the whole, this study unbosomed and explained the comparatively higher binding affinity of entinostat for HDAC2 and its wide spectrum cytotoxicity than SFN-NAC.