Search

Vaccination is one of the most significant inventions in medicine. Reverse vaccinology (RV) is a state-of-the-art technique to predict vaccine candidates from pathogen's genome(s). To promote vaccine development, we updated Vaxign2, the first web-based vaccine design program using reverse vaccinology with machine learning. Vaxign2 is a comprehensive web server for rational vaccine design, consisting of predictive and computational workflow components. The predictive part includes the original Vaxign filtering-based method and a new machine learning-based method, Vaxign-ML. The benchmarking results using a validation dataset showed that Vaxign-ML had superior prediction performance compared to other RV tools. Besides the prediction component, Vaxign2 implemented various post-prediction analyses to significantly enhance users’ capability to refine the prediction results based on different vaccine design rationales and considerably reduce user time to analyze the Vaxign/Vaxign-ML prediction results. Users provide proteome sequences as input data, select candidates based on Vaxign outputs and Vaxign-ML scores, and perform post-prediction analysis. Vaxign2 also includes precomputed results from approximately 1 million proteins in 398 proteomes of 36 pathogens. As a demonstration, Vaxign2 was used to effectively analyse SARS-CoV-2, the coronavirus causing COVID-19. The comprehensive framework of Vaxign2 can support better and more rational vaccine design. Vaxign2 is publicly accessible at http://www.violinet.org/vaxign2., Graphical Abstract Graphical AbstractVaxign2 is the first comprehensive web server application for vaccine design using reverse vaccinology and machine learning, showing superior prediction performance compared to other similar tools. Vaxign2 consists of (i) protective antigen prediction from a user-provided protein(s) or proteome of the target pathogen and (ii) post-prediction analyses of protein functions and immunogenicity.

Graphical abstract, Highlights • Evolutionary design of new SARS-CoV-2 spike proteins by perturbing the core protein sequence without changing the surface conformation. • Computationally designed SARS-CoV-2 spike proteins contain additional MHC-II T cell promiscuous epitopes with enhanced T cell immunity. • Structurally designed SARS-CoV-2 spike proteins preserve function and B cell immunity. • The newly designed SARS-CoV-2 spike proteins can be explored as new COVID-19 vaccine candidates., The development of effective and safe vaccines is the ultimate way to efficiently stop the ongoing COVID-19 pandemic, which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Built on the fact that SARS-CoV-2 utilizes the association of its Spike (S) protein with the human Angiotensin-converting enzyme 2 (ACE2) receptor to invade host cells, we computationally redesigned the S protein sequence to improve its immunogenicity and antigenicity. Toward this purpose, we extended an evolutionary protein design algorithm, EvoDesign, to create thousands of stable S protein variants that perturb the core protein sequence but keep the surface conformation and B cell epitopes. The T cell epitope content and similarity scores of the perturbed sequences were calculated and evaluated. Out of 22,914 designs with favorable stability energy, 301 candidates contained at least two pre-existing immunity-related epitopes and had promising immunogenic potential. The benchmark tests showed that, although the epitope restraints were not included in the scoring function of EvoDesign, the top S protein design successfully recovered 31 out of the 32 major histocompatibility complex (MHC) -II T cell promiscuous epitopes in the native S protein, where two epitopes were present in all seven human coronaviruses. Moreover, the newly designed S protein introduced nine new MHC-II T cell promiscuous epitopes that do not exist in the wildtype SARS-CoV-2. These results demonstrated a new and effective avenue to enhance a target protein’s immunogenicity using rational protein design, which could be applied for new vaccine design against COVID-19 and other human viruses.

Motivation Reverse vaccinology (RV) is a milestone in rational vaccine design, and machine learning (ML) has been applied to enhance the accuracy of RV prediction. However, ML-based RV still faces challenges in prediction accuracy and program accessibility. Results This study presents Vaxign-ML, a supervised ML classification to predict bacterial protective antigens (BPAgs). To identify the best ML method with optimized conditions, five ML methods were tested with biological and physiochemical features extracted from well-defined training data. Nested 5-fold cross-validation and leave-one-pathogen-out validation were used to ensure unbiased performance assessment and the capability to predict vaccine candidates against a new emerging pathogen. The best performing model (eXtreme Gradient Boosting) was compared to three publicly available programs (Vaxign, VaxiJen, and Antigenic), one SVM-based method, and one epitope-based method using a high-quality benchmark dataset. Vaxign-ML showed superior performance in predicting BPAgs. Vaxign-ML is hosted in a publicly accessible web server and a standalone version is also available. Availability and implementation Vaxign-ML website at http://www.violinet.org/vaxign/vaxign-ml, Docker standalone Vaxign-ML available at https://hub.docker.com/r/e4ong1031/vaxign-ml and source code is available at https://github.com/VIOLINet/Vaxign-ML-docker. Supplementary information Supplementary data are available at Bioinformatics online.

Rational vaccine design, especially vaccine antigen identification and optimization, is critical to successful and efficient vaccine development against various infectious diseases including coronavirus disease 2019 (COVID-19). In general, computational vaccine design includes three major stages: (i) identification and annotation of experimentally verified gold standard protective antigens through literature mining, (ii) rational vaccine design using reverse vaccinology (RV) and structural vaccinology (SV) and (iii) post-licensure vaccine success and adverse event surveillance and its usage for vaccine design. Protegen is a database of experimentally verified protective antigens, which can be used as gold standard data for rational vaccine design. RV predicts protective antigen targets primarily from genome sequence analysis. SV refines antigens through structural engineering. Recently, RV and SV approaches, with the support of various machine learning methods, have been applied to COVID-19 vaccine design. The analysis of post-licensure vaccine adverse event report data also provides valuable results in terms of vaccine safety and how vaccines should be used or paused. Ontology standardizes and incorporates heterogeneous data and knowledge in a human- and computer-interpretable manner, further supporting machine learning and vaccine design. Future directions on rational vaccine design are discussed.

Tuberculosis (TB) is the leading cause of death of any single infectious agent, having led to 1.4 million deaths in 2019 alone. Moreover, an estimated one-quarter of the global population is latently infected with Mycobacterium tuberculosis (MTB), presenting a huge pool of potential future disease. Nonetheless, the only currently licensed TB vaccine fails to prevent the activation of latent TB infections (LTBI). These facts together illustrate the desperate need for a more effective TB vaccine strategy that can prevent both primary infection and the activation of LTBI. In this study, we employed a machine learning-based reverse vaccinology approach to predict the likelihood that each protein within the proteome of MTB laboratory reference strain H37Rv would be a protective antigen (PAg). The proteins predicted most likely to be a PAg were assessed for their belonging to a protein family of previously established PAgs, the relevance of their biological processes to MTB virulence and latency, and finally the immunogenic potential that they may provide in terms of the number of promiscuous epitopes within each. This study led to the identification of 16 proteins with the greatest vaccine potential for further in vitro and in vivo studies. It also demonstrates the value of computational methods in vaccine development.

Background Drug adverse events (AEs), or called adverse drug events (ADEs), are ranked one of the leading causes of mortality. The Ontology of Adverse Events (OAE) has been widely used for adverse event AE representation, standardization, and analysis. OAE-based ADE-specific ontologies, including ODNAE for drug-associated neuropathy-inducing AEs and OCVDAE for cardiovascular drug AEs, have also been developed and used. However, these ADE-specific ontologies do not consider the effects of other factors (e.g., age and drug-treated disease) on the outcomes of ADEs. With more ontological studies of ADEs, it is also critical to develop a general purpose ontology for representing ADEs for various types of drugs. Results Our survey of FDA drug package insert documents and other resources for 224 neuropathy-inducing drugs discovered that many drugs (e.g., sirolimus and linezolid) cause different AEs given patients’ age or the diseases treated by the drugs. To logically represent the complex relations among drug, drug ingredient and mechanism of action, AE, age, disease, and other related factors, an ontology design pattern was developed and applied to generate a community-driven open-source Ontology of Drug Adverse Events (ODAE). The ODAE development follows the OBO Foundry ontology development principles (e.g., openness and collaboration). Built on a generalizable ODAE design pattern and extending the OAE and NDF-RT ontology, ODAE has represented various AEs associated with the over 200 neuropathy-inducing drugs given different age and disease conditions. ODAE is now deposited in the Ontobee for browsing and queries. As a demonstration of usage, a SPARQL query of the ODAE knowledge base was developed to identify all the drugs having the mechanisms of ion channel interactions, the diseases treated with the drugs, and AEs after the treatment in adult patients. AE-specific drug class effects were also explored using ODAE and SPARQL. Conclusion ODAE provides a general representation of ADEs given different conditions and can be used for querying scientific questions. ODAE is also a robust knowledge base and platform for semantic and logic representation and study of ADEs of more drugs in the future. Electronic supplementary material The online version of this article (10.1186/s12859-019-2729-1) contains supplementary material, which is available to authorized users.

An important need exists to better understand and stratify kidney disease according to its underlying pathophysiology in order to develop more precise and effective therapeutic agents. National collaborative efforts such as the Kidney Precision Medicine Project are working towards this goal through the collection and integration of large, disparate clinical, biological and imaging data from patients with kidney disease. Ontologies are powerful tools that facilitate these efforts by enabling researchers to organize and make sense of different data elements and the relationships between them. Ontologies are critical to support the types of big data analysis necessary for kidney precision medicine, where heterogeneous clinical, imaging and biopsy data from diverse sources must be combined to define a patient's phenotype. The development of two new ontologies - the Kidney Tissue Atlas Ontology and the Ontology of Precision Medicine and Investigation - will support the creation of the Kidney Tissue Atlas, which aims to provide a comprehensive molecular, cellular and anatomical map of the kidney. These ontologies will improve the annotation of kidney-relevant data, and eventually lead to new definitions of kidney disease in support of precision medicine.

The Coronavirus Infectious Disease Ontology (CIDO) is a community-based ontology that supports coronavirus disease knowledge and data standardization, integration, sharing, and analysis.

Ontologies are critical to data/metadata and knowledge standardization, sharing, and analysis. With hundreds of biological and biomedical ontologies developed, it has become critical to ensure ontology interoperability and the usage of interoperable ontologies for standardized data representation and integration. The suite of web-based Ontoanimal tools (e.g., Ontofox, Ontorat, and Ontobee) support different aspects of extensible ontology development. By summarizing the common features of Ontoanimal and other similar tools, we identified and proposed an “eXtensible Ontology Development” (XOD) strategy and its associated four principles. These XOD principles reuse existing terms and semantic relations from reliable ontologies, develop and apply well-established ontology design patterns (ODPs), and involve community efforts to support new ontology development, promoting standardized and interoperable data and knowledge representation and integration. The adoption of the XOD strategy, together with robust XOD tool development, will greatly support ontology interoperability and robust ontology applications to support data to be Findable, Accessible, Interoperable and Reusable (i.e., FAIR).

Tuberculosis (TB) is the leading infectious cause of death worldwide and claimed over 1.6 million lives in 2017. Furthermore, one-third of the world population is estimated to be latently infected with Mycobacterium tuberculosis (MTB). A safe and effective MTB vaccine that can prevent both the primary infection and the reactivation of latent tuberculosis infection (LTBI), and that can protect against all forms of TB in adults and adolescents is urgently needed. In this study, using computational approaches, we predicted the capacity of the epitopes to be presented by the HLA molecules for ten MTB protein antigens (Mtb39a, Mtb32a, Ag85B, ESAT-6, TB10.4, Rv2660, Rv2608, Rv3619, Rv3620, and Rv1813) constituting five MTB subunit vaccines (M72, H1, H4, H56, and ID93) that are currently in clinical trials. We also assessed the promiscuity of the predicted epitopes based on a reference set of alleles and supertype alleles, and estimated the population coverage of the ten antigens in three high TB burden countries (China, India, and South Africa). Among the ten antigens evaluated, Rv2608 was found to have the highest number of promiscuous epitopes predicted to bind the most MHC-I and MHC-II supertype alleles, highest predicted immunogenicity, and the broadest population coverage in three high burden countries. Between the two latency-related antigens (Rv1813 and Rv2660), Rv1813 was predicted to have a better epitope diversity and promiscuity, immunogenicity, and population coverage. As a result, the ID93 vaccine consisted of Rv2608, Rv1813, Rv3619, and Rv3620 was predicted to have the best potential for preventing both active and latent TB infection. Our results highlighted the importance and usefulness of a systematic and comprehensive assessment of protein antigens using computational approaches in MTB vaccine development.

Background Stem cells and stem cell lines are widely used in biomedical research. The Cell Ontology (CL) and Cell Line Ontology (CLO) are two community-based OBO Foundry ontologies in the domains of in vivo cells and in vitro cell line cells, respectively. Results To support standardized stem cell investigations, we have developed an Ontology for Stem Cell Investigations (OSCI). OSCI imports stem cell and cell line terms from CL and CLO, and investigation-related terms from existing ontologies. A novel focus of OSCI is its application in representing metadata types associated with various stem cell investigations. We also applied OSCI to systematically categorize experimental variables in an induced pluripotent stem cell line cell study related to bipolar disorder. In addition, we used a semi-automated literature mining approach to identify over 200 stem cell gene markers. The relations between these genes and stem cells are modeled and represented in OSCI. Conclusions OSCI standardizes stem cells found in vivo and in vitro and in various stem cell investigation processes and entities. The presented use cases demonstrate the utility of OSCI in iPSC studies and literature mining related to bipolar disorder.

© The Author(s) 2016. Linked Data (LD) aims to achieve interconnected data by representing entities using Unified Resource Identifiers (URIs), and sharing information using Resource Description Frameworks (RDFs) and HTTP. Ontologies, which logically represent entities and relations in specific domains, are the basis of LD. Ontobee (http://www.ontobee.org/) is a linked ontology data server that stores ontology information using RDF triple store technology and supports query, visualization and linkage of ontology terms. Ontobee is also the default linked data server for publishing and browsing biomedical ontologies in the Open Biological Ontology (OBO) Foundry (http://obofoundry.org) library. Ontobee currently hosts more than 180 ontologies (including 131 OBO Foundry Library ontologies) with over four million terms. Ontobee provides a user-friendly web interface for querying and visualizing the details and hierarchy of a specific ontology term. Using the eXtensible Stylesheet Language Transformation (XSLT) technology, Ontobee is able to dereference a single ontology term URI, and then output RDF/eXtensible Markup Language (XML) for computer processing or display the HTML information on a web browser for human users. Statistics and detailed information are generated and displayed for each ontology listed in Ontobee. In addition, a SPARQL web interface is provided for custom advanced SPARQL queries of one or multiple ontologies.

Background Aiming to understand cellular responses to different perturbations, the NIH Common Fund Library of Integrated Network-based Cellular Signatures (LINCS) program involves many institutes and laboratories working on over a thousand cell lines. The community-based Cell Line Ontology (CLO) is selected as the default ontology for LINCS cell line representation and integration. Results CLO has consistently represented all 1097 LINCS cell lines and included information extracted from the LINCS Data Portal and ChEMBL. Using MCF 10A cell line cells as an example, we demonstrated how to ontologically model LINCS cellular signatures such as their non-tumorigenic epithelial cell type, three-dimensional growth, latrunculin-A-induced actin depolymerization and apoptosis, and cell line transfection. A CLO subset view of LINCS cell lines, named LINCS-CLOview, was generated to support systematic LINCS cell line analysis and queries. In summary, LINCS cell lines are currently associated with 43 cell types, 131 tissues and organs, and 121 cancer types. The LINCS-CLO view information can be queried using SPARQL scripts. Conclusions CLO was used to support ontological representation, integration, and analysis of over a thousand LINCS cell line cells and their cellular responses.

A key aim of systems biology is the reconstruction of molecular networks. We do not yet, however, have networks that integrate information from all datasets available for a particular clinical condition. This is in part due to the limited scalability, in terms of required computational time and power, of existing algorithms. Network reconstruction methods should also be scalable in the sense of allowing scientists from different backgrounds to efficiently integrate additional data. We present a network model of acute myeloid leukemia (AML). In the current version (AML 2.1), we have used gene expression data (both microarray and RNA-seq) from 5 different studies comprising a total of 771 AML samples and a protein-protein interactions dataset. Our scalable network reconstruction method is in part based on the well-known property of gene expression correlation among interacting molecules. The difficulty of distinguishing between direct and indirect interactions is addressed by optimizing the coefficient of variation of gene expression, using a validated gold-standard dataset of direct interactions. Computational time is much reduced compared to other network reconstruction methods. A key feature is the study of the reproducibility of interactions found in independent clinical datasets. An analysis of the most significant clusters, and of the network properties (intraset efficiency, degree, betweenness centrality, and PageRank) of common AML mutations demonstrated the biological significance of the network. A statistical analysis of the response of blast cells from 11 AML patients to a library of kinase inhibitors provided an experimental validation of the network. A combination of network and experimental data identified CDK1, CDK2, CDK4, and CDK6 and other kinases as potential therapeutic targets in AML.

While effective preventive vaccines induce intended protective immunity, they also induce unintended adverse events (AEs). Generally speaking, compared to killed, inactivated vaccines and protein vaccines, live attenuated vaccines induce more protective immune responses. However, live attenuated vaccines are also associated with more AEs and even more serious AEs. For example, while live attenuated smallpox vaccines were critical to the eradication of smallpox, approximately 20–30% of smallpox vaccine recipients also experienced with various AEs that range in prevalence and severity [1]. Inter-individual variations in cytokine and AE response after smallpox vaccinations are in part due to genetic variation. For another example, the attenuated oral poliovirus vaccine (OPV) efficiently induces intestinal immunity and durable humoral immunity. However, OPV has the disadvantage of genetic instability, contributing to rare and sporadic cases of vaccine-associated paralytic poliomyelitis and the emergence of genetically divergent vaccine-derived polioviruses [2]. These AEs are worsened in patients with primary immunodeficiencies. These results suggest that the intended protective immune responses and unintended adverse events are correlated and deserve being studied simultaneously.

H56/AERAS-456+IC31 (H56), composed of two early secretion proteins, Ag85B and ESAT-6, and a latency associated protein, Rv2660, and the IC31 Intercell adjuvant, is a new fusion subunit vaccine candidate designed to induce immunity against both new infection and reactivation of latent tuberculosis infection. Efficacy of subunit vaccines may be affected by the diversity of vaccine antigens among clinical strains and the extent of recognition by the diverse HLA molecules in the recipient population. Although a previous study showed the conservative nature of Ag85B- and ESAT-6-encoding genes, genetic diversity of Rv2660c that encodes RV2660 is largely unknown. The population coverage of H56 as a whole yet remains to be assessed. The present study was conducted to address these important knowledge gaps. DNA sequence analysis of Rv2660c found no variation among 83 of the 84 investigated clinical strains belonging to four genetic lineages. H56 was predicted to have as high as 99.6% population coverage in the South Africa population using the Immune Epitope Database (IEDB) Population Coverage Tool. Further comparison of H56 population coverage between South African Blacks and Caucasians based on the phenotypic frequencies of binding MHC Class I and Class II supertype alleles found that all of the nine MHC-I and six of eight MHC-II human leukocyte antigen (HLA) supertype alleles analyzed were significantly differentially expressed between the two subpopulations. This finding suggests the presence of race-specific functional binding motifs of MHC-I and MHC-II HLA alleles, which, in turn, highlights the importance of including diverse populations in vaccine clinical evaluation. In conclusion, H56 vaccine is predicted to have a promising population coverage in South Africa; this study demonstrates the utility of integrating comparative genomics and bioinformatics in bridging animal and clinical studies of novel TB vaccines.

Hundreds of biological and biomedical ontologies have been developed to support data standardization, integration and analysis. Although ontologies are typically developed for community usage, community efforts in ontology development are limited. To support ontology visualization, distribution, and community-based annotation and development, we have developed Ontokiwi, an ontology extension to the MediaWiki software. Ontokiwi displays hierarchical classes and ontological axioms. Ontology classes and axioms can be edited and added using Ontokiwi form or MediaWiki source editor. Ontokiwi also inherits MediaWiki features such as Wikitext editing and version control. Based on the Ontokiwi/MediaWiki software package, we have developed Ontobedia, which targets to support community-based development and annotations of biological and biomedical ontologies. As demonstrations, we have loaded the Ontology of Adverse Events (OAE) and the Cell Line Ontology (CLO) into Ontobedia. Our studies showed that Ontobedia was able to achieve expected Ontokiwi features.

Cell-based therapies to treat skeletal muscle disease are limited by the poor survival of donor myoblasts, due in part to acute hypoxic stress. After confirming that the microenvironment of transplanted myoblasts is hypoxic, we screened a kinase inhibitor library in vitro and identified five kinase inhibitors that protected myoblasts from cell death or growth arrest in hypoxic conditions. A systematic, combinatorial study of these compounds further improved myoblast viability, showing both synergistic and additive effects. Pathway and target analysis revealed CDK5, CDK2, CDC2, WEE1, and GSK3β as the main target kinases. In particular, CDK5 was the center of the target kinase network. Using our recently developed statistical method based on elastic net regression we computationally validated the key role of CDK5 in cell protection against hypoxia. This method provided a list of potential kinase targets with a quantitative measure of their optimal amount of relative inhibition. A modified version of the method was also able to predict the effect of combinations using single-drug response data. This work is the first step towards a broadly applicable system-level strategy for the pharmacology of hypoxic damage.

As one of the most influential and troublesome human pathogens, Acinetobacter baumannii (A. baumannii) has emerged with many multidrug-resistant strains. After collecting 33 complete A. baumannii genomes and 84 representative antibiotic resistance determinants, we used the Vaxign reverse vaccinology approach to predict classical type vaccine candidates against A. baumannii infections and new type vaccine candidates against antibiotic resistance. Our genome analysis identified 35 outer membrane or extracellular adhesins that are conserved among all 33 genomes, have no human protein homology, and have less than 2 transmembrane helices. These 35 antigens include 11 TonB dependent receptors, 8 porins, 7 efflux pump proteins, and 2 fimbrial proteins (FilF and CAM87009.1). CAM86003.1 was predicted to be an adhesin outer membrane protein absent from 3 antibiotic-sensitive strains and conserved in 21 antibiotic-resistant strains. Feasible anti-resistance vaccine candidates also include one extracellular protein (QnrA), 3 RND type outer membrane efflux pump proteins, and 3 CTX-M type β-lactamases. Among 39 β-lactamases, A. baumannii CTX-M-2, -5, and -43 enzymes are predicted as adhesins and better vaccine candidates than other β-lactamases to induce preventive immunity and enhance antibiotic treatments. This report represents the first reverse vaccinology study to systematically predict vaccine antigen candidates against antibiotic resistance for a microbial pathogen.

Gene regulatory network inference uses genome-wide transcriptome measurements in response to genetic, environmental or dynamic perturbations to predict causal regulatory influences between genes. We hypothesized that evolution also acts as a suitable network perturbation and that integration of data from multiple closely related species can lead to improved reconstruction of gene regulatory networks. To test this hypothesis, we predicted networks from temporal gene expression data for 3,610 genes measured during early embryonic development in six Drosophila species and compared predicted networks to gold standard networks of ChIP-chip and ChIP-seq interactions for developmental transcription factors in five species. We found that (i) the performance of single-species networks was independent of the species where the gold standard was measured; (ii) differences between predicted networks reflected the known phylogeny and differences in biology between the species; (iii) an integrative consensus network which minimized the total number of edge gains and losses with respect to all single-species networks performed better than any individual network. Our results show that in an evolutionarily conserved system, integration of data from comparable experiments in multiple species improves the inference of gene regulatory networks. They provide a basis for future studies on the numerous multi-species gene expression datasets for other biological processes available in the literature., Comment: 10 pages text + 3 figures + 1 table + 2 supplementary figures + 3 supplementary tables

The BCR-ABL translocation is found in chronic myeloid leukemia (CML) and in Ph+ acute lymphoblastic leukemia (ALL) patients. Although imatinib and its analogues have been used as front-line therapy to target this mutation and control the disease for over a decade, resistance to the therapy is still observed and most patients are not cured but need to continue the therapy indefinitely. It is therefore of great importance to find new therapies, possibly as drug combinations, which can overcome drug resistance. In this study, we identified eleven candidate anti-leukemic drugs that might be combined with imatinib, using three approaches: a kinase inhibitor library screen, a gene expression correlation analysis, and literature analysis. We then used an experimental search algorithm to efficiently explore the large space of possible drug and dose combinations and identified drug combinations that selectively kill a BCR-ABL+ leukemic cell line (K562) over a normal fibroblast cell line (IMR-90). Only six iterations of the algorithm were needed to identify very selective drug combinations. The efficacy of the top forty-nine combinations was further confirmed using Ph+ and Ph- ALL patient cells, including imatinib-resistant cells. Collectively, the drug combinations and methods we describe might be a first step towards more effective interventions for leukemia patients, especially those with the BCR-ABL translocation.