Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol. , 27 February 2025

Sec. Experimental Pharmacology and Drug Discovery

Volume 16 - 2025 | https://doi.org/10.3389/fphar.2025.1510319

This article is part of the Research Topic Drug Discovery in Cancer Research: Success Stories and Open Challenges View all 26 articles

Discovery of highly potent and novel LSD1 inhibitors for the treatment of acute myeloid leukemia: structure-based virtual screening, molecular dynamics simulation, and biological evaluation

Ye Hong,&#x;Ye Hong1,2Yuting Wang&#x;Yuting Wang3Ziyi Hao&#x;Ziyi Hao1Xingxia ZhangXingxia Zhang1Yejun SiYejun Si1Guoqiang LinGuoqiang Lin1Shurong ZhangShurong Zhang1Miao-Miao NiuMiao-Miao Niu3Xiaotian Yang
Xiaotian Yang1*Yanming Zhang
Yanming Zhang1*
  • 1Department of Hematology, The Affiliated Huai’an Hospital of Xuzhou Medical University, Huai’an, Jiangsu, China
  • 2Department of Hematology, Binhai Couty People’s Hospital, Yancheng, Jiangsu, China
  • 3Department of Pharmaceutical Analysis, China Pharmaceutical University, Nanjing, Jiangsu, China

Acute myeloid leukemia (AML) is a highly aggressive hematological malignancy with a significant unmet clinical need for new therapeutic agents. Lysine-specific demethylase 1 (LSD1), a key regulator of leukemia stem cell self-renewal, has emerged as a promising epigenetic target for AML treatment. Herein, we employed an innovative multi-step integrated screening protocol, encompassing pharmacophore modeling, docking screening, molecular dynamics simulation, and biological evaluation, to identify novel LSD1 inhibitors. This comprehensive approach led to the discovery of six potent LSD1 inhibitors (we named these inhibitors LTMs 1–6), with LTM-1 exhibiting the most pronounced inhibitory effects on LSD1 (IC50 = 2.11 ± 0.14 nM) and the highest selectivity for LSD1 over LSD2 (>2370-fold). Notably, LTM-1 demonstrated outstanding antitumor activity both in vitro and in vivo. In vitro, LTM-1 showed potent anti-proliferative effects against LSD1-addicted MV-4-11 leukemia cells (IC50 = 0.16 ± 0.01 μM). In vivo, LTM-1 treatment significantly reduced tumor growth in MV-4-11 xenografted mice. Moreover, LTM-1 did not induce significant changes in liver and kidney function indices, suggesting a favorable safety profile. These results indicate that LTM-1 is a highly promising preclinical candidate for AML treatment, offering a new strategy for the development of more effective and selective LSD1 inhibitors.

1 Introduction

Acute myeloid leukemia (AML) accounts for about 80% of all adult acute leukemia patients and is the most common type of acute leukemia in adults (Dong et al., 2020). AML is a heterogeneous hematological malignancy characterized by abnormal clonal expansion of hematopoietic progenitor cells with differentiation defects (Miles et al., 2020; Rubnitz et al., 2010; Prada-Arismendy et al., 2017). With a mortality rate of up to 60%, AML is one of the deadliest leukemias (Jani et al., 2023). In the past 30 years, despite significant progress in traditional cytotoxic chemotherapy, targeted therapy, and immunotherapy for AML, the number of newly diagnosed leukemia patients worldwide has increased by 46% annually due to aging and the widespread use of chemotherapy leading to secondary leukemia (Wei et al., 2023; Martínez-Cuadrón et al., 2022). Furthermore, the 5-year overall survival rate of patients with AML is less than 50% (Shallis et al., 2019; Tang et al., 2024). With the rapid development of molecular cancer biology, more precise and personalized targeted therapy strategies have shown great potential (Newell and Cook, 2021). However, considerable challenges persist in overcoming treatment resistance and improving the efficacy of AML therapies. There is an urgent and unmet clinical need for the development of novel therapeutic agents that can significantly enhance the prognosis and survival of AML patients.

As the first discovered histone demethylase in 2004, LSD1 is a member of the flavin adenine dinucleotide (FAD)-dependent amine oxidase demethylase family, and another isoform LSD2 was subsequently discovered in 2009 (Shi et al., 2004; Ciccone et al., 2009). LSD1 catalyzes the demethylation of histone 3 lysine 4 methyl 1/2 (H3K4me1/2) and histone 3 lysine 9 methyl 1/2 (H3K9me1/2) (Metzger et al., 2005). Additionally, LSD1 has been found to catalyze the demethylation of some non-histone lysines, including p53, STAT3, E2F transcription factor 1, and DNA methyltransferase 1 (DNMT1), thereby regulating their downstream cellular functions (Xie et al., 2011; Huang et al., 2007; Wang et al., 2008). However, abnormal overexpression of LSD1 has been found in various hematological diseases, including AML, acute lymphocytic leukemia (ALL), myeloproliferative tumors, and chronic myelomonocytic leukemia. Particularly in AML, LSD1 is overexpressed in approximately 60% of cases (Niebel et al., 2014). Moreover, Harri et al. have elucidated that LSD1 is a critical regulator for the self-renewal of leukemia stem cells using the human MLL-AF9 leukemia model (Harris et al., 2012). Furthermore, previous studies have shown that the compound effectively inhibiting LSD1 exhibits synergistic activity with anti-leukemic drugs (Binda et al., 2010). In addition, LSD1 inhibition can enhance the sensitivity of AML cells to all-trans retinoic acid (ATRA) treatment, which was first used to treat AML by inducing cell differentiation (Lokken and Zeleznik-Le, 2012; Schenk et al., 2012). These findings collectively underscore the promise of LSD1 as a highly promising epigenetic target for the treatment of AML (Zhang et al., 2021), providing a compelling biological basis for our investigation.

Given the central role of LSD1 in AML pathogenesis and its association with poor prognosis, targeting LSD1 presents a unique opportunity to develop a new class of therapeutics that could potentially overcome the limitations of current treatments. By inhibiting LSD1, we may be able to disrupt the self-renewal capacity of leukemia stem cells, a key factor in disease relapse, and also modulate the expression of critical genes involved in leukemia progression. This could lead to more effective and durable responses in patients, improving both survival rates and quality of life. Therefore, the motivation for our study stems from the need to develop more effective and selective LSD1 inhibitors. Tranylcypromine (TCP) has been identified as an irreversible and non-selective LSD1 inhibitor by forming a covalent TCP-FAD adduct (Figure 1) (Lee et al., 2006; Fioravanti et al., 2020). Based on the importance of the TCP scaffold in inhibiting LSD1 activity and the clinical demand for developing more effective and selective LSD1 inhibitors, various TCP derivatives have been synthesized (Ganesan et al., 2017; Schulz-Fincke et al., 2018; Koda et al., 2022). Currently, irreversible LSD1 inhibitors based on the TCP scaffold, including TCP, GSK2879552, and ORY-1001 are undergoing clinical trials for AML (Zheng et al., 2016; Huang et al., 2017; Somervaille et al., 2016). In addition, reversible LSD1 inhibitors such as CC-90011 and SP-2577 have also been evaluated in clinical trials (Kanouni et al., 2020; Kurmasheva et al., 2021). However, Sacilotto et al. conducted a comprehensive evaluation of the in vitro inhibitory potential of LSD1 inhibitors at clinical stages and found that SP-2577 had minimal effects on AML markers and did not exhibit satisfactory targeting of LSD1, thereby raising doubts about the efficacy of reversible LSD1 inhibitors in the treatment of AML (Sacilotto et al., 2021; Baby et al., 2023). In summary, there is still an urgent need for the development of highly selective LSD1 inhibitors inhibitors that can potently inhibit cancer cell proliferation without side effects for AML treatment.

Figure 1
www.frontiersin.org

Figure 1. Structures of published LSD1 inhibitors.

The novelty of our approach lies in the application of a multi-step integrated screening protocol to identify novel LSD1 inhibitors. At present, structure-based virtual screening is widely utilized for identifying lead compounds and optimizing drug candidates, attributed to its high efficiency and cost-effectiveness (Aziz et al., 2022; Ballante et al., 2021). Pharmacophore-based screening, utilizing key pharmacophore elements, rapidly screens large databases to identify potential active molecules (Lanka et al., 2023). Molecular docking-based screening analysis evaluates the binding modes and affinities of target proteins with ligands for identification of candidate drugs (Raval and Ganatra, 2022). In prior studies conducted by our group, we successfully employed an integrated screening strategy combining pharmacophore modeling and molecular docking to discover effective inhibitors, including dual KRASG12D single targeted inhibitors and tubulin/NRP1 dual targeted inhibitors (Wang et al., 2022; Zheng et al., 2023). In this study, a multi-step integrated screening protocol encompassing pharmacophore modeling, docking screening, molecular dynamics simulation, and biological evaluation led to the identification of six novel LSD1 inhibitors for the treatment of myeloid leukemia (we named these inhibitors LTMs 1–6). Notably, LTM-1 and LTM-3 exhibited strong in vitro enzymatic inhibitory activity. In particular, LTM-1demonstrated outstanding antitumor activity both in vitro and in vivo. Overall, these results indicate that LTM-1 may be a promising preclinical candidate for further investigation.

2 Materials and methods

2.1 General

The compounds LTMs 1–6, identified in our study, and the positive control GSK2879552 were purchased from WuXi AppTec (Shanghai, China). These were utilized in both in vitro and in vivo assays. Iscove’s Modified Dulbecco’s Medium (IMDM), Roswell Park Memorial Institute (RPMI)-1640 medium, fetal bovine serum (FBS) and penicillin/streptomycin were purchased from Gibco (Grand Island, New Nork, United States). Recombinant human LSD1 (172–852 aa) and full-length LSD2 were obtained from Sino Biological Inc. and ActiveMotif.

2.2 Cell lines and culture conditions

The myeloid leukemia cell lines MV-4-11, MOLT-4, MOLM-16, HL-60, and the human B lymphocyte cell line RPMI-226 were purchased from the American Type Culture Collection (ATCC, Manassas, VA, United States). The myeloid leukemia cell line HAL-01 was obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH (DSMZ, Braunschweig, Germany). All cells were cultured in IMDM or RPMI -1640 medium supplemented with 10% FBS, 100 U/mL penicillin, and 100 μg/mL streptomycin. The cells were maintained in a humidified incubator with an atmosphere of 5% CO2 at 37°C.

2.3 Establishment of a structure-based pharmacophore model for LSD1

The X-ray crystal structure of LSD1 complexed with the inhibitor CC-90011 (PDB ID: 6W4K) was obtained from the Protein Data Bank (PDB) and loaded into the molecular operating environment (MOE, Chemical Computing Group Inc, Montreal, Quebec, Canada). 6W4K was selected due to its high-resolution structure and relevance in previous studies related to LSD1 inhibitors (Kanouni et al., 2020). Within MOE, the QuickPrep tool was employed to execute a series of preparatory steps essential for refining the crystal structure. These steps included the elimination of unbound water molecules, the addition of polar hydrogen atoms, the computation of partial charges, and energy minimization. Subsequently, the Ligand Interaction tool in MOE was engaged to dissect the key key interactions between the LSD1 protein and its bound ligand. Finally, leveraging the detailed understanding of the LSD1 crystal structure, pharmacophore features were meticulously constructed using Pharmacophore Query Editor. These features encapsulate the critical interaction points necessary for the binding of inhibitors to LSD1, including hydrogen-bond donors (Don), hydrophobic features (Hyd), and aromatic centers (Aro).

2.4 Virtual screening

As described in the previously reported methods, virtual screening was performed using the MOE software (Liang et al., 2023). Initially, a database of 121,000 two-dimensional compounds assembled through combinatorial chemistry was converted into three-dimensional structures using an energy minimization algorithm. Subsequently, the prepared LSD1 pharmacophore features were used as a template for preliminary pharmacophore-based screening. Finally, the docking score of −8.53 kcal/mol, obtained from the LSD1 inhibitor GSK2879552, which is in the clinical trial phase, was used as a reference value for docking-based screening. The MOE docking scoring function was employed to evaluate the binding free energy of the ligand with the target protein. Compounds with scores lower than the set reasonable reference value were selected to obtain potential LSD1-targeting candidate inhibitors.

2.5 Molecular dynamics simulation

The crystal structure of LSD1 (PDB ID: 6W4K) was retrieved from the PDB. The crystal complex structures of the ligand LTM-1 with the protein LSD1 were constructed in the MOE and used as the initial coordinates for MD simulations. The ligand LTM-1 and the protein LSD1 were individually subjected to MD simulations using GROMACS (version 2021.5). LSD1 was topologically modeled under periodic boundary conditions, employing the AMBER99SB-ILDN force field for accurate representation. The Acpype Server (www.bio2byte.be) was engaged to generate the topology file of LSD1, which was then integrated with the ligand file to create a comprehensive complex system. Subsequently, the ligand file was added to the protein file to form a complex system. This complex system was immersed in a 1.0 nm cubic simulation box, employing the SPC water model to ensure solvation and equilibrium. To ensure system neutrality, Na+ and Cl were introduced. The solvated system underwent energy minimization, applying the steepest descent algorithm for 50,000 steps. A 100 ps NVT simulation was then performed on the system using the V-rescale thermostat to maintain the system temperature at 300 K. This was followed by a 100 ps NPT simulation, where the Parrinello-Rahman barostat was used to maintain a pressure of 1 bar. Finally, a 500 ns molecular dynamics simulation is conducted for each system, with trajectory data recorded at intervals of 0.1 ns. Data processing is performed using GraphPad Prism 10 software (GraphPad Software Inc., San Diego, CA), and stability is analyzed in conjunction with the root mean square deviation (RMSD), root mean square fluctuation (RMSF), binding-induced conformational changes and protein’s radius of gyration (Rg).

2.6 Binding free energy calculation using MM/PBSA method

The Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method was exploited through the GMX - PBSA tool for calculating the free binding energy of the LSD1 - LTM - 1 complex (Al-Khafaji and Tok, 2020). This computational tool included the PCM and SGB models for solvation energy calculations, along with consideration of van der Waals and electrostatic interaction energies to evaluate the overall energy within the complex. Finally, from the calculation results, we obtained the total binding Gibbs free energy (ΔGtotal). ΔGtotal comprises Van der Waals, Electrostatic, Polar solvation, Non - polar solvation, Net gas phase, and Net solvation. Van der Waals is the contribution of attractive or repulsive forces between molecules or within a single molecule. Electrostatic is the contribution of electrostatic interactions between charged particles within molecules. Polar solvation is the polar contribution to solvation free energy. Non - polar solvation is the non - polar contribution to solvation free energy. Net gas phase is the sum of internal energies of the isolated molecules in vacuum. Net solvation is the energy change when molecules interact with a solvent environment.

2.7 Enzyme inhibition assay

The inhibitory effects of the tested compounds on LSD1 and LSD2 enzymes were methodically assessed using the Lance Ultra LSD1 Histone H3-Lysine 4 Demethylase Assay Kit from PerkinElmer, adhering to previously published methodologies (Li et al., 2022). In brief, the reaction was initiated by adding 4 μL of enzyme solution, with concentrations specified at 2 nM for LSD1 or 172 nM for LSD2, 4 μL of substrate solution [2.5 μM Bio-H3K4me2 (124 amino acids)], and 2 μL of tested compound to a Tris buffer (50 mM Tris-HCl, 50 mM NaCl, 0.01% Tween 20, 1 mM DTT, 10 μM FAD, 10% glycerol, pH 9.0) at room temperature and incubated for 1 h. Following the incubation, 5 μL of a detection mixture, comprising an Eu-labeled H3K4me0 antibody and ULight Streptavidin, was added to the reaction. The fluorescence intensity was then measured employing TR-FRET mode with excitation and emission wavelengths set at 320 nm and 665 nm, respectively, utilizing the Envision (PerkinElmer) system. Each experiment was performed in triplicate.

2.8 Kinase selectivity assays

The kinase selectivity analysis of LTM-1 was performed utilizing the SelectScreen Kinase Profiling Service offered by Thermo Fisher Scientific. In a dose-response assay featuring a ten-point dilution series from 0.25 μM to 128 μM, the inhibitory effects of LTM-1 on substrate phosphorylation reactions catalyzed by a panel of kinases were assessed. Consequently, the half-maximal inhibitory concentrations (IC50 values) for each kinase were determined.

2.9 Cell proliferation assay

The in vitro anti-proliferative activity of the compounds was determined using the CellTiter 96® AQueous NonRadioactive Cell Proliferation Assay Kit (Promega, Madison, United States) as previously described (Ji et al., 2017). In brief, 5–10 × 103 cells and 180 μL of medium were first seeded into each well of a 96-well plate. Subsequently, the cells were exposed to either 0.2% Dimethyl sulfoxide (DMSO, Aladdin Reagent, Shanghai, China) or a serial dilution of the compounds, prepared from a 10 mM stock solution in DMSO, with the final concentration of DMSO being 0.2% for 72 h. 20 μL of MTS reagent was added to each well, and the culture was incubated in the dark at 37°C for 4 h. Finally, the optical density (OD) was then recorded at 490 nm and 690 nm using a microplate reader from Molecular Devices. The proliferation inhibition rate was calculated according to the following formula: inhibition ratio = (OD of DMSO treated wells - OD of compound treated wells)/(OD of DMSO treated wells - OD of blank wells) × 100%. The IC50 values of the compounds inhibiting cell proliferation were calculated using GraphPad Prism 10 (GraphPad Software Inc., San Diego, CA).

2.10 In vivo experiments

Given the potent inhibitory effect of LTM-1 on MV-4-11 cell proliferation in vitro, a previously reported xenograft animal model was established to measure the anti-tumor growth and development effects of LTM-1 (Zhou et al., 2022). The animal experimental design of this study was reviewed and confirmed by the Ethics Committee of China Pharmaceutical University. The experimental procedures complied with the ARRIVE guidelines (https://arriveguidelines.org/). Female BALB/c nude mice aged 4–6 weeks were purchased from Changzhou Cave Animal Co., Ltd. (Changzhou, China). Mice were housed in a specific pathogen-free (SPF) environment, where the temperature was precisely regulated to 25°C ± 2°C, the relative humidity was kept at 50% ± 5%, and a strict 12-h light-dark cycle was maintained. A 200 μL suspension of MV-4-11 cells (1 × 106 cells) was injected subcutaneously into the mice. When the tumor volume was observed to reach 80–100 mm3, the mice were randomly divided into two groups: the vehicle group and the LTM-1 treatment group (10 mg/kg). The mice were administered drugs via intraperitoneal injection daily, with a total treatment period of 15 days. Mouse weight and tumor volume were measured every 3 days. The tumor volume was calculated using the formula (c × c × d)/2 (where c represents the smallest diameter, and d represents the largest diameter). When the animals could no longer access food and/or the average tumor diameter exceeded 15 mm, the mice were anesthetized with 2.5% isoflurane and euthanized by cervical dislocation to minimize suffering. In addition, the liver and kidney index levels were measured using an automatic biochemical analyzer to evaluate the toxic effects of LTM-1 on mice.

2.11 Statistical analysis

The entire statistical analysis was performed with GraphPad Prism 10. The t-tests were employed to evaluate whether there was a significant difference between the two sets of data. P values less than 0.05 were considered as significant. The data are presented as the mean ± standard deviation (SD), n = 3.

3 Results

3.1 Establishment of pharmacophore model

Based on the crystal structure of LSD1 (PDB ID: 6W4K), a pharmacophore model was constructed using MOE to provide clues for identifying novel LSD1 inhibitors. Firstly, the LSD1 protein was preprocessed by removing water molecules, adding polar hydrogens, and calculating partial charges. Subsequently, the binding details of the LSD1 protein and its ligand were analyzed using the ligand interaction tool of MOE to yield key information for constructing the pharmacophore model. As depicted in Figure 2A, the ligand formed two hydrogen bonds with Asp555. Moreover, the ligand engaged in hydrophobic interactions with key amino acid residues Val333, Ile356, Phe538, Ala539, Tyr761, and Ala809 of LSD1. Based on the structure-activity relationship analysis of the LSD1 protein-ligand complex, three representative pharmacophore features were constructed using the pharmacophore query editor: a hydrogen-bond donor feature (F1) corresponding to the Asp555 residue; a hydrophobic interaction feature (F2) corresponding to the Phe538, Ala539, and Tyr761 residues; and an aromatic feature (F3) corresponding to the Trp695 residue. Additionally, the spatial constraints of the pharmacophore model within the LSD1 active site demonstrated that the established pharmacophore features were properly embedded into the depressions of the LSD1 binding pocket (Figure 2B). The constructed pharmacophore model offered critical chemical features for the subsequent virtual screening of novel LSD1 inhibitors.

Figure 2
www.frontiersin.org

Figure 2. Pharmacophore model based on LSD1. (A) Details of pharmacophore model features. Hydrogen-bond donor (F1) is represented by magenta sphere, hydrophobic feature (F2) is depicted by yellow sphere, and aromatic center (F3) is illustrated by brown sphere. Each critical amino acid is represented by a pink stick model and annotated with a three-letter amino acid code. Hydrogen bonds are depicted with black dashed lines. (B) Spatial constraints of the pharmacophore model within the LSD1 active site.

3.2 Virtual screening

As illustrated in Figure 3, this study employed a comprehensive screening strategy that integrated pharmacophore-based screening, molecular docking-based screening, molecular dynamics (MD) simulation, and biological evaluation to identify novel and potent LSD1 inhibitors. Initially, an energy optimization algorithm was applied to convert 121,000 compounds from a combinatorial chemistry-derived database from two-dimensional (2D) to three-dimensional (3D) structures. Subsequently, the LSD1 pharmacophore model was utilized as a probe for pharmacophore-based docking. The 207 compounds obtained through preliminary screening were docked to the LSD1 binding site for further docking-based screening and the binding free energy of each compound was calculated. The lower the binding free energy, the stronger the binding affinity between the compound and the target protein. Based on the docking score of the positive control GSK2879552 (−8.53 kcal/mol), which served as the threshold, the top six compounds (LTMs 1–6) were ultimately identified (Figure 4). The structures of these compounds are displayed in Figure 5, and according to SciFinder and Reaxys databases, LTMs 1–6 have not yet been reported internationally. Subsequently, these hit compounds were prepared for subsequent in vitro enzyme inhibition assays.

Figure 3
www.frontiersin.org

Figure 3. A flowchart of the identification of LSD1 inhibitors via an integrated strategy of in silico screening and biological evaluation.

Figure 4
www.frontiersin.org

Figure 4. The docking scores of LTMs 1–6 and GSK2879552. aBinding free energy between the compounds and the targets (lower binding free energies show stronger binding affinities). The data are presented as the mean ± SD, n = 3. ***p < 0.001, as compared to GSK2879552.

Figure 5
www.frontiersin.org

Figure 5. The chemical structures of LTMs 1–6.

3.3 In vitro LSD1 inhibitory activity

The enzyme inhibitory activity of the six candidate inhibitors obtained was further determined. The LSD1 inhibitor GSK2879552, which has entered the clinical trial stage, was used as a positive control. As shown in Table 1, LTMs 1–6 exhibited nanomolar inhibition of LSD1 (IC50 = 2.11–16.54 nM), and LTMs 1–6 had stronger enzyme inhibition of LSD1 than GSK2879552 (IC50 = 24.53 ± 2.26 nM). Particularly, the enzyme inhibitory activities of LTM-1 (IC50 = 2.11 ± 0.14 nM) and LTM-3 (IC50 = 4.57 ± 0.32 nM) were the most prominent. Moreover, the enzyme inhibitory capacities of LTMs 1–6 against LSD2 were evaluated to determine the selectivity of their inhibition. The selectivity of LTMs 1–6 for LSD1 over LSD2 was superior to that of GSK2879552, with LTM-1 (>2370-fold) and LTM-3 (668-fold) showing the most significant selectivity. In addition, the excellent enzyme inhibitory activity and selectivity of LTM-1 and LTM-3 are consistent with the lowest docking scores of LTM-1and LTM-3 with LSD1 mentioned above. Additionally, the compound with the most potent inhibitory activity and selectivity against LSD1, LTM-1, underwent comprehensive kinase panel analysis to assess its broad-spectrum inhibition potential against 60 different kinases. The data presented in Supplementary Table S1 demonstrate that the IC50 values of LTM-1 against each of the 60 evaluated kinases were all greater than 10 μM. Overall, these results indicate the significant selectivity of LTM-1 for LSD1, reducing the risk of off-target effects and associated toxicity.

Table 1
www.frontiersin.org

Table 1. The inhibitory activities of LTMs 1–6 against LSD1 and LSD2.

3.4 Interaction analysis

Given that LTM-1 and LTM-3 were identified as having the most favorable docking scores for binding to LSD1 and exhibited the highest potency in in vitro enzyme inhibition assays among the six candidates evaluated. Consequently, we conducted an analysis to elucidate the potential binding modes of LTM-1 and LTM-3 with LSD1. As illustrated in Figures 6A, C, LTM-1 and LTM-3 engaged in hydrogen bonding with the critical amino acid residues Asp555, Pro808, and His564 of LSD1. Additionally, they both formed hydrophobic interactions with Val333, Phe538, Ala539, Tyr761, and Trp695 of LSD1. It is noteworthy that LTM-1 established an additional hydrogen bond with Gln358, potentially contributing to its superior binding affinity and inhibitory efficacy. Additionally, the surface topography of the LSD1’s active site, as depicted in Figures 6B, D, revealed a pronounced geometric congruence with the molecular contours of LTM-1 and LTM-3, indicative of an optimal shape complementarity within the binding pocket.

Figure 6
www.frontiersin.org

Figure 6. The predicted docking poses and binding surface of LTM-1 and LTM-3 in the active site of LSD1, respectively. (A) and (B) represent LTM-1; (C, D) represent LTM-3. The structure of LTM-1 is represented by cyan sticks. The structure of LTM-3 is depicted by yellow sticks. Key amino acids in the active site are represented by pink sticks and annotated with three-letter amino acid codes. The surface of the LSD1 protein is rendered with hydrogen bond regions (purple), hydrophobic regions (green), and moderately polar regions (blue).

3.5 MD simulation

To determine the binding affinity and persistence of LTM-1 with LSD1, comprehensive molecular dynamics (MD) simulations were conducted over a 500 ns timescale using GROMACS software (version 2021.5). Trajectory analyses indicated that for the LSD1-LTM-1 complex, the root mean square deviation (RMSD) displayed a moderate initial rise, subsequently stabilizing within the range of 0.14 nm–0.20 nm after approximately 7 ns (Figure 7A). This observation suggests that LTM-1 exhibits stable binding with LSD1. Additionally, Figure 7B presents the root mean square fluctuation (RMSF) of LSD1 residues over the 500 ns simulation period. The RMSF values for key residues, including Val333, Phe538, Ala539, Asp555, His564, Trp695, Tyr761, and Pro808, were all below 0.15 nm, which play a crucial role in binding, exhibit minimal fluctuation and maintain stable binding. Furthermore, the radius of gyration (Rg) values for the LSD1-LTM-1 complex exhibited fluctuations of less than 0.1 nm, suggesting that the protein maintained its structural compactness throughout the simulation (Figure 7C). Based on the obtained structural analysis and RMSD values, the conformational stability of the LSD1-LTM1 complex was also analyzed. As shown in Figure 7D, compared with the initial pose in the LSD1-LTM-1 complex, the superposition of the initial pose and the binding-induced pose did not show any conformational changes, and the RMSD value of the LSD1-LTM-1 complex before and after dynamic simulation was less than 0.5 Å, indicating that the LSD1-LTM-1 complex is in a stable conformation. Collectively, these findings substantiate the stable binding of LTM-1 to the LSD1 active site throughout the entire simulation timeframe.

Figure 7
www.frontiersin.org

Figure 7. MD simulation of LSD1 in complex with LTM-1. (A) RMSD of LSD1-LTM-1 complex. (B) RMSF of LSD1 residues in LSD1-LTM-1 complex. (C) Rg of LSD1 in the complex of LSD1-LTM-1. (D) The conformational changes induced upon binding of the LTM-1 (represented in cyan) to the LSD1 protein (represented in green and gray). The data are presented as the mean ± SD, n = 3. Error bars representing the standard deviation are indicated by red dashed lines.

3.6 Binding free energy calculation using MM/PBSA method

In this study, the MM/GBSA method was employed to evaluate the strength and stability of interactions between the LSD-1 protein and the LTM-1, performing free energy calculations. Detailed energy components are provided in Table 2. The results showed that the LSD-1-LTM-1 complex had a ΔGtotal energy of −8,995.02 ± 210.65 kcal/mol, and the electrostatic energy of the LSD-1-LTM-1 complex was −29,188.78 ± 146.72 kcal/mol. The lower the electrostatic energy of the LSD-1-LTM-1 complex, the stronger its binding stability. Further analysis revealed that the electrostatic energy plays a crucial role in providing maximum stability during the binding process, as its contribution to the total binding free energy is significantly greater than the other components. This confirms its importance in the affinity of LTM-1 with the target protein LSD-1.

Table 2
www.frontiersin.org

Table 2. MM/GBSA energy components of the LSD-1-LTM-1 complex expressed in kcal/mol.

3.7 In vitro cancer cell assays

To further evaluate the in vitro anti-proliferative effects of LTMs 1–6, we conducted an experiment utilizing the MTS assay to assess the cytotoxicity of these compounds against LSD1-addicted MV-4-11 leukemia cells and LSD1-non-sensitive RPMI-8226 cells. The results, as presented in Table 3, indicated that LTMs 1–6 demonstrated significant anti-proliferative activity against the LSD1-addicted MV-4-11 cells compared to GSK2879552 (IC50 = 1.17 ± 0.28 μM), with LTM-1 showing particularly potent effects (IC50 = 0.16 ± 0.01 μM). It was noteworthy that the impact of LTMs 1–6 on the LSD1-non-sensitive RPMI-8226 cells was markedly weaker (IC50 > 100 μM). Taken together, the potent inhibition of LSD1 by LTM-1 shows its potential as a therapeutic candidate for AML. Additionally, we assessed the antiproliferative activity of LTM-1 against other AML cell lines, including MOLT-4, MOLM-16, HAL-01, and HL-60, to confirm its broader applicability. As shown in Supplementary Table S2, LTM-1 demonstrated significant antiproliferative effects across these AML cell lines (IC50 = 0.19–0.39 μM). These findings indicate that LTM-1 not only exerts potent antiproliferative activity against the MV-4-11 cell line but also shows a broad spectrum of anticancer activity against other AML cell lines.

Table 3
www.frontiersin.org

Table 3. Antiproliferative activities of LTMs 1–6 against human MV-4-11 and RPMI-8226 cells by MTS assay.

3.8 Antitumor activity of LTM-1 in vivo

Considering the outstanding inhibitory effect of LTM-1 on LSD1 protein and MV-4-11 cell proliferation, the MV-4-11 xenograft model was used to further evaluate the in vivo anti-tumor activity of LTM-1. Tumor-bearing mice were randomly treated with the vehicle group and 10 mg/kg of LTM-1. As depicted in Figure 8A; Supplementary Figure S1, the tumor growth in mice was significantly reduced upon treatment with LTM-1, compared to the control group. Furthermore, throughout the experimental process, the nude mice in two groups were slowly gaining weight, with no significant differences (Figure 8B), suggesting that treatment with LTM-1 did not induce severe systemic side effects.

Figure 8
www.frontiersin.org

Figure 8. Evaluation of in vivo anti-tumor activity of LTM-1. (A) Changes in tumor volume. (B) Body weight of mice. Results are expressed as means ± SD (n = 3). ***p < 0.001, as compared to the vehicle group.

We conducted an additional assessment of the potential toxicity of LTM-1. The levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), blood urea nitrogen (BUN), and creatinine (CRE) in mouse serum were measured to evaluate the functional status of the liver and kidneys, respectively. As illustrated in Supplementary Figure S2, compared to the vehicle control group, there were no significant changes in the levels of ALT, AST, BUN and CRE among mice treated with various doses of LTM-1. This indicates that LTM-1 does not exert notable toxic effects on the liver or kidneys in mice. In summary, while effectively suppressing tumor growth, LTM-1 demonstrates no apparent toxic impact on major organs such as the liver and kidneys in mice. These findings suggest that LTM-1 possesses a favorable safety profile and potential clinical application value.

4 Discussion

Acute myeloid leukemia (AML) is a formidable heterogeneous malignant hematologic neoplasm. An increasing body of research has confirmed that lysine-specific demethylase 1 (LSD1) is a key regulatory factor in the self-renewal of leukemic stem cells. Although several LSD1 inhibitors have been reported, they have exhibited unsatisfactory targeting specificity and concerning side effects. Consequently, there is an urgent clinical demand for the development of highly selective and potent LSD1 inhibitors. In this study, an integrated strategy combining pharmacophore modeling, docking screening, MD simulation, and biological validation was employed to identify six LSD1 inhibitor (LTMs 1–6). The results of in vitro enzymatic inhibition assays revealed that all six compounds demonstrated potent nanomolar-level inhibitory activity against LSD1. Moreover, the weak inhibitory capacity of LTMs 1–6 against LSD2 underscored their excellent target selectivity. Particularly, LTM-1 and LTM-3 exhibited the most pronounced inhibitory effects on LSD1 and the highest selectivity for LSD1 over LSD2. Then, subsequent interaction analyses indicated that both LTM-1 and LTM-3 formed hydrogen bond interactions with key residues within the LSD1 active site, including Asp555, Pro808, and His564, and also engaged in hydrophobic interactions with Val333, Phe538, Ala539, Tyr761, and Trp695. Additionally, LTM-1 maintained a stable binding state with LSD1 throughout the entire 500 ns MD simulation process. Furthermore, the results of the cell proliferation activity assessment demonstrated that LTMs 1–6 possess strong cytotoxicity against the LSD1-addicted MV-4-11 cells, while having negligible effects on the LSD1-insensitive RPMI-8226 cells, with LTM-1 exhibiting the most potent ability to inhibit cancer cell proliferation. Notably, in the MV-4-11 xenograft model, LTM-1 displayed excellent in vivo antitumor activity without significant toxic side effect. In comparison to other LSD1 inhibitors in clinical trials, such as SP-2577 and GSK2879552, LTM-1 has shown superior enzyme inhibition and selectivity against LSD1 (Liu et al., 2023). This comprehensive comparison underscores the potential of LTM-1 as a promising lead compound for preclinical research of AML and reinforces the significance of our research in the context of ongoing clinical trials. While LTM-1 shows promising results in inhibiting AML progression, we acknowledge the need for further structural optimization to improve its pharmacokinetics and minimize toxicity. Therefore, we propose a detailed plan for future structural modifications, which involves substituting the trifluoromethyl group in LTM-1 with other halogens like chlorine, bromine, or iodine. This modification aims to alter the compound’s lipophilicity, potentially enhancing its pharmacokinetic properties and reducing toxicity, thereby providing valuable insights for designing more effective and safer LSD1 inhibitors. Additionally, the virtual screening methodologies utilized in this study, comprising pharmacophore modeling, docking screening, and molecular dynamics (MD) simulations, are conducive to the expeditious identification of novel therapeutics. This methodological framework can be adapted to recognize potential inhibitors from other natural databases and FDA marketed drug databases, so that discovering other novel and efficient natural active small molecules and marketed drugs targeting LSD1. Consequently, it constitutes a potent instrument in the nascent stages of drug discovery, facilitating the identification of target compounds and hastening the lead optimization process.

5 Conclusion

In summary, we have designed and identified a series of new potent LSD1 inhibitors (LTMs 1–6). In particular, LTM-1 exhibited the most prominent nanomolar enzyme inhibition ability and strongest inhibition selectivity against LSD1. Our in vitro and in vivo assays have shown that LTM-1 significantly inhibits the progression of AML without notable toxic side effects, positioning it as a highly promising lead compound for preclinical AML research. While these findings are promising, we recognize the necessity for further experimental validation to substantiate the safety and efficacy of LTM-1 before advancing to preclinical development. As we continue to optimize the structure of LTM-1 in our laboratory, we are also looking ahead to further enhance its potency and selectivity through additional structural modifications. We plan to conduct extensive pharmacokinetic studies to fully characterize LTM-1, ensuring its safety and efficacy. Moreover, we are committed to moving our preclinical findings towards clinical applications and will explore the potential of LTM-1 and its analogs as dual-target inhibitors to improve therapeutic outcomes. We believe that the integration of computational modeling and experimental validation, as demonstrated in this study, will streamline the drug discovery process and significantly contribute to the development of new treatments for AML. The structure-based virtual screening, molecular dynamics simulation, and biological evaluation methods outlined here offer a novel strategy for designing and optimizing more effective and selective LSD1 inhibitors, setting the stage for future research directions that will build upon our current findings.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics statement

The animal study was approved by the Ethics Committee of China Pharmaceutical University. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

YH: Writing–review and editing, Investigation, Methodology, Validation, Software. YW: Investigation, Supervision, Validation, Writing–original draft, Writing–review and editing. ZH: Conceptualization, Investigation, Methodology, Visualization, Writing–original draft. XZ: Data curation, Methodology, Writing–review and editing. YS: Data curation, Investigation, Validation, Writing–review and editing. GL: Investigation, Supervision, Validation, Writing–review and editing. SZ: Project administration, Resources, Software, Writing–review and editing. M-MN: Formal Analysis, Supervision, Visualization, Writing–review and editing. XY: Conceptualization, Investigation, Methodology, Writing–original draft. YZ: Data curation, Funding acquisition, Investigation, Methodology, Validation, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by Jiangsu Commission of Health scientific research project (Geriatrics) (LKM2023045), Campus Collaborative Innovation Research Project of Jiangsu Medical Vocational College (202491004), Huai’an Health Commission Research General Project (HAWJ2024017), and the research project of Huai’an Science and Technology Bureau (HAB2023).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1510319/full#supplementary-material

References

Al-Khafaji, K., and Tok, T. T. (2020). Molecular dynamics simulation, free energy landscape and binding free energy computations in exploration the anti-invasive activity of amygdalin against metastasis. Comput. Meth Prog. Bio 195, 105660. doi:10.1016/j.cmpb.2020.105660

CrossRef Full Text | Google Scholar

Aziz, M., Ejaz, S. A., Zargar, S., Akhtar, N., Aborode, A. T., Wani, T. A., et al. (2022). Deep learning and structure-based virtual screening for drug discovery against NEK7: a novel target for the treatment of cancer. Molecules 27, 4098. doi:10.3390/molecules27134098

PubMed Abstract | CrossRef Full Text | Google Scholar

Baby, S., Shinde, S. D., Kulkarni, N., and Sahu, B. (2023). Lysine-specific demethylase 1 (LSD1) inhibitors: peptides as an emerging class of therapeutics. ACS Chem. Biol. 18, 2144–2155. doi:10.1021/acschembio.3c00386

PubMed Abstract | CrossRef Full Text | Google Scholar

Ballante, F., Kooistra, A. J., Kampen, S., de Graaf, C., Carlsson, J., and Schulte, G. (2021). Structure-based virtual screening for ligands of G protein–coupled receptors: what can molecular docking do for you? Pharmacol. Rev. 73, 527–565. doi:10.1124/pharmrev.120.000246

PubMed Abstract | CrossRef Full Text | Google Scholar

Binda, C., Valente, S., Romanenghi, M., Pilotto, S., Cirilli, R., Karytinos, A., et al. (2010). Biochemical, structural, and biological evaluation of tranylcypromine derivatives as inhibitors of histone demethylases LSD1 and LSD2. J. Am. Chem. Soc. 132, 6827–6833. doi:10.1021/ja101557k

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciccone, D. N., Su, H., Hevi, S., Gay, F., Lei, H., Bajko, J., et al. (2009). KDM1B is a histone H3K4 demethylase required to establish maternal genomic imprints. Nature 461, 415–418. doi:10.1038/nature08315

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, Y., Shi, O., Zeng, Q., Lu, X., Wang, W., Li, Y., et al. (2020). Leukemia incidence trends at the global, regional, and national level between 1990 and 2017. Exp. Hematol. and Oncol. 9, 14. doi:10.1186/s40164-020-00170-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fioravanti, R., Romanelli, A., Mautone, N., Di Bello, E., Rovere, A., Corinti, D., et al. (2020). Tranylcypromine-based LSD1 inhibitors: structure-activity relationships, antiproliferative effects in leukemia, and gene target modulation. ChemMedChem 15, 643–658. doi:10.1002/cmdc.201900730

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganesan, A., Zhou, C., Wu, F., Lu, L., Wei, L., Pai, E., et al. (2017). Structure activity relationship and modeling studies of inhibitors of lysine specific demethylase 1. Plos One 12, e0170301. doi:10.1371/journal.pone.0170301

PubMed Abstract | CrossRef Full Text | Google Scholar

Harris, W. J., Huang, X., Lynch, J. T., Spencer, G. J., Hitchin, J. R., Li, Y., et al. (2012). The histone demethylase KDM1A sustains the oncogenic potential of MLL-AF9 leukemia stem cells. Cancer Cell 21, 473–487. doi:10.1016/j.ccr.2012.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Sengupta, R., Espejo, A. B., Lee, M. G., Dorsey, J. A., Richter, M., et al. (2007). p53 is regulated by the lysine demethylase LSD1. Nature 449, 105–108. doi:10.1038/nature06092

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Chen, C., Geng, J., Han, D., Wang, T., Xie, T., et al. (2017). Targeting KDM1A attenuates Wnt/β-catenin signaling pathway to eliminate sorafenib-resistant stem-like cells in hepatocellular carcinoma. Cancer Lett. 398, 12–21. doi:10.1016/j.canlet.2017.03.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Jani, C. T., Ahmed, A., Singh, H., Mouchati, C., Al Omari, O., Bhatt, P. S., et al. (2023). Burden of AML, 1990-2019: estimates from the global burden of disease study. JCO Glob. Oncol. 9, e2300229. doi:10.1200/GO.23.00229

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, Y.-Y., Lin, S.-D., Wang, Y.-J., Su, M.-B., Zhang, W., Gunosewoyo, H., et al. (2017). Tying up tranylcypromine: novel selective histone lysine specific demethylase 1 (LSD1) inhibitors. Eur. J. Med. Chem. 141, 101–112. doi:10.1016/j.ejmech.2017.09.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanouni, T., Severin, C., Cho, R. W., Yuen, N. Y. Y., Xu, J., Shi, L., et al. (2020). Discovery of CC-90011: a potent and selective reversible inhibitor of lysine specific demethylase 1 (LSD1). J. Med. Chem. 63, 14522–14529. doi:10.1021/acs.jmedchem.0c00978

PubMed Abstract | CrossRef Full Text | Google Scholar

Koda, Y., Sato, S., Yamamoto, H., Niwa, H., Watanabe, H., Watanabe, C., et al. (2022). Design and synthesis of tranylcypromine-derived LSD1 inhibitors with improved hERG and microsomal stability profiles. ACS Med. Chem. Lett. 13, 848–854. doi:10.1021/acsmedchemlett.2c00120

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurmasheva, R. T., Erickson, S. W., Han, R., Teicher, B. A., Smith, M. A., Roth, M., et al. (2021). In vivo evaluation of the lysine-specific demethylase (KDM1A/LSD1) inhibitor SP-2577 (Seclidemstat) against pediatric sarcoma preclinical models: a report from the Pediatric Preclinical Testing Consortium (PPTC). Pediatr. Blood and Cancer 68, 68. doi:10.1002/pbc.29304

PubMed Abstract | CrossRef Full Text | Google Scholar

Lanka, G., Begum, D., Banerjee, S., Adhikari, N., P, Y., and Ghosh, B. (2023). Pharmacophore-based virtual screening, 3D QSAR, Docking, ADMET, and MD simulation studies: an in silico perspective for the identification of new potential HDAC3 inhibitors. Comput. Biol. Med. 166, 107481. doi:10.1016/j.compbiomed.2023.107481

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. G., Wynder, C., Schmidt, D. M., McCafferty, D. G., and Shiekhattar, R. (2006). Histone H3 lysine 4 demethylation is a target of nonselective antidepressive medications. Chem. and Biol. 13, 563–567. doi:10.1016/j.chembiol.2006.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Su, M., Zhu, W., Kan, W., Ge, T., Xu, G., et al. (2022). Structure–activity relationship study of indolin-5-yl-cyclopropanamine derivatives as selective lysine specific demethylase 1 (LSD1) inhibitors. J. Med. Chem. 65, 4335–4349. doi:10.1021/acs.jmedchem.1c02156

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, S., Geng, Y., Niu, M.-M., Zhang, Y., He, W., Li, J., et al. (2023). Discovery and biological evaluation of novel CARM1/HDAC2 dual-targeting inhibitors with anti-prostate cancer agents. J. Enzyme Inhibition Med. Chem. 38, 2241118. doi:10.1080/14756366.2023.2241118

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H. M., Zhou, Y., Chen, H. X., Wu, J. W., Ji, S. K., Shen, L., et al. (2023). LSD1 in drug discovery: from biological function to clinical application. Med. Res. Rev. 44, 833–866. doi:10.1002/med.22000

PubMed Abstract | CrossRef Full Text | Google Scholar

Lokken, A. A., and Zeleznik-Le, N. J. (2012). Breaking the LSD1/kdm1a addiction: therapeutic targeting of the epigenetic modifier in AML. Cancer Cell 21, 451–453. doi:10.1016/j.ccr.2012.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Cuadrón, D., Megías-Vericat, J. E., Serrano, J., Martínez-Sánchez, P., Rodríguez-Arbolí, E., Gil, C., et al. (2022). Treatment patterns and outcomes of 2310 patients with secondary acute myeloid leukemia: a PETHEMA registry study. Blood Adv. 6, 1278–1295. doi:10.1182/bloodadvances.2021005335

PubMed Abstract | CrossRef Full Text | Google Scholar

Metzger, E., Wissmann, M., Yin, N., Müller, J. M., Schneider, R., Peters, A. H. F. M., et al. (2005). LSD1 demethylates repressive histone marks to promote androgen-receptor-dependent transcription. Nature 437, 436–439. doi:10.1038/nature04020

PubMed Abstract | CrossRef Full Text | Google Scholar

Miles, L. A., Bowman, R. L., Merlinsky, T. R., Csete, I. S., Ooi, A. T., Durruthy-Durruthy, R., et al. (2020). Single-cell mutation analysis of clonal evolution in myeloid malignancies. Nature 587, 477–482. doi:10.1038/s41586-020-2864-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Newell, L. F., and Cook, R. J. (2021). Advances in acute myeloid leukemia. Bmj 375, n2026. doi:10.1136/bmj.n2026

PubMed Abstract | CrossRef Full Text | Google Scholar

Niebel, D., Kirfel, J., Janzen, V., Höller, T., Majores, M., and Gütgemann, I. (2014). Lysine-specific demethylase 1 (LSD1) in hematopoietic and lymphoid neoplasms. Blood 124, 151–152. doi:10.1182/blood-2014-04-569525

PubMed Abstract | CrossRef Full Text | Google Scholar

Prada-Arismendy, J., Arroyave, J. C., and Röthlisberger, S. (2017). Molecular biomarkers in acute myeloid leukemia. Blood Rev. 31, 63–76. doi:10.1016/j.blre.2016.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Raval, K., and Ganatra, T. (2022). Basics, types and applications of molecular docking: a review. IP Int. J. Compr. Adv. Pharmacol. 7, 12–16. doi:10.18231/j.ijcaap.2022.003

CrossRef Full Text | Google Scholar

Rubnitz, J. E., Gibson, B., and Smith, F. O. (2010). Acute myeloid leukemia. Hematology/Oncology Clin. N. Am. 24, 35–63. doi:10.1016/j.hoc.2009.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sacilotto, N., Dessanti, P., Lufino, M. M. P., Ortega, A., Rodríguez-Gimeno, A., Salas, J., et al. (2021). Comprehensive in vitro characterization of the LSD1 small molecule inhibitor class in oncology. ACS Pharmacol. and Transl. Sci. 4, 1818–1834. doi:10.1021/acsptsci.1c00223

PubMed Abstract | CrossRef Full Text | Google Scholar

Schenk, T., Chen, W. C., Göllner, S., Howell, L., Jin, L., Hebestreit, K., et al. (2012). Inhibition of the LSD1 (KDM1A) demethylase reactivates the all-trans-retinoic acid differentiation pathway in acute myeloid leukemia. Nat. Med. 18, 605–611. doi:10.1038/nm.2661

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulz-Fincke, J., Hau, M., Barth, J., Robaa, D., Willmann, D., Kürner, A., et al. (2018). Structure-activity studies on N -Substituted tranylcypromine derivatives lead to selective inhibitors of lysine specific demethylase 1 (LSD1) and potent inducers of leukemic cell differentiation. Eur. J. Med. Chem. 144, 52–67. doi:10.1016/j.ejmech.2017.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shallis, R. M., Wang, R., Davidoff, A., Ma, X., and Zeidan, A. M. (2019). Epidemiology of acute myeloid leukemia: recent progress and enduring challenges. Blood Rev. 36, 70–87. doi:10.1016/j.blre.2019.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Lan, F., Matson, C., Mulligan, P., Whetstine, J. R., Cole, P. A., et al. (2004). Histone demethylation mediated by the nuclear amine oxidase homolog LSD1. Cell 119, 941–953. doi:10.1016/j.cell.2004.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Somervaille, T., Salamero, O., Montesinos, P., Willekens, C., Perez Simon, J. A., Pigneux, A., et al. (2016). Safety, phamacokinetics (PK), pharmacodynamics (PD) and preliminary activity in acute leukemia of ory-1001, a first-in-class inhibitor of lysine-specific histone demethylase 1A (LSD1/KDM1A): initial results from a first-in-human phase 1 study. Blood 128, 4060. doi:10.1182/blood.v128.22.4060.4060

CrossRef Full Text | Google Scholar

Tang, Y., Zheng, Y., Hu, X., Zhao, H., and Cui, S. (2024). Discovery of potent and selective PI3Kδ inhibitors for the treatment of acute myeloid leukemia. J. Med. Chem. 67, 6638–6657. doi:10.1021/acs.jmedchem.4c00094

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Hevi, S., Kurash, J. K., Lei, H., Gay, F., Bajko, J., et al. (2008). The lysine demethylase LSD1 (KDM1) is required for maintenance of global DNA methylation. Nat. Genet. 41, 125–129. doi:10.1038/ng.268

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhang, H., Li, J., Niu, M.-M., Zhou, Y., and Qu, Y. (2022). Discovery of potent and noncovalent KRASG12D inhibitors: structure-based virtual screening and biological evaluation. Front. Pharmacol. 13, 1094887. doi:10.3389/fphar.2022.1094887

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Q., Yao, K., Yang, J., Zhou, Q., Liu, P., Chen, J., et al. (2023). Structure-based drug design of novel triaminotriazine derivatives as orally bioavailable IDH2R140Q inhibitors with high selectivity and reduced hERG inhibitory activity for the treatment of acute myeloid leukemia. J. Med. Chem. 66, 12894–12910. doi:10.1021/acs.jmedchem.3c00835

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Q., Bai, Y., Wu, J., Sun, Y., Wang, Y., Zhang, Y., et al. (2011). Methylation-mediated regulation of E2F1 in DNA damage-induced cell death. J. Recept. Signal Transduct. 31, 139–146. doi:10.3109/10799893.2011.552914

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Liu, M., Yao, Y., Yu, B., and Liu, H. (2021). Targeting LSD1 for acute myeloid leukemia (AML) treatment. Pharmacol. Res., 164. doi:10.1016/j.phrs.2020.105335

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, L., Zou, Y., Xie, T., Wu, X., Tan, Y., Mei, S., et al. (2023). Discovery of a dual tubulin and neuropilin-1 (NRP1) inhibitor with potent in vivo anti-tumor activity via pharmacophore-based docking screening, structure optimization, and biological evaluation. J. Med. Chem. 66, 16187–16200. doi:10.1021/acs.jmedchem.3c01572

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y.-C., Yu, B., Chen, Z.-S., Liu, Y., and Liu, H.-M. (2016). TCPs: privileged scaffolds for identifying potent LSD1 inhibitors for cancer therapy. Epigenomics 8, 651–666. doi:10.2217/epi-2015-0002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Zou, Y., Yang, M., Mei, S., Liu, X., Han, H., et al. (2022). Highly potent, selective, biostable, and cell-permeable cyclic d-peptide for dual-targeting therapy of lung cancer. J. Am. Chem. Soc. 144, 7117–7128. doi:10.1021/jacs.1c12075

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: acute myeloid leukemia, LSD1 inhibitors, structure-based virtual screening, molecular dynamics simulation, biological evaluation

Citation: Hong Y, Wang Y, Hao Z, Zhang X, Si Y, Lin G, Zhang S, Niu M-M, Yang X and Zhang Y (2025) Discovery of highly potent and novel LSD1 inhibitors for the treatment of acute myeloid leukemia: structure-based virtual screening, molecular dynamics simulation, and biological evaluation. Front. Pharmacol. 16:1510319. doi: 10.3389/fphar.2025.1510319

Received: 12 October 2024; Accepted: 03 February 2025;
Published: 27 February 2025.

Edited by:

Viviana di Giacomo, University of Studies G. d’Annunzio Chieti and Pescara, Italy

Reviewed by:

Sirajudheen Anwar, University of Hail, Saudi Arabia
Yuwei Wang, Shaanxi University of Chinese Medicine, China

Copyright © 2025 Hong, Wang, Hao, Zhang, Si, Lin, Zhang, Niu, Yang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaotian Yang, MTU2OTUxMTA0OTJAMTYzLmNvbQ==; Yanming Zhang, emhhbmd5YW5taW5nMjAwNUAxMjYuY29t

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

94% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more