An expanded protein-protein interaction network in Bacillus subtilis reveals a group of hubs: Exploration by an integrative approach

Structural insights into RNA cleavage by a novel family of bacterial RNases

Processing of RNA is a key regulatory mechanism for all living systems. We recently discovered a novel family of endoribonucleases that is conserved across all bacteria. Here, using crystallography, cryo-EM microscopy, biochemical, biophysical, and mass spectrometry techniques, we are able to shed light on a novel RNA cleavage mechanism in bacteria. We show that YicC, the prototypical member of this family, forms a hexameric channel that closes down on a 26-mer RNA substrate, and find that it cleaves across an RNA hairpin to generate several short fragments.

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.

Development of a Dynamic Network Model to Identify Temporal Patterns of Structural Malformations in Zebrafish Embryos Exposed to a Model Toxicant, Tris(4-chlorophenyl)methanol

Embryogenesis is a well-coordinated process relying on precise cues and environmental signals that direct spatiotemporal embryonic patterning. Quite often, when one error in this process occurs, others tend to co-occur. We posit that investigating the co-occurrence of these abnormalities over time would yield additional information about the mode of toxicity for chemicals. Here, we use the environmental contaminant tris(4-chlorophenyl)methanol (TCPMOH) as a model toxicant to assess the relationship between exposures and co-occurrence of developmental abnormalities in zebrafish embryos. We propose a dynamic network modeling approach to study the co-occurrence of abnormalities, including pericardial edema, yolk sac edema, cranial malformation, spinal deformity, delayed/failed swim bladder inflation, and mortality induced by TCPMOH exposure. TCPMOH-exposed samples revealed increased abnormality co-occurrence when compared to controls. The abnormalities were represented as nodes in the dynamic network model. Abnormalities with high co-occurrence over time were identified using network centrality scores. We found that the temporal patterns of abnormality co-occurrence varied between exposure groups. In particular, the high TCPMOH exposure group experienced abnormality co-occurrence earlier than the low exposure group. The network model also revealed that pericardial and yolk sac edema are the most common critical nodes among all TCPMOH exposure levels, preceding further abnormalities. Overall, this study introduces a dynamic network model as a tool for assessing developmental toxicology, integrating structural and temporal features with a concentration response.

Protein complexes in cells by AI-assisted structural proteomics

Accurately modeling the structures of proteins and their complexes using artificial intelligence is revolutionizing molecular biology. Experimental data enable a candidate-based approach to systematically model novel protein assemblies. Here, we use a combination of in-cell crosslinking mass spectrometry and co-fractionation mass spectrometry (CoFrac-MS) to identify protein-protein interactions in the model Gram-positive bacterium Bacillus subtilis. We show that crosslinking interactions prior to cell lysis reveals protein interactions that are often lost upon cell lysis. We predict the structures of these protein interactions and others in the SubtiWiki database with AlphaFold-Multimer and, after controlling for the false-positive rate of the predictions, we propose novel structural models of 153 dimeric and 14 trimeric protein assemblies. Crosslinking MS data independently validates the AlphaFold predictions and scoring. We report and validate novel interactors of central cellular machineries that include the ribosome, RNA polymerase, and pyruvate dehydrogenase, assigning function to several uncharacterized proteins. Our approach uncovers protein-protein interactions inside intact cells, provides structural insight into their interaction interfaces, and is applicable to genetically intractable organisms, including pathogenic bacteria.

Prophage-encoded small protein YqaH counteracts the activities of the replication initiator DnaA in Bacillus subtilis

Bacterial genomes harbour cryptic prophages that are mostly transcriptionally silent with many unannotated genes. Still, cryptic prophages may contribute to their host fitness and phenotypes. In Bacillus subtilis, the yqaF-yqaN operon belongs to the prophage element skin, and is tightly repressed by the Xre-like repressor SknR. This operon contains several small ORFs (smORFs) potentially encoding small-sized proteins. The smORF-encoded peptide YqaH was previously reported to bind to the replication initiator DnaA. Here, using a yeast two-hybrid assay, we found that YqaH binds to the DNA binding domain IV of DnaA and interacts with Spo0A, a master regulator of sporulation. We isolated single amino acid substitutions in YqaH that abolished the interaction with DnaA but not with Spo0A. Then, using a plasmid-based inducible system to overexpress yqaH WT and mutant derivatives, we studied in B. subtilis the phenotypes associated with the specific loss-of-interaction with DnaA (DnaA_LOI). We found that expression of yqaH carrying DnaA_LOI mutations abolished the deleterious effects of yqaH WT expression on chromosome segregation, replication initiation and DnaA-regulated transcription. When YqaH was induced after vegetative growth, DnaA_LOI mutations abolished the drastic effects of YqaH WT on sporulation and biofilm formation. Thus, YqaH inhibits replication, sporulation and biofilm formation mainly by antagonizing DnaA in a manner that is independent of the cell cycle checkpoint Sda.

Interrelation between Stress Management and Secretion Systems of Ralstonia solanacearum: An In Silico Assessment

Ralstonia solanacearum (Rs), the causative agent of devastating wilt disease in several major and minor economic crops, is considered one of the most destructive bacterial plant pathogens. However, the mechanism(s) by which Rs counteracts host-associated environmental stress is still not clearly elucidated. To investigate possible stress management mechanisms, orthologs of stress-responsive genes in the Rs genome were searched using a reference set of known genes. The genome BLAST approach was used to find the distributions of these orthologs within different Rs strains. BLAST results were first confirmed from the KEGG Genome database and then reconfirmed at the protein level from the UniProt database. The distribution pattern of these stress-responsive factors was explored through multivariate analysis and STRING analysis. STRING analysis of stress-responsive genes in connection with different secretion systems of Rs was also performed. Initially, a total of 28 stress-responsive genes of Rs were confirmed in this study. STRING analysis revealed an additional 7 stress-responsive factors of Rs, leading to the discovery of a total of 35 stress-responsive genes. The segregation pattern of these 35 genes across 110 Rs genomes was found to be almost homogeneous. Increasing interactions of Rs stress factors were observed in six distinct clusters, suggesting six different types of stress responses: membrane stress response (MSR), osmotic stress response (OSR), oxidative stress response (OxSR), nitrosative stress response (NxSR), and DNA damage stress response (DdSR). Moreover, a strong network of these stress responses was observed with type 3 secretion system (T3SS), general secretory proteins (GSPs), and different types of pili (T4P, Tad, and Tat). To the best of our knowledge, this is the first report on overall stress response management by Rs and the potential connection with secretion systems.

SepF supports the recruitment of the DNA translocase SftA to the Z-ring

In many bacteria, cell division begins before the sister chromosomes are fully segregated. Specific DNA translocases ensure that the chromosome is removed from the closing septum, such as the transmembrane protein FtsK in Escherichia coli. Bacillus subtilis contains two FtsK homologues, SpoIIIE and SftA. SftA is active during vegetative growth whereas SpoIIIE is primarily active during sporulation and pumps the chromosome into the spore compartment. FtsK and SpoIIIE contain several transmembrane helices, however SftA is assumed to be a cytoplasmic protein. It is unknown how SftA is recruited to the cell division site. Here we show that SftA is a peripheral membrane protein, containing an N-terminal amphipathic helix that reversibly anchors the protein to the cell membrane. Using a yeast two-hybrid screen we found that SftA interacts with the conserved cell division protein SepF. Based on extensive genetic analyses and previous data we propose that the septal localization of SftA depends on either SepF or the cell division protein FtsA. Since SftA seems to interfere with the activity of SepF, and since inactivation of SepF mitigates the sensitivity of a ∆sftA mutant for ciprofloxacin, we speculate that SftA might delay septum synthesis when chromosomal DNA is in the vicinity.

Computational Network Inference for Bacterial Interactomics

Since the large-scale experimental characterization of protein–protein interactions (PPIs) is not possible for all species, several computational PPI prediction methods have been developed that harness existing data from other species. While PPI network prediction has been extensively used in eukaryotes, microbial network inference has lagged behind. However, bacterial interactomes can be built using the same principles and techniques; in fact, several methods are better suited to bacterial genomes. These predicted networks allow systems-level analyses in species that lack experimental interaction data. This review describes the current network inference and analysis techniques and summarizes the use of computationally-predicted microbial interactomes to date.

Evolutionary Analysis of the Bacillus subtilis Genome Reveals New Genes Involved in Sporulation

Bacilli can form dormant, highly resistant, and metabolically inactive spores to cope with extreme environmental challenges. In this study, we examined the evolutionary age of Bacillus subtilis sporulation genes using the approach known as genomic phylostratigraphy. We found that B. subtilis sporulation genes cluster in several groups that emerged at distant evolutionary time-points, suggesting that the sporulation process underwent several stages of expansion. Next, we asked whether such evolutionary stratification of the genome could be used to predict involvement in sporulation of presently uncharacterized genes (y-genes). We individually inactivated a representative sample of uncharacterized genes that arose during the same evolutionary periods as the known sporulation genes and tested the resulting strains for sporulation phenotypes. Sporulation was significantly affected in 16 out of 37 (43%) tested strains. In addition to expanding the knowledge base on B. subtilis sporulation, our findings suggest that evolutionary age could be used to help with genome mining.

Analysis of the RNA and Protein Complexome by Grad-seq

The complexome of a cell is the entirety of its complexes. Complexome capture studies have mostly focused on protein-protein interactions, which has left a gap in our knowledge of the global interactions of RNAs. To overcome these limitations, we recently introduced gradient profiling by sequencing (Grad-seq), which analyzes in a high-throughput fashion soluble cellular complexes after their separation in a glycerol gradient by fraction-wise RNA-seq and mass spectrometry. Here, we describe a detailed Grad-seq protocol for Streptococcus pneumoniae, which should also be applicable to other bacterial species.

Applicability of Chromatographic Co-Elution for Antibiotic Target Identification

Identification of the molecular target is a crucial step in evaluating novel antibiotics. To support target identification, a label-free method based on chromatographic co-elution has previously been developed. Target identification by chromatographic coelution (TICC) exploits the alteration of the elution profile of target-bound drug versus free drug in ion exchange (IEX) chromatography to identify potential target proteins from elution fractions. The applicability of TICC for antibiotic research is investigated by evaluating which proteins, that is, putative targets, can be monitored in Bacillus subtilis. Coelution of components of known protein complexes provides a read-out for how well the native state of proteins is conserved during chromatography. Rifampicin, which targets RNA polymerase, is used in a proof-of-concept study.

Comparative analysis of differential proteome-wide protein-protein interaction network of Methanobrevibacter ruminantium M1

A proteome-wide protein-protein interaction (PPI) network of Methanobrevibacter ruminantium M1 (MRU), a predominant rumen methanogen, was constructed from its metabolic genes using a gene neighborhood algorithm and then compared with closely related rumen methanogens Using proteome-wide PPI approach, we constructed network encompassed 2194 edges and 637 nodes interacting with 634 genes. Network quality and robustness of functional modules were assessed with gene ontology terms. A structure-function-metabolism mapping for each protein has been carried out with efforts to extract experimental PPI concomitant information from the literature. The results of our study revealed that some topological properties of its network were robust for sharing homologous protein interactions across heterotrophic and hydrogenotrophic methanogens. MRU proteome has shown to establish many PPI sub-networks for associated metabolic subsystems required to survive in the rumen environment. MRU genome found to share interacting proteins from its PPI network involved in specific metabolic subsystems distinct to heterotrophic and hydrogenotrophic methanogens. Across these proteomes, the interacting proteins from differential PPI networks were shared in common for the biosynthesis of amino acids, nucleosides, and nucleotides and energy metabolism in which more fractions of protein pairs shared with Methanosarcina acetivorans. Our comparative study expedites our knowledge to understand a complex proteome network associated with typical metabolic subsystems of MRU and to improve its genome-scale reconstruction in the future.

RocS drives chromosome segregation and nucleoid protection in Streptococcus pneumoniae

Chromosome segregation in bacteria is poorly understood outside some prominent model strains^1-5 and even less is known about how it is coordinated with other cellular processes. This is the case for the opportunistic human pathogen Streptococcus pneumoniae (the pneumococcus)⁶, which lacks the Min and the nucleoid occlusion systems⁷, and possesses only an incomplete chromosome partitioning Par(A)BS system, in which ParA is absent⁸. The bacterial tyrosine kinase⁹ CpsD, which is required for capsule production, was previously found to interfere with chromosome segregation¹⁰. Here, we identify a protein of unknown function that interacts with CpsD and drives chromosome segregation. RocS (Regulator of Chromosome Segregation) is a membrane-bound protein that interacts with both DNA and the chromosome partitioning protein ParB to properly segregate the origin of replication region to new daughter cells. In addition, we show that RocS interacts with the cell division protein FtsZ and hinders cell division. Altogether, this work reveals that RocS is the cornerstone of a nucleoid protection system ensuring proper chromosome segregation and cell division in coordination with the biogenesis of the protective capsular layer.

The Protein-Protein Interaction Network of Litopenaeus vannamei Haemocytes

Protein–protein interaction networks (PINs) have been constructed in various organisms and utilized to conduct evolutionary analyses and functional predictions. Litopenaeus vannamei is a high-valued commercial aquaculture species with an uncharacterized interactome. With the development of RNA-seq techniques and systems biology, it is possible to obtain genome-wide transcriptional information for L. vannamei and construct a systematic network based on these data. In this work, based on the RNA-seq of haemocytes we constructed the first L. vannamei PIN including 4,858 proteins and 104,187 interactions. The PIN constructed here is the first large-scale PIN for shrimp. The confidence scores of interactions in the PIN were evaluated on the basis of sequence homology and genetic relationships. The immune-specific sub-network was extracted from global PIN, and more than a third of proteins were found in signaling pathways in the sub-network, which indicates an inseparable relationship between signaling processes and immune mechanisms. Six selected signaling pathways were constructed at different age groups based on evolutionary analyses. Furthermore, we showed that the functions of the pathways’ proteins were associated with their evolutionary history based on the evolutionary analyses combining with protein functional analyses. In addition, the functions of 1,955 unclassified proteins which were associated with 3,191 unigenes were assigned using the PIN, which account for approximately 70.3 and 44.9% of the previously unclassified proteins and unigenes in the network, respectively. The annotation of unclassified proteins and unigenes based on the PIN provides new candidates for further functional studies. The immune-specific sub-network and the pathways extracted from the PIN provide a novel information source for studying of immune mechanisms and disease resistances in shrimp.

Protein-protein interactions in bacteria: a promising and challenging avenue towards the discovery of new antibiotics

Antibiotics are potent pharmacological weapons against bacterial infections; however, the growing antibiotic resistance of microorganisms is compromising the efficacy of the currently available pharmacotherapies. Even though antimicrobial resistance is not a new problem, antibiotic development has failed to match the growth of resistant pathogens and hence, it is highly critical to discover new anti-infective drugs with novel mechanisms of action which will help reducing the burden of multidrug-resistant microorganisms. Protein–protein interactions (PPIs) are involved in a myriad of vital cellular processes and have become an attractive target to treat diseases. Therefore, targeting PPI networks in bacteria may offer a new and unconventional point of intervention to develop novel anti-infective drugs which can combat the ever-increasing rate of multidrug-resistant bacteria. This review describes the progress achieved towards the discovery of molecules that disrupt PPI systems in bacteria for which inhibitors have been identified and whose targets could represent an alternative lead discovery strategy to obtain new anti-infective molecules.

Proteome Data Improves Protein Function Prediction in the Interactome of Helicobacter pylori

Helicobacter pylori is a common pathogen that is estimated to infect half of the human population, causing several diseases such as duodenal ulcer. Despite one of the first pathogens to be sequenced, its proteome remains poorly characterized as about one-third of its proteins have no functional annotation. Here, we integrate and analyze known protein interactions with proteomic and genomic data from different sources. We find that proteins with similar abundances tend to interact. Such an observation is accompanied by a trend of interactions to appear between proteins of similar functions, although some show marked cross-talk to others. Protein function prediction with protein interactions is significantly improved when interactions from other bacteria are included in our network, allowing us to obtain putative functions of more than 300 poorly or previously uncharacterized proteins. Proteins that are critical for the topological controllability of the underlying network are significantly enriched with genes that are up-regulated in the spiral compared with the coccoid form of H. pylori Determining their evolutionary conservation, we present evidence that 80 protein complexes are identical in composition with their counterparts in Escherichia coli, while 85 are partially conserved and 120 complexes are completely absent. Furthermore, we determine network clusters that coincide with related functions, gene essentiality, genetic context, cellular localization, and gene expression in different cellular states.

Metabolism Shapes the Cell

More than 5 decades of work support the idea that cell envelope synthesis, including the inward growth of cell division, is tightly coordinated with DNA replication and protein synthesis through central metabolism. Remarkably, no unifying model exists to account for how these fundamentally disparate processes are functionally coupled. Recent studies demonstrate that proteins involved in carbohydrate and nitrogen metabolism can moonlight as direct regulators of cell division, coordinate cell division and DNA replication, and even suppress defects in DNA replication. In this minireview, we focus on studies illustrating the intimate link between metabolism and regulation of peptidoglycan (PG) synthesis during growth and division, and we identify the following three recurring themes. (i) Nutrient availability, not growth rate, is the primary determinant of cell size. (ii) The degree of gluconeogenic flux is likely to have a profound impact on the metabolites available for cell envelope synthesis, so growth medium selection is a critical consideration when designing and interpreting experiments related to morphogenesis. (iii) Perturbations in pathways relying on commonly shared and limiting metabolites, like undecaprenyl phosphate (Und-P), can lead to pleotropic phenotypes in unrelated pathways.

The Protein Interactome of Streptococcus pneumoniae and Bacterial Meta-interactomes Improve Function Predictions

The functions of roughly a third of all proteins in Streptococcus pneumoniae, a significant human-pathogenic bacterium, are unknown. Using a yeast two-hybrid approach, we have determined more than 2,000 novel protein interactions in this organism. We augmented this network with meta-interactome data that we defined as the pool of all interactions between evolutionarily conserved proteins in other bacteria. We found that such interactions significantly improved our ability to predict a protein's function, allowing us to provide functional predictions for 299 S. pneumoniae proteins with previously unknown functions. IMPORTANCE Identification of protein interactions in bacterial species can help define the individual roles that proteins play in cellular pathways and pathogenesis. Very few protein interactions have been identified for the important human pathogen S. pneumoniae. We used an experimental approach to identify over 2,000 new protein interactions for S. pneumoniae, the most extensive interactome data for this bacterium to date. To predict protein function, we used our interactome data augmented with interactions from other closely related bacteria. The combination of the experimental data and meta-interactome data significantly improved the prediction results, allowing us to assign possible functions to a large number of poorly characterized proteins.

Reconstruction and Application of Protein-Protein Interaction Network

The protein-protein interaction network (PIN) is a useful tool for systematic investigation of the complex biological activities in the cell. With the increasing interests on the proteome-wide interaction networks, PINs have been reconstructed for many species, including virus, bacteria, plants, animals, and humans. With the development of biological techniques, the reconstruction methods of PIN are further improved. PIN has gradually penetrated many fields in biological research. In this work we systematically reviewed the development of PIN in the past fifteen years, with respect to its reconstruction and application of function annotation, subsystem investigation, evolution analysis, hub protein analysis, and regulation mechanism analysis. Due to the significant role of PIN in the in-depth exploration of biological process mechanisms, PIN will be preferred by more and more researchers for the systematic study of the protein systems in various kinds of organisms.

Quantitative Tagless Copurification: A Method to Validate and Identify Protein-Protein Interactions

Identifying protein-protein interactions (PPIs) at an acceptable false discovery rate (FDR) is challenging. Previously we identified several hundred PPIs from affinity purification - mass spectrometry (AP-MS) data for the bacteria Escherichia coli and Desulfovibrio vulgaris. These two interactomes have lower FDRs than any of the nine interactomes proposed previously for bacteria and are more enriched in PPIs validated by other data than the nine earlier interactomes. To more thoroughly determine the accuracy of ours or other interactomes and to discover further PPIs de novo, here we present a quantitative tagless method that employs iTRAQ MS to measure the copurification of endogenous proteins through orthogonal chromatography steps. 5273 fractions from a four-step fractionation of a D. vulgaris protein extract were assayed, resulting in the detection of 1242 proteins. Protein partners from our D. vulgaris and E. coli AP-MS interactomes copurify as frequently as pairs belonging to three benchmark data sets of well-characterized PPIs. In contrast, the protein pairs from the nine other bacterial interactomes copurify two- to 20-fold less often. We also identify 200 high confidence D. vulgaris PPIs based on tagless copurification and colocalization in the genome. These PPIs are as strongly validated by other data as our AP-MS interactomes and overlap with our AP-MS interactome for D.vulgaris within 3% of expectation, once FDRs and false negative rates are taken into account. Finally, we reanalyzed data from two quantitative tagless screens of human cell extracts. We estimate that the novel PPIs reported in these studies have an FDR of at least 85% and find that less than 7% of the novel PPIs identified in each screen overlap. Our results establish that a quantitative tagless method can be used to validate and identify PPIs, but that such data must be analyzed carefully to minimize the FDR.

Role of Protein Phosphorylation in the Regulation of Cell Cycle and DNA-Related Processes in Bacteria

In all living organisms, the phosphorylation of proteins modulates various aspects of their functionalities. In eukaryotes, protein phosphorylation plays a key role in cell signaling, gene expression, and differentiation. Protein phosphorylation is also involved in the global control of DNA replication during the cell cycle, as well as in the mechanisms that cope with stress-induced replication blocks. Similar to eukaryotes, bacteria use Hanks-type kinases and phosphatases for signal transduction, and protein phosphorylation is involved in numerous cellular processes. However, it remains unclear whether protein phosphorylation in bacteria can also regulate the activity of proteins involved in DNA-mediated processes such as DNA replication or repair. Accumulating evidence supported by functional and biochemical studies suggests that phospho-regulatory mechanisms also take place during the bacterial cell cycle. Recent phosphoproteomics and interactomics studies identified numerous phosphoproteins involved in various aspect of DNA metabolism strongly supporting the existence of such level of regulation in bacteria. Similar to eukaryotes, bacterial scaffolding-like proteins emerged as platforms for kinase activation and signaling. This review reports the current knowledge on the phosphorylation of proteins involved in the maintenance of genome integrity and the regulation of cell cycle in bacteria that reveals surprising similarities to eukaryotes.

Bacterial Interactomes: Interacting Protein Partners Share Similar Function and Are Validated in Independent Assays More Frequently Than Previously Reported

Numerous affinity purification-mass spectrometry (AP-MS) and yeast two-hybrid screens have each defined thousands of pairwise protein-protein interactions (PPIs), most of which are between functionally unrelated proteins. The accuracy of these networks, however, is under debate. Here, we present an AP-MS survey of the bacterium Desulfovibrio vulgaris together with a critical reanalysis of nine published bacterial yeast two-hybrid and AP-MS screens. We have identified 459 high confidence PPIs from D. vulgaris and 391 from Escherichia coli Compared with the nine published interactomes, our two networks are smaller, are much less highly connected, and have significantly lower false discovery rates. In addition, our interactomes are much more enriched in protein pairs that are encoded in the same operon, have similar functions, and are reproducibly detected in other physical interaction assays than the pairs reported in prior studies. Our work establishes more stringent benchmarks for the properties of protein interactomes and suggests that bona fide PPIs much more frequently involve protein partners that are annotated with similar functions or that can be validated in independent assays than earlier studies suggested.

SubtiWiki 2.0--an integrated database for the model organism Bacillus subtilis

To understand living cells, we need knowledge of each of their parts as well as about the interactions of these parts. To gain rapid and comprehensive access to this information, annotation databases are required. Here, we present SubtiWiki 2.0, the integrated database for the model bacterium Bacillus subtilis (http://subtiwiki.uni-goettingen.de/). SubtiWiki provides text-based access to published information about the genes and proteins of B. subtilis as well as presentations of metabolic and regulatory pathways. Moreover, manually curated protein-protein interactions diagrams are linked to the protein pages. Finally, expression data are shown with respect to gene expression under 104 different conditions as well as absolute protein quantification for cytoplasmic proteins. To facilitate the mobile use of SubtiWiki, we have now expanded it by Apps that are available for iOS and Android devices. Importantly, the App allows to link private notes and pictures to the gene/protein pages. Today, SubtiWiki has become one of the most complete collections of knowledge on a living organism in one single resource.

MreB-Dependent Inhibition of Cell Elongation during the Escape from Competence in Bacillus subtilis

During bacterial exponential growth, the morphogenetic actin-like MreB proteins form membrane-associated assemblies that move processively following trajectories perpendicular to the long axis of the cell. Such MreB structures are thought to scaffold and restrict the movement of peptidoglycan synthesizing machineries, thereby coordinating sidewall elongation. In Bacillus subtilis, this function is performed by the redundant action of three MreB isoforms, namely MreB, Mbl and MreBH. mreB and mbl are highly transcribed from vegetative promoters. We have found that their expression is maximal at the end of exponential phase, and rapidly decreases to a low basal level upon entering stationary phase. However, in cells developing genetic competence, a stationary phase physiological adaptation, expression of mreB was specifically reactivated by the central competence regulator ComK. In competent cells, MreB was found in complex with several competence proteins by in vitro pull-down assays. In addition, it co-localized with the polar clusters formed by the late competence peripheral protein ComGA, in a ComGA-dependent manner. ComGA has been shown to be essential for the inhibition of cell elongation characteristic of cells escaping the competence state. We show here that the pathway controlling this elongation inhibition also involves MreB. Our findings suggest that ComGA sequesters MreB to prevent cell elongation and therefore the escape from competence.

In bacterial cells, like in their eukaryotic counterparts, precise spatiotemporal localization of proteins is critical for their cellular function. This study shows that the expression and the localization of the bacterial actin-like MreB protein are growth phase-dependent. During exponential growth, we previously showed that MreB, together with other morphogenetic factors, forms discrete assemblies that move in a directed manner along peripheral tracks. Here, we demonstrate that in cells that develop genetic competence during stationary phase, transcription of mreB is specifically activated and MreB relocalizes to the cell poles. Our findings suggest a model in which MreB sequestration by the late competence protein ComGA prevents cell elongation during the escape from competence.

Single cell FRET analysis for the identification of optimal FRET-pairs in Bacillus subtilis using a prototype MEM-FLIM system

Protein-protein interactions can be studied in vitro, e.g. with bacterial or yeast two-hybrid systems or surface plasmon resonance. In contrast to in vitro techniques, in vivo studies of protein-protein interactions allow examination of spatial and temporal behavior of such interactions in their native environment. One approach to study protein-protein interactions in vivo is via Förster Resonance Energy Transfer (FRET). Here, FRET efficiency of selected FRET-pairs was studied at the single cell level using sensitized emission and Frequency Domain-Fluorescence Lifetime Imaging Microscopy (FD-FLIM). For FRET-FLIM, a prototype Modulated Electron-Multiplied FLIM system was used, which is, to the best of our knowledge, the first account of Frequency Domain FLIM to analyze FRET in single bacterial cells. To perform FRET-FLIM, we first determined and benchmarked the best fluorescent protein-pair for FRET in Bacillus subtilis using a novel BglBrick-compatible integration vector. We show that GFP-tagRFP is an excellent donor-acceptor pair for B. subtilis in vivo FRET studies. As a proof of concept, selected donor and acceptor fluorescent proteins were fused using a linker that contained a tobacco etch virus (TEV)-protease recognition sequence. Induction of TEV-protease results in loss of FRET efficiency and increase in fluorescence lifetime. The loss of FRET efficiency after TEV induction can be followed in time in single cells via time-lapse microscopy. This work will facilitate future studies of in vivo dynamics of protein complexes in single B. subtilis cells.

Network representations of immune system complexity

The mammalian immune system is a dynamic multiscale system composed of a hierarchically organized set of molecular, cellular, and organismal networks that act in concert to promote effective host defense. These networks range from those involving gene regulatory and protein–protein interactions underlying intracellular signaling pathways and single‐cell responses to increasingly complex networks of in vivo cellular interaction, positioning, and migration that determine the overall immune response of an organism. Immunity is thus not the product of simple signaling events but rather nonlinear behaviors arising from dynamic, feedback‐regulated interactions among many components. One of the major goals of systems immunology is to quantitatively measure these complex multiscale spatial and temporal interactions, permitting development of computational models that can be used to predict responses to perturbation. Recent technological advances permit collection of comprehensive datasets at multiple molecular and cellular levels, while advances in network biology support representation of the relationships of components at each level as physical or functional interaction networks. The latter facilitate effective visualization of patterns and recognition of emergent properties arising from the many interactions of genes, molecules, and cells of the immune system. We illustrate the power of integrating ‘omics’ and network modeling approaches for unbiased reconstruction of signaling and transcriptional networks with a focus on applications involving the innate immune system. We further discuss future possibilities for reconstruction of increasingly complex cellular‐ and organism‐level networks and development of sophisticated computational tools for prediction of emergent immune behavior arising from the concerted action of these networks. WIREs Syst Biol Med 2015, 7:13–38. doi: 10.1002/wsbm.1288

This article is categorized under: 1

Analytical and Computational Methods > Computational Methods

2

Laboratory Methods and Technologies > Macromolecular Interactions, Methods

Applications of omics approaches to the development of microbiological risk assessment using RNA virus dose-response models as a case study

T e in the amount of ‘omics’ data available and in our ability to interpret those data. The aim of this paper was to consider how omics techniques can be used to improve and refine microbiological risk assessment, using dose–response models for RNA viruses, with particular reference to norovirus through the oral route as the case study. The dose–response model for initial infection in the gastrointestinal tract is broken down into the component steps at the molecular level and the feasibility of assigning probabilities to each step assessed. The molecular mechanisms are not sufficiently well understood at present to enable quantitative estimation of probabilities on the basis of omics data. At present, the great strength of gene sequence data appears to be in giving information on the distribution and proportion of susceptible genotypes (for example due to the presence of the appropriate pathogen‐binding receptor) in the host population rather than in predicting specificities from the amino acid sequences concurrently obtained. The nature of the mutant spectrum in RNA viruses greatly complicates the application of omics approaches to the development of mechanistic dose–response models and prevents prediction of risks of disease progression (given infection has occurred) at the level of the individual host. However, molecular markers in the host and virus may enable more broad predictions to be made about the consequences of exposure in a population. In an alternative approach, comparing the results of deep sequencing of RNA viruses in the faeces/vomitus from donor humans with those from their infected recipients may enable direct estimates of the average probability of infection per virion to be made.

Protein-tyrosine phosphorylation interaction network in Bacillus subtilis reveals new substrates, kinase activators and kinase cross-talk

Signal transduction in eukaryotes is generally transmitted through phosphorylation cascades that involve a complex interplay of transmembrane receptors, protein kinases, phosphatases and their targets. Our previous work indicated that bacterial protein-tyrosine kinases and phosphatases may exhibit similar properties, since they act on many different substrates. To capture the complexity of this phosphorylation-based network, we performed a comprehensive interactome study focused on the protein-tyrosine kinases and phosphatases in the model bacterium Bacillus subtilis. The resulting network identified many potential new substrates of kinases and phosphatases, some of which were experimentally validated. Our study highlighted the role of tyrosine and serine/threonine kinases and phosphatases in DNA metabolism, transcriptional control and cell division. This interaction network reveals significant crosstalk among different classes of kinases. We found that tyrosine kinases can bind to several modulators, transmembrane or cytosolic, consistent with a branching of signaling pathways. Most particularly, we found that the division site regulator MinD can form a complex with the tyrosine kinase PtkA and modulate its activity in vitro. In vivo, it acts as a scaffold protein which anchors the kinase at the cell pole. This network highlighted a role of tyrosine phosphorylation in the spatial regulation of the Z-ring during cytokinesis.

The ESX system in Bacillus subtilis mediates protein secretion

Esat-6 protein secretion systems (ESX or Ess) are required for the virulence of several human pathogens, most notably Mycobacterium tuberculosis and Staphylococcus aureus. These secretion systems are defined by a conserved FtsK/SpoIIIE family ATPase and one or more WXG100 family secreted substrates. Gene clusters coding for ESX systems have been identified amongst many organisms including the highly tractable model system, Bacillus subtilis. In this study, we demonstrate that the B. subtilis yuk/yue locus codes for a nonessential ESX secretion system. We develop a functional secretion assay to demonstrate that each of the locus gene products is specifically required for secretion of the WXG100 virulence factor homolog, YukE. We then employ an unbiased approach to search for additional secreted substrates. By quantitative profiling of culture supernatants, we find that YukE may be the sole substrate that depends on the FtsK/SpoIIIE family ATPase for secretion. We discuss potential functional implications for secretion of a unique substrate.

The protein-protein interaction network of eyestalk, Y-organ and hepatopancreas in Chinese mitten crab Eriocheir sinensis

Background

The protein-protein interaction network (PIN) is an effective information tool for understanding the complex biological processes inside the cell and solving many biological problems such as signaling pathway identification and prediction of protein functions. Eriocheir sinensis is a highly-commercial aquaculture species with an unclear proteome background which hinders the construction and development of PIN for E. sinensis. However, in recent years, the development of next-generation deep-sequencing techniques makes it possible to get high throughput data of E. sinensis tanscriptome and subsequently obtain a systematic overview of the protein-protein interaction system.

Results

In this work we sequenced the transcriptional RNA of eyestalk, Y-organ and hepatopancreas in E. sinensis and generated a PIN of E. sinensis which included 3,223 proteins and 35,787 interactions. Each protein-protein interaction in the network was scored according to the homology and genetic relationship. The signaling sub-network, representing the signal transduction pathways in E. sinensis, was extracted from the global network, which depicted a global view of the signaling systems in E. sinensis. Seven basic signal transduction pathways were identified in E. sinensis. By investigating the evolution paths of the seven pathways, we found that these pathways got mature in different evolutionary stages. Moreover, the functions of unclassified proteins and unigenes in the PIN of E. sinensis were predicted. Specifically, the functions of 549 unclassified proteins related to 864 unclassified unigenes were assigned, which respectively covered 76% and 73% of all the unclassified proteins and unigenes in the network.

Conclusions

The PIN generated in this work is the first large-scale PIN of aquatic crustacean, thereby providing a paradigmatic blueprint of the aquatic crustacean interactome. Signaling sub-network extracted from the global PIN depicts the interaction of different signaling proteins and the evolutionary paths of the identified signal transduction pathways. Furthermore, the function assignment of unclassified proteins based on the PIN offers a new reference in protein function exploration. More importantly, the construction of the E. sinensis PIN provides necessary experience for the exploration of PINs in other aquatic crustacean species.

Bacillus subtilis Systems Biology: Applications of -Omics Techniques to the Study of Endospore Formation

Endospore-forming bacteria, with Bacillus subtilis being the prevalent model organism, belong to the phylum Firmicutes. Although the last common ancestor of all Firmicutes is likely to have been an endospore-forming species, not every lineage in the phylum has maintained the ability to produce endospores (hereafter, spores). In 1997, the release of the full genome sequence for B. subtilis strain 168 marked the beginning of the genomic era for the study of spore formation (sporulation). In this original genome sequence, 139 of the 4,100 protein-coding genes were annotated as sporulation genes. By the time a revised genome sequence with updated annotations was published in 2009, that number had increased significantly, especially since transcriptional profiling studies (transcriptomics) led to the identification of several genes expressed under the control of known sporulation transcription factors. Over the past decade, genome sequences for multiple spore-forming species have been released (including several strains in the Bacillus anthracis / Bacillus cereus group and many Clostridium species), and phylogenomic analyses have revealed many conserved sporulation genes. Parallel advances in transcriptomics led to the identification of small untranslated regulatory RNAs (sRNAs), including some that are expressed during sporulation. An extended array of -omics techniques, i.e., techniques designed to probe gene function on a genome-wide scale, such as proteomics, metabolomics, and high-throughput protein localization studies, have been implemented in microbiology. Combined with the use of new computational methods for predicting gene function and inferring regulatory relationships on a global scale, these -omics approaches are uncovering novel information about sporulation and a variety of other bacterial cell processes.

A second-generation protein-protein interaction network of Helicobacter pylori

Helicobacter pylori infections cause gastric ulcers and play a major role in the development of gastric cancer. In 2001, the first protein interactome was published for this species, revealing over 1500 binary protein interactions resulting from 261 yeast two-hybrid screens. Here we roughly double the number of previously published interactions using an ORFeome-based, proteome-wide yeast two-hybrid screening strategy. We identified a total of 1515 protein-protein interactions, of which 1461 are new. The integration of all the interactions reported in H. pylori results in 3004 unique interactions that connect about 70% of its proteome. Excluding interactions of promiscuous proteins we derived from our new data a core network consisting of 908 interactions. We compared our data set to several other bacterial interactomes and experimentally benchmarked the conservation of interactions using 365 protein pairs (interologs) of E. coli of which one third turned out to be conserved in both species.

DNA double strand break end-processing and RecA induce RecN expression levels in Bacillus subtilis

Bacillus subtilis cells respond to double strand breaks (DSBs) with an ordered recruitment of repair proteins to the site lesion, being RecN one of the first responders. In B. subtilis, one of the responses to DSBs is to increase RecN expression rather than modifying its turnover rate. End-processing activities and the RecA protein itself contribute to increase RecN levels after DNA DSBs. RecO is required for RecA filament formation and full SOS induction, but its absence did not significantly affect RecN expression. Neither the absence of LexA nor the phosphorylation state of RecA or SsbA significantly affect RecN expression levels. These findings identify two major mechanisms (SOS and DSB response) used to respond to DSBs, with LexA required for one of them (SOS response). The DSB response, which requires end-processing and RecA or short RecO-independent RecA filaments, highlights the importance of guarding genome stability by modulating the DNA damage responses.

An early cytoplasmic step of peptidoglycan synthesis is associated to MreB in Bacillus subtilis

MreB proteins play a major role during morphogenesis of rod-shaped bacteria by organizing biosynthesis of the peptidoglycan cell wall. However, the mechanisms underlying this process are not well understood. In Bacillus subtilis, membrane-associated MreB polymers have been shown to be associated to elongation-specific complexes containing transmembrane morphogenetic factors and extracellular cell wall assembly proteins. We have now found that an early intracellular step of cell wall synthesis is also associated to MreB. We show that the previously uncharacterized protein YkuR (renamed DapI) is required for synthesis of meso-diaminopimelate (m-DAP), an essential constituent of the peptidoglycan precursor, and that it physically interacts with MreB. Highly inclined laminated optical sheet microscopy revealed that YkuR forms uniformly distributed foci that exhibit fast motion in the cytoplasm, and are not detected in cells lacking MreB. We propose a model in which soluble MreB organizes intracellular steps of peptidoglycan synthesis in the cytoplasm to feed the membrane-associated cell wall synthesizing machineries.

Interaction of bacterial fatty-acid-displaced regulators with DNA is interrupted by tyrosine phosphorylation in the helix-turn-helix domain

Bacteria possess transcription regulators (of the TetR family) specifically dedicated to repressing genes for cytochrome P450, involved in oxidation of polyunsaturated fatty acids. Interaction of these repressors with operator sequences is disrupted in the presence of fatty acids, and they are therefore known as fatty-acid-displaced regulators. Here, we describe a novel mechanism of inactivating the interaction of these proteins with DNA, illustrated by the example of Bacillus subtilis regulator FatR. FatR was found to interact in a two-hybrid assay with TkmA, an activator of the protein-tyrosine kinase PtkA. We show that FatR is phosphorylated specifically at the residue tyrosine 45 in its helix-turn-helix domain by the kinase PtkA. Structural modelling reveals that the hydroxyl group of tyrosine 45 interacts with DNA, and we show that this phosphorylation reduces FatR DNA binding capacity. Point mutants mimicking phosphorylation of FatR in vivo lead to a strong derepression of the fatR operon, indicating that this regulatory mechanism works independently of derepression by polyunsaturated fatty acids. Tyrosine 45 is a highly conserved residue, and PtkA from B. subtilis can phosphorylate FatR homologues from other bacteria. This indicates that phosphorylation of tyrosine 45 may be a general mechanism of switching off bacterial fatty-acid-displaced regulators.

Mapping the twin-arginine protein translocation network of Bacillus subtilis

Bacteria employ twin-arginine translocation (Tat) pathways for the transport of folded proteins to extracytoplasmic destinations. In recent years, most studies on bacterial Tat pathways addressed the membrane-bound TatA(B)C subunits of the Tat translocase, and the specific interactions between this translocase and its substrate proteins. In contrast, relatively few studies investigated possible coactors in the TatA(B)C-dependent protein translocation process. The present studies were aimed at identifying interaction partners of the Tat pathway of Bacillus subtilis, which is a paradigm for studies on protein secretion by Gram-positive bacteria. Specifically, 36 interaction partners of the TatA and TatC subunits were identified by rigorous application of the yeast two-hybrid (Y2H) approach. Our Y2H analyses revealed that the three TatA isoforms of B. subtilis can form homo- and heterodimers. Subsequently, the secretion of the Tat substrates YwbN and PhoD was tested in mutant strains lacking genes for the TatAC interaction partners identified in our genome-wide Y2H screens. Our results show that the cell wall-bound protease WprA is important for YwbN secretion, and that the HemAT and CsbC proteins are required for PhoD secretion under phosphate starvation conditions. Taken together, our findings imply that the Bacillus Tat pathway is embedded in an intricate protein-protein interaction network.

Systems biology of pathogen-host interaction: networks of protein-protein interaction within pathogens and pathogen-human interactions in the post-genomic era

Infectious diseases comprise some of the leading causes of death and disability worldwide. Interactions between pathogen and host proteins underlie the process of infection. Improved understanding of pathogen-host molecular interactions will increase our knowledge of the mechanisms involved in infection, and allow novel therapeutic solutions to be devised. Complete genome sequences for a number of pathogenic microorganisms, as well as the human host, has led to the revelation of their protein-protein interaction (PPI) networks. In this post-genomic era, pathogen-host interactions (PHIs) operating during infection can also be mapped. Detailed systematic analyses of PPI and PHI data together are required for a complete understanding of pathogenesis of infections. Here we review the striking results recently obtained during the construction and investigation of these networks. Emphasis is placed on studies producing large-scale interaction data by high-throughput experimental techniques.

Diversity in genetic in vivo methods for protein-protein interaction studies: from the yeast two-hybrid system to the mammalian split-luciferase system

The yeast two-hybrid system pioneered the field of in vivo protein-protein interaction methods and undisputedly gave rise to a palette of ingenious techniques that are constantly pushing further the limits of the original method. Sensitivity and selectivity have improved because of various technical tricks and experimental designs. Here we present an exhaustive overview of the genetic approaches available to study in vivo binary protein interactions, based on two-hybrid and protein fragment complementation assays. These methods have been engineered and employed successfully in microorganisms such as Saccharomyces cerevisiae and Escherichia coli, but also in higher eukaryotes. From single binary pairwise interactions to whole-genome interactome mapping, the self-reassembly concept has been employed widely. Innovative studies report the use of proteins such as ubiquitin, dihydrofolate reductase, and adenylate cyclase as reconstituted reporters. Protein fragment complementation assays have extended the possibilities in protein-protein interaction studies, with technologies that enable spatial and temporal analyses of protein complexes. In addition, one-hybrid and three-hybrid systems have broadened the types of interactions that can be studied and the findings that can be obtained. Applications of these technologies are discussed, together with the advantages and limitations of the available assays.

Physical interaction between coat morphogenetic proteins SpoVID and CotE is necessary for spore encasement in Bacillus subtilis

Endospore formation by Bacillus subtilis is a complex and dynamic process. One of the major challenges of sporulation is the assembly of a protective, multilayered, proteinaceous spore coat, composed of at least 70 different proteins. Spore coat formation can be divided into two distinct stages. The first is the recruitment of proteins to the spore surface, dependent on the morphogenetic protein SpoIVA. The second step, known as encasement, involves the migration of the coat proteins around the circumference of the spore in successive waves, a process dependent on the morphogenetic protein SpoVID and the transcriptional regulation of individual coat genes. We provide genetic and biochemical evidence supporting the hypothesis that SpoVID promotes encasement of the spore by establishing direct protein-protein interactions with other coat morphogenetic proteins. It was previously demonstrated that SpoVID directly interacts with SpoIVA and the inner coat morphogenetic protein, SafA. Here, we show by yeast two-hybrid and pulldown assays that SpoVID also interacts directly with the outer coat morphogenetic protein, CotE. Furthermore, by mutational analysis, we identified a specific residue in the N-terminal domain of SpoVID that is essential for the interaction with CotE but dispensable for the interaction with SafA. We propose an updated model of coat assembly and spore encasement that incorporates several physical interactions between the principal coat morphogenetic proteins.

SubtiWiki--a comprehensive community resource for the model organism Bacillus subtilis

In the post-genomic era, most components of a cell are known and they can be quantified by large-scale functional genomics approaches. However, genome annotation is the bottleneck that hampers our understanding of living cells and organisms. Up-to-date functional annotation is of special importance for model organisms that provide a frame of reference for studies with other relevant organisms. We have generated a Wiki-type database for the Gram-positive model bacterium Bacillus subtilis, SubtiWiki (http://subtiwiki.uni-goettingen.de/). This Wiki is centered around the individual genes and gene products of B. subtilis and provides information on each aspect of gene function and expression as well as protein activity and its control. SubtiWiki is accompanied by two companion databases SubtiPathways and SubtInteract that provide graphical representations of B. subtilis metabolism and its regulation and of protein-protein interactions, respectively. The diagrams of both databases are easily navigatable using the popular Google maps API, and they are extensively linked with the SubtiWiki gene pages. Moreover, each gene/gene product was assigned to one or more functional categories and transcription factor regulons. Pages for the specific categories and regulons provide a rapid overview of functionally related genes/proteins. Today, SubtiWiki can be regarded as one of the most complete inventories of knowledge on a living organism in one single resource.