Mutaib M. Mashraqi1
Ahmad Alzamami2
Norah A. Alturki3
Hassan H. Almasaudi1
Ibrar Ahmed4,5Saleh Alshamrani1*
Zarrin Basharat4*- 1Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Najran University, Najran, Saudi Arabia
- 2Clinical Laboratory Science Department, College of Applied Medical Science, Shaqra University, AlQuwayiyah, Saudi Arabia
- 3Clinical Laboratory Science Department, College of Applied Medical Science, King Saud University, Riyadh, Saudi Arabia
- 4Alpha Genomics Private Limited, Islamabad, Pakistan
- 5Group for Biometrology, Korea Research Institute of Standards and Science (KRISS), Daejeon, Republic of Korea
Kingella kingae is a Gram-negative bacterium that primarily causes pediatric infections such as septicemia, endocarditis, and osteoarticular infections. Its virulence is attributed to the outer membrane proteins having implications in bacterial adhesion, invasion, nutrition, and host tissue damage. TonB-dependent receptors (TBDRs) play an important role in nutrition and were previously implicated as vaccine targets in other bacteria. Therefore, we targeted the conserved β-barrel TBDR domain of these proteins for designing a vaccine construct that could elicit humoral and cellular immune responses. We used bioinformatic tools to mine TBDR-containing proteins from K. kingae ATCC 23330 and then predict B- and T-cell epitopes from their conserved β-barrel TDR domain. A chimeric vaccine construct was designed using three antigenic epitopes, covering >98% of the world population and capable of inciting humoral and adaptive immune responses. The final construct elicited a robust immune response. Docking and dynamics simulation showed good binding affinity of the vaccine construct to various receptors of the immune system. Additionally, the vaccine was predicted to be safe and non-allergenic, making it a promising candidate for further development. In conclusion, our study demonstrates the potential of immunoinformatics approaches in designing chimeric vaccines against K. kingae infections. The chimeric vaccine we designed can serve as a blueprint for future experimental studies to develop an effective vaccine against this pathogen, which can serve as a potential strategy to prevent K. kingae infections.
1 Introduction
Kingella kingae is a Gram-negative bacterium that was first identified in 1960 (Lebel et al., 2006). It was initially considered a commensal organism of the upper respiratory tract (Al-Qwbani et al., 2016). However, subsequent studies have shown that it can cause a range of infections, particularly in children under the age of 5 years (Williams et al., 2014). Currently, it is considered an emerging pathogen in the field of pediatrics (Dubnov-Raz et al., 2010). The bacterium’s virulence is attributed to its ability to adhere to host cells, invade host tissues, and damage host cells (Porsch, 2022). Infections caused by it are usually localized and can involve various organs, such as the joints, bones, and heart (Ceroni et al., 2013). The most common manifestation of its infection is osteoarticular infection, which can present as septic arthritis or osteomyelitis (Lundy and Kehl, 1998). The infection can also affect the heart valves and cause endocarditis, which is a serious and potentially life-threatening condition (Taubert and Dajani, 2001). Prompt diagnosis and appropriate treatment are crucial to prevent complications and ensure a good clinical outcome, while prevention through vaccine is a better strategy. However, to date, no vaccine exists against this bacterium.
TonB-dependent receptors (TBDRs) are outer membrane proteins found in many Gram-negative bacteria that play a key role in the uptake of essential nutrients, such as iron and vitamin B12, from the host environment (Blanvillain et al., 2007; Moynié et al., 2017). TBDRs are highly conserved among bacteria and are essential for bacterial survival and growth, making them attractive targets for the development of vaccines (Grassmann et al., 2017; Bettin et al., 2022). One advantage of targeting TBDRs is that they are expressed on the outer surface of the bacterial cell, making them accessible to antibodies and other immune system components. Several studies have investigated the potential of TBDRs as vaccine candidates, with promising results. For example, a vaccine based on the TBDR of the pathogenic bacterium Haemophilus influenzae was shown to be highly effective in preventing infection in a mouse model (Webb and Cripps, 1999). Other studies have focused on the development of vaccines targeting TBDRs in pathogenic bacteria such as Aeromonas hydrophila (Abdelhamed et al., 2017), Mycobacterium bovis (Bettin et al., 2022), Pseudomonas fluorescens (Hu et al., 2012), and Acinetobacter baumannii (Yang et al., 2017).
Traditional vaccine development approaches can be time-consuming and costly, often taking several years and requiring large amounts of resources (Mahoney and Maynard, 1999). Immunoinformatics approaches can accelerate the vaccine design process and reduce the cost by predicting epitopes that are likely to be effective (Kazi et al., 2018). The utilization of immunoinformatics allows researchers to analyze vast amounts of genomic and proteomic data to identify antigenic epitopes with desirable characteristics, such as high immunogenicity, conservation across strains, and binding affinity to immune receptors. By narrowing down the pool of potential epitopes, immunoinformatics streamlines the selection process, enabling more focused and efficient experimental validation. Chimeric vaccines, which are composed of multiple antigenic epitopes, have been widely used in infectious disease research and developed using immunoinformatics (Rahmani et al., 2019; Mishra et al., 2022). Using bioinformatics tools, researchers can predict B-cell and T-cell epitopes against the protein sequence of choice, to design a chimeric vaccine and predict the safety and immunogenicity in silico, reducing the cost and time required for experimental testing (Enayatkhani et al., 2021). Several studies have demonstrated the potential of immunoinformatics approaches in designing chimeric vaccines against infectious agents. For instance, Kota et al. (2021) identified B- and T-cell epitopes from the toxoid and hemolysin protein of Staphylococcus aureus to design a chimeric vaccine. The vaccine was predicted to be highly immunogenic and provided 83% direct and 50% passive immunization in a murine model. Another study identified MHC-binding epitopes from a Lom-like protein, a putative pilin subunit, and a section of the type III secretion structural protein EscC from the Escherichia coli O157:H7 (García-Angulo et al., 2014). The chimeric vaccine was designed and tested in mice, as well as through other experimental assays. The vaccine induced cytokine response and reduced colonization of the pathogen in a mouse model. Similarly, immunoinformatics-based predictions for the dengue virus vaccine led to good immune response mapping in a rabbit model (Kaushik et al., 2022).
Hence, designing a chimeric construct targeting the K. kingae TDR domain is assumed to be a promising approach for developing a vaccine against strains of this pathogen. Therefore, the objective of this study was to design a vaccine construct against K. kingae using a computational approach. Prediction of antigenic epitopes that are likely to be recognized by the immune system was followed by construct design. The best immunogenic construct was docked with immune receptors and immune response simulation. This bioinformatics-aided vaccine design approach can accelerate the vaccine design process to provide protection against K. kingae and reduce the cost of development.
2 Materials and methods
2.1 Target sequence retrieval
The whole proteome of K. kingae ATCC 23330 (GenBank accession: AFHS01000057.1; comprising 3,995 sequences) was scanned against the TonB-dependent receptor domain (Pfam ID: PF00593; length = 553 amino acids) using BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi; accessed 18 March 2023) (Altschul et al., 1997). Three protein sequences, comprising two homologs of TbpA (accession: EGK07547.1; WP_257003592.1) and an hpuB (accession: WP_003785727.1), were obtained as significant hits. Both consisted of a TonB-dependent receptor plug domain and a TonB-dependent receptor-like β-barrel domain (Figure 1). The sequence of the larger β-barrel-like domain (Pfam ID: PF00593; Figure 1) was taken and used for creating a vaccine construct against K. kingae ATCC 23330.
FIGURE 1. 3D structures of (A) TbpA (https://alphafold.ebi.ac.uk/entry/F5S962; retrieved 18 March 2023) and (B) hpuB (https://alphafold.ebi.ac.uk/entry/F5S5I2; retrieved 18 March 2023) of K. kingae, with the larger TBDR β-barrel domain shown in bluish purple and the smaller TonB plug domain shown in teal. The rest of the sequence is shown in yellowish green.
2.2 Epitope mapping and construct preparation
Immunoinformatics analysis was carried out using several public servers. Antigenicity of the sequence was predicted using the VaxiJen server (http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html; retrieved 18 March 2023) (Doytchinova and Flower, 2007). The B-cell epitopes were predicted using the two servers ABCpred (https://webs.iiitd.edu.in/raghava/abcpred; retrieved 19 March 2023) (Saha and Raghava, 2006) and BcePred (https://webs.iiitd.edu.in/raghava/bcepred/; retrieved 19 March 2023) (Saha and Raghava, 2004). ABCpred uses an artificial neural network to predict continuous B-cell epitopes based on a dataset of 700 experimental epitopes. We took a threshold of 0.5 and a window of 16 for prediction. BcePred uses data on 1,029 B-cell epitopes and takes into account parameters like hydrophobicity and flexibility. We used a default threshold ranging from −3 to 3 for these parameters. A consensus of results from both servers was taken. T-cell epitopes were predicted using IEDB-AR (https://www.iedb.org/; retrieved 20 March 2023) (Dhanda et al., 2019) and were filtered based on their antigenicity, conservancy, and accessibility to the immune system. Overlap with B-cell epitopes was studied, and the selected ones that fulfilled all criteria were arranged in a linear order to form a chimeric vaccine.
Accurately assessing and maximizing world population coverage is essential in designing vaccines that have broad efficacy and can effectively combat infectious diseases on a global scale. It helps identify the extent to which a vaccine is effective in protecting a specific pathogen across diverse populations worldwide, using HLA allele frequencies and MHC binding data. World population coverage of these epitopes was mapped using the IEDB server (http://tools.iedb.org/population/; retrieved 21 March 2023). This server maps coverage for 115 countries, covering 21 ethnicities.
A flexible linker sequence was introduced between the epitopes to ensure that each epitope maintained its conformation and function. In total, nine constructs were created and checked for antigenicity using VaxiJen and allergenicity using AllerTOP (http://www.ddg-pharmfac.net/AllerTOP/; retrieved 22 March 2023) (Dimitrov et al., 2014) and AllerCatPro (https://allercatpro.bii.a-star.edu.sg/; retrieved 22 March 2023) (Maurer-Stroh et al., 2019). Construct no. 4 (referred to as C4 hereafter), having the best antigenicity score and a non-allergen, was taken for structure modeling.
2.3 Structure modeling and dynamics simulation
The 3D structure was modeled using AlphaFold v2.3.1 (Jumper et al., 2021). AlphaFold is a deep learning algorithm developed by DeepMind that predicts the 3D structure of a protein from its amino acid sequence. The final predicted structure was selected based on the conformation having the lowest energy. A per-residue estimate of its confidence on a scale from 0 to 100 was measured using a scale called pLDDT, with a value >90 referring to high accuracy, and 70 or above expected to be modeled well (Zifruddin et al., 2023). The structure was minimized using the OPLS3e force field in the software package Maestro (van Zundert et al., 2021). Bond orders were assigned using the CCD database, and Het states were generated using Epik at pH 7.0+/-2.0. Hydrogen bond assignment was carried out using PROPKA at pH 7. Heavy metals were converged to 0.3 Å. This minimization process enhances quality by removing any small stereochemical violations present and generally results in minimal differences between the predicted and relaxed structures. Molecular dynamics simulation was also conducted for 100 ns in Desmond, Schrodinger Inc. (United States). The construct was prepared by adding hydrogens, assigning bond orders, and adding disulfide bonds. Refinement was carried out at pH 7 in the PROPKA module. Restrained minimization was carried out by converging heavy atoms to 0.30 Å. The solvation model was TIP3P, box=orthorhombic, box size calculation method: Buffer and neutralization was carried out by adding counter ions (Na+ or Cl-according to requirement). Simulation parameters were force field: OPLS2005 force field; Ensemble: NVT; system temperature: 300 K; pressure: 1.01 (bar); recording interval: 100 ps. The rest of the options were default. RMSD and RMSF were obtained along with the secondary structure over a time course of 100 ns.
2.4 Molecular docking
The binding affinity of the chimeric vaccine to the immune receptors (HLA alleles and TLRs), specifically known to interact with the adjuvant used in C4, was carried out using ClusPro (Kozakov et al., 2017). TLRs are involved in recognizing pathogens, while HLAs help present antigens to T cells. Among TLRs, TLR-1, TLR-1/2, and TLR-4 were included, while HLA receptors included HLA-A*0201, HLA-B*5301, HLA-CW3, HLA-DRA1, HLA-DRB1, HLA-DP1, HLA-DP2, HLA-DQA1, and HLA-DQB1. ClusPro is based on a fast Fourier transform (FFT) algorithm that uses rigid-body docking to generate a large number of candidate models for protein–protein complexes. It then uses a clustering algorithm to group these candidate models into clusters based on their similarity, and it ranks the top models by their calculated binding energy (Comeau et al., 2004). Interaction prediction could help infer binding between C4 and human immune receptors, and the subsequent initiation of immune response. Post-docking interactions were visualized using Protein Interactions Calculator (PIC) (http://pic.mbu.iisc.ernet.in/; accessed 22 March 2023) (Tina et al., 2007).
Normal mode analysis is a computational method used to study the collective vibrational motions of a protein or other biomolecules (Case, 1994). It can provide insights into protein dynamics and conformational changes. WEBnm@ (http://apps.cbu.uib.no/webnma3/; accessed 23 March 2023) (Tiwari et al., 2014) was used to study conformational changes and for comparing unbound and bound C4 via the lowest-energy collective vibrational modes, such as global twisting or bending motions. Atomic displacement graphs in several modes and a correlation matrix were used to observe the changes in bound and unbound C4. MD simulation was carried out for all the receptors bound with C4. Parameters were the same as described in Section 2.3, except for the time duration being reduced to 50 ns due to the complexity of the molecular system.
2.5 Immune response simulation
C-IMMSIM (Castiglione and Bernaschi, 2004) was used to simulate the interactions between the vaccine construct, immune cells, antigens, and cytokines, and to study the dynamics of the immune response under different conditions. It is a C++-based agent-based modeling framework that integrates both cellular and molecular components of the immune system. Hence, the response of various immune cells, such as T cells, B cells, and dendritic cells, to C4 and K. kingae TbpA and hpuB challenge was attempted. For most of the vaccines currently in use, 4 weeks is the minimum recommended time between the first and second doses (Saha et al., 2021). We utilized the same values. The entire simulation ran for 1,000 time steps, which are ∼11 months (1 timestep=8 h). Two vaccine injections were given on days 1 and 30 (i.e., after 1 month), followed by the bacterial challenge on day 240 (i.e., after 8 months since vaccine initiation and 7 months after the end of vaccination) to check the efficacy of the vaccination process. Host HLA selection was A0101, A0102, B0702, B0704, DRB101, and DRB102. The rest of the parameters were as described previously (Rahman et al., 2020).
2.6 Cloning
In silico cloning is a computational method used to predict the optimal restriction enzyme sites and sequence overlap for the insertion of a DNA sequence of interest into a cloning vector (Celie et al., 2016). It is carried out to design the cloning strategy and is helpful in guiding wet laboratory experiments. JCat (Grote et al., 2005) was used for optimization of the codon usage of C4 to enhance its expression in E. coli. Codon optimization was attempted because it can lead to higher yields of the antigen protein and a stronger immune response in the host (Ramakrishna et al., 2004). The C4 protein sequence was first reverse-translated to the DNA sequence, and then, the codon usage was adjusted to match the preferred codon usage of E. coli. A suitable pET vector was then selected for cloning this sequence in SnapGene software (available at https://www.snapgene.com/) based on optimal restriction enzyme sites for C4. It was then cloned, and the image was saved in .jpeg format.
3 Results
3.1 Construct assembly
In total, 538 16-mer B-cell epitopes were obtained from the ABCpred server, and 11 epitopes were obtained from the BcePred server. Additionally, 29,403 epitopes were predicted for MHC-I and 14,553 for MHC-II. Finalization was carried out based on various properties like immunogenicity and sequence coverage. The final MHC-I and MHC-II epitopes were checked for cumulative population coverage. This was carried out to identify epitopes that are likely to be recognized by a large proportion of the population to design a vaccine effective in a broad range of individuals. By targeting epitopes that are highly conserved and commonly presented by MHC molecules, we can increase the likelihood of generating a robust and broadly protective immune response. According to the IEDB server, the cumulative MHC-I and MHC-II epitope population coverage of the world was predicted to be 98.55% (Supplementary Figure S1). This means that more than 98% of the world’s population is likely to have at least one MHC-I or MHC-II allele that can bind to this set of epitopes. Achieving high world population coverage is a desirable goal in vaccine design as it aims to ensure that a significant majority of individuals, regardless of geographic location, genetic variations, or other demographic factors, can mount an effective immune response against the targeted pathogen. A vaccine with high world population coverage shows potential to reduce the overall burden of the disease, prevent transmission, and contribute to global public health.
Finally, three overlapping validated epitopes (Table 1) were selected for designing the vaccine constructs (Supplementary Table S1). The overlap with B-cell epitopes was then studied to select those that fulfilled all criteria and arranged in a linear order to form a chimeric vaccine. To ensure the stability and functionality of the epitopes, a flexible linker sequence was introduced between them, along with adjuvants like β-defensin and flagellin. Adjuvants are added to enhance the immunogenicity of the vaccine and stimulate the immune response (Bastola et al., 2017). Nine different vaccine constructs were prepared and checked for antigenicity and allergenicity using VaxiJen, AllerTOP, and AllerCatPro servers. C4, which had the best antigenicity score and was non-allergenic, was selected for structure modeling.
TABLE 1. Final epitopes used for preparing the vaccine construct.
3.2 Structure modeling and simulation
The 3D structure of a vaccine construct is important because it can provide insights into its function, interactions with other molecules, and potential immunogenicity. The 3D structure (Figure 2A) comprised 2 sheets, 5 β-hairpins, 3 β-bulges, 8 strands (Figure 2B), 3 helices, 14 β-turns, 15 gamma turns, and 4 disulfides. Disulfides are covalent bonds between two cysteine residues that help stabilize the protein structure and prevent degradation (Trivedi et al., 2009).
FIGURE 2. (A) 3D model of C4 colored according to model confidence. (B) 2D depiction of C4 topology. (C) LDDT plot of the predicted structure. (D). Ramachandran plot of the predicted structure.
The first 50 residues modeled had a high confidence level (Figure 2C), while the overall Ramachandran plot (Figure 2D) showed 74% residues in highly favored regions, 24% in allowed regions, and the rest in disallowed regions.
The construct stabilized after 30 ns of simulation, and its RMSD remained approximately 15 Å after 60 ns (Figure 3A). RMSF remained less than 12 Å for most of the residues, and there was very little variation between heavy atoms, side chains, and the backbone (Figure 3B). The percentage of helices was 8.89 and that of strands was 25.58%, making SSE>30% of the residues over the course of simulation (Figure 3C). A major portion of the structure retained their composition of helices and sheets throughout the simulation (Figure 3D).
FIGURE 3. (A) RMSD plot of construct C4. (B) RMSF plot of C4. (C) Secondary structural elements (SSEs) like α-helices and β-strands throughout the simulation. (D) Residue scale SSE (α-helices shown in orange; β-strands shown in sky blue) assignment over 100 ns.
3.3 Docking and dynamics simulation analysis
ClusPro is a widely used docking program that predicts the binding affinity and orientation of two protein structures (Kozakov et al., 2017). Here, it was used to predict the interaction between the vaccine construct and the target proteins TLR1, TLR2, and TLR4, and various HLA molecules. The results showed that C4 has a high predicted binding affinity with all the target proteins, and docking scores range from −1550.0 to −1017.4, indicating strong interactions (Table 2). The lowest docking score was obtained with the protein TLR2, followed by TLR1 and TLR4, which suggests that C4 may have a high binding affinity for these receptors. The HLA molecules tested also showed a good interaction with C4, with docking scores ranging from −1709.4 to −1162.9, and binding is likely to potentially increase C4 efficacy as a vaccine. Previously, the interaction of TLR4 with the vaccine construct designed against Haemophilus influenzae (ClusPro score: −1035.2 kcal/mol) and Streptococcus pneumoniae (ClusPro score: −1150.5 kcal/mol) was reported to be robust (Mazumder et al., 2023). Our score was higher than the score in these studies, indicating a highly stable binding affinity (AlChalabi et al., 2022). Similarly, the binding interaction score of the vaccine construct designed against Candida tropicalis and HLA-A*0201 has been reported as − 1178.4 kcal/mol (Akhtar et al., 2022), whereas it was −1218.9 kcal/mol in this study, indicating good binding affinity. It is important to ensure the reliability and accuracy of computational methods in determining binding affinities. For this purpose, we incorporated cross-validation of docking calculations by introducing a control (TLR1–TLR2 interaction; PDB ID: 2Z7X) with a strong experimentally known binding affinity (Jin et al., 2007) and compared against the docking prediction for our construct. The TLR1–TLR2 interaction cluster with the largest number of members had the lowest energy of −936.5 kcal/mol. Energy values of C4 binding with immune receptors were much lower than this value, indicating more stable and possibly robust interactions. The stability of the interactions between C4 and immune receptors is an important factor in determining the effectiveness of the immune response. Lower energy values imply stronger binding between C4 and the receptors, suggesting a higher likelihood of successful recognition and activation of the immune system. Robust interactions can trigger a more vigorous immune response, which is desirable in combating pathogens.
TABLE 2. HLA and TLR interaction statistics with the designed vaccine construct. 3D depiction of interactions is shown in Supplementary Figure S3.
Delving into the details of interactions, C4 formed a significant number of hydrophobic interactions and hydrogen bonds with all interactors, which are crucial for stabilizing the interaction. A smaller number of ionic bonds and aromatic interactions were also observed, which contributed to the stability of the interaction. The highest number of hydrophobic interactions were observed for HLA-DRB1 (n= 32), the largest number of hydrogen bonds was observed for HLA-DP2 (n = 23), and maximum ionic interactions were detected for HLA-DRB1 (n = 3). C4 also formed aromatic interactions with some interactors, with the highest number observed for HLA-DRB1 (n = 10). Aromatic interactions are relatively weaker than hydrophobic interactions and hydrogen bonds, but they can add to the strength of the interaction. Normal mode analysis for C4 showed major fluctuations around residues 30 and 40, and between residues 160 and 170 (Supplementary Figure S2), while binding appeared restrained for the TLR4 receptor. Displacements were also observed for HLA-DRA1 and HLA-DQB1.
Comparative MD simulation (Figure 4) of the C4 vaccine construct and bound immune receptors revealed that the RMSD of the complexes (Supplementary Figure S4) was much lower than that of the vaccine construct alone. This suggests that binding to the immune receptors stabilizes the structure of the vaccine construct. RMSD values varied among different immune receptors. HLA-DP1, HLA-DQA1, and TLR-4 had comparatively low average RMSD values (less than 6), indicating that the complexed vaccine remained relatively stable when bound to these receptors, while TLR1/2, HLA-B*5301, HLA-DRA1, HLA-DRB1, HLA-DP2, and HLA-CW3 complexed with the vaccine had higher average RMSD values (less than 8), suggesting that the complex with these receptors experienced more structural fluctuations. HLA-A*0201 and HLA-DQB1 exhibited a more dynamic behavior than the rest of the immune receptors. RMSD of HLA-A*0201-C4 reached up to 15 Å at approximately 10 ns and then stabilized for some time for less than 8 Å until 40 ns. It again increased to more than 8, and the trend continued upward. HLA-DQB1 reached up to 12 Å after 20 ns and then stabilized but remained higher than 10.5 Å, peaking toward 12 Å after 35 ns. The initial increase in RMSD indicates structural adjustments or fluctuations as they accommodate the vaccine construct. The subsequent stabilization suggests a period of relative structural stability, but this stability is at a higher RMSD level, indicating that these receptors undergo some level of structural change even in their stable states and as time progresses. The results are in alliance with the findings of normal mode analysis, where HLA-DQB1 showed higher displacement.
FIGURE 4. RMSD plot of 50-ns simulation showing C4=vaccine construct; C4-6NIH=vaccine–TLR1 complex; C4-2Z80=vaccine–TLR1/2 complex; C4-3FXI=vaccine–TLR4 complex; C4-1AKJ=vaccine–HLA-A*0201 complex; C4-1A1M=vaccine–HLA-B*5301 complex; C4-1EFX=vaccine–HLA-CW3 complex; C4-1A6AA=vaccine–HLA-DRA1 complex; C4-1A6AB=vaccine–HLA-DRB1 complex; C4-3LZQA=vaccine–HLA-DP1 complex; C4-3LZQB=vaccine–HLA-DP2 complex; C4-1JK8A=vaccine–HLA-DQA1 complex; and C4-1JK8B=vaccine–HLA-DQB1 complex.
Similarly, it was noted that the RMSF varied for complexes (Figure 5) compared to the vaccine construct alone (Figure 3B). Findings provide insights into the dynamic nature of the interaction between these immune receptors and the vaccine construct, suggesting that the receptors undergo structural adaptations during the simulation. The initial residue portion of the graph shows the immune receptor RMSF, while the last 200 residues depict the vaccine construct C4. The RMSF of TLR1 was less than 3 Å while that of the construct C4 was less than 5 Å on average. It exceeded 4.5 Å on average in the unbound form. It was even lower for TLR1/2, TLR4, HLA-DRA1, HLA-DRB1, and HLA-DP1. HLA-B*5301, HLA-CW3, and HLA-DQB1 showed higher fluctuations. These observations highlight the dynamic nature of the interaction between the immune receptors and the vaccine construct. Different receptors exhibit varying levels of stability and flexibility during the simulation, which suggests that the dynamic nature of receptor–vaccine interactions observed in the simulation may have implications for the effectiveness of a vaccine. A stable interaction with immune receptors could be important for eliciting a strong and targeted immune response, which was depicted by the receptors TLR1/2, TLR4, HLA-DRA1, HLA-DRB1, and HLA-DP1.
FIGURE 5. RMSF plots of the vaccine construct bound with (A) TLR1. (B) TLR1/2. (C) TLR4. (D) HLA-A*0201. (E) HLA-B*5301. (F) HLA-CW3. (G) HLA-DRA1. (H) HLA-DRB1. (I) HLA-DP1. (J) HLA-DP2. (K) HLA-DQA1. (L) HLA-DQB1.
3.4 Immune response simulation
The C-IMMSIM results describe the efficacy of the vaccination process against K. kingae infection, attempted by injecting tbpA and hpuB proteins after 7 months of vaccination. Memory B cells kept decreasing after 3 months, while non-memory B cells kept increasing after 5 months (Figure 6A). IgM (first antibody to be produced in response to an initial exposure to a pathogen) and the total B-cell population plateaued (at ∼500 cells per mm3) 4 months after vaccination. IgM is the first antibody that is produced in response to an initial exposure to a pathogen, and the plateauing of its levels and that of the total B-cell population suggests that the immune system mounted a robust initial response to the vaccine, but this response might not continue to increase or remain at a high level over an extended period.
FIGURE 6. Immune cell population post-C4 injection. (A) Graph depicting the change in the B-cell population. (B) Varied antibody titer. (C) Altered TH cell population. (D) TC cell population for active, resting, duplicating, and anergic states.
The immunoglobulin population stabilized after 4 months, with the antibody titer count of IgG1+ IgG2 (antibodies that are involved in long-term immunity) and IgM+IgG being approximately 20,000 for 100,000 antigens per ml (Figure 6B), indicating that the immune system has produced a significant amount of antibodies against the pathogen, which should provide strong protection. T-helper (TH) cells reached up to 8,000 cells per mm3 after 20 days of the first injection (Figure 6C), with a doubling of TH1-active cells. Resting and active TH cell population stabilized at 1,000 cells per mm3 after 300 days. TH cells have a relatively short lifespan and are continuously produced and replaced by the immune system. Therefore, the stabilization of the TH cell population may not necessarily imply a lack of long-term protection. Immune response to the vaccine was not uniform across all T-cell populations. The TC memory cell population remained stable, while non-memory cells fluctuated, with active cells decreasing and resting cells increasing after 4 months of vaccination (Figure 6D).
3.5 Cloning of C4
In silico cloning allows researchers to test different cloning strategies and find the best combination of restriction enzymes and vector sequences to use without the need to perform multiple trial-and-error experiments in the laboratory. This can save time and resources in the cloning process. C4 was cloned in the pET-30a (+) vector in E. coli K12 (Figure 7).
FIGURE 7. Cloned C4 in the pET-30a+ vector, shown in red.
The pET-30a (+) vector is commonly used in molecular biology as it contains several useful features for protein expression in E. coli, such as strong promoters, a His-tag for protein purification, and antibiotic resistance markers for selection. The choice of HindIII and BstEII enzymes was based on their ability to generate compatible overhangs or “sticky ends” that can be easily ligated to the corresponding ends of the pET-30a (+) vector. C4 was then inserted in the vector and cloned in E. coli. This approach of in silico cloning could help save time and money by allowing the design and test of virtual constructs before actually synthesizing or cloning them in the laboratory.
4 Discussion
Proteins known to be involved in iron acquisition in bacteria are potential targets for vaccine development (Schryvers and Stojiljkovic, 1999; Glanfield et al., 2007; Ismail et al., 2022). We aimed to utilize the conserved domain of TBDRs involved in iron and nutrient acquisition for vaccine design against the K. kingae reference strain ATCC 23330. For this purpose, we scanned the proteome of K. kingae and identified two homologs of the TBDR protein (TbpA and hpuB) comprising the TBDR β-barrel domain. Previous studies identified TbpA and hpuB proteins as important virulence factors in other bacteria, such as Actinobacillus pleuropneumoniae (Klitgaard et al., 2010) and Neisseria gonorrhoeae (Lin et al., 2022). Finney et al. showed that antibodies against these proteins can protect against Neisseria lactamica meningococcal infection in rabbits (Finney et al., 2008). The TBDR transporter of Neisseria gonorrhoeae, i.e., TbpB, has also been targeted for vaccine design, in addition to TbpA, and shown to invoke an immune response in mice (Greenawalt, 2021). Additionally, several vaccine candidates targeting TBDRs have been developed for various bacterial pathogens, including Acinetobacter baumannii, Pseudomonas aeruginosa, and Klebsiella pneumonia (Wang et al., 2021). Therefore, we used the beta-barrel domain of TBDR from the K. kingae reference strain ATCC 23330 for creating a vaccine construct by predicting B- and T-cell epitopes. The selected ones fulfilling all criteria were manually arranged with adjuvants, linkers, etc., to assemble a vaccine construct. The generation and use of B- and T-cell epitopes through computational servers is an established strategy for fabricating a vaccine construct capable of inducing protective immune responses against a pathogen (Wang et al., 2020). B-cell epitopes are regions on the surface of a pathogen that can be recognized by antibodies (Ponomarenko and Van Regenmortel, 2009), while T-cell epitopes are fragments of a pathogen’s proteins that can be recognized by T cells (Khatoon et al., 2017). By including both types of epitopes in a vaccine construct, the immune system is activated to generate a robust and specific immune response. Population coverage analysis (>98%) of the prioritized epitopes suggests that the antigens are likely to be effective against a majority of individuals within the world population.
Adjuvants enhance the immune response (Bastola et al., 2017), and the best construct in this study (C4), with the highest antigenicity, had a β-defensin adjuvant. β-defensin is an antimicrobial peptide that can stimulate the innate immune system (dendritic cells and macrophages) and enhance antigen-specific immune responses (Shelley et al., 2020). β-defensin may be induced by inflammatory mediators in the body and has broad-spectrum antimicrobial activity against bacteria, viruses, and fungi (Gao et al., 2021). After assembling the construct, it was subjected to 3D structure modeling to analyze its capability to bind the immune system components like HLA alleles and TLRs, and invoke an immune response. Recent advances in deep learning techniques have led to the development of AlphaFold (Jumper et al., 2021), an AI-based tool for predicting protein structures with remarkable accuracy. We attempted 3D structure modeling of C4 using AlphaFold, and it was noted that the structure consisted of four natural disulfide bonds. Disulfide bonds are formed between two cysteine residues, where the thiol groups (-SH) on each cysteine react to form a covalent bond (-S-S-) under oxidizing conditions (Saito et al., 2003; Barford, 2004). These are crucial for maintaining the stability and conformation of the protein in the context of vaccine design and can be engineered into the construct to enhance its stability and immunogenicity (Scheiblhofer et al., 2017). However, in C4, these were present naturally, depicting a stable construct. C4 was then docked with critical components of the immune system using ClusPro (Comeau et al., 2004). ClusPro is a widely used tool for predicting protein–protein interactions, and it has been applied in several studies to predict the interaction between vaccine constructs and immune system components. One such example is a study between TLR9, HLA class I/II alleles, and a vaccine construct against Kaposi’s sarcoma (Chauhan et al., 2019), while another study on Candida auris vaccine design reported docking between TLR5 and MHC class-II HLA DRB_0101 with scores in a similar range (Akhtar et al., 2021). The docked results indicated that C4 has a high binding affinity with all the target proteins, with hydrophobic interactions and hydrogen bonds stabilizing the interactions. However, it is important to consider that binding affinity alone does not guarantee an effective immune response or protection against a specific pathogen. Other factors, such as the overall immunogenicity of the antigen, the activation of appropriate immune pathways, and the generation of memory responses, are also critical in achieving a robust and protective immune reaction.
Normal mode analysis was then performed to investigate the flexibility and dynamics of C4 with immune components. Low-frequency vibrations, representing the collective movements of the atoms in the protein, were calculated, and the results showed that C4 had major fluctuations around residues Pro30 and Thr40, and between the residues at position 160–170 (DCSGRVIKGSG). These residues were in the loop regions and, therefore, more flexible than the rest of the structure (Supplementary Figure S5). Loop regions are segments of the protein structure that connect secondary structural elements such as ɑ-helices and β-sheets (Chavan, 2014). These regions are often more flexible and exhibit more conformational changes than the other parts of the protein (Papaleo et al., 2016). Loop regions have immunogenic characteristics and are used in vaccine design for HIV (Montero et al., 2008), Leptotricha buccalis (Alshammari et al., 2022), and dengue virus (Urakami et al., 2017). Previously, Ismail et al. described fluctuations in the loop regions of the vaccine construct against Enterobacteriaceae, with loops inferred to adopt a stable conformation with increasing time, during the simulation (Ismail et al., 2020). They suggested that this might help in flexible binding with host immune components. Previously, the CDR H3 loop-dominated recognition of CD4 epitopes was also reported for HIV (Pancera et al., 2017). Thus, it is inferred that these regions may become displaced and help in flexible binding with HLAs and TLRs of the immune system. The varying levels of structural stability and flexibility observed through the RMSD and RMSF plot of different receptors via MD simulation also have implications for understanding immune responses and guiding the development of vaccines with improved efficacy. The stability and specificity of these interactions may be optimized to enhance the C4-mediated immune response.
Immune response simulation analysis revealed that the vaccine construct induced a robust initial immune response against K. kingae infection through the introduction of tbpA and hpuB proteins in the body. The production of immunoglobulin antibodies and their subsequent stabilization, along with the doubling of TH1-active cells, indicates strong protection against the pathogen. However, the immune response was not uniform across all T-cell populations. Although the long-term protection provided by the TC memory cells remains stable, the short-term immune response decreased after a few months of vaccination. Fluctuations in some cells, e.g., the effector T cells, may indicate a natural reaction of the immune system to the vaccine and not a lack of efficacy. The immune system is complex and dynamic, and its response to a vaccine can vary among individuals and over time (Shen-Orr and Furman, 2013; Laserson et al., 2014; Brodin and Davis, 2017). Therefore, it is important to consider multiple parameters of the immune response when assessing the efficacy of a vaccine. In silico immune simulations can provide valuable insights into immune responses, but they cannot fully capture the complexity of the immune response (Russo et al., 2020) and, therefore, must still be validated through in vitro and in vivo experiments. Hence, further studies are needed to fully understand the dynamics of the immune response over an extended period of time for this vaccine construct against K. kingae.
There were several limitations to our work. In silico protein structure prediction is challenging, and inaccuracies in the 3D structure of the vaccine construct might have affected the reliability of docked interaction findings. Mutations in the proteins used for designing vaccines over time may alter their antigenic properties, hence compromising the reliability of our findings over time. Apart from this, biological systems are highly complex, and many factors contribute to the immune response. Vaccine efficacy depends on multiple factors, including innate and adaptive immune responses, cellular and humoral immunity, and long-term memory. Integrating all these factors accurately is difficult, and computational models may oversimplify certain reactions, leading to inaccuracies in predicting actual immune responses. Our understanding of how the immune system recognizes and responds to the computationally predicted epitopes based on molecular docking analysis might not always be accurate, and information about the induction of the desired immune response is uncertain. Despite these limitations, computational vaccine design remains a valuable tool that can significantly accelerate the vaccine development process and guide experimental research.
5 Conclusion
The K. kingae species hosts TBDRs, which acquire and transport important substances like iron and vitamins. These can be used for the design of a chimeric vaccine construct against the conserved region (existing as a motif or domain). Chimeric vaccines can reduce the risk of adverse reactions as they typically contain only selected antigenic epitopes rather than whole pathogens or toxins. Therefore, a multi-epitope vaccine construct was designed and tested for immune response elicitation through docking with immune system components and simulation. The use of computational modeling and docking studies can provide valuable insights into the interactions between vaccine components and immune system components, which can inform the design and development of more effective vaccines. RMSD and RMSF plots after MD simulation revealed that different receptors exhibited varying levels of stability and flexibility during the simulation, which depicted that C4 adapts and interacts differently with different receptors. Immune response simulation through agent-based modeling showed good binding affinity along with a robust immune response elicitation for C4 against K. kingae infection. However, the decrease in memory B cells may indicate that the immunity provided by the vaccine could wane over time, potentially leading to a higher risk of reinfection and the need for booster doses. Additional experimental validation through in vitro and in vivo studies is suggested to confirm the efficacy of the vaccine construct in inducing an immune response and improving the longevity of immunity.
Data availability statement
The proteome data used in this study is present in the NCBI database (GenBank accession number: AFHS01000057.1). Protein sequences for vaccine design exist in NCBI as TbpA (Accession: EGK07547.1; WP_257003592.1) and an hpuB (Accession: WP_003785727.1).
Author contributions
MM: conceptualization, data curation, formal analysis, funding acquisition, investigation, and writing–original draft. AA: data curation, methodology, and writing–review and editing. NA: data curation, formal analysis, visualization, and writing–review and editing. HA: data curation, funding acquisition, and writing–review and editing. IA: resources, and formal analysis. SA: conceptualization, funding acquisition, methodology, project administration, validation, visualization, and writing–review and editing. ZB: formal analysis, funding acquisition, supervision, visualization, and writing–original draft.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Deanship of Scientific Research at the Najran University Research Support Program under the grant program code NU/DRP/MRC/12/42. IA was funded by the National Research Foundation (NRF) of Korea, Grant No. RS-2023-00223245.
Acknowledgments
The authors acknowledge the Deanship of Scientific Research at Najran University for funding this work, under the General Research Funding program grant code NU/DRP/MRC/12/42.
Conflict of interest
Authors IA and ZB were employed by Alpha Genomics Private Limited.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1258834/full#supplementary-material
References
Abdelhamed, H., Ibrahim, I., Nho, S. W., Banes, M. M., Wills, R. W., Karsi, A., et al. (2017). Evaluation of three recombinant outer membrane proteins, OmpA1, Tdr, and TbpA, as potential vaccine antigens against virulent Aeromonas hydrophila infection in channel catfish (Ictalurus punctatus). Fish. Shellfish Immunol. 66, 480–486. doi:10.1016/j.fsi.2017.05.043
Akhtar, N., Joshi, A., Kaushik, V., Kumar, M., and Mannan, M. A. U. (2021). In-silico design of a multivalent epitope-based vaccine against Candida auris. Microb. Pathog. 155, 104879. doi:10.1016/j.micpath.2021.104879
Akhtar, N., Singh, A., Upadhyay, A. K., and Mannan, M. (2022). Design of a multi-epitope vaccine against the pathogenic fungi Candida tropicalis using an in silico approach. J. Genet. Eng. Biotechnol. 20, 140–215. doi:10.1186/s43141-022-00415-3
Al-Qwbani, M., Jiang, N., and Yu, B. (2016). Kingella kingae–associated pediatric osteoarticular infections: an overview of 566 reported cases. Clin. Pediatr. 55, 1328–1337. doi:10.1177/0009922816629620
Alchalabi, R., Al-Rahim, A., Omer, D., and Suleiman, A. A. (2022). Immunoinformatics design of multi-epitope peptide-based vaccine against Haemophilus influenzae strain using cell division protein. Netw. Model. Analysis Health Inf. Bioinforma. 12, 1. doi:10.1007/s13721-022-00395-x
Alshammari, A., Alasmari, A. F., Alharbi, M., Ali, N., Muhseen, Z. T., Ashfaq, U. A., et al. (2022). Novel chimeric vaccine candidate development against leptotrichia buccalis. Int. J. Environ. Res. Public Health 19, 10742. doi:10.3390/ijerph191710742
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi:10.1093/nar/25.17.3389
Barford, D. (2004). The role of cysteine residues as redox-sensitive regulatory switches. Curr. Opin. Struct. Biol. 14, 679–686. doi:10.1016/j.sbi.2004.09.012
Bastola, R., Noh, G., Keum, T., Bashyal, S., Seo, J. E., Choi, J., et al. (2017). Vaccine adjuvants: smart components to boost the immune system. Archives Pharmacal Res. 40, 1238–1248. doi:10.1007/s12272-017-0969-z
Bettin, E. B., Dorneles, J., Hecktheuer, A. S., Madruga, A. B., Seixas Neto, A. C., Mcbride, A. J., et al. (2022). TonB-dependent receptor epitopes expressed in M. bovis BCG induced significant protection in the hamster model of leptospirosis. Appl. Microbiol. 106, 173–184. doi:10.1007/s00253-021-11726-9
Blanvillain, S., Meyer, D., Boulanger, A., Lautier, M., Guynet, C., Denancé, N., et al. (2007). Plant carbohydrate scavenging through TonB-dependent receptors: a feature shared by phytopathogenic and aquatic bacteria. PLoS One 2, e224. doi:10.1371/journal.pone.0000224
Brodin, P., and Davis, M. M. (2017). Human immune system variation. Nat. Rev. Immunol. 17, 21–29. doi:10.1038/nri.2016.125
Case, D. A. (1994). Normal mode analysis of protein dynamics. Curr. Opin. Struct. Biol. 4, 285–290. doi:10.1016/s0959-440x(94)90321-2
Castiglione, F., and Bernaschi, M. (2004). “C-immsim: playing with the immune response,” in Proceedings of the sixteenth international symposium on mathematical theory of networks and systems (MTNS2004) (Belgium: Katholieke Universiteit Leuven Belgium).
Celie, P. H., Parret, A. H., and Perrakis, A. J. C. O. I. S. B. (2016). Recombinant cloning strategies for protein expression. Curr. Opin. Struct. Biol. 38, 145–154. doi:10.1016/j.sbi.2016.06.010
Ceroni, D., Dubois-Ferrière, V., Cherkaoui, A., Lamah, L., Renzi, G., Lascombes, P., et al. (2013). 30 years of study of Kingella kingae: post tenebras, lux. Future Microbiol. 8, 233–245. doi:10.2217/fmb.12.144
Chauhan, V., Rungta, T., Goyal, K., and Singh, M. P. (2019). Designing a multi-epitope based vaccine to combat Kaposi Sarcoma utilizing immunoinformatics approach. Sci. Rep. 9, 2517. doi:10.1038/s41598-019-39299-8
Chavan, A. G. (2014). Exploring the molecular architecture of proteins: method developments in structure prediction and design. University of the Pacific, California: University of the Pacific.
Comeau, S. R., Gatchell, D. W., Vajda, S., and Camacho, C. J. (2004). ClusPro: a fully automated algorithm for protein–protein docking. Nucleic Acids Res. 32, W96–W99. doi:10.1093/nar/gkh354
Dhanda, S. K., Mahajan, S., Paul, S., Yan, Z., Kim, H., Jespersen, M. C., et al. (2019). IEDB-AR: immune epitope database—analysis resource in 2019. Nucleic Acids Res. 47, W502–W506. doi:10.1093/nar/gkz452
Dimitrov, I., Bangov, I., Flower, D. R., and Doytchinova, I. (2014). AllerTOP v. 2—a server for in silico prediction of allergens. J. Mol. Model. 20, 2278–2286. doi:10.1007/s00894-014-2278-5
Doytchinova, I. A., and Flower, D. R. (2007). VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinforma. 8, 4–7. doi:10.1186/1471-2105-8-4
Dubnov-Raz, G., Ephros, M., Garty, B. Z., Schlesinger, Y., Maayan-Metzger, A., Hasson, J., et al. (2010). Invasive pediatric Kingella kingae infections: a nationwide collaborative study. Pediatr. Infect. Dis. J. 29, 639–643. doi:10.1097/INF.0b013e3181d57a6c
Enayatkhani, M., Hasaniazad, M., Faezi, S., Gouklani, H., Davoodian, P., Ahmadi, N., et al. (2021). Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: an in silico study. J. Biomol. Struct. Dyn. 39, 2857–2872. doi:10.1080/07391102.2020.1756411
Finney, M., Vaughan, T., Taylor, S., Hudson, M. J., Pratt, C., Wheeler, J. X., et al. (2008). Characterization of the key antigenic components and pre-clinical immune responses to a meningococcal disease vaccine based on Neisseria lactamica outer membrane vesicles. Hum. Vaccines 4, 23–30. doi:10.4161/hv.4.1.4806
Gao, X., Ding, J., Liao, C., Xu, J., Liu, X., and Lu, W. (2021). Defensins: the natural peptide antibiotic. Adv. Drug Deliv. Rev. 179, 114008. doi:10.1016/j.addr.2021.114008
García-Angulo, V. A., Kalita, A., Kalita, M., Lozano, L., and Torres, A. G. (2014). Comparative genomics and immunoinformatics approach for the identification of vaccine candidates for enterohemorrhagic Escherichia coli O157: H7. Infect. Immun. 82, 2016–2026. doi:10.1128/IAI.01437-13
Glanfield, A., Mcmanus, D. P., Anderson, G. J., and Jones, M. K. (2007). Pumping iron: a potential target for novel therapeutics against schistosomes. Trends Parasitol. 23, 583–588. doi:10.1016/j.pt.2007.08.018
Grassmann, A. A., Kremer, F. S., Dos Santos, J. C., Souza, J. D., Pinto, L. D. S., and Mcbride, A. J. A. (2017). Discovery of novel leptospirosis vaccine candidates using reverse and structural vaccinology. Front. Immunol. 8, 463. doi:10.3389/fimmu.2017.00463
Greenawalt, A. (2021). Targeting transferrin binding protein A for vaccine development in Neisseria gonorrhoeae. Dissertation. Virginia Commonwealth University, Richmond: VCU University Archives.
Grote, A., Hiller, K., Scheer, M., Münch, R., Nörtemann, B., Hempel, D. C., et al. (2005). JCat: a novel tool to adapt codon usage of a target gene to its potential expression host. Nucleic Acids Res. 33, W526–W531. doi:10.1093/nar/gki376
Hu, Y. H., Dang, W., and Sun, L. (2012). A TonB-dependent outer membrane receptor of Pseudomonas fluorescens: virulence and vaccine potential. Archives Microbiol. 194, 795–802. doi:10.1007/s00203-012-0812-3
Ismail, C. M. K. H. B., Shabani, N. R. M., Chuah, C., Hassan, Z., Majeed, A. B. A., Leow, C. H., et al. (2022). Shigella iron-binding proteins: an insight into molecular physiology, pathogenesis, and potential target vaccine development. Vaccine 40, 3991–3998. doi:10.1016/j.vaccine.2022.05.061
Ismail, S., Ahmad, S., and Azam, S. S. (2020). Vaccinomics to design a novel single chimeric subunit vaccine for broad-spectrum immunological applications targeting nosocomial Enterobacteriaceae pathogens. Eur. J. Pharm. Sci. 146, 105258. doi:10.1016/j.ejps.2020.105258
Jin, M. S., Kim, S. E., Heo, J. Y., Lee, M. E., Kim, H. M., Paik, S. G., et al. (2007). Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell. 130, 1071–1082. doi:10.1016/j.cell.2007.09.008
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nat. Rev. Immunol. 596, 583–589. doi:10.1038/s41586-021-03819-2
Kaushik, V., Krishnan, S., Gupta, L. R., Kalra, U., Shaikh, A. R., Cavallo, L., et al. (2022). Immunoinformatics aided design and in-vivo validation of a cross-reactive peptide based multi-epitope vaccine targeting multiple serotypes of dengue virus. Front. Immunol. 13, 865180. doi:10.3389/fimmu.2022.865180
Kazi, A., Chuah, C., Majeed, A. B. A., Leow, C. H., Lim, B. H., and Leow, C. Y. (2018). Current progress of immunoinformatics approach harnessed for cellular-and antibody-dependent vaccine design. Pathog. Glob. Health 112, 123–131. doi:10.1080/20477724.2018.1446773
Khatoon, N., Pandey, R. K., and Prajapati, V. K. (2017). Exploring Leishmania secretory proteins to design B and T cell multi-epitope subunit vaccine using immunoinformatics approach. Sci. Rep. 7, 8285–8312. doi:10.1038/s41598-017-08842-w
Klitgaard, K., Friis, C., Angen, Ø., and Boye, M. (2010). Comparative profiling of the transcriptional response to iron restriction in six serotypes of Actinobacillus pleuropneumoniae with different virulence potential. BMC Genomics 11, 698–717. doi:10.1186/1471-2164-11-698
Kota, R. K., Kolla, H. B., Reddy, P. N., Kalagatur, N. K., and Samudrala, S. K. (2021). Immunoinformatics analysis and evaluation of recombinant chimeric triple antigen toxoid (r-HAB) against Staphylococcus aureus toxaemia in mouse model. Appl. Microbiol. Biotechnol. 105, 8297–8311. doi:10.1007/s00253-021-11609-z
Kozakov, D., Hall, D. R., Xia, B., Porter, K. A., Padhorny, D., Yueh, C., et al. (2017). The ClusPro web server for protein–protein docking. Nat. Protoc. 12, 255–278. doi:10.1038/nprot.2016.169
Laserson, U., Vigneault, F., Gadala-Maria, D., Yaari, G., Uduman, M., Vander Heiden, J. A., et al. (2014). High-resolution antibody dynamics of vaccine-induced immune responses. Proc. Natl. Acad. Sci. 111, 4928–4933. doi:10.1073/pnas.1323862111
Lebel, E., Rudensky, B., Karasik, M., Itzchaki, M., and Schlesinger, Y. (2006). Kingella kingae infections in children. J. Pediatr. Orthop. B 15, 289–292. doi:10.1097/01202412-200607000-00011
Lin, X., Chen, W., Xie, Q., Yu, Y., Liao, Y., Feng, Z., et al. (2022). Dissemination and genome analysis of high-level ceftriaxone-resistant penA 60.001 Neisseria gonorrhoeae strains from the Guangdong Gonococcal antibiotics susceptibility Programme (GD-GASP), 2016–2019. Emerg. Microbes Infect. 11, 344–350. doi:10.1080/22221751.2021.2011618
Lundy, D. W., and Kehl, D. K. (1998). Increasing prevalence of Kingella kingae in osteoarticular infections in young children. J. Pediatr. Orthop. B 18, 262–267. doi:10.1097/01241398-199803000-00025
Mahoney, R. T., and Maynard, J. E. (1999). The introduction of new vaccines into developing countries. Vaccine 17, 646–652. doi:10.1016/s0264-410x(98)00246-1
Maurer-Stroh, S., Krutz, N. L., Kern, P. S., Gunalan, V., Nguyen, M. N., Limviphuvadh, V., et al. (2019). AllerCatPro—prediction of protein allergenicity potential from the protein sequence. Bioinformatics 35, 3020–3027. doi:10.1093/bioinformatics/btz029
Mazumder, L., Shahab, M., Islam, S., Begum, M., Oliveira, J. I. N., Begum, S., et al. (2023). An immunoinformatics approach to epitope-based vaccine design against PspA in Streptococcus pneumoniae. J. Genet. Eng. Biotechnol. 21, 57–15. doi:10.1186/s43141-023-00506-9
Mishra, S., Rout, M., Panda, S., Singh, S. K., Sinha, R., Dehury, B., et al. (2022). An immunoinformatic approach towards development of a potent and effective multi-epitope vaccine against monkeypox virus (MPXV). J. Biomol. Struct. Dyn. 2022, 1–14. doi:10.1080/07391102.2022.2163426
Montero, M., Van Houten, N. E., Wang, X., and Scott, J. K. (2008). The membrane-proximal external region of the human immunodeficiency virus type 1 envelope: dominant site of antibody neutralization and target for vaccine design. Microbiol. Mol. Biol. Rev. 72, 54–84. doi:10.1128/MMBR.00020-07
Moynié, L., Luscher, A., Rolo, D., Pletzer, D., Tortajada, A., Weingart, H., et al. (2017). Structure and function of the PiuA and PirA siderophore-drug receptors from Pseudomonas aeruginosa and Acinetobacter baumannii. Antimicrob. Agents Chemother. 61, e02531–e02516. doi:10.1128/AAC.02531-16
Pancera, M., Changela, A., and Kwong, P. D. (2017). How HIV-1 entry mechanism and broadly neutralizing antibodies guide structure-based vaccine design. Curr. Opin. HIV AIDS 12, 229–240. doi:10.1097/COH.0000000000000360
Papaleo, E., Saladino, G., Lambrughi, M., Lindorff-Larsen, K., Gervasio, F. L., and Nussinov, R. (2016). The role of protein loops and linkers in conformational dynamics and allostery. Chem. Rev. 116, 6391–6423. doi:10.1021/acs.chemrev.5b00623
Ponomarenko, J. V., and Van Regenmortel, M. H. (2009). B cell epitope prediction. Struct. Bioinforma. 2, 849–879.
Porsch, E. A. (2022). Kingella kingae virulence factors and insights into pathogenicity. Microorganisms 10, 997. doi:10.3390/microorganisms10050997
Rahman, N., Ali, F., Basharat, Z., Shehroz, M., Khan, M. K., Jeandet, P., et al. (2020). Vaccine design from the ensemble of surface glycoprotein epitopes of SARS-CoV-2: an immunoinformatics approach. Vaccines 8, 423. doi:10.3390/vaccines8030423
Rahmani, A., Baee, M., Rostamtabar, M., Karkhah, A., Alizadeh, S., Tourani, M., et al. (2019). Development of a conserved chimeric vaccine based on helper T-cell and CTL epitopes for induction of strong immune response against Schistosoma mansoni using immunoinformatics approaches. Int. J. Biol. Macromol. 141, 125–136. doi:10.1016/j.ijbiomac.2019.08.259
Ramakrishna, L., Anand, K. K., Mohankumar, K. M., and Ranga, U. (2004). Codon optimization of the tat antigen of human immunodeficiency virus type 1 generates strong immune responses in mice following genetic immunization. J. Virology 78, 9174–9189. doi:10.1128/JVI.78.17.9174-9189.2004
Russo, G., Pennisi, M., Fichera, E., Motta, S., Raciti, G., Viceconti, M., et al. (2020). In silico trial to test COVID-19 candidate vaccines: a case study with UISS platform. BMC Bioinforma. 21, 527–616. doi:10.1186/s12859-020-03872-0
Saha, R., Ghosh, P., and Burra, V. P. (2021). Designing a next generation multi-epitope based peptide vaccine candidate against SARS-CoV-2 using computational approaches. 3 Biotech. 11, 47. doi:10.1007/s13205-020-02574-x
Saha, S., and Raghava, G. P. S. (2004). “BcePred: prediction of continuous B-cell epitopes in antigenic sequences using physico-chemical properties,” in Artificial immune systems: third international conference, ICARIS 2004 (Catania, Sicily, Italy: Springer), 197–204.
Saha, S., and Raghava, G. P. S. (2006). Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins Struct. Funct. Bioinforma. 65, 40–48. doi:10.1002/prot.21078
Saito, G., Swanson, J. A., and Lee, K. D. (2003). Drug delivery strategy utilizing conjugation via reversible disulfide linkages: role and site of cellular reducing activities. Adv. Drug Deliv. Rev. 55, 199–215. doi:10.1016/s0169-409x(02)00179-5
Scheiblhofer, S., Laimer, J., Machado, Y., Weiss, R., and Thalhamer, J. (2017). Influence of protein fold stability on immunogenicity and its implications for vaccine design. Expert Rev. Vaccines 16, 479–489. doi:10.1080/14760584.2017.1306441
Schryvers, A. B., and Stojiljkovic, I. (1999). Iron acquisition systems in the pathogenic Neisseria. Mol. Microbiol. 32, 1117–1123. doi:10.1046/j.1365-2958.1999.01411.x
Shelley, J. R., Davidson, D. J., and Dorin, J. R. (2020). The dichotomous responses driven by beta-defensins. Front. Immunol. 11, 1176. doi:10.3389/fimmu.2020.01176
Shen-Orr, S. S., and Furman, D. (2013). Variability in the immune system: of vaccine responses and immune states. Curr. Opin. Immunol. 25, 542–547. doi:10.1016/j.coi.2013.07.009
Taubert, K. A., and Dajani, A. S. (2001). Optimisation of the prevention and treatment of bacterial endocarditis. Drugs Aging 18, 415–424. doi:10.2165/00002512-200118060-00004
Tina, K., Bhadra, R., and Srinivasan, N. (2007). PIC: protein interactions calculator. Nucleic Acids Res. 35, W473–W476. doi:10.1093/nar/gkm423
Tiwari, S. P., Fuglebakk, E., Hollup, S. M., Skjærven, L., Cragnolini, T., Grindhaug, S. H., et al. (2014). WEBnm@ v2. 0: web server and services for comparing protein flexibility. BMC Bioinforma. 15, 1–12. doi:10.1186/s12859-014-0427-6
Trivedi, M. V., Laurence, J. S., and Siahaan, T. J. (2009). The role of thiols and disulfides on protein stability. Curr. Protein Peptide Sci. 10, 614–625. doi:10.2174/138920309789630534
Urakami, A., Ngwe Tun, M. M., Moi, M. L., Sakurai, A., Ishikawa, M., Kuno, S., et al. (2017). An envelope-modified tetravalent dengue virus-like-particle vaccine has implications for flavivirus vaccine design. J. Virology 91, e01181–e01117. doi:10.1128/JVI.01181-17
Van Zundert, G. C., Moriarty, N. W., Sobolev, O. V., Adams, P. D., and Borrelli, K. W. (2021). Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality. Structure 29, 913–921.e4. doi:10.1016/j.str.2021.03.011
Wang, D., Mai, J., Zhou, W., Yu, W., Zhan, Y., Wang, N., et al. (2020). Immunoinformatic analysis of T-and B-cell epitopes for SARS-CoV-2 vaccine design. Vaccines 8, 355. doi:10.3390/vaccines8030355
Wang, J., Xiong, K., Pan, Q., He, W., and Cong, Y. (2021). Application of TonB-dependent transporters in vaccine development of gram-negative bacteria. Front. Cell. Infect. Microbiol. 10, 589115. doi:10.3389/fcimb.2020.589115
Webb, D. C., and Cripps, A. W. (1999). Immunization with recombinant transferrin binding protein B enhances clearance of nontypeable Haemophilus influenzae from the rat lung. Infect. Immun. 67, 2138–2144. doi:10.1128/IAI.67.5.2138-2144.1999
Williams, N., Cooper, C., and Cundy, P. (2014). Kingella kingae septic arthritis in children: recognising an elusive pathogen. J. Children's Orthop. 8, 91–95. doi:10.1007/s11832-014-0549-4
Yang, F. L., Lou, T. C., Kuo, S. C., Wu, W. L., Chern, J., Lee, Y. T., et al. (2017). A medically relevant capsular polysaccharide in Acinetobacter baumannii is a potential vaccine candidate. Vaccine 35, 1440–1447. doi:10.1016/j.vaccine.2017.01.060
Zifruddin, A. N., Yusoff, M. A. M., Abd Ghani, N. S., Muhammad, N. A. N., Lam, K. W., and Hassan, M. (2023). Ensemble-based, high-throughput virtual screening of potential inhibitor targeting putative farnesol dehydrogenase of Metisa plana (Lepidoptera: psychidae). Comput. Biol. Chem. 103, 107811. doi:10.1016/j.compbiolchem.2023.107811
Keywords: Kingella kingae, immunoinformatics, chimeric vaccine, epitope mapping, bioinformatics
Citation: Mashraqi MM, Alzamami A, Alturki NA, Almasaudi HH, Ahmed I, Alshamrani S and Basharat Z (2023) Chimeric vaccine design against the conserved TonB-dependent receptor-like β-barrel domain from the outer membrane tbpA and hpuB proteins of Kingella kingae ATCC 23330. Front. Mol. Biosci. 10:1258834. doi: 10.3389/fmolb.2023.1258834
Received: 14 July 2023; Accepted: 12 October 2023;
Published: 20 November 2023.
Edited by:
Parimal Kar, Indian Institute of Technology Indore, IndiaReviewed by:
Kalyani Dhusia, University of Arkansas for Medical Sciences, United StatesBudheswar Dehury, Regional Medical Research Center (ICMR), India
Copyright © 2023 Mashraqi, Alzamami, Alturki, Almasaudi, Ahmed, Alshamrani and Basharat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Saleh Alshamrani, saalshamrani@nu.edu.sa; Zarrin Basharat, zarrin.iiui@gmail.com