Discovery and Characterization of Linear Motif Mediated Protein-Protein Complexes.

There are myriads of protein-protein complexes that form within the cell. In addition to classical binding events between globular domains, many protein-protein interactions involve short disordered protein regions. The latter contain so-called linear motifs binding specifically to ordered protein domain surfaces. Linear binding motifs are classified based on their consensus sequence, where only a few amino acids are conserved. In this chapter we will review experimental and in silico techniques that can be used for the discovery and characterization of linear motif mediated protein-protein complexes involved in cellular signaling, protein level and gene expression regulation.

Phosphorylation-Assisted Luciferase Complementation Assay Designed to Monitor Kinase Activity and Kinase-Domain-Mediated Protein-Protein Binding.

Protein kinases are key regulators of cell signaling and have been important therapeutic targets for three decades. ATP-competitive drugs directly inhibit the activity of kinases but these enzymes work as part of complex protein networks in which protein-protein interactions (often referred to as kinase docking) may govern a more complex activation pattern. Kinase docking is indispensable for many signaling disease-relevant Ser/Thr kinases and it is mediated by a dedicated surface groove on the kinase domain which is distinct from the substrate-binding pocket. Thus, interfering with kinase docking provides an alternative strategy to control kinases. We describe activity sensors developed for p90 ribosomal S6 kinase (RSK) and mitogen-activated protein kinases (MAPKs: ERK, p38, and JNK) whose substrate phosphorylation is known to depend on kinase-docking-groove-mediated protein-protein binding. The in vitro assays were based on fragment complementation of the NanoBit luciferase, which is facilitated upon substrate motif phosphorylation. The new phosphorylation-assisted luciferase complementation (PhALC) sensors are highly selective and the PhALC assay is a useful tool for the quantitative analysis of kinase activity or kinase docking, and even for high-throughput screening of academic compound collections.

Structural insights into the pSer/pThr dependent regulation of the SHP2 tyrosine phosphatase in insulin and CD28 signaling.

Serine/threonine phosphorylation of insulin receptor substrate (IRS) proteins is well known to modulate insulin signaling. However, the molecular details of this process have mostly been elusive. While exploring the role of phosphoserines, we have detected a direct link between Tyr-flanking Ser/Thr phosphorylation sites and regulation of specific phosphotyrosine phosphatases. Here we present a concise structural study on how the activity of SHP2 phosphatase is controlled by an asymmetric, dual phosphorylation of its substrates. The structure of SHP2 has been determined with three different substrate peptides, unveiling the versatile and highly dynamic nature of substrate recruitment. What is more, the relatively stable pre-catalytic state of SHP2 could potentially be useful for inhibitor design. Our findings not only show an unusual dependence of SHP2 catalytic activity on Ser/Thr phosphorylation sites in IRS1 and CD28, but also suggest a negative regulatory mechanism that may also apply to other tyrosine kinase pathways as well.

Incorporation of Oxidized Phenylalanine Derivatives into Insulin Signaling Relevant Proteins May Link Oxidative Stress to Signaling Conditions Underlying Chronic Insulin Resistance.

A link between oxidative stress and insulin resistance has been suggested. Hydroxyl free radicals are known to be able to convert phenylalanine (Phe) into the non-physiological tyrosine isoforms ortho- and meta-tyrosine (o-Tyr, m-Tyr). The aim of our study was to examine the role of o-Tyr and m-Tyr in the development of insulin resistance. We found that insulin-induced uptake of glucose was blunted in cultures of 3T3-L1 grown on media containing o- or m-Tyr. We show that these modified amino acids are incorporated into cellular proteins. We focused on insulin receptor substrate 1 (IRS-1), which plays a role in insulin signaling. The activating phosphorylation of IRS-1 was increased by insulin, the effect of which was abolished in cells grown in m-Tyr or o-Tyr media. We found that phosphorylation of m- or o-Tyr containing IRS-1 segments by insulin receptor (IR) kinase was greatly reduced, PTP-1B phosphatase was incapable of dephosphorylating phosphorylated m- or o-Tyr IRS-1 peptides, and the SH2 domains of phosphoinositide 3-kinase (PI3K) bound the o-Tyr IRS-1 peptides with greatly reduced affinity. According to our data, m- or o-Tyr incorporation into IRS-1 modifies its protein-protein interactions with regulating enzymes and effectors, thus IRS-1 eventually loses its capacity to play its role in insulin signaling, leading to insulin resistance.

Systematic Discovery of FBXW7-Binding Phosphodegrons Highlights Mitogen-Activated Protein Kinases as Important Regulators of Intracellular Protein Levels.

A FBXW7 is an F-box E3 ubiquitin-ligase affecting cell growth by controlling protein degradation. Mechanistically, its effect on its substrates depends on the phosphorylation of degron motifs, but the abundance of these phosphodegrons has not been systematically explored. We used a ratiometric protein degradation assay geared towards the identification of FBXW7-binding degron motifs phosphorylated by mitogen-activated protein kinases (MAPKs). Most of the known FBXW7 targets are localized in the nucleus and function as transcription factors. Here, in addition to more transcription affecting factors (ETV5, KLF4, SP5, JAZF1, and ZMIZ1 CAMTA2), we identified phosphodegrons located in proteins involved in chromatin regulation (ARID4B, KMT2E, KMT2D, and KAT6B) or cytoskeletal regulation (MAP2, Myozenin-2, SMTL2, and AKAP11), and some other proteins with miscellaneous functions (EIF4G3, CDT1, and CCAR2). We show that the protein level of full-length ARID4B, ETV5, JAZF1, and ZMIZ1 are affected by different MAPKs since their FBXW7-mediated degradation was diminished in the presence of MAPK-specific inhibitors. Our results suggest that MAPK and FBXW7 partnership plays an important cellular role by directly affecting the level of key regulatory proteins. The data also suggest that the p38α-controlled phosphodegron in JAZF1 may be responsible for the pathological regulation of the cancer-related JAZF1-SUZ12 fusion construct implicated in endometrial stromal sarcoma.

A non-catalytic herpesviral protein reconfigures ERK-RSK signaling by targeting kinase docking systems in the host.

The Kaposi's sarcoma associated herpesvirus protein ORF45 binds the extracellular signal-regulated kinase (ERK) and the p90 Ribosomal S6 kinase (RSK). ORF45 was shown to be a kinase activator in cells but a kinase inhibitor in vitro, and its effects on the ERK-RSK complex are unknown. Here, we demonstrate that ORF45 binds ERK and RSK using optimized linear binding motifs. The crystal structure of the ORF45-ERK2 complex shows how kinase docking motifs recognize the activated form of ERK. The crystal structure of the ORF45-RSK2 complex reveals an AGC kinase docking system, for which we provide evidence that it is functional in the host. We find that ORF45 manipulates ERK-RSK signaling by favoring the formation of a complex, in which activated kinases are better protected from phosphatases and docking motif-independent RSK substrate phosphorylation is selectively up-regulated. As such, our data suggest that ORF45 interferes with the natural design of kinase docking systems in the host.

Characterization of the Intramolecular Interactions and Regulatory Mechanisms of the Scaffold Protein Tks4.

The scaffold protein Tks4 is a member of the p47phox-related organizer superfamily. It plays a key role in cell motility by being essential for the formation of podosomes and invadopodia. In addition, Tks4 is involved in the epidermal growth factor (EGF) signaling pathway, in which EGF induces the translocation of Tks4 from the cytoplasm to the plasma membrane. The evolutionarily-related protein p47phox and Tks4 share many similarities in their N-terminal region: a phosphoinositide-binding PX domain is followed by two SH3 domains (so called "tandem SH3") and a proline-rich region (PRR). In p47phox, the PRR is followed by a relatively short, disordered C-terminal tail region containing multiple phosphorylation sites. These play a key role in the regulation of the protein. In Tks4, the PRR is followed by a third and a fourth SH3 domain connected by a long (~420 residues) unstructured region. In p47phox, the tandem SH3 domain binds the PRR while the first SH3 domain interacts with the PX domain, thereby preventing its binding to the membrane. Based on the conserved structural features of p47phox and Tks4 and the fact that an intramolecular interaction between the third SH3 and the PX domains of Tks4 has already been reported, we hypothesized that Tks4 is similarly regulated by autoinhibition. In this study, we showed, via fluorescence-based titrations, MST, ITC, and SAXS measurements, that the tandem SH3 domain of Tks4 binds the PRR and that the PX domain interacts with the third SH3 domain. We also investigated a phosphomimicking Thr-to-Glu point mutation in the PRR as a possible regulator of intramolecular interactions. Phosphatidylinositol-3-phosphate (PtdIns(3)P) was identified as the main binding partner of the PX domain via lipid-binding assays. In truncated Tks4 fragments, the presence of the tandem SH3, together with the PRR, reduced PtdIns(3)P binding, while the presence of the third SH3 domain led to complete inhibition.

Peptide Based Inhibitors of Protein Binding to the Mitogen-Activated Protein Kinase Docking Groove.

Mitogen-activated protein kinases (MAPK) are important regulatory units in cells and they take part in the regulation of many cellular functions such as cell division, differentiation or apoptosis. All MAPKs have a shallow docking groove that interacts with linear binding motifs of their substrate proteins and their regulatory proteins such as kinases, phosphatases, scaffolds. Inhibition of these protein-protein interactions may reduce or abolish the activity of the targeted kinase. Based on the wide range of their biological activity, this kind of inhibition can be useful in the treatment of many disorders like tumors, inflammation or undesired cell apoptosis. In this study a linear binding motif from the RHDF1 protein-a 15 amino acids long peptide-was selected for optimization to increase its cellular uptake but retaining its low micromolar binding affinity. First, we synthesized an octaarginine conjugate that showed efficient cellular uptake. Next, we set out to reduce the size of this construct. We were able to decrease the length of the original peptide, and to increase its cellular uptake with specific chemical modifications. These new constructs bound better to ERK2 and p38 kinases than the original peptide and they showed markedly increased cellular uptake. The new octaarginine conjugate and one of the minimized bicyclic derivatives could inhibit the phosphorylation of intracellular ERK or p38. However, the modulation of MAPK phosphorylation levels by these cell-penetrating peptides were complex, despite that in biochemical assays they all inhibited MAPK-substrate binding as well as phosphorylation. The optimized peptides depending on the applied concentration caused an expected decrease, but also some unexpected increase in MAPK phosphorylation patterns in the cell. This possibly reflects the complexity of MAPK docking groove mediated protein-protein interactions including bone fide MAPK clients such activator kinases, deactivating phosphatases or regulatory scaffolds. Thus, our findings with optimized cell-penetrating "inhibitory" peptides highlight the opportunities but also the pitfalls of docking peptide based MAPK activity regulation and call for a better quantitative understanding of MAPK mediated protein-protein interactions in cells.

Co-regulation of the transcription controlling ATF2 phosphoswitch by JNK and p38.

Transcription factor phosphorylation at specific sites often activates gene expression, but how environmental cues quantitatively control transcription is not well-understood. Activating protein 1 transcription factors are phosphorylated by mitogen-activated protein kinases (MAPK) in their transactivation domains (TAD) at so-called phosphoswitches, which are a hallmark in response to growth factors, cytokines or stress. We show that the ATF2 TAD is controlled by functionally distinct signaling pathways (JNK and p38) through structurally different MAPK binding sites. Moreover, JNK mediated phosphorylation at an evolutionarily more recent site diminishes p38 binding and made the phosphoswitch differently sensitive to JNK and p38 in vertebrates. Structures of MAPK-TAD complexes and mechanistic modeling of ATF2 TAD phosphorylation in cells suggest that kinase binding motifs and phosphorylation sites line up to maximize MAPK based co-regulation. This study shows how the activity of an ancient transcription controlling phosphoswitch became dependent on the relative flux of upstream signals.

MAP Kinase-Mediated Activation of RSK1 and MK2 Substrate Kinases.

Mitogen-activated protein kinases (MAPKs) control essential eukaryotic signaling pathways. While much has been learned about MAPK activation, much less is known about substrate recruitment and specificity. MAPK substrates may be other kinases that are crucial to promote a further diversification of the signaling outcomes. Here, we used a variety of molecular and cellular tools to investigate the recruitment of two substrate kinases, RSK1 and MK2, to three MAPKs (ERK2, p38α, and ERK5). Unexpectedly, we identified that kinase heterodimers form structurally and functionally distinct complexes depending on the activation state of the MAPK. These may be incompatible with downstream signaling, but naturally they may also form structures that are compatible with the phosphorylation of the downstream kinase at the activation loop, or alternatively at other allosteric sites. Furthermore, we show that small-molecule inhibitors may affect the quaternary arrangement of kinase heterodimers and thus influence downstream signaling in a specific manner.

P38 MAPK Promotes Migration and Metastatic Activity of BRAF Mutant Melanoma Cells by Inducing Degradation of PMCA4b.

Metastatic melanoma is the most aggressive type of skin cancer. Previously, we identified the plasma membrane Ca2+ pump isoform 4b (PMCA4b or ATP2B4) as a putative metastasis suppressor in BRAF mutant melanoma cells. Metastasis suppressors are often downregulated in cancer, therefore, it is important to identify the pathways involved in their degradation. Here, we studied the role of p38 MAPK in PMCA4b degradation and its effect on melanoma metastasis. We found that activation of p38 MAPK induces internalization and subsequent degradation of PMCA4b through the endo/lysosomal system that contributes to the low PMCA4b steady-state protein level of BRAF mutant melanoma cells. Moreover, BRAF wild type cell models including a doxycycline-inducible HEK cell system revealed that p38 MAPK is a universal modulator of PMCA4b endocytosis. Inhibition of the p38 MAPK pathway markedly reduced migration, colony formation and metastatic activity of BRAF mutant cells in vitro partially through an increase in PMCA4b and a decrease in ß4 integrin abundance. In conclusion, our data suggest that the p38 MAPK pathway plays a key role in PMCA4b degradation and inhibition of this pathway-by increasing the stability of PMCA4b-may provide a potential therapeutic target for inhibition of melanoma progression and metastasis.

Ndr/Lats Kinases Bind Specific Mob-Family Coactivators through a Conserved and Modular Interface.

Ndr/Lats kinases bind Mob coactivator proteins to form complexes that are essential and evolutionarily conserved components of "Hippo" signaling pathways, which control cell proliferation and morphogenesis in eukaryotes. All Ndr/Lats kinases have a characteristic N-terminal regulatory (NTR) region that binds a specific Mob cofactor: Lats kinases associate with Mob1 proteins, and Ndr kinases associate with Mob2 proteins. To better understand the significance of the association of Mob protein with Ndr/Lats kinases and selective binding of Ndr and Lats to distinct Mob cofactors, we determined crystal structures of Saccharomyces cerevisiae Cbk1NTR-Mob2 and Dbf2NTR-Mob1 and experimentally assessed determinants of Mob cofactor binding and specificity. This allowed a significant improvement in the previously determined structure of Cbk1 kinase bound to Mob2, presently the only crystallographic model of a full length Ndr/Lats kinase complexed with a Mob cofactor. Our analysis indicates that the Ndr/LatsNTR-Mob interface provides a distinctive kinase regulation mechanism, in which the Mob cofactor organizes the Ndr/Lats NTR to interact with the AGC kinase C-terminal hydrophobic motif (HM), which is involved in allosteric regulation. The Mob-organized NTR appears to mediate association of the HM with an allosteric site on the N-terminal kinase lobe. We also found that Cbk1 and Dbf2 associated specifically with Mob2 and Mob1, respectively. Alteration of residues in the Cbk1 NTR allows association of the noncognate Mob cofactor, indicating that cofactor specificity is restricted by discrete sites rather than being broadly distributed. Overall, our analysis provides a new picture of the functional role of Mob association and indicates that the Ndr/LatsNTR-Mob interface is largely a common structural platform that mediates kinase-cofactor binding.

High-throughput competitive fluorescence polarization assay reveals functional redundancy in the S100 protein family.

The calcium-binding, vertebrate-specific S100 protein family consists of 20 paralogs in humans (referred as the S100ome), with several clinically important members. To explore their protein-protein interactions (PPIs) quantitatively, we have chosen an unbiased, high-throughput, competitive fluorescence polarization (FP) assay that revealed a partial functional redundancy when the complete S100ome (n = 20) was tested against numerous model partners (n = 13). Based on their specificity, the S100ome can be grouped into two distinct classes: promiscuous and orphan. In the first group, members bound to several ligands (> 4-5) with comparable high affinity, while in the second one, the paralogs bound only one partner weakly, or no ligand was identified. Our results demonstrate that FP assays are highly suitable for quantitative interaction profiling of selected protein families. Moreover, we provide evidence that PPI-based phenotypic characterization can complement or even exceed the information obtained from the sequence-based phylogenetic analysis of the S100ome, an evolutionary young protein family.

Rewiring of RSK-PDZ Interactome by Linear Motif Phosphorylation.

Phosphorylation of short linear peptide motifs is a widespread process for the dynamic regulation of protein-protein interactions. However, the global impact of phosphorylation events on the protein-protein interactome is rarely addressed. The disordered C-terminal tail of ribosomal S6 kinase 1 (RSK1) binds to PDZ domain-containing scaffold proteins, and it harbors a phosphorylatable PDZ-binding motif (PBM) responsive to epidermal growth factor stimulation. Here, we examined binding of two versions of the RSK1 PBM, either phosphorylated or unphosphorylated at position -3, to almost all (95%) of the 266 PDZ domains of the human proteome. PBM phosphorylation dramatically altered the PDZ domain-binding landscape of RSK1, by strengthening or weakening numerous interactions to various degrees. The RSK-PDZome interactome analyzed in this study reveals how linear motif-based phospho-switches convey stimulus-dependent changes in the context of related network components.

Disordered Protein Kinase Regions in Regulation of Kinase Domain Cores.

Since publication of the crystal structure of protein kinase (PK)A three decades ago, a structural portrait of the conserved kinase core has been drawn. The next challenge is to elucidate structures of full-length kinases and to address the intrinsically disordered regions (IDRs) that typically flank the core as well as the small linear motifs (SLiMs) that are embedded within the IDRs. It is increasingly apparent that unstructured regions integrate the kinase catalytic chassis into multienzyme-based regulatory networks. The extracellular signal-regulated kinase-ribosomal S6 PK-phosphoinositide-dependent kinase (ERK-RSK-PDK) complex is an excellent example to demonstrate how IDRs and SLiMs govern communication between four different kinase catalytic cores to mediate activation and how in molecular terms these promote the formation of kinase heterodimers in a context dependent fashion.

Theileria highjacks JNK2 into a complex with the macroschizont GPI (GlycosylPhosphatidylInositol)-anchored surface protein p104.

Constitutive c-Jun N-terminal kinase (JNK) activity characterizes bovine T and B cells infected with Theileria parva, and B cells and macrophages infected with Theileria annulata. Here, we show that T. annulata infection of macrophages manipulates JNK activation by recruiting JNK2 and not JNK1 to the parasite surface, whereas JNK1 is found predominantly in the host cell nucleus. At the parasite's surface, JNK2 forms a complex with p104, a GPI-(GlycosylPhosphatidylInositol)-anchor T. annulata plasma membrane protein. Sequestration of JNK2 depended on Protein Kinase-A (PKA)-mediated phosphorylation of a JNK-binding motif common to T. parva and a cell penetrating peptide harbouring the conserved p104 JNK-binding motif competitively ablated binding, whereupon liberated JNK2 became ubiquitinated and degraded. Cytosolic sequestration of JNK2 suppressed small mitochondrial ARF-mediated autophagy, whereas it sustained nuclear JNK1 levels, c-Jun phosphorylation, and matrigel traversal. Therefore, T. annulata sequestration of JNK2 contributes to both survival and dissemination of Theileria-transformed macrophages.

Dynamic control of RSK complexes by phosphoswitch-based regulation.

Assembly and disassembly of protein-protein complexes needs to be dynamically controlled and phosphoswitches based on linear motifs are crucial in this process. Extracellular signal-regulated kinase 2 (ERK2) recognizes a linear-binding motif at the C-terminal tail (CTT) of ribosomal S6 kinase 1 (RSK1), leading to phosphorylation and subsequent activation of RSK1. The CTT also contains a classical PDZ domain-binding motif which binds RSK substrates (e.g. MAGI-1). We show that autophosphorylation of the disordered CTT promotes the formation of an intramolecular charge clamp, which efficiently masks critical residues and indirectly hinders ERK binding. Thus, RSK1 CTT operates as an autoregulated phosphoswitch: its phosphorylation at specific sites affects its protein-binding capacity and its conformational dynamics. These biochemical feedbacks, which form the structural basis for the rapid dissociation of ERK2-RSK1 and RSK1-PDZ substrate complexes under sustained epidermal growth factor (EGF) stimulation, were structurally characterized and validated in living cells. Overall, conformational changes induced by phosphorylation in disordered regions of protein kinases, coupled to allosteric events occurring in the kinase domain cores, may provide mechanisms that contribute to the emergence of complex signaling activities. In addition, we show that phosphoswitches based on linear motifs can be functionally classified as ON and OFF protein-protein interaction switches or dimmers, depending on the specific positioning of phosphorylation target sites in relation to functional linear-binding motifs. Moreover, interaction of phosphorylated residues with positively charged residues in disordered regions is likely to be a common mechanism of phosphoregulation. DATABASE: Structural data are available in the PDB database under the accession numbers 5N7D, 5N7F and 5N7G. NMR spectral assignation data are available in the BMRB database under the accession numbers 27213 and 27214.

The eukaryotic linear motif resource - 2018 update.

Short linear motifs (SLiMs) are protein binding modules that play major roles in almost all cellular processes. SLiMs are short, often highly degenerate, difficult to characterize and hard to detect. The eukaryotic linear motif (ELM) resource (elm.eu.org) is dedicated to SLiMs, consisting of a manually curated database of over 275 motif classes and over 3000 motif instances, and a pipeline to discover candidate SLiMs in protein sequences. For 15 years, ELM has been one of the major resources for motif research. In this database update, we present the latest additions to the database including 32 new motif classes, and new features including Uniprot and Reactome integration. Finally, to help provide cellular context, we present some biological insights about SLiMs in the cell cycle, as targets for bacterial pathogenicity and their functionality in the human kinome.

JNK Signaling: Regulation and Functions Based on Complex Protein-Protein Partnerships.

The c-Jun N-terminal kinases (JNKs), as members of the mitogen-activated protein kinase (MAPK) family, mediate eukaryotic cell responses to a wide range of abiotic and biotic stress insults. JNKs also regulate important physiological processes, including neuronal functions, immunological actions, and embryonic development, via their impact on gene expression, cytoskeletal protein dynamics, and cell death/survival pathways. Although the JNK pathway has been under study for >20 years, its complexity is still perplexing, with multiple protein partners of JNKs underlying the diversity of actions. Here we review the current knowledge of JNK structure and isoforms as well as the partnerships of JNKs with a range of intracellular proteins. Many of these proteins are direct substrates of the JNKs. We analyzed almost 100 of these target proteins in detail within a framework of their classification based on their regulation by JNKs. Examples of these JNK substrates include a diverse assortment of nuclear transcription factors (Jun, ATF2, Myc, Elk1), cytoplasmic proteins involved in cytoskeleton regulation (DCX, Tau, WDR62) or vesicular transport (JIP1, JIP3), cell membrane receptors (BMPR2), and mitochondrial proteins (Mcl1, Bim). In addition, because upstream signaling components impact JNK activity, we critically assessed the involvement of signaling scaffolds and the roles of feedback mechanisms in the JNK pathway. Despite a clarification of many regulatory events in JNK-dependent signaling during the past decade, many other structural and mechanistic insights are just beginning to be revealed. These advances open new opportunities to understand the role of JNK signaling in diverse physiological and pathophysiological states.

Systematic analysis of somatic mutations driving cancer: uncovering functional protein regions in disease development.

BACKGROUND: Recent advances in sequencing technologies enable the large-scale identification of genes that are affected by various genetic alterations in cancer. However, understanding tumor development requires insights into how these changes cause altered protein function and impaired network regulation in general and/or in specific cancer types. RESULTS: In this work we present a novel method called iSiMPRe that identifies regions that are significantly enriched in somatic mutations and short in-frame insertions or deletions (indels). Applying this unbiased method to the complete human proteome, by using data enriched through various cancer genome projects, we identified around 500 protein regions which could be linked to one or more of 27 distinct cancer types. These regions covered the majority of known cancer genes, surprisingly even tumor suppressors. Additionally, iSiMPRe also identified novel genes and regions that have not yet been associated with cancer. CONCLUSIONS: While local somatic mutations correspond to only a subset of genetic variations that can lead to cancer, our systematic analyses revealed that they represent an accompanying feature of most cancer driver genes regardless of the primary mechanism by which they are perturbed during tumorigenesis. These results indicate that the accumulation of local somatic mutations can be used to pinpoint genes responsible for cancer formation and can also help to understand the effect of cancer mutations at the level of functional modules in a broad range of cancer driver genes. REVIEWERS: This article was reviewed by Sándor Pongor, Michael Gromiha and Zoltán Gáspári.