A model of estrogen-related gene expression reveals non-linear effects in transcriptional response to tamoxifen

Estrogen receptors alpha (ER) are implicated in many types of female cancers, and are the common target for anti-cancer therapy using selective estrogen receptor modulators (SERMs, such as tamoxifen). However, cell-type specific and patient-to-patient variability in response to SERMs (from suppression to stimulation of cancer growth), as well as frequent emergence of drug resistance, represents a serious problem. The molecular processes behind mixed effects of SERMs remain poorly understood, and this strongly motivates application of systems approaches. In this work, we aimed to establish a mathematical model of ER-dependent gene expression to explore potential mechanisms underlying the variable actions of SERMs. We developed an equilibrium model of ER binding with 17β-estradiol, tamoxifen and DNA, and linked it to a simple ODE model of ER-induced gene expression. The model was parameterised on the broad range of literature available experimental data, and provided a plausible mechanistic explanation for the dual agonism/antagonism action of tamoxifen in the reference cell line used for model calibration. To extend our conclusions to other cell types we ran global sensitivity analysis and explored model behaviour in the wide range of biologically plausible parameter values, including those found in cancer cells. Our findings suggest that transcriptional response to tamoxifen is controlled in a complex non-linear way by several key parameters, including ER expression level, hormone concentration, amount of ER-responsive genes and the capacity of ER-tamoxifen complexes to stimulate transcription (e.g. by recruiting co-regulators of transcription). The model revealed non-monotonic dependence of ER-induced transcriptional response on the expression level of ER, that was confirmed experimentally in four variants of the MCF-7 breast cancer cell line. We established a minimal mechanistic model of ER-dependent gene expression, that predicts complex non-linear effects in transcriptional response to tamoxifen in the broad range of biologically plausible parameter values. Our findings suggest that the outcome of a SERM's action is defined by several key components of cellular micro-environment, that may contribute to cell-type-specific effects of SERMs and justify the need for the development of combinatorial biomarkers for more accurate prediction of the efficacy of SERMs in specific cell types.