Centrality in the host-pathogen interactome is associated with pathogen fitness during infection

High throughput methods to study protein-protein interactions during host-pathogen interactions

The ability of a pathogen to survive and cause an infection is often determined by specific interactions between the host and pathogen proteins. Such interactions can be both intra- and extracellular and may define the outcome of an infection. There are a range of innovative biochemical, biophysical and bioinformatic techniques currently available to identify protein-protein interactions (PPI) between the host and the pathogen. However, the complexity and the diversity of host-pathogen PPIs has led to the development of several high throughput (HT) techniques that enable the study of multiple interactions at once and/or screen multiple samples at the same time, in an unbiased manner. We review here the major HT laboratory-based technologies employed for host-bacterial interaction studies.

Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications

SARS-CoV-2 infections, commonly referred to as COVID-19, remain a critical risk to both human life and global economies. Particularly, COVID-19 patients with weak immunity may suffer from different complications due to the bacterial co-infections/super-infections/secondary infections. Therefore, different variants of alternative antibacterial therapeutic agents are required to inhibit those infection-causing drug-resistant pathogenic bacteria. This study attempted to explore these bacterial pathogens and their inhibitors by using integrated statistical and bioinformatics approaches. By analyzing bacterial 16S rRNA sequence profiles, at first, we detected five bacterial genera and taxa (Bacteroides, Parabacteroides, Prevotella Clostridium, Atopobium, and Peptostreptococcus) based on differentially abundant bacteria between SARS-CoV-2 infection and control samples that are significantly enriched in 23 metabolic pathways. A total of 183 bacterial genes were found in the enriched pathways. Then, the top-ranked 10 bacterial genes (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) were selected as the pathogenic bacterial key genes (bKGs) by their protein-protein interaction (PPI) network analysis. Then, we detected bKG-guided top-ranked eight drug molecules (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Entreatinib, Digitoxin, and Elbasvir) by molecular docking. Finally, the binding stability of the top-ranked three drug molecules (Bemcentinib, Ledipasvir, and Velpatasvir) against three receptors (hldD, mlaA, and lptD) was investigated by computing their binding free energies with molecular dynamic (MD) simulation-based MM-PBSA techniques, respectively, and was found to be stable. Therefore, the findings of this study could be useful resources for developing a proper treatment plan against bacterial co-/super-/secondary-infection in SARS-CoV-2 infections.

Identification of host proteins interacting with the E protein of porcine epidemic diarrhea virus

Introduction

Porcine epidemic diarrhea (PED) is an acute, highly contagious, and high-mortality enterophilic infectious disease caused by the porcine epidemic diarrhea virus (PEDV). PEDV is globally endemic and causes substantial economic losses in the swine industry. The PEDV E protein is the smallest structural protein with high expression levels that interacts with the M protein and participates in virus assembly. However, how the host proteins interact with E proteins in PEDV replication remains unknown.

Methods

We identified host proteins that interact with the PEDV E protein using a combination of PEDV E protein-labeled antibody co-immunoprecipitation and tandem liquid-chromatography mass-spectroscopy (LC-MS/MS).

Results

Bioinformatical analysis showed that in eukaryotes, ribosome biogenesis, RNA transport, and amino acid biosynthesis represent the three main pathways that are associated with the E protein. The interaction between the E protein and isocitrate dehydrogenase [NAD] β-subunit (NAD-IDH-β), DNA-directed RNA polymerase II subunit RPB9, and mRNA-associated protein MRNP 41 was validated using co-immunoprecipitation and confocal assays. NAD-IDH-β overexpression significantly inhibited viral replication.

Discussion

The antiviral effect of NAD-IDH-β suggesting that the E protein may regulate host metabolism by interacting with NAD-IDH-β, thereby reducing the available energy for viral replication. Elucidating the interaction between the PEDV E protein and host proteins may clarify its role in viral replication. These results provide a theoretical basis for the study of PEDV infection mechanism and antiviral targets.

Structural assembly of the bacterial essential interactome

The study of protein interactions in living organisms is fundamental for understanding biological processes and central metabolic pathways. Yet, our knowledge of the bacterial interactome remains limited. Here, we combined gene deletion mutant analysis with deep-learning protein folding using AlphaFold2 to predict the core bacterial essential interactome. We predicted and modeled 1402 interactions between essential proteins in bacteria and generated 146 high-accuracy models. Our analysis reveals previously unknown details about the assembly mechanisms of these complexes, highlighting the importance of specific structural features in their stability and function. Our work provides a framework for predicting the essential interactomes of bacteria and highlight the potential of deep-learning algorithms in advancing our understanding of the complex biology of living organisms. Also, the results presented here offer a promising approach to identify novel antibiotic targets.

DomiRank Centrality reveals structural fragility of complex networks via node dominance

Determining the key elements of interconnected infrastructure and complex systems is paramount to ensure system functionality and integrity. This work quantifies the dominance of the networks' nodes in their respective neighborhoods, introducing a centrality metric, DomiRank, that integrates local and global topological information via a tunable parameter. We present an analytical formula and an efficient parallelizable algorithm for DomiRank centrality, making it applicable to massive networks. From the networks' structure and function perspective, nodes with high values of DomiRank highlight fragile neighborhoods whose integrity and functionality are highly dependent on those dominant nodes. Underscoring this relation between dominance and fragility, we show that DomiRank systematically outperforms other centrality metrics in generating targeted attacks that effectively compromise network structure and disrupt its functionality for synthetic and real-world topologies. Moreover, we show that DomiRank-based attacks inflict more enduring damage in the network, hindering its ability to rebound and, thus, impairing system resilience. DomiRank centrality capitalizes on the competition mechanism embedded in its definition to expose the fragility of networks, paving the way to design strategies to mitigate vulnerability and enhance the resilience of critical infrastructures.

The protein interactome of the citrus Huanglongbing pathogen Candidatus Liberibacter asiaticus

The bacterium Candidatus Liberibacter asiaticus (CLas) causes citrus Huanglongbing disease. Our understanding of the pathogenicity and biology of this microorganism remains limited because CLas has not yet been cultivated in artificial media. Its genome is relatively small and encodes approximately 1136 proteins, of which 415 have unknown functions. Here, we use a high-throughput yeast-two-hybrid (Y2H) screen to identify interactions between CLas proteins, thus providing insights into their potential functions. We identify 4245 interactions between 542 proteins, after screening 916 bait and 936 prey proteins. The false positive rate of the Y2H assay is estimated to be 2.9%. Pull-down assays for nine protein-protein interactions (PPIs) likely involved in flagellar function support the robustness of the Y2H results. The average number of PPIs per node in the CLas interactome is 15.6, which is higher than the numbers previously reported for interactomes of free-living bacteria, suggesting that CLas genome reduction has been accompanied by increased protein multi-functionality. We propose potential functions for 171 uncharacterized proteins, based on the PPI results, guilt-by-association analyses, and comparison with data from other bacterial species. We identify 40 hub-node proteins, including quinone oxidoreductase and LysR, which are known to protect other bacteria against oxidative stress and might be important for CLas survival in the phloem. We expect our PPI database to facilitate research on CLas biology and pathogenicity mechanisms.

Role of molecular mimicry in the SARS-CoV-2-human interactome for pathogenesis of cardiovascular diseases: An update to ImitateDB

Mimicry of host proteins is a strategy employed by pathogens to hijack host functions. Domain and motif mimicry was explored in the experimental and predicted SARS-CoV-2-human interactome. The host first interactor proteins were also added to capture the continuum of the interactions. The domains and motifs of the proteins were annotated using NCBI CD Search and ScanProsite, respectively. Host and pathogen proteins with a common host interactor and similar domain/motif constitute a mimicry pair indicating global structural similarity (domain mimicry pair; DMP) or local sequence similarity (motif mimicry pair; MMP). 593 DMPs and 7,02,472 MMPs were determined. AAA, DEXDc and Macro domains were frequent among DMPs whereas glycosylation, myristoylation and RGD motifs were abundant among MMP. The proteins involved in mimicry were visualised as a SARS-CoV-2 mimicry interaction network. The host proteins were enriched in multiple CVD pathways indicating the role of mimicry in COVID-19 associated CVDs. Bridging nodes were identified as potential drug targets. Approved antihypertensive and anti-inflammatory drugs are proposed for repurposing against COVID-19 associated CVDs. The SARS-CoV-2 mimicry data has been updated in ImitateDB (http://imitatedb.sblab-nsit.net/SARSCoV2Mimicry). Determination of key mechanisms, proteins, pathways, drug targets and repurposing candidates is critical for developing therapeutics for SARS CoV-2 associated CVDs.

A binary interaction map between turnip mosaic virus and Arabidopsis thaliana proteomes

Viruses are obligate intracellular parasites that have co-evolved with their hosts to establish an intricate network of protein-protein interactions. Here, we followed a high-throughput yeast two-hybrid screening to identify 378 novel protein-protein interactions between turnip mosaic virus (TuMV) and its natural host Arabidopsis thaliana. We identified the RNA-dependent RNA polymerase NIb as the viral protein with the largest number of contacts, including key salicylic acid-dependent transcription regulators. We verified a subset of 25 interactions in planta by bimolecular fluorescence complementation assays. We then constructed and analyzed a network comprising 399 TuMV-A. thaliana interactions together with intravirus and intrahost connections. In particular, we found that the host proteins targeted by TuMV are enriched in different aspects of plant responses to infections, are more connected and have an increased capacity to spread information throughout the cell proteome, display higher expression levels, and have been subject to stronger purifying selection than expected by chance. The proviral or antiviral role of ten host proteins was validated by characterizing the infection dynamics in the corresponding mutant plants, supporting a proviral role for the transcriptional regulator TGA1. Comparison with similar studies with animal viruses, highlights shared fundamental features in their mode of action.

HPIPred: Host-pathogen interactome prediction with phenotypic scoring

Protein-protein interactions (PPIs) are involved in most cellular processes. Unfortunately, current knowledge of host-pathogen interactomes is still very limited. Experimental methods used to detect PPIs have several limitations, including increasing complexity and economic cost in large-scale screenings. Hence, computational methods are commonly used to support experimental data, although they generally suffer from high false-positive rates. To address this issue, we have created HPIPred, a host-pathogen PPI prediction tool based on numerical encoding of physicochemical properties. Unlike other available methods, HPIPred integrates phenotypic data to prioritize biologically meaningful results. We used HPIPred to screen the entire Homo sapiens and Pseudomonas aeruginosa PAO1 proteomes to generate a host-pathogen interactome with 763 interactions displaying a highly connected network topology. Our predictive model can be used to prioritize protein-protein interactions as potential targets for antibacterial drug development. Available at: https://github.com/SysBioUAB/hpi_predictor.

Fungal Effectoromics: A World in Constant Evolution

Effectors are small, secreted molecules that mediate the establishment of interactions in nature. While some concepts of effector biology have stood the test of time, this area of study is ever-evolving as new effectors and associated characteristics are being revealed. In the present review, the different characteristics that underly effector classifications are discussed, contrasting past and present knowledge regarding these molecules to foster a more comprehensive understanding of effectors for the reader. Research gaps in effector identification and perspectives for effector application in plant disease management are also presented, with a focus on fungal effectors in the plant-microbe interaction and interactions beyond the plant host. In summary, the review provides an amenable yet thorough introduction to fungal effector biology, presenting noteworthy examples of effectors and effector studies that have shaped our present understanding of the field.

The Network of Interactions between the Porcine Epidemic Diarrhea Virus Nucleocapsid and Host Cellular Proteins

Host–virus protein interactions are critical for intracellular viral propagation. Understanding the interactions between cellular and viral proteins may help us develop new antiviral strategies. Porcine epidemic diarrhea virus (PEDV) is a highly contagious coronavirus that causes severe damage to the global swine industry. Here, we employed co-immunoprecipitation and liquid chromatography-mass spectrometry to characterize 426 unique PEDV nucleocapsid (N) protein-binding proteins in infected Vero cells. A protein–protein interaction network (PPI) was created, and gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) database analyses revealed that the PEDV N-bound proteins belong to different cellular pathways, such as nucleic acid binding, ribonucleoprotein complex binding, RNA methyltransferase, and polymerase activities. Interactions of the PEDV N protein with 11 putative proteins: tripartite motif containing 21, DEAD-box RNA helicase 24, G3BP stress granule assembly factor 1, heat shock protein family A member 8, heat shock protein 90 alpha family class B member 1, YTH domain containing 1, nucleolin, Y-box binding protein 1, vimentin, heterogeneous nuclear ribonucleoprotein A2/B1, and karyopherin subunit alpha 1, were further confirmed by in vitro co-immunoprecipitation assay. In summary, studying an interaction network can facilitate the identification of antiviral therapeutic strategies and novel targets for PEDV infection.

Analysis of Host–Bacteria Protein Interactions Reveals Conserved Domains and Motifs That Mediate Fundamental Infection Pathways

Adhesion and colonization of host cells by pathogenic bacteria depend on protein–protein interactions (PPIs). These interactions are interesting from the pharmacological point of view since new molecules that inhibit host-pathogen PPIs would act as new antimicrobials. Most of these interactions are discovered using high-throughput methods that may display a high false positive rate. The absence of curation of these databases can make the available data unreliable. To address this issue, a comprehensive filtering process was developed to obtain a reliable list of domains and motifs that participate in PPIs between bacteria and human cells. From a structural point of view, our analysis revealed that human proteins involved in the interactions are rich in alpha helix and disordered regions and poorer in beta structure. Disordered regions in human proteins harbor short sequence motifs that are specifically recognized by certain domains in pathogenic proteins. The most relevant domain–domain interactions were validated by AlphaFold, showing that a proper analysis of host-pathogen PPI databases can reveal structural conserved patterns. Domain–motif interactions, on the contrary, were more difficult to validate, since unstructured regions were involved, where AlphaFold could not make a good prediction. Moreover, these interactions are also likely accommodated by post-translational modifications, especially phosphorylation, which can potentially occur in 25–50% of host proteins. Hence, while common structural patterns are involved in host–pathogen PPIs and can be retrieved from available databases, more information is required to properly infer the full interactome. By resolving these issues, and in combination with new prediction tools like Alphafold, new classes of antimicrobials could be discovered from a more detailed understanding of these interactions.

dRNASb: a systems biology approach to decipher dynamics of host-pathogen interactions using temporal dual RNA-seq data

Infection triggers a dynamic cascade of reciprocal events between host and pathogen wherein the host activates complex mechanisms to recognise and kill pathogens while the pathogen often adjusts its virulence and fitness to avoid eradication by the host. The interaction between the pathogen and the host results in large-scale changes in gene expression in both organisms. Dual RNA-seq, the simultaneous detection of host and pathogen transcripts, has become a leading approach to unravelling complex molecular interactions between the host and the pathogen and is particularly informative for intracellular organisms. The amount of in vitro and in vivo dual RNA-seq data is rapidly growing, which demands computational pipelines to effectively analyse such data. In particular, holistic, systems-level, and temporal analyses of dual RNA-seq data are essential to enable further insights into the host-pathogen transcriptional dynamics and potential interactions. Here, we developed an integrative network-driven bioinformatics pipeline, dRNASb, a systems biology-based computational pipeline to analyse temporal transcriptional clusters, incorporate molecular interaction networks (e.g. protein-protein interactions), identify topologically and functionally key transcripts in host and pathogen, and associate host and pathogen temporal transcriptome to decipher potential between-species interactions. The pipeline is applicable to various dual RNA-seq data from different species and experimental conditions. As a case study, we applied dRNASb to analyse temporal dual RNA-seq data of Salmonella-infected human cells, which enabled us to uncover genes contributing to the infection process and their potential functions and to identify putative associations between host and pathogen genes during infection. Overall, dRNASb has the potential to identify key genes involved in bacterial growth or host defence mechanisms for future uses as therapeutic targets.

Transcriptomic profiling of Escherichia coli K-12 in response to a compendium of stressors

Environmental perturbations impact multiple cellular traits, including gene expression. Bacteria respond to these stressful situations through complex gene interaction networks, thereby inducing stress tolerance and survival of cells. In this paper, we study the response mechanisms of E. coli when exposed to different environmental stressors via differential expression and co-expression analysis. Gene co-expression networks were generated and analyzed via Weighted Gene Co-expression Network Analysis (WGCNA). Based on the gene co-expression networks, genes with similar expression profiles were clustered into modules. The modules were analysed for identification of hub genes, enrichment of biological processes and transcription factors. In addition, we also studied the link between transcription factors and their differentially regulated targets to understand the regulatory mechanisms involved. These networks validate known gene interactions and provide new insights into genes mediating transcriptional regulation in specific stress environments, thus allowing for in silico hypothesis generation.

ImitateDB: A database for domain and motif mimicry incorporating host and pathogen protein interactions

Molecular mimicry of host proteins by pathogens constitutes a strategy to hijack the host pathways. At present, there is no dedicated resource for mimicked domains and motifs in the host-pathogen interactome. In this work, the experimental host-pathogen (HP) and host-host (HH) protein-protein interactions (PPIs) were collated. The domains and motifs of these proteins were annotated using CD Search and ScanProsite, respectively. Host and pathogen proteins with a shared host interactor and similar domain/motif constitute a mimicry pair exhibiting global structural similarity (domain mimicry pair; DMP) or local sequence motif similarity (motif mimicry pair; MMP). Mimicry pairs are likely to be co-expressed and co-localized. 1,97,607 DMPs and 32,67,568 MMPs were identified in 49,265 experimental HP-PPIs and organized in a web-based resource, ImitateDB ( http://imitatedb.sblab-nsit.net ) that can be easily queried. The results are externally integrated using hyperlinked domain PSSM ID, motif ID, protein ID and PubMed ID. Kinase, UL36, Smc and DEXDc were frequent DMP domains whereas protein kinase C phosphorylation, casein kinase 2 phosphorylation, glycosylation and myristoylation sites were frequent MMP motifs. Novel DMP domains SANT, Tudor, PhoX and MMP motif microbody C-terminal targeting signal, cornichon signature and lipocalin signature were proposed. ImitateDB is a novel resource for identifying mimicry in interacting host and pathogen proteins.

Interaction Network of Porcine Circovirus Type 3 and 4 Capsids with Host Proteins

An extensive understanding of the interactions between host cellular and viral proteins provides clues for studying novel antiviral strategies. Porcine circovirus type 3 (PCV3) and type 4 (PCV4) have recently been identified as viruses that can potentially damage the swine industry. Herein, 401 putative PCV3 Cap-binding and 484 putative PCV4 Cap-binding proteins were characterized using co-immunoprecipitation and liquid chromatography-mass spectrometry. Both PCV3 and PCV4 Caps shared 278 identical interacting proteins, but some putative interacting proteins (123 for PCV3 Cap and 206 for PCV4 Cap) differed. A protein-protein interaction network was constructed, and according to gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) database analyses, both PCV3 Cap- and PCV4 Cap-binding proteins participated mainly in ribosome biogenesis, nucleic acid binding, and ATP-dependent RNA helicase activities. Verification assays of eight putative interacting proteins indicated that nucleophosmin-1, nucleolin, DEAD-box RNA helicase 21, heterogeneous nuclear ribonucleoprotein A2/B1, YTH N6-methyladenosine RNA binding protein 1, and Y-box binding protein 1 bound directly to both PCV3 and PCV4 Caps, but ring finger protein 2 and signal transducer and activator of transcription 6 did not. Therefore, the interaction network provided helpful information to support further research into the underlying mechanisms of PCV3 and PCV4 infection.

Graphlet eigencentralities capture novel central roles of genes in pathways

MOTIVATION

Graphlet adjacency extends regular node adjacency in a network by considering a pair of nodes being adjacent if they participate in a given graphlet (small, connected, induced subgraph). Graphlet adjacencies captured by different graphlets were shown to contain complementary biological functions and cancer mechanisms. To further investigate the relationships between the topological features of genes participating in molecular networks, as captured by graphlet adjacencies, and their biological functions, we build more descriptive pathway-based approaches.

CONTRIBUTION

We introduce a new graphlet-based definition of eigencentrality of genes in a pathway, graphlet eigencentrality, to identify pathways and cancer mechanisms described by a given graphlet adjacency. We compute the centrality of genes in a pathway either from the local perspective of the pathway or from the global perspective of the entire network.

RESULTS

We show that in molecular networks of human and yeast, different local graphlet adjacencies describe different pathways (i.e., all the genes that are functionally important in a pathway are also considered topologically important by their local graphlet eigencentrality). Pathways described by the same graphlet adjacency are functionally similar, suggesting that each graphlet adjacency captures different pathway topology and function relationships. Additionally, we show that different graphlet eigencentralities describe different cancer driver genes that play central roles in pathways, or in the crosstalk between them (i.e. we can predict cancer driver genes participating in a pathway by their local or global graphlet eigencentrality). This result suggests that by considering different graphlet eigencentralities, we can capture different functional roles of genes in and between pathways.

A determination of pan-pathogen antimicrobials?

While antimicrobial drug development has historically mitigated infectious diseases that are known, COVID-19 revealed a dearth of ‘in-advance’ therapeutics suitable for infections by pathogens that have not yet emerged. Such drugs must exhibit a property that is antithetical to the classical paradigm of antimicrobial development: the ability to treat infections by any pathogen. Characterisation of such ‘pan-pathogen’ antimicrobials requires consolidation of drug repositioning studies, a new and growing field of drug discovery. In this review, a previously-established system for evaluating repositioning studies is used to highlight 4 therapeutics which exhibit pan-pathogen properties, namely azithromycin, ivermectin, niclosamide, and nitazoxanide. Recognition of the pan-pathogen nature of these antimicrobials is the cornerstone of a novel paradigm of antimicrobial development that is not only anticipatory of pandemics and bioterrorist attacks, but cognisant of conserved anti-infective mechanisms within the host-pathogen interactome which are only now beginning to emerge. Ultimately, the discovery of pan-pathogen antimicrobials is concomitantly the discovery of a new class of antivirals, and begets significant implications for pandemic preparedness research in a world after COVID-19.

Machine Learning for Prediction of Drug Targets in Microbe Associated Cardiovascular Diseases by Incorporating Host-pathogen Interaction Network Parameters

Host-pathogen interactions play a crucial role in invasion, infection, and induction of immune response in humans. In this work, four machine learning algorithms, namely Logistic regression, K-nearest neighbor, Support Vector Machine, and Random Forest were implemented for the classification of drug targets. The algorithms were trained using 3400 hosts and 3800 pathogen drug and non-drug target proteins as learning instances. For each protein, 68 pathogen and 73 host features were computed that included sequence, structure, biological and host-pathogen network centrality characteristics. The Random Forest classifier model achieved the best accuracy after 10-fold cross-validation. 99 % accuracy was achieved with a ROC-AUC score of 0.99±0.01 for both pathogen and host training sets. The Eigenvector Centrality of host-pathogen interactions and host-host interactions was the top feature in performing classification of pathogen and host targets respectively. Other features important for classification were the presence of catalytic and binding sites, low instability/aliphatic index, and cellular location. The Random Forest classifier was then used for prediction of drug targets involved in Microbe Associated Cardiovascular Diseases. 331 host and 743 pathogen proteins were predicted as drug targets by the random forest model and can be validated experimentally for therapeutic intervention in Microbe Associated Cardiovascular Diseases.

Aspergillus fumigatus versus Genus Aspergillus: Conservation, Adaptive Evolution and Specific Virulence Genes

Aspergillus is an important fungal genus containing economically important species, as well as pathogenic species of animals and plants. Using eighteen fungal species of the genus Aspergillus, we conducted a comprehensive investigation of conserved genes and their evolution. This also allows us to investigate the selection pressure driving the adaptive evolution in the pathogenic species A. fumigatus. Among single-copy orthologs (SCOs) for A. fumigatus and the closely related species A. fischeri, we identified 122 versus 50 positively selected genes (PSGs), respectively. Moreover, twenty conserved genes of unknown function were established to be positively selected and thus important for adaption. A. fumigatus PSGs interacting with human host proteins show over-representation of adaptive, symbiosis-related, immunomodulatory and virulence-related pathways, such as the TGF-β pathway, insulin receptor signaling, IL1 pathway and interfering with phagosomal GTPase signaling. Additionally, among the virulence factor coding genes, secretory and membrane protein-coding genes in multi-copy gene families, 212 genes underwent positive selection and also suggest increased adaptation, such as fungal immune evasion mechanisms (aspf2), siderophore biosynthesis (sidD), fumarylalanine production (sidE), stress tolerance (atfA) and thermotolerance (sodA). These genes presumably contribute to host adaptation strategies. Genes for the biosynthesis of gliotoxin are shared among all the close relatives of A. fumigatus as an ancient defense mechanism. Positive selection plays a crucial role in the adaptive evolution of A. fumigatus. The genome-wide profile of PSGs provides valuable targets for further research on the mechanisms of immune evasion, antimycotic targeting and understanding fundamental virulence processes.

An Ebola, Neisseria and Trypanosoma human protein interaction census reveals a conserved human protein cluster targeted by various human pathogens

Filovirus ebolavirus (ZE; Zaire ebolavirus, Bundibugyo ebolavirus), Neisseria meningitidis (NM), and Trypanosoma brucei (Tb) are serious infectious pathogens, spanning viruses, bacteria and protists and all may target the blood and central nervous system during their life cycle. NM and Tb are extracellular pathogens while ZE is obligatory intracellular, targetting immune privileged sites. By using interactomics and comparative evolutionary analysis we studied whether conserved human proteins are targeted by these pathogens. We examined 2797 unique pathogen-targeted human proteins. The information derived from orthology searches of experimentally validated protein-protein interactions (PPIs) resulted both in unique and shared PPIs for each pathogen. Comparing and analyzing conserved and pathogen-specific infection pathways for NM, TB and ZE, we identified human proteins predicted to be targeted in at least two of the compared host-pathogen networks. However, four proteins were common to all three host-pathogen interactomes: the elongation factor 1-alpha 1 (EEF1A1), the SWI/SNF complex subunit SMARCC2 (matrix-associated actin-dependent regulator of chromatin subfamily C), the dolichyl-diphosphooligosaccharide--protein glycosyltransferase subunit 1 (RPN1), and the tubulin beta-5 chain (TUBB). These four human proteins all are also involved in cytoskeleton and its regulation and are often addressed by various human pathogens. Specifically, we found (i) 56 human pathogenic bacteria and viruses that target these four proteins, (ii) the well researched new pandemic pathogen SARS-CoV-2 targets two of these four human proteins and (iii) nine human pathogenic fungi (yet another evolutionary distant organism group) target three of the conserved proteins by 130 high confidence interactions.

A Molecular Interaction Map of Klebsiella pneumoniae and Its Human Host Reveals Potential Mechanisms of Host Cell Subversion

Klebsiella pneumoniae is a leading cause of pneumonia and septicemia across the world. The rapid emergence of multidrug-resistant K. pneumoniae strains necessitates the discovery of effective drugs against this notorious pathogen. However, there is a dearth of knowledge on the mechanisms by which this deadly pathogen subverts host cellular machinery. To fill this knowledge gap, our study attempts to identify the potential mechanisms of host cell subversion by building a K. pneumoniae-human interactome based on rigorous computational methodology. The putative host targets inferred from the predicted interactome were found to be functionally enriched in the host's immune surveillance system and allied functions like apoptosis, hypoxia, etc. A multifunctionality-based scoring system revealed P53 as the most multifunctional protein among host targets accompanied by HIF1A and STAT1. Moreover, mining of host protein-protein interaction (PPI) network revealed that host targets interact among themselves to form a network (TTPPI), where P53 and CDC5L occupy a central position. The TTPPI is composed of several inter complex interactions which indicate that K. pneumoniae might disrupt functional coordination between these protein complexes through targeting of P53 and CDC5L. Furthermore, we identified four pivotal K. pneumoniae-targeted transcription factors (TTFs) that are part of TTPPI and are involved in generating host's transcriptional response to K. pneumoniae-mediated sepsis. In a nutshell, our study identifies some of the pivotal molecular targets of K. pneumoniae which primarily correlate to the physiological response of host during K. pneumoniae-mediated sepsis.

Time-Resolved Transcriptional Profiling of Epithelial Cells Infected by Intracellular Acinetobacter baumannii

The rise in the number of antibiotic-resistant bacteria has become a serious threat to health, making it important to identify, characterize and optimize new molecules to help us to overcome the infections they cause. It is well known that Acinetobacter baumannii has a significant capacity to evade the actions of antibacterial drugs, leading to its emergence as one of the bacteria responsible for hospital and community-acquired infections. Nonetheless, how this pathogen infects and survives inside the host cell is unclear. In this study, we analyze the time-resolved transcriptional profile changes observed in human epithelial HeLa cells after infection by A. baumannii, demonstrating how it survives in host cells and starts to replicate 4 h post infection. These findings were achieved by sequencing RNA to obtain a set of Differentially Expressed Genes (DEGs) to understand how bacteria alter the host cells’ environment for their own benefit. We also determine common features observed in this set of genes and identify the protein–protein networks that reveal highly-interacted proteins. The combination of these findings paves the way for the discovery of new antimicrobial candidates for the treatment of multidrug-resistant bacteria.

DualSeqDB: the host-pathogen dual RNA sequencing database for infection processes

Despite antibiotic resistance being a matter of growing concern worldwide, the bacterial mechanisms of pathogenesis remain underexplored, restraining our ability to develop new antimicrobials. The rise of high-throughput sequencing technology has made available a massive amount of transcriptomic data that could help elucidate the mechanisms underlying bacterial infection. Here, we introduce the DualSeqDB database, a resource that helps the identification of gene transcriptional changes in both pathogenic bacteria and their natural hosts upon infection. DualSeqDB comprises nearly 300 000 entries from eight different studies, with information on bacterial and host differential gene expression under in vivo and in vitro conditions. Expression data values were calculated entirely from raw data and analyzed through a standardized pipeline to ensure consistency between different studies. It includes information on seven different strains of pathogenic bacteria and a variety of cell types and tissues in Homo sapiens, Mus musculus and Macaca fascicularis at different time points. We envisage that DualSeqDB can help the research community in the systematic characterization of genes involved in host infection and help the development and tailoring of new molecules against infectious diseases. DualSeqDB is freely available at http://www.tartaglialab.com/dualseq.

Bacteria use structural imperfect mimicry to hijack the host interactome

Bacteria use protein-protein interactions to infect their hosts and hijack fundamental pathways, which ensures their survival and proliferation. Hence, the infectious capacity of the pathogen is closely related to its ability to interact with host proteins. Here, we show that hubs in the host-pathogen interactome are isolated in the pathogen network by adapting the geometry of the interacting interfaces. An imperfect mimicry of the eukaryotic interfaces allows pathogen proteins to actively bind to the host's target while preventing deleterious effects on the pathogen interactome. Understanding how bacteria recognize eukaryotic proteins may pave the way for the rational design of new antibiotic molecules.

A systematic strategy for the investigation of vaccines and drugs targeting bacteria

Infectious and epidemic diseases induced by bacteria have historically caused great distress to people, and have even resulted in a large number of deaths worldwide. At present, many researchers are working on the discovery of viable drug and vaccine targets for bacteria through multiple methods, including the analyses of comparative subtractive genome, core genome, replication-related proteins, transcriptomics and riboswitches, which plays a significant part in the treatment of infectious and pandemic diseases. The 3D structures of the desired target proteins, drugs and epitopes can be predicted and modeled through target analysis. Meanwhile, molecular dynamics (MD) analysis of the constructed drug/epitope-protein complexes is an important standard for testing the suitability of these screened drugs and vaccines. Currently, target discovery, target analysis and MD analysis are integrated into a systematic set of drug and vaccine analysis strategy for bacteria. We hope that this comprehensive strategy will help in the design of high-performance vaccines and drugs.

Identification of host cell proteins that interact with the M protein of porcine epidemic diarrhea virus

•

Interaction of PEDV M protein with host cellular proteins eIF3L, CDC42 and Rab11A was confirmed.

•

PEDV replication may be regulated by eIF3L expression.

•

218 host cell proteins were designated putative PEDV M protein interacting proteins.

Porcine epidemic diarrhea virus (PEDV) is a coronavirus that causes severe diarrhea in pigs of all ages and a high fatality rate in neonates. The PEDV membrane protein (M) plays crucial roles in viral assembly, viral budding and host immune regulation, most likely by interacting with host cell proteins that have yet to be identified. In this study, co-immunoprecipitation (Co-IP) using an M-specific monoclonal antibody, coupled with LC-MS/MS, was employed to identify M protein-interacting proteins in PEDV-infected cells. Three viral proteins (S, E and ORF3) and 218 host cell proteins were identified as putative M-interacting partners. Bioinformatic analysis showed that the identified host cell proteins were related to 131 signal pathways and 10 biological processes. In addition, interaction between translation initiation factor 3(eIF3L) and M protein was validated by Co-IP. Down-regulation of eIF3L expression significantly increased viral production, which suggests that eIF3L could be a negative regulator in PEDV replication. This interactome study of the PEDV M protein will serve to clarify its function during viral replication.

BacFITBase: a database to assess the relevance of bacterial genes during host infection

Bacterial infections have been on the rise world-wide in recent years and have a considerable impact on human well-being in terms of attributable deaths and disability-adjusted life years. Yet many mechanisms underlying bacterial pathogenesis are still poorly understood. Here, we introduce the BacFITBase database for the systematic characterization of bacterial proteins relevant for host infection aimed to enable the identification of new antibiotic targets. BacFITBase is manually curated and contains more than 90 000 entries with information on the contribution of individual genes to bacterial fitness under in vivo infection conditions in a range of host species. The data were collected from 15 different studies in which transposon mutagenesis was performed, including top-priority pathogens such as Acinetobacter baumannii and Campylobacter jejuni, for both of which increasing antibiotic resistance has been reported. Overall, BacFITBase includes information on 15 pathogenic bacteria and 5 host vertebrates across 10 different tissues. It is freely available at www.tartaglialab.com/bacfitbase.

Network analysis of host-pathogen protein interactions in microbe induced cardiovascular diseases

Large-scale visualization and analysis of HPIs involved in microbial CVDs can provide crucial insights into the mechanisms of pathogenicity. The comparison of CVD associated HPIs with the entire set of HPIs can identify the pathways specific to CVDs. Therefore, topological properties of HPI networks in CVDs and all pathogens was studied using Cytoscape3.5.1. Ontology and pathway analysis were done using KOBAS 3.0. HPIs of Papilloma, Herpes, Influenza A virus as well as Yersinia pestis and Bacillus anthracis among bacteria were predominant in the whole (wHPI) and the CVD specific (cHPI) network. The central viral and secretory bacterial proteins were predicted virulent. The central viral proteins had higher number of interactions with host proteins in comparison with bacteria. Major fraction of central and essential host proteins interacts with central viral proteins. Alpha-synuclein, Ubiquitin ribosomal proteins, TATA-box-binding protein, and Polyubiquitin-C &B proteins were the top interacting proteins specific to CVDs. Signaling by NGF, Fc epsilon receptor, EGFR and ubiquitin mediated proteolysis were among the top enriched CVD specific pathways. DEXDc and HELICc were enriched host mimicry domains that may help in hijacking of cellular machinery by pathogens. This study provides a system level understanding of cardiac damage in microbe induced CVDs.

Combination of SAXS and Protein Painting Discloses the Three-Dimensional Organization of the Bacterial Cysteine Synthase Complex, a Potential Target for Enhancers of Antibiotic Action

The formation of multienzymatic complexes allows for the fine tuning of many aspects of enzymatic functions, such as efficiency, localization, stability, and moonlighting. Here, we investigated, in solution, the structure of bacterial cysteine synthase (CS) complex. CS is formed by serine acetyltransferase (CysE) and O-acetylserine sulfhydrylase isozyme A (CysK), the enzymes that catalyze the last two steps of cysteine biosynthesis in bacteria. CysK and CysE have been proposed as potential targets for antibiotics, since cysteine and related metabolites are intimately linked to protection of bacterial cells against redox damage and to antibiotic resistance. We applied a combined approach of small-angle X-ray scattering (SAXS) spectroscopy and protein painting to obtain a model for the solution structure of CS. Protein painting allowed the identification of protein–protein interaction hotspots that were then used as constrains to model the CS quaternary assembly inside the SAXS envelope. We demonstrate that the active site entrance of CysK is involved in complex formation, as suggested by site-directed mutagenesis and functional studies. Furthermore, complex formation involves a conformational change in one CysK subunit that is likely transmitted through the dimer interface to the other subunit, with a regulatory effect. Finally, SAXS data indicate that only one active site of CysK is involved in direct interaction with CysE and unambiguously unveil the quaternary arrangement of CS.

A Coordinated Response at The Transcriptome and Interactome Level is Required to Ensure Uropathogenic Escherichia coli Survival during Bacteremia

Localized infections or disruption of the skin barrier can enable the entry of bacteria into the bloodstream, possibly leading to acute inflammation and sepsis. There is currently no holistic view on how bacteria can survive and spread in the bloodstream. In this context, we combined transposon mutagenesis, gene-expression profiling and a protein interaction network analysis to examine how uropathogenic Escherichia coli can proliferate in blood. Our results indicate that, upon migration from the urea to serum, E. coli reacts to the osmolarity difference, triggering a transcriptomic response in order to express survival genes. The proteins codified by these genes are precisely organized at the interactome level and specifically target short linear motifs located in disordered regions of host proteins. Such a coordinated response helps to explain how bacteria can adapt to and survive environmental changes within the host. Overall, our results provide a general framework for the study of bacteremia and reveal new targets for potential study as novel antimicrobials.

Analysis of protein targets in pathogen-host interaction in infectious diseases: a case study on Plasmodium falciparum and Homo sapiens interaction network

Infection and disease progression is the outcome of protein interactions between pathogen and host. Pathogen, the role player of Infection, is becoming a severe threat to life as because of its adaptability toward drugs and evolutionary dynamism in nature. Identifying protein targets by analyzing protein interactions between host and pathogen is the key point. Proteins with higher degree and possessing some topologically significant graph theoretical measures are found to be drug targets. On the other hand, exceptional nodes may be involved in infection mechanism because of some pathway process and biologically unknown factors. In this article, we attempt to investigate characteristics of host-pathogen protein interactions by presenting a comprehensive review of computational approaches applied on different infectious diseases. As an illustration, we have analyzed a case study on infectious disease malaria, with its causative agent Plasmodium falciparum acting as 'Bait' and host, Homo sapiens/human acting as 'Prey'. In this pathogen-host interaction network based on some interconnectivity and centrality properties, proteins are viewed as central, peripheral, hub and non-hub nodes and their significance on infection process. Besides, it is observed that because of sparseness of the pathogen and host interaction network, there may be some topologically unimportant but biologically significant proteins, which can also act as Bait/Prey. So, functional similarity or gene ontology mapping can help us in this case to identify these proteins.

Systems Biology of Plant-Microbiome Interactions

In natural environments, plants are exposed to diverse microbiota that they interact with in complex ways. While plant-pathogen interactions have been intensely studied to understand defense mechanisms in plants, many microbes and microbial communities can have substantial beneficial effects on their plant host. Such beneficial effects include improved acquisition of nutrients, accelerated growth, resilience against pathogens, and improved resistance against abiotic stress conditions such as heat, drought, and salinity. However, the beneficial effects of bacterial strains or consortia on their host are often cultivar and species specific, posing an obstacle to their general application. Remarkably, many of the signals that trigger plant immune responses are molecularly highly similar and often identical in pathogenic and beneficial microbes. Thus, it is unclear what determines the outcome of a particular microbe-host interaction and which factors enable plants to distinguish beneficials from pathogens. To unravel the complex network of genetic, microbial, and metabolic interactions, including the signaling events mediating microbe-host interactions, comprehensive quantitative systems biology approaches will be needed.

Analysis of Predicted Host-Parasite Interactomes Reveals Commonalities and Specificities Related to Parasitic Lifestyle and Tissues Tropism

The study of molecular host–parasite interactions is essential to understand parasitic infection and adaptation within the host system. As well, prevention and treatment of infectious diseases require a clear understanding of the molecular crosstalk between parasites and their hosts. Yet, large-scale experimental identification of host–parasite molecular interactions remains challenging, and the use of computational predictions becomes then necessary. Here, we propose a computational integrative approach to predict host—parasite protein—protein interaction (PPI) networks resulting from the human infection by 15 different eukaryotic parasites. We used an orthology-based approach to transfer high-confidence intraspecies interactions obtained from the STRING database to the corresponding interspecies homolog protein pairs in the host–parasite system. Our approach uses either the parasites predicted secretome and membrane proteins, or only the secretome, depending on whether they are uni- or multi-cellular, respectively, to reduce the number of false predictions. Moreover, the host proteome is filtered for proteins expressed in selected cellular localizations and tissues supporting the parasite growth. We evaluated the inferred interactions by analyzing the enriched biological processes and pathways in the predicted networks and their association with known parasitic invasion and evasion mechanisms. The resulting PPI networks were compared across parasites to identify common mechanisms that may define a global pathogenic hallmark. We also provided a study case focusing on a closer examination of the human–S. mansoni predicted interactome, detecting central proteins that have relevant roles in the human–S. mansoni network, and identifying tissue-specific interactions with key roles in the life cycle of the parasite. The predicted PPI networks can be visualized and downloaded at http://orthohpi.jensenlab.org.

Interaction paths promote module integration and network-level robustness of spliceosome to cascading effects

The functionality of distinct types of protein networks depends on the patterns of protein-protein interactions. A problem to solve is understanding the fragility of protein networks to predict system malfunctioning due to mutations and other errors. Spectral graph theory provides tools to understand the structural and dynamical properties of a system based on the mathematical properties of matrices associated with the networks. We combined two of such tools to explore the fragility to cascading effects of the network describing protein interactions within a key macromolecular complex, the spliceosome. Using S. cerevisiae as a model system we show that the spliceosome network has more indirect paths connecting proteins than random networks. Such multiplicity of paths may promote routes to cascading effects to propagate across the network. However, the modular network structure concentrates paths within modules, thus constraining the propagation of such cascading effects, as indicated by analytical results from the spectral graph theory and by numerical simulations of a minimal mathematical model parameterized with the spliceosome network. We hypothesize that the concentration of paths within modules favors robustness of the spliceosome against failure, but may lead to a higher vulnerability of functional subunits, which may affect the temporal assembly of the spliceosome. Our results illustrate the utility of spectral graph theory for identifying fragile spots in biological systems and predicting their implications.

Identification of Antifungal Targets Based on Computer Modeling

Aspergillus fumigatus is a saprophytic, cosmopolitan fungus that attacks patients with a weak immune system. A rational solution against fungal infection aims to manipulate fungal metabolism or to block enzymes essential for Aspergillus survival. Here we discuss and compare different bioinformatics approaches to analyze possible targeting strategies on fungal-unique pathways. For instance, phylogenetic analysis reveals fungal targets, while domain analysis allows us to spot minor differences in protein composition between the host and fungi. Moreover, protein networks between host and fungi can be systematically compared by looking at orthologs and exploiting information from host⁻pathogen interaction databases. Further data—such as knowledge of a three-dimensional structure, gene expression data, or information from calculated metabolic fluxes—refine the search and rapidly put a focus on the best targets for antimycotics. We analyzed several of the best targets for application to structure-based drug design. Finally, we discuss general advantages and limitations in identification of unique fungal pathways and protein targets when applying bioinformatics tools.

Network biology discovers pathogen contact points in host protein-protein interactomes

In all organisms, major biological processes are controlled by complex protein-protein interactions networks (interactomes), yet their structural complexity presents major analytical challenges. Here, we integrate a compendium of over 4300 phenotypes with Arabidopsis interactome (AI-1MAIN). We show that nodes with high connectivity and betweenness are enriched and depleted in conditional and essential phenotypes, respectively. Such nodes are located in the innermost layers of AI-1MAIN and are preferential targets of pathogen effectors. We extend these network-centric analyses to Cell Surface Interactome (CSILRR) and predict its 35 most influential nodes. To determine their biological relevance, we show that these proteins physically interact with pathogen effectors and modulate plant immunity. Overall, our findings contrast with centrality-lethality rule, discover fast information spreading nodes, and highlight the structural properties of pathogen targets in two different interactomes. Finally, this theoretical framework could possibly be applicable to other inter-species interactomes to reveal pathogen contact points.

Network Analysis of MERS Coronavirus within Households, Communities, and Hospitals to Identify Most Centralized and Super-Spreading in the Arabian Peninsula, 2012 to 2016

Contact history is crucial during an infectious disease outbreak and vital when seeking to understand and predict the spread of infectious diseases in human populations. The transmission connectivity networks of people infected with highly contagious Middle East respiratory syndrome coronavirus (MERS-CoV) in Saudi Arabia were assessed to identify super-spreading events among the infected patients between 2012 and 2016. Of the 1379 MERS cases recorded during the study period, 321 (23.3%) cases were linked to hospital infection, out of which 203 (14.7%) cases occurred among healthcare workers. There were 1113 isolated cases while the number of recorded contacts per MERS patient is between 1 (n=210) and 17 (n=1), with a mean of 0.27 (SD = 0.76). Five super-important nodes were identified based on their high number of connected contacts worthy of prioritization (at least degree of 5). The number of secondary cases in each SSE varies (range, 5–17). The eigenvector centrality was significantly (p < 0.05) associated with place of exposure, with hospitals having on average significantly higher eigenvector centrality than other places of exposure. Results suggested that being a healthcare worker has a higher eigenvector centrality score on average than being nonhealthcare workers. Pathogenic droplets are easily transmitted within a confined area of hospitals; therefore, control measures should be put in place to curtail the number of hospital visitors and movements of nonessential staff within the healthcare facility with MERS cases.

Elucidation of host-pathogen protein-protein interactions to uncover mechanisms of host cell rewiring

Infectious diseases are the result of molecular cross-talks between hosts and their pathogens. These cross-talks are in part mediated by host-pathogen protein-protein interactions (HP-PPI). HP-PPI play crucial roles in infections, as they may tilt the balance either in favor of the pathogens' spread or their clearance. The identification of host proteins targeted by viral or bacterial pathogenic proteins necessary for the infection can provide insights into their underlying molecular mechanisms of pathogenicity, and potentially even single out pharmacological intervention targets. Here, we review the available methods to study HP-PPI, with a focus on recent mass spectrometry based methods to decipher bacterial-human infectious diseases and examine their relevance in uncovering host cell rewiring by pathogens.