key: cord-0075037-k7kitg4a authors: Jahangirian, Ehsan; Jamal, Ghadir A.; Nouroozi, MohammadReza; Mohammadpour, Alemeh title: A Novel Multiepitope Vaccine Against Bladder Cancer Based on CTL and HTL Epitopes for Induction of Strong Immune Using Immunoinformatics Approaches date: 2022-02-24 journal: Int J Pept Res Ther DOI: 10.1007/s10989-022-10380-7 sha: a4613ef08ecb3c5e04f74c2c1310372b611df288 doc_id: 75037 cord_uid: k7kitg4a Bladder cancer is well-known cancer in two forms of muscle-invasive and non-muscle-invasive bladder cancer which is responsible for annual deaths worldwide. Common therapies methods are somewhat successful; however, these methods have the limitations such as the side effects of chemotherapy which necessitate the requirement for new preventive methods against bladder cancer. Hence, we explain a novel designed multi-epitope vaccine against bladder cancer using the immunoinformatics tool. Three well-known BLCAP, PRAM, and BAGE4 antigens were evaluated due to most repetitive CTL and HTL epitopes binding. IFNγ and IL10 inducer potential of selected epitopes were investigated, as well as liner and conformational B-cell epitopes. Human beta-defensin 3 and PADRE sequence were added to construct as adjuvants, along with EAAAK, AAY, and GGGS linkers to fuse CTL and HTL epitopes. Results showed this construct encodes a soluble, non-toxic, and non-allergic protein with 70 kDa molecular weight. Modeled 3D structure of vaccine was docked whit Toll-Like Receptors (TLR) of 7/8. Docking, molecular dynamics simulation and MMBPSA analysis confirmed stability of vaccine-TLR complexes. The immunogenicity showed this construct could elicit humoral and cellular immune responses. In silico and immunoinformatics evaluations suggest that this construct is a recombinant candidate vaccine against bladder cancer. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10989-022-10380-7. Bladder cancer is common cancer worldwide known for two forms of muscle-invasive and non-muscle-invasive bladder cancer (NMIBC) (Jiang et al. 2020; Lenis et al. 2020; Slovacek et al. 2021) . Although the type of muscle-invasive bladder cancer is more dangerous, NMIBC type requires more preventative treatment because of higher recurrence rate and longer and more costly care (Shore et al. 2021 ). These preventative treatments include cystoscopies and transurethral resection (TUR), intravesical chemotherapy, and immunotherapy (Melekos and Moutzouris 2000; Fang and Huang 2009) . One successful way is using Bacillus Calmette-Guérin (BCG) vaccine intravesical immunotherapy, which decreased the risk of recurrence and progression of NMIBCs (Shelley et al. 2000) . The intrauterine injection of BCG causes extensive inflammation in the bladder wall which helps kill tumor cells, but BCG intravesical immunotherapy may have short-term immune-stimulating effects (Bevers et al. 2004 ). The recurrence of bladder cancer increases the cost of treatment, as well as the side effects of using chemotherapy drugs in the process of treating bladder cancer, today clinical studies is given in priority to discover effective, cost-effective, and accurate treatment strategies such as utilization of cytotoxicity properties of nanoparticles in the direction of suppressing the expression of cancer antigens and metastasis of cancer cells (Yin et al. 2016; DeGeorge et al. 2017; Hosseini et al. 2019 ). Ehsan Jahangirian and Ghadir A. Jamal are joint first authors (equally contributed). Page 2 of 25 Antibody therapy is one of the well-known methods against cancerous tumors. However, administration of monoclonal antibodies may cause adverse side effects due to their accumulation in non-target organs. Nikpoor et al. (2015) succeeded in using encapsulated monoclonal antibodies in Nano-liposomes as a model antibody for intravenous immunoglobulin (IVIG) in mice bearing C-26 colon carcinoma tumors, to make the efficiency of PEGylated liposomes more efficient in delivering antibodies to the tumor site than non-PEGylated liposomes (Nikpoor et al. 2015) . Non-invasive bladder cancer is dangerous in which it may not be of interest to the patient because it may not cause pain or early detection in the early stages and may progress easily, as well as the risk of developing secondary primary tumors also follows (Lenis et al. 2020) . Chemotherapy drugs are commonly used as antibiotics or anti-cancer drugs, but they often affect normal human cells which have serious side effects (Koch et al. 2021) . Moreover, cancer cells can escape chemotherapeutic agents with cellular changes (Schirrmacher 2019) . The first barrier to bladder infection and cancerous tissue cancer is asymmetric umbrella cells around the urinary tract (Romih et al. 2005 ). Due to their apical membrane, these umbrella cells form an impermeable shell covered with layers of glycan's against the penetration of bacteria and other pathogens and ultimately cause the urothelium to benefit from a highly paradigmatic mucosal defense system (Veranic et al. 2004; Khandelwal et al. 2009 ). However, a very important cellular immune response dependent on TLRs plays a key role in combating pathogens and carcinogens (Kawai and Akira 2007; Ohadian Moghadam and Nowroozi 2019) . Today, it was shown that activating TLRs by their agonists triggers immune response mediators such as cytokines, chemokines that help eliminate the cellular infection (Adams 2009; Hennessy et al. 2010; Urban-Wojciuk et al. 2019) . In general, TLR5 is poorly expressed in normal bladder cells. TLR2, TLR3, and TLR7 are moderately expressed, and finally, TLR9 and TLR4 are strongly expressed (Ohadian Moghadam and Nowroozi 2019) . Although TLRs have a courtly role in fighting cancer cells, including bladder cancer, but they act dually due to their ability to regulate immune responses (LaRue et al. 2013a) . It was shown that TLR2, TLR3, TLR4, TLR5, TLR7, and TLR9 levels of expression in normal bladder cells are very high, while their levels of expression in bladder cancer cells are reduced (So and Ouchi 2010) . But research has shown that cytokines and chemokines involved in the activation of some TLRs may contribute to bladder cancer progression (Vasekar et al. 2016; Wołącewicz et al. 2020) . The low-grade bladder cancer cell line (MGH-U3) shows similar inflammatory responses in the imitation of grade 1 NMIBC bladder cancer in decreased expression with TLR (LaRue et al. 2013b ; Ohadian Moghadam and Nowroozi 2019) . Studies showed BCG treatment of superficial urothelial cancer (UC) leaded to activation of apoptosis in UC (Yu et al. 2015) . Moreover, phagocytosis-related markers and interleukin-1 receptor-associated kinases 2 and 4 were increased when TSGH8301 cells, low-grade UC cells, were treated by BCG (Yu et al. 2015) . Increasing the expression of TLR7 by activating protein kinases activated by mitogen ERK and JNK provides the conditions that lead to tumor escape via the inhibition of T cell immunity during bladder cancer and contributes to the progression of cancer (Chi et al. 2017 ). However, the stimulation of TLR2 agonists triggers the release of pro-inflammatory cytokines that impair bladder cancer progression. Although TLRs are much more pronounced in normal urinary cells, their stable expression and activity in tumor cells make it possible to use TLR ligands for treatment. For example, the success of TLR-mediated BCG immunotherapy for NMIBCs suggests that TLR-based alternative immunotherapy may also be successful for these cancers. In many studies, the role of TLR agonists to stimulate cellular immunity was indicated. Recently, a new agonist called Beta-defensin 3 was introduced for TLRs (Funderburg et al. 2007) . Studies showed that Beta-defensin 3 can effectively activate dendritic cells (DCs) and helper T cells (Th)-1, increase the expression levels of CD40, CD80, and CD86 molecules, MHC I and II classes, and IL-preinflammatory cytokines (Tani et al. 2000; Ferris et al. 2013 ). IL-12, IL-1β, TNF-α, as well as stimulation of DC migration capacity in vitro and in vivo (Semple et al. 2010; Jiang et al. 2012; Bian et al. 2017) . This ultimately improves DC-based anti-tumor immune system and secretes CD4 ( +) and CD8 to secrete IFN-γ and mediate T cellular toxicity and help promote immunosuppressive T-cell therapy vaccines (Joly et al. 2005; Kim et al. 2018) . Several specific antigens are strongly expressed during cancers associated to the urinary system, such as Bladder Cancer (BLCA), and their introduction into the immune system can trigger a specific cellular immune system for bladder cancer (Santoni et al. 2012) . BLCA antigens are members of six bladder-specific nuclear matrix proteins (EMPs), which are known as specific urinary markers of bladder cancer (Getzenberg et al. 1996; Santoni et al. 2012) . Studies showed BLCA has a role in DNA replication and RNA synthesis of cellular functions, as well as in nuclear morphology, which leads is associated to tumor cell proliferation, survival, and angiogenesis (Berezney and Buchholtz 1981; van Eekelen and van Venrooij 1981; Seddighzadeh et al. 2003) . BLCA is expressed in carcinogenesis and tumor, and because of their high sensitivity, they are considered as markers in the diagnosis of bladder cancer (Feng et al. 2011) . PReferentially expressed Antigen of MElanoma (PRAME) is known as an atypical cancer-testis gene which is the high expression in several cancers such as melanoma, breast cancer, and leukemia. Besides, several studies indicated that testicular cancer-related antigens, such as PRAME, are also expressed in bladder cancer (Dyrskjøt et al. 2012; Kulkarni et al. 2012; Lerut et al. 2015; ). The capability of triggering an anti-tumor immune response in a melanoma, along with inducing specific killing of PRAME expressing leukemia cells, has caused PRAME to have great immunogenicity. Several studies confirmed the expression of PRAME in urothelial carcinoma and nonmuscle-invasive tumors (Dyrskjøt et al. 2012) . It has shown that cancer-testis genes are expressed in bladder cancer, particularly associated with high-grade tumors (Simpson et al. 2005) . Utilizing CT antigens such as PRAME as bladder antigen cancer vaccine candidates may ultimately induce an immune response in specific immunotherapy against bladder cancer. We used antigen, expressed in non-muscle-invasive bladder carcinomas to stimulate and identify a specific immune response to bladder cancer. BAGE is an antigen associated with the human B melanoma antigen and encodes a putative protein of 43 amino acids that are found to be expressed in different tumors (Boël et al. 1995) . BAGEL has a pattern of expression similar to the MAGE, GAGE, and LAGE/NY-ESO-1 families, which are highly expressed in urinary cancers (De Plaen et al. 1994; Lethé et al. 1998; De Backer et al. 1999) . Immunohistochemistry has identified the antigens in cancers of histological origin, especially bladder cancer and melanoma (Amin 2009; Bellizzi 2020) . Moreover, two studies have shown that BAGE is involved in cancer/testis expression profile and is also expressed in bladder carcinoma simultaneously (Boël et al. 1995; Ruault et al. 2002) . Using traditional methods of vaccine production is very costly and time-consuming, in addition to the need for repeated tests and errors in the laboratory. Today, it is possible to predict epithelial B-cells and T cells using immunoinformatics tools (Safavi et al. 2019 ). These tools help achieve a recombinant vaccine that accurately stimulates the cell's immune response (Validi et al. 2018; Beg et al. 2021) . The use of molecular binding studies to predict or confirm experimental results plays an important role in understanding the function of protein-protein, protein-peptide, protein-DNA or protein-ligand such as Drug-protein complexes and also provides a better basis for the production of more efficient and newer drugs and vaccines as well as experimental and in silico approaches for researchers (Dareini et al. 2020; Sadeghzadeh et al. 2020) . Molecular modeling approaches and determining the binding site and binding affinity in order to monitor and identify the binding sites between two protein or peptide complexes, can help to consciously and accurately select the desired epitopes to be used in the vaccine structure or design of peptide inhibitors (Danesh et al. 2018; Mokaberi et al. 2021 ). We designed a recombinant multi-epitope vaccine based on the potential of cellular immune response stimulation of three proprietary bladder cancer antigens, and TLR agonists, using immunoinformatics tools and then evaluated vaccine structure via protein-protein docking and molecular dynamics (MD) simulation. This study aimed to design a multiepitope recombinant vaccine based on bladder cancer antigens to investigate the specific stimulation of cellular immune system against bladder cancer. Uniprot database-base was used to retrieve HUMAN B melanoma antigen 4 (BAGE4), HUMAN Bladder cancerassociated protein (BLCAP), and HUMAN Melanoma antigen preferentially expressed in tumors (PRAME) proteins based on anti-cancerous properties. The NetMHC 4.0 server (http:// www. cbs. dtu. dk/ servi ces/ NetMHC/) was used to prediction of peptide-MHC class I binding molecules. This server is harnessing Artificial Neural Networks (ANN) capable of performing sensitive, quantitative predictions of peptide binding to the MHC class I molecule. The HLA-A, HLA-B and HLA-C MHC class I molecules were select for binding to antigens. We used the MHC-II-peptide binding affinity prediction methods to predict of MHC class II alleles. This server (http:// www. cbs. dtu. dk/ servi ces/ NetMH CII/) uses an extended data set of quantitative MHC-peptide binding affinity data obtained from the Immune Epitope Database. The Human alleles including DR, DP, and DQ were used for the prediction of MHC class II alleles. The AAY linker bonded CTL epitopes. The GGGS linker also fused HTL epitopes. The EAAAK linker was used to bind Human Beta-defensin 3 and BAGE4. Also, Beta-defensin and PADRE adjuvants were added to the N-terminal and C-terminal region of the construct, respectively. Finally, the amino acid structure of the bladder cancer vaccine was formed. The theoretical properties of the construct were estimated using ProtParam (http:// web. expasy. org/ protp aram). Furthermore, PsiPred online software (http:// bioinf. cs. ucl. ac. uk/ psipr ed/) was used to predicts the secondary structure of the protein. PsiPred utilizes the self-optimized prediction method (SOPM) that defined based on the amino acid sequence. The rate of construct allergenicity was determined by Alg-Pred (http:// www. imtech. res. in/ ragha va/ algpr ed/). This server predicts the allergenicity of a protein with a sensitivity accuracy of about 85% (Saha and Raghava 2006) . The predicting of the antigenic value of the construct was determined by the VaxiJen v2.0 server on (http:// www. ddg-pharm fac. net/ vaxij en/ VaxiJ en. html). VaxiJen v2.0 proposes a new alignment-free approach for antigen prediction to overcome the limitations of alignment-dependent methods, which is based on auto cross-covariance (ACC) transformation of protein sequences into uniform vectors of principal amino acid properties (Doytchinova and Flower 2007) . The solubility of protein was evaluated via the SOLpro server (http:// prote omics. ics. uci. edu). SOLpro is a server for protein secondary structure prediction based on protein evolutionary information (sequence homology) and homologous protein's secondary structure (structure homology) which achieves a performance exceeding 79% correctly classified residues on proteins with no homologs in the PDB and exceeding 92% correctly classified residues on proteins where homologs can be found in the PDB (Magnan et al. 2009 ). Java Utility Server (JCat) (http:// www. prodo ric. de/ JCat) was utilized for codon optimizing the recombinant cancer vaccine construct. The pET-26b(+) was selected as a proteinexpressing vector. The sequence was optimized for expression in a negative gram bacteria E. coli K12. The Nde I and Xho I restriction enzymes were selected for restriction sites in the SnapGene tool. The B-Cells epitopes were predicted using the IEDB server (https:// www. iedb. org/). This server has collected methods to predict linear B-cell epitopes based on sequence characteristics of the antigen using amino acid scales and HMMs. B-Cells epitopes are one of the most important factors contributing to the vaccine design because they are characterized by the immune system. All discontinuous B cell epitopes were determined using the ElliPro server. The importance of these epitopes is due to their decisive role in the spatial form of antigenic structures (Ponomarenko et al. 2008) . We developed a model with a prediction accuracy of more than 85% and Area Under the Curve (AUC) 0.9 via the ElliPro which uses the amino acid composition as an input feature for the Support Vector Machine (SVM). GalaxyPepDock server was used to evaluate the interaction between the selected epitopes and HLA class I (Ko et al. 2012 ). Since HLA-A6823 was repeated for 12 of the total HLA class I epitopes of both BCAP and PRAME antigens (in Table 1 ), HLA-A6823-epitope was selected for docking (Crooke et al. 2020) . The 3D structure of HLA-A6823 was modeled with 67% similarity to the crystal structure of HLA-B*1501 (RCSB code: 1xr9) and then docked with SMFMGFYLL, CTICALVFL, LWFSHSMFM, HSMFMG-FYL, CTICALVFL, TLAKFSPYL, HVMNPLETL, YIAQFTSQF, KAMVQAWPF, FPPLFMAAF, RTFYDPEPI and LPRELFPPL epitopes. Results was evaluated based on protein structure similarity (TM-score) and Interaction similarity score. Ten models of each peptide-HLA complex were generated on the basis of minimized energy scores, and the top model for each complex was selected for comparative analysis. I-Tasser modeled the 3D structure of the construct (Yang and Zhang 2015) . I-TASSER server is an online platform that implements algorithms for protein structuring and function predictions. It allows us to automatically generate high-quality model predictions of the 3D structure and biological function of protein molecules from their amino acid sequences. The modeled construct of the cancer vaccine was refined by https:// zhang lab. ccmb. med. umich. edu/ ModRe finer. The Ramachandran plot graph and Z-Score were used to validate of modeled construct in http:// molpr obity. bioch em. duke. edu/ and https:// www. came. sbg. ac. at/ prosa. php, respectively. We used the crystallographed structures of TLR2 and TLR7 which were selected based on their interaction with the bladder cancer antigens. Also, we selected the binding site TLR's receptor based on the binding site belonging to I-TASSER server information. Each of the TLR2 (6NIG PDB code) and TLR7 (7CYN PDB code) PDBs structures were retrieved from the RCSB server and saved. Then, structures of the TLR PDBs were minimized energy in the MVD software. All of the additional ligands and water molecules were deleted. In the next step, the cancer vaccine construct was introduced to MVD software and all the steps for docking preparation were performed like TLRs structures. Protein-protein docking was performed by the HADDOCK server. HADDOCK distinguishes itself from ab-initio docking methods in that it encodes information from identified or predicted protein interfaces in Ambiguous Interaction Restraints (AIRs) to drive the docking process. We defined amino acids identified in the previous step for interaction HADDOCK. Then, we selected the best cluster which had the most negative HADDOCK score after docking. MD simulation was utilized as an important technique for analyzing and investigating the protein-protein complex validation. GROMACS v.5.1.4 (http:// manual. groma cs. org/ docum entat ion/5. 1.4/ index. html) was used to simulation of the construct-TLRs complexes. The system was neutralized by the addition of an appropriate concentration of sodium or chloride ions. Energy minimization was performed based on which this process was a descent steeped algorithm. System equilibration was performed for NVT and NPT ensemble that the temperature was then brought up to 300 K and 1 bar pressure. Simulation of complexes were carried out for the protein-protein by 40 ns time scale to analyze the trajectories after the equilibration process. Contributions for the binding free energy of complexes were calculated by MM-PBSA analysis. Solvation free energies was determined via solving the linearised Poisson Boltzman equation for reaching the electrostatic contribution in the solvation free energy. Average interaction energy between receptor, ligand and taking the entropy change upon binding into account was used to calculates ΔGvacuum. Simulation of the immunological response excitation for constructed bladder vaccine was validated by the C-ImmSim server (http:// 150. 146.2. 1/C-IMMSIM/ index. php). Three major functional mammal system components (bone marrow, thymus, and lymph node) are simulated in this server. The ability to simulate various types of immune cells such as CTL, HTL, B-cells, dendritic cells and cytokines was considered for this candidate bladder vaccine when it is made. The immune simulation was performed taking into account the minimum clinical interval between doses of two (1) ΔGbind, solv = ΔGbind, vaccum + ΔGsolv, complex − (ΔGsolv, ligand + ΔGsolv, complex) (2) ΔGsol = Gelecterostatic = 80 − Gelecterostatic = 1 + ΔGhydrophobic (3) ΔGvaccum = ΔEmolecular mechanics − T. ΔSnormal mode analysis All three proteins were stored in the FASTA form after retrieval in the Uniprot. NetMHC 4.0 server was used to predict MHC-I epitopes. As can be seen in Table 1 , 20 Cytotoxic T-cell Lymphocyte epitopes (9-mer) were selected for HLA-A, HLA-B, and HLA-C of MHC Class-I alleles molecules based on most repeatability among strong binds. Those T-cell epitopes that had recurrence and overlap were identified and selected along with other T-cell epitopes for the construct. All MHC-I epitopes are included in supplementary Table S1 . The helper T lymphocyte (HTL) epitopes were predicted using NetMHC II 2.3 server. During this step, 10 Human MHC-II allele molecules, including DR, DP, and DQ were selected based on the most repeatability among strong binds. Each of the HTL epitopes molecules has defaulted to 15mer (with 9mer core) ( Table 2 ). All investigated epitopes related to helper T lymphocyte are shown in Supplementary Table S2 . We used complete BAGE4 protein sequence. The design of multiepitope bladder vaccine used 5 important factors. These items including 3 proteins, along with Beta-defensin 3 and PADRE adjuvants. Different linkers were used to fuse each of peptides. Beta-defensin 3 and PADRE were bonded using EAAAK linker. The connection among CTL peptides was performed using the AAY linker. Furthermore, the connection among HTL peptides was created via the GGGS linker (Fig. 1 ). The results of ProtParam server showed that the bladder construct vaccine is protein, including 743 amino acids, the molecular weight of 82.18 kDa, theoretical isoelectric (pI) of about 9.44, Aliphatic Index of 83.15, and an Instability Index of 50.5. These results confirmed that the construct is a stable protein. Moreover, the Grand average of hydropathicity of this bladder peptide vaccine is 0.341. The construct has a half-life of about 30 h in mammalian reticulocytes (in vitro), > 20 h in yeast (in vivo), and > 10 h in Escherichia coli (in vivo). The secondary structure of bladder cancer construct is included in Supplementary Fig. S1 . Two Xho I and Nde I restriction enzymes were selected for in silico cloning nucleotide sequence with the length of 2229 bp. The mean GC of sequences and codon compatibility index (CAI) of the nucleotide sequences was 52.26% and 1.0, respectively. Finally, because E. coli bacterium K12 will be amused as an expression host, the nucleotide sequence was cloned inside the pET26b (+) vector. HisTag was considered in the C-terminal sequence, and the final vaccine sequence was designed by SnapGene software (Fig. 2 ). I-Tasser server modeled the vaccine construct. Then, the modeled construct was submitted to the ModRefiner server for refining. TM score difference between the raw and refined structure was 0.46 angstrom, respectively. Using the ProSA server, the Z-score of construct was determined about −6.35 before refinement, while this value was obtained −6.9 after refinement (Fig. 3A, B) . Ramachandran graph details indicated most of the amino acids which make up the model are in favored and allowed areas (%86.7 in favored regions, and 95.9% in allowed regions) (Fig. 3C, D) . Favored regions allowed regions and outlier areas were calculated 91.9% and 98.2%, after refinement, respectively. Superimposition of bladder cancer vaccine structure, before and after the refinement, respectively (Fig. 3E ). B-cells epitopes of final construct were predicted using the IEDB server. As shown in Table 3 , all predicted epitopes had an antigenicity greater than 0.4. In the next step, conformational B-cells epitopes for each of the peptides were predicted using the ElliPro server. Conformational B-cells The HLA class I candidate peptides were predicted to provide the highest coverage in the HLA-A6823 allele. Because the HLA-A6823 was the only candidate to show the most connection to class 1 epitopes (with 12 repetitions), with a relatively high affinity. Therefore, we performed the molecular binding studies of each peptide with the molecular structure of HLA-A6823. The results showed that all peptides bind in the peptide binding groove (Fig. 4 and Table 4 ). Binding of epitopes in the hydrophobic binding cavity of the HLA groove supports the hypothesis that these peptides Fig. 3 Modeling, refinement and validation of bladder cancer vaccine 3D structure. A and B Z-Score of bladder cancer vaccine 3D structure, before and after structure refinement. C and D Ramachandran graph related to bladder cancer vaccine 3D structure, before and after structure refinement. E Superimposition of modeled to bladder cancer vaccine 3D structure. Blue structure is referred to before refinement and yellow structure is referred to after refinement bind strongly to HLA-A6823 and could be promising epitope candidates for vaccine development studies and if used in the structure of the vaccine, they can induce an appropriate immune response. In the first step, the raw 3D structure of bladder cancer vaccine was minimized by MVD software. Then, the bladder cancer vaccine and TLRs were introduced to HADDOCK server for docking. Each protein-protein docking was separately performed (vaccine-TLR2 and vaccine-TLR7). HAD-DOCK score in the bladder vaccine-TLR2 complex was obtained − 40 (Table 5 ) ( Fig. 5A-C) , moreover, the docking results of HADDOCK server showed the best complex with a HADDOCK score of -104 there are in the vaccine-TLR7 complex (Table 6 ) ( Fig. 6A-C) . All of the interactions between vaccines and TLRs are shown as 2D in Fig. 7A , B. The MD simulation method performs interactions between the protein/protein, such as vaccine/receptor analysis. The structural integrity of vaccine-TLR2 and vaccine-TLR7 docked complexes were simulated via 35 ns by MD. Then, MD's results, such as root-mean square-deviations (RMSD) was used followed for the analysis of hydrogen bonds interaction between bladder vaccine structure and TLRs. As shown in Fig. 8A , RMSD graph relates to the vaccine-TLR7 complex indicated an upward trend; however, after 5 ns, it has almost stabilized. Vaccine-TLR2 and vaccine-TLR7 complexes had an RMSD average value of about 2.5 and 2.2 nm after the complexes stabilized, respectively. RMSD results is inferred that the vaccine-TLR2 and vaccine-TLR7 complexes remained stable via 35 ns simulations. Root mean-square-fluctuations (RMSF) were calculated to further calculations of residual and side-chain flexibility of complexes (Figs. 8B and 9B ). Almost no obvious fluctuations in complexes, nonetheless fluctuations at the amino acid interaction regions up to 0.4 nm, but the other complexes residues remained highly stable in their fluctuations. The results of MD simulation endorse that referring to the HADDOCK docking, this vaccine could strongly evoke a significant immune response against TLRs related to bladder cancer. Furthermore, the Radius of gyration related to vaccine-TLR2 and vaccine-TLR7 complexes (~ 1.5-1.4) showed the potential of stability maintenance in complexes (Figs. 8C and 9C). Also, H-bond graphs in vaccine-TLR2 and vaccine-TLR7 complexes are remained unchanged via the MD, indicating that the complexes remain stable (Figs. 8D and 9D). MD results were more evaluated by MM/BPSA analysis to determine the energy contribution of involved amino acids in each complex. As shown in Table 8 , the energy contribution of vaccine-TLR7 complex is higher than vaccine-TLR2 (Table 7) . Besides, further analysis indicates van der Waal energy is negative in both complexes, but vaccine-TLR7 complex shows higher negative electrostatic energy compare to vaccine-TLR2. In general, MM/BPSA analysis showed a better interaction in vaccine-TLR7 complex, although overall results of docking, MD simulation, and MM/BPSA analysis confirm high stable interaction for vaccine-TLR2 and vaccine-TLR7 complexes. Real-life phenomena adequation was achieved based on the response of in silico mediated immune. As shown in Fig. 10A , the primary reaction was lower than both secondary and tertiary reactions, which follows rapid antigen clearance leads to the distinguishing of the immunoglobulin movement greater such as IgM, IgG1 + IgG2, and IgG + IgM antibodies. Significant development of memory cells is shown that is relevant to a higher activation level of B-cells with IgM and IgG1 biotype (Fig. 10B, C) . In the next context, the number of elicited active T-cells by the bladder vaccine was enhanced throughout the secondary or tertiary reactions, but gradually decreased after that (Fig. 10D-G) . Bladder cancer vaccine simulations lead to evoke immune responses as well as significantly increased levels of Despite significant advances in cancer treatment today, cancer is still one of the leading causes of death worldwide. Most cancer treatment regimens for the treatment of choice for advanced or metastatic disease steps have originated focus primarily on chemotherapy (Schirrmacher 2019) . The overall quality of life in the patients undergoing chemotherapy is decreased under affecting the wide range of side effects generated by high toxicity and induce from Current chemotherapeutic drugs (Yin et al. 2016) . Furthermore, a wide range from chemotherapeutic agents is identified and expelled by those cancer cells which display a high tendency to develop resistance against chemotherapeutic compounds (Schirrmacher 2019) . Traditional methods of designing a cancer vaccine are costly and time-consuming tasks requiring many experimental studies (Fatoba et al. 2022 ). Today, due to the advent of computational biology methods such as genomics, proteomics, and immunoinformatics tools, the ability to construct, predict, and interact a virtual cancer vaccine with immune system molecules has been achieved. It is not possible to analyze these benefits in traditional vaccine production methods. The rapid development of computational methods and models, along with a vast body of experimental data to analyze the immune system response using immunological data, has led to the simulation of immune system response to the cancer vaccine (Aminnezhad et al. 2020) . Using cancer vaccines is more effective than current treatments for cancer, including chemotherapy and radiation, especially for treating melanoma-related cancers such as bladder cancer. Activating and stimulating the immune system against specific antigens to elicit an anti-tumor immune response using a vaccine is the great purpose of immunization against cancer. Numerous studies showed that using different types of antigens associated with melanoma-related tumors, has yielded excellent results for vaccination (Boël et al. 1995; Ruault et al. 2002 Ruault et al. , 2003 Weber et al. 2011; Al-Khadairi and Decock 2019; Safavi et al. 2019; Ellingsen et al. 2021; Ouyang et al. 2021; Richard et al. 2021) . We used BLCAP, PRAM, and BAGE4 antigens as key factors in our construct. Regarding the expression of BLCAP, PRAM, and BAGE4 in malignant tissues, these antigens could potentially be used as appropriate targets for anticancer immunization (Moreira et al. 2010; Gordeeva 2018; Lezcano et al. 2018; Lohman et al. 2021) . Although so far there is no report that simultaneously examines the interaction between the three antigens BLCAP, PRAM and BAGE4 in physiological conditions and the interaction behavior among them in vivo; However, a number of studies have shown that each of the three antigens BLCAP, PRAM, and BAGE4 may have an interaction partner in bladder or testis cancer separately (Cho et al. 2006; Chen et al. 2017; Gromova et al. 2017) . For example, Gromova et al. (2017) showed that BLCAP and Stat3 are in close physical proximity of each other in bladder tissue, and that BLCAP physically interacts with Stat3, so that BLCAP is a novel Stat3 interaction partner and suggested a role for BLCAP in the Stat3-mediated progression of precancerous lesions to invasive tumors of the bladder. BLCAP plays a role as a novel STAT3 interaction partner in bladder cancer. Functional association for each of the antigens can be associated with concomitant and same direction responses to elicit immune responses (Gromova et al. 2017) . The selected region was fused to the CTL epitopes by GGGS linker to facilitate the delivery of epitopes via MHC processing and effective fragment separation. CTLs play a vital role in the cell's immune response to cancer by inducing apoptosis in tumor cells as well as producing interferongamma (Yong et al. 2012 ). The stimulation of CTL epitopes was the most important parameter to select BLCAP, PRAM, and BAGE4 antigens. An appropriate vaccine should be able to stimulate helper T lymphocyte (HTL) responses to achieve maximum efficacy. HTL-related responses can play a key role to regulate a wide range from different immune system responses. Stimulation of HTL can trigger the production of memory T 20 CTL epitopes and 20 HTL epitopes derived from BLCAP and PRAM antigens were fused by AAY and GGGS linkers, respectively. Besides, we considered and BAGE4 sequence as an independent antigen. In addition to activating HTL, these linkers maintain vaccine-dependent immunogenicity for the CTL response. Beta-defensin 3 and PADER sequence selected in this study are considered the detectable adjuvants for the Toll-like receptor (TLR) that activate the cellular immune response (Funderburg et al. 2007; Ferris et al. 2013) . The integration of these adjuvants into the structure of cancer vaccine could increase its effectiveness, as TLRs are important in identifying pathogens by various immune cells, including B cells, dendritic cells, and macrophages (Adams 2009 ). The activation of TLRs by their specific adjuvants can increase the secretion of inflammatory cytokines and the over-regulation of MHC molecules. The activation of TLRs, in addition to differentiating simple HTL cells from T-helper1, leads to the activation of CTLs (Kawai and Akira 2007) . Therefore, choosing the right agonist to stimulate TLR can greatly enhance the success of the immune response elicited by the cancer vaccine. Beta-defensin 3 is a well-known agonist for stimulating TLR complex. The interaction of Beta-defensin with TLRs leads dendritic cells to T-helper1 cells and subsequently stimulates cellular immune responses (Safavi et al. 2019 ). The Beta-defensin was fused to the BAGE4 via EAAAK linkers. EAAAK, as a helical structure, was placed in the C-terminal of Betadefensin to reduce its interaction with other peptide vaccine domains. The presence of an EAAAK linker can spatially separate B and T cell epitopes and ultimately lead to proper vaccine interaction with TLRs. B cells affect dendritic cells, natural killer cells, and myeloid-suppressing cells that express Fc receptors by producing antibodies and forming antigen-antibody complexes (Xu-Monette and Young 2020) . Furthermore, the mechanism of antibody-dependent cell-mediated cytotoxicity is one of the most important strengths of B cells in the elimination of cancer cells (Zappasodi et al. 2015; Xu-Monette and Young 2020) . Since the selected MHC-I and MHC-II epitopes provided maximum antigen coverage for the vaccine, we selected B-cell epitopes from the final sequence of the vaccine to prevent the vaccine structure from over-enlarging. As expected, the overall antigenicity of B cell epitopes from the final structure of the vaccine was very significant. Considering the importance of the fact that B cells are very important in responding to the multi-peptide peptide vaccine, 12 linear B cell epitopes were selected from the final construct that showed the highest antigenicity. Moreover, conformational B cell epitopes were identified on the amino acid sequence of the vaccine. Considering this strategy for the cancer vaccine by determining a large number of B cell epitopes leads to an effective humoral immune response in the body. The high solubility of structure after protein expression increases stability and immunogenicity. The final structure is codon-optimized to maximize expression in the prokaryotic host. The final DNA structure was integrated into the pET26b vector to in silico cloning of bladder cancer vaccine. In addition, studies have shown that protein folding is influenced by physicochemical properties, so understanding intermolecular relationships, such as protein complexes, can improve the prediction of molecular structures (Chamani and Moosavi-Movahedi 2006) . Regarding bladder cancer immunotherapy, Galluzzi et al. (2012) showed that TLR2 signaling is strongly stimulated through BCG therapy (Galluzzi et al. 2012) . This finding confirms that appropriate stimulation of TLR2 may occur if antigenspecific antigens are used. Phase II studies of the CDX-1307 bladder cancer vaccine have also demonstrated the activation of TLR signaling pathways in patients with newly diagnosed bladder cancer (Morse et al. 2011) . Alignment of laboratory results and bioinformatics simulations can increase the reliability of studies. Application of molecular modeling in the study of Zare-Feizabadi et al. (2021) in the interaction of calf thymus DNA (ct DNA) with anastrozole (which is acknowledged as an antineoplastic drug), supported and confirmed anti-cancer results of anastrozole use on MCF 7 cell line (Zare-Feizabadi et al. 2021 ). TLR2 represents a double-edged sword whose role in cancer must be carefully understood in order to launch effective therapies. Because it is expressed in various malignancies and promotes their proliferation through activation of the myeloid differentiation primary response protein 88 (MyD88)/nuclear factor kappalight-chain-enhancer of activated B cell (NF-κB) pathway. Hence, to activate immune responses, TLR2 is considered as a player in anti-cancer immunity. Therefore, TLR2 agonists have been used as adjuncts to anti-cancer immunotherapies. Anticancer inhibitors or agonists induce apoptosis or death of cancer cells by activating immune pathways (Moosavi-Movahedi et al. 2003; Di Lorenzo et al. 2020) . Docking results showed a high affinity between Bladder cancer vaccine structure and both of TLR2/TLR7. Then docking results were evaluated by MD simulation to confirmation of structure interactions. As it was identified, MD simulation confirmed the stability of both bladder vaccine-TLR2 and bladder vaccine-TLR7 by RMSD, RMSF, Rg graphs. Furthermore, MM/PBSA analysis indicated high negative energy contributions of amino acids in complexes, as well as electrostatics and Van der Waals energy. Using docking techniques and MD simulations, the results of immune responses at the cellular and molecular levels can be used before laboratory tests to evaluate cancer vaccines. Using prevention methods, such as multi-epitope vaccines containing antigens expressed in bladder cancer can provide effective immunity against bladder cancer. The MHC-I and MHC-II epitopes selected from BLCAP and PRAM antigens had the highest coverage in the final structure of the vaccine. B cell epitopes were selected and identified from the final structure of the vaccine that had high antigenicity to induce hemorrhagic immunity. In addition, the interaction of HLA class I molecules with MHC-I epitopes indicated that highaffinity peptides were expected to be used to develop specific vaccines against bladder cancer. As shown in our study, the bladder cancer vaccine candidate can establish strong interactions with prognostic molecules such as TLR2 and TLR7. Therefore, studies simulating the immune responses given to this vaccine showed that both humoral and cellular immune responses are evoked. The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s10989-022-10380-7. Toll-like receptor agonists in cancer therapy Cancer testis antigens and immunotherapy: where do we stand in the targeting of PRAME? Histological variants of urothelial carcinoma: diagnostic, therapeutic and prognostic implications Immunoinformatics design of multivalent chimeric vaccine for modulation of the immune system in Pseudomonas aeruginosa infection Designing multi-epitope vaccine candidates against functional amyloids in Pseudomonas aeruginosa through immunoinformatic and structural bioinformatics approach An algorithmic immunohistochemical approach to define tumor type and assign site of origin Dynamic association of replicating DNA fragments with the nuclear matrix of regenerating liver Role of urothelial cells in BCG immunotherapy for superficial bladder cancer Human β-defensin 3 reduces TNF-α-induced inflammation and monocyte adhesion in human umbilical vein endothelial cells BAGE: a new gene encoding an antigen recognized on human melanomas by cytolytic T lymphocytes Effect of n-alkyl trimethylammonium bromides on folding and stability of alkaline and acid-denatured cytochrome c: a spectroscopic approach A-to-I RNA editing of BLCAP lost the inhibition to STAT3 activation in cervical cancer Antitumor activity of toll-like receptor 7 agonists Physical interaction of two cancer-testis antigens, MAGE-C1 (CT7) and NY-ESO-1 (CT6) Immunoinformatic identification of B cell and T cell epitopes in the SARS-CoV-2 proteome Determining the binding site and binding affinity of estradiol to human serum albumin and holo-transferrin: fluorescence spectroscopic, isothermal titration calorimetry and molecular modeling approaches A novel view of the separate and simultaneous binding effects of docetaxel and anastrozole with calf thymus DNA: experimental and in silico approaches Characterization of the GAGE genes that are expressed in various human cancers and in normal testis Structure, chromosomal localization, and expression of 12 genes of the MAGE family Bladder cancer: diagnosis and treatment Toll-like receptor 2 at the crossroad between cancer cells, the immune system, and the microbiota VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines Expression of MAGE-A3, NY-ESO-1, LAGE-1 and PRAME in urothelial carcinoma Telomerase as a target for therapeutic cancer vaccines and considerations for optimizing their clinical potential. Front Immunol Intravesical drug delivery into the bladder to treat cancers Transmission of antibiotic-resistant Escherichia coli from chicken litter to agricultural soil Urinary BLCA-4 is highly specific for detection of bladder cancer in Chinese Han population and is related to tumour invasiveness Human beta-defensin 3 induces maturation of human langerhans cell-like dendritic cells: an antimicrobial peptide that functions as an endogenous adjuvant Human-defensin-3 activates professional antigen-presenting cells via Toll-like receptors 1 and 2 Trial watch: experimental toll-like receptor agonists for cancer therapy Bladder cancer-associated nuclear matrix proteins Cancer-testis antigens: Unique cancer stem cell biomarkers and targets for cancer therapy Identification of BLCAP as a novel STAT3 interaction partner in bladder cancer Targeting Toll-like receptors: emerging therapeutics? Nano-curcumin's suppression of breast cancer cells (MCF7) through the inhibition of cyclinD1 expression Antiviral activity of recombinant mouse β-defensin 3 against influenza A virus in vitro and in vivo Trimodality therapy for muscle-invasive bladder cancer: recent advances and unanswered questions Correlation between beta-defensin expression and induction profiles in gingival keratinocytes Cell biology and physiology of the uroepithelium Human β-defensin 2 plays a regulatory role in innate antiviral immunity and is capable of potentiating the induction of antigen-specific immunity GalaxyWEB server for protein structure prediction and refinement Side effects of intravesical BCG and chemotherapy for bladder cancer: what they are and how to manage them Cancer/testis antigens and urological malignancies Toll-like receptors in urothelial cells-targets for cancer immunotherapy Toll-like receptors in urothelial cells-targets for cancer immunotherapy Bladder cancer: a review Rates of MAGE-A3 and PRAME expressing tumors in FFPE tissue specimens from bladder cancer patients: potential targets for antigen-specific cancer immunotherapeutics LAGE-1, a new gene with tumor specificity PRAME expression in melanocytic tumors The utility of PRAME staining in identifying malignant transformation of melanocytic nevi SOLpro: accurate sequence-based prediction of protein solubility Intravesical therapy of superficial bladder cancer Analysis of the interaction behavior between Nano-Curcumin and two human serum proteins: combining spectroscopy and molecular stimulation to understand protein-protein interaction Design, synthesis, and anticancer activity of phosphonic acid diphosphate derivative of adenine-containing butenolide and its water-soluble derivatives of paclitaxel with high antitumor activity Bladder cancer-associated protein, a potential prognostic biomarker in human bladder cancer CDX-1307: a novel vaccine under study as treatment for muscle-invasive bladder cancer A novel multi-epitope peptide vaccine against cancer: an in silico approach Nanoliposome-mediated targeting of antibodies to tumors: IVIG antibodies as a model Toll-like receptors: the role in bladder cancer development, progression and immunotherapy Antitumor effects of iPSC-based cancer vaccine in pancreatic cancer ElliPro: a new structure-based tool for the prediction of antibody epitopes Multi-step screening of neoantigens' HLA-and TCR-interfaces improves prediction of survival Differentiation of epithelial cells in the urinary tract New BAGE (B melanoma antigen) genes mapping to the juxtacentromeric regions of human chromosomes 13 and 21 have a cancer/testis expression profile BAGE genes generated by juxtacentromeric reshuffling in the Hominidae lineage are under selective pressure Characterizing the binding of angiotensin converting enzyme I inhibitory peptide to human hemoglobin: influence of electromagnetic fields In silico analysis of transmembrane protein 31 (TMEM31) antigen to design novel multiepitope peptide and DNA cancer vaccines against melanoma AlgPred: prediction of allergenic proteins and mapping of IgE epitopes Pathogenic and diagnostic potential of BLCA-1 and BLCA-4 nuclear proteins in urothelial cell carcinoma of human bladder From chemotherapy to biological therapy: a review of novel concepts to reduce the side effects of systemic cancer treatment (review) Low IL-1alpha expression in bladder cancer tissue and survival Human beta-defensin 3 has immunosuppressive activity in vitro and in vivo Intravesical bacillus calmette-guerin in Ta and T1 bladder cancer Non-muscle-invasive bladder cancer: an overview of potential new treatment options Cancer/testis antigens, gametogenesis and cancer Approaches to non-muscleinvasive bladder The application of Toll like receptors for cancer therapy Defensins act as potent adjuvants that promote cellular and humoral immune responses in mice to a lymphoma idiotype and carrier antigens The role of TLRs in anti-cancer immunity and tumor rejection. Front Immunol Immuno-informatics based approaches to design a novel multi epitope-based vaccine for immune response reinforcement against Leptospirosis Van Venrooij WJ (1981) hnRNA and its attachment to a nuclear protein matrix Immunotherapy in bladder cancer What determines differentiation of urothelial umbrella cells? A phase 1 study of a vaccine targeting preferentially expressed antigen in melanoma and prostate-specific membrane antigen in patients with advanced solid tumors Immunotherapy in bladder cancer: current methods and future perspectives The role of the cancer testis antigen PRAME in tumorigenesis and immunotherapy in human cancer Therapeutic Vaccines for Aggressive B-Cell Lymphoma I-TASSER server: new development for protein structure and function predictions Neoadjuvant chemotherapy for muscle-invasive bladder cancer: a systematic review and two-step meta-analysis Strategies for enhancing vaccine-induced CTL antitumor immune responses Bacille Calmette-Guerin can induce cellular apoptosis of urothelial cancer directly through toll-like receptor 7 activation Lymphoma immunotherapy: current status. Front Immunol Determining the interaction behavior of calf thymus DNA with anastrozole in the presence of histone H1: spectroscopies and cell viability of MCF-7 cell line investigations Author Contributions Conceptualization, EJ and GhAJ; methodology, EJ, GhAJ, MN and AM; formal analysis, EJ and GhAJ; resources, MN and AM; data curation, EJ, GhAJ, MN and AM; writing-original draft preparation, EJ, GhAJ, MN and AM; writing-review and editing, EJ and GhAJ; project administration, EJ and GhAJ; funding acquisition, EJ and GhAJ. All authors have read and agreed to the published version of the manuscript. Conflict of interest None.