Identifying the hotspots on the top faces of WD40-repeat proteins from their primary sequences by β-bulges and DHSW tetrads

Characterization of the RACK1 gene of Aips cerana cerana and its role in adverse environmental stresses

Receptors for Activated C Kinase 1 (RACK1s) are a kind of multifunction scaffold protein that plays an important role in cell signal transductions and animal development. However, the function of RACK1 in the Chinese honeybee Apis cerana cerana is little known. Here, we isolated and identified a RACK1 gene from Apis cerana cerana, named AccRACK1. By bioinformatic analysis, we revealed a high nucleic acid homology between AccRACK1 and RACK1 of Apis cerana. RT-qPCR analyses demonstrated AccRACK1 was mostly expressed in 3rd instar larvae, darked-eyed pupae and adults (one and thirty days post-emergence), suggesting it might participate in the development of A. cerana cerana. Moreover, the expression of AccRACK1 was highest in the thorax, followed by the venom gland. Compared to the blank control group, AccRACK1 was induced by 24 and 44 °C, HgCl2 and pesticides (paraquat, pyridaben and methomyl) but inhibited by 14 °C, H2O2, UV light and cyhalothrin. Additionally, 0.05, 0.1, 1, 5 and 10 mg/ml PPN (juvenile hormone analogue pyriproxyfen) could promote the expression of AccRACK1, with 1 mg/ml showing the highest upregulation, suggesting it was regulated by hormones. Further study found that after knockdown of AccRACK1 by RNAi, the expression of the eukaryotic initiation factor 6 of A. cerana cerana (AcceIF6), an initiation factor regulating the initiation of translation, was inhibited, indicating AccRACK1 might affect cellular responses by translation. These findings, taken together, suggest AccRACK1 is involved in the development and responses to abiotic stresses of A. cerana cerana, and therefore, it may be of critical importance to the survival of A. cerana cerana.

A VPS13D spastic ataxia mutation disrupts the conserved adaptor-binding site in yeast Vps13

Mutations in each of the four human VPS13 (VPS13A-D) proteins are associated with distinct neurological disorders: chorea-acanthocytosis, Cohen syndrome, early-onset Parkinson's disease and spastic ataxia. Recent evidence suggests that the different VPS13 paralogs transport lipids between organelles at different membrane contact sites. How each VPS13 isoform is targeted to organelles is not known. We have shown that the localization of yeast Vps13 protein to membranes requires a conserved six-repeat region, the Vps13 Adaptor Binding (VAB) domain, which binds to organelle-specific adaptors. Here, we use a systematic mutagenesis strategy to determine the role of each repeat in recognizing each known adaptor. Our results show that mutation of invariant asparagines in repeats 1 and 6 strongly impacts the binding of all adaptors and blocks Vps13 membrane recruitment. However, we find that repeats 5-6 are sufficient for localization and interaction with adaptors. This supports a model where a single adaptor-binding site is found in the last two repeats of the VAB domain, while VAB domain repeat 1 may influence domain conformation. Importantly, a disease-causing mutation in VPS13D, which maps to the highly conserved asparagine residue in repeat 6, blocks adaptor binding and Vps13 membrane recruitment when modeled in yeast. Our findings are consistent with a conserved adaptor binding role for the VAB domain and suggest the presence of as-yet-unidentified adaptors in both yeast and humans.

WD40 domain of RqkA regulates its kinase activity and role in extraordinary radioresistance of D. radiodurans

RqkA, a DNA damage responsive serine/threonine kinase, is characterized for its role in DNA repair and cell division in D. radiodurans. It has a unique combination of a kinase domain at N-terminus and a WD40 type domain at C-terminus joined through a linker. WD40 domain is comprised of eight β-propeller repeats held together via 'tryptophan-docking motifs' and forming a typical 'velcro' closure structure. RqkA mutants lacking the WD40 region (hereafter referred to as WD mutant) could not complement RqkA loss in γ radiation resistance in D. radiodurans and lacked γ radiation-mediated activation of kinase activity in vivo. WD mutants failed to phosphorylate its cognate substrate (e.g. DrRecA) in surrogate E. coli cells. Unlike wild-type enzyme, the kinase activity of its WD40 mutants was not stimulated by pyrroloquinoline quinine (PQQ) indicating the role of the WD motifs in PQQ interaction and stimulation of its kinase activity. Together, results highlighted the importance of the WD40 domain in the regulation of RqkA kinase signaling functions in vivo, and thus, the role of WD40 domain in the regulation of any STPK is first time demonstrated in bacteria.Communicated by Ramaswamy H. Sarma.

PcG Proteins MSI1 and BMI1 Function Upstream of miR156 to Regulate Aerial Tuber Formation in Potato

Polycomb Repressive Complexes (PRC1 and PRC2) regulate developmental transitions in plants. AtBMI1, a PRC1 member, represses micro RNA156 (miR156) to trigger the onset of adult phase in Arabidopsis (Arabidopsis thaliana). miR156 overexpression (OE) reduces below-ground tuber yield, but stimulates aerial tubers in potato (Solanum tuberosum ssp andigena) under short-day (SD) photoperiodic conditions. Whether PRC members could govern tuber development through photoperiod-mediated regulation of miR156 is unknown. Here, we investigated the role of two PRC proteins, StMSI1 (PRC2 member) and StBMI1-1, in potato development. In wild-type andigena plants, StMSI1 and miR156 levels increased in stolon, whereas StBMI1-1 decreased under SD conditions. StMSI1-OE and StBMI1-1-antisense (AS) lines produced pleiotropic effects, including altered leaf architecture/compounding and reduced below-ground tuber yield. Notably, these lines showed enhanced miR156 accumulation accompanied by aerial stolons and tubers from axillary nodes, similar to miR156-OE lines. Further, grafting of StMSI1-OE or StBMI1-1-AS on wild-type stock resulted in reduced root biomass and showed increased accumulation of miR156a/b and -c precursors in the roots of wild-type stocks. RNA-sequencing of axillary nodes from StMSI1-OE and StBMI1-1-AS lines revealed downregulation of auxin and brassinosteroid genes, and upregulation of cytokinin transport/signaling genes, from 1,023 differentially expressed genes shared between the two lines. Moreover, we observed downregulation of genes encoding H2A-ubiquitin ligase and StBMI1-1/3, and upregulation of Trithorax group H3K4-methyl-transferases in StMSI1-OE Chromatin immunoprecipitation-quantitative PCR confirmed H3K27me3-mediated suppression of StBMI1-1/3, and H3K4me3-mediated activation of miR156 in StMSI1-OE plants. In summary, we show that cross talk between histone modifiers regulates miR156 and alters hormonal response during aerial tuber formation in potato under SD conditions.

The Cajal Body Protein WRAP53β Prepares the Scene for Repair of DNA Double-Strand Breaks by Regulating Local Ubiquitination

Proper repair of DNA double-strand breaks is critical for maintaining genome integrity and avoiding disease. Modification of damaged chromatin has profound consequences for the initial signaling and regulation of repair. One such modification involves ubiquitination by E3 ligases RNF8 and RNF168 within minutes after DNA double-strand break formation, altering chromatin structure and recruiting factors such as 53BP1 and BRCA1 for repair via non-homologous end-joining (NHEJ) and homologous recombination (HR), respectively. The WD40 protein WRAP53β plays an essential role in localizing RNF8 to DNA breaks by scaffolding its interaction with the upstream factor MDC1. Loss of WRAP53β impairs ubiquitination at DNA lesions and reduces downstream repair by both NHEJ and HR. Intriguingly, WRAP53β depletion attenuates repair of DNA double-strand breaks more than depletion of RNF8, indicating functions other than RNF8-mediated ubiquitination. WRAP53β plays key roles with respect to the nuclear organelles Cajal bodies, including organizing the genome to promote associated transcription and collecting factors involved in maturation of the spliceosome and telomere elongation within these organelles. It is possible that similar functions may aid also in DNA repair. Here we describe the involvement of WRAP53β in Cajal bodies and DNA double-strand break repair in detail and explore whether and how these processes may be linked. We also discuss the possibility that the overexpression of WRAP53β detected in several cancer types may reflect its normal participation in the DNA damage response rather than oncogenic properties.

Genome Wide Analysis of WD40 Proteins in Saccharomyces cerevisiae and Their Orthologs in Candida albicans

The WD40 domain containing proteins are present in the lower organisms (Monera) to higher complex metazoans with involvement in diverse cellular processes. The WD40 repeats fold into β propeller structure due to which the proteins harbouring WD40 domains function as scaffold by offering platform for interactions, bring together diverse cellular proteins to form a single complex for mediating downstream effects. Multiple functions of WD40 domain containing proteins in lower eukaryote as in Fungi have been reported with involvement in vegetative and reproductive growth, virulence etc. In this article insilico analysis of the WDR proteins in the budding yeast Saccharomyces cerevisiae was performed. By WDSP software 83 proteins in S. cerevisiae were identified with at least one WD40 motif. WD40 proteins with 6 or more WD40 motifs were considered for further studies. The WD40 proteins in yeast which are involved in various biological processes show distribution on all chromosomes (16 chromosomes in yeast) except chromosome 1. Besides the WD40 domain some of these proteins also contain other protein domains which might be responsible for the diversity in the functions of WD40 proteins in the budding yeast. These proteins in budding yeast were analysed by DAVID and Blast2Go software for functional and domains categorization. Candida albicans, an opportunistic fungal pathogen also have orthologs of these WD40 proteins with possible similar functions. This is the first time genome wide analysis of WD40 proteins in lower eukaryote i.e. budding yeast. This data may be useful in further study of the functional diversity of yeast proteomes.

RAMOSA1 ENHANCER LOCUS2-Mediated Transcriptional Repression Regulates Vegetative and Reproductive Architecture

Transcriptional repression in multicellular organisms orchestrates dynamic and precise gene expression changes that enable complex developmental patterns. Here, we present phenotypic and molecular characterization of the maize (Zea mays) transcriptional corepressor RAMOSA1 ENHANCER LOCUS2 (REL2), a unique member of the highly conserved TOPLESS (TPL) family. Analysis of single recessive mutations in rel2 revealed an array of vegetative and reproductive phenotypes, many related to defects in meristem initiation and maintenance. To better understand how REL2-mediated transcriptional complexes relate to rel2 phenotypes, we performed protein interaction assays and transcriptional profiling of mutant inflorescences, leading to the identification of different maize transcription factors and regulatory pathways that employ REL2 repression to control traits directly impacting maize yields. In addition, we used our REL2 interaction data to catalog conserved repression motifs present on REL2 interactors and showed that two of these, RLFGV- and DLN-type motifs, interact with the C-terminal WD40 domain of REL2 rather than the N terminus, which is known to bind LxLxL EAR motifs. These findings establish that the WD40 domain of TPL family proteins is an independent protein interaction surface that may work together with the N-terminal domain to allow the formation of large macromolecular complexes of functionally related transcription factors.

Ciliary proteins Fap43 and Fap44 interact with each other and are essential for proper cilia and flagella beating

Cilia beating is powered by the inner and outer dynein arms (IDAs and ODAs). These multi-subunit macrocomplexes are arranged in two rows on each outer doublet along the entire cilium length, except its distal end. To generate cilia beating, the activity of ODAs and IDAs must be strictly regulated locally by interactions with the dynein arm-associated structures within each ciliary unit and coordinated globally in time and space between doublets and along the axoneme. Here, we provide evidence of a novel ciliary complex composed of two conserved WD-repeat proteins, Fap43p and Fap44p. This complex is adjacent to another WD-repeat protein, Fap57p, and most likely the two-headed inner dynein arm, IDA I1. Loss of either protein results in altered waveform, beat stroke and reduced swimming speed. The ciliary localization of Fap43p and Fap44p is interdependent in the ciliate Tetrahymena thermophila.

PPI network analyses of human WD40 protein family systematically reveal their tendency to assemble complexes and facilitate the complex predictions

Background

WD40 repeat proteins constitute one of the largest families in eukaryotes, and widely participate in various fundamental cellular processes by interacting with other molecules. Based on individual WD40 proteins, previous work has demonstrated that their structural characteristics should confer great potential of interaction and complex formation, and has speculated that they may serve as hubs in the protein-protein interaction (PPI) network. However, what roles the whole family plays in organizing the PPI network, and whether this information can be utilized in complex prediction remain unclear. To address these issues, quantitative and systematic analyses of WD40 proteins from the perspective of PPI networks are highly required.

Results

In this work, we built two human PPI networks by using data sets with different confidence levels, and studied the network properties of the whole human WD40 protein family systematically. Our analyses have quantitatively confirmed that the human WD40 protein family, as a whole, tends to be hubs with an odds ratio of about 1.8 or greater, and the network decomposition has revealed that they are prone to enrich near the global center of the whole network with a fold change of two in the median k-values. By integrating expression profiles, we have further shown that WD40 hub proteins are inclined to be intramodular, which is indicative of complex assembling. Based on this information, we have further predicted 1674 potential WD40-associated complexes by choosing a clique-based method, which is more sensitive than others, and an indirect evaluation by co-expression scores has demonstrated its reliability.

Conclusions

At the systems level but not sporadic examples’ level, this work has provided rich knowledge for better understanding WD40 proteins’ roles in organizing the PPI network. These findings and predicted complexes can offer valuable clues for prioritizing candidates for further studies.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0567-9) contains supplementary material, which is available to authorized users.

The WD40 domain of ATG16L1 is required for its non-canonical role in lipidation of LC3 at single membranes

A hallmark of macroautophagy is the covalent lipidation of LC3 and insertion into the double-membrane phagophore, which is driven by the ATG16L1/ATG5-ATG12 complex. In contrast, non-canonical autophagy is a pathway through which LC3 is lipidated and inserted into single membranes, particularly endolysosomal vacuoles during cell engulfment events such as LC3-associated phagocytosis. Factors controlling the targeting of ATG16L1 to phagophores are dispensable for non-canonical autophagy, for which the mechanism of ATG16L1 recruitment is unknown. Here we show that the WD repeat-containing C-terminal domain (WD40 CTD) of ATG16L1 is essential for LC3 recruitment to endolysosomal membranes during non-canonical autophagy, but dispensable for canonical autophagy. Using this strategy to inhibit non-canonical autophagy specifically, we show a reduction of MHC class II antigen presentation in dendritic cells from mice lacking the WD40 CTD Further, we demonstrate activation of non-canonical autophagy dependent on the WD40 CTD during influenza A virus infection. This suggests dependence on WD40 CTD distinguishes between macroautophagy and non-canonical use of autophagy machinery.

Targeted copy number screening highlights an intragenic deletion of WDR63 as the likely cause of human occipital encephalocele and abnormal CNS development in zebrafish

Congenital malformations affecting the neural tube can present as isolated malformations or occur in association with other developmental abnormalities and syndromes. Using high-resolution copy number screening in 66 fetuses with neural tube defects, we identified six fetuses with likely pathogenic mutations, three aneuploidies (one trisomy 13 and two trisomy 18) and three deletions previously reported in NTDs (one 22q11.2 deletion and two 1p36 deletions) corresponding to 9% of the cohort. In addition, we identified five rare deletions and two duplications of uncertain significance including a rare intragenic heterozygous in-frame WDR63 deletion in a fetus with occipital encephalocele. Whole genome sequencing verified the deletion and excluded known pathogenic variants. The deletion spans exons 14-17 resulting in the expression of a protein missing the third and fourth WD-repeat domains. These findings were supported by CRISPR/Cas9-mediated somatic deletions in zebrafish. Injection of two different sgRNA-pairs targeting relevant intronic regions resulted in a deletion mimicking the human deletion and a concomitant increase of abnormal embryos with body and brain malformations (41%, n = 161 and 62%, n = 224, respectively), including a sac-like brain protrusion (7% and 9%, P < 0.01). Similar results were seen with overexpression of RNA encoding the deleted variant in zebrafish (total abnormal; 46%, n = 255, P < 0.001) compared with the overexpression of an equivalent amount of wild-type RNA (total abnormal; 3%, n = 177). We predict the in-frame WDR63 deletion to result in a dominant negative or gain-of-function form of WDR63. These are the first findings supporting a role for WDR63 in encephalocele formation.

Identifying novel small molecule antagonists for mLST8 protein using computational approaches

Mammalian lethal with SEC13 protein 8 (mLST8), is an indispensable protein subunit of mammalian target of rapamycin (mTOR) signaling pathway that interacts with the kinase domain of mTOR protein, thereby stabilizing its active site. Experimental studies reported the over expression of mLST8 in human colon and prostate cancers by activation of both mTORC1/2 complexes and subsequent downstream substrates leading to tumor progression. Considering its role, targeting mLST8 protein would be a therapeutic approach against tumor progression in colon and prostate cancers. Hence, using in silico structure based drug design approach, the comparative binding patterns of 1,1'-binapthyl-2,2'diol (BINOL), 1-(2-carboxynaphth-1yl)-2-naphthoic acid (SCF-12) and their analogs in the cavity of mLST8 were explored. ADME and binding energy calculations led to the identification of five compounds with favorable Glide (G) scores and implicated the importance of Asn132 and Gln225 as key binding residues. Molecular dynamics (MD) simulations and free energy landscape (FEL) approaches helped in elucidating the binding mechanism and suggested the possibility of ligands 1-3 namely, ZINC01765622, ZINC62723702 and ZINC02576980 to be promising antagonists for mLST8. Thus, this study substantiates the prospect of targeting mLST8 protein using potent hits which could hinder tumor progression in colon and prostate cancers.

Prokaryotic and Highly-Repetitive WD40 Proteins: A Systematic Study

As an ancient protein family, the WD40 repeat proteins often play essential roles in fundamental cellular processes in eukaryotes. Although investigations of eukaryotic WD40 proteins have been frequently reported, prokaryotic ones remain largely uncharacterized. In this paper, we report a systematic analysis of prokaryotic WD40 proteins and detailed comparisons with eukaryotic ones. About 4,000 prokaryotic WD40 proteins have been identified, accounting for 6.5% of all WD40s. While their abundances are less than 0.1% in most prokaryotes, they are enriched in certain species from Cyanobacteria and Planctomycetes, and participate in various functions such as prokaryotic signal transduction and nutrient synthesis. Comparisons show that a higher proportion of prokaryotic WD40s tend to contain multiple WD40 domains and a large number of hydrogen bond networks. The observation that prokaryotic WD40 proteins tend to show high internal sequence identity suggests that a substantial proportion of them (~20%) should be formed by recent or young repeat duplication events. Further studies demonstrate that the very young WD40 proteins, i.e., Highly-Repetitive WD40s, should be of higher stability. Our results have presented a catalogue of prokaryotic WD40 proteins, and have shed light on their evolutionary origins.

Remarkable Evolutionary Conservation of Antiobesity ADIPOSE/WDTC1 Homologs in Animals and Plants

ASG2 (Altered Seed Germination 2) is a prenylated protein in Arabidopsis thaliana that participates to abscisic acid signaling and is proposed to act as a substrate adaptor for the DDB1 (DNA damage-binding protein 1)-CUL4 (Cullin 4) E3 ubiquitin ligase complex. ASG2 harbors WD40 and TetratricoPeptide Repeat (TPR) domains, and resembles the well-conserved animal gene called ADP (antiobesity factor ADIPOSE) in fly and WDTC1 (WD40 and TPR 1) in humans. Loss of function of WDTC1 results in an increase in adipocytes, fat accumulation, and obesity. Antiadipogenic functions of WDTC1 involve regulation of fat-related gene transcription, notably through its binding to histone deacetylases (HDACs). Our sequence and phylogenetic analysis reveals that ASG2 belongs to the ADP/WDTC1 cluster. ASG2 and WDTC1 share a highly conserved organization that encompasses structural and functional motifs: seven WD40 domains and WD40 hotspot-related residues, three TPR protein-protein interaction domains, DDB1-binding elements [H-box and DWD (DDB1-binding WD40 protein)-box], and a prenylatable C-terminus. Furthermore, ASG2 involvement in fat metabolism was confirmed by reverse genetic approaches using asg2 knockout Arabidopsis plants. Under limited irradiance, asg2 mutants produce "obese" seeds characterized by increased weight, oil body density, and higher fatty acid contents. In addition, considering some ASG2- and WDTC1-peculiar properties, we show that the WDTC1 C-terminus is prenylated in vitro and HDAC-binding capability is conserved in ASG2, suggesting that the regulation mechanism and targets of ADP/WDTC1-like proteins may be conserved features. Our findings reveal the remarkable evolutionary conservation of the structure and the physiological role of ADIPOSE homologs in animals and plants.

Missense mutations in the WD40 domain of AHI1 cause non-syndromic retinitis pigmentosa

Background

Recent findings suggesting that Abelson helper integration site 1 (AHI1) is involved in non-syndromic retinal disease have been debated, as the functional significance of identified missense variants was uncertain. We assessed whether AHI1 variants cause non-syndromic retinitis pigmentosa (RP).

Methods

Exome sequencing was performed in three probands with RP. The effects of the identified missense variants in AHI1 were predicted by three-dimensional structure homology modelling. Ciliary parameters were evaluated in patient’s fibroblasts, and recombinant mutant proteins were expressed in ciliated retinal pigmented epithelium cells.

Results

In the three patients with RP, three sets of compound heterozygous variants were detected in AHI1 (c.2174G>A; p.Trp725* and c.2258A>T; p.Asp753Val, c.660delC; p.Ser221Glnfs*10 and c.2090C>T; p.Pro697Leu, c.2087A>G; p.His696Arg and c.2429C>T; p.Pro810Leu). All four missense variants were present in the conserved WD40 domain of Jouberin, the ciliary protein encoded by AHI1, with variable predicted implications for the domain structure. No significant changes in the percentage of ciliated cells, nor in cilium length or intraflagellar transport were detected. However, expression of mutant recombinant Jouberin in ciliated cells showed a significantly decreased enrichment at the ciliary base.

Conclusions

This report confirms that mutations in AHI1 can underlie autosomal recessive RP. Moreover, it structurally and functionally validates the effect of the RP-associated AHI1 variants on protein function, thus proposing a new genotype–phenotype correlation for AHI1 mutation associated retinal ciliopathies.

Planar Cell Polarity Effector Fritz Interacts with Dishevelled and Has Multiple Functions in Regulating PCP

The Planar cell Polarity Effector (PPE) genes inturned, fuzzy, and fritz are downstream components in the frizzled/starry night signaling pathway, and their function is instructed by upstream Planar Cell Polarity (PCP) core genes such as frizzled and dishevelled. PPE proteins accumulate asymmetrically in wing cells and function in a protein complex mediated by direct interactions between In and Frtz and In and Fy. How the PCP proteins instruct the accumulation of PPE protein is unknown. We found a likely direct interaction between Dishevelled and Fritz and Dishevelled and Fuzzy that could play a role in this. We previously found that mild overexpression of frtz rescued a weak in allele. To determine if this was due to extra Frtz stabilizing mutant In or due to Frtz being able to bypass the need for In we generate a precise deletion of the inturned gene (in^PD). We found that mild overexpression of Fritz partially rescued in^PD, indicating that fritz has In independent activity in PCP. Previous studies of PPE proteins used fixed tissues, and did not provide any insights into the dynamic properties of PPE proteins. We used CRISPR/Cas9 genome editing technology to edit the fritz gene to add a green fluorescent protein tag. fritz^mNeonGreen provides complete rescue activity and works well for in vivo imaging. Our data showed that Fritz is very dynamic in epidermal cells and preferentially distributed to discrete membrane subdomains (“puncta”). Surprisingly, we found it in stripes in developing bristles.

Structure of the MeCP2-TBLR1 complex reveals a molecular basis for Rett syndrome and related disorders

Rett syndrome (RTT) is an X-linked neurological disorder caused by mutations in the methyl-CpG-binding protein 2 (MeCP2) gene. The majority of RTT missense mutations disrupt the interaction of the MeCP2 with DNA or the nuclear receptor corepressor (NCoR)/silencing mediator of retinoic acid and thyroid receptors (SMRT) corepressor complex. Here, we show that the "NCoR/SMRT interaction domain" (NID) of MeCP2 directly contacts transducin beta-like 1 (TBL1) and TBL1 related (TBLR1), two paralogs that are core components of NCoR/SMRT. We determine the cocrystal structure of the MeCP2 NID in complex with the WD40 domain of TBLR1 and confirm by in vitro and ex vivo assays that mutation of interacting residues of TBLR1 and TBL1 disrupts binding to MeCP2. Strikingly, the four MeCP2-NID residues mutated in RTT are those residues that make the most extensive contacts with TBLR1. Moreover, missense mutations in the gene for TBLR1 that are associated with intellectual disability also prevent MeCP2 binding. Our study therefore reveals the molecular basis of an interaction that is crucial for optimal brain function.

Direct Calculation of Protein Fitness Landscapes through Computational Protein Design

Naturally selected amino-acid sequences or experimentally derived ones are often the basis for understanding how protein three-dimensional conformation and function are determined by primary structure. Such sequences for a protein family comprise only a small fraction of all possible variants, however, representing the fitness landscape with limited scope. Explicitly sampling and characterizing alternative, unexplored protein sequences would directly identify fundamental reasons for sequence robustness (or variability), and we demonstrate that computational methods offer an efficient mechanism toward this end, on a large scale. The dead-end elimination and A(∗) search algorithms were used here to find all low-energy single mutant variants, and corresponding structures of a G-protein heterotrimer, to measure changes in structural stability and binding interactions to define a protein fitness landscape. We established consistency between these algorithms with known biophysical and evolutionary trends for amino-acid substitutions, and could thus recapitulate known protein side-chain interactions and predict novel ones.

Recurrent somatic mutations in POLR2A define a distinct subset of meningiomas

RNA polymerase II mediates the transcription of all protein-coding genes in eukaryotic cells, a process that is fundamental to life. Genomic mutations altering this enzyme have not previously been linked to any pathology in humans, which is a testament to its indispensable role in cell biology. On the basis of a combination of next-generation genomic analyses of 775 meningiomas, we report that recurrent somatic p.Gln403Lys or p.Leu438_His439del mutations in POLR2A, which encodes the catalytic subunit of RNA polymerase II (ref. 1), hijack this essential enzyme and drive neoplasia. POLR2A mutant tumors show dysregulation of key meningeal identity genes, including WNT6 and ZIC1/ZIC4. In addition to mutations in POLR2A, NF2, SMARCB1, TRAF7, KLF4, AKT1, PIK3CA, and SMO, we also report somatic mutations in AKT3, PIK3R1, PRKAR1A, and SUFU in meningiomas. Our results identify a role for essential transcriptional machinery in driving tumorigenesis and define mutually exclusive meningioma subgroups with distinct clinical and pathological features.

A WDR Gene Is a Conserved Member of a Chitin Synthase Gene Cluster and Influences the Cell Wall in Aspergillus nidulans

WD40 repeat (WDR) proteins are pleiotropic molecular hubs. We identify a WDR gene that is a conserved genomic neighbor of a chitin synthase gene in Ascomycetes. The WDR gene is unique to fungi and plants, and was called Fungal Plant WD (FPWD). FPWD is within a cell wall metabolism gene cluster in the Ascomycetes (Pezizomycotina) comprising chsD, a Chs activator and a GH17 glucanase. The FPWD, AN1556.2 locus was deleted in Aspergillus nidulans strain SAA.111 by gene replacement and only heterokaryon transformants were obtained. The re-annotation of Aspergilli genomes shows that AN1556.2 consists of two tightly linked separate genes, i.e., the WDR gene and a putative beta-flanking gene of unknown function. The WDR and the beta-flanking genes are conserved genomic neighbors localized within a recently identified metabolic cell wall gene cluster in genomes of Aspergilli. The heterokaryons displayed increased susceptibility to drugs affecting the cell wall, and their phenotypes, observed by optical, confocal, scanning electron and atomic force microscopy, suggest cell wall alterations. Quantitative real-time PCR shows altered expression of some cell wall-related genes. The possible implications on cell wall biosynthesis are discussed.

The adaptor protein DCAF7 mediates the interaction of the adenovirus E1A oncoprotein with the protein kinases DYRK1A and HIPK2

DYRK1A is a constitutively active protein kinase that has a critical role in growth and development which functions by regulating cell proliferation, differentiation and survival. DCAF7 (also termed WDR68 or HAN11) is a cellular binding partner of DYRK1A and also regulates signalling by the protein kinase HIPK2. DCAF7 is an evolutionarily conserved protein with a single WD40 repeat domain and has no catalytic activity. We have defined a DCAF7 binding motif of 12 amino acids in the N-terminal domain of class 1 DYRKs that is functionally conserved in DYRK1 orthologs from Xenopus, Danio rerio and the slime mold Dictyostelium discoideum. A similar sequence was essential for DCAF7 binding to HIPK2, whereas the closely related HIPK1 family member did not bind DCAF7. Immunoprecipitation and pulldown experiments identified DCAF7 as an adaptor for the association of the adenovirus E1A protein with DYRK1A and HIPK2. Furthermore, DCAF7 was required for the hyperphosphorylation of E1A in DYRK1A or HIPK2 overexpressing cells. Our results characterize DCAF7 as a substrate recruiting subunit of DYRK1A and HIPK2 and suggest that it is required for the negative effect of DYRK1A on E1A-induced oncogenic transformation.

Integrating population variation and protein structural analysis to improve clinical interpretation of missense variation: application to the WD40 domain

We present a generic, multidisciplinary approach for improving our understanding of novel missense variants in recently discovered disease genes exhibiting genetic heterogeneity, by combining clinical and population genetics with protein structural analysis. Using six new de novo missense diagnoses in TBL1XR1 from the Deciphering Developmental Disorders study, together with population variation data, we show that the β-propeller structure of the ubiquitous WD40 domain provides a convincing way to discriminate between pathogenic and benign variation. Children with likely pathogenic mutations in this gene have severely delayed language development, often accompanied by intellectual disability, autism, dysmorphology and gastrointestinal problems. Amino acids affected by likely pathogenic missense mutations are either crucial for the stability of the fold, forming part of a highly conserved symmetrically repeating hydrogen-bonded tetrad, or located at the top face of the β-propeller, where 'hotspot' residues affect the binding of β-catenin to the TBLR1 protein. In contrast, those altered by population variation are significantly less likely to be spatially clustered towards the top face or to be at buried or highly conserved residues. This result is useful not only for interpreting benign and pathogenic missense variants in this gene, but also in other WD40 domains, many of which are associated with disease.

Recruitment of Saccharomyces cerevisiae Cmr1/Ydl156w to Coding Regions Promotes Transcription Genome Wide

Cmr1 (changed mutation rate 1) is a largely uncharacterized nuclear protein that has recently emerged in several global genetic interaction and protein localization studies. It clusters with proteins involved in DNA damage and replication stress response, suggesting a role in maintaining genome integrity. Under conditions of proteasome inhibition or replication stress, this protein localizes to distinct sub-nuclear foci termed as intranuclear quality control (INQ) compartments, which sequester proteins for their subsequent degradation. Interestingly, it also interacts with histones, chromatin remodelers and modifiers, as well as with proteins involved in transcription including subunits of RNA Pol I and Pol III, but not with those of Pol II. It is not known whether Cmr1 plays a role in regulating transcription of Pol II target genes. Here, we show that Cmr1 is recruited to the coding regions of transcribed genes of S. cerevisiae. Cmr1 occupancy correlates with the Pol II occupancy genome-wide, indicating that it is recruited to coding sequences in a transcription-dependent manner. Cmr1-enriched genes include Gcn4 targets and ribosomal protein genes. Furthermore, our results show that Cmr1 recruitment to coding sequences is stimulated by Pol II CTD kinase, Kin28, and the histone deacetylases, Rpd3 and Hos2. Finally, our genome-wide analyses implicate Cmr1 in regulating Pol II occupancy at transcribed coding sequences. However, it is dispensable for maintaining co-transcriptional histone occupancy and histone modification (acetylation and methylation). Collectively, our results show that Cmr1 facilitates transcription by directly engaging with transcribed coding regions.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

Genome-wide identification, sequence characterization, and protein-protein interaction properties of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat family members in Solanum lycopersicum

MAIN CONCLUSIONS

One hundred DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family genes were identified in the S. lycopersicum genome. The DWD genes encode proteins presumably functioning as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. These findings provide candidate genes and a research platform for further gene functionality and molecular breeding study. A subclass of DDB1 (damaged DNA binding protein-1)-binding WD40-repeat domain (DWD) family proteins has been demonstrated to function as the substrate recognition subunits of the cullin4-ring ubiquitin E3 ligase complex. However, little information is available about the cognate subfamily genes in tomato (S. lycopersicum). In this study, based on the recently released tomato genome sequences, 100 tomato genes encoding DWD proteins that potentially interact with DDB1 were identified and characterized, including analyses of the detailed annotations, chromosome locations and compositions of conserved amino acid domains. In addition, a phylogenetic tree, which comprises of three main groups, of the subfamily genes was constructed. The physical interaction between tomato DDB1 and 14 representative DWD proteins was determined by yeast two-hybrid and co-immunoprecipitation assays. The subcellular localization of these 14 representative DWD proteins was determined. Six of them were localized in both nucleus and cytoplasm, seven proteins exclusively in cytoplasm, and one protein either in nucleus and cytoplasm, or exclusively in cytoplasm. Comparative genomic analysis demonstrated that the expansion of these subfamily members in tomato predominantly resulted from two whole-genome triplication events in the evolution history.

The carboxy-terminal domain of Erb1 is a seven-bladed ß-propeller that binds RNA

Erb1 (Eukaryotic Ribosome Biogenesis 1) protein is essential for the maturation of the ribosomal 60S subunit. Functional studies in yeast and mammalian cells showed that altogether with Nop7 and Ytm1 it forms a stable subcomplex called PeBoW that is crucial for a correct rRNA processing. The exact function of the protein within the process remains unknown. The N-terminal region of the protein includes a well conserved region shown to be involved in PeBoW complex formation whereas the carboxy-terminal half was predicted to contain seven WD40 repeats. This first structural report on Erb1 from yeast describes the architecture of a seven-bladed β-propeller domain that revealed a characteristic extra motif formed by two α-helices and a β-strand that insert within the second WD repeat. We performed analysis of molecular surface and crystal packing, together with multiple sequence alignment and comparison of the structure with other β-propellers, in order to identify areas that are more likely to mediate protein-protein interactions. The abundance of many positively charged residues on the surface of the domain led us to investigate whether the propeller of Erb1 might be involved in RNA binding. Three independent assays confirmed that the protein interacted in vitro with polyuridilic acid (polyU), thus suggesting a possible role of the domain in rRNA rearrangement during ribosome biogenesis.

Cmr1/WDR76 defines a nuclear genotoxic stress body linking genome integrity and protein quality control

DNA replication stress is a source of genomic instability. Here we identify changed mutation rate 1 (Cmr1) as a factor involved in the response to DNA replication stress in Saccharomyces cerevisiae and show that Cmr1--together with Mrc1/Claspin, Pph3, the chaperonin containing TCP1 (CCT) and 25 other proteins--define a novel intranuclear quality control compartment (INQ) that sequesters misfolded, ubiquitylated and sumoylated proteins in response to genotoxic stress. The diversity of proteins that localize to INQ indicates that other biological processes such as cell cycle progression, chromatin and mitotic spindle organization may also be regulated through INQ. Similar to Cmr1, its human orthologue WDR76 responds to proteasome inhibition and DNA damage by relocalizing to nuclear foci and physically associating with CCT, suggesting an evolutionarily conserved biological function. We propose that Cmr1/WDR76 plays a role in the recovery from genotoxic stress through regulation of the turnover of sumoylated and phosphorylated proteins.

Loss-of-function mutations in WDR73 are responsible for microcephaly and steroid-resistant nephrotic syndrome: Galloway-Mowat syndrome

Galloway-Mowat syndrome is a rare autosomal-recessive condition characterized by nephrotic syndrome associated with microcephaly and neurological impairment. Through a combination of autozygosity mapping and whole-exome sequencing, we identified WDR73 as a gene in which mutations cause Galloway-Mowat syndrome in two unrelated families. WDR73 encodes a WD40-repeat-containing protein of unknown function. Here, we show that WDR73 was present in the brain and kidney and was located diffusely in the cytoplasm during interphase but relocalized to spindle poles and astral microtubules during mitosis. Fibroblasts from one affected child and WDR73-depleted podocytes displayed abnormal nuclear morphology, low cell viability, and alterations of the microtubule network. These data suggest that WDR73 plays a crucial role in the maintenance of cell architecture and cell survival. Altogether, WDR73 mutations cause Galloway-Mowat syndrome in a particular subset of individuals presenting with late-onset nephrotic syndrome, postnatal microcephaly, severe intellectual disability, and homogenous brain MRI features. WDR73 is another example of a gene involved in a disease affecting both the kidney glomerulus and the CNS.

WDSPdb: a database for WD40-repeat proteins

WD40-repeat proteins, as one of the largest protein families, often serve as platforms to assemble functional complexes through the hotspot residues on their domain surfaces, and thus play vital roles in many biological processes. Consequently, it is highly required for researchers who study WD40 proteins and protein-protein interactions to obtain structural information of WD40 domains. Systematic identification of WD40-repeat proteins, including prediction of their secondary structures, tertiary structures and potential hotspot residues responsible for protein-protein interactions, may constitute a valuable resource upon this request. To achieve this goal, we developed a specialized database WDSPdb (http://wu.scbb.pkusz.edu.cn/wdsp/) to provide these details of WD40-repeat proteins based on our recently published method WDSP. The WDSPdb contains 63,211 WD40-repeat proteins identified from 3383 species, including most well-known model organisms. To better serve the community, we implemented a user-friendly interactive web interface to browse, search and download the secondary structures, 3D structure models and potential hotspot residues provided by WDSPdb.

Mutation of POC1B in a severe syndromic retinal ciliopathy

We describe a consanguineous Iraqi family with Leber congenital amaurosis (LCA), Joubert syndrome (JBTS), and polycystic kidney disease (PKD). Targeted next-generation sequencing for excluding mutations in known LCA and JBTS genes, homozygosity mapping, and whole-exome sequencing identified a homozygous missense variant, c.317G>C (p.Arg106Pro), in POC1B, a gene essential for ciliogenesis, basal body, and centrosome integrity. In silico modeling suggested a requirement of p.Arg106 for the formation of the third WD40 repeat and a protein interaction interface. In human and mouse retina, POC1B localized to the basal body and centriole adjacent to the connecting cilium of photoreceptors and in synapses of the outer plexiform layer. Knockdown of Poc1b in zebrafish caused cystic kidneys and retinal degeneration with shortened and reduced photoreceptor connecting cilia, compatible with the human syndromic ciliopathy. A recent study describes homozygosity for p.Arg106ProPOC1B in a family with nonsyndromic cone-rod dystrophy. The phenotype associated with homozygous p.Arg106ProPOC1B may thus be highly variable, analogous to homozygous p.Leu710Ser in WDR19 causing either isolated retinitis pigmentosa or Jeune syndrome. Our study indicates that POC1B is required for retinal integrity, and we propose POC1B mutations as a probable cause for JBTS with severe PKD.

Direct interaction between the WD40 repeat protein WDR-23 and SKN-1/Nrf inhibits binding to target DNA

SKN-1/Nrf transcription factors activate cytoprotective genes in response to reactive small molecules and strongly influence stress resistance, longevity, and development. The molecular mechanisms of SKN-1/Nrf regulation are poorly defined. We previously identified the WD40 repeat protein WDR-23 as a repressor of Caenorhabditis elegans SKN-1 that functions with a ubiquitin ligase to presumably target the factor for degradation. However, SKN-1 activity and nuclear accumulation are not always correlated, suggesting that there could be additional regulatory mechanisms. Here, we integrate forward genetics and biochemistry to gain insights into how WDR-23 interacts with and regulates SKN-1. We provide evidence that WDR-23 preferentially regulates one of three SKN-1 variants through a direct interaction that is required for normal stress resistance and development. Homology modeling predicts that WDR-23 folds into a β-propeller, and we identify the top of this structure and four motifs at the termini of SKN-1c as essential for the interaction. Two of these SKN-1 motifs are highly conserved in human Nrf1 and Nrf2 and two directly interact with target DNA. Lastly, we demonstrate that WDR-23 can block the ability of SKN-1c to interact with DNA sequences of target promoters identifying a new mechanism of regulation that is independent of the ubiquitin proteasome system, which can become occupied with damaged proteins during stress.

β-Bulges: extensive structural analyses of β-sheets irregularities

β-Sheets are quite frequent in protein structures and are stabilized by regular main-chain hydrogen bond patterns. Irregularities in β-sheets, named β-bulges, are distorted regions between two consecutive hydrogen bonds. They disrupt the classical alternation of side chain direction and can alter the directionality of β-strands. They are implicated in protein-protein interactions and are introduced to avoid β-strand aggregation. Five different types of β-bulges are defined. Previous studies on β-bulges were performed on a limited number of protein structures or one specific family. These studies evoked a potential conservation during evolution. In this work, we analyze the β-bulge distribution and conservation in terms of local backbone conformations and amino acid composition. Our dataset consists of 66 times more β-bulges than the last systematic study (Chan et al. Protein Science 1993, 2:1574-1590). Novel amino acid preferences are underlined and local structure conformations are highlighted by the use of a structural alphabet. We observed that β-bulges are preferably localized at the N- and C-termini of β-strands, but contrary to the earlier studies, no significant conservation of β-bulges was observed among structural homologues. Displacement of β-bulges along the sequence was also investigated by Molecular Dynamics simulations.

A method for WD40 repeat detection and secondary structure prediction

WD40-repeat proteins (WD40s), as one of the largest protein families in eukaryotes, play vital roles in assembling protein-protein/DNA/RNA complexes. WD40s fold into similar β-propeller structures despite diversified sequences. A program WDSP (WD40 repeat protein Structure Predictor) has been developed to accurately identify WD40 repeats and predict their secondary structures. The method is designed specifically for WD40 proteins by incorporating both local residue information and non-local family-specific structural features. It overcomes the problem of highly diversified protein sequences and variable loops. In addition, WDSP achieves a better prediction in identifying multiple WD40-domain proteins by taking the global combination of repeats into consideration. In secondary structure prediction, the average Q3 accuracy of WDSP in jack-knife test reaches 93.7%. A disease related protein LRRK2 was used as a representive example to demonstrate the structure prediction.