Abstract
Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.
Highlights
Protein kinases act as dynamic molecular switches in cellular signaling and their functional activity is essential for the integrity and viability of signaling pathways involved in cell cycle control, organism development, and stress response [1,2,3,4,5,6,7,8,9,10,11,12]
Using Molecular dynamics (MD) simulations, we first characterized conformational dynamics of the ABL kinase complexes with the type 2 inhibitors (Nilotinib, Ponatinib) and type 1 inhibitor Bosutinib
We have introduced a computational framework for explaining binding preferences of the ABL kinase inhibitors by unraveling a relationship between energetic hot spots of ligand binding and organization of the residue interaction networks in various kinase states
Summary
Protein kinases act as dynamic molecular switches in cellular signaling and their functional activity is essential for the integrity and viability of signaling pathways involved in cell cycle control, organism development, and stress response [1,2,3,4,5,6,7,8,9,10,11,12]. The human protein kinases represent one of the largest protein families that orchestrate functional processes in cellular networks and comprise an important class of therapeutic targets, owing to the presence of a highly conserved ATP binding pocket that can be exploited by small molecule inhibitors [13,14,15,16,17]. A diverse repertoire of crystallographic conformations has indicated that molecular switching mechanism of protein kinases may not necessarily imply an on–off binary operation (from inactive to active), but could rather represent a continuous dynamic process in which kinases may adopt a wide spectrum of inactive-like and active-like conformations exhibiting a range of activity levels. Conformational transitions between kinase states are orchestrated by three conserved structural motifs in the catalytic domain: the αC-helix, the DFG-Asp motif (DFG-Asp in, active; DFG-Asp out, inactive), and the activation loop (A-loop open, active; A-loop closed, inactive). The conserved His-Arg-Asp (HRD) motif in the catalytic loop and the DFG motif are coupled with the αC-helix to form conserved intramolecular networks termed regulatory spine (R-spine) and catalytic spine (C-spine) whose assembly and stabilization are intimately linked with the conformational transformations and kinase activation [25,26]
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have