UISS, an in silico platform for the human immune system simulation
Agent-Based Models (ABMs) belong to the class of mechanistic models, a family of models that, differently from data-driven models, uses a description of the underlying mechanisms of a given phenomenon to reproduce it. Such a description is usually based on different observational data, previous knowledge and/or hypotheses, and is usually aggregated and rationalized into a conceptual map (i.e., a flow chart and/or a schematic disease model) that reassumes the cascade of events of the phenomenon under investigation. The conceptual map is then translated into mathematical/computational terms and then executed by computers to observe, in silico, the evolution of the phenomenon over time. Besides ABMs, other modeling techniques based on the mechanistic approach can be used. Among these, we recall, for example, ordinary and partial differential equations [40,41,42] and Petri nets [43, 44].
As the name suggests, agent-based models are based on the paradigm of ‘agents’, autonomous entities that behave individually according to established rules. Such entities can be heterogeneous in nature, and are usually represented on a simulation space where they are free to move, interact each-other and change their internal state as a consequence of interactions. From a computer science perspective, agents can be seen as stochastic finite-state machines, capable of assuming a limited number of discrete states. Using ABMs, the global evolution of the phenomena is observed by taking into account the sum of the individual behaviors of all agents, and sometimes unexpected “emergent” behaviors may be observed.
ABMs have been successfully applied in many research fields, from social sciences to ecology, from epidemiology to biology. In the field of immunology, we developed the Universal Immune System Simulator (UISS), an agent-based framework that has been extended through the last decades to simulate the behavior of the immune system response when challenged against many diseases.
In UISS agents are used to describe cells and molecules of the immune system, as well as external actors that can destabilize (i.e., pathogens such as viruses and bacteria) or restore (i.e. prophylactic and therapeutic treatments) the normal health of the host.
One of the main features of UISS is its ability to mimic the adaptive immune response mechanisms. Mammals have in fact developed an advanced immune system machinery capable to specifically recognize pathogens in order to better react against them. This advanced response is based on the ability to exactly recognize foreign proteins (i.e., epitopes) on pathogens surface by means of receptors, through a key-to-lock mechanism. While an explicit implementation would be both unfeasible and partially inaccurate from a computational point of view, in UISS we mimic such a process through the use of binary strings. Binary strings are used for both representing epitopes and immune system cells’ receptors, and the probability that an immune system cell recognizes a pathogen is proportional to the Hamming distance (the number of mismatching bits) between the two strings involved into the interaction. Although this abstraction may seem binding, millions of interactions can be simulated quickly on modern computers, making easier the reproduction of many features of the immune system such as memory, specificity, tolerance and homeostasis. For example, this abstraction demonstrated able to allow the selection of the best adjuvant among a series of candidates for an influenza vaccine when properly coupled with results coming from existing binding prediction tools [32]. This suggests how such an abstraction is able to capture the complexity of the problem.
Besides of receptors, UISS implements many other immune system mechanisms, as thymus selection, haematopoiesis, cell maturation, Hayflick limit, aging, immunological memory, antibody hyper-mutation, bystander effect, cell anergy, antigen processing and presentation [45, 46].
Up to now, UISS in silico platform has been successfully applied to the design and verification of novel treatments for many diseases in both preclinical and clinical environments, including pathologies such as mammary carcinoma [47] and derived lung metastases [48], melanoma [49], atherosclerosis [50], multiple sclerosis [35] and influenza [32].
More recently, UISS has been used as the centerpiece of the StriTuVaD H2020 project with the aim to create an in silico trial for tuberculosis. In this context, observations from virtual patients will be coupled with results from a real clinical trial to obtain an in silico augmented clinical trial, with greater accuracy and more statistical power [34].
SARS-CoV-2 disease model
The SARS-CoV-2 disease model has been implemented in UISS computational framework starting by identifying a question of interest. The question of interest describes the specific question, decision or concern that is being addressed with a computational model. In other words, the question of interest lays out the engineering question that is to be answered (at least in part) through a model. The next step is to define the context of use (CoU), which provides a detailed and complete explanation of how the computational model output will be used to answer the question of interest.
In this specific study, the question of interest is how potential prophylactic or therapeutic vaccines could cure COVID-19, building or stimulating an effective immune response against SARS-CoV-2 virus. UISS must then represent and reproduce the fundamental SARS-CoV-2—immune system competition and dynamics. To this end, we first selected all the players that have a role in the viral infection both at cellular and molecular scale and then we categorized all the interactions among entities that play a relevant role in this biological scenario. Finally compartment assumptions have to be done to let the entities move and interact each other. In our case, we considered the lung compartment that models the main organ target of the virus and the generic lymph node that allows immune system entities to be activated and selected. Figure 1 gives a detailed sketch on the main compartments, entities and interactions.
SARS-CoV-2 first entry is located in the upper respiratory tract. Then it proceeds to bronchial and finally to lungs in which it reaches its main cellular target i.e., the epithelial lung cells (LEP) [2]. The virus is eventually captured by dendritic cells (DC) and macrophages (M).
DC are the main antigen processing cells of the immune system [51] that are able to present the peptides antigen complexed in both major histocompatibility class I and class II (MHC-I and MHC-II, respectively). If a DC encounters the native virus form, it can be able to process it and present its peptides complexed with MHC-II to CD4 T cells for further actions. DC, upon virus activation, release interferon type A and B (IFN-A and IFN-B) and interleukin-12 (IL-12) that are important cytokines in fighting intracellular pathogens. Also, M are able to capture the native form of the SARS-CoV-2 and, if properly activated by pro-inflammatory cytokines, be able to internally destroy it. After their successful activation, macrophages release a pro-inflammatory cytokine that is tumor necrosis factor alpha (TNF-alpha).
A fraction of SARS-CoV-2 viruses reach LEP and through the envelope spike glycoprotein binds to their cellular receptor, angiotensin-converting enzyme 2 (ACE2). Doing that, the viral RNA genome starts to be released into the cytoplasm and is translated into two polyproteins and structural proteins, after which the viral genome begins to replicate inside the cell [52].
Following the flux of the conceptual disease model represented in Fig. 1, after a certain amount of time (that we tuned with available data, as described in the next sections), new copies of the virus are released from the infected LEP that eventually dies. New released copies of functional SARS-CoV-2 infect new cells, spreading further the infection in the lungs.
When a cell is infected by a virus, it can be susceptible of different destinies. One of them is the shutting down of MHC-I expression to avoid immune system recognition from specific CD8 T cells. In this case, a population of innate immunity cells, natural killer cells (NK) may identify them and proceed to kill them through specific actions. The other one is represented by a different MHC-I presentation on the cell surface, as the virus has modified the normal behavior of cell to let the host to make functioning virus copies. In this circumstance, (that we supposed to happen during SARS-CoV-2 infection) cell MHC-I presentation is different from the normal case.
DC are able (through a mechanism known as “nibbling” process [53]) to cross present the antigen complexed with MHC-I proteins to let adaptive immune response to recognize and kill virus infected cells. Activated and antigen presenting cells (both DC and/or M) migrate into the proximal lymph nodes to present their content to adaptive immune cells i.e., T cells and B cells. We implemented nibbling process in a specific UISS interaction in which DC capture Ag from live cells through intimate cell contact, presenting in MHC class I complex to T cells for further actions.
Also, a portion of viruses could eventually migrate to the lymph nodes. Here, B cells can be activated by virus if specific immunoglobulin receptor in B cell surface binds to it. In this context, B cell is activated, and it immediately releases immunoglobulins of M class (IgM) that are the first antibody response that can be measured. Further, APC cells activate CD4 T cells (helper T cells, Th) that under the influence of specific cytokines released before, differentiate into helper T cell type 1 (Th1). Th1 migrate under chemokines gradient to the site of infection. There, they release interferon gamma (IFN-G) that makes macrophages able to destroy captured viral particles and allow them to release IL-12 that promotes immune system activation against the virus. Th1 cells allow the differentiation and the iso-switching B cells into immunoglobulins class G (IgG) producing plasma cells. IgG are specific antibodies that bind against virus receptors, eventually inhibiting its capacity to infect cells. MHC-I/peptides DC presenting cells are also able to activate CD8 cytotoxic T cells (Tc) to destroy SARS-CoV-2 infected cells and then eliminate the reservoir of infection.
Eventually, Tc migrate into the site of infection and recognize and kill infected LEP. Tc killed infected LEP release chemokines and interleukin 1 and 6 (IL-1 and IL-6). IL-1 is the main cytokine that induces several systemic effects in the host, for example fever. IL-6 is a proinflammatory cytokine that can change the severity of COVID-19 disease as reported in very recent literature [54]. Our disease model takes good account of the cytokines storm in the prognosis of the severity of the disease.
Entities (both cellular and molecular) move and diffuse in a simulation space represented as a L X L lattice (L is set depending on the dimension of the compartment one intends to reproduce), with periodic boundary conditions. There is no correlation between entities residing on different sites at a fixed time as the interactions among cells and molecules take place within a lattice-site in a single time step.
All entities are allowed to move with a uniform probability between neighboring lattices in the grid and with an equal diffusion coefficient (Brownian motion).