Identification of HLA-I restricted epitopes in six vaccine candidates of Leishmania tropica using immunoinformatics and molecular dynamics simulation approaches

Engineering of methionine-auxotroph Escherichia coli via parallel evolution of two enzymes from Corynebacterium glutamicum's direct-sulfurylation pathway enables its recovery in minimal medium

Methionine biosynthesis relies on the sequential catalysis of multiple enzymes. Escherichia coli, the main bacteria used in research and industry for protein production and engineering, utilizes the three-step trans-sulfurylation pathway catalyzed by L-homoserine O-succinyl transferase, cystathionine gamma synthase and cystathionine beta lyase to convert L-homoserine to L-homocysteine. However, most bacteria employ the two-step direct-sulfurylation pathway involving L-homoserine O-acetyltransferases and O-acetyl homoserine sulfhydrylase. We previously showed that a methionine-auxotroph Escherichiacoli strain (MG1655) with deletion of metA, encoding for L-homoserine O-succinyl transferase, and metB, encoding for cystathionine gamma synthase, could be complemented by introducing the genes metX, encoding for L-homoserine O-acetyltransferases and metY, encoding for O-acetyl homoserine sulfhydrylase, from various sources, thus altering the Escherichia coli methionine biosynthesis metabolic pathway to direct-sulfurylation. However, introducing metX and metY from Corynebacterium glutamicum failed to complement methionine auxotrophy. Herein, we generated a randomized genetic library based on the metX and metY of Corynebacterium glutamicum and transformed it into a methionine-auxotrophic Escherichia coli strain lacking the metA and metB genes. Through multiple enrichment cycles, we successfully isolated active clones capable of growing in M9 minimal media. The dominant metX mutations in the evolved methionine-autotrophs Escherichia coli were L315P and H46R. Interestingly, we found that a metY gene encoding only the N-terminus 106 out of 438 amino acids of the wild-type MetY enzyme is functional and supports the growth of the methionine auxotroph. Recloning the new genes into the original plasmid and transforming them to methionine auxotroph Escherichia coli validated their functionality. These results show that directed enzyme-evolution enables fast and simultaneous engineering of new active variants within the Escherichia coli methionine direct-sulfurylation pathway, leading to efficient complementation.

Advancements and hurdles in the development of a vaccine for triple-negative breast cancer: A comprehensive review of multi-omics and immunomics strategies

Triple-Negative Breast Cancer (TNBC) presents a significant challenge in oncology due to its aggressive behavior and limited therapeutic options. This review explores the potential of immunotherapy, particularly vaccine-based approaches, in addressing TNBC. It delves into the role of immunoinformatics in creating effective vaccines against TNBC. The review first underscores the distinct attributes of TNBC and the importance of tumor antigens in vaccine development. It then elaborates on antigen detection techniques such as exome sequencing, HLA typing, and RNA sequencing, which are instrumental in identifying TNBC-specific antigens and selecting vaccine candidates. The discussion then shifts to the in-silico vaccine development process, encompassing antigen selection, epitope prediction, and rational vaccine design. This process merges computational simulations with immunological insights. The role of Artificial Intelligence (AI) in expediting the prediction of antigens and epitopes is also emphasized. The review concludes by encapsulating how Immunoinformatics can augment the design of TNBC vaccines, integrating tumor antigens, advanced detection methods, in-silico strategies, and AI-driven insights to advance TNBC immunotherapy. This could potentially pave the way for more targeted and efficacious treatments.

CRM197-conjugated peptides vaccine of HCMV pp65 and gH induce maturation of DC and effective viral-specific T cell responses

Human cytomegalovirus (HCMV) infection is prevalent worldwide, and there is currently no licenced HCMV vaccine to control it. Therefore, developing an effective HCMV vaccine is a significant priority. Because of their excellent immunogenicity, the crucial components of HCMV, phosphoprotein 65 (pp65) and glycoproteins H (gH) are potential target proteins for HCMV vaccine design. In this study, we predicted and screened the dominant antigenic epitopes of B and T cells from pp65 and gH conjugated with the carrier protein cross-reacting material 197 (CRM197) to form three peptide-CRM197 vaccines (pp65-CRM197, gH-CRM197, and pp65-CRM197+gH-CRM197). Furthermore, the immunogenicity of the peptide-CRM197 vaccines and their effects on dendritic cells (DCs) were explored. The results showed that three peptide-CRM197 vaccines could induce maturation of DCs through the p38 MAPK signalling pathway and promote the release of proinflammatory factors, such as TNF-α and interleukin (IL) -6. Meanwhile, the peptide-CRM197 vaccines could effectively activate T cell and humoral immunity, which were far better than the inactivated HCMV vaccine. In conclusion, we constructed three peptide-CRM197 vaccines, which could induce multiple immune effects, providing a novel approach for HCMV vaccine design.

Multi-epitope vaccine candidates based on mycobacterial membrane protein large (MmpL) proteins against Mycobacterium ulcerans

Buruli ulcer (BU) is a neglected tropical disease. It is caused by the bacterium Mycobacterium ulcerans and is characterized by skin lesions. Several studies were performed testing the Bacillus Calmette-Guérin (BCG) vaccine in human and animal models and M. ulcerans-specific vaccines in animal models. However, there are currently no clinically accepted vaccines to prevent M. ulcerans infection. The aim of this study was to identify T-cell and B-cell epitopes from the mycobacterial membrane protein large (MmpL) proteins of M. ulcerans. These epitopes were analysed for properties including antigenicity, immunogenicity, non-allergenicity, non-toxicity, population coverage and the potential to induce cytokines. The final 8 CD8+, 12 CD4+ T-cell and 5 B-cell epitopes were antigenic, non-allergenic and non-toxic. The estimated global population coverage of the CD8+ and CD4+ epitopes was 97.71%. These epitopes were used to construct five multi-epitope vaccine constructs with different adjuvants and linker combinations. The constructs underwent further structural analyses and refinement. The constructs were then docked with Toll-like receptors. Three of the successfully docked complexes were structurally analysed. Two of the docked complexes successfully underwent molecular dynamics simulations (MDS) and post-MDS analysis. The complexes generated were found to be stable. However, experimental validation of the complexes is required.

Two peptides derivate from Acinetobacter baumannii outer membrane protein K as vaccine candidates: a comprehensive in silico study

BACKGROUND

The lack of appropriate vaccines is an obstacle to the effective management of A. baumannii infections. Peptide vaccines offer an attractive and promising preventive strategy against A. baumannii.

OBJECTIVE

In this study, we identified specific T cell epitopes of A. baumannii outer membrane protein K (OMPK) using comprehensive bioinformatics and detailed molecular docking analysis.

METHODS

Both class-I and class-II T cell epitopes of A. baumannii OMPK were predicted by three tools namely IEDB, SYFPEITHI, and ProPred. The predicted epitopes were shortlisted based on several analyses including prediction scoring, clustering, exclusion of human similarity, considering immunogenicity and cytokine production, and removal of toxic and/or allergen epitopes. The epitopic peptides with high prediction scores and appropriate properties containing both class-I and class-II T cell epitopes were selected. Two of these class I/II epitopic peptides were chosen for molecular docking studies and assessing their physicochemical properties as vaccine candidates.

RESULTS

The results showed many T-cell epitopes of OMPK that could be evaluated for possible immunogenicity. Two of these epitopes (containing both class-I and II epitopes) had high prediction scores, were predicted by several tools, attached to several HLAs, and had the best docking score. They had different physicochemical properties and were conserved among Acinetobacter species.

DISCUSSION

We identified the A. baumannii OMPK high immunogenic class-I and class-II T cell epitopes and introduced two promising high immunogenic peptides as vaccine candidates. It is recommended to perform in vitro/in vivo investigation of these peptides to determine their true efficacy and efficiency.

A New Strategy for Mapping Epitopes of LACK and PEPCK Proteins of Leishmania amazonensis Specific for Major Histocompatibility Complex Class I

Leishmaniasis represents a complex of diseases with a broad clinical spectrum and epidemiological diversity, considered a major public health problem. Although there is treatment, there are still no vaccines for cutaneous leishmaniasis. Because Leishmania spp. is an intracellular protozoan with several escape mechanisms, a vaccine must provoke cellular and humoral immune responses. Previously, we identified the Leishmania homolog of receptors for activated C kinase (LACK) and phosphoenolpyruvate carboxykinase (PEPCK) proteins as strong immunogens and candidates for the development of a vaccine strategy. The present work focuses on the in silico prediction and characterization of antigenic epitopes that might interact with mice or human major histocompatibility complex class I. After immunogenicity prediction on the Immune Epitope Database (IEDB) and the Database of MHC Ligands and Peptide Motifs (SYFPEITHI), 26 peptides were selected for interaction assays with infected mouse lymphocytes by flow cytometry and ELISpot. This strategy identified nine antigenic peptides (pL1-H2, pPL3-H2, pL10-HLA, pP13-H2, pP14-H2, pP15-H2, pP16-H2, pP17-H2, pP18-H2, pP26-HLA), which are strong candidates for developing a peptide vaccine against leishmaniasis.

Finding epitopes of Klebsiella pneumoniae outer membrane protein-K17 (OMPK17) and introducing a 25-mer peptide of it as a vaccine candidate

No approved vaccine exists for Klebsiella pneumoniae yet. Outer membrane protein-K17 (OMPK17) is involved in K. pneumoniae pathogenesis. No information has been found about OMPK17 dominant epitopes in the literature. Therefore, this study aimed to predict both T cell and B cell epitopes of K. pneumoniae OMPK17 via immunoinformatics approaches. Both T cell (class-I and II) and B cell (linear and discontinuous) epitopes of OMPK17 were predicted. Several screening analyses were performed including clustering, immunogenicity, human similarity, toxicity, allergenicity, conservancy, docking, and structural/physicochemical suitability. The results showed that some regions of OMPK17 have more potential as epitopes. The most possible epitopes were found via several analyses including the selection of higher-scoring epitopes, the epitopes predicted with more tools, more immunogenic epitopes, the epitopes capable of producing interferon-gamma, the epitopes with more dissimilarity to human peptides, and non-toxic and non-allergenic epitopes. By comparing the best T cell and B cell epitopes, we reached a 25-mer peptide containing both T cell (class-I and class-II) and B cell (linear) epitopes and comprising appropriate physicochemical characteristics that are required for K. pneumoniae vaccine development. The in vitro/in vivo study of this peptide is recommended to clarify its actual efficiency and efficacy.

Supplementary information

The online version contains supplementary material available at 10.1007/s11756-023-01371-0.

Biomolecular interactions and binding dynamics of inhibitor arachidonic acid, with tyrosinase enzyme

Hyperpigmentation is a disorder caused by increased melanin deposition and changes in skin pigmentation. Inhibition of tyrosinase activity contributes to the control of food browning and skin pigmentation diseases. The effects of arachidonic acid (AA) on tyrosinase activity were examined using different spectroscopy methods including UV-VIS spectrophotometry, fluorescence spectroscopy, circular dichroism (CD) differential scanning calorimetry, and molecular dynamics (MD) simulations. Based on the kinetic results, arachidonic acid showed mixed-type of inhibition with Ki = 4.7 µM. Fluorescence and CD studies showed changes of secondary and tertiary structures of enzyme and a reduction of α-helix* amino acids after its incubation with different concentrations of AA, which is also confirmed by DSSP analysis. In addition, differential scanning calorimetry (DSC) studies showed a decrease in thermodynamic stability of enzyme from Tm = 338.65k for sole enzyme after incubation with AA in comparison with complex enzyme with Tm= 334.26k, ΔH =7.52 kJ/mol, and ΔS = 0.15 kJ/mol k. Based on the theoretical methods, it was found that the interaction between enzyme and AA follows an electrostatic manner with ΔG = -8.314 kJ/mol and ΔH = -12.9 kJ/mol. The MD results showed the lowest flexibility in the complex amino acids and minimal fluctuations in AA interaction with tyrosinase in Residue 240 to 260 and 66 to 80. Thus, AA inhibitory and structural and thermodynamic instability of tyrosinase supported advantages of this fatty acid for prevention of medical hyperpigmentation. Therefore, it is a good candidate for cosmetic applications. Communicated by Ramaswamy H. Sarma.

Affinity enhancement of CR3022 binding to RBD; in silico site directed mutagenesis using molecular dynamics simulation approaches

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a disease which caused by a novel beta coronavirus. Structural and non-structural proteins are expressed by the virus gene fragments. The RBD of the S1 protein of the virus has the ability to interact with potent antibodies including CR3022, which was characterized to target the S protein of the virus which can efficiently neutralize the SARS-CoV in vitro and in vivo. In current study, we aimed to design CR3022 based antibody with high affinity compared with wild-type CR3022 using MD simulation method. Two variants were designed based on the amino acid binding conformation and the free binding energy of the critical amino acids which involved in CR3022-RBD interactions were evaluated. In this study three complexes were evaluated; CR3022-RBD, V1-RBD and V2-RBD using molecular dynamics simulations carried out for 100 ns in each case. Then, all the complexes were simulated for 100 ns. In the next step, to calculate the free binding affinity of the wild CR3022 and mutant antibody (V1 and V2) with RBD, the PMF method was performed. The RMSD profile demonstrated that all three complexes were equilibrated after 85 ns. Furthermore, the free binding energy results indicated that the V2-RBD complex has the higher binding affinity than V1-RBD and CR3022-RBD complexes. It should be noted that in above variants, the electrostatic energy and the number of H-bonds between the antibody and RBD increased. Thus, it is suggested that both designed antibodies could be considered as appropriate candidates for covid-19 disease treatment.Communicated by Ramaswamy H. Sarma.

Identification of potential candidate vaccines against Mycobacterium ulcerans based on the major facilitator superfamily transporter protein

Buruli ulcer is a neglected tropical disease that is characterized by non-fatal lesion development. The causative agent is Mycobacterium ulcerans (M. ulcerans). There are no known vectors or transmission methods, preventing the development of control methods. There are effective diagnostic techniques and treatment routines; however, several socioeconomic factors may limit patients' abilities to receive these treatments. The Bacillus Calmette-Guérin vaccine developed against tuberculosis has shown limited efficacy, and no conventionally designed vaccines have passed clinical trials. This study aimed to generate a multi-epitope vaccine against M. ulcerans from the major facilitator superfamily transporter protein using an immunoinformatics approach. Twelve M. ulcerans genome assemblies were analyzed, resulting in the identification of 11 CD8+ and 7 CD4+ T-cell epitopes and 2 B-cell epitopes. These conserved epitopes were computationally predicted to be antigenic, immunogenic, non-allergenic, and non-toxic. The CD4+ T-cell epitopes were capable of inducing interferon-gamma and interleukin-4. They successfully bound to their respective human leukocyte antigens alleles in in silico docking studies. The expected global population coverage of the T-cell epitopes and their restricted human leukocyte antigens alleles was 99.90%. The population coverage of endemic regions ranged from 99.99% (Papua New Guinea) to 21.81% (Liberia). Two vaccine constructs were generated using the Toll-like receptors 2 and 4 agonists, LprG and RpfE, respectively. Both constructs were antigenic, non-allergenic, non-toxic, thermostable, basic, and hydrophilic. The DNA sequences of the vaccine constructs underwent optimization and were successfully in-silico cloned with the pET-28a(+) plasmid. The vaccine constructs were successfully docked to their respective toll-like receptors. Molecular dynamics simulations were carried out to analyze the binding interactions within the complex. The generated binding energies indicate the stability of both complexes. The constructs generated in this study display severable favorable properties, with construct one displaying a greater range of favorable properties. However, further analysis and laboratory validation are required.

Immunoinformatics-Aided Analysis of RSV Fusion and Attachment Glycoproteins to Design a Potent Multi-Epitope Vaccine

Respiratory syncytial virus (RSV) usually causes respiratory tract infections of upper airways in infants and young children. Despite recent medical advances, no approved vaccine is available to control RSV infections. Therefore, we conducted an immunoinformatics study to design and evaluate a potential multi-epitope vaccine against RSV. Sequence-based analyses of the glycoproteins F and G revealed a total of eight CD8 T-cell and three CD4 T-cell epitopes after considering antigenicity, binding affinity and other parameters. Molecular docking analysis confirmed that these T-cell epitopes developed strong structural associations with HLA allele(s). By integrating these prioritized epitopes with linkers and a cholera toxin-derived adjuvant, a multi-epitope vaccine was designed. The developed vaccine was found to be stable, non-allergenic, flexible and antigenic. Molecular docking analysis revealed a striking mean HADDOCK score (-143.3) of top-ranked vaccine-TLR cluster and a Gibbs free energy change (ΔG) value of -11.3 kcal mol-1. As per computational immune simulation results, the vaccine generated a high titer of antibodies (especially IgM) and effector T-cells. Also, codon optimization and in silico cloning ensured the increased expression of vaccine in Escherichia coli. Altogether, we anticipate that the multi-epitope vaccine reported in this study will stimulate humoral and cellular responses against RSV infection, subject to follow-up experimental validation.

Revisiting the Principles of Designing a Vaccine

Immune principles formulated by Jenner, Pasteur, and early immunologists served as fundamental propositions for vaccine discovery against many dreadful pathogens. However, decisive success in the form of an efficacious vaccine still eludes for diseases such as tuberculosis, leishmaniasis, and trypanosomiasis. Several antileishmanial vaccine trials have been undertaken in past decades incorporating live, attenuated, killed, or subunit vaccination, but the goal remains unmet. In light of the above facts, we have to reassess the principles of vaccination by dissecting factors associated with the hosts' immune response. This chapter discusses the pathogen-associated perturbations at various junctures during the generation of the immune response which inhibits antigenic processing, presentation, or remodels memory T cell repertoire. This can lead to ineffective priming or inappropriate activation of memory T cells during challenge infection. Thus, despite a protective primary response, vaccine failure can occur due to altered immune environments in the presence of pathogens.

Detecting the Dominant T and B Epitopes of Klebsiella pneumoniae Ferric Enterobactin Protein (FepA) and Introducing a Single Epitopic Peptide as Vaccine Candidate

Klebsiella pneumoniae causes various human infections. Ferric enterobactin protein (FepA) is a conserved protein of K. pneumoniae with high immunogenicity. In the present study, using comprehensive in silico approaches the T and B cell-specific epitopes of K. pneumoniae FepA were identified. The T (both class I and class II) and B (both linear and conformational) epitopes of FepA were predicted using prediction tools. The predicted epitopes were screened for human similarity, immunogenicity, antigenicity, allergenicity, toxicity, conservancy, structural and physicochemical suitability, and in case of T epitopes binding to HLA alleles, using numerous immune-informatics, homology modeling, and molecular docking approaches. These analyses led to introduce the most dominant FepA epitopes that are appropriate for vaccine development. Furthermore, we introduced an antigenic peptide containing both T and B epitopes which comprises suitable structural and physiochemical properties needed for vaccine development and it is conserved in many bacteria. Altogether, here the highly immunogenic T and B epitopes of FepA as well as a final epitopic peptide containing both T and B epitopes were found and introduced for future vaccine development studies. It is suggested that the actual efficiency and efficacy of our final epitopic peptide be investigated by in vitro/in vivo testing.

HLA-DRB1 Alleles Associated with Lower Leishmaniasis Susceptibility Share Common Amino Acid Polymorphisms and Epitope Binding Repertoires

Susceptibility for leishmaniasis is largely dependent on host genetic and immune factors. Despite the previously described association of human leukocyte antigen (HLA) gene cluster variants as genetic susceptibility factors for leishmaniasis, little is known regarding the mechanisms that underpin these associations. To better understand this underlying functionality, we first collected all known leishmaniasis-associated HLA variants in a thorough literature review. Next, we aligned and compared the protection- and risk-associated HLA-DRB1 allele sequences. This identified several amino acid polymorphisms that distinguish protection- from risk-associated HLA-DRB1 alleles. Subsequently, T cell epitope binding predictions were carried out across these alleles to map the impact of these polymorphisms on the epitope binding repertoires. For these predictions, we used epitopes derived from entire proteomes of multiple Leishmania species. Epitopes binding to protection-associated HLA-DRB1 alleles shared common binding core motifs, mapping to the identified HLA-DRB1 amino acid polymorphisms. These results strongly suggest that HLA polymorphism, resulting in differential antigen presentation, affects the association between HLA and leishmaniasis disease development. Finally, we established a valuable open-access resource of putative epitopes. A set of 14 HLA-unrestricted strong-binding epitopes, conserved across species, was prioritized for further epitope discovery in the search for novel subunit-based vaccines.

Prediction and identification of CD4+ T cell epitope for the protective antigens of Mycobacterium tuberculosis

CD4+T cell epitopes plays a key role in anti-tuberculosis (TB) immunity, CD4+T cell epitopes suitable for the domestic population are lacking. Therefore, we predicted and identified novel CD4+T cell epitopes.The bioinformatics software, namely, DNAStar (DNASTAR of the United States), SYFPEITHI (INTERFACTORS INSTITUT Für ZELL Biologie of Germany), RANKPEP, and NetMHC IIpan (National Cancer Institute, United States of America), were used to comprehensively predict the CD4+T cell immune epitope of Mycobacterium TB, and the predicted epitope polypeptide was synthesized by the standard Fmoc scheme. The proliferation of PBMC and CD4+T cells stimulated by peptides was preliminarily detected by the CCK8 method. Then, the candidate polypeptides screened out by the CCK8 method were verified again by the BrdU assay, and flow cytometry was performed to analyze further the extent of their stimulation on the proliferation of CD4+T cells. The changes in the secreted cytokines IFN-γ, TNF-α, IL-2, and IL-10 before and after the candidate polypeptide stimulation of CD4+T lymphocytes were detected by ELISA. The preliminary humoral immunity test was conducted by indirect ELISA to evaluate the serological diagnostic value of the CD4+T cell epitope polypeptide.In this study, 5 novel candidate CD4+T cell epitope polypeptides with the amino acid sequences of LQGQWRGAAGTAAQA, PVTLAETGSTLLYPL, AAAWGGSGSEAYQGV, QFVYAGAMSGLLDPS, and KAALTRTASNMNAAA and others that have not been reported in the research were predicted. For convenience, the 5 candidates were successively named as P39, P50, P40, P185, and P62. P39, P62, and the mixed peptide P39+P62 could effectively induce the proliferation of CD4+T cells and increase the secretion of IFN-γ, TNF-α, and IL-2 from the CD4+T cells, while reducing the content of IL-10. The serological test showed that the sensitivity, specificity, and area under the receiver operating characteristic curve (AUC) of P39 were 75%, 67.71%, and 0.844, respectively. The sensitivity, specificity, and AUC of P62 were 91.66%, 46.87%, and 0.649, respectively. The sensitivity of the mixed peptide P39+P62 was 95.83%, the specificity was 97.91%, and the AUC was 0.793.The P39 and P62 polypeptides were predicted and identified as potential CD4+T cell immune epitope polypeptides of M. TB. The polypeptide had better diagnosis effect, which provided potential candidate epitope polypeptides for the development of TB-specific diagnosis reagents and novel TB epitope vaccines.

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma

Immunoinformatics and molecular dynamics studies to predict T-cell-specific epitopes of four Klebsiella pneumoniae fimbriae antigens

Klebsiella pneumoniae (K. pneumoniae) is a causative agent of severe infections in humans. There is no publically available vaccine for K. pneumoniae infections yet. Here, using comprehensive immunoinformatics methods, T-cell-specific epitopes of four type 1 fimbriae antigens of K. pneumoniae were predicted and evaluated as potential vaccine candidates. Both CD8+ (class I) and CD4+ (class II) T-cell-specific epitopes were predicted and the epitopes similar to human proteome were excluded. Subsequently, the windows of class-II epitopes containing class-I epitopes were determined. The immunogenicity, IFN-γ production and population coverage were also estimated. Using the 3D structure of HLA and epitopes, molecular docking was carried out. Two best epitopes were selected for molecular dynamics studies. Our prediction and analyses resulted in the several dominant epitopes for each antigen. The docking results showed that all selected epitopes can bind to their restricted HLA molecules with high affinity. The molecular dynamics results indicated the stability of system with minimum possible deviation, suggesting the selected epitopes can be promising candidates for stably binding to HLA molecules. Altogether, our results suggest that the selected T-cell-specific epitopes of K. pneumoniae fimbriae antigens, particularly the two epitopes confirmed by molecular dynamics, can be applied for vaccine development. However, the in vitro and in vivo studies are required to authenticate the results of the present study.Communicated by Ramaswamy H. Sarma.

Effects of unsaturated fatty acids (Arachidonic/Oleic Acids) on stability and structural properties of Calprotectin using molecular docking and molecular dynamics simulation approach

Calprotectin is a heterodimeric protein complex with two subunits called S100A8/A9. The protein has an essential role in inflammation process and various human diseases. It has the ability to bind to unsaturated fatty acids including Arachidonic acid, Oleic acid and etc., which could be considered as a major carrier for fatty acids. In this study we aimed to appraise the thermodynamics and structural changes of Calprotectin in presence of Arachidonic acid/Oleic acid) using docking and molecular dynami simulation method. To create the best conformation of Calprotectin-Oleic acid/Arachidonic acid complexes, the docking process was performed. The complexes with the best binding energy were selected as the models for molecular dynamics simulation process. Furthermore, the structural and thermodynamics properties of the complexes were evaluated too. The Root Mean Square Deviation and Root Mean Square Fluctuation results showed that the binding of Arachidonic acid/Oleic acid to Calprotectin can cause the protein structural changes which was confirmed by Define Secondary Structure of Proteins results. Accordingly, the binding free energy results verified that binding of Oleic acid to Calprotectin leads to instability of S100A8/A9 subunits in the protein. Moreover, the electrostatic energy contribution of the complexes (Calprotectin-Oleic acid/Arachidonic acid) was remarkably higher than van der Waals energy. Thus, the outcome of this study confirm that Oleic acid has a stronger interaction with Calprotectin in comparison with Arachidonic acid. Our findings indicated that binding of unsaturated fatty acids to Calprotectin leads to structural changes of the S100A8/A9 subunits which could be beneficial to play a biological role in inflammation process.

In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach

Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.