BACKGROUND: Determining the impact of somatic mutations requires understanding the functional relationship of genes acquiring mutations; however, it is largely unknown how mutations in functionally related genes influence each other.
METHODS: We employed non-synonymous-to-synonymous or dNdS ratios to evaluate the evolutionary dependency (ED) of gene pairs, assuming a mutation in one gene of a gene pair can affect the evolutionary fitness of mutations in its partner genes as mutation context. We employed PanCancer- and tumor type-specific mutational profiles to infer the ED of gene pairs and evaluated their biological relevance with respect to gene dependency and drug sensitivity.
RESULTS: We propose that dNdS ratios of gene pairs and their derived cdNS (context-dependent dNdS) scores as measure of ED distinguishing gene pairs either as synergistic (SYN) or antagonistic (ANT). Mutation contexts can induce substantial changes in the evolutionary fitness of mutations in the paired genes, e.g., IDH1 and IDH2 mutation contexts lead to substantial increase and decrease of dNdS ratios of ATRX indels and IDH1 missense mutations corresponding to SYN and ANT relationship with positive and negative cdNS scores, respectively. The impact of gene silencing or knock-outs on cell viability (genetic dependencies) often depends on ED, suggesting that ED can guide the selection of candidates for synthetic lethality such as TCF7L2-KRAS mutations. Using cell line-based drug sensitivity data, the effects of targeted agents on cell lines are often associated with mutations of genes exhibiting ED with the target genes, informing drug sensitizing or resistant mutations for targeted inhibitors, e.g., PRSS1 and CTCF mutations as resistant mutations to EGFR and BRAF inhibitors for lung adenocarcinomas and melanomas, respectively.
CONCLUSIONS: We propose that the ED of gene pairs evaluated by dNdS ratios can advance our understanding of the functional relationship of genes with potential biological and clinical implications.

The development of cancer involves the accumulation of somatic mutations in several essential biological pathways. Delineating the temporal order of pathway mutations during tumorigenesis is crucial for comprehending the biological mechanisms underlying cancer development and identifying potential targets for therapeutic intervention. Several computational and statistical methods have been introduced for estimating the order of somatic mutations based on mutation profile data from a cohort of patients. However, one major issue of current methods is that they do not take into account intra-tumor heterogeneity (ITH), which limits their ability to accurately discern the order of pathway mutations. To address this problem, we propose PATOPAI, a probabilistic approach to estimate the temporal order of mutations at the pathway level by incorporating ITH information as well as pathway and functional annotation information of mutations. PATOPAI uses a maximum likelihood approach to estimate the probability of pathway mutational events occurring in a specific sequence, wherein it focuses on the orders that are consistent with the phylogenetic structure of the tumors. Applications to whole exome sequencing data from The Cancer Genome Atlas (TCGA) illustrate our method\'s ability to recover the temporal order of pathway mutations in several cancer types.

κB-Ras (NF-κB inhibitor-interacting Ras-like protein) GTPases are small Ras-like GTPases but harbor interesting differences in important sequence motifs. They act in a tumor-suppressive manner as negative regulators of Ral (Ras-like) GTPase and NF-κB signaling, but little is known about their mode of function. Here, we demonstrate that, in contrast to predictions based on primary structure, κB-Ras GTPases possess hydrolytic activity. Combined with low nucleotide affinity, this renders them fast-cycling GTPases that are predominantly GTP-bound in cells. We characterize the impact of κB-Ras mutations occurring in tumors and demonstrate that nucleotide binding affects κB-Ras stability but is not strictly required for RalGAP (Ral GTPase-activating protein) binding. This demonstrates that κB-Ras control of RalGAP/Ral signaling occurs in a nucleotide-binding- and switch-independent fashion.

mRNAs containing premature stop codons are responsible for various genetic diseases as well as cancers. The truncated proteins synthesized from these aberrant mRNAs are seldom detected due to the nonsense-mediated mRNA decay (NMD) pathway. Such a surveillance mechanism detects most of these aberrant mRNAs and rapidly destroys them from the pool of mRNAs. Here, we implemented chemical cross-linking mass spectrometry (CLMS) techniques to trace novel biology consisting of protein-protein interactions (PPIs) within the NMD machinery. A set of novel complex networks between UPF2 (Regulator of nonsense transcripts 2), SMG1 (Serine/threonine-protein kinase SMG1), and SMG7 from the NMD pathway were identified, among which UPF2 was found as a connection bridge between SMG1 and SMG7. The UPF2 N-terminal formed most interactions with SMG7, and a set of residues emerged from the MIF4G-I, II, and III domains docked with SMG1 or SMG7. SMG1 mediated interactions with initial residues of UPF2, whereas SMG7 formed very few interactions in this region. Modelled structures highlighted that PPIs for UPF2 and SMG1 emerged from the well-defined secondary structures, whereas SMG7 appeared from the connecting loops. Comparing the influence of cancer-derived mutations over different CLMS sites revealed that variants in the PPIs for UPF2 or SMG1 have significant structural stability effects. Our data highlights the protein-protein interface of the SMG1, UPF2, and SMG7 genes that can be used for potential therapeutic approaches. Blocking the NMD pathway could enhance the production of neoantigens or internal cancer vaccines, which could provide a platform to design potential peptide-based vaccines.

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

TYK2 is a nonreceptor tyrosine kinase, member of the Janus kinases (JAK), with a central role in several diseases, including cancer. The JAKs\' catalytic domains (KD) are highly conserved, yet the isolated TYK2-KD exhibits unique specificities. In a previous work, using molecular dynamics (MD) simulations of a catalytically impaired TYK2-KD variant (P1104A) we found that this amino acid change of its JAK-characteristic insert (αFG), acts at the dynamics level. Given that structural dynamics is key to the allosteric activation of protein kinases, in this study we applied a long-scale MD simulation and investigated an active TYK2-KD form in the presence of adenosine 5\'-triphosphate and one magnesium ion that represents a dynamic and crucial step of the catalytic cycle, in other protein kinases. Community analysis of the MD trajectory shed light, for the first time, on the dynamic profile and dynamics-driven allosteric communications within the TYK2-KD during activation and revealed that αFG and amino acids P1104, P1105, and I1112 in particular, hold a pivotal role and act synergistically with a dynamically coupled communication network of amino acids serving intra-KD signaling for allosteric regulation of TYK2 activity. Corroborating our findings, most of the identified amino acids are associated with cancer-related missense/splice-site mutations of the Tyk2 gene. We propose that the conformational dynamics at this step of the catalytic cycle, coordinated by αFG, underlie TYK2-unique substrate recognition and account for its distinct specificity. In total, this work adds to knowledge towards an in-depth understanding of TYK2 activation and may be valuable towards a rational design of allosteric TYK2-specific inhibitors.

ROS1 is the largest receptor tyrosine kinase in the human genome. Rearrangements of the ROS1 gene result in oncogenic ROS1 kinase fusion proteins that are currently the only validated biomarkers for targeted therapy with ROS1 TKIs in patients. While numerous somatic missense mutations in ROS1 exist in the cancer genome, their impact on catalytic activity and pathogenic potential is unknown. We interrogated the AACR Genie database and identified 34 missense mutations in the ROS1 tyrosine kinase domain for further analysis. Our experiments revealed that these mutations have varying effects on ROS1 kinase function, ranging from complete loss to significantly increased catalytic activity. Notably, Asn and Gly substitutions at Asp2113 in the ROS1 kinase domain were found to be TKI-sensitive oncogenic variants in cell-based model systems. In vivo experiments showed that ROS1 D2113N induced tumor formation that was sensitive to crizotinib and lorlatinib, FDA-approved ROS1-TKIs. Collectively, these findings highlight the tumorigenic potential of specific point mutations within the ROS1 kinase domain and their potential as therapeutic targets with FDA-approved ROS1-TKIs.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3\' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3\' splice site recognition of long introns, which represent 80% of all human introns.

Tumour suppressor p53 plays a key role in the development of cancer and has therefore been widely studied in recent decades. While it is well known that p53 is biologically active as a tetramer, the tetramerisation mechanism is still not completely understood. p53 is mutated in nearly 50% of cancers, and mutations can alter the oligomeric state of the protein, having an impact on the biological function of the protein and on cell fate decisions. Here, we describe the effects of a number of representative cancer-related mutations on tetramerisation domain (TD) oligomerisation defining a peptide length that permits having a folded and structured domain, thus avoiding the effect of the flanking regions and the net charges at the N- and C-terminus. These peptides have been studied under different experimental conditions. We have applied a variety of techniques, including circular dichroism (CD), native mass spectrometry (MS) and high-field solution NMR. Native MS allows us to detect the native state of complexes maintaining the peptide complexes intact in the gas phase; the secondary and quaternary structures were analysed in solution by NMR, and the oligomeric forms were assigned by diffusion NMR experiments. A significant destabilising effect and a variable monomer population were observed for all the mutants studied.

In this autobiographical article, I reflect on my Swedish background. Then I discuss endogenous DNA alterations and the base excision repair pathway and alternative repair strategies for some unusual DNA lesions. Endogenous DNA damage, such as loss of purine bases and cytosine deamination, is proposed as a major source of cancer-causing mutations.