Mathematical modelling of promoter occupancies in MYC-dependent gene regulation

The human MYC proto-oncogene protein (MYC) is a transcription factor that plays a major role in the regulation of cell proliferation. Deregulation of MYC expression is often found in cancer. In the last years, several hypotheses have been proposed to explain cell type specific MYC target gene expression patterns despite genome wide DNA binding of MYC. In a recent publication, a mathematical modelling approach in combination with experimental data demonstrated that differences in MYC-DNA-binding affinity are sufficient to explain distinct promoter occupancies and allow stratification of distinct MYC-regulated biological processes at different MYC concentrations. Here, we extend the analysis of the published mathematical model of DNA-binding behaviour of MYC to demonstrate that the insights gained in the investigation of the human osteosarcoma cell line U2OS can be generalized to other human cell types

Paracrine and autocrine regulation of gene expression by Wnt-inhibitor Dickkopf in wild-type and mutant hepatocytes

BACKGROUND: Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/β-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/β-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells. RESULTS: Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/β-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength. CONCLUSIONS: Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter- and intracellular processes, such as Dkk feedback, diffusion and Wnt/β-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression

Modeling Wnt/β-catenin target gene expression in APC and Wnt gradients under wild type and mutant conditions

The Wnt/β-catenin pathway is involved in the regulation of a multitude of physiological processes by controlling the differential expression of target genes. In certain tissues such as the adult liver, the Wnt/β-catenin pathway can attain different levels of activity due to gradients of Wnt ligands and/or intracellular pathway components like APC. How graded pathway activity is converted into regionally distinct patterns of Wnt/β-catenin target gene expression is largely unknown. Here, we apply a mathematical modeling approach to investigate the impact of different regulatory mechanisms on target gene expression within Wnt or APC concentration gradients. We develop a minimal model of Wnt/beta-catenin signal transduction and combine it with various mechanisms of target gene regulation. In particular, the effects of activation, inhibition, and an incoherent feedforward loop (iFFL) are compared. To specify activation kinetics, we analyze experimental data that quantify the response of β-catenin/TCF reporter constructs to different Wnt concentrations, and demonstrate that the induction of these constructs occurs in a cooperative manner with Hill coefficients between 2 and 5. In summary, our study shows that the combination of specific gene regulatory mechanisms with a time-independent gradient of Wnt or APC is sufficient to generate distinct target gene expression patterns as have been experimentally observed in liver. We find that cooperative gene activation in combination with a TCF feedback can establish sharp borders of target gene expression in Wnt or APC gradients. In contrast, the iFFL renders gene expression independent of gradients of the upstream signaling components. Our subsequent analysis of carcinogenic pathway mutations reveals that their impact on gene expression is determined by the gene regulatory mechanism and the APC concentration of the cell in which the mutation occurs

Different promoter affinities account for specificity in MYC-dependent gene regulation

Enhanced expression of the MYC transcription factor is observed in the majority of tumors. Two seemingly conflicting models have been proposed for its function: one proposes that MYC enhances expression of all genes, while the other model suggests gene-specific regulation. Here, we have explored the hypothesis that specific gene expression profiles arise since promoters differ in affinity for MYC and high-affinity promoters are fully occupied by physiological levels of MYC. We determined cellular MYC levels and used RNA- and ChIP-sequencing to correlate promoter occupancy with gene expression at different concentrations of MYC. Mathematical modeling showed that binding affinities for interactions of MYC with DNA and with core promoter-bound factors, such as WDR5, are sufficient to explain promoter occupancies observed in vivo. Importantly, promoter affinity stratifies different biological processes that are regulated by MYC, explaining why tumor-specific MYC levels induce specific gene expression programs and alter defined biological properties of cells

A quantitative modular modeling approach reveals the effects of different A20 feedback implementations for the NF-κB signaling dynamics

Signaling pathways involve complex molecular interactions and are controled by non-linear regulatory mechanisms. If details of regulatory mechanisms are not fully elucidated, they can be implemented by different, equally reasonable mathematical representations in computational models. The study presented here focusses on NF-κB signaling, which is regulated by negative feedbacks via IκBα and A20. A20 inhibits NF-κB activation indirectly through interference with proteins that transduce the signal from the TNF receptor complex to activate the IκB kinase (IKK) complex. A number of pathway models has been developed implementing the A20 effect in different ways. We here focus on the question how different A20 feedback implementations impact the dynamics of NF-κB. To this end, we develop a modular modeling approach that allows combining previously published A20 modules with a common pathway core module. The resulting models are fitted to a published comprehensive experimental data set and therefore show quantitatively comparable NF-κB dynamics. Based on defined measures for the initial and long-term behavior we analyze the effects of a wide range of changes in the A20 feedback strength, the IκBα feedback strength and the TNFα stimulation strength on NF-κB dynamics. This shows similarities between the models but also model-specific differences. In particular, the A20 feedback strength and the TNFα stimulation strength affect initial and long-term NF-κB concentrations differently in the analyzed models. We validated our model predictions experimentally by varying TNFα concentrations applied to HeLa cells. These time course data indicate that only one of the A20 feedback models appropriately describes the impact of A20 on the NF-κB dynamics in this cell type

N-Myc-induced metabolic rewiring creates novel therapeutic vulnerabilities in neuroblastoma

N-Myc is a transcription factor that is aberrantly expressed in many tumor types and is often correlated with poor patient prognosis. Recently, several lines of evidence pointed to the fact that oncogenic activation of Myc family proteins is concomitant with reprogramming of tumor cells to cope with an enhanced need for metabolites during cell growth. These adaptions are driven by the ability of Myc proteins to act as transcriptional amplifiers in a tissue-of-origin specific manner. Here, we describe the effects of N-Myc overexpression on metabolic reprogramming in neuroblastoma cells. Ectopic expression of N-Myc induced a glycolytic switch that was concomitant with enhanced sensitivity towards 2-deoxyglucose, an inhibitor of glycolysis. Moreover, global metabolic profiling revealed extensive alterations in the cellular metabolome resulting from overexpression of N-Myc. Limited supply with either of the two main carbon sources, glucose or glutamine, resulted in distinct shifts in steady-state metabolite levels and significant changes in glutathione metabolism. Interestingly, interference with glutamine-glutamate conversion preferentially blocked proliferation of N-Myc overexpressing cells, when glutamine levels were reduced. Thus, our study uncovered N-Myc induction and nutrient levels as important metabolic master switches in neuroblastoma cells and identified critical nodes that restrict tumor cell proliferation

Leveraging large language models for decision support in personalized oncology

IMPORTANCE: Clinical interpretation of complex biomarkers for precision oncology currently requires manual investigations of previous studies and databases. Conversational large language models (LLMs) might be beneficial as automated tools for assisting clinical decision-making. OBJECTIVE: To assess performance and define their role using 4 recent LLMs as support tools for precision oncology. DESIGN, SETTING, AND PARTICIPANTS: This diagnostic study examined 10 fictional cases of patients with advanced cancer with genetic alterations. Each case was submitted to 4 different LLMs (ChatGPT, Galactica, Perplexity, and BioMedLM) and 1 expert physician to identify personalized treatment options in 2023. Treatment options were masked and presented to a molecular tumor board (MTB), whose members rated the likelihood of a treatment option coming from an LLM on a scale from 0 to 10 (0, extremely unlikely; 10, extremely likely) and decided whether the treatment option was clinically useful. MAIN OUTCOMES AND MEASURES: Number of treatment options, precision, recall, F1 score of LLMs compared with human experts, recognizability, and usefulness of recommendations. RESULTS: For 10 fictional cancer patients (4 with lung cancer, 6 with other; median [IQR] 3.5 [3.0-4.8] molecular alterations per patient), a median (IQR) number of 4.0 (4.0-4.0) compared with 3.0 (3.0-5.0), 7.5 (4.3-9.8), 11.5 (7.8-13.0), and 13.0 (11.3-21.5) treatment options each was identified by the human expert and 4 LLMs, respectively. When considering the expert as a criterion standard, LLM-proposed treatment options reached F1 scores of 0.04, 0.17, 0.14, and 0.19 across all patients combined. Combining treatment options from different LLMs allowed a precision of 0.29 and a recall of 0.29 for an F1 score of 0.29. LLM-generated treatment options were recognized as AI-generated with a median (IQR) 7.5 (5.3-9.0) points in contrast to 2.0 (1.0-3.0) points for manually annotated cases. A crucial reason for identifying AI-generated treatment options was insufficient accompanying evidence. For each patient, at least 1 LLM generated a treatment option that was considered helpful by MTB members. Two unique useful treatment options (including 1 unique treatment strategy) were identified only by LLM. CONCLUSIONS AND RELEVANCE: In this diagnostic study, treatment options of LLMs in precision oncology did not reach the quality and credibility of human experts; however, they generated helpful ideas that might have complemented established procedures. Considering technological progress, LLMs could play an increasingly important role in assisting with screening and selecting relevant biomedical literature to support evidence-based, personalized treatment decisions

Hunt for new phenomena using large jet multiplicities and missing transverse momentum with ATLAS in 4.7 fb−1 of s√=7TeV proton-proton collisions

Results are presented of a search for new particles decaying to large numbers of jets in association with missing transverse momentum, using 4.7 fb−1 of pp collision data at s√=7TeV collected by the ATLAS experiment at the Large Hadron Collider in 2011. The event selection requires missing transverse momentum, no isolated electrons or muons, and from ≥6 to ≥9 jets. No evidence is found for physics beyond the Standard Model. The results are interpreted in the context of a MSUGRA/CMSSM supersymmetric model, where, for large universal scalar mass m 0, gluino masses smaller than 840 GeV are excluded at the 95% confidence level, extending previously published limits. Within a simplified model containing only a gluino octet and a neutralino, gluino masses smaller than 870 GeV are similarly excluded for neutralino masses below 100 GeV

Measurements of fiducial and differential cross sections for Higgs boson production in the diphoton decay channel at s√=8 TeV with ATLAS

Measurements of fiducial and differential cross sections are presented for Higgs boson production in proton-proton collisions at a centre-of-mass energy of s√=8 TeV. The analysis is performed in the H → γγ decay channel using 20.3 fb−1 of data recorded by the ATLAS experiment at the CERN Large Hadron Collider. The signal is extracted using a fit to the diphoton invariant mass spectrum assuming that the width of the resonance is much smaller than the experimental resolution. The signal yields are corrected for the effects of detector inefficiency and resolution. The pp → H → γγ fiducial cross section is measured to be 43.2 ±9.4(stat.) − 2.9 + 3.2 (syst.) ±1.2(lumi)fb for a Higgs boson of mass 125.4GeV decaying to two isolated photons that have transverse momentum greater than 35% and 25% of the diphoton invariant mass and each with absolute pseudorapidity less than 2.37. Four additional fiducial cross sections and two cross-section limits are presented in phase space regions that test the theoretical modelling of different Higgs boson production mechanisms, or are sensitive to physics beyond the Standard Model. Differential cross sections are also presented, as a function of variables related to the diphoton kinematics and the jet activity produced in the Higgs boson events. The observed spectra are statistically limited but broadly in line with the theoretical expectations