Skip to main content

In silico trial to test COVID-19 candidate vaccines: a case study with UISS platform



SARS-CoV-2 is a severe respiratory infection that infects humans. Its outburst entitled it as a pandemic emergence. To get a grip on this outbreak, specific preventive and therapeutic interventions are urgently needed. It must be said that, until now, there are no existing vaccines for coronaviruses. To promptly and rapidly respond to pandemic events, the application of in silico trials can be used for designing and testing medicines against SARS-CoV-2 and speed-up the vaccine discovery pipeline, predicting any therapeutic failure and minimizing undesired effects.


We present an in silico platform that showed to be in very good agreement with the latest literature in predicting SARS-CoV-2 dynamics and related immune system host response. Moreover, it has been used to predict the outcome of one of the latest suggested approach to design an effective vaccine, based on monoclonal antibody. Universal Immune System Simulator (UISS) in silico platform is potentially ready to be used as an in silico trial platform to predict the outcome of vaccination strategy against SARS-CoV-2.


In silico trials are showing to be powerful weapons in predicting immune responses of potential candidate vaccines. Here, UISS has been extended to be used as an in silico trial platform to speed-up and drive the discovery pipeline of vaccine against SARS-CoV-2.


As the epicenter of Coronavirus disease 2019 (COVID-19) and emerging severe acute respiratory syndrome (SARS) caused by novel Coronavirus (2019-nCoV) spread is making its way across the world, global healthcare finds itself facing tremendous challenges. According to the World Health Organization (WHO) situation report (91st), updated on 20 April 2020, there have been globally 72,846 confirmed cases of 2019-nCoV and 5296 cases of death caused by the virus itself [1].

2019-nCoV (also referred to as SARS-CoV-2 or HCoV-19) [2], is the seventh coronavirus known to infect humans along with SARS-CoV, MERS-CoV, HKU1, NL63, OC43 and 229E [3]. While these last four coronaviruses are associated with mild symptoms, SARS-CoV, MERS-CoV and SARS-CoV-2 can cause severe acute respiratory syndrome [4], especially in elderlies, of which men, and those individuals with comorbidities and immunocompromised conditions [5]). Although it is similar to SARS-CoV, SARS-CoV-2 has an improved ability for pathogenicity [6]. In particular, latest evidences during the ongoing pandemic reveal that patients affected by SARS-CoV-2 can progress their clinical picture from fever, cough, ageusia and anosmia, sore throat, breathlessness, fatigue, or malaise to pneumonia, acute respiratory distress syndrome (ARDS) and multi organ dysfunction illness [7]. Significantly, in most critically ill patients, SARS-CoV-2 infection is also associated with a severe clinical inflammatory picture based on a serious cytokine storm that is mainly characterized by elevated plasma concentrations of interleukins 6 (IL-6) [8]. In this scenario, it seems that IL-6 owns an important driving role on the cytokine storm, leading to lung damage and reduced survival [9].

Recently, a growing body of evidence has demonstrated a plethora of symptoms related to COVID-19 infection, ranging from cardiovascular to neurological clinical manifestations, and a different severity in young and adult patients as well as in fragile patients, including diabetic, cancer and immunodeficient patients [10,11,12].

To get a grip on this outbreak and flatten the curve of infection, a specific therapeutic intervention to prevent the severity of the disease is urgently needed to reduce morbidity and mortality because, until now, there are no existing vaccines for coronaviruses.

The ideal profile for a targeted SARS-CoV-2 vaccine must address the need of vaccinating human population, with particular regard of those individuals classified as at high risk, comprising, for example, frontline healthcare workers, individuals over the age of 60 and those that show debilitating chronic diseases.

Recently, specific findings about the genome sequencing of SARS-CoV-2 in different countries where cases of infection were registered, revealed its relative intrinsic genomic variability, its virus dynamics and the related host response mechanisms, unveiling interesting knowledge useful for the formulation of innovative strategies for preventive vaccination.

Specifically, SARS-CoV-2 sequencing along with its relative intrinsic genomic variability [13], the presence of minority variants generated during SARS-CoV-2 replication [14], the involved cellular factors that favors SARS-CoV-2 cell entry [15], the timing in which viral load peaks (during the first week of illness), its gradual decline (over the second week) and the increasing of both IgG and IgM antibodies (around day 10 after symptom onset) represent some of the relevant insights so far delineated and considered by research community about SARS-CoV-2 virus [16].

Even though these findings are having several practice consequences and suggest valuable conclusions, SARS-CoV-2 dynamics has not been yet fully understood. Information about which parts of SARS-CoV-2 sequence are recognized by the human immune system is still limited and scarcely available. Such knowledge would be of immediate relevance and great help for the design of new vaccines, facilitating the evaluation of potential immunogenic candidates, as well as monitoring the virus mutation events that would be transmitted through the human population.

Currently, there are at least 42 vaccine candidates around the world under development and evaluation at different stages against COVID-19 [17], also accordingly from what reported by WHO through its continuously undergoing landscapes documents concerning the COVID-19 candidate vaccines. These promising vaccine candidates deal with several vaccine technologies based on recombinant protein subunits [18], nucleic acids [19], non-replicating and replicating viral vectors [20, 21], protein constructs [22], virus-like particles [23], live-attenuated virus strains [24], inactivated virus [17], or human monoclonal antibodies (mAbs) [25]. Very recently, it has been shows that DTP vaccinations could protect against COVID-19 through potential cross-reactive immunity [26].

Today, challenges of continuing development of solutions for COVID-19 pandemic are a mandatory need. As never before, the application of modeling and simulation can actively design better vaccine prototypes, support decision making, decrease experimental costs and time, and eventually improve success rates of the trials. To this aim, in silico trials (ISTs) for design and testing medicines [27,28,29] can accelerate and speed-up the vaccine discovery pipeline, predicting any therapeutic failure and minimizing undesired effects.

Beyond traditional modeling techniques or applications, Agent-Based Models (ABMs) represent a paradigm that can cover the entire spectrum of the vaccine development process [30], especially for the quantification and prediction of the humoral and cellular response of a specific candidate vaccine as well as its efficacy [31].

The simulation platform we use from 15 years, named Universal Immune System Simulator (UISS), is based on agent-based methodology, which is able to brilliantly simulate each single entity of the immune system (and consequently its dynamics), along with the significant immune responses induced by a specific pathogen or stimulus. Recently, UISS provided different success stories in immunology field as it is most widely reported in the literature [32,33,34,35].

As the actual diagnostic strategies based on RT-qPCR [36] often fail in diagnose correctly COVID-19 patients (including also asymptomatic or false-negative ones) [37], in silico trials applied to the development of an effective vaccine are desirable.

We chose to analyze, within the wide landscape of potential candidate vaccines against SARS-CoV-2, a specific cross-neutralizing antibody that Wang et al. [38] suggest to be promising in targeting and binding a communal conserved epitope of SARS-CoV-2 and SARS-CoV on the spike receptor binding domain [39], through an independent mechanism of receptor binding inhibition.

As a case study, here we report a first application of UISS in silico platform to provide predictions of the efficacy of a potential therapy against COVID-19 based on a mAb strategy intervention like the one proposed by Wang et al.


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.

Fig. 1

SARS-CoV-2 disease model implemented in UISS. Main compartments (lung, and lymph-nodes) are delimited with dashed lines. Peripheral blood compartment is seen as connecting duct, not explicitly represented. The starting point is the SARS-CoV-2 droplets entrance in the upper respiratory tract (not shown). Then, all the main infection dynamics is described. The immune system cascade is shown as it was implemented, based on the latest research results published in specialized literature. For each entity, the localization (i.e., the biological compartment in which the entities are present) and the status (i.e., the differentiation states that an entity can own) are defined. The results of the immune system mounting process is the killing of the infected lung epithelial cells by the cytotoxic T lymphocyte and the local release of both chemokine factors and cytokines. At the humoral level, specific IgM (first) and IgG (after) directed against SARS-CoV-2 virus are released by plasma B cells. Regulatory system is also involved in the process. If the immune system machinery works correctly, regulatory arm shutdowns excessive cytokines storm, avoiding the severe prognosis of COVID-19. All entities are allowed to move with a uniform probability between neighboring lattices in the grid with an equal diffusion coefficient (Brownian motion). If a chemokines gradient is present, then to mimic short-range chemotaxis effects, higher probabilities of being chosen are given to sites containing chemokines

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).

Results and discussion

Tuning and validation of SARS-CoV-2 disease model

Scientific knowledge about SARS-CoV-2 is still not complete and research contributions appear every day. Apart from this, we used all the available literature data to compare the dynamics predicted by the UISS platform with all findings we were able to fetch. All the simulations we run represent the mean patient for the two different scenario we considered i.e., mild to moderate and the severe one. The mean patient was calculated from 100 different in silico simulations.

The first task we accomplished with success was the evaluation of the replication kinetic of SARS-CoV-2. To this end, we set a first use case simulation considering a digital patient in which a virus challenge dose of 0.1 multiplicity of infection (MOI) was administered at day 0. Simulation space was 5 cubic millimeters of lung tissue, 5 cubic millimeters of lymph tissue and 5 µl of peripheral blood. Figure 2 (right panel) shows that peak viral titers are reached by 48 h post-inoculation. We also plotted IL-6 dynamics: as reported in [9] the levels of IL-6 could be provide a prognosis on the severity of the infection.

Fig. 2

In silico SARS-CoV-2 viral dynamics and related CPE in a mild to moderate “mean in silico patient” scenario. In the left panel, one can observe the mild digital patient case in which a virus challenge dose of 0.1 multiplicity of infection (MOI) was administered at day 0 (green line). Peak viral titers are reached by 48 h post-inoculation. IL-6 dynamics and its related plasma levels (fg/μL) are also shown in the inner panel (purple line). In the right one, the dynamics of CPE on the lung infected cells is measured: they started at day 3.5 and peak around day 5. After 21 days, the simulated digital patient almost recovers from the infection. One can notice how UISS is capable to simulate, accordingly to the recent literature, the early viral clearance by day 10 post-onset in mild cases

We also measured cytopathic effects (CPE) on the lung compartment. CPE are defined as changes occurred in the infected cell that eventually lead its lysis or inability to reproduce. Figure 2 (left panel) highlights the dynamics of CPE: they started at day 3.5 and peak around day 5. After 21 days, the simulated digital patient almost recovers from the infections. These findings are in good agreement with actual literature [16, 55].

In a recent work, Liu et al. [56] reported that mild cases were found to have an early viral clearance, with 90% of these patients repeatedly testing negative by day 10 post-onset. At the same time, they found that all severe cases still tested positive at or beyond day 10 post-onset. Moreover, severe cases tended to have a higher viral load both at the beginning and later. In contrast, mild cases had early viral clearance 10 days post on-set. UISS was also able to reproduce this scenario. As one can see, Fig. 2 is in very good agreement for viral clearance.

We also were able to reproduce severe conditions acting on the immune system aging parameters obtaining results showed in Fig. 3. In this case, Fig. 3 (right panel) shows virus presence after day 10, until day 15, and its complete clearance about day 19. Moreover, CPE are much more severe and the recover from infection was clearly delayed (left panel). IL-6 dynamics shows a much more prominent peak of values. This is in very good agreement with latest literature data, as explained before.

Fig. 3

In silico SARS-CoV-2 viral dynamics and related CPE in a “mean in silico patient” severe scenario. In the left panel, one can observe the severe digital patient case in which a virus challenge dose of 0.1 multiplicity of infection (MOI) was administered at day 0 (green line). Peak viral titers are reached by 48 h post-inoculation. In addition, it is wort to note that virus persists after day 10, until day 15, and its complete clearance is around day 19. In the inner panel (purple line), IL-6 dynamics and its related plasma levels (fg/μL) are shown. IL-6 dynamics shows a much more prominent peak of values. This is in very good agreement with latest literature data, as explained within the manuscript. In the right panel, the dynamics of CPE on the lung infected cells is measured: in this case, CPE are much more severe and the recover from infection is clearly delayed. UISS is capable to simulate, accordingly to the recent literature, how the severe cases tend to have a higher viral load both at the beginning and later on

To validate the main immune system response of mild-to-moderate COVID-19 patients, we used the results available in [57]. In this work, the authors report the kinetics of immune responses in terms of activated CD4 + T cells, CD8 + T cells, IgM and IgG antibodies, detected in blood before symptomatic recovery. As one can see from Fig. 4, the kinetics of activated Th1 cells (panel A), activated CD8 T cells (panel B) and the IgM and IgG (panel C) predicted by the simulator are in good agreement with their findings.

Fig. 4

Cellular and humoral response mounted by the host immune system against SARS-CoV-2. Panel A shows the dynamics of CD4 + T cells, subtype 1 (Th1). Th1 are primed by dendritic cells that present the viral particles complexed with MHC-II of the host. Th1 cells help the activation of B cells, eventually favoring their iso-type switching to IgG producing plasma cell. B cells dynamics is depicted in panel B. Antigen activated B cells initially releases IgM; then, after interacting with Th1 and their released pro-inflammatory cytokines, they start to release specific IgG directed against SARS-CoV-2 virus

UISS IST to predict mAb efficacy against SARS-CoV-2

UISS is an immune system simulation platform that was designed to be applied to several and different scenarios, especially to carry on in silico trials to predict the efficacy of a specific prophylactic or therapeutic vaccine against a particular disease. In silico trials aim to strongly reduce the time to develop new effective therapeutics: this is particularly crucial in situation like the one we are facing with. As soon as a disease model incorporated into UISS is tuned and validated against available data, it can be used as an in silico lab to test new vaccines. In the previous section, we demonstrated that UISS-SARS-CoV-2 is able to reproduce and predict the main aspects of the viral infection.

As a working example, here we show how the platform can be immediately used to predict the efficacy of a human monoclonal antibody that neutralizes SARS-CoV-2 developed by Wang et al. [38]. In this work the authors suppose that the developed antibody (named 47D11) neutralizes SARS-CoV-2 through a yet unknown mechanism that is different from receptor binding interference. Hence, in implementing the mechanism of action of 47D11 into UISS computational framework we used the alternative mechanisms of coronavirus neutralization by receptor-binding domain (RBD) targeting antibodies that have been reported, including spike inactivation through antibody-induced destabilization of its prefusion structure, which Wang et al. indicated also applicable for 47D11. We then modeled the receptor interaction to trigger irreversible conformational changes in coronavirus spike proteins enabling membrane fusion, as described in [58].

The validation in silico trial consists in simulating the in vitro experiment conducted by Wang et al. where they showed that the monoclonal antibody was effective in contrasting SARS-CoV-2 to infect the target cells. To mimic the in vitro system that is an isolated system, we disabled both the lung and the lymph node compartments and we turned off all the immune system interactions in UISS. We then injected virus particles in the peripheral blood compartment of the simulator along with different concentrations of mAb. Only LEP-SARS-CoV-2 interaction has been allowed to happen. Figure 5 shows the obtained results.

Fig. 5

Antibody-mediated neutralization of SARS-CoV-2 infection on simulated lung epithelial cells. 10 ng/ml revealed the best concentration to obtain maximum clearance of the virus

As one can envisage from the figure, we reproduced with great accuracy the in vitro results reported by Wang et al. In particular, the computational framework was able to correctly predict the efficacy of 47D11 mAb simulating its mechanism of action that induces the entering of the virus inside LEP.

The most effective concentration is 10 μg/ml. This step makes the simulation platform ready for usage as a in silico trial to predict now the effects of a mAb therapy. For this purpose, we designed two kind of in silico experiments. The first one dealt with mAb vaccination strategy used to prevent the onset of infection. The second one involved mAb as interventional drug to treat already infected hosts.

mAb as a preventive vaccine

Figure 6 depicts the dynamics of LEP while we injected the 47D11 mAb at day 1 at a concentration of 10 μg/ml. Subsequently, we injected SARS-CoV-2 virus particles at day 7 (panel A), day 14 (panel B), month 1 (panel C), month 3 (panel D), month 6 (panel E) and after 1 year (panel F). Solid lines refer to the no-treated digital patient while dashed lines refer to mAb treated one. Blue lines depict actively infected LEP while red lines represent LEP to show the CPE.

Fig. 6

In silico trial of 47D11 to predict preventive efficacy. The overall prediction dynamics of LEP at the injection time of 47D11 mAb (day 1) at a concentration of 10 ng/ml after different SARS-CoV-2 challenge in time (af) is depicted. Specifically, one can observe the exposures at virus particles at day 7 (a), day 14 (b), month 1 (c), month 3 (d), month 6 (e) and after 1 year (f). Solid lines refer to the no-treated digital patient while dashed lines refer to mAb treated one. Blue lines depict actively infected LEP while red lines represent LEP to show the CPE. Panel A highlights that if a patient enters in contact with the virus just after 7 days post vaccination, he is fully protected from the infection (SARS-CoV-2 actively infected values of LEP are low). This is also true if a potential subject is being infected after 2 weeks (b), 1 month (c) or 3 months (d) after vaccination. Oppositely, in e (subject infected after 6 months) and f (subject infected after 1 year) mAb vaccination is practically ineffective in protecting the onset of the disease

Panel A of Fig. 6 clearly demonstrates that if a patient enters in contact with the virus just after 7 days post vaccination, he is fully protected from the infection: as one can see, SARS-CoV-2 actively infected values of LEP are low. This is also true if a potential subject is being infected after 2 weeks (panel B), 1 month (panel C) or 3 months (panel D) after vaccination. Things change for the other two cases i.e., panel E (subject infected after 6 months) or panel F (subject infected after 1 year). In this circumstance, the computational framework predicts that the mAb vaccination is practically ineffective in protecting the onset of the disease.

A second injection of the mAb is suggested around month 4 to extend the protection of the host for 1 year.

mAb as a therapeutic vaccine

The second experiment we designed is to use our in silico trial platform to predict the efficacy of a mAb-based vaccine in therapeutic settings.

We envisaged both mild-moderate case (the same digital patient type shown in Fig. 2) and severe case (the same digital patient shown in Fig. 3). For both cases, we administered the mAb vaccine 1 day and 2 days after the onset of infection.

Figure 7 shows different behavior for both mild-moderate case (panel A, mAb injected 1 day after the onset of infection, and panel B, mAb injected 2 days after the onset of infection) and severe case (panels C, mAb injected 1 day after the onset of infection, and panel D, mAb injected 2 days after the onset of infection). As Fig. 7 depicts, mAb vaccine is effective in preventing or strongly limiting the CPE. If the mAb vaccine is injected after 2 days, it is not able to protect the LPE of the host to be infected by the virus and consequently from COVID-19 pathology.

Fig. 7

In silico trial of 47D11 to predict therapeutic efficacy. Different behaviors for both mild-moderate case (a, b) and severe case (c, d) are shown. In a mAb is injected 1 day after the onset of infection. In b mAb is injected 2 days after the onset of infection. In c mAb is injected 1 day after the onset of infection. In d mAb is injected 2 days after the onset of infection. Blue lines depict actively infected LEP while red lines represent LEP to show the CPE. As one can notice, mAb vaccine is effective in preventing or strongly limiting the CPE only by 2 days of the infection


In this paper, we present an in silico platform that was demonstrated able to reproduce the main dynamics of SARS-CoV-2 virus and the elicited host immune response against it.

The disease model was implemented inside UISS computational framework, an in silico trial platform that has been applied to several biological scenario. UISS shows that the simulated SARS-CoV-2 dynamics is in very good agreement with the one described in the latest literature; also, the immune system response predicted by UISS against the virus mirrored the one observed in state of the art research data.

This validation step entitled UISS-SARS-CoV-2 to be used as a in silico lab to test the efficacy of potential vaccines for COVID-19, knowing a priori their mechanism of action. Hence, we set an in silico trial to test a recent vaccination strategy based on the employment of monoclonal antibodies directed against a specific target protein of the virus. The simulator is in good agreement in predicting the in vitro experiment outcome performed by the inventors of 47D11 mAb.

Finally, we designed two experimental settings to predict the efficacy of mAb vaccination when used in both preventive and therapeutic cases. We predicted that mAb is an effective therapy when used as a preventive vaccine (granting up 1-year protection when injected with two times schedule). Moreover, we envisaged that 47D11 mAb if effective only if administered in a very stringent time-frame if employed as a therapeutic strategy.

UISS-SARS-CoV-2 could potentially be used also to predict the effects of other kind of therapeutic strategies. However, it is designed to be focused on immune system machinery. This is not preventing to test antivirals compounds or other strategies that directly affect the virus dynamics or, in some way, stimulates the immune system to react in a more powerful way.

The web interface to access the simulator is freely available at Users can register or alternatively can access anonymously. UISS is implemented in ANSI C and it is supported on all platforms.

Availability and requirements

Project name UISS-COVID19.

Project home page

Operating system(s) Platform independent.

Programming language: C and Python

Other requirements: none

Any restrictions to use by non-academics: not applicable.

Availability of data and materials

The main computational framework is fully described in the paper.



Coronavirus disease 2019


Severe acute respiratory syndrome


World Health Organization

SARS-CoV-2, HCoV-19, 2019-nCoV:



Respiratory distress syndrome




Human monoclonal antibodies


In silico trials


Agent Based Models


Reverse transcription Polymerase Chain Reaction


Universal Immune System Simulator


In Silico TRIal for tuberculosis VAccine Development


Context of Use


Epithelial lung cells


Dendritic cells




Major histocompatibility class I


Major histocompatibility class II


Interferon type A


Interferon type B




Tumor necrosis factor alpha


Angiotensin-converting enzyme 2


Natural killer cells


Class M immunoglobulins


CD4 T cells


CD8 T cells


Interferon gamma


Class G immunoglobulins




Multiplicity of infection


Cytopathic effects


Receptor-binding domain


  1. 1.

    World Health Organization. Novel Coronavirus (2019-nCoV) situation reports. 2020.

  2. 2.

    Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG, et al. A new coronavirus associated with human respiratory disease in China. Nature. 2020;579:265–9.

    CAS  Article  Google Scholar 

  3. 3.

    Raoult D, Zumla A, Locatelli F, Ippolito G, Kroemer G. Coronavirus infections: epidemiological, clinical and immunological features and hypotheses. Cell Stress. 2020;4:66–75.

    CAS  Article  Google Scholar 

  4. 4.

    Peiris JSM, Guan Y, Yuen KY. Severe acute respiratory syndrome. Nat Med. 2004;10:S88–97.

    CAS  Article  Google Scholar 

  5. 5.

    Wu J, Xu F, Zhou W, Feikin DR, Lin C-Y, He X, et al. Risk factors for SARS among persons without known contact with SARS patients, Beijing, China. Emerg Infect Dis. 2004;10:210–6.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Xu J, Zhao S, Teng T, Abdalla AE, Zhu W, Xie L, et al. Systematic comparison of two animal-to-human transmitted human coronaviruses: SARS-CoV-2 and SARS-CoV. Viruses. 2020;12:244.

    Article  Google Scholar 

  7. 7.

    Lai CC, Shih TP, Ko WC, Tang HJ, Hsueh PR. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges. Int J Antimicrob Agents. 2020;55:105924.

    CAS  Article  Google Scholar 

  8. 8.

    Chen X, Zhao B, Qu Y, Chen Y, Xiong J, Feng Y, et al. Detectable serum SARS-CoV-2 viral load (RNAaemia) is closely associated with drastically elevated interleukin 6 (IL-6) level in critically ill COVID-19 patients. medRxiv. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Pedersen SF, Ho Y-C. SARS-CoV-2: a storm is raging. J Clin Investig. 2020;130:2202–5.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Pennisi M, Lanza G, Falzone L, Fisicaro F, Ferri R, Bella R. SARS-CoV-2 and the nervous system: from clinical features to molecular mechanisms. Int J Mol Sci. 2020;21:5475.

    CAS  Article  Google Scholar 

  11. 11.

    Aryal MR, Gosain R, Donato A, Pathak R, Bhatt VR, Katel A, et al. Venous thromboembolism in COVID-19: towards an ideal approach to thromboprophylaxis, screening, and treatment. Curr Cardiol Rep. 2020;22:52.

    Article  Google Scholar 

  12. 12.

    Vivarelli S, Falzone L, Grillo CM, Scandurra G, Torino F, Libra M. Cancer management during covid-19 pandemic: Is immune checkpoint inhibitors-based immunotherapy harmful or beneficial? Cancers. 2020;12:2237.

    Article  Google Scholar 

  13. 13.

    Capobianchi MR, Rueca M, Messina F, Giombini E, Carletti F, Colavita F, et al. Molecular characterization of SARS-CoV-2 from the first case of COVID-19 in Italy. Clin Microbiol Infect. 2020;26:954–6.

    CAS  Article  Google Scholar 

  14. 14.

    Smith EC, Blanc H, Vignuzzi M, Denison MR. Coronaviruses lacking exoribonuclease activity are susceptible to lethal mutagenesis: evidence for proofreading and potential therapeutics. PLoS Pathog. 2013;9:e1003565.

    CAS  Article  Google Scholar 

  15. 15.

    Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell. 2020;181:271.e8-280.e8.

    Article  Google Scholar 

  16. 16.

    Zou L, Ruan F, Huang M, Liang L, Huang H, Hong Z, et al. SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N Engl J Med. 2020;382:1177–9.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Amanat Fatima KF. SARS-CoV-2 vaccines: status report. Cell. 2020;52:1–7.

    Google Scholar 

  18. 18.

    Kim E, Erdos G, Huang S, Kenniston TW, Balmert SC, Carey CD, et al. Microneedle array delivered recombinant coronavirus vaccines: immunogenicity and rapid translational development. EBioMedicine. 2020;55:102743.

    Article  Google Scholar 

  19. 19.

    Pang J, Wang MX, Ang IYH, Tan SHX, Lewis RF, Chen JI-P, et al. Potential rapid diagnostics, vaccine and therapeutics for 2019 novel coronavirus (2019-nCoV): a systematic review. J Clin Med. 2020;9:623.

    CAS  Article  Google Scholar 

  20. 20.

    Tu Y-F, Chien C-S, Yarmishyn AA, Lin Y-Y, Luo Y-H, Lin Y-T, et al. A review of SARS-CoV-2 and the ongoing clinical trials. Int J Mol Sci. 2020;21:2657.

    CAS  Article  Google Scholar 

  21. 21.

    Rabaan AA, Al-ahmed SH, Sah R, Tiwari R, Iqbal M, Patel SK, et al. SARS-CoV-2/COVID-19 and advances in developing potential therapeutics and vaccines to counter this emerging pandemic virus—a review. Ann Clin Microbiol Antimicrob. 2020;4:1–46.

    Google Scholar 

  22. 22.

    Zakhartchouk AN, Viswanathan S, Moshynskyy I, Petric M, Babiuk LA. Optimization of a DNA vaccine against SARS. DNA Cell Biol. 2007;26:721–6.

    CAS  Article  Google Scholar 

  23. 23.

    Zhang J, Zeng H, Gu J, Li H, Zheng L, Zou Q. Progress and prospects on vaccine development against SARS-CoV-2. Vaccines. 2020;8:153.

    CAS  Article  Google Scholar 

  24. 24.

    Shang W, Yang Y, Rao Y, Rao X. The outbreak of SARS-CoV-2 pneumonia calls for viral vaccines. NPJ Vaccines. 2020;5:18.

    CAS  Article  Google Scholar 

  25. 25.

    Chen X, Li R, Pan Z, Qian C, Yang Y, You R, et al. Human monoclonal antibodies block the binding of SARS-CoV-2 spike protein to angiotensin converting enzyme 2 receptor; 2020.

  26. 26.

    Reche PA. Cross-reactive immunity from combination DTP vaccines could protect against COVID-19. OSF Prepr. 2020.

    Article  Google Scholar 

  27. 27.

    Pappalardo F, Russo G, Tshinanu FM, Viceconti M. In silico clinical trials: concepts and early adoptions. Brief Bioinform. 2019;20:1699–708.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Viceconti M, Pappalardo F, Rodriguez B, Horner M, Bischoff J, Musuamba Tshinanu F. In silico trials: Verification, validation and uncertainty quantification of predictive models used in the regulatory evaluation of biomedical products. Methods. 2020.

    Article  PubMed  Google Scholar 

  29. 29.

    Carlier A, Vasilevich A, Marechal M, de Boer J, Geris L. In silico clinical trials for pediatric orphan diseases. Sci Rep. 2018;8:2465.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Cosgrove J, Butler J, Alden K, Read M, Kumar V, Cucurull-Sanchez L, et al. Agent-based modeling in systems pharmacology. CPT Pharmacometr Syst Pharmacol. 2015;4:615–29.

    CAS  Article  Google Scholar 

  31. 31.

    Bonabeau E. Agent-based modeling: methods and techniques for simulating human systems. Proc Natl Acad Sci USA. 2002;99:7280–7.

    CAS  Article  Google Scholar 

  32. 32.

    Pappalardo F, Fichera E, Paparone N, Lombardo A, Pennisi M, Russo G, et al. A computational model to predict the immune system activation by citrus-derived vaccine adjuvants. Bioinformatics. 2016;32:2672–80.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Pennisi M, Russo G, Ravalli S, Pappalardo F. Combining agent based-models and virtual screening techniques to predict the best citrus-derived vaccine adjuvants against human papilloma virus. BMC Bioinform. 2017;18:544.

    CAS  Article  Google Scholar 

  34. 34.

    Pennisi M, Russo G, Sgroi G, Bonaccorso A, Parasiliti Palumbo GA, Fichera E, et al. Predicting the artificial immunity induced by RUTI® vaccine against tuberculosis using universal immune system simulator (UISS). BMC Bioinform. 2019;20:504.

    CAS  Article  Google Scholar 

  35. 35.

    Pappalardo F, Russo G, Pennisi M, Parasiliti Palumbo GA, Sgroi G, Motta S, et al. The potential of computational modeling to predict disease course and treatment response in patients with relapsing multiple sclerosis. Cells. 2020;9:586.

    CAS  Article  PubMed Central  Google Scholar 

  36. 36.

    Tahamtan A, Ardebili A. Real-time RT-PCR in COVID-19 detection: issues affecting the results. Expert Rev Mol Diagn. 2020;20:453–4.

    CAS  Article  Google Scholar 

  37. 37.

    Falzone L, Musso N, Gattuso G, Bongiorno D, Palermo CI, Scalia G, et al. Sensitivity assessment of droplet digital PCR for SARS-CoV-2 detection. Int J Mol Med. 2020;46:957–64.

    CAS  Article  Google Scholar 

  38. 38.

    Wang C, Li W, Drabek D, Okba NMA, van Haperen R, Osterhaus ADME, et al. A human monoclonal antibody blocking SARS-CoV-2 infection. Nat Commun. 2020;11:2251.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Lan J, Ge J, Yu J, Shan S, Zhou H, Fan S, et al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature. 2020;581:215–20.

    CAS  Article  Google Scholar 

  40. 40.

    Gullo F, van der Garde M, Russo G, Pennisi M, Motta S, Pappalardo F, et al. Computational modeling of the expansion of human cord blood CD133+ hematopoietic stem/progenitor cells with different cytokine combinations. Bioinformatics. 2015;31:2514–22.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Ragusa MA, Russo G. ODEs approaches in modeling fibrosis. Phys Life Rev. 2016;17:112–3.

    Article  PubMed  Google Scholar 

  42. 42.

    Komarova NL, Wodarz D. ODE models for oncolytic virus dynamics. J Theor Biol. 2010;263:530–43.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Russo G, Pennisi M, Boscarino R, Pappalardo F. Continuous petri nets and microRNA analysis in melanoma. IEEE/ACM Trans Comput Biol Bioinform. 2017;15:1492–9.

    Article  PubMed  Google Scholar 

  44. 44.

    Liu F, Heiner M, Gilbert D. Fuzzy petri nets for modelling of uncertain biological systems. Brief Bioinform. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Bianca C, Brézin L. Modeling the antigen recognition by B-cell and T-cell receptors through thermostatted kinetic theory methods. Int J Biomath. 2017;10:1750072.

    Article  Google Scholar 

  46. 46.

    Bianca C, Riposo J. Mimic therapeutic actions against keloid by thermostatted kinetic theory methods. Eur Phys J Plus. 2015;130:159.

    Article  Google Scholar 

  47. 47.

    Palladini A, Nicoletti G, Pappalardo F, Murgo A, Grosso V, Stivani V, et al. In silico modeling and in vivo efficacy of cancer-preventive vaccinations. Cancer Res. 2010;70:7755–63.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Pennisi M, Pappalardo F, Palladini A, Nicoletti G, Nanni P, Lollini P-L, et al. Modeling the competition between lung metastases and the immune system using agents. BMC Bioinform. 2010;11(Suppl 7):S13.

    Article  Google Scholar 

  49. 49.

    Pappalardo F, Forero IM, Pennisi M, Palazon A, Melero I, Motta S. SimB16: modeling induced immune system response against B16-melanoma. PLoS ONE. 2011;6:e26523.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Pappalardo F, Musumeci S, Motta S. Modeling immune system control of atherogenesis. Bioinformatics. 2008;24:1715–21.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Patente TA, Pinho MP, Oliveira AA, Evangelista GCM, Bergami-Santos PC, Barbuto JAM. Human dendritic cells: their heterogeneity and clinical application potential in cancer immunotherapy. Front Immunol. 2019;9:3176.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Chen Y, Liu Q, Guo D. Emerging coronaviruses: genome structure, replication, and pathogenesis. J Med Virol. 2020;92:418–23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Matheoud D, Perie L, Hoeffel G, Vimeux L, Parent I, Maranon C, et al. Cross-presentation by dendritic cells from live cells induces protective immune responses in vivo. Blood. 2010;115:4412–20.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395:497–506.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Harcourt J, Tamin A, Lu X, Kamili S, Sakthivel SK, Murray J, et al. Severe acute respiratory syndrome coronavirus 2 from patient with 2019 novel coronavirus disease, United States. Emerg Infect Dis. 2020;26:1266–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Liu Y, Yan L-M, Wan L, Xiang T-X, Le A, Liu J-M, et al. Viral dynamics in mild and severe cases of COVID-19. Lancet Infect Dis. 2020;2019:2019–20.

    Article  Google Scholar 

  57. 57.

    Thevarajan I, Nguyen THO, Koutsakos M, Druce J, Caly L, van de Sandt CE, et al. Breadth of concomitant immune responses prior to patient recovery: a case report of non-severe COVID-19. Nat Med. 2020;26:453–5.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Walls AC, Xiong X, Park Y-J, Tortorici MA, Snijder J, Quispe J, et al. Unexpected receptor functional mimicry elucidates activation of coronavirus fusion. Cell. 2019;176(1026–1039):e15.

    CAS  Article  Google Scholar 

Download references


Not applicable.

About this supplement

This article has been published as part of BMC Bioinformatics Volume 21 Supplement 17 2020: Selected papers from the 3rd International Workshop on Computational Methods for the Immune System Function (CMISF 2019). The full contents of the supplement are available at


Publication costs are funded by the University of Catania internal grant.

Author information




GR: developed the conceptual model, checked the biological adherence and meaning, analyzed data, wrote the manuscript. MP: performed numerical simulations, analyzed data and wrote the manuscript. EF: checked the biological adherence and meaning, analyzed data, wrote the manuscript. SM: gave useful insights about numerical simulations, analyzed data, wrote the manuscript. GRAC: checked the biological adherence and meaning, analyzed data, wrote the manuscript. MV: gave useful insights about in silico clinical trial development, curated the validation of the simulator. FP: supervised, conceived and developed the model, gave computational immunological knowledge, drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Giuseppina Raciti or Francesco Pappalardo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Russo, G., Pennisi, M., Fichera, E. et al. In silico trial to test COVID-19 candidate vaccines: a case study with UISS platform. BMC Bioinformatics 21, 527 (2020).

Download citation


  • Agent-based model
  • Human monoclonal antibodies
  • In silico trials
  • SARS-CoV-2
  • Vaccines