Drug development in Systemic Lupus Erythematosus (SLE) has been hindered by poor translation from successful preclinical experiments to clinical efficacy. This lack of success has been attributed to the high heterogeneity of SLE patients and to the lack of understanding of disease physiopathology. Modelling approaches could be useful for supporting the identification of targets, biomarkers and patient subpopulations with differential response to drugs. However, the use of traditional quantitative models based on differential equations is not justifiable in a sparse data situation. Boolean networks models are less demanding on the required data to be implemented and can provide insights into the dynamics of biological networks. This methodology allows the integration of all the available knowledge into a single framework to evaluate the behavior of the system under different conditions and test hypotheses about unknown aspects of the disease. In this proof-of-concept study, we explored the potential of a systems pharmacology model based on Boolean networks to support drug development in SLE. We focused the analysis on the antigen presentation by the antigen presenting cells (APC) to the T-cells to evaluate the scope of this methodology in a medium size network before full implementation of the whole SLE pathway. The heterogeneity of SLE patients was replicated using this methodology simulating subjects with distinct pathway alterations. A perturbation analysis of the network coupled with clustering analysis showed potential to identify drug targets, optimal combinatorial regimens and subpopulations of responders and non-responders to drug treatment. We propose this approach as a first step towards the development of more quantitative platforms to address the current challenges in drug development for complex diseases.