AliSim: A Fast and Versatile Phylogenetic Sequence Simulator for the Genomic Era

Sequence simulators play an important role in phylogenetics. Simulated data has many applications, such as evaluating the performance of different methods, hypothesis testing with parametric bootstraps, and, more recently, generating data for training machine-learning applications. Many sequence simulation programmes exist, but the most feature-rich programmes tend to be rather slow, and the fastest programmes tend to be feature-poor. Here, we introduce AliSim, a new tool that can efficiently simulate biologically realistic alignments under a large range of complex evolutionary models. To achieve high performance across a wide range of simulation conditions, AliSim implements an adaptive approach that combines the commonly used rate matrix and probability matrix approaches. AliSim takes 1.4 h and 1.3 GB RAM to simulate alignments with one million sequences or sites, whereas popular software Seq-Gen, Dawg, and INDELible require 2–5 h and 50–500 GB of RAM. We provide AliSim as an extension of the IQ-TREE software version 2.2, freely available at www.iqtree.org, and a comprehensive user tutorial at http://www.iqtree.org/doc/AliSim.     

The Abundant Histone Chaperones Spt6 and FACT Collaborate to Assemble,Inspect, and Maintain Chromatin Structure in Saccharomyces cerevisiae

Saccharomyces cerevisiae Spt6 protein is a conserved chromatin factor with several distinct functional domains, including a natively unstructured 30-residue N-terminal region that binds competitively with Spn1 or nucleosomes. To uncover physiological roles of these interactions, we isolated histone mutations that suppress defects caused by weakening Spt6:Spn1 binding with the spt6-F249K mutation. The strongest suppressor was H2A-N39K, which perturbs the point of contact between the two H2A-H2B dimers in an assembled nucleosome. Substantial suppression also was observed when the H2A-H2B interface with H3-H4 was altered, and many members of this class of mutations also suppressed a defect in another essential histone chaperone, FACT. Spt6 is best known as an H3-H4 chaperone, but we found that it binds with similar affinity to H2A-H2B or H3-H4. Like FACT, Spt6 is therefore capable of binding each of the individual components of a nucleosome, but unlike FACT, Spt6 did not produce endonuclease-sensitive reorganized nucleosomes and did not displace H2A-H2B dimers from nucleosomes. Spt6 and FACT therefore have distinct activities, but defects can be suppressed by overlapping histone mutations. We also found that Spt6 and FACT together are nearly as abundant as nucleosomes, with ∼24,000 Spt6 molecules, ∼42,000 FACT molecules, and ∼75,000 nucleosomes per cell. Histone mutations that destabilize interfaces within nucleosomes therefore reveal multiple spatial regions that have both common and distinct roles in the functions of these two essential and abundant histone chaperones. We discuss these observations in terms of different potential roles for chaperones in both promoting the assembly of nucleosomes and monitoring their quality.     

Chromatin Remodeling Factors Isw2 and Ino80 Regulate Checkpoint Activity and Chromatin Structure in S Phase

When cells undergo replication stress, proper checkpoint activation and deactivation are critical for genomic stability and cell survival and therefore must be highly regulated. Although mechanisms of checkpoint activation are well studied, mechanisms of checkpoint deactivation are far less understood. Previously, we reported that chromatin remodeling factors Isw2 and Ino80 attenuate the S-phase checkpoint activity in Saccharomyces cerevisiae, especially during recovery from hydroxyurea. In this study, we found that Isw2 and Ino80 have a more pronounced role in attenuating checkpoint activity during late S phase in the presence of methyl methanesulfonate (MMS). We therefore screened for checkpoint factors required for Isw2 and Ino80 checkpoint attenuation in the presence of MMS. Here we demonstrate that Isw2 and Ino80 antagonize checkpoint activators and attenuate checkpoint activity in S phase in MMS either through a currently unknown pathway or through RPA. Unexpectedly, we found that Isw2 and Ino80 increase chromatin accessibility around replicating regions in the presence of MMS through a novel mechanism. Furthermore, through growth assays, we provide additional evidence that Isw2 and Ino80 partially counteract checkpoint activators specifically in the presence of MMS. Based on these results, we propose that Isw2 and Ino80 attenuate S-phase checkpoint activity through a novel mechanism.     

GRIMD: distributed computing for chemists and biologists

Motivation: Biologists and chemists are facing problems of high computational complexity that require the use of several computers organized in clusters or in specialized grids. Examples of such problems can be found in molecular dynamics (MD), in silico screening, and genome analysis. Grid Computing and Cloud Computing are becoming prevalent mainly because of their competitive performance/cost ratio. Regrettably, the diffusion of Grid Computing is strongly limited because two main limitations: it is confined to scientists with strong Computer Science background and the analyses of the large amount of data produced can be cumbersome it. We have developed a package named GRIMD to provide an easy and flexible implementation of distributed computing for the Bioinformatics community. GRIMD is very easy to install and maintain, and it does not require any specific Computer Science skill. Moreover, permits preliminary analysis on the distributed machines to reduce the amount of data to transfer. GRIMD is very flexible because it shields the typical computational biologist from the need to write specific code for tasks such as molecular dynamics or docking calculations. Furthermore, it permits an efficient use of GPU cards whenever is possible. GRIMD calculations scale almost linearly and, therefore, permits to exploit efficiently each machine in the network. Here, we provide few examples of grid computing in computational biology (MD and docking) and bioinformatics (proteome analysis).

Availability

GRIMD is available for free for noncommercial research at www.yadamp.unisa.it/grimd

Supplementary information

www.yadamp.unisa.it/grimd/howto.aspx     

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.     

Systematic Analyses of rpm-1 Suppressors Reveal Roles for ESS-2 in mRNA Splicing in Caenorhabditis elegans

The PHR (Pam/Highwire/RPM-1) family of ubiquitin E3 ligases plays conserved roles in axon patterning and synaptic development. Genetic modifier analysis has greatly aided the discovery of the signal transduction cascades regulated by these proteins. In Caenorhabditis elegans, loss of function in rpm-1 causes axon overgrowth and aberrant presynaptic morphology, yet the mutant animals exhibit little behavioral deficits. Strikingly, rpm-1 mutations strongly synergize with loss of function in the presynaptic active zone assembly factors, syd-1 and syd-2, resulting in severe locomotor deficits. Here, we provide ultrastructural evidence that double mutants, between rpm-1 and syd-1 or syd-2, dramatically impair synapse formation. Taking advantage of the synthetic locomotor defects to select for genetic suppressors, previous studies have identified the DLK-1 MAP kinase cascade negatively regulated by RPM-1. We now report a comprehensive analysis of a large number of suppressor mutations of this screen. Our results highlight the functional specificity of the DLK-1 cascade in synaptogenesis. We also identified two previously uncharacterized genes. One encodes a novel protein, SUPR-1, that acts cell autonomously to antagonize RPM-1. The other affects a conserved protein ESS-2, the homolog of human ES2 or DGCR14. Loss of function in ess-2 suppresses rpm-1 only in the presence of a dlk-1 splice acceptor mutation. We show that ESS-2 acts to promote accurate mRNA splicing when the splice site is compromised. The human DGCR14/ES2 resides in a deleted chromosomal region implicated in DiGeorge syndrome, and its mutation has shown high probability as a risk factor for schizophrenia. Our findings provide the first functional evidence that this family of proteins regulate mRNA splicing in a context-specific manner.     

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.     

The Role of Dbf4-Dependent Protein Kinase in DNA Polymerase ζ-Dependent Mutagenesis in Saccharomyces cerevisiae

The yeast Dbf4-dependent kinase (DDK) (composed of Dbf4 and Cdc7 subunits) is an essential, conserved Ser/Thr protein kinase that regulates multiple processes in the cell, including DNA replication, recombination and induced mutagenesis. Only DDK substrates important for replication and recombination have been identified. Consequently, the mechanism by which DDK regulates mutagenesis is unknown. The yeast mcm5-bob1 mutation that bypasses DDK’s essential role in DNA replication was used here to examine whether loss of DDK affects spontaneous as well as induced mutagenesis. Using the sensitive lys2ΔA746 frameshift reversion assay, we show DDK is required to generate “complex” spontaneous mutations, which are a hallmark of the Polζ translesion synthesis DNA polymerase. DDK co-immunoprecipitated with the Rev7 regulatory, but not with the Rev3 polymerase subunit of Polζ. Conversely, Rev7 bound mainly to the Cdc7 kinase subunit and not to Dbf4. The Rev7 subunit of Polζ may be regulated by DDK phosphorylation as immunoprecipitates of yeast Cdc7 and also recombinant Xenopus DDK phosphorylated GST-Rev7 in vitro. In addition to promoting Polζ-dependent mutagenesis, DDK was also important for generating Polζ-independent large deletions that revert the lys2ΔA746 allele. The decrease in large deletions observed in the absence of DDK likely results from an increase in the rate of replication fork restart after an encounter with spontaneous DNA damage. Finally, nonepistatic, additive/synergistic UV sensitivity was observed in cdc7Δ pol32Δ and cdc7Δ pol30-K127R,K164R double mutants, suggesting that DDK may regulate Rev7 protein during postreplication “gap filling” rather than during “polymerase switching” by ubiquitinated and sumoylated modified Pol30 (PCNA) and Pol32.     

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.     

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.     

The Mps1 Kinase Modulates the Recruitment and Activity of Cnn1CENP-T at Saccharomyces cerevisiae Kinetochores

Kinetochores are conserved protein complexes that bind the replicated chromosomes to the mitotic spindle and then direct their segregation. To better comprehend Saccharomyces cerevisiae kinetochore function, we dissected the phospho-regulated dynamic interaction between conserved kinetochore protein Cnn1^CENP-T, the centromere region, and the Ndc80 complex through the cell cycle. Cnn1 localizes to kinetochores at basal levels from G1 through metaphase but accumulates abruptly at anaphase onset. How Cnn1 is recruited and which activities regulate its dynamic localization are unclear. We show that Cnn1 harbors two kinetochore-localization activities: a C-terminal histone-fold domain (HFD) that associates with the centromere region and a N-terminal Spc24/Spc25 interaction sequence that mediates linkage to the microtubule-binding Ndc80 complex. We demonstrate that the established Ndc80 binding site in the N terminus of Cnn1, Cnn1^60–84, should be extended with flanking residues, Cnn1^25–91, to allow near maximal binding affinity to Ndc80. Cnn1 localization was proposed to depend on Mps1 kinase activity at Cnn1–S74, based on in vitro experiments demonstrating the Cnn1–Ndc80 complex interaction. We demonstrate that from G1 through metaphase, Cnn1 localizes via both its HFD and N-terminal Spc24/Spc25 interaction sequence, and deletion or mutation of either region results in anomalous Cnn1 kinetochore levels. At anaphase onset (when Mps1 activity decreases) Cnn1 becomes enriched mainly via the N-terminal Spc24/Spc25 interaction sequence. In sum, we provide the first in vivo evidence of Cnn1 preanaphase linkages with the kinetochore and enrichment of the linkages during anaphase.     

New Regulators of Clathrin-Mediated Endocytosis Identified in Saccharomyces cerevisiae by Systematic Quantitative Fluorescence Microscopy

Despite the importance of clathrin-mediated endocytosis (CME) for cell biology, it is unclear if all components of the machinery have been discovered and many regulatory aspects remain poorly understood. Here, using Saccharomyces cerevisiae and a fluorescence microscopy screening approach we identify previously unknown regulatory factors of the endocytic machinery. We further studied the top scoring protein identified in the screen, Ubx3, a member of the conserved ubiquitin regulatory X (UBX) protein family. In vivo and in vitro approaches demonstrate that Ubx3 is a new coat component. Ubx3-GFP has typical endocytic coat protein dynamics with a patch lifetime of 45 ± 3 sec. Ubx3 contains a W-box that mediates physical interaction with clathrin and Ubx3-GFP patch lifetime depends on clathrin. Deletion of the UBX3 gene caused defects in the uptake of Lucifer Yellow and the methionine transporter Mup1 demonstrating that Ubx3 is needed for efficient endocytosis. Further, the UBX domain is required both for localization and function of Ubx3 at endocytic sites. Mechanistically, Ubx3 regulates dynamics and patch lifetime of the early arriving protein Ede1 but not later arriving coat proteins or actin assembly. Conversely, Ede1 regulates the patch lifetime of Ubx3. Ubx3 likely regulates CME via the AAA-ATPase Cdc48, a ubiquitin-editing complex. Our results uncovered new components of the CME machinery that regulate this fundamental process.     

The Jaw of the Worm: GTPase-activating Protein EAT-17 Regulates Grinder Formation in Caenorhabditis elegans

Constitutive transport of cellular materials is essential for cell survival. Although multiple small GTPase Rab proteins are required for the process, few regulators of Rabs are known. Here we report that EAT-17, a novel GTPase-activating protein (GAP), regulates RAB-6.2 function in grinder formation in Caenorhabditis elegans. We identified EAT-17 as a novel RabGAP that interacts with RAB-6.2, a protein that presumably regulates vesicle trafficking between Golgi, the endoplasmic reticulum, and plasma membrane to form a functional grinder. EAT-17 has a canonical GAP domain that is critical for its function. RNA interference against 25 confirmed and/or predicted RABs in C. elegans shows that RNAi against rab-6.2 produces a phenotype identical to eat-17. A directed yeast two-hybrid screen using EAT-17 as bait and each of the 25 RAB proteins as prey identifies RAB-6.2 as the interacting partner of EAT-17, confirming that RAB-6.2 is a specific substrate of EAT-17. Additionally, deletion mutants of rab-6.2 show grinder defects identical to those of eat-17 loss-of-function mutants, and both RAB-6.2 and EAT-17 are expressed in the terminal bulb of the pharynx where the grinder is located. Collectively, these results suggest that EAT-17 is a specific GTPase-activating protein for RAB-6.2. Based on the conserved function of Rab6 in vesicular transport, we propose that EAT-17 regulates the turnover rate of RAB-6.2 activity in cargo trafficking for grinder formation.