Next Article in Journal
Gastrointestinal Mechanisms Underlying the Cardiovascular Effect of Metformin
Next Article in Special Issue
Tree-Based QSAR Model for Drug Repurposing in the Discovery of New Antibacterial Compounds against Escherichia coli
Previous Article in Journal
In Vitro and Ex Vivo Evaluation of Penetratin as a Non-invasive Permeation Enhancer in the Penetration of Salmon Calcitonin through TR146 Buccal Cells and Porcine Buccal Tissues
Previous Article in Special Issue
The Treatment of Impaired Wound Healing in Diabetes: Looking among Old Drugs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Perspective

A Multi-Objective Approach for Anti-Osteosarcoma Cancer Agents Discovery through Drug Repurposing

1
Grupo de Bio-Quimioinformática, Universidad de Las Américas, Quito 170125, Ecuador
2
Carrera de Enfermería, Facultad de Ciencias de la Salud, Universidad de Las Américas, Quito 170125, Ecuador
3
Department of Computer Science and Information Technologies, Faculty of Computer Science, University of A Coruña, CITIC, Campus Elviña s/n, 15071 A Coruña, Spain
4
Centro de Investigación Genética y Genómica, Facultad de Ciencias de la Salud Eugenio Espejo, Universidad UTE, Quito 170129, Ecuador
5
Latin American Network for Implementation and Validation of Clinical Pharmacogenomics Guidelines (RELIVAF-CYTED), 28029 Madrid, Spain
6
Laboratorio de Biología Molecular, Subproceso de Anatomía Patológica, Hospital de Especialidades Eugenio Espejo, Quito 170403, Ecuador
7
Department of Organic and Inorganic Chemistry, and Basque Center for Biophysics CSIC-UPV/EHU, University of the Basque Country UPV/EHU, 48940 Leioa, Spain
8
IKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Spain
9
Biomedical Research Institute of A Coruña (INIBIC), University Hospital Complex of A Coruña (CHUAC), 15006 A Coruña, Spain
10
Escuela de Ciencias Físicas y Matemáticas, Universidad de Las Américas, Quito 170125, Ecuador
11
Facultad de Ingeniería y Ciencias Agropecuarias, Universidad de Las Américas, Quito 170125, Ecuador
*
Authors to whom correspondence should be addressed.
Pharmaceuticals 2020, 13(11), 409; https://0-doi-org.brum.beds.ac.uk/10.3390/ph13110409
Submission received: 2 October 2020 / Revised: 11 November 2020 / Accepted: 12 November 2020 / Published: 22 November 2020
(This article belongs to the Collection Old Pharmaceuticals with New Applications)

Abstract

:
Osteosarcoma is the most common type of primary malignant bone tumor. Although nowadays 5-year survival rates can reach up to 60–70%, acute complications and late effects of osteosarcoma therapy are two of the limiting factors in treatments. We developed a multi-objective algorithm for the repurposing of new anti-osteosarcoma drugs, based on the modeling of molecules with described activity for HOS, MG63, SAOS2, and U2OS cell lines in the ChEMBL database. Several predictive models were obtained for each cell line and those with accuracy greater than 0.8 were integrated into a desirability function for the final multi-objective model. An exhaustive exploration of model combinations was carried out to obtain the best multi-objective model in virtual screening. For the top 1% of the screened list, the final model showed a BEDROC = 0.562, EF = 27.6, and AUC = 0.653. The repositioning was performed on 2218 molecules described in DrugBank. Within the top-ranked drugs, we found: temsirolimus, paclitaxel, sirolimus, everolimus, and cabazitaxel, which are antineoplastic drugs described in clinical trials for cancer in general. Interestingly, we found several broad-spectrum antibiotics and antiretroviral agents. This powerful model predicts several drugs that should be studied in depth to find new chemotherapy regimens and to propose new strategies for osteosarcoma treatment.

1. Introduction

Osteosarcoma (OS) is the most common primary bone tumor in children, adolescents and young adults, representing approximately 3.5% of all childhood cancers and 56% of malignant bone tumors in children. Its incidence rate ranges between 1 and 5 cases per million people and it is usually diagnosed in patients who are 10 to 19 years old. OS follows a bimodal distribution, with an initial peak in late adolescence and young adulthood and a second peak in old age [1].
The management of patients diagnosed with OS has not changed in recent decades. Current systemic OS first-line therapy includes cycles of cisplatin, doxorubicin, and high-dose methotrexate (MAP). Second-line therapy can integrate some tyrosine kinase inhibitors such as sorafenib and everolimus, plus antineoplastic agents like etoposide, topotecan and cyclophosphamide [2]. Neoadjuvant chemotherapy is generally administered for a period of 10 weeks, followed by the surgical resection of the compromised tumor area and radiotherapy. If 90% or more of the tumor area shows necrosis, additional cycles of postoperative therapy are applied to reject micrometastasis [3,4].
The prognosis of this disease is highly variable, possibly due to its high rate of tumor mutations, which leads to widespread dysregulation in cell signaling pathways and genomic instability [5]. Patients with localized disease show a 5-year survival rate of 65 to 70%, while for those who develop metastases the rate drops to 19–30% [6]. These metastatic events involve the lung parenchyma in 75% of the cases and distant skeletal sites [7,8], hindering treatment efficacy [9]. In this scenario, current therapy shows little response sensitivity and the survival rate decreases considerably.
Despite current chemotherapy regimens being the most effective strategy for OS treatment, patients’ sensitivity to these agents regarding the toxic side-effects and antitumor effects varies considerably [10,11]. Several clinical trials have developed experimental designs to improve survival rates by testing dose intensification, and also adding or combining various chemotherapeutic agents [12]. There is a dose effect on treatment response, but several studies have shown that high-dose chemotherapy may not increase survival rates any more than less toxic moderate doses [13]. Due to the lack of tumor specificity or metastasis events or the complex etiology of these bone tumors, the anti-OS compounds currently used have a narrow therapeutic index and no increase in survival rates have been achieved in the last three decades [14], thus the therapeutic strategies need to be optimized.
The development and validation of novel therapeutic compounds is a time-consuming and labor-intensive process. Drug repositioning, which explores potential novel uses for known molecules based on prediction algorithms has become an effective and innovative approach [15,16,17]. One approach is based on multi-objective computational models where the repositioning process is addressed from a set of potentially desirable solutions. To do this, some computational techniques have been applied and these include Quantitative Structure-Activity Relationship (QSAR) and Ligand-Based Virtual Screening, which aid the identification of hit structures [18,19,20,21]. These QSAR studies are used to perform virtual drug screening that has been integrated into the drug discovery pipeline and could save both time and money, especially in the early phase of drug discovery [22,23].
In this sense, the application of these models is of high interest to researchers specializing in cancer. Several studies have focused on the description of new therapeutic agents, especially for treating carcinomas [24,25,26,27,28,29,30,31,32]. However, very few have concentrated on tumors of mesenchymal origin [33,34]. Thus, we developed a multi-objective model for the prediction of drugs with potential biological activity towards OS, one of the most prevalent cancers in pediatric populations where current chemotherapy treatments have not varied in the last decades.

2. Results

2.1. Datasets and Molecular Descriptors

The ChEMBL database reports a total of 1250 compounds with biological activity for HOS, MG63, U2OS, and SAOS2 cell lines. Of these, 1036 shows complete information on their biological activities evaluated by IC50, GI50, and C50 assays (Table S1). Before constructing the prediction algorithms for each cell line, we inspected all those compounds reported in the DrugBank and separated them from this list for later use in virtual screening (VS).
Of the 1036 compounds, 28 drugs are reported in DrugBank for the HOS cell line, 30 for MG63, 31 for SAOS2, and 32 for U2OS. In this way, we removed these 121 drugs from the 1036, and the prediction models were built on the remaining 915 compounds. Thus, the calculation of the molecular descriptors (MDs) was performed on 277 compounds described for HOS, 124 for MG63, 173 for U2OS and 341 for SAOS2 (Figure 1A) and we obtained 500 ISIDA variables for each cell line (Table S2).
Since inactive compounds are described in greater numbers than active ones in all cell lines, we evaluated the chemical diversity in the inactive series by applying a hierarchical clustering. Thus, we calculated the degree of similarity in the inactive compounds to balance the data through stratified random sampling instead of a solely random partitioning. As a result, hierarchical representations were generated in which the clusters at each level of the hierarchy were created by merging clusters at the next level down [35]. In our case, we chose a strict cut-off to show all the possible groups and to make sure we had a wide representation of the chemical space within the data in the inactive series. From this, we identified four clusters in the HOS cell line, six in MG63, 14 in SAOS2 and 15 in U2OS (Figure 1B).
Therefore, we separated 24 inactive compounds for all cell lines (24 in each) reported in the DrugBank, as mentioned above, and then 7, 14, 19 and 29 compounds from each list. The balanced datasets resulted in 246 molecules for HOS, 86 for MG63, 130 for SAOS2, and 288 for U2OS, with a ratio of 1:1 between active and inactive compounds.

2.2. Construction of Models

The prediction algorithms used were: support vector machine (SVM), random forest (RF), neural networks (NN), decision tree (DTREE), k-Nearest Neighbors (KNN), and a scalable end-to-end tree boosting system (XGBoost) [36].
As seen in Figure 2, each of the six trained models demonstrated different performance metrics on the external data. HOS had similar achievements in the six learning techniques, but only the SVM, RF and XGBoost models showed optimal accuracy (AC) for subsequent assembly (0.836, 0.828 and 0.833 respectively). For MG63, the best strategies were SVM (0.882) and KNN (0.833). Prediction values when using RF, NN, XGBoost DTREE were less than 80%, and the true positive rate was lower than 0.7.
SAOS2 modeling resulted in only one algorithm with AC greater than 0.8, namely KNN (0.833). Lastly, the prediction models for U2OS with ACs > 0.8 were SVM, RF, NN, and KNN (0.886, 0.857, 0.812, and 0.857 respectively), all of whose SN and SP rates were higher than 0.8. XGBoost and DTREE were not taken into account for the assembly.
It is interesting to note that the feature selection by genetic algorithm (GA) reduced from 500 to 19–101, depending on each dataset. This strategy allowed us to generate the models described above and obtain the desired performance measurements for the final model (Table S3).

2.3. Multi-Objective Model Assessment and Virtual Screening

The AC, SN and specificity (SP) evaluated the performance of a model based on its training and data described as external, but these metrics do not always describe a desirable recovery rate at the time of performing screening for drug repositioning [37]. Therefore, the multi-objective model was evaluated in a virtual screening setting, where we mainly took into account the Area Under the Accumulation Curve (AUC), the Boltzmann-Enhanced Discrimination of ROC (BEDROC) and recovery efficiency (EF) at 1% of the screened list. The VS was developed on a dataset of 772 compounds. Of these, 14 corresponded to drugs used for previously described OS treatment, 653 were decoy molecules calculated from these 14 drugs, and 105 compounds were described as inactive and removed during the data balancing process.
We used all base-models with an AC greater than 0.8 in previously described external validation, and tested all possible combinations. Based on the VS results, the best multi-objective model was made up of the desirability values of the algorithms HOS-SVM, HOS-RF, MG63-SVM, SAOS2-KNN, U2OS-NN, and U2OS-KNN.
As seen in Figure 3A, our strategy generated a prediction method with better early recognition rates than the individual models. BEDROC is a metric that assigns more weight to early ranked molecules than late ranked molecules, therefore the initial enrichment was weighted. This enrichment was higher in our multi-objective model (BEDROC = 0.562) when calculated with an α = 160.9. This means that our algorithm turned out to be the best strategy for recognizing “active” anti-sarcoma molecules in 1% of the list of therapeutic drugs for OS. Likewise, the EF value was higher in our algorithm when analyzing the recovery rate at 0.01. The EF values calculated for a recovery efficiency of 1% in the base models HOS-RF, SAOS2-KNN, and HOS-SVM were 20.68, 13.8, and 20.68, respectively, while in our method they reached 27.57. This indicates that with our protocol, it is possible to retrieve almost 27 times more the number of multi-targeted compounds in the first 1% of the ranked list than what is expected from a uniform distribution of the active ones in the virtual screening database.
When analyzing the active retrieved fractions of all the models (Figure 3B), one notices that all protocols have similar recovery rates within 30% of the data screened. However, a closer inspection of the screened data shows that the multi-objective has the highest recovery rate of anti-OS compounds at 1% (or less) of the data screened (Figure 3C). This algorithm recognized four compounds within the first six positions. These retrieve rates suggest that a strategy made up of several methods is capable of predicting those molecules described as active within 1% of a screened list. In a drug repositioning scenario, this is important since the compounds ranked in the first positions have a high probability of presenting biological activity in vitro. Using our prediction algorithm, a prediction rate of 59.9% would be expected in 1% of screening for drugs with anti-OS activity.

2.4. Analysis of Repurposed Drugs

The screening weighted four principal drug classes in the first 1% of the screened list: anti-infectives for systemic use (antimycobacterials, macrolides, protease inhibitors and tetracyclines; which represent 55%); antineoplastic/immunomodulating agents (immunosuppressants, protein kinase inhibitors and taxanes; 32%); dermatological/immunosuppressant (agents for dermatitis, excluding corticosteroids; 4%); and antiparasitics (an antinematodal agent and a broad-spectrum endectocide; 9%). The first two groups represent more than 85% of all repositioned drugs (Figure 4A). All the desirability values for each repositioned drug are detailed in Table S4.
The action mechanism of antineoplastic and immunomodulating agent mainly inhibits the mTOR pathway and microtubules polymerization. Among these, temsirolimus, paclitaxel, sirolimus (rapamicine), everolimus, cabazitaxel and docetaxel are ranked at the top. On the other hand, broad-spectrum antibacterials described as drugs that bind to the 30S/50s subunit of bacterial ribosome, HIV-1 protease inhibitors and antimycobacterials, which inhibit DNA-dependent RNA bacterial polymerase, were weighted in the screening. We also found two molecules used for the treatment of HIV, described as inhibitors of HIV-1 protease (tipranavir and fosamprenavir) (Figure 4B). It is interesting to note that several repositioned drugs have been found within clinical trials for cancer patients. Out of the antineoplastic and immunomodulating agents, only cabazitaxel has not yet been studied in trials related to bone sarcomas. Moreover, broad-spectrum antibacterial compounds such as clarithromycin, erythromycin, doxycycline and tetracycline are top-ranked drugs that are registered in clinical trials for carcinomas.

3. Discussion

Several reports have shown that multi-objective models have a better prediction rate during screening time since they approach the problem with a particular perspective from a set of potentially desirable solutions [18,21,38,39]. In our case, each of these possible desirable solutions was made up of each algorithm constructed from the described compounds with activity for the OS cell lines HOS, MG63, SAOS2 and U2OS. One of the major outcomes of this study is the improvement of the AUC and BEDROC values obtained in the VS, especially the EF of the multi-objective model when comparing with the performance of the base models (Figure 3A). This suggests that our algorithm improves the recognition rate of molecules described as therapeutic for OS treatment, especially within 1% of the data screened (Figure 3B). Specifically, the EF obtained indicates that it is possible to retrieve in the first 1% of the ranked list almost 27 times more multi-targeted compounds than what is expected from a uniform distribution of the active ones in the virtual screening database, something that is not obtained from the algorithms generated by each cell line.
Drug repositioning is an effective strategy for finding novel drug-disease relationships for existing molecules. The development of these strategies has gained considerable interest in recent years compared with de novo drug discovery pipeline, which demands more research time and experimental hours in the case of new drug development, and requires a greater financial investment. On the other hand, the use of already proven drugs is a highly efficient, low-cost and low-risk strategy since screening is carried out on molecules that have passed all clinical safety tests at Phase I, Phase II, and Phase III [40,41]. We used our multi-objective approach in order to propose new agents with chemotherapeutic activity for osteosarcoma treatment. Given the high recovery rate of active compounds obtained in our model (EF 0.01 = 27,571), we considered the first 22 highest-ranking compounds belonging to 1% of the 2218 approved FDA drugs reported in the DrugBank.
Of these 22 drugs, 13 (59.1%) are enrolled in clinical trials for cancer patients (reviewed at https://clinicaltrials.gov/) (Table S5): temsirolimus, paclitaxel, sirolimus/rapamycin, everolimus, cabazitaxel, docetaxel, rifampicin, tacrolimus, clarithromycin, erylinethromycin, doxycycline, tetracycline and minocycin. Interestingly, only five of these drugs are included in trials of patients with OS: temsirolimus, paclitaxel, sirolimus/rapamicin, everolimus, and docetaxel. The remaining 10 drugs (ivermectin, pimecrolimus, roxithromycin, tipranavir, oxytetracycline, fosamprenavir, rifapentine, troleandomycin and moxidectin) are not registered in any clinical trial for cancer patients, however, their action mechanisms are similar to various chemotherapeutic agents used in oncology practice.
Cancer cells are characterized by unregulated proliferation, which leads to cellular undifferentiation and disruption on the function of tissues. Cell proliferation can be caused by a checkpoint failure in cell cycle or a disruption in the cell death pathway. In this sense, any agent that affects the metabolism of cancer cells by reducing or inhibiting cell proliferation and promotes apoptosis is a potential target for cancer treatment [42]. Several agents used as first-line treatment for OS, such as methotrexate, doxorubicin, etoposide, cisplatin and ifosfamide, induce a disruption in these cellular functions, either by interrupting nucleotide synthesis, by DNA synthesis by inhibiting topoisomerase II, or by binding to double-strand DNA to promote apoptosis. On the other hand, several second-line drugs act on mTOR, a pathway considered pathogenic within the development and progression of OS [43,44,45], and on the formation of microtubules, inhibiting the progression from the G1 to the S phase of the cell cycle. In the top six positions of our screening, we found chemotherapeutic drugs described as therapeutic agents for various types of cancer (temsirolimus, paclitaxel, sirolimus, everolimus, cabazitaxel and docetaxel). Indeed, these compounds belong to one of the four principal drug classes found in our repositioning called antineoplastic and immune-modulatory agents (Figure 4A). Their action mechanism resembles those previously described as second-line drugs, which mainly inhibit mTOR and interfere with the microtubule depolymerization (Figure 4B). Interestingly, cabazitaxel is the only one of these six top-ranked compounds that is not reported in clinical trials of OS patients. This molecule is a semi-synthetic derivative of a natural taxoid that considerably increases overall survival versus mitoxantrone after prior docetaxel treatment in patients with metastatic castration-resistant prostate cancer [46,47,48]. Cabazitaxel induces cell cycle arresting by interacting with the microtubule depolymerization by what is defined as a microtubule destabilizing agent. These types of agents show high antineoplastic activity and have been reported in previous studies into drug repositioning [49]. Although they commonly used in pediatric oncology [50], the microtubule-stabilizing taxanes are not often used to treat childhood cancers due to limited activity, even if safety is observed [51]. In this sense, cabazitaxel can be an important therapeutic agent for the treatment of OS, especially in patients who can progress onto it after docetaxel.
It is interesting to note that 54.5% of the total predicted compounds (12 out of 22) are classified as anti-infectives for systemic use. More specifically, taking into account the Anatomical Therapeutic Chemical (ATC) classification system, our protocol weighted several macrolides (roxithromycin, clarithromycin, erythromycin and troleandomycin), tetracyclines (oxytetracycline, doxycycline, tetracycline and minocycline), protease inhibitors (tipranavir and fosamprenavir) and antimycobacterial (rifampicin and rifapentine) as possible anti-OS agents (Figure 4B). On the one hand, prior studies into cancer therapy have noted the importance of macrolide and tetracycline compounds in cancer treatment [52,53]. Some authors have suggested that these groups of compounds inhibit the action of matrix metalloproteinases (MMPs) in order to reduce the degree of tumor invasion and metastases [54]. Others have observed that these drugs act on mitochondrial biogenesis [55,56], disrupting this process and thus increasing the effectiveness of chemotherapy or radiotherapy on tumor cells. On the other hand, the therapeutic action of HIV-protease inhibitors for the treatment of cancer has been reported. Although these molecules are not expected to cross-react with human peptides, preclinical data suggest that their antitumor activity may be linked in part to the inhibition of endopeptidases, such as metalloproteases and proteasomes [57]. Of our repositioned drugs, clarithromycin, erythromycin and doxycycline are currently under study as possible therapeutic agents for leukemia, colorectal, prostate and lung cancer, among others [58,59,60,61], and are involved in clinical trials of cancer patients (Table S5). Based on our findings, these agents could demonstrate antitumor activity in bone tumors.
These results may be promising for future preclinical and clinical studies. The lack of therapeutic options for OS should be the basis of searches for new agents as potential treatments. The discovery of molecular targets in OS will be part of the development of new molecules that could give these patients more options.

4. Materials and Methods

4.1. Preprocessing Datasets and Molecular Descriptors

Prediction models were developed from compounds described in the Chemical database (Version 25) of the European Molecular Biology Laboratory (ChEMBL) [62,63] with biological activity against the OS cell lines HOS (ChEMBL614736), MG63 (ChEMBL614347), SAOS (ChEMBL614894) and U2OS (ChEMBL615023). We considered all standard values evaluated by IC50 (half-maximal inhibitory concentration), GI50 (percentage of cell grow inhibition at a fixed concentration), and EC50 (a concentration that inhibited half the cell culture growth), and from these scores, we defined a class for each compound.
Compounds with standard values > 10 µM were classified as inactive (0), and those with values < 10 µM as active (1). In those drugs where two or more assays are reported, and the standard values classify these compounds in different classes, the final criteria were assigned by most of the set of classes obtained. If more than 75% of the tests obtained classify a compound in the same class, this drug was included in the study, otherwise it was rejected. On the other hand, compounds that did not show information about their biological activity, inconclusive data about their activity, or incomplete information regarding ChEMBL ID or canonical SMILES were removed from the analysis.
We used the ChemAxon’s JChem for Excel (18.8.0.253) [64] software to code the chemical structures in SMILES format. This information was converted to SD files (SDFs) and the structure of each compound was standardized using ChemAxon’s Standardizer [65]. Explicit hydrogen atoms were removed. Then we normalized specific chemotypes, such as nitro to one unique representation, the rings aromatization, the curation of tautomeric forms, the striping of salts and small fragments. Furthermore, all duplicate structures were identified using the EdiSDF tool within the ISIDA/QSPR package and subsequently withdrawn from the list [66].
Two-dimensional molecular descriptors were computed with ISIDA Fragmentor 2017 [67,68]. The types of descriptors calculated were: Sequences of atoms and bonds; Atom-centered fragments based on sequences of atoms and bonds; Atom-centered fragments based on sequences of atoms and bonds of fixed length; and Triplets. For these calculations, the minimum and maximum length of fragments as sequences were set to 2 and 8, respectively. Molecular descriptors were calculated separately for each dataset.
The computed descriptors for each dataset were first filtered to remove those present in less than 1% of the compounds. Next, the Minimal Redundancy Maximal Relevance (mRMR) algorithm [69] was employed to keep the top 500 features in each dataset. For mRMR, the Mutual Information Quotient (MIQ) score was used as a features-ranking metric. This subset of 500 selected molecular descriptors was employed for QSAR modeling.

4.2. Machine Learning Models and Quality Evaluation

In each cell line under study, the amount of active and inactive compounds varied considerably. This unbalance in the dataset is not desirable for modeling. In order to balance the classes, the following procedure was executed in all cell lines. (1) Using the previously computed molecular descriptors, we carried out a hierarchical clustering. We applied the interval measure, the Euclidean distance and Ward’s method for clustering, both in active and inactive compounds for each cell line. The IBM SPSS Statistics software v.25 (IBM Corp., Armonk, NY, USA) was employed to generate dendrograms and define all the clusters within the data. (2) Once the number of clusters had been defined, we continued with a random stratify extraction of the same amount of compound in both classes. This procedure had previously been used by other authors in order to obtain a balanced data representative of the chemical diversity space [70]. The training set consisted of 75% of randomly chosen compounds from the balanced dataset and the remaining percentage was utilized as external data. The external dataset was used to evaluate the model performance metrics.
To obtain each model, we applied genetic algorithms as feature selection by considering an initial population of 50 chromosomes and 30 generations. For validation of the fitness function in GA, we performed a cross-validation strategy using the average balanced classification rate (BRC) across 100 random splits (bootstrap sampling). This means that in each generation, 100 models were evaluated and the average AC was extracted. The models used together with genetic algorithm were: the support vector machine, random forest, neural networks, decision tree, k-Nearest Neighbors, and a scalable end-to-end tree boosting system [36]. The SVM kernel was fixed to RBF. For performance metrics of models, we calculated the total accuracy (AC), sensitivity (SN), specificity (SP) and the balanced classification rate (BCR) as follows:
A C = N u m b e r   o f   c o r r e c t l y   c l a s s i f i e d   c o m p o u n d s T o t a l   n u m b e r   o f   c o m p o u n d s
S N = N u m b e r   o f   c o r r e c t l y   c l a s s i f i e d   a c t i v e   c o m p o u n d s T o t a l   n u m b e r   a c t i v e   c o m p o u n d s
S P = N u m b e r   o f   c o r r e c t l y   c l a s s i f i e d   i n a c t i v e   c o m p o u n d s T o t a l   n u m b e r   i n a c t i v e   c o m p o u n d s
B C R = S N + S P 2 × ( 1 | S N S P | )

4.3. Multi-Objective Model Assembly and Virtual Screening

The construction of the multi-objective model was performed by computing the global desirability as:
D 1 = ( d ( y 1 ) d ( y 2 ) ,   ,   d ( y k ) ) 1 / k
where y k corresponds with the desirability scores of each cell line (k = 1, …, 4). For each of the cell lines, several possible models are available. The resulting prediction of each model for a given compound resulted in a score linked to the class membership- either a prediction score for an active class, and/or score for an inactive class. In all cases, we established the score for which the compound was active against the cell line. However, the calculated class membership for some machine learning algorithm occurred in the range of positives and negatives, where active compounds showed positive values and vice versa. For all the cases, the geometrical mean of all scores of the compound to be active in a particular cell line (or to be 0–1, normalized transformation) was used as a desirability score for each model ( y k ). Since there are several possible combinations, we performed an exhaustive exploration to obtain the best possible model. Hence, we explored the combination of all possible models in the computation of each d ( y k ) and consequently D1 in order to obtain the best performance in early recognition metrics for virtual screening.
For VS, we developed a dataset with those antitumor compounds used in the current management of osteosarcoma, not included in either the training or external sets for any cell line and with compounds validated in clinical studies for OS, published on the US government’s Clinical Trials website (https://clinicaltrials.gov/). As first- and second-line therapy drugs, we included [4,12,71,72]: doxorubicin (ChEMBL53463), methotrexate (ChEMBL34259), ifosfamide (ChEMBL1024), etoposide (ChEMBL44657), sorafenib (ChEMBL1336), cyclophosphamide (ChEMBL88), docetaxel (ChEMBL92), gemcitabine (ChEMBL888), dactinomycin (ChEMBL1554) and vincristine (ChEMBL90555). Additionally, we incorporated as validated drugs in clinical trials: temsirolimus (ChEMBL1201182) [73,74], ridaforolimus (ChEMBL2103839) [75,76], sirolimus (ChEMBL413) [77] and pazopanib (ChEMBL477772) [78,79].
As inactive compounds for screening, we considered those molecules withdrawn in the data balancing process (previously described, thus not employed for model training and selection), and common ChEMBL compounds for the four cell lines that showed no biological activity (standard values > 10 µM). Additionally, we generated Decoy molecules based on the selected active compounds by employing the DUD-E server 5 [80]. We incorporated around 50 inactive molecules for each active compound, which is the proportion used in the DUD-E database that is widely employed to validate virtual screening workflows [80].
The performance of our models within this VS scenario was evaluated by computing the AUAC, BEDROC and EF [81,82]:
A U C = 1 1 n i = 1 n x i
E F = i = 1 n δ i χ n , where   δ i = { 1   r i χ N 0   r i > χ N }
B E D R O C = R I E R I E m i n R I E m a x R I E m i n , with   R I E m i n = 1 e α R a R a ( 1 e α ) , R I E m a x = 1 e α R a R a ( 1 e α )   and R I E = 1 n i = 1 n e α x i 1 N ( 1 e α e α / N 1 )
In the above equations, n represents the number of active compounds, N the total number of compounds in the dataset, xi the relative ranking of active compound i in the ranked list, χ the fraction of data for which EF will be computed, Ra the rate of active compounds in the dataset (n/N), and α is the α parameter which ensures that active compounds ranked at the beginning of the ordered list result in higher weights than those at the tail. The α parameter is computed using the following equation:
θ ( 1 e α ) 1 + e α z = 0
where z represents the fraction of the ranked list at which enrichment is important and θ is the expected contribution of the enrichment at this z% fraction to the overall enrichment.

5. Conclusions

In conclusion, this study presents a multi-objective prediction algorithm developed from compounds described with biological activity for the osteosarcoma cell lines HOS, MG63, SAOS2 and U2OS. The performance of this multi-objective model considerably improves the recognition rate in a virtual screening scenario, developed on drugs used as first- and second-line treatment for OS. Specifically, a high level of performance was observed for the recognition of molecules with biological activity within 1%. Using this ML algorithm on 2218 compounds described in the DrugBank, we found several antineoplastic agents currently being studied in clinical trials for the treatment of OS. Interestingly, Cabazitaxel is a compound with chemotherapeutic activity that is being studied in several clinical trials for different types of carcinomas and not in sarcomas, therefore it can be taken into account for clinical validations in patients with OS. Furthermore, several broad-spectrum antibiotics, for instance clarithromycin, erythromycin and doxycycline, were top-ranked drugs in our screening. These compounds have already been studied in various types of carcinomas, so they comprise an interesting group of drugs for developing therapeutic validation studies in bone cancers. One of the main limitations was the lack of experimental validation of the drugs proposed for repositioning because this was an initial study. Although it is true that several compounds have already been studied in clinical trials, the validation process of their biological activities on bone tumor cells is an indispensable step in order to proceed with a validation strategy in patients, hence this will be the next procedure in further studies.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-8247/13/11/409/s1. Table S1. ChEMBL dataset of compounds with biological activity for HOS, MG63, U2OS and SAOS2 cell lines. Table S2. ISIDA Molecular descriptors calculated for each cell line’s dataset. Table S3. Performance on machine learning models for HOS, MG63, SAOS2 and U2OS cell lines. Table S4. Ranking and desirability values of drugs repositioned by the multi-objective model. Table S5. Clinical trials reported for the 22 top-ranked compounds repositioned by the multi-objective model.

Author Contributions

E.T., A.C.-A. and Y.P.-C. conceived the project and wrote the manuscript. E.T. and Y.P.-C. designed the algorithms. A.C.-A. and A.L.-C. implemented the algorithm and performed the data analysis. G.J.-K. made substantial contributions to the discussion section. C.R.M., H.G.-D. and A.P. helped with study design and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universidad de Las Américas, Quito, Ecuador, grant number ENF.RCA.18.01, by Ministry of Competitiveness and Economy (CTQ2016-74881-P), Ministry of Science and Innovation (PID2019-104148GB-I00), and Basque Government (IT1045-16)-2016–2021.

Acknowledgments

This work was supported by Universidad de Las Américas (Quito, Ecuador) and the Competitive Reference Groups (Ref. ED431C 2018/49) in Galicia, Spain.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Ottaviani, G.; Jaffe, N. The epidemiology of osteosarcoma. Cancer Treat. Res. 2009, 152, 3–13. [Google Scholar]
  2. National Comprehensive Cancer Network. NCCN Clinical Practice Guidelines in Oncology (NCC Guidelines)—Bone Cancer. Available online: https://www.nccn.org/ (accessed on 4 March 2020).
  3. Misaghi, A.; Goldin, A.; Awad, M.; Kulidjian, A.A. Osteosarcoma: A comprehensive review. SICOT J. 2018, 4, 12. [Google Scholar] [CrossRef] [Green Version]
  4. Biermann, J.S.; Chow, W.; Reed, D.R.; Lucas, D.; Adkins, D.R.; Agulnik, M.; Benjamin, R.S.; Brigman, B.; Budd, G.T.; Curry, W.T.; et al. NCCN Guidelines Insights: Bone Cancer, Version 2.2017. J. Natl. Compr. Cancer Netw. 2017, 15, 155–167. [Google Scholar] [CrossRef] [Green Version]
  5. de Azevedo, J.W.V.; de Medeiros Fernandes, T.A.A.; Fernandes, J.V., Jr.; de Azevedo, J.C.V.; Lanza, D.C.F.; Bezerra, C.M.; Andrade, V.S.; de Araujo, J.M.G.; Fernandes, J.V. Biology and pathogenesis of human osteosarcoma. Oncol. Lett. 2020, 19, 1099–1116. [Google Scholar]
  6. Xin, S.; Wei, G. Prognostic factors in osteosarcoma: A study level meta-analysis and systematic review of current practice. J. Bone Oncol. 2020, 21, 100281. [Google Scholar] [CrossRef]
  7. Marko, T.A.; Diessner, B.J.; Spector, L.G. Prevalence of Metastasis at Diagnosis of Osteosarcoma: An International Comparison. Pediatr. Blood Cancer 2016, 63, 1006–1011. [Google Scholar] [CrossRef] [Green Version]
  8. Duchman, K.R.; Gao, Y.; Miller, B.J. Prognostic factors for survival in patients with high-grade osteosarcoma using the Surveillance, Epidemiology, and End Results (SEER) Program database. Cancer Epidemiol. 2015, 39, 593–599. [Google Scholar] [CrossRef]
  9. Song, K.; Song, J.; Lin, K.; Chen, F.; Ma, X.; Jiang, J.; Li, F. Survival analysis of patients with metastatic osteosarcoma: A Surveillance, Epidemiology, and End Results population-based study. Int. Orthop. 2019, 43, 1983–1991. [Google Scholar] [CrossRef]
  10. Taran, S.J.; Taran, R.; Malipatil, N.B. Pediatric Osteosarcoma: An Updated Review. Indian J. Med. Paediatr. Oncol. 2017, 38, 33–43. [Google Scholar] [CrossRef] [Green Version]
  11. Vos, H.I.; Coenen, M.J.; Guchelaar, H.J.; Te Loo, D.M. The role of pharmacogenetics in the treatment of osteosarcoma. Drug Discov. Today 2016, 21, 1775–1786. [Google Scholar] [CrossRef]
  12. Durfee, R.A.; Mohammed, M.; Luu, H.H. Review of Osteosarcoma and Current Management. Rheumatol. Ther. 2016, 3, 221–243. [Google Scholar] [CrossRef] [Green Version]
  13. Omer, N.; Le Deley, M.C.; Piperno-Neumann, S.; Marec-Berard, P.; Italiano, A.; Corradini, N.; Bellera, C.; Brugieres, L.; Gaspar, N. Phase-II trials in osteosarcoma recurrences: A systematic review of past experience. Eur. J. Cancer 2017, 75, 98–108. [Google Scholar] [CrossRef]
  14. Harrison, D.J.; Geller, D.S.; Gill, J.D.; Lewis, V.O.; Gorlick, R. Current and future therapeutic approaches for osteosarcoma. Expert Rev. Anticancer Ther. 2018, 18, 39–50. [Google Scholar] [CrossRef]
  15. Lo, Y.C.; Rensi, S.E.; Torng, W.; Altman, R.B. Machine learning in chemoinformatics and drug discovery. Drug Discov. Today 2018, 23, 1538–1546. [Google Scholar] [CrossRef]
  16. Brown, A.S.; Patel, C.J. A review of validation strategies for computational drug repositioning. Brief. Bioinform. 2018, 19, 174–177. [Google Scholar] [CrossRef]
  17. Pushpakom, S.; Iorio, F.; Eyers, P.A.; Escott, K.J.; Hopper, S.; Wells, A.; Doig, A.; Guilliams, T.; Latimer, J.; McNamee, C.; et al. Drug repurposing: Progress, challenges and recommendations. Nat. Rev. Drug Discov. 2019, 18, 41–58. [Google Scholar] [CrossRef]
  18. Ma, X.H.; Shi, Z.; Tan, C.; Jiang, Y.; Go, M.L.; Low, B.C.; Chen, Y.Z. In-silico approaches to multi-target drug discovery: Computer aided multi-target drug design, multi-target virtual screening. Pharm. Res. 2010, 27, 739–749. [Google Scholar] [CrossRef]
  19. Huang, G.; Li, J.; Wang, P.; Li, W. A Review of Computational Drug Repositioning Approaches. Comb. Chem. High Throughput Screen 2017, 20, 831–838. [Google Scholar] [CrossRef]
  20. Park, K. A review of computational drug repurposing. Transl. Clin. Pharmacol. 2019, 27, 59–63. [Google Scholar] [CrossRef] [Green Version]
  21. Cruz-Monteagudo, M.; Schurer, S.; Tejera, E.; Perez-Castillo, Y.; Medina-Franco, J.L.; Sanchez-Rodriguez, A.; Borges, F. Systemic QSAR and phenotypic virtual screening: Chasing butterflies in drug discovery. Drug Discov. Today 2017, 22, 994–1007. [Google Scholar] [CrossRef]
  22. Murphy, R.F. An active role for machine learning in drug development. Nat. Chem. Biol. 2011, 7, 327–330. [Google Scholar] [CrossRef] [Green Version]
  23. Ramsay, R.R.; Popovic-Nikolic, M.R.; Nikolic, K.; Uliassi, E.; Bolognesi, M.L. A perspective on multi-target drug discovery and design for complex diseases. Clin. Transl. Med. 2018, 7, 3. [Google Scholar] [CrossRef] [Green Version]
  24. Nagamalla, L.; Kumar, J.V.S. In silico screening of FDA approved drugs on AXL kinase and validation for breast cancer cell line. J. Biomol. Struct. Dyn. 2020. [Google Scholar] [CrossRef]
  25. Kumar, R.; Chaudhary, K.; Singla, D.; Gautam, A.; Raghava, G.P. Designing of promiscuous inhibitors against pancreatic cancer cell lines. Sci. Rep. 2014, 4, 4668. [Google Scholar] [CrossRef] [Green Version]
  26. Issa, N.T.; Stathias, V.; Schurer, S.; Dakshanamurthy, S. Machine and deep learning approaches for cancer drug repurposing. Semin. Cancer Biol. 2020. [Google Scholar] [CrossRef]
  27. Koudijs, K.K.M.; Terwisscha van Scheltinga, A.G.T.; Bohringer, S.; Schimmel, K.J.M.; Guchelaar, H.J. Personalised drug repositioning for Clear Cell Renal Cell Carcinoma using gene expression. Sci. Rep. 2018, 8, 5250. [Google Scholar] [CrossRef] [Green Version]
  28. Wei, G.G.; Gao, L.; Tang, Z.Y.; Lin, P.; Liang, L.B.; Zeng, J.J.; Chen, G.; Zhang, L.C. Drug repositioning in head and neck squamous cell carcinoma: An integrated pathway analysis based on connectivity map and differential gene expression. Pathol. Res. Pract. 2019, 215, 152378. [Google Scholar] [CrossRef]
  29. Singh, H.; Kumar, R.; Singh, S.; Chaudhary, K.; Gautam, A.; Raghava, G.P. Prediction of anticancer molecules using hybrid model developed on molecules screened against NCI-60 cancer cell lines. BMC Cancer 2016, 16, 77. [Google Scholar] [CrossRef] [Green Version]
  30. Speck-Planche, A.; Cordeiro, M. Fragment-based in silico modeling of multi-target inhibitors against breast cancer-related proteins. Mol. Divers. 2017, 21, 511–523. [Google Scholar] [CrossRef]
  31. Bediaga, H.; Arrasate, S.; Gonzalez-Diaz, H. PTML Combinatorial Model of ChEMBL Compounds Assays for Multiple Types of Cancer. ACS Comb. Sci 2018, 20, 621–632. [Google Scholar] [CrossRef]
  32. Lopez-Cortes, A.; Paz, Y.M.C.; Guerrero, S.; Cabrera-Andrade, A.; Barigye, S.J.; Munteanu, C.R.; Gonzalez-Diaz, H.; Pazos, A.; Perez-Castillo, Y.; Tejera, E. OncoOmics approaches to reveal essential genes in breast cancer: A panoramic view from pathogenesis to precision medicine. Sci. Rep. 2020, 10, 5285. [Google Scholar] [CrossRef] [Green Version]
  33. Li, X.; Yan, M.L.; Yu, Q. Identification of candidate drugs for the treatment of metastatic osteosarcoma through a subpathway analysis method. Oncol. Lett. 2017, 13, 4378–4384. [Google Scholar] [CrossRef] [Green Version]
  34. Speck-Planche, A.; Kleandrova, V.V.; Luan, F.; Cordeiro, M.N. Fragment-based QSAR model toward the selection of versatile anti-sarcoma leads. Eur. J. Med. Chem. 2011, 46, 5910–5916. [Google Scholar] [CrossRef]
  35. Chalise, P.; Koestler, D.C.; Bimali, M.; Yu, Q.; Fridley, B.L. Integrative clustering methods for high-dimensional molecular data. Transl. Cancer Res. 2014, 3, 202–216. [Google Scholar]
  36. Chen, T.; Guestrin, C. Xgboost: A scalable tree boosting system. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 785–794. [Google Scholar]
  37. Braga, R.C.; Alves, V.M.; Silva, A.C.; Nascimento, M.N.; Silva, F.C.; Liao, L.M.; Andrade, C.H. Virtual screening strategies in medicinal chemistry: The state of the art and current challenges. Curr. Top. Med. Chem. 2014, 14, 1899–1912. [Google Scholar] [CrossRef]
  38. Bolognesi, M.L.; Cavalli, A. Multitarget Drug Discovery and Polypharmacology. Chem. Med. Chem. 2016, 11, 1190–1192. [Google Scholar] [CrossRef] [Green Version]
  39. Perez-Castillo, Y.; Sanchez-Rodriguez, A.; Tejera, E.; Cruz-Monteagudo, M.; Borges, F.; Cordeiro, M.; Le-Thi-Thu, H.; Pham-The, H. A desirability-based multi objective approach for the virtual screening discovery of broad-spectrum anti-gastric cancer agents. PLoS ONE 2018, 13, e0192176. [Google Scholar] [CrossRef] [Green Version]
  40. Xue, H.; Li, J.; Xie, H.; Wang, Y. Review of Drug Repositioning Approaches and Resources. Int. J. Biol Sci 2018, 14, 1232–1244. [Google Scholar] [CrossRef] [Green Version]
  41. Langedijk, J.; Mantel-Teeuwisse, A.K.; Slijkerman, D.S.; Schutjens, M.H. Drug repositioning and repurposing: Terminology and definitions in literature. Drug Discov. Today 2015, 20, 1027–1034. [Google Scholar] [CrossRef]
  42. Parvathaneni, V.; Kulkarni, N.S.; Muth, A.; Gupta, V. Drug repurposing: A promising tool to accelerate the drug discovery process. Drug Discov. Today 2019, 24, 2076–2085. [Google Scholar] [CrossRef]
  43. Ding, L.; Congwei, L.; Bei, Q.; Tao, Y.; Ruiguo, W.; Heze, Y.; Bo, D.; Zhihong, L. mTOR: An attractive therapeutic target for osteosarcoma? Oncotarget 2016, 7, 50805–50813. [Google Scholar] [CrossRef] [Green Version]
  44. Bishop, M.W.; Janeway, K.A. Emerging concepts for PI3K/mTOR inhibition as a potential treatment for osteosarcoma. F1000Research 2016, 5, 1590. [Google Scholar] [CrossRef]
  45. Cabrera-Andrade, A.; Lopez-Cortes, A.; Jaramillo-Koupermann, G.; Paz, Y.M.C.; Perez-Castillo, Y.; Munteanu, C.R.; Gonzalez-Diaz, H.; Pazos, A.; Tejera, E. Gene Prioritization through Consensus Strategy, Enrichment Methodologies Analysis, and Networking for Osteosarcoma Pathogenesis. Int. J. Mol. Sci 2020, 21, 1053. [Google Scholar] [CrossRef] [Green Version]
  46. de Bono, J.S.; Oudard, S.; Ozguroglu, M.; Hansen, S.; Machiels, J.P.; Kocak, I.; Gravis, G.; Bodrogi, I.; Mackenzie, M.J.; Shen, L.; et al. Prednisone plus cabazitaxel or mitoxantrone for metastatic castration-resistant prostate cancer progressing after docetaxel treatment: A randomised open-label trial. Lancet 2010, 376, 1147–1154. [Google Scholar] [CrossRef]
  47. de Wit, R.; de Bono, J.; Sternberg, C.N.; Fizazi, K.; Tombal, B.; Wulfing, C.; Kramer, G.; Eymard, J.C.; Bamias, A.; Carles, J.; et al. Cabazitaxel versus Abiraterone or Enzalutamide in Metastatic Prostate Cancer. N. Engl. J. Med. 2019, 381, 2506–2518. [Google Scholar] [CrossRef]
  48. Oudard, S.; Fizazi, K.; Sengelov, L.; Daugaard, G.; Saad, F.; Hansen, S.; Hjalm-Eriksson, M.; Jassem, J.; Thiery-Vuillemin, A.; Caffo, O.; et al. Cabazitaxel Versus Docetaxel As First-Line Therapy for Patients With Metastatic Castration-Resistant Prostate Cancer: A Randomized Phase III Trial-FIRSTANA. J. Clin. Oncol. 2017, 35, 3189–3197. [Google Scholar] [CrossRef]
  49. Lo, Y.C.; Senese, S.; France, B.; Gholkar, A.A.; Damoiseaux, R.; Torres, J.Z. Computational Cell Cycle Profiling of Cancer Cells for Prioritizing FDA-Approved Drugs with Repurposing Potential. Sci. Rep. 2017, 7, 11261. [Google Scholar] [CrossRef] [Green Version]
  50. Reynolds, C.P.; Kang, M.H.; Maris, J.M.; Kolb, E.A.; Gorlick, R.; Wu, J.; Kurmasheva, R.T.; Houghton, P.J.; Smith, M.A. Initial testing (stage 1) of the anti-microtubule agents cabazitaxel and docetaxel, by the pediatric preclinical testing program. Pediatr. Blood Cancer 2015, 62, 1897–1905. [Google Scholar] [CrossRef] [Green Version]
  51. Amoroso, L.; Castel, V.; Bisogno, G.; Casanova, M.; Marquez-Vega, C.; Chisholm, J.C.; Doz, F.; Moreno, L.; Ruggiero, A.; Gerber, N.U.; et al. Phase II results from a phase I/II study to assess the safety and efficacy of weekly nab-paclitaxel in paediatric patients with recurrent or refractory solid tumours: A collaboration with the European Innovative Therapies for Children with Cancer Network. Eur. J. Cancer 2020, 135, 89–97. [Google Scholar] [CrossRef]
  52. Hussain, A.; Dar, M.S.; Bano, N.; Hossain, M.M.; Basit, R.; Bhat, A.Q.; Aga, M.A.; Ali, S.; Hassan, Q.P.; Dar, M.J. Identification of dinactin, a macrolide antibiotic, as a natural product-based small molecule targeting Wnt/beta-catenin signaling pathway in cancer cells. Cancer Chemother. Pharmacol. 2019, 84, 551–559. [Google Scholar] [CrossRef]
  53. Gupta, A.; Okesli-Armlovich, A.; Morgens, D.; Bassik, M.C.; Khosla, C. A genome-wide analysis of targets of macrolide antibiotics in mammalian cells. J. Biol. Chem. 2020, 295, 2057–2067. [Google Scholar] [CrossRef]
  54. Bahrami, F.; Morris, D.L.; Pourgholami, M.H. Tetracyclines: Drugs with huge therapeutic potential. Mini Rev. Med. Chem 2012, 12, 44–52. [Google Scholar] [CrossRef]
  55. Fiorillo, M.; Toth, F.; Sotgia, F.; Lisanti, M.P. Doxycycline, Azithromycin and Vitamin C (DAV): A potent combination therapy for targeting mitochondria and eradicating cancer stem cells (CSCs). Aging 2019, 11, 2202–2216. [Google Scholar] [CrossRef]
  56. Lamb, R.; Ozsvari, B.; Lisanti, C.L.; Tanowitz, H.B.; Howell, A.; Martinez-Outschoorn, U.E.; Sotgia, F.; Lisanti, M.P. Antibiotics that target mitochondria effectively eradicate cancer stem cells, across multiple tumor types: Treating cancer like an infectious disease. Oncotarget 2015, 6, 4569–4584. [Google Scholar] [CrossRef] [Green Version]
  57. Maksimovic-Ivanic, D.; Fagone, P.; McCubrey, J.; Bendtzen, K.; Mijatovic, S.; Nicoletti, F. HIV-protease inhibitors for the treatment of cancer: Repositioning HIV protease inhibitors while developing more potent NO-hybridized derivatives? Int. J. Cancer 2017, 140, 1713–1726. [Google Scholar] [CrossRef] [Green Version]
  58. Petroni, G.; Stefanini, M.; Pillozzi, S.; Crociani, O.; Becchetti, A.; Arcangeli, A. Data describing the effects of the Macrolide Antibiotic Clarithromycin on preclinical mouse models of Colorectal Cancer. Data Brief. 2019, 26, 104406. [Google Scholar] [CrossRef]
  59. Van Nuffel, A.M.; Sukhatme, V.; Pantziarka, P.; Meheus, L.; Sukhatme, V.P.; Bouche, G. Repurposing Drugs in Oncology (ReDO)-clarithromycin as an anti-cancer agent. Ecancermedicalscience 2015, 9, 513. [Google Scholar] [CrossRef] [Green Version]
  60. de Jong, J.; Hellemans, P.; De Wilde, S.; Patricia, D.; Masterson, T.; Manikhas, G.; Myasnikov, A.; Osmanov, D.; Cordoba, R.; Panizo, C.; et al. A drug-drug interaction study of ibrutinib with moderate/strong CYP3A inhibitors in patients with B-cell malignancies. Leuk. Lymphoma 2018, 59, 2888–2895. [Google Scholar] [CrossRef]
  61. Markowska, A.; Kaysiewicz, J.; Markowska, J.; Huczynski, A. Doxycycline, salinomycin, monensin and ivermectin repositioned as cancer drugs. Bioorg. Med. Chem. Lett. 2019, 29, 1549–1554. [Google Scholar] [CrossRef]
  62. Mendez, D.; Gaulton, A.; Bento, A.P.; Chambers, J.; De Veij, M.; Felix, E.; Magarinos, M.P.; Mosquera, J.F.; Mutowo, P.; Nowotka, M.; et al. ChEMBL: Towards direct deposition of bioassay data. Nucleic Acids Res. 2019, 47, D930–D940. [Google Scholar] [CrossRef]
  63. Gaulton, A.; Hersey, A.; Nowotka, M.; Bento, A.P.; Chambers, J.; Mendez, D.; Mutowo, P.; Atkinson, F.; Bellis, L.J.; Cibrian-Uhalte, E.; et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017, 45, D945–D954. [Google Scholar] [CrossRef]
  64. Chem Axon J. Chem for Office. Available online: https://chemaxon.com (accessed on 12 March 2020).
  65. Chem Axon Chemaxon Standardizer. Available online: http://www.chemaxon.com (accessed on 24 March 2020).
  66. Varnek, A.; Fourches, D.; Horvath, D.; Klimchuk, O.; Gaudin, C.; Vayer, P.; Hoonakker, F.; Tetko, I.V.; Marcou, G. ISIDA-Platform for virtual screening based on fragment and pharmacophoric descriptors. Curr. Comput. Aided Drug Des. 2008, 4, 191. [Google Scholar] [CrossRef]
  67. Ruggiu, F.; Marcou, G.; Varnek, A.; Horvath, D. ISIDA Property-Labelled Fragment Descriptors. Mol. Inform. 2010, 29, 855–868. [Google Scholar] [CrossRef]
  68. Varnek, A.; Fourches, D.; Hoonakker, F.; Solovev, V.P. Substructural fragments: An universal language to encode reactions, molecular and supramolecular structures. J. Comput. Aided Mol. Des. 2005, 19, 693–703. [Google Scholar] [CrossRef]
  69. Peng, H.; Long, F.; Ding, C. Feature selection based on mutual information: Criteria of max-dependency, max-relevance, and min-redundancy. IEEE Trans. Pattern Anal. Mach. Intell. 2005, 27, 1226–1238. [Google Scholar] [CrossRef]
  70. Potter, T.; Matter, H. Random or rational design? Evaluation of diverse compound subsets from chemical structure databases. J. Med. Chem 1998, 41, 478–488. [Google Scholar] [CrossRef]
  71. Serra, M.; Hattinger, C.M. The pharmacogenomics of osteosarcoma. Pharmacogenom. J. 2017, 17, 11–20. [Google Scholar] [CrossRef]
  72. Hattinger, C.M.; Vella, S.; Tavanti, E.; Fanelli, M.; Picci, P.; Serra, M. Pharmacogenomics of second-line drugs used for treatment of unresponsive or relapsed osteosarcoma patients. Pharmacogenomics 2016, 17, 2097–2114. [Google Scholar] [CrossRef] [Green Version]
  73. Schwartz, G.K.; Tap, W.D.; Qin, L.X.; Livingston, M.B.; Undevia, S.D.; Chmielowski, B.; Agulnik, M.; Schuetze, S.M.; Reed, D.R.; Okuno, S.H.; et al. Cixutumumab and temsirolimus for patients with bone and soft-tissue sarcoma: A multicentre, open-label, phase 2 trial. Lancet Oncol. 2013, 14, 371–382. [Google Scholar] [CrossRef] [Green Version]
  74. Trucco, M.M.; Meyer, C.F.; Thornton, K.A.; Shah, P.; Chen, A.R.; Wilky, B.A.; Carrera-Haro, M.A.; Boyer, L.C.; Ferreira, M.F.; Shafique, U.; et al. A phase II study of temsirolimus and liposomal doxorubicin for patients with recurrent and refractory bone and soft tissue sarcomas. Clin. Sarcoma Res. 2018, 8, 21. [Google Scholar] [CrossRef]
  75. Demetri, G.D.; Chawla, S.P.; Ray-Coquard, I.; Le Cesne, A.; Staddon, A.P.; Milhem, M.M.; Penel, N.; Riedel, R.F.; Bui-Nguyen, B.; Cranmer, L.D.; et al. Results of an international randomized phase III trial of the mammalian target of rapamycin inhibitor ridaforolimus versus placebo to control metastatic sarcomas in patients after benefit from prior chemotherapy. J. Clin. Oncol. 2013, 31, 2485–2492. [Google Scholar] [CrossRef]
  76. Chawla, S.P.; Staddon, A.P.; Baker, L.H.; Schuetze, S.M.; Tolcher, A.W.; D’Amato, G.Z.; Blay, J.Y.; Mita, M.M.; Sankhala, K.K.; Berk, L.; et al. Phase II study of the mammalian target of rapamycin inhibitor ridaforolimus in patients with advanced bone and soft tissue sarcomas. J. Clin. Oncol. 2012, 30, 78–84. [Google Scholar] [CrossRef]
  77. Qayed, M.; Cash, T.; Tighiouart, M.; MacDonald, T.J.; Goldsmith, K.C.; Tanos, R.; Kean, L.; Watkins, B.; Suessmuth, Y.; Wetmore, C.; et al. A phase I study of sirolimus in combination with metronomic therapy (CHOAnome) in children with recurrent or refractory solid and brain tumors. Pediatr. Blood Cancer 2020, 67, e28134. [Google Scholar] [CrossRef]
  78. van der Graaf, W.T.; Blay, J.Y.; Chawla, S.P.; Kim, D.W.; Bui-Nguyen, B.; Casali, P.G.; Schoffski, P.; Aglietta, M.; Staddon, A.P.; Beppu, Y.; et al. Pazopanib for metastatic soft-tissue sarcoma (PALETTE): A randomised, double-blind, placebo-controlled phase 3 trial. Lancet 2012, 379, 1879–1886. [Google Scholar] [CrossRef]
  79. Longhi, A.; Paioli, A.; Palmerini, E.; Cesari, M.; Abate, M.E.; Setola, E.; Spinnato, P.; Donati, D.; Hompland, I.; Boye, K. Pazopanib in relapsed osteosarcoma patients: Report on 15 cases. Acta Oncol. 2019, 58, 124–128. [Google Scholar] [CrossRef] [Green Version]
  80. Mysinger, M.M.; Carchia, M.; Irwin, J.J.; Shoichet, B.K. Directory of useful decoys, enhanced (DUD-E): Better ligands and decoys for better benchmarking. J. Med. Chem. 2012, 55, 6582–6594. [Google Scholar] [CrossRef]
  81. Truchon, J.F.; Bayly, C.I. Evaluating virtual screening methods: Good and bad metrics for the “early recognition” problem. J. Chem. Inf. Model. 2007, 47, 488–508. [Google Scholar] [CrossRef]
  82. Kirchmair, J.; Markt, P.; Distinto, S.; Wolber, G.; Langer, T. Evaluation of the performance of 3D virtual screening protocols: RMSD comparisons, enrichment assessments, and decoy selection—What can we learn from earlier mistakes? J. Comput. Aided. Mol. Des. 2008, 22, 213–228. [Google Scholar] [CrossRef]
Figure 1. The chemical diversity of inactive compounds in OS cell lines. (A) Compounds with biological activity reported in ChEMBL for OS cell lines. (B) Dendrograms calculated for inactive compounds in the HOS, MG63, SAOS2 and U2OS cell lines.
Figure 1. The chemical diversity of inactive compounds in OS cell lines. (A) Compounds with biological activity reported in ChEMBL for OS cell lines. (B) Dendrograms calculated for inactive compounds in the HOS, MG63, SAOS2 and U2OS cell lines.
Pharmaceuticals 13 00409 g001
Figure 2. Performance of machine learning models constructed from compounds with biological activity for the HOS, MG63, SAOS2 and U2OS cell lines. Accuracy, sensitivity and specificity values correspond to the external data.
Figure 2. Performance of machine learning models constructed from compounds with biological activity for the HOS, MG63, SAOS2 and U2OS cell lines. Accuracy, sensitivity and specificity values correspond to the external data.
Pharmaceuticals 13 00409 g002
Figure 3. Results of the performance of base models and multi-objective models for the VS. (A) Comparison of AUC values (black bars) and BEDROC with a = 160.9 of base models and the multi-objective algorithm. (B) Accumulative curves for the four top-performing VS protocols. The comparison includes the best 3 simple models and the multi-objective algorithm. Results are presented for the whole screening, and (C) for the top 5% of screened data.
Figure 3. Results of the performance of base models and multi-objective models for the VS. (A) Comparison of AUC values (black bars) and BEDROC with a = 160.9 of base models and the multi-objective algorithm. (B) Accumulative curves for the four top-performing VS protocols. The comparison includes the best 3 simple models and the multi-objective algorithm. Results are presented for the whole screening, and (C) for the top 5% of screened data.
Pharmaceuticals 13 00409 g003
Figure 4. Repositioned drugs for OS treatment. (A) The central pie chart (black) shows the distribution of the 4 main drug classes repositioned at the first 1% of the screened list, whereas the outer pie chart shows the groups they represent. Each color represents a specific group obtained from the Anatomical Therapeutic Chemical (ATC) classification system. The action mechanisms of the screened drugs are also included in italics. (B) Correlation between the top-ranked drugs using the multi-objective model and its action mechanism. Listed are the first 22 positions (1%) of the 2218 DrugBank compounds screened. Drugs and their desirability values obtained using the prediction algorithm are described in the left column, while their action mechanism is in the right column. The colors represent the drug groups described in the graph above.
Figure 4. Repositioned drugs for OS treatment. (A) The central pie chart (black) shows the distribution of the 4 main drug classes repositioned at the first 1% of the screened list, whereas the outer pie chart shows the groups they represent. Each color represents a specific group obtained from the Anatomical Therapeutic Chemical (ATC) classification system. The action mechanisms of the screened drugs are also included in italics. (B) Correlation between the top-ranked drugs using the multi-objective model and its action mechanism. Listed are the first 22 positions (1%) of the 2218 DrugBank compounds screened. Drugs and their desirability values obtained using the prediction algorithm are described in the left column, while their action mechanism is in the right column. The colors represent the drug groups described in the graph above.
Pharmaceuticals 13 00409 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cabrera-Andrade, A.; López-Cortés, A.; Jaramillo-Koupermann, G.; González-Díaz, H.; Pazos, A.; Munteanu, C.R.; Pérez-Castillo, Y.; Tejera, E. A Multi-Objective Approach for Anti-Osteosarcoma Cancer Agents Discovery through Drug Repurposing. Pharmaceuticals 2020, 13, 409. https://0-doi-org.brum.beds.ac.uk/10.3390/ph13110409

AMA Style

Cabrera-Andrade A, López-Cortés A, Jaramillo-Koupermann G, González-Díaz H, Pazos A, Munteanu CR, Pérez-Castillo Y, Tejera E. A Multi-Objective Approach for Anti-Osteosarcoma Cancer Agents Discovery through Drug Repurposing. Pharmaceuticals. 2020; 13(11):409. https://0-doi-org.brum.beds.ac.uk/10.3390/ph13110409

Chicago/Turabian Style

Cabrera-Andrade, Alejandro, Andrés López-Cortés, Gabriela Jaramillo-Koupermann, Humberto González-Díaz, Alejandro Pazos, Cristian R. Munteanu, Yunierkis Pérez-Castillo, and Eduardo Tejera. 2020. "A Multi-Objective Approach for Anti-Osteosarcoma Cancer Agents Discovery through Drug Repurposing" Pharmaceuticals 13, no. 11: 409. https://0-doi-org.brum.beds.ac.uk/10.3390/ph13110409

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