Search

The transcriptional regulatory network (TRN) of Pseudomonas aeruginosa plays a critical role in coordinating numerous cellular processes. We extracted and quality controlled all publicly available RNA-sequencing datasets for P. aeruginosa to find 281 high-quality transcriptomes. We produced 83 new RNAseq data sets under critical conditions to generate a comprehensive compendium of 364 transcriptomes. We used this compendium to reconstruct the TRN of P. aeruginosa using independent component analysis (ICA). We identified 104 independently modulated sets of genes (called iModulons), among which 81 (78%) reflect the effects of known transcriptional regulators. We show that iModulons: 1) play an important role in defining the genomic boundaries of biosynthetic gene clusters (BGCs); 2) show increased expression of the BGCs and associated secretion systems in conditions that emulate cystic fibrosis (CF); 3) show the presence of a novel BGC named RiPP (bacteriocin producer) which might have a role in worsening CF outcomes; 4) exhibit the interplay of amino acid metabolism regulation and central metabolism across carbon sources, and 5) clustered according to their activity changes to define iron and sulfur stimulons. Finally, we compare the iModulons of P. aeruginosa with those of E. coli to observe conserved regulons across two gram negative species. This comprehensive TRN framework covers almost every aspect of the transcriptional regulatory machinery in P. aeruginosa, and thus could prove foundational for future research of its physiological functions.

The ESKAPEE pathogens are the leading cause of health care-associated infections worldwide. Two-component systems (TCSs) can be used as effective targets against pathogenic bacteria since they are ubiquitous and manage various vital functions such as antibiotic resistance, virulence, biofilms, quorum sensing, and pH balance, among others., The two-component system (TCS) helps bacteria sense and respond to environmental stimuli through histidine kinases and response regulators. TCSs are the largest family of multistep signal transduction processes, and they are involved in many important cellular processes such as antibiotic resistance, pathogenicity, quorum sensing, osmotic stress, and biofilms. Here, we perform the first comprehensive study to highlight the role of TCSs as potential drug targets against ESKAPEE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter spp., and Escherichia coli) pathogens through annotation, mapping, pangenomic status, gene orientation, and sequence variation analysis. The distribution of the TCSs is group specific with regard to Gram-positive and Gram-negative bacteria, except for KdpDE. The TCSs among ESKAPEE pathogens form closed pangenomes, except for Pseudomonas aeruginosa. Furthermore, their conserved nature due to closed pangenomes might make them good drug targets. Fitness score analysis suggests that any mutation in some TCSs such as BaeSR, ArcBA, EvgSA, and AtoSC, etc., might be lethal to the cell. Taken together, the results of this pangenomic assessment of TCSs reveal a range of strategies deployed by the ESKAPEE pathogens to manifest pathogenicity and antibiotic resistance. This study further suggests that the conserved features of TCSs might make them an attractive group of potential targets with which to address antibiotic resistance. IMPORTANCE The ESKAPEE pathogens are the leading cause of health care-associated infections worldwide. Two-component systems (TCSs) can be used as effective targets against pathogenic bacteria since they are ubiquitous and manage various vital functions such as antibiotic resistance, virulence, biofilms, quorum sensing, and pH balance, among others. This study provides a comprehensive overview of the pangenomic status of the TCSs among ESKAPEE pathogens. The annotation and pangenomic analysis of TCSs show that they are significantly distributed and conserved among the pathogens, as most of them form closed pangenomes. Furthermore, our analysis also reveals that the removal of the TCSs significantly affects the fitness of the cell. Hence, they may be used as promising drug targets against bacteria.

A number of anti-retroviral drugs are being used for treating Human Immunodeficiency Virus (HIV) infection. Due to emergence of drug resistant strains, there is a constant quest to discover more effective anti-HIV compounds. In this endeavor, computational tools have proven useful in accelerating drug discovery. Although methods were published to design a class of compounds against a specific HIV protein, but an integrated web server for the same is lacking. Therefore, we have developed support vector machine based regression models using experimentally validated data from ChEMBL repository. Quantitative structure activity relationship based features were selected for predicting inhibition activity of a compound against HIV proteins namely protease (PR), reverse transcriptase (RT) and integrase (IN). The models presented a maximum Pearson correlation coefficient of 0.78, 0.76, 0.74 and 0.76, 0.68, 0.72 during tenfold cross-validation on IC50 and percent inhibition datasets of PR, RT, IN respectively. These models performed equally well on the independent datasets. Chemical space mapping, applicability domain analyses and other statistical tests further support robustness of the predictive models. Currently, we have identified a number of chemical descriptors that are imperative in predicting the compound inhibition potential. HIVprotI platform (http://bioinfo.imtech.res.in/manojk/hivproti) would be useful in virtual screening of inhibitors as well as designing of new molecules against the important HIV proteins for therapeutics development. Electronic supplementary material The online version of this article (10.1186/s13321-018-0266-y) contains supplementary material, which is available to authorized users.

Viral epidemics and pandemics are considered public health emergencies. However, traditional and novel antiviral discovery approaches are unable to mitigate them in a timely manner. Notably, drug repurposing emerged as an alternative strategy to provide antiviral solutions in a timely and cost-effective manner. In the literature, many FDA-approved drugs have been repurposed to inhibit viruses, while a few among them have also entered clinical trials. Using experimental data, we identified repurposed drugs against 14 viruses responsible for causing epidemics and pandemics such as SARS-CoV-2, SARS, Middle East respiratory syndrome, influenza H1N1, Ebola, Zika, Nipah, chikungunya, and others. We developed a novel computational "drug-target-drug" approach that uses the drug-targets extracted for specific drugs, which are experimentally validated in vitro or in vivo for antiviral activity. Furthermore, these extracted drug-targets were used to fetch the novel FDA-approved drugs for each virus and prioritize them by calculating their confidence scores. Pathway analysis showed that the majority of the extracted targets are involved in cancer and signaling pathways. For SARS-CoV-2, our method identified 21 potential repurposed drugs, of which 7 (e.g., baricitinib, ramipril, chlorpromazine, enalaprilat, etc.) have already entered clinical trials. The prioritized drug candidates were further validated using a molecular docking approach. Therefore, we anticipate success during the experimental validation of our predicted FDA-approved repurposed drugs against 14 viruses. This study will assist the scientific community in hastening research aimed at the development of antiviral therapeutics.

Flaviviruses are arboviruses, which comprised of more than 70 viruses, covering broad geographic ranges, and responsible for significant mortality and morbidity globally. Due to the lack of efficient inhibitors targeting flaviviruses, the designing of novel and efficient anti-flavi agents is an important problem. Therefore, in the current study we have developed a dedicated prediction algorithm anti-flavi, to identify inhibition ability of chemicals and peptides against flaviviruses through Quantitative structure–activity relationship (QSAR) based method. We extracted the non-redundant 2168 chemicals and 117 peptides from ChEMBL and AVPpred databases respectively with reported IC50 values. The regression based model developed on training/testing datasets of 1952 chemicals and 105 peptides displayed the Pearson’s correlation coefficient (PCC) of 0.87, 0.84 and 0.87, 0.83 using Support Vector Machine and Random Forest techniques correspondingly. We also explored the peptidomimetics approach, in which the most contributing descriptors of peptides were used to identify chemicals having anti-flavi potential. Conversely, the selected descriptors of chemicals performed well to predict anti-flavi peptides. Moreover, the developed model proved to be highly robust while checked through various approaches like independent validation and decoy datasets. We hope that our web server would prove a useful tool to predict and design the efficient anti-flavi agents. The anti-flavi webserver is freely available at URL http://bioinfo.imtech.res.in/manojk/antiflavi.

The quorum sensing (QS) is a well-characterized phenomenon in the microbial world for cellular communication and being exploited for intraspecies, interspecies and even interkingdom interactions. This robust and widely distributed social networking cascade assists the microbes to emerge as strong parasites for the hosts. Moreover, the colonization with more efficient cross talk among the microbes further intensifies their infections. Here, we will be focusing to decipher the evolutionary status of the QS regulators (LuxI and LuxR) in the prokaryotic world. LuxI is a signal synthase, while LuxR is the recipient for sensing internal (cognate LuxR/solos LuxR) and external (solos LuxR) signals. These regulators are reported to evolve vertically as well as borrowed through horizontal gene transfer w.r.t. ecological niche. Their universal distribution in the microbial world further corroborates the need for targeting multiple signaling system regulators.

Quorum sensing is a widespread phenomenon in prokaryotes that helps them to communicate among themselves and with eukaryotes. It is driven through quorum sensing signaling molecules (QSSMs) in a density dependent manner that assists in numerous biological functions like biofilm formation, virulence factors secretion, swarming motility, bioluminescence, etc. Despite immense implications, dedicated resources of QSSMs are lacking. Therefore, we have developed SigMol (http://bioinfo.imtech.res.in/manojk/sigmol), a specialized repository of these molecules in prokaryotes. SigMol harbors information on QSSMs pertaining to different quorum sensing signaling systems namely acylated homoserine lactones (AHLs), diketopiperazines (DKPs), 4-hydroxy-2-alkylquinolines (HAQs), diffusible signal factors (DSFs), autoinducer-2 (AI-2) and others. Database contains 1382: entries of 182: unique signaling molecules from 215: organisms. It encompasses biological as well as chemical aspects of signaling molecules. Biological information includes genes, preliminary bioassays, identification assays and applications, while chemical detail comprises of IUPAC name, SMILES and structure. We have provided user-friendly browsing and searching facilities for easy data retrieval and comparison. We have gleaned information of diverse QSSMs reported in literature at a single platform 'SigMol'. This comprehensive resource will assist the scientific community in understanding intraspecies, interspecies or interkingdom networking and further help to unfold different facets of quorum sensing and related therapeutics.

Biofilms play an important role in the antibiotic drug resistance, which is threatening public health globally. Almost, all microbes mimic multicellular lifestyle to form biofilm by undergoing phenotypic changes to adapt adverse environmental conditions. Many anti-biofilm agents have been experimentally validated to disrupt the biofilms during last three decades. To organize this data, we developed the ‘aBiofilm’ resource (http://bioinfo.imtech.res.in/manojk/abiofilm/) that harbors a database, a predictor, and the data visualization modules. The database contains biological, chemical, and structural details of 5027 anti-biofilm agents (1720 unique) reported from 1988–2017. These agents target over 140 organisms including Gram-negative, Gram-positive bacteria, and fungus. They are mainly chemicals, peptides, phages, secondary metabolites, antibodies, nanoparticles and extracts. They show the diverse mode of actions by attacking mainly signaling molecules, biofilm matrix, genes, extracellular polymeric substances, and many more. The QSAR based predictor identifies the anti-biofilm potential of an unknown chemical with an accuracy of ∼80.00%. The data visualization section summarized the biofilm stages targeted (Circos plot); interaction maps (Cytoscape) and chemicals diversification (CheS-Mapper) of the agents. This comprehensive platform would help the researchers to understand the multilevel communication in the microbial consortium. It may aid in developing anti-biofilm therapeutics to deal with antibiotic drug resistance menace.

Genome editing by sgRNA a component of CRISPR/Cas system emerged as a preferred technology for genome editing in recent years. However, activity and stability of sgRNA in genome targeting is greatly influenced by its sequence features. In this endeavor, a few prediction tools have been developed to design effective sgRNAs but these methods have their own limitations. Therefore, we have developed “ge-CRISPR” using high throughput data for the prediction and analysis of sgRNAs genome editing efficiency. Predictive models were employed using SVM for developing pipeline-1 (classification) and pipeline-2 (regression) using 2090 and 4139 experimentally verified sgRNAs respectively from Homo sapiens, Mus musculus, Danio rerio and Xenopus tropicalis. During 10-fold cross validation we have achieved accuracy and Matthew’s correlation coefficient of 87.70% and 0.75 for pipeline-1 on training dataset (T1840) while it performed equally well on independent dataset (V250). In pipeline-2 we attained Pearson correlation coefficient of 0.68 and 0.69 using best models on training (T3169) and independent dataset (V520) correspondingly. ge-CRISPR (http://bioinfo.imtech.res.in/manojk/gecrispr/) for a given genomic region will identify potent sgRNAs, their qualitative as well as quantitative efficiencies along with potential off-targets. It will be useful to scientific community engaged in CRISPR research and therapeutics development.

Knowledge of the subcellular location (SCL) of viral proteins in the host cell is important for understanding their function in depth. Therefore, we have developed "MSLVP", a two-tier prediction algorithm for predicting multiple SCLs of viral proteins. For this study, data sets of comprehensive viral proteins with experimentally validated SCL annotation were collected from UniProt. Non-redundant (90%) data sets of 3480 viral proteins that belonged to single (2715), double (391) and multiple (374) sites were employed. Additionally, 1687 (30% sequence identity) viral proteins were categorised into single (1366), double (167) and multiple (154) sites. Single, double and multiple locations further comprised of eight, four and six categories, respectively. Viral protein locations include the nucleus, cytoplasm, endoplasmic reticulum, extracellular, single-pass membrane, multi-pass membrane, capsid, remaining others and combinations thereof. Support vector machine based models were developed using sequence features like amino acid composition, dipeptide composition, physicochemical properties and their hybrids. We have employed "one-versus-one" as well as "one-versus-other" strategies for multiclass classification. The performance of "one-versus-one" is better than the "one-versus-other" approach during 10-fold cross-validation. For the 90% data set, we achieved an accuracy, a Matthew's correlation coefficient (MCC) and a receiver operating characteristic (ROC) of 99.99%, 1.00, 1.00; 100.00%, 1.00, 1.00 and 99.90%; 1.00, 1.00 for single, double and multiple locations, respectively. Similar results were achieved for a 30% sequence identity data set. Predictive models for each SCL performed equally well on the independent dataset. The MSLVP web server () can predict subcellular locations i.e. single (8; including single and multi-pass membrane), double (4) and multiple (6). This would be helpful for elucidating the functional annotation of viral proteins and potential drug targets.