A Comprehensive Machine Learning Benchmark Study for Radiomics-Based Survival Analysis of CT Imaging Data in Patients With Hepatic Metastases of CRC.

OBJECTIVES: Optimizing a machine learning (ML) pipeline for radiomics analysis involves numerous choices in data set composition, preprocessing, and model selection. Objective identification of the optimal setup is complicated by correlated features, interdependency structures, and a multitude of available ML algorithms. Therefore, we present a radiomics-based benchmarking framework to optimize a comprehensive ML pipeline for the prediction of overall survival. This study is conducted on an image set of patients with hepatic metastases of colorectal cancer, for which radiomics features of the whole liver and of metastases from computed tomography images were calculated. A mixed model approach was used to find the optimal pipeline configuration and to identify the added prognostic value of radiomics features. MATERIALS AND METHODS: In this study, a large-scale ML benchmark pipeline consisting of preprocessing, feature selection, dimensionality reduction, hyperparameter optimization, and training of different models was developed for radiomics-based survival analysis. Portal-venous computed tomography imaging data from a previous prospective randomized trial evaluating radioembolization of liver metastases of colorectal cancer were quantitatively accessible through a radiomics approach. One thousand two hundred eighteen radiomics features of hepatic metastases and the whole liver were calculated, and 19 clinical parameters (age, sex, laboratory values, and treatment) were available for each patient. Three ML algorithms-a regression model with elastic net regularization (glmnet), a random survival forest (RSF), and a gradient tree-boosting technique (xgboost)-were evaluated for 5 combinations of clinical data, tumor radiomics, and whole-liver features. Hyperparameter optimization and model evaluation were optimized toward the performance metric integrated Brier score via nested cross-validation. To address dependency structures in the benchmark setup, a mixed-model approach was developed to compare ML and data configurations and to identify the best-performing model. RESULTS: Within our radiomics-based benchmark experiment, 60 ML pipeline variations were evaluated on clinical data and radiomics features from 491 patients. Descriptive analysis of the benchmark results showed a preference for RSF-based pipelines, especially for the combination of clinical data with radiomics features. This observation was supported by the quantitative analysis via a linear mixed model approach, computed to differentiate the effect of data sets and pipeline configurations on the resulting performance. This revealed the RSF pipelines to consistently perform similar or better than glmnet and xgboost. Further, for the RSF, there was no significantly better-performing pipeline composition regarding the sort of preprocessing or hyperparameter optimization. CONCLUSIONS: Our study introduces a benchmark framework for radiomics-based survival analysis, aimed at identifying the optimal settings with respect to different radiomics data sources and various ML pipeline variations, including preprocessing techniques and learning algorithms. A suitable analysis tool for the benchmark results is provided via a mixed model approach, which showed for our study on patients with intrahepatic liver metastases, that radiomics features captured the patients' clinical situation in a manner comparable to the provided information solely from clinical parameters. However, we did not observe a relevant additional prognostic value obtained by these radiomics features.

CXCR4 expression is elevated in TNBC patient derived samples and Z-guggulsterone abrogates tumor progression by targeting CXCL12/CXCR4 signaling axis in preclinical breast cancer model.

Environmental factors such as exposure to ionizing radiations, certain environmental pollutants, and toxic chemicals are considered as risk factors in the development of breast cancer. Triple-negative breast cancer (TNBC) is a molecular variant of breast cancer that lacks therapeutic targets such as progesterone receptor, estrogen receptor, and human epidermal growth factor receptor-2 which makes the targeted therapy ineffective in TNBC patients. Therefore, identification of new therapeutic targets for the treatment of TNBC and the discovery of new therapeutic agents is the need of the hour. In this study, CXCR4 was found to be highly expressed in majority of breast cancer tissues and metastatic lymph nodes derived from TNBC patients. CXCR4 expression is positively correlated with breast cancer metastasis and poor prognosis of TNBC patients suggesting that suppression of CXCR4 expression could be a good strategy in the treatment of TNBC patients. Therefore, the effect of Z-guggulsterone (ZGA) on the expression of CXCR4 in TNBC cells was examined. ZGA downregulated protein and mRNA expression of CXCR4 in TNBC cells and proteasome inhibition or lysosomal stabilization had no effect on the ZGA-induced CXCR4 reduction. CXCR4 is under the transcriptional control of NF-κB, whereas ZGA was found to downregulate transcriptional activity of NF-κB. Functionally, ZGA downmodulated the CXCL12-driven migration/invasion in TNBC cells. Additionally, the effect of ZGA on growth of tumor was investigated in the orthotopic TNBC mice model. ZGA presented good inhibition of tumor growth and liver/lung metastasis in this model. Western blotting and immunohistochemical analysis indicated a reduction of CXCR4, NF-κB, and Ki67 in tumor tissues. Computational analysis suggested PXR agonism and FXR antagonism as targets of ZGA. In conclusion, CXCR4 was found to be overexpressed in majority of patient-derived TNBC tissues and ZGA abrogated the growth of TNBC tumors by partly targeting the CXCL12/CXCR4 signaling axis.

Explaining Blood-Brain Barrier Permeability of Small Molecules by Integrated Analysis of Different Transport Mechanisms.

The blood-brain barrier (BBB) represents a major obstacle to delivering drugs to the central nervous system (CNS), resulting in the lack of effective treatment for many CNS diseases including brain cancer. To accelerate CNS drug development, computational prediction models could save the time and effort needed for experimental evaluation. Here, we studied BBB permeability focusing on active transport (influx and efflux) as well as passive diffusion using previously published and self-curated data sets. We created prediction models based on physicochemical properties, molecular substructures, or their combination to understand which mechanisms contribute to BBB permeability. Our results show that features that predicted passive diffusion over membranes overlap with features that explain endothelial permeation of approved CNS-active drugs. We also identified physical properties and molecular substructures that positively or negatively predicted BBB transport. These findings provide guidance toward identifying BBB-permeable compounds by optimally matching physicochemical and molecular properties to BBB transport mechanisms.

In silico prediction and biological assessment of novel angiogenesis modulators from traditional Chinese medicine.

Uncontrolled angiogenesis is a common denominator underlying many deadly and debilitating diseases such as myocardial infarction, chronic wounds, cancer, and age-related macular degeneration. As the current range of FDA-approved angiogenesis-based medicines are far from meeting clinical demands, the vast reserve of natural products from traditional Chinese medicine (TCM) offers an alternative source for developing pro-angiogenic or anti-angiogenic modulators. Here, we investigated 100 traditional Chinese medicine-derived individual metabolites which had reported gene expression in MCF7 cell lines in the Gene Expression Omnibus (GSE85871). We extracted literature angiogenic activities for 51 individual metabolites, and subsequently analysed their predicted targets and differentially expressed genes to understand their mechanisms of action. The angiogenesis phenotype was used to generate decision trees for rationalising the poly-pharmacology of known angiogenesis modulators such as ferulic acid and curculigoside and validated by an in vitro endothelial tube formation assay and a zebrafish model of angiogenesis. Moreover, using an in silico model we prospectively examined the angiogenesis-modulating activities of the remaining 49 individual metabolites. In vitro, tetrahydropalmatine and 1 beta-hydroxyalantolactone stimulated, while cinobufotalin and isoalantolactone inhibited endothelial tube formation. In vivo, ginsenosides Rb3 and Rc, 1 beta-hydroxyalantolactone and surprisingly cinobufotalin, restored angiogenesis against PTK787-induced impairment in zebrafish. In the absence of PTK787, deoxycholic acid and ursodeoxycholic acid did not affect angiogenesis. Despite some limitations, these results suggest further refinements of in silico prediction combined with biological assessment will be a valuable platform for accelerating the research and development of natural products from traditional Chinese medicine and understanding their mechanisms of action, and also for other traditional medicines for the prevention and treatment of angiogenic diseases.

EMT-Related Genes Have No Prognostic Relevance in Metastatic Colorectal Cancer as Opposed to Stage II/III: Analysis of the Randomised, Phase III Trial FIRE-3 (AIO KRK 0306; FIRE-3).

INTRODUCTION: There is no standard treatment after resection of colorectal liver metastases and the role of systemic therapy remains controversial. To avoid over- or undertreatment, proper risk stratification with regard to postoperative treatment strategy is highly needed. We recently demonstrated the prognostic relevance of EMT-related (epithelial-mesenchymal transition) genes in stage II/III CRC. As EMT is a major step in CRC progression, we now aimed to analyse the prognostic relevance of EMT-related genes in stage IV CRC using the study cohort of the FIRE-3 trial, an open-label multi-centre randomised controlled phase III trial of patients with metastatic CRC. METHODS: Overall and progression free survival were considered as endpoints (n = 350). To investigate the prognostic relevance of EMT-related genes on either endpoint, we compared predictive performance of different models using clinical data only to models using gene data in addition to clinical data, expecting better predictive performance if EMT-related genes have prognostic value. In addition to baseline models (Kaplan Meier (KM), (regularised) Cox), Random Survival Forest (RSF), and gradient boosted trees (GBT) were fit to the data. Repeated, nested five-fold cross-validation was used for hyperparameter optimisation and performance evaluation. Predictive performance was measured by the integrated Brier score (IBS). RESULTS: The baseline KM model showed the best performance (OS: 0.250, PFS: 0.251). None of the other models were able to outperform the KM when using clinical data only according to the IBS scores (OS: 0.253 (Cox), 0.256 (RSF), 0.284 (GBT); PFS: 0.254 (Cox), 0.256 (RSF), 0.276 (GBT)). When adding gene data, performance of GBT improved slightly (OS: 0.262 vs. 0.284; PFS: 0.268 vs. 0.276), however, none of the models performed better than the KM baseline. CONCLUSION: Overall, the results suggest that the prognostic relevance of EMT-related genes may be stage-dependent and that EMT-related genes have no prognostic relevance in stage IV CRC.

Monomerization of Homodimeric Trefoil Factor 3 (TFF3) by an Aminonitrile Compound Inhibits TFF3-Dependent Cancer Cell Survival.

Trefoil factor 3 (TFF3) is a secreted protein with an established oncogenic function and a highly significant association with clinical progression of various human malignancies. Herein, a novel small molecule that specifically targets TFF3 homodimeric functions was identified. Utilizing the concept of reversible covalent interaction, 2-amino-4-(4-(6-fluoro-5-methylpyridin-3-yl)phenyl)-5-oxo-4H,5H-pyrano[3,2-c]chromene-3-carbonitrile (AMPC) was identified as a molecule that interacted with TFF3. AMPC monomerized the cellular and secreted TFF3 homodimer at the cysteine (Cys)57-Cys57 residue with subsequent more rapid degradation of the generated TFF3 monomers. Hence, AMPC treatment also resulted in cellular depletion of TFF3 with consequent decreased cell viability in various human carcinoma-derived TFF3 expressing cell lines, including estrogen receptor positive (ER+) mammary carcinoma (MC). AMPC treatment of TFF3 expressing ER+ MC cells significantly suppressed total cell number in a dose-dependent manner. Consistently, exposure of TFF3 expressing ER+ MC cells to AMPC decreased soft agar colony formation, foci formation, and growth in suspension culture and inhibited growth of preformed colonies in 3D Matrigel. AMPC increased apoptosis in TFF3 expressing ER+ MC cells associated with decreased activity of EGFR, p38, STAT3, AKT, and ERK, decreased protein levels of CCND1, CCNE1, BCL2, and BCL-XL, and increased protein levels of TP53, CDKN1A, CASP7, and CASP9. siRNA-mediated depletion of TFF3 expression in ER+ MC cells efficiently abrogated AMPC-stimulated loss of cell viability and CASPASE 3/7 activities. Furthermore, in mice bearing ER+ MC cell-generated xenografts, AMPC treatment significantly impeded xenograft growth. Hence, AMPC exemplifies a novel mechanism by which small molecule drugs may inhibit a dimeric oncogenic protein and provides a strategy to impede TFF3-dependent cancer progression.

Targeting Cell Death Mechanism Specifically in Triple Negative Breast Cancer Cell Lines.

Triple negative breast cancer (TNBC) is currently associated with a lack of treatment options. Arsenic derivatives have shown antitumoral activity both in vitro and in vivo; however, their mode of action is not completely understood. In this work we evaluate the response to arsenate of the double positive MCF-7 breast cancer cell line as well as of two different TNBC cell lines, Hs578T and MDA-MB-231. Multimodal experiments were conducted to this end, using functional assays and microarrays. Arsenate was found to induce cytoskeletal alteration, autophagy and apoptosis in TNBC cells, and moderate effects in MCF-7 cells. Gene expression analysis showed that the TNBC cell lines' response to arsenate was more prominent in the G2M checkpoint, autophagy and apoptosis compared to the Human Mammary Epithelial Cells (HMEC) and MCF-7 cell lines. We confirmed the downregulation of anti-apoptotic genes (MCL1, BCL2, TGFß1 and CCND1) by qRT-PCR, and on the protein level, for TGFß2, by ELISA. Insight into the mode of action of arsenate in TNBC cell lines it is provided, and we concluded that TNBC and non-TNBC cell lines reacted differently to arsenate treatment in this particular experimental setup. We suggest the future research of arsenate as a treatment strategy against TNBC.

Latent Variables Capture Pathway-Level Points of Departure in High-Throughput Toxicogenomic Data.

Estimation of points of departure (PoDs) from high-throughput transcriptomic data (HTTr) represents a key step in the development of next-generation risk assessment (NGRA). Current approaches mainly rely on single key gene targets, which are constrained by the information currently available in the knowledge base and make interpretation challenging as scientists need to interpret PoDs for thousands of genes or hundreds of pathways. In this work, we aimed to address these issues by developing a computational workflow to investigate the pathway concentration-response relationships in a way that is not fully constrained by known biology and also facilitates interpretation. We employed the Pathway-Level Information ExtractoR (PLIER) to identify latent variables (LVs) describing biological activity and then investigated in vitro LVs' concentration-response relationships using the ToxCast pipeline. We applied this methodology to a published transcriptomic concentration-response data set for 44 chemicals in MCF-7 cells and showed that our workflow can capture known biological activity and discriminate between estrogenic and antiestrogenic compounds as well as activity not aligning with the existing knowledge base, which may be relevant in a risk assessment scenario. Moreover, we were able to identify the known estrogen activity in compounds that are not well-established ER agonists/antagonists supporting the use of the workflow in read-across. Next, we transferred its application to chemical compounds tested in HepG2, HepaRG, and MCF-7 cells and showed that PoD estimates are in strong agreement with those estimated using a recently developed Bayesian approach (cor = 0.89) and in weak agreement with those estimated using a well-established approach such as BMDExpress2 (cor = 0.57). These results demonstrate the effectiveness of using PLIER in a concentration-response scenario to investigate pathway activity in a way that is not fully constrained by the knowledge base and to ease the biological interpretation and support the development of an NGRA framework with the ability to improve current risk assessment strategies for chemicals using new approach methodologies.

Prediction and identification of synergistic compound combinations against pancreatic cancer cells.

Resistance to current therapies is common for pancreatic cancer and hence novel treatment options are urgently needed. In this work, we developed and validated a computational method to select synergistic compound combinations based on transcriptomic profiles from both the disease and compound side, combined with a pathway scoring system, which was then validated prospectively by testing 30 compounds (and their combinations) on PANC-1 cells. Some compounds selected as single agents showed lower GI50 values than the standard of care, gemcitabine. Compounds suggested as combination agents with standard therapy gemcitabine based on the best performing scoring system showed on average 2.82-5.18 times higher synergies compared to compounds that were predicted to be active as single agents. Examples of highly synergistic in vitro validated compound pairs include gemcitabine combined with Entinostat, thioridazine, loperamide, scriptaid and Saracatinib. Hence, the computational approach presented here was able to identify synergistic compound combinations against pancreatic cancer cells.

Bi-Centric Independent Validation of Outcome Prediction after Radioembolization of Primary and Secondary Liver Cancer.

BACKGROUND: Yttrium-90 radioembolization (RE) plays an important role in the treatment of liver malignancies. Optimal patient selection is crucial for an effective and safe treatment. In this study, we aim to validate the prognostic performance of a previously established random survival forest (RSF) with an external validation cohort from a different national center. Furthermore, we compare outcome prediction models with different established metrics. METHODS: A previously established RSF model, trained on a consecutive cohort of 366 patients who had received RE due to primary or secondary liver tumor at a national center (center 1), was used to predict the outcome of an independent consecutive cohort of 202 patients from a different national center (center 2) and vice versa. Prognostic performance was evaluated using the concordance index (C-index) and the integrated Brier score (IBS). The prognostic importance of designated baseline parameters was measured with the minimal depth concept, and the influence on the predicted outcome was analyzed with accumulated local effects plots. RSF values were compared to conventional cox proportional hazards models in terms of C-index and IBS. RESULTS: The established RSF model achieved a C-index of 0.67 for center 2, comparable to the results obtained for center 1, which it was trained on (0.66). The RSF model trained on center 2 achieved a C-index of 0.68 on center 2 data and 0.66 on center 1 data. CPH models showed comparable results on both cohorts, with C-index ranging from 0.68 to 0.72. IBS validation showed more differentiated results depending on which cohort was trained on and which cohort was predicted (range: 0.08 to 0.20). Baseline cholinesterase was the most important variable for survival prediction. CONCLUSION: The previously developed predictive RSF model was successfully validated with an independent external cohort. C-index and IBS are suitable metrics to compare outcome prediction models, with IBS showing more differentiated results. The findings corroborate that survival after RE is critically determined by functional hepatic reserve and thus baseline liver function should play a key role in patient selection.

Transcriptional drug repositioning and cheminformatics approach for differentiation therapy of leukaemia cells.

Differentiation therapy is attracting increasing interest in cancer as it can be more specific than conventional chemotherapy approaches, and it has offered new treatment options for some cancer types, such as treating acute promyelocytic leukaemia (APL) by retinoic acid. However, there is a pressing need to identify additional molecules which act in this way, both in leukaemia and other cancer types. In this work, we hence developed a novel transcriptional drug repositioning approach, based on both bioinformatics and cheminformatics components, that enables selecting such compounds in a more informed manner. We have validated the approach for leukaemia cells, and retrospectively retinoic acid was successfully identified using our method. Prospectively, the anti-parasitic compound fenbendazole was tested in leukaemia cells, and we were able to show that it can induce the differentiation of leukaemia cells to granulocytes in low concentrations of 0.1 µM and within as short a time period as 3 days. This work hence provides a systematic and validated approach for identifying small molecules for differentiation therapy in cancer.

Structure-based identification of dual ligands at the A2AR and PDE10A with anti-proliferative effects in lung cancer cell-lines.

Enhanced/prolonged cAMP signalling has been suggested as a suppressor of cancer proliferation. Interestingly, two key modulators that elevate cAMP, the A2A receptor (A2AR) and phosphodiesterase 10A (PDE10A), are differentially co-expressed in various types of non-small lung cancer (NSCLC) cell-lines. Thus, finding dual-target compounds, which are simultaneously agonists at the A2AR whilst also inhibiting PDE10A, could be a novel anti-proliferative approach. Using ligand- and structure-based modelling combined with MD simulations (which identified Val84 displacement as a novel conformational descriptor of A2AR activation), a series of known PDE10A inhibitors were shown to dock to the orthosteric site of the A2AR. Subsequent in-vitro analysis confirmed that these compounds bind to the A2AR and exhibit dual-activity at both the A2AR and PDE10A. Furthermore, many of the compounds exhibited promising anti-proliferative effects upon NSCLC cell-lines, which directly correlated with the expression of both PDE10A and the A2AR. Thus, we propose a structure-based methodology, which has been validated in in-vitro binding and functional assays, and demonstrated a promising therapeutic value.

Comparison of Cellular Morphological Descriptors and Molecular Fingerprints for the Prediction of Cytotoxicity- and Proliferation-Related Assays.

Cell morphology features, such as those from the Cell Painting assay, can be generated at relatively low costs and represent versatile biological descriptors of a system and thereby compound response. In this study, we explored cell morphology descriptors and molecular fingerprints, separately and in combination, for the prediction of cytotoxicity- and proliferation-related in vitro assay endpoints. We selected 135 compounds from the MoleculeNet ToxCast benchmark data set which were annotated with Cell Painting readouts, where the relatively small size of the data set is due to the overlap of required annotations. We trained Random Forest classification models using nested cross-validation and Cell Painting descriptors, Morgan and ErG fingerprints, and their combinations. While using leave-one-cluster-out cross-validation (with clusters based on physicochemical descriptors), models using Cell Painting descriptors achieved higher average performance over all assays (Balanced Accuracy of 0.65, Matthews Correlation Coefficient of 0.28, and AUC-ROC of 0.71) compared to models using ErG fingerprints (BA 0.55, MCC 0.09, and AUC-ROC 0.60) and Morgan fingerprints alone (BA 0.54, MCC 0.06, and AUC-ROC 0.56). While using random shuffle splits, the combination of Cell Painting descriptors with ErG and Morgan fingerprints further improved balanced accuracy on average by 8.9% (in 9 out of 12 assays) and 23.4% (in 8 out of 12 assays) compared to using only ErG and Morgan fingerprints, respectively. Regarding feature importance, Cell Painting descriptors related to nuclei texture, granularity of cells, and cytoplasm as well as cell neighbors and radial distributions were identified to be most contributing, which is plausible given the endpoint considered. We conclude that cell morphological descriptors contain complementary information to molecular fingerprints which can be used to improve the performance of predictive cytotoxicity models, in particular in areas of novel structural space.

Towards understanding antimicrobial activity, cytotoxicity and the mode of action of dichapetalins A and M using in silico and in vitro studies.

Dichapetalum madagascariense Poir (Dichapetalaceae) is traditionally used to treat bacterial infections, jaundice, urethritis and viral hepatitis in Africa. Its root contains a broad spectrum of biologically active dichapetalins. To evaluate the plant's effect on human MCF-7 cells and its' antibacterial and antiparasitic potentials, we isolated and identified the known dichapetalins A and M from the roots. Both dichapetalins were tested on six bacterial strains (Shigella flexneri, Bacillus cereus, Salmonella paratyphi B, Listeria monocytogenes, Escherichia coli, Staphylococcus aureus) and two parasite strains; Trypanosoma brucei brucei, and Leishmania donovani using the Alamar Blue assay system. Dichapetalins A and M were more potent against B. cereus with IC50 values of 11.15 and 3.15 µg/ml, respectively, compared to the positive control ampicillin (IC50 = 19.50 µg/ml). Dichapetalins A (IC50 = 74.22 µg/ml) and M (IC50 = 72.34 µg/ml) were less active against T. b. brucei, compared to the standard Suramin (IC50 = 4.96 µg/ml). Dichapetalin M showed moderate activity against L. donovani (Amphotericin B: IC50 = 0.21 µg/ml) with an IC50 of 16.80 µg/ml. In human MCF-7 cells expressing the NR1I2 receptor, the activity of dichapetalin M was higher (IC50 = 4.71 µM and 3.95 µM) for 48 and 72 h of treatment, respectively compared to Curcumin with IC50 of 17.49 µM and 12.53 µM for 48 and 72 h of treatment, respectively. Results from in vitro expression studies with qPCR confirmed an antagonistic effect of dichapetalin M on PXR (NR1I2) signaling; supporting the PXR signaling pathway as a possible mode of action of dichapetalin M as predicted by in silico results. These findings confirm previous studies that D. madagascariense can be a source of potential lead compounds for development of novel antibiotic, antiparasitic and anticancer medicines, and provide further insights into the mechanism of action of the dichapetalins.

Identification of Intrinsic Drug Resistance and Its Biomarkers in High-Throughput Pharmacogenomic and CRISPR Screens.

High-throughput drug screens in cancer cell lines test compounds at low concentrations, thereby enabling the identification of drug-sensitivity biomarkers, while resistance biomarkers remain underexplored. Dissecting meaningful drug responses at high concentrations is challenging due to cytotoxicity, i.e., off-target effects, thus limiting resistance biomarker discovery to frequently mutated cancer genes. To address this, we interrogate subpopulations carrying sensitivity biomarkers and consecutively investigate unexpectedly resistant (UNRES) cell lines for unique genetic alterations that may drive resistance. By analyzing the GDSC and CTRP datasets, we find 53 and 35 UNRES cases, respectively. For 24 and 28 of them, we highlight putative resistance biomarkers. We find clinically relevant cases such as EGFRT790M mutation in NCI-H1975 or PTEN loss in NCI-H1650 cells, in lung adenocarcinoma treated with EGFR inhibitors. Interrogating the underpinnings of drug resistance with publicly available CRISPR phenotypic assays assists in prioritizing resistance drivers, offering hypotheses for drug combinations.

Transcriptomics predicts compound synergy in drug and natural product treated glioblastoma cells.

Pathway analysis is an informative method for comparing and contrasting drug-induced gene expression in cellular systems. Here, we define the effects of the marine natural product fucoxanthin, separately and in combination with the prototypic phosphatidylinositol 3-kinase (PI3K) inhibitor LY-294002, on gene expression in a well-established human glioblastoma cell system, U87MG. Under conditions which inhibit cell proliferation, LY-294002 and fucoxanthin modulate many pathways in common, including the retinoblastoma, DNA damage, DNA replication and cell cycle pathways. In sharp contrast, we see profound differences in the expression of genes characteristic of pathways such as apoptosis and lipid metabolism, contributing to the development of a differentiated and distinctive drug-induced gene expression signature for each compound. Furthermore, in combination, fucoxanthin synergizes with LY-294002 in inhibiting the growth of U87MG cells, suggesting complementarity in their molecular modes of action and pointing to further treatment combinations. The synergy we observe between the dietary nutraceutical fucoxanthin and the synthetic chemical LY-294002 in producing growth arrest in glioblastoma, illustrates the potential of nutri-pharmaceutical combinations in targeting this challenging disease.

QSAR-derived affinity fingerprints (part 2): modeling performance for potency prediction.

Affinity fingerprints report the activity of small molecules across a set of assays, and thus permit to gather information about the bioactivities of structurally dissimilar compounds, where models based on chemical structure alone are often limited, and model complex biological endpoints, such as human toxicity and in vitro cancer cell line sensitivity. Here, we propose to model in vitro compound activity using computationally predicted bioactivity profiles as compound descriptors. To this aim, we apply and validate a framework for the calculation of QSAR-derived affinity fingerprints (QAFFP) using a set of 1360 QSAR models generated using Ki, Kd, IC50 and EC50 data from ChEMBL database. QAFFP thus represent a method to encode and relate compounds on the basis of their similarity in bioactivity space. To benchmark the predictive power of QAFFP we assembled IC50 data from ChEMBL database for 18 diverse cancer cell lines widely used in preclinical drug discovery, and 25 diverse protein target data sets. This study complements part 1 where the performance of QAFFP in similarity searching, scaffold hopping, and bioactivity classification is evaluated. Despite being inherently noisy, we show that using QAFFP as descriptors leads to errors in prediction on the test set in the ~ 0.65-0.95 pIC50 units range, which are comparable to the estimated uncertainty of bioactivity data in ChEMBL (0.76-1.00 pIC50 units). We find that the predictive power of QAFFP is slightly worse than that of Morgan2 fingerprints and 1D and 2D physicochemical descriptors, with an effect size in the 0.02-0.08 pIC50 units range. Including QSAR models with low predictive power in the generation of QAFFP does not lead to improved predictive power. Given that the QSAR models we used to compute the QAFFP were selected on the basis of data availability alone, we anticipate better modeling results for QAFFP generated using more diverse and biologically meaningful targets. Data sets and Python code are publicly available at https://github.com/isidroc/QAFFP_regression .

Applying synergy metrics to combination screening data: agreements, disagreements and pitfalls.

Synergistic drug combinations are commonly sought to overcome monotherapy resistance in cancer treatment. To identify such combinations, high-throughput cancer cell line combination screens are performed; and synergy is quantified using competing models based on fundamentally different assumptions. Here, we compare the behaviour of four synergy models, namely Loewe additivity, Bliss independence, highest single agent and zero interaction potency, using the Merck oncology combination screen. We evaluate agreements and disagreements between the models and investigate putative artefacts of each model's assumptions. Despite at least moderate concordance between scores (Pearson's r >0.32, Spearman's ρ>0.34), multiple instances of strong disagreement were observed. Those disagreements are driven by, among others, large differences in tested concentrations, maximum response values and median effective concentrations.

Elucidating Compound Mechanism of Action and Predicting Cytotoxicity Using Machine Learning Approaches, Taking Prediction Confidence into Account.

The modes of action (MoAs) of drugs frequently are unknown, because many are small molecules initially identified from phenotypic screens, giving rise to the need to elucidate their MoAs. In addition, the high attrition rate for candidate drugs in preclinical studies due to intolerable toxicity has motivated the development of computational approaches to predict drug candidate (cyto)toxicity as early as possible in the drug-discovery process. Here, we provide detailed instructions for capitalizing on bioactivity predictions to elucidate the MoAs of small molecules and infer their underlying phenotypic effects. We illustrate how these predictions can be used to infer the underlying antidepressive effects of marketed drugs. We also provide the necessary functionalities to model cytotoxicity data using single and ensemble machine-learning algorithms. Finally, we give detailed instructions on how to calculate confidence intervals for individual predictions using the conformal prediction framework. © 2019 by John Wiley & Sons, Inc.

Triazole-Pyridine Dicarbonitrile Targets Phosphodiesterase 4 to Induce Cytotoxicity in Lung Carcinoma Cells.

Phosphodiesterase 4 (PDE4) is a key enzyme involved in the hydrolysis of cyclic adenosine monophosphate (cAMP) and widely expressed in several types of cancers. The inhibition of PDE4 results in an increased concentration of intracellular cAMP levels that imparts the anti-inflammatory response in the target cells. In the present report, two series of triazolo-pyridine dicarbonitriles and substituted dihydropyridine dicarbonitriles were synthesized using green protocol (TBAB in refluxed water). We next evaluated the title compounds for their cytotoxicity towards lung cancer (A549) cells and identified 7'-[4-(methylsulfonyl)phenyl]-5'-oxo-1',5'-dihydrospiro[cyclohexane-1,2'-[1,2,4]triazolo[1,5-a]pyridine]-6',8'-dicarbonitrile (5h) and 7'-(1-methyl-1H-imidazol-2-yl)-5'-oxo-1',5'-dihydrospiro[cyclohexane-1,2'-[1,2,4]triazolo[1,5-a]pyridine]-6',8'-dicarbonitrile (5j) as lead analogs with the IC50 values of 15.2 and 24.1 µm, respectively. Furthermore, all the new compounds were tested for PDE4 inhibitory activity and 5j showed relatively good inhibitory activity towards PDE4 with inhibition of 50.9 % at 10 µm. In silico analysis demonstrated the favorable interaction of the title compounds with the target enzyme. Taken together, the present study introduces a new scaffold for the development of novel PDE4 inhibitors to fight against inflammatory diseases.