Predicting DNA-Binding Proteins and Binding Residues by Complex Structure Prediction and Application to Human Proteome

As more and more protein sequences are uncovered from increasingly inexpensive sequencing techniques, an urgent task is to find their functions. This work presents a highly reliable computational technique for predicting DNA-binding function at the level of protein-DNA complex structures, rather than low-resolution two-state prediction of DNA-binding as most existing techniques do. The method first predicts protein-DNA complex structure by utilizing the template-based structure prediction technique HHblits, followed by binding affinity prediction based on a knowledge-based energy function (Distance-scaled finite ideal-gas reference state for protein-DNA interactions). A leave-one-out cross validation of the method based on 179 DNA-binding and 3797 non-binding protein domains achieves a Matthews correlation coefficient (MCC) of 0.77 with high precision (94%) and high sensitivity (65%). We further found 51% sensitivity for 82 newly determined structures of DNA-binding proteins and 56% sensitivity for the human proteome. In addition, the method provides a reasonably accurate prediction of DNA-binding residues in proteins based on predicted DNA-binding complex structures. Its application to human proteome leads to more than 300 novel DNA-binding proteins; some of these predicted structures were validated by known structures of homologous proteins in APO forms. The method [SPOT-Seq (DNA)] is available as an on-line server at http://sparks-lab.org.     

Probing the Boundaries of Orthology: The Unanticipated Rapid Evolution of Drosophila centrosomin

The rapid evolution of essential developmental genes and their protein products is both intriguing and problematic. The rapid evolution of gene products with simple protein folds and a lack of well-characterized functional domains typically result in a low discovery rate of orthologous genes. Additionally, in the absence of orthologs it is difficult to study the processes and mechanisms underlying rapid evolution. In this study, we have investigated the rapid evolution of centrosomin (cnn), an essential gene encoding centrosomal protein isoforms required during syncytial development in Drosophila melanogaster. Until recently the rapid divergence of cnn made identification of orthologs difficult and questionable because Cnn violates many of the assumptions underlying models for protein evolution. To overcome these limitations, we have identified a group of insect orthologs and present conserved features likely to be required for the functions attributed to cnn in D. melanogaster. We also show that the rapid divergence of Cnn isoforms is apparently due to frequent coding sequence indels and an accelerated rate of intronic additions and eliminations. These changes appear to be buffered by multi-exon and multi-reading frame maximum potential ORFs, simple protein folds, and the splicing machinery. These buffering features also occur in other genes in Drosophila and may help prevent potentially deleterious mutations due to indels in genes with large coding exons and exon-dense regions separated by small introns. This work promises to be useful for future investigations of cnn and potentially other rapidly evolving genes and proteins.     

A New Player in the Spermiogenesis Pathway of Caenorhabditis elegans

Precise timing of sperm activation ensures the greatest likelihood of fertilization. Precision in Caenorhabditis elegans sperm activation is ensured by external signaling, which induces the spherical spermatid to reorganize and extend a pseudopod for motility. Spermatid activation, also called spermiogenesis, is prevented from occurring prematurely by the activity of SPE-6 and perhaps other proteins, termed “the brake model.” Here, we identify the spe-47 gene from the hc198 mutation that causes premature spermiogenesis. The mutation was isolated in a suppressor screen of spe-27(it132ts), which normally renders worms sterile, due to defective transduction of the activation signal. In a spe-27(+) background, spe-47(hc198) causes a temperature-sensitive reduction of fertility, and in addition to premature spermiogenesis, many mutant sperm fail to activate altogether. The hc198 mutation is semidominant, inducing a more severe loss of fertility than do null alleles generated by CRISPR-associated protein 9 (Cas9) technology. The hc198 mutation affects an major sperm protein (MSP) domain, altering a conserved amino acid residue in a β-strand that mediates MSP–MSP dimerization. Both N- and C-terminal SPE-47 reporters associate with the forming fibrous body (FB)-membranous organelle, a specialized sperm organelle that packages MSP and other components during spermatogenesis. Once the FB is fully formed, the SPE-47 reporters dissociate and disappear. SPE-47 reporter localization is not altered by either the hc198 mutation or a C-terminal truncation deleting the MSP domain. The disappearance of SPE-47 reporters prior to the formation of spermatids requires a reevaluation of the brake model for prevention of premature spermatid activation.     

DNA Replication Error-Induced Extinction of Diploid Yeast

Genetic defects in DNA polymerase accuracy, proofreading, or mismatch repair (MMR) induce mutator phenotypes that accelerate adaptation of microbes and tumor cells. Certain combinations of mutator alleles synergistically increase mutation rates to levels that drive extinction of haploid cells. The maximum tolerated mutation rate of diploid cells is unknown. Here, we define the threshold for replication error-induced extinction (EEX) of diploid Saccharomyces cerevisiae. Double-mutant pol3 alleles that carry mutations for defective DNA polymerase-δ proofreading (pol3-01) and accuracy (pol3-L612M or pol3-L612G) induce strong mutator phenotypes in heterozygous diploids (POL3/pol3-01,L612M or POL3/pol3-01,L612G). Both pol3-01,L612M and pol3-01,L612G alleles are lethal in the homozygous state; cells with pol3-01,L612M divide up to 10 times before arresting at random stages in the cell cycle. Antimutator eex mutations in the pol3 alleles suppress this lethality (pol3-01,L612M,eex or pol3-01,L612G,eex). MMR defects synergize with pol3-01,L612M,eex and pol3-01,L612G,eex alleles, increasing mutation rates and impairing growth. Conversely, inactivation of the Dun1 S-phase checkpoint kinase suppresses strong pol3-01,L612M,eex and pol3-01,L612G,eex mutator phenotypes as well as the lethal pol3-01,L612M phenotype. Our results reveal that the lethal error threshold in diploids is 10 times higher than in haploids and likely determined by homozygous inactivation of essential genes. Pronounced loss of fitness occurs at mutation rates well below the lethal threshold, suggesting that mutator-driven cancers may be susceptible to drugs that exacerbate replication errors.     

The Drosophila wings apart Gene Anchors a Novel,Evolutionarily Conserved Pathway of Neuromuscular Development

wings apart (wap) is a recessive, semilethal gene located on the X chromosome in Drosophila melanogaster, which is required for normal wing-vein patterning. We show that the wap mutation also results in loss of the adult jump muscle. We use complementation mapping and gene-specific RNA interference to localize the wap locus to the proximal X chromosome. We identify the annotated gene CG14614 as the gene affected by the wap mutation, since one wap allele contains a non-sense mutation in CG14614, and a genomic fragment containing only CG14614 rescues the jump-muscle phenotypes of two wap mutant alleles. The wap gene lies centromere-proximal to touch-insensitive larva B and centromere-distal to CG14619, which is tentatively assigned as the gene affected in introverted mutants. In mutant wap animals, founder cell precursors for the jump muscle are specified early in development, but are later lost. Through tissue-specific knockdowns, we demonstrate that wap function is required in both the musculature and the nervous system for normal jump-muscle formation. wap/CG14614 is homologous to vertebrate wdr68, DDB1 and CUL4 associated factor 7, which also are expressed in neuromuscular tissues. Thus, our findings provide insight into mechanisms of neuromuscular development in higher animals and facilitate the understanding of neuromuscular diseases that may result from mis-expression of muscle-specific or neuron-specific genes.     

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.     

The Aggregation-Prone Intracellular Serpin SRP-2 Fails to Transit the ER in Caenorhabditis elegans

Familial encephalopathy with neuroserpin inclusions bodies (FENIB) is a serpinopathy that induces a rare form of presenile dementia. Neuroserpin contains a classical signal peptide and like all extracellular serine proteinase inhibitors (serpins) is secreted via the endoplasmic reticulum (ER)–Golgi pathway. The disease phenotype is due to gain-of-function missense mutations that cause neuroserpin to misfold and aggregate within the ER. In a previous study, nematodes expressing a homologous mutation in the endogenous Caenorhabditis elegans serpin, srp-2, were reported to model the ER proteotoxicity induced by an allele of mutant neuroserpin. Our results suggest that SRP-2 lacks a classical N-terminal signal peptide and is a member of the intracellular serpin family. Using confocal imaging and an ER colocalization marker, we confirmed that GFP-tagged wild-type SRP-2 localized to the cytosol and not the ER. Similarly, the aggregation-prone SRP-2 mutant formed intracellular inclusions that localized to the cytosol. Interestingly, wild-type SRP-2, targeted to the ER by fusion to a cleavable N-terminal signal peptide, failed to be secreted and accumulated within the ER lumen. This ER retention phenotype is typical of other obligate intracellular serpins forced to translocate across the ER membrane. Neuroserpin is a secreted protein that inhibits trypsin-like proteinase. SRP-2 is a cytosolic serpin that inhibits lysosomal cysteine peptidases. We concluded that SRP-2 is neither an ortholog nor a functional homolog of neuroserpin. Furthermore, animals expressing an aggregation-prone mutation in SRP-2 do not model the ER proteotoxicity associated with FENIB.     

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.     

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.     

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.     

A Novel Gene Controlling the Timing of Courtship Initiation in Drosophila melanogaster

Over the past 35 years, developmental geneticists have made impressive progress toward an understanding of how genes specify morphology and function, particularly as they relate to the specification of each physical component of an organism. In the last 20 years, male courtship behavior in Drosophila melanogaster has emerged as a robust model system for the study of genetic specification of behavior. Courtship behavior is both complex and innate, and a single gene, fruitless (fru), is both necessary and sufficient for all aspects of the courtship ritual. Typically, loss of male-specific Fruitless protein function results in male flies that perform the courtship ritual incorrectly, slowly, or not at all. Here we describe a novel requirement for fru: we have identified a group of cells in which male Fru proteins are required to reduce the speed of courtship initiation. In addition, we have identified a gene, Trapped in endoderm 1 (Tre1), which is required in these cells for normal courtship and mating behavior. Tre1 encodes a G-protein-coupled receptor required for establishment of cell polarity and cell migration and has previously not been shown to be involved in courtship behavior. We describe the results of feminization of the Tre1-expressing neurons, as well as the effects on courtship behavior of mutation of Tre1. In addition, we show that Tre1 is expressed in a sexually dimorphic pattern in the central and peripheral nervous systems and investigate the role of the Tre1 cells in mate identification.     

Analysis of PHA-1 Reveals a Limited Role in Pharyngeal Development and Novel Functions in Other Tissues

PHA-1 encodes a cytoplasmic protein that is required for embryonic morphogenesis and attachment of the foregut (pharynx) to the mouth (buccal capsule). Previous reports have in some cases suggested that PHA-1 is essential for the differentiation of most or all pharyngeal cell types. By performing mosaic analysis with a recently acquired pha-1 null mutation (tm3671), we found that PHA-1 is not required within most or all pharyngeal cells for their proper specification, differentiation, or function. Rather, our evidence suggests that PHA-1 acts in the arcade or anterior epithelial cells of the pharynx to promote attachment of the pharynx to the future buccal capsule. In addition, PHA-1 appears to be required in the epidermis for embryonic morphogenesis, in the excretory system for osmoregulation, and in the somatic gonad for normal ovulation and fertility. PHA-1 activity is also required within at least a subset of intestinal cells for viability. To better understand the role of PHA-1 in the epidermis, we analyzed several apical junction markers in pha-1(tm3671) homozygous embryos. PHA-1 regulates the expression of several components of two apical junction complexes including AJM-1–DLG-1/discs large complex and the classical cadherin–catenin complex, which may account for the role of PHA-1 in embryonic morphogenesis.     

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.     

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.

     

Emergence of DNA Polymerase ε Antimutators That Escape Error-Induced Extinction in Yeast

DNA polymerases (Pols) ε and δ perform the bulk of yeast leading- and lagging-strand DNA synthesis. Both Pols possess intrinsic proofreading exonucleases that edit errors during polymerization. Rare errors that elude proofreading are extended into duplex DNA and excised by the mismatch repair (MMR) system. Strains that lack Pol proofreading or MMR exhibit a 10- to 100-fold increase in spontaneous mutation rate (mutator phenotype), and inactivation of both Pol δ proofreading (pol3-01) and MMR is lethal due to replication error-induced extinction (EEX). It is unclear whether a similar synthetic lethal relationship exists between defects in Pol ε proofreading (pol2-4) and MMR. Using a plasmid-shuffling strategy in haploid Saccharomyces cerevisiae, we observed synthetic lethality of pol2-4 with alleles that completely abrogate MMR (msh2Δ, mlh1Δ, msh3Δ msh6Δ, or pms1Δ mlh3Δ) but not with partial MMR loss (msh3Δ, msh6Δ, pms1Δ, or mlh3Δ), indicating that high levels of unrepaired Pol ε errors drive extinction. However, variants that escape this error-induced extinction (eex mutants) frequently emerged. Five percent of pol2-4 msh2Δ eex mutants encoded second-site changes in Pol ε that reduced the pol2-4 mutator phenotype between 3- and 23-fold. The remaining eex alleles were extragenic to pol2-4. The locations of antimutator amino-acid changes in Pol ε and their effects on mutation spectra suggest multiple mechanisms of mutator suppression. Our data indicate that unrepaired leading- and lagging-strand polymerase errors drive extinction within a few cell divisions and suggest that there are polymerase-specific pathways of mutator suppression. The prevalence of suppressors extragenic to the Pol ε gene suggests that factors in addition to proofreading and MMR influence leading-strand DNA replication fidelity.     

Genetic Interactions Between UNC-17/VAChT and a Novel Transmembrane Protein in Caenorhabditis elegans

The unc-17 gene encodes the vesicular acetylcholine transporter (VAChT) in Caenorhabditis elegans. unc-17 reduction-of-function mutants are small, slow growing, and uncoordinated. Several independent unc-17 alleles are associated with a glycine-to-arginine substitution (G347R), which introduces a positive charge in the ninth transmembrane domain (TMD) of UNC-17. To identify proteins that interact with UNC-17/VAChT, we screened for mutations that suppress the uncoordinated phenotype of UNC-17(G347R) mutants. We identified several dominant allele-specific suppressors, including mutations in the sup-1 locus. The sup-1 gene encodes a single-pass transmembrane protein that is expressed in a subset of neurons and in body muscles. Two independent suppressor alleles of sup-1 are associated with a glycine-to-glutamic acid substitution (G84E), resulting in a negative charge in the SUP-1 TMD. A sup-1 null mutant has no obvious deficits in cholinergic neurotransmission and does not suppress unc-17 mutant phenotypes. Bimolecular fluorescence complementation (BiFC) analysis demonstrated close association of SUP-1 and UNC-17 in synapse-rich regions of the cholinergic nervous system, including the nerve ring and dorsal nerve cords. These observations suggest that UNC-17 and SUP-1 are in close proximity at synapses. We propose that electrostatic interactions between the UNC-17(G347R) and SUP-1(G84E) TMDs alter the conformation of the mutant UNC-17 protein, thereby restoring UNC-17 function; this is similar to the interaction between UNC-17/VAChT and synaptobrevin.