Bivalent chromatin as a therapeutic target in cancer: Anin silicopredictive approach for combining epigenetic drugs
Authors: Alarcón Tomás , Sardanyés Josep , Guillamon Antoni , Menendez Javier A. 
Published: October 2020

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.

Author summary

Heterogeneity in cancer cell populations is one of the main engines of resistance to targeted therapies, as it induces nonuniform responses within the population that clears the sensitive subpopulation, whilst leaving unaffected the non-responsive cells. Although this is a well-known fact, few successful approaches have been proposed aimed at both quantifying the variability associated to cell heterogeneity, and characterising strategies that circumvent its drug-resistance inducing effects. Here we present a computational approach that addresses these issues in the particular context of targeting epigenetic regulators (specifically, chromatin modifiers), which have been proposed as therapeutic targets in several types of cancer and also in ageing-related diseases. Our model predicts that the more successful combinations involve modulators of demethylase activity (specifically, KDM5/6 and UTX). By contrast, strategies involving EZH2 activity are predicted to be less effective. Our results support the use of our framework as a platform forin silicodrug 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.