Prediction of factor Xa inhibitors by machine learning methods

Machine Learning Enabled Structure-Based Drug Repurposing Approach to Identify Potential CYP1B1 Inhibitors

Drug-metabolizing enzyme (DME)-mediated pharmacokinetic resistance of some clinically approved anticancer agents is one of the main reasons for cancer treatment failure. In particular, some commonly used anticancer medicines, including docetaxel, tamoxifen, imatinib, cisplatin, and paclitaxel, are inactivated by CYP1B1. Currently, no approved drugs are available to treat this CYP1B1-mediated inactivation, making the pharmaceutical industries strive to discover new anticancer agents. Because of the extreme complexity and high risk in drug discovery and development, it is worthwhile to come up with a drug repurposing strategy that may solve the resistance problem of existing chemotherapeutics. Therefore, in the current study, a drug repurposing strategy was implemented to find the possible CYP1B1 inhibitors using machine learning (ML) and structure-based virtual screening (SB-VS) approaches. Initially, three different ML models were developed such as support vector machines (SVMs), random forest (RF), and artificial neural network (ANN); subsequently, the best-selected ML model was employed for virtual screening of the selleckchem database to identify potential CYP1B1 inhibitors. The inhibition potency of the obtained hits was judged by analyzing the crucial active site amino acid interactions against CYP1B1. After a thorough assessment of docking scores, binding affinities, as well as binding modes, four compounds were selected and further subjected to in vitro analysis. From the in vitro analysis, it was observed that chlorprothixene, nadifloxacin, and ticagrelor showed promising inhibitory activity toward CYP1B1 in the IC₅₀ range of 0.07-3.00 μM. These new chemical scaffolds can be explored as adjuvant therapies to address CYP1B1-mediated drug-resistance problems.

Identification of potential matrix metalloproteinase-2 inhibitors from natural products through advanced machine learning-based cheminformatics approaches

Matrix metalloproteinase-2 (MMP-2) is capable of degrading Collage TypeIV in the vascular basement membrane and extracellular matrix. Studies have shown that MMP-2 is tightly associated with the biological behavior of malignant tumors. Therefore, the identification of inhibitors targeting MMP-2 could be effective in treating the disease by maintaining extracellular matrix homeostasis. In the pharmaceutical and biomedical fields, many computational tools are widely used, which improve the efficiency of the whole process to some extent. Apart from the conventional cheminformatics approaches (e.g., pharmacophore model and molecular docking), virtual screening strategies based on machine learning also have promising applications. In this study, we collected 2871 compound activity data against MMP-2 from the ChEMBL database and divided the training and test sets in a 3:1 ratio. Four machine learning algorithms were then selected to construct the classification models, and the best-performing model, i.e., the stacking-based fusion model with the highest AUC value in both training and test datasets, was used for the virtual screening of ZINC database. Next, we screened 17 potential MMP-2 inhibitors from the results predicted by the machine learning model via ADME/T analysis. The interactions between these compounds and the target protein were explored through molecular docking calculations, and the results showed that ZINC712249, ZINC4270723, and ZINC15858504 had lower binding free energies than the co-crystal ligand. To further examine the binding stability of the complexes, we performed molecular dynamics simulations and finally identified these three hits as the most promising natural products for MMP-2 inhibitors.

Identification of kinase inhibitors that rule out the CYP27B1-mediated activation of vitamin D: an integrated machine learning and structure-based drug designing approach

CYP27B1, a cytochrome P450-containing hydroxylase enzyme, converts vitamin D precursor calcidiol (25-hydroxycholecalciferol) to its active form calcitriol (1α,25(OH)₂D₃). Tyrosine kinase inhibitor such as imatinib is reported to interfere with the activation of vitamin D₃ by inhibiting CYP27B1 enzyme. Consequently, there is a decrease in the serum levels of active vitamin D that in turn may increase the relapse risk among the cancer patients treated with imatinib. Within this framework, the current study focuses on identifying other possible kinase inhibitors that may affect the calcitriol level in the body by inhibiting CYP27B1. To achieve this, we explored multiple machine learning approaches including support vector machine (SVM), random forest (RF), and artificial neural network (ANN) to identify possible CYP27B1 inhibitors from a pool of kinase inhibitors database. The most reliable classification model was obtained from the SVM approach with Matthews correlation coefficient of 0.82 for the external test set. This model was further employed for the virtual screening of kinase inhibitors from the binding database (DB), which tend to interfere with the CYP27B1-mediated activation of vitamin D. This screening yielded around 4646 kinase inhibitors that were further subjected to structure-based analyses using the homology model of CYP27B1, as the 3D structure of CYP27B1 complexed with heme was not available. Overall, five kinase inhibitors including two well-known drugs, i.e., AT7867 (Compound-2) and amitriptyline N-oxide (Compound-3), were found to interact with CYP27B1 in such a way that may preclude the conversion of vitamin D to its active form and hence testify the impairment of vitamin D activation pathway.

Multiple machine learning models combined with virtual screening and molecular docking to identify selective human ALDH1A1 inhibitors

Aldehyde dehydrogenases (ALDHs) are the enzymes of oxidoreductase family that are responsible for the aldehyde metabolism. The unbalanced expression of these enzymes may be associated with a variety of disease conditions including cancers. ALDH1A1 is one of the isoform of ALDHs majorly overexpressed in a variety of tumors and responsible for the anti-cancer drug resistance. This makes ALDH1A1 as a specific target to develop small molecule ALDH1A1 inhibitors for resistant cancer condition. Number of ALDH1A1 inhibitors have been developed and reported in the literature, but because of non-selectivity and inappropriate pharmacokinetic properties till now none of these have reached in the market for clinical use. Therefore, multiple machine learning models of different isoforms of ALDHs are integrated with in-silico techniques including virtual screening, docking, ADMET profiling, and MD simulation to identify selective ALDH1A1 inhibitors. Total ten selective ALDH1A1 inhibitors with diverse scaffolds and appropriate ADMET were identified that can be further developed as adjuvant therapy in cyclophosphamide and cisplatin resistance cancer.

Multiple machine learning, molecular docking, and ADMET screening approach for identification of selective inhibitors of CYP1B1

Cytochrome P4501B1 is a ubiquitous family protein that is majorly overexpressed in tumors and is responsible for biotransformation-based inactivation of anti-cancer drugs. This inactivation marks the cause of resistance to chemotherapeutics. In the present study, integrated in-silico approaches were utilized to identify selective CYP1B1 inhibitors. To achieve this objective, we initially developed different machine learning models corresponding to two isoforms of the CYP1 family i.e. CYP1A1 and CYP1B1. Subsequently, small molecule databases including ChemBridge, Maybridge, and natural compound library were screened from the selected models of CYP1B1 and CYP1A1. The obtained CYP1B1 inhibitors were further subjected to molecular docking and ADMET analysis. The selectivity of the obtained hits for CYP1B1 over the other isoforms was also judged with molecular docking analysis. Finally, two hits were found to be the most stable which retained key interactions within the active site of CYP1B1 after the molecular dynamics simulations. Novel compound with CYP-D9 and CYP-14 IDs were found to be the most selective CYP1B1 inhibitors which may address the issue of resistance. Moreover, these compounds can be considered as safe agents for further cell-based and animal model studies.Communicated by Ramaswamy H. Sarma.

Drug-target interaction prediction: A Bayesian ranking approach

BACKGROUND AND OBJECTIVE

In silico prediction of drug-target interactions (DTI) could provide valuable information and speed-up the process of drug repositioning - finding novel usage for existing drugs. In our work, we focus on machine learning algorithms supporting drug-centric repositioning approach, which aims to find novel usage for existing or abandoned drugs. We aim at proposing a per-drug ranking-based method, which reflects the needs of drug-centric repositioning research better than conventional drug-target prediction approaches.

METHODS

We propose Bayesian Ranking Prediction of Drug-Target Interactions (BRDTI). The method is based on Bayesian Personalized Ranking matrix factorization (BPR) which has been shown to be an excellent approach for various preference learning tasks, however, it has not been used for DTI prediction previously. In order to successfully deal with DTI challenges, we extended BPR by proposing: (i) the incorporation of target bias, (ii) a technique to handle new drugs and (iii) content alignment to take structural similarities of drugs and targets into account.

RESULTS

Evaluation on five benchmark datasets shows that BRDTI outperforms several state-of-the-art approaches in terms of per-drug nDCG and AUC. BRDTI results w.r.t. nDCG are 0.929, 0.953, 0.948, 0.897 and 0.690 for G-Protein Coupled Receptors (GPCR), Ion Channels (IC), Nuclear Receptors (NR), Enzymes (E) and Kinase (K) datasets respectively. Additionally, BRDTI significantly outperformed other methods (BLM-NII, WNN-GIP, NetLapRLS and CMF) w.r.t. nDCG in 17 out of 20 cases. Furthermore, BRDTI was also shown to be able to predict novel drug-target interactions not contained in the original datasets. The average recall at top-10 predicted targets for each drug was 0.762, 0.560, 1.000 and 0.404 for GPCR, IC, NR, and E datasets respectively.

CONCLUSIONS

Based on the evaluation, we can conclude that BRDTI is an appropriate choice for researchers looking for an in silico DTI prediction technique to be used in drug-centric repositioning scenarios. BRDTI Software and supplementary materials are available online at www.ksi.mff.cuni.cz/∼peska/BRDTI.

Classification models for safe drug molecules

Frequent failure of drug candidates during development stages remains the major deterrent for an early introduction of new drug molecules. The drug toxicity is the major cause of expensive late-stage development failures. An early identification/optimization of the most favorable molecule will naturally save considerable cost, time, human efforts and minimize animal sacrifice. (Quantitative) Structure Activity Relationships [(Q)SARs] represent statistically derived predictive models correlating biological activity (including desirable therapeutic effect and undesirable side effects) of chemicals (drugs/toxicants/environmental pollutants) with molecular descriptors and/or properties. (Q)SAR models which categorize the available data into two or more groups/classes are known as classification models. Numerous techniques of diverse nature are being presently employed for development of classification models. Though there is an increasing use of classification models for prediction of either biological activity or toxicity, the future trend will naturally be towards the development of classification models capable of simultaneous prediction of biological activity, toxicity, and pharmacokinetic parameters so as to accelerate development of bioavailable safe drug molecules.

Models for anti-tumor activity of bisphosphonates using refined topochemical descriptors

An in silico approach comprising of decision tree (DT), random forest (RF) and moving average analysis (MAA) was successfully employed for development of models for prediction of anti-tumor activity of bisphosphonates. A dataset consisting of 65 analogues of both nitrogen-containing and non-nitrogen-containing bisphosphonates was selected for the present study. Four refinements of eccentric distance sum topochemical index termed as augmented eccentric distance sum topochemical indices 1-4 [formula: see text] have been proposed so as to significantly augment discriminating power. Proposed topological indices (TIs) along with the exiting TIs (>1,400) were subsequently utilized for development of models for prediction of anti-tumor activity of bisphosphonates. A total of 43 descriptors of diverse nature, from a large pool of molecular descriptors, calculated through E-Dragon software (version 1.0) and an in-house computer program were selected for development of suitable models by employing DT, RF and MAA. DT identified two TIs as most important and classified the analogues of the dataset with an accuracy of 97% in training set and 90.7% in tenfold cross-validated set. Random forest correctly classified the analogues with an accuracy of 89.2%. Four independent models developed through MAA predicted the activity of analogues of the dataset with an accuracy of 87.6% to 89%. The statistical significance of proposed models was assessed through intercorrelation analysis, specificity, sensitivity and Matthew's correlation coefficient. The proposed models offer a vast potential for providing lead structures for development of potent anti-tumor agents for treatment of cancer that has spread to the bone.