Interconversion of inactive to active conformation of MARK2: Insights from molecular modeling and molecular dynamics simulation

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.

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.

MARK2 regulates chemotherapeutic responses through class IIa HDAC-YAP axis in pancreatic cancer

Despite paclitaxel's wide use in cancer treatment, patient response rate is still low and drug resistance is a major clinical obstacle. Through a Phos-tag-based kinome-wide screen, we identified MARK2 as a critical regulator for paclitaxel chemosensitivity in PDAC. We show that MARK2 is phosphorylated by CDK1 in response to antitubulin chemotherapeutics and in unperturbed mitosis. Phosphorylation is essential for MARK2 in regulating mitotic progression and paclitaxel cytotoxicity in PDAC cells. Mechanistically, our findings also suggest that MARK2 controls paclitaxel chemosensitivity by regulating class IIa HDACs. MARK2 directly phosphorylates HDAC4 specifically during antitubulin treatment. Phosphorylated HDAC4 promotes YAP activation and controls expression of YAP target genes induced by paclitaxel. Importantly, combination of HDAC inhibition and paclitaxel overcomes chemoresistance in organoid culture and preclinical PDAC animal models. The expression levels of MARK2, HDACs, and YAP are upregulated and positively correlated in PDAC patients. Inhibition of MARK2 or class IIa HDACs potentiates paclitaxel cytotoxicity by inducing mitotic abnormalities in PDAC cells. Together, our findings identify the MARK2-HDAC axis as a druggable target for overcoming chemoresistance in PDAC.

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.

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

Salt inducible kinase 2 (SIK2) is a calcium/calmodulin-dependent protein kinase-like kinase that is implicated in a variety of biological phenomena, including cellular metabolism, growth, and apoptosis. SIK2 is the key target for various cancers, including ovarian, breast, prostate, and lung cancers. Although potent inhibitors of SIK2 are being developed, their binding stability and functional role are not presently known. In this work, we studied the detailed interactions between SIK2 and four of its inhibitors, HG-9-91-01, KIN112, MRT67307, and MRT199665, using molecular docking, molecular dynamics simulation, binding free energy calculation, and interaction fingerprint analysis. Intermolecular interactions revealed that HG-9-91-01 and KIN112 have stronger interactions with SIK2 than those of MRT199665 and MRT67307. The key residues involved in binding with SIK2 are conserved among all four inhibitors. Our results explain the detailed interaction of SIK2 with its inhibitors at the molecular level, thus paving the way for the development of targeted efficient anti-cancer drugs.

Dasatinib-SIK2 Binding Elucidated by Homology Modeling, Molecular Docking, and Dynamics Simulations

Salt-inducible kinases (SIKs) are calcium/calmodulin-dependent protein kinase (CAMK)-like (CAMKL) family members implicated in insulin signal transduction, metabolic regulation, inflammatory response, and other processes. Here, we focused on SIK2, which is a target of the Food and Drug Administration (FDA)-approved pan inhibitor N-(2-chloro-6-methylphenyl)-2-(6-(4-(2-hydroxyethyl)piperazin-1-yl)-2-methylpyrimidin-4-ylamino)thiazole-5-carboxamide (dasatinib), and constructed four representative SIK2 structures by homology modeling. We investigated the interactions between dasatinib and SIK2 via molecular docking, molecular dynamics simulation, and binding free energy calculation and found that dasatinib showed strong binding affinity for SIK2. Binding free energy calculations suggested that the modification of various dasatinib regions may provide useful information for drug design and to guide the discovery of novel dasatinib-based SIK2 inhibitors.

MARK4 protein can explore the active-like conformations in its non-phosphorylated state

Microtubule affinity-regulating kinase 4 (MARK4) is a Ser/Thr protein kinase, best known for its role in phosphorylating microtubule associated proteins, causing their detachment from microtubules. In the current study, the non-phosphorylated conformation of the activation loop was modeled in a structure representing the enzymatically inactive form of this protein, and its dynamics were evaluated through a 100 ns initial all-atom simulation, which was prolonged by another 2 μs. Although the activation loop was folding on itself and was leaning toward ATP site in the initial modeled structure, soon after the initiating the simulation, this loop stretched away from the ATP binding site and stably settled in its new position for the rest of simulation time. A network of hydrogen bonds, mainly between the activation segment residues, αC-helix and the catalytic loop reinforced this conformation. Interestingly, several features of active kinase conformation such as formation of R-spine, Glu106-Lys88 salt-bridge, and DFG-In motif were observed during a considerable number of trajectory frames. However, they were not sustainably established during the simulation time, except for the DFG-In motif. Consequently, this study introduces a stable conformation of the non-phosphorylated form of MARK4 protein with a partially stretched activation loop conformation as well as partial formation of R-spine, closely resembling the active kinase.