Phosphorylation-dependent activity-based conformational changes in P21-activated kinase family members and screening of novel ATP competitive inhibitors

Structural basis of constitutive c-Src kinase activity due to R175L and W118A mutations

Cellular Src (c-Src) belongs to a non-receptor membrane-associated tyrosine kinase family that plays essential roles in cellular processes. Growing evidence suggests that R175L and W118A mutations in SH2/SH3 domains of c-Src functionally inactivate these domains leading to constitutive activation of kinase domain (KD). Here we modeled c-Src^R175L, c-Src^W118A and c-Src^W118A+R175L structures by inducing phosphorylation at Y416 or Y527, respectively to characterize the comparative dynamics in the active versus inactive states through molecular dynamics simulation assay. We observed more conformational readjustments in c-Src^open than its close variants. In particular, C-terminal tail residues of c-Src^W118A-open and c-Src^{W118A+R175L-open} demonstrate significantly higher transitions. The cross-correlation analysis revealed an anticorrelation behavior in the motion of KD with respect to SH2, SH3 and the linker region of Src^{W118A+R175L-open}, while in c-Src^WT-open, SH2 and SH3 domains were anticorrelated, while KD and C-terminal tail motions were correlated. Due to these conformational differences, c-Src open forms exhibited lower interaction between pY527 and SH2 domain. Through detailed structural analysis, we observed a uniform myristate binding cavity in c-Src^WT-open, while the myristoyl pockets of mutant forms were deformed. We propose that constitutive activation of mutant Src forms may presumably be achieved by the prolonged membrane binding due to unusual conformations of C-terminal and myristoyl switch residues that may result in a higher dephosphorylation rate at pY527 in the myristoylated c-Src. Thus, our study establishes novel clues to decipher the constitutive activation status of c-Src in response to known mutations that may help in devising novel therapeutic strategies for cancer metastasis treatment.Communicated by Ramaswamy H. Sarma.

The Inhibitory Mechanism of 7H-Pyrrolo[2,3-d]pyrimidine Derivatives as Inhibitors of P21-Activated Kinase 4 through Molecular Dynamics Simulation

The overexpression of p21-activated kinase 4 (PAK4) is associated with a variety of cancers. In this paper, the binding modes and inhibitory mechanisms of four 7H-pyrrolo[2,3-d]pyrimidine competitive inhibitors of PAK4 were investigated at the molecular level, mainly using molecular dynamics simulations and binding free energy calculations. The results show that the inhibitors had strong interactions with the hinge region, the β-sheets, and the residues with charged side chains around the 4-substituent. The terminal amino group of the inhibitor 5n was different from the other three, which could cause the enhancement of hydrogen bonds or electrostatic interactions formed with the surrounding residues. Thus, inhibitor 5n had the strongest inhibition capacity. The different halogen atoms on the 2-substituents of the inhibitors 5h, 5g, and 5e caused differences in the positions of the 2-benzene rings and affected the interactions of the hinge region. It also affected to some extent the orientations of the 4-imino groups and consequently their affinities for the surrounding charged residues. The combined results lead to the weakest inhibitory capacity of inhibitor 5e.

SHP-1 tyrosine phosphatase binding to c-Src kinase phosphor-dependent conformations: A comparative structural framework

SHP-1 is a cytosolic tyrosine phosphatase that is primarily expressed in hematopoietic cells. It acts as a negative regulator of numerous signaling pathways and controls multiple cellular functions involved in cancer pathogenesis. This study describes the binding preferences of SHP-1 (pY536) to c-Srcopen (pY416) and c-Srcclose (pY527) through in silico approaches. Molecular dynamics simulation analysis revealed more conformational changes in c-Srcclose upon binding to SHP-1, as compared to its active/open conformation that is stabilized by the cooperative binding of the C-SH2 domain and C-terminal tail of SHP-1 to c-Src SH2 and KD. In contrast, c-Srcclose and SHP-1 interaction is mediated by PTP domain-specific WPD-loop (WPDXGXP) and Q-loop (QTXXQYXF) binding to c-Srcclose C-terminal tail residues. The dynamic correlation analysis demonstrated a positive correlation for SHP-1 PTP with KD, SH3, and the C-terminal tail of c-Srcclose. In the case of the c-Srcopen-SHP-1 complex, SH3 and SH2 domains of c-Srcopen were correlated to C-SH2 and the C-terminal tail of SHP-1. Our findings reveal that SHP1-dependent c-Src activation through dephosphorylation relies on the conformational shift in the inhibitory C-terminal tail that may ease the recruitment of the N-SH2 domain to phosphotyrosine residue, resulting in the relieving of the PTP domain. Collectively, this study delineates the intermolecular interaction paradigm and underlying conformational readjustments in SHP-1 due to binding with the c-Src active and inactive state. This study will largely help in devising novel therapeutic strategies for targeting cancer development.

WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches

p21-Activated kinase 1 (PAK1) is positioned at the nexus of several oncogenic signaling pathways. Currently, there are no approved inhibitors for disabling the transfer of phosphate in the active site directly, as they are limited by lower affinity, and poor kinase selectivity. In this work, a repurposing study utilizing FDA-approved drugs from the DrugBank database was pursued with an initial selection of 27 molecules out of ∼2162 drug molecules, based on their docking energies and molecular interaction patterns. From the molecules that were considered for WaterMap analysis, seven molecules, namely, Mitoxantrone, Labetalol, Acalabrutinib, Sacubitril, Flubendazole, Trazodone, and Niraparib, ascertained the ability to overlap with high-energy hydration sites. Considering many other displaced unfavorable water molecules, only Acalabrutinib, Flubendazole, and Trazodone molecules highlighted their prominence in terms of binding affinity gains through ΔΔG that ranges between 6.44 and 2.59 kcal/mol. Even if Mitoxantrone exhibited the highest docking score and greater interaction strength, it did not comply with the WaterMap and molecular dynamics simulation results. Moreover, detailed MD simulation trajectory analyses suggested that the drug molecules Flubendazole, Niraparib, and Acalabrutinib were highly stable, observed from their RMSD values and consistent interaction pattern with Glu315, Glu345, Leu347, and Asp407 including the hydrophobic interactions maintained in the three replicates. However, the drug molecule Trazodone displayed a loss of crucial interaction with Leu347, which was essential to inhibit the kinase activity of PAK1. The molecular orbital and electrostatic potential analyses elucidated the reactivity and strong complementarity potentials of the drug molecules in the binding pocket of PAK1. Therefore, the CADD-based reposition efforts, reported in this work, helped in the successful identification of new PAK1 inhibitors that requires further investigation by in vitro analysis.