Yanting Zhou1†
Xiandeng Li2†Peifang Luo3Huiting Chen2Yan Zhou4Xueting Zheng4Yuan Yin4
Haoche Wei5
Hongji Liu6
Wen Xia1*
Mingsong Shi4*Xiaoan Li4- 1Key Laboratory of Basic Pharmacology of Ministry of Education and Joint International Research Laboratory of Ethnocentric of Ministry of Education, Zunyi Medical University, Zunyi, Guizhou, China
- 2College of Pharmacy, Chongqing Medical University, Chongqing, China
- 3Department of Cardiovascular Surgery, Affiliated Hospital of Zunyi Medical University, Zunyi, Guizhou, China
- 4NHC Key Laboratory of Nuclear Technology Medical Transformation, Mianyang Central Hospital, School of Medicine, University of Electronic Science and Technology of China, Mianyang, Sichuan, China
- 5State Key Laboratory of Biotherapy/Collaborative Innovation Center of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, China
- 6Department of Ophthalmology, Mianyang Central Hospital, School of Medicine, University of Electronic Science and Technology of China, Mianyang, Sichuan, China
CDK4/6 plays a crucial role in various cancers and is an effective anticancer drug target. However, the gap between clinical requirements and approved CDK4/6 drugs is unresolved. Thus, there is an urgent need to develop selective and oral CDK4/6 inhibitors, particularly for monotherapy. Here, we studied the interaction between abemaciclib and human CDK6 using molecular dynamics simulations, binding free energy calculations, and energy decomposition. V101 and H100 formed stable hydrogen bonds with the amine-pyrimidine group, and K43 interacted with the imidazole ring via an unstable hydrogen bond. Meanwhile, I19, V27, A41, and L152 interacted with abemaciclib through π-alkyl interactions. Based on the binding model, abemaciclib was divided into four regions. With one region modification, 43 compounds were designed and evaluated using molecular docking. From each region, three favorable groups were selected and combined with each other to obtain 81 compounds. Among them, C2231-A, which was obtained by removing the methylene group from C2231, showed better inhibition than C2231. Kinase profiling revealed that C2231-A showed inhibitory activity similar to that of abemaciclib; additionally, C2231-A inhibited the growth of MDA-MB-231 cells to a greater extent than did abemaciclib. Based on molecular dynamics simulation, C2231-A was identified as a promising candidate compound with considerable inhibitory effects on human breast cancer cell lines.
1 Introduction
Cyclin-dependent kinases (CDKs) belong to the serine/threonine protein kinase family and have a classical protein kinase structure (Figure 1). They are activated by cyclin binding to form cyclin-CDK heterodimers. The cyclin-CDK heterodimer regulates cell cycle progression (Yousuf et al., 2022; Zabihi et al., 2022; Hope et al., 2023) and phosphorylates the upstream signals, including the retinoblastoma (Rb) protein (Macaluso et al., 2006; Choi and Anders, 2014; Spring et al., 2020). The phosphorylated Rb induces the release of the critical transcription factor early 2 factor (E2F) to initiate the transcription of target genes, thereby facilitating the activation of various signaling pathways, including the G1/S transition, mitochondrial dynamics, and metabolism (Sanchez-Martinez et al., 2015; Ingham and Schwartz, 2017; Matson and Cook, 2017; Sanchez-Martinez et al., 2019). When CDK4/6 is inhibited, the G1 phase of proliferating cells is blocked (Yuan et al., 2021; Chen et al., 2022; Watt and Goel, 2022; Yousuf et al., 2022). Additionally, CDK4/6 is disordered in numerous cancers (Wesierska-Gadek et al., 2009; Sanchez-Martinez et al., 2019; Yuan et al., 2021; Ghafouri-Fard et al., 2022; Kargbo, 2022; Yousuf et al., 2022; Zabihi et al., 2022), such as breast cancer (Gelbert et al., 2014; Chen et al., 2022; Watt and Goel, 2022), osteosarcoma (Wang and Bao, 2022), and acute megakaryoblastic leukemia (Qi et al., 2022). Thus, CDK4/6 is an effective anticancer drug target, especially in hormone receptor-positive and human epidermal growth factor receptor 2-negative (HR+/HER2-) breast cancer.
FIGURE 1. Structure of human CDK6 and abemaciclib. (A) Structure of human CDK6 labeled as different domains. (B) Different regions labeled with human CDK6. (C) Structure of abemaciclib, and biochemical activity of abemaciclib with human CDK6.
Several CDK4/6 inhibitors have been recently developed and employed to treat cancer (Supplementary Figure S1) (Bhurta and Bharate, 2022; Shi Z. F. et al., 2022; Ettl et al., 2022; Ji et al., 2022; Jorda et al., 2022; Manohar, 2022; Watt and Goel, 2022; Yousuf et al., 2022; Yu et al., 2022; Zeng et al., 2022). Three adenosine triphosphate (ATP)-competitive CDK4/6 inhibitors have been approved by the U.S. Food and Drug Administration (FDA); these are palbociclib (PD-0332991), ribociclib (LEE011), and abemaciclib (LY2835219) (Dhillon, 2015; Kim, 2017; Syed, 2017). Palbociclib is an oral, reversible, and selective CDK4/6 inhibitor with inhibition activity (Ki) of 0.26 ± 0.03 nM (CDK4) and 0.26 ± 0.07 nM (CDK6) (Chen et al., 2016). The FDA approved palbociclib in combination with letrozole as the first-line treatment for advanced breast cancer (Dhillon, 2015). Meanwhile, the selective CDK4/6 inhibitor ribociclib inhibits the activity, with Ki of 0.53 ± 0.08 nM and 2.3 ± 0.3 nM for CDK4 and CDK6, respectively (Chen et al., 2016). Furthermore, the combination of ribociclib and an aromatase inhibitor was approved as a treatment for advanced breast cancer in the USA in March 2017 (Syed, 2017). Conversely, abemaciclib was approved as a combination with fulvestrant or monotherapy (Kim, 2017). Abemaciclib is the only CDK4/6 inhibitor approved for breast cancer treatment as monotherapy (Abraham, 2018). Nevertheless, these three inhibitors have several side effects, such as neutropenia, gastrointestinal toxicities, anemia, and thrombopenia (Syed, 2017; Braal et al., 2021). The reasons for dose modification were myelosuppression (ribociclib and palbociclib) and diarrhea (abemaciclib). Therefore, developing novel CDK6 inhibitors is necessary. Additionally, there is a gap between the clinical requirements and approved CDK4/6 drugs. Accordingly, there is an urgent need to develop selective and oral CDK4/6 inhibitors, particularly for monotherapy.
N-{5-[(4-ethylpiperazin-1-yl) methyl] pyridin-2-yl}-5-fluoro-4-[4-fluoro-1-isopropyl-2-methyl-1H-benzo (d) imidazol-6-yl] pyrimidin-2-amine (abemaciclib; LY2835219; Verzenio; Figure 1) is a selective, administered, ATP-competitive, and reversible kinase inhibitor of CDK4/6 (Gelbert et al., 2014). Abemaciclib has several distinct chemical [2-anilino-2, 4-pyrimidine-(5-benzimidazole) scaffold], pharmacological (continuous dosing schedule), and clinical properties (single agent therapy) compared to palbociclib and ribociclib (Yuan et al., 2021; Wander et al., 2022). Abemaciclib inhibits Rb phosphorylation and induces G1 cell cycle arrest, resulting in antitumor activity (Gelbert et al., 2014; Wright and Md, 2021; Nardone et al., 2022; Nawa et al., 2022; Sun et al., 2022; Wander et al., 2022; Wang and Bao, 2022). Abemaciclib is a mainstay of HR + breast cancer treatment. It is currently under clinical trials for several other tumors, such as non-small cell lung cancer, Ewing sarcoma, and renal cell carcinoma (Kim, 2017; Kosovec et al., 2017; Small et al., 2017; Dowless et al., 2018; Naz et al., 2018; Wander et al., 2022; Wang and Bao, 2022). Meanwhile, the adverse effects of abemaciclib mainly include neutropenia, leukopenia, anemia, diarrhea, and eosinophilic pneumonia (Xie Q. et al., 2022; Lim et al., 2022; Mitarai et al., 2022). Notably, a black-box warning for abemaciclib resulted from venous thromboembolism (Watson et al., 2022). In vitro, abemaciclib inhibits CDK6 functions with Ki of 2.4 ± 1.2 nM (or 0.07 ± 0.01 nM or 8.2 ± 1.1 nM) and half-maximal inhibitory concentration (IC50) of 9.9 nM (Gelbert et al., 2014; Chen et al., 2016; Torres-Guzman et al., 2017). Abemaciclib potently inhibits osteosarcoma cell lines, with IC50 ranging from 90 nM to >20 μM (Wang and Bao, 2022). Some other potential targets have been identified from kinase screens in vitro, specifically CDK9, proviral insertion in murine malignancies serine/threonine kinases (PIM1), glycogen synthase kinase-3 beta (GSK3β), homeodomain-interacting protein kinase 2 (HIPK2), dual-specificity tyrosine-phosphorylation-regulated kinase 2 (DYRK2), and casein kinase 2 (CK2) with 57 ± 42 nM, 50 nM, 192 nM, 31 nM, 61 nM, and 117 nM, respectively (Gelbert et al., 2014; Kaltheuner et al., 2021). These studies indicate that abemaciclib has remarkable potential for use in developing novel inhibitors with improved inhibition activity and selectivity.
In this study, we describe a binding model and the key interactions between abemaciclib and human CDK6. The binding model of abemaciclib with human CDK6 was optimized using molecular dynamics (MD) simulations. The binding free energy was calculated using the molecular mechanics/generalized Born (GB) surface area (MM/GBSA) method and decomposed to identify the key residues that mainly contribute to the binding of abemaciclib with CDK6. Based on the binding models, some abemaciclib derivatives were designed, and their binding affinities were tested using molecular docking. Candidate inhibitors were synthesized, and their inhibitory activity against CDK6 was evaluated in vitro using MDA-MB-231 and MCF-7 cell lines.
2 Materials and methods
2.1 Molecular dynamics simulations
The abemaciclib/CDK6 complex structure with a resolution of 2.27 Å was obtained from the Protein Database Bank (PDB ID: 5L2S) (Chen et al., 2016). The restrained electrostatic potential (RESP) protocol (Bayly et al., 1993) was employed to generate the atom charge of abemaciclib. The force field parameters were produced based on the general amber force field (GAFF, version 2) (Wang et al., 2004) in AMBERTools21 (Case et al., 2020). The AMBER ff19SB force field (Tian et al., 2020) was applied to construct the topology parameters of human CDK6. Afterward, the abemaciclib/CDK6 system was neutralized by adding sodium chloride ions and solvated in a cuboid box with TIP3P water (Jorgensen et al., 1983). Finally, the solvated abemaciclib/CDK6 complex system included human CDK6 (291 residues), one small molecule (abemaciclib), and solvent water (19,723). The complex was subjected to the steepest descent method (9,000 steps) and conjugate gradient (1,000 steps), while abemaciclib and human CDK6 were fixed. Next, the entire abemaciclib/CDK6 system was optimized using the 10,000-step conjugate gradient method. Subsequently, Langevin dynamics (Martyna et al., 1994; Feller et al., 1995) was performed to increase the temperature of the solvated abemaciclib/CDK6 system from 0 K to 300 K for 200 ps. Subsequently, isotropic position scaling (Berendsen et al., 1984) was employed to maintain the pressure at 1 bar and perform a 200 ps simulation. Then, the abemaciclib/CDK6 system was equilibrated with an isothermal-isobaric ensemble at 300 K and 1 bar, and then a 200 ps simulation was run. Finally, the solvated abemaciclib/CDK6 system was subjected to 500 ns MD simulations to collect and analyze data for abemaciclib binding with human CDK6. In this study, MD simulations were performed using the AMBER20 software (Case et al., 2020). To analyze the trajectories of abemaciclib/CDK6, the CPPTRAJ module (Roe and Cheatham, 2018; Cheatham et al., 2019) was employed for calculating the root mean square deviation (RMSD), root-mean-square fluctuation, hydrogen bonds, distance, and angle. Detailed information regarding MD simulations and cluster analysis is provided in the Supplementary Material.
2.2 Binding free energy calculation
The MM/GBSA approach has been widely used to evaluate ligand and enzyme systems; it was used to calculate the binding free energies between abemaciclib and human CDK6 (Srinivasan et al., 1998; Lee et al., 2004; Tse and Verkhivker, 2015; Wang et al., 2017; Shi and Xu, 2019). The framework of MM/GBSA has been extensively discussed (Honig and Nicholls, 1995; Genheden and Ryde, 2015; Onufriev and Case, 2019). Different energy terms of abemaciclib/CDK6 complex system were obtained via a statistical average from the last 100 ns of the MD trajectory over 1,000 frames. The entropy was averaged over 100 frames. The contribution of each residue was obtained by decomposing the enthalpy energy between human CDK6 and abemaciclib with MM/GBSA energy decomposition (Gaillard and Simonson, 2014). The energy terms were assessed using the MMPBA.py program in AMBERTools21 (Miller et al., 2012). Additional information about the calculation of binding free energy is provided in the Supplementary Material.
2.3 Chemistry
Compounds 1 and 2 underwent N coupling reactions in the presence of potassium carbonate to generate compound 3 (Scheme 1). Intermediate 3 was reduced to 4 in a hydrogen atmosphere under palladium catalysis. Compound 5 was reduced to compound 6 under the action of iron powder and ammonium chloride. Then, subsequent amide condensation of 6 with acetic anhydride afforded compound 7. Intermediate 7 underwent cyclization and hydrolysis to yield compounds 8 and 9 successively. Intermediate 9 and isopropylboronic acid underwent a C-N coupling reaction to generate compound 10. Intermediate 10 then underwent a Suzuki coupling reaction to generate compound 12. Finally, intermediates 12 and 4 underwent Buchwald coupling to yield the target product C2213-A. Detailed chemical information is provided in the Supplementary Material.
SCHEME 1. Synthetic route of compound C2213-A. Reagents and conditions: (A) K2CO3, MeCN, 80°C; (B) Pd/C, H2, rt; (C) Fe/NH4Cl, MeOH/H2O, 80°C; (D) AcOK, acetic anhydride, CHCl3, reflux; (E) isopentyl nitrite, 18-Crown-6, K2CO3, CHCl3, reflux; (F) KOH, H2O, rt; (G) cyclopropylboronic acid, Na2CO3, Cu (AcO)2, 2,2′-Bipyridine, DCM, 80°C; (H) bis (pinacolato) diboron, PdCl2 (dppf), AcOK, dioxane, 90°C; (I) 2,4-Dichloro-5-fluoropyrimidine, PdCl2 (dppf), K2CO3, dioxane/water (V/V, 10/1), 80°C; (J) compound 4, Pd2 (dba)3, Xantphos, Cs2CO3, dioxane, reflux.
2.4 Pharmacological evaluation
The inhibitory activity of C2213-A against CDK6 was validated using the KinaseProfiler radiometric protein kinase assay and IC50 values, which were obtained from Eurofins Pharma Discovery Services UK Limited (Wolverhampton, UK). The ATP concentrations used here represent the Km values of human CDK6. Cell proliferation assays for MDA-MB-231 and MCF-7 cells, colony formation assays for MDA-MB-231 cells, and cell cycle analysis based on MDA-MB-231 cells are described in the Supplementary Material.
3 Results and discussion
3.1 System stability
The stability of the abemaciclib/CDK6 complex system was determined from the RMSD values of CDK6 and abemaciclib. The structure of the receptor protein CDK6 was stable after 100 ns (Figure 2). However, the conformation of abemaciclib fluctuated during the simulation. The average RMSD value of CDK6 (2.57 ± 0.28 Å) was larger than that of abemaciclib (2.03 ± 0.59 Å) for the abemaciclib/CDK6 complex system. Meanwhile, the RMSD value of CDK6 was centered at 2.57 Å and that of abemaciclib was between 1 and 4 Å (Supplementary Figure S3). This result suggested that the fluctuation in the conformation of abemaciclib was greater than that in CDK6, which is more stable than abemaciclib. However, the RMSD value was increased to 3.33 ± 0.43 Å when abemaciclib was not used for human CDK6 (Supplementary Figure S4). This stable CDK6 and fluctuating abemaciclib conformation can also be found in the other two replicated abemaciclib/CDK6 complex systems (Supplementary Figure S5). In addition, the RMSD value of the C-lobe region (residues: 104–300, 3.18 ± 0.41 Å) was larger than that of the N-lobe region (residues: 11–98, 2.21 ± 0.35 Å), as shown in Supplementary Figure S6. The structure of human CDK6 was stable in the simulations after 100 ns, whereas that of abemaciclib was unstable throughout the simulation. The stability of this system suggests that CDK6 with abemaciclib is more stable than apo-CDK6, and abemaciclib conformation fluctuates in the ATP-binding pocket of human CDK6.
FIGURE 2. Stability analysis of the abemaciclib/CDK6 system. (A) Root mean square deviation (RMSD) value of heavy backbone atoms for human CDK6 and that of heavy atoms of abemaciclib along 500 ns molecular dynamics (MD) simulation for the abemaciclib/CDK6 system. (B) Root-mean-square fluctuation variations for Cα atom of human CDK6 for the abemaciclib/CDK6 system from the 500 ns MD simulation.
The root-mean-square fluctuation in the conformation of human CDK6 residues was analyzed based on the frames in the 500 ns simulations. Notably, there were four regions of human CDK6 with significant fluctuations (Figure 2)—the regions at positions 47–55 (loop region between β3 and αC), 84–92 (loop region between β4 and β5), 167–181 (T-loop), and 245–260 (loop region between αG and αH). The T-loop is in an open conformation from the crystal structures of human CDK6, such as fisetin (PDB ID: 1XO2) (Lu et al., 2005), palbociclib (PDB ID: 2EUF) (Lu and Schulze-Gahmen, 2006), aminopurvalanol (PDB ID: 2F2C) (Lu and Schulze-Gahmen, 2006), and pyrido-[4′, 3′: 4, 5] pyrrolo [2,3-d] pyrimidine compound (PDB ID: 4TTH) (Li et al., 2014) binding with human CDK6 (Supplementary Figure S7). Meanwhile, some published crystal structures, such as 4-(pyrazol-4-yl)-pyrimidines (PDB ID:3NUX and 3NUP) (Cho et al., 2010), 7-azabenzimidazoles (PDB ID:4AUA and 4EZ5) (Cho et al., 2012), palbociclib (PDB ID: 5L2I) (Chen et al., 2016), abemaciclib (PDB ID: 5L2S) (Chen et al., 2016), ribociclib (PDB ID: 5L2T) (Chen et al., 2016), pyrazolopyrimidine compounds (PDB ID: 6OQL and 6OQO) (Bronner et al., 2019), and human CDK6, were missing from the T-loop. In addition, the T-loop for protein kinase has diversified conformations, with open, closed, and intermediate conformations. For instance, human FAK may exhibit an open or closed T-loop conformation while binding with FAK inhibitors (Shi et al., 2022a; Guo et al., 2022). However, some inhibitors, such as the pan-kinase inhibitor bosutinib-bound SIK2, can also bind to the protein kinase with conformational plasticity in the active pocket of the T-loop (Shi et al., 2022c). Thus, the large fluctuation in the conformation of the T-loop of CDK6 binding with abemaciclib is a classical phenomenon. Additionally, the T-loop of CDK6 is phosphorylated to activate kinase activity and regulate downstream signaling.
These observations indicate that different conformations of the T-loop are extremely important for studying the interactions between abemaciclib and human CDK6. The loop region (residues 47–55) between β3 and αC was also found to be flexible. This region directly or indirectly regulates the salt bridge between K43 in β3 and E61 in αC; it also affects the kinase activity of CDK6. From the crystal inhibitor/human CDK6 complex structures, the salt bridge can be formed with fisetin, palbociclib, and aminopurvalanol compounds binding to human CDK6 (Supplementary Figure S8). Meanwhile, a salt bridge cannot form when palbociclib, abemaciclib, and ribociclib bind to human CDK6. Thus, the flexibility of this loop region suggested that the conformation of αC fluctuated during the simulation. The RMSD value of αC was calculated after the superposition of the entire human CDK6 structure, indicating that the conformation of αC changed during the simulation time (Supplementary Figure S9). The distance between K43 and E61 was also determined; it was 11.51 ± 1.22 Å and 11.46 ± 1.18 Å for K43_NZ−E61_OE1 and K43_NZ−E61_OE2, respectively (Supplementary Figure S10). This result implied that the salt bridge was not formed in the simulation, consistent with the experimental results (9.88 Å and 10.82 Å for K43_NZ−E61_OE1 and K43_NZ−E61_OE2, respectively) (Chen et al., 2016). The flexible region (residues: 84–92) was located between β4 and β5 in the N-lobe of the kinase domain, far from the ATP-binding site. The last flexible region in the C-lobe was also far from the ATP-binding pocket, which may not directly affect inhibitor binding with CDK6. These flexible regions, especially the loop region (between β3 and αC) and the T-loop, can be induced by the binding of abemaciclib to the ATP-binding site of human CDK6. Furthermore, the overall structure of CDK6 was stable near the T-loop region in the simulation.
Snapshots at 100 ns, 200 ns, 300 ns, 400 ns, and 500 ns, obtained from the MD simulation, were drawn to demonstrate the conformational change in abemaciclib in the ATP-binding site of human CDK6 (Supplementary Figure S11). However, abemaciclib remained as a stable binding complex during the simulation period, which indicated that abemaciclib can directly bind to the ATP-binding pocket of human CDK6. The piperazine group of abemaciclib, which was exposed to the solvent, fluctuated during the simulation. This result also suggests that the piperazine group of abemaciclib can be replaced with other solvent-friendly groups, such as piperidine, morpholine, cyclohexane, and cyclopentane. The overall structure of abemaciclib fluctuated in the ATP-binding site, which was consistent with the fluctuation in the RMSD value of abemaciclib. This surface-exposed piperazine group of abemaciclib can be considered as a modification.
3.2 Hydrogen bond analysis
Hydrogen bonds play an important role as kinase inhibitors bound to protein kinases and are commonly required for the actin of potential kinase inhibitors (Stavenger, 2008; Xing et al., 2015; Pang et al., 2021). Moreover, hydrogen bonds may significantly contribute to the binding of inhibitors with human CDK6 as a protein kinase (Xie Z. L. et al., 2022). Accordingly, the number of hydrogen bonds formed between abemaciclib and human CDK6 was first examined for the abemaciclib/CDK6 complex system from the simulation trajectories (Supplementary Figure S12). The wide range of hydrogen bond numbers (0–7) for the abemaciclib/CDK6 system indicated that some unstable hydrogen bonds were formed during the simulation. The number of hydrogen bonds formed between human CDK6 and abemaciclib was 2.59 ± 0.86. Hydrogen bonds were mainly occupied by 2 or 3 number, indicating that there are disordered hydrogen bonds between abemaciclib and the ATP-binding site of human CDK6.
Hydrogen bond occupancy analysis was performed to determine stable hydrogen bonds between abemaciclib and CDK6. In this study, a hydrogen bond was formed between the acceptor and donor atoms <3.5 Å, and the internal angle of the acceptor···H-donor was >120°. The results of hydrogen bond occupancy analysis are summarized in Figure 3. A high occupancy was observed for abemaciclib with the hinge loop of human CDK6. The hydrogen bond occupancy between atom N of V101 and atom N3 of abemaciclib was 86.40%, indicating that this hydrogen bond can be stably formed during the simulation. The other hydrogen bond formed between abemaciclib and the hinge loop of CDK6 was transformed between the ND1 atom of H100 (54.31%) and the O atom of V101 (32.27%) to form a hydrogen bond with N8 of abemaciclib. However, in the last 200 ns of the simulation, hydrogen bonds were mainly formed with H100 for 80.27% occupancy (Supplementary Figure S13). The formation of these two hydrogen bonds, which formed between abemaciclib and the hinge loop of CDK6, was consistent with the results of experimental and theoretical studies (Yuan et al., 2021; Chen et al., 2022; Sokolsky et al., 2022; Wang et al., 2023). In addition, the hydrogen bond between the atom N4 of abemaciclib and the NZ atom of K43 was not unstable from the occupancy with 22.07% for the 500 ns simulation and 20.74% for the last 200 ns simulation, respectively. Therefore, two stable hydrogen bonds were formed between abemaciclib and the hinge loop of CDK6.
FIGURE 3. Hydrogen analysis between abemaciclib and human CDK6.
Occupancy was expressed as a percentage of the period (500 ns, 50,000 frames) during which specific hydrogen bonds were formed. A hydrogen bond is defined as the distance between the acceptor and donor atoms (<3.5 Å), with an internal angle between the H-acceptor, hydrogen atom, and H-donor (>120°).
The distance between the acceptor and donor atoms for the hydrogen bond and the angles between the acceptor, hydrogen, and donor atoms were also determined (Supplementary Figures S14, S15). The distances for Abe@N3 … V101@N (abemaciclib was labeled as Abe), H100@ND1 … Abe@N8, V101@O … Abe@N8, and Abe@N4 … K43@NZ were 3.12 ± 0.18 Å, 4.01 ± 1.10 Å, 3.69 ± 0.53 Å, and 3.21 ± 0.32 Å from the 500 ns simulation, respectively. Meanwhile, those distance were 3.13 ± 0.18 Å, 3.32 ± 0.35 Å, 3.94 ± 0.42 Å, and 3.19 ± 0.25 Å for Abe@N3 … V101@N, H100@ND1 … Abe@N8, V101@O … Abe@N8, and Abe@N4 … K43@NZ in the last 200 ns simulation, respectively. This result indicates that the hydrogen bond between the N3 of V101 and the amine-pyrimidine group has the same strength as that between ND1 of H100 or the O of V101 and the amine-pyrimidine group. These results were consistent with those of the occupancy analysis. Meanwhile, the distance between ND1 of H100 and the amine-pyrimidine group of abemaciclib (3.32 ± 0.35 Å) was maintained to form a hydrogen bond. They also indicate that this hydrogen bond is stable during the last 200 ns of the simulation. However, these two hydrogen bonds should also be maintained when designing novel CDK6 inhibitors. Thus, other scaffolds, which also formed hydrogen bonds with the hinge loop, such as 1H-pyrrolo [2,3-b]pyridine, 1H-indazole, adenine, 7H-pyrrolo [2,3-d] pyrimidin-4-amine, 1H-pyrazole-3-amine, 7H-pyrrolo [2,3-d] pyrimidine pyrimidine-4, 6-diamine, and pyridine, can be replaced with the 2-amine-pyrimidine group.
3.3 Binding free energy
The information regarding the binding of abemaciclib to CDK6 was obtained from simulations. However, determining the binding affinity of abemaciclib and CDK6 is challenging. Thus, the binding free energy for abemaciclib and CDK6 was calculated using the MM/GBSA method as a quantitative measurement. This study used frames from the last 100 ns to calculate the binding free energy. Detailed information regarding the binding free energies of the abemaciclib/CDK6 system is summarized in Supplementary Table S1. Abemaciclib bound to the ATP-binding pocket of human CDK6, with
TABLE 1. Binding free energies, decomposition, electrostatic interactions (
Van der Waals (vdW) interactions (−51.70 kcal/mol) mainly contributed to the binding between abemaciclib and CDK6. The electronic interaction (−27.31 kcal/mol) is also favorable for the binding process. vdW interactions contribute more energy (−24.39 kcal/mol) than electronic interactions for the binding of abemaciclib with CDK6. Generally, the binding energy is divided into two major categories—polar binding energy (
3.4 Free energy decomposition
The binding of abemaciclib to human CDK6 has been explained by MD simulation. The key residues in the abemaciclib/CDK6 complex system remain unclear. Thus, the per-residue energy decomposition approach was employed to identify the key residues, which contributed more than 0.50 kcal/mol to the abemaciclib/CDK6 system (Supplementary Table S6).
H100 and V101 in the hinge loop of CDK6 contributed −1.17 and −0.91 kcal/mol from the hydrogen bonds formed with the amine-pyrimidine group of abemaciclib, which agreed with hydrogen bond analysis results (Table 2). H100 formed a stable hydrogen bond and contributed more than V101 did. Thus, the hydrogen bonds formed between the hinge loop of CDK6 and the amine-pyrimidine group of abemaciclib must remain for novel inhibitor binding with CDK6.
TABLE 2. Free energy decomposition for the abemaciclib/CDK6 complex on the individual residue basis.
However, the largest contribution residue, K43 (−3.14 kcal/mol), forms an unstable hydrogen bond with benzo [d]imidazole group of abemaciclib. Meanwhile, K43 also formed a hydrogen bond with D163 (the distance between NZ of K43 and OD1 of D163 was 2.90 ± 0.30 Å, Supplementary Figure S17). This hydrogen bond between K43 and D163 reduced the strength of the hydrogen bond between abemaciclib and K43. Thus, the contribution of D163 (3.94 kcal/mol) was unfavorable for the binding of abemaciclib to human CDK6. Additionally, the OD2 of D163 formed a hydrogen bond with the ND2 of N150 (distance between the ND2 of N150 and the OD2 of D163 was 3.58 ± 0.57 Å). This weakened the hydrogen bond between D163 and K43 to contribute to the binding of abemaciclib with K43; N150 contributed −0.65 kcal/mol. The orientation of R168 must be considered to disturb the interactions among abemaciclib, D163, K43, and N150. For example, the distance between the NH2 of R168 and the OD2 of D163 was 14.49 ± 1.83 Å for the first 10 ns and 5.20 ± 0.32 Å for the last 10 ns. In the simulation, the location of D163 remained, and R168 was transformed into the ATP-binding site (Supplementary Figure S18) and contributed −1.32 kcal/mol. Meanwhile, the region (163DFGLAR168) at the beginning of the T-loop plays a key role in the kinase activity of CDK6. The conserved DFG motif also regulates the activity of other protein kinases, such as SIKs, FAK, and ROCKs (Roskoski, 2016; Bronner et al., 2019; Tesch et al., 2021; Arter et al., 2022; Shi et al., 2022a; Shi et al., 2022b; Ladduwahetty et al., 2022). Furthermore, the binding of palbociclib with human CDK6 showed that R168 was far from the ATP-binding site (Supplementary Figure S19) (Lu and Schulze-Gahmen, 2006; Chen et al., 2016). This orientation of R168 means that R168 can be employed to form an interaction with the other groups by substituting the benzo [d]imidazole group of abemaciclib.
The P-loop also plays an important role in the binding of inhibitors to protein kinases, such as TYK2 and ABL1 (Yang et al., 2014; Wang et al., 2017; Fensome et al., 2018). The I19 in P-loop region also contributes the largest energy (−3.27 kcal/mol) for the binding of abemaciclib with human CDK6. In particular, the side chain of I19 contributed −3.59 kcal/mol, and the backbone contributed 0.31 kcal/mol. The side chain of I19 interacted with the π-alkyl interaction of the pyridine ring of abemaciclib, which also resulted from the distance between the pyridine ring of abemaciclib and CD1 of I19 (3.96 ± 0.38 Å for the last 200 ns, Supplementary Figure S20). Experimental and computational studies have also shown that I19 plays a key role in the binding of other inhibitors with CDK6 (Chen et al., 2022; Yuan et al., 2022; Bhattacharya et al., 2023; Wang et al., 2023).
V27, located in β2, contributed approximately −1.95 kcal/mol to abemaciclib binding with CDK6; the side chain of V27 interacted with the imidazole ring of abemaciclib. The distances between the CG1/CG2 of V27 and the imidazole ring of abemaciclib were 4.36 ± 0.50 Å and 4.18 ± 0.34 Å, respectively (Supplementary Figure S21). The distance of the V27_CG1-imidazole ring decreased by approximately 170 ns. However, the distances for the V27-imidazole ring were well maintained in the last 300 ns of the simulation. This result indicated that the π-alkyl interaction for V27-abemaciclib was mainly attributed to the abemaciclib bound to human CDK6. In other words, the imidazole ring must balance the π-alkyl from V27 and the hydrogen bond from K43. Thus, the imidazole ring can be substituted by the 1, 2, 3-triazole ring.
The ATP-binding site is formed between the N-lobe and C-lobe for protein kinases. β3, which includes A41, is the main domain of the N-lobe that forms the ATP-binding pocket of human CDK6. A41 forms a π-alkyl interaction with the pyrimidine ring of abemaciclib. However, L152 of the C-lobe also interacted with the pyrimidine ring via π-alkyl interactions. A41 and L152 formed clips to stabilize the interaction between abemaciclib and CDK6. Those π-alkyl interactions can be verified from the distance between CB of A41 or CD1 of L152 and the pyrimidine ring (A41_CB-Abe_pyrimidine: 3.70 ± 0.26 Å and L152_CD1-Abe_pyrimidine: 3.49 ± 0.25 Å, Supplementary Figure S22). The angle between CB of A41, the pyrimidine ring of abemaciclib, and CD1 of L152 was approximately 154.96° ± 9.97°. This suggests that the three points were related to the line structure. These π-alkyl interactions resulted in A41 and L152 contributing to −1.40 and −2.36 kcal/mol for abemaciclib, forming a complex system with human CDK6. Hence, the pyrimidine ring needs to remain as an aromatic ring and maintain the hydrogen bond with the hinge loop of CDK6, such as 1H-pyrrolo [2, 3-b] pyridine, 1H-indazole, adenine, 7H-pyrrolo [2, 3-d] pyrimidin-4-amine, 1H-pyrazole-3-amine, 7H-pyrrolo [2, 3-d] pyrimidine pyrimidine-4,6-diamine, and pyridine.
Additionally, some hydrophobic residues possessed substantial subtotal binding free energies, such as V77, F98, E99, and A162. These residues form a hydrophobic cavity to accommodate fluorine on the pyrimidine ring of abemaciclib (Supplementary Figure S23). Meanwhile, the cavity was not large enough for large groups and only accommodated small groups, such as halogen, methyl, ethyl, propyl, and cyclopropyl groups.
The conformation of the piperazine ring of abemaciclib fluctuated during the simulation and was exposed to the solvent. This indicates that the piperazine group may interact with the N-lobe (K29) or C-lobe (D104) of human CDK6 (Supplementary Figure S24). K29 contributed to the binding of abemaciclib with CDK6 via the van der Waals interactions (−0.84 kcal/mol) and electrostatic interactions (−1.71 kcal/mol). Those interactions also formed with the side chain (−0.19 kcal/mol) and backbone (−0.37 kcal/mol) of K29. However, D104 was unfavorable (1.64 kcal/mol. These disadvantageous interactions were mainly obtained from polar solvation interaction (3.90 kcal/mol). In other words, if D104 forms a hydrogen bond with the piperazine group of abemaciclib, the solvated water around D104 will first repush, possibly increasing the energy barrier. Thus, water molecules around the ATP-binding site of human CDK6 should also be considered when designing novel CDK6 inhibitors.
The main contributions of H100 and V101 to human CDK6 were obtained from hydrogen bonds. K43, N150, D163, R168, and the imidazole ring of abemaciclib form a hydrogen bond network. I19, V27, K29, A41, V77, F98, L152, and A162 were favorable for abemaciclib binding with CDK6 while E99 and D104 were unfavorable. The ATP-binding pocket binds to the inhibitor mainly via hydrogen bonds and π-alkyl interactions.
3.5 Design of novel CDK6 inhibitors
Based on the MD simulation and binding free energy calculations, a binding model and key interactions for abemaciclib binding with human CDK6 were obtained (Figure 4). The benzo [d] imidazole ring of abemaciclib, defined as R1, formed hydrogen bonds with K43 and π-alkyl interactions with V27. V101 and H100 in the hinge loop formed two hydrogen bonds with amine-pyrimidine (named R2), constituting flips with A41 and L152. Meanwhile, the pyridine ring of abemaciclib, which interacted with I19 of CDK6 via π-alkyl interactions, was probed as R3. The piperazine group was also modified to R4 and exposed to the solvent.
FIGURE 4. Binding model for abemaciclib binding with human CDK6.
Cluster analysis for the 20,000 frames obtained from the last 200 ns simulation was performed to find the representative frame for the binding of abemaciclib with CDK6. There were two representative clusters of abemaciclib binding CDK6 with 82.4% and 17.6% for Clusters 1 and 2, respectively (Supplementary Figure S25). The centroid frame for every cluster was then selected as the representative complex structure and labeled Cluster1 and Cluster2. To test the binding affinity of the designed inhibitors, we extracted the CDK6 structures of Cluster1 and Cluster2 as receptor structures for molecular docking. Ensemble docking (Vilar and Costanzi, 2013; Amaro et al., 2018) was performed to assess the effects of these modifications on inhibitor binding affinity. The ensemble docking method for the designed inhibitors considered Cluster1 and Cluster2 conformations. The RMSD of 2.01 Å between Cluster1 and Cluster2 indicates that the conformational difference was insignificant.
The selected representative complex structures for abemaciclib/CDK6 were used to evaluate the power of the docking process. The native conformation in the complex structures and the docking conformation were superimposed well (Supplementary Figure S26). This result confirmed that the docking process was robust for molecular docking in the design of inhibitors.
3.5.1 One-region modification
Based on the binding model and different regions, 42 novel compounds with one-region modification (Supplementary Figure S27) were designed to inhibit the activity of CDK6. The lowest binding energy of Cluster1 and Cluster2 docking were defined as the docking scores for the designed compound (Supplementary Table S7). The score for abemaciclib (−9.46 kcal/mol) was employed as one standard value to determine whether the designed compound has a higher affinity than abemaciclib for targeting human CDK6. Most compounds had a docking score similar to that of abemaciclib, and a few designed compounds had better docking scores than abemaciclib.
As indicated, the docking score for the R1 region translation from benzene ring (abemaciclib) to pyridine ring (R1a), 1H-1, 2, 3-triazole ring (R1b), and 1H-pyrazole (R1i) decreased the binding affinity from −9.46 to −9.31, −9.31, and −8.58 kcal/mol, respectively, indicating that benzene and pyridine rings showed a similar binding affinity with V27 of CDK6. Conversely, 1-cyclopentyl-1H-[1, 2, 3] triazolo [4,5-b] pyridine (R1f) and 1-cyclopentyl-1H-pyrazolo [4,3-b] pyridine (R1m), which replaced 4-fluoro-1-isopropyl-2-methyl-1H-benzo [d] imidazole (abemaciclib), improved the activity by −10.03, −9.78, and −9.46 kcal/mol, respectively. Thus, the R1 region will be considered an isopropyl, cyclopropyl, and cyclopentyl group for the remaining similarity in binding affinity with human CDK6.
The amine-pyrimidine group plays a key role in abemaciclib binding to human CDK6 via hydrogen bonding with the hinge loop. Thus, this amine-pyrimidine group was reserved, and fluorine atoms were substituted with other groups to locate the hydrophobic cavity. Methyl (R2a), cyclopropyl (R2b), formyl (R2e), cyclopropylethan (R2f), ethyl (R2h), and hydrogen groups were used to occupy the hydrophobic cavity. Thus, the hydrogen atom, fluorine atom, and methyl group were selected as the favorable groups.
In addition, certain modifications have been made in the R3 region to improve binding affinity. R3e (fluorine group, −10.39 kcal/mol), R3f (methyl, −10.37 kcal/mol), and methoxy group (R3i, −10.28 kcal/mol) showed a higher dock score than abemaciclib (−9.46 kcal/mol). This region has been previously used to improve the inhibitory activity of protein kinases, for example, the addition of a methyl group in dasatinib. Thus, hydrogen, fluorine, and methoxy groups were considered better groups for the R3 region.
The R4 region of abemaciclib was unstable during the simulation. In molecular docking, the exposed amino group (R4a) increased the binding affinity to −10.31 kcal/mol. Isopropyl (R4b), methyl (R4g), cyclopropyl (R4i), and cyclopentyl (R4j) groups were also favorable for binding with ROCK6. Thus, the hydrogen, methyl, and isopropyl groups were considered suitable for the R4 region.
3.5.2 Combination four regions
From the single-region modification, some better groups were selected for R1 (isopropyl, cyclopropyl, and cyclopentyl groups), R2 (hydrogen, fluorine, and methyl groups), R3 (hydrogen, fluorine, and methoxy groups), and R4 (hydrogen, methyl, and isopropyl groups). In this step, a combination of R1, R2, R3, and R4 was employed to identify additional potential inhibitors. The strategies for this combination are shown in Figure 5; 81 compounds were identified in the combination strategies. The combination compounds were named as R1, R2, R3, and R4. For example, C1333 with an isopropyl group was named as R1 (1), that with a methyl group as R2 (3), that with a methoxy group as R3 (3), and that with an isopropyl group as R4 (3). These combination inhibitors have also been subjected to molecular docking using Cluster1 and Cluster2 receptors (Supplementary Table S8).
FIGURE 5. Combination strategies for designing novel inhibitors. The combination compound was named with the group numbers R1, R2, R3, and R4. For example, C1333 with an isopropyl group was named as R1, with a methyl group as R2, with a methoxy group as R3, and with an isopropyl group as R4.
When R1 remained in the isopropyl group, four compounds (C1231, C1311, C1321, and C1331) had docking scores less than −10.00 kcal/mol. Meanwhile, three compounds had dock scores < −10.00 kcal/mol (C2121, C2213, and C2331) for R1 with the cyclopropyl group. R1 with the cyclopentyl group showed a stronger binding affinity than that with the isopropyl and cyclopropyl groups. Most (12/27) designed compounds had higher binding affinities than −10.00 kcal/mol. Thus, the cyclopentyl group was the best among the isopropyl, cyclopropyl, and cyclopentyl groups. However, these three groups had similar binding conformations in the ATP-binding site.
Six compounds had better binding energy than −10.00 kcal/mol with hydrogen, fluorine, and methyl groups for R2. The trifluoromethyl group for R2 was also found for roniciclib and SY-5609, and the halogen group for abemaciclib, QHRD107, THZ1, NVP2, and THZ531 (Supplementary Figure S28). Thus, the hydrogen, halogen, and methyl groups are suitable for the R2 site and occupy the hydrophobic cavity of human CDK6. Eight, six, and four designing inhibitors had calculated binding energy above −10.00 kcal/mol for R3 with hydrogen, fluorine, and methoxy groups, respectively. Thus, the benzene ring retains a π-alkyl interaction with I19.
Hydrogen atoms with R4 can be used to form hydrogen bonds with CDK6, which contributes to inhibitor binding. Thus, in this study, the secondary amine group exposed one hydrogen atom as a hydrogen bond donor and a tertiary amine without a hydrogen atom for hydrogen bonding. The fluctuation in the conformation of the piperazine ring of abemaciclib also indicates that the flexibility of this region is insufficient for abemaciclib binding to human CDK6.
By combining R1, R2, R3, and R4, the C2213 compound was selected as the superior inhibitor. Furthermore, the variant of C2213 without the methylene group between piperazine and benzene groups also decreased the fluctuation of the piperazine group, similar to compound C2213-A (Figure 5).
3.5.3 Candidate inhibitors
Representative inhibitors, including C2213 and C2213-A, have a binding model with human CDK6 similar to that of abemaciclib (Supplementary Figure S29). These two complex systems (C2213/CDK6 and C2213-A/CDK6) were subjected to 100-ns MD simulations to calculate the binding free energies between C2213 and C2213-A and human CDK6. The RMSD vs time plot showed that the C2213/CDK6 and C213-A/CDK6 systems remained stable after 40 ns (Supplementary Figure S30). Frames 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 ns for the C2213/CDK6 and C2213-A/CDK6 systems were extracted. They revealed that C2213-A was tightly bound to the ATP-binding site of human CDK6, and C2213 conformation fluctuated in the pocket between the N-lobe and C-lobe of CDK6 (Supplementary Figures S31, S32). The side chain of the piperazine ring with C2213, which has a methylene group between the piperazine and benzene groups, exposed the piperazine group to the solvent. It disordered the hydrogen bond between the hinge loop of CDK6 and the amine-pyrimidine group of abemaciclib. In addition, the binding free energies were calculated based on the last 20 ns of the simulation trajectories (Figure 6; Supplementary Tables S9, S10). The binding free energy for C2213-A binding with CDK6 was −9.14 kcal/mol, which was 0.46 kcal/mol lower than that of abemaciclib (
FIGURE 6. Binding free energy for C2213 and C2213-A. (A) Structures of C2213 and C2213-A. (B) Binding free energy for C2213 and C2213-A binding with human CDK6. (C) The interaction model for C2213-A from 100 ns frame.
3.6 Pharmacological evaluation
C2213-A was selected as the candidate drug for the calculation studies. C2213-A {4-[1-cyclopropyl-1H-pyrazolo (4, 3-b) pyridin-6-yl]-5-fluoro-N-(4-(4-isopropylpiperazin-1-yl) phenyl) pyrimidin-2-amine} was synthesized as shown in Scheme 1. Detailed information can also be found in the Supplementary Material.
The IC50 values of abemaciclib and C2213-A against human CDK6 were evaluated using the commercial KinaseProfiler Service. The IC50 value for C2213-A was 290 nM, similar to the estimate of 238 nM for abemaciclib targeting human CDK6/cyclin D3 (Figure 7 and Supplementary Figure S33). This agreed with the binding free energy calculation with −9.14 and −9.74 kcal/mol for C2213-A and abemaciclib, respectively. Thus, C2213-A was identified as an abemaciclib derivative for continued optimization to target CDK6.
FIGURE 7. Pharmacological evaluation of C2213-A and abemaciclib. (A) Structures of abemaciclib and C2231-A. (B) Biochemical and cellular profiling of abemaciclib and C2231-A. (C) Colony-formation assays for MDA-MB-231 with different concentrations of C2231-A. (D) Percentage of apoptotic cells after treatment with different concentrations of compound C2213-A for 48 h (E) MDA-MB-231 cell lines exhibited obvious G1 arrest and a decrease in S phase cells after incubation with different doses of compound C2213-A or abemaciclib. (F) Compound C2213-A inhibited the phosphorylation of Rb at S780 in a dose-dependent manner in MDA-MB-231.
In addition, the antiproliferative effects of abemaciclib and C2213-A were evaluated in MDA-MB-231 (a human breast cancer cell line). The antiproliferative activity of C2213-A was significantly higher than that of abemaciclib alone, with IC50 values of 2.95 ± 0.15 µM and 3.51 ± 0.28 µM, respectively (Figure 7). In MCF-7 cells, the IC50 values of abemaciclib and C2213-A were 0.69 ± 0.04 µM and 0.77 ± 0.05 µM, respectively (Supplementary Figure S33). To further demonstrate the broad applications of our findings, the inhibitory activities of abemaciclib and C2213-A were tested using T-47D, MDA-MB-453, and MDA-MB-468 cell lines. C2231-A showed a better inhibitory effect than that of abemaciclib in T-47D (IC50 values were 1.13 ± 0.15 µM and 3.70 ± 0.18 µM, respectively). However, C2231-A was less effective than abemaciclib in inhibiting the growth of MDA-MB-468 cells (3.61 ± 0.11 µM and 2.20 ± 0.11 µM, respectively). Moreover, C2231-A showed activity similar to that of abemaciclib in MDA-MB-453 cells (1.28 ± 0.16 µM and 2.34 ± 0.08 µM, respectively). These observations indicated that the inhibitory activities of C2231-A and abemaciclib differ among cell lines and suggest that C2213-A inhibits the proliferation of MDA-MB-231, MCF-7, T-47D, MDA-MB-453, and MDA-MB-468 cells.
A clone formation assay was conducted using MDA-MB-231 cells to investigate the inhibitory activity of C2213-A, with abemaciclib for comparison (Supplementary Figure S34). The proliferation of MDA-MB-231 cells was impaired (inhibition rate >45%) by abemaciclib (62.5 nM) and C2231-A (125 nM) compared with that in the control group (Figure 7). These results suggest that C2213-A inhibits the activity of MDA-MB-231 cells.
An annexin V-FITC/propidium iodide staining assay using human breast tumor MDA-MB-231 cells showed that approximately 83% of treated cells underwent apoptosis after a 48-h treatment with compound C2213-A at a concentration of 5 μM (Supplementary Figure S35). These findings confirmed that the anti-proliferative effect of C2213-A is most likely due to the induction of apoptosis.
We also compared the effects of compound C2213-A and abemaciclib on the cell cycle profile of MDA-MB-231 cells by flow cytometry (Supplementary Figure S36). These findings revealed that both compounds elicited a dose-dependent effect, leading to the accumulation of cells in the G1 phase, as expected.
The effects of compound C2213-A on CDK4/6 and the phosphorylation of Rb (S780) were assessed by a Western blot assay on MDA-MB-231 cells. Compound C2213-A blocked the CDK4/6/Rb/E2F signaling pathway in a dose-dependent manner after 24 h of incubation (Supplementary Figure S37). Treatment with compound C2213-A resulted in the dose-dependent blockage of Rb phosphorylation, consistent with its targeted mechanism of action on CDK4/6 in cells.
4 Conclusion
CDK4/6 is disordered in numerous cancers, such as breast cancer, osteosarcoma, and acute megakaryoblastic leukemia. CDK4/6 is considered an effective anticancer drug target. Three CDK4/6 inhibitors have been approved; however, there is still a gap between the clinical requirements and approved CDK4/6 drugs. Thus, selective and oral CDK4/6 inhibitors urgently need to be developed, particularly for monotherapy. This study investigated the interaction between abemaciclib and human CDK6 using MD simulations. Simultaneously, the binding free energy and energy decomposition were calculated using the MM/GBSA approach. V101 and H100 formed stable hydrogen bonds with the amine pyrimidine group of abemaciclib. K43 interacted with an imidazole ring via an unstable hydrogen bond. Meanwhile, I19 (pyridine ring), V27 (imidazole ring), A41 (pyrimidine ring), and L152 (pyrimidine ring) interacted with abemaciclib via an π-alkyl force. Based on the binding model of abemaciclib/CDK6. Abemaciclib treatment was divided into four regions; 43 compounds were designed with only one region modification. Three favorable groups for each region were selected from the docking score for single-region modification. Subsequently, a combination of the favorable groups within the four regions was performed, and 81 compounds were obtained. Among them, C2231 emerged as a better compound from the four-region combination strategy. C2231-A, which lacks a methylene group compared to C2213, showed superior flexibility reduction. Additionally, kinase profiling and cell profiling of MDA-MB-231 and MCF-7 cells were performed after synthesizing C2231-A. C2231-A caused greater inhibition than did abemaciclib in MDA-MB-231 cells. Based on MD simulation, one candidate compound was designed, which showed a significant inhibitory effect on a human breast cancer cell line. However, computational tools and experimental methods must be used to identify more compounds based on C2231-A. In the future, C2213-A derivatives will be studied using molecular docking, MD simulation, binding free energy calculation, kinase, and cellular profiling.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number (s) can be found in the article/Supplementary Material.
Author contributions
Conceptualization: YtZ, XdL, and MS; methodology: YtZ, XdL, and MS; software: MS, HC, and PL validation: YZ and XZ; formal analysis: YY and MS; investigation: YtZ, XdL, and MS; resources: HW and WX; data curation: WX and HL; writing—original draft preparation: MS; writing—review and editing: WX, MS, and XaL; visualization: HL, HC, and MS; supervision: WX, MS, and XaL; project administration: MS and WX; funding acquisition: XZ, HW, and MS. All the authors have read and agreed to the published version of the manuscript. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication. All authors contributed to the article and approved the submitted version.
Funding
The present work was supported by Natural Science Foundation of Sichuan, China (No. 2022NSFSC1467), the National Natural Science Foundation of China (No. 22203056), the Science and Technology Planning Project of Sichuan Province (No. 2021YJ0483), NHC Key Laboratory of Nuclear Technology Medical Transformation (Mianyang Central Hospital) (Nos 2022HYX001 and 2022HYX005), GHfund B (ghfund202302022385), and the Postdoctoral Interdisciplinary Innovation Fund of Sichuan University.
Acknowledgments
Certain data were obtained from the National Supercomputing Center of Guangzhou, the Supercomputing Center of Sichuan University, and the Chengdu Supercomputing Center. The authors thank Editage (www.editage.cn) for English language editing.
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.
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.2023.1154654/full#supplementary-material
References
Abraham, J. (2018). Abemaciclib becomes first cdk inhibitor to clinch single-agent approval for breast cancer. J. Community Support. Oncol. 16 (1), E2–E4. doi:10.12788/jcso.0378
Amaro, R. E., Baudry, J., Chodera, J., Demir, O., McCammon, J. A., Miao, Y. L., et al. (2018). Ensemble docking in drug discovery. Biophys. J. 114 (10), 2271–2278. doi:10.1016/j.bpj.2018.02.038
Arter, C., Trask, L., Ward, S., Yeoh, S., and Bayliss, R. (2022). Structural features of the protein kinase domain and targeted binding by small-molecule inhibitors. J. Biol. Chem. 298 (8), 102247. doi:10.1016/j.jbc.2022.102247
Bayly, C. I., Cieplak, P., Cornell, W. D., and Kollman, P. A. (1993). A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: The resp model. J. Phys. Chem. 97 (40), 10269–10280. doi:10.1021/j100142a004
Berendsen, H. J. C., Postma, J. P. M., Vangunsteren, W. F., Dinola, A., and Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81 (8), 3684–3690. doi:10.1063/1.448118
Bhattacharya, A., Sen Guha, P., Chowdhury, N., Bagchi, A., and Guha, D. (2023). Virtual screening and molecular docking of flavone derivatives as a potential anticancer drug in the presence of dexamethasone. Biointerface Res. Appl. Chem. 13 (3), 19. doi:10.33263/briac1333.23.215
Bhurta, D., and Bharate, S. B. (2022). Analyzing the scaffold diversity of cyclin-dependent kinase inhibitors and revisiting the clinical and preclinical pipeline. Med. Res. Rev. 42 (2), 654–709. doi:10.1002/med.21856
Braal, C. L., Jongbloed, E. M., Wilting, S. M., Mathijssen, R. H. J., Koolen, S. L. W., and Jager, A. (2021). Inhibiting cdk4/6 in breast cancer with palbociclib, ribociclib, and abemaciclib: Similarities and differences. Drugs 81 (3), 317–331. doi:10.1007/s40265-020-01461-2
Bronner, S. M., Merrick, K. A., Murray, J., Salphati, L., Moffat, J. G., Pang, J., et al. (2019). Design of a brain-penetrant cdk4/6 inhibitor for glioblastoma. Bioorg. Med. Chem. Lett. 29 (16), 2294–2301. doi:10.1016/j.bmcl.2019.06.021
Case, D. A., Ben-Shalom, I. Y., Brozell, S. R., Cerutti, D. S., Cheatham, T. E., Cruzeiro, V. W. D., et al. (2020). Amber 2020. san francisco: University of California.
Cheatham, T., Galindo, R., and Roe, D. (2019). Parallel analysis of large ensembles of molecular dynamics simulation derived trajectories with the open-source cpptraj tools. Abstr. Pap. Am. Chem. Soc. 257, 1.
Chen, P., Lee, N. V., Hu, W. Y., Xu, M. R., Ferre, R. A., Lam, H., et al. (2016). Spectrum and degree of cdk drug interactions predicts clinical performance. Mol. Cancer Ther. 15 (10), 2273–2281. doi:10.1158/1535-7163.Mct-16-0300
Chen, W. J., Ji, M. H., Cheng, H., Zheng, M. M., Xia, F., Min, W. J., et al. (2022). Discovery, optimization, and evaluation of selective cdk4/6 inhibitors for the treatment of breast cancer. J. Med. Chem. 65 (22), 15102–15122. doi:10.1021/acs.jmedchem.2c00947
Cho, Y. S., Angove, H., Brain, C., Chen, C. H. T., Cheng, H., Cheng, R., et al. (2012). Fragment-based discovery of 7-azabenzimidazoles as potent, highly selective, and orally active cdk4/6 inhibitors. ACS Med. Chem. Lett. 3 (6), 445–449. doi:10.1021/ml200241a
Cho, Y. S., Borland, M., Brain, C., Chen, C. H. T., Cheng, H., Chopra, R., et al. (2010). 4-(pyrazol-4-yl)-pyrimidines as selective inhibitors of cyclin-dependent kinase 4/6. J. Med. Chem. 53 (22), 7938–7957. doi:10.1021/jm100571n
Choi, Y. J., and Anders, L. (2014). Signaling through cyclin d-dependent kinases. Oncogene 33 (15), 1890–1903. doi:10.1038/onc.2013.137
Dhillon, S. (2015). Palbociclib: First global approval. Drugs 75 (5), 543–551. doi:10.1007/s40265-015-0379-9
Dowless, M., Lowery, C. D., Shackleford, T., Renschler, M., Stephens, J., Flack, R., et al. (2018). Abemaciclib is active in preclinical models of ewing sarcoma via multipronged regulation of cell cycle, DNA methylation, and interferon pathway signaling. Clin. Cancer Res. 24 (23), 6028–6039. doi:10.1158/1078-0432.Ccr-18-1256
Ettl, T., Schulz, D., and Bauer, R. J. (2022). The renaissance of cyclin dependent kinase inhibitors. Cancers 14 (2), 293. doi:10.3390/cancers14020293
Feller, S. E., Zhang, Y. H., Pastor, R. W., and Brooks, B. R. (1995). Constant pressure molecular dynamics simulation: The Langevin piston method. J. Chem. Phys. 103 (11), 4613–4621. doi:10.1063/1.470648
Fensome, A., Ambler, C. M., Arnold, E., Banker, M. E., Brown, M. F., Chrencik, J., et al. (2018). Dual inhibition of tyk2 and jak1 for the treatment of autoimmune diseases: Discovery of ((s)-2,2-difluorocyclopropyl)((1r,5s)-3-(2-(1-methyl-1h-pyrazol-4-yl)ami no)pyrimidin-4-yl)-3,8-diazabicyclo 3.2.1 octan-8-yl)methanone (pf-06700841). J. Med. Chem. 61 (19), 8597–8612. doi:10.1021/acs.jmedchem.8b00917
Gaillard, T., and Simonson, T. (2014). Pairwise decomposition of an mmgbsa energy function for computational protein design. J. Comput. Chem. 35 (18), 1371–1387. doi:10.1002/jcc.23637
Gelbert, L. M., Cai, S. F., Lin, X., Sanchez-Martinez, C., del Prado, M., Lallena, M. J., et al. (2014). Preclinical characterization of the cdk4/6 inhibitor ly2835219: In-vivo cell cycle-dependent/independent anti-tumor activities alone/in combination with gemcitabine. Invest. New Drugs 32 (5), 825–837. doi:10.1007/s10637-014-0120-7
Genheden, S., and Ryde, U. (2015). The mm/pbsa and mm/gbsa methods to estimate ligand-binding affinities. Expert. Opin. Drug Discov. 10 (5), 449–461. doi:10.1517/17460441.2015.1032936
Ghafouri-Fard, S., Khoshbakht, T., Hussen, B. M., Dong, P. X., Gassler, N., Taheri, M., et al. (2022). A review on the role of cyclin dependent kinases in cancers. Cancer Cell Int. 22 (1), 325. doi:10.1186/s12935-022-02747-z
Guo, C., Li, Q. X., Xiao, J. J., Ma, F., Xia, X., and Shi, M. S. (2022). Identification of defactinib derivatives targeting focal adhesion kinase using ensemble docking, molecular dynamics simulations and binding free energy calculations. J. Biomol. Struct. Dyn. 17, 1–17. doi:10.1080/07391102.2022.2135601
Honig, B., and Nicholls, A. (1995). Classical electrostatics in biology and chemistry. Science 268 (5214), 1144–1149. doi:10.1126/science.7761829
Hope, I., Endicott, J. A., and Watt, J. E. (2023). Emerging approaches to cdk inhibitor development, a structural perspective. RSC Chem. Biol. 19, 146–164. doi:10.1039/d2cb00201a
Ingham, M., and Schwartz, G. K. (2017). Cell-cycle therapeutics come of age. J. Clin. Oncol. 35(25), 2949–2959. doi:10.1200/jco.2016.69.0032
Ji, J., Feng, J., Pan, G., Lv, M. X., Lv, J. Y., Liu, W. W., et al. (2022). A novel cdk4/6 inhibitor wxj-111, against breast cancer through mediated cell-cycle arrest and apoptosis. Acta Pol. Pharm. 79 (3), 353–365. doi:10.32383/appdr/152215
Jorda, R., Havlicek, L., Perina, M., Vojackova, V., Pospisil, T., Djukic, S., et al. (2022). 3,5,7-substituted pyrazolo 4,3-d pyrimidine inhibitors of cyclin-dependent kinases and cyclin k degraders. J. Med. Chem. 65 (13), 8881–8896. doi:10.1021/acs.jmedchem.1c02184
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79 (2), 926–935. doi:10.1063/1.445869
Kaltheuner, I. H., Anand, K., Moecking, J., Duster, R., Wang, J. H., Gray, N. S., et al. (2021). Abemaciclib is a potent inhibitor of dyrk1a and hip kinases involved in transcriptional regulation. Nat. Commun. 12 (1), 6607. doi:10.1038/s41467-021-26935-z
Kargbo, R. B. (2022). Selective cyclin-dependent kinase inhibitors and their application in cancer therapy. ACS Med. Chem. Lett. 13 (10), 1561–1563. doi:10.1021/acsmedchemlett.2c00417
Kim, E. S. (2017). Abemaciclib: First global approval. Drugs 77 (18), 2063–2070. doi:10.1007/s40265-017-0840-z
Kosovec, J. E., Zaidi, A. H., Omstead, A. N., Matsui, D., Biedka, M. J., Cox, E. J., et al. (2017). Cdk4/6 dual inhibitor abemaciclib demonstrates compelling preclinical activity against esophageal adenocarcinoma: A novel therapeutic option for a deadly disease. Oncotarget 8 (59), 100421–100432. doi:10.18632/oncotarget.22244
Ladduwahetty, T., Lee, M. R., Maillard, M. C., Cachope, R., Todd, D., Barnes, M., et al. (2022). Identification of a potent, selective, and brain-penetrant rho kinase inhibitor and its activity in a mouse model of huntington's disease. J. Med. Chem. 27, 9819–9845. doi:10.1021/acs.jmedchem.2c00474
Lee, M. S., Salsbury, F. R., and Olson, M. A. (2004). An efficient hybrid explicit/implicit solvent method for biomolecular simulations. J. Comput. Chem. 25 (16), 1967–1978. doi:10.1002/jcc.20119
Li, Z. H., Wang, X. H., Eksterowicz, J., Gribble, M. W., Alba, G. Q., Ayres, M., et al. (2014). Discovery of amg 925, a flt3 and cdk4 dual kinase inhibitor with preferential affinity for the activated state of flt3. J. Med. Chem. 57 (8), 3430–3449. doi:10.1021/jm500118j
Lim, E., Boyle, F., Okera, M., Loi, S., Goksu, S. S., van Hal, G., et al. (2022). An open label, randomized phase 2 trial assessing the impact of food on the tolerability of abemaciclib in patients with advanced breast cancer. Breast Cancer Res. Treat. 195 (3), 275–287. doi:10.1007/s10549-022-06690-5
Lu, H. S., Chang, D. J., Baratte, B., Meijer, L., and Schulze-Gahmen, U. (2005). Crystal structure of a human cyclin-dependent kinase 6 complex with a flavonol inhibitor, fisetin. J. Med. Chem. 48 (3), 737–743. doi:10.1021/jm049353p
Lu, H., and Schulze-Gahmen, U. (2006). Toward understanding the structural basis of cyclin-dependent kinase 6 specific inhibition. J. Med. Chem. 49 (13), 3826–3831. doi:10.1021/jm0600388
Macaluso, M., Montanari, M., and Giordano, A. (2006). Rb family proteins as modulators of gene expression and new aspects regarding the interaction with chromatin remodeling enzymes. Oncogene 25 (38), 5263–5267. doi:10.1038/sj.onc.1209680
Manohar, S. M. (2022). Cyclin-dependent kinases as potential targets for colorectal cancer: Past, present and future. Future Med. Chem. 14 (14), 1087–1105. doi:10.4155/fmc-2022-0064
Martyna, G. J., Tobias, D. J., and Klein, M. L. (1994). Constant pressure molecular dynamics algorithms. J. Chem. Phys. 101 (5), 4177–4189. doi:10.1063/1.467468
Matson, J. P., and Cook, J. G. (2017). Cell cycle proliferation decisions: The impact of single cell analyses. Febs J. 284 (3), 362–375. doi:10.1111/febs.13898
Miller, B. R., McGee, T. D., Swails, J. M., Homeyer, N., Gohlke, H., and Roitberg, A. E. (2012). Mmpbsa.Py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 8(9), 3314–3321.doi:10.1021/Ct300418h
Mitarai, Y., Tsubata, Y., Hyakudomi, M., and Isobe, T. (2022). Drug-induced eosinophilic pneumonia as an adverse event of abemaciclib. Cureus J. Med. Sci. 14 (1), e21741. doi:10.7759/cureus.21741
Nardone, V., Porta, C., Giannicola, R., Correale, P., and Mutti, L. (2022). Abemaciclib for malignant pleural mesothelioma. Lancet Oncol. 23 (6), E237. doi:10.1016/S1470-2045(22)00276-5
Nawa, H., Niimura, T., Yagi, K., Goda, M., Zamami, Y., and Ishizawa, K. (2022). Evaluation of potential complication of interstitial lung disease with abemaciclib and palbociclib treatments. Cancer Rep.-US 5 (1), e1402. doi:10.1002/cnr2.1402
Naz, S., Sowers, A., Choudhuri, R., Wissler, M., Gamson, J., Mathias, A., et al. (2018). Abemaciclib, a selective cdk4/6 inhibitor, enhances the radiosensitivity of non-small cell lung cancer in vitro and in vivo. Clin. Cancer Res. 24 (16), 3994–4005. doi:10.1158/1078-0432.Ccr-17-3575
Onufriev, A. V., and Case, D. A. (2019). “Generalized born implicit solvent models for biomolecules,”. Annual review of biophysics. 48, 275–296. doi:10.1146/annurev-biophys-052118-115325
Pang, X. J., Liu, X. J., Liu, Y., Liu, W. B., Li, Y. R., Yu, G. X., et al. (2021). Drug discovery targeting focal adhesion kinase (fak) as a promising cancer therapy. Molecules 26 (14), 4250. doi:10.3390/molecules26144250
Qi, K. M., Hu, X. T., Yu, X. R., Cheng, H., Wang, C. Q., Wang, S. J., et al. (2022). Targeting cyclin-dependent kinases 4/6 inhibits survival of megakaryoblasts in acute megakaryoblastic leukaemia. Leuk. Res. 120, 106920. doi:10.1016/j.leukres.2022.106920
Roe, D. R., and Cheatham, T. E. (2018). Parallelization of cpptraj enables large scale analysis of molecular dynamics trajectory data. J. Comput. Chem. 39 (25), 2110–2117. doi:10.1002/jcc.25382
Roskoski, R. (2016). Classification of small molecule protein kinase inhibitors based upon the structures of their drug-enzyme complexes. Pharmacol. Res. 103, 26–48. doi:10.1016/j.phrs.2015.10.021
Sanchez-Martinez, C., Gelbert, L. M., Lallena, M. J., and de Dios, A. (2015). Cyclin dependent kinase (cdk) inhibitors as anticancer drugs. Bioorg. Med. Chem. Lett. 25 (17), 3420–3435. doi:10.1016/j.bmcl.2015.05.100
Sanchez-Martinez, C., Lallena, M. J., Sanfeliciano, S. G., and de Dios, A. (2019). Cyclin dependent kinase (cdk) inhibitors as anticancer drugs: Recent advances (2015-2019). Bioorg. Med. Chem. Lett. 29 (20), 126637. doi:10.1016/j.bmcl.2019.126637
Shi, M., Chen, T., Wei, S., Zhao, C., Zhang, X., Li, X., et al. (2022a). Molecular docking, molecular dynamics simulations, and free energy calculation insights into the binding mechanism between vs-4718 and focal adhesion kinase. ACS Omega 7 (36), 32442–32456. doi:10.1021/acsomega.2c03951
Shi, M., Liu, J., Fu, S., Pei, H., Peng, B., Wen, Y., et al. (2022b). Structural insights into the interactions of belumosudil with rho-associated coiled-coil containing protein kinases 1 and 2 based on molecular docking, molecular dynamics simulations, and free energy calculations. J. Comput. Biophysics Chem., 1–22. doi:10.1142/S2737416523500163
Shi, M., Wang, L., Liu, K., Chen, Y., Hu, M., Yang, L., et al. (2022c). Molecular dynamics simulations of the conformational plasticity in the active pocket of salt-inducible kinase 2 (sik2) multi-state binding with bosutinib. Comput. Struct. Biotechnol. J. 20, 2574–2586. doi:10.1016/j.csbj.2022.05.039
Shi, M., and Xu, D. (2019). Molecular dynamics investigations suggest a non-specific recognition strategy of 14-3-3σ protein by tweezer: Implication for the inhibition mechanism. Front. Chem. 7 (237), 237. doi:10.3389/fchem.2019.00237
Shi, Z. F., Tian, L., Qiang, T. T., Li, J. Y., Xing, Y., Ren, X. D., et al. (2022d). From structure modification to drug launch: A systematic reviewof the ongoing development of cyclin-dependent kinase inhibitorsfor multiple cancer therapy. J. Med. Chem. 65 (9), 6390–6418. doi:10.1021/acs.jmedchem.1c02064
Small, J., Washburn, E., Millington, K., Zhu, J. J., and Holder, S. L. (2017). The addition of abemaciclib to sunitinib induces regression of renal cell carcinoma xenograft tumors. Oncotarget 8 (56), 95116–95134. doi:10.18632/oncotarget.19618
Sokolsky, A., Winterton, S., Kennedy, K., Drake, K., Stump, K., Huo, L., et al. (2022). Discovery of 5,7-dihydro-6h-pyrrolo 2,3-d pyrimidin-6-ones as highly selective cdk2 inhibitors. ACS Med. Chem. Lett. 13 (11), 1797–1804. doi:10.1021/acsmedchemlett.2c00408
Spring, L. M., Wander, S. A., Andre, F., Moy, B., Turner, N. C., and Bardia, A. (2020). Cyclin-dependent kinase 4 and 6 inhibitors for hormone receptor-positive breast cancer: Past, present, and future. Lancet 395 (10226), 817–827. doi:10.1016/S0140-6736(20)30165-3
Srinivasan, J., Cheatham, T. E., Cieplak, P., Kollman, P. A., and Case, D. A. (1998). Continuum solvent studies of the stability of DNA, rna, and phosphoramidate - DNA helices. J. Am. Chem. Soc. 120 (37), 9401–9409. doi:10.1021/ja981844+
Stavenger, R. A. (2008). “Rho-kinase inhibitors for cardiovascular disease,”. Annual reports in medicinal chemistry. 43, 87–102. doi:10.1161/CIRCRESAHA.115.306532
Sun, S., Liu, L., Song, H. M., and Li, H. (2022). Pharmacokinetic study on the co-administration of abemaciclib and astragaloside iv in rats. Pharm. Biol. 60 (1), 1944–1948. doi:10.1080/13880209.2022.2125539
Syed, Y. Y. (2017). Ribociclib: First global approval. Drugs 77 (7), 799–807. doi:10.1007/s40265-017-0742-0
Tesch, R., Rak, M., Raab, M., Berger, L. M., Kronenberger, T., Joerger, A. C., et al. (2021). Structure-based design of selective salt-inducible kinase inhibitors. J. Med. Chem. 64 (12), 8142–8160. doi:10.1021/acs.jmedchem.0c02144
Tian, C., Kasavajhala, K., Belfon, K. A. A., Raguette, L., Huang, H., Migues, A. N., et al. (2020). Ff19sb: Amino-acid-specific protein backbone parameters trained against quantum mechanics energy surfaces in solution. J. Chem. Theory Comput. 16 (1), 528–552. doi:10.1021/acs.jctc.9b00591
Torres-Guzman, R., Calsina, B., Hermoso, A., Baquero, C., Alvarez, B., Amat, J., et al. (2017). Preclinical characterization of abemaciclib in hormone receptor positive breast cancer. Oncotarget 8 (41), 69493–69507. doi:10.18632/oncotarget.17778
Tse, A., and Verkhivker, G. M. (2015). Molecular dynamics simulations and structural network analysis of c-abl and c-src kinase core proteins: Capturing allosteric mechanisms and communication pathways from residue centrality. J. Chem. Inf. Model. 55 (8), 1645–1662. doi:10.1021/acs.jcim.5b00240
Vilar, S., and Costanzi, S. (2013). “Application of Monte Carlo-based receptor ensemble docking to virtual screening for gpcr ligands,” Methods Enzymol. 522, 263–278. doi:10.1016/B978-0-12-407865-9.00014-5
Wander, S. A., O'Brien, N., Litchfield, L. M., O'Dea, D., Guimaraes, C. M., Slamon, D. J., et al. (2022). Targeting cdk4 and 6 in cancer therapy: Emerging preclinical insights related to abemaciclib. Oncologist 27 (10), 811–821. doi:10.1093/oncolo/oyac138
Wang, D. L., and Bao, H. Q. (2022). Abemaciclib is synergistic with doxorubicin in osteosarcoma pre-clinical models via inhibition of cdk4/6-cyclin d-rb pathway. Cancer Chemother. Pharmacol. 89 (1), 31–40. doi:10.1007/s00280-021-04363-6
Wang, J. M., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2004). Development and testing of a general amber force field. J. Comput. Chem. 25 (9), 1157–1174. doi:10.1002/jcc.20035
Wang, J. Y., Chen, Q., Wang, M., and Zhong, C. (2017). The opening/closure of the p-loop and hinge of bcr-abl1 decodes the low/high bioactivities of dasatinib and axitinib. Phys. Chem. Chem. Phys. 19 (33), 22444–22453. doi:10.1039/c7cp03443a
Wang, L. F., Lu, D., Wang, Y., Xu, X. Y., Zhong, P. H., and Yang, Z. Y. (2023). Binding selectivity-dependent molecular mechanism of inhibitors towards cdk2 and cdk6 investigated by multiple short molecular dynamics and free energy landscapes. J. Enzym. Inhib. Med. Chem. 38 (1), 84–99. doi:10.1080/14756366.2022.2135511
Watson, N. W., Wander, S. A., Shatzel, J. J., and Al-Samkari, H. (2022). Venous and arterial thrombosis associated with abemaciclib therapy for metastatic breast cancer. Cancer 128 (17), 3224–3232. doi:10.1002/cncr.34367
Watt, A. C., and Goel, S. (2022). Cellular mechanisms underlying response and resistance to cdk4/6 inhibitors in the treatment of hormone receptor-positive breast cancer. Breast Cancer Res. 24 (1), 17. doi:10.1186/s13058-022-01510-6
Wesierska-Gadek, J., Chamrad, I., and Krystof, V. (2009). Novel potent pharmacological cyclin-dependent kinase inhibitors. Future Med. Chem. 1 (9), 1561–1581. doi:10.4155/fmc.09.110
Wright, M. D., and Md, J. A. (2021). Preclinical discovery and development of abemaciclib used to treat breast cancer. Expert. Opin. Drug Discov. 16 (5), 485–496. doi:10.1080/17460441.2021.1853097
Xie, Q., Zheng, H. R., and Li, Q. (2022a). Cost-effectiveness analysis of abemaciclib plus fulvestrant versus placebo plus fulvestrant in patients with hormone receptor-positive, erbb2-negative breast cancer. Breast Care 17 (3), 237–243. doi:10.1159/000518551
Xie, Z. L., Hou, S. Z., Yang, X. X., Duan, Y. J., Han, J. H., Wang, Q., et al. (2022b). Lessons learned from past cyclin-dependent kinase drug discovery efforts. J. Med. Chem. 65 (9), 6356–6389. doi:10.1021/acs.jmedchem.1c02190
Xing, L., Klug-Mcleod, J., Rai, B., and Lunney, E. A. (2015). Kinase hinge binding scaffolds and their hydrogen bond patterns. Bioorg. Med. Chem. 23 (19), 6520–6527. doi:10.1016/j.bmc.2015.08.006
Yang, L. J., Zhao, W. H., and Liu, Q. (2014). Detailed conformation dynamics and activation process of wild type c-Abl and T315I mutant. J. Mol. Struct. 1075, 292–298. doi:10.1016/j.molstruc.2014.06.087
Yousuf, M., Alam, M., Shamsi, A., Khan, P., Hasan, G. M., Haque, Q. M. R., et al. (2022). Structure-guided design and development of cyclin-dependent kinase 4/6 inhibitors: A review on therapeutic implications. Int. J. Biol. Macromol. 218, 394–408. doi:10.1016/j.ijbiomac.2022.07.156
Yu, B., Du, Z. K., Zhang, Y. M., Li, Z. Y., and Bian, J. L. (2022). Small-molecule degraders of cyclin-dependent kinase protein: A review. Future Med. Chem. 14 (3), 167–185. doi:10.4155/fmc-2021-0154
Yuan, K., Kuang, W. B., Chen, W. J., Ji, M. H., Min, W. J., Zhu, Y. S., et al. (2022). Discovery of novel and orally bioavailable cdk 4/6 inhibitors with high kinome selectivity, low toxicity and long-acting stability for the treatment of multiple myeloma. Eur. J. Med. Chem. 228, 114024. doi:10.1016/j.ejmech.2021.114024
Yuan, K., Wang, X., Dong, H. J., Min, W. J., Hao, H. P., and Yang, P. (2021). Selective inhibition of cdk4/6: A safe and effective strategy for developing anticancer drugs. Acta Pharm. Sin. B 11 (1), 30–54. doi:10.1016/j.apsb.2020.05.001
Zabihi, M., Lotfi, R., Yousefi, A. M., and Bashash, D. (2022). Cyclins and cyclin-dependent kinases: From biology to tumorigenesis and therapeutic opportunities. J. Cancer Res. Clin. Oncol. 22, 1585–1606. doi:10.1007/s00432-022-04135-6
Keywords: abemaciclib, cyclin-dependent kinase 4, cyclin-dependent kinase 6, inhibitor, molecular dynamics simulation, binding free energy
Citation: Zhou Y, Li X, Luo P, Chen H, Zhou Y, Zheng X, Yin Y, Wei H, Liu H, Xia W, Shi M and Li X (2023) Identification of abemaciclib derivatives targeting cyclin-dependent kinase 4 and 6 using molecular dynamics, binding free energy calculation, synthesis, and pharmacological evaluation. Front. Pharmacol. 14:1154654. doi: 10.3389/fphar.2023.1154654
Received: 13 February 2023; Accepted: 02 May 2023;
Published: 10 May 2023.
Edited by:
Kehao Zhao, Yantai University, ChinaCopyright © 2023 Zhou, Li, Luo, Chen, Zhou, Zheng, Yin, Wei, Liu, Xia, Shi and Li. 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: Wen Xia, eHh4eXhra0AxNjMuY29t; Mingsong Shi, dGhlcm90eW9udGhAMTYzLmNvbQ==
†These authors have contributed equally to this work