Molecular modeling approaches to address drug-metabolizing enzymes (DMEs) mediated chemoresistance: a review

DPP9 regulates NQO1 and ROS to promote resistance to chemotherapy in liver cancer cells

Chemotherapy has been the standard treatment for liver cancer. However, intrinsic or acquired drug resistance remains a major barrier to successful treatment. At present, the underlying molecular mechanisms of chemoresistance in liver cancer have not been elucidated. Dipeptidyl peptidase 9 (DPP9) is a member of the dipeptidyl peptidase IV family that has been found to be highly expressed in a variety of tumors, including liver cancer. It is unclear whether DPP9 affects chemoresistance in liver cancer. In this study, we find that DPP9 weakens the responses of liver cancer cells to chemotherapy drugs by up-regulating NQO1 and inhibiting intracellular ROS levels. In terms of mechanism, DPP9 inhibits ubiquitin-mediated degradation of NRF2 protein by binding to KEAP1, up-regulates NRF2 protein levels, promotes mRNA transcription of NQO1, and inhibits intracellular ROS levels. In addition, the NQO1 inhibitor dicoumarol can enhance the efficacy of chemotherapy drugs in liver cancer cells. Collectively, our findings suggest that inhibiting DPP9/NQO1 signaling can serve as a potential therapeutic strategy for liver cancer.

Breaking the Barriers of Therapy Resistance: Harnessing Ferroptosis for Effective Hepatocellular Carcinoma Therapy

Ferroptosis is a type of cell death that relies on iron and is distinguished by the occurrence of lipid peroxidation and the buildup of reactive oxygen species. Ferroptosis has been demonstrated to have a significant impact on the advancement and resistance to treatment of hepatocellular carcinoma (HCC), thereby highlighting its potential as a viable therapeutic target. Ferroptosis was observed in HCC tissues in contrast to normal liver tissue. The inhibition of ferroptosis has been found to increase the viability of HCC cells and decrease their susceptibility to various anticancer therapies, including chemotherapy, radiotherapy, and immune checkpoint blockade. The administration of drugs that directly modulate ferroptosis regulators or induce excessive production of lipid-reactive oxygen species has demonstrated the potential to enhance the responsiveness of drug-resistant HCC cells to treatment. However, the precise mechanism underlying this phenomenon remains ambiguous. This review presents a comprehensive overview of the crucial role played by ferroptosis in enhancing the efficacy of treatment for hepatocellular carcinoma (HCC). The main aim of this study is to examine the feasibility of utilizing ferroptosis as a therapeutic approach to improve the efficacy of HCC treatment and overcome drug resistance.

Deciphering the Molecular Mechanisms of Reactive Metabolite Formation in the Mechanism-Based Inactivation of Cytochrome p450 1B1 by 8-Methoxypsoralen and Assessing the Driving Effect of phe268

This study provides a comprehensive computational exploration of the inhibitory activity and metabolic pathways of 8-methoxypsoralen (8-MP), a furocoumarin derivative used for treating various skin disorders, on cytochrome P450 (P450). Employing quantum chemical DFT calculations, molecular docking, and molecular dynamics (MD) simulations analyses, the biotransformation mechanisms and the active site binding profile of 8-MP in CYP1B1 were investigated. Three plausible inactivation mechanisms were minutely scrutinized. Further analysis explored the formation of reactive metabolites in subsequent P450 metabolic processes, including covalent adduct formation through nucleophilic addition to the epoxide, 8-MP epoxide hydrolysis, and non-CYP-catalyzed epoxide ring opening. Special attention was paid to the catalytic effect of residue Phe268 on the mechanism-based inactivation (MBI) of P450 by 8-MP. Energetic profiles and facilitating conditions revealed a slight preference for the C4'=C5' epoxidation pathway, while recognizing a potential kinetic competition with the 8-OMe demethylation pathway due to comparable energy demands. The formation of covalent adducts via nucleophilic addition, particularly by phenylalanine, and the generation of potentially harmful reactive metabolites through autocatalyzed ring cleavage are likely to contribute significantly to P450 metabolism of 8-MP. Our findings highlight the key role of Phe268 in retaining 8-MP within the active site of CYP1B1, thereby facilitating initial oxygen addition transition states. This research offers crucial molecular-level insights that may guide the early stages of drug discovery and risk assessment related to the use of 8-MP.

Scaffold hopping for designing of potent and selective CYP1B1 inhibitors to overcome docetaxel resistance: synthesis and evaluation

Cytochrome P450 1B1, a tumor-specific overexpressed enzyme, significantly impairs the pharmacokinetics of several commonly used anticancer drugs including docetaxel, paclitaxel and cisplatin, leading to the problem of resistance to these drugs. Currently, there is no CYP1B1 inhibition-based adjuvant therapy available to treat this resistance problem. Hence, in the current study, exhaustive in-silico studies including scaffold hopping followed by molecular docking, three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular dynamics and free energy perturbation studies were carried out to identify potent and selective CYP1B1 inhibitors. Initially, scaffold hopping analysis was performed against a well-reported potent and selective CYP1B1 inhibitor (i.e. compound 3n). A total of 200 scaffolds were identified along with their shape and field similarity scores. The top three scaffolds were further selected on the basis of these scores and their synthesis feasibility to design some potent and selective CYP1B1 inhibitors using the aforementioned in-silico techniques. Designed molecules were further synthesized to evaluate their CYP1B1 inhibitory activity and docetaxel resistance reversal potential against CYP1B1 overexpressed drug resistance MCF-7 cell line. In-vitro results indicated that compounds 2a, 2c and 2d manifested IC50 values for CYP1B1 ranging from 0.075, 0.092 to 0.088 μM with at least 10-fold selectivity. At low micromolar concentrations, compounds 1e, 1f, 2a and 2d exhibited promising cytotoxic effects in the docetaxel-resistant CYP1B1 overexpressed MCF-7 cell line. In particular, compound 2a is most effective in reversing the resistance with IC50 of 29.0 ± 3.6 μM. All of these discoveries could pave the way for the development of adjuvant therapy capable of overcoming CYP1B1-mediated resistance.Communicated by Ramaswamy H. Sarma.

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.

CYP1A inhibitors: Recent progress, current challenges, and future perspectives

Mammalian cytochrome P450 1A (CYP1A) are key phase I xenobiotic-metabolizing enzymes that play a distinctive role in metabolic activation or metabolic clearance of a variety of procarcinogens, drugs, and endogenous substances. Human CYP1A subfamily contains two members (hCYP1A1 and hCYP1A2), which are known to catalyze the oxidative activation of some environmental procarcinogens into carcinogenic species. Increasing evidence has demonstrated that CYP1A inhibitor therapies are promising strategies for cancer chemoprevention or overcoming CYP1A-associated drug toxicity and resistance. Herein, we reviewed recent advances in the discovery and characterization of hCYP1A inhibitors, from the discovery approaches to structural features and biomedical applications of hCYP1A inhibitors. The inhibition potentials, inhibition modes, and inhibition constants of all reported hCYP1A inhibitors are comprehensively summarized. Meanwhile, the structural features and structure-activity relationships of different classes of hCYP1A1 and hCYP1A2 inhibitors are analyzed and discussed in depth. Furthermore, the major challenges and future directions for this field are presented and highlighted. Collectively, the information and knowledge presented here will strongly facilitate the researchers to discover and develop more efficacious CYP1A inhibitors for specific purposes, such as chemo-preventive agents or as tool molecules in hCYP1A-related fundamental studies.

3D-QSAR assisted identification of selective CYP1B1 inhibitors: an effective bioisosteric replacement/molecular docking/electrostatic complementarity analysis

Cytochrome P450-1B1 is a majorly overexpressed drug-metabolizing enzyme in tumors and is responsible for inactivation and subsequent resistance to a variety of anti-cancer drugs, i.e., docetaxel, tamoxifen, and cisplatin. In the present study, a 3D quantitative structure-activity relationship (3D-QSAR) model has been constructed for the identification, design, and optimization of novel CYP1B1 inhibitors. The model has been built using a set of 148 selective CYP1B1 inhibitors. The developed model was evaluated based on certain statistical parameters including q2 and r2 which showed the acceptable predictive and descriptive capability of the generated model. The developed 3D-QSAR model assisted in understanding the key molecular fields which were firmly related to the selective CYP1B1 inhibition. A theoretic approach for the generation of new lead compounds with optimized CYP1B1 receptor affinity has been performed utilizing bioisosteric replacement analysis. These generated molecules were subjected to a developed 3D-QSAR model to predict the inhibitory activity potentials. Furthermore, these compounds were scrutinized through the activity atlas model, molecular docking, electrostatic complementarity, molecular dynamics, and waterswap analysis. The final hits might act as selective CYP1B1 inhibitors which could address the issue of resistance. This 3D-QSAR includes several chemically diverse selective CYP1B1 receptor ligands and well accounts for the individual ligand's inhibition affinities. These features of the developed 3D-QSAR model will ensure future prospective applications of the model to speed up the identification of new potent and selective CYP1B1 receptor ligands.

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.

In silico guided designing of optimized benzochalcones derivatives as potent CYP1B1 inhibitors: An integrated in vitro and ONIOM study

Cytochrome P4501B1 (CYP1B1) is reported to be overexpressed in various malignancies including ovarian, lung, lymph, and breast cancers. The overexpression of this enzyme is accountable for the biotransformation-based inactivation of some anti-cancer drugs i.e. Docetaxel, Paclitaxel, and Cisplatin. To circumvent solutions to this issue, the current study reports some optimized derivatives of benzochalcone as selective CYP1B1 inhibitors. The optimized derivatives were screened using some structure-based drug-designing approaches including molecular docking and molecular dynamics. The implemented approaches revealed that all the designed molecules demonstrated not only essential interactions with key amino acid residues but also maintained stability within the active site of CYP1B1. Furthermore, to validate the in-silico results and develop a SAR, the designed molecules were subsequently synthesized and tested for their ability to selectively inhibit CYP1B1 over CYP1A1 using well established EROD assay. This assay results suggested that compounds 1(c), 1(d), and 1(e) are eightfold more selective CYP1B1 inhibitors over CYP1A1 with IC₅₀ values ranging from 0.06 to 0.09 μM respectively. Among these, compound 1(d) manifested potent inhibitory activity i.e. IC₅₀ of 0.06 μM with 24 folds selectivity over 1A1. To have a better insight into the binding pattern of 1(d) within CYP1B1 and precisely compute binding affinity for 1(d)-CYP1B1 complex, one of the advanced QM/MM approaches i.e. ONIOM has been implemented. Where 1(d)-CYP1B1 complex conferred comparable binding affinity in terms of ΔG (kcal/mol) with that of ANF-CYP1B1 complex. This research could provide a suitable starting point for the development of more potent multi-functional compounds with CYP1B1 inhibitory activity.

Dihydropyrimidine Dehydrogenase-Mediated Resistance to 5-Fluorouracil: Mechanistic Investigation and Solution

5-Fluorouracil (5-FU) is one of the most widely used chemotherapeutics for the treatment of cancers associated with the aerodigestive tract, breast, and colorectal system. The efficacy of 5-FU is majorly affected by dihydropyrimidine dehydrogenase (DPD) as it degrades more than 80% of administered 5-FU into an inactive metabolite, dihydrofluorouracil. Herein we discuss the molecular mechanism of this inactivation by analyzing the interaction pattern and electrostatic complementarity of the DPD-5-FU complex. The basis of DPD overexpression in cancer cell lines due to significantly distinct levels of the miRNAs (miR-134, miR-27b, and miR-27a) compared to normal cells has also been outlined. Additionally, some kinases including sphingosine kinase 2 (SphK2) have been reported to correlate with DPD expression. Currently, to address this problem various strategies are reported in the literature, including 5-FU analogues (bypass the DPD-mediated inactivation), DPD downregulators (regulate the DPD expression levels in tumors), inhibitors (as promising adjuvants), and formulation development loaded with 5-FU (liposomes, nanoparticles, nanogels, etc.), which are briefly discussed in this Review.

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.

3D-QSAR and scaffold hopping based designing of benzo[d]ox-azol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one derivatives as selective aldehyde dehydrogenase 1A1 inhibitors: Synthesis and biological evaluation

Aldehyde dehydrogenase 1 (ALDH1A1), an oxidoreductase class of enzymes, is overexpressed in various types of cancer cell lines and is the major cause of resistance to the Food and Drug Administration (FDA)-approved drug, cyclophosphamide (CP). In cancer conditions, CP undergoes a sequence of biotransformations to form an active metabolite, aldophosphamide, which further biotransforms to its putative cytotoxic metabolite, phosphoramide mustard. However, in resistant cancer conditions, aldophosphamide is converted into its inactive metabolite, carboxyphosphamide, via oxidation with ALDH1A1. Herein, to address the issue of ALDH1A1 mediated CP resistance, we report a series of benzo[d]oxazol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one derivatives as selective ALDH1A1 inhibitors. These inhibitors were designed using a validated 3D-quantitative structure activity relationship (3D-QSAR) model coupled with scaffold hopping. The 3D-QSAR model was developed using reported indole-2,3-diones based ALDH1A1 inhibitors, which provided field points in terms of electrostatic, van der Waals and hydrophobic potentials required for selectively inhibiting ALDH1A1. The most selective indole-2,3-diones-based compound, that is, cmp 3, was further considered for scaffold hopping. Two top-ranked bioisosteres, that is, benzo[d]oxazol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one, were selected for designing new inhibitors by considering the field pattern of 3D-QSAR. All designed molecules were mapped perfectly on the 3D-QSAR model and found to be predictive with good inhibitory potency (pIC₅₀ range: 7.5-6.8). Molecular docking was carried out for each designed molecule to identify key interactions that are required for ALDH1A1 inhibition and to authenticate the 3D-QSAR result. The top five inhibitor-ALDH1A1 complexes were also submitted for molecular dynamics simulations to access their stability. In vitro enzyme assays of 21 compounds suggested that these compounds are selective toward ALDH1A1 over the other two isoforms, that is, ALDH2 and ALDH3A1. All the compounds were found to be at least three and two times more selective toward ALDH1A1 over ALDH2 and ALDH3A1, respectively. All the compounds showed an IC₅₀ value in the range of 0.02-0.80 μM, which indicates the potential for these to be developed as adjuvant therapy for CP resistance.