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

Scaffold hopping based designing of selective ALDH1A1 inhibitors to overcome cyclophosphamide resistance: synthesis and biological evaluation

Aldehyde dehydrogenase 1A1 (ALDH1A1) is an isoenzyme that catalyzes the conversion of aldehydes to acids. However, the overexpression of ALDH1A1 in a variety of malignancies is the major cause of resistance to an anti-cancer drug, cyclophosphamide (CP). CP is a prodrug that is initially converted into 4-hydroxycyclophosphamide and its tautomer aldophosphamide, in the liver. These compounds permeate into the cell and are converted as active metabolites, i.e., phosphoramide mustard (PM), through spontaneous beta-elimination. On the other hand, the conversion of CP to PM is diverted at the level of aldophosphamide by converting it into inactive carboxyphosphamide using ALDH1A1, which ultimately leads to high drug inactivation and CP resistance. Hence, in combination with our earlier work on the target of resistance, i.e., ALDH1A1, we hereby report selective ALDH1A1 inhibitors. Herein, we selected a lead molecule from our previous virtual screening and implemented scaffold hopping analysis to identify a novel scaffold that can act as an ALDH1A1 inhibitor. This results in the identification of various novel scaffolds. Among these, on the basis of synthetic feasibility, the benzimidazole scaffold was selected for the design of novel ALDH1A1 inhibitors, followed by machine learning-assisted structure-based virtual screening. Finally, the five best compounds were selected and synthesized. All synthesized compounds were evaluated using in vitro enzymatic assay against ALDH1A1, ALDH2, and ALDH3A1. The results disclosed that three molecules A1, A2, and A3 showed significant selective ALDH1A1 inhibitory potential with an IC50 value of 0.32 μM, 0.55 μM, and 1.63 μM, respectively, and none of the compounds exhibits potency towards the other two ALDH isoforms i.e. ALDH2 and ALDH3A1. Besides, the potent compounds (A1, A2, and A3) have been tested for in vitro cell line assay in combination with mafosfamide (analogue of CP) on two cell lines i.e. A549 and MIA-PaCa-2. All three compounds show significant potency to reverse mafosfamide resistance by inhibiting ALDH1A1 against these cell lines.

Exploring Prognostic Biomarkers of Acute Myeloid Leukemia to Determine Its Most Effective Drugs from the FDA-Approved List through Molecular Docking and Dynamic Simulation

Acute myeloid leukemia (AML) is a blood cancer caused by the abnormal proliferation and differentiation of hematopoietic stem cells in the bone marrow. The actual genetic markers and molecular mechanisms of AML prognosis are unclear till today. This study used bioinformatics approaches for identifying hub genes and pathways associated with AML development to uncover potential molecular mechanisms. The expression profiles of RNA-Seq datasets, GSE68925 and GSE183817, were retrieved from the Gene Expression Omnibus (GEO) database. These two datasets were analyzed by GREIN to obtain differentially expressed genes (DEGs), which were used for performing the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, protein-protein interaction (PPI), and survival analysis. The molecular docking and dynamic simulation were performed to identify the most effective drug/s for AML from the drug list approved by the Food and Drug Administration (FDA). By integrating the two datasets, 238 DEGs were identified as likely to be affected by AML progression. GO enrichment analyses exhibited that the upregulated genes were mainly associated with inflammatory response (BP) and extracellular region (CC). The downregulated DEGs were involved in the T-cell receptor signalling pathway (BP), an integral component of the lumenal side of the endoplasmic reticulum membrane (CC) and peptide antigen binding (MF). The pathway enrichment analysis showed that the upregulated DEGs were mainly associated with the T-cell receptor signalling pathway. Among the top 15 hub genes, the expression levels of ALDH1A1 and CFD were associated with the prognosis of AML. Four FDA-approved drugs were selected, and a top-ranked drug was identified for each biomarker through molecular docking studies. The top-ranked drugs were further confirmed by molecular dynamic simulation that revealed their binding stability and confirmed their stable performance. Therefore, the drug compounds, enasidenib and gilteritinib, can be recommended as the most effective drugs against the ALDH1A1 and CFD proteins, respectively.

Raloxifene and bazedoxifene as selective ALDH1A1 inhibitors to ameliorate cyclophosphamide resistance: A drug repurposing approach

Cyclophosphamide (CP) is one of the most widely used anticancer drugs for various malignancies. However, its long-term use leads to ALDH1A1-mediated inactivation and subsequent resistance which necessitates the development of potential ALDH1A1 inhibitors. Currently, ALDH1A1 inhibitors from different chemical classes have been reported, but these failed to reach the market due to safety and efficacy problems. Developing a new treatment from the ground requires a huge amount of time, effort, and money, therefore it is worthwhile to improve CP efficacy by proposing better adjuvants as ALDH1A1 inhibitors. Herein, the database constituting the FDA-approved drugs with well-established safety and toxicity profiles was screened through already reported machine learning models by our research group. This model is validated for discriminating the ALDH1A1 inhibitors and non-inhibitors. Virtual screening protocol (VS) from this model identified four FDA-approved drugs, raloxifene, bazedoxifene, avanafil, and betrixaban as selective ALDH1A1 inhibitors. The molecular docking, dynamics, and water swap analysis also suggested these drugs to be promising ALDH1A1 inhibitors which were further validated for their CP resistance reversal potential by in-vitro analysis. The in-vitro enzymatic assay results indicated that raloxifene and bazedoxifene selectively inhibited the ALDH1A1 enzyme with IC₅₀ values of 2.35 and 4.41 μM respectively, whereas IC₅₀ values of both the drugs against ALDH2 and ALDH3A1 was >100 μM. Additional in-vitro stu = dies with well-reported ALDH1A1 overexpressing A549 and MIA paCa-2 cell lines suggested that mafosfamide sensitivity was further ameliorated by the combination of both raloxifene and bazedoxifene. Collectively, in-silico and in-vitro studies indicate raloxifene and bazedoxifene act as promising adjuvants with CP that may improve the quality of treatment for cancer patients with minimal toxicities.

TargIDe: a machine-learning workflow for target identification of molecules with antibiofilm activity against Pseudomonas aeruginosa

Bacterial biofilms are a source of infectious human diseases and are heavily linked to antibiotic resistance. Pseudomonas aeruginosa is a multidrug-resistant bacterium widely present and implicated in several hospital-acquired infections. Over the last years, the development of new drugs able to inhibit Pseudomonas aeruginosa by interfering with its ability to form biofilms has become a promising strategy in drug discovery. Identifying molecules able to interfere with biofilm formation is difficult, but further developing these molecules by rationally improving their activity is particularly challenging, as it requires knowledge of the specific protein target that is inhibited. This work describes the development of a machine learning multitechnique consensus workflow to predict the protein targets of molecules with confirmed inhibitory activity against biofilm formation by Pseudomonas aeruginosa. It uses a specialized database containing all the known targets implicated in biofilm formation by Pseudomonas aeruginosa. The experimentally confirmed inhibitors available on ChEMBL, together with chemical descriptors, were used as the input features for a combination of nine different classification models, yielding a consensus method to predict the most likely target of a ligand. The implemented algorithm is freely available at https://github.com/BioSIM-Research-Group/TargIDe under licence GNU General Public Licence (GPL) version 3 and can easily be improved as more data become available.

Targeting protein conformations with small molecules to control protein complexes

Dynamic protein complexes function in all cellular processes, from signaling to transcription, using distinct conformations that regulate their activity. Conformational switching of proteins can turn on or off their activity through protein-protein interactions, catalytic function, cellular localization, or membrane interaction. Recent advances in structural, computational, and chemical methodologies have enabled the discovery of small-molecule activators and inhibitors of conformationally dynamic proteins by using a more rational design than a serendipitous screening approach. Here, we discuss such recent examples, focusing on the mechanism of protein conformational switching and its regulation by small molecules. We emphasize the rational approaches to control protein oligomerization with small molecules that offer exciting opportunities for investigation of novel biological mechanisms and drug discovery.

Role of Genetic Polymorphisms in Drug-Metabolizing Enzyme-Mediated Toxicity and Pharmacokinetic Resistance to Anti-Cancer Agents: A Review on the Pharmacogenomics Aspect

The inter-individual differences in cancer susceptibility are somehow correlated with the genetic differences that are caused by the polymorphisms. These genetic variations in drug-metabolizing enzymes/drug-inactivating enzymes may negatively or positively affect the pharmacokinetic profile of chemotherapeutic agents that eventually lead to pharmacokinetic resistance and toxicity against anti-cancer drugs. For instance, the CYP1B1*3 allele is associated with CYP1B1 overexpression and consequent resistance to a variety of taxanes and platins, while 496T>G is associated with lower levels of dihydropyrimidine dehydrogenase, which results in severe toxicities related to 5-fluorouracil. In this context, a pharmacogenomics approach can be applied to ascertain the role of the genetic make-up in a person's response to any drug. This approach collectively utilizes pharmacology and genomics to develop effective and safe medications that are devoid of resistance problems. In addition, recently reported genomics studies revealed the impact of many single nucleotide polymorphisms in tumors. These studies emphasized the importance of single nucleotide polymorphisms in drug-metabolizing enzymes on the effect of anti-tumor drugs. In this review, we discuss the pharmacogenomics aspect of polymorphisms in detail to provide an insight into the genetic manipulations in drug-metabolizing enzymes that are responsible for pharmacokinetic resistance or toxicity against well-known anti-cancer drugs. Special emphasis is placed on different deleterious single nucleotide polymorphisms and their effect on pharmacokinetic resistance. The information provided in this report may be beneficial to researchers, especially those who are working in the field of biotechnology and human genetics, in rationally manipulating the genetic information of patients with cancer who are undergoing chemotherapy to avoid the problem of pharmacokinetic resistance/toxicity associated with drug-metabolizing enzymes.

Multi-stage virtual screening of natural products against p38α mitogen-activated protein kinase: predictive modeling by machine learning, docking study and molecular dynamics simulation

p38α is a mitogen-activated protein kinase (MAPK), and the signaling pathways involved are closely related to the inflammation, apoptosis and differentiation of cells, which also makes it an attractive target for drug discovery. With the high efficiency and low cost, virtual screening technology is becoming an indispensable part of drug development. In this study, a novel multi-stage virtual screening method based on machine learning, molecular docking and molecular dynamics simulation was developed to identify p38α MAPK inhibitors from natural products in ZINC database, which improves the prediction accuracy by considering and utilizing both ligand and receptor information compared to any individual approach. Ultimately, we screened out two candidate inhibitors with acceptable ADMET properties (ZINC4260400 and ZINC8300300). Among the generated machine learning models, Random Forest (RF) and Support Vector Machine (SVM) performed better, with the area under the receiver operating characteristic curve (AUC) values of 0.932 and 0.931 on the test set, as well as 0.834 and 0.850 on the external validation set. In addition, the results of molecular docking and ADMET prediction showed that two compounds with appropriate pharmacokinetic properties had binding free energies less than −8.0 kcal/mol for the target protein, and the results of molecular dynamics simulations further confirmed that they were stable during the process of inhibition.

ALDH1A1 in Cancers: Bidirectional Function, Drug Resistance, and Regulatory Mechanism

Aldehyde dehydrogenases 1 family member A1(ALDH1A1) gene codes a cytoplasmic enzyme and shows vital physiological and pathophysiological functions in many areas. ALDH1A1 plays important roles in various diseases, especially in cancers. We reviewed and summarized representative correlative studies and found that ALDH1A1 could induce cancers via the maintenance of cancer stem cell properties, modification of metabolism, promotion of DNA repair. ALDH1A1 expression is regulated by several epigenetic processes. ALDH1A1 also acted as a tumor suppressor in certain cancers. The detoxification of ALDH1A1 often causes chemotherapy failure. Currently, ALDH1A1-targeted therapy is widely used in cancer treatment, but the mechanism by which ALDH1A1 regulates cancer development is not fully understood. This review will provide insight into the status of ALDH1A1 research and new viewpoint for cancer therapy.