Deciphering the conformational transitions of LIMK2 active and inactive states to ponder specific druggable states through microsecond scale molecular dynamics simulation

LIMK2 inhibitors are one of the potential therapeutic modalities for treating various diseases. In the current scenario, there is a paucity of effective LIMK inhibitors that are highly specific with minimal off-target effects. To date, the conformational transitions of LIMK2 from DFG_inαC_in (CIDI) (active) to DFG_outαC_out (CODO) (inactive) states are yet to be probed and are essential for capturing the unique, druggable conformations. Therefore, this study was intended to capture the diverse conformational states of LIMK2 for accelerating the rational identification of conformation specific inhibitors through high-end structural bioinformatics protocols. Hence, in this study, molecular modelling followed by an extensive microsecond timescale of molecular dynamics simulation was performed encompassing perturbation response scanning, metapath, and community analysis towards the conformational sampling of LIMK2. Overall this study precisely identifies the conformational ensemble of LIMK2 the intermediate inactive states namely, CIDO, CinterDinter, CIDinter, CinterDI, CinterDO, CODI, CODinter apart from CIDI and CODO. This also facilitated observing that β8 preceding XDFG, αC (F373, L374), and αD (L413) as the major effectors that may facilitate the regulation of varying conformational transitions among the states. Additionally, the conserved β sheets and the loops namely, C.l, b.l, and G/P.loop were observed to be involved in the metapath for allosteric communication among the intermediates with CIDI and CODO state. Moreover, only the CODO state was observed to have closed type A.l, while the CIDI and other intermediate states except for CIDO were observed to have open-DFG out type A.l, thereby enabling the binding of substrate. Apart from these, the druggable site analysis inferred that the CIDI and CODO states harbor prominent druggable sites spanning the conserved N-lobe, while the intermediates were observed to have unraveled allosteric druggable sites distal from the ATP binding site, majorly spanning the C-lobe of LIMK2. Thus, this study provides potential insights into the intermediate conformational druggable states of LIMK2 and also the druggable conformations, especially the inactive states of LIMK2, as a specific therapeutic targeting mode. Thus, providing a widened avenue to ponder the allosteric sites or the isoform selectivity conformations for targeting LIMK2 in various disease conditions.

The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning

BACKGROUND

Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest.

METHOD

A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition.

RESULTS

The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space.

CONCLUSIONS

The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements.

GENERAL SIGNIFICANCE

The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.

Activation loop targeting strategy for design of receptor-interacting protein kinase 2 (RIPK2) inhibitors

Development of selective kinase inhibitors remains a challenge due to considerable amino acid sequence similarity among family members particularly in the ATP binding site. Targeting the activation loop might offer improved inhibitor selectivity since this region of kinases is less conserved. However, the strategy presents difficulties due to activation loop flexibility. Herein, we report the design of receptor-interacting protein kinase 2 (RIPK2) inhibitors based on pan-kinase inhibitor regorafenib that aim to engage basic activation loop residues Lys169 or Arg171. We report development of CSR35 that displayed >10-fold selective inhibition of RIPK2 versus VEGFR2, the target of regorafenib. A co-crystal structure of CSR35 with RIPK2 revealed a resolved activation loop with an ionic interaction between the carboxylic acid installed in the inhibitor and the side-chain of Lys169. Our data provides principle feasibility of developing activation loop targeting type II inhibitors as a complementary strategy for achieving improved selectivity.

Tyrosine 842 in the activation loop is required for full transformation by the oncogenic mutant FLT3-ITD

The type III receptor tyrosine kinase FLT3 is frequently mutated in acute myeloid leukemia. Oncogenic FLT3 mutants display constitutive activity leading to aberrant cell proliferation and survival. Phosphorylation on several critical tyrosine residues is known to be essential for FLT3 signaling. Among these tyrosine residues, Y842 is located in the so-called activation loop. The position of this tyrosine residue is well conserved in all receptor tyrosine kinases. It has been reported that phosphorylation of the activation loop tyrosine is critical for catalytic activity for some but not all receptor tyrosine kinases. The role of Y842 residue in FLT3 signaling has not yet been studied. In this report, we show that Y842 is not important for FLT3 activation or ubiquitination but plays a critical role in regulating signaling downstream of the receptor as well as controlling receptor stability. We found that mutation of Y842 in the FLT3-ITD oncogenic mutant background reduced cell viability and increased apoptosis. Furthermore, the introduction of the Y842 mutation in the FLT3-ITD background led to a dramatic reduction in in vitro colony forming capacity. Additionally, mice injected with cells expressing FLT3-ITD/Y842F displayed a significant delay in tumor formation, compared to FLT3-ITD expressing cells. Microarray analysis comparing gene expression regulated by FLT3-ITD versus FLT3-ITD/Y842F demonstrated that mutation of Y842 causes suppression of anti-apoptotic genes. Furthermore, we showed that cells expressing FLT3-ITD/Y842F display impaired activity of the RAS/ERK pathway due to reduced interaction between FLT3 and SHP2 leading to reduced SHP2 activation. Thus, we suggest that Y842 is critical for FLT3-mediated RAS/ERK signaling and cellular transformation.

Structural Studies of ERK2 Protein Complexes

ERK1 and ERK2 (ERK1/2) are the primary effector kinases of the RAS-RAF-MEK-ERK signaling pathway. A variety of substrates and regulatory partners associate with ERK1/2 through distinct D-peptide- and DEF-docking sites on their kinase domains. While understanding of D-peptides that bind to ERK1/2 has become increasingly clear over the last decade, only more recently have structures of proteins interacting with other binding sites on ERK1/2 become available. PEA-15 is a 130-residue ERK1/2 regulator that engages both the D-peptide- and DEF-docking sites of ERK kinases, and directly sequesters the ERK2 activation loop in various different phosphorylation states. Here we describe the methods used to derive crystallization-grade complexes of ERK2-PEA-15, which may also be adapted for other regulators that associate with the activation loop of ERK1/2.

Activation and inhibition of Snf1 kinase activity by phosphorylation within the activation loop

The AMP-activated protein kinase is a metabolic regulator that transduces information about energy and nutrient availability. In yeast, the AMP-activated protein kinase, called Snf1, is activated when energy and nutrients are scarce. Earlier studies have demonstrated that activation of Snf1 requires the phosphorylation of the activation loop on threonine 210. Here we examined the regulation of Snf1 kinase activity in response to phosphorylation at other sites. Phosphoproteomic studies have identified numerous phosphorylation sites within the Snf1 kinase enzyme. We made amino acid substitutions in the Snf1 protein that were either non-phosphorylatable (serine to alanine) or phospho-mimetic (serine to glutamate) and examined the effects of these changes on Snf1 kinase function in vivo and on its catalytic activity in vitro. We found that changes to most of the phosphorylation sites had no effect on Snf1 kinase function. However, changes to serine 214, a site within the kinase activation loop, inhibited Snf1 kinase activity. Snf1-activating kinase 1 still phosphorylates Snf1-S214E on threonine 210 but the S214E enzyme is non-functional in vivo and catalytically inactive in vitro. We conclude that yeast have developed two distinct pathways for down-regulating Snf1 activity. The first is through direct dephosphorylation of the conserved activation loop threonine. The second is through phosphorylation of serine 214.

Differential regulation of FGFR3 by PTPN1 and PTPN2

Aberrant expression and activation of FGFR3 is associated with disease states including bone dysplasia and malignancies of bladder, cervix, and bone marrow. MS analysis of protein-phosphotyrosine in multiple myeloma cells revealed a prevalent phosphorylated motif, D/EYYR/K, derived from the kinase domain activation loops of tyrosine kinases including FGFR3 corresponding to a recognition sequence of protein-tyrosine phosphatase PTPN1. Knockdown of PTPN1 or the related enzyme PTPN2 by RNAi resulted in ligand-independent activation of FGFR3. Modulation of FGFR3 activation loop phosphorylation by both PTPN1 and PTPN2 was a function of receptor trafficking and phosphotyrosine phosphatase (PTP) compartmentalization. The FGFR3 activation loop motif DYYKK(650) is altered to DYYKE(650) in the oncogenic variant FGFR3(K650E) , and consequently it is constitutively fully activated and unaffected by activation loop phosphorylation. FGFR3(K650E) was nevertheless remarkably sensitive to negative regulation by PTPN1 and PTPN2. This suggests that in addition to modulating FGFR3 phosphorylation, PTPN1 and PTPN2 constrain the kinase domain by fostering an inactive-state. Loss of this constraint in response to ligand or impaired PTPN1/N2 may initiate FGFR3 activation. These results suggest a model wherein PTP expression levels may define conditions that select for ectopic FGFR3 expression and activation during tumorigenesis.