Search

Protein kinase activity forms the backbone of cellular information transfer, acting both individually and as part of a broader network, the kinome. Correspondingly, their central role in signaling implicates kinome dysfunction as a common driver of cancer, where numerous kinases have been identified as having a causal or modulating role in cancer development and progression. Driven by their importance, the development of therapies targeting kinases has rapidly grown, with over 70 kinase inhibitors approved for use in the clinic and over double this number currently in clinical trials. Given the growing importance of kinase-targeted therapies, linking the relationship between kinase inhibitor treatment and their effects on downstream cellular phenotype is of clear importance for understanding treatment mechanisms and streamlining compound screening in therapy development. In this work, we combine two large-scale kinome profiling data sets and use them to link inhibitor-kinome interactions with cell line treatment responses (AUC/IC50). We then built computational models on this data set that achieve a high degree of prediction accuracy (R2of 0.7 and RMSE of 0.9), and were able to identify a set of well-characterized and understudied kinases that significantly affect cell responses. Further, we validated these models experimentally by testing predicted effects in breast cancer cell lines, and extended the model scope by performing additional validation in patient-derived pancreatic cancer cell lines. Overall, these results demonstrate that broad quantification of kinome inhibition state is highly predictive of downstream cellular phenotypes.

Numerous aspects of cellular signaling are regulated by the kinome – the network of over 500 protein kinases that guides and modulates information transfer throughout the cell. The key role played by both individual kinases and assemblies of kinases organized into functional subnetworks leads to kinome dysregulation being a key driver of many diseases, particularly cancer. In the case of pancreatic ductal adenocarcinoma (PDAC), a variety of kinases and associated signaling pathways have been identified for their key role in the establishment of disease as well as its progression. However, the identification of additional relevant therapeutic targets has been slow and is further confounded by interactions between the tumor and the surrounding tumor microenvironment. Fundamentally, it is an open question as to the degree to which knowledge of the state of the kinome at the protein level is able to provide insight into the downstream phenotype of the cell.In this work, we attempt to link the state of the kinome, or kinotype, with cell viability in representative PDAC tumor and stroma cell lines. Through the application of both regression and classification models to independent kinome perturbation and kinase inhibitor cell screen data, we find that the inferred kinotype of a cell has a significant and predictive relationship with cell viability. While regression models perform poorly, we find that classification approaches are able to predict drug viability effects. We further find that models are able to identify a set of kinases whose behavior in response to perturbation drive the majority of viability responses in these cell lines. Using the models to predict new compounds with cell viability effects and not in the initial data set, we conducted a validation screen that confirmed the accuracy of the models. These results suggest that characterizing the state of the protein kinome provides significant opportunity for better understanding signaling behavior and downstream cell phenotypes, as well as providing insight into the broader design of potential therapy design for PDAC.

From DNA base pairs to drug-receptor binding, hydrogen (H-)bonding and aromaticity are common features of heterocycles. Herein, the interplay of these bonding aspects is explored; H-bond strength modulation due to enhancement or disruption of heterocycles' aromaticity is revealed by comparing homodimer H-bond energies of aromatic heterocycles with analogues that have the same H-bonding moieties but lack cyclic π-conjugation. NMR studies of dimerization in C6D6 find aromaticity-modulated H-bonding (AMHB) energy effects of ca. ±30%, depending on whether they enhance or weaken aromatic delocalization. The attendant ring current perturbations expected from such modulation are confirmed by chemical shift changes in both observed ring C-H and calculated nucleus-independent sites. In silico modeling confirms that AMHB effects outweigh those of hybridization or dipole-dipole interaction.