The complex differentiation process of the CD4+ T helper lymphocytes shapes the form and the range of the immune response to different antigenic challenges. Few mathematical and computational models have addressed this key phenomenon. We here present a multiscale approach in which two different levels of description, i.e. a gene regulatory network model and an agent-based simulator for cell population dynamics, are integrated into a single immune system model.