Functional annotation of hypothetical proteins from the Exiguobacterium antarcticum strain B7 reveals proteins involved in adaptation to extreme environments, including high arsenic resistance

Targeting Essential Hypothetical Proteins of Pseudomonas aeruginosa PAO1 for Mining of Novel Therapeutics: An In Silico Approach

As an omnipresent opportunistic bacterium, Pseudomonas aeruginosa PAO1 is responsible for acute and chronic infection in immunocompromised individuals. Currently, this bacterium is on WHO's red list where new antibiotics are urgently required for the treatment. Finding essential genes and essential hypothetical proteins (EHP) can be crucial in identifying novel druggable targets and therapeutics. This study is aimed at characterizing these EHPs and analyzing subcellular and physiochemical properties, PPI network, nonhomologous analysis against humans, virulence factor and novel drug target prediction, and finally structural analysis of the identified target employing around 42 robust bioinformatics tools/databases, the output of which was evaluated using the ROC analysis. The study discovered 18 EHPs from 336 essential genes, with domain and functional annotation revealing that 50% of these proteins belong to the enzyme category. The majority are cytoplasmic and cytoplasmic membrane proteins, with half being stable proteins subjected to PPIs network analysis. The network contains 261 nodes and 269 edges for 9 proteins of interest, with 11 hubs containing at least three nodes each. Finally, a pipeline builder predicts 7 proteins with novel drug targets, 5 nonhomologous proteins against human proteome, human antitargets, and human gut flora, and 3 virulent proteins. Among these, homology modeling of NP_249450 and NP_251676 was done, and the Ramachandran plot analysis revealed that more than 94% of the residues were in the preferred region. By analyzing functional attributes and virulence characteristics, the findings of this study may facilitate the development of innovative antibacterial drug targets and drugs of Pseudomonas aeruginosa PAO1.

Secretome Analysis of the Plant Biostimulant Bacteria Strains Bacillus subtilis (EB2004S) and Lactobacillus helveticus (EL2006H) in Response to pH Changes

It is well-known that there is a high frequency of plant-growth-promoting strains in Bacillus subtilis and that these can be effective under both stressful and stress-free conditions. There are very few studies of this activity in the case of Lactobacillus helveticus. In this study, the effects of pH on the secretome (proteins) in the cell-free supernatants of two bacterial strains were evaluated. The bacteria were cultured at pH 5, 7 and 8, and their secretome profiles were analyzed, with pH 7 (optimal growth pH) considered as the "control". The results showed that acidity (lower pH 5) diminishes the detectable production of most of the secretome proteins, whereas alkalinity (higher pH 8) increases the detectable protein production. At pH 5, five (5) new proteins were produced by L. helveticus, including class A sortase, fucose-binding lectin II, MucBP-domain-containing protein, SLAP-domain-containing protein and hypothetical protein LHEJCM1006_11110, whereas for B. subtilis, four (4) types of proteins were uniquely produced (p ≤ 0.05), including helicase-exonuclease AddAB subunit AddB, 5-methyltetrahydropteroyltriglutamate-homocysteine S-methyltransferase, a cluster of ABC-F family ATP-binding-cassette-domain-containing proteins and a cluster of excinuclease ABC (subunit B). At pH 8, Bacillus subtilis produced 56 unique proteins. Many of the detected proteins were involved in metabolic processes, whereas the others had unknown functions. The unique and new proteins with known and unknown functions suggest potential the acclimatization of the microbes to pH stress.

In-silico functional and structural annotation of hypothetical protein from Klebsiella pneumonia: A potential drug target

Klebsiella pneumonia is known to cause several nosocomial infections in immunocompromised patients. It has developed resistance against a broad range of presently available antibiotics, resulting in high mortality rates in patients and declared an urgent threat. Therefore, exploration of possible novel drug targets against this opportunistic bacteria needs to be undertaken. In the present study, we performed an extensive in-silico analysis for functional and structural annotation and characterized HP CP995_08280 from K. pneumonia as a drug target and aimed to identify potent drug candidates. The functional and structural studies using several bioinformatics tools and databases predicted that HP CP995_08280 is a cytosolic protein that belongs to the β-lactamase family and shares structural similarity with FmtA protein from Staphylococcus aureus (PDB ID: 5ZH8). The structure of HP CP995_08280 was successfully modeled followed by structure-based virtual screening, docking, molecular dynamics, and Molecular mechanic/Poisson-Boltzmann surface area (MMPBSA) were performed to identify the potential compounds. We have found five potent antibacterial molecules, namely BDD 24083171, BDD 24085737, BDE 25098678, BDE 33638819, and BDE 33672484, which exhibited high binding affinity (>-7.5 kcal/mol) and were stabilized by hydrogen bonding and hydrophobic interactions with active site residues (Ser42, Lys45, Tyr126, and Asp128) of protein. Molecular dynamics and MMPBSA revealed that HP CP995_08280 - ligand(s) complexes were less dynamic and more stable than native HP CP995_08280. Hence, the present study may serve as a potential lead for developing inhibitors against drug-resistant Klebsiella pneumonia.

In silico functional annotation of hypothetical proteins from the Bacillus paralicheniformis strain Bac84 reveals proteins with biotechnological potentials and adaptational functions to extreme environments

Members of the Bacillus genus are industrial cell factories due to their capacity to secrete significant quantities of biomolecules with industrial applications. The Bacillus paralicheniformis strain Bac84 was isolated from the Red Sea and it shares a close evolutionary relationship with Bacillus licheniformis. However, a significant number of proteins in its genome are annotated as functionally uncharacterized hypothetical proteins. Investigating these proteins' functions may help us better understand how bacteria survive extreme environmental conditions and to find novel targets for biotechnological applications. Therefore, the purpose of our research was to functionally annotate the hypothetical proteins from the genome of B. paralicheniformis strain Bac84. We employed a structured in-silico approach incorporating numerous bioinformatics tools and databases for functional annotation, physicochemical characterization, subcellular localization, protein-protein interactions, and three-dimensional structure determination. Sequences of 414 hypothetical proteins were evaluated and we were able to successfully attribute a function to 37 hypothetical proteins. Moreover, we performed receiver operating characteristic analysis to assess the performance of various tools used in this present study. We identified 12 proteins having significant adaptational roles to unfavorable environments such as sporulation, formation of biofilm, motility, regulation of transcription, etc. Additionally, 8 proteins were predicted with biotechnological potentials such as coenzyme A biosynthesis, phenylalanine biosynthesis, rare-sugars biosynthesis, antibiotic biosynthesis, bioremediation, and others. Evaluation of the performance of the tools showed an accuracy of 98% which represented the rationality of the tools used. This work shows that this annotation strategy will make the functional characterization of unknown proteins easier and can find the target for further investigation. The knowledge of these hypothetical proteins' potential functions aids B. paralicheniformis strain Bac84 in effectively creating a new biotechnological target. In addition, the results may also facilitate a better understanding of the survival mechanisms in harsh environmental conditions.

Exploration of the presence and abundance of multidrug resistance efflux genes in oil and gas environments

As sequencing technology improves and the cost of metagenome sequencing decreases, the number of sequenced environments increases. These metagenomes provide a wealth of data in the form of annotated and unannotated genes. The role of multidrug resistance efflux pumps (MDREPs) is the removal of antibiotics, biocides and toxic metabolites created during aromatic hydrocarbon metabolism. Due to their naturally occurring role in hydrocarbon metabolism and their role in biocide tolerance, MDREP genes are of particular importance for the protection of pipeline assets. However, the heterogeneity of MDREP genes creates a challenge during annotation and detection. Here we use a selection of primers designed to target MDREPs in six pure species and apply them to publicly available metagenomes associated with oil and gas environments. Using in silico PCR with relaxed primer binding conditions we probed the metagenomes of a shale reservoir, a heavy oil tailings pond, a civil wastewater treatment, two marine sediments exposed to hydrocarbons following the Deepwater Horizon oil spill and a non-exposed marine sediment to assess the presence and abundance of MDREP genes. Through relaxed primer binding conditions during in silico PCR, the prevalence of MDREPs was determined. The percentage of nucleotide sequences identified by the MDREP primers was partially augmented by exposure to hydrocarbons in marine sediment and in shale reservoir compared to hydrocarbon-free marine sediments while tailings ponds and wastewater had the highest percentages. We believe this approach lays the groundwork for a supervised method of identifying poorly conserved genes within metagenomes.

Functional Annotation of Hypothetical Proteins From the Enterobacter cloacae B13 Strain and Its Association With Pathogenicity

Enterobacter cloacae B13 strain is a rod-shaped gram-negative bacterium that belongs to the Enterobacteriaceae family. It can cause respiratory and urinary tract infections, and is responsible for several outbreaks in hospitals. E. cloacae has become an important pathogen and an emerging global threat because of its opportunistic and multidrug resistant ability. However, little knowledge is present about a large portion of its proteins and functions. Therefore, functional annotation of the hypothetical proteins (HPs) can provide an improved understanding of this organism and its virulence activity. The workflow in the study included several bioinformatic tools which were utilized to characterize functions, family and domains, subcellular localization, physiochemical properties, and protein-protein interactions. The E. cloacae B13 strain has overall 604 HPs, among which 78 were functionally annotated with high confidence. Several proteins were identified as enzymes, regulatory, binding, and transmembrane proteins with essential functions. Furthermore, 23 HPs were predicted to be virulent factors. These virulent proteins are linked to pathogenesis with their contribution to biofilm formation, quorum sensing, 2-component signal transduction or secretion. Better knowledge about the HPs’ characteristics and functions will provide a greater overview of the proteome. Moreover, it will help against E. cloacae in neonatal intensive care unit (NICU) outbreaks and nosocomial infections.

In Silico Characterization of Uncharacterized Proteins From Multiple Strains of Clostridium Difficile

Clostridium difficile (C. difficile) is a multi-strain, spore-forming, Gram-positive, opportunistic enteropathogen bacteria, majorly associated with nosocomial infections, resulting in severe diarrhoea and colon inflammation. Several antibiotics including penicillin, tetracycline, and clindamycin have been employed to control C. difficile infection, but studies have suggested that injudicious use of antibiotics has led to the development of resistance in C. difficile strains. However, many proteins from its genome are still considered uncharacterized proteins that might serve crucial functions and assist in the biological understanding of the organism. In this study, we aimed to annotate and characterise the 6 C. difficile strains using in silico approaches. We first analysed the complete genome of 6 C. difficile strains using standardised approaches and analysed hypothetical proteins (HPs) employing various bioinformatics approaches coalescing, including identifying contigs, coding sequences, phage sequences, CRISPR-Cas9 systems, antimicrobial resistance determination, membrane helices, instability index, secretory nature, conserved domain, and vaccine target properties like comparative homology analysis, allergenicity, antigenicity determination along with structure prediction and binding-site analysis. This study provides crucial supporting information about the functional characterization of the HPs involved in the pathophysiology of the disease. Moreover, this information also aims to assist in mechanisms associated with bacterial pathogenesis and further design candidate inhibitors and bona fide pharmaceutical targets.

In silico functional and structural characterization revealed virulent proteins of Francisella tularensis strain SCHU4

Francisella tularensis is a pathogenic, aerobic gram-negative coccobacillus bacterium. It is the causative agent of tularemia, a rare infectious disease that can attack skin, lungs, eyes, and lymph nodes. The genome of F. tularensis has been sequenced, and ~16% of the proteome is still uncharacterized. Characterizations of these proteins are essential to find new drug targets for better therapeutics. In silico characterization of proteins has become an extremely important approach to determine the functionality of proteins as experimental functional elucidation is unable to keep pace with the current growth of the sequence database. Initially, we have annotated 577 Hypothetical Proteins (HPs) of F. tularensis strain SCHU4 with seven bioinformatics tools which characterized them based on the family, domain and motif. Out of 577 HPs, 119 HPs were annotated by five or more tools and are further screened to predict their virulence properties, subcellular localization, transmembrane helices as well as physicochemical parameters. VirulentPred predicted 66 HPs out of 119 as virulent. These virulent proteins were annotated to find the interacting partner using STRING, and proteins with high confidence interaction scores were used to predict their 3D structures using Phyre2. The three virulent proteins Q5NH99 (phosphoserine phosphatase), Q5NG42 (Cystathionine beta-synthase) and Q5NG83 (Rrf2-type helix turn helix domain) were predicted to involve in modulation of cytoskeletal and innate immunity of host, H2S (hydrogen sulfide) based antibiotic tolerance and nitrite and iron metabolism of bacteria. The above predicted virulent proteins can serve as novel drug targets in the era of antibiotic resistance.

The Virulent Hypothetical Proteins: The Potential Drug Target Involved in Bacterial Pathogenesis

Abstract:

Hypothetical proteins (HPs) are non-predicted sequences that are identified only by open reading frames in sequenced genomes but their protein products remain uncharacterized by any experimental means. The genome of every species consists of HPs that are involved in various cellular processes and signaling pathways. Annotation of HPs is important as they play a key role in disease mechanisms, drug designing, vaccine production, antibiotic production, and host adaptation. In the case of bacteria, 25-50% of the genome comprises of HPs, which are involved in metabolic pathways and pathogenesis. The characterization of bacterial HPs helps to identify virulent proteins that are involved in pathogenesis. This can be done using in-silico studies, which provide sequence analogs, physiochemical properties, cellular or subcellular localization, structure and function validation, and protein-protein interactions. The most diverse types of virulent proteins are exotoxins, endotoxins, and adherent virulent factors that are encoded by virulent genes present on the chromosomal DNA of the bacteria. This review evaluates virulent HPs of pathogenic bacteria, such as Staphylococcus aureus, Chlamydia trachomatis, Fusobacterium nucleatum, and Yersinia pestis. The potential of these HPs as a drug target in bacteria-caused infectious diseases along with the mode of action and treatment approaches have been discussed.

Investigating the Functional Role of Hypothetical Proteins From an Antarctic Bacterium Pseudomonas sp. Lz4W: Emphasis on Identifying Proteins Involved in Cold Adaptation

Exploring the molecular mechanisms behind bacterial adaptation to extreme temperatures has potential biotechnological applications. In the present study, Pseudomonas sp. Lz4W, a Gram-negative psychrophilic bacterium adapted to survive in Antarctica, was selected to decipher the molecular mechanism underlying the cold adaptation. Proteome analysis of the isolates grown at 4°C was performed to identify the proteins and pathways that are responsible for the adaptation. However, many proteins from the expressed proteome were found to be hypothetical proteins (HPs), whose function is unknown. Investigating the functional roles of these proteins may provide additional information in the biological understanding of the bacterial cold adaptation. Thus, our study aimed to assign functions to these HPs and understand their role at the molecular level. We used a structured insilico workflow combining different bioinformatics tools and databases for functional annotation. Pseudomonas sp. Lz4W genome (CP017432, version 1) contains 4493 genes and 4412 coding sequences (CDS), of which 743 CDS were annotated as HPs. Of these, from the proteome analysis, 61 HPs were found to be expressed consistently at the protein level. The amino acid sequences of these 61 HPs were submitted to our workflow and we could successfully assign a function to 18 HPs. Most of these proteins were predicted to be involved in biological mechanisms of cold adaptations such as peptidoglycan metabolism, cell wall organization, ATP hydrolysis, outer membrane fluidity, catalysis, and others. This study provided a better understanding of the functional significance of HPs in cold adaptation of Pseudomonas sp. Lz4W. Our approach emphasizes the importance of addressing the “hypothetical protein problem” for a thorough understanding of mechanisms at the cellular level, as well as, provided the assessment of integrating proteomics methods with various annotation and curation approaches to characterize hypothetical or uncharacterized protein data. The MS proteomics data generated from this study has been deposited to the ProteomeXchange through PRIDE with the dataset identifier–PXD029741.

The Chironomid Microbiome Plays a Role in Protecting Its Host From Toxicants

Organisms are assemblages of the host and their endogenous bacteria, which are defined as microbiomes. The host and its microbiome undergo a mutual evolutionary process to adapt to changes in the environment. Chironomids (Diptera; Chironomidae), are aquatic insects that grow and survive in polluted environments; however, the mechanisms that protect them under these conditions are not fully understood. Here we present evidence that the chironomids’ microbiome enables them to survival in polluted environments. It has been demonstrated that about 40% of the microbiota that inhabit Chironomus transvaalensis egg masses and larvae has the potential to detoxify different toxicants. Metagenomic analysis of Chironomus ramosus larvae demonstrated the presence of genes in the insects’ microbiome that can help the insects to survive in hostile environments. A set of experiments demonstrated that short exposure of C. transvaalensis larvae to metals significantly changed their microbiota composition in comparison to unexposed larvae. Another experiment, that followed Koch’s postulates, demonstrated that disinfected C. transvaalensis larvae can survive toxic lead and chromium exposure when they are recolonized with bacteria that can detoxify these toxic metals. This accumulating research, points to the conclusion that the chironomid microbiome plays a role in protecting its host from toxicants.

Integrated bioinformatics based subtractive genomics approach to decipher the therapeutic function of hypothetical proteins from Salmonella typhi XDR H-58 strain

Purpose

The efficacy of drugs against Salmonella infection have compromised due to emerging XDR H58 strain. There is a dire need to find novel antimicrobial drug targets as well as drug candidates to cure by the XDR strain of Salmonella. It is observed that the complete genome sequence of the XDR H58 strain contains a large number of hypothetical proteins with unknown cellular and biological functions. Hence, it is indispensable to annotate these proteins functionally as well as structurally to identify novel drug targets.

Methods

In the current study, a comparative genomics and proteomics based approach was applied to find the novel drug targets in XDR strain while comparing the MDR and NR strains of Salmonella typhi.

Results

The characterization of ~ 350 hypothetical proteins were performed through determination of their physio-chemical properties, sub-cellular localization, functional annotation, and structure-based studies. As a result, only five proteins were prioritized as essential, druggable, and virulent proteins. Moreover, only one protein i.e. WP_000916613.1 was functionally annotated with high confidence and subjected to further structure-based analysis.

Conclusion

The current study presents a hypothetical protein from the XDR S. typhi proteome as a potential pharmacological target against which novel therapeutic candidates may be predicted. The outcome of the current study may lead to formulate a general set of pipelines for better understanding of the role of hypothetical proteins in pathogenesis of not only Salmonella but also for other pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10529-021-03219-6.

Functional Analysis of Hypothetical Proteins of Vibrio parahaemolyticus Reveals the Presence of Virulence Factors and Growth-Related Enzymes With Therapeutic Potential

Vibrio parahaemolyticus, an aquatic pathogen, is a major concern in the shrimp aquaculture industry. Several strains of this pathogen are responsible for causing acute hepatopancreatic necrosis disease as well as other serious illness, both of which result in severe economic losses. The genome sequence of two pathogenic strains of V. parahaemolyticus, MSR16 and MSR17, isolated from Bangladesh, have been reported to gain a better understanding of their diversity and virulence. However, the prevalence of hypothetical proteins (HPs) makes it challenging to obtain a comprehensive understanding of the pathogenesis of V. parahaemolyticus. The aim of the present study is to provide a functional annotation of the HPs to elucidate their role in pathogenesis employing several in silico tools. The exploration of protein domains and families, similarity searches against proteins with known function, gene ontology enrichment, along with protein-protein interaction analysis of the HPs led to the functional assignment with a high level of confidence for 656 proteins out of a pool of 2631 proteins. The in silico approach used in this study was important for accurately assigning function to HPs and inferring interactions with proteins with previously described functions. The HPs with function predicted were categorized into various groups such as enzymes involved in small-compound biosynthesis pathway, iron binding proteins, antibiotics resistance proteins, and other proteins. Several proteins with potential druggability were identified among them. In addition, the HPs were investigated in search of virulent factors, which led to the identification of proteins that have the potential to be exploited as vaccine candidate. The findings of the study will be effective in gaining a better understanding of the molecular mechanisms of bacterial pathogenesis. They may also provide an insight into the process of evaluating promising targets for the development of drugs and vaccines against V. parahaemolyticus.

Living to the High Extreme: Unraveling the Composition, Structure, and Functional Insights of Bacterial Communities Thriving in the Arsenic-Rich Salar de Huasco Altiplanic Ecosystem

Microbial communities inhabiting extreme environments such as Salar de Huasco (SH) in northern Chile are adapted to thrive while exposed to several abiotic pressures and the presence of toxic elements such as arsenic (As). Hence, we aimed to uncover the role of As in shaping bacterial composition, structure, and functional potential in five different sites in this altiplanic wetland using a shotgun metagenomic approach. The sites exhibit wide gradients of As (9 to 321 mg/kg), and our results showed highly diverse communities and a clear dominance exerted by the Proteobacteria and Bacteroidetes phyla. Functional potential analyses show broadly convergent patterns, contrasting with their great taxonomic variability. As-related metabolism, as well as other functional categories such as those related to the CH4 and S cycles, differs among the five communities. Particularly, we found that the distribution and abundance of As-related genes increase as the As concentration rises. Approximately 75% of the detected genes for As metabolism belong to expulsion mechanisms; arsJ and arsP pumps are related to sites with higher As concentrations and are present almost exclusively in Proteobacteria. Furthermore, taxonomic diversity and functional potential are reflected in the 12 reconstructed high-quality metagenome assembled genomes (MAGs) belonging to the Bacteroidetes (5), Proteobacteria (5), Cyanobacteria (1), and Gemmatimonadetes (1) phyla. We conclude that SH microbial communities are diverse and possess a broad genetic repertoire to thrive under extreme conditions, including increasing concentrations of highly toxic As. Finally, this environment represents a reservoir of unknown and undescribed microorganisms, with great metabolic versatility, which needs further study. IMPORTANCE As microbial communities inhabiting extreme environments are fundamental for maintaining ecosystems, many studies concerning composition, functionality, and interactions have been carried out. However, much is still unknown. Here, we sampled microbial communities in the Salar de Huasco, an extreme environment subjected to several abiotic stresses (high UV radiation, salinity and arsenic; low pressure and temperatures). We found that although microbes are taxonomically diverse, functional potential seems to have an important degree of convergence, suggesting high levels of adaptation. Particularly, arsenic metabolism showed differences associated with increasing concentrations of the metalloid throughout the area, and it effectively exerts a significant pressure over these organisms. Thus, the significance of this research is that we describe highly specialized communities thriving in little-explored environments subjected to several pressures, considered analogous of early Earth and other planets, that have the potential for unraveling technologies to face the repercussions of climate change in many areas of interest.

Screening and characterization of hypothetical proteins of Plasmodium falciparum as novel vaccine candidates in the fight against malaria using reverse vaccinology

BACKGROUND

Plasmodium falciparum is the most deadly and leading cause of morbidity and mortality in Africa. About 90% of all malaria deaths in the world today occur in Sub-Saharan Africa especially in children aged < 5 years. In 2018, it was reported that there were 228 million malaria cases that resulted in 405,000 deaths from 91 countries. Currently, a fully effective and long-lasting preventive malaria vaccine is still elusive therefore more effort is needed to identify better effective vaccine candidates. The aim of this study was to identify and characterize hypothetical proteins as vaccine candidates derived from Plasmodium falciparum 3D7 genome by reverse vaccinology.

RESULTS

Of the 23 selected hypothetical proteins, 5 were predicted on the extracellular localization by WoLFPSORTv.2.0 program and all the 5 had less than 2 transmembrane regions that were predicted by TMHMMv2.0 and HMMTOP programs at default settings. Two out of the five proteins lacked secretory signal peptides as predicted by SignalP program. Among the 5 extracellular proteins, 3 were predicted to be antigenic by VaxiJen (score ≥ 0.5) and had negative GRAVY values ranging from - 1.156 to - 0.440. B cell epitope prediction by ABCpred and BCpred programs revealed a total of 15 antigenic epitopes. A total of 13 cytotoxic T cells were predicted from the 3 proteins using CTLPred online server. Only 2 out of the 13 CTL were antigenic, immunogenic, non-allergenic, and non-toxic using VaxiJen, IEDB, AllergenFp, and Toxinpred servers respectively in that order. Five HTL peptides from XP_001351030.1 protein are predicted inducers of all the three cytokines. STRING protein-protein network analysis of HPs revealed XP_001350955.1 closely interacts with nucleoside diphosphate kinase (PF13-0349) at 0.704, XP_001351030.1 interacts with male development protein1 (Mdv-1) at 0.645, and XP_001351047.1 with an uncharacterized protein (MAL8P1.53) at 0.400.

CONCLUSION

Reverse vaccinology is a promising strategy for the screening and identification of antigenic antigens with potential capacity to elicit cellular and humoral immune responses against P. falciparum infection. In this study, potential vaccine candidates of Plasmodium falciparum were identified and screened using standard bioinformatics tools. The vaccine candidates contained antigenic and immunogenic epitopes which could be considered for novel and effective vaccine targets. However, we strongly recommend in vivo and in vitro experiments to validate their immunogenicity and protective efficacy to completely decipher the vaccine targets against malaria.

Questing functions and structures of hypothetical proteins from Campylobacter jejuni: a computer-aided approach

Campylobacter jejuni (C. jejuni) is considered to be one of the most frequent causes of bacterial gastroenteritis globally, especially in young children. The genome of C. jejuni contains many proteins with unknown functions termed as hypothetical proteins (HPs). These proteins might have essential biological role to show the full spectrum of this bacterium. Hence, our study aimed to determine the functions of HPs, pertaining to the genome of C. jejuni. An in-silico work flow integrating various tools were performed for functional assignment, three-dimensional structure determination, domain architecture predictors, subcellular localization, physicochemical characterization, and protein-protein interactions (PPIs). Sequences of 267 HPs of C. jejuni were analyzed and successfully attributed the function of 49 HPs with higher confidence. Here, we found proteins with enzymatic activity, transporters, binding and regulatory proteins as well as proteins with biotechnological interest. Assessment of the performance of various tools used in this analysis revealed an accuracy of 95% using receiver operating characteristic (ROC) curve analysis. Functional and structural predictions and the results from ROC analyses provided the validity of in-silico tools used in the present study. The approach used for this analysis leads us to assign the function of unknown proteins and relate them with the functions that have already been described in previous literature.

Genomics Insights into Pseudomonas sp. CG01: An Antarctic Cadmium-Resistant Strain Capable of Biosynthesizing CdS Nanoparticles Using Methionine as S-Source

Here, we present the draft genome sequence of Pseudomonas sp. GC01, a cadmium-resistant Antarctic bacterium capable of biosynthesizing CdS fluorescent nanoparticles (quantum dots, QDs) employing a unique mechanism involving the production of methanethiol (MeSH) from methionine (Met). To explore the molecular/metabolic components involved in QDs biosynthesis, we conducted a comparative genomic analysis, searching for the genes related to cadmium resistance and sulfur metabolic pathways. The genome of Pseudomonas sp. GC01 has a 4,706,645 bp size with a 58.61% G+C content. Pseudomonas sp. GC01 possesses five genes related to cadmium transport/resistance, with three P-type ATPases (cadA, zntA, and pbrA) involved in Cd-secretion that could contribute to the extracellular biosynthesis of CdS QDs. Furthermore, it exhibits genes involved in sulfate assimilation, cysteine/methionine synthesis, and volatile sulfur compounds catabolic pathways. Regarding MeSH production from Met, Pseudomonas sp. GC01 lacks the genes E4.4.1.11 and megL for MeSH generation. Interestingly, despite the absence of these genes, Pseudomonas sp. GC01 produces high levels of MeSH. This is probably associated with the metC gene that also produces MeSH from Met in bacteria. This work is the first report of the potential genes involved in Cd resistance, sulfur metabolism, and the process of MeSH-dependent CdS QDs bioproduction in Pseudomonas spp. strains.

An Educational Bioinformatics Project to Improve Genome Annotation

Scientific advancement is hindered without proper genome annotation because biologists lack a complete understanding of cellular protein functions. In bacterial cells, hypothetical proteins (HPs) are open reading frames with unknown functions. HPs result from either an outdated database or insufficient experimental evidence (i.e., indeterminate annotation). While automated annotation reviews help keep genome annotation up to date, often manual reviews are needed to verify proper annotation. Students can provide the manual review necessary to improve genome annotation. This paper outlines an innovative classroom project that determines if HPs have outdated or indeterminate annotation. The Hypothetical Protein Characterization Project uses multiple well-documented, freely available, web-based, bioinformatics resources that analyze an amino acid sequence to (1) detect sequence similarities to other proteins, (2) identify domains, (3) predict tertiary structure including active site characterization and potential binding ligands, and (4) determine cellular location. Enough evidence can be generated from these analyses to support re-annotation of HPs or prioritize HPs for experimental examinations such as structural determination via X-ray crystallography. Additionally, this paper details several approaches for selecting HPs to characterize using the Hypothetical Protein Characterization Project. These approaches include student- and instructor-directed random selection, selection using differential gene expression from mRNA expression data, and selection based on phylogenetic relations. This paper also provides additional resources to support instructional use of the Hypothetical Protein Characterization Project, such as example assignment instructions with grading rubrics, links to training videos in YouTube, and several step-by-step example projects to demonstrate and interpret the range of achievable results that students might encounter. Educational use of the Hypothetical Protein Characterization Project provides students with an opportunity to learn and apply knowledge of bioinformatic programs to address scientific questions. The project is highly customizable in that HP selection and analysis can be specifically formulated based on the scope and purpose of each student's investigations. Programs used for HP analysis can be easily adapted to course learning objectives. The project can be used in both online and in-seat instruction for a wide variety of undergraduate and graduate classes as well as undergraduate capstone, honor's, and experiential learning projects.

Deciphering the functional role of hypothetical proteins from Chloroflexus aurantiacs J-10-f1 using bioinformatics approach

Chloroflexus aurantiacus J-10-f1 is an anoxygenic, photosynthetic, facultative autotrophic gram negative bacterium found from hot spring at a temperature range of 50-60°C. It can sustain itself in dark only if oxygen is available thereby exhibiting a dark orange color, however display a dark green color when grown in sunlight. Genome of the organism contains total of 3853 proteins out of which 785 (~20%) proteins are uncharacterised or hypothetical proteins (HPs). Therefore in this work we have characterized the 785 hypothetical proteins of Chloroflexus aurantiacus J-10-f1 using bioinformatics tools and databases. HPs annotated by more than five domain prediction tools were filtered and named high confidence-hypothetical proteins (HC-HPs). These HC-HPs were further annotated by calculating their physiochemical properties, homologous, subcellular locations, signal peptides and transmembrane regions. We found most of the HC-HPs were involved in photosynthesis, carbohydrate metabolism, biofuel production and cellulose synthesis processes. Furthermore, few of these HC-HPs could provide resistance to bacteria at high temperature due to their thermophilic nature. Hence these HC-HPs have the potential to be used in industrial as well as in biomedical needs. To conclude, the bioinformatics approach used in this study provides an insight to better understand the nature and role of Chloroflexus aurantiacus J-10-f1 hypothetical proteins.

In silico functional prediction of hypothetical proteins from the core genome of Corynebacterium pseudotuberculosis biovar ovis

Corynebacterium pseudotuberculosis is a pathogen of veterinary relevance diseases, being divided into two biovars: equi and ovis; causing ulcerative lymphangitis and caseous lymphadenitis, respectively. The isolation and sequencing of C. pseudotuberculosis biovar ovis strains in the Northern and Northeastern regions of Brazil exhibited the emergence of this pathogen, which causes economic losses to small ruminant producers, and condemnation of carcasses and skins of animals. Through the pan-genomic approach, it is possible to determine and analyze genes that are shared by all strains of a species—the core genome. However, many of these genes do not have any predicted function, being characterized as hypothetical proteins (HP). In this study, we considered 32 C. pseudotuberculosis biovar ovis genomes for the pan-genomic analysis, where were identified 172 HP present in a core genome composed by 1255 genes. We are able to functionally annotate 80 sequences previously characterized as HP through the identification of structural features as conserved domains and families. Furthermore, we analyzed the physicochemical properties, subcellular localization and molecular function. Additionally, through RNA-seq data, we investigated the differential gene expression of the annotated HP. Genes inserted in pathogenicity islands had their virulence potential evaluated. Also, we have analyzed the existence of functional associations for their products based on protein–protein interaction networks, and perform the structural prediction of three targets. Due to the integration of different strategies, this study can underlie deeper in vitro researches in the characterization of these HP and the search for new solutions for combat this pathogen.

In silico Functional Annotation and Characterization of Hypothetical Proteins from Serratia marcescens FGI94

Serratia marcescens, rod-shaped Gram-negative bacteria is classified as an opportunistic pathogen in the family Enterobacteriaceae. It causes a wide variety of infections in humans, including urinary, respiratory, ocular lens and ear infections, osteomyelitis, endocarditis, meningitis and septicemia. Unfortunately, over the past decade, antibiotic resistance has become a serious health care issue; the effective means to control and dissemination of S. marcescens resistance is the need of hour. The whole genome sequencing of S. marcescens FGI94 strain contains 4434 functional proteins, among which 690 (15.56%) proteins were classified under hypothetical. In the present study, we applied the power of various bioinformatics tools on the basis of protein family comparison, motifs, functional properties of amino acids and genome context to assign the possible functions for the HPs. The pseudo sequences (protein sequence that contain ≤100 amino acid residues) are eliminated from the study. Although we have successfully predicted the function for 483 proteins, we were able to infer the high level of confidence only for 108 proteins. The predicted HPs were classified into various classes such as enzymes, transporters, binding proteins, cell division, cell regulatory and other proteins. The outcome of the study could be helpful to understand the molecular mechanism in bacterial pathogenesis and also provide an insight into the identification of potential targets for drug and vaccine development.

Genomic Variation and Arsenic Tolerance Emerged as Niche Specific Adaptations by Different Exiguobacterium Strains Isolated From the Extreme Salar de Huasco Environment in Chilean - Altiplano

Polyextremophilic bacteria can thrive in environments with multiple stressors such as the Salar de Huasco (SH). Microbial communities in SH are exposed to low atmospheric pressure, high UV radiation, wide temperature ranges, salinity gradient and the presence of toxic compounds such as arsenic (As). In this work we focus on arsenic stress as one of the main adverse factors in SH and bacteria that belong to the Exiguobacterium genus due to their plasticity and ubiquity. Therefore, our aim was to shed light on the effect of niche conditions pressure (particularly arsenic), on the adaptation and divergence (at genotypic and phenotypic levels) of Exiguobacterium strains from five different SH sites. Also, to capture greater diversity in this genus, we use as outgroup five As(III) sensitive strains isolated from Easter Island (Chile) and The Great Salt Lake (United States). For this, samples were obtained from five different SH sites under an arsenic gradient (9 to 321 mg/kg: sediment) and isolated and sequenced the genomes of 14 Exiguobacterium strains, which had different arsenic tolerance levels. Then, we used comparative genomic analysis to assess the genomic divergence of these strains and their association with phenotypic differences such as arsenic tolerance levels and the ability to resist poly-stress. Phylogenetic analysis showed that SH strains share a common ancestor. Consequently, populations were separated and structured in different SH microenvironments, giving rise to multiple coexisting lineages. Hence, this genotypic variability is also evidenced by the COG (Clusters of Orthologous Groups) composition and the size of their accessory genomes. Interestingly, these observations correlate with physiological traits such as growth patterns, gene expression, and enzyme activity related to arsenic response and/or tolerance. Therefore, Exiguobacterium strains from SH are adapted to physiologically overcome the contrasting environmental conditions, like the arsenic present in their habitat.

Toward a chimeric vaccine against multiple isolates of Mycobacteroides - An integrative approach

AIM

Nontuberculous mycobacterial (NTM) infection such as endophthalmitis, dacryocystitis, and canaliculitis are pervasive across the globe and are currently managed by antibiotics. However, the recent cases of Mycobacteroides developing drug resistance reported along with the improper practice of medicine intrigued us to explore its genomic and proteomic canvas at a global scale and develop a chimeric vaccine against Mycobacteroides.

MAIN METHODS

We carried out a vivid genomic study on five recently sequenced strains of Mycobacteroides and explored their Pan-core genome/proteome in three different phases. The promiscuous antigenic proteins were identified via a subtractive proteomics approach that qualified for virulence causation, resistance and essentiality factors for this notorious bacterium. An integrated pipeline was developed for the identification of B-Cell, MHC (Major histocompatibility complex) class I and II epitopes.

KEY FINDINGS

Phase I identified the shreds of evidence of reductive evolution and propensity of the Pan-genome of Mycobacteroides getting closed soon. Phase II and Phase III produced 8 vaccine constructs. Our final vaccine construct, V6 qualified for all tests such as absence for allergenicity, presence of antigenicity, etc. V6 contains β-defensin as an adjuvant, linkers, Lysosomal-associated membrane protein 1 (LAMP1) signal peptide, and PADRE (Pan HLA-DR epitopes) amino acid sequence. Besides, V6 also interacts with a maximum number of MHC molecules and the TLR4/MD2 (Toll-like receptor 4/Myeloid differentiation factor 2) complex confirmed by docking and molecular dynamics simulation studies.

SIGNIFICANCE

The knowledge harnessed from the current study can help improve the current treatment regimens or in an event of an outbreak and propel further related studies.

Seasonal dynamics of Chironomus transvaalensis populations and the microbial community composition of their egg masses

Chironomids (Diptera; Chironomidae) are the most abundant insects in freshwater environments and are considered natural reservoirs of Vibrio cholerae. We monitored the annual dynamics of chironomid populations along with their microbiota in order to better understand host–microbiota interactions. Chironomus transvaalensis populations peaked biannually in August and May–June. The composition of the endogenous bacterial communities of their egg masses clustered in two groups according to the sampling periods August–November and May–July. Nevertheless, a core bacterial community (43%) was present in all egg-mass samples. The most abundant phyla were: Proteobacteria, Firmicutes, Cyanobacteria and Bacteroidetes. The abundance of several genera (e.g. Rheinheimera and Pseudomonas) was positively correlated with C. transvaalensis population dynamics, while a predator–prey interaction was observed between the relative abundance of Vibrio OTUs and C. transvaalensis population size. Chironomids are known to tolerate toxic and stress conditions, and our results demonstrated that bacterial genera that may protect the insect under these conditions are present in the egg masses. After hatching, the first larval meal is the gelatinous matrix that surrounds the eggs. This meal contains a probiotic consortium that may protect the larva during its metamorphosis. The results provide important insights into the host–microbe interactions of chironomids.

Comparative Genomics Reveals Metabolic Specificity of Endozoicomonas Isolated from a Marine Sponge and the Genomic Repertoire for Host-Bacteria Symbioses

The most recently described bacterial members of the genus Endozoicomonas have been found in association with a wide variety of marine invertebrates. Despite their ubiquity in the host holobiont, limited information is available on the molecular genomic signatures of the symbiotic association of Endozoicomonas with marine sponges. Here, we generated a draft genome of Endozoicomonas sp. OPT23 isolated from the intertidal marine sponge Ophlitaspongia papilla and performed comprehensive comparative genomics analyses. Genome-specific analysis and metabolic pathway comparison of the members of the genus Endozoicomonas revealed the presence of gene clusters encoding for unique metabolic features, such as the utilization of carbon sources through lactate, L-rhamnose metabolism, and a phenylacetic acid degradation pathway in Endozoicomonas sp. OPT23. Moreover, the genome harbors genes encoding for eukaryotic-like proteins, such as ankyrin repeats, tetratricopeptide repeats, and Sel1 repeats, which likely facilitate sponge-bacterium attachment. The genome also encodes major secretion systems and homologs of effector molecules that seem to enable the sponge-associated bacterium to interact with the sponge and deliver the virulence factors for successful colonization. In conclusion, the genome analysis of Endozoicomonas sp. OPT23 revealed the presence of adaptive genomic signatures that might favor their symbiotic lifestyle within the sponge host.

Current Status and Potential Applications of Underexplored Prokaryotes

Thousands of prokaryotic genera have been published, but methodological bias in the study of prokaryotes is noted. Prokaryotes that are relatively easy to isolate have been well-studied from multiple aspects. Massive quantities of experimental findings and knowledge generated from the well-known prokaryotic strains are inundating scientific publications. However, researchers may neglect or pay little attention to the uncommon prokaryotes and hard-to-cultivate microorganisms. In this review, we provide a systematic update on the discovery of underexplored culturable and unculturable prokaryotes and discuss the insights accumulated from various research efforts. Examining these neglected prokaryotes may elucidate their novelties and functions and pave the way for their industrial applications. In addition, we hope that this review will prompt the scientific community to reconsider these untapped pragmatic resources.

Arsenic Response of Three Altiplanic Exiguobacterium Strains With Different Tolerance Levels Against the Metalloid Species: A Proteomics Study

Exiguobacterium is a polyextremophile bacterial genus with a physiology that allows it to develop in different adverse environments. The Salar de Huasco is one of these environments due to its altitude, atmospheric pressure, solar radiation, temperature variations, pH, salinity, and the presence of toxic compounds such as arsenic. However, the physiological and/or molecular mechanisms that enable them to prosper in these environments have not yet been described. Our research group has isolated several strains of Exiguobacterium genus from different sites of Salar de Huasco, which show different resistance levels to As(III) and As(V). In this work, we compare the protein expression patterns of the three strains in response to arsenic by a proteomic approach; strains were grown in absence of the metalloid and in presence of As(III) and As(V) sublethal concentrations and the protein separation was carried out in 2D electrophoresis gels (2D-GE). In total, 999 spots were detected, between 77 and 173 of which showed significant changes for As(III) among the three strains, and between 90 and 143 for As(V), respectively, compared to the corresponding control condition. Twenty-seven of those were identified by mass spectrometry (MS). Among these identified proteins, the ArsA [ATPase from the As(III) efflux pump] was found to be up-regulated in response to both arsenic conditions in the three strains, as well as the Co-enzyme A disulfide reductase (Cdr) in the two more resistant strains. Interestingly, in this genus the gene that codifies for Cdr is found within the genic context of the ars operon. We suggest that this protein could be restoring antioxidants molecules, necessary for the As(V) reduction. Additionally, among the proteins that change their expression against As, we found several with functions relevant to stress response, e.g., Hpf, LuxS, GLpX, GlnE, and Fur. This study allowed us to shed light into the physiology necessary for these bacteria to be able to tolerate the toxicity and stress generated by the presence of arsenic in their niche.

Functional prediction, characterization, and categorization of operome from Acetoanaerobium sticklandii DSM 519

Acetoanaerobium sticklandii DSM 519 is a hyper-ammonia producing anaerobic bacterium that can be able utilizes amino acids as sole carbon and energy sources for its growth and energetic metabolism. A lack of knowledge on its molecular machinery and 30.5% conserved hypothetical proteins (HPs; operome) hinders the successful utility in biofuel applications. In this study, we have predicted, characterized and categorized its operome whose functions are still not determined accurately using a combined bioinformatics approach. The functions of 64 of the 359 predicted HPs are involved in diverse metabolic subsystems. A. sticklandii operome has consisted of 16% Rossmann fold and 46% miscellaneous folds. Subsystems-based technology has classified 51 HPs contributing to the small-molecular reactions, 26 in macromolecular reactions and 12 in the biosynthesis of cofactors, prosthetic groups and electron carriers. A generality of functions predicted from its operome contributed to the cell cycle, amino acid metabolism, membrane transport, and regulatory processes. Many of them have duplicated functions as paralogs in this genome. A. sticklandii has the ability to compete with invading microorganisms and tolerate abiotic stresses, which can be overwhelmed by the predicted functions of its operome. Results of this study revealed that it has specialized systems for amino acid catabolism-directed solventogenesis and acidogenesis but the level of gene expression may determine the metabolic function in amino acid fermenting niches in the rumina of cattle. As shown by our analysis, the predicted functions of its operome allow us for a better understanding of the growth and physiology at systems-scale.

Identifying miltefosine-resistant key genes in protein-protein interactions network and experimental verification in Iranian Leishmania major

Drug resistance is a complex phenomenon during leishmaniasis chemotherapy. In this study, the genes and pathways involved in miltefosine (MIL)-resistant Leishmania were identified using microarray data and in silico approaches. GSE30685 and GSE45496 were obtained from GEO database and analyzed with GEO2R tool to identify genes involved in MIL-resistant Leishmania. 177 differentially expressed genes (DEGs) were selected from these GSEs, which about half of them were uncharacterized/hypothetical proteins. The interactions between DEGs were investigated using STRING database and protein-protein interaction (PPI) networks. Five hub nodes were found in the PPI network. The gene ontology (GO) analysis of the resulting network revealed that DNA replication (GO:0006260) and ATP hydrolysis coupled proton transport (GO:0015991) were the most enriched GO term. Iranian MIL-resistant Leishmania major (L. major) parasites were generated by exposure of wild-type isolates to the increasing concentrations of MIL over a period of 5 months. Proof of mRNA expression levels of the obtained hub genes was assessed in Iranian wild-type and acquired resistant L. major parasites by real-time PCR. A significant higher expression level of LDBPK_150170 (encoding protein phosphatase 2C, PP2C), was only observed in Iranian L. major parasites resistance to MIL. Moreover, the RT-PCR results showed that the expression of metacyclic marker (small hydrophilic endoplasmic reticulum-associated protein, SHERP) and MIL-resistant marker (Leishmania MIL-transporter, LMT) was significantly increased and decreased, respectively, in Iranian MIL-resistant L. major parasites. Taken together, these data suggested that PP2C as well as SHERP and LMT genes may be prospective targets for the treatment of MIL-resistant Leishmania.