Search

Abstract When it was first introduced in 2000, reverse vaccinology was defined as an in silico approach that begins with the pathogen's genomic sequence. It concludes with a list of potential proteins with a possible, but not necessarily, list of peptide candidates that need to be experimentally confirmed for vaccine production. During the subsequent years, reverse vaccinology has dramatically changed: now it consists of a large number of bioinformatics tools and processes, namely subtractive proteomics, computational vaccinology, immunoinformatics, and in silico related procedures. However, the state of the art of reverse vaccinology still misses the ability to predict the efficacy of the proposed vaccine formulation. Here, we describe how to fill the gap by introducing an advanced immune system simulator that tests the efficacy of a vaccine formulation against the disease for which it has been designed. As a working example, we entirely apply this advanced reverse vaccinology approach to design and predict the efficacy of a potential vaccine formulation against influenza H5N1. Climate change and melting glaciers are critical due to reactivating frozen viruses and emerging new pandemics. H5N1 is one of the potential strains present in icy lakes that can raise a pandemic. Investigating structural antigen protein is the most profitable therapeutic pipeline to generate an effective vaccine against H5N1. In particular, we designed a multi-epitope vaccine based on predicted epitopes of hemagglutinin and neuraminidase proteins that potentially trigger B-cells, CD4, and CD8 T-cell immune responses. Antigenicity and toxicity of all predicted CTL, Helper T-lymphocytes, and B-cells epitopes were evaluated, and both antigenic and non-allergenic epitopes were selected. From the perspective of advanced reverse vaccinology, the Universal Immune System Simulator, an in silico trial computational framework, was applied to estimate vaccine efficacy using a cohort of 100 digital patients.

Multiple sclerosis is an autoimmune inflammatory disease that affects the central nervous system through chronic demyelination and loss of oligodendrocytes. Since the relapsing-remitting form is the most prevalent, relapse-reducing therapies are a primary choice for specialists. Universal Immune System Simulator is an agent-based model that simulates the human immune system dynamics under physiological conditions and during several diseases, including multiple sclerosis. In this work, we extended the UISS-MS disease layer by adding two new treatments, i.e., cladribine and ocrelizumab, to show that UISS-MS can be potentially used to predict the effects of any existing or newly designed treatment against multiple sclerosis. To retrospectively validate UISS-MS with ocrelizumab and cladribine, we extracted the clinical and MRI data from patients included in two clinical trials, thus creating specific cohorts of digital patients for predicting and validating the effects of the considered drugs. The obtained results mirror those of the clinical trials, demonstrating that UISS-MS can correctly simulate the mechanisms of action and outcomes of the treatments. The successful retrospective validation concurred to confirm that UISS-MS can be considered a digital twin solution to be used as a support system to inform clinical decisions and predict disease course and therapeutic response at a single patient level.

COVID-19 was declared a pandemic in March 2020, and since then, it has not stopped spreading like wildfire in almost every corner of the world, despite the many efforts made to stem its spread. SARS-CoV-2 has one of the biggest genomes among RNA viruses and presents unique characteristics that differentiate it from other coronaviruses, making it even more challenging to find a cure or vaccine that is efficient enough. This work aims, using RNA sequencing (RNA-Seq) data, to evaluate whether the expression of specific human genes in the host can vary in different grades of disease severity and to determine the molecular origins of the differences in response to SARS-CoV-2 infection in different patients. In addition to quantifying gene expression, data coming from RNA-Seq allow for the discovery of new transcripts, the identification of alternative splicing events, the detection of allele-specific expression, and the detection of post-transcriptional alterations. For this reason, we performed differential expression analysis on different expression profiles of COVID-19 patients, using RNA-Seq data coming from NCBI public repository, and we obtained the lists of all differentially expressed genes (DEGs) emerging from 7 experimental conditions. We performed a Gene Set Enrichment Analysis (GSEA) on these genes to find possible correlations between DEGs and known disease phenotypes. We mainly focused on DEGs coming out from the analysis of the contrasts involving severe conditions to infer any possible relation between a worsening of the clinical picture and an over-representation of specific genes. Based on the obtained results, this study indicates a small group of genes that result up-regulated in the severe form of the disease. EXOSC5, MESD, REXO2, and TRMT2A genes are not differentially expressed or not present in the other conditions, being for that reason, good biomarkers candidates for the severe form of COVID-19 disease. The use of specific over-expressed genes, whether up-regulated or down-regulated, which have an individual role in each different condition of COVID-19 as a biomarker, can assist in early diagnosis.

Abstract Background Influenza A virus is one of the leading causes of annual mortality. The emerging of novel escape variants of the influenza A virus is still a considerable challenge in the annual process of vaccine production. The evolution of vaccines ranks among the most critical successes in medicine and has eradicated numerous infectious diseases. Recently, multi-epitope vaccines, which are based on the selection of epitopes, have been increasingly investigated. Results This study utilized an immunoinformatic approach to design a recombinant multi-epitope vaccine based on a highly conserved epitope of hemagglutinin, neuraminidase, and membrane matrix proteins with fewer changes or mutate over time. The potential B cells, cytotoxic T lymphocytes (CTL), and CD4 T cell epitopes were identified. The recombinant multi-epitope vaccine was designed using specific linkers and a proper adjuvant. Moreover, some bioinformatics online servers and datasets were used to evaluate the immunogenicity and chemical properties of selected epitopes. In addition, Universal Immune System Simulator (UISS) in silico trial computational framework was run after influenza exposure and recombinant multi-epitope vaccine administration, showing a good immune response in terms of immunoglobulins of class G (IgG), T Helper 1 cells (TH1), epithelial cells (EP) and interferon gamma (IFN-g) levels. Furthermore, after a reverse translation (i.e., convertion of amino acid sequence to nucleotide one) and codon optimization phase, the optimized sequence was placed between the two EcoRV/MscI restriction sites in the PET32a+ vector. Conclusions The proposed “Recombinant multi-epitope vaccine” was predicted with unique and acceptable immunological properties. This recombinant multi-epitope vaccine can be successfully expressed in the prokaryotic system and accepted for immunogenicity studies against the influenza virus at the in silico level. The multi-epitope vaccine was then tested with the Universal Immune System Simulator (UISS) in silico trial platform. It revealed slight immune protection against the influenza virus, shedding the light that a multistep bioinformatics approach including molecular and cellular level is mandatory to avoid inappropriate vaccine efficacy predictions.

Background and Aims: The rate of variation in various genes of a bacterial species is different during evolution. Therefore, in systematic bacterial studies many researchers compare the phylogenetic tree of a particular gene to the standard tree of an rRNA gene. Regarding the importance of 16SrRNA gene and the evolutional process of RecA protein family, we investigated the changes in the selected phylum of bacterial RecA family in comparing with the 16SrRNA gene. Materials and Methods: For this purpose, sequences of the RecA protein family were extracted from Uniprot database and categorized by using CD-hit algorithm. One species was selected from each category. Then we found 16SrRNA complete sequences for same species. After determining the index, based on the Average Alignment Score (AAS), the 16s-taxonomic tree was obtained. Furthermore, Similar calculations were considered for corresponding RecA proteins phylums. Results: By comparing amount of AAS in 16srRNA phylums and RecA phylums, we observed that the Actinobacteria phylum is the closest to the header phylum in the 16s-taxonomic tree, but the same phylum in the RecA is the most distant to the header phylum; and the position of the cyanobacteria phylum remains the same in both trees, which indicates the least amount of changes in the genus and species of this phylum. Conclusion: The 16s-taxonomic tree which is presented in this study for the first time is different from the available bioinformatics algorithms for tree drawing. Finding the species with the highest and lowest rates of changes, can be a type of prediction method for indicating the reasons why bacteria become resistant to drugs over a long period of time.