Subtractive proteomics-guided vaccine targets identification and designing of multi-epitopes vaccine for immune response instigation against Burkholderia pseudomallei.

In this study, a methodical workflow using subtractive proteomics, vaccine designing, molecular simulation, and agent-based modeling approaches were used to annotate the whole proteome of Burkholderia pseudomallei (strain K96243) for vaccine designing. Among the total 5717 proteins in the whole proteome, 505 were observed to be essential for the pathogen's survival and pathogenesis predicted by the Database of Essential Genes. Among these, 23 vaccine targets were identified, of which fimbrial assembly chaperone (Q63UH5), Outer membrane protein (Q63UH1), and Hemolysin-like protein (Q63UE4) were selected for the subsequent analysis based on the systematic approaches. Using immunoinformatic approaches CTL (cytotoxic T lymphocytes), HTL (helper T lymphocytes), IFN-positive, and B cell epitopes were predicted for these targets. A total of 9 CTL epitopes were added using the GSS linker, 6 HTL epitopes using the GPGPG linker, and 6 B cell epitopes using the KK linker. An adjuvant was added for enhanced antigenicity, an HIV-TAT peptide for improved delivery, and a PADRE sequence was added to form a 466 amino acids long vaccine construct. The construct was classified as non-allergenic, highly antigenic, and experimentally feasible. Molecular docking results validated the robust interaction of MEVC with immune receptors such as TLR2/4. Furthermore, molecular simulation revealed stable dynamics and compact nature of the complexes. The binding free energy results further validated the robust binding. In silico cloning, results revealed GC contents of 50.73 % and a CIA value of 0.978 which shows proper downstream processing. Immune simulation results reported that after the three injections of the vaccine a robust secondary immune response, improved antigen clearance, and effective immune memory generation were observed highlighting its potential for effective and sustained immunity. Future directions should encompass experimental validations, animal model studies, and clinical trials to substantiate the vaccine's efficacy, safety, and immunogenicity.

Influence of Maqian essential oil on gut microbiota and immunoresponses in type 1 diabetes: In silico study.

Diversity and homeostasis of gut bacterial composition is highly associated with the pathogenesis of insulin dysfunction and type 1 diabetes melittus (T1D), hence emerged in parallel with the activation of autoimmunity. We aimed to study the bioactive potential of essential oil from Zanthoxylum myriacanthum var. pubescens Huang (Maqian) through computational approaches. Twelve chemical constituents derived from Maqian essential oil were docked with selected proteins (i.e., 3pig, 1kho, 7dmq, 4m4d, 2z65, 4glp, and 3fxi) in which are involved in gut microbiota modulation in T1D. Subsequently, the prediction of bioavailability properties of the small molecules were evaluated. Among all chemical constituents, the post-docking interaction analysis demonstrated that α-phellandrene exhibits the strongest binding affinity and induces gut microbiota modulation with ß-fructofuranosidase from Bifidobacterium longum. The current result revealed the potential of 3-Carene and α-Pinene in inducing specific changes in gut microbiota downregulating Clostridium perfringens and quenching Leptotrichia shahii respectively. ß-Pinene possess exceptionally strong binding affinity that effectively disrupt the interaction between lipopolysaccharide and its cognate receptors, while α-Phellandrene was exhibited the uppermost binding affinity with TLR4/MD2 and could likely target TLR4 stimulating lipopolysaccharide. Our results are the first to report on the gut microbiota modulation effects of α-Phellandrene and ß-Phellandrene via actions on LPS binding to CD14 and the TLR4 co-receptor signaling. In conclusion, our findings based on computational approaches, small molecules from Maqian present as promising agents which could regulate inflammatory response and modulate gut microbiota in type 1 diabetes mellitus.

Structural and molecular investigation of the impact of S30L and D88N substitutions in G9R protein on coupling with E4R from Monkeypox virus (MPXV).

Understanding the pathogenesis mechanism of the Monkeypox virus (MPXV) is essential to guide therapeutic development against the Monkeypox virus. In the current study, we investigated the impact of the only two reported substitutions, S30L, D88N, and S30L-D88N on the G9R of the replication complex in 2022 with E4R using structural modeling, simulation, and free energy calculation methods. From the molecular docking and dissociation constant (KD) results, it was observed that the binding affinity did not increase in the mutants, but the interaction paradigm was altered by these substitutions. Molecular simulation data revealed that these mutations are responsible for destabilization, changes in protein packing, and internal residue fluctuations, which can cause functional variance. Additionally, hydrogen bonding analysis revealed that the estimated number of hydrogen bonds are almost equal among the wild-type G9R and each mutant. The total binding free energy for the wild-type G9R with E4R was -85.00 kcal/mol while for the mutants the TBE was -42.75 kcal/mol, -43.68 kcal/mol, and -48.65 kcal/mol respectively. This shows that there is no direct impact of these two reported mutations on the binding with E4R, or it may affect the whole replication complex or any other mechanism involved in pathogenesis. To explore these variations further, we conducted PCA and FEL analyses. Based on our findings, we speculate that within the context of interaction with E4R, the mutations in the G9R protein might be benign, potentially leading to functional diversity associated with other proteins.Communicated by Ramaswamy H. Sarma.

Unraveling the mechanisms of Cefoxitin resistance in methicillin-resistant Staphylococcus aureus (MRSA): structural and molecular simulation-based insights.

Methicillin-resistant Staphylococcus aureus (MRSA) severely affects human health, including the skin glands, nasal cavity, wound infections, bone infections, and pneumonia. Among the most effective MRSA drugs, Cefoxitin also develops resistance due to mutations in the mecA gene. Four mutations at positions E229K, E239R, G246K, and E447K are classified as high-level resistance mutations. However, the resistance mechanism of MRSA towards Cefoxitin caused by these mutations is still unclear, as there is less information available regarding the structural and functional effects of the mutations against Cefoxitin. Therefore, our present study was designed to explore the mechanisms of binding interactions between wild-type and mutated PBP2a against Cefoxitin using molecular docking and MD simulations. Subsequently, we identified that the mutant form of PBP2a affects the activity of Cefoxitin. Interestingly, the binding of Cefoxitin with G246K and E239R mutants demonstrates unstable behavior compared to E447K-Cefoxitin and E229K-Cefoxitin. In this study, we propose the resistance mechanism of Cefoxitin at the atomic level. The proposed drug-resistance mechanism will provide valuable guidance for the design of MRSA drugs. This research might provide a new framework for designing new agents against the mutated form of PBP2a.Communicated by Ramaswamy H. Sarma.

Deciphering the mechanism of resistance by novel double mutations in pncA in Mycobacterium tuberculosis using protein structural graphs (PSG) and structural bioinformatic approaches.

The evolution of MDR and XDR-TB is a growing concern and public health safety threat around the world. Gene mutations are the prime cause of drug resistance in tuberculosis, however the reports of double mutations further aggravated the situation. Despite the large-scale genomic sequencing and identification of novel mutations, structure investigation of the protein is still required to structurally and functionally characterize these novel mutations to design novel drugs for improved clinical outcome. Hence, we used structural bioinformatics approaches i.e. molecular modeling, residues communication and molecular simulation to understand the impact of novel double S59Y-L85P, D86G-V180F and S104G-V130 M mutation on the structure, function of pncA encoded Pyrazinamidase (PZase) and resistance of Pyrazinamide (PZA). Our results revealed that these mutations alter the binding paradigm and destabilize the protein to release the drug. Protein commination network (PCN) revealed variations in the hub residues and sub-networks which consequently alter the internal communication and signaling. The region 1-75 demonstrated higher flexibility in the mutant structures and minimal by the wild type which destabilize of the internally arranged beta-sheets which consequently reduce the binding of PZA and potentially Fe ion in the mutants. Hydrogen bonding analysis further validated the findings. The total binding free energy (ΔG) for each complex i.e. wild type -7.46 kcal/mol, S59Y-L85P -5.21 kcal/mol, S104G-V130 M -5.33 kcal/mol while for the D86G-V180F mutant the TBE was calculated to be -6.26 kcal/mol. This further confirms that these mutations reduce the binding energy of PZA for PZase and causes resistance in the effective therapy for TB. The trajectories motion was also observed to be affected by these mutations. In conclusion, these mutations use destabilizing approach to reduce the binding of PZA and causes resistance. These features can be used to design novel structure-based drugs against Tuberculosis.

Pan-Genome Analysis of Oral Bacterial Pathogens to Predict a Potential Novel Multi-Epitopes Vaccine Candidate.

Porphyromonas gingivalis is a Gram-negative anaerobic bacterium, mainly present in the oral cavity and causes periodontal infections. Currently, no licensed vaccine is available against P. gingivalis and other oral bacterial pathogens. To develop a vaccine against P. gingivalis, herein, we applied a bacterial pan-genome analysis (BPGA) on the bacterial genomes that retrieved a total number of 4908 core proteins, which were further utilized for the identification of good vaccine candidates. After several vaccine candidacy analyses, three proteins, namely lytic transglycosylase domain-containing protein, FKBP-type peptidyl-propyl cis-trans isomerase and superoxide dismutase, were shortlisted for epitopes prediction. In the epitopes prediction phase, different types of B and T-cell epitopes were predicted and only those with an antigenic, immunogenic, non-allergenic, and non-toxic profile were selected. Moreover, all the predicted epitopes were joined with each other to make a multi-epitopes vaccine construct, which was linked further to the cholera toxin B-subunit to enhance the antigenicity of the vaccine. For downward analysis, a three dimensional structure of the designed vaccine was modeled. The modeled structure was checked for binding potency with major histocompatibility complex I (MHC-I), major histocompatibility complex II (MHC-II), and Toll-like receptor 4 (TLR-4) immune cell receptors which revealed that the designed vaccine performed proper binding with respect to immune cell receptors. Additionally, the binding efficacy of the vaccine was validated through a molecular dynamic simulation that interpreted strong intermolecular vaccine-receptor binding and confirmed the exposed situation of vaccine epitopes to the host immune system. In conclusion, the study suggested that the model vaccine construct has the potency to generate protective host immune responses and that it might be a good vaccine candidate for experimental in vivo and in vitro studies.

Discovery of Rift Valley fever virus natural pan-inhibitors by targeting its multiple key proteins through computational approaches.

The Rift Valley fever virus (RVFV) is a zoonotic arbovirus and pathogenic to both humans and animals. Currently, no proven effective RVFV drugs or licensed vaccine are available for human or animal use. Hence, there is an urgent need to develop effective treatment options to control this viral infection. RVFV glycoprotein N (GN), glycoprotein C (GC), and nucleocapsid (N) proteins are attractive antiviral drug targets due to their critical roles in RVFV replication. In present study, an integrated docking-based virtual screening of more than 6000 phytochemicals with known antiviral activities against these conserved RVFV proteins was conducted. The top five hit compounds, calyxin C, calyxin D, calyxin J, gericudranins A, and blepharocalyxin C displayed optimal binding against all three target proteins. Moreover, multiple parameters from the molecular dynamics (MD) simulations and MM/GBSA analysis confirmed the stability of protein-ligand complexes and revealed that these compounds may act as potential pan-inhibitors of RVFV replication. Our computational analyses may contribute toward the development of promising effective drugs against RVFV infection.

Structural Elucidation of Rift Valley Fever Virus L Protein towards the Discovery of Its Potential Inhibitors.

Rift valley fever virus (RVFV) is the causative agent of a viral zoonosis that causes a significant clinical burden in domestic and wild ruminants. Major outbreaks of the virus occur in livestock, and contaminated animal products or arthropod vectors can transmit the virus to humans. The viral RNA-dependent RNA polymerase (RdRp; L protein) of the RVFV is responsible for viral replication and is thus an appealing drug target because no effective and specific vaccine against this virus is available. The current study reported the structural elucidation of the RVFV-L protein by in-depth homology modeling since no crystal structure is available yet. The inhibitory binding modes of known potent L protein inhibitors were analyzed. Based on the results, further molecular docking-based virtual screening of Selleckchem Nucleoside Analogue Library (156 compounds) was performed to find potential new inhibitors against the RVFV L protein. ADME (Absorption, Distribution, Metabolism, and Excretion) and toxicity analysis of these compounds was also performed. Besides, the binding mechanism and stability of identified compounds were confirmed by a 50 ns molecular dynamic (MD) simulation followed by MM/PBSA binding free energy calculations. Homology modeling determined a stable multi-domain structure of L protein. An analysis of known L protein inhibitors, including Monensin, Mycophenolic acid, and Ribavirin, provide insights into the binding mechanism and reveals key residues of the L protein binding pocket. The screening results revealed that the top three compounds, A-317491, Khasianine, and VER155008, exhibited a high affinity at the L protein binding pocket. ADME analysis revealed good pharmacodynamics and pharmacokinetic profiles of these compounds. Furthermore, MD simulation and binding free energy analysis endorsed the binding stability of potential compounds with L protein. In a nutshell, the present study determined potential compounds that may aid in the rational design of novel inhibitors of the RVFV L protein as anti-RVFV drugs.

Immunoinformatics and Biophysics Approaches to Design a Novel Multi-Epitopes Vaccine Design against Staphylococcus auricularis.

Due to the misuse of antibiotics in our daily lives, antimicrobial resistance (AMR) has become a major health problem. Penicillin, the first antibiotic, was used in the 1930s and led to the emergence of AMR. Due to alterations in the microbe's genome and the evolution of new resistance mechanisms, antibiotics are losing efficacy against microbes. There are high rates of mortality and morbidity due to antibiotic resistance, so addressing this major health issue requires new approaches. Staphylococcus auricularis is a Gram-positive cocci and is capable of causing opportunistic infections and sepsis. S. auricularis is resistant to several antibiotics and does not currently have a licensed vaccine. In this study, we used bacterial pan-genome analysis (BPGA) to study S. auricularis pan-genome and applied a reverse immunology approach to prioritize vaccine targets against S. auricularis. A total of 15,444 core proteins were identified by BPGA analysis, which were then used to identify good vaccine candidates considering potential vaccine filters. Two vaccine candidates were evaluated for epitope prediction including the superoxide dismutase and gamma-glutamyl transferase protein. The epitope prediction phase involved the prediction of a variety of B-Cell and T-cell epitopes, and the epitopes that met certain criteria, such as antigenicity, immunogenicity, non-allergenicity, and non-toxicity were chosen. A multi-epitopes vaccine construct was then constructed from all the predicted epitopes, and a cholera toxin B-subunit adjuvant was also added to increase vaccine antigenicity. Three-dimensional models of the vaccine were used for downward analyses. Using the best-modeled structure, binding potency was tested with MHC-I, MHC-II and TLR-4 immune cells receptors, proving that the vaccine binds strongly with the receptors. Further, molecular dynamics simulations interpreted strong intermolecular binding between the vaccine and receptors and confirmed the vaccine epitopes exposed to the host immune system. The results support that the vaccine candidate may be capable of eliciting a protective immune response against S. auricularis and may be a promising candidate for experimental in vitro and in vivo studies.

Structural and Dynamic Insights into the W68L, L85P, and T87A Mutations of Mycobacterium tuberculosis Inducing Resistance to Pyrazinamide.

Tuberculosis (TB), the most frequent bacterium-mediated infectious disease caused by Mycobacterium tuberculosis, has been known to infect humans since ancient times. Although TB is common worldwide, the most recent report by the WHO (World Health Organization) listed the three countries of India, China, and Russia with 27%, 14%, and 8% of the global burden of TB, respectively. It has been reported that resistance to TB drugs, particularly by the pncA gene to the pyrazinamide drug due to mutations, significantly affects the effective treatment of TB. Understanding the mechanism of drug resistance using computational methods is of great interest to design effective TB treatment, exploring the structural features with these tools. Thus, keeping in view the importance of these methods, we employed state-of-the-art computational methods to study the mechanism of resistance caused by the W68L, L85P, and T87A mutations recently reported in 2021. We employed a molecular docking approach to predict the binding conformation and studied the dynamic properties of each complex using molecular dynamics simulation approaches. Our analysis revealed that compared to the wildtype, these three mutations altered the binding pattern and reduced the binding affinity. Moreover, the structural dynamic features also showed that these mutations significantly reduced the structural stability and packing, particularly by the W68L and L85P mutations. Moreover, principal component analysis, free energy landscape, and the binding free energy results revealed variation in the protein's motion and the binding energy. The total binding free energy was for the wildtype -9.61 kcal/mol, W68L -7.57 kcal/mol, L85P -6.99 kcal/mol, and T87A -7.77 kcal/mol. Our findings can help to design a structure-based drug against the MDR (multiple drug-resistant) TB.

Insights into the Binding of Receptor-Binding Domain (RBD) of SARS-CoV-2 Wild Type and B.1.620 Variant with hACE2 Using Molecular Docking and Simulation Approaches.

Recently, a new variant, B.1620, with mutations (S477N-E484K) in the spike protein's receptor-binding domain (RBD) has been reported in Europe. In order to design therapeutic strategies suitable for B.1.620, further studies are required. A detailed investigation of the structural features and variations caused by these substitutions, that is, a molecular level investigation, is essential to uncover the role of these changes. To determine whether and how the binding affinity of ACE2-RBD is affected, we used protein-protein docking and all-atom simulation approaches. Our analysis revealed that B.1.620 binds more strongly than the wild type and alters the hydrogen bonding network. The docking score for the wild type was reported to be -122.6 +/- 0.7 kcal/mol, while for B.1.620, the docking score was -124.9 +/- 3.8 kcal/mol. A comparative binding investigation showed that the wild-type complex has 11 hydrogen bonds and one salt bridge, while the B.1.620 complex has 14 hydrogen bonds and one salt bridge, among which most of the interactions are preserved between the wild type and B.1.620. A dynamic analysis of the two complexes revealed stable dynamics, which corroborated the global stability trend, compactness, and flexibility of the three essential loops, providing a better conformational optimization opportunity and binding. Furthermore, binding free energy revealed that the wild type had a total binding energy of -51.14 kcal/mol, while for B.1.628, the total binding energy was -68.25 kcal/mol. The current findings based on protein complex modeling and bio-simulation methods revealed the atomic features of the B.1.620 variant harboring S477N and E484K mutations in the RBD and the basis for infectivity. In conclusion, the current study presents distinguishing features of B.1.620, which can be used to design structure-based drugs against the B.1.620 variant.

Designing a Multi-Epitope Vaccine against Chlamydia trachomatis by Employing Integrated Core Proteomics, Immuno-Informatics and In Silico Approaches.

Chlamydia trachomatis, a Gram-negative bacterium that infects the rectum, urethra, congenital sites, and columnar epithelium of the cervix. It is a major cause of preventable blindness, ectopic pregnancy, and bacterial sexually transmitted infections worldwide. There is currently no licensed multi-epitope vaccination available for this pathogen. This study used core proteomics, immuno-informatics, and subtractive proteomics approaches to identify the best antigenic candidates for the development of a multi-epitope-based vaccine (MEBV). These approaches resulted in six vaccine candidates: Type III secretion system translocon subunit CopD2, SctW family type III secretion system gatekeeper subunit CopN, SycD/LcrH family type III secretion system chaperone Scc2, CT847 family type III secretion system effector, hypothetical protein CTDEC_0668, and CHLPN 76kDa-like protein. A variety of immuno-informatics tools were used to predict B and T cell epitopes from vaccine candidate proteins. An in silico vaccine was developed using carefully selected epitopes (11 CTL, 2 HTL & 10 LBL) and then docked with the MHC molecules (MHC I & MHC II) and human TLR4. The vaccine was coupled with Cholera toxin subunit B (CTB) adjuvant to boost the immune response. Molecular dynamics (MD) simulations, molecular docking, and MMGBSA analysis were carried out to analyze the molecular interactions and binding affinity of MEBV with TLR4 and MHC molecules. To achieve the highest level of vaccine protein expression, the MEBV was cloned and reverse-translated in Escherichia coli. The highest level of expression was achieved, and a CAI score of 0.97 was reported. Further experimental validation of the MEBV is required to prove its efficacy. The vaccine developed will be useful in preventing infections caused by C. trachomatis.

Core-Proteomics-Based Annotation of Antigenic Targets and Reverse-Vaccinology-Assisted Design of Ensemble Immunogen against the Emerging Nosocomial Infection-Causing Bacterium Elizabethkingia meningoseptica.

Elizabethkingia meningoseptica is a ubiquitous Gram-negative emerging pathogen that causes hospital-acquired infection in both immunocompromised and immunocompetent patients. It is a multi-drug-resistant bacterium; therefore, an effective subunit immunogenic candidate is of great interest to encounter the pathogenesis of this pathogen. A protein-wide annotation of immunogenic targets was performed to fast-track the vaccine development against this pathogen, and structural-vaccinology-assisted epitopes were predicted. Among the total proteins, only three, A0A1T3FLU2, A0A1T3INK9, and A0A1V3U124, were shortlisted, which are the essential vaccine targets and were subjected to immune epitope mapping. The linkers EAAK, AAY, and GPGPG were used to link CTL, HTL, and B-cell epitopes and an adjuvant was also added at the N-terminal to design a multi-epitope immunogenic construct (MEIC). The computationally predicted physiochemical properties of the ensemble immunogen reported a highly antigenic nature and produced multiple interactions with immune receptors. In addition, the molecular dynamics simulation confirmed stable binding and good dynamic properties. Furthermore, the computationally modeled immune response proposed that the immunogen triggered a strong immune response after several doses at different intervals. Neutralization of the antigen was observed on the 3rd day of injection. Conclusively, the immunogenic construct produces protection against Elizabethkingia meningoseptica; however, further immunological testing is needed to unveil its real efficacy.