A Real-Time All-Atom Structural Search Engine for Proteins

Protein designers use a wide variety of software tools for de novo design, yet their repertoire still lacks a fast and interactive all-atom search engine. To solve this, we have built the Suns program: a real-time, atomic search engine integrated into the PyMOL molecular visualization system. Users build atomic-level structural search queries within PyMOL and receive a stream of search results aligned to their query within a few seconds. This instant feedback cycle enables a new “designability”-inspired approach to protein design where the designer searches for and interactively incorporates native-like fragments from proven protein structures. We demonstrate the use of Suns to interactively build protein motifs, tertiary interactions, and to identify scaffolds compatible with hot-spot residues. The official web site and installer are located at http://www.degradolab.org/suns/ and the source code is hosted at https://github.com/godotgildor/Suns (PyMOL plugin, BSD license), https://github.com/Gabriel439/suns-cmd (command line client, BSD license), and https://github.com/Gabriel439/suns-search (search engine server, GPLv2 license).

This is a PLOS Computational Biology Software Article

     

Approximate Bayesian Computation

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

This is a “Topic Page” article for PLOS Computational Biology.

     

Flow Cytometry Bioinformatics

Flow cytometry bioinformatics is the application of bioinformatics to flow cytometry data, which involves storing, retrieving, organizing, and analyzing flow cytometry data using extensive computational resources and tools. Flow cytometry bioinformatics requires extensive use of and contributes to the development of techniques from computational statistics and machine learning. Flow cytometry and related methods allow the quantification of multiple independent biomarkers on large numbers of single cells. The rapid growth in the multidimensionality and throughput of flow cytometry data, particularly in the 2000s, has led to the creation of a variety of computational analysis methods, data standards, and public databases for the sharing of results. Computational methods exist to assist in the preprocessing of flow cytometry data, identifying cell populations within it, matching those cell populations across samples, and performing diagnosis and discovery using the results of previous steps. For preprocessing, this includes compensating for spectral overlap, transforming data onto scales conducive to visualization and analysis, assessing data for quality, and normalizing data across samples and experiments. For population identification, tools are available to aid traditional manual identification of populations in two-dimensional scatter plots (gating), to use dimensionality reduction to aid gating, and to find populations automatically in higher dimensional space in a variety of ways. It is also possible to characterize data in more comprehensive ways, such as the density-guided binary space partitioning technique known as probability binning, or by combinatorial gating. Finally, diagnosis using flow cytometry data can be aided by supervised learning techniques, and discovery of new cell types of biological importance by high-throughput statistical methods, as part of pipelines incorporating all of the aforementioned methods. Open standards, data, and software are also key parts of flow cytometry bioinformatics. Data standards include the widely adopted Flow Cytometry Standard (FCS) defining how data from cytometers should be stored, but also several new standards under development by the International Society for Advancement of Cytometry (ISAC) to aid in storing more detailed information about experimental design and analytical steps. Open data is slowly growing with the opening of the CytoBank database in 2010 and FlowRepository in 2012, both of which allow users to freely distribute their data, and the latter of which has been recommended as the preferred repository for MIFlowCyt-compliant data by ISAC. Open software is most widely available in the form of a suite of Bioconductor packages, but is also available for web execution on the GenePattern platform.

This is a “Topic Page” article for PLOS Computational Biology.

     

Remodeling of the Rad51 DNA Strand-Exchange Protein by the Srs2 Helicase

Homologous recombination is associated with the dynamic assembly and disassembly of DNA–protein complexes. Assembly of a nucleoprotein filament comprising ssDNA and the RecA homolog, Rad51, is a key step required for homology search during recombination. The budding yeast Srs2 DNA translocase is known to dismantle Rad51 filament in vitro. However, there is limited evidence to support the dismantling activity of Srs2 in vivo. Here, we show that Srs2 indeed disrupts Rad51-containing complexes from chromosomes during meiosis. Overexpression of Srs2 during the meiotic prophase impairs meiotic recombination and removes Rad51 from meiotic chromosomes. This dismantling activity is specific for Rad51, as Srs2 Overexpression does not remove Dmc1 (a meiosis-specific Rad51 homolog), Rad52 (a Rad51 mediator), or replication protein A (RPA; a single-stranded DNA-binding protein). Rather, RPA replaces Rad51 under these conditions. A mutant Srs2 lacking helicase activity cannot remove Rad51 from meiotic chromosomes. Interestingly, the Rad51-binding domain of Srs2, which is critical for Rad51-dismantling activity in vitro, is not essential for this activity in vivo. Our results suggest that a precise level of Srs2, in the form of the Srs2 translocase, is required to appropriately regulate the Rad51 nucleoprotein filament dynamics during meiosis.     

PathVisio 3: An Extendable Pathway Analysis Toolbox

PathVisio is a commonly used pathway editor, visualization and analysis software. Biological pathways have been used by biologists for many years to describe the detailed steps in biological processes. Those powerful, visual representations help researchers to better understand, share and discuss knowledge. Since the first publication of PathVisio in 2008, the original paper was cited more than 170 times and PathVisio was used in many different biological studies. As an online editor PathVisio is also integrated in the community curated pathway database WikiPathways.Here we present the third version of PathVisio with the newest additions and improvements of the application. The core features of PathVisio are pathway drawing, advanced data visualization and pathway statistics. Additionally, PathVisio 3 introduces a new powerful extension systems that allows other developers to contribute additional functionality in form of plugins without changing the core application.PathVisio can be downloaded from http://www.pathvisio.org and in 2014 PathVisio 3 has been downloaded over 5,500 times. There are already more than 15 plugins available in the central plugin repository. PathVisio is a freely available, open-source tool published under the Apache 2.0 license (http://www.apache.org/licenses/LICENSE-2.0). It is implemented in Java and thus runs on all major operating systems. The code repository is available at http://svn.bigcat.unimaas.nl/pathvisio. The support mailing list for users is available on https://groups.google.com/forum/#!forum/wikipathways-discuss and for developers on https://groups.google.com/forum/#!forum/wikipathways-devel.

This is a PLOS Computational Biology software article.

     

High quality draft genome sequence of the slightly halophilic bacterium Halomonas zhanjiangensis type strain JSM 078169T (DSM 21076T) from a sea urchin in southern China

Halomonas zhanjiangensis Chen et al. 2009 is a member of the genus Halomonas, family Halomonadaceae, class Gammaproteobacteria. Representatives of the genus Halomonas are a group of halophilic bacteria often isolated from salty environments. The type strain H. zhanjiangensis JSM 078169^T was isolated from a sea urchin (Hemicentrotus pulcherrimus) collected from the South China Sea. The genome of strain JSM 078169^T is the fourteenth sequenced genome in the genus Halomonas and the fifteenth in the family Halomonadaceae. The other thirteen genomes from the genus Halomonas are H. halocynthiae, H. venusta, H. alkaliphila, H. lutea, H. anticariensis, H. jeotgali, H. titanicae, H. desiderata, H. smyrnensis, H. salifodinae, H. boliviensis, H. elongata and H stevensii. Here, we describe the features of strain JSM 078169^T, together with the complete genome sequence and annotation from a culture of DSM 21076^T. The 4,060,520 bp long draft genome consists of 17 scaffolds with the 3,659 protein-coding and 80 RNA genes and is a part of Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG) project.     

DNA Replication Checkpoint Signaling Depends on a Rad53–Dbf4 N-Terminal Interaction in Saccharomyces cerevisiae

Dbf4-dependent kinase (DDK) and cyclin-dependent kinase (CDK) are essential to initiate DNA replication at individual origins. During replication stress, the S-phase checkpoint inhibits the DDK- and CDK-dependent activation of late replication origins. Rad53 kinase is a central effector of the replication checkpoint and both binds to and phosphorylates Dbf4 to prevent late-origin firing. The molecular basis for the Rad53–Dbf4 physical interaction is not clear but occurs through the Dbf4 N terminus. Here we found that both Rad53 FHA1 and FHA2 domains, which specifically recognize phospho-threonine (pT), interacted with Dbf4 through an N-terminal sequence and an adjacent BRCT domain. Purified Rad53 FHA1 domain (but not FHA2) bound to a pT Dbf4 peptide in vitro, suggesting a possible phospho-threonine-dependent interaction between FHA1 and Dbf4. The Dbf4–Rad53 interaction is governed by multiple contacts that are separable from the Cdc5- and Msa1-binding sites in the Dbf4 N terminus. Importantly, abrogation of the Rad53–Dbf4 physical interaction blocked Dbf4 phosphorylation and allowed late-origin firing during replication checkpoint activation. This indicated that Rad53 must stably bind to Dbf4 to regulate its activity.     

Wide-Ranging Effects of the Yeast Ptc1 Protein Phosphatase Acting Through the MAPK Kinase Mkk1

The Saccharomyces cerevisiae type 2C protein phosphatase Ptc1 is required for a wide variety of cellular functions, although only a few cellular targets have been identified. A genetic screen in search of mutations in protein kinase–encoding genes able to suppress multiple phenotypic traits caused by the ptc1 deletion yielded a single gene, MKK1, coding for a MAPK kinase (MAPKK) known to activate the cell-wall integrity (CWI) Slt2 MAPK. In contrast, mutation of the MKK1 paralog, MKK2, had a less significant effect. Deletion of MKK1 abolished the increased phosphorylation of Slt2 induced by the absence of Ptc1 both under basal and CWI pathway stimulatory conditions. We demonstrate that Ptc1 acts at the level of the MAPKKs of the CWI pathway, but only the Mkk1 kinase activity is essential for ptc1 mutants to display high Slt2 activation. We also show that Ptc1 is able to dephosphorylate Mkk1 in vitro. Our results reveal the preeminent role of Mkk1 in signaling through the CWI pathway and strongly suggest that hyperactivation of Slt2 caused by upregulation of Mkk1 is at the basis of most of the phenotypic defects associated with lack of Ptc1 function.     

Nucleoporin FG Domains Facilitate mRNP Remodeling at the Cytoplasmic Face of the Nuclear Pore Complex

Directional export of messenger RNA (mRNA) protein particles (mRNPs) through nuclear pore complexes (NPCs) requires multiple factors. In Saccharomyces cerevisiae, the NPC proteins Nup159 and Nup42 are asymmetrically localized to the cytoplasmic face and have distinct functional domains: a phenylalanine-glycine (FG) repeat domain that docks mRNP transport receptors and domains that bind the DEAD-box ATPase Dbp5 and its activating cofactor Gle1, respectively. We speculated that the Nup42 and Nup159 FG domains play a role in positioning mRNPs for the terminal mRNP-remodeling steps carried out by Dbp5. Here we find that deletion (Δ) of both the Nup42 and Nup159 FG domains results in a cold-sensitive poly(A)+ mRNA export defect. The nup42ΔFG nup159ΔFG mutant also has synthetic lethal genetic interactions with dbp5 and gle1 mutants. RNA cross-linking experiments further indicate that the nup42ΔFG nup159ΔFG mutant has a reduced capacity for mRNP remodeling during export. To further analyze the role of these FG domains, we replaced the Nup159 or Nup42 FG domains with FG domains from other Nups. These FG “swaps” demonstrate that only certain FG domains are functional at the NPC cytoplasmic face. Strikingly, fusing the Nup42 FG domain to the carboxy-terminus of Gle1 bypasses the need for the endogenous Nup42 FG domain, highlighting the importance of proximal positioning for these factors. We conclude that the Nup42 and Nup159 FG domains target the mRNP to Gle1 and Dbp5 for mRNP remodeling at the NPC. Moreover, these results provide key evidence that character and context play a direct role in FG domain function and mRNA export.     

Non-contiguous genome sequence of Mycobacterium simiae strain DSM 44165T

Mycobacterium simiae is a non-tuberculosis mycobacterium causing pulmonary infections in both immunocompetent and imunocompromized patients. We announce the draft genome sequence of M. simiae DSM 44165^T. The 5,782,968-bp long genome with 65.15% GC content (one chromosome, no plasmid) contains 5,727 open reading frames (33% with unknown function and 11 ORFs sizing more than 5000 -bp), three rRNA operons, 52 tRNA, one 66-bp tmRNA matching with tmRNA tags from Mycobacterium avium, Mycobacterium tuberculosis, Mycobacterium bovis, Mycobacterium microti, Mycobacterium marinum, and Mycobacterium africanum and 389 DNA repetitive sequences. Comparing ORFs and size distribution between M. simiae and five other Mycobacterium species M. simiae clustered with M. abscessus and M. smegmatis. A 40-kb prophage was predicted in addition to two prophage-like elements, 7-kb and 18-kb in size, but no mycobacteriophage was seen after the observation of 10⁶ M. simiae cells. Fifteen putative CRISPRs were found. Three genes were predicted to encode resistance to aminoglycosides, betalactams and macrolide-lincosamide-streptogramin B. A total of 163 CAZYmes were annotated. M. simiae contains ESX-1 to ESX-5 genes encoding for a type-VII secretion system. Availability of the genome sequence may help depict the unique properties of this environmental, opportunistic pathogen.     

The Principal Role of Ku in Telomere Length Maintenance Is Promotion of Est1 Association with Telomeres

Telomere length is tightly regulated in cells that express telomerase. The Saccharomyces cerevisiae Ku heterodimer, a DNA end-binding complex, positively regulates telomere length in a telomerase-dependent manner. Ku associates with the telomerase RNA subunit TLC1, and this association is required for TLC1 nuclear retention. Ku–TLC1 interaction also impacts the cell-cycle-regulated association of the telomerase catalytic subunit Est2 to telomeres. The promotion of TLC1 nuclear localization and Est2 recruitment have been proposed to be the principal role of Ku in telomere length maintenance, but neither model has been directly tested. Here we study the impact of forced recruitment of Est2 to telomeres on telomere length in the absence of Ku’s ability to bind TLC1 or DNA ends. We show that tethering Est2 to telomeres does not promote efficient telomere elongation in the absence of Ku–TLC1 interaction or DNA end binding. Moreover, restoration of TLC1 nuclear localization, even when combined with Est2 recruitment, does not bypass the role of Ku. In contrast, forced recruitment of Est1, which has roles in telomerase recruitment and activation, to telomeres promotes efficient and progressive telomere elongation in the absence of Ku–TLC1 interaction, Ku DNA end binding, or Ku altogether. Ku associates with Est1 and Est2 in a TLC1-dependent manner and enhances Est1 recruitment to telomeres independently of Est2. Together, our results unexpectedly demonstrate that the principal role of Ku in telomere length maintenance is to promote the association of Est1 with telomeres, which may in turn allow for efficient recruitment and activation of the telomerase holoenzyme.     

Integrity and Function of the Saccharomyces cerevisiae Spindle Pole Body Depends on Connections Between the Membrane Proteins Ndc1, Rtn1, and Yop1

The nuclear envelope in Saccharomyces cerevisiae harbors two essential macromolecular protein assemblies: the nuclear pore complexes (NPCs) that enable nucleocytoplasmic transport, and the spindle pole bodies (SPBs) that mediate chromosome segregation. Previously, based on metazoan and budding yeast studies, we reported that reticulons and Yop1/DP1 play a role in the early steps of de novo NPC assembly. Here, we examined if Rtn1 and Yop1 are required for SPB function in S. cerevisiae. Electron microscopy of rtn1Δ yop1Δ cells revealed lobular abnormalities in SPB structure. Using an assay that monitors lateral expansion of the SPB central layer, we found that rtn1Δ yop1Δ SPBs had decreased connections to the NE compared to wild type, suggesting that SPBs are less stable in the NE. Furthermore, large budded rtn1Δ yop1Δ cells exhibited a high incidence of short mitotic spindles, which were frequently misoriented with respect to the mother–daughter axis. This correlated with cytoplasmic microtubule defects. We found that overexpression of the SPB insertion factors NDC1, MPS2, or BBP1 rescued the SPB defects observed in rtn1Δ yop1Δ cells. However, only overexpression of NDC1, which is also required for NPC biogenesis, rescued both the SPB and NPC associated defects. Rtn1 and Yop1 also physically interacted with Ndc1 and other NPC membrane proteins. We propose that NPC and SPB biogenesis are altered in cells lacking Rtn1 and Yop1 due to competition between these complexes for Ndc1, an essential common component of both NPCs and SPBs.     

Translational Control of the Oogenic Program by Components of OMA Ribonucleoprotein Particles in Caenorhabditis elegans

The oocytes of most sexually reproducing animals arrest in meiotic prophase I. Oocyte growth, which occurs during this period of arrest, enables oocytes to acquire the cytoplasmic components needed to produce healthy progeny and to gain competence to complete meiosis. In the nematode Caenorhabditis elegans, the major sperm protein hormone promotes meiotic resumption (also called meiotic maturation) and the cytoplasmic flows that drive oocyte growth. Prior work established that two related TIS11 zinc-finger RNA-binding proteins, OMA-1 and OMA-2, are redundantly required for normal oocyte growth and meiotic maturation. We affinity purified OMA-1 and identified associated mRNAs and proteins using genome-wide expression data and mass spectrometry, respectively. As a class, mRNAs enriched in OMA-1 ribonucleoprotein particles (OMA RNPs) have reproductive functions. Several of these mRNAs were tested and found to be targets of OMA-1/2-mediated translational repression, dependent on sequences in their 3′-untranslated regions (3′-UTRs). Consistent with a major role for OMA-1 and OMA-2 in regulating translation, OMA-1-associated proteins include translational repressors and activators, and some of these proteins bind directly to OMA-1 in yeast two-hybrid assays, including OMA-2. We show that the highly conserved TRIM-NHL protein LIN-41 is an OMA-1-associated protein, which also represses the translation of several OMA-1/2 target mRNAs. In the accompanying article in this issue, we show that LIN-41 prevents meiotic maturation and promotes oocyte growth in opposition to OMA-1/2. Taken together, these data support a model in which the conserved regulators of mRNA translation LIN-41 and OMA-1/2 coordinately control oocyte growth and the proper spatial and temporal execution of the meiotic maturation decision.     

Evolving Digital Ecological Networks

“It is hard to realize that the living world as we know it is just one among many possibilities” [1]. Evolving digital ecological networks are webs of interacting, self-replicating, and evolving computer programs (i.e., digital organisms) that experience the same major ecological interactions as biological organisms (e.g., competition, predation, parasitism, and mutualism). Despite being computational, these programs evolve quickly in an open-ended way, and starting from only one or two ancestral organisms, the formation of ecological networks can be observed in real-time by tracking interactions between the constantly evolving organism phenotypes. These phenotypes may be defined by combinations of logical computations (hereafter tasks) that digital organisms perform and by expressed behaviors that have evolved. The types and outcomes of interactions between phenotypes are determined by task overlap for logic-defined phenotypes and by responses to encounters in the case of behavioral phenotypes. Biologists use these evolving networks to study active and fundamental topics within evolutionary ecology (e.g., the extent to which the architecture of multispecies networks shape coevolutionary outcomes, and the processes involved).

This is a “Topic Page” article for PLOS Computational Biology.