key: cord-0047951-x2ehhe0s
authors: Muhammad, Syed Aun; Ashfaq, Hiba; Zafar, Sidra; Munir, Fahad; Jamshed, Muhammad Babar; Chen, Jake; Zhang, Qiyu
title: Polyvalent therapeutic vaccine for type 2 diabetes mellitus: Immunoinformatics approach to study co-stimulation of cytokines and GLUT1 receptors
date: 2020-07-23
journal: BMC Mol Cell Biol
DOI: 10.1186/s12860-020-00279-w
sha: cffb835ae23577b52b077d14dc448557cf952495
doc_id: 47951
cord_uid: x2ehhe0s

BACKGROUND: Type 2 diabetes mellitus (T2DM) is a worldwide disease that have an impact on individuals of all ages causing micro and macro vascular impairments due to hyperglycemic internal environment. For ultimate treatment to cure T2DM, association of diabetes with immune components provides a strong basis for immunotherapies and vaccines developments that could stimulate the immune cells to minimize the insulin resistance and initiate gluconeogenesis through an insulin independent route. METHODOLOGY: Immunoinformatics based approach was used to design a polyvalent vaccine for T2DM that involved data accession, antigenicity analysis, T-cell epitopes prediction, conservation and proteasomal evaluation, functional annotation, interactomic and in silico binding affinity analysis. RESULTS: We found the binding affinity of antigenic peptides with major histocompatibility complex (MHC) Class-I molecules for immune activation to control T2DM. We found 13-epitopes of 9 amino acid residues for multiple alleles of MHC class-I bears significant binding affinity. The downstream signaling resulted by T-cell activation is directly regulated by the molecular weight, amino acid properties and affinity of these epitopes. Each epitope has important percentile rank with significant ANN IC(50) values. These high score potential epitopes were linked using AAY, EAAAK linkers and HBHA adjuvant to generate T-cell polyvalent vaccine with a molecular weight of 35.6 kDa containing 322 amino acids residues. In silico analysis of polyvalent construct showed the significant binding affinity (− 15.34 Kcal/mol) with MHC Class-I. This interaction would help to understand our hypothesis, potential activation of T-cells and stimulatory factor of cytokines and GLUT1 receptors. CONCLUSION: Our system-level immunoinformatics approach is suitable for designing potential polyvalent therapeutic vaccine candidates for T2DM by reducing hyperglycemia and enhancing metabolic activities through the immune system.

Type 2 diabetes mellitus (T2DM), a non-insulin dependent metabolic disorder, is a pandemic disease affecting large set population of the world [1] . It is responsible for 90% of total diabetic population and sixth prime cause of disability. T2DM is characterized by the inability of pancreatic β-cells to produce enough insulin resulting hyperglycemia and the inability of insulin to bind with its receptors restrict the absorption of glucose (insulin resistance) into the cells [2, 3] . Disease prevalence is increasing for sure due to unknown causes and the lack of therapeutic strategies [4, 5] .

Normally, insulin regulates the absorption of glucose through glucose transporter type 4 (GLUT-4) protein channels [6] present in cell membranes. In case of insulin impairment, the absorption of glucose by GLUT-4 doesn't occur, causing T2DM [7] . The glucose transporters (GLUT1 and GLUT4) facilitate glucose transport into cells. GLUT1 is insulinindependent and is widely distributed in different tissues [8, 9] . Cells need growth factors to facilitate glucose absorption for subsistence and development. Tcell stimulation leads to fast proliferation and differentiation into effector cells that release cytokines and mediate the immunity [10] [11] [12] .

The non-insulin growth factors such as cytokines including interleukin IL3 and IL7 may absorb the glucose through glucose transporter type 1 (GLUT-1) proteins. These cytokines can trigger the cascade of important signals to promote glucose uptake via different pathways. GLUT1 proteins in response to cytokines released by active immune components facilitate the constitutive, insulin-independent glucose uptake in most of the cells including hematopoietic and muscle cells [13] [14] [15] .

The progression of T2DM is also associated to abnormal immune responses [16] , and therefore the cytokinemediated regulation of GLUT-1 can be thought of playing some significant role in this respect. The metabolic reprogramming is shaped to help definite cell functions [17] and glucose uptake delivers a key metabolic control point through the GLUT family of facilitative glucose transporters.

In this study, we hypothesize the development of potential immunotherapeutic vaccine candidates for the activation and secretion of cytokines (IL-1, IL-3, and IL-7) to facilitate glucose absorption and cure T2DM. The issues related to insulin resistance could be minimized through alternative non-insulin dependent GLUT1 pathway. We predicted the T-cell epitopes and analyzed the in-silico binding affinity with MHC class-I molecules. Our vaccine would target the T-cells resulting the secretion of interleukins. Instead of insulin dependent GLUT-4 channels, these interleukins open up the GLUT-1 proteins and regulate the glucose absorption. Our hypothesis has been illustrated in Fig. 1 . This study would modulate the therapeutic strategies to manage type 2 diabetes mellitus.

The tissue specific (bearing GLUT-1 receptors) protein sequences were retrieved from NCBI and Uniprot databases (Supplementary Table 1 ). The data covered all information including protein names, gene symbols, Uniprot accession numbers, protein description and sequences. The currently available proteins sequences associated to type 2 diabetes mellitus (Supplementary  Table 2 ) were accessed from diabetic databases. We carried out this study using integrated framework (Fig. 2 ) by computational tools, databases, online servers and software ( Table 1) .

We mapped the non-tissue specific list of diabetic proteins with the entire list of tissue specific proteins using "Compare Two Lists" tool [23] to shortlist the T2DM associated proteins in the relevant tissues. The antigenicity of these shortlisted proteins were determined (based on the significant threshold level > 0.5) using "VaxiJen v2.0" server [47] . To select the robust and effective antigenic proteins, these proteins were further screened and filtered based on the molecular weight (threshold of > 85 kDa) using Expasy Compute pI/MW tool [24, 50] , hydrophilicity using SPpred server, amino acid composition using Expasy ProtParam tool [24] and subcellular localization using CELLO v.2.5: Subcellular Localization Predictor tool [25] .

We predicted multi-allelic T-Cells epitopes of selected proteins using ProPred-I online server [24] . This server identified the MHC Class-I regions in the input sequences of selected proteins by applying matrices of all 47 alleles of this class. These epitopes have the potential to target several alleles of MHC Class-I. The multiallelic epitopes-based vaccine has more worth and chance of success as compared to the vaccine that target only one type of allele in the whole population. The proteasomal cleavage sites of these antigenic proteins were predicted using NetChop predictor of IEDB tool [27] . The proteasomal prediction is important to find out the potential immunogenic regions in the selected proteins. We verified the immunogenicity of MHC Class-I epitopes using IEDB server. This tool uses amino acid properties and their position within the peptide to find the immunogenicity of a peptide MHC (pMHC) complex [27] . 

The conservational analysis of selected epitopes was analyzed using IEDB server [27] to observe the conservancy between epitopes and the proteins. We predicted the physicochemical properties of candidate peptides including half-life, instability index, net charge, peptide toxicity, and hydropathicity using ToxinPred and PredSTP servers [20, 21] . The half-life of peptide varies among species and it is important to determine the half-life of peptides in humans. The stability profile of peptide is indicated by the value of the instability index. The charge of a peptide is pH-dependent. At the isoelectric point, the net charge on protein is zero. The solubility of protein is lowest at its isoelectric point. The hydropathicity analysis provides detail information about individual amino acid that enables us to predict the overall threedimensional structure of a protein from its amino acid sequence. The Support Vector Machine (SVM) model of amino acid composition was used to predict the toxicity profile of peptides.

The three-dimensional (3D) structure of peptides is usually helpful to understand its topological description, biological activity and function. The protein structure prediction often offers a suitable alternative to facilitate structure-based studies. 3D models of selected epitopes were generated by using PEP-FOLD [29] , I-Tasser [34] and Chimera [37] tools to observe epitopic regions and underlying pattern of amino acids. These tools require amino acid input sequences to build the peptide folds.

Gene enrichment and annotation analysis of these peptides were performed using FunRich [38] and David tool [30] . The annotation profiles including functional and cluster details were studied.

We studied PPI network of selected antigenic proteins to carryout system level investigations. The target interactions of selected source proteins were retrieved using STRING database [32] . This database provides the comprehensive detail about the interactions, functions and pathways of sample proteins. The PPI network was constructed by using Cytoscape v3.6.0 software [33] .

To study the physiological role of these potential vaccinating proteins, we designed and constructed the integrated, interactive and metabolic network of T2DMrelated proteins and observed the correlation between these pathways. Cellular and signaling pathways were reconstructed from the combined gene signatures using PathVisio3 tool. These proteins were mapped and curated using KEGG (Kyoto Encyclopedia of Genes and Genomes) and WIKI pathways on the basis of literature and database evidence.

Polyvalent vaccine was designed and constructed by linking selected epitopes of MHC Class-I using potential linkers. We used AAY and EAAAK linkers to link 13-epitopes to minimize the undesirable attachment of their ends that might cause change in amino acid arrangement and even the functionality of proteins [51] . To enhance the immunogenicity, the amino acid sequence (159 residues) of the heparin-binding hemagglutinin (HBHA) was used as adjuvant. On both ends of HBHA, EAAAK linkers were used to make it non-reactive [51] . Epitopes cannot be bluntly linked with each other to design a polyvalent vaccine. The arrangement of amino acid sequences is critical and their order is based on their affinity and compatibility with each other. To design the best and compatible construct of selected epitopes, we used the HADDOCK 2.2 server [26] to analyze the binding affinity of these epitopes. Initially, the binding affinity of each epitope was evaluated followed by the combined sequences to shortlist the best construct. The refinement and residues determination were carried out using HADDOCK: Refinement Interface and Cport modules [26] respectively.

The polyvalent vaccine construct was modelled using I-Tasser server [34] . The input FASTA format of amino acid sequence was used to generate 3D model. The model was refined using 3D Refine tool. The quality of the model was estimated by Ramachandran plot and quality model energy analysis (QMEAN) score. Ramachandran plot accurately describes the protein conformation and illustrates the favorable regions for backbone dihedral angles against amino acid residues in protein structure. Similarly, QMean score describes the protein quality on the basis of different geometrical aspects of its structure.

The binding affinity of polyvalent vaccine model with MHC class-I molecule was analyzed using Molecular Operating Environment (MOE) software. The MOE software was set at default parameters and the interaction between molecules was visualized. The active binding sites of the MHC molecule was observed and the binding interaction between the amino acid residues was assessed based on binding energy.

From the list of 2601 diabetic proteins, we identified 13antigenic tissue specific proteins associated with T2DM based on antigenicity, molecular weight, subcellular localization, amino acid composition, length of protein sequence and solubility using successive screening tools (Fig. 3a) . These selected extracellular and membrane bound proteins are associated to blood brain barrier, pancreas, muscles, lymphocytes and intestines (Fig. 3b) . These proteins like interleukin-32 (IL-32), insulin like growth factor 1 (IGF1), transforming growth factor beta-1 (TGFB1), toll like receptors-3 (TLR3) and ras-related C3 botulinum toxin substrate 1 (RAC1), with their important biological functions including immune modulation, insulin signaling, cell survival, immune components activation and phosphorylation during glucose metabolism ( Table 2) .

T-cell epitopes (9-mers in length) of 13-selected proteins were predicted. These MHC class-I specific epitopes are multi-allelic in nature and could target several alleles of MHC Class-I of human population. We determined 13immunogenic T-cell epitopes with significant percentile rank (% rank: 0.1), MHC Class-I immunogenicity score (> 0.5) and multiple alleles hit. These selected T-cell epitopes of 13-proteins as potential candidate antigens have been listed in Table 3 .

We applied NetChop tool to predict proteasomal cleavage that depends upon a neural network. This method identified the C-terminal at cleavage sites with the threshold value of 0.5 to categorize the cleavage and non-cleavage sites. We observed the cleavage (positive predictions) and non-cleavage (negative predictions) sites of antigenic proteins by proteasomes indicating their significant role in antigen presentation to MHC class-I molecule (Fig. 4) . The scores of combined predictions of proteasomal cleavage, TAP translocations, and MHC binding shows each peptide intrinsic capacity of being a T-cell epitope. Epitope conservation analysis of antigenic peptides showed broad coverage and 100% protein sequence matches. Conservation analysis was performed by IEDB and epitopes were designed by ProPred-I identified the MHC Class-I regions in the input sequences of selected proteins by applying matrices of all 47 alleles of this class. Those epitopes were selected that are present in maximum alleles out of 47 reported alleles of MHC class I. These epitopes have the potential to target several alleles of MHC Class-I. The multi-allelic epitopes-based vaccine has more worth and the chance of success as compared to the vaccine that targets only one type of allele in the whole population. The conserved epitopes have significant ability to provide the effective immune response ( Table 4 ). Most of the predicted epitopes have instability index < 40 indicating their stability. A good epitope must be stable and stay in the body to activate the immune components. Most of our predicted peptide has a half-life of more than 5 h and all are non-toxic. SVM score is used to classify either the epitope is toxic or non-toxic based on the threshold value (non-toxic: < 5). Almost all predicted epitopes have an SVM score of less than 5 indicating non-toxicity. Amino acid chains are either hydrophobic or hydrophilic. The lesser value of hydropathicity means the more hydrophilic character of epitopes (Table 5) . 

We predicted the 3D models of 13-epitopes from input 9-mers amino acid sequences using PEP-FOLD server.

In local structure analysis, the probabilities of each Structural Alphabet (SA) on vertical and horizontal axis has been shown (Fig. 5) . In heat map Figure, red color codes indicate the helical form of the structure, green presents extended and blue indicates coil conformations. This server found lowest energy conformations with an average RMSD of 2.1 Å.

The functional annotation and enrichment of the T2DM associated proteins were analyzed using DAVID and Fun-Rich tools. The functional annotation identifies protein functional domains, disease associations, protein-protein interactions and biological pathways. This cluster analysis showed that these proteins are significantly associated with cell-cell communication, glycosylation, cell surface responses and cell proliferation regulation (Fig. 6a) . The biological role of these potential antigenic proteins has been importantly seen in immune cell migration, wound healing, cell mobility, cell communication, anti-apoptosis and lipid metabolism. Their transcription factors involve NKX2-1, ELF1. ZNF513, ZNF238, GLI1 and IRF1 with significant pvalues (p < 0.05) expressed in different tissues (Fig. 6b) .

To understand and analyze the topology and functional annotation of protein-protein interaction (PPI) of product of many signal transduction pathways stimulated by inflammation, apoptosis and cell growth) and LRP6 (LRP6_HUMAN) interacts with GSK3B (regulator of glucose homeostasis) (Fig. 7) . Such interaction of antigenic proteins makes them effective to regulate the immune responses and metabolic pathways.

The pathways modeling showed that several pathways are involved in T2DM-pathophysiology and they are connected with our core vaccinating proteins. This analysis indicated the association of antigenic protein with immune system and cellular signaling mechanism. Besides insulin-signaling and insulin-resistance pathways, the others are AMPK signaling, JAK/STAT pathway, FOXO signaling, P13K-AKT signaling, and WNT pathways are associated to our potential antigenic proteins and responsible for the activation of immune components and GLUT-1 receptors to regulate T2DM. It was found that some proteins like RAC1_Human, IGF1_Human, ITB1_Human were associated with the activation of P13K-AKT pathway which is the main regulator of many cellular processes like cell growth, survival and glucose metabolism. In insulin resistance, the FOXO signaling pathway activates GLUT-1 and inhibition of GLUT-4 by IKKB. The immune-regulatory mechanism and increased production of cytokines is associated with the activation of JAK/STAT pathway and GLUT-1 receptors to control glucose metabolism and T2DM (Fig. 8) . 

For the construction of final vaccine construct, predicted T-cell epitopes were linked to each other with help of suitable and flexible AAY linkers. HBHA (1-159 residues) as adjuvant was added at N-terminals of final constructs with the help of EAAAK linker (Fig. 9a) . After adding linkers and adjuvant, the final construct contained 322 amino acid residues with the molecular weight of 35.6 kDa and 3D-model was generated (Fig. 9b) . Our polyvalent epitope model showed significant antigenicity (> 0.5) determined by Vaxijen server. Analysis of Ramachandran plots indicated that 68.8% residues of polyvalent protein model (alpha-helices) were in the most favored regions with the 98% expectation while 2% in the allowed region. The Ramachandran plot is demonstrating low energy conformations for ϕ (phi) and (psi) angles of the model. The graphical representation shows the local backbone conformation of each residue and more than 68% residues are in favorable region (Fig. 9c) . The model was refined and the quality of the construct was estimated by QMEAN score (− 15.9) and C-score (− 3.62) (Fig. 9d ).

We analyzed the binding affinity of polyvalent vaccine construct with MHC class-I molecule using MOE software. The lowest binding energies confirmed the optimum binding patterns between these molecules. The interaction of poly-epitopes (ligand) and MHC class-I (target) have been shown in Fig. 10a . We found the following binding energy of polyvalent construct with MHC class-I molecule: − 15.34 Kcal/mole (Illustration Table of Fig. 10a ). It has been observed that polyvalent interacted with ARG, PRO, GLU, TRP, ASP, GLY, THR, LYS, ALA, PHE, TYR, GLN, MET, HIS and LEU amino acid residues of MHC class-I molecule (Fig. 10b) . This interaction would prove that our polyvalent would be effective therapeutic vaccine candidate for T2DM.

Type 2 diabetes mellitus is a metabolic complex disease caused by several factors that is characterized by relative lack of obesity, insulin, insulin resistance, and high blood sugar [52, 53] . Although new therapies and management protocols to cure T2DM are being introduced, still the exact causes and effective treatment is imprecise. Currently, we are using different drug molecules to maintain blood glucose level and cure T2DM [54, 55] . These medications only minimize the symptoms of hyperglycemia and provide short term control. As a result, the drug resistance is developed in most of the diabetic patients and in the while, disease progress so fast that drugs alone cannot encounter it [56] . The current study signifies the importance of immunotherapies and vaccine development to modify the therapeutic strategies. Therapeutic vaccines could be effective to manage T2DM [57] as immune components have been found interactive with glucose transporter proteins [58] . It has been observed that the blood serum cytokines level elevates in diabetic conditions and in turn cytokines opens up the insulin independent glucose transporter channels and cells absorb the glucose [59] . The significant association of interleukin-1β (IL-1β) with T2DM is reported [53] and the role of antibodies has been studied in diabetic conditions [60] .

The evidences about the association of GLUT-1 with immune components are available. It has been revealed that macrophages critically require GLUT-1 for their inflammatory response [61] . Also, the lymphocytes functioning demands high levels of energy that is assisted by GLUT-1 and in turn the expression of GLUT-1 remarkably increases with the increase of immune cells activation [17] . CD4 T-cells and T-effector cells also depends on GLUT-1 for their survival and programming regarding glucose metabolism and the GLUT-1 deficiency causes impairment of both glycolysis and T cells survival [15] . IL3 has reported to be involved in the trafficking and recycling of GLUT-1 intracellularly and controlling its expression on the cell surface, thus, maintaining the glucose uptake [62] . Cytokines like IL-1β enhances the expression of GLUT-1 and increase the glucose uptake in human articular chondrocytes [63] .

Cytokines have also been known to affect and upregulate the GLUT-1 expression significantly by increasing the nitric oxide production levels [64] . Cytokines inhibit the tyrosine kinase activity which downregulates the insulin dependent GLUT-4 and upregulates the expression of insulin independent GLUT-1 [65] . CD28 alone as well as by secreting other immune mediators increases the GLUT-1 induction considerably [66] . TNFA (Tumor Necrosis Factor Alpha) has a cause and effect bond with GLUT-1 and plays a key role in the stabilization of GLUT-1 [67] . Similarly, stimulation of CD46 in CD4 T cells of humans causes the elevated expression of GLUT-1 on cells surfaces [68] . Numerous methods and databases for developing vaccines or immunotherapy against various pathogenic diseases have been developed over the last 20 years. T-cell epitope prediction methods that include indirect techniques such as Major Histocompatibility Complex prediction and transporterassociated protein binders [69] . Endogenous proteins may promote the hypercreativity and autoimmune reactions, however few studies showed that predicted epitopes of endogenous proteins exhibited significant, safe and effective results in animals. Recently, the vaccines for metabolic diseases have made considerable progress, particularly in the treatment of dyslipidemia, atherosclerosis, diabetes mellitus and hypertension, but comprehensive studies are required before any clinical applications [70] . Since T2DM is a multifactorial disorder, T2DM therapeutic-vaccine has been designed to predict obesity protein antigens [71] . Such vaccine targets include adipose tissue antigens, somatostatin, glucose-dependent insulinotropic polypeptide (GIP), and ghrelin [71] . In the studies of Zhang et al., 2018 and Zha et al., 2016 reported that cytokine L-1ß is a key proinflammatory substance in T2DM pathogenesis and has shown a reduction in weight gain, improved glucose tolerance and insulin sensitivity, and a lower β-cell loss in their vaccine with predicted peptide epitopes [55, 72] . Similarly, in phase I/II clinical trials vaccine Hil1bQb targeting IL-1β was found safe and well-tolerant [73] . Dipeptidyl peptidase 4 (DPP4) has been studied as an inhibitor of the glucagon-like peptide-1 (GLP-1) glucosedependent insulinotropic peptide (GLP) [74] , a therapeutic vaccine against DPP4 that has shown important, safe and successful results in the control of glucose levels in the mouse using GLP-1 [75] . In-silico method for the prediction of antigenic peptides and polyvalent construction could be used to develop vaccines for other human diseases. Molecular docking analysis of the refined vaccine structure with different MHC molecules and human immune TLR-2 receptor proved significant interaction [76] . In this case, the sequences of those proteins would be required against which we have to design a vaccine. This strategy could be effectively applied for the development of polyvalent vaccine candidates for other diseases specifically cancer and autoimmune diseases.

This study focused to design potential T-cell polyepitopes to simulate the immune components to provoke glucose transporters to control hyperglycemia. Our 

Pathophysiology of type 1 and type 2 diabetes mellitus: a 90-year perspective

The relative contributions of insulin resistance and beta-cell dysfunction to the pathophysiology of type 2 diabetes

Prevalence and risk factors of microalbuminuria among type 2 diabetes mellitus: a hospital-based study from

A review of clinical efficacy and safety of canagliflozin 300 mg in the management of patients with type 2 diabetes mellitus

Examining trends in type 2 diabetes incidence, prevalence and mortality in the UK between

Regulation of insulin resistance by multiple MiRNAs via targeting the GLUT4 signalling pathway

Akt/PKB activation and insulin signaling: a novel insulin signaling pathway in the treatment of type 2 diabetes

Insulin resistance in type 2 (non-insulin-dependent) diabetic patients and their relatives is not associated with a defect in the expression of the insulin-responsive glucose transporter (GLUT-4) gene in human skeletal muscle

Astrocytic insulin signaling couples brain glucose uptake with nutrient availability

T cell metabolism drives immunity

Transcriptional control of effector and memory CD8+ T cell differentiation

Production of transforming growth factor beta by human T lymphocytes and its potential role in the regulation of T cell growth

effects of varying degrees of intermittent hypoxia on proinflammatory cytokines and adipokines in rats and 3T3-L1 adipocytes

Placental glucose transporter (GLUT)-1 is down-regulated in preeclampsia

The glucose transporter Glut1 is selectively essential for CD4 T cell activation and effector function

Selfish brain and selfish immune system interplay: a theoretical framework for metabolic comorbidities of mood disorders

Glucose metabolism in lymphocytes is a regulated process with significant effects on immune cell function and survival

T2D@ ZJU: a knowledgebase integrating heterogeneous connections associated with type 2 diabetes mellitus

DAPD: a knowledgebase for diabetes associated proteins

Open Source Drug Discovery Consortium, et al. In Silico approach for predicting toxicity of peptides and proteins

PredSTP: a highly accurate SVM based model to predict sequential cystine stabilized peptides

UniProt: the universal protein knowledgebase

Translocon Declogger Ste24 protects against IAPP oligomer-induced Proteotoxicity

Protein identification and analysis tools on the ExPASy server

Predicting subcellular localization of proteins for gram-negative bacteria by support vector machines based on n-peptide compositions

The HADDOCK2. 2 web server: user-friendly integrative modeling of biomolecular complexes

The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage

Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines

A hidden markov model derived structural alphabet for proteins

AIM2 recognizes cytosolic dsDNA and forms a caspase-1-activating inflammasome with ASC

HAPPI-2: a comprehensive and high-quality map of human annotated and predicted protein interactions

The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible

Cytoscape: a software environment for integrated models of biomolecular interaction networks

I-TASSER server for protein 3D structure prediction

I-TASSER: a unified platform for automated protein structure and function prediction

The I-TASSER suite: protein structure and function prediction

UCSF chimera-a visualization system for exploratory research and analysis

A novel community driven software for functional enrichment analysis of extracellular vesicles data

Verification of protein structures: patterns of nonbonded atomic interactions

QMEAN: a comprehensive scoring function for model quality assessment

Assessing the local structural quality of transmembrane protein models using statistical potentials (QMEANBrane)

Structure validation by Cα geometry: ϕ, ψ and Cβ deviation

3Drefine: an interactive web server for efficient protein structure refinement

i3Drefine software for protein 3D structure refinement and its assessment in CASP10

3Drefine: consistent protein structure refinement by optimizing hydrogen bonding network and atomic-level energy minimization

SCRATCH: a protein structure and structural feature prediction server

VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines

Identifying candidate subunit vaccines using an alignment-independent method based on principal amino acid properties

Bioinformatic approach for identifying parasite and fungal candidate subunit vaccines

The focusing positions of polypeptides in immobilized pH gradients can be predicted from their amino acid sequences

A novel multiepitope peptide vaccine against cancer: an in silico approach

Mechanisms of pancreatic β-cell death in type 1 and type 2 diabetes: many differences, few similarities

Preclinical efficacy and safety of an anti-IL-1β vaccine for the treatment of type 2 diabetes

Pathophysiology and treatment of type 2 diabetes: perspectives on the past, present, and future

Therapeutic vaccine against IL-1β improved glucose control in a mouse model of type 2 diabetes

The IL-1 pathway in type 2 diabetes and cardiovascular complications

GAD treatment and insulin secretion in recent-onset type 1 diabetes

Prospects of a type 1 diabetes vaccine

IL-7 promotes Glut1 trafficking and glucose uptake via STAT5-mediated activation of Akt to support T-cell survival

Studies on glycoxidatively modified human IgG: implications in immuno-pathology of type 2 diabetes mellitus

Metabolic reprogramming of macrophages glucose transporter 1 (GLUT1)-mediated glucose metabolism drives a proinflammatory phenotype

Cytokine stimulation promotes glucose uptake via phosphatidylinositol-3 kinase/Akt regulation of Glut1 activity and trafficking

Cytokine regulation of facilitated glucose transport in human articular chondrocytes

Cytokines modulate glucose transport in skeletal muscle by inducing the expression of inducible nitric oxide synthase

Plasma interleukin-6, tumour necrosis factor α and blood cytokine production in type 2 diabetes

The CD28 signaling pathway regulates glucose metabolism

Tumor necrosis factor alphainduced glucose transporter (GLUT-1) mRNA stabilization in 3T3-L1 preadipocytes. Regulation by the adenosine-uridine binding factor

CD46 activation regulates miR-150-mediated control of GLUT1 expression and cytokine secretion in human CD4+ T cells

Novel in silico tools for designing peptide-based subunit vaccines and Immunotherapeutics

Recent advances in the development of vaccines for diabetes, hypertension, and atherosclerosis

Obesity vaccines

Interleukin-1β-targeted vaccine improves glucose control and β-cell function in a diabetic KK-ay mouse model

Development of an interleukin-1β vaccine in patients with type 2 diabetes

Differentiating among incretin therapies: a multiple-target approach to type 2 diabetes

Therapeutic vaccine against DPP4 improves glucose metabolism in mice

Contriving a chimeric polyvalent vaccine to prevent infections caused by herpes simplex virus (Type-1 and Type-2): an exploratory immunoinformatic approach

Publisher's Note

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

Supplementary information accompanies this paper at https://doi.org/10. 1186/s12860-020-00279-w.Additional file 1.

Authors' contributions AH, ZS and MS designed the study. AH, and ZS performed the experiments. CJ, JM, QZ, and MS analyzed and interpreted the data. AH, MS and MF wrote the manuscript. All the authors read and approved the final article.

Not applicable.

The datasets generated during the current study have been included in manuscript and available as additional files.Ethics approval and consent to participate Not applicable.

Not applicable.

The authors declare that they have no competing interests.Author details 1 Institute of Molecular Biology and Biotechnology, Bahauddin Zakariya