An automated method for finding molecular complexes in large protein interaction networks

BACKGROUND: Recent advances in proteomics technologies such as two-hybrid, phage display and mass spectrometry have enabled us to create a detailed map of biomolecular interaction networks. Initial mapping efforts have already produced a wealth of data. As the size of the interaction set increases, databases and computational methods will be required to store, visualize and analyze the information in order to effectively aid in knowledge discovery. RESULTS: This paper describes a novel graph theoretic clustering algorithm, "Molecular Complex Detection" (MCODE), that detects densely connected regions in large protein-protein interaction networks that may represent molecular complexes. The method is based on vertex weighting by local neighborhood density and outward traversal from a locally dense seed protein to isolate the dense regions according to given parameters. The algorithm has the advantage over other graph clustering methods of having a directed mode that allows fine-tuning of clusters of interest without considering the rest of the network and allows examination of cluster interconnectivity, which is relevant for protein networks. Protein interaction and complex information from the yeast Saccharomyces cerevisiae was used for evaluation. CONCLUSION: Dense regions of protein interaction networks can be found, based solely on connectivity data, many of which correspond to known protein complexes. The algorithm is not affected by a known high rate of false positives in data from high-throughput interaction techniques. The program is available from

SeqHound: biological sequence and structure database as a platform for bioinformatics research

BACKGROUND: SeqHound has been developed as an integrated biological sequence, taxonomy, annotation and 3-D structure database system. It provides a high-performance server platform for bioinformatics research in a locally-hosted environment. RESULTS: SeqHound is based on the National Center for Biotechnology Information data model and programming tools. It offers daily updated contents of all Entrez sequence databases in addition to 3-D structural data and information about sequence redundancies, sequence neighbours, taxonomy, complete genomes, functional annotation including Gene Ontology terms and literature links to PubMed. SeqHound is accessible via a web server through a Perl, C or C++ remote API or an optimized local API. It provides functionality necessary to retrieve specialized subsets of sequences, structures and structural domains. Sequences may be retrieved in FASTA, GenBank, ASN.1 and XML formats. Structures are available in ASN.1, XML and PDB formats. Emphasis has been placed on complete genomes, taxonomy, domain and functional annotation as well as 3-D structural functionality in the API, while fielded text indexing functionality remains under development. SeqHound also offers a streamlined WWW interface for simple web-user queries. CONCLUSIONS: The system has proven useful in several published bioinformatics projects such as the BIND database and offers a cost-effective infrastructure for research. SeqHound will continue to develop and be provided as a service of the Blueprint Initiative at the Samuel Lunenfeld Research Institute. The source code and examples are available under the terms of the GNU public license at the Sourceforge site http://sourceforge.net/projects/slritools/ in the SLRI Toolkit

Structural analysis of the starfish SALMFamide neuropeptides S1 and S2: The N-terminal region of S2 facilitates self-association

The neuropeptides S1 (GFNSALMFamide) and S2 (SGPYSFNSGLTFamide), which share sequence similarity, were discovered in the starfish Asterias rubens and are prototypical members of the SALMFamide family of neuropeptides in echinoderms. SALMFamide neuropeptides act as muscle relaxants and both S1 and S2 cause relaxation of cardiac stomach and tube foot preparations in vitro but S2 is an order of magnitude more potent than S1. Here we investigated a structural basis for this difference in potency using spectroscopic techniques. Circular dichroism spectroscopy showed that S1 does not have a defined structure in aqueous solution and this was supported by 2D nuclear magnetic resonance experiments. In contrast, we found that S2 has a well-defined conformation in aqueous solution. However, the conformation of S2 was concentration dependent, with increasing concentration inducing a transition from an unstructured to a structured conformation. Interestingly, this property of S2 was not observed in an N-terminally truncated analogue of S2 (short S2 or SS2; SFNSGLTFamide). Collectively, the data obtained indicate that the N-terminal region of S2 facilitates peptide self-association at high concentrations, which may have relevance to the biosynthesis and/or bioactivity of S2 in vivo

miR-34a Silences c-SRC to Attenuate Tumor Growth in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is an aggressive subtype with no clinically proven biologically targeted treatment options. The molecular heterogeneity of TNBC and lack of high frequency driver mutations other than TP53 have hindered the development of new and effective therapies that significantly improve patient outcomes. miRNAs, global regulators of survival and proliferation pathways important in tumor development and maintenance, are becoming promising therapeutic agents. We performed miRNA-profiling studies in different TNBC subtypes to identify miRNAs that significantly contribute to disease progression. We found that miR-34a was lost in TNBC, specifically within mesenchymal and mesenchymal stem cell-like subtypes, whereas expression of miR-34a targets was significantly enriched. Furthermore, restoration of miR-34a in cell lines representing these subtypes inhibited proliferation and invasion, activated senescence, and promoted sensitivity to dasatinib by targeting the proto-oncogene c-SRC. Notably, SRC depletion in TNBC cell lines phenocopied the effects of miR-34a reintroduction, whereas SRC overexpression rescued the antitumorigenic properties mediated by miR-34a. miR-34a levels also increased when cells were treated with c-SRC inhibitors, suggesting a negative feedback exists between miR-34a and c-SRC. Moreover, miR-34a administration significantly delayed tumor growth of subcutaneously and orthotopically implanted tumors in nude mice, and was accompanied by c-SRC downregulation. Finally, we found that miR-34a and SRC levels were inversely correlated in human tumor specimens. Together, our results demonstrate that miR-34a exerts potent antitumorigenic effects in vitro and in vivo and suggests that miR-34a replacement therapy, which is currently being tested in human clinical trials, represents a promising therapeutic strategy for TNBC. Cancer Res; 76(4); 1-13. (c)2015 AACR

Alpha-Amino-Beta-Carboxy-Muconate-Semialdehyde Decarboxylase Controls Dietary Niacin Requirements for NAD+ Synthesis

NAD+ is essential for redox reactions in energy metabolism and necessary for DNA repair and epigenetic modification. Humans require sufficient amounts of dietary niacin (nicotinic acid, nicotinamide, and nicotinamide riboside) for adequate NAD+ synthesis. In contrast, mice easily generate sufficient NAD+ solely from tryptophan through the kynurenine pathway. We show that transgenic mice with inducible expression of human alpha-amino-beta-carboxy-muconate-semialdehyde decarboxylase (ACMSD) become niacin dependent similar to humans when ACMSD expression is high. On niacin-free diets, these acquired niacin dependency (ANDY) mice developed reversible, mild-to-severe NAD+ deficiency, depending on the nutrient composition of the diet. NAD deficiency in mice contributed to behavioral and health changes that are reminiscent of human niacin deficiency. This study shows that ACMSD is a key regulator of mammalian dietary niacin requirements and NAD+ metabolism and that the ANDY mouse represents a versatile platform for investigating pathologies linked to low NAD+ levels in aging and neurodegenerative diseases

A human MAP kinase interactome.

Mitogen-activated protein kinase (MAPK) pathways form the backbone of signal transduction in the mammalian cell. Here we applied a systematic experimental and computational approach to map 2,269 interactions between human MAPK-related proteins and other cellular machinery and to assemble these data into functional modules. Multiple lines of evidence including conservation with yeast supported a core network of 641 interactions. Using small interfering RNA knockdowns, we observed that approximately one-third of MAPK-interacting proteins modulated MAPK-mediated signaling. We uncovered the Na-H exchanger NHE1 as a potential MAPK scaffold, found links between HSP90 chaperones and MAPK pathways and identified MUC12 as the human analog to the yeast signaling mucin Msb2. This study makes available a large resource of MAPK interactions and clone libraries, and it illustrates a methodology for probing signaling networks based on functional refinement of experimentally derived protein-interaction maps

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization.

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal mendelian long-QT syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals, we identified 35 common variant loci associated with QT interval that collectively explain ∼8-10% of QT-interval variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 new QT interval-associated loci in 298 unrelated probands with LQTS identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies new candidate genes for ventricular arrhythmias, LQTS and SCD

Deep functional analysis of synII, a 770-kilobase synthetic yeast chromosome

INTRODUCTION Although much effort has been devoted to studying yeast in the past few decades, our understanding of this model organism is still limited. Rapidly developing DNA synthesis techniques have made a “build-to-understand” approach feasible to reengineer on the genome scale. Here, we report on the completion of a 770-kilobase synthetic yeast chromosome II (synII). SynII was characterized using extensive Trans-Omics tests. Despite considerable sequence alterations, synII is virtually indistinguishable from wild type. However, an up-regulation of translational machinery was observed and can be reversed by restoring the transfer RNA (tRNA) gene copy number. RATIONALE Following the “design-build-test-debug” working loop, synII was successfully designed and constructed in vivo. Extensive Trans-Omics tests were conducted, including phenomics, transcriptomics, proteomics, metabolomics, chromosome segregation, and replication analyses. By both complementation assays and SCRaMbLE (synthetic chromosome rearrangement and modification by loxP -mediated evolution), we targeted and debugged the origin of a growth defect at 37°C in glycerol medium. RESULTS To efficiently construct megabase-long chromosomes, we developed an I- Sce I–mediated strategy, which enables parallel integration of synthetic chromosome arms and reduced the overall integration time by 50% for synII. An I- Sce I site is introduced for generating a double-strand break to promote targeted homologous recombination during mitotic growth. Despite hundreds of modifications introduced, there are still regions sharing substantial sequence similarity that might lead to undesirable meiotic recombinations when intercrossing the two semisynthetic chromosome arm strains. Induction of the I- Sce I–mediated double-strand break is otherwise lethal and thus introduced a strong selective pressure for targeted homologous recombination. Since our strategy is designed to generate a markerless synII and leave the URA3 marker on the wild-type chromosome, we observed a tenfold increase in URA3 -deficient colonies upon I- Sce I induction, meaning that our strategy can greatly bias the crossover events toward the designated regions. By incorporating comprehensive phenotyping approaches at multiple levels, we demonstrated that synII was capable of powering the growth of yeast indistinguishably from wild-type cells (see the figure), showing highly consistent biological processes comparable to the native strain. Meanwhile, we also noticed modest but potentially significant up-regulation of the translational machinery. The main alteration underlying this change in expression is the deletion of 13 tRNA genes. A growth defect was observed in one very specific condition—high temperature (37°C) in medium with glycerol as a carbon source—where colony size was reduced significantly. We targeted and debugged this defect by two distinct approaches. The first approach involved phenotype screening of all intermediate strains followed by a complementation assay with wild-type sequences in the synthetic strain. By doing so, we identified a modification resulting from PCRTag recoding in TSC10 , which is involved in regulation of the yeast high-osmolarity glycerol (HOG) response pathway. After replacement with wild-type TSC10 , the defect was greatly mitigated. The other approach, debugging by SCRaMbLE, showed rearrangements in regions containing HOG regulation genes. Both approaches indicated that the defect is related to HOG response dysregulation. Thus, the phenotypic defect can be pinpointed and debugged through multiple alternative routes in the complex cellular interactome network. CONCLUSION We have demonstrated that synII segregates, replicates, and functions in a highly similar fashion compared with its wild-type counterpart. Furthermore, we believe that the iterative “design-build-test-debug” cycle methodology, established here, will facilitate progression of the Sc2.0 project in the face of the increasing synthetic genome complexity. SynII characterization. ( A ) Cell cycle comparison between synII and BY4741 revealed by the percentage of cells with separated CEN2-GFP dots, metaphase spindles, and anaphase spindles. ( B ) Replication profiling of synII (red) and BY4741 (black) expressed as relative copy number by deep sequencing. ( C ) RNA sequencing analysis revealed that the significant up-regulation of translational machinery in synII is induced by the deletion of tRNA genes in synII. </jats:sec

Assessing the functional coherence of modules found in multiple-evidence networks from Arabidopsis

<p>Abstract</p> <p>Background</p> <p>Combining multiple evidence-types from different information sources has the potential to reveal new relationships in biological systems. The integrated information can be represented as a relationship network, and clustering the network can suggest possible functional modules. The value of such modules for gaining insight into the underlying biological processes depends on their functional coherence. The challenges that we wish to address are to define and quantify the functional coherence of modules in relationship networks, so that they can be used to infer function of as yet unannotated proteins, to discover previously unknown roles of proteins in diseases as well as for better understanding of the regulation and interrelationship between different elements of complex biological systems.</p> <p>Results</p> <p>We have defined the functional coherence of modules with respect to the Gene Ontology (GO) by considering two complementary aspects: (i) the fragmentation of the GO functional categories into the different modules and (ii) the most representative functions of the modules. We have proposed a set of metrics to evaluate these two aspects and demonstrated their utility in <it>Arabidopsis thaliana</it>. We selected 2355 proteins for which experimentally established protein-protein interaction (PPI) data were available. From these we have constructed five relationship networks, four based on single types of data: PPI, co-expression, co-occurrence of protein names in scientific literature abstracts and sequence similarity and a fifth one combining these four evidence types. The ability of these networks to suggest biologically meaningful grouping of proteins was explored by applying Markov clustering and then by measuring the functional coherence of the clusters.</p> <p>Conclusions</p> <p>Relationship networks integrating multiple evidence-types are biologically informative and allow more proteins to be assigned to a putative functional module. Using additional evidence types concentrates the functional annotations in a smaller number of modules without unduly compromising their consistency. These results indicate that integration of more data sources improves the ability to uncover functional association between proteins, both by allowing more proteins to be linked and producing a network where modular structure more closely reflects the hierarchy in the gene ontology.</p

Lifetime measurement of the 2₁⁺ state in ⁷⁴Rb and isospin properties of quadrupole transition strengths at N = Z

Self-conjugate nuclei in the A ≈ 70–80 region have attracted a great deal of attention due to phenomena such as shape coexistence and increasing collectivity along the N=Zline. We investigate the structure of nuclei in this region through lifetime measurements using the GRETINA array. The first implementation of the Differential Recoil Distance Doppler Shift technique with fast radioactive beams is demonstrated and verified through a measurement of the well-known B(E2; 2₁⁺→ 0₁⁺)transition strength in ⁷⁴Kr. The method is then applied to determine the B(E2; 2₁⁺→ 0₁⁺) transition strength in ⁷⁴Rb, the heaviest odd–odd N=Z nucleus for which this quantity has been determined. This result and extended systematics along N=Z suggest the dominance of the isoscalar part of the quadrupole transition strengths in self-conjugate nuclei, as well as the possible presence of shape coexistence in ⁷⁴Rb