Abstract
As immune sentinel cells, macrophages are required to respond specifically to diverse immune threats and initiate appropriate immune responses. This stimulus-response specificity (SRS) is in part encoded in the signaling dynamics of the NFκB transcription factor. While experimental stimulus-response studies have typically focused on single defined ligands, in physiological contexts cells are exposed to multi-ligand mixtures. It remains unclear how macrophages combine multi-ligand information and whether they are able to maintain SRS in such complex exposure conditions. Here, we leveraged an established mathematical model that captures the heterogeneous single-cell NFκB responses of macrophage populations to extend experimental studies with systematic simulations of complex mixtures containing up to five ligands. Live-cell microscopy experiments for some conditions validated model predictions but revealed a discrepancy when TLR3 and TLR9 are stimulated. Refining the model suggested that the observed but unexpected ligand antagonism arises from a limited capacity for endosomal transport which is required for responses to CpG and pIC. With the updated model, we systematically analyzed SRS across all combinatorial-ligand conditions and employed three ways of quantifying SRS involving trajectory decomposition into informative trajectory features or machine learning. Our findings show that macrophages most effectively distinguish single-ligand stimuli, and distinguishability declines as more ligands are combined. However, even in complex combinatorial conditions, macrophages still maintain statistically significant distinguishability. These results indicate a robustness of innate immune response specificity: even in the context of complex exposure conditions, the NFκB temporal signaling code of macrophages can still classify immune threats to direct an appropriate response. SIGNIFICANCE ≤120: Macrophages sense diverse pathogens within complex environments and respond appropriately. Experimental studies have found that the NFκB pathway responds with stimulus-specific dynamics when macrophages are exposed to single ligand stimuli. However, it remains unclear complex contexts might erode this stimulus-specificity. Here we systematically studies NFκB responses using a mathematical model that provides simulations of the heterogeneous population of single cell responses. We show that although the model is parameterized to single ligand data it can predict the responses to multi-ligand mixtures. Indeed, model validation uncovered signaling antagonism between two ligands and the underlying mechanism. Importantly, we found that NFκB signaling dynamics distinguish ligands within multi-ligand mixtures indicating a robustness of the NFκB temporal code that was not previously appreciated.