- Open Access
In silico trial to test COVID-19 candidate vaccines: a case study with UISS platform
BMC Bioinformatics volume 21, Article number: 527 (2020)
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 .
2019-nCoV (also referred to as SARS-CoV-2 or HCoV-19) , is the seventh coronavirus known to infect humans along with SARS-CoV, MERS-CoV, HKU1, NL63, OC43 and 229E . While these last four coronaviruses are associated with mild symptoms, SARS-CoV, MERS-CoV and SARS-CoV-2 can cause severe acute respiratory syndrome , especially in elderlies, of which men, and those individuals with comorbidities and immunocompromised conditions ). Although it is similar to SARS-CoV, SARS-CoV-2 has an improved ability for pathogenicity . 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 . 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) . In this scenario, it seems that IL-6 owns an important driving role on the cytokine storm, leading to lung damage and reduced survival .
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 , the presence of minority variants generated during SARS-CoV-2 replication , the involved cellular factors that favors SARS-CoV-2 cell entry , 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 .
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 , 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 , nucleic acids , non-replicating and replicating viral vectors [20, 21], protein constructs , virus-like particles , live-attenuated virus strains , inactivated virus , or human monoclonal antibodies (mAbs) . Very recently, it has been shows that DTP vaccinations could protect against COVID-19 through potential cross-reactive immunity .
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 , especially for the quantification and prediction of the humoral and cellular response of a specific candidate vaccine as well as its efficacy .
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  often fail in diagnose correctly COVID-19 patients (including also asymptomatic or false-negative ones) , 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.  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 , 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 . 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  and derived lung metastases , melanoma , atherosclerosis , multiple sclerosis  and influenza .
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 .
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) . The virus is eventually captured by dendritic cells (DC) and macrophages (M).
DC are the main antigen processing cells of the immune system  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 .
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 ) 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 . 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  the levels of IL-6 could be provide a prognosis on the severity of the infection.
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.  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.
To validate the main immune system response of mild-to-moderate COVID-19 patients, we used the results available in . 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.
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. . 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 .
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.
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.
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.
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 https://combine.dmi.unict.it/UISS-COVID19. 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 https://combine.dmi.unict.it/UISS-COVID19/.
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
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
Class G immunoglobulins
Multiplicity of infection
World Health Organization. Novel Coronavirus (2019-nCoV) situation reports. 2020.
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.
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.
Peiris JSM, Guan Y, Yuen KY. Severe acute respiratory syndrome. Nat Med. 2004;10:S88–97.
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. https://doi.org/10.3201/eid1002.030730.
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.
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.
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. https://doi.org/10.1093/cid/ciaa449.
Pedersen SF, Ho Y-C. SARS-CoV-2: a storm is raging. J Clin Investig. 2020;130:2202–5. https://doi.org/10.1172/JCI137647.
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.
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.
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.
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.
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.
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.
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. https://doi.org/10.1056/NEJMc2001737.
Amanat Fatima KF. SARS-CoV-2 vaccines: status report. Cell. 2020;52:1–7.
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.
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.
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.
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.
Zakhartchouk AN, Viswanathan S, Moshynskyy I, Petric M, Babiuk LA. Optimization of a DNA vaccine against SARS. DNA Cell Biol. 2007;26:721–6.
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.
Shang W, Yang Y, Rao Y, Rao X. The outbreak of SARS-CoV-2 pneumonia calls for viral vaccines. NPJ Vaccines. 2020;5:18.
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.
Reche PA. Cross-reactive immunity from combination DTP vaccines could protect against COVID-19. OSF Prepr. 2020. https://doi.org/10.31219/osf.io/sbgy3.
Pappalardo F, Russo G, Tshinanu FM, Viceconti M. In silico clinical trials: concepts and early adoptions. Brief Bioinform. 2019;20:1699–708. https://doi.org/10.1093/bib/bby043.
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. https://doi.org/10.1016/j.ymeth.2020.01.011.
Carlier A, Vasilevich A, Marechal M, de Boer J, Geris L. In silico clinical trials for pediatric orphan diseases. Sci Rep. 2018;8:2465. https://doi.org/10.1038/s41598-018-20737-y.
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.
Bonabeau E. Agent-based modeling: methods and techniques for simulating human systems. Proc Natl Acad Sci USA. 2002;99:7280–7.
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. https://doi.org/10.1093/bioinformatics/btw293.
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. https://doi.org/10.1186/s12859-017-1961-9.
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. https://doi.org/10.1186/s12859-019-3045-5.
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. https://doi.org/10.3390/cells9030586.
Tahamtan A, Ardebili A. Real-time RT-PCR in COVID-19 detection: issues affecting the results. Expert Rev Mol Diagn. 2020;20:453–4.
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.
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. https://doi.org/10.1038/s41467-020-16256-y.
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.
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. https://doi.org/10.1093/bioinformatics/btv172.
Ragusa MA, Russo G. ODEs approaches in modeling fibrosis. Phys Life Rev. 2016;17:112–3. https://doi.org/10.1016/j.plrev.2016.05.012.
Komarova NL, Wodarz D. ODE models for oncolytic virus dynamics. J Theor Biol. 2010;263:530–43. https://doi.org/10.1016/j.jtbi.2010.01.009.
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. https://doi.org/10.1109/TCBB.2017.2733529.
Liu F, Heiner M, Gilbert D. Fuzzy petri nets for modelling of uncertain biological systems. Brief Bioinform. 2018. https://doi.org/10.1093/bib/bby118.
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.
Bianca C, Riposo J. Mimic therapeutic actions against keloid by thermostatted kinetic theory methods. Eur Phys J Plus. 2015;130:159.
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. https://doi.org/10.1158/0008-5472.CAN-10-0701.
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. https://doi.org/10.1186/1471-2105-11-S7-S13.
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. https://doi.org/10.1371/journal.pone.0026523.
Pappalardo F, Musumeci S, Motta S. Modeling immune system control of atherogenesis. Bioinformatics. 2008;24:1715–21. https://doi.org/10.1093/bioinformatics/btn306.
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. https://doi.org/10.3389/fimmu.2018.03176.
Chen Y, Liu Q, Guo D. Emerging coronaviruses: genome structure, replication, and pathogenesis. J Med Virol. 2020;92:418–23. https://doi.org/10.1002/jmv.25681.
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. https://doi.org/10.1182/blood-2009-11-255935.
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. https://doi.org/10.1016/S0140-6736(20)30183-5.
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. https://doi.org/10.3201/eid2606.200516.
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. https://doi.org/10.1016/S1473-3099(20)30232-2.
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. https://doi.org/10.1038/s41591-020-0819-2.
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. https://doi.org/10.1016/j.cell.2018.12.028.
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 https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-21-supplement-17.
Publication costs are funded by the University of Catania internal grant.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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). https://doi.org/10.1186/s12859-020-03872-0
- Agent-based model
- Human monoclonal antibodies
- In silico trials