Personalized logical models to investigate cancer response to BRAF treatments in melanomas and colorectal cancers

Abstract

The study of response to cancer treatments has benefited greatly from the contribution of different omics data but their interpretation is sometimes difficult. Some mathematical models based on prior biological knowledge of signalling pathways, facilitate this interpretation but often require fitting of their parameters using perturbation data. We propose a more qualitative mechanistic approach, based on logical formalism and on the sole mapping and interpretation of omics data, and able to recover differences in sensitivity to gene inhibition without model training. This approach is showcased by the study of BRAF inhibition in patients with melanomas and colorectal cancers who experience significant differences in sensitivity despite similar omics profiles. We first gather information from literature and build a logical model summarizing the regulatory network of the mitogen-activated protein kinase (MAPK) pathway surrounding BRAF, with factors involved in the BRAF inhibition resistance mechanisms. The relevance of this model is verified by automatically assessing that it qualitatively reproduces response or resistance behaviours identified in the literature. Data from over 100 melanoma and colorectal cancer cell lines are then used to validate the models ability to explain differences in sensitivity. This generic model is transformed into personalized cell line-specific logical models by integrating the omics information of the cell lines as constraints of the model. The use of mutations alone allows personalized models to correlate significantly with experimental sensitivities to BRAF inhibition, both from drug and CRISPR targeting, and even better with the joint use of mutations and RNA, supporting multi-omics mechanistic models. A comparison of these untrained models with learning approaches highlights similarities in interpretation and complementarity depending on the size of the datasets. This parsimonious pipeline, which can easily be extended to other biological questions, makes it possible to explore the mechanistic causes of the response to treatment, on an individualized basis. Author summaryWe constructed a logical model to study, from a dynamical perspective, the differences between melanomas and colorectal cancers that share the same BRAF mutations but exhibit different sensitivities to anti-BRAF treatments. The model was built from the literature and completed from existing pathway databases. The model encompasses the key proteins of the MAPK pathway and was made specific to each cancer cell line (100 melanoma and colorectal cell lines from public database) using available omics data, including mutations and RNAseq data. It can simulate the effect of drugs and show high correlation with experimental results. Moreover, the structure of the network confirms both the importance of the reactivation of the MAPK pathway through CRAF and the involvement of PI3K/AKT pathway in the mechanisms of resistance to BRAF inhibition. The study shows that, because of the low number of samples, the mechanistic approach that we propose provides different insights than powerful standard machine learning methodologies would, showing the complementarity between the two approaches. An important aspect to mention is that the mechanistic approach presented here does not rely on training datasets but directly interprets and maps data on the model to simulate drug responses.