The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins

Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum)

BACKGROUND

WD40 proteins, which are highly prevalent in eukaryotes, play important roles in plant development and stress responses. However, systematic identification and exploration of WD40 proteins in tobacco have not yet been conducted.

RESULTS

In this study, a total of 399 WD40 regulatory genes were identified in common tobacco (Nicotiana tabacum). Gene structure and motif analysis revealed structural and functional diversity among different clades of tobacco WD40 regulatory genes. The expansion of tobacco WD40 regulatory genes was mainly driven by segmental duplication and purifying selection. A potential regulatory network of NtWD40s suggested that NtWD40s might be regulated by miRNAs and transcription factors in various biological processes. Expression pattern analysis via transcriptome analysis and qRT-PCR revealed that many NtWD40s exhibited tissue-specific expression patterns and might be involved in various biotic and abiotic stresses. Furthermore, we have validated the critical role of NtTTG1, which was located in the nuclei of trichome cells, in enhancing the drought tolerance of tobacco plants.

CONCLUSIONS

Our study provides comprehensive information to better understand the evolution of WD40 regulatory genes and their roles in different stress responses in tobacco.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

Structure and stability of the designer protein WRAP-T and its permutants

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-Propeller proteins are common natural disc-like pseudo-symmetric proteins that contain multiple repeats (‘blades’) each consisting of a 4-stranded anti-parallel \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-sheet. So far, 4- to 12-bladed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-propellers have been discovered in nature showing large functional and sequential variation. Using computational design approaches, we created perfectly symmetric \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-propellers out of natural pseudo-symmetric templates. These proteins are useful tools to study protein evolution of this very diverse fold. While the 7-bladed architecture is the most common, no symmetric 7-bladed monomer has been created and characterized so far. Here we describe such a engineered protein, based on a highly symmetric natural template, and test the effects of circular permutation on its stability. Geometrical analysis of this protein and other artificial symmetrical proteins reveals no systematic constraint that could be used to help in engineering of this fold, and suggests sequence constraints unique to each \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-propeller sub-family.

Kinetics of the multitasking high-affinity Win binding site of WDR5 in restricted and unrestricted conditions

Recent advances in quantitative proteomics show that WD40 proteins play a pivotal role in numerous cellular networks. Yet, they have been fairly unexplored and their physical associations with other proteins are ambiguous. A quantitative understanding of these interactions has wide-ranging significance. WD40 repeat protein 5 (WDR5) interacts with all members of human SET1/MLL methyltransferases, which regulate methylation of the histone 3 lysine 4 (H3K4). Here, using real-time binding measurements in a high-throughput setting, we identified the kinetic fingerprint of transient associations between WDR5 and 14-residue WDR5 interaction (Win) motif peptides of each SET1 protein (SET1Win). Our results reveal that the high-affinity WDR5-SET1Win interactions feature slow association kinetics. This finding is likely due to the requirement of SET1Win to insert into the narrow WDR5 cavity, also named the Win binding site. Furthermore, our explorations indicate fairly slow dissociation kinetics. This conclusion is in accordance with the primary role of WDR5 in maintaining the functional integrity of a large multisubunit complex, which regulates the histone methylation. Because the Win binding site is considered a key therapeutic target, the immediate outcomes of this study could form the basis for accelerated developments in medical biotechnology.

Increased WD40 motifs in Planctomycete bacteria and their evolutionary relevance

Species of the family Planctomycetes have a complex intracellular structure, which is distinct from that of the majority of non-Planctomycetes bacteria. At present, genomic evidence of the evolution of intracellular complexity is lacking, cognitions of Planctomycetes's intracellular structure mainly rely on electron microscope observation. As the presence of WD40 motifs in eukaryotic proteins probably links to intracellular complexity, bioinformatic studies were conducted to detect and enumerate WD40 motifs, WD40 domains, and WD40 motif-bearing proteins in the genomes of 11 Planctomycetes species, 2775 non-Planctomycetes bacteria, and 63 representative eukaryotes. Compared to non-Planctomycetes bacteria (average 5 WD40 motifs and 1 WD40 motif-bearing protein per genome), a large increase in the number of WD40 motifs in Planctomycetes species (average 116 WD40 motifs and 26 WD40 motif-bearing proteins per genome) was observed. However, the average number of WD40 motifs in Planctomycetes species was significantly lower than that of eukaryotes (average 584 WD40 motifs and 193 WD40 motif-bearing proteins per genome). The number of WD40 motif-bearing proteins was found to correlate with genome size and gene number. Most WD40 motif-bearing proteins of Planctomycetes species belonged to the categories of 'ribosome assembly protein 4' and 'eukaryotic-like serine/threonine protein kinase.' Collinearity analysis of amino acid compositions of Planctomycetes and eukaryotic WD40 motifs revealed that the sequences of the four anti-parallel β-sheets of WD40 motifs were conserved. However, a number of Planctomycetes WD40 motifs had increased size of the interval region of β-sheets D and A. Taken together, results of this study suggest a positive correlation between the number of WD40 motif-bearing proteins and the evolution of Planctomycetes species toward a complex intracellular structure similar to that of eukaryotes.

Cloning and characterization of the target protein subunit lst8 of rapamycin in Apostichopus japonicus

Autophagy plays an important role in the immune defense systems of vertebrates through the interaction between the lethal with SEC13 protein 8 (lst8) and the mechanistic target of rapamycin. In the present study, a novel invertebrate lst8 homologue is identified from Apostichopus japonicus (designated as Ajlst8) via polymerase chain reaction. The full-length complementary DNA of Ajlst8 comprises a 5'-untranslated region (UTR) of 78 base pair (bp), a 3'-UTR of 479 bp, and a putative open reading frame of 951 bp; hence, 316 amino acids are encoded. Structural analysis shows that the deduced amino acid of Ajlst8 shares six typical WD40 domains (28 aa-248 aa). Spatial expression analysis indicates that Ajlst8 is ubiquitously expressed in all the examined tissues, with a larger magnitude in coelomocytes. Vibrio splendidus infection in vivo and lipopolysaccharide exposure in vitro can significantly upregulate the messenger RNA expression of Ajlst8 by 2.39-fold and 1.93-fold compared with the control group, respectively. LPS exposure could also significantly induced the protein level of Ajlst8 to 2.38-fold and the autophagy level was markedly increased by 3.08-fold under same condition. The RNA interference of Ajlst8 in primary coelomocytes also reduces the relative expression of autophagy with a 0.71-fold decrease in the ratio of LC3-II/LC3-I compared with that in the control group. These results indicate that Ajlst8 is a novel immune regulator that may be involved in the antibacterial response process of sea cucumber by regulating autophagy.

TRIP-1 in the extracellular matrix promotes nucleation of calcium phosphate polymorphs

In search for bone and dentin extracellular matrix (ECM) proteins, transforming growth factor beta receptor II interacting protein 1 (TRIP-1) was identified as a novel protein synthesized by osteoblasts and odontoblasts and exported to the ECM. TRIP-1 is a WD-40 (WD is Tryptophan-Aspartic acid dipeptide) protein that has been well recognized for its physiological role in the endoplasmic reticulum (ER). In the ER, TRIP-1 functions as an essential subunit of eukaryotic elongation initiation factor 3 and is involved in the protein translational machinery. Recently, we reported that TRIP-1 is localized in the ECM of bone and dentin. In this study, we demonstrate that varying concentrations of TRIP-1 can participate in the nucleation of calcium phosphate polymorphs. Nucleation studies performed with high calcium and phosphate concentration demonstrated that recombinant TRIP-1 could orchestrate the formation of hydroxyapatite crystals. Nucleation experiments performed on demineralized and deproteinized dentin wafer under physiological conditions and subsequent transmission electron microscope analysis of the deposits at the end of 7 and 14 days showed that TRIP-1 promoted the deposition of calcium phosphate mineral aggregates in the gap-overlap region of type I collagen. Taken together, we provide mechanistic insight into the role of this intracellular protein in matrix mineralization.

Genome-wide identification and analysis of WD40 proteins in wheat (Triticum aestivum L.)

BACKGROUND

WD40 domains are abundant in eukaryotes, and they are essential subunits of large multiprotein complexes, which serve as scaffolds. WD40 proteins participate in various cellular processes, such as histone modification, transcription regulation, and signal transduction. WD40 proteins are regarded as crucial regulators of plant development processes. However, the systematic identification and analysis of WD40 proteins have yet to be reported in wheat.

RESULTS

In this study, a total of 743 WD40 proteins were identified in wheat, and they were grouped into 5 clusters and 11 subfamilies. Their gene structures, chromosomal locations, and evolutionary relationships were analyzed. Among them, 39 and 46 pairs of TaWD40s were distinguished as tandem duplication and segmental duplication genes. The 123 OsWD40s were identified to exhibit synteny with TaWD40s. TaWD40s showed the specific characteristics at the reproductive developmental stage, and numerous TaWD40s were involved in responses to stresses, including cold, heat, drought, and powdery mildew infection pathogen, based on the result of RNA-seq data analysis. The expression profiles of some TaWD40s in wheat seed development were confirmed through qRT-PCR technique.

CONCLUSION

In this study, 743 TaWD40s were identified from the wheat genome. As the main driving force of evolution, duplication events were observed, and homologous recombination was another driving force of evolution. The expression profiles of TaWD40s revealed their importance for the growth and development of wheat and their response to biotic and abiotic stresses. Our study also provided important information for further functional characterization of some WD40 proteins in wheat.

Biophysical and structural characterization of the thermostable WD40 domain of a prokaryotic protein, Thermomonospora curvata PkwA

WD40 proteins belong to a big protein family with members identified in every eukaryotic proteome. However, WD40 proteins were only reported in a few prokaryotic proteomes. Using WDSP (http://wu.scbb.pkusz.edu.cn/wdsp/), a prediction tool, we identified thousands of prokaryotic WD40 proteins, among which few proteins have been biochemically characterized. As shown in our previous bioinformatics study, a large proportion of prokaryotic WD40 proteins have higher intramolecular sequence identity among repeats and more hydrogen networks, which may indicate better stability than eukaryotic WD40s. Here we report our biophysical and structural study on the WD40 domain of PkwA from Thermomonospora curvata (referred as tPkwA-C). We demonstrated that the stability of thermophilic tPkwA-C correlated to ionic strength and tPkwA-C exhibited fully reversible unfolding under different denaturing conditions. Therefore, the folding kinetics was also studied through stopped-flow circular dichroism spectra. The crystal structure of tPkwA-C was further resolved and shed light on the key factors that stabilize its beta-propeller structure. Like other WD40 proteins, DHSW tetrad has a significant impact on the stability of tPkwA-C. Considering its unique features, we proposed that tPkwA-C should be a great structural template for protein engineering to study key residues involved in protein-protein interaction of a WD40 protein.

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.

The Role of PALB2 in the DNA Damage Response and Cancer Predisposition

The deoxyribonucleic acid (DNA) damage response (DDR) is a major feature in the maintenance of genome integrity and in the suppression of tumorigenesis. PALB2 (Partner and Localizer of Breast Cancer 2 (BRCA2)) plays an important role in maintaining genome integrity through its role in the Fanconi anemia (FA) and homologous recombination (HR) DNA repair pathways. Since its identification as a BRCA2 interacting partner, PALB2 has emerged as a pivotal tumor suppressor protein associated to hereditary cancer susceptibility to breast and pancreatic cancers. In this review, we discuss how other DDR proteins (such as the kinases Ataxia Telangiectasia Mutated (ATM) and ATM- and Rad3-Related (ATR), mediators BRCA1 (Breast Cancer 1)/BRCA2 and effectors RAD51/DNA Polymerase η (Polη) interact with PALB2 to orchestrate DNA repair. We also examine the involvement of PALB2 mutations in the predisposition to cancer and the role of PALB2 in stimulating error-free DNA repair through the FA/HR pathway.

RNAi-mediated silencing of the Arabidopsis thaliana ULCS1 gene, encoding a WDR protein, results in cell wall modification impairment and plant infertility

Ubiquitin mediated protein degradation constitutes one of the most complex post translational gene regulation mechanisms in eukaryotes. This fine-tuned proteolytic machinery is based on a vast number of E3 ubiquitin ligase complexes that mark target proteins with ubiquitin. The specificity is accomplished by a number of adaptor proteins that contain functional binding domains, including the WD40 repeat motif (WDRs). To date, only few of these proteins have been identified in plants. An RNAi mediated silencing approach was used here to functionally characterize the Arabidopsis thaliana ULCS1 gene, which encodes for a small molecular weight WDR protein. AtULCS1 interacts with the E3Cullin Ring Ligase subunit DDB1a, regulating most likely the degradation of specific proteins involved in the manifestation of diverse developmental events. Silencing of AtULCS1 results in sterile plants with pleiotropic phenotypes. Detailed analysis revealed that infertility is the outcome of anther indehiscence, which in turn is due to the impairment of the plants to accomplish secondary wall modifications. Furthermore, IREGULAR XYLEM gene expression and lignification is diminished in anther endothecium and the stem vascular tissue of the silenced plants. These data underline the importance of AtULCS1 in plant development and reproduction.

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.

The architecture of the spliceosomal U4/U6.U5 tri-snRNP

U4/U6.U5 tri-snRNP is a 1.5-megadalton pre-assembled spliceosomal complex comprising U5 small nuclear RNA (snRNA), extensively base-paired U4/U6 snRNAs and more than 30 proteins, including the key components Prp8, Brr2 and Snu114. The tri-snRNP combines with a precursor messenger RNA substrate bound to U1 and U2 small nuclear ribonucleoprotein particles (snRNPs), and transforms into a catalytically active spliceosome after extensive compositional and conformational changes triggered by unwinding of the U4 and U6 (U4/U6) snRNAs. Here we use cryo-electron microscopy single-particle reconstruction of Saccharomyces cerevisiae tri-snRNP at 5.9 Å resolution to reveal the essentially complete organization of its RNA and protein components. The single-stranded region of U4 snRNA between its 3' stem-loop and the U4/U6 snRNA stem I is loaded into the Brr2 helicase active site ready for unwinding. Snu114 and the amino-terminal domain of Prp8 position U5 snRNA to insert its loop I, which aligns the exons for splicing, into the Prp8 active site cavity. The structure provides crucial insights into the activation process and the active site of the spliceosome.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly

BACKGROUND

Binding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure.

RESULTS

We investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical.

CONCLUSIONS

The reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly

Background

Binding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure.

Results

We investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical.

Conclusions

The reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.

Reviewers

This article was reviewed by Andrei L. Osterman, Narayanaswamy Srinivasan, Igor N. Berezovsky, and Gerrit Vriend (nominated by Martijn Huynen).

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0059-4) contains supplementary material, which is available to authorized users.

Non-manifesting AHI1 truncations indicate localized loss-of-function tolerance in a severe Mendelian disease gene

Determination of variant pathogenicity represents a major challenge in the era of high-throughput sequencing. Erroneous categorization may result if variants affect genes that are in fact dispensable. We demonstrate that this also applies to rare, apparently unambiguous truncating mutations of an established disease gene. By whole-exome sequencing (WES) in a consanguineous family with congenital non-syndromic deafness, we unexpectedly identified a homozygous nonsense variant, p.Arg1066*, in AHI1, a gene associated with Joubert syndrome (JBTS), a severe recessive ciliopathy. None of four homozygotes expressed any signs of JBTS, and one of them had normal hearing, which also ruled out p.Arg1066* as the cause of deafness. Homozygosity mapping and WES in the only other reported JBTS family with a homozygous C-terminal truncation (p.Trp1088Leufs*16) confirmed AHI1 as disease gene, but based on a more N-terminal missense mutation impairing WD40-repeat formation. Morpholinos against N-terminal zebrafish Ahi1, orthologous to where human mutations cluster, produced a ciliopathy, but targeting near human p.Arg1066 and p.Trp1088 did not. Most AHI1 mutations in JBTS patients result in truncated protein lacking WD40-repeats and the SH3 domain; disease was hitherto attributed to loss of these protein interaction modules. Our findings indicate that normal development does not require the C-terminal SH3 domain. This has far-reaching implications, considering that variants like p.Glu984* identified by preconception screening ('Kingsmore panel') do not necessarily indicate JBTS carriership. Genomes of individuals with consanguineous background are enriched for homozygous variants that may unmask dispensable regions of disease genes and unrecognized false positives in diagnostic large-scale sequencing and preconception carrier screening.

A general computational approach for repeat protein design

Repeat proteins have considerable potential for use as modular binding reagents or biomaterials in biomedical and nanotechnology applications. Here we describe a general computational method for building idealized repeats that integrates available family sequences and structural information with Rosetta de novo protein design calculations. Idealized designs from six different repeat families were generated and experimentally characterized; 80% of the proteins were expressed and soluble and more than 40% were folded and monomeric with high thermal stability. Crystal structures determined for members of three families are within 1Å root-mean-square deviation to the design models. The method provides a general approach for fast and reliable generation of stable modular repeat protein scaffolds.

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.

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.

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

The analysis of 36 available crystal structures of WD40 repeat proteins reveals widespread existence of a beta-bulge formed at the beginning of strand a and the end of strand b, termed as WD_b–a bulge: among a total of 259 WD40 blades, there are 243 such β-bulges. The R₁ positions in these WD_b–a bulges have fair distributions of Arg, His, Ile, Leu, Lys, Met, Phe, Trp, Tyr and Val residues. These residues protrude on the top face of the WD40 proteins and can serve as hotspots for protein-protein interactions. An analysis of 29 protein complexes formed by 17 WD proteins reveals that these R₁ residues, along with two other residues (R₁-2 and D-1), are indeed widely involved in protein-protein interactions. Interestingly, these WD_b–a bulges can be easily identified by the 4-amino acid sequences of (V, L, I), R₁, R₂, (V, L, I), along with some other significant amino acids. Thus, the hotspots of WD40 proteins on the top face can be readily predicted based on the primary sequences of the proteins. The literature-reported mutagenesis studies for Met30, MDV1, Tup11, COP1 and SPA1, which crystal structures are not available, can be readily understood based on the feature-based method. Applying the method, the twelve potential hotspots on the top face of Tup11 from S. japonicas have been identified. Our ITC measurements confirm seven of them, Tyr382, Arg284, Tyr426, Tyr508, Leu559, Lys575 and Ile601, are essential for recognizing Fep1. The ITC measurements further convinced that the feature-based method provides accurate prediction of hotspots on the top face.

Structure and function of WD40 domain proteins

The WD40 domain exhibits a β-propeller architecture, often comprising seven blades. The WD40 domain is one of the most abundant domains and also among the top interacting domains in eukaryotic genomes. In this review, we will discuss the identification, definition and architecture of the WD40 domains. WD40 domain proteins are involved in a large variety of cellular processes, in which WD40 domains function as a protein-protein or protein-DNA interaction platform. WD40 domain mediates molecular recognition events mainly through the smaller top surface, but also through the bottom surface and sides. So far, no WD40 domain has been found to display enzymatic activity. We will also discuss the different binding modes exhibited by the large versatile family of WD40 domain proteins. In the last part of this review, we will discuss how post-translational modifications are recognized by WD40 domain proteins.