AVIA 3.0: interactive portal for genomic variant and sample level analysis.

SUMMARY: The Annotation, Visualization and Impact Analysis (AVIA) is a web application combining multiple features to annotate and visualize genomic variant data. Users can investigate functional significance of their genetic alterations across samples, genes and pathways. Version 3.0 of AVIA offers filtering options through interactive charts and by linking disease relevant data sources. Newly incorporated services include gene, variant and sample level reporting, literature and functional correlations among impacted genes, comparative analysis across samples and against data sources such as TCGA and ClinVar, and cohort building. Sample and data management is now feasible through the application, which allows greater flexibility with sharing, reannotating and organizing data. Most importantly, AVIA's utility stems from its convenience for allowing users to upload and explore results without any a priori knowledge or the need to install, update and maintain software or databases. Together, these enhancements strengthen AVIA as a comprehensive, user-driven variant analysis portal. AVAILABILITYAND IMPLEMENTATION: AVIA is accessible online at https://avia-abcc.ncifcrf.gov.

Mannose receptor (CD206) activation in tumor-associated macrophages enhances adaptive and innate antitumor immune responses.

Solid tumors elicit a detectable immune response including the infiltration of tumor-associated macrophages (TAMs). Unfortunately, this immune response is co-opted into contributing toward tumor growth instead of preventing its progression. We seek to reestablish an antitumor immune response by selectively targeting surface receptors and endogenous signaling processes of the macrophage subtypes driving cancer progression. RP-182 is a synthetic 10-mer amphipathic analog of host defense peptides that selectively induces a conformational switch of the mannose receptor CD206 expressed on TAMs displaying an M2-like phenotype. RP-182-mediated activation of this receptor in human and murine M2-like macrophages elicits a program of endocytosis, phagosome-lysosome formation, and autophagy and reprograms M2-like TAMs to an antitumor M1-like phenotype. In syngeneic and autochthonous murine cancer models, RP-182 suppressed tumor growth, extended survival, and was an effective combination partner with chemo- or immune checkpoint therapy. Antitumor activity of RP-182 was also observed in CD206high patient-derived xenotransplantation models. Mechanistically, via selective reduction of immunosuppressive M2-like TAMs, RP-182 improved adaptive and innate antitumor immune responses, including increased cancer cell phagocytosis by reprogrammed TAMs.

Associations of CDH1 germline variant location and cancer phenotype in families with hereditary diffuse gastric cancer (HDGC).

INTRODUCTION: Hereditary diffuse gastric cancer (HDGC) is a cancer syndrome associated with variants in E-cadherin (CDH1), diffuse gastric cancer and lobular breast cancer. There is considerable heterogeneity in its clinical manifestations. This study aimed to determine associations between CDH1 germline variant status and clinical phenotypes of HDGC. METHODS: One hundred and fifty-two HDGC families, including six previously unreported families, were identified. CDH1 gene-specific guidelines released by the Clinical Genome Resource (ClinGen) CDH1 Variant Curation Expert Panel were applied for pathogenicity classification of truncating, missense and splice site CDH1 germline variants. We evaluated ORs between location of truncating variants of CDH1 and incidence of colorectal cancer, breast cancer and cancer at young age (gastric cancer at <40 or breast cancer <50 years of age). RESULTS: Frequency of truncating germline CDH1 variants varied across functional domains of the E-cadherin receptor gene and was highest in linker (0.05785 counts/base pair; p=0.0111) and PRE regions (0.10000; p=0.0059). Families with truncating CDH1 germline variants located in the PRE-PRO region were six times more likely to have family members affected by colorectal cancer (OR 6.20, 95% CI 1.79 to 21.48; p=0.004) compared with germline variants in other regions. Variants in the intracellular E-cadherin region were protective for cancer at young age (OR 0.2, 95% CI 0.06 to 0.64; p=0.0071) and in the linker regions for breast cancer (OR 0.35, 95% CI 0.12 to 0.99; p=0.0493). Different CDH1 genotypes were associated with different intracellular signalling activation levels including different p-ERK, p-mTOR and ß-catenin levels in early submucosal T1a lesions of HDGC families with different CDH1 variants. CONCLUSION: Type and location of CDH1 germline variants may help to identify families at increased risk for concomitant cancers that might benefit from individualised surveillance and intervention strategies.

Ca2+-Sensor Neurocalcin δ and Hormone ANF Modulate ANF-RGC Activity by Diverse Pathways: Role of the Signaling Helix Domain.

Prototype member of the membrane guanylate cyclase family, ANF-RGC (Atrial Natriuretic Factor Receptor Guanylate Cyclase), is the physiological signal transducer of two most hypotensive hormones ANF and BNP, and of the intracellular free Ca2+. Both the hormonal and the Ca2+-modulated signals operate through a common second messenger, cyclic GMP; yet, their operational modes are divergent. The hormonal pathways originate at the extracellular domain of the guanylate cyclase; and through a cascade of structural changes in its successive domains activate the C-terminal catalytic domain (CCD). In contrast, the Ca2+ signal operating via its sensor, myristoylated neurocalcin δ both originates and is translated directly at the CCD. Through a detailed sequential deletion and expression analyses, the present study examines the role of the signaling helix domain (SHD) in these two transduction pathways. SHD is a conserved 35-amino acid helical region of the guanylate cyclase, composed of five heptads, each meant to tune and transmit the hormonal signals to the CCD for their translation and generation of cyclic GMP. Its structure is homo-dimeric and the molecular docking analyses point out to the possibility of antiparallel arrangement of the helices. Contrary to the hormonal signaling, SHD has no role in regulation of the Ca2+- modulated pathway. The findings establish and define in molecular terms the presence of two distinct non-overlapping transduction modes of ANF-RGC, and for the first time demonstrate how differently they operate, and, yet generate cyclic GMP utilizing common CCD machinery.

Plasma metabolomic profiling of a ketamine and placebo crossover trial of major depressive disorder and healthy control subjects.

(R,S)-Ketamine produces rapid, robust, and sustained antidepressant effects in major depressive disorder. Specifically, its pharmacological efficacy in treatment refractory depression is considered a major breakthrough in the field. However, the mechanism of action of ketamine's rapid effect remains to be determined. In order to identify pathways that are responsible for ketamine's effect, a targeted metabolomic approach was carried out using a double-blind, placebo-controlled crossover design, with infusion order randomized with medication-free patients with treatment-resistant major depressive disorder (29 subjects) and healthy controls (25 subjects). The metabolomic profile of these subjects was characterized at multiple time points, and a comprehensive analysis was investigated between the following: MDD and healthy controls, treatment and placebo in both groups and the corresponding response to ketamine treatment. Ketamine treatment resulted in a general increase in circulating sphingomyelins, levels which were not correlated with response. Ketamine response resulted in more pronounced effects in the kynurenine pathway and the arginine pathway at 4 h post-infusion, where a larger decrease in circulating kynurenine levels and a larger increase in the bioavailability of arginine were observed in responders to ketamine treatment, suggesting possible mechanisms for response to ketamine treatment.

Natural polyphenols as sirtuin 6 modulators.

Flavonoids are polyphenolic secondary metabolites synthesized by plants and fungus with various pharmacological effects. Due to their plethora of biological activities, they have been studied extensively in drug development. They have been shown to modulate the activity of a NAD+-dependent histone deacetylase, SIRT6. Because SIRT6 has been implicated in longevity, metabolism, DNA-repair, and inflammatory response reduction, it is an interesting target in inflammatory and metabolic diseases as well as in cancer. Here we show, that flavonoids can alter SIRT6 activity in a structure dependent manner. Catechin derivatives with galloyl moiety displayed significant inhibition potency against SIRT6 at 10 µM concentration. The most potent SIRT6 activator, cyanidin, belonged to anthocyanidins, and produced a 55-fold increase in SIRT6 activity compared to the 3-10 fold increase for the others. Cyanidin also significantly increased SIRT6 expression in Caco-2 cells. Results from the docking studies indicated possible binding sites for the inhibitors and activators. Inhibitors likely bind in a manner that could disturb NAD+ binding. The putative activator binding site was found next to a loop near the acetylated peptide substrate binding site. In some cases, the activators changed the conformation of this loop suggesting that it may play a role in SIRT6 activation.

Synthesis and Biological Evaluation of Fluorescent Bryostatin Analogues.

To investigate the cellular distribution of tumor-promoting vs. non-tumor-promoting bryostatin analogues, we synthesized fluorescently labeled variants of two bryostatin derivatives that have previously shown either phorbol ester-like or bryostatin-like biological activity in U937 leukemia cells. These new fluorescent analogues both displayed high affinity for protein kinase C (PKC) binding and retained the basic properties of the parent unlabeled compounds in U937 assays. The fluorescent compounds showed similar patterns of intracellular distribution in cells, however; this argues against an existing hypothesis that various patterns of intracellular distribution are responsible for differences in biological activity. Upon further characterization, the fluorescent compounds revealed a slow rate of cellular uptake; correspondingly, they showed reduced activity for cellular responses that were only transient upon treatment with phorbol ester or bryostatin 1.

The evolutionary landscape of chronic lymphocytic leukemia treated with ibrutinib targeted therapy.

Treatment of chronic lymphocytic leukemia (CLL) has shifted from chemo-immunotherapy to targeted agents. To define the evolutionary dynamics induced by targeted therapy in CLL, we perform serial exome and transcriptome sequencing for 61 ibrutinib-treated CLLs. Here, we report clonal shifts (change >0.1 in clonal cancer cell fraction, Q < 0.1) in 31% of patients during the first year of therapy, associated with adverse outcome. We also observe transcriptional downregulation of pathways mediating energy metabolism, cell cycle, and B cell receptor signaling. Known and previously undescribed mutations in BTK and PLCG2, or uncommonly, other candidate alterations are present in seventeen subjects at the time of progression. Thus, the frequently observed clonal shifts during the early treatment period and its potential association with adverse outcome may reflect greater evolutionary capacity, heralding the emergence of drug-resistant clones.

Membrane Guanylate Cyclase catalytic Subdomain: Structure and Linkage with Calcium Sensors and Bicarbonate.

Membrane guanylate cyclase (MGC) is a ubiquitous multi-switching cyclic GMP generating signaling machine linked with countless physiological processes. In mammals it is encoded by seven distinct homologous genes. It is a single transmembrane spanning multi-modular protein; composed of integrated blocks and existing in homo-dimeric form. Its core catalytic domain (CCD) module is a common transduction center where all incoming signals are translated into the production of cyclic GMP, a cellular signal second messenger. Crystal structure of the MGC's CCD does not exist and its precise identity is ill-defined. Here, we define it at a sub-molecular level for the phototransduction-linked MGC, the rod outer segment guanylate cyclase type 1, ROS-GC1. (1) The CCD is a conserved 145-residue structural unit, represented by the segment V820-P964. (2) It exists as a homo-dimer and contains seven conserved catalytic elements (CEs) wedged into seven conserved motifs. (3) It also contains a conserved 21-residue neurocalcin δ-modulated structural domain, V836-L857. (4) Site-directed mutagenesis documents that each of the seven CEs governs the cyclase's catalytic activity. (5) In contrast to the soluble and the bacterium MGC which use Mn2+-GTP substrate for catalysis, MGC CCD uses the natural Mg2+-GTP substrate. (6) Strikingly, the MGC CCD requires anchoring by the Transmembrane Domain (TMD) to exhibit its major (∼92%) catalytic activity; in isolated form the activity is only marginal. This feature is not linked with any unique sequence of the TMD; there is minimal conservation in TMD. Finally, (7) the seven CEs control each of four phototransduction pathways- -two Ca2+-sensor GCAPs-, one Ca2+-sensor, S100B-, and one bicarbonate-modulated. The findings disclose that the CCD of ROS-GC1 has built-in regulatory elements that control its signal translational activity. Due to conservation of these regulatory elements, it is proposed that these elements also control the physiological activity of other members of MGC family.

Transcriptomic profiling and quantitative high-throughput (qHTS) drug screening of CDH1 deficient hereditary diffuse gastric cancer (HDGC) cells identify treatment leads for familial gastric cancer.

BACKGROUND: Patients with hereditary diffuse gastric cancer (HDGC), a cancer predisposition syndrome associated with germline mutations of the CDH1 (E-cadherin) gene, have few effective treatment options. Despite marked differences in natural history, histopathology, and genetic profile to patients afflicted by sporadic gastric cancer, patients with HDGC receive, in large, identical systemic regimens. The lack of a robust preclinical in vitro system suitable for effective drug screening has been one of the obstacles to date which has hampered therapeutic advances in this rare disease. METHODS: In order to identify therapeutic leads selective for the HDGC subtype of gastric cancer, we compared gene expression profiles and drug phenotype derived from an oncology library of 1912 compounds between gastric cancer cells established from a patient with metastatic HDGC harboring a c.1380delA CDH1 germline variant and sporadic gastric cancer cells. RESULTS: Unsupervised hierarchical cluster analysis shows select gene expression alterations in c.1380delA CDH1 SB.mhdgc-1 cells compared to a panel of sporadic gastric cancer cell lines with enrichment of ERK1-ERK2 (extracellular signal regulated kinase) and IP3 (inositol trisphosphate)/DAG (diacylglycerol) signaling as the top networks in c.1380delA SB.mhdgc-1 cells. Intracellular phosphatidylinositol intermediaries were increased upon direct measure in c.1380delA CDH1 SB.mhdgc-1 cells. Differential high-throughput drug screening of c.1380delA CDH1 SB.mhdgc-1 versus sporadic gastric cancer cells identified several compound classes with enriched activity in c.1380 CDH1 SB.mhdgc-1 cells including mTOR (Mammalian Target Of Rapamycin), MEK (Mitogen-Activated Protein Kinase), c-Src kinase, FAK (Focal Adhesion Kinase), PKC (Protein Kinase C), or TOPO2 (Topoisomerase II) inhibitors. Upon additional drug response testing, dual PI3K (Phosphatidylinositol 3-Kinase)/mTOR and topoisomerase 2A inhibitors displayed up to >100-fold increased activity in hereditary c.1380delA CDH1 gastric cancer cells inducing apoptosis most effectively in cells with deficient CDH1 function. CONCLUSION: Integrated pharmacological and transcriptomic profiling of hereditary diffuse gastric cancer cells with a loss-of-function c.1380delA CDH1 mutation implies various pharmacological vulnerabilities selective to CDH1-deficient familial gastric cancer cells and suggests novel treatment leads for future preclinical and clinical treatment studies of familial gastric cancer.

Somatic Genomics and Clinical Features of Lung Adenocarcinoma: A Retrospective Study.

BACKGROUND: Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer and has a high risk of distant metastasis at every disease stage. We aimed to characterize the genomic landscape of LUAD and identify mutation signatures associated with tumor progression. METHODS AND FINDINGS: We performed an integrative genomic analysis, incorporating whole exome sequencing (WES), determination of DNA copy number and DNA methylation, and transcriptome sequencing for 101 LUAD samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study. We detected driver genes by testing whether the nonsynonymous mutation rate was significantly higher than the background mutation rate and replicated our findings in public datasets with 724 samples. We performed subclonality analysis for mutations based on mutant allele data and copy number alteration data. We also tested the association between mutation signatures and clinical outcomes, including distant metastasis, survival, and tumor grade. We identified and replicated two novel candidate driver genes, POU class 4 homeobox 2 (POU4F2) (mutated in 9 [8.9%] samples) and ZKSCAN1 (mutated in 6 [5.9%] samples), and characterized their major deleterious mutations. ZKSCAN1 was part of a mutually exclusive gene set that included the RTK/RAS/RAF pathway genes BRAF, EGFR, KRAS, MET, and NF1, indicating an important driver role for this gene. Moreover, we observed strong associations between methylation in specific genomic regions and somatic mutation patterns. In the tumor evolution analysis, four driver genes had a significantly lower fraction of subclonal mutations (FSM), including TP53 (p = 0.007), KEAP1 (p = 0.012), STK11 (p = 0.0076), and EGFR (p = 0.0078), suggesting a tumor initiation role for these genes. Subclonal mutations were significantly enriched in APOBEC-related signatures (p < 2.5×10-50). The total number of somatic mutations (p = 0.0039) and the fraction of transitions (p = 5.5×10-4) were associated with increased risk of distant metastasis. Our study's limitations include a small number of LUAD patients for subgroup analyses and a single-sample design for investigation of subclonality. CONCLUSIONS: These data provide a genomic characterization of LUAD pathogenesis and progression. The distinct clonal and subclonal mutation signatures suggest possible diverse carcinogenesis pathways for endogenous and exogenous exposures, and may serve as a foundation for more effective treatments for this lethal disease. LUAD's high heterogeneity emphasizes the need to further study this tumor type and to associate genomic findings with clinical outcomes.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Identification of IL-23p19 as an endothelial proinflammatory peptide that promotes gp130-STAT3 signaling.

Interleukin-23 (IL-23), a heterodimeric cytokine composed of the unique p19 peptide (IL-23p19) and a peptide called IL-12p40, which is shared with IL-12, is implicated in Crohn's disease, rheumatoid arthritis, psoriasis, and other immune-mediated inflammatory diseases. Endothelial cells produce the IL-23p19 peptide in the absence of the IL-12p40 chain and thus do not make heterodimeric IL-23. We found that intercellular IL-23p19 increased the cell surface abundances of intercellular adhesion molecule-1 (ICAM-1) and vascular cell adhesion molecule-1 (VCAM-1) on endothelial cells, which enhanced the attachment of leukocytes and increased their transendothelial migration. Intracellular p19 associated with the cytokine receptor subunit gp130 and stimulated the gp130-dependent activation of signal transducer and activator of transcription 3 (STAT3) signaling. Proinflammatory factors promoted the generation of IL-23p19 in endothelial cells. The adventitial capillaries of inflamed temporal arteries in patients with giant-cell arteritis (GCA) had endothelial p19 protein associated with gp130, but did not contain the IL-12p40 chain. Because adventitial capillaries are essential for the entry of inflammatory cells into arterial walls, these data suggest that p19 may contribute to GCA disease and could represent a therapeutic target. Our results provide evidence that IL-23p19 is a previously unrecognized endothelial proinflammatory peptide that promotes leukocyte transendothelial migration, advancing our current understanding of the complexities of inflammatory responses.

Whole exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We used whole exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mu-tant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

Deep sequencing of HPV16 genomes: A new high-throughput tool for exploring the carcinogenicity and natural history of HPV16 infection.

For unknown reasons, there is huge variability in risk conferred by different HPV types and, remarkably, strong differences even between closely related variant lineages within each type. HPV16 is a uniquely powerful carcinogenic type, causing approximately half of cervical cancer and most other HPV-related cancers. To permit the large-scale study of HPV genome variability and precancer/cancer, starting with HPV16 and cervical cancer, we developed a high-throughput next-generation sequencing (NGS) whole-genome method. We designed a custom HPV16 AmpliSeq™ panel that generated 47 overlapping amplicons covering 99% of the genome sequenced on the Ion Torrent Proton platform. After validating with Sanger, the current "gold standard" of sequencing, in 89 specimens with concordance of 99.9%, we used our NGS method and custom annotation pipeline to sequence 796 HPV16-positive exfoliated cervical cell specimens. The median completion rate per sample was 98.0%. Our method enabled us to discover novel SNPs, large contiguous deletions suggestive of viral integration (OR of 27.3, 95% CI 3.3-222, P=0.002), and the sensitive detection of variant lineage coinfections. This method represents an innovative high-throughput, ultra-deep coverage technique for HPV genomic sequencing, which, in turn, enables the investigation of the role of genetic variation in HPV epidemiology and carcinogenesis.

Juvenile myelomonocytic leukemia due to a germline CBL Y371C mutation: 35-year follow-up of a large family.

Juvenile myelomonocytic leukemia (JMML) is a pediatric myeloproliferative neoplasm that arises from malignant transformation of the stem cell compartment and results in increased production of myeloid cells. Somatic and germline variants in CBL (Casitas B-lineage lymphoma proto-oncogene) have been associated with JMML. We report an incompletely penetrant CBL Y371C mutation discovered by whole-exome sequencing in three individuals with JMML in a large pedigree with 35 years of follow-up. The Y371 residue is highly evolutionarily conserved among CBL orthologs and paralogs. In silico bioinformatics prediction programs suggested that the Y371C mutation is highly deleterious. Protein structural modeling revealed that the Y371C mutation abrogated the ability of the CBL protein to adopt a conformation that is required for ubiquitination. Clinically, the three mutation-positive JMML individuals exhibited variable clinical courses; in two out of three, primary hematologic abnormalities persisted into adulthood with minimal clinical symptoms. The penetrance of the CBL Y371C mutation was 30% for JMML and 40% for all leukemia. Of the 8 mutation carriers in the family with available photographs, only one had significant dysmorphic features; we found no evidence of a clinical phenotype consistent with a "CBL syndrome". Although CBL Y371C has been previously reported in familial JMML, we are the first group to follow a complete pedigree harboring this mutation for an extended period, revealing additional information about this variant's penetrance, function and natural history.

AVIA v2.0: annotation, visualization and impact analysis of genomic variants and genes.

UNLABELLED: As sequencing becomes cheaper and more widely available, there is a greater need to quickly and effectively analyze large-scale genomic data. While the functionality of AVIA v1.0, whose implementation was based on ANNOVAR, was comparable with other annotation web servers, AVIA v2.0 represents an enhanced web-based server that extends genomic annotations to cell-specific transcripts and protein-level functional annotations. With AVIA's improved interface, users can better visualize their data, perform comprehensive searches and categorize both coding and non-coding variants. AVAILABILITY AND IMPLEMENTATION: AVIA is freely available through the web at http://avia.abcc.ncifcrf.gov. CONTACT: Hue.Vuong@fnlcr.nih.gov SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Can structural features of kinase receptors provide clues on selectivity and inhibition? A molecular modeling study.

Cancer is a complex disease resulting from the uncontrolled proliferation of cell signaling events. Protein kinases have been identified as central molecules that participate overwhelmingly in oncogenic events, thus becoming key targets for anticancer drugs. A majority of studies converged on the idea that ligand-binding pockets of kinases retain clues to the inhibiting abilities and cross-reacting tendencies of inhibitor drugs. Even though these ideas are critical for drug discovery, validating them using experiments is not only difficult, but also in some cases infeasible. To overcome these limitations and to test these ideas at the molecular level, we present here the results of receptor-focused in-silico docking of nine marketed drugs to 19 different wild-type and mutated kinases chosen from a wide range of families. This investigation highlights the need for using relevant models to explain the correct inhibition trends and the results are used to make predictions that might be able to influence future experiments. Our simulation studies are able to correctly predict the primary targets for each drug studied in majority of cases and our results agree with the existing findings. Our study shows that the conformations a given receptor acquires during kinase activation, and their micro-environment, defines the ligand partners. Type II drugs display high compatibility and selectivity for DFG-out kinase conformations. On the other hand Type I drugs are less selective and show binding preferences for both the open and closed forms of selected kinases. Using this receptor-focused approach, it is possible to capture the observed fold change in binding affinities between the wild-type and disease-centric mutations in ABL kinase for Imatinib and the second-generation ABL drugs. The effects of mutation are also investigated for two other systems, EGFR and B-Raf. Finally, by including pathway information in the design it is possible to model kinase inhibitors with potentially fewer side-effects.