Exploring the stability of inhibitor binding to SIK2 using molecular dynamics simulation and binding free energy calculation

SIK2-positive tumor epithelial cells in breast cancer maybe potential anti-cancer messengers: A systematic analysis from a single-cell perspective

Breast cancer is the most common malignancy in the world and one of the leading causes of cancer death, which is a heterogeneous disease involving genetic and environmental factors. Breast cancer stem cells (BCSCs) are the main players in the aggressiveness of different tumors, at the same time, these cells are the main challenge for cancer treatment. There are multiple treatment options for breast cancer (BC) patients and the lack of understanding of prognostic and predictive biomarkers for breast cancer is a potential research direction for us to develop better treatments in the future. In this paper, we conducted a correlation analysis between SIK2 and clinical traits by searching numerous BRCA datasets in the GEO database. The model was constructed and validated by incorporating tumor samples from the TCGA-BRCA cohort. Surprisingly, we found differential expression of SIK2 gene in individual tumor samples from the UCSC database. Subsequently, we found significantly high expression of SIK2 in epithelial cells by comparing the differential expression of SIK2 in different cell subpopulations and performed subsequent immune infiltration and pathway correlation analysis. Differential genes in SIK2+ epithelial cells, which may be potential therapeutic targets for breast cancer. In conclusion, our results suggest that SIK2 may be a potential prognostic and predictive biomarker that could serve as an oncogenic messenger for breast cancer. This discovery of SIK2 may provide more valuable references for potential therapeutic tools for breast cancer.

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.

Discovery of pyrimidine-5-carboxamide derivatives as novel salt-inducible kinases (SIKs) inhibitors for inflammatory bowel disease (IBD) treatment

Salt-inducible kinases (SIKs) play a crucial role in inflammation process, acting as molecular switches that regulate the transformation of M1/M2 macrophages. HG-9-91-01 is a SIKs inhibitor with potent inhibitory activity against SIKs in the nanomolar range. However, its poor drug-like properties, including a rapid elimination rate, low in vivo exposure and high plasma protein binding rate, have hindered further research and clinical application. To improve the drug-like properties of HG-9-91-01, a series of pyrimidine-5-carboxamide derivatives were designed and synthesized through a molecular hybridization strategy. The most promising compound 8h was obtained with favorable activity and selectivity on SIK1/2, excellent metabolic stability in human liver microsome, enhanced in vivo exposure and suitable plasma protein binding rate. Mechanism research showed that compound 8h significantly up-regulated the expression of anti-inflammatory cytokine IL-10 and reduced the expression of pro-inflammatory cytokine IL-12 in bone marrow-derived macrophages. Furthermore, it significantly elevated expression of cAMP response element-binding protein (CREB) target genes IL-10, c-FOS and Nurr77. Compound 8h also induced the translocation of CREB-regulated transcriptional coactivator 3 (CRTC3) and elevated the expression of LIGHT, SPHK1 and Arginase 1. Additionally, compound 8h demonstrated excellent anti-inflammatory effects in a DSS-induced colitis model. Generally, this research indicated that compound 8h has the potential to be developed as an anti-inflammatory drug candidate.

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.

Molecular Docking, Molecular Dynamics Simulations, and Free Energy Calculation Insights into the Binding Mechanism between VS-4718 and Focal Adhesion Kinase

Focal adhesion kinase (FAK) is a 125 kDa nonreceptor tyrosine kinase that plays an important role in many carcinomas. Thus, the targeting of FAK by small molecules is considered to be promising for cancer therapy. Some FAK inhibitors have been reported as potential anticancer drugs and have entered into clinical development; for example, VS-4718 is currently undergoing clinical trials. However, the lack of crystal structural data for the binding of VS-4718 with FAK has hindered the optimization of this anticancer agent. In this work, the VS-4718/FAK interaction model was obtained by molecular docking and molecular dynamics simulations. The binding free energies of VS-4718/FAK were also calculated using the molecular mechanics generalized Born surface area method. It was found that the aminopyrimidine group formed hydrogen bonds with the C502 residue of the hinge loop, while the D564 residue of the T-loop interacted with the amide group. In addition, I428, A452, V484, M499, G505, and L553 residues formed hydrophobic interactions with VS-4718. The obtained results therefore provide an improved understanding of the interaction between human FAK and VS-4718. Based on the obtained binding mechanism, 47 novel compounds were designed to target the adenosine 5'-triphosphate-binding pocket of human FAK, and ensemble docking was performed to assess the effects of these modifications on the inhibitor binding affinity. This work is also expected to provide additional insights into potential future target design strategies based on VS-4718.

Designing of potent anti-diabetic molecules by targeting SIK2 using computational approaches

Diabetes mellitus (DM) is one of the major health problems worldwide. WHO have estimated that 439 million people may have DM by the year 2030. Several classes of drugs such as sulfonylureas, meglitinides, thiazolidinediones etc. are available to manage this disease, however, there is no cure for this disease. Salt inducible kinase 2 (SIK2) is expressed several folds in adipose tissue than in normal tissues and thus SIK2 is one of the attractive targets for DM treatment. SIK2 inhibition improves glucose homeostasis. Several analogues have been reported and experimentally proven against SIK for DM treatment. But, identifying potential SIK2 inhibitors with improved efficacy and good pharmacokinetic profiles will be helpful for the effective treatment of DM. The objective of the present study is to identify selective SIK2 inhibitors with good pharmacokinetic profiles. Due to the unavailability of SIK2 structure, the modeled structure of SIK2 will be an important to understand the atomic level of SIK2 inhibitors in the binding site pocket. In this study, different molecular modeling studies such as Homology Modeling, Molecular Docking, Pharmacophore-based virtual screening, MD simulations, Density Functional Theory calculations and WaterMap analysis were performed to identify potential SIK2 inhibitors. Five molecules from different databases such as Binding_4067, TosLab_837067, NCI_349155, Life chemicals_ F2565-0113, Enamine_7623111186 molecules were identified as possible SIK2 inhibitors.

Molecular Dynamics Simulations of the Conformational Plasticity in the Active Pocket of Salt-Inducible Kinase 2 (SIK2) Multi-State Binding with Bosutinib

The kinase domain is highly conserved among protein kinases 'in terms of both sequence and structure. Conformational rearrangements of the kinase domain are affected by the phosphorylation of residues and the binding of kinase inhibitors. Interestingly, the conformational rearrangement of the active pocket plays an important role in kinase activity and can be used to design novel kinase inhibitors. We characterized the conformational plasticity of the active pocket when bosutinib was bound to salt-inducible kinase 2 (SIK2) using homology modeling and molecular dynamics simulations. Ten different initial complex models were constructed using the Morph server, ranging from open to closed conformations of SIK2 binding with bosutinib. Our simulation showed that bosutinib binds SIK2 with up or down conformations of the P-loop and with all the conformations of the activation loop. In addition, the αC-helix conformation was induced by the conformation of the activation loop, and the salt bridge formed only with its open conformation. The binding affinity of the models was also determined using the molecular mechanics generalized Born surface area method. Bosutinib was found to form a strong binding model with SIK2 and hydrophobic interactions were the dominant factor. This discovery may help guide the design of novel SIK2 inhibitors.

The binding mechanism of NHWD-870 to bromodomain-containing protein 4 based on molecular dynamics simulations and free energy calculation

Bromodomain and extra-terminal (BET) proteins (BRD2, BRD3, BRD4, and BRDT) are epigenetic readers with tandem bromodomains.