Lysosomal cystine export regulates mTORC1 signaling to guide kidney epithelial cell fate specialization.

Differentiation is critical for cell fate decisions, but the signals involved remain unclear. The kidney proximal tubule (PT) cells reabsorb disulphide-rich proteins through endocytosis, generating cystine via lysosomal proteolysis. Here we report that defective cystine mobilization from lysosomes through cystinosin (CTNS), which is mutated in cystinosis, diverts PT cells towards growth and proliferation, disrupting their functions. Mechanistically, cystine storage stimulates Ragulator-Rag GTPase-dependent recruitment of mechanistic target of rapamycin complex 1 (mTORC1) and its constitutive activation. Re-introduction of CTNS restores nutrient-dependent regulation of mTORC1 in knockout cells, whereas cell-permeant analogues of L-cystine, accumulating within lysosomes, render wild-type cells resistant to nutrient withdrawal. Therapeutic mTORC1 inhibition corrects lysosome and differentiation downstream of cystine storage, and phenotypes in preclinical models of cystinosis. Thus, cystine serves as a lysosomal signal that tailors mTORC1 and metabolism to direct epithelial cell fate decisions. These results identify mechanisms and therapeutic targets for dysregulated homeostasis in cystinosis.

Identification of dual-purpose therapeutic targets implicated in aging and glioblastoma multiforme using PandaOmics - an AI-enabled biological target discovery platform.

Glioblastoma Multiforme (GBM) is the most aggressive and most common primary malignant brain tumor. The age of GBM patients is considered as one of the disease's negative prognostic factors and the mean age of diagnosis is 62 years. A promising approach to preventing both GBM and aging is to identify new potential therapeutic targets that are associated with both conditions as concurrent drivers. In this work, we present a multi-angled approach of identifying targets, which takes into account not only the disease-related genes but also the ones important in aging. For this purpose, we developed three strategies of target identification using the results of correlation analysis augmented with survival data, differences in expression levels and previously published information of aging-related genes. Several studies have recently validated the robustness and applicability of AI-driven computational methods for target identification in both cancer and aging-related diseases. Therefore, we leveraged the AI predictive power of the PandaOmics TargetID engine in order to rank the resulting target hypotheses and prioritize the most promising therapeutic gene targets. We propose cyclic nucleotide gated channel subunit alpha 3 (CNGA3), glutamate dehydrogenase 1 (GLUD1) and sirtuin 1 (SIRT1) as potential novel dual-purpose therapeutic targets to treat aging and GBM.

AlphaFold accelerates artificial intelligence powered drug discovery: efficient discovery of a novel CDK20 small molecule inhibitor.

The application of artificial intelligence (AI) has been considered a revolutionary change in drug discovery and development. In 2020, the AlphaFold computer program predicted protein structures for the whole human genome, which has been considered a remarkable breakthrough in both AI applications and structural biology. Despite the varying confidence levels, these predicted structures could still significantly contribute to structure-based drug design of novel targets, especially the ones with no or limited structural information. In this work, we successfully applied AlphaFold to our end-to-end AI-powered drug discovery engines, including a biocomputational platform PandaOmics and a generative chemistry platform Chemistry42. A novel hit molecule against a novel target without an experimental structure was identified, starting from target selection towards hit identification, in a cost- and time-efficient manner. PandaOmics provided the protein of interest for the treatment of hepatocellular carcinoma (HCC) and Chemistry42 generated the molecules based on the structure predicted by AlphaFold, and the selected molecules were synthesized and tested in biological assays. Through this approach, we identified a small molecule hit compound for cyclin-dependent kinase 20 (CDK20) with a binding constant Kd value of 9.2 ± 0.5 µM (n = 3) within 30 days from target selection and after only synthesizing 7 compounds. Based on the available data, a second round of AI-powered compound generation was conducted and through this, a more potent hit molecule, ISM042-2-048, was discovered with an average Kd value of 566.7 ± 256.2 nM (n = 3). Compound ISM042-2-048 also showed good CDK20 inhibitory activity with an IC50 value of 33.4 ± 22.6 nM (n = 3). In addition, ISM042-2-048 demonstrated selective anti-proliferation activity in an HCC cell line with CDK20 overexpression, Huh7, with an IC50 of 208.7 ± 3.3 nM, compared to a counter screen cell line HEK293 (IC50 = 1706.7 ± 670.0 nM). This work is the first demonstration of applying AlphaFold to the hit identification process in drug discovery.

Radioprotectors.org: an open database of known and predicted radioprotectors.

The search for radioprotectors is an ambitious goal with many practical applications. Particularly, the improvement of human radioresistance for space is an important task, which comes into view with the recent successes in the space industry. Currently, all radioprotective drugs can be divided into two large groups differing in their effectiveness depending on the type of exposure. The first of these is radioprotectors, highly effective for pulsed, and some types of relatively short exposure to irradiation. The second group consists of long-acting radioprotectors. These drugs are effective for prolonged and fractionated irradiation. They also protect against impulse exposure to ionizing radiation, but to a lesser extent than short-acting radioprotectors. Creating a database on radioprotectors is a necessity dictated by the modern development of science and technology. We have created an open database, Radioprotectors.org, containing an up-to-date list of substances with proven radioprotective properties. All radioprotectors are annotated with relevant chemical and biological information, including transcriptomic data, and can be filtered according to their properties. Additionally, the performed transcriptomics analysis has revealed specific transcriptomic profiles of radioprotectors, which should facilitate the search for potent radioprotectors.

PIM1 kinase promotes gallbladder cancer cell proliferation via inhibition of proline-rich Akt substrate of 40 kDa (PRAS40).

Gallbladder cancer (GBC) is a rare malignancy, associated with poor disease prognosis with a 5-year survival of only 20%. This has been attributed to late presentation of the disease, lack of early diagnostic markers and limited efficacy of therapeutic interventions. Elucidation of molecular events in GBC can contribute to better management of the disease by aiding in the identification of therapeutic targets. To identify aberrantly activated signaling events in GBC, tandem mass tag-based quantitative phosphoproteomic analysis of five GBC cell lines was carried out. Proline-rich Akt substrate 40 kDa (PRAS40) was one of the proteins found to be hyperphosphorylated in all the invasive GBC cell lines. Tissue microarray-based immunohistochemical labeling of phospho-PRAS40 (T246) revealed moderate to strong staining in 77% of the primary gallbladder adenocarcinoma cases. Regulation of PRAS40 activity by inhibiting its upstream kinase PIM1 resulted in a significant decrease in cell proliferation, colony forming and invasive ability of GBC cells. Our results support the role of PRAS40 phosphorylation in GBC cell survival and aggressiveness. This study also elucidates phospho-PRAS40 as a clinical marker in GBC and the role of PIM1 as a therapeutic target in GBC.

Use of deep neural network ensembles to identify embryonic-fetal transition markers: repression of COX7A1 in embryonic and cancer cells.

Here we present the application of deep neural network (DNN) ensembles trained on transcriptomic data to identify the novel markers associated with the mammalian embryonic-fetal transition (EFT). Molecular markers of this process could provide important insights into regulatory mechanisms of normal development, epimorphic tissue regeneration and cancer. Subsequent analysis of the most significant genes behind the DNNs classifier on an independent dataset of adult-derived and human embryonic stem cell (hESC)-derived progenitor cell lines led to the identification of COX7A1 gene as a potential EFT marker. COX7A1, encoding a cytochrome C oxidase subunit, was up-regulated in post-EFT murine and human cells including adult stem cells, but was not expressed in pre-EFT pluripotent embryonic stem cells or their in vitro-derived progeny. COX7A1 expression level was observed to be undetectable or low in multiple sarcoma and carcinoma cell lines as compared to normal controls. The knockout of the gene in mice led to a marked glycolytic shift reminiscent of the Warburg effect that occurs in cancer cells. The DNN approach facilitated the elucidation of a potentially new biomarker of cancer and pre-EFT cells, the embryo-onco phenotype, which may potentially be used as a target for controlling the embryonic-fetal transition.

Data aggregation at the level of molecular pathways improves stability of experimental transcriptomic and proteomic data.

High throughput technologies opened a new era in biomedicine by enabling massive analysis of gene expression at both RNA and protein levels. Unfortunately, expression data obtained in different experiments are often poorly compatible, even for the same biologic samples. Here, using experimental and bioinformatic investigation of major experimental platforms, we show that aggregation of gene expression data at the level of molecular pathways helps to diminish cross- and intra-platform bias otherwise clearly seen at the level of individual genes. We created a mathematical model of cumulative suppression of data variation that predicts the ideal parameters and the optimal size of a molecular pathway. We compared the abilities to aggregate experimental molecular data for the 5 alternative methods, also evaluated by their capacity to retain meaningful features of biologic samples. The bioinformatic method OncoFinder showed optimal performance in both tests and should be very useful for future cross-platform data analyses.

MiRImpact, a new bioinformatic method using complete microRNA expression profiles to assess their overall influence on the activity of intracellular molecular pathways.

MicroRNAs (miRs) are short noncoding RNA molecules that regulate expression of target mRNAs. Many published sources provide information about miRs and their targets. However, bioinformatic tools elucidating higher level impact of the established total miR profiles, are still largely missing. Recently, we developed a method termed OncoFinder enabling quantification of the activities of intracellular molecular pathways basing on gene expression data. Here we propose a new technique, MiRImpact, which enables to link miR expression data with its estimated outcome on the regulation of molecular pathways, like signaling, metabolic, cytoskeleton rearrangement, and DNA repair pathways. MiRImpact uses OncoFinder rationale for pathway activity calculations, with the major distinctions that (i) it deals with the concentrations of miRs--known regulators of gene products participating in molecular pathways, and (ii) miRs are considered as negative regulators of target molecules, if other is not specified. MiRImpact operates with 2 types of databases: for molecular targets of miRs and for gene products participating in molecular pathways. We applied MiRImpact to compare regulation of human bladder cancer-specific signaling pathways at the levels of mRNA and miR expression. We took 2 most complete alternative databases of experimentally validated miR targets--miRTarBase and DianaTarBase, and an OncoFinder database featuring 2725 gene products and 271 signaling pathways. We showed that the impact of miRs is orthogonal to pathway regulation at the mRNA level, which stresses the importance of studying posttranscriptional regulation of gene expression. We also report characteristic set of miR and mRNA regulation features linked with bladder cancer.

Molecular pathway activation features linked with transition from normal skin to primary and metastatic melanomas in human.

Melanoma is the most aggressive and dangerous type of skin cancer, but its molecular mechanisms remain largely unclear. For transcriptomic data of 478 primary and metastatic melanoma, nevi and normal skin samples, we performed high-throughput analysis of intracellular molecular networks including 592 signaling and metabolic pathways. We showed that at the molecular pathway level, the formation of nevi largely resembles transition from normal skin to primary melanoma. Using a combination of bioinformatic machine learning algorithms, we identified 44 characteristic signaling and metabolic pathways connected with the formation of nevi, development of primary melanoma, and its metastases. We created a model describing formation and progression of melanoma at the level of molecular pathway activation. We discovered six novel associations between activation of metabolic molecular pathways and progression of melanoma: for allopregnanolone biosynthesis, L-carnitine biosynthesis, zymosterol biosynthesis (inhibited in melanoma), fructose 2, 6-bisphosphate synthesis and dephosphorylation, resolvin D biosynthesis (activated in melanoma), D-myo-inositol hexakisphosphate biosynthesis (activated in primary, inhibited in metastatic melanoma). Finally, we discovered fourteen tightly coordinated functional clusters of molecular pathways. This study helps to decode molecular mechanisms underlying the development of melanoma.

Signaling pathways activation profiles make better markers of cancer than expression of individual genes.

Identification of reliable and accurate molecular markers remains one of the major challenges of contemporary biomedicine. We developed a new bioinformatic technique termed OncoFinder that for the first time enables to quantatively measure activation of intracellular signaling pathways basing on transcriptomic data. Signaling pathways regulate all major cellular events in health and disease. Here, we showed that the Pathway Activation Strength (PAS) value itself may serve as the biomarker for cancer, and compared it with the "traditional" molecular markers based on the expression of individual genes. We applied OncoFinder to profile gene expression datasets for the nine human cancer types including bladder cancer, basal cell carcinoma, glioblastoma, hepatocellular carcinoma, lung adenocarcinoma, oral tongue squamous cell carcinoma, primary melanoma, prostate cancer and renal cancer, totally 292 cancer and 128 normal tissue samples taken from the Gene expression omnibus (GEO) repository. We profiled activation of 82 signaling pathways that involve ~2700 gene products. For 9/9 of the cancer types tested, the PAS values showed better area-under-the-curve (AUC) scores compared to the individual genes enclosing each of the pathways. These results evidence that the PAS values can be used as a new type of cancer biomarkers, superior to the traditional gene expression biomarkers.

Novel robust biomarkers for human bladder cancer based on activation of intracellular signaling pathways.

We recently proposed a new bioinformatic algorithm called OncoFinder for quantifying the activation of intracellular signaling pathways. It was proved advantageous for minimizing errors of high-throughput gene expression analyses and showed strong potential for identifying new biomarkers. Here, for the first time, we applied OncoFinder for normal and cancerous tissues of the human bladder to identify biomarkers of bladder cancer. Using Illumina HT12v4 microarrays, we profiled gene expression in 17 cancer and seven non-cancerous bladder tissue samples. These experiments were done in two independent laboratories located in Russia and Canada. We calculated pathway activation strength values for the investigated transcriptomes and identified signaling pathways that were regulated differently in bladder cancer (BC) tissues compared with normal controls. We found, for both experimental datasets, 44 signaling pathways that serve as excellent new biomarkers of BC, supported by high area under the curve (AUC) values. We conclude that the OncoFinder approach is highly efficient in finding new biomarkers for cancer. These markers are mathematical functions involving multiple gene products, which distinguishes them from "traditional" expression biomarkers that only assess concentrations of single genes.

Oncofinder, a new method for the analysis of intracellular signaling pathway activation using transcriptomic data.

We propose a new biomathematical method, OncoFinder, for both quantitative and qualitative analysis of the intracellular signaling pathway activation (SPA). This method is universal and may be used for the analysis of any physiological, stress, malignancy and other perturbed conditions at the molecular level. In contrast to the other existing techniques for aggregation and generalization of the gene expression data for individual samples, we suggest to distinguish the positive/activator and negative/repressor role of every gene product in each pathway. We show that the relative importance of each gene product in a pathway can be assessed using kinetic models for "low-level" protein interactions. Although the importance factors for the pathway members cannot be so far established for most of the signaling pathways due to the lack of the required experimental data, we showed that ignoring these factors can be sometimes acceptable and that the simplified formula for SPA evaluation may be applied for many cases. We hope that due to its universal applicability, the method OncoFinder will be widely used by the researcher community.