Molecular dynamics simulations of the conformational plasticity in the active pocket of salt-inducible kinase 2 (SIK2) multi-state binding with bosutinib

Advanced modeling of salt-inducible kinase (SIK) inhibitors incorporating protein flexibility through molecular dynamics and cross-docking

Salt-inducible kinases (SIK1, SIK2, and SIK3) regulate metabolism and immune responses, making them promising targets for inflammatory and autoimmune diseases. Understanding inhibitor selectivity among isoforms is crucial for therapeutic development. In this study, 44 compounds were investigated as SIK inhibitors using molecular modeling. A flexible treatment of the kinases via molecular dynamics (MD) simulations captured binding site conformational changes, followed by molecular docking to generate protein kinase (PK)-ligand complex models. Ligand orientations were validated against crystallographic data using LigRMSD and interaction fingerprints (IFPs). A genetic algorithm was applied to select conformations that maximize correlation between docking energies and biological activities, yielding R² values of 0.821, 0.646, and 0.620 for SIK1, SIK2, and SIK3, respectively. Our results highlight the importance of protein flexibility in achieving accurate correlations between docking energies and experimental pIC50 values, enhancing inhibitor selectivity predictions.

Identification of defactinib derivatives targeting focal adhesion kinase using ensemble docking, molecular dynamics simulations and binding free energy calculations

Focal adhesion kinase (FAK) belongs to the nonreceptor tyrosine kinases, which selectively phosphorylate tyrosine residues on substrate proteins. FAK is associated with bladder, esophageal, gastric, neck, breast, ovarian and lung cancers. Thus, FAK has been considered as a potential target for tumor treatment. Currently, there are six adenosine triphosphate (ATP)-competitive FAK inhibitors tested in clinical trials but no approved inhibitors targeting FAK. Defactinib (VS-6063) is a second-generation FAK inhibitor with an IC50 of 0.6 nM. The binding model of VS-6063 with FAK may provide a reference model for developing new antitumor FAK-targeting drugs. In this study, the VS-6063/FAK binding model was constructed using ensemble docking and molecular dynamics simulations. Furthermore, the molecular mechanics/generalized Born (GB) surface area (MM/GBSA) method was employed to estimate the binding free energy between VS-6063 and FAK. The key residues involved in VS-6063/FAK binding were also determined using per-residue energy decomposition analysis. Based on the binding model, VS-6063 could be separated into seven regions to enhance its binding affinity with FAK. Meanwhile, 60 novel defactinib-based compounds were designed and verified using ensemble docking. Overall, the present study improves our understanding of the binding mechanism of human FAK with VS-6063 and provides new insights into future drug designs targeting FAK.Communicated by Ramaswamy H. Sarma.

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

Identification of abemaciclib derivatives targeting cyclin-dependent kinase 4 and 6 using molecular dynamics, binding free energy calculation, synthesis, and pharmacological evaluation

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.

Interactions between curcumin and human salt-induced kinase 3 elucidated from computational tools and experimental methods

Natural products are widely used for treating mitochondrial dysfunction-related diseases and cancers. Curcumin, a well-known natural product, can be potentially used to treat cancer. Human salt-induced kinase 3 (SIK3) is one of the target proteins for curcumin. However, the interactions between curcumin and human SIK3 have not yet been investigated in detail. In this study, we studied the binding models for the interactions between curcumin and human SIK3 using computational tools such as homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations. The open activity loop conformation of SIK3 with the ketoenol form of curcumin was the optimal binding model. The I72, V80, A93, Y144, A145, and L195 residues played a key role for curcumin binding with human SIK3. The interactions between curcumin and human SIK3 were also investigated using the kinase assay. Moreover, curcumin exhibited an IC₅₀ (half-maximal inhibitory concentration) value of 131 nM, and it showed significant antiproliferative activities of 9.62 ± 0.33 µM and 72.37 ± 0.37 µM against the MCF-7 and MDA-MB-23 cell lines, respectively. This study provides detailed information on the binding of curcumin with human SIK3 and may facilitate the design of novel salt-inducible kinases inhibitors.