Biological network models offer a framework for understanding disease by describing the relationships between the mechanisms involved in the regulation of biological processes. Crowdsourcing can efficiently gather feedback from a wide audience with varying expertise. In the Network Verification Challenge, scientists verified and enhanced a set of 46 biological networks relevant to lung and chronic obstructive pulmonary disease. The networks were built using Biological Expression Language and contain detailed information for each node and edge, including supporting evidence from the literature. Network scoring of public transcriptomics data inferred perturbation of a subset of mechanisms and networks that matched the measured outcomes. These results, based on a computable network approach, can be used to identify novel mechanisms activated in disease, quantitatively compare different treatments and time points, and allow for assessment of data with low signal. These networks are periodically verified by the crowd to maintain an up-to-date suite of networks for toxicology and drug discovery applications.

Epithelial to mesenchymal transition (EMT) is an important event during development and cancer metastasis. There is limited understanding of the metabolic alterations that give rise to and take place during EMT. Dysregulation of signalling pathways that impact metabolism, including epidermal growth factor receptor (EGFR), are however a hallmark of EMT and metastasis. In this study, we report the investigation into EGFR signalling and metabolic crosstalk of EMT through constraint-based modelling and analysis of the breast epithelial EMT cell model D492 and its mesenchymal counterpart D492M. We built an EGFR signalling network for EMT based on stoichiometric coefficients and constrained the network with gene expression data to build epithelial (EGFR_E) and mesenchymal (EGFR_M) networks. Metabolic alterations arising from differential expression of EGFR genes was derived from a literature review of AKT regulated metabolic genes. Signaling flux differences between EGFR_E and EGFR_M models subsequently allowed metabolism in D492 and D492M cells to be assessed. Higher flux within AKT pathway in the D492 cells compared to D492M suggested higher glycolytic activity in D492 that we confirmed experimentally through measurements of glucose uptake and lactate secretion rates. The signaling genes from the AKT, RAS/MAPK and CaM pathways were predicted to revert D492M to D492 phenotype. Follow-up analysis of EGFR signaling metabolic crosstalk in three additional breast epithelial cell lines highlighted variability in in vitro cell models of EMT. This study shows that the metabolic phenotype may be predicted by in silico analyses of gene expression data of EGFR signaling genes, but this phenomenon is cell-specific and does not follow a simple trend.

The metabolism of gluconate is well characterized in prokaryotes where it is known to be degraded following phosphorylation by gluconokinase. Less is known of gluconate metabolism in humans. Human gluconokinase activity was recently identified proposing questions about the metabolic role of gluconate in humans. Here we report the recombinant expression, purification and biochemical characterization of isoform I of human gluconokinase alongside substrate specificity and kinetic assays of the enzyme catalyzed reaction. The enzyme, shown to be a dimer, had ATP dependent phosphorylation activity and strict specificity towards gluconate out of 122 substrates tested. In order to evaluate the metabolic impact of gluconate in humans we modeled gluconate metabolism using steady state metabolic network analysis. The results indicate that significant metabolic flux changes in anabolic pathways linked to the hexose monophosphate shunt (HMS) are induced through a small increase in gluconate concentration. We argue that the enzyme takes part in a context specific carbon flux route into the HMS that, in humans, remains incompletely explored. Apart from the biochemical description of human gluconokinase, the results highlight that little is known of the mechanism of gluconate metabolism in humans despite its widespread use in medicine and consumer products.

Oncogene-induced senescence provides a barrier against malignant transformation. However, it can also promote cancer through the secretion of a plethora of factors released by senescent cells, called the senescence associated secretory phenotype (SASP). We have previously shown that in proliferating cells, nuclear lncRNA MIR31HG inhibits p16/CDKN2A expression through interaction with polycomb repressor complexes and that during BRAF-induced senescence, MIR31HG is overexpressed and translocates to the cytoplasm. Here, we show that MIR31HG regulates the expression and secretion of a subset of SASP components during BRAF-induced senescence. The SASP secreted from senescent cells depleted for MIR31HG fails to induce paracrine invasion without affecting the growth inhibitory effect. Mechanistically, MIR31HG interacts with YBX1 facilitating its phosphorylation at serine 102 (p-YBX1S102) by the kinase RSK. p-YBX1S102 induces IL1A translation which activates the transcription of the other SASP mRNAs. Our results suggest a dual role for MIR31HG in senescence depending on its localization and points to the lncRNA as a potential therapeutic target in the treatment of senescence-related pathologies.

Overexpression of epidermal growth factor receptor (EGFR), and downstream pathway activation appears to be a common oncogenic driver in the majority of head and neck squamous cell cancers (HNSCCs); yet targeting EGFR for the treatment of HNSCC has met with limited success. Apart from the anti-EGFR antibody cetuximab, no small molecule EGFR/tyrosine kinase inhibitors (TKIs) have progressed to routine clinical use. The aim of this study was to determine factors contributing to the lack of response to TKIs and identify alternative therapeutic vulnerabilities.

Cullin3 (Cul3)-based ubiquitin E3 ligase complexes catalyze the transfer of ubiquitin from an E2 enzyme to target substrate proteins. In these assemblies, the C-terminal region of Cul3 binds Rbx1/E2-ubiquitin, while the N-terminal region interacts with various BTB (bric-à-brac, tramtrack, broad complex) domain proteins that serve as substrate adaptors. Previous crystal structures of the homodimeric BTB proteins KLHL3, KLHL11 and SPOP in complex with the N-terminal domain of Cul3 revealed the features required for Cul3 recognition in these proteins. A second class of BTB-domain-containing proteins, the KCTD proteins, is also Cul3 substrate adaptors, but these do not share many of the previously identified determinants for Cul3 binding. We report the pentameric crystal structures of the KCTD1 and KCTD9 BTB domains and identify plasticity in the KCTD1 rings. We find that the KCTD proteins 5, 6, 9 and 17 bind to Cul3 with high affinity, while the KCTD proteins 1 and 16 do not have detectable binding. Finally, we confirm the 5:5 assembly of KCTD9/Cul3 complexes by cryo-electron microscopy and provide a molecular rationale for BTB-mediated Cul3 binding specificity in the KCTD family.

Despite the established role of EGFR tyrosine kinase inhibitors (TKIs) in EGFR-mutated NSCLC, drug resistance inevitably ensues, with a paucity of treatment options especially in EGFR T790M-negative resistance.