Cheminformatics models based on machine learning approaches for design of USP1/UAF1 abrogators as anticancer agents

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.

Ubiquitin-specific protease 1 inhibition sensitizes hepatocellular carcinoma cells to doxorubicin by ubiquitinated proliferating cell nuclear antigen-mediated attenuation of stemness

Currently, resistance to the chemotherapeutic agent doxorubicin (Dox) in hepatocellular carcinoma (HCC) cells is an obstacle in developing effective Dox-targeted clinical therapies. Ubiquitin-specific protease 1 (USP1) plays a crucial role in the progression of multiple cancers. In this study, the purpose was to investigate the effect of USP1 depletion with chemotherapeutant Dox on the HCC cells. Flow cytometry was used to detect the ratio of apoptosis. The expression levels of selected proteins were evaluated by western blotting. In addition, the expression of genes was quantitated by quantitative real-time PCR assay. Coimmunoprecipitation was performed to confirm the interaction between USP1 and proliferating cell nuclear antigen (PCNA). Sphere formation assay was carried out to investigate the cancer stemness. Subcutaneous xenograft and orthotopic liver tumor models were established to examine the growth of tumor. Knockdown of USP1 increased the rate of Dox-induced apoptosis in stem-like and nonstem-like HCC cells. The combination of Dox and the USP1 inhibitor SJB3-019A (SJB3) markedly enhanced apoptosis in the primary liver carcinoma/PRF/5 and MHCC-97H cell lines. Notably, Dox/SJB3-induced tumor inhibition was further determined in vivo using a xenograft and orthotopic liver tumor model. Mechanically, USP1 inhibition via SJB3 or short hairpin RNA significantly decreased cancer stemness, including sphere formation ability and the expression of Nanog, Sox2, and c-Myc. The sensitization of HCC to Dox by SJB3 is attributed to the upregulation of PCNA ubiquitylation. Thus, genetic or pharmacological inhibition of USP1 restored the sensitivity of HCC cells to Dox in vitro and in vivo , representing a new potential therapeutic strategy for HCC.

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.

Machine Learning From Molecular Dynamics Trajectories to Predict Caspase-8 Inhibitors Against Alzheimer's Disease

Alzheimer’s disease (AD) is a neurodegenerative disorder in which the death of brain cells takes place leading to loss of memory and decreased cognitive ability. AD is a leading cause of death worldwide and is progressive in nature with symptoms worsening over time. Machine learning–based computational predictive models based on 2D and 3D descriptors have been effective in identifying potential active compounds. However, the use of data from molecular dynamics (MD) trajectories for training machine learning models still needs to be explored. In the present study, descriptors have been extracted from the MD trajectories of caspase-8 ligand complexes to train models using artificial neural networks and random forest algorithms. Caspase-8 plays a key role in causing AD by cleaving amyloid precursor proteins during apoptosis leading to increased formation of the amyloid-beta peptide. A total of 43 ligands were docked using the glide module of Schrodinger software, and short MD simulations of 10 ns were performed for the calculation of MD descriptors. The MD descriptors were also combined with the 2D and 3D descriptors of chemical compounds, and individual descriptor based as well as combination models were generated. This study demonstrated that MD descriptors could be effectively used for the characterization of bioactive compounds along with lead prioritization and optimization.

Computational models for the prediction of adverse cardiovascular drug reactions

Background

Predicting adverse drug reactions (ADRs) has become very important owing to the huge global health burden and failure of drugs. This indicates a need for prior prediction of probable ADRs in preclinical stages which can improve drug failures and reduce the time and cost of development thus providing efficient and safer therapeutic options for patients. Though several approaches have been put forward for in silico ADR prediction, there is still room for improvement.

Methods

In the present work, we have used machine learning based approach for cardiovascular (CV) ADRs prediction by integrating different features of drugs, biological (drug transporters, targets and enzymes), chemical (substructure fingerprints) and phenotypic (therapeutic indications and other identified ADRs), and their two and three level combinations. To recognize quality and important features, we used minimum redundancy maximum relevance approach while synthetic minority over-sampling technique balancing method was used to introduce a balance in the training sets.

Results

This is a rigorous and comprehensive study which involved the generation of a total of 504 computational models for 36 CV ADRs using two state-of-the-art machine-learning algorithms: random forest and sequential minimization optimization. All the models had an accuracy of around 90% and the biological and chemical features models were more informative as compared to the models generated using chemical features.

Conclusions

The results obtained demonstrated that the predictive models generated in the present study were highly accurate, and the phenotypic information of the drugs played the most important role in drug ADRs prediction. Furthermore, the results also showed that using the proposed method, different drugs properties can be combined to build computational predictive models which can effectively predict potential ADRs during early stages of drug development.

Electronic supplementary material

The online version of this article (10.1186/s12967-019-1918-z) contains supplementary material, which is available to authorized users.

Design of novel amyloid β aggregation inhibitors using QSAR, pharmacophore modeling, molecular docking and ADME prediction

The inhibition of abnormal amyloid β (Aβ) aggregation has been regarded as a good target to control Alzheimer's disease. The present study adopted 2D-QSAR, HQSAR and 3D QSAR (CoMFA & CoMSIA) modeling approaches to identify the structural and physicochemical requirements for the potential Aβ aggregation inhibition. A structure-based molecular docking technique is utilized to approve the features that are obtained from the ligand-based techniques on 30 curcumin derivatives. The combined outputs were then used to screen the modified 10 compounds. The 2D QSAR model on curcumin derivatives gave statistical values R2 = 0.9086 and SEE = 0.1837. The model was further confirmed by Y-randomization test and Applicability domain analysis by the standardization approach. The HQSAR study (Q2 = 0.615, Rncv 2 = 0.931, Rpred 2 = 0.956) illustrated the important molecular fingerprints for inhibition. Contour maps of 3D QSAR models, CoMFA (Q2 = 0.687, Rncv 2 = 0.787, Rpred 2 = 0.731) and CoMSIA (Q2 = 0.743, Rncv 2 = 0.972, Rpred 2 = 0.713), depict that the models are robust and provide explanation of the important features, like steric, electrostatic and hydrogen bond acceptor, which play important role for interaction with the receptor site cavity. The molecular docking study of the curcumin derivatives elucidates the important interactions between the amino acid residues at the catalytic site of the receptor and the ligands, indicating the structural requirements of the inhibitors. The ligand-receptor interactions of top hits were analyzed to explore the pharmacophore features of Aβ aggregation inhibition. The Aβ aggregation inhibitory activities of novel chemical entities were then obtained through inverse QSAR. The newly designed molecules were further screened through machine learning, prediction of toxicity and nature of metabolism to get the proposed six lead compounds.

Gene silencing of USP1 by lentivirus effectively inhibits proliferation and invasion of human osteosarcoma cells

Osteosarcoma is the most frequent malignant bone tumor, affecting the extremities of adolescents and young adults. Ubiquitin-specific protease 1 (USP1) plays a critical role in many cellular processes including proteasome degradation, chromatin remodeling and cell cycle regulation. In the present study, we discovered that USP1 was overexpressed in 26 out of 30 osteosarcoma tissues compared to cartilage tumor tissues and normal bone tissues. We then constructed a lentiviral vector mediating RNA interference (RNAi) targeting USP1 and demonstrated that it significantly suppressed the mRNA and protein expression of the USP1 gene in U2OS cells. Knockdown of USP1 inhibited the growth and colony-forming, as well as significantly reduced the invasiveness of U2OS cells. Western blot analysis indicated that suppression of USP1 downregulated the expression of many proteins including SIK2, MMP-2, GSK-3β, Bcl-2, Stat3, cyclin E1, Notch1, Wnt-1 and cyclin A1. Most of these proteins are associated with tumor genesis and development. RNAi of SIK2 significantly decreased SIK2 protein expression and inhibited the ability of forming colonies, as well as induced apoptosis and reduced the invasiveness of U2OS cells. Collectively, our results suggest that silencing USP1 inhibits cell proliferation and invasion in U2OS cells. Therefore, USP1 may provide a novel therapeutic target for the treatment of osteosarcoma.

Integrating network, sequence and functional features using machine learning approaches towards identification of novel Alzheimer genes

Background

Alzheimer’s disease (AD) is a complex progressive neurodegenerative disorder commonly characterized by short term memory loss. Presently no effective therapeutic treatments exist that can completely cure this disease. The cause of Alzheimer’s is still unclear, however one of the other major factors involved in AD pathogenesis are the genetic factors and around 70 % risk of the disease is assumed to be due to the large number of genes involved. Although genetic association studies have revealed a number of potential AD susceptibility genes, there still exists a need for identification of unidentified AD-associated genes and therapeutic targets to have better understanding of the disease-causing mechanisms of Alzheimer’s towards development of effective AD therapeutics.

Results

In the present study, we have used machine learning approach to identify candidate AD associated genes by integrating topological properties of the genes from the protein-protein interaction networks, sequence features and functional annotations. We also used molecular docking approach and screened already known anti-Alzheimer drugs against the novel predicted probable targets of AD and observed that an investigational drug, AL-108, had high affinity for majority of the possible therapeutic targets. Furthermore, we performed molecular dynamics simulations and MM/GBSA calculations on the docked complexes to validate our preliminary findings.

Conclusions

To the best of our knowledge, this is the first comprehensive study of its kind for identification of putative Alzheimer-associated genes using machine learning approaches and we propose that such computational studies can improve our understanding on the core etiology of AD which could lead to the development of effective anti-Alzheimer drugs.

Electronic supplementary material

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

Integrated machine learning, molecular docking and 3D-QSAR based approach for identification of potential inhibitors of trypanosomal N-myristoyltransferase

Integrated in silico approaches for the identification of antitrypanosomal inhibitors.