Is Asp-His-Ser/Thr-Trp tetrad hydrogen-bond network important to WD40-repeat proteins: a statistical and theoretical study

GORI, encoding the WD40 domain protein, is required for pollen tube germination and elongation in rice

Successful delivery of sperm cells to the embryo sac in higher plants is mediated by pollen tube growth. The molecular mechanisms underlying pollen germination and tube growth in crop plants remain rather unclear, although these mechanisms are crucial to plant reproduction and seed formation. By screening pollen-specific gene mutants in rice (Oryza sativa), we identified a T-DNA insertional mutant of Germinating modulator of rice pollen (GORI) that showed a one-to-one segregation ratio for wild type (WT) to heterozygous. GORI encodes a seven-WD40-motif protein that is homologous to JINGUBANG/REN4 in Arabidopsis. GORI is specifically expressed in rice pollen, and its protein is localized in the nucleus, cytosol and plasma membrane. Furthermore, a homozygous mutant, gori-2, created through CRISPR-Cas9 clearly exhibited male sterility with disruption of pollen tube germination and elongation. The germinated pollen tube of gori-2 exhibited decreased actin filaments and altered pectin distribution. Transcriptome analysis revealed that 852 pollen-specific genes were downregulated in gori-2 compared with the WT, and Gene Ontology enrichment analysis indicated that these genes are strongly associated with cell wall modification and clathrin coat assembly. Based on the molecular features of GORI, phenotypical observation of the gori mutant and its interaction with endocytic proteins and Rac GTPase, we propose that GORI plays key roles in forming endo-/exocytosis complexes that could mediate pollen tube growth in rice.

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.

Structural modules of the stress-induced protein HflX: an outlook on its evolution and biological role

Multifunctional proteins often show modular structures. A functional domain and the structural modules within the domain show evolutionary conservation of their spatial arrangement since that gives the protein its functionality. However, the question remains as to how members of different domains of life (Archaea, Bacteria, Eukarya), polish and perfect these modules within conserved multidomain proteins, to tailor functional proteins according to their specific requirements. In the quest for plausible answers to this question, we studied the bacterial protein HflX. HflX is a universally conserved member of the Obg-GTPase superfamily but its functional role in Archaea and Eukarya is barely known. It is a multidomain protein and possesses, in addition to its conserved GTPase domain, an ATP-binding N-terminal domain. It is involved in heat stress response in Escherichia coli and our laboratory recently identified an ATP-dependent RNA helicase activity of E. coli HflX, which is likely instrumental in rescuing ribosomes during heat stress. Because perception and response to stress is expected to be different in different life forms, the question is whether this activity is preserved in higher organisms or not. Thus, we explored the evolution pattern of different structural modules of HflX, with particular emphasis on the ATP-binding domain, to understand plausible biological role of HflX in other forms of life. Our analyses indicate that, while the evolutionary pattern of the GTPase domain follows a conserved phylogeny, conservation of the ATP-binding domain shows a complicated pattern. The limited analysis described here hints towards possible evolutionary adaptations and modifications of the domain, something which needs to be investigated in more depth in homologs from other life forms. Deciphering how nature 'tweaks' such modules, both structurally and functionally, may help in understanding the evolution of such proteins, and, on a large-scale, of stress-related proteins in general as well.

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.

Towards the construction of an interactome for Human WD40 protein family

WD40 proteins are involved in a variety of protein-protein interactions as part of a multi-protein assembly modulating diverse and critical cellular process. It is known that several proteins of this family have been implicated in different disorders such as developmental abnormalities and cancer. However, molecular functions of many proteins in this family are yet unknown and it is of clinical interest. Therefore, it is of interest to define, construct, understand, analyze, evaluate, redefine and refine an interactome for WD40 protein family. We used data from literature mining using Cytoscape followed by linear regression analysis between Betweenness centrality and stress scores to define a model to filter the nodes in a representative WD40 interactome construction. We identified 10 ranked nodes in this analysis and subsequent microarray data selected three of them in insulin resistance that is further demonstrated in HepG2 cell culture models. We also observed the expression of GRWD1, RBBP5 and WDR5 genes during perturbation. Thus, we report hub nodes of WD40 interactome in insulin resistance. It should be noted that the pipeline using protein interaction network help find new proteins of clinical importance.

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.

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-function relationships of the G domain, a canonical switch motif

GTP-binding (G) proteins constitute a class of P-loop (phosphate-binding loop) proteins that work as molecular switches between the GDP-bound OFF and the GTP-bound ON state. The common principle is the 160-180-residue G domain with an α,β topology that is responsible for nucleotide-dependent conformational changes and drives many biological functions. Although the G domain uses a universally conserved switching mechanism, its structure, function, and GTPase reaction are modified for many different pathways and processes.

Functional analysis of archaeal MBF1 by complementation studies in yeast

Background

Multiprotein-bridging factor 1 (MBF1) is a transcriptional co-activator that bridges a sequence-specific activator (basic-leucine zipper (bZIP) like proteins (e.g. Gcn4 in yeast) or steroid/nuclear-hormone receptor family (e.g. FTZ-F1 in insect)) and the TATA-box binding protein (TBP) in Eukaryotes. MBF1 is absent in Bacteria, but is well- conserved in Eukaryotes and Archaea and harbors a C-terminal Cro-like Helix Turn Helix (HTH) domain, which is the only highly conserved, classical HTH domain that is vertically inherited in all Eukaryotes and Archaea. The main structural difference between archaeal MBF1 (aMBF1) and eukaryotic MBF1 is the presence of a Zn ribbon motif in aMBF1. In addition MBF1 interacting activators are absent in the archaeal domain. To study the function and therefore the evolutionary conservation of MBF1 and its single domains complementation studies in yeast (mbf1Δ) as well as domain swap experiments between aMBF1 and yMbf1 were performed.

Results

In contrast to previous reports for eukaryotic MBF1 (i.e. Arabidopsis thaliana, insect and human) the two archaeal MBF1 orthologs, TMBF1 from the hyperthermophile Thermoproteus tenax and MMBF1 from the mesophile Methanosarcina mazei were not functional for complementation of an Saccharomyces cerevisiae mutant lacking Mbf1 (mbf1Δ). Of twelve chimeric proteins representing different combinations of the N-terminal, core domain, and the C-terminal extension from yeast and aMBF1, only the chimeric MBF1 comprising the yeast N-terminal and core domain fused to the archaeal C-terminal part was able to restore full wild-type activity of MBF1.

However, as reported previously for Bombyx mori, the C-terminal part of yeast Mbf1 was shown to be not essential for function. In addition phylogenetic analyses revealed a common distribution of MBF1 in all Archaea with available genome sequence, except of two of the three Thaumarchaeota; Cenarchaeum symbiosum A and Nitrosopumilus maritimus SCM1.

Conclusions

The absence of MBF1-interacting activators in the archaeal domain, the presence of a Zn ribbon motif in the divergent N-terminal domain of aMBF1 and the complementation experiments using archaeal- yeast chimeric proteins presented here suggests that archaeal MBF1 is not able to functionally interact with the transcription machinery and/or Gcn4 of S. cerevisiae. Based on modeling and structural prediction it is tempting to speculate that aMBF1 might act as a single regulator or non-essential transcription factor, which directly interacts with DNA via the positive charged linker or the basal transcription machinery via its Zn ribbon motif and the HTH domain. However, also alternative functions in ribosome biosynthesis and/or functionality have been discussed and therefore further experiments are required to unravel the function of MBF1 in Archaea.

Reviewers

This article was reviewed by William Martin, Patrick Forterre, John van der Oost and Fabian Blombach (nominated by Eugene V Koonin (United States)). For the full reviews, please go to the Reviewer's Reports section.

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

We recently found that Asp-His-Ser/Thr-Trp hydrogen-bonded tetrads are widely and uniquely present in the WD40-repeat proteins. WDR5 protein is a seven WD40-repeat propeller with five such tetrads. To explore the effect of the tetrad on the structure and stability of WD40-repeat proteins, the wild-type WDR5 and its seven mutants involving the substitutions of tetrad residues have been isolated. The crystal structures of the wild-type WDR5 and its three WDR5 mutants have been determined by X-ray diffraction method. The mutations of the tetrad residues are found not to change the basic structural features. The denaturing profiles of the wild type and the seven mutants with the use of denaturant guanidine hydrochloride have been studied by circular dichroism spectroscopy to determine the folding free energies of these proteins. The folding free energies of the wild type and the S62A, S146A, S188A, D192E, W330F, W330Y, and D324E mutants are measured to be about -11.6, -2.7, -3.1, -2.9, -3.6, -7.1, -7.0, and -7.5 kcal/mol, respectively. These suggest that (1) the hydrogen bonds in these hydrogen bond networks are unusually strong; (2) each hydrogen-bonded tetrad provides over 12 kcal/mol stability to the protein; thus, the removal of any single tetrad would cause unfolding of the protein; (3) since there are five tetrads, the protein must be in a highly unstable state without the tetrads, which might be related to its biological functions.