Computational modeling reveals MAP3K8 as mediator of resistance to vemurafenib in thyroid cancer stem cells

Motivation: Val600Glu (V600E) mutation is the most common BRAF mutation detected in thyroid cancer. Hence, recent research efforts have been performed trying to explore several inhibitors of the V600E mutation-containing BRAF kinase as potential therapeutic options in thyroid cancer refractory to standard interventions. Among them, vemurafenib is a selective BRAF inhibitor approved by Food and Drug Administration for clinical practice. Unfortunately, vemurafenib often displays limited efficacy in poorly differentiated and anaplastic thyroid carcinomas probably because of intrinsic and/or acquired resistance mechanisms. In this view, cancer stem cells (CSCs) may represent a possible mechanism of resistance to vemurafenib, due to their self-renewal and chemo resistance properties. Results: We present a computational framework to suggest new potential targets to overcome drug resistance. It has been validated with an in vitro model based upon a spheroid-forming method able to isolate thyroid CSCs that may mimic resistance to vemurafenib. Indeed, vemurafenib did not inhibit cell proliferation of BRAF V600E thyroid CSCs, but rather stimulated cell proliferation along with a paradoxical over-activation of ERK and AKT pathways. The computational model identified a fundamental role of mitogen-activated protein kinase 8 (MAP3K8), a serine/threonine kinase expressed in thyroid CSCs, in mediating this drug resistance. To confirm model prediction, we set a suitable in vitro experiment revealing that the treatment with MAP3K8 inhibitor restored the effect of vemurafenib in terms of both DNA fragmentation and poly (ADP-ribose) polymerase cleavage (apoptosis) in thyroid CSCs. Moreover, MAP3K8 expression levels may be a useful marker to predict the response to vemurafenib.