Computational drug repurposing for primary hyperparathyroidism.

In hyperparathyroidism (hyperPTH), excessive amounts of PTH are secreted, interfering with calcium regulation in the body. Several drugs can control the disease's side effects, but none of them is an alternative treatment to surgery. Therefore, new drug candidates are necessary. In this study, three computationally repositioned drugs, DG 041, IMD 0354, and cucurbitacin I, are evaluated in an in vitro model of hyperPTH. First, we integrated publicly available transcriptomics datasets to propose drug candidates. Using 3D spheroids derived from a single primary hyperPTH patient, we assessed their in vitro efficacy. None of the proposed drugs affected the viability of healthy cell control (HEK293) or overactive parathyroid cells at the level of toxicity. This behavior was attributed to the non-cancerous nature of the parathyroid cells, establishing the hyperPTH disease model. Cucurbitacin I and IMD 0354 exhibited a slight inverse relationship between increased drug concentrations and cell viability, whereas DG 041 increased viability. Based on these results, further studies are needed on the mechanism of action of the repurposed drugs, including determining the effects of these drugs on cellular PTH synthesis and secretion and on the metabolic pathways that regulate PTH secretion.

Three candidate anticancer drugs were repositioned by integrative analysis of the transcriptomes of species with different regenerative abilities after injury.

Regeneration is a homeostatic process that involves the restoration of cells and body parts. Most of the molecular mechanisms and signalling pathways involved in wound healing, such as proliferation, have also been associated with cancer cell growth, suggesting that cancer is an over/unhealed wound. In this study, we examined differentially expressed genes in spinal cord samples from regenerative organisms (axolotl and zebrafish) and nonregenerative organisms (mouse and rat) compared to intact control spinal cord samples using publicly available transcriptomics data and bioinformatics analyses. Based on these gene signatures, we investigated 3 small compounds, namely cucurbitacin I, BMS-754807, and PHA-793887 as potential candidates for the treatment of cancer. The predicted target genes of the repositioned compounds were mainly enriched with the greatest number of genes in cancer pathways. The molecular docking results on the binding affinity between the repositioned compounds and their target genes are also reported. The repositioned 3 small compounds showed anticancer effect both in 2D and 3D cell cultures using the prostate cancer cell line as a model. We propose cucurbitacin I, BMS-754807, and PHA-793887 as potential anticancer drug candidates. Future studies on the mechanisms associated with the revealed gene signatures and anticancer effects of these three small compunds would allow scientists to develop therapeutic approaches to combat cancer. This research contributes to the evaluation of mechanisms and gene signatures that either limit or cause cancer, and to the development of new cancer therapies by establishing a link between regeneration and carcinogenesis.

Drug Repositioning Identifies Six Drug Candidates for Systemic Autoimmune Diseases by Integrative Analyses of Transcriptomes from Scleroderma, Systemic Lupus Erythematosus, and Sjogren's Syndrome.

The mechanisms of systemic autoimmune diseases (ADs) are still not clearly understood. Understanding the etiology of systemic ADs and identifying new therapeutic targets require a systems science approach. Using publicly available transcriptome data and bioinformatic analysis, we investigated the differential gene expression profiles of patients with scleroderma, systemic lupus erythematosus, and Sjogren's syndrome. Of these common differentially expressed gene signatures, 208 were regulated in the same direction (either upregulated or downregulated in all datasets) and used for drug repositioning. Six small molecule drug candidates (KU-0063794, YM-155 [sepantronium bromide], MST-312 [telomerase inhibitor IX], PLX-4720, ZM 336372, and 528116.cdx [PIK-75]) were discovered by drug repositioning as potential therapeutics for systemic ADs. The Search Tool for Chemical Interactions was used to find the anticipated target genes of the repositioned molecules. The PI3K/AKT pathway topped the list of common enriched pathways with the most anticipated target genes of the six repositioned small molecules. We also report here the molecular docking findings on the binding affinity between the repositioned drug candidates and genes from the protein-protein interaction network modules of anticipated target genes. Taken together, this study provides new insights and opens up new possibilities on both pathogenesis and treatment of systemic ADs through drug repositioning.

Integrative Analysis of Motor Neuron and Microglial Transcriptomes from SOD1G93A Mice Models Uncover Potential Drug Treatments for ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal disease of motor neurons that mainly affects the motor cortex, brainstem, and spinal cord. Under disease conditions, microglia could possess two distinct profiles, M1 (toxic) and M2 (protective), with the M2 profile observed at disease onset. SOD1 (superoxide dismutase 1) gene mutations account for up to 20% of familial ALS cases. Comparative gene expression differences in M2-protective (early) stage SOD1G93A microglia and age-matched SOD1G93A motor neurons are poorly understood. We evaluated the differential gene expression profiles in SOD1G93A microglia and SOD1G93A motor neurons utilizing publicly available transcriptomics data and bioinformatics analyses, constructed biomolecular networks around them, and identified gene clusters as potential drug targets. Following a drug repositioning strategy, 5 small compounds (belinostat, auranofin, BRD-K78930611, AZD-8055, and COT-10b) were repositioned as potential ALS therapeutic candidates that mimic the protective state of microglia and reverse the toxic state of motor neurons. We anticipate that this study will provide new insights into the ALS pathophysiology linking the M2 state of microglia and drug repositioning.

Assessment of 13 in silico pathogenicity methods on cancer-related variants.

Single nucleotide variants (SNVs) are single base substitutions that could influence many biological functions in the cell including gene expression, protein folding, and protein-protein interactions among many others. Thus, predictions of functional effects of cancer-related variants are crucial for drug responses and treatment options in clinical oncology. Experimental identification of these effects could be slow, inefficient, and inconvenient, hence in silico methods are gaining popularity in predicting the variants' effects. There are many studies on the cancer variants, however, up to date, none of these have been aimed to assess the performance metrics of in silico pathogenicity methods on functional relevance of cancer variants obtained from ClinVar. To this end, we examined the pathogenicity predictions of cancer-related variant datasets of 8 cancer types (bladder, breast, colon, colorectal, kidney, liver, lung, and pancreas cancer) retrieved from ClinVar using 13 different in silico methods including SIFT, CADD, FATHMM-weighted, FATHMM-unweighted, GERP++, MetaSVM, Mutation Assessor, MutationTaster, MutPred, PolyPhen-2, Provean, Revel and VEST4. A combination of statistical performance metric analysis, prediction distribution frequency data and ROC curve analysis results have suggested that; among all in silico prediction tools, top three tools with the highest discriminatory power were found to be MutPred (AUC = 0.677), MetaSVM (AUC = 0.645) and Revel (AUC = 0.637).

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review.

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Revisiting allostery in CREB-binding protein (CBP) using residue-based interaction energy.

CREB-binding protein (CBP) is a multi-subunit scaffold protein complex in transcription regulation process, binding and interacting with ligands such as mixed-lineage leukemia (MLL) and c-Myb allosterically. Here in this study, we have revisited the concept of allostery in CBP via residue-based interaction energy calculation based on molecular dynamics (MD) simulations. To this end, we conducted MD simulations of KIX:MLL:c-Myb ternary complex, its binary components and kinase-inducible domain (KID) interacting domain (KIX) backbone. Interaction energy profiles and cross correlation analysis were performed and the results indicated that KIX:MLL and KIX:c-Myb:MLL complexes demonstrate significant similarities according to both analysis methods. Two regions in the KIX backbone were apparent from the interaction energy and cross correlation maps that hold a key to allostery phenomena observed in CBP. While one of these regions are related to the ligand binding residues, the other comprises of L12-G2 loop and α3 helix regions that have been found to have a significant role in allosteric signal propagation. All in all, residue-based interaction energy calculation method is demonstrated to be a valuable calculation technique for the detection of allosteric signal propagation and ligand interaction regions.

Cancer Drug Repositioning by Comparison of Gene Expression in Humans and Axolotl (Ambystoma mexicanum) During Wound Healing.

Urodele amphibians such as the axolotl (Ambystoma mexicanum) display a large capacity for tissue regeneration and remarkable resistance to cancer. As a model organism, axolotl thus offers a unique opportunity for cancer research and anticancer drug discovery, not to mention the discerning mechanisms that underpin controlled cellular growth and regeneration versus cancer. To the best of our knowledge, little is known on comparative gene expression changes during regeneration events such as wound healing in axolotl and humans. Using publicly available transcriptomics data and bioinformatics analyses, we examined the differential gene expression signatures in skin wound samples from axolotl and humans after skin biopsy punch injury, in comparison with intact (uninjured) control skin samples. We identified 95 genes exhibiting a reversal expression pattern between humans and axolotl during the wound healing/regeneration period. These genes were significantly associated with collagen biosynthesis, extracellular matrix organization, PI3K-Akt signaling pathway, immune system response, and apoptotic process. Furthermore, this new gene set exhibited high prognostic performance in discriminating the survival risk in skin-related cancers, including melanoma (hazard ratio [HR] = 8.14, p < 10-30), oral cancer (HR >100, p < 10-12), and head and neck carcinoma (HR = 5.29, p < 10-30). Moreover, considering these gene signatures, we repositioned 11 small molecules as potential anticancer drug candidates indicating reversal effects on upregulated human genes and downregulated axolotl genes or mimicking downregulated human genes and upregulated axolotl genes. We anticipate that this study offers new insights on gene signatures bridging regeneration mechanisms with tumorigenesis and cancer drug repositioning.

BLK pathway-associated rs13277113 GA genotype is more frequent in SLE patients and associated with low gene expression and increased flares.

We aimed to evaluate the relationship between some important genetic variations and expressions of these genes in our SLE population. We also determined their association with clinical parameters. Eighty-four SLE patients (79 F, 5 M) and 105 healthy controls (98 F, 7 M) were included in the study. rs13277113, rs2736340, rs7829816, rs6983130, rs2613310, and rs704853 polymorphisms, gene expressions of Src family kinases (Blk, Hck, Lck, and Lyn), and Syk kinases (Syk, ZAP70) were studied by real-time PCR. The heterozygous genotypic pattern (GA) for rs13277113 polymorphism was more frequent in patients with SLE when compared to that in controls (48.8 vs. 31.4%, p = 0.035). Other genotype variants were similar in SLE patients and controls. In the SLE group, the heterozygous genotype for rs13277113 was significantly less frequent in active SLE patients (58.8 vs. 26.7%, p = 0.01). SLE flares according to the SELENA-SLEDAI flare index were significantly more frequent in GA (rs13277113) (70 vs. 37%) and CT (rs2736340) genotypes (66.7 vs. 35.2%) than those in other genotypes (p values <0.01). The relative expression of Blk gene was significantly decreased in the SLE group as compared to that in controls (0.52 times, 95%CI 0.19-0.85). The gene expressions of Blk and ZAP70 were significantly lower in SLE patients who had flares according to the SELENA-SLEDAI flare index when compared to those in others (p values 0.01 and 0.017). We observed more frequent heterozygous GA genotypic pattern (rs13277113) in our SLE patients compared to that in controls; and it was associated with disease flares. Blk gene expression in SLE was lower, especially in relapsing patients.

Increased frequency of class I and II anti-human leukocyte antigen antibodies in systemic lupus erythematosus and scleroderma and associated factors: a comparative study.

AIM: There is significant autoantibody production in systemic lupus erythematosus (SLE) and scleroderma (SSc); microchimerism is also thought to play a role in pathogenesis. We determined the frequency of anti-HLA antibodies in SLE and SSc patients and evaluated associated clinical factors. METHODS: We included 77 SLE patients, 46 SSc patients and 53 healthy controls into the study. Clinical data about the patients were obtained from hospital records. Anti-human leukocyte (anti-HLA) antigen antibody analysis of sera was performed by applying Lifecodes anti-HLA Class I and Class II Screening kits based on xMAP technology. RESULTS: The frequencies of class I and II anti-HLA antibodies were significantly higher in SLE (27.3% and 41.6%) and SSc (26.1% and 41.3%) groups than in healthy controls (1.9% and 5.7%) (all P < 0.001). Frequencies of thrombocytopenia (P = 0.021), anti-ribonucleoprotein (P = 0.037) and anti-Ro (P = 0.027) were significantly higher in the class I antibody-positive SLE group; however, pericarditis was less frequent (P = 0.05). On the other hand, the class II antibody-positive SLE group had more frequent anti-ribosomal P antibody (P = 0.038), but less frequent active disease (P = 0.038). In the SSc group, class I antibody-positive patients had more frequent digital ulcers (P = 0.048) and anti-centromere antibodies (P = 0.01). There was no association of anti-HLA antibodies with pulmonary hypertension and interstitial fibrosis in SSc patients. CONCLUSIONS: Both class I and class II antibodies were found to be significantly increased in SLE and SSc. Rather than major organ involvement, anti-HLA antibodies were associated with the presence of other antibodies in both diseases.

PECAM-1 gene polymorphisms and soluble PECAM-1 level in rheumatoid arthritis and systemic lupus erythematosus patients: any link with clinical atherosclerotic events?

Genetic polymorphisms of platelet endothelial cell adhesion molecule-1 (PECAM-1) were found to play roles in atherosclerotic events. We determined PECAM-1 polymorphisms, soluble PECAM-1, and CD40L levels in rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) and evaluated their associations with clinical atherosclerotic complications. We included 100 RA patients, 81 SLE patients, and 94 healthy controls. The clinical features about the patients were obtained from medical records. Past cardiovascular complications were recorded. The most frequent gene polymorphisms of PECAM-1 were studied in our genetics laboratory. Soluble PECAM-1 and CD40L levels in serum were determined with ELISA. The frequencies of 373C (rs668) and 1688A (rs12953) alleles were higher in RA patients when compared to controls (p values, 0.028 and 0.016). RA and SLE patients had significantly higher allele frequencies for 2008A (rs1131012) when compared to controls (p values, 0.016 and 0.001). SLE patients had significantly more frequent AA genotype for rs1131012 polymorphism than RA patients and controls (p values, 0.007 and <0.001). Soluble PECAM-1 level was significantly higher in RA patients than in SLE patients and healthy controls (p values <0.001). Atherosclerotic complications were more frequent in SLE patients with AG genotype (rs12953) than those with AA genotype (p = 0.021). SLE patients with CC genotype (rs668) had a significantly lower frequency of atherosclerotic complications than those with CG genotype (p = 0.045). Nevertheless, in multivariate analysis, there was no association between genotype and atherosclerotic complications. We found associations between various PECAM-1 polymorphisms in RA and SLE; PECAM-1 and soluble CD40 ligand (sCD40L) levels were significantly higher in RA patients than in SLE and control groups. PECAM-1 polymorphisms in SLE were protective against atherosclerotic complications.