Next Article in Journal
MdPP2C24/37, Protein Phosphatase Type 2Cs from Apple, Interact with MdPYL2/12 to Negatively Regulate ABA Signaling in Transgenic Arabidopsis
Next Article in Special Issue
hgtseq: A Standard Pipeline to Study Horizontal Gene Transfer
Previous Article in Journal
Eukaryotic CRFK Cells Motion Characterized with Atomic Force Microscopy
Previous Article in Special Issue
What Is a Digital Twin? Experimental Design for a Data-Centric Machine Learning Perspective in Health
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antiproliferative Activity Predictor: A New Reliable In Silico Tool for Drug Response Prediction against NCI60 Panel

by
Annamaria Martorana
1,†,
Gabriele La Monica
1,†,
Alessia Bono
1,
Salvatore Mannino
1,
Silvestre Buscemi
1,
Antonio Palumbo Piccionello
1,
Carla Gentile
1,
Antonino Lauria
1,* and
Daniele Peri
2
1
Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche “STEBICEF”, Università Degli Studi di Palermo, Viale Delle Scienze Ed. 17, I-90128 Palermo, Italy
2
Dipartimento di Ingegneria Dell’innovazione Industriale e Digitale, Università Degli Studi di Palermo, Viale Delle Scienze Ed. 6, I-90128 Palermo, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(22), 14374; https://doi.org/10.3390/ijms232214374
Submission received: 22 September 2022 / Revised: 13 November 2022 / Accepted: 16 November 2022 / Published: 19 November 2022
(This article belongs to the Special Issue Data Mining and Bioinformatic Tools for Health)

Abstract

:
In vitro antiproliferative assays still represent one of the most important tools in the anticancer drug discovery field, especially to gain insights into the mechanisms of action of anticancer small molecules. The NCI-DTP (National Cancer Institute Developmental Therapeutics Program) undoubtedly represents the most famous project aimed at rapidly testing thousands of compounds against multiple tumor cell lines (NCI60). The large amount of biological data stored in the National Cancer Institute (NCI) database and many other databases has led researchers in the fields of computational biology and medicinal chemistry to develop tools to predict the anticancer properties of new agents in advance. In this work, based on the available antiproliferative data collected by the NCI and the manipulation of molecular descriptors, we propose the new in silico Antiproliferative Activity Predictor (AAP) tool to calculate the GI50 values of input structures against the NCI60 panel. This ligand-based protocol, validated by both internal and external sets of structures, has proven to be highly reliable and robust. The obtained GI50 values of a test set of 99 structures present an error of less than ±1 unit. The AAP is more powerful for GI50 calculation in the range of 4–6, showing that the results strictly correlate with the experimental data. The encouraging results were further supported by the examination of an in-house database of curcumin analogues that have already been studied as antiproliferative agents. The AAP tool identified several potentially active compounds, and a subsequent evaluation of a set of molecules selected by the NCI for the one-dose/five-dose antiproliferative assays confirmed the great potential of our protocol for the development of new anticancer small molecules. The integration of the AAP tool in the free web service DRUDIT provides an interesting device for the discovery and/or optimization of anticancer drugs to the medicinal chemistry community. The training set will be updated with new NCI-tested compounds to cover more chemical spaces, activities, and cell lines. Currently, the same protocol is being developed for predicting the TGI (total growth inhibition) and LC50 (median lethal concentration) parameters to estimate toxicity profiles of small molecules.

1. Introduction

In the search for new chemical compounds endowed with anticancer properties, in vitro antiproliferative screening remains one of the most commonly used approaches to identify new biologically active compounds. As demonstrated by the numerous research projects focused on the characterization of tumor cells (including National Cancer Institute Human Tumor Cell Lines Screen, NCI60 [1]; Cancer Cell Line Encyclopedia, CCLE [2]; Genomics of Drug Sensitivity in Cancer, GDSC [3]; The Cancer Genome Atlas Program, TCGA [4]; Cancer Therapeutics Response Portal, CTRP [5]), the study of cancer cell lines has led to gaining insights into the metabolic pathways involved in uncontrolled growth and the mechanisms of action of the tested anticancer compounds [6,7,8].
The NCI60 Human Tumor Cell Lines Screen is the pioneering and best-known project based on the large-scale screening of chemical compounds and cancer cell phenotypes. For more than 20 years, it has been of great benefit to the global research community and has proven to be a rich and reliable source of information. The NCI60 high-throughput screening program, in its current version, allows the screening of up to 3000 compounds per year and consists of a standardized assay performed on approximately sixty cancer cell lines belonging to nine different subpanels (each one related to a specific tumor: leukemia, non-small-cell lung (NSCLC), colon, central nervous system (CNS), melanoma, ovarian, renal, prostate, and breast cancer cells) characterized at the genomic, transcriptomic, and proteomic levels (Molecular Target Characterization Program) [1,9,10,11,12,13].
Considering the high failure rate in anticancer drug discovery research and the significant amount of time and resources required for in vitro screening, in silico predictions, especially in the early stages of the drug development, could represent a critical aid. To this end, several computational protocols and machine learning approaches have been proposed, as recently reviewed by Firoozbakht et al. [14].
In particular, COMPARE Analysis was the first tool that used the NCI database to correlate biological data with chemical entries (compound–compound; compound–target; target–target) permitting researchers to hypothesize the mechanism of action/resistance of a compound and the putative targets involved in its antiproliferative activity [15,16]. The great potential of this tool encouraged the development of other high-performing protocols, such as CellMiner, a web-accessible correlation application exploiting the NCI database [17,18,19]. CellMiner tools enable the rapid data retrieval of transcripts for genes and microRNAs, along with activity reports for thousands of chemical compounds, including several drugs approved by the Food and Drug Administration (FDA).
The elaboration of these data into quantitative patterns against the NCI60 dataset clarified their cross-correlational organization using a novel pattern matching tool. To this end, several case studies of the in silico discovery process afforded by CellMiner have been reported in the literature (i.e., analyses of multidrug resistance and doxorubicin activity; the identification of colon-specific genes, microRNAs, and drugs; etc.).
Recently, Lind et al. developed a regression random forest model by the integration of mutational status data of 145 selected oncogenes from the GDSC database and more than 1200 molecular descriptors to reliably calculate the GI50 values of selected compounds against cancer cells [20]. Similarly, Zhang et al. designed a dual-layer integrated cell–drug network capable of predicting drug responses using the gene expression profiles of numerous cancer cell types available in the CCLE and the comprehensive genomic profiling (CGP), providing the chemical proprieties of drugs captured by molecular descriptors [21]. In addition, numerous web-based tools have been described, such as the PaccMann web service, which can estimate the chemosensitivity of a cancer cell line by integrating both cancer cell and chemical structure features [22], and CDRscan [23], a deep learning model that assesses anticancer drug responsiveness based on large-scale drug screening assay data. It employs a two-step convolution architecture, where the genomic mutational fingerprints of cell lines and the molecular fingerprints of drugs are processed individually then merged by ‘virtual docking’. An analysis of the extrapolation capability revealed a high accuracy (R2 > 0.84; AUROC > 0.98). The tool was applied to a large set of approved drugs, allowing the identification of 14 oncology and 23 non-oncology drugs. DeepIC50 [24], a one-dimensional convolution neural network model designed to predict drug responsiveness classes, was based on a large set of features. As a result, it showed better cell viability (IC50) prediction accuracy in pancancer cell lines over two independent cancer cell line datasets. More recently, pdCSM-cancer, which uses a graph-based signature representation, has been used to estimate the antiproliferative activity against multiple cancer cell lines [25].
In this work, considering our expertise in using molecular descriptors for biological purposes [26,27,28,29,30], we propose a new in silico Antiproliferative Activity Predictor (AAP) tool that can predict the antiproliferative activity of compounds against the NCI60 panel. The protocol, which is freely available on the DRUDIT web service (https://www.drudit.com (accessed on 15 November 2022) [28]), was based on both the molecular descriptors and in vitro antiproliferative data of tested compounds belonging to the NCI database.
In this article, the computational functions and the validation of the AAP tool are discussed in detail. Since our research group has focused on the synthesis of new small molecules with anticancer properties for many years [31,32,33], we decided to evaluate the performance of the AAP protocol by predicting the anticancer potential of an in-house small molecule database. The anticancer properties of several curcumin-like derivatives, some of which had already been evaluated as neuroprotective agents, were investigated [34,35,36]. To validate and confirm the AAP in silico data, in vitro antiproliferative assays of selected compounds were performed as a part of the NCI60 DTP screening program.

2. Results and Discussion

2.1. NCI60 Antiproliferative Activity Predictor Tool

The Antiproliferative Activity Predictor (AAP) tool is a new molecular-descriptor-based protocol for predicting the anticancer activities, expressed as GI50, of small molecules against the NCI60 panel. The AAP has been included as a module in DRUDIT, a web service that has already proven to be a reliable and valuable support for the development of new small molecules with biological activities [28,37].

2.1.1. Description of the Tool Learning Process

The AAP tool consists of sequential steps, as reported in the flowchart (Figure 1).
First, the NCI60 database, which contains in vitro data on the antiproliferative effects, expressed as GI50, of many compounds, was selected [10] (Figure 2).
From the thousands of structures tested by the NCI (Figure 2A), those that were screened in a five-dose assay were selected by considering the experimental GI50 values (Figure 2C). Then, according to publication dates, they were split into a training set (Figure 2C, data published until September 2014, referred to as NCI2014DB), which was used to build the model, and a test set (Figure 2E, data published until June 2016, referred to as NCI2016DB (Figure 2B)), which was used to validate the tool [38].
The AAP protocol results were obtained from the weighed contributions of two “modules”, called Finger-Print (FP) and Cell-Lines (CL); these contributed through a series of well-considered steps (Figure 1) to assign -logGI50 values (indicated as GI50 in the following sections) to input structures against the NCI60 panel cell lines (Ln).
Preliminarily, a set of molecular descriptors (D1, D2, D3, and Dn) was calculated for the training set (38 k structures, N, from NCI2014DB, Figure 2C). The molecular descriptor calculation was performed by MOLDESTO software (Supporting Information S1 contains a list of molecular descriptors implemented in MOLDESTO; see the Methods Section), which calculates more than 1000 1D, 2D, and 3D molecular descriptors (Di).
Then, using the molecular descriptors matrix (structures versus molecular descriptors) the above-described modules (FP and CL) were built as described below.
The FP module relies on structure similarity. It matches the molecular descriptor sequence of the input structure Di(Xj) (Figure 1) with the Di sequences of the structures belonging to the training set and assigns the S score as follows:
S = N(Di(Xj))/N(Di);
where N(Di(Xj)) is the number of the molecular descriptors, Di(Xj), which has a value in the range Di ± Di × 0.05, and N(Di) is the number of the molecular descriptors used.
By ranking according to the S score, the protocol assigns to the input structure Xj the experimental GI50 values of the best-scored structure of the training set (GI50(FP)). If experimental GI50 data of the training set structure are missing for one or more cell lines, a GI50 value is not assigned.
The CL module, on the other hand, is based on the cell lines. For each of the sixty NCI cell lines, 42 templates (CLi) were built: 40 for the GI50 values in the range of 4–8 (0.1 units each: 4.0–4.1; 4.1–4.2;….;7.9–8.0) and 2 for GI50 values <4 and >8.
With this aim, the structures of the training set were assigned to each template according to the experimental GI50 values. In detail, all the structures that, for the specific cell line CLi, had a GI50 value in the appropriate range were assigned to the related template (Figure 3).
Then, the mean (μi) and standard deviation (σi) were computed for all molecular descriptors, considering the structures belonging to each template.
Thus, for each of the sixty cell lines, the molecular descriptor values of the input structure Xj were matched with the 42 cell line templates. The input structure Xj was assigned the GI50 value (GI50(CL)) of the corresponding template with the higher matching score (Figure 4).
This score was assigned by considering the sequence of Di values of the input structure Xj matched to the sequence of the µi(Di) ± σi(Di) range of the cell line template. If a Di value was in the range, a value of 1 was assigned; otherwise, a value of 0 was assigned. The sum of these binary scores, normalized to all molecular descriptors, gave the CL scores (Figure 4) and consequently the higher value, allowing GI50(CL) assignation (Figure 4).
Once an input structure, Xj, was uploaded into the DRUDIT tools interface, MOLDESTO optimized the geometry in vacuo and calculated the molecular descriptors described above for the training set. Then, the molecular descriptor values were submitted to the FP and CL modules.
The output data from these modules were weighted as shown below:
GI50i = GI50i(FP) × S + GI50i(CL) × (1 − S)
where GI50i is the GI50 value for that cell line, S is the fingerprint score, GI50i(FP) is the GI50 value predicted by the FP module, and GI50i(CL) is the GI50 value assigned by the CL module.
From this formula, if the structural similarity between the input structure and the best-scored structure of the training set was high (S score close to 1), it was assumed that the biological activity of the input structure was very similar to that of the compound from the training set (a similar structure could correspond to a similar biological activity). When S = 1, the input structure was included in the training set. Thus, the predicted GI50 corresponded to the experimental values. Instead, if S was not close to 1, GI50(CL) contributed more to the overall result, according to the S value. When the GI50(FP) was not available (unavailable data from the NCI screening), the GI50 corresponded to the GI50(CL).

2.1.2. Validation of the AAP Tool

The predictive ability of the AAP was validated by internal and external validation. Internal validation: First, 5‰ of the training database structures (193 molecules), randomly selected from NCI2014DB (Figure 2D, Supporting Information S2), were used to validate the CL module by matching the calculated GI50(CL) values with the experimental GI50 data. Because these structures were used to generate the AAP protocol, their experimental GI50 values were indicated by the FP protocol, except for those that were not available (the experimental GI50 values are listed in Supporting Information S3; an empty box in the matrix indicates unavailable experimental data). The 193 structures were clustered into three groups according to their GI50 values: the most active compounds (more than 40/60 GI50 values equal to 8); the structures with GI50 values in the range of 4–8; and the cluster of less active/inactive compounds with GI50 values close to 4.
Therefore, GI50(CL) values were first calculated for the selected structures by setting the DRUDIT parameters (see the Methods Section for the meaning of the DRUDIT parameters). This step had two aims: it allowed us to verify the predictive capability of the CL module and, more importantly, to tune the DRUDIT parameters for the best prediction of antiproliferative activity for new compounds. Thus, runs 1–18 were performed by modulating the values of the parameters N, Z, and G, as reported in Table 1. The 18 outputs for the CL module are reported in Supporting Information S4).
The eighteen matrices from CL were matched with the experimental GI50 values to obtain 18 new matrices in which the |DTV (GI50)|, i.e., the absolute deviation of the calculated GI50(CL) value from the experimental GI50 value, was reported for each structure. Furthermore, for each entry, the average |DTV(GI50)| for all runs was examined.
From the analysis of these data, it appears that the protocol allows the identification of potentially inactive or moderately active structures (below <4 or in the range of 4–7) with a remarkable degree of reliability, while it is less effective for structures with high activity (in the range of 7–8) but with acceptable errors.
Moreover, it was demonstrated that the quality of the prediction was closely related to the amount of available biological data used to build the model. Since the number of highly active compounds (7 < GI50 < 8) was very low with respect to inactive or moderately active compounds, the prediction was negatively affected (higher error). Thus, the M19-MEL melanoma cancer cell line was kept out of the analysis due to a lack of sufficient biological data.
The matrix of GI50 values provided by each run was further elaborated to calculate the overall average |DTV(GI50)| for each cancer cell line. SK-MEL-28, belonging to the melanoma panel, gave the best predictions (|DTV(GI50)| of 1.14).
To select the optimized DRUDIT parameters, |DTV(GI50)| was calculated for each GI50(CL) matrix (Table 1). The parameters of run 1 (N = 240, Z = 50, G = a) were identified as the best, with an overall average error of 1.22 (Table 1). In this run, the renal cancer panel was the best prediction, with an average error of 1.18. The full results are reported in Supporting Information S5.
External test validation: A set of 99 molecules was collected from NCI2016DB (Figure 2B, see Methods Section for database selection and Supporting Information S6). Their known GI50 values were compared with the predicted ones, which were calculated using the optimized DRUDIT parameters (run 1, see above). The output matrix showed an interesting scenario with a |DTV(GI50)| of 0.87 and excellent predictions for structures with low activity (GI50 > 4 for at least 40 cell lines). On the other hand, significant errors were recorded for structures with high experimental GI50 values (GI50 above or close to 8). This evidence confirmed the capability of the protocol to better predict GI50 values for low-activity molecules. The average errors for each cancer cell line were also calculated, and they showed excellent prediction for the breast cancer panel (average error of 0.77 for the panel), especially against MDA-MB-231-ATCC (average error for the cancer line of 0.64).
Analyzing the |DTV(GI50)| for all selected databases, considered by ranges, it was found that the protocol was able to assign the correct value returning ÷1 for 65% of the data (Figure 5).
A matrix that compares the experimental GI50 values with the calculated ones for all structures is reported in Supporting Information S7, in addition to the complete data analysis.
With the aim to demonstrate the relevant contribution of the CL module in the prediction, we also compared the experimental GI50 values with the calculated GI50(FP) values. The total average |DTV(GI50)| of 0.95, which is higher than that obtained by combining the predicted GI50 values given by both modules, shows that the CL module improves accuracy and leads to better predictions (Supporting Information S8).

2.1.3. Parameter Optimization for Cell Line/Subpanel Activity Prediction

Tuning the DRUDIT parameters (N, Z, and G) could also allow the optimization of the prediction for a specific cell line or subpanel (all the following results are shown in Supporting Information S9).
With this aim, the analysis was addressed to determinate the best combination of parameters (1–18) for the nine NCI subpanels and for each cancer cell line, as described above.
Then, the average |DTV(GI50)| values obtained for the cancer cell lines in all runs were analyzed (Supporting Information S5).
Regarding the optimization of parameters for each cancer cell line, SK-MEL-28 (melanoma panel) gave the best prediction, with an average error of 0.97 for parameter combinations 2 and 3. The full results are shown in Table 2.
To identify the best parameter combination for each panel, the average |DTV(GI50)| values of cell lines obtained in each run were grouped by the panel, and the average |DTV(GI50)| for the entire panel was calculated for runs 1–18. The results of the best combinations of parameters for each panel are given in Table 3.
The renal cancer panel resulted the best prediction of all, with an average error of 1.15 using parameter combination 10.

2.1.4. Application of the AAP Tool for the Virtual Screening of an In-House Structure Database

The predictive capability of the AAP tool was exploited in the analysis of an in-house small molecule database to select compounds to be submitted to the NCI-DTP screening program. Recently, a few curcumin-like compounds were designed and biologically evaluated for neuroprotective and anti-Alzheimer properties, showing interesting results (Figure 6) [34,35,36].
Since curcumin has also been extensively studied in recent decades for its significant anticancer activity against various malignant cell types [39], the evaluation of curcumin-like compounds as antiproliferative agents against the NCI60 panel may be remarkably interesting. Indeed, many curcumin analogues have already been studied as antiproliferative agents, and some of them (e.g., EF24, UBS109, and CDF) showed higher activities and improved drug-like properties compared to curcumin (Figure 6) [40,41,42].
In detail, the selected in-house database included three different subclasses of curcumin-like derivatives, as reported in Figure 7: 1,2-diones (1ao); 1,2,4-oxadiazoles (2ah); and 1,3,4-oxadiazoles (3ah). They were developed by replacing the symmetrical β-diketone core of curcumin, which is responsible for unfavorable physicochemical properties and a weak pharmacokinetic profile [43,44,45,46,47,48], with stacked moieties such as heterocycles or α-diketones; these replacements have been shown to improve stability, solubility, oral absorption, and bioavailability [49].
After selecting the molecule database, the next step was to tune the parameters of the AAP protocol. As mentioned earlier, the GI50 calculation can be targeted to a specific cell line or class of compounds by optimizing the parameters N, G, and Z. In this light, curcumin, tested by NCI (NSC code 32982) and included in the training set (NCI2014DB), was selected as a reference compound to determinate the best combination of parameters for the CL module. The tuning was performed with the parameters in the ranges 250 < N < 800 and 50 < Z < 100 while considering a, b, or c for the G function. Eighteen runs were started, following the procedure described above for the internal validation (the combinations are reported in Table 1). Consequently, the total absolute deviations from the experimental GI50 values were calculated for each run, and the set of N = 760, Z = 50, and G = c was identified (total |DTV(GI50)| = 0.44) and applied to the AAP tool for the selected database (see Supporting Information S10).
The output of the AAP tool with the calculated GI50 values of the 39 in-house compounds is listed in Supporting Information S11. The AAP GI50 values of the curcumin-like analogues were compared with the experimental ones determined by the NCI for the curcumin lead compound (Supporting Information S11). The analysis of the GI50 mean of each compound for the full panel highlighted several curcumin-like molecules with predicted antiproliferative activities better than that of the reference curcumin (average GI50 value of 5.17), such as 1a (5.47), 1e (5.65), and 1m (5.82) for the diones class; 2a (5.49), 2d (5.57), and 2g (5.44) for the 1,2,4-oxadiazole class; and 3a (5.72), 3e (5.49), 3g (6.27), and 3h (5.47) for the 1,3,4-oxadiazole class. To test the consistency of the AAP protocol, both compounds classified as active and inactive (with a GI50 mean of less than 4.5) were synthesized and subsequently proposed to the NCI for the in vitro evaluation of the antiproliferative activity against the NCI60 human tumor cell lines, assuming that a reliable protocol must be able to identify both active and inactive compounds.

2.2. Chemistry

The three classes of compounds of types 1, 2, and 3 were synthetized as previously described in the literature [35,36,50]. Cinnamils (1,6-diarylhexa-1,5-diene-3,4-diones) 1ao were achieved by the aldol condensation of aromatic aldehydes 4ao with diacetyl 5, leading to the formation of the double bonds, both with E geometry (Scheme 1) [51].
The 1,2,4-oxadiazole derivatives 2ae were synthesized following the conventional amidoxime synthetic strategy [52], starting from the esters 6 and the amidoximes 7 (Scheme 2).
The 1,3,4-oxadiazoles 3ah were achieved through the one-pot synthesis described in Scheme 3. Diacylhydrazine intermediates were obtained by the reaction of the cinnamic acid analogues 8 and hydrazine. The subsequent cyclization led to the isolation of the regio-isomers (E) 3ah in good overall yields [34] (synthetic details and spectroscopic characterization for all compounds are reported in Supporting Information S12) [34,35,36].

2.3. Biological Assays: NCI60 Human Tumor Cell Lines Screen Selected Compounds

All synthesized curcumin-like compounds were submitted to NCI cell-line-based in vitro screening for anticancer drugs. As described in the Methods Section (compound selection guidelines paragraph), the NCI applied specific criteria for compound selection. In the case of analogues, the selected compounds were those that were the most representative of the series and had significant structural novelty compared to the NCI collection.
From the three series of curcumin-like derivatives, five molecules were selected by the NCI for the in vitro biological screening: 1a (NSC785541), 1b (NSC785539), and 1c (NSC785540) for the dione series; 2a (NSC785542) for the 1,2,4-oxadiazole series; and 3e (NSC785543) for the 1,3,4-oxadiazole series (Figure 8).

2.3.1. One-Dose Antiproliferative Assay

The NCI screening protocol consisted of a preliminary one-dose assay (concentration of 10 μM) against the full NCI60 panel. Compounds that met the NCI selection criteria and had a significant growth-inhibitory effect on a minimum number of cell lines proceed to the five-dose screening (experimental details are described in the Methods Section). The results are expressed as the percent of growth (G%) of the treated cells when compared to the untreated control cells. This parameter accurately expresses the anticancer potential of the drug. At G% > 100, the compound has no effect on cancer cell proliferation (inactive). In the range, the compound inhibits cell proliferation by a percentage indicated by 100-G%. When G% is <0, the compound is cytotoxic and lethal to the cancer cells. To graphically represent the most sensitive panels/cell lines, a mean growth percentage is also provided.
The mean G% values of the five selected compounds for each of the nine subgroups of cancer cell lines are shown in Table 4. The full results and the mean graphs from the one-dose screening are reported in Supporting Information S13.
Consistent with the AAP-predicted GI50 values for these compounds, the biological data confirmed the dimethoxy-dione 1a (NSC785541) and the dimethoxy-1,3,4-oxadiazole 3e (NSC785543) to be the most active curcumin-like derivatives. They showed remarkable overall average G% values (26.29 and 26.93, respectively) and an average inhibition of cell growth of about 75% against the full NCI60.
Compound 1a demonstrated high lethality against the colon cancer panel (average G% of −16.06), with the highest cytotoxic effect on cell lines HCT-116 (G% = −74.85) and HT29 (G% = −35.20) and against leukemia panel, with an average G% of 14.84 and an almost complete blockade of cell growth (G% ~0) for cell lines K562, RPMI-8226, and SR (see Supporting Information S13). Notable data were obtained for the LOX IMVI cell line (melanoma panel), with a G% of −66.41, and against the RXF 393 (renal cancer panel), with a G% of −63.27.
Similarly, dimethoxy-1,3,4-oxadiazole 3e showed a strong antiproliferative effect against the colon cancer panel but with lower toxicity; it exhibited a G% close to 0, implying an arrest of cell growth with low/no lethality against the most sensitive cell lines, HCT116 and HT-29. Furthermore, the leukemia, melanoma, and CNS cancer panels also showed interesting sensitivity to the tested compound. Therefore, according to the selection criteria of the DTP NCI protocol, these two molecules progressed to the full five-dose assay (see the next section).
The other three compounds tested in the NCI one-dose protocol, 1b, 1c, and 2a, generally exhibited less inhibitory activity, with an average G% close to 100.
These results confirm the AAP in silico data, according to which the compounds 1b and 1c were predicted to be almost inactive, with mean GI50 values of 4.38 and 4.54 (low millimolar range), respectively.
On the other hand, when it was analyzed for its antiproliferative effect on specific human cancer cell lines, dione 1b reduced the growth of the RPMI-8226 and MCF-7 cell lines by 55% and 76%, respectively (G% = 45.15 and 23.84) and induced a remarkable death in the HCT-116 colon cancer cell line (G% = −36.84). Compound 1c, instead, exhibited a G% of 69.01 against HCT-116 colon cancer cells.
The 1,2,4-oxadiazole 2a, which was predicted to be more active than curcumin, did not exhibit any appreciable anticancer activity against the NCI60 database, with the exception of HT29 colon cancer cells (G% of 56.32). This result was unexpected, considering that the isomer 3e was selected for investigation in the five-dose screening. It could probably be suggested that the change from the 1,3,4-oxadizole to the 1,2,4-oxadiazole core affects, in particular, the ability of the compound to interfere with biological targets.
Therefore, both compounds 1a and 3e were selected for the five-dose screening to measure the GI50s, which permitted us to further validate our tool.

2.3.2. Five-Dose Antiproliferative Assay for the Most Active Derivatives, 1a and 3e

The two selected compounds, 1a and 3e, were tested with the five-dose assay by measuring the percentage of cell growth at five different concentrations (from 10−8 to 10−4 M), as described in detail in the Methods Section. For each selected compound, NCI provided the measured GI50, TGI, and LC50 values against the NCI60 cell lines, with the corresponding mean graphs (see Supporting Information S14).
In the first part of this section, attention is focused exclusively on the GI50 values, as these data allowed us to further assess the predictive ability of the proposed AAP protocol.
The comparison of the average predicted GI50 with the experimental values obtained by NCI for the two compounds confirmed that the protocol was able to predict with high accuracy the range of activity of both compounds against the full NCI60 (for 1a, the average predicted GI50 was 5.39, whereas the average experimental GI50 was 5.49; for 3e, the average predicted GI50 was 5.41, whereas the average experimental GI50 was 5.28).
To further analyze the performance of the AAP protocol, the predicted GI50 values were matched to the experimental values; moreover, the |DTV(GI50)| was computed for the two tested compounds. This allowed us to calculate the average absolute error values for the compounds, which were 0.39 and 0.40 for 1a and 3e, respectively (see Supporting Information S15), indicating the capability to assign a GI50 value with an error of less than one order of magnitude.
In general, considering the mean |DTV(GI50)| for each panel, the AAP returned very low errors in activity prediction against specific panels, such as, for example, the prostate cancer panel for 1a (average |DTV(GI50)| for the panel of 0.07) and the colon cancer panel for 3e (average |DTV(GI50)| for the panel of 0.25). Furthermore, a detailed analysis of specific cell lines revealed that the protocol was able to predict GI50 values for some specific cell lines with remarkable precision: for 1a, the |DTV(GI50)| against HL-60TB and NCI-H322M was only 0.02; for 2a, a |DTV(GI50)| of only 0.03 was computed for HT-29 and OVCAR-5.
On the other hand, BT-549 (breast cancer panel) gave the worst predictions for both compounds, with |DTV(GI50)| values of 2.21 and 2.41; this evidence is consistent with the results presented in the previous section (tool validation), where this cell line showed the highest error. As previously demonstrated, high prediction error can be attributed to the lack of biological data for selected cell lines. This was confirmed by looking at the GI50(FP) values assigned to both compounds: the structure selected with the best score in the FP module was not tested against BT-549. Thus, the final GI50 relied solely on the CL protocol rather than combining outputs from both modules, which may have drastically affected the quality of the prediction.
In Figure 9a, the two bar graphs depicting the comparisons between predicted GI50 and experimental GI50 are reported to graphically appreciate the excellent performance of our protocol; in Figure 9b, instead, the mean error graphs are shown to highlight the cell lines for which the highest/lowest errors were recorded.
In Table 5, the full NCI output data, (GI50s, TGI, and LC50 values) of compounds 1a, 3e, and curcumin are shown (the mean graphs and the full NCI schedules are reported in Supporting Information S14). The analysis of the values allowed us to highlight the noteworthy antiproliferative activity of the two selected curcumin-like compounds compared to curcumin; indeed, the protocol predicted a higher activity with respect to the parent compound.
Through the evaluation of the GI50 values (the most diagnostic parameter used to compare antiproliferative potential), it emerged that the average GI50 values were higher for both the tested compounds, in the high micromolar range, than for curcumin (5.59 for 1a and 5.37 for 3e vs. 5.16 for curcumin), confirming the in silico predictions.
With respect to the tumor subpanels, the most active compound, the dione 1a, proved to be particularly effective against leukemia, colon cancer, and breast cancer. In fact, the calculated average GI50 values for these subpanels (5.87, 6.00, and 5.79, respectively) were always higher than the overall average GI50 value (5.59). In detail, among these subpanels, several cell lines showed remarkable sensitivity to the compound, with excellent GI50 values in the low micromolar range: RPMI-8226 (6.41), HCT-116 (6.5, the most sensitive cell line), HCT-15 (6.11), HT-29 (6.12), and MCF-7 (6.29). Moreover, the analysis of TGI values showed that the dione 1a was the most active compound of the series (overall average TGI of 4.81), confirming its selectivity against the aforementioned subpanels and cell lines, especially against colon cancer (average TGI of 5.48). This trend was also confirmed at the LC50 level, with a strong cytotoxic effect against colon cancer cell lines (average GI50 of 4.91 for colon cancer). Interestingly, it is important to highlight the very low toxicity against RPMI-8226 (LC50 < 4), which proved to be one of the most sensitive (GI50 in the sub-micromolar range), demonstrating the high potency and low toxicity of the compound against this cell line.
The curcumin-like 3e, although less active than the previous one, was more effective than curcumin. Remarkable results were obtained against the leukemia subpanel, with a panel average of 5.57, much higher than the average value for the full NCI60 (5.37). Moreover, the oxadiazole derivative exhibited high potency against two cell lines, MDA-MB-435 (melanoma) and A498 (renal cancer), with excellent GI50 values in the low micromolar range (6.03 and 6.61, respectively). In terms of the TGI level, it was slightly less effective than curcumin, but the average LC50 value of 4.01, even against the most susceptible cells, indicated high potency with low cytotoxicity, even at high concentrations.
To further evaluate the prediction capability of the AAP, the pdCSM-cancer tool was selected as a comparative approach. Thus, structures 1a and 3a were submitted to pdCSM-cancer and compared with the results obtained by the AAP. The obtained data are reported in Figure 10.
It is noteworthy that the AAP tool showed |DTV(GI50)| values of 0.39 and 0.41 for the compounds 1a and 3a, while the pdCSM-cancer tool showed values of 0.60 and 0.63, respectively. The cell line that showed a low grade of reliability was BT-549, not only in these cases in the study but also in the external test screening. A change in dataset could solve this problem.

3. Materials and Methods

3.1. Computational Studies

3.1.1. Hardware

The DRUDIT web service runs on four servers that are automatically selected according to the number of jobs and online availability. Each server can support up to 10 simultaneous jobs, while the exceeding jobs are placed in a queue.

3.1.2. Software

DRUDIT consists of several software modules implemented in C and JAVA and running on MacOS Mojave.

3.1.3. Database Selection and Dataset Building

The NCI60 database, containing both antiproliferative and chemical data of thousands of compounds, was selected as a reliable source for the building of the protocol.
In detail, since the presented tool is based on molecular descriptors, the 2D chemical structures of the NCI-tested compounds (.mol files, only available until the June 2016 release) and the corresponding growth inhibition data were retrieved from the NCI website (284,176 chemical structures) [38,53]. Among these thousands of compounds, only those tested with the five-dose assay, which provided GI50 data, were selected to build and validate the model. In particular, the structures were split into two sets: a training set containing more than 38 k compounds released until 2014 (NCI2014DB) was used to build the protocol, and a test set containing about 100 compounds that were first released in 2016 (NCI2016DB) was used to validate the AAP tool.

3.1.4. MOLDESTO: A New Software for Molecular Descriptor Calculations

MOLDESTO (molecular descriptor tool), as described previously [28], is a software tool implemented in DRUDIT that represents the evolution of our expertise in the calculation/manipulation of molecular descriptors [30]. It is currently able to calculate more than 1000 molecular descriptors (1D, 2D, and 3D) for each input structure (the full list of molecular descriptors calculated by MOLDESTO is reported in Supporting Information S1). The input structures can be drawn directly in the web interface or uploaded as commonly used molecule file formats (e.g., SMILES, SDF, Inchi, Mdl, and Mol2). The software is provided with a caching system to boost the calculation speed of previously submitted structures.

3.1.5. DRUDIT Settings for Antiproliferative Activity Predictor (AAP) Tool

The AAP tool comprises the fingerprint (FP) and cell line (CL) modules, which cooperate simultaneously to assign the predicted GI50 values to an input structure. In each module, the performed calculation is dynamic; indeed, it can be modulated by appropriately tuning the values of the available parameters (three for each module, see below).
The FP module parameters are a choice of biological activity, such as GI50, TGI, LC50, or G% (in this work only the first choice was considered); N (-b), the best number of the dynamically selected molecular descriptors; Z (-m), the number of descriptors for which |v-m|/m < <value> applies (v: descriptor value, m: target mean); and G (-c), the max number of zero percentage values per descriptor. The DRUDIT parameters for the CL module are a choice of biological activity as GI50, TGI, LC50, or G% (in this work only the first parameter was considered); N (-b), the best number of dynamically selected molecular descriptors; Z (-m), the max number of zero percentage values per descriptor; and G (-f), the Gaussian smoothing function to be used (a, b, or c mode).

3.2. Chemistry

All solvents and reagents were used as received unless otherwise stated. Melting points were determined on a hot-stage apparatus. The 1H-NMR and 13C-NMR spectra were recorded at the indicated frequencies; the residual solvent peak was used as a reference. Chromatography was performed using silica gel (0.040–0.063 mm) and mixtures of ethyl acetate and petroleum ether (fraction boiling in the range of 40–60 °C) in various ratios (v/v). Compounds 1do [36], 2ac [36], 2d [35], 3ag [36], and 3h [35] were prepared as previously reported. Compounds 1ac and 2eh were prepared by adapting previously reported methods. The synthetic details and spectroscopic characterizations of all compounds are reported in Supporting Information S12.

3.3. NCI60 Antiproliferative Screenings

3.3.1. Compound Selection Guidelines

The compounds to be screened were selected according to precise and rigorous guidelines; in general, submission was encouraged for molecules that bring some novelties (novel heterocyclic ring systems and privileged scaffolds) to the NCI collection and compounds that emerged from computer-aided drug design. In addition, in the case of a series of analogues, it was preferred to select only the one that was expected to provide the greatest information. On the other hand, the submission of compounds with the following features was discouraged: excessive flexibility; the presence of non-drug-like functional groups (nitro, nitroso, diazo, imine, etc.); and the presence of chemical portions that could affect the reliability of the assays (PAINS) [54].

3.3.2. One-Dose Assay

All compounds submitted to NCI were first assayed against the NCI60 DB in a one-dose screen (concentration of 105 M); this kind of assay aims to determine the G% (growth inhibition percent) of the compounds against the considered cells. The results were plotted in a one-dose graph showing the G% of the single compound against the 60 cell lines. This first assay was considered passed only for the most promising compounds (satisfaction of predetermined threshold criteria); in this case, the compound passed to the five-dose screen (for further experimental details about the standardized assay procedures, see [55,56]).

3.3.3. Five-Dose Assay

The most active compounds were submitted to a multiple-dose screen using five different concentrations (ranging from 10−8 to 10−4 M). The dose–response curves obtained from this assay permitted the extrapolation of the GI50 (the molar concentration of the compound that inhibits 50% of cell growth), TGI (the molar concentration of the compound leading to total inhibition of cell growth), and LC50 (the molar concentration of the compound that induces 50% cell death) values of the selected compounds against each cancer cell line. For each of the mentioned parameters, a mean graph midpoint (MG_MID) was calculated, providing an average activity parameter over all cell lines (for further experimental details about the standardized assay procedures, see [55,56]).

4. Conclusions

In the field of antitumor drug discovery, in vitro antiproliferative assays still represent one of the most important tools for identifying new small molecules against cancer. In recent years, numerous projects and online databases have been launched aiming to collect both drug response and cellular data, with NCI60 undoubtedly being the best-known and most complete [1].
To avoid the high failure rate and the enormous number of resources invested to perform such intensive in vitro screenings, computational chemistry and biology, in the last few years, have sought to develop in silico techniques that are able to predict the antiproliferative potential of new small molecules in the early stage of the drug discovery pipeline.
In this light, we have presented the antiproliferative activity predictor (AAP), a new molecular-descriptor-based tool capable of reliably predicting the anticancer potential (expressed as GI50, the most diagnostic parameter to measure anticancer drug responses) of small molecule libraries against the full NCI60 DB.
Using both the structural and antiproliferative data (GI50) of thousands of compounds stored in the NCI DB, we applied our expertise in manipulating molecular descriptors to build two convergent modules, the FP (fingerprint) and the CL (cell line); as shown, these operate synergistically to assign GI50s values to the input structures.
Both internal and external validation were performed to validate the CL module and the entire tool, demonstrating the reliability and the robustness of the proposed protocol. Interestingly, the possibility to appropriately tune the available parameters (N, Z, and G) allows researchers to address the screening to a specific subpanel/cell line or class of compounds. An important aspect that should be highlighted is the deep correlation between the quality of the prediction and the availability of biological data. In the case of the lack of sufficient GI50 data, as for a cell line (M19-MEL) or for some ranges of activity (structures with GI50 values in the range of 7–8), the prediction is negatively affected.
Moreover, the application of the AAP tool to quickly screen an in-house structure database of curcumin-like compounds permitted us to further corroborate the already obtained encouraging results: compounds 1a and 3e, predicted to be highly active by the protocol, were found to be significantly antiproliferative for several human cancer cell lines in the five-dose NCI screen. This further goal confirmed that the AAP could be an invaluable help in the in silico design of new effective anticancer small molecules, permitting the selection of the most promising molecules in the first phases of the drug design process.
It is worth noting that the integration of the AAP tool into our free and easy-to-use DRUDIT web service (available online at https://www.drudit.com (accessed on 15 November 2022)) provides an interesting tool for the entire medicinal chemistry community.
In the future, we plan to continuously update the training set with new NCI-tested compounds to cover more and more chemical spaces, activities, and cell lines. Finally, the extension of the AAP’s potential for the calculation of more parameters, such as TGI or LC50, could also be interesting. The possibility to calculate not only the GI50 but also these parameters (evaluated in vitro by the NCI) could allow the identification and eventually the elimination of small molecules with toxicity profiles that are unacceptable for further development.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232214374/s1.

Author Contributions

Conceptualization, A.L. and D.P.; Data curation, A.M., G.L.M., C.G., A.L. and D.P.; Formal analysis, C.G.; Investigation, A.P.P.; Methodology, S.M., A.L. and D.P.; Project administration, A.L.; Resources, A.M. and A.P.P.; Supervision, A.L.; Validation, A.L.; Visualization, S.B. and C.G.; Writing—original draft, A.M., G.L.M., A.B., S.B. and A.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by PJ_MIN_SALUTE_METODI_SOSTITUTIVI—CUP: B75F21002350001 and PJ_RIC_FFABR_2022_005832 Grant—University of Palermo.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shoemaker, R.H. The NCI60 human tumour cell line anticancer drug screen. Nat. Rev. Cancer 2006, 6, 813–823. [Google Scholar] [CrossRef]
  2. Barretina, J.; Caponigro, G.; Stransky, N.; Venkatesan, K.; Margolin, A.A.; Kim, S.; Wilson, C.J.; Lehár, J.; Kryukov, G.V.; Sonkin, D.; et al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 2012, 483, 603–607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yang, W.; Soares, J.; Greninger, P.; Edelman, E.J.; Lightfoot, H.; Forbes, S.; Bindal, N.; Beare, D.; Smith, J.A.; Thompson, I.R.; et al. Genomics of Drug Sensitivity in Cancer (GDSC): A resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res. 2013, 41, D955–D961. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tomczak, K.; Czerwińska, P.; Wiznerowicz, M. The Cancer Genome Atlas (TCGA): An immeasurable source of knowledge. Contemp. Oncol. 2015, 19, A68–A77. [Google Scholar] [CrossRef] [PubMed]
  5. Basu, A.; Bodycombe, N.E.; Cheah, J.H.; Price, E.V.; Liu, K.; Schaefer, G.I.; Ebright, R.Y.; Stewart, M.L.; Ito, D.; Wang, S.; et al. An interactive resource to identify cancer genetic and lineage dependencies targeted by small molecules. Cell 2013, 154, 1151–1161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Gillet, J.P.; Varma, S.; Gottesman, M.M. The clinical relevance of cancer cell lines. J. Natl. Cancer Inst. 2013, 105, 452–458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Mirabelli, P.; Coppola, L.; Salvatore, M. Cancer cell lines are useful model systems for medical research. Cancers 2019, 11, 1098. [Google Scholar] [CrossRef] [Green Version]
  8. Goodspeed, A.; Heiser, L.M.; Gray, J.W.; Costello, J.C. Tumor-derived cell lines as molecular models of cancer pharmacogenomics. Mol. Cancer Res. 2016, 14, 3–13. [Google Scholar] [CrossRef] [Green Version]
  9. Takimoto, C.H. Anticancer drug development at the US National Cancer Institute. Cancer Chemother. Pharmacol. 2003, 52 (Suppl. 1), S29–S33. [Google Scholar] [CrossRef] [PubMed]
  10. NCI-60 Human Tumor Cell Lines Screen—Introduction. Available online: https://dtp.cancer.gov/discovery_development/nci-60/ (accessed on 15 November 2022).
  11. Molecular Characterization of the NCI-60. Available online: https://dtp.cancer.gov/discovery_development/nci-60/characterization.htm (accessed on 15 November 2022).
  12. Covell, D.G.; Huang, R.; Wallqvist, A. Anticancer medicines in development: Assessment of bioactivity profiles within the National Cancer Institute anticancer screening data. Mol. Cancer Ther. 2007, 6, 2261–2270. [Google Scholar] [CrossRef]
  13. NCI60 Human Timor Cell Lines Screen—Cell Lines in the In Vitro Screen. Available online: https://dtp.cancer.gov/discovery_development/nci-60/cell_list.htm (accessed on 15 November 2022).
  14. Firoozbakht, F.; Yousefi, B.; Schwikowski, B. An overview of machine learning methods for monotherapy drug response prediction. Brief. Bioinform. 2021, 23, bbab408. [Google Scholar] [CrossRef]
  15. Paull, K.D.; Shoemaker, R.H.; Hodes, L.; Monks, A.; Scudiero, D.A.; Rubinstein, L.; Plowman, J.; Boyd, M.R. Display and analysis of patterns of differential activity of drugs against human tumor cell lines: Development of mean graph and COMPARE algorithm. J. Natl. Cancer Inst. 1989, 81, 1088–1092. [Google Scholar] [CrossRef] [PubMed]
  16. Zaharevitz, D.W.; Holbeck, S.L.; Bowerman, C.; Svetlik, P.A. COMPARE: A web accessible tool for investigating mechanisms of cell growth inhibition. J. Mol. Graph. Model. 2002, 20, 297–303. [Google Scholar] [CrossRef]
  17. NCI-60 Analysis Tools—CellMiner. Available online: http://discover.nci.nih.gov/cellminer/ (accessed on 15 November 2022).
  18. Reinhold, W.C.; Sunshine, M.; Liu, H.; Varma, S.; Kohn, K.W.; Morris, J.; Doroshow, J.; Pommier, Y. CellMiner: A web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Cancer Res. 2012, 72, 3499–3511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Luna, A.; Elloumi, F.; Varma, S.; Wang, Y.; Rajapakse, V.N.; Aladjem, M.I.; Robert, J.; Sander, C.; Pommier, Y.; Reinhold, W.C. CellMiner Cross-Database (CellMinerCDB) version 1.2: Exploration of patient-derived cancer cell line pharmacogenomics. Nucleic Acids Res. 2021, 49, D1083–D1093. [Google Scholar] [CrossRef]
  20. Lind, A.P.; Anderson, P.C. Predicting drug activity against cancer cells by random forest models based on minimal genomic information and chemical properties. PLoS ONE 2019, 14, e0219774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Zhang, N.; Wang, H.; Fang, Y.; Wang, J.; Zheng, X.; Liu, X.S. Predicting anticancer drug responses using a dual-layer integrated cell line-drug network model. PLoS Comput. Biol. 2015, 11, e1004498. [Google Scholar] [CrossRef] [PubMed]
  22. Cadow, J.; Born, J.; Manica, M.; Oskooei, A.; Rodríguez Martínez, M. PaccMann: A web service for interpretable anticancer compound sensitivity prediction. Nucleic Acids Res. 2020, 48, W502–W508. [Google Scholar] [CrossRef] [PubMed]
  23. Chang, Y.; Park, H.; Yang, H.J.; Lee, S.; Lee, K.Y.; Kim, T.S.; Jung, J.; Shin, J.M. Cancer Drug Response Profile scan (CDRscan): A deep learning model that predicts drug effectiveness from cancer genomic signature. Sci. Rep. 2018, 8, 8857. [Google Scholar] [CrossRef] [Green Version]
  24. Joo, M.; Park, A.; Kim, K.; Son, W.J.; Lee, H.S.; Lim, G.; Lee, J.; Lee, D.H.; An, J.; Kim, J.H.; et al. A deep learning model for cell growth inhibition IC50 prediction and its application for gastric cancer patients. Int. J. Mol. Sci. 2019, 20, 6276. [Google Scholar] [CrossRef]
  25. Al-Jarf, R.; de Sá, A.G.C.; Pires, D.E.V.; Ascher, D.B. pdCSM-cancer: Using graph-based signatures to identify small molecules with anticancer properties. J. Chem. Inf. Model. 2021, 61, 3314–3322. [Google Scholar] [CrossRef] [PubMed]
  26. Lauria, A.; Tutone, M.; Almerico, A.M. Virtual lock-and-key approach: The in silico revival of Fischer model by means of molecular descriptors. Eur. J. Med. Chem. 2011, 46, 4274–4280. [Google Scholar] [CrossRef] [PubMed]
  27. Lauria, A.; Tutone, M.; Barone, G.; Almerico, A.M. Multivariate analysis in the identification of biological targets for designed molecular structures: The BIOTA protocol. Eur. J. Med. Chem. 2014, 75, 106–110. [Google Scholar] [CrossRef] [PubMed]
  28. Lauria, A.; Mannino, S.; Gentile, C.; Mannino, G.; Martorana, A.; Peri, D. DRUDIT: Web-based DRUgs DIscovery Tools to design small molecules as modulators of biological targets. Bioinformatics 2020, 36, 1562–1569. [Google Scholar] [CrossRef] [PubMed]
  29. Lauria, A.; Abbate, I.; Patella, C.; Martorana, A.; Dattolo, G.; Almerico, A.M. New annelated thieno[2,3-e][1,2,3]triazolo[1,5-a]pyrimidines, with potent anticancer activity, designed through VLAK protocol. Eur. J. Med. Chem. 2013, 62, 416–424. [Google Scholar] [CrossRef] [PubMed]
  30. Lauria, A.; Patella, C.; Abbate, I.; Martorana, A.; Almerico, A.M. Lead optimization through VLAK protocol: New annelated pyrrolo-pyrimidine derivatives as antitumor agents. Eur. J. Med. Chem. 2012, 55, 375–383. [Google Scholar] [CrossRef] [PubMed]
  31. Lauria, A.; Abbate, I.; Gentile, C.; Angileri, F.; Martorana, A.; Almerico, A.M. Synthesis and biological activities of a new class of heat shock protein 90 inhibitors, designed by energy-based pharmacophore virtual screening. J. Med. Chem. 2013, 56, 3424–3428. [Google Scholar] [CrossRef] [PubMed]
  32. Diana, P.; Martorana, A.; Barraja, P.; Montalbano, A.; Carbone, A.; Cirrincione, G. Nucleophilic substitutions in the isoindole series as a valuable tool to synthesize derivatives with antitumor activity. Tetrahedron 2011, 67, 2072–2080. [Google Scholar] [CrossRef]
  33. Mingoia, F.; Di Sano, C.; Di Blasi, F.; Fazzari, M.; Martorana, A.; Almerico, A.M.; Lauria, A. Exploring the anticancer potential of pyrazolo[1,2-a]benzo[1,2,3,4]tetrazin-3-one derivatives: The effect on apoptosis induction, cell cycle and proliferation. Eur. J. Med. Chem. 2013, 64, 345–356. [Google Scholar] [CrossRef] [PubMed]
  34. Kayed, R.; Lo Cascio, F.; Piccionello Palumbo, A.; Pace, A. Novel Small Molecules That Bind and/or Modulate Different Forms of Tau Oligomers. Patent WO2020/219714 A1, 23 April 2020. [Google Scholar]
  35. Battisti, A.; Palumbo Piccionello, A.; Sgarbossa, A.; Vilasi, S.; Ricci, C.; Ghetti, F.; Spinozzi, F.; Marino Gammazza, A.; Giacalone, V.; Martorana, A.; et al. Curcumin-like compounds designed to modify amyloid beta peptide aggregation patterns. RSC Adv. 2017, 7, 31714–31724. [Google Scholar] [CrossRef]
  36. Lo Cascio, F.; Puangmalai, N.; Ellsworth, A.; Bucchieri, F.; Pace, A.; Palumbo Piccionello, A.; Kayed, R. Toxic tau oligomers modulated by novel curcumin derivatives. Sci. Rep. 2019, 9, 19011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lauria, A.; Martorana, A.; La Monica, G.; Mannino, S.; Mannino, G.; Peri, D.; Gentile, C. In silico identification of small molecules as new Cdc25 inhibitors through the correlation between chemosensitivity and protein expression pattern. Int. J. Mol. Sci. 2021, 22, 3714. [Google Scholar] [CrossRef] [PubMed]
  38. NCI DTP Chemical Data. Available online: https://wiki.nci.nih.gov/display/NCIDTPdata/Chemical+Data (accessed on 15 November 2022).
  39. Sharma, R.A.; Gescher, A.J.; Steward, W.P. Curcumin: The story so far. Eur. J. Cancer 2005, 41, 1955–1968. [Google Scholar] [CrossRef] [PubMed]
  40. Adeluola, A.; Zulfiker, A.H.M.; Brazeau, D.; Amin, A.R.M.R. Perspectives for synthetic curcumins in chemoprevention and treatment of cancer: An update with promising analogues. Eur. J. Pharmacol. 2021, 906, 174266. [Google Scholar] [CrossRef]
  41. Ahsan, M.J.; Choudhary, K.; Jadav, S.S.; Yasmin, S.; Ansari, M.Y.; Sreenivasulu, R. Synthesis, antiproliferative activity, and molecular docking studies of curcumin analogues bearing pyrazole ring. Med. Chem. Res. 2015, 24, 4166–4180. [Google Scholar] [CrossRef]
  42. Ahsan, M.J. Evaluation of anticancer activity of curcumin analogues bearing a heterocyclic nucleus. Asian Pac. J. Cancer Prev. 2016, 17, 1739–1744. [Google Scholar] [CrossRef] [Green Version]
  43. Anand, P.; Kunnumakkara, A.B.; Newman, R.A.; Aggarwal, B.B. Bioavailability of curcumin: Problems and promises. Mol. Pharm. 2007, 4, 807–818. [Google Scholar] [CrossRef]
  44. Sabet, S.; Rashidinejad, A.; Melton, L.D.; McGillivray, D.J. Recent advances to improve curcumin oral bioavailability. Trends Food Sci. Technol. 2021, 110, 253–266. [Google Scholar] [CrossRef]
  45. Sanidad, K.Z.; Sukamtoh, E.; Xiao, H.; McClements, D.J.; Zhang, G. Curcumin: Recent advances in the development of strategies to improve oral bioavailability. Annu. Rev. Food Sci. Technol. 2019, 10, 597–617. [Google Scholar] [CrossRef] [Green Version]
  46. Kotha, R.R.; Luthria, D.L. Curcumin: Biological, pharmaceutical, nutraceutical, and analytical aspects. Molecules 2019, 24, 2930. [Google Scholar] [CrossRef]
  47. Nelson, K.M.; Dahlin, J.L.; Bisson, J.; Graham, J.; Pauli, G.F.; Walters, M.A. The essential medicinal chemistry of curcumin. J. Med. Chem. 2017, 60, 1620–1637. [Google Scholar] [CrossRef] [PubMed]
  48. Olotu, F.; Agoni, C.; Soremekun, O.; Soliman, M.E.S. An update on the pharmacological usage of curcumin: Has it failed in the drug discovery pipeline? Cell Biochem. Biophys. 2020, 78, 267–289. [Google Scholar] [CrossRef] [PubMed]
  49. Rodrigues, F.C.; Kumar, N.A.; Thakur, G. The potency of heterocyclic curcumin analogues: An evidence-based review. Pharmacol. Res. 2021, 166, 105489. [Google Scholar] [CrossRef] [PubMed]
  50. Lo Cascio, F.; Garcia, S.; Montalbano, M.; Puangmalai, N.; McAllen, S.; Pace, A.; Palumbo Piccionello, A.; Kayed, R. Modulating disease-relevant tau oligomeric strains by small molecules. J. Biol. Chem. 2020, 295, 14807–14825. [Google Scholar] [CrossRef] [PubMed]
  51. Sinu, C.R.; Padmaja, D.V.; Ranjini, U.P.; Seetha Lakshmi, K.C.; Suresh, E.; Nair, V. A cascade reaction actuated by nucleophilic heterocyclic carbene catalyzed intramolecular addition of enals via homoenolate to α,β-unsaturated esters: Efficient synthesis of coumarin derivatives. Org. Lett. 2013, 15, 68–71. [Google Scholar] [CrossRef] [PubMed]
  52. Pace, A.; Buscemi, S.; Piccionello, A.P.; Pibiri, I. Recent advances in the chemistry of 1,2,4-oxadiazoles. Adv. Heterocycl. Chem. 2015, 116, 85–136. [Google Scholar] [CrossRef]
  53. NCI60 Growth Inhibition Data—Download NCI Cell Line Data. Available online: https://wiki.nci.nih.gov/display/NCIDTPdata/NCI-60+Growth+Inhibition+Data (accessed on 15 November 2022).
  54. Compound Submission for NCI60 Testing—Selection Guidelines. Available online: https://dtp.cancer.gov/organization/dscb/compoundSubmission/structureSelection.htm (accessed on 15 November 2022).
  55. NCI60 Screening Methodology—NCI60 Cell One/Five Doses Screen. Available online: https://dtp.cancer.gov/discovery_development/nci-60/methodology.htm (accessed on 15 November 2022).
  56. NCI-60 Human Cancer Cell Line Screen—Standard Operating Procedures for Sample Preparation for NCI60 Screen. Available online: https://dtp.cancer.gov/discovery_development/nci-60/handling.htm (accessed on 15 November 2022).
Figure 1. Flowchart of the antiproliferative predictor protocol: GI50i is the GI50 value for cell line, S is the fingerprint score, GI50i(FP) is the GI50 value predicted by the FP module, and GI50i(CL) is the GI50 value assigned by the CL module.
Figure 1. Flowchart of the antiproliferative predictor protocol: GI50i is the GI50 value for cell line, S is the fingerprint score, GI50i(FP) is the GI50 value predicted by the FP module, and GI50i(CL) is the GI50 value assigned by the CL module.
Ijms 23 14374 g001
Figure 2. NCI selection data used in the AAP tool: compounds screened in five-dose assay and published until 2014 were used as training set; the 5‰ of these structures was used for the internal validation of the protocol (panels A,C,D); those structures screened in five-dose assay and published until 2016 were used as test set to evaluate the predictive performance of the entire protocol (panels B,E).
Figure 2. NCI selection data used in the AAP tool: compounds screened in five-dose assay and published until 2014 were used as training set; the 5‰ of these structures was used for the internal validation of the protocol (panels A,C,D); those structures screened in five-dose assay and published until 2016 were used as test set to evaluate the predictive performance of the entire protocol (panels B,E).
Ijms 23 14374 g002
Figure 3. Template building process.
Figure 3. Template building process.
Ijms 23 14374 g003
Figure 4. Graphical representation of GI50 prediction by the CL module for a cell line. D1, D2…Dn: molecular descriptors values for the input structure.
Figure 4. Graphical representation of GI50 prediction by the CL module for a cell line. D1, D2…Dn: molecular descriptors values for the input structure.
Ijms 23 14374 g004
Figure 5. |DTV(GI50)| vs. frequencies of the data for the external test validation.
Figure 5. |DTV(GI50)| vs. frequencies of the data for the external test validation.
Ijms 23 14374 g005
Figure 6. Curcumin and curcumin-like biologically active compounds.
Figure 6. Curcumin and curcumin-like biologically active compounds.
Ijms 23 14374 g006
Figure 7. In-house structure database of curcumin-like compounds investigated by the AAP tool.
Figure 7. In-house structure database of curcumin-like compounds investigated by the AAP tool.
Ijms 23 14374 g007
Scheme 1. Synthesis of cinnamils 1ao.
Scheme 1. Synthesis of cinnamils 1ao.
Ijms 23 14374 sch001
Scheme 2. Synthesis of 1,2,4-oxadiazole derivatives 2ae.
Scheme 2. Synthesis of 1,2,4-oxadiazole derivatives 2ae.
Ijms 23 14374 sch002
Scheme 3. Synthesis of 1,3,4-oxadiazole derivatives 3ae.
Scheme 3. Synthesis of 1,3,4-oxadiazole derivatives 3ae.
Ijms 23 14374 sch003
Figure 8. Chemical structures of the five curcumin-like compounds selected by NCI for the one-dose antiproliferative assay.
Figure 8. Chemical structures of the five curcumin-like compounds selected by NCI for the one-dose antiproliferative assay.
Ijms 23 14374 g008
Figure 9. (a) Comparison between AAP-predicted GI50 values and the corresponding experimental GI50 values measured by NCI for the two selected compounds, 1a and 3e (inside each bar, the corresponding GI50 value is indicated). (b) Mean error graphs for the two selected compounds, 1a and 3e.
Figure 9. (a) Comparison between AAP-predicted GI50 values and the corresponding experimental GI50 values measured by NCI for the two selected compounds, 1a and 3e (inside each bar, the corresponding GI50 value is indicated). (b) Mean error graphs for the two selected compounds, 1a and 3e.
Ijms 23 14374 g009aIjms 23 14374 g009b
Figure 10. Performance of AAP vs. pdCSM-cancer tools: blue and red vertical lines indicate the |DTV(GI50)| for the AAP and pdCSM tools, respectively.
Figure 10. Performance of AAP vs. pdCSM-cancer tools: blue and red vertical lines indicate the |DTV(GI50)| for the AAP and pdCSM tools, respectively.
Ijms 23 14374 g010
Table 1. Overview of DRUDIT runs for parameter tuning. Absolute deviation from true values (|DTV(GI50)|) is reported, and run codes (1–18) are given in parentheses.
Table 1. Overview of DRUDIT runs for parameter tuning. Absolute deviation from true values (|DTV(GI50)|) is reported, and run codes (1–18) are given in parentheses.
ZNG
abc
502401.22 (1)1.23 (2)1.23 (3)
5001.22 (7)1.30 (8)1.31 (9)
7601.32 (13)1.44 (14)1.42 (15)
1002401.31 (4)1.64(5)1.72 (6)
5001.23 (10)1.51 (11)1.53 (12)
7601.28 (16)1.42 (17)1.44 (18)
Table 2. Tuning of DRUDIT parameters for each cancer cell line.
Table 2. Tuning of DRUDIT parameters for each cancer cell line.
PANELSCELL LINESRUNAVERAGE |DTV(GI50)|
Breast CancerBT-54941.35
HS-578T31.30
MCF71/71.30
MDA-MB-231-ATCC71.22
T-47D71.16
CNS CancerSF-26871.17
SF-29511.25
SF-53941.18
SNB-1971.15
SNB-7511.16
U25111.24
Colon CancerCOLO-205101.13
HCC-299811.09
HCT-1162/71.13
HCT-1521.21
HT2911.14
KM1211.19
SW-62021.14
LeukemiaCCRF-CEM71.13
HL-60TB71.22
K-56221.27
MOLT-431.12
RPMI-8226101.12
SR21.28
MelanomaLOX-IMVI31.16
M141/31.20
MALME-3M101.19
MDA-MB-43531.22
SK-MEL-231.03
SK-MEL-282/30.97
SK-MEL-521.26
UACC-25721.07
UACC-62101.31
Non-Small-Cell Lung CancerA549-ATCC31.18
EKVX11.02
HOP-621/81.19
HOP-92101.21
NCI-H22671.07
NCI-H2341.16
NCI-H322M11.10
NCI-H46021.26
NCI-H52271.09
Ovarian CancerIGROV11/31.25
NCI-ADR-RES41.31
OVCAR-31/41.22
OVCAR-471.00
OVCAR-5161.02
OVCAR-811.14
SK-OV-371.18
Prostate CancerDU-145101.19
PC-321.19
Renal Cancer786-0101.16
A49811.21
ACHN71.19
CAKI-111.11
RXF-393101.12
SN12C1/101.16
TK-10100.99
UO-3121.16
Table 3. Best parameter combinations for each panel.
Table 3. Best parameter combinations for each panel.
PANELSRUNAVERAGE|DTV(GI50)|
Breast Cancer1/31.37
CNS Cancer11.23
Colon Cancer11.19
Leukemia21.23
Melanoma31.18
Non-Small-Cell Lung Cancer2/71.20
Ovarian Cancer11.21
Prostate Cancer11.23
Renal Cancer101.15
Table 4. G% values determined for the five selected compounds against the NCI60 panel with the one-dose assay.
Table 4. G% values determined for the five selected compounds against the NCI60 panel with the one-dose assay.
PANEL 11a1b1c2a3e
Leukemia14.8478.7196.3777.4718.53
Non-Small-Cell Lung Cancer61.7998.0295.4784.3929.86
Colon Cancer−16.0680.6595.8486.5221.07
CNS Cancer32.9898.55101.4098.1918.66
Melanoma24.2697.49100.5996.7422.40
Ovarian Cancer46.48103.35101.4595.2932.68
Renal Cancer21.07100.02100.3194.3634.17
Prostate Cancer33.06102.04101.6288.8940.34
Breast Cancer18.2186.7799.6288.9324.71
Overall average26.2993.9699.1990.0926.93
1 For each compound, the average G% values against the nine subpanels are reported.
Table 5. DTP NCI five-dose screening for compounds 1a, 3e, and curcumin.
Table 5. DTP NCI five-dose screening for compounds 1a, 3e, and curcumin.
PANELCELL LINE 11a3eCurcumin
GI50TGILC50GI50TGILC50GI50TGILC50
LeukemiaCCRF-CEM5.654.745.4445.524.814
HL-60(TB)5.634.7745.635.0745.144.604.04
K-5625.97445.74445.514.264
MOLT-45.694.7245.49445.334.754.12
RPMI-82266.415.6345.58445.685.204
Panel average5.874.7645.574.2145.434.724.03
Non-Small-Cell Lung CancerA549/ATCC4445.32444.894.504.11
EKVX5.3445.21444.824.454.10
HOP-625.17445.28445.444.724.24
HOP-924.824.0845.64.634NTNTNT
NCI-H2265.63NT 144.83444.734.274
NCI-H235.52445.36445.254.504
NCI-H322M5.17444.84444.784.494.21
NCI-H4605.514.9545.37445.094.644.22
NCI-H5225.434.7245.725.174.025.274.784.07
Panel average5.174.224.005.284.2045.034.544.12
Colon CancerCOLO-2055.725.314.45.084.0144.874.544.21
HCC-29985.775.495.224.71445.525.094.53
HCT-1166.55.885.395.594.7545.535.034.28
HCT-156.115.214.255.56445.394.734.14
HT-296.125.585.095.584.9745.294.494
KM125.85.435.075.16445.274.714.19
SW-6205.975.484.975.48445.384.674.07
Panel average6.005.484.915.314.2545.324.754.20
CNS CancerSF-2685.65.0745.11445.154.444
SF-2955.524.5145.534.7545.104.684.32
SF-5395.675.294.225.575.034.095.555.054.48
SNB-195.514.6345.394.0445.054.614.20
SNB-755.64.4745.384.3445.174.744.35
U2515.815.445.075.424.7345.334.784.31
Panel average5.624.904.225.404.484.025.224.724.28
MelanomaLOX IMVI5.845.495.155.514.5445.575.074
MALME-3M5.244.145.21444.854.564.27
M145.775.364.245.514.5845.424.804.35
MDA-MB-4355.945.535.116.035.414.195.534.924.40
SK-MEL-24.51445.544.7344.784.394.06
SK-MEL-285.745.4NT5.3445.354.804.30
SK-MEL-55.675.2145.664.9945.064.654.28
UACC-2575.615.1344.97444.944.624.31
UACC-625.725.314.525.62545.194.694.26
Panel average5.565.064.385.484.584.025.194.724.25
Ovarian CancerIGROV-15.57NT45.37445.104.574.09
OVCAR-35.55545.34.0945.184.614.17
OVCAR-45.39445.04445.034.444
OVCAR-55.65.0945.314.2744.784.454.12
OVCAR-85.53445.06445.134.554.08
NCI/ADR-RES5.38445.45445.144.124
SK-OV-34.44445.064.0345.054.684.33
Panel average5.354.354.005.234.0645.064.494.11
Renal Cancer786-05.645.145.25445.484.974.42
A-4985.654.746.625.0444.804.484.16
ACHN5.695.2744.86444.914.544.17
CAKI-15.5444.94444.924.604.30
RXF-3935.835.525.25.14.1145.524.954.27
SN12C5.514.7545.26445.084.604.20
TK-105.424.645.424.5644.854.514.18
UO-315.795.475.145.32444.954.614.27
Panel average5.634.934.295.354.2145.064.664.25
Prostate CancerPC-35.48445.39445.064.594.15
DU-1455.514.8544.82444.814.534.25
Panel average5.504.434.005.114.0044.934.564.20
Breast CancerMCF76.295.0245.46445.484.464
MDA-MB-231/ATCC5.655.2245.444.3644.754.254
HS 578T5.57445.434.3244.964.234
BT-5495.745.354.665.544.4545.304.864.37
T-47D5.65445.464.0945.084.334
MDA-MB-4685.815.4345.544.564NTNTNT
Panel average5.794.844.115.484.3045.114.434.07
Overall average5.594.814.245.374.284.015.164.634.17
Range4–6.54–5.884–5.394.71–6.624–5.414–4.194.73–5.684–5.204–4.53
1 NT = not tested against the cell line.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martorana, A.; La Monica, G.; Bono, A.; Mannino, S.; Buscemi, S.; Palumbo Piccionello, A.; Gentile, C.; Lauria, A.; Peri, D. Antiproliferative Activity Predictor: A New Reliable In Silico Tool for Drug Response Prediction against NCI60 Panel. Int. J. Mol. Sci. 2022, 23, 14374. https://doi.org/10.3390/ijms232214374

AMA Style

Martorana A, La Monica G, Bono A, Mannino S, Buscemi S, Palumbo Piccionello A, Gentile C, Lauria A, Peri D. Antiproliferative Activity Predictor: A New Reliable In Silico Tool for Drug Response Prediction against NCI60 Panel. International Journal of Molecular Sciences. 2022; 23(22):14374. https://doi.org/10.3390/ijms232214374

Chicago/Turabian Style

Martorana, Annamaria, Gabriele La Monica, Alessia Bono, Salvatore Mannino, Silvestre Buscemi, Antonio Palumbo Piccionello, Carla Gentile, Antonino Lauria, and Daniele Peri. 2022. "Antiproliferative Activity Predictor: A New Reliable In Silico Tool for Drug Response Prediction against NCI60 Panel" International Journal of Molecular Sciences 23, no. 22: 14374. https://doi.org/10.3390/ijms232214374

APA Style

Martorana, A., La Monica, G., Bono, A., Mannino, S., Buscemi, S., Palumbo Piccionello, A., Gentile, C., Lauria, A., & Peri, D. (2022). Antiproliferative Activity Predictor: A New Reliable In Silico Tool for Drug Response Prediction against NCI60 Panel. International Journal of Molecular Sciences, 23(22), 14374. https://doi.org/10.3390/ijms232214374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop