MyxoPortal: a database of myxobacterial genomic features.

Myxobacteria are predatory bacteria with antimicrobial activity, utilizing complex mechanisms to kill their prey and assimilate their macromolecules. Having large genomes encoding hundreds of secondary metabolites, hydrolytic enzymes and antimicrobial peptides, these organisms are widely studied for their antibiotic potential. MyxoPortal is a comprehensive genomic database hosting 262 genomes of myxobacterial strains. Datasets included provide genome annotations with gene locations, functions, amino acids and nucleotide sequences, allowing analysis of evolutionary and taxonomical relationships between strains and genes. Biosynthetic gene clusters are identified by AntiSMASH, and dbAMP-generated antimicrobial peptide sequences are included as a resource for novel antimicrobial discoveries, while curated datasets of CRISPR/Cas genes, regulatory protein sequences, and phage associated genes give useful insights into each strain's biological properties. MyxoPortal is an intuitive open-source database that brings together application-oriented genomic features that can be used in taxonomy, evolution, predation and antimicrobial research. MyxoPortal can be accessed at http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Database URL:  http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Graphical Abstract.

Designing multi-epitope vaccine against human cytomegalovirus integrating pan-genome and reverse vaccinology pipelines.

Human cytomegalovirus (HCMV) is accountable for high morbidity in neonates and immunosuppressed individuals. Due to the high genetic variability of HCMV, current prophylactic measures are insufficient. In this study, we employed a pan-genome and reverse vaccinology approach to screen the target for efficient vaccine candidates. Four proteins, envelope glycoprotein M, UL41A, US23, and US28, were shortlisted based on cellular localization, high solubility, antigenicity, and immunogenicity. A total of 29 B-cell and 44 T-cell highly immunogenic and antigenic epitopes with high global population coverage were finalized using immunoinformatics tools and algorithms. Further, the epitopes that were overlapping among the finalized B-cell and T-cell epitopes were linked with suitable linkers to form various combinations of multi-epitopic vaccine constructs. Among 16 vaccine constructs, Vc12 was selected based on physicochemical and structural properties. The docking and molecular simulations of VC12 were performed, which showed its high binding affinity (-23.35 kcal/mol) towards TLR4 due to intermolecular hydrogen bonds, salt bridges, and hydrophobic interactions, and there were only minimal fluctuations. Furthermore, Vc12 eliciting a good response was checked for its expression in Escherichia coli through in silico cloning and codon optimization, suggesting it to be a potent vaccine candidate.

Integrating pan-genome and reverse vaccinology to design multi-epitope vaccine against Herpes simplex virus type-1.

Herpes simplex virus type-1 (HSV-1), the etiological agent of sporadic encephalitis and recurring oral (sometimes genital) infections in humans, affects millions each year. The evolving viral genome reduces susceptibility to existing antivirals and, thus, necessitates new therapeutic strategies. Immunoinformatics strategies have shown promise in designing novel vaccine candidates in the absence of a clinically licensed vaccine to prevent HSV-1. However, to encourage clinical translation, the HSV-1 pan-genome was integrated with the reverse-vaccinology pipeline for rigorous screening of universal vaccine candidates. Viral targets were screened from 104 available complete genomes. Among 364 proteins, envelope glycoprotein D being an outer membrane protein with a high antigenicity score (> 0.4) and solubility (> 0.6) was selected for epitope screening. A total of 17 T-cell and 4 B-cell epitopes with highly antigenic, immunogenic, non-toxic properties and high global population coverage were identified. Furthermore, 8 vaccine constructs were designed using different combinations of epitopes and suitable linkers. VC-8 was identified as the most potential vaccine candidate regarding chemical and structural stability. Molecular docking revealed high interactive affinity (low binding energy: - 56.25 kcal/mol) of VC-8 with the target elicited by firm intermolecular H-bonds, salt-bridges, and hydrophobic interactions, which was validated with simulations. Compatibility of the vaccine candidate to be expressed in pET-29(a) + plasmid was established by in silico cloning studies. Immune simulations confirmed the potential of VC-8 to trigger robust B-cell, T-cell, cytokine, and antibody-mediated responses, thereby suggesting a promising candidate for the future of HSV-1 prevention. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-04022-6.

Genomic investigation unveils colistin resistance mechanism in carbapenem-resistant Acinetobacter baumannii clinical isolates.

Colistin resistance in Acinetobacter baumannii is mediated by multiple mechanisms. Recently, mutations within pmrABC two-component system and overexpression of eptA gene due to upstream insertion of ISAba1 have been shown to play a major role. Thus, the aim of our study is to characterize colistin resistance mechanisms among the clinical isolates of A. baumannii in India. A total of 207 clinical isolates of A. baumannii collected from 2016 to 2019 were included in this study. Mutations within lipid A biosynthesis and pmrABC genes were characterized by whole-genome shotgun sequencing. Twenty-eight complete genomes were further characterized by hybrid assembly approach to study insertional inactivation of lpx genes and the association of ISAba1-eptA. Several single point mutations (SNPs), like M12I in pmrA, A138T and A444V in pmrB, and E117K in lpxD, were identified. We are the first to report two novel SNPs (T7I and V383I) in the pmrC gene. Among the five colistin-resistant A. baumannii isolates where complete genome was available, the analysis showed that three of the five isolates had ISAba1 insertion upstream of eptA. No mcr genes were identified among the isolates. We mapped the SNPs on the respective protein structures to understand the effect on the protein activity. We found that majority of the SNPs had little effect on the putative protein function; however, some SNPs might destabilize the local structure. Our study highlights the diversity of colistin resistance mechanisms occurring in A. baumannii, and ISAba1-driven eptA overexpression is responsible for colistin resistance among the Indian isolates.IMPORTANCEAcinetobacter baumannii is a Gram-negative, emerging and opportunistic bacterial pathogen that is often associated with a wide range of nosocomial infections. The treatment of these infections is hindered by increase in the occurrence of A. baumannii strains that are resistant to most of the existing antibiotics. The current drug of choice to treat the infection caused by A. baumannii is colistin, but unfortunately, the bacteria started to show resistance to the last-resort antibiotic. The loss of lipopolysaccharides and mutations in lipid A biosynthesis genes are the main reasons for the colistin resistance. The present study characterized 207 A. baumannii clinical isolates and constructed complete genomes of 28 isolates to recognize the mechanisms of colistin resistance. We showed the mutations in the colistin-resistant variants within genes essential for lipid A biosynthesis and that cause these isolates to lose the ability to produce lipopolysaccharides.

Multi-Epitope Vaccine for Monkeypox Using Pan-Genome and Reverse Vaccinology Approaches.

Outbreaks of monkeypox virus infections have imposed major health concerns worldwide, with high morbidity threats to children and immunocompromised adults. Although repurposed drugs and vaccines are being used to curb the disease, the evolving traits of the virus, exhibiting considerable genetic dynamicity, challenge the limits of a targeted treatment. A pan-genome-based reverse vaccinology approach can provide fast and efficient solutions to resolve persistent inconveniences in experimental vaccine design during an outbreak-exigency. The approach encompassed screening of available monkeypox whole genomes (n = 910) to identify viral targets. From 102 screened viral targets, viral proteins L5L, A28, and L5 were finalized based on their location, solubility, and antigenicity. The potential T-cell and B-cell epitopes were extracted from the proteins using immunoinformatics tools and algorithms. Multiple vaccine constructs were designed by combining the epitopes. Based on immunological properties, chemical stability, and structural quality, a novel multi-epitopic vaccine construct, V4, was finalized. Flexible-docking and coarse-dynamics simulation portrayed that the V4 had high binding affinity towards human HLA-proteins (binding energy < -15.0 kcal/mol) with low conformational fluctuations (<1 Å). Thus, the vaccine construct (V4) may act as an efficient vaccine to induce immunity against monkeypox, which encourages experimental validation and similar approaches against emerging viral infections.

R521C and R521H mutations in FUS result in weak binding with Karyopherinß2 leading to Amyotrophic lateral sclerosis: a molecular docking and dynamics study.

Fused in sarcoma (FUS) gene encodes the RNA binding protein FUS. This gene is mapped to chromosome 16p11.2. The FUS protein binds with karyopherineß2 (Kapß2) through its proline/tyrosine nuclear localization signal (PY-NLS) that helps in the localization of FUS protein within the nucleus. Arginine residue in 521 position (R521) of PY-NLS plays a vital role in the binding of FUS protein with Kapß2. Mutations in this position (R521C and R521H) are the most predominant mutations associated with amyotrophic lateral sclerosis (ALS). However, the mechanism by which these mutations lead to ALS is poorly understood. We examined the binding behaviour of the mutants FUS (R521C) and FUS (R521H) with Kapß2 through protein-protein docking and molecular dynamics simulation. The binding patterns of mutants were compared with the binding behaviour of wild FUS-Kapß2. Our results suggest that these mutants have relatively weak binding affinity with Kapß2 when compared with wild FUS-Kapß2 as indicated by the lesser number of interactions found between the mutant FUS and Kapß2. Hence, these mutations weakens the binding and this results in the cytoplasmic mislocalization of mutant FUS; and thereby it increases the severity of ALS.

Ebolavirus Database: Gene and Protein Information Resource for Ebolaviruses.

Ebola Virus Disease (EVD) is a life-threatening haemorrhagic fever in humans. Even though there are many reports on EVD, the protein precursor functions and virulent factors of ebolaviruses remain poorly understood. Comparative analyses of Ebolavirus genomes will help in the identification of these important features. This prompted us to develop the Ebolavirus Database (EDB) and we have provided links to various tools that will aid researchers to locate important regions in both the genomes and proteomes of Ebolavirus. The genomic analyses of ebolaviruses will provide important clues for locating the essential and core functional genes. The aim of EDB is to act as an integrated resource for ebolaviruses and we strongly believe that the database will be a useful tool for clinicians, microbiologists, health care workers, and bioscience researchers.

Identification of CD4+ T-cell epitope and investigation of HLA distribution for the immunogenic proteins of Burkholderia pseudomallei using in silico approaches - A key vaccine development strategy for melioidosis.

Melioidosis is a serious infectious diseases affecting multi-organ system in humans with high mortality rate. The disease is caused by the bacterium, Burkholderia pseudomallei and it is intrinsically resistant to many antibiotics. Thus, there is an urgent need for protective vaccine against B. pseudomallei; which may reduce morbidity and mortality in endemic areas. The identification of peptides that bind to major histocompatibility complex II class helps in understanding the nature of immune response and identifying T-cell epitopes for the design of new vaccines. Previous studies indicate that, ompA, bipB, fliC and groEL proteins of B. pseudomallei stimulate CD4+ T-cell immune response and act as protective immunogens. However, the data for CD4+ T-cell epitopes of these immunogenic proteins are very limited. Hence, in this present study we attempted to identify CD4+ T-cell epitopes in B. pseudomallei immunogenic proteins using in silico approaches. We did population coverage analysis for these identified epitopic core sequences to identify individuals in endemic areas expected to respond to a given set of these epitopes on the basis of HLA genotype frequencies. We observed that eight epitopic core sequences, two from each immunogenic protein, were associated with the maximum number of HLA-DR binding alleles. These eight peptides are found to be immunogenic in more than 90% of population in endemic areas considered. Thus, these eight peptides containing epitopic core sequences may act as probable vaccine candidates and they may be considered for the development of epitope-based vaccines for melioidosis.

Essentiality Assessment of Cysteinyl and Lysyl-tRNA Synthetases of Mycobacterium smegmatis.

Discovery of mupirocin, an antibiotic that targets isoleucyl-tRNA synthetase, established aminoacyl-tRNA synthetase as an attractive target for the discovery of novel antibacterial agents. Despite a high degree of similarity between the bacterial and human aminoacyl-tRNA synthetases, the selectivity observed with mupirocin triggered the possibility of targeting other aminoacyl-tRNA synthetases as potential drug targets. These enzymes catalyse the condensation of a specific amino acid to its cognate tRNA in an energy-dependent reaction. Therefore, each organism is expected to encode at least twenty aminoacyl-tRNA synthetases, one for each amino acid. However, a bioinformatics search for genes encoding aminoacyl-tRNA synthetases from Mycobacterium smegmatis returned multiple genes for glutamyl (GluRS), cysteinyl (CysRS), prolyl (ProRS) and lysyl (LysRS) tRNA synthetases. The pathogenic mycobacteria, namely, Mycobacterium tuberculosis and Mycobacterium leprae, were also found to possess two genes each for CysRS and LysRS. A similar search indicated the presence of additional genes for LysRS in gram negative bacteria as well. Herein, we describe sequence and structural analysis of the additional aminoacyl-tRNA synthetase genes found in M. smegmatis. Characterization of conditional expression strains of Cysteinyl and Lysyl-tRNA synthetases generated in M. smegmatis revealed that the canonical aminoacyl-tRNA synthetase are essential, while the additional ones are not essential for the growth of M. smegmatis.

Molecular Dynamics Studies on D835N Mutation in FLT3-Its Impact on FLT3 Protein Structure.

Mutations in Fetal Liver Tyrosine Kinase 3 (FLT3) genes are implicated in the constitutive activation and development of Acute Myeloid Leukaemia (AML). They are involved in signalling pathway of autonomous proliferation and block differentiation in leukaemia cells. FLT3 is considered as a promising target for the therapeutic intervention of AML. There are a few missense mutations associated with FLT3 that are found in AML patients. The D835N mutation is the most frequently observed and the aspartic acid in this position acts as a key residue for the receptor activation. The present study aims to understand the structural effect of D835N mutation in FLT3. We carried out the molecular dynamics (MD) simulation for a period of 120 ns at 300 K. Root-mean square deviation, root-mean square fluctuations, surface accessibility, radius of gyration, hydrogen bond, eigenvector projection analysis, trace of covariance matrix, and density analysis revealed the instability of mutant (D835N) protein. Our study provides new insights on the conformational changes in the mutant (D835N) structure of FLT3 protein. Our observations will be useful for researchers exploring AML and for the development of FLT3 inhibitors.

Tryptophan to Glycine mutation in the position 116 leads to protein aggregation and decreases the stability of the LITAF protein.

Mutations in the gene-encoding vesicle lipopolysaccharide-induced tumor necrosis factor (LITAF) protein cause Charcot-Marie-Tooth type 1C (CMT1C) disease, a neurological disorder. The LITAF gene is mapped to chromosome number 16 and can be found at cytogenetic location 16p13 of the chromosome. CMT1C-linked small integral membrane protein of lysosome/late endosome mutants are loss-of-function mutants that act in a dominant negative manner to impair endosomal trafficking, leading to prolonged extracellular signal-regulated kinases 1/2 signaling downstream of ErbB activation. Mutation W116G in the LITAF decreases the stability of the protein and also interrupts the functioning of gene. We have analyzed the single nucleotide polymorphism (SNP) results of 28 nsSNPs obtained from dbSNP. We also carried out multiple molecular dynamics simulations of 200 ns and obtained results of root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent-accessible surface area, H-bond, and principal component analysis to check and prove the stability of both the wild type and the mutant. The protein was then checked for its aggregation and the results showed loss of helix. The loss of helix leads to the instability of the protein.

Haemophilus influenzae Genome Database (HIGDB): a single point web resource for Haemophilus influenzae.

BACKGROUND: Haemophilus influenzae (H. Influenzae) is the causative agent of pneumonia, bacteraemia and meningitis. The organism is responsible for large number of deaths in both developed and developing countries. Even-though the first bacterial genome to be sequenced was that of H. Influenzae, there is no exclusive database dedicated for H. Influenzae. This prompted us to develop the Haemophilus influenzae Genome Database (HIGDB). METHODS: All data of HIGDB are stored and managed in MySQL database. The HIGDB is hosted on Solaris server and developed using PERL modules. Ajax and JavaScript are used for the interface development. RESULTS: The HIGDB contains detailed information on 42,741 proteins, 18,077 genes including 10 whole genome sequences and also 284 three dimensional structures of proteins of H. influenzae. In addition, the database provides "Motif search" and "GBrowse". The HIGDB is freely accessible through the URL: http://bioserver1.physics.iisc.ernet.in/HIGDB/. DISCUSSION: The HIGDB will be a single point access for bacteriological, clinical, genomic and proteomic information of H. influenzae. The database can also be used to identify DNA motifs within H. influenzae genomes and to compare gene or protein sequences of a particular strain with other strains of H. influenzae.

Computational Analysis Reveals the Association of Threonine 118 Methionine Mutation in PMP22 Resulting in CMT-1A.

The T118M mutation in PMP22 gene is associated with Charcot Marie Tooth, type 1A (CMT1A). CMT1A is a form of Charcot-Marie-Tooth disease, the most common inherited disorder of the peripheral nervous system. Mutations in CMT related disorder are seen to increase the stability of the protein resulting in the diseased state. We performed SNP analysis for all the nsSNPs of PMP22 protein and carried out molecular dynamics simulation for T118M mutation to compare the stability difference between the wild type protein structure and the mutant protein structure. The mutation T118M resulted in the overall increase in the stability of the mutant protein. The superimposed structure shows marked structural variation between the wild type and the mutant protein structures.

Streptococcus pneumoniae Genome Database (SPGDB): a database for strain specific comparative analysis of Streptococcus pneumoniae genes and proteins.

Streptococcus pneumoniae causes pneumonia, septicemia and meningitis. S. pneumoniae is responsible for significant mortality both in children and in the elderly. In recent years, the whole genome sequencing of various S. pneumoniae strains have increased manifold and there is an urgent need to provide organism specific annotations to the scientific community. This prompted us to develop the Streptococcus pneumoniae Genome Database (SPGDB) to integrate and analyze the completely sequenced and available S. pneumoniae genome sequences. Further, links to several tools are provided to compare the pool of gene and protein sequences, and proteins structure across different strains of S. pneumoniae. SPGDB aids in the analysis of phenotypic variations as well as to perform extensive genomics and evolutionary studies with reference to S. pneumoniae. The database will be updated at regular intervals and is freely accessible through the URL: http://pranag.physics.iisc.ernet.in/SPGDB/.