Abstract
In this model we use a dynamic and multistable Boolean regulatory network to provide a mechanistic explanation of the lymphopenia and dysregulation of CD4+ T cell subsets in COVID-19 and provide therapeutic targets. Using a previous model, the cytokine micro-environments found in mild, moderate, and severe COVID-19 with and without TGF-β and IL-10 was we simulated. It shows that as the severity of the disease increases, the number of antiviral Th1 cells decreases, while the the number of Th1-like regulatory and exhausted cells and the proportion between Th1 and Th1R cells increases. The addition of the regulatory cytokines TFG-β and IL-10 makes the Th1 attractor unstable and favors the Th17 and regulatory subsets. This is associated with the contradictory signals in the micro-environment that activate SOCS proteins that block the signaling pathways. Furthermore, it determined four possible therapeutic targets that increase the Th1 compartment in severe COVID-19: the activation of the IFN-γ pathway, or the inhibition of TGF-β or IL-10 pathways or SOCS1 protein; from these, inhibiting SOCS1 has the lowest number of predicted collateral effects. Finally, a tool is provided that allows simulations of specific cytokine environments and predictions of CD4 T cell subsets and possible interventions, as well as associated secondary effects.