The role of protein interactions in mediating essentiality and synthetic lethality

MAGICAL: A multi-class classifier to predict synthetic lethal and viable interactions using protein-protein interaction network

Synthetic lethality (SL) and synthetic viability (SV) are commonly studied genetic interactions in the targeted therapy approach in cancer. In SL, inhibiting either of the genes does not affect the cancer cell survival, but inhibiting both leads to a lethal phenotype. In SV, inhibiting the vulnerable gene makes the cancer cell sick; inhibiting the partner gene rescues and promotes cell viability. Many low and high-throughput experimental approaches have been employed to identify SLs and SVs, but they are time-consuming and expensive. The computational tools for SL prediction involve statistical and machine-learning approaches. Almost all machine learning tools are binary classifiers and involve only identifying SL pairs. Most importantly, there are limited properties known that best describe and discriminate SL from SV. We developed MAGICAL (Multi-class Approach for Genetic Interaction in Cancer via Algorithm Learning), a multi-class random forest based machine learning model for genetic interaction prediction. Network properties of protein derived from physical protein-protein interactions are used as features to classify SL and SV. The model results in an accuracy of ~80% for the training dataset (CGIdb, BioGRID, and SynLethDB) and performs well on DepMap and other experimentally derived reported datasets. Amongst all the network properties, the shortest path, average neighbor2, average betweenness, average triangle, and adhesion have significant discriminatory power. MAGICAL is the first multi-class model to identify discriminatory features of synthetic lethal and viable interactions. MAGICAL can predict SL and SV interactions with better accuracy and precision than any existing binary classifier.

Synthetic negative genome screen of the GPN-loop GTPase NPA3 in Saccharomyces cerevisiae

The GPN-loop GTPase Npa3 is encoded by an essential gene in the yeast Saccharomyces cerevisiae. Npa3 plays a critical role in the assembly and nuclear accumulation of RNA polymerase II (RNAPII), a function that may explain its essentiality. Genetic interactions describe the extent to which a mutation in a particular gene affects a specific phenotype when co-occurring with an alteration in a second gene. Discovering synthetic negative genetic interactions has long been used as a tool to delineate the functional relatedness between pairs of genes participating in common or compensatory biological pathways. Previously, our group showed that nuclear targeting and transcriptional activity of RNAPII were unaffected in cells expressing exclusively a C-terminal truncated mutant version of Npa3 (npa3∆C) lacking the last 106 residues naturally absent from the single GPN protein in Archaea, but universally conserved in all Npa3 orthologs of eukaryotes. To gain insight into novel cellular functions for Npa3, we performed here a genome-wide Synthetic Genetic Array (SGA) study coupled to bulk fluorescence monitoring to identify negative genetic interactions of NPA3 by crossing an npa3∆C strain with a 4,389 nonessential gene-deletion collection. This genetic screen revealed previously unknown synthetic negative interactions between NPA3 and 15 genes. Our results revealed that the Npa3 C-terminal tail extension regulates the participation of this essential GTPase in previously unknown biological processes related to mitochondrial homeostasis and ribosome biogenesis.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Archetypal transcriptional blocks underpin yeast gene regulation in response to changes in growth conditions

The transcriptional responses of yeast cells to diverse stresses typically include gene activation and repression. Specific stress defense, citric acid cycle and oxidative phosphorylation genes are activated, whereas protein synthesis genes are coordinately repressed. This view was achieved from comparative transcriptomic experiments delineating sets of genes whose expression greatly changed with specific stresses. Less attention has been paid to the biological significance of 1) consistent, albeit modest, changes in RNA levels across multiple conditions, and 2) the global gene expression correlations observed when comparing numerous genome-wide studies. To address this, we performed a meta-analysis of 1379 microarray-based experiments in yeast, and identified 1388 blocks of RNAs whose expression changes correlate across multiple and diverse conditions. Many of these blocks represent sets of functionally-related RNAs that act in a coordinated fashion under normal and stress conditions, and map to global cell defense and growth responses. Subsequently, we used the blocks to analyze novel RNA-seq experiments, demonstrating their utility and confirming the conclusions drawn from the meta-analysis. Our results provide a new framework for understanding the biological significance of changes in gene expression: 'archetypal' transcriptional blocks that are regulated in a concerted fashion in response to external stimuli.

Binding interface change and cryptic variation in the evolution of protein-protein interactions

Background

Physical interactions between proteins are essential for almost all biological functions and systems. To understand the evolution of function it is therefore important to understand the evolution of molecular interactions. Of key importance is the evolution of binding specificity, the set of interactions made by a protein, since change in specificity can lead to “rewiring” of interaction networks. Unfortunately, the interfaces through which proteins interact are complex, typically containing many amino-acid residues that collectively must contribute to binding specificity as well as binding affinity, structural integrity of the interface and solubility in the unbound state.

Results

In order to study the relationship between interface composition and binding specificity, we make use of paralogous pairs of yeast proteins. Immediately after duplication these paralogues will have identical sequences and protein products that make an identical set of interactions. As the sequences diverge, we can correlate amino-acid change in the interface with any change in the specificity of binding. We show that change in interface regions correlates only weakly with change in specificity, and many variants in interfaces are functionally equivalent. We show that many of the residue replacements within interfaces are silent with respect to their contribution to binding specificity.

Conclusions

We conclude that such functionally-equivalent change has the potential to contribute to evolutionary plasticity in interfaces by creating cryptic variation, which in turn may provide the raw material for functional innovation and coevolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0608-1) contains supplementary material, which is available to authorized users.

UDoNC: An Algorithm for Identifying Essential Proteins Based on Protein Domains and Protein-Protein Interaction Networks

Prediction of essential proteins which are crucial to an organism's survival is important for disease analysis and drug design, as well as the understanding of cellular life. The majority of prediction methods infer the possibility of proteins to be essential by using the network topology. However, these methods are limited to the completeness of available protein-protein interaction (PPI) data and depend on the network accuracy. To overcome these limitations, some computational methods have been proposed. However, seldom of them solve this problem by taking consideration of protein domains. In this work, we first analyze the correlation between the essentiality of proteins and their domain features based on data of 13 species. We find that the proteins containing more protein domain types which rarely occur in other proteins tend to be essential. Accordingly, we propose a new prediction method, named UDoNC, by combining the domain features of proteins with their topological properties in PPI network. In UDoNC, the essentiality of proteins is decided by the number and the frequency of their protein domain types, as well as the essentiality of their adjacent edges measured by edge clustering coefficient. The experimental results on S. cerevisiae data show that UDoNC outperforms other existing methods in terms of area under the curve (AUC). Additionally, UDoNC can also perform well in predicting essential proteins on data of E. coli.