Integration of a Randomized Sequence Scanning Approach in AlphaFold2 and Local Frustration Profiling of Conformational States Enable Interpretable Atomistic Characterization of Conformational Ensembles and Detection of Hidden Allosteric States in the ABL1 Protein Kinase.

Despite the success of AlphaFold methods in predicting single protein structures, these methods showed intrinsic limitations in the characterization of multiple functional conformations of allosteric proteins. The recent NMR-based structural determination of the unbound ABL kinase in the active state and discovery of the inactive low-populated functional conformations that are unique for ABL kinase present an ideal challenge for the AlphaFold2 approaches. In the current study, we employ several adaptations of the AlphaFold2 methodology to predict protein conformational ensembles and allosteric states of the ABL kinase including randomized alanine sequence scanning combined with the multiple sequence alignment subsampling proposed in this study. We show that the proposed new AlphaFold2 adaptation combined with local frustration profiling of conformational states enables accurate prediction of the protein kinase structures and conformational ensembles, also offering a robust approach for interpretable characterization of the AlphaFold2 predictions and detection of hidden allosteric states. We found that the large high frustration residue clusters are uniquely characteristic of the low-populated, fully inactive ABL form and can define energetically frustrated cracking sites of conformational transitions, presenting difficult targets for AlphaFold2. The results of this study uncovered previously unappreciated fundamental connections between local frustration profiles of the functional allosteric states and the ability of AlphaFold2 methods to predict protein structural ensembles of the active and inactive states. This study showed that integration of the randomized sequence scanning adaptation of AlphaFold2 with a robust landscape-based analysis allows for interpretable atomistic predictions and characterization of protein conformational ensembles, providing a physical basis for the successes and limitations of current AlphaFold2 methods in detecting functional allosteric states that play a significant role in protein kinase regulation.

The molecular basis of Abelson kinase regulation by its αI-helix.

Abelson tyrosine kinase (Abl) is regulated by the arrangement of its regulatory core, consisting sequentially of the SH3, SH2, and kinase (KD) domains, where an assembled or disassembled core corresponds to low or high kinase activity, respectively. It was recently established that binding of type II ATP site inhibitors, such as imatinib, generates a force from the KD N-lobe onto the SH3 domain and in consequence disassembles the core. Here, we demonstrate that the C-terminal αI-helix exerts an additional force toward the SH2 domain, which correlates both with kinase activity and type II inhibitor-induced disassembly. The αI-helix mutation E528K, which is responsible for the ABL1 malformation syndrome, strongly activates Abl by breaking a salt bridge with the KD C-lobe and thereby increasing the force onto the SH2 domain. In contrast, the allosteric inhibitor asciminib strongly reduces Abl's activity by fixating the αI-helix and reducing the force onto the SH2 domain. These observations are explained by a simple mechanical model of Abl activation involving forces from the KD N-lobe and the αI-helix onto the KD/SH2SH3 interface.

Myristoyl's dual role in allosterically regulating and localizing Abl kinase.

c-Abl kinase, a key signaling hub in many biological processes ranging from cell development to proliferation, is tightly regulated by two inhibitory Src homology domains. An N-terminal myristoyl modification can bind to a hydrophobic pocket in the kinase C-lobe, which stabilizes the autoinhibitory assembly. Activation is triggered by myristoyl release. We used molecular dynamics simulations to show how both myristoyl and the Src homology domains are required to impose the full inhibitory effect on the kinase domain and reveal the allosteric transmission pathway at residue-level resolution. Importantly, we find myristoyl insertion into a membrane to thermodynamically compete with binding to c-Abl. Myristoyl thus not only localizes the protein to the cellular membrane, but membrane attachment at the same time enhances activation of c-Abl by stabilizing its preactivated state. Our data put forward a model in which lipidation tightly couples kinase localization and regulation, a scheme that currently appears to be unique for this non-receptor tyrosine kinase.

A biophysical framework for double-drugging kinases.

Selective orthosteric inhibition of kinases has been challenging due to the conserved active site architecture of kinases and emergence of resistance mutants. Simultaneous inhibition of distant orthosteric and allosteric sites, which we refer to as "double-drugging", has recently been shown to be effective in overcoming drug resistance. However, detailed biophysical characterization of the cooperative nature between orthosteric and allosteric modulators has not been undertaken. Here, we provide a quantitative framework for double-drugging of kinases employing isothermal titration calorimetry, Förster resonance energy transfer, coupled-enzyme assays, and X-ray crystallography. We discern positive and negative cooperativity for Aurora A kinase (AurA) and Abelson kinase (Abl) with different combinations of orthosteric and allosteric modulators. We find that a conformational equilibrium shift is the main principle governing cooperativity. Notably, for both kinases, we find a synergistic decrease of the required orthosteric and allosteric drug dosages when used in combination to inhibit kinase activities to clinically relevant inhibition levels. X-ray crystal structures of the double-drugged kinase complexes reveal the molecular principles underlying the cooperative nature of double-drugging AurA and Abl with orthosteric and allosteric inhibitors. Finally, we observe a fully closed conformation of Abl when bound to a pair of positively cooperative orthosteric and allosteric modulators, shedding light on the puzzling abnormality of previously solved closed Abl structures. Collectively, our data provide mechanistic and structural insights into rational design and evaluation of double-drugging strategies.

Correspondence on "Synergy and Antagonism between Allosteric and Active-Site Inhibitors of Abl Tyrosine Kinase".

Soellner published on the interplay between allosteric and adenosine triphosphate (ATP)-competitive inhibitors of ABL kinase, showing that the latter preferably binds to different conformational states of ABL compared to allosteric agents that specifically target the ABL myristate pocket (STAMP) and deducing that asciminib cannot bind to ABL simultaneously with ATP-competitive drugs. These results are to some extent in line with ours, although our analyses of dose-response matrices from combinations of asciminib with imatinib, nilotinib or dasatinib, show neither synergy nor antagonism, but suggest additive antiproliferative effects on BCR-ABL-dependent KCL22 cells. Furthermore, our X-ray crystallographic, solution nuclear magnetic resonance (NMR), and isothermal titration calorimetry studies show that asciminib can bind ABL concomitantly with type-1 or -2 ATP-competitive inhibitors to form ternary complexes. Concomitant binding of asciminib with imatinib, nilotinib, or dasatinib might translate to benefit some chronic myeloid leukaemia patients.

Betulin, a Newly Characterized Compound in Acacia auriculiformis Bark, Is a Multi-Target Protein Kinase Inhibitor.

The purpose of this work is to investigate the protein kinase inhibitory activity of constituents from Acacia auriculiformis stem bark. Column chromatography and NMR spectroscopy were used to purify and characterize betulin from an ethyl acetate soluble fraction of acacia bark. Betulin, a known inducer of apoptosis, was screened against a panel of 16 disease-related protein kinases. Betulin was shown to inhibit Abelson murine leukemia viral oncogene homolog 1 (ABL1) kinase, casein kinase 1ε (CK1ε), glycogen synthase kinase 3α/ß (GSK-3 α/ß), Janus kinase 3 (JAK3), NIMA Related Kinase 6 (NEK6), and vascular endothelial growth factor receptor 2 kinase (VEGFR2) with activities in the micromolar range for each. The effect of betulin on the cell viability of doxorubicin-resistant K562R chronic myelogenous leukemia cells was then verified to investigate its putative use as an anti-cancer compound. Betulin was shown to modulate the mitogen-activated protein (MAP) kinase pathway, with activity similar to that of imatinib mesylate, a known ABL1 kinase inhibitor. The interaction of betulin and ABL1 was studied by molecular docking, revealing an interaction of the inhibitor with the ABL1 ATP binding pocket. Together, these data demonstrate that betulin is a multi-target inhibitor of protein kinases, an activity that can contribute to the anticancer properties of the natural compound and to potential treatments for leukemia.

Molecular insights on ABL kinase activation using tree-based machine learning models and molecular docking.

Abelson kinase (c-Abl) is a non-receptor tyrosine kinase involved in several biological processes essential for cell differentiation, migration, proliferation, and survival. This enzyme's activation might be an alternative strategy for treating diseases such as neutropenia induced by chemotherapy, prostate, and breast cancer. Recently, a series of compounds that promote the activation of c-Abl has been identified, opening a promising ground for c-Abl drug development. Structure-based drug design (SBDD) and ligand-based drug design (LBDD) methodologies have significantly impacted recent drug development initiatives. Here, we combined SBDD and LBDD approaches to characterize critical chemical properties and interactions of identified c-Abl's activators. We used molecular docking simulations combined with tree-based machine learning models-decision tree, AdaBoost, and random forest to understand the c-Abl activators' structural features required for binding to myristoyl pocket, and consequently, to promote enzyme and cellular activation. We obtained predictive and robust models with Matthews correlation coefficient values higher than 0.4 for all endpoints and identified characteristics that led to constructing a structure-activity relationship model (SAR).

Role of the ABL tyrosine kinases in the epithelial-mesenchymal transition and the metastatic cascade.

The ABL kinases, ABL1 and ABL2, promote tumor progression and metastasis in various solid tumors. Recent reports have shown that ABL kinases have increased expression and/or activity in solid tumors and that ABL inactivation impairs metastasis. The therapeutic effects of ABL inactivation are due in part to ABL-dependent regulation of diverse cellular processes related to the epithelial to mesenchymal transition and subsequent steps in the metastatic cascade. ABL kinases target multiple signaling pathways required for promoting one or more steps in the metastatic cascade. These findings highlight the potential utility of specific ABL kinase inhibitors as a novel treatment paradigm for patients with advanced metastatic disease. Video abstract.

Design, synthesis, and anticancer activity of novel 4-thiazolidinone-phenylaminopyrimidine hybrids.

4-Thiazolidinones and phenylaminopyrimidines are known as anticancer agents. Imatinib is the pioneer phenylaminopyrimidine derivative kinase inhibitor, which is used for the treatment of chronic myeloid leukemia. With a hybrid approach, a novel series of 5-benzylidene-2-arylimino-4-thiazolidinone derivatives containing phenylaminopyrimidine core were designed, synthesized, and tested for their anticancer activity on K562 (chronic myeloid leukemia), PC3 (prostat cancer), and SHSY-5Y (neuroblastoma) cells. Since superior anticancer activity was observed on K562 cells, further biological studies of selected compounds (8, 15, and 34) were performed on K562 cells. For the synthesis of designed compounds, thiourea compounds were converted to 2-imino-1,3-thiazolidin-4-ones with α-chloroacetic acid in the presence of sodium acetate. 5-Benzylidene-2-imino-1,3-thiazolidin-4-one derivatives were obtained by Knoevenagel condensation of 2-imino-1,3-thiazolidin-4-ones with related aldehydes. Compounds 8, 15, and 34 were evaluated for cell viability, apoptosis studies, cell cycle experiments, and DNA damage assays. IC50 values of compounds 8, 15, and 34 were found as 5.26 ± 1.03, 3.52 ± 0.91, and 8.16 ± 1.27 µM, respectively, in K562 cells. Preferably, these compounds showed less toxicity towards L929 cells compared to imatinib. Furthermore, compounds 8 and 15 significantly induced early and late apoptosis in a time-dependent manner. Compounds 15 and 34 induced cell cycle arrest at G0/G1 phase and compound 8 caused cell cycle arrest at G2/M phase. Based on DNA damage assay, compounds 8 and 15 were found to be more genotoxic than imatinib towards K562 cells. To put more molecular insight, possible Abl inhibition mechanisms of most active compounds were predicted by molecular docking studies. In conclusion, a novel series of 5-benzylidene-2-arylimino-4-thiazolidinone derivatives and their promising anticancer activities were reported herein.

Pathogenic variants causing ABL1 malformation syndrome cluster in a myristoyl-binding pocket and increase tyrosine kinase activity.

ABL1 is a proto-oncogene encoding a nonreceptor tyrosine kinase, best known in the somatic BCR-ABL fusion gene associated with chronic myeloid leukaemia. Recently, germline missense variants in ABL1 have been found to cause an autosomal dominant developmental syndrome with congenital heart disease, skeletal malformations and characteristic facies. Here, we describe a series of six new unrelated individuals with heterozygous missense variants in ABL1 (including four novel variants) identified via whole exome sequencing. All the affected individuals in this series recapitulate the phenotype of the ABL1 developmental syndrome and additionally we affirm that hearing impairment is a common feature of the condition. Four of the variants cluster in the myristoyl-binding pocket of ABL1, a region critical for auto-inhibitory regulation of the kinase domain. Bio-informatic analysis of transcript-wide conservation and germline/somatic variation reveals that this pocket region is subject to high missense constraint and evolutionary conservation. Functional work to investigate ABL1 kinase activity in vitro by transient transfection of HEK293T cells with variant ABL1 plasmid constructs revealed increased phosphorylation of ABL1-specific substrates compared to wild-type. The increased tyrosine kinase activity was suppressed by imatinib treatment. This case series of six new patients with germline heterozygous ABL1 missense variants further delineates the phenotypic spectrum of this condition and recognises microcephaly as a common finding. Our analysis supports an ABL1 gain-of-function mechanism due to loss of auto-inhibition, and demonstrates the potential for pharmacological inhibition using imatinib.

Multiplexable fluorescence lifetime imaging (FLIM) probes for Abl and Src-family kinases.

Many commonly employed strategies to map kinase activities in live cells require expression of genetically encoded proteins (e.g. FRET sensors). In this work, we describe the development and preliminary application of a set of cell-penetrating, fluorophore labelled peptide substrates for fluorescence lifetime imaging (FLIM) of Abl and Src-family kinase activities. These probes do not rely on FRET pairs or genetically-encoded protein expression. We further demonstrate probe multiplexing and pixel-by-pixel quantification to estimate the relative proportion of modified probe, suggesting that this strategy will be useful for detailed mapping of single cell and subcellular dynamics of multiple kinases concurrently in live cells.

c­Abl­mediated tyrosine phosphorylation of DNA damage response proteins and implications in important cellular functions (Review).

Tyrosine phosphorylation is an essential post­translational protein modification catalyzed by tyrosine kinases. c­Abl is a crucial non­receptor tyrosine kinase, which is most commonly activated by auto­phosphorylation, DNA damage and by interacting with other protein kinases. DNA damage response (DDR) proteins stimulated by DNA lesions or chromatin alterations recruit the DNA repair and cell cycle checkpoint machinery to restore genome integrity and cellular homeostasis. The fundamental roles of activated c­Abl tyrosine kinase in cellular response pathways have been intensively and extensively investigated and in recent years, a number of c­Abl protein binding partners have been determined; however, the functional roles of these molecules remain to be determined. The present review aimed to summarize the DDR proteins phosphorylated by c­Abl tyrosine kinase that have been identified to date, in addition to the functional outcomes of these phosphotyrosine events. Notably, it has been discovered that c­Abl tyrosine kinase can bind with and phosphorylate DDR proteins at different tyrosine sites, which serve distinct roles in various cellular contexts.

Ultra-accurate Duplex Sequencing for the assessment of pretreatment ABL1 kinase domain mutations in Ph+ ALL.

Mutations of ABL1 are the dominant mechanism of relapse in Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph + ALL). We performed highly accurate Duplex Sequencing of exons 4-10 of ABL1 on bone marrow or peripheral blood samples from 63 adult patients with previously untreated Ph + ALL who received induction with intensive chemotherapy plus a BCR-ABL1 TKI. We identified ABL1 mutations prior to BCR-ABL1 TKI exposure in 78% of patients. However, these mutations were generally present at extremely low levels (median variant allelic frequency 0.008% [range, 0.004%-3.71%] and did not clonally expand and lead to relapse in any patient, even when the pretreatment mutation was known to confer resistance to the TKI received. In relapse samples harboring a TKI-resistant ABL1 mutation, the corresponding mutation could not be detected pretreatment, despite validated sequencing sensitivity of Duplex Sequencing down to 0.005%. In samples under the selective pressure of ongoing TKI therapy, we detected low-level, emerging resistance mutations up to 5 months prior to relapse. These findings suggest that pretreatment ABL1 mutation assessment should not guide upfront TKI selection in Ph + ALL, although serial testing while on TKI therapy may allow for early detection of clinically actionable resistant clones.

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.

Design, Synthesis and Biological Evaluation of Pentacyclic Triterpene Derivatives: Optimization of Anti-ABL Kinase Activity.

Imatinib, an Abelson (ABL) tyrosine kinase inhibitor, is a lead molecular-targeted drug against chronic myelogenous leukemia (CML). To overcome its resistance and adverse effects, new inhibitors of ABL kinase are needed. Our previous study showed that the benzyl ester of gypsogenin (1c), a pentacyclic triterpene, has anti-ABL kinase and a subsequent anti-CML activity. To optimize its activities, benzyl esters of carefully selected triterpenes (PT1-PT6), from different classes comprising oleanane, ursane and lupane, and new substituted benzyl esters of gypsogenin (GP1-GP5) were synthesized. All of the synthesized compounds were purified and charachterized by different spectroscopic methods. Cytotoxicity of the parent triterpenes and the synthesized compounds against CML cell line K562 was examined; revealing three promising compounds PT5, GP2 and GP5 (IC50 5.46, 4.78 and 3.19 µM, respectively). These compounds were shown to inhibit extracellular signal-regulated kinase (ERK) downstream signaling, and induce apoptosis in K562 cells. Among them, PT5 was identified to have in vitro activity (IC50 = 1.44 µM) against ABL1 kinase, about sixfold of 1c, which was justified by molecular docking. The in vitro activities of GP2 and GP5 are less than PT5, hence they were supposed to possess other more mechanisms of cytotoxicity. In general, our design and derivatizations resulted in enhancing the activity against ABL1 kinase and CML cells.

Mapping RNAPII CTD Phosphorylation Reveals That the Identity and Modification of Seventh Heptad Residues Direct Tyr1 Phosphorylation.

The C-terminal domain (CTD) of the largest subunit in eukaryotic RNA polymerase II has a repetitive heptad sequence of Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7 which is responsible for recruiting transcriptional regulatory factors. The seventh heptad residues in mammals are less conserved and subject to various post-translational modifications, but the consequences of such variations are not well understood. In this study, we use ultraviolet photodissociation mass spectrometry, kinetic assays, and structural analyses to dissect how different residues or modifications at the seventh heptad position alter Tyr1 phosphorylation. We found that negatively charged residues in this position promote phosphorylation of adjacent Tyr1 sites, whereas positively charged residues discriminate against it. Modifications that alter the charges on seventh heptad residues such as arginine citrullination negate such distinctions. Such specificity can be explained by conserved, positively charged pockets near the active sites of ABL1 and its homologues. Our results reveal a novel mechanism for variations or modifications in the seventh heptad position directing subsequent phosphorylation of other CTD sites, which can contribute to the formation of various modification combinations that likely impact transcriptional regulation.

Rapid Screen for Tyrosine Kinase Inhibitor Resistance Mutations and Substrate Specificity.

We present a rapid and high-throughput yeast and flow cytometry based method for predicting kinase inhibitor resistance mutations and determining kinase peptide substrate specificity. Despite the widespread success of targeted kinase inhibitors as cancer therapeutics, resistance mutations arising within the kinase domain of an oncogenic target present a major impediment to sustained treatment efficacy. Our method, which is based on the previously reported YESS system, recapitulated all validated BCR-ABL1 mutations leading to clinical resistance to the second-generation inhibitor dasatinib, in addition to identifying numerous other mutations which have been previously observed in patients, but not yet validated as drivers of resistance. Further, we were able to demonstrate that the newer inhibitor ponatinib is effective against the majority of known single resistance mutations, but ineffective at inhibiting many compound mutants. These results are consistent with preliminary clinical and in vitro reports, indicating that mutations providing resistance to ponatinib are significantly less common; therefore, predicting ponatinib will be less susceptible to clinical resistance relative to dasatinib. Using the same yeast-based method, but with random substrate libraries, we were able to identify consensus peptide substrate preferences for the SRC and LYN kinases. ABL1 lacked an obvious consensus sequence, so a machine learning algorithm utilizing amino acid covariances was developed which accurately predicts ABL1 kinase peptide substrates.

Anticancer properties of N-alkyl-2, 4-diphenylimidazo [1, 2-a] quinoxalin-1-amine derivatives; kinase inhibitors.

Structure activity correlation revealed that the quinoxaline ring is a satisfactory backbone for anticancer activity and a specific functional group at position 1 and 2 can improve the activity. In this basis, besides quinoxaline, imidazoles as potential anticancer agents were used as a supplementary agents for cancer treatment. In this paper, a new series of N-alkyl-2, 4-diphenylimidazo [1, 2-a] quinoxalin-1-amine derivatives were synthesized in a simple and efficient step. The products are fully characterized by 1H NMR, 13C NMR, FT-IR, HRMS, and CHN elemental analysis. Several starting materials with different functionalities have been used for the synthesis of the final products with high isolated yields. The biological activities of the synthesized compounds were evaluated in kinase inhibition and cytotoxic activity in several cancerous cell lines. All compounds (6) were evaluated for inhibition of the cell proliferation using 4 cancerous cell lines. Five of the more active compounds were studied for determination of IC50%. Compounds 6(32-34) showed good activity on some of cancerous cell lines. The results showed that compound 6-32 has the highest biological activity (IC50% 9.77 for K562 cell line). An IC50% value of 15.84 µM was observed for 6-34. Furthermore 6-34 exhibited inhibition of ABL1 and c-Src kinases with an IC50% value of 5.25 µM and 3.94 µM respectively. Docking simulation was performed to position active synthesized compounds 6-32, 6-33, and 6-34 over the ABL1 active site in two different wild-type (DFG-in and DFG-out motif conformer) and T315I mutant to determine the probable binding orientation, conformation and mode of interaction. According to docking study, the docked location in wild type forms is similar and can be found near the P-loop region while in the case of T315I mutant form, the compounds have a distinct docked location which is close to the αC helix and activation loop. Also, it concluded the role of R1 substituent on phenyl ring produced higher interaction energy. Additionally, the detailed inter-molecular energy and types of non-bonding interaction of these compounds over the wild-type and mutant form of ABL1.

New Heteroleptic Ruthenium(II) Complexes with Sulfamethoxypyridazine and Diimines as Potential Antitumor Agents.

Two new complexes of Ru(II) with mixed ligands were prepared: [Ru(bpy)2smp](PF6) (1) and [Ru(phen)2smp](PF6) (2), in which smp = sulfamethoxypyridazine; bpy = 2,2'-bipyridine; phen = 1,10-phenanthroline. The complexes have been characterized by elemental and conductivity analyses; infrared, NMR, and electrospray ionization mass spectroscopies; and X-ray diffraction of single crystal. Structural analyses reveal a distorted octahedral geometry around Ru(II) that is bound to two bpy (in 1) or two phen (in 2) via their two heterocyclic nitrogens and to two nitrogen atoms from sulfamethoxypyridazine-one of the methoxypyridazine ring and the sulfonamidic nitrogen, which is deprotonated. Both complexes inhibit the growth of chronic myelogenous leukemia cells. The interaction of the complexes with bovine serum albumin and DNA is described. DNA footprinting using an oligonucleotide as substrate showed the complexes' preference for thymine base rich sites. It is worth notifying that the complexes interact with the Src homology SH3 domain of the Abl tyrosine kinase protein. Abl protein is involved in signal transduction and implicated in the development of chronic myelogenous leukemia. Nuclear magnetic resonance (NMR) studies of the interaction of complex 2 with the Abl-SH3 domain showed that the most affected residues were T79, G97, W99, and Y115.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications.

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.