High-quality binary protein interaction map of the yeast interactome network

Structural and mechanistic basis for preferential deadenylation of U6 snRNA by Usb1

Post-transcriptional modification of snRNA is central to spliceosome function. Usb1 is an exoribonuclease that shortens the oligo-uridine tail of U6 snRNA, resulting in a terminal 2',3' cyclic phosphate group in most eukaryotes, including humans. Loss of function mutations in human Usb1 cause the rare disorder poikiloderma with neutropenia (PN), and result in U6 snRNAs with elongated 3' ends that are aberrantly adenylated. Here, we show that human Usb1 removes 3' adenosines with 20-fold greater efficiency than uridines, which explains the presence of adenylated U6 snRNAs in cells lacking Usb1. We determined three high-resolution co-crystal structures of Usb1: wild-type Usb1 bound to the substrate analog adenosine 5'-monophosphate, and an inactive mutant bound to RNAs with a 3' terminal adenosine and uridine. These structures, along with QM/MM MD simulations of the catalytic mechanism, illuminate the molecular basis for preferential deadenylation of U6 snRNA. The extent of Usb1 processing is influenced by the secondary structure of U6 snRNA.

Metabolic network-based stratification of hepatocellular carcinoma reveals three distinct tumor subtypes

Hepatocellular carcinoma (HCC) is one of the most frequent forms of liver cancer, and effective treatment methods are limited due to tumor heterogeneity. There is a great need for comprehensive approaches to stratify HCC patients, gain biological insights into subtypes, and ultimately identify effective therapeutic targets. We stratified HCC patients and characterized each subtype using transcriptomics data, genome-scale metabolic networks and network topology/controllability analysis. This comprehensive systems-level analysis identified three distinct subtypes with substantial differences in metabolic and signaling pathways reflecting at genomic, transcriptomic, and proteomic levels. These subtypes showed large differences in clinical survival associated with altered kynurenine metabolism, WNT/β-catenin-associated lipid metabolism, and PI3K/AKT/mTOR signaling. Integrative analyses indicated that the three subtypes rely on alternative enzymes (e.g., ACSS1/ACSS2/ACSS3, PKM/PKLR, ALDOB/ALDOA, MTHFD1L/MTHFD2/MTHFD1) to catalyze the same reactions. Based on systems-level analysis, we identified 8 to 28 subtype-specific genes with pivotal roles in controlling the metabolic network and predicted that these genes may be targeted for development of treatment strategies for HCC subtypes by performing in silico analysis. To validate our predictions, we performed experiments using HepG2 cells under normoxic and hypoxic conditions and observed opposite expression patterns between genes expressed in high/moderate/low-survival tumor groups in response to hypoxia, reflecting activated hypoxic behavior in patients with poor survival. In conclusion, our analyses showed that the heterogeneous HCC tumors can be stratified using a metabolic network-driven approach, which may also be applied to other cancer types, and this stratification may have clinical implications to drive the development of precision medicine.

Similarity-based future common neighbors model for link prediction in complex networks

Link prediction aims to predict the existence of unknown links via the network information. However, most similarity-based algorithms only utilize the current common neighbor information and cannot get high enough prediction accuracy in evolving networks. So this paper firstly defines the future common neighbors that can turn into the common neighbors in the future. To analyse whether the future common neighbors contribute to the current link prediction, we propose the similarity-based future common neighbors (SFCN) model for link prediction, which accurately locate all the future common neighbors besides the current common neighbors in networks and effectively measure their contributions. We also design and observe three MATLAB simulation experiments. The first experiment, which adjusts two parameter weights in the SFCN model, reveals that the future common neighbors make more contributions than the current common neighbors in complex networks. And two more experiments, which compares the SFCN model with eight algorithms in five networks, demonstrate that the SFCN model has higher accuracy and better performance robustness.

Changes of Cell Biochemical States Are Revealed in Protein Homomeric Complex Dynamics

We report here a simple and global strategy to map out gene functions and target pathways of drugs, toxins, or other small molecules based on "homomer dynamics" protein-fragment complementation assays (hdPCA). hdPCA measures changes in self-association (homomerization) of over 3,500 yeast proteins in yeast grown under different conditions. hdPCA complements genetic interaction measurements while eliminating the confounding effects of gene ablation. We demonstrate that hdPCA accurately predicts the effects of two longevity and health span-affecting drugs, the immunosuppressant rapamycin and the type 2 diabetes drug metformin, on cellular pathways. We also discovered an unsuspected global cellular response to metformin that resembles iron deficiency and includes a change in protein-bound iron levels. This discovery opens a new avenue to investigate molecular mechanisms for the prevention or treatment of diabetes, cancers, and other chronic diseases of aging.

Global Identification of Protein Complexes within the Membrane Proteome of Arabidopsis Roots Using a SEC-MS Approach

Biological processes consist of several consecutive and interacting steps as, for example, in signal transduction cascades or metabolic reaction chains. These processes are regulated by protein-protein interactions and the formation of larger protein complexes, which also occur within biological membranes. To gain a large-scale overview of complex-forming proteins and the composition of such complexes within the cellular membranes of Arabidopsis roots, we use the combination of size-exclusion chromatography and mass spectrometry. First, we identified complex-forming proteins by a retention shift analysis relative to expected retention times of monomeric proteins during size-exclusion chromatography. In a second step we predicted complex composition through pairwise correlation of elution profiles. As result we present an interactome of 963 proteins within cellular membranes of Arabidopsis roots. Identification of complex-forming proteins was highly robust between two independently grown root proteomes. The protein complex composition derived from pairwise correlations of coeluting proteins reproducibly identified stable protein complexes (ribosomes, proteasome, mitochondrial respiratory chain supercomplexes) but showed higher variance between replicates regarding transient interactions (e.g., interactions with kinases) within membrane protein complexes.

A genome-scale metabolic network alignment method within a hypergraph-based framework using a rotational tensor-vector product

Biological network alignment aims to discover important similarities and differences and thus find a mapping between topological and/or functional components of different biological molecular networks. Then, the mapped components can be considered to correspond to both their places in the network topology and their biological attributes. Development and evolution of biological network alignment methods has been accelerated by the rapidly increasing availability of such biological networks, yielding a repertoire of tens of methods based upon graph theory. However, most biological processes, especially the metabolic reactions, are more sophisticated than simple pairwise interactions and contain three or more participating components. Such multi-lateral relations are not captured by graphs, and computational methods to overcome this limitation are currently lacking. This paper introduces hypergraphs and association hypergraphs to describe metabolic networks and their potential alignments, respectively. Within this framework, metabolic networks are aligned by identifying the maximal Z-eigenvalue of a symmetric tensor. A shifted higher-order power method was utilized to identify a solution. A rotational strategy has been introduced to accelerate the tensor-vector product by 250-fold on average and reduce the storage cost by up to 1,000-fold. The algorithm was implemented on a spark-based distributed computation cluster to significantly increase the convergence rate further by 50- to 80-fold. The parameters have been explored to understand their impact on alignment accuracy and speed. In particular, the influence of initial value selection on the stationary point has been simulated to ensure an accurate approximation of the global optimum. This framework was demonstrated by alignments among the genome-wide metabolic networks of Escherichia coli MG-1655 and Halophilic archaeon DL31. To our knowledge, this is the first genome-wide metabolic network alignment at both the metabolite level and the enzyme level. These results demonstrate that it can supply quite a few valuable insights into metabolic networks. First, this method can access the driving force of organic reactions through the chemical evolution of metabolic network. Second, this method can incorporate the chemical information of enzymes and structural changes of compounds to offer new way defining reaction class and module, such as those in KEGG. Third, as a vertex-focused treatment, this method can supply novel structural and functional annotation for ill-defined molecules. The related source code is available on request.

Resurrected Protein Interaction Networks Reveal the Innovation Potential of Ancient Whole-Genome Duplication

The evolution of plants is characterized by whole-genome duplications, sometimes closely associated with the origin of large groups of species. The gamma (γ) genome triplication occurred at the origin of the core eudicots, which comprise ∼75% of flowering plants. To better understand the impact of whole-genome duplication, we studied the protein interaction network of MADS domain transcription factors, which are key regulators of reproductive development. We reconstructed, synthesized, and tested the interactions of ancestral proteins immediately before and closely after the triplication and directly compared these ancestral networks to the extant networks of Arabidopsis thaliana and tomato (Solanum lycopersicum). We found that gamma expanded the MADS domain interaction network more strongly than subsequent genomic events. This event strongly rewired MADS domain interactions and allowed for the evolution of new functions and installed robustness through new redundancy. Despite extensive rewiring, the organization of the network was maintained through gamma. New interactions and protein retention compensated for its potentially destructive impact on network organization. Post gamma, the network evolved from an organization around the single hub SEP3 to a network organized around multiple hubs and well-connected proteins lost, rather than gained, interactions. The data provide a resource for comparative developmental biology in flowering plants.

Multiscale effects of heating and cooling on genes and gene networks

Fluctuating environments such as changes in ambient temperature represent a fundamental challenge to life. Cells must protect gene networks that protect them from such stresses, making it difficult to understand how temperature affects gene network function in general. Here, we focus on single genes and small synthetic network modules to reveal four key effects of nonoptimal temperatures at different biological scales: (i) a cell fate choice between arrest and resistance, (ii) slower growth rates, (iii) Arrhenius reaction rates, and (iv) protein structure changes. We develop a multiscale computational modeling framework that captures and predicts all of these effects. These findings promote our understanding of how temperature affects living systems and enables more robust cellular engineering for real-world applications.

Most organisms must cope with temperature changes. This involves genes and gene networks both as subjects and agents of cellular protection, creating difficulties in understanding. Here, we study how heating and cooling affect expression of single genes and synthetic gene circuits in Saccharomyces cerevisiae. We discovered that nonoptimal temperatures induce a cell fate choice between stress resistance and growth arrest. This creates dramatic gene expression bimodality in isogenic cell populations, as arrest abolishes gene expression. Multiscale models incorporating population dynamics, temperature-dependent growth rates, and Arrhenius scaling of reaction rates captured the effects of cooling, but not those of heating in resistant cells. Molecular-dynamics simulations revealed how heating alters the conformational dynamics of the TetR repressor, fully explaining the experimental observations. Overall, nonoptimal temperatures induce a cell fate decision and corrupt gene and gene network function in computationally predictable ways, which may aid future applications of engineered microbes in nonstandard temperatures.

A type 2 diabetes disease module with a high collective influence for Cdk2 and PTPLAD1 is localized in endosomes

Despite the identification of many susceptibility genes our knowledge of the underlying mechanisms responsible for complex disease remains limited. Here, we identified a type 2 diabetes disease module in endosomes, and validate it for functional relevance on selected nodes. Using hepatic Golgi/endosomes fractions, we established a proteome of insulin receptor-containing endosomes that allowed the study of physical protein interaction networks on a type 2 diabetes background. The resulting collated network is formed by 313 nodes and 1147 edges with a topology organized around a few major hubs with Cdk2 displaying the highest collective influence. Overall, 88% of the nodes are associated with the type 2 diabetes genetic risk, including 101 new candidates. The Type 2 diabetes module is enriched with cytoskeleton and luminal acidification–dependent processes that are shared with secretion-related mechanisms. We identified new signaling pathways driven by Cdk2 and PTPLAD1 whose expression affects the association of the insulin receptor with TUBA, TUBB, the actin component ACTB and the endosomal sorting markers Rab5c and Rab11a. Therefore, the interactome of internalized insulin receptors reveals the presence of a type 2 diabetes disease module enriched in new layers of feedback loops required for insulin signaling, clearance and islet biology.

The budding-yeast RWD protein Csm1 scaffolds diverse protein complexes through a conserved structural mechanism

RWD domains mediate protein-protein interactions in a variety of pathways in eukaryotes. In budding yeast, the RWD domain protein Csm1 is particularly versatile, assembling key complexes in the nucleolus and at meiotic kinetochores through multiple protein interaction surfaces. Here, we reveal a third functional context for Csm1 by identifying a new Csm1-interacting protein, Dse3. We show that Dse3 interacts with Csm1 in a structurally equivalent manner to its known binding partners Mam1 and Ulp2, despite these three proteins' lack of overall sequence homology. We theorize that the unique "clamp" structure of Csm1 and the loose sequence requirements for Csm1 binding have led to its incorporation into at least three different structural/signaling pathways in budding yeast.

Modulating transcription factor activity: Interfering with protein-protein interaction networks

Biophysical parameters that govern transcription factors activity are binding locations across the genome, dwelling time at these regulatory elements and specific protein-protein interactions. Most molecular strategies used to develop small compounds that block transcription factors activity have been based on biochemistry and cell biology methods that that do not take into consideration these key biophysical features. Here, we review the advance in the field of transcription factor biology and describe how their interactome and transcriptional regulation on a genome wide scale have been deciphered. We suggest that this new knowledge has the potential to be used to implement innovative research drug discovery program.

Link Prediction based on Quantum-Inspired Ant Colony Optimization

Incomplete or partial observations of network structures pose a serious challenge to theoretical and engineering studies of real networks. To remedy the missing links in real datasets, topology-based link prediction is introduced into the studies of various networks. Due to the complexity of network structures, the accuracy and robustness of most link prediction algorithms are not satisfying enough. In this paper, we propose a quantum-inspired ant colony optimization algorithm that integrates ant colony optimization and quantum computing to predict links in networks. Extensive experiments on both synthetic and real networks show that the accuracy and robustness of the new algorithm is competitive in respect to most of the state of the art algorithms. This result suggests that the application of intelligent optimization to link prediction is promising for boosting its accuracy and robustness.

Whi2 is a conserved negative regulator of TORC1 in response to low amino acids

Yeast WHI2 was originally identified in a genetic screen for regulators of cell cycle arrest and later suggested to function in general stress responses. However, the function of Whi2 is unknown. Whi2 has predicted structure and sequence similarity to human KCTD family proteins, which have been implicated in several cancers and are causally associated with neurological disorders but are largely uncharacterized. The identification of conserved functions between these yeast and human proteins may provide insight into disease mechanisms. We report that yeast WHI2 is a new negative regulator of TORC1 required to suppress TORC1 activity and cell growth specifically in response to low amino acids. In contrast to current opinion, WHI2 is dispensable for TORC1 inhibition in low glucose. The only widely conserved mechanism that actively suppresses both yeast and mammalian TORC1 specifically in response to low amino acids is the conserved SEACIT/GATOR1 complex that inactivates the TORC1-activating RAG-like GTPases. Unexpectedly, Whi2 acts independently and simultaneously with these established GATOR1-like Npr2-Npr3-Iml1 and RAG-like Gtr1-Gtr2 complexes, and also acts independently of the PKA pathway. Instead, Whi2 inhibits TORC1 activity through its binding partners, protein phosphatases Psr1 and Psr2, which were previously thought to only regulate amino acid levels downstream of TORC1. Furthermore, the ability to suppress TORC1 is conserved in the SKP1/BTB/POZ domain-containing, Whi2-like human protein KCTD11 but not other KCTD family members tested.

Integrating genetic and protein-protein interaction networks maps a functional wiring diagram of a cell

Systematic experimental approaches have led to construction of comprehensive genetic and protein-protein interaction networks for the budding yeast, Saccharomyces cerevisiae. Genetic interactions capture functional relationships between genes using phenotypic readouts, while protein-protein interactions identify physical connections between gene products. These complementary, and largely non-overlapping, networks provide a global view of the functional architecture of a cell, revealing general organizing principles, many of which appear to be evolutionarily conserved. Here, we focus on insights derived from the integration of large-scale genetic and protein-protein interaction networks, highlighting principles that apply to both unicellular and more complex systems, including human cells. Network integration reveals fundamental connections involving key functional modules of eukaryotic cells, defining a core network of cellular function, which could be elaborated to explore cell-type specificity in metazoans.

Metabolic Network-Based Identification and Prioritization of Anticancer Targets Based on Expression Data in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a deadly form of liver cancer with high mortality worldwide. Unfortunately, the large heterogeneity of this disease makes it difficult to develop effective treatment strategies. Cellular network analyses have been employed to study heterogeneity in cancer, and to identify potential therapeutic targets. However, the existing approaches do not consider metabolic growth requirements, i.e., biological network functionality, to rank candidate targets while preventing toxicity to non-cancerous tissues. Here, we developed an algorithm to overcome these issues based on integration of gene expression data, genome-scale metabolic models, network controllability, and dispensability, as well as toxicity analysis. This method thus predicts and ranks potential anticancer non-toxic controlling metabolite and gene targets. Our algorithm encompasses both objective-driven and-independent tasks, and uses network topology to finally rank the predicted therapeutic targets. We employed this algorithm to the analysis of transcriptomic data for 50 HCC patients with both cancerous and non-cancerous samples. We identified several potential targets that would prevent cell growth, including 74 anticancer metabolites, and 3 gene targets (PRKACA, PGS1, and CRLS1). The predicted anticancer metabolites showed good agreement with existing FDA-approved cancer drugs, and the 3 genes were experimentally validated by performing experiments in HepG2 and Hep3B liver cancer cell lines. Our observations indicate that our novel approach successfully identifies therapeutic targets for effective treatment of cancer. This approach may also be applied to any cancer type that has tumor and non-tumor gene or protein expression data.

An optimization framework for network annotation

Motivation

A chief goal of systems biology is the reconstruction of large-scale executable models of cellular processes of interest. While accurate continuous models are still beyond reach, a powerful alternative is to learn a logical model of the processes under study, which predicts the logical state of any node of the model as a Boolean function of its incoming nodes. Key to learning such models is the functional annotation of the underlying physical interactions with activation/repression (sign) effects. Such annotations are pretty common for a few well-studied biological pathways.

Results

Here we present a novel optimization framework for large-scale sign annotation that employs different plausible models of signaling and combines them in a rigorous manner. We apply our framework to two large-scale knockout datasets in yeast and evaluate its different components as well as the combined model to predict signs of different subsets of physical interactions. Overall, we obtain an accurate predictor that outperforms previous work by a considerable margin.

Availability and implementation

The code is publicly available at https://github.com/spatkar94/NetworkAnnotation.git.

High-Quality Yeast-2-Hybrid Interaction Network Mapping

In this article, we describe a Y2H interaction mapping protocol for systematically screening medium-to-large collections of open reading frames (ORFs) against each other. The protocol is well suited for analysis of focused networks, for modules of interest, assembling genome-scale interactome maps, and investigating host-microbe interactions. The four-step high-throughput protocol has a demonstrated low false-discovery rate that is essential for producing meaningful network maps. Following the assembly of yeast expression clones from an existing ORFeome collection, we describe in detail the four-step procedure that begins with the primary screen using minipools, followed by secondary verification of primary positives, identification of candidate interaction pairs by sequencing, and a final verification step using fresh inoculants. In addition to this core protocol, aspects of network mapping and quality control will be discussed briefly. © 2018 by John Wiley & Sons, Inc.

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.

An interactome perturbation framework prioritizes damaging missense mutations for developmental disorders

Identifying disease-associated missense mutations remains a challenge, especially in large-scale sequencing studies. Here we establish an experimentally and computationally integrated approach to investigate the functional impact of missense mutations in the context of the human interactome network and test our approach by analyzing ~2,000 de novo missense mutations found in autism subjects and their unaffected siblings. Interaction-disrupting de novo missense mutations are more common in autism probands, principally affect hub proteins, and disrupt a significantly higher fraction of hub interactions than in unaffected siblings. Moreover, they tend to disrupt interactions involving genes previously implicated in autism, providing complementary evidence that strengthens previously identified associations and enhances the discovery of new ones. Importantly, by analyzing de novo missense mutation data from six disorders, we demonstrate that our interactome perturbation approach offers a generalizable framework for identifying and prioritizing missense mutations that contribute to the risk of human disease.

Modeling and simulating networks of interdependent protein interactions

Protein interactions are fundamental building blocks of biochemical reaction systems underlying cellular functions. The complexity and functionality of these systems emerge not only from the protein interactions themselves but also from the dependencies between these interactions, as generated by allosteric effects or mutual exclusion due to steric hindrance. Therefore, formal models for integrating and utilizing information about interaction dependencies are of high interest. Here, we describe an approach for endowing protein networks with interaction dependencies using propositional logic, thereby obtaining constrained protein interaction networks ("constrained networks"). The construction of these networks is based on public interaction databases as well as text-mined information about interaction dependencies. We present an efficient data structure and algorithm to simulate protein complex formation in constrained networks. The efficiency of the model allows fast simulation and facilitates the analysis of many proteins in large networks. In addition, this approach enables the simulation of perturbation effects, such as knockout of single or multiple proteins and changes of protein concentrations. We illustrate how our model can be used to analyze a constrained human adhesome protein network, which is responsible for the formation of diverse and dynamic cell-matrix adhesion sites. By comparing protein complex formation under known interaction dependencies versus without dependencies, we investigate how these dependencies shape the resulting repertoire of protein complexes. Furthermore, our model enables investigating how the interplay of network topology with interaction dependencies influences the propagation of perturbation effects across a large biochemical system. Our simulation software CPINSim (for Constrained Protein Interaction Network Simulator) is available under the MIT license at http://github.com/BiancaStoecker/cpinsim and as a Bioconda package (https://bioconda.github.io).

Recent Trends in Plant Protein Complex Analysis in a Developmental Context

Because virtually all proteins interact with other proteins, studying protein-protein interactions (PPIs) is fundamental in understanding protein function. This is especially true when studying specific developmental processes, in which proteins often make developmental stage- or tissue specific interactions. However, studying these specific PPIs in planta can be challenging. One of the most widely adopted methods to study PPIs in planta is affinity purification coupled to mass spectrometry (AP/MS). Recent developments in the field of mass spectrometry have boosted applications of AP/MS in a developmental context. This review covers two main advancements in the field of affinity purification to study plant developmental processes: increasing the developmental resolution of the harvested tissues and moving from affinity purification to affinity enrichment. Furthermore, we discuss some new affinity purification approaches that have recently emerged and could have a profound impact on the future of protein interactome analysis in plants.

Network inference reveals novel connections in pathways regulating growth and defense in the yeast salt response

Cells respond to stressful conditions by coordinating a complex, multi-faceted response that spans many levels of physiology. Much of the response is coordinated by changes in protein phosphorylation. Although the regulators of transcriptome changes during stress are well characterized in Saccharomyces cerevisiae, the upstream regulatory network controlling protein phosphorylation is less well dissected. Here, we developed a computational approach to infer the signaling network that regulates phosphorylation changes in response to salt stress. We developed an approach to link predicted regulators to groups of likely co-regulated phospho-peptides responding to stress, thereby creating new edges in a background protein interaction network. We then use integer linear programming (ILP) to integrate wild type and mutant phospho-proteomic data and predict the network controlling stress-activated phospho-proteomic changes. The network we inferred predicted new regulatory connections between stress-activated and growth-regulating pathways and suggested mechanisms coordinating metabolism, cell-cycle progression, and growth during stress. We confirmed several network predictions with co-immunoprecipitations coupled with mass-spectrometry protein identification and mutant phospho-proteomic analysis. Results show that the cAMP-phosphodiesterase Pde2 physically interacts with many stress-regulated transcription factors targeted by PKA, and that reduced phosphorylation of those factors during stress requires the Rck2 kinase that we show physically interacts with Pde2. Together, our work shows how a high-quality computational network model can facilitate discovery of new pathway interactions during osmotic stress.

Cells sense and respond to stressful environments by utilizing complex signaling networks that integrate diverse signals to coordinate a multi-faceted physiological response. Much of this response is controlled by post-translational protein phosphorylation. Although many regulators that mediate changes in protein phosphorylation are known, how these regulators inter-connect in a single regulatory network that can transmit cellular signals is not known. It is also unclear how regulators that promote growth and regulators that activate the stress response interconnect to reorganize resource allocation during stress. Here, we developed an integrated experimental and computational workflow to infer the signaling network that regulates phosphorylation changes during osmotic stress in the budding yeast Saccharomyces cerevisiae. The workflow integrates data measuring protein phosphorylation changes in response to osmotic stress with known physical interactions between yeast proteins from large-scale datasets, along with other information about how regulators recognize their targets. The resulting network suggested new signaling connections between regulators and pathways, including those involved in regulating growth and defense, and predicted new regulators involved in stress defense. Our work highlights the power of using network inference to deliver new insight on how cells coordinate a diverse adaptive strategy to stress.

Predicted Arabidopsis Interactome Resource and Gene Set Linkage Analysis: A Transcriptomic Analysis Resource

An advanced functional understanding of omics data is important for elucidating the design logic of physiological processes in plants and effectively controlling desired traits in plants. We present the latest versions of the Predicted Arabidopsis Interactome Resource (PAIR) and of the gene set linkage analysis (GSLA) tool, which enable the interpretation of an observed transcriptomic change (differentially expressed genes [DEGs]) in Arabidopsis (Arabidopsis thaliana) with respect to its functional impact for biological processes. PAIR version 5.0 integrates functional association data between genes in multiple forms and infers 335,301 putative functional interactions. GSLA relies on this high-confidence inferred functional association network to expand our perception of the functional impacts of an observed transcriptomic change. GSLA then interprets the biological significance of the observed DEGs using established biological concepts (annotation terms), describing not only the DEGs themselves but also their potential functional impacts. This unique analytical capability can help researchers gain deeper insights into their experimental results and highlight prospective directions for further investigation. We demonstrate the utility of GSLA with two case studies in which GSLA uncovered how molecular events may have caused physiological changes through their collective functional influence on biological processes. Furthermore, we showed that typical annotation-enrichment tools were unable to produce similar insights to PAIR/GSLA. The PAIR version 5.0-inferred interactome and GSLA Web tool both can be accessed at http://public.synergylab.cn/pair/.

Determining the minimum number of protein-protein interactions required to support known protein complexes

The prediction of protein complexes from protein-protein interactions (PPIs) is a well-studied problem in bioinformatics. However, the currently available PPI data is not enough to describe all known protein complexes. In this paper, we express the problem of determining the minimum number of (additional) required protein-protein interactions as a graph theoretic problem under the constraint that each complex constitutes a connected component in a PPI network. For this problem, we develop two computational methods: one is based on integer linear programming (ILPMinPPI) and the other one is based on an existing greedy-type approximation algorithm (GreedyMinPPI) originally developed in the context of communication and social networks. Since the former method is only applicable to datasets of small size, we apply the latter method to a combination of the CYC2008 protein complex dataset and each of eight PPI datasets (STRING, MINT, BioGRID, IntAct, DIP, BIND, WI-PHI, iRefIndex). The results show that the minimum number of additional required PPIs ranges from 51 (STRING) to 964 (BIND), and that even the four best PPI databases, STRING (51), BioGRID (67), WI-PHI (93) and iRefIndex (85), do not include enough PPIs to form all CYC2008 protein complexes. We also demonstrate that the proposed problem framework and our solutions can enhance the prediction accuracy of existing PPI prediction methods. ILPMinPPI can be freely downloaded from http://sunflower.kuicr.kyoto-u.ac.jp/~nakajima/.

PPI network analyses of human WD40 protein family systematically reveal their tendency to assemble complexes and facilitate the complex predictions

Background

WD40 repeat proteins constitute one of the largest families in eukaryotes, and widely participate in various fundamental cellular processes by interacting with other molecules. Based on individual WD40 proteins, previous work has demonstrated that their structural characteristics should confer great potential of interaction and complex formation, and has speculated that they may serve as hubs in the protein-protein interaction (PPI) network. However, what roles the whole family plays in organizing the PPI network, and whether this information can be utilized in complex prediction remain unclear. To address these issues, quantitative and systematic analyses of WD40 proteins from the perspective of PPI networks are highly required.

Results

In this work, we built two human PPI networks by using data sets with different confidence levels, and studied the network properties of the whole human WD40 protein family systematically. Our analyses have quantitatively confirmed that the human WD40 protein family, as a whole, tends to be hubs with an odds ratio of about 1.8 or greater, and the network decomposition has revealed that they are prone to enrich near the global center of the whole network with a fold change of two in the median k-values. By integrating expression profiles, we have further shown that WD40 hub proteins are inclined to be intramodular, which is indicative of complex assembling. Based on this information, we have further predicted 1674 potential WD40-associated complexes by choosing a clique-based method, which is more sensitive than others, and an indirect evaluation by co-expression scores has demonstrated its reliability.

Conclusions

At the systems level but not sporadic examples’ level, this work has provided rich knowledge for better understanding WD40 proteins’ roles in organizing the PPI network. These findings and predicted complexes can offer valuable clues for prioritizing candidates for further studies.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0567-9) contains supplementary material, which is available to authorized users.

Systematic analysis of complex genetic interactions

To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double-mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes, and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and displayed a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance, including the genotype-to-phenotype relationship.

The Interactomic Analysis Reveals Pathogenic Protein Networks in Phomopsis longicolla Underlying Seed Decay of Soybean

Phomopsis longicolla T. W. Hobbs (syn. Diaporthe longicolla) is the primary cause of Phomopsis seed decay (PSD) in soybean, Glycine max (L.) Merrill. This disease results in poor seed quality and is one of the most economically important seed diseases in soybean. The objectives of this study were to infer protein-protein interactions (PPI) and to identify conserved global networks and pathogenicity subnetworks in P. longicolla including orthologous pathways for cell signaling and pathogenesis. The interlog method used in the study identified 215,255 unique PPIs among 3,868 proteins. There were 1,414 pathogenicity related genes in P. longicolla identified using the pathogen host interaction (PHI) database. Additionally, 149 plant cell wall degrading enzymes (PCWDE) were detected. The network captured five different classes of carbohydrate degrading enzymes, including the auxiliary activities, carbohydrate esterases, glycoside hydrolases, glycosyl transferases, and carbohydrate binding molecules. From the PPI analysis, novel interacting partners were determined for each of the PCWDE classes. The most predominant class of PCWDE was a group of 60 glycoside hydrolases proteins. The glycoside hydrolase subnetwork was found to be interacting with 1,442 proteins within the network and was among the largest clusters. The orthologous proteins FUS3, HOG, CYP1, SGE1, and the g5566t.1 gene identified in this study could play an important role in pathogenicity. Therefore, the P. longicolla protein interactome (PiPhom) generated in this study can lead to a better understanding of PPIs in soybean pathogens. Furthermore, the PPI may aid in targeting of genes and proteins for further studies of the pathogenicity mechanisms.

The life of U6 small nuclear RNA, from cradle to grave

Removal of introns from precursor messenger RNA (pre-mRNA) and some noncoding transcripts is an essential step in eukaryotic gene expression. In the nucleus, this process of RNA splicing is carried out by the spliceosome, a multi-megaDalton macromolecular machine whose core components are conserved from yeast to humans. In addition to many proteins, the spliceosome contains five uridine-rich small nuclear RNAs (snRNAs) that undergo an elaborate series of conformational changes to correctly recognize the splice sites and catalyze intron removal. Decades of biochemical and genetic data, along with recent cryo-EM structures, unequivocally demonstrate that U6 snRNA forms much of the catalytic core of the spliceosome and is highly dynamic, interacting with three snRNAs, the pre-mRNA substrate, and >25 protein partners throughout the splicing cycle. This review summarizes the current state of knowledge on how U6 snRNA is synthesized, modified, incorporated into snRNPs and spliceosomes, recycled, and degraded.

SMILE: a novel procedure for subcellular module identification with localisation expansion

Computational clustering methods help identify functional modules in protein–protein interaction (PPI) network, in which proteins participate in the same biological pathways or specific functions. Subcellular localisation is crucial for proteins to implement biological functions and each compartment accommodates specific portions of the protein interaction structure. However, the importance of protein subcellular localisation is often neglected in the studies of module identification. In this study, the authors propose a novel procedure, subcellular module identification with localisation expansion (SMILE), to identify super modules that consist of several subcellular modules performing specific biological functions among cell compartments. These super modules identified by SMILE are more functionally diverse and have been verified to be more associated with known protein complexes and biological pathways compared with the modules identified from the global PPI networks in both the compartmentalised PPI and InWeb_InBioMap datasets. The authors’ results reveal that subcellular localisation is a principal feature of functional modules and offers important guidance in detecting biologically meaningful results.

Yeast Aim21/Tda2 both regulates free actin by reducing barbed end assembly and forms a complex with Cap1/Cap2 to balance actin assembly between patches and cables

Yeast Aim21 is recruited by the SH3-containing proteins Bbc1 and Abp1 to patches and, with Tda2, reduces barbed end assembly to balance the distribution of actin between patches and cables. Aim21/Tda2 also interacts with Cap1/Cap2, revealing a complex interplay between actin assembly regulators.

How cells balance the incorporation of actin into diverse structures is poorly understood. In budding yeast, a single actin monomer pool is used to build both actin cables involved in polarized growth and actin cortical patches involved in endocytosis. Here we report how Aim21/Tda2 is recruited to the cortical region of actin patches, where it negatively regulates actin assembly to elevate the available actin monomer pool. Aim21 has four polyproline regions and is recruited by two SH3-containing patch proteins, Bbc1 and Abp1. The C-terminal region, which is required for its function, binds Tda2. Cell biological and biochemical data reveal that Aim21/Tda2 is a negative regulator of barbed end filamentous actin (F-actin) assembly, and this activity is necessary for efficient endocytosis and plays a pivotal role in balancing the distribution of actin between cables and patches. Aim21/Tda2 also forms a complex with the F-actin barbed end capping protein Cap1/Cap2, revealing an interplay between regulators and showing the complexity of regulation of barbed end assembly.

Hub Protein Controversy: Taking a Closer Look at Plant Stress Response Hubs

Plant stress responses involve numerous changes at the molecular and cellular level and are regulated by highly complex signaling pathways. Studying protein-protein interactions (PPIs) and the resulting networks is therefore becoming increasingly important in understanding these responses. Crucial in PPI networks are the so-called hubs or hub proteins, commonly defined as the most highly connected central proteins in scale-free PPI networks. However, despite their importance, a growing amount of confusion and controversy seems to exist regarding hub protein identification, characterization and classification. In order to highlight these inconsistencies and stimulate further clarification, this review critically analyses the current knowledge on hub proteins in the plant interactome field. We focus on current hub protein definitions, including the properties generally seen as hub-defining, and the challenges and approaches associated with hub protein identification. Furthermore, we give an overview of the most important large-scale plant PPI studies of the last decade that identified hub proteins, pointing out the lack of overlap between different studies. As such, it appears that although major advances are being made in the plant interactome field, defining hub proteins is still heavily dependent on the quality, origin and interpretation of the acquired PPI data. Nevertheless, many hub proteins seem to have a reported role in the plant stress response, including transcription factors, protein kinases and phosphatases, ubiquitin proteasome system related proteins, (co-)chaperones and redox signaling proteins. A significant number of identified plant stress hubs are however still functionally uncharacterized, making them interesting targets for future research. This review clearly shows the ongoing improvements in the plant interactome field but also calls attention to the need for a more comprehensive and precise identification of hub proteins, allowing a more efficient systems biology driven unraveling of complex processes, including those involved in stress responses.

Interactome INSIDER: a structural interactome browser for genomic studies

We present Interactome INSIDER, a tool to link genomic variant information with structural protein-protein interactomes. Underlying this tool is the application of machine learning to predict protein interaction interfaces for 185,957 protein interactions with previously unresolved interfaces, in human and 7 model organisms, including the entire experimentally determined human binary interactome. Predicted interfaces exhibit similar functional properties as known interfaces, including enrichment for disease mutations and recurrent cancer mutations. Through 2,164 de novo mutagenesis experiments, we show that mutations of predicted and known interface residues disrupt interactions at a similar rate, and much more frequently than mutations outside of predicted interfaces. To spur functional genomic studies, Interactome INSIDER (http://interactomeinsider.yulab.org) enables users to identify whether variants or disease mutations are enriched in known and predicted interaction interfaces at various resolutions. Users may explore known population variants, disease mutations, and somatic cancer mutations, or upload their own set of mutations for this purpose.

Meloidogyne incognita PASSE-MURAILLE (MiPM) Gene Encodes a Cell-Penetrating Protein That Interacts With the CSN5 Subunit of the COP9 Signalosome

The pathogenicity of phytonematodes relies on secreted virulence factors to rewire host cellular pathways for the benefits of the nematode. In the root-knot nematode (RKN) Meloidogyne incognita, thousands of predicted secreted proteins have been identified and are expected to interact with host proteins at different developmental stages of the parasite. Identifying the host targets will provide compelling evidence about the biological significance and molecular function of the predicted proteins. Here, we have focused on the hub protein CSN5, the fifth subunit of the pleiotropic and eukaryotic conserved COP9 signalosome (CSN), which is a regulatory component of the ubiquitin/proteasome system. We used affinity purification-mass spectrometry (AP-MS) to generate the interaction network of CSN5 in M. incognita-infected roots. We identified the complete CSN complex and other known CSN5 interaction partners in addition to unknown plant and M. incognita proteins. Among these, we described M. incognita PASSE-MURAILLE (MiPM), a small pioneer protein predicted to contain a secretory peptide that is up-regulated mostly in the J2 parasitic stage. We confirmed the CSN5-MiPM interaction, which occurs in the nucleus, by bimolecular fluorescence complementation (BiFC). Using MiPM as bait, a GST pull-down assay coupled with MS revealed some common protein partners between CSN5 and MiPM. We further showed by in silico and microscopic analyses that the recombinant purified MiPM protein enters the cells of Arabidopsis root tips in a non-infectious context. In further detail, the supercharged N-terminal tail of MiPM (NTT-MiPM) triggers an unknown host endocytosis pathway to penetrate the cell. The functional meaning of the CSN5-MiPM interaction in the M. incognita parasitism is discussed. Moreover, we propose that the cell-penetrating properties of some M. incognita secreted proteins might be a non-negligible mechanism for cell uptake, especially during the steps preceding the sedentary parasitic phase.

The Notch Interactome: Complexity in Signaling Circuitry

The Notch pathway controls a very broad spectrum of cell fates in metazoans during development, influencing proliferation, differentiation and cell death. Given its central role in normal development and homeostasis, misregulation of Notch signals can lead to various disorders including cancer. How the Notch pathway mediates such pleiotropic and differential effects is of fundamental importance. It is becoming increasingly clear through a number of large-scale genetic and proteomic studies that Notch interacts with a staggeringly large number of other genes and pathways in a context-dependent, complex, and highly regulated network, which determines the ultimate biological outcome. How best to interpret and analyze the continuously increasing wealth of data on Notch interactors remains a challenge. Here we review the current state of genetic and proteomic data related to the Notch interactome.

Dynamic patterns of information flow in complex networks

Although networks are extensively used to visualize information flow in biological, social and technological systems, translating topology into dynamic flow continues to challenge us, as similar networks exhibit fundamentally different flow patterns, driven by different interaction mechanisms. To uncover a network’s actual flow patterns, here we use a perturbative formalism, analytically tracking the contribution of all nodes/paths to the flow of information, exposing the rules that link structure and dynamic information flow for a broad range of nonlinear systems. We find that the diversity of flow patterns can be mapped into a single universal function, characterizing the interplay between the system’s topology and its dynamics, ultimately allowing us to identify the network’s main arteries of information flow. Counter-intuitively, our formalism predicts a family of frequently encountered dynamics where the flow of information avoids the hubs, favoring the network’s peripheral pathways, a striking disparity between structure and dynamics.

Complex networks are a useful tool to investigate spreading processes but topology alone is insufficient to predict information flow. Here the authors propose a measure of information flow and predict its behavior from the interplay between structure and dynamics.

RedNemo: topology-based PPI network reconstruction via repeated diffusion with neighborhood modifications

Motivation

Analysis of protein-protein interaction (PPI) networks provides invaluable insight into several systems biology problems. High-throughput experimental techniques together with computational methods provide large-scale PPI networks. However, a major issue with these networks is their erroneous nature; they contain false-positive interactions and usually many more false-negatives. Recently, several computational methods have been proposed for network reconstruction based on topology, where given an input PPI network the goal is to reconstruct the network by identifying false-positives/-negatives as correctly as possible.

Results

We observe that the existing topology-based network reconstruction algorithms suffer several shortcomings. An important issue is regarding the scalability of their computational requirements, especially in terms of execution times, with the network sizes. They have only been tested on small-scale networks thus far and when applied on large-scale networks of popular PPI databases, the executions require unreasonable amounts of time, or may even crash without producing any output for some instances even after several months of execution. We provide an algorithm, RedNemo, for the topology-based network reconstruction problem. It provides more accurate networks than the alternatives as far as biological qualities measured in terms of most metrics based on gene ontology annotations. The recovery of a high-confidence network modified via random edge removals and rewirings is also better with RedNemo than with the alternatives under most of the experimented removal/rewiring ratios. Furthermore, through extensive tests on databases of varying sizes, we show that RedNemo achieves these results with much better running time performances.

Availability and Implementation

Supplementary material including source code, useful scripts, experimental data and the results are available at http://webprs.khas.edu.tr/~cesim/RedNemo.tar.gz.

Contact

cesim@khas.edu.tr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mechanisms of action of Coxiella burnetii effectors inferred from host-pathogen protein interactions

Coxiella burnetii is an obligate Gram-negative intracellular pathogen and the etiological agent of Q fever. Successful infection requires a functional Type IV secretion system, which translocates more than 100 effector proteins into the host cytosol to establish the infection, restructure the intracellular host environment, and create a parasitophorous vacuole where the replicating bacteria reside. We used yeast two-hybrid (Y2H) screening of 33 selected C. burnetii effectors against whole genome human and murine proteome libraries to generate a map of potential host-pathogen protein-protein interactions (PPIs). We detected 273 unique interactions between 20 pathogen and 247 human proteins, and 157 between 17 pathogen and 137 murine proteins. We used orthology to combine the data and create a single host-pathogen interaction network containing 415 unique interactions between 25 C. burnetii and 363 human proteins. We further performed complementary pairwise Y2H testing of 43 out of 91 C. burnetii-human interactions involving five pathogen proteins. We used the combined data to 1) perform enrichment analyses of target host cellular processes and pathways, 2) examine effectors with known infection phenotypes, and 3) infer potential mechanisms of action for four effectors with uncharacterized functions. The host-pathogen interaction profiles supported known Coxiella phenotypes, such as adapting cell morphology through cytoskeletal re-arrangements, protein processing and trafficking, organelle generation, cholesterol processing, innate immune modulation, and interactions with the ubiquitin and proteasome pathways. The generated dataset of PPIs-the largest collection of unbiased Coxiella host-pathogen interactions to date-represents a rich source of information with respect to secreted pathogen effector proteins and their interactions with human host proteins.

Identification of genes and critical control proteins associated with inflammatory breast cancer using network controllability

One of the most aggressive forms of breast cancer is inflammatory breast cancer (IBC), whose lack of tumour mass also makes a prompt diagnosis difficult. Moreover, genomic differences between common breast cancers and IBC have not been completely assessed, thus substantially limiting the identification of biomarkers unique to IBC. Here, we developed a novel statistical analysis of gene expression profiles corresponding to microdissected IBC, non-IBC (nIBC) and normal samples that enabled us to identify a set of genes significantly associated with a specific disease state. Second, by using advanced methods based on controllability network theory, we identified a set of critical control proteins that uniquely and structurally control the entire proteome. By mapping high change variance genes in protein interaction networks, we found that a large statistically significant fraction of genes whose variance changed significantly between normal and IBC and nIBC disease states were among the set of critical control proteins. Moreover, this analysis identified the overlapping genes with the highest statistical significance; these genes may assist in developing future biomarkers and determining drug targets to disrupt the molecular pathways driving carcinogenesis in IBC.

Interactome evolution: insights from genome-wide analyses of protein-protein interactions

We highlight new evolutionary insights enabled by recent genome-wide studies on protein-protein interaction (PPI) networks ('interactomes'). While most PPIs are mediated by a single sequence region promoting or inhibiting interactions, many PPIs are mediated by multiple sequence regions acting cooperatively. Most PPIs perform important functions maintained by negative selection: we estimate that less than ∼10% of the human interactome is effectively neutral upon perturbation (i.e. 'junk' PPIs), and the rest are deleterious upon perturbation; interfacial sites evolve more slowly than other sites; many conserved PPIs show signatures of co-evolution at the interface; PPIs evolve more slowly than protein sequence. At the same time, many PPIs undergo rewiring during evolution for lineage-specific adaptation. Finally, chaperone-protein and host-pathogen interactomes are governed by distinct evolutionary principles.

MCC/Eisosomes Regulate Cell Wall Synthesis and Stress Responses in Fungi

The fungal plasma membrane is critical for cell wall synthesis and other important processes including nutrient uptake, secretion, endocytosis, morphogenesis, and response to stress. To coordinate these diverse functions, the plasma membrane is organized into specialized compartments that vary in size, stability, and composition. One recently identified domain known as the Membrane Compartment of Can1 (MCC)/eisosome is distinctive in that it corresponds to a furrow-like invagination in the plasma membrane. MCC/eisosomes have been shown to be formed by the Bin/Amphiphysin/Rvs (BAR) domain proteins Lsp1 and Pil1 in a range of fungi. MCC/eisosome domains influence multiple cellular functions; but a very pronounced defect in cell wall synthesis has been observed for mutants with defects in MCC/eisosomes in some yeast species. For example, Candida albicans MCC/eisosome mutants display abnormal spatial regulation of cell wall synthesis, including large invaginations and altered chemical composition of the walls. Recent studies indicate that MCC/eisosomes affect cell wall synthesis in part by regulating the levels of the key regulatory lipid phosphatidylinositol 4,5-bisphosphate (PI_4,5P₂) in the plasma membrane. One general way MCC/eisosomes function is by acting as protected islands in the plasma membrane, since these domains are very stable. They also act as scaffolds to recruit >20 proteins. Genetic studies aimed at defining the function of the MCC/eisosome proteins have identified important roles in resistance to stress, such as resistance to oxidative stress mediated by the flavodoxin-like proteins Pst1, Pst2, Pst3 and Ycp4. Thus, MCC/eisosomes play multiple roles in plasma membrane organization that protect fungal cells from the environment.

Numerous interactions act redundantly to assemble a tunable size of P bodies in Saccharomyces cerevisiae

Eukaryotic cells contain multiple RNA-protein assemblies referred to as RNP granules, which are thought to form through multiple protein-protein interactions analogous to a liquid-liquid phase separation. One class of RNP granules consists of P bodies, which consist of nontranslating mRNAs and the general translation repression and mRNA degradation machinery. P bodies have been suggested to form predominantly through interactions of Edc3 and a prion-like domain on Lsm4. In this work, we provide evidence that P-body assembly can be driven by multiple different protein-protein and/or protein-RNA interactions, including interactions involving Dhh1, Psp2, and Pby1. Moreover, the relative importance of specific interactions can vary with different growth conditions. Based on these observations, we develop a summative model wherein the P-body assembly phenotype of a given mutant can be predicted from the number of currently known protein-protein interactions between P-body components.

Emerging and evolving concepts in gene essentiality

Gene essentiality is a founding concept of genetics with important implications in both fundamental and applied research. Multiple screens have been performed over the years in bacteria, yeasts, animals and more recently in human cells to identify essential genes. A mounting body of evidence suggests that gene essentiality, rather than being a static and binary property, is both context dependent and evolvable in all kingdoms of life. This concept of a non-absolute nature of gene essentiality changes our fundamental understanding of essential biological processes and could directly affect future treatment strategies for cancer and infectious diseases.

Systematic Identification of Mycobacterium tuberculosis Effectors Reveals that BfrB Suppresses Innate Immunity

Mycobacterium tuberculosis (Mtb) has evolved multiple strategies to counter the human immune system. The effectors of Mtb play important roles in the interactions with the host. However, because of the lack of highly efficient strategies, there are only a handful of known Mtb effectors, thus hampering our understanding of Mtb pathogenesis. In this study, we probed Mtb proteome microarray with biotinylated whole-cell lysates of human macrophages, identifying 26 Mtb membrane proteins and secreted proteins that bind to macrophage proteins. Combining GST pull-down with mass spectroscopy then enabled the specific identification of all binders. We refer to this proteome microarray-based strategy as SOPHIE (Systematic unlOcking of Pathogen and Host Interacting Effectors). Detailed investigation of a novel effector identified here, the iron storage protein BfrB (Rv3841), revealed that BfrB inhibits NF-κB-dependent transcription through binding and reducing the nuclear abundance of the ribosomal protein S3 (RPS3), which is a functional subunit of NF- κB. The importance of this interaction was evidenced by the promotion of survival in macrophages of the mycobacteria, Mycobacterium smegmatis, by overexpression of BfrB. Thus, beyond demonstrating the power of SOPHIE in the discovery of novel effectors of human pathogens, we expect that the set of Mtb effectors identified in this work will greatly facilitate the understanding of the pathogenesis of Mtb, possibly leading to additional potential molecular targets in the battle against tuberculosis.

A DNA nanoscope via auto-cycling proximity recording

Analysis of the spatial arrangement of molecular features enables the engineering of synthetic nanostructures and the understanding of natural ones. The ability to acquire a comprehensive set of pairwise proximities between components would satisfy an increasing interest in investigating individual macromolecules and their interactions, but current biochemical techniques detect only a single proximity partner per probe. Here, we present a biochemical DNA nanoscopy method that records nanostructure features in situ and in detail for later readout. Based on a conceptually novel auto-cycling proximity recording (APR) mechanism, it continuously and repeatedly produces proximity records of any nearby pairs of DNA-barcoded probes, at physiological temperature, without altering the probes themselves. We demonstrate the production of dozens of records per probe, decode the spatial arrangements of 7 unique probes in a homogeneous sample, and repeatedly sample the same probes in different states.

The spatial organisation of nanostructures is fundamental to their function. Here, the authors develop a non-destructive, proximity-based method to record extensive spatial organization information in DNA molecules for later readout.