Identifying Pathway Proteins in Networks using Convergence

One of the key goals of systems biology concerns the analysis of experimental biological data available to the scientific public. New technologies are rapidly developed to observe and report whole-scale biological phenomena; however, few methods exist with the ability to produce specific, testable hypotheses from this noisy ‘big’ data. In this work, we propose an approach that combines the power of data-driven network theory along with knowledge-based ontology to tackle this problem. Network models are especially powerful due to their ability to display elements of interest and their relationships as internetwork structures. Additionally, ontological data actually supplements the confidence of relationships within the model without clouding critical structure identification. As such, we postulate that given a (gene/protein) marker set of interest, we can systematically identify the core of their interactions (if they are indeed working together toward a biological function), via elimination of original markers and addition of additional necessary markers. This concept, which we refer to as “convergence,” harnesses the idea of “guilt-by-association” and recursion to identify whether a core of relationships exists between markers. In this study, we test graph theoretic concepts such as shortest-path, k-Nearest- Neighbor and clustering) to identify cores iteratively in data- and knowledge-based networks in the canonical yeast Pheromone Mating Response pathway. Additionally, we provide results for convergence application in virus infection, hearing loss, and Parkinson’s disease. Our results indicate that if a marker set has common discrete function, this approach is able to identify that function, its interacting markers, and any new elements necessary to complete the structural core of that function. The result below find that the shortest path function is the best approach of those used, finding small target sets that contain a majority or all of the markers in the gold standard pathway. The power of this approach lies in its ability to be used in investigative studies to inform decisions concerning target selection

Functional Characteristics of a Highly Specific Integrase Encoded by an LTR-Retrotransposon

Background: The retroviral Integrase protein catalyzes the insertion of linear viral DNA into host cell DNA. Although different retroviruses have been shown to target distinctive chromosomal regions, few of them display a site-specific integration. ZAM, a retroelement from Drosophila melanogaster very similar in structure and replication cycle to mammalian retroviruses is highly site-specific. Indeed, ZAM copies target the genomic 59-CGCGCg-39 consensus-sequences. To enlighten the determinants of this high integration specificity, we investigated the functional properties of its integrase protein denoted ZAM-IN. Principal Findings: Here we show that ZAM-IN displays the property to nick DNA molecules in vitro. This endonuclease activity targets specific sequences that are present in a 388 bp fragment taken from the white locus and known to be a genomic ZAM integration site in vivo. Furthermore, ZAM-IN displays the unusual property to directly bind specific genomic DNA sequences. Two specific and independent sites are recognized within the 388 bp fragment of the white locus: the CGCGCg sequence and a closely apposed site different in sequence. Conclusion: This study strongly argues that the intrinsic properties of ZAM-IN, ie its binding properties and its endonuclease activity, play an important part in ZAM integration specificity. Its ability to select two binding sites and to nick the DNA molecule reminds the strategy used by some site-specific recombination enzymes and forms the basis for site-specifi

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome

Retroviral insertions in the VISION database identify molecular pathways in mouse lymphoid leukemia and lymphoma

AKXD recombinant inbred (RI) strains develop a variety of leukemias and lymphomas due to somatically acquired insertions of retroviral DNA into the genome of hematopoetic cells that can mutate cellular proto-oncogenes and tumor suppressor genes. We generated a new set of tumors from nine AKXD RI strains selected for their propensity to develop B-cell tumors, the most common type of human hematopoietic cancers. We employed a PCR technique called viral insertion site amplification (VISA) to rapidly isolate genomic sequence at the site of provirus insertion. Here we describe 550 VISA sequence tags (VSTs) that identify 74 common insertion sites (CISs), of which 21 have not been identified previously. Several suspected proto-oncogenes and tumor suppressor genes lie near CISs, providing supportive evidence for their roles in cancer. Furthermore, numerous previously uncharacterized genes lie near CISs, providing a pool of candidate disease genes for future research. Pathway analysis of candidate genes identified several signaling pathways as common and powerful routes to blood cancer, including Notch, E-protein, NFκB, and Ras signaling. Misregulation of several Notch signaling genes was confirmed by quantitative RT-PCR. Our data suggest that analyses of insertional mutagenesis on a single genetic background are biased toward the identification of cooperating mutations. This tumor collection represents the most comprehensive study of the genetics of B-cell leukemia and lymphoma development in mice. We have deposited the VST sequences, CISs in a genome viewer, histopathology, and molecular tumor typing data in a public web database called VISION (Viral Insertion Sites Identifying Oncogenes), which is located at http://www.mouse-genome.bcm.tmc.edu/vision

Synthetic Morphology Using Alternative Inputs

Designing the shape and size of a cell is an interesting challenge for synthetic biology. Prolonged exposure to the mating pheromone α-factor induces an unusual morphology in yeast cells: multiple mating projections. The goal of this work was to reproduce the multiple projections phenotype in the absence of α-factor using a gain-of-function approach termed “Alternative Inputs (AIs)”. An alternative input is defined as any genetic manipulation that can activate the signaling pathway instead of the natural input. Interestingly, none of the alternative inputs were sufficient to produce multiple projections although some produced a single projection. Then, we extended our search by creating all combinations of alternative inputs and deletions that were summarized in an AIs-Deletions matrix. We found a genetic manipulation (AI-Ste5p ste2Δ) that enhanced the formation of multiple projections. Following up this lead, we demonstrated that AI-Ste4p and AI-Ste5p were sufficient to produce multiple projections when combined. Further, we showed that overexpression of a membrane-targeted form of Ste5p alone could also induce multiple projections. Thus, we successfully re-engineered the multiple projections mating morphology using alternative inputs without α-factor

The Ankyrin Repeats and DHHC S-acyl Transferase Domain of AKR1 Act Independently to Regulate Switching from Vegetative to Mating States in Yeast

Signal transduction from G-protein coupled receptors to MAPK cascades through heterotrimeric G-proteins has been described for many eukaryotic systems. One of the best-characterised examples is the yeast pheromone response pathway, which is negatively regulated by AKR1. AKR1-like proteins are present in all eukaryotes and contain a DHHC domain and six ankyrin repeats. Whilst the DHHC domain dependant S-acyl transferase (palmitoyl transferase) function of AKR1 is well documented it is not known whether the ankyrin repeats are also required for this activity. Here we show that the ankyrin repeats of AKR1 are required for full suppression of the yeast pheromone response pathway, by sequestration of the Gβγ dimer, and act independently of AKR1 S-acylation function. Importantly, the functions provided by the AKR1 ankyrin repeats and DHHC domain are not required on the same molecule to fully restore WT phenotypes and function. We also show that AKR1 molecules are S-acylated at locations other than the DHHC cysteine, increasing the abundance of AKR1 in the cell. Our results have important consequences for studies of AKR1 function, including recent attempts to characterise S-acylation enzymology and kinetics. Proteins similar to AKR1 are found in all eukaryotes and our results have broad implications for future work on these proteins and the control of switching between Gβγ regulated pathways

A Man-Made ATP-Binding Protein Evolved Independent of Nature Causes Abnormal Growth in Bacterial Cells

Recent advances in de novo protein evolution have made it possible to create synthetic proteins from unbiased libraries that fold into stable tertiary structures with predefined functions. However, it is not known whether such proteins will be functional when expressed inside living cells or how a host organism would respond to an encounter with a non-biological protein. Here, we examine the physiology and morphology of Escherichia coli cells engineered to express a synthetic ATP-binding protein evolved entirely from non-biological origins. We show that this man-made protein disrupts the normal energetic balance of the cell by altering the levels of intracellular ATP. This disruption cascades into a series of events that ultimately limit reproductive competency by inhibiting cell division. We now describe a detailed investigation into the synthetic biology of this man-made protein in a living bacterial organism, and the effect that this protein has on normal cell physiology

Deciphering the Code for Retroviral Integration Target Site Selection

Upon cell invasion, retroviruses generate a DNA copy of their RNA genome and integrate retroviral cDNA within host chromosomal DNA. Integration occurs throughout the host cell genome, but target site selection is not random. Each subgroup of retrovirus is distinguished from the others by attraction to particular features on chromosomes. Despite extensive efforts to identify host factors that interact with retrovirion components or chromosome features predictive of integration, little is known about how integration sites are selected. We attempted to identify markers predictive of retroviral integration by exploiting Precision-Recall methods for extracting information from highly skewed datasets to derive robust and discriminating measures of association. ChIPSeq datasets for more than 60 factors were compared with 14 retroviral integration datasets. When compared with MLV, PERV or XMRV integration sites, strong association was observed with STAT1, acetylation of H3 and H4 at several positions, and methylation of H2AZ, H3K4, and K9. By combining peaks from ChIPSeq datasets, a supermarker was identified that localized within 2 kB of 75% of MLV proviruses and detected differences in integration preferences among different cell types. The supermarker predicted the likelihood of integration within specific chromosomal regions in a cell-type specific manner, yielding probabilities for integration into proto-oncogene LMO2 identical to experimentally determined values. The supermarker thus identifies chromosomal features highly favored for retroviral integration, provides clues to the mechanism by which retrovirus integration sites are selected, and offers a tool for predicting cell-type specific proto-oncogene activation by retroviruses

Novel Escape Mutants Suggest an Extensive TRIM5α Binding Site Spanning the Entire Outer Surface of the Murine Leukemia Virus Capsid Protein

After entry into target cells, retroviruses encounter the host restriction factors such as Fv1 and TRIM5α. While it is clear that these factors target retrovirus capsid proteins (CA), recognition remains poorly defined in the absence of structural information. To better understand the binding interaction between TRIM5α and CA, we selected a panel of novel N-tropic murine leukaemia virus (N-MLV) escape mutants by a serial passage of replication competent N-MLV in rhesus macaque TRIM5α (rhTRIM5α)-positive cells using a small percentage of unrestricted cells to allow multiple rounds of virus replication. The newly identified mutations, many of which involve changes in charge, are distributed over the outer ‘top’ surface of N-MLV CA, including the N-terminal β-hairpin, and map up to 29 Ao apart. Biological characterisation with a number of restriction factors revealed that only one of the new mutations affects restriction by human TRIM5α, indicating significant differences in the binding interaction between N-MLV and the two TRIM5αs, whereas three of the mutations result in dual sensitivity to Fv1n and Fv1b. Structural studies of two mutants show that no major changes in the overall CA conformation are associated with escape from restriction. We conclude that interactions involving much, if not all, of the surface of CA are vital for TRIM5α binding

Checkpoints in a Yeast Differentiation Pathway Coordinate Signaling during Hyperosmotic Stress

All eukaryotes have the ability to detect and respond to environmental and hormonal signals. In many cases these signals evoke cellular changes that are incompatible and must therefore be orchestrated by the responding cell. In the yeast Saccharomyces cerevisiae, hyperosmotic stress and mating pheromones initiate signaling cascades that each terminate with a MAP kinase, Hog1 and Fus3, respectively. Despite sharing components, these pathways are initiated by distinct inputs and produce distinct cellular behaviors. To understand how these responses are coordinated, we monitored the pheromone response during hyperosmotic conditions. We show that hyperosmotic stress limits pheromone signaling in at least three ways. First, stress delays the expression of pheromone-induced genes. Second, stress promotes the phosphorylation of a protein kinase, Rck2, and thereby inhibits pheromone-induced protein translation. Third, stress promotes the phosphorylation of a shared pathway component, Ste50, and thereby dampens pheromone-induced MAPK activation. Whereas all three mechanisms are dependent on an increase in osmolarity, only the phosphorylation events require Hog1. These findings reveal how an environmental stress signal is able to postpone responsiveness to a competing differentiation signal, by acting on multiple pathway components, in a coordinated manner