Entry into mitosis is driven by the phosphorylation of thousands of substrates, under the master control of Cdk1. During entry into mitosis, Cdk1, in collaboration with MASTL kinase, represses the ...activity of the major mitotic protein phosphatases, PP1 and PP2A, thereby ensuring mitotic substrates remain phosphorylated. For cells to complete and exit mitosis, these phosphorylation events must be removed, and hence, phosphatase activity must be reactivated. This reactivation of phosphatase activity presumably requires the inhibition of MASTL; however, it is not currently understood what deactivates MASTL and how this is achieved. In this study, we identified that PP1 is associated with, and capable of partially dephosphorylating and deactivating, MASTL during mitotic exit. Using mathematical modelling, we were able to confirm that deactivation of MASTL is essential for mitotic exit. Furthermore, small decreases in Cdk1 activity during metaphase are sufficient to initiate the reactivation of PP1, which in turn partially deactivates MASTL to release inhibition of PP2A and, hence, create a feedback loop. This feedback loop drives complete deactivation of MASTL, ensuring a strong switch-like activation of phosphatase activity during mitotic exit.
Wnt signalling is involved in the formation, metastasis and relapse of a wide array of cancers. However, there is ongoing debate as to whether activation or inhibition of the pathway holds the most ...promise as a therapeutic treatment for cancer, with conflicting evidence from a variety of tumour types. We show that Wnt/β-catenin signalling is a bi-directional vulnerability of neuroblastoma, malignant melanoma and colorectal cancer, with hyper-activation or repression of the pathway both representing a promising therapeutic strategy, even within the same cancer type. Hyper-activation directs cancer cells to undergo apoptosis, even in cells oncogenically driven by β-catenin. Wnt inhibition blocks proliferation of cancer cells and promotes neuroblastoma differentiation. Wnt and retinoic acid co-treatments synergise, representing a promising combination treatment for MYCN-amplified neuroblastoma. Additionally, we report novel cross-talks between MYCN and β-catenin signalling, which repress normal β-catenin mediated transcriptional regulation. A β-catenin target gene signature could predict patient outcome, as could the expression level of its DNA binding partners, the TCF/LEFs. This β-catenin signature provides a tool to identify neuroblastoma patients likely to benefit from Wnt-directed therapy. Taken together, we show that Wnt/β-catenin signalling is a bi-directional vulnerability of a number of cancer entities, and potentially a more broadly conserved feature of malignant cells.
Despite intensive study, many mysteries remain about the MYCN oncogene's functions. Here we focus on MYCN's role in neuroblastoma, the most common extracranial childhood cancer. MYCN gene ...amplification occurs in 20% of cases, but other recurrent somatic mutations are rare. This scarcity of tractable targets has hampered efforts to develop new therapeutic options. We employed a multi-level omics approach to examine MYCN functioning and identify novel therapeutic targets for this largely un-druggable oncogene. We used systems medicine based computational network reconstruction and analysis to integrate a range of omic techniques: sequencing-based transcriptomics, genome-wide chromatin immunoprecipitation, siRNA screening and interaction proteomics, revealing that MYCN controls highly connected networks, with MYCN primarily supressing the activity of network components. MYCN's oncogenic functions are likely independent of its classical heterodimerisation partner, MAX. In particular, MYCN controls its own protein interaction network by transcriptionally regulating its binding partners.Our network-based approach identified vulnerable therapeutically targetable nodes that function as critical regulators or effectors of MYCN in neuroblastoma. These were validated by siRNA knockdown screens, functional studies and patient data. We identified β-estradiol and MAPK/ERK as having functional cross-talk with MYCN and being novel targetable vulnerabilities of MYCN-amplified neuroblastoma. These results reveal surprising differences between the functioning of endogenous, overexpressed and amplified MYCN, and rationalise how different MYCN dosages can orchestrate cell fate decisions and cancerous outcomes. Importantly, this work describes a systems-level approach to systematically uncovering network based vulnerabilities and therapeutic targets for multifactorial diseases by integrating disparate omic data types.
Early-stage lung cancer is crucial clinically due to its insidious nature and rapid progression. Most of the prediction models designed to predict tumour recurrence in the early stage of lung cancer ...rely on the clinical or medical history of the patient. However, their performance could likely be improved if the input patient data contained genomic information. Unfortunately, such data is not always collected. This is the main motivation of our work, in which we have imputed and integrated specific type of genomic data with clinical data to increase the accuracy of machine learning models for prediction of relapse in early-stage, non-small cell lung cancer patients. Using a publicly available TCGA lung adenocarcinoma cohort of 501 patients, their aneuploidy scores were imputed into similar records in the Spanish Lung Cancer Group (SLCG) data, more specifically a cohort of 1348 early-stage patients. First, the tumor recurrence in those patients was predicted without the imputed aneuploidy scores. Then, the SLCG data were enriched with the aneuploidy scores imputed from TCGA. This integrative approach improved the prediction of the relapse risk, achieving area under the precision-recall curve (PR-AUC) score of 0.74, and area under the ROC (ROC-AUC) score of 0.79. Using the prediction explanation model SHAP (SHapley Additive exPlanations), we further explained the predictions performed by the machine learning model. We conclude that our explainable predictive model is a promising tool for oncologists that addresses an unmet clinical need of post-treatment patient stratification based on the relapse risk, while also improving the predictive power by incorporating proxy genomic data not available for the actual specific patients.
Protein degradation via ubiquitination is a major proteolytic mechanism in cells. Once a protein is destined for degradation, it is tagged by multiple ubiquitin (Ub) molecules. The synthesized ...polyubiquitin chains can be recognized by the 26S proteosome where proteins are degraded. These chains form through multiple ubiquitination cycles that are similar to multi-site phosphorylation cycles. As kinases and phosphatases, two opposing enzymes (E3 ligases and deubiquitinases DUBs) catalyze (de)ubiquitination cycles. Although multi-ubiquitination cycles are fundamental mechanisms of controlling protein concentrations within a cell, their dynamics have never been explored. Here, we fill this knowledge gap. We show that under permissive physiological conditions, the formation of polyubiquitin chain of length greater than two and subsequent degradation of the ubiquitinated protein, which is balanced by protein synthesis, can display bistable, switch-like responses. Interestingly, the occurrence of bistability becomes pronounced, as the chain grows, giving rise to "all-or-none" regulation at the protein levels. We give predictions of protein distributions under bistable regime awaiting experimental verification. Importantly, we show for the first time that sustained oscillations can robustly arise in the process of formation of ubiquitin chain, largely due to the degradation of the target protein. This new feature is opposite to the properties of multi-site phosphorylation cycles, which are incapable of generating oscillation if the total abundance of interconverted protein forms is conserved. We derive structural and kinetic constraints for the emergence of oscillations, indicating that a competition between different substrate forms and the E3 and DUB is critical for oscillation. Our work provides the first detailed elucidation of the dynamical features brought about by different molecular setups of the polyubiquitin chain assembly process responsible for protein degradation.
A network of the Rho family GTPases, which cycle between inactive GDP-bound and active GTP-bound states, controls key cellular processes, including proliferation and migration. Activating and ...deactivating GTPase transitions are controlled by guanine nucleotide exchange factors (GEFs), GTPase activating proteins (GAPs) and GDP dissociation inhibitors (GDIs) that sequester GTPases from the membrane to the cytoplasm. Here we show that a cascade of two Rho family GTPases, RhoA and Rac1, regulated by RhoGDI1, exhibits distinct modes of the dynamic behavior, including abrupt, bistable switches, excitable overshoot transitions and oscillations. The RhoGDI1 abundance and signal-induced changes in the RhoGDI1 affinity for GTPases control these different dynamics, enabling transitions from a single stable steady state to bistability, to excitable pulses and to sustained oscillations of GTPase activities. These RhoGDI1-controlled dynamic modes of RhoA and Rac1 activities form the basis of cell migration behaviors, including protrusion-retraction cycles at the leading edge of migrating cells.
The Morris water maze is an experimental procedure in which animals learn to escape swimming in a pool using environmental cues. Despite its success in neuroscience and psychology for studying ...spatial learning and memory, the exact mnemonic and navigational demands of the task are not well understood. Here, we provide a mathematical model of rat swimming dynamics on a behavioural level. The model consists of a random walk, a heading change and a feedback control component in which learning is reflected in parameter changes of the feedback mechanism. The simplicity of the model renders it accessible and useful for analysis of experiments in which swimming paths are recorded. Here, we used the model to analyse an experiment in which rats were trained to find the platform with either three or one extramaze cue. Results indicate that the 3-cues group employs stronger feedback relying only on the actual visual input, whereas the 1-cue group employs weaker feedback relying to some extent on memory. Because the model parameters are linked to neurological processes, identifying different parameter values suggests the activation of different neuronal pathways.
Early-stage lung cancer is crucial clinically due to its insidious nature and rapid progression. Most of the prediction models designed to predict tumour recurrence in the early stage of lung cancer ...rely on the clinical or medical history of the patient. However, their performance could likely be improved if the input patient data contained genomic information. Unfortunately, such data is not always collected. This is the main motivation of our work, in which we have imputed and integrated specific type of genomic data with clinical data to increase the accuracy of machine learning models for prediction of relapse in early-stage, non-small cell lung cancer patients. Using a publicly available TCGA lung adenocarcinoma cohort of 501 patients, their aneuploidy scores were imputed into similar records in the Spanish Lung Cancer Group (SLCG) data, more specifically a cohort of 1348 early-stage patients. First, the tumor recurrence in those patients was predicted without the imputed aneuploidy scores. Then, the SLCG data were enriched with the aneuploidy scores imputed from TCGA. This integrative approach improved the prediction of the relapse risk, achieving area under the precision-recall curve (PR-AUC) score of 0.74, and area under the ROC (ROC-AUC) score of 0.79. Using the prediction explanation model SHAP (SHapley Additive exPlanations), we further explained the predictions performed by the machine learning model. We conclude that our explainable predictive model is a promising tool for oncologists that addresses an unmet clinical need of post-treatment patient stratification based on the relapse risk, while also improving the predictive power by incorporating proxy genomic data not available for the actual specific patients.
We developed a multiscale model to bridge neuropeptide receptor-activated signaling pathway activity with membrane electrophysiology. Typically, the neuromodulation of biochemical signaling and ...biophysics have been investigated separately in modeling studies. We studied the effects of Angiotensin II (AngII) on neuronal excitability changes mediated by signaling dynamics and downstream phosphorylation of ion channels. Experiments have shown that AngII binding to the AngII receptor type-1 elicits baseline-dependent regulation of cytosolic Ca2+ signaling. Our model simulations revealed a baseline Ca2+-dependent response to AngII receptor type-1 activation by AngII. Consistent with experimental observations, AngII evoked a rise in Ca2+ when starting at a low baseline Ca2+ level, and a decrease in Ca2+ when starting at a higher baseline. Our analysis predicted that the kinetics of Ca2+ transport into the endoplasmic reticulum play a critical role in shaping the Ca2+ response. The Ca2+ baseline also influenced the AngII-induced excitability changes such that lower Ca2+ levels were associated with a larger firing rate increase. We examined the relative contributions of signaling kinases protein kinase C and Ca2+/Calmodulin-dependent protein kinase II to AngII-mediated excitability changes by simulating activity blockade individually and in combination. We found that protein kinase C selectively controlled firing rate adaptation whereas Ca2+/Calmodulin-dependent protein kinase II induced a delayed effect on the firing rate increase. We tested whether signaling kinetics were necessary for the dynamic effects of AngII on excitability by simulating three scenarios of AngII-mediated KDR channel phosphorylation: (1), an increased steady state; (2), a step-change increase; and (3), dynamic modulation. Our results revealed that the kinetics emerging from neuromodulatory activation of the signaling network were required to account for the dynamical changes in excitability. In summary, our integrated multiscale model provides, to our knowledge, a new approach for quantitative investigation of neuromodulatory effects on signaling and electrophysiology.