Cluster Headache Genomewide Association Study and Meta-Analysis Identifies Eight Loci and Implicates Smoking as Causal Risk Factor.

OBJECTIVE: The objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological insights. METHODS: A total of 4,777 cases (3,348 men and 1,429 women) with clinically diagnosed cluster headache were recruited from 10 European and 1 East Asian cohorts. We first performed an inverse-variance genomewide association meta-analysis of 4,043 cases and 21,729 controls of European ancestry. In a secondary trans-ancestry meta-analysis, we included 734 cases and 9,846 controls of East Asian ancestry. Candidate causal genes were prioritized by 5 complementary methods: expression quantitative trait loci, transcriptome-wide association, fine-mapping of causal gene sets, genetically driven DNA methylation, and effects on protein structure. Gene set and tissue enrichment analyses, genetic correlation, genetic risk score analysis, and Mendelian randomization were part of the downstream analyses. RESULTS: The estimated single nucleotide polymorphism (SNP)-based heritability of cluster headache was 14.5%. We identified 9 independent signals in 7 genomewide significant loci in the primary meta-analysis, and one additional locus in the trans-ethnic meta-analysis. Five of the loci were previously known. The 20 genes prioritized as potentially causal for cluster headache showed enrichment to artery and brain tissue. Cluster headache was genetically correlated with cigarette smoking, risk-taking behavior, attention deficit hyperactivity disorder (ADHD), depression, and musculoskeletal pain. Mendelian randomization analysis indicated a causal effect of cigarette smoking intensity on cluster headache. Three of the identified loci were shared with migraine. INTERPRETATION: This first genomewide association study meta-analysis gives clues to the biological basis of cluster headache and indicates that smoking is a causal risk factor. ANN NEUROL 2023;94:713-726.

Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET.

Potentilla nepalensis Hook is a perennial Himalayan medicinal herb of the Rosaceae family. The present study aimed to evaluate biological activities such as the antioxidant, antibacterial, and anticancer activities of roots and shoots of P. nepalensis and its synergistic antibacterial activity with antibacterial drugs. Folin-Ciocalteau and aluminium chloride methods were used for the calculation of total phenolic (TPC) and flavonoid content (TFC). A DPPH radical scavenging assay and broth dilution method were used for the determination of the antioxidant and antibacterial activity of the root and shoot extracts of P. nepalensis. Cytotoxic activity was determined using a colorimetric MTT assay. Further, phytochemical characterization of the root and shoot extracts was performed using the Gas chromatography-mass spectrophotometry (GC-MS) method. The TPC and TFC were found to be higher in the methanolic root extract of P. nepalensis. The methanolic shoot extract of P. nepalensis showed good antioxidant activity, while then-hexane root extract of P. nepalensis showed strong cytotoxic activity against tested SK-MEL-28 cells. Subsequently, in silico molecular docking studies of the identified bioactive compounds predicted potential anticancer properties. This study can lead to the production of new herbal medicines for various diseases employing P. nepalensis, leading to the creation of new medications.

Targeting leukemic stem cell subpopulation in AML using phytochemicals: An in-silico and in-vitro approach.

It has been indicated that leukemic stem cells (LSCs), a subset of leukaemia cells, are responsible for therapy resistance and relapse in acute myeloid leukaemia (AML). Therefore, the current study aimed to discover an LSC biomarker in AML patients and identify a natural compound that may target the same. By performing the different gene expression analyses, we identified 12 up-regulated and 192 down-regulated genes in LSCs of AML compared to normal bone marrow-derived HSCs. Further STRING interaction, GO enrichment and KEGG pathway analysis were carried out to top hub genes. Wilms' tumour-1 (WT1) transcription factor was pointed out as the top hub gene and a potential biomarker for LSCs in AML. For the targeted inhibition of WT1, we performed screening and stimulation of potential natural compounds. The results revealed Gallic acid (GA) and Chlorogenic acid (CA) as promising WT1 inhibitors. In-vitro validation of cytotoxic effects of both GA and CA on THP-1 and HL-60 cell lines suggested that both these compounds inhibited cell proliferation. Still, GA has a more cytotoxic effect compared to CA. Next, we performed cell cycle analysis and apoptosis analysis and found that both compounds arrested cells in G0/G1 phase and induced apoptosis in both cell lines. Surprisingly, a significant decrease in colony formation and cell migration was also observed. However, GA gave more promising results in all cellular assays than CA. Furthermore, we studied the mRNA expression of WT1 and BCL2, which are transcriptionally activated by it. We found that GA significantly downregulated both these genes compared to CA. Our results suggested that GA is a potential inhibitor of WT1 and might be an excellent anti-LSCs natural drug for AML patients.

MHC2AffyPred: A machine-learning approach to estimate affinity of MHC class II peptides based on structural interaction fingerprints.

Understanding how MHC class II (MHC-II) binding peptides with differing lengths exhibit specific interaction at the core and extended sites within the large MHC-II pocket is a very important aspect of immunological research for designing peptides. Certain efforts were made to generate peptide conformations amenable for MHC-II binding and calculate the binding energy of such complex formation but not directed toward developing a relationship between the peptide conformation in MHC-II structures and the binding affinity (BA) (IC50 ). We present here a machine-learning approach to calculate the BA of the peptides within the MHC-II pocket for HLA-DRA1, HLA-DRB1, HLA-DP, and HLA-DQ allotypes. Instead of generating ensembles of peptide conformations conventionally, the biased mode of conformations was created by considering the peptides in the crystal structures of pMHC-II complexes as the templates, followed by site-directed peptide docking. The structural interaction fingerprints generated from such docked pMHC-II structures along with the Moran autocorrelation descriptors were trained using a random forest regressor specific to each MHC-II peptide lengths (9-19). The entire workflow is automated using Linux shell and Perl scripts to promote the utilization of MHC2AffyPred program to any characterized MHC-II allotypes and is made for free access at https://github.com/SiddhiJani/MHC2AffyPred. The MHC2AffyPred attained better performance (correlation coefficient [CC] of .612-.898) than MHCII3D (.03-.594) and NetMHCIIpan-3.2 (.289-.692) programs in the HLA-DRA1, HLA-DRB1 types. Similarly, the MHC2AffyPred program achieved CC between .91 and .98 for HLA-DP and HLA-DQ peptides (13-mer to 17-mer). Further, a case study on MHC-II binding 15-mer peptides of severe acute respiratory syndrome coronavirus-2 showed very close competency in computing the IC50 values compared to the sequence-based NetMHCIIpan v3.2 and v4.0 programs with a correlation of .998 and .570, respectively.

Theaflavin-3-gallate, a natural antagonist for Hsp90: In-silico and in-vitro approach.

Poor prognosis and metastasis have been recognized as the major cause of breast cancer related deaths worldwide. Recent experimental evidence has shown that Hsp90, the prime chaperone, is overexpressed in many cancers and is responsible if reducing the 5-year survival rate of cancer patients. Therefore, targeted inhibition of Hsp90 may be a new and effective way to target cancer as well as enhancing therapeutic outcomes. In the present study, screening and simulation of potential natural compounds result in the identification of theaflavin-3-gallate as a promising inhibitory compound of Hsp90. Further in-vitro validation of the cytotoxic effect of theaflavin-3-gallate in human breast carcinoma cell line MCF7 and normal cell line MCF10A revealed that theaflavin-3-gallate significantly inhibited the cell proliferation of MCF7 cells whereas no cytotoxic effect was observed on MCF10A cells. We also found that theaflavin-3-gallate significantly induced programmed cell death by arresting cells in the G2/M phase of the cell cycle. A significant decrease in cell migration and colony formation by theaflavin-3-gallate treatment was also observed in MCF7 cells. Furthermore, theaflavin-3-gallate significantly downregulated the mRNA expression patterns of the HSP90, MMP9, VEGFA, and SPP1 genes. Collectively, our results demonstrated theaflavin-3-gallate as a potential natural Hsp90 inhibitor that can be used to enhance the therapeutic efficacy of existing breast cancer therapies and improve overall survival of breast cancer patients.

Repurposing of anticancer phytochemicals for identifying potential fusion inhibitor for SARS-CoV-2 using molecular docking and molecular dynamics (MD) simulations.

The viral particle, SARS-CoV-2 is responsible for causing the epidemic of Coronavirus disease 2019 (COVID-19). To combat this situation, numerous strategies are being thought for either creating its antidote, vaccine, or agents that can prevent its infection. For enabling research on these strategies, several target proteins are identified where, Spike (S) protein is of great potential. S-protein interacts with human angiotensin-converting-enzyme-2 (ACE2) for entering the cell. S-protein is a large protein and a portion of it designated as a receptor-binding domain (RBD) is the key region that interacts with ACE2, following to which the viral membrane fuses with the alveolar membrane to enter the human cell. The hypothesis is to identify molecules from the pool of anticancer phytochemicals as a lead possessing the ability to interact and mask the amino acids of RBD, making them unavailable to form associations with ACE2. Such a molecule is termed as 'fusion inhibitor'. We hypothesized to identify fusion inhibitors from the NPACT library of anticancer phytochemicals. For this, all the molecules from the NPACT were screened using molecular docking, the five top hits (Theaflavin, Ginkgetin, Ursolic acid, Silymarin and Spirosolane) were analyzed for essential Pharmacophore features and their ADMET profiles were studied following to which the best two hits were further analyzed for their interaction with RBD using Molecular Dynamics (MD) simulation. Binding free energy calculations were performed using MM/GBSA, proving these phytochemicals containing anticancer properties to serve as fusion inhibitors.Communicated by Ramaswamy H. Sarma.

Computational investigation of ginkgetin and theaflavin as potential inhibitors of heat shock protein 90 (Hsp90).

Heat shock protein 90 (Hsp90) is the prime molecular chaperone found to be overexpressed in cancer cells and pose as an anti-cancer therapeutic drug target for cancer chemotherapy. Even drugs are available which inhibit Hsp90, the associated side effects along with multi-drug regimen necessitate the identification of natural molecules to block the activity of Hsp90. In this present investigation, we performed virtual screening of Hsp90 inhibitors from a curated collection of natural molecules with proven pharmacological effects. This process helped in the identification of the top two scoring ligands, ginkgetin and theaflavin with favorable as well as crucial interactions with the Hsp90 ligand-binding pocket. Molecular dynamics simulations of these two natural molecules exhibited minimal fluctuations in the binding pattern of ginkgetin and theaflavin to Hsp90 which retained crucial contacts throughout the simulation time. We anticipate that ginkgetin and theaflavin could act as potent Hsp90 inhibitors which are under current investigation in our laboratory.Communicated by Ramaswamy H. Sarma.

Identification of potential inhibitors of coronavirus hemagglutinin-esterase using molecular docking, molecular dynamics simulation and binding free energy calculation.

The pandemic outbreak of the Corona viral infection has become a critical global health issue. Biophysical and structural evidence shows that spike protein possesses a high binding affinity towards host angiotensin-converting enzyme 2 and viral hemagglutinin-acetylesterase (HE) glycoprotein receptor. We selected HE as a target in this study to identify potential inhibitors using a combination of various computational approaches such as molecular docking, ADMET analysis, dynamics simulations and binding free energy calculations. Virtual screening of NPACT compounds identified 3,4,5-Trihydroxy-1,8-bis[(2R,3R)-3,5,7-trihydroxy-3,4-dihydro-2H-chromen-2-yl]benzo[7]annulen-6-one, Silymarin, Withanolide D, Spirosolane and Oridonin as potential HE inhibitors with better binding energy. Furthermore, molecular dynamics simulations for 100 ns time scale revealed that most of the key HE contacts were retained throughout the simulations trajectories. Binding free energy calculations using MM/PBSA approach ranked the top-five potential NPACT compounds which can act as effective HE inhibitors.

Truncated stathmin-2 is a marker of TDP-43 pathology in frontotemporal dementia.

No treatment for frontotemporal dementia (FTD), the second most common type of early-onset dementia, is available, but therapeutics are being investigated to target the 2 main proteins associated with FTD pathological subtypes: TDP-43 (FTLD-TDP) and tau (FTLD-tau). Testing potential therapies in clinical trials is hampered by our inability to distinguish between patients with FTLD-TDP and FTLD-tau. Therefore, we evaluated truncated stathmin-2 (STMN2) as a proxy of TDP-43 pathology, given the reports that TDP-43 dysfunction causes truncated STMN2 accumulation. Truncated STMN2 accumulated in human induced pluripotent stem cell-derived neurons depleted of TDP-43, but not in those with pathogenic TARDBP mutations in the absence of TDP-43 aggregation or loss of nuclear protein. In RNA-Seq analyses of human brain samples from the NYGC ALS cohort, truncated STMN2 RNA was confined to tissues and disease subtypes marked by TDP-43 inclusions. Last, we validated that truncated STMN2 RNA was elevated in the frontal cortex of a cohort of patients with FTLD-TDP but not in controls or patients with progressive supranuclear palsy, a type of FTLD-tau. Further, in patients with FTLD-TDP, we observed significant associations of truncated STMN2 RNA with phosphorylated TDP-43 levels and an earlier age of disease onset. Overall, our data uncovered truncated STMN2 as a marker for TDP-43 dysfunction in FTD.

Cardiotonic steroids as potential Na+/K+-ATPase inhibitors - a computational study.

Cardiotonic steroids (CTS) are steroidal drugs, processed from the seeds and dried leaves of the genus Digitalis as well as from the skin and parotid gland of amphibians. The most commonly known CTS are ouabain, digoxin, digoxigenin and bufalin. CTS can be used for safer medication of congestive heart failure and other related conditions due to promising pharmacological and medicinal properties. Ouabain isolated from plants is widely utilized in in vitro studies to specifically block the sodium potassium (Na+/K+-ATPase) pump. For checking, whether ouabain derivatives are robust inhibitors of Na+/K+-ATPase pump, molecular docking simulation was performed between ouabain and its derivatives using YASARA software. The docking energy falls within the range of 8.470 kcal/mol to 7.234 kcal/mol, in which digoxigenin was found to be the potential ligand with the best docking energy of 8.470 kcal/mol. Furthermore, pharmacophore modeling was applied to decipher the electronic features of CTS. Molecular dynamics simulation was also employed to determine the conformational properties of Na+/K+-ATPase-ouabain and Na+/K+-ATPase-digoxigenin complexes with the plausible structural integrity through conformational ensembles for 100 ns which promoted digoxigenin as the most promising CTS for treating conditions of congestive heart failure patients.

Molecular dynamics-assisted pharmacophore modeling of caspase-3-isatin sulfonamide complex: Recognizing essential intermolecular contacts and features of sulfonamide inhibitor class for caspase-3 binding.

The identification of isatin sulfonamide as a potent small molecule inhibitor of caspase-3 had fuelled the synthesis and characterization of the numerous sulfonamide class of inhibitors to optimize for potency. Recent works that relied on the ligand-based approaches have successfully shown the regions of optimizations for sulfonamide scaffold. We present here molecular dynamics-based pharmacophore modeling of caspase-3-isatin sulfonamide crystal structure, to elucidate the essential non-covalent contacts and its associated pharmacophore features necessary to ensure caspase-3 optimal binding. We performed 20ns long dynamics of this crystal structure to extract global conformation states and converted into structure-based pharmacophore hypotheses which were rigorously validated using an exclusive focussed library of experimental actives and inactives of sulfonamide class by Receiver Operating Characteristic (ROC) statistic. Eighteen structure-based pharmacophore hypotheses with better sensitivity and specificity measures (>0.6) were chosen which collectively showed the role of pocket residues viz. Cys163 (S1 sub-site; required for covalent and H bonding with Michael acceptor of inhibitors), His121 (S1; π stack with bicyclic isatin moiety), Gly122 (S1; H bond with carbonyl oxygen) and Tyr204 (S2; π stack with phenyl group of the isatin sulfonamide molecule) as stringent binding entities for enabling caspase-3 optimal binding. The introduction of spatial pharmacophore site points obtained from dynamics-based pharmacophore models in a virtual screening strategy will be helpful to screen and optimize molecules belonging to sulfonamide class of caspase-3 inhibitors.

Humanized mutant FUS drives progressive motor neuron degeneration without aggregation in 'FUSDelta14' knockin mice.

Mutations in FUS are causative for amyotrophic lateral sclerosis with a dominant mode of inheritance. In trying to model FUS-amyotrophic lateral sclerosis (ALS) in mouse it is clear that FUS is dosage-sensitive and effects arise from overexpression per se in transgenic strains. Novel models are required that maintain physiological levels of FUS expression and that recapitulate the human disease-with progressive loss of motor neurons in heterozygous animals. Here, we describe a new humanized FUS-ALS mouse with a frameshift mutation, which fulfils both criteria: the FUS Delta14 mouse. Heterozygous animals express mutant humanized FUS protein at physiological levels and have adult onset progressive motor neuron loss and denervation of neuromuscular junctions. Additionally, we generated a novel antibody to the unique human frameshift peptide epitope, allowing specific identification of mutant FUS only. Using our new FUSDelta14 ALS mouse-antibody system we show that neurodegeneration occurs in the absence of FUS protein aggregation. FUS mislocalization increases as disease progresses, and mutant FUS accumulates at the rough endoplasmic reticulum. Further, transcriptomic analyses show progressive changes in ribosomal protein levels and mitochondrial function as early disease stages are initiated. Thus, our new physiological mouse model has provided novel insight into the early pathogenesis of FUS-ALS.

Receptor-Guided De Novo Design of Dengue Envelope Protein Inhibitors.

Inhibitor design associated with the dynamics of dengue envelope protein at pre-fusion stage is a prominent strategy to interfere fusion transition of dengue virus with the host cell membrane. Receptor-guided de novo inhibitors were designed based on the knowledge of co-crystallized detergent, ß-octyl glucoside. Pharmacophore features distribution showed the preference of aromatic groups with H bonding features connected to aliphatic bulky group as the skeleton for inhibitor design. Molecular dynamic simulations revealed (2R)-2-[(6-amino-1-oxohexan-2-yl)amino]-4-[6-(4-phenylpiperidine-1-yl)-1,2-benzoxazol-3-yl]butanoate as the probable binder which developed extensive conservative interactions despite the local pocket residues movements especially from kl ß-hairpin, the key structural unit for initiating conformational changes required for fusion transition. The electronic and hydrophobic potentials also indicated that butanoate molecule as the initial lead for envelope protein inhibitors.

Structural insights into the theoretical model of Plasmodium falciparum NADH dehydrogenase and its interaction with artemisinin and derivatives: towards global health therapeutics.

It is a continuing quest to uncover the principal molecular targets of malarial parasites to understand the antimalarial activity and mechanism of action of artemisinin, a potent antimalarial. A series of parasite proteins are experimentally validated as potential targets, such as translationally controlled tumor protein (TCTP) and sarco/endoplasmic reticulum membrane calcium ATP-ase (SERCA). The present study addressed the development of a theoretical model of Plasmodium falciparum NADH dehydrogenase with inference from artemisinin in vivo inhibitory activity. We report here the predicted binding modes of artemisinin and its derivatives. The modeled protein resembled the structural architecture of flavoproteins and oxidoreductases, consisting of two Rossmann folds and dedicated binding sites for its cofactors. Docked poses of the ligand dataset revealed its interactions at or near the si face, indicating being activated. This may aid in generation of reactive oxygen species, thereby disrupting the membrane potential of parasite mitochondria and leading to the clearance from the blood. These observations open up new strategies for development of novel therapeutics, or improvement of existing pharmacotherapies against malaria, a major burden for global health.

Epitope-based immunoinformatics and molecular docking studies of nucleocapsid protein and ovarian tumor domain of crimean-congo hemorrhagic Fever virus.

Crimean-Congo hemorrhagic fever virus (CCHFV), the fatal human pathogen is transmitted to humans by tick bite, or exposure to infected blood or tissues of infected livestock. The CCHFV genome consists of three RNA segments namely, S, M, and L. The unusual large viral L protein has an ovarian tumor (OTU) protease domain located in the N terminus. It is likely that the protein may be autoproteolytically cleaved to generate the active virus L polymerase with additional functions. Identification of the epitope regions of the virus is important for the diagnosis, phylogeny studies, and drug discovery. Early diagnosis and treatment of CCHF infection is critical to the survival of patients and the control of the disease. In this study, we undertook different in silico approaches using molecular docking and immunoinformatics tools to predict epitopes which can be helpful for vaccine designing. Small molecule ligands against OTU domain and protein-protein interaction between a viral and a host protein have been studied using docking tools.