Bivalent chromatin as a therapeutic target in cancer: An in silico predictive approach for combining epigenetic drugs.

Tomás Alarcón,Josep Sardanyés,Javier A Menendez,Antoni Guillamon,Ilya Ioshikhes

doi:10.1371/journal.pcbi.1008408

Abstract

Tumour cell heterogeneity is a major barrier for efficient design of targeted anti-cancer therapies. A diverse distribution of phenotypically distinct tumour-cell subpopulations prior to drug treatment predisposes to non-uniform responses, leading to the elimination of sensitive cancer cells whilst leaving resistant subpopulations unharmed. Few strategies have been proposed for quantifying the variability associated to individual cancer-cell heterogeneity and minimizing its undesirable impact on clinical outcomes. Here, we report a computational approach that allows the rational design of combinatorial therapies involving epigenetic drugs against chromatin modifiers. We have formulated a stochastic model of a bivalent transcription factor that allows us to characterise three different qualitative behaviours, namely: bistable, high- and low-gene expression. Comparison between analytical results and experimental data determined that the so-called bistable and high-gene expression behaviours can be identified with undifferentiated and differentiated cell types, respectively. Since undifferentiated cells with an aberrant self-renewing potential might exhibit a cancer/metastasis-initiating phenotype, we analysed the efficiency of combining epigenetic drugs against the background of heterogeneity within the bistable sub-ensemble. Whereas single-targeted approaches mostly failed to circumvent the therapeutic problems represented by tumour heterogeneity, combinatorial strategies fared much better. Specifically, the more successful combinations were predicted to involve modulators of the histone H3K4 and H3K27 demethylases KDM5 and KDM6A/UTX. Those strategies involving the H3K4 and H3K27 methyltransferases MLL2 and EZH2, however, were predicted to be less effective. Our theoretical framework provides a coherent basis for the development of an in silico platform capable of identifying the epigenetic drugs combinations best-suited to therapeutically manage non-uniform responses of heterogenous cancer cell populations.

Highlights

Heterogeneity is a primary cause of de novo and acquired resistance to targeted treatments in cancer therapeutics [1,2,3]
We present a computational approach that addresses these issues in the particular context of targeting epigenetic regulators, which have been proposed as therapeutic targets in several types of cancer and in ageing-related diseases
Our results support the use of our framework as a platform for in silico drug trials, as it accounts for non-homogeneous response of cell populations to drugs as well as identifying which subpopulations are more likely to respond to specific strategies

Summary

Introduction

Heterogeneity is a primary cause of de novo and acquired resistance to targeted treatments in cancer therapeutics [1,2,3]. Non-genetic heterogeneity has two broad sources, epigenetic and stochastic, the latter being due to intrinsic factors such as noise in gene expression [10,11,12,13] and asymmetric cell division [14, 15]. These sources of heterogeneity can produce phenotype diversity even in genetically identical cells [8, 9, 16]. The analysis of the consequences of non-genetic cell-to-cell variability in the cellular response to drugs and its potential impact for the treatment of human diseases including cancer has become a crucial issue to understanding drug resistance phenomena and developing new targeted agents [1, 3, 17]

Methods

Results

Discussion

Conclusion