Comparative Network Reconstruction using Mixed Integer Programming

New anti-cancer drugs that specifically target oncogenes involved in signalling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. Addressing this problem requires an understanding of these networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells. Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on incomplete perturbation data, and to identify which edges differ quantitatively between two or more signalling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re-establishing wildtype MAPK signaling, possibly through an alternative RAF-isoform.