Rational drug repurposing for cancer by inclusion of the unbiased molecular dynamics simulation in the structure-based virtual screening approach: Challenges and breakthroughs

Pyruvate kinase deficiency and PKLR gene mutations: Insights from molecular dynamics simulation analysis

Pyruvate kinase deficiency is a rare hereditary erythrocyte enzyme disease caused by mutations in the pyruvate kinase liver and red blood cell gene. The clinical presentations of pyruvate kinase deficiency are significantly heterogeneous, ranging from just mild anemia to hemolytic crisis or even death. The proband in our study was a 2-year-old girl for severe skin and scleral icterus with progressive aggravation. We collected the family's data for further analysis. Whole exome genome sequencing of the pedigree revealed a novel compound heterozygous mutation, c.1097del (p.P366Lfs*12) and c.1493G > A (p.R498H), in the pyruvate kinase liver and red blood cell gene. Furthermore, molecular dynamics simulations were employed to uncover differences between the wild type and mutant pyruvate kinase liver and red blood cell proteins, focusing on structural stability, protein flexibility, secondary structure, and overall conformation. The combined bioinformatic tools were also utilised to assess the effects of the missense mutation on protein function. Thereafter, wild type and mutant plasmids were constructed and transfected into 293T cells, and Western blot assay was conducted to validate the impact of the mutations on the expression of pyruvate kinase liver and red blood cell protein. The data presented in our study enriches the genotype database and provides evidence for genetic counseling and molecular diagnosis of pyruvate kinase deficiency.

Structural Dynamics-Driven Discovery of Anticancer and Antimetastatic Effects of Diltiazem and Glibenclamide Targeting Urokinase Receptor

Diltiazem and glibenclamide are commonly used hypotensive and antidiabetic drugs. This study reports the discovery of the potential antitumor and antimetastatic effects of these two drugs using a structural dynamics-driven virtual screening targeting urokinase receptor (uPAR). Owing to uPAR's high flexibility, currently resolved crystal structures of uPAR, all in ligand-bound states, provide limited representations of its physiological conformation. To improve the accuracy of screening, we performed a long-timescale molecular dynamics simulation and obtained the representative conformations of apo-uPAR as the targets for our screening. Experimentally, we demonstrated that diltiazem and glibenclamide bound uPAR with K_D values in the micromolar range. In addition, both compounds effectively suppressed tumor growth and metastasis in a uPAR-dependent manner in vitro and in vivo. This work not only provides two potent uPAR inhibitors but also reports a proof-of-concept study on the potential off-label antitumor and antimetastatic uses of diltiazem and glibenclamide.

Complete reconstruction of dasatinib unbinding pathway from c-Src kinase by supervised molecular dynamics simulation method; assessing efficiency and trustworthiness of the method

Over the past years, rational drug design has gained lots of attention since employing it gave the world targeted therapy and more effective treatment solutions. Structure-based drug design (SBDD) is an excellent tool in rational drug design that takes advantage of accurate methods such as unbiased molecular dynamics (UMD) simulation for designing and optimizing molecular entities by understanding the binding and unbinding pathways of the binders. Supervised molecular dynamics (SuMD) simulation is a branch of UMD in which long-duration simulations are turned into short simulations, called replica, and a specific parameter is monitored throughout the simulation. In this work, we utilized this strategy to reconstruct the unbinding pathway of the anticancer drug dasatinib from its target protein, the c-Src kinase. Several unbinding events with valuable details were achieved. Then, to assess the efficiency and trustworthiness of the SuMD method, the unbinding pathway was also reconstructed by conventional UMD simulation, which uncovered some of the limitations of this method, such as limited sampling of the active site and finding the metastable states in the unbinding pathway. Furthermore, in times like these, when the world is desperate to find treatments for the Covid-19 disease, we think these methods are of exceptional value.Communicated by Ramaswamy H. Sarma.

Mitochondrial adaptation in cancer drug resistance: prevalence, mechanisms, and management

Drug resistance represents a major obstacle in cancer management, and the mechanisms underlying stress adaptation of cancer cells in response to therapy-induced hostile environment are largely unknown. As the central organelle for cellular energy supply, mitochondria can rapidly undergo dynamic changes and integrate cellular signaling pathways to provide bioenergetic and biosynthetic flexibility for cancer cells, which contributes to multiple aspects of tumor characteristics, including drug resistance. Therefore, targeting mitochondria for cancer therapy and overcoming drug resistance has attracted increasing attention for various types of cancer. Multiple mitochondrial adaptation processes, including mitochondrial dynamics, mitochondrial metabolism, and mitochondrial apoptotic regulatory machinery, have been demonstrated to be potential targets. However, recent increasing insights into mitochondria have revealed the complexity of mitochondrial structure and functions, the elusive functions of mitochondria in tumor biology, and the targeting inaccessibility of mitochondria, which have posed challenges for the clinical application of mitochondrial-based cancer therapeutic strategies. Therefore, discovery of both novel mitochondria-targeting agents and innovative mitochondria-targeting approaches is urgently required. Here, we review the most recent literature to summarize the molecular mechanisms underlying mitochondrial stress adaptation and their intricate connection with cancer drug resistance. In addition, an overview of the emerging strategies to target mitochondria for effectively overcoming chemoresistance is highlighted, with an emphasis on drug repositioning and mitochondrial drug delivery approaches, which may accelerate the application of mitochondria-targeting compounds for cancer therapy.

Pharmacological inhibition of KDM5A for cancer treatment

Lysine-specific demethylase 5A (KDM5A, also named RBP2 or JARID1A) is a demethylase that can remove methyl groups from histones H3K4me1/2/3. It is aberrantly expressed in many cancers, where it impedes differentiation and contributes to cancer cell proliferation, cell metastasis and invasiveness, drug resistance, and is associated with poor prognosis. Pharmacological inhibition of KDM5A has been reported to significantly attenuate tumor progression in vitro and in vivo in a range of solid tumors and acute myeloid leukemia. This review will present the structural aspects of KDM5A, its role in carcinogenesis, a comparison of currently available approaches for screening KDM5A inhibitors, a classification of KDM5A inhibitors, and its potential as a drug target in cancer therapy.

Reconstruction of the binding pathway of an anti-HIV drug, Indinavir, in complex with the HTLV-1 protease using unaggregated unbiased molecular dynamics simulation

Retroviruses are a growing concern for the health of human beings, and one of the dangerous members of this family is the Human T-cell Leukemia Virus 1 (HTLV-1) virus. It has affected more than 20 million people so far, and since there are no registered treatments against it yet, urgent treatment solutions are needed. One of the most promising drug targets to fight this virus is the protease enzyme of the virus's protein machinery. In this study, by utilizing a computational method called Unaggregated Unbiased Molecular Dynamics (UUMD), we reconstructed the binding pathway of a HTLV-1 protease inhibitor, Indinavir, to find the details of the binding pathway, the influential residues, and also the stable states of the binding pathway. We achieved the native conformation of the inhibitor in 6 rounds, 360 replicas by performing over 4 micro-seconds of UMD simulations. We found 3 Intermediate states between the solvated state and the native conformation state in the binding pathway. We also discovered that aromatic residues such as Trp98 and Trp98', catalytic residues Asp32 and Asp32', and the flap region's residues have the most influential roles in the binding pathway and also have the most contribution to the total interaction energies. We believe that the details found in this study would be a great guide for developing new treatment solutions against the HTLV-1 virus by inhibiting the HTLV-1 protease.

Comparative analysis of the unbinding pathways of antiviral drug Indinavir from HIV and HTLV1 proteases by supervised molecular dynamics simulation

Determining the unbinding pathways of potential small molecule compounds from their target proteins is of great significance for designing efficacious treatment solutions. One of these potential compounds is the approved HIV-1 protease inhibitor, Indinavir, which has a weak effect on the HTLV-1 protease. In this work, by employing the SuMD method, we reconstructed the unbinding pathways of Indinavir from HIV and HTLV-1 proteases to compare and understand the mechanism of the unbinding and to discover the reasons for the lack of inhibitory activity of Indinavir against the HTLV-1 protease. We achieved multiple unbinding events from both HIV and HTLV-1 proteases in which the RMSD values of Indinavir reached over 40 Å. Also, we found that the mobility and fluctuations of the flap region are higher in the HTLV-1 protease, making the drug less stable. We realized that critically positioned aromatic residues such as Trp98/Trp98' and Phe67/Phe67' in the HTLV-1 protease could make strong π-Stacking interactions with Indinavir in the unbinding pathway, which are unfavorable for the stability of Indinavir in the active site. The details found in this study can make a reasonable explanation for the lack of inhibitory activity of this drug against HTLV-1 protease. We believe the details discovered in this work can help design more effective and selective inhibitors for the HTLV-1 protease.

Unraveling the unbinding pathways of SARS-CoV-2 Papain-like proteinase known inhibitors by Supervised Molecular Dynamics simulation

The COVID-19 disease has infected and killed countless people all over the world since its emergence at the end of 2019. No specific therapy for COVID-19 is not currently available, and urgent treatment solutions are needed. Recent studies have found several potential molecular targets, and one of the most critical proteins of the SARS-CoV-2 virus work machine is the Papain-like protease (Plpro). Potential inhibitors are available, and their X-ray crystallographic structures in complex with this enzyme have been determined recently. However, their activities against this enzyme are insufficient and need to be characterized and improved to be of clinical values. Therefore, in this work, by utilizing the Supervised Molecular Dynamics (SuMD) simulation method, we achieved multiple unbinding events of Plpro inhibitors, GRL0617, and its derivates, and captured and understood the details of the unbinding pathway. We found that residues of the BL2 loop, such as Tyr268 and Gln269, play major roles in the unbinding pathways, but the most important contributing factor is the natural movements and behavior of the BL2 loop, which can control the entire process. We believe that the details found in this study can be used to refine and optimize potential inhibitors like GRL0617 and design more efficacious inhibitors as a treatment for the SARS-CoV-2 virus.

CRISPR-dCas9-Based Artificial Transcription Factors to Improve Efficacy of Cancer Treatment With Drug Repurposing: Proposal for Future Research

Due to the high resistance that cancer has shown to conventional therapies, it is difficult to treat this disease, particularly in advanced stages. In recent decades, treatments have been improved, being more specific according to the characteristics of the tumor, becoming more effective, less toxic, and invasive. Cancer can be treated by the combination of surgery, radiation therapy, and/or drug administration, but therapies based on anticancer drugs are the main cancer treatment. Cancer drug development requires long-time preclinical and clinical studies and is not cost-effective. Drug repurposing is an alternative for cancer therapies development since it is faster, safer, easier, cheaper, and repurposed drugs do not have serious side effects. However, cancer is a complex, heterogeneous, and highly dynamic disease with multiple evolving molecular constituents. This tumor heterogeneity causes several resistance mechanisms in cancer therapies, mainly the target mutation. The CRISPR-dCas9-based artificial transcription factors (ATFs) could be used in cancer therapy due to their possibility to manipulate DNA to modify target genes, activate tumor suppressor genes, silence oncogenes, and tumor resistance mechanisms for targeted therapy. In addition, drug repurposing combined with the use of CRISPR-dCas9-based ATFs could be an alternative cancer treatment to reduce cancer mortality. The aim of this review is to describe the potential of the repurposed drugs combined with CRISPR-dCas9-based ATFs to improve the efficacy of cancer treatment, discussing the possible advantages and disadvantages.