Chemoinformatics in anti-cancer chemotherapy: Multi-target QSAR model for the in silico discovery of anti-breast cancer agents

Artificial Intelligence Models and Tools for the Assessment of Drug-Herb Interactions

Artificial intelligence (AI) has emerged as a powerful tool in medical sciences that is revolutionizing various fields of drug research. AI algorithms can analyze large-scale biological data and identify molecular targets and pathways advancing pharmacological knowledge. An especially promising area is the assessment of drug interactions. The AI analysis of large datasets, such as drugs' chemical structure, pharmacological properties, molecular pathways, and known interaction patterns, can provide mechanistic insights and identify potential associations by integrating all this complex information and returning potential risks associated with these interactions. In this context, an area where AI may prove valuable is in the assessment of the underlying mechanisms of drug interactions with natural products (i.e., herbs) that are used as dietary supplements. These products pose a challenging problem since they are complex mixtures of constituents with diverse and limited information regarding their pharmacological properties, especially their pharmacokinetic data. As the use of herbal products and supplements continues to grow, it becomes increasingly important to understand the potential interactions between them and conventional drugs and the associated adverse drug reactions. This review will discuss AI approaches and how they can be exploited in providing valuable mechanistic insights regarding the prediction of interactions between drugs and herbs, and their potential exploitation in experimental validation or clinical utilization.

Comprehending the pharmacological mechanism of marine phenolic acids in bladder cancer therapy against matrix metalloproteinase 9 protein by integrated network pharmacology and in-silico approaches

Bladder cancer (BC) is the 10th most common tumour with a high incidence and recurrence rate worldwide; however, the current therapies present limitations as, regularly, not all patients benefit from treatment. Therefore, the search for new, active marine phenolic acids with anti-tumour properties is imperative. In this study, we subjected marine phenolic acids to in silico investigations such as network pharmacology, molecular docking, and molecular dynamics simulation (MD) to identify a plausible pathway and the lead compound that inhibits BC. According to the network pharmacology analysis, eight hub genes (PLAU, MMP2, ITGB3, MAPK1, PTPN11, ESR1, TLR4, MMP9) were found and linked to the enrichment of hsa05205: proteoglycans in cancer, and four hub genes (MMP1, MMP2, MAPK1, MMP9) were involved in the enrichment of hsa05219: BC. Subsequently, molecular docking studies showed that the marine phenolic acids exhibit a strong binding affinity for the target protein, matrix metalloproteinase-9 (MPP9). Among these 14 marine phenolic acids, chicoric acid showed the highest binding affinity of -67.1445 kcal/mol and formed hydrogen bonds with the residues of Ala189, Gln227, Leu188, His226, Ala242, Arg249, Ala191, and Gly186 in the active site of the MPP9 protein. Then, molecular dynamics simulation revealed that chicoric acid formed a stable protein-ligand complex with RMSD and RMSF values of 0.72 nm and 0.53 nm, respectively. Furthermore, the PCA method was employed to understand the dynamical behaviour in the conformational space of MPP9 protein bound to chicoric acid, and the results showed the good conformational space behaviour of MPP9 protein. Moreover, chicoric acid showed a free binding energy value of -32.62 kcal/mol, which indicated it could be a BC inhibitor. Overall, chicoric acid demonstrated potential anti-BC activity through MPP9 protein inhibition.

Support Vector Machine-Based Prediction Models for Drug Repurposing and Designing Novel Drugs for Colorectal Cancer

Colorectal cancer (CRC) has witnessed a concerning increase in incidence and poses a significant therapeutic challenge due to its poor prognosis. There is a pressing demand to identify novel drug therapies to combat CRC. In this study, we addressed this need by utilizing the pharmacological profiles of anticancer drugs from the Genomics of Drug Sensitivity in Cancer (GDSC) database and developed QSAR models using the Support Vector Machine (SVM) algorithm for prediction of alternative and promiscuous anticancer compounds for CRC treatment. Our QSAR models demonstrated their robustness by achieving a high correlation of determination (R2) after 10-fold cross-validation. For 12 CRC cell lines, R2 ranged from 0.609 to 0.827. The highest performance was achieved for SW1417 and GP5d cell lines with R2 values of 0.827 and 0.786, respectively. Further, we listed the most common chemical descriptors in the drug profiles of the CRC cell lines and we also further reported the correlation of these descriptors with drug activity. The KRFP314 fingerprint was the predominantly occurring descriptor, with the KRFPC314 fingerprint following closely in prevalence within the drug profiles of the CRC cell lines. Beyond predictive modeling, we also confirmed the applicability of our developed QSAR models via in silico methods by conducting descriptor-drug analyses and recapitulating drug-to-oncogene relationships. We also identified two potential anti-CRC FDA-approved drugs, viomycin and diamorphine, using QSAR models. To ensure the easy accessibility and utility of our research findings, we have incorporated these models into a user-friendly prediction Web server named "ColoRecPred", available at https://project.iith.ac.in/cgntlab/colorecpred. We anticipate that this Web server can be used for screening of chemical libraries to identify potential anti-CRC drugs.

Quantitative multi-species toxicity modeling: Does a multi-species, machine learning model provide better performance than a single-species model for the evaluation of acute aquatic toxicity by organic pollutants?

The toxicological profile of any chemical is defined by multiple endpoints and testing procedures, including representative test species from different trophic levels. While computer-aided methods play an increasingly important role in supporting ecotoxicology research and chemical hazard assessment, most of the recently developed machine learning models are directed towards a single, specific endpoint. To overcome this limitation and accelerate the process of identifying potentially hazardous environmental pollutants, we are introducing an effective approach for quantitative, multi-species modeling. The proposed approach is based on canonical correlation analysis that finds a pair(s) of uncorrelated, linear combinations of the original variables that best defines the overall variability within and between multiple biological responses and predictor variables. Its effectiveness was confirmed by the machine learning model for estimating acute toxicity of diverse organic pollutants in aquatic species from three trophic levels: algae (Pseudokirchneriella subcapitata), daphnia (Daphnia magna), and fish (Oryzias latipes). The multi-species model achieved a favorable predictive performance that were in line with predictive models derived for the aquatic organisms individually. The chemical bioavailability and reactivity parameters (n-octanol/water partition coefficient, chemical potential, and molecular size and volume) were important to accurately predict acute ecotoxicity to the three aquatic organisms. To facilitate the use of this approach, an open-source, Python-based script, named qMTM (quantitative Multi-species Toxicity Modeling) has been provided.

Moving Average-Based Multitasking In Silico Classification Modeling: Where Do We Stand and What Is Next?

Conventional in silico modeling is often viewed as 'one-target' or 'single-task' computer-aided modeling since it mainly relies on forecasting an endpoint of interest from similar input data. Multitasking or multitarget in silico modeling, in contrast, embraces a set of computational techniques that efficiently integrate multiple types of input data for setting up unique in silico models able to predict the outcome(s) relating to various experimental and/or theoretical conditions. The latter, specifically, based upon the Box-Jenkins moving average approach, has been applied in the last decade to several research fields including drug and materials design, environmental sciences, and nanotechnology. The present review discusses the current status of multitasking computer-aided modeling efforts, meanwhile describing both the existing challenges and future opportunities of its underlying techniques. Some important applications are also discussed to exemplify the ability of multitasking modeling in deriving holistic and reliable in silico classification-based models as well as in designing new chemical entities, either through fragment-based design or virtual screening. Focus will also be given to some software recently developed to automate and accelerate such types of modeling. Overall, this review may serve as a guideline for researchers to grasp the scope of multitasking computer-aided modeling as a promising in silico tool.

Using Stochastic Computing for Virtual Screening Acceleration

Stochastic computing is an emerging scientific field pushed by the need for developing high-performance artificial intelligence systems in hardware to quickly solve complex data processing problems. This is the case of virtual screening, a computational task aimed at searching across huge molecular databases for new drug leads. In this work, we show a classification framework in which molecules are described by an energy-based vector. This vector is then processed by an ultra-fast artificial neural network implemented through FPGA by using stochastic computing techniques. Compared to other previously published virtual screening methods, this proposal provides similar or higher accuracy, while it improves processing speed by about two or three orders of magnitude.

Antiproliferative activity of cardenolides on cell line A549: structure-activity relationship analysis

Since the beginning, natural products have represented an important source of bioactive molecules for cancer treatment. Among them, cardenolides attract the attention of different research groups due to their cardiotonic and antitumor activity. The observed biological activity is closely related to their Na+/K+-ATPase inhibition potency. Currently, the discovery of new compounds against cancer is an urgent need in modern pharmaceutical research. Thus, the aim of this work is to determine the physicochemical properties and substituent effects that module the antiproliferative activity of cardenolides on the human lung cancer cell line A549. We build and curate a library with results obtained from literature; molecular descriptors were calculated in PaDEL software, and SAR/QSAR analysis was performed. The SAR results showed that cardenolides were sensitive to modifications in C and D steroidal ring and required substituent groups with the function of hydrogen bond acceptor at the C3 position. QSAR models to doubly linked-type cardenolides indicated that properties as lipoaffinity and atoms with the capacity to be hydrogen bond acceptors are involved in the increment of antiproliferative activity on A549 cell line. In contrast, the presence and position of very electro-negative atoms on the molecule decreased the antiproliferative effect on A549 cells. These results suggest that the antiproliferative capacity of cardenolides on the cell line A549 is strongly related to substituent groups on the C3 position, which must not be carbohydrate. Additionally, the steroidal rings C and D must remain without modifications.

Structure based designing of thiazolidinone-pyrimidine derivatives as ERK2 inhibitors: Synthesis and in vitro evaluation

Breast cancer has been associated with an overexpression of various molecular targets; accordingly, various target-specific chemotherapeutic agents have been developed. Inhibition of ERK2, a member of MAPK pathway, is an important target involved in the treatment of both oestrogen receptor-positive and triple-negative breast cancer. Thus, in continuation of our previous work on the ERK2 target, we here report novel inhibitors of this kinase. Out of three lead molecules reported in our previous study, we selected the thiazolidinone-pyrimidine scaffold for further development of small molecule inhibitors of ERK2. Analogues of the lead molecule were docked in the target kinase, followed by molecular dynamic simulations and MM-GBSA calculations. Analogues maintaining key interactions with amino acid residues in the ATP-binding domain of ERK2 were selected and duly synthesized. In vitro biochemical evaluation of these molecules against ERK2 kinase disclosed that two molecules possess significant kinase inhibitory potential with IC₅₀ values ≤ 0.5 µM.

The QSAR Paradigm in Fragment-Based Drug Discovery: From the Virtual Generation of Target Inhibitors to Multi-Scale Modeling

Fragment-Based Drug Design (FBDD) has established itself as a promising approach in modern drug discovery, accelerating and improving lead optimization, while playing a crucial role in diminishing the high attrition rates at all stages in the drug development process. On the other hand, FBDD has benefited from the application of computational methodologies, where the models derived from the Quantitative Structure-Activity Relationships (QSAR) have become consolidated tools. This mini-review focuses on the evolution and main applications of the QSAR paradigm in the context of FBDD in the last five years. This report places particular emphasis on the QSAR models derived from fragment-based topological approaches to extract physicochemical and/or structural information, allowing to design potentially novel mono- or multi-target inhibitors from relatively large and heterogeneous databases. Here, we also discuss the need to apply multi-scale modeling, to exemplify how different datasets based on target inhibition can be simultaneously integrated and predicted together with other relevant endpoints such as the biological activity against non-biomolecular targets, as well as in vitro and in vivo toxicity and pharmacokinetic properties. In this context, seminal papers are briefly analyzed. As huge amounts of data continue to accumulate in the domains of the chemical, biological and biomedical sciences, it has become clear that drug discovery must be viewed as a multi-scale optimization process. An ideal multi-scale approach should integrate diverse chemical and biological data and also serve as a knowledge generator, enabling the design of potentially optimal chemicals that may become therapeutic agents.

Predicting Metabolic Reaction Networks with Perturbation-Theory Machine Learning (PTML) Models

BACKGROUND

Checking the connectivity (structure) of complex Metabolic Reaction Networks (MRNs) models proposed for new microorganisms with promising properties is an important goal for chemical biology.

OBJECTIVE

In principle, we can perform a hand-on checking (Manual Curation). However, this is a challenging task due to the high number of combinations of pairs of nodes (possible metabolic reactions).

RESULTS

The CPTML linear model obtained using the LDA algorithm is able to discriminate nodes (metabolites) with the correct assignation of reactions from incorrect nodes with values of accuracy, specificity, and sensitivity in the range of 85-100% in both training and external validation data series.

METHODS

In this work, we used Combinatorial Perturbation Theory and Machine Learning techniques to seek a CPTML model for MRNs >40 organisms compiled by Barabasis' group. First, we quantified the local structure of a very large set of nodes in each MRN using a new class of node index called Markov linear indices fk. Next, we calculated CPT operators for 150000 combinations of query and reference nodes of MRNs. Last, we used these CPT operators as inputs of different ML algorithms.

CONCLUSION

Meanwhile, PTML models based on Bayesian network, J48-Decision Tree and Random Forest algorithms were identified as the three best non-linear models with accuracy greater than 97.5%. The present work opens the door to the study of MRNs of multiple organisms using PTML models.

Synthesis, Pharmacological, and Biological Evaluation of 2-Furoyl-Based MIF-1 Peptidomimetics and the Development of a General-Purpose Model for Allosteric Modulators (ALLOPTML)

This work describes the synthesis and pharmacological evaluation of 2-furoyl-based Melanostatin (MIF-1) peptidomimetics as dopamine D₂ modulating agents. Eight novel peptidomimetics were tested for their ability to enhance the maximal effect of tritiated N-propylapomorphine ([³H]-NPA) at D₂ receptors (D₂R). In this series, 2-furoyl-l-leucylglycinamide (6a) produced a statistically significant increase in the maximal [³H]-NPA response at 10 pM (11 ± 1%), comparable to the effect of MIF-1 (18 ± 9%) at the same concentration. This result supports previous evidence that the replacement of proline residue by heteroaromatic scaffolds are tolerated at the allosteric binding site of MIF-1. Biological assays performed for peptidomimetic 6a using cortex neurons from 19-day-old Wistar-Kyoto rat embryos suggest that 6a displays no neurotoxicity up to 100 μM. Overall, the pharmacological and toxicological profile and the structural simplicity of 6a makes this peptidomimetic a potential lead compound for further development and optimization, paving the way for the development of novel modulating agents of D₂R suitable for the treatment of CNS-related diseases. Additionally, the pharmacological and biological data herein reported, along with >20 000 outcomes of preclinical assays, was used to seek a general model to predict the allosteric modulatory potential of molecular candidates for a myriad of target receptors, organisms, cell lines, and biological activity parameters based on perturbation theory (PT) ideas and machine learning (ML) techniques, abbreviated as ALLOPTML. By doing so, ALLOPTML shows high specificity Sp = 89.2/89.4%, sensitivity Sn = 71.3/72.2%, and accuracy Ac = 86.1%/86.4% in training/validation series, respectively. To the best of our knowledge, ALLOPTML is the first general-purpose chemoinformatic tool using a PTML-based model for the multioutput and multicondition prediction of allosteric compounds, which is expected to save both time and resources during the early drug discovery of allosteric modulators.

Cell-based multi-target QSAR model for design of virtual versatile inhibitors of liver cancer cell lines

Liver cancers are one of the leading fatal diseases among malignant neoplasms. Current chemotherapeutic treatments used to fight these illnesses have become less efficient in terms of both efficacy and safety. Therefore, there is a great need of search for new anti-liver cancer agents and this can be accelerated by using computer-aided drug discovery approaches. In this work, we report the development of the first cell-based multi-target model based on quantitative structure-activity relationships (CBMT-QSAR) for the design and prediction of chemicals as anticancer agents against 17 liver cancer cell lines. While having a good quality and predictive power (accuracy higher than 80%) in the training and test sets, respectively, the CBMT-QSAR model was employed as a tool to directly extract suitable fragments from the physicochemical and structural interpretations of the molecular descriptors. Some of these desirable fragments were assembled, leading to the virtual design of eight molecules with drug-like properties, with six of them being predicted as versatile anticancer agents against the 17 liver cancer cell lines reported here.

Perturbation-Theory Machine Learning (PTML) Multilabel Model of the ChEMBL Dataset of Preclinical Assays for Antisarcoma Compounds

Sarcomas are a group of malignant neoplasms of connective tissue with a different etiology than carcinomas. The efforts to discover new drugs with antisarcoma activity have generated large datasets of multiple preclinical assays with different experimental conditions. For instance, the ChEMBL database contains outcomes of 37,919 different antisarcoma assays with 34,955 different chemical compounds. Furthermore, the experimental conditions reported in this dataset include 157 types of biological activity parameters, 36 drug targets, 43 cell lines, and 17 assay organisms. Considering this information, we propose combining perturbation theory (PT) principles with machine learning (ML) to develop a PTML model to predict antisarcoma compounds. PTML models use one function of reference that measures the probability of a drug being active under certain conditions (protein, cell line, organism, etc.). In this paper, we used a linear discriminant analysis and neural network to train and compare PT and non-PT models. All the explored models have an accuracy of 89.19-95.25% for training and 89.22-95.46% in validation sets. PTML-based strategies have similar accuracy but generate simplest models. Therefore, they may become a versatile tool for predicting antisarcoma compounds.

Thidiazuron suppresses breast cancer via targeting miR-132 and dysregulation of the PI3K-Akt signaling pathway mediated by the miR-202-5p-PTEN axis

Chemo-resistance and metastasis are the most common causes of breast cancer recurrence and death. Thidiazuron (TDZ) is a plant growth regulator (phytohormone) whose biological effects on humans and animals has not yet been determined. In this study, we investigated the anticancer activity of this phytohormone on the drug resistant-triple negative breast cancer cell line MDA-MB-231. Treatment of the breast cancer cells with TDZ (1-50 μmol/L) caused more stressful environment and induced a significant increase in active caspase-positive cells. In addition, TDZ treatment (5 and 10 μmol/L) significantly attenuated the migration and the invasiveness of these highly metastatic cancer cells. Mechanistically, TDZ reduces cancer progression and invasiveness by targeting miR-202-5p, which stimulates the expression of phosphatase and tensin homolog (PTEN), the tumor suppressor that downregulates the PI3K-Akt signaling pathway. Treatment with TDZ significantly upregulates miRNA-132, the suppressor of breast cancer proliferation, which is also implicated in dysregulation of the TEN-Akt-NFκB signaling pathway. Interestingly, our molecular docking analysis revealed a potential non-covalent interaction between TDZ and Akt, PTEN, and PI3K. These findings suggest that TDZ suppresses breast cancer metastasis by targeting miRNA-132, the miR-202-5p-PTEN axis, and the PI3K-Akt signaling pathway downstream.

Advances in Nanoparticles as Anticancer Drug Delivery Vector: Need of this Century

Background:

Nanotechnology has contributed a great deal to the field of medical science. Smart drugdelivery vectors, combined with stimuli-based characteristics, are becoming increasingly important. The use of external and internal stimulating factors can have enormous benefits and increase the targeting efficiency of nanotechnology platforms. The pH values of tumor vascular tissues are acidic in nature, allowing the improved targeting of anticancer drug payloads using drug-delivery vectors. Nanopolymers are smart drug-delivery vectors that have recently been developed and recommended for use by scientists because of their potential targeting capabilities, non-toxicity and biocompatibility, and make them ideal nanocarriers for personalized drug delivery.

Method:

The present review article provides an overview of current advances in the use of nanoparticles (NPs) as anticancer drug-delivery vectors.

Results:

This article reviews the molecular basis for the use of NPs in medicine, including personalized medicine, personalized therapy, emerging vistas in anticancer therapy, nanopolymer targeting, passive and active targeting transports, pH-responsive drug carriers, biological barriers, computer-aided drug design, future challenges and perspectives, biodegradability and safety.

Conclusions:

This article will benefit academia, researchers, clinicians, and government authorities by providing a basis for further research advancements.

A Review on Applications of Computational Methods in Drug Screening and Design

Drug development is one of the most significant processes in the pharmaceutical industry. Various computational methods have dramatically reduced the time and cost of drug discovery. In this review, we firstly discussed roles of multiscale biomolecular simulations in identifying drug binding sites on the target macromolecule and elucidating drug action mechanisms. Then, virtual screening methods (e.g., molecular docking, pharmacophore modeling, and QSAR) as well as structure- and ligand-based classical/de novo drug design were introduced and discussed. Last, we explored the development of machine learning methods and their applications in aforementioned computational methods to speed up the drug discovery process. Also, several application examples of combining various methods was discussed. A combination of different methods to jointly solve the tough problem at different scales and dimensions will be an inevitable trend in drug screening and design.

MCDCalc: Markov Chain Molecular Descriptors Calculator for Medicinal Chemistry

AIMS

Cheminformatics models are able to predict different outputs (activity, property, chemical reactivity) in single molecules or complex molecular systems (catalyzed organic synthesis, metabolic reactions, nanoparticles, etc.).

BACKGROUND

Cheminformatics models are able to predict different outputs (activity, property, chemical reactivity) in single molecules or complex molecular systems (catalyzed organic synthesis, metabolic reactions, nanoparticles, etc.).

OBJECTIVE

Cheminformatics prediction of complex catalytic enantioselective reactions is a major goal in organic synthesis research and chemical industry. Markov Chain Molecular Descriptors (MCDs) have been largely used to solve Cheminformatics problems. There are different types of Markov chain descriptors such as Markov-Shannon entropies (Shk), Markov Means (Mk), Markov Moments (πk), etc. However, there are other possible MCDs that have not been used before. In addition, the calculation of MCDs is done very often using specific software not always available for general users and there is not an R library public available for the calculation of MCDs. This fact, limits the availability of MCMDbased Cheminformatics procedures.

METHODS

We studied the enantiomeric excess ee(%)[Rcat] for 324 α-amidoalkylation reactions. These reactions have a complex mechanism depending on various factors. The model includes MCDs of the substrate, solvent, chiral catalyst, product along with values of time of reaction, temperature, load of catalyst, etc. We tested several Machine Learning regression algorithms. The Random Forest regression model has R2 > 0.90 in training and test. Secondly, the biological activity of 5644 compounds against colorectal cancer was studied.

RESULTS

We developed very interesting model able to predict with Specificity and Sensitivity 70-82% the cases of preclinical assays in both training and validation series.

CONCLUSION

The work shows the potential of the new tool for computational studies in organic and medicinal chemistry.

Rational Drug Design of Antineoplastic Agents Using 3D-QSAR, Cheminformatic, and Virtual Screening Approaches

Background:

Computer-Aided Drug Design has strongly accelerated the development of novel antineoplastic agents by helping in the hit identification, optimization, and evaluation.

Results:

Computational approaches such as cheminformatic search, virtual screening, pharmacophore modeling, molecular docking and dynamics have been developed and applied to explain the activity of bioactive molecules, design novel agents, increase the success rate of drug research, and decrease the total costs of drug discovery. Similarity, searches and virtual screening are used to identify molecules with an increased probability to interact with drug targets of interest, while the other computational approaches are applied for the design and evaluation of molecules with enhanced activity and improved safety profile.

Conclusion:

In this review are described the main in silico techniques used in rational drug design of antineoplastic agents and presented optimal combinations of computational methods for design of more efficient antineoplastic drugs.

Predicting Synergism of Cancer Drug Combinations Using NCI-ALMANAC Data

Drug combinations are of great interest for cancer treatment. Unfortunately, the discovery of synergistic combinations by purely experimental means is only feasible on small sets of drugs. In silico modeling methods can substantially widen this search by providing tools able to predict which of all possible combinations in a large compound library are synergistic. Here we investigate to which extent drug combination synergy can be predicted by exploiting the largest available dataset to date (NCI-ALMANAC, with over 290,000 synergy determinations). Each cell line is modeled using primarily two machine learning techniques, Random Forest (RF) and Extreme Gradient Boosting (XGBoost), on the datasets provided by NCI-ALMANAC. This large-scale predictive modeling study comprises more than 5,000 pair-wise drug combinations, 60 cell lines, 4 types of models, and 5 types of chemical features. The application of a powerful, yet uncommonly used, RF-specific technique for reliability prediction is also investigated. The evaluation of these models shows that it is possible to predict the synergy of unseen drug combinations with high accuracy (Pearson correlations between 0.43 and 0.86 depending on the considered cell line, with XGBoost providing slightly better predictions than RF). We have also found that restricting to the most reliable synergy predictions results in at least 2-fold error decrease with respect to employing the best learning algorithm without any reliability estimation. Alkylating agents, tyrosine kinase inhibitors and topoisomerase inhibitors are the drugs whose synergy with other partner drugs are better predicted by the models. Despite its leading size, NCI-ALMANAC comprises an extremely small part of all conceivable combinations. Given their accuracy and reliability estimation, the developed models should drastically reduce the number of required in vitro tests by predicting in silico which of the considered combinations are likely to be synergistic.

Toward steroidal anticancer drugs: Non-parametric and 3D-QSAR modeling of 17-picolyl and 17-picolinylidene androstanes with antiproliferative activity on breast adenocarcinoma cells

The present study is aimed to analyze lipophilicity and ADMET profiles, and to develop field based 3D-QSAR and ligand-based pharmacophore hypothesis for a series of 17α-picolyl and 17(E)-picolinylidene androstane derivatives in order to give detailed structural insights and to highlight important binding features of novel androstane derivatives, as compounds with antiproliferative activity toward breast adenocarcinoma cells. This study can provide guidelines for the rational design of novel potent compounds. Sum of ranking differences (SRD), as a non-parametric method, was applied for compounds ranking. 3D-QSAR methods, including comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), were applied to predict the antiproliferative effect on breast adenocarcinoma cells and provide the regions in space where interactive fields may influence the activity. The compounds are ranked so the compounds with the most favorable ADME and lipophilicity features together with significant anticancer activity can be distinguished. The established 3D-QSAR model could be used for design of new compounds with antiproliferative activity on the human ER- breast adenocarcinoma cells. The pharmacophore model is able to accurately predict antiproliferative activity. Generally, the present study provides significant guidelines for further selection, synthesis and rational design of new highly potential androstane derivatives as anticancer drugs.

PTML Combinatorial Model of ChEMBL Compounds Assays for Multiple Types of Cancer

Determining the target proteins of new anticancer compounds is a very important task in Medicinal Chemistry. In this sense, chemists carry out preclinical assays with a high number of combinations of experimental conditions (c _j). In fact, ChEMBL database contains outcomes of 65 534 different anticancer activity preclinical assays for 35 565 different chemical compounds (1.84 assays per compound). These assays cover different combinations of c _j formed from >70 different biological activity parameters ( c₀), >300 different drug targets ( c₁), >230 cell lines ( c₂), and 5 organisms of assay ( c₃) or organisms of the target ( c₄). It include a total of 45 833 assays in leukemia, 6227 assays in breast cancer, 2499 assays in ovarian cancer, 3499 in colon cancer, 3159 in lung cancer, 2750 in prostate cancer, 601 in melanoma, etc. This is a very complex data set with multiple Big Data features. This data is hard to be rationalized by researchers to extract useful relationships and predict new compounds. In this context, we propose to combine perturbation theory (PT) ideas and machine learning (ML) modeling to solve this combinatorial-like problem. In this work, we report a PTML (PT + ML) model for ChEMBL data set of preclinical assays of anticancer compounds. This is a simple linear model with only three variables. The model presented values of area under receiver operating curve = AUROC = 0.872, specificity = Sp(%) = 90.2, sensitivity = Sn(%) = 70.6, and overall accuracy = Ac(%) = 87.7 in training series. The model also have Sp(%) = 90.1, Sn(%) = 71.4, and Ac(%) = 87.8 in external validation series. The model use PT operators based on multicondition moving averages to capture all the complexity of the data set. We also compared the model with nonlinear artificial neural network (ANN) models obtaining similar results. This confirms the hypothesis of a linear relationship between the PT operators and the classification as anticancer compounds in different combinations of assay conditions. Last, we compared the model with other PTML models reported in the literature concluding that this is the only one PTML model able to predict activity against multiple types of cancer. This model is a simple but versatile tool for the prediction of the targets of anticancer compounds taking into consideration multiple combinations of experimental conditions in preclinical assays.

Antiproliferative and proapoptotic activities of aza-annulated naphthoquinone analogs

1,4-Naphthoquinone derivatives have been widely documented with regard to their biological properties, and particularly their anticancer activities. In the 9,10-anthraquinone family, aza-annulation involving one of the carbonyl oxygen atoms has afforded more potent, possibly less toxic analogues. We recently carried out different modifications on the naphthoquinone skeleton to generate 3-chloro-2-amino- and 3-chloro-2-(N-acetamido)-1,4-naphthoquinone and 3,4-dihydrobenzo[f]quinoxalin-6(2H)-one derivatives. These three series of compounds were now tested against normal human fibroblasts and six human cancer cell lines. Some of the dihydrobenzoquinoxalinone derivatives were not only more potent than their 1,4-naphthoquinone counterparts, but also exhibited 10- to 14-fold selectivity between bladder carcinoma and normal cells and were equipotent with the non-selective reference drug used (etoposide). The fusion of an additional azaheterocycle to the 1,4-naphthoquinone nucleus modulates both the activity, selectivity and mechanism of action of the compounds. The electrochemical properties of selected compounds were evaluated in an attempt to correlate them with cytotoxic activity and mechanism of action. Finally, 3D-QSAR CoMFA and CoMSIA models were built on the AGS, J82, and HL-60 cell lines. The best models had values of r²_pred = 0.815; 0.823 and 0.925. The main structural relationships found, suggest that acetylation and alkylation of the amino group with large groups would be beneficial for cytotoxic activity.

Machine learning and data science in soft materials engineering

In many branches of materials science it is now routine to generate data sets of such large size and dimensionality that conventional methods of analysis fail. Paradigms and tools from data science and machine learning can provide scalable approaches to identify and extract trends and patterns within voluminous data sets, perform guided traversals of high-dimensional phase spaces, and furnish data-driven strategies for inverse materials design. This topical review provides an accessible introduction to machine learning tools in the context of soft and biological materials by 'de-jargonizing' data science terminology, presenting a taxonomy of machine learning techniques, and surveying the mathematical underpinnings and software implementations of popular tools, including principal component analysis, independent component analysis, diffusion maps, support vector machines, and relative entropy. We present illustrative examples of machine learning applications in soft matter, including inverse design of self-assembling materials, nonlinear learning of protein folding landscapes, high-throughput antimicrobial peptide design, and data-driven materials design engines. We close with an outlook on the challenges and opportunities for the field.

The advancement of multidimensional QSAR for novel drug discovery - where are we headed?

INTRODUCTION

The Multidimensional quantitative structure-activity relationship (multidimensional-QSAR) method is one of the most popular computational methods employed to predict interesting biochemical properties of existing or hypothetical molecules. With continuous progress, the QSAR method has made remarkable success in various fields, such as medicinal chemistry, material science and predictive toxicology. Areas covered: In this review, the authors cover the basic elements of multidimensional -QSAR including model construction, validation and application. It includes and emphasizes the very recent developments of multidimensional -QSAR such as: HQSAR, G-QSAR, MIA-QSAR, multi-target QSAR. The advantages and disadvantages of each method are also discussed and typical examples of their application are detailed. Expert opinion: Although there are defects in multidimensional-QSAR modeling, it is still of enormous help to chemists, biologists and other researchers in various fields. In the authors' opinion, the latest more precise and feasible QSAR models should be further developed by integrating new descriptors, algorithms and other relevant computational techniques. Apart from being applied in traditional fields (e.g. lead optimization and predictive risk assessment), QSAR should be used more widely as a routine method in other emerging research fields including the modeling of nanoparticles(NPs), mixture toxicity and peptides.

Designing multi-targeted agents: An emerging anticancer drug discovery paradigm

The dominant paradigm in drug discovery is to design ligands with maximum selectivity to act on individual drug targets. With the target-based approach, many new chemical entities have been discovered, developed, and further approved as drugs. However, there are a large number of complex diseases such as cancer that cannot be effectively treated or cured only with one medicine to modulate the biological function of a single target. As simultaneous intervention of two (or multiple) cancer progression relevant targets has shown improved therapeutic efficacy, the innovation of multi-targeted drugs has become a promising and prevailing research topic and numerous multi-targeted anticancer agents are currently at various developmental stages. However, most multi-pharmacophore scaffolds are usually discovered by serendipity or screening, while rational design by combining existing pharmacophore scaffolds remains an enormous challenge. In this review, four types of multi-pharmacophore modes are discussed, and the examples from literature will be used to introduce attractive lead compounds with the capability of simultaneously interfering with different enzyme or signaling pathway of cancer progression, which will reveal the trends and insights to help the design of the next generation multi-targeted anticancer agents.

Biological Evaluation in Vitro and in Silico of Azetidin-2-one Derivatives as Potential Anticancer Agents

Potential anticancer activity of 16 azetidin-2-one derivatives was evaluated showing that compound 6 [N-(p-methoxy-phenyl)-2-(p-methyl-phenyl)-3-phenoxy-azetidin-2-one] presented cytotoxic activity in SiHa cells and B16F10 cells. The caspase-3 assay in B16F10 cells displayed that azetidin-2-one derivatives induce apoptosis. Microarray and molecular analysis showed that compound 6 was involved on specific gene overexpression of cytoskeleton regulation and apoptosis due to the inhibition of some cell cycle genes. From the 16 derivatives, compound 6 showed the highest selectivity to neoplastic cells, it was an inducer of apoptosis, and according to an in silico analysis of chemical interactions with colchicine binding site of human α/β-tubulin, the mechanism of action could be a molecular interaction involving the amino acids outlining such binding site.

Enabling the Discovery and Virtual Screening of Potent and Safe Antimicrobial Peptides. Simultaneous Prediction of Antibacterial Activity and Cytotoxicity

Antimicrobial peptides (AMPs) represent promising alternatives to fight against bacterial pathogens. However, cellular toxicity remains one of the main concerns in the early development of peptide-based drugs. This work introduces the first multitasking (mtk) computational model focused on performing simultaneous predictions of antibacterial activities, and cytotoxicities of peptides. The model was created from a data set containing 3592 cases, and it displayed accuracy higher than 96% for classifying/predicting peptides in both training and prediction (test) sets. The technique known as alanine scanning was computationally applied to illustrate the calculation of the quantitative contributions of the amino acids (in their respective positions of the sequence) to the biological effects of a defined peptide. A small library formed by 10 peptides was generated, where peptides were designed by considering the interpretations of the different descriptors in the mtk-computational model. All the peptides were predicted to exhibit high antibacterial activities against multiple bacterial strains, and low cytotoxicity against various cell types. The present mtk-computational model can be considered a very useful tool to support high throughput research for the discovery of potent and safe AMPs.

A web server for analysis, comparison and prediction of protein ligand binding sites

BACKGROUND

One of the major challenges in the field of system biology is to understand the interaction between a wide range of proteins and ligands. In the past, methods have been developed for predicting binding sites in a protein for a limited number of ligands.

RESULTS

In order to address this problem, we developed a web server named 'LPIcom' to facilitate users in understanding protein-ligand interaction. Analysis, comparison and prediction modules are available in the "LPIcom' server to predict protein-ligand interacting residues for 824 ligands. Each ligand must have at least 30 protein binding sites in PDB. Analysis module of the server can identify residues preferred in interaction and binding motif for a given ligand; for example residues glycine, lysine and arginine are preferred in ATP binding sites. Comparison module of the server allows comparing protein-binding sites of multiple ligands to understand the similarity between ligands based on their binding site. This module indicates that ATP, ADP and GTP ligands are in the same cluster and thus their binding sites or interacting residues exhibit a high level of similarity. Propensity-based prediction module has been developed for predicting ligand-interacting residues in a protein for more than 800 ligands. In addition, a number of web-based tools have been integrated to facilitate users in creating web logo and two-sample between ligand interacting and non-interacting residues.

CONCLUSIONS

In summary, this manuscript presents a web-server for analysis of ligand interacting residue. This server is available for public use from URL http://crdd.osdd.net/raghava/lpicom .

Prediction of anticancer molecules using hybrid model developed on molecules screened against NCI-60 cancer cell lines

Background

In past, numerous quantitative structure-activity relationship (QSAR) based models have been developed for predicting anticancer activity for a specific class of molecules against different cancer drug targets. In contrast, limited attempt have been made to predict the anticancer activity of a diverse class of chemicals against a wide variety of cancer cell lines. In this study, we described a hybrid method developed on thousands of anticancer and non-anticancer molecules tested against National Cancer Institute (NCI) 60 cancer cell lines.

Results

Our analysis of anticancer molecules revealed that majority of anticancer molecules contains 18–24 carbon atoms and are dominated by functional groups like R₂NH, R₃N, ROH, RCOR, and ROR. It was also observed that certain substructures (e.g., 1-methoxy-4-methylbenzene, 1-methoxy benzene, Nitrobenzene, Indole, Propenyl benzene) are more abundant in anticancer molecules. Next, we developed anticancer molecule prediction models using various machine-learning techniques and achieved maximum matthews correlation coefficient (MCC) of 0.81 with 90.40 % accuracy using support vector machine (SVM) based models. In another approach, a novel similarity or potency score based method has been developed using selected fragments/fingerprints and achieved maximum MCC of 0.82 with 90.65 % accuracy. Finally, we combined the strength of above methods and developed a hybrid method with maximum MCC of 0.85 with 92.47 % accuracy.

Conclusions

We developed a hybrid method utilizing the best of machine learning and potency score based method. The highly accurate hybrid method can be used for classification of anticancer and non-anticancer molecules. In order to facilitate scientific community working in the field of anticancer drug discovery, we integrate hybrid and potency method in a web server CancerIN. This server provides various facilities that includes; virtual screening of anticancer molecules, analog based drug design, and similarity with known anticancer molecules (http://crdd.osdd.net/oscadd/cancerin).

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2082-y) contains supplementary material, which is available to authorized users.

Virtual screening of chemical compounds active against breast cancer cell lines based on cell cycle modelling, prediction of cytotoxicity and interaction with targets

Bio- and chemoinformatics methods are widely used for the detection of mechanisms of cancer, to search for potential drug targets and their ligands. Regulatory network analysis based on signalling pathways, and cell cycle regulation provides better understanding of diseases with multiple mechanisms of pathogenesis. We developed an approach for in silico prediction of the cytotoxic effect of chemical compounds in non-transformed and breast cancer cell lines. This approach combines the prediction of the interaction between chemical compounds and human proteins, cytotoxicity and regulatory network modelling taking into account gene expression. Application of our approach to virtual screening of libraries of commercially available compounds allowed selection of dozens of promising hits. These molecules are predicted to interact with the identified targets and exhibit cytotoxicity against breast cancer cell lines but not non-tumour human cell lines. Experimental testing of 49 selected compounds against MDA-MB-231 and MCF7 breast cancer cell lines confirmed the activity of eight compounds with IC50 values ranged from 0.8 to 50 μM. Thus, the developed approach may be applied for virtual screening for cytotoxic compounds against tumour cell lines.

Hologram quantitative structure–activity relationship and comparative molecular interaction field analysis of aminothiazole and thiazolesulfonamide as reversible LSD1 inhibitors

Background: LSD-1 is an enzyme that removes methyl groups from lysine residues of histone proteins. LSD-1 inhibition decreases cellular proliferation and therefore represents a therapeutic target for cancer treatment. MAO and LSD-1 are both flavin adenine dinucleotide-dependent MAOs, and the MAO inhibitor, tranylcypromine, is currently undergoing clinical trials for cancer treatment because it acts as an irreversible LSD-1 inhibitor. Materials & methods: The present study investigated new reversible LSD-1 inhibitors, in order to develop novel selective anticancer agents. We constructed 2 and 3D quantitative structure–activity relationship models by using a series of 54 aminothiazole and thiazolesulfonamide derivatives. Results: The models were validated internally and externally (q2 , 0.691 and 0.701; r2 , 0.894 and 0.937; r2 test , 0.785 and 0.644, for 2 and 3D models, respectively). Fragment contribution maps, as well as steric and electrostatic contour maps were generated in order to obtain chemical information related to LSD-1 inhibition. Conclusion: The thiazolesulfonamide group was fundamental to the inhibition of LSD-1 by these compounds and that bulky and aromatic substituents at the thiazole ring were important for their steric and electrostatic interactions with the active site of LSD-1.

Experimental and computational studies of fatty acid distribution networks

A new PT-LFER model is useful for predicting a distribution network in terms of specific fatty acid distribution.

A general ANN-based multitasking model for the discovery of potent and safer antibacterial agents

Bacteria have been one of the world's most dangerous and deadliest pathogens for mankind, nowadays giving rise to significant public health concerns. Given the prevalence of these microbial pathogens and their increasing resistance to existing antibiotics, there is a pressing need for new antibacterial drugs. However, development of a successful drug is a complex, costly, and time-consuming process. Quantitative Structure-Activity Relationships (QSAR)-based approaches are valuable tools for shortening the time of lead compound identification but also for focusing and limiting time-costly synthetic activities and in vitro/vivo evaluations. QSAR-based approaches, supported by powerful statistical techniques such as artificial neural networks (ANNs), have evolved to the point of integrating dissimilar types of chemical and biological data. This chapter reports an overview of the current research and potential applications of QSAR modeling tools toward the rational design of more efficient antibacterial agents. Particular emphasis is given to the setup of multitasking models along with ANNs aimed at jointly predicting different antibacterial activities and safety profiles of drugs/chemicals under diverse experimental conditions.

Prediction of multi-target networks of neuroprotective compounds with entropy indices and synthesis, assay, and theoretical study of new asymmetric 1,2-rasagiline carbamates

In a multi-target complex network, the links (L_ij) represent the interactions between the drug (d_i) and the target (t_j), characterized by different experimental measures (K_i, K_m, IC₅₀, etc.) obtained in pharmacological assays under diverse boundary conditions (c_j). In this work, we handle Shannon entropy measures for developing a model encompassing a multi-target network of neuroprotective/neurotoxic compounds reported in the CHEMBL database. The model predicts correctly >8300 experimental outcomes with Accuracy, Specificity, and Sensitivity above 80%–90% on training and external validation series. Indeed, the model can calculate different outcomes for >30 experimental measures in >400 different experimental protocolsin relation with >150 molecular and cellular targets on 11 different organisms (including human). Hereafter, we reported by the first time the synthesis, characterization, and experimental assays of a new series of chiral 1,2-rasagiline carbamate derivatives not reported in previous works. The experimental tests included: (1) assay in absence of neurotoxic agents; (2) in the presence of glutamate; and (3) in the presence of H₂O₂. Lastly, we used the new Assessing Links with Moving Averages (ALMA)-entropy model to predict possible outcomes for the new compounds in a high number of pharmacological tests not carried out experimentally.

Simultaneous virtual prediction of anti-Escherichia coli activities and ADMET profiles: A chemoinformatic complementary approach for high-throughput screening

Escherichia coli remains one of the principal pathogens that cause nosocomial infections, medical conditions that are increasingly common in healthcare facilities. E. coli is intrinsically resistant to many antibiotics, and multidrug-resistant strains have emerged recently. Chemoinformatics has been a great ally of experimental methodologies such as high-throughput screening, playing an important role in the discovery of effective antibacterial agents. However, there is no approach that can design safer anti-E. coli agents, because of the multifactorial nature and complexity of bacterial diseases and the lack of desirable ADMET (absorption, distribution, metabolism, elimination, and toxicity) profiles as a major cause of disapproval of drugs. In this work, we introduce the first multitasking model based on quantitative-structure biological effect relationships (mtk-QSBER) for simultaneous virtual prediction of anti-E. coli activities and ADMET properties of drugs and/or chemicals under many experimental conditions. The mtk-QSBER model was developed from a large and heterogeneous data set of more than 37800 cases, exhibiting overall accuracies of >95% in both training and prediction (validation) sets. The utility of our mtk-QSBER model was demonstrated by performing virtual prediction of properties for the investigational drug avarofloxacin (AVX) under 260 different experimental conditions. Results converged with the experimental evidence, confirming the remarkable anti-E. coli activities and safety of AVX. Predictions also showed that our mtk-QSBER model can be a promising computational tool for virtual screening of desirable anti-E. coli agents, and this chemoinformatic approach could be extended to the search for safer drugs with defined pharmacological activities.

Model for vaccine design by prediction of B-epitopes of IEDB given perturbations in peptide sequence, in vivo process, experimental techniques, and source or host organisms

Perturbation methods add variation terms to a known experimental solution of one problem to approach a solution for a related problem without known exact solution. One problem of this type in immunology is the prediction of the possible action of epitope of one peptide after a perturbation or variation in the structure of a known peptide and/or other boundary conditions (host organism, biological process, and experimental assay). However, to the best of our knowledge, there are no reports of general-purpose perturbation models to solve this problem. In a recent work, we introduced a new quantitative structure-property relationship theory for the study of perturbations in complex biomolecular systems. In this work, we developed the first model able to classify more than 200,000 cases of perturbations with accuracy, sensitivity, and specificity >90% both in training and validation series. The perturbations include structural changes in >50000 peptides determined in experimental assays with boundary conditions involving >500 source organisms, >50 host organisms, >10 biological process, and >30 experimental techniques. The model may be useful for the prediction of new epitopes or the optimization of known peptides towards computational vaccine design.

S. Cordeiro MN, Yañez M, González-Díaz H, García-Mera X. Model for high-throughput screening of multitarget drugs in chemical neurosciences: synthesis, assay, and theoretic study of rasagiline carbamates

The disappointing results obtained in recent clinical trials renew the interest in experimental/computational techniques for the discovery of neuroprotective drugs. In this context, multitarget or multiplexing QSAR models (mt-QSAR/mx-QSAR) may help to predict neurotoxicity/neuroprotective effects of drugs in multiple assays, on drug targets, and in model organisms. In this work, we study a data set downloaded from CHEMBL; each data point (>8000) contains the values of one out of 37 possible measures of activity, 493 assays, 169 molecular or cellular targets, and 11 different organisms (including human) for a given compound. In this work, we introduce the first mx-QSAR model for neurotoxicity/neuroprotective effects of drugs based on the MARCH-INSIDE (MI) method. First, we used MI to calculate the stochastic spectral moments (structural descriptors) of all compounds. Next, we found a model that classified correctly 2955 out of 3548 total cases in the training and validation series with Accuracy, Sensitivity, and Specificity values>80%. The model also showed excellent results in Computational-Chemistry simulations of High-Throughput Screening (CCHTS) experiments, with accuracy=90.6% for 4671 positive cases. Next, we reported the synthesis, characterization, and experimental assays of new rasagiline derivatives. We carried out three different experimental tests: assay (1) in the absence of neurotoxic agents, assay (2) in the presence of glutamate, and assay (3) in the presence of H2O2. Compounds 11 with 27.4%, 8 with 11.6%, and 9 with 15.4% showed the highest neuroprotective effects in assays (1), (2), and (3), respectively. After that, we used the mx-QSAR model to carry out a CCHTS of the new compounds in >400 unique pharmacological tests not carried out experimentally. Consequently, this model may become a promising auxiliary tool for the discovery of new drugs for the treatment of neurodegenerative diseases.

Fragment-based optimization of small molecule CXCL12 inhibitors for antagonizing the CXCL12/CXCR4 interaction

The chemokine CXCL12 and its G protein-coupled receptor (GPCR) CXCR4 are high-priority clinical targets because of their involvement in metastatic cancers (also implicated in autoimmune disease and cardiovascular disease). Because chemokines interact with two distinct sites to bind and activate their receptors, both the GPCRs and chemokines are potential targets for small molecule inhibition. A number of chemokines have been validated as targets for drug development, but virtually all drug discovery efforts focus on the GPCRs. However, all CXCR4 receptor antagonists with the exception of MSX-122 have failed in clinical trials due to unmanageable toxicities, emphasizing the need for alternative strategies to interfere with CXCL12/CXCR4-guided metastatic homing. Although targeting the relatively featureless surface of CXCL12 was presumed to be challenging, focusing efforts at the sulfotyrosine (sY) binding pockets proved successful for procuring initial hits. Using a hybrid structure-based in silico/NMR screening strategy, we recently identified a ligand that occludes the receptor recognition site. From this initial hit, we designed a small fragment library containing only nine tetrazole derivatives using a fragment-based and bioisostere approach to target the sY binding sites of CXCL12. Compound binding modes and affinities were studied by 2D NMR spectroscopy, X-ray crystallography, molecular docking and cell-based functional assays. Our results demonstrate that the sY binding sites are conducive to the development of high affinity inhibitors with better ligand efficiency (LE) than typical protein-protein interaction inhibitors (LE ≤ 0.24). Our novel tetrazole-based fragment 18 was identified to bind the sY21 site with a K(d) of 24 μM (LE = 0.30). Optimization of 18 yielded compound 25 which specifically inhibits CXCL12-induced migration with an improvement in potency over the initial hit 9. The fragment from this library that exhibited the highest affinity and ligand efficiency (11: K(d) = 13 μM, LE = 0.33) may serve as a starting point for development of inhibitors targeting the sY12 site.