GIGANTUS1 (GTS1), a member of Transducin/WD40 protein superfamily, controls seed germination, growth and biomass accumulation through ribosome-biogenesis protein interactions in Arabidopsis thaliana
- Emma W Gachomo†1, 2,
- Jose C Jimenez-Lopez†3, 4,
- Lyla Jno Baptiste1 and
- Simeon O Kotchoni1, 2Email author
© Gachomo et al.; licensee BioMed Central Ltd. 2014
Received: 30 November 2013
Accepted: 24 January 2014
Published: 27 January 2014
WD40 domains have been found in a plethora of eukaryotic proteins, acting as scaffolding molecules assisting proper activity of other proteins, and are involved in multi-cellular processes. They comprise several stretches of 44-60 amino acid residues often terminating with a WD di-peptide. They act as a site of protein-protein interactions or multi-interacting platforms, driving the assembly of protein complexes or as mediators of transient interplay among other proteins. In Arabidopsis, members of WD40 protein superfamily are known as key regulators of plant-specific events, biologically playing important roles in development and also during stress signaling.
Using reverse genetic and protein modeling approaches, we characterize GIGANTUS1 (GTS1), a new member of WD40 repeat protein in Arabidopsis thaliana and provide evidence of its role in controlling plant growth development. GTS1 is highly expressed during embryo development and negatively regulates seed germination, biomass yield and growth improvement in plants. Structural modeling analysis suggests that GTS1 folds into a β-propeller with seven pseudo symmetrically arranged blades around a central axis. Molecular docking analysis shows that GTS1 physically interacts with two ribosomal protein partners, a component of ribosome Nop16, and a ribosome-biogenesis factor L19e through β-propeller blade 4 to regulate cell growth development.
Our results indicate that GTS1 might function in plant developmental processes by regulating ribosomal structural features, activities and biogenesis in plant cells. Our results suggest that GIGANTUS1 might be a promising target to engineer transgenic plants with higher biomass and improved growth development for plant-based bioenergy production.
KeywordsArabidopsis Gigantus1 Gene expression Homology modeling Docking
In plants, growth development, cell patterning, yield, and biomass accumulation are controlled by functional genetic networks of several orders of magnitude that will require various complementary research approaches to uncover the basis of this functional genetic complexity. Using a combination of reverse genetic and computational protein modeling, we identified a protein, termed GIGANTUS1 (GTS1), a member of transducin/WD40 protein superfamily that regulates growth development in plants.
Transducin/WD40 repeat proteins are prominent features within proteins that mediate diverse protein-protein interactions, including those involved in scaffolding and the cooperative assembly and regulation of dynamic multi-subunit complexes [1, 2]. The common and defined feature of these proteins are short ~40 amino acid motifs, typically ending by Trp-Asp sequence, and usually composed of a 7-8 bladed beta-propeller fold. However, a diversity of proteins has been found with 4 to 16 repeated units . The low level of sequence conservation and functional diversity of WD40 domains make their identification difficult.
Repeated WD40 domains play central roles in biological processes such as cell division and cytokinesis, apoptosis, light signaling and vision, cell motility, flowering, floral development, meristem organization, protein trafficking, cytoskeleton dynamics, chemotaxis, nuclear export to RNA processing, chromatin modification, and transcriptional mechanism . They act as a site for protein-protein interaction, where the specificity of the proteins is determined by the sequences outside the repeats themselves. Their functional importance resides largely on the protein surfaces. They serve as multi-interacting platforms in cellular networks for the assembly of protein complexes or mediators of transient interplay among other proteins . Structural studies suggest that this property stems from their ability to interact with diverse proteins, peptides or nucleic acids using multiple surfaces, where the most common peptide interaction site of WD40 proteins is located on the top surface of the propeller close to the central channel .
Although WD40 proteins also are present in bacteria, e.g. Thermomonospora curvata and Cyanobacterium synechocystis, WD40 domains are among the ten most abundant domain types across eukaryotic proteomes, and interactome studies suggest that they are among the most promiscuous interactors. Despite several WD40-containing proteins acting as key regulators of plant-specific developmental events, WD40 domains have been given less research attention compared to other common domains, for example kinase domains. Up to now, there have been no comprehensive 3D structural analyses of WD40 protein revealing the interacting partners and highlighting the relevance of WD40 domains mediating different biological functions . Furthermore, in contrast to other members of the β-propeller family, and despite being crucial for and residing in enzymatic complexes, no WD40 protein has been reported to possess catalytic activity.
In order to fully understand the molecular- and structural-based regulatory mechanisms of GTS1 and the transducin/WD40 protein family, we carried out a comprehensive expression profiling, a mutational-based phenotypic characterization and a structural-based protein modeling of functional features of GIGANTUS1 in Arabidopsis. This is the first study in which the interacting partners of Arabidosis GTS1 protein have been identified. A comprehensive molecular and structural analysis based on protein homology modeling, and a docking interaction study to elucidate the functional mechanisms of GTS1 protein in regulating growth development in Arabidopsis thaliana were carried out. We described the phenotypic characterization of a gts1 knockout mutant during growth development and assessed the GTS1-3D molecular structure and its docking features to uncover the regulatory relationship of GTS1 with other proteins.
Plant material and GTS1expression profiling analysis
Sequences of oligonucleotide primers used in this study
For TDNA insertion
For complementation test
For complementation test
GTS1database search and phylogenetic analysis
The Arabidopsis thaliana GTS1 cDNA sequence (AT2G47790) was obtained from the Arabidopsis-TAIR website and used to perform a nucleotide BLAST (BLASTn) search of the Oryza sativa (rice) genomic sequence on the SALK Institute RiceGE2 web interface. The rice sequence identified as having the highest degree of homology to the Arabidopsis cDNA was downloaded and translated. The translated rice sequence was then aligned to the Arabidopsis protein sequence to validate the identification of the gene. The rice cDNA sequences were then used to perform a (BLASTn) search against the sequenced Zea mays genome on the MaizeSequence.org website. BLAST searches against the maize genome produced a list of BAC sequences that aligned to the query sequence. The BAC with the highest level of similarity was indicated on a genome map. The complete BAC sequence was downloaded and aligned to the Oryza sativa cDNA sequence using the NCBI BLAST (bl2seq) algorithm. Once the putative exons had been identified for a specific gene homolog, the exon start and end positions were manually corrected based on the canonical splice site donor/acceptor sequences and the overlapping sequence from one putative exon to the next. Each putative maize cDNA was then translated into a protein sequence. The GTS1 protein sequences from both Oryza sativa and Zea mays were aligned to the Arabidopsis thaliana sequence using the EMBOSS global pairwise alignment algorithm  to get the percent identity between the proteins. Next a BLASTp search of the Homo sapiens proteome in the NCBI database was performed using the Arabidopsis thaliana GTS1 protein sequence. From this search the Wdr89 sequence of Homo sapiens was identified as the most homologous protein. The Wdr89 sequence was used to retrieve (BLASTp) the Rattus norvegicus and the Mus musculis sequence homologs. Finally, a ClustalW2 alignment  was performed aligning all of the GTS1 protein sequences from Arabidopsis thaliana, Oryza sativa, and Zea mays with the Wdr89 protein sequences from Homo sapiens, Rattus norvegicus, and Mus musculis.
The conserved WD40 protein sequence from Arabidopsis thaliana was identified using ScanProsite and the Swiss-Prot/TrEMBLE databases. The identified WD40 sequence was used to perform a BLASTp search against plant proteomes on NCBI. From the list of potential WD40 containing proteins, the cDNA sequence of several GTS1-like proteins were downloaded from Nicotiana benthamiana, Nicotiana tabacum, Pisum sativum, Phaseolus vulgaris, Medicago truncatula, Gossypium hirsutum, Lycopersicon esculentum, Solanum chacoense, Solanum lycopersicum, and Physcomitrella patens. The collected plant WD40 containing cDNA sequences were combined with the Arabidopsis thaliana, Oryza sativa, and Zea mays cDNA sequences previously obtained and aligned using RevTrans version 1.4 . Using MrAIC [15–20] model test software the best AIC model to use in constructing phylogenetic relationships was determined to be the General Transition Rate (GTR) with Gamma.
Phylogenetic trees were generated in MrBayes using four Markov Chain Monte Carlo runs with three cold and one hot chain each for 1,000,000 generations sampling every 100 generations. The burnin was determined to be within the first 4000 generations for each phylogeny. The first 10000 generations sampled were removed and a 50% consensus majority tree was constructed from the remaining trees. Trees were then drawn using TreeView version 1.6.6 and rooted using the Physcomitrella patens sequence as an out group.
Sequences interaction database search
Arabidopsis thaliana GTS1, a Transducin/WD40 repeat protein (NCBI accession number AEC10888) was used as query to retrieve WD40 protein sequences from Uniprot (http://www.uniprot.org), and NCBI (http://www.ncbi.nlm.nih.gov) databases using BLAST tools (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
GTS1 protein interaction network was obtained using the STRING v9.0 (http://string-db.org) database. STRING outcome gave the two most possible interacting protein counterparts for GTS1 protein, a 60S ribosome structural protein L19e (NCBI accession number AEE75864), and Nop16 (NCBI accession number AAP21378), a protein involved in 60S subunit ribosomal biogenesis.
Functional domains analyses
GTS1 protein functional domains were studied by querying different structure-functional motifs and/or patterns databases such as Pfam v25.0 (http://pfam.sanger.ac.uk), Prosite (http://prosite.expasy.org/scanprosite), SMART v6.0 (http://smart.embl-heidelberg.de), Conserved Domain Database (CDD) v3.02, CDART (Conserved Domain Architecture Retrieval Tool) and CD-Search tools (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml), InterPRO v35.0 (http://www.ebi.ac.uk/interpro), ProDom release 2010.1 (http://prodom.prabi.fr/prodom/current/html/home.php), CATH v3.4 (http://www.cathdb.info), Superfamily v1.75 (http://supfam.cs.bris.ac.uk/SUPERFAMILY), PIRSF (http://pir.georgetown.edu), and functional searched by PANTHER (http://www.pantherdb.org). Similar analysis was performed for both interacting ribosomal protein counterparts, Nop16 and L19e.
Secondary structure prediction
Secondary structural elements of the GTS1 protein were initially assessed for substructure conserved motifs by threading the sequences through the Protein Data Bank (PDB) (http://www.pdb.org) library using threading algorithm the Segmer . These elements of the secondary structure were also confirmed by comparison with the results obtained with other additional 2-D structure servers: SSpro8 (Scratch Protein Predictor, http://scratch.proteomics.ics.uci.edu), NetSurfP ver. 1.1 (http://www.cbs.dtu.dk), and PSIPRED (http://bioinf.cs.ucl.ac.uk/psipred) fold servers. These secondary structure predictions were also performed for both interacting ribosomal protein counterparts, Nop16 and L19e.
Structural templates searching
Protein sequences of GTS1, Nop16 and L19e were searched for homology in the Protein Data Bank (PDB). Homologous templates suitable for these three proteins were selected by BLASTp from the BLAST server (http://ncbi.nlm.nih.gov). The BioInfoBank Metaserver (http://meta.bioinfo.pl), which employs fold recognition for homology search, was also used for the selection of templates. The results obtained by previous methods were compared with the results obtained by Swiss-model server for template identification (http://swissmodel.expasy.org).
Homology modeling was performed by the I-TASSER server . An initial structural model was generated and checked for recognition of errors in 3D structures using ProSA (http://prosa.services.came.sbg.ac.at/prosa.php), and for a first overall quality estimation of the model with QMEAN (http://swissmodel.expasy.org/qmean/cgi/index.cgi).
Energy minimization of the final structures was performed using GROMOS96 force field energy implemented in DeepView/Swiss-PDBViewer v3.7 (http://spdbv.vital-it.ch) in order to improve the van der Waals contacts and correct the stereochemistry of the model.
Each structure was assessed using the following softwares: QMEAN for quality, PROCHECK (http://www.ebi.ac.uk/thornton-srv/software/PROCHECK) for stereological corrections, ProSA, and ANOLEA (http://protein.bio.puc.cl/cardex/servers/anolea) for protein energy. The number of protein residues in the favored regions for each structure were calculated and visualized by The Ramachandran plot.
Ligand-binding domains and conservational analysis
Ligand-binding sites in the 3D protein structures were analyzed using Cofactor software (http://zhanglab.ccmb.med.umich.edu/COFACTOR), to identify functional homology. Gene Ontology (GO) terms (The Gene Ontology project) were used to identify functional analogs based on the 3D built models, indicating the possible functions and biological pathway in which the proteins might be involved (http://www.geneontology.org).
Conservational analyses of the proteins were made by generating evolutionary related conservation scores using ConSurf server (http://consurf.tau.ac.il). Structural function conservation and key residues in the query proteins were identified by ConSeq server (http://conseq.tau.ac.il).
Surface electrostatic potential analysis
The electrostatic Poisson-Boltzmann (PB) potentials for the surface amino acids of the structures were obtained using APBS software implemented in PyMol 0.99 (http://www.pymol.org) with AMBER99, and optimized with the Python software package PDB2PQR. Fine grid spaces of 0.35 Å were used to solve the linearized PB equation in sequential focusing multigrid calculations in a mesh of 130 points per dimension at 310.00 K. The dielectric constants were two for the proteins and 80.00 for water. The output mesh was processed in the scalar OpenDX format to render the isocontours and maps on the surfaces with PyMOL 0.99. Potential values are given in units of kT per unit charge (k Boltzmann’s constant; T temperature).
Molecuolar docking analysis
The analysis of interactions between GTS1 with each of its ribosomal counterparts (Nop16 and L19e,) was performed using CLUSpro server . During the workflow, backbone flexibility analysis was done using rigid-body ensemble docking with multiple structures derived from NMR, while the ZDOCK option for sampling at 6-degree rotational steps was used to obtain the decoys. The energies of docked conformations from protein-protein docking were evaluated by applying the Fast Fourier Transform (FFT) correlation approach.Docking scores were calculated by considering shape complementarity, desolvation, and electrostatics potential. The top docked conformations, along with their ZDOCK scores, were used as candidates of near-native structures, and for clustering of binding sites where the ligand was within 10 Å of its receptor.
After clustering, the ranked complexes were subjected to van der Waals minimization using CHARMM, and the protein-inhibitor structure with the best score was chosen as the most fitting model for GTS1-ribosomal interacting proteins.
Molecular and structural analysis of GTS1
Expression profiling and phylogenetic analysis of GTS1
Transducin/WD40 repeat-like superfamily protein
Evolutionarily conserved C-terminal region 2
Kow domain-containing transcription factor 1
VIRE2-interacting protein 1
Protein of unknown function DUF1446 (InterPro:IPR010839)
Protein of unknown function (DUF2361)
2-isopropylmalate synthase 1
Protein kinase superfamily protein
RNA-binding (RRM/RBD/RNP motifs) family protein
SWAP (Suppressor-of-White-APricot)/surp domain-containing protein/D111/G-patch domain-containing protein
Ribosomal protein L19e family protein
Pale cress protein (PAC)
Ribosomal protein S25 family protein
NAC domain containing protein 16
Argonaute family protein
KNOTTED-like homeobox of Arabidopsis thaliana 7
Surfeit locus protein 5 subunit 22 of mediator complex
GTS1 regulates seed germination and plant growth development
Searching for structural templates
Since we confirmed that GTS1 regulates seed germination and growth development (biomass yield, and flowering time) in plants (Figures 2 and 3), we next examined the protein structure scaffold and interacting partners of GTS1 in accomplishing its functions. In order to study the physical interactions of GTS1 with other proteins in regulating growth development in plants, we performed a Protein Data Bank (PDB) search for GTS1 protein with known tertiary structure in PDB. The search yielded the crystal structures/PDB accession numbers 2h9l, 3iz6, 3ow8, 2gnq, and 1tbg, showing the highest sequence identity (22, 17, 16, 24, and 15%, respectively). The suitability of selected model was checked by BioInfoBank Metaserver, which returned 3D Jury scores (J-score) of 208.4 (2h9l), 200.2 (3iz6), 198.5 (3ow8), 205.7 (2gnq), and 201.9 (1tbg) for GTS1, respectively. In order to confirm the best possible templates to use for building the GTS1 structure, Swiss-Model server was used, finding high scores (64, 61, 63, 60, and 69) and very low E-values (1.5E-37, 2.1E-30, 4.9E-29, 2E-10, and 1.8E-43) for the templates 2h9l, (3iz6), 3ow8, 2gnq, and 1tbg respectively.
The same workflow was followed to obtain the best crystal templates in order 3D-structural protein model for the GTS1 ribosomal interacting partners. The search in the Protein Data Bank (PDB) for the protein L19e and Nop16 yielded the crystal structures of 3iz5, 3jyw, and 3u5e for L19e, and 2aje, 1w0t, 2juh, 2ckx, and 2roh for Nop16, showing comparable values in identity and suitability by BioInforMank Metaserver and Swiss-Model server analysis.
Quality of threading models
To assess the quality of the protein models I-TASSER and Procheck analysis were performed. The I-TASSER analysis gave the accuracy parameters such as a C-score of -0.9, 0.60 ± 0.14 TM-score with 1848 decoys and 0.1467 of cluster density for GTS1, while Procheck analysis revealed that the main chain conformations of GTS1 protein model were located in the acceptable regions of the Ramachandran plot. A majority of residues (77.6%) were in the most favorable regions, whereas 14.5% of the residues were placed in the allowed regions, and 6.5% were in generously allowed regions. Only 1.4% of the residues were present in the disallowed regions, respectively. The plot of ×1 versus ×2 torsion angles for each residue showed that most of the rotamers in GTS1 model was localized in low energy regions. iii) The ProSa analysis gave Z-scores of -5.71 for GTS1. The scores were within the range usually found for native proteins of similar size, i.e., -7.31, -4.02, -6.63, -7.93, and -7.33 for the templates 2h9l, 3iz6, 3ow8, 2gnq, and 1tbg crystal structures, respectively. iv) QMEAN analysis of GTS1 model revealed Q-values of 0.67. A quality factor of 0.793, 0.315, 0.71, 0.786, and 0.849 was estimated for the crystal structures of the templates 2h9l, 3iz6, 3ow8, 2gnq, and 1tbg, respectively, indicating that the GTS1 model is within the range of accuracy of the templates crystallographic structures. v) Root mean square deviation (RMSD) between GTS1 model and the crystal templates Cα backbones of the closed templates were 2.408 Å for 2gnq and 3.192 Å for 2h9l.
All of the above parameters were also determined for L19e and Nop16 protein models resulting in comparable structural quality values of our modeled L19e and Nop16 proteins.
3D structure of Arabidopsis GTS1
Surface electrostatic potential analysis (Figure 5B) reveals several prominent positively charged residues (blue regions), predominantly in the walls of the cylinder (tunnel) and C-terminal arm. The environment of this protein is essentially negatively charged (red regions) (Figure 5B), as highlighted by the Poisson - Boltzmann electrostatic potential. By assigning a value of +1 to basic residues (R, K) and -1 to the acidic residues (D, E), net charge of protein was calculated to be -22 for GTS1. The central tunnel of the GTS1 structure exhibited a predominantly positive charge in the top view. Comparison between GTS1 and other WD40 repeat proteins, such as templates 2gnq and 2h9l did not exhibited large differences in the general topology as it was further confirmed by the RMSD value of 2.408 Å and 3.192 Å, respectively, whereas significant differences were found in particular regions of the proteins such as the N-terminal, and C-terminal regions.
Conservational and functional/ligand-binding site analysis
Figure 2B shows a general view of the amino acids in the tunnel-like structure holding up the interacting peptides. GTS1 exhibit the putative active (most representative ligand-binding) site located in the center of the structure (Figure 6B, top view), which normally is another interacting polypeptide, containing several conserved but also few variable amino acids. A detailed view, showing the spatial distribution of residues responsible for the conformation of the ligand binding domain surrounding the ligand and directly implicated in this interaction are Y43, V44, F45, S61, N87, S107, F134, V194, S267, R329, and the peptide chain substrate bound to the tunnel-like structure (Figure 6B, detailed view at the right side). Conservational analysis of WD40 repeat proteins with significant close identity to Arabidopsis GTS1 among other species returned a large number of highly variable (bold) and small number of conserved (italic) residues as written above. Conformational predictions indicate that the peptide (ligand) is projected to and partially located inside the tunnel-like structure of GTS1.
Consurf conservation analysis showed in the Figure 7C highlights an equally distributed number of conserved and variable residues along the surface of L19e protein. The area of interaction with GTS1 seems not to be as well conserved as its counterpart Nop16. Few variable residues (blue color) are located in the area of interaction with the N-terminal tail of GTS1 (red arrows in Figure 7C). In addition, among the surface directly implicated in the interaction with the ribonucleic acid (blue arrows in Figures. 7C, D), S3, K5, I6, R9, L10, N36 seems to be well conserved, since only one residue, I6, was found to be variable (italic), and the rest exhibit an average index of conservation or highly conserved like R9 (bold).
GTS1 interaction mechanism: molecular docking analysis with ribosomal counterparts
In order to get insights into the GTS1 regulatory/multi-interacting mode with other proteins, we analyzed the conformational interaction between GTS1 and two interacting partners involved in the structure and biogenesis of ribosomes in Arabidopsis. This analysis was carried out by molecular docking, using newly modeled structures of the two ribosomal proteins.
With regard to energy, the electrostatic potential analysis for the contact surface of both structures exhibited a highly compatible fingerprint distribution of opposite charges, i.e., positive (in the contact area of Nop16 protein), and negative (mainly for GTS1) (Figure 8C). There are not large areas with hydrophilic character in the contact surface between both proteins, and the formation of the complex may be mediated by a high number of direct and water-mediated H-bonds.
GIGANTUS1 is here described to be very important in regulating plant growth development (seed germination, faster growth, flowering time, and biomass accumulation) (Figures 1, 2 and 3). This is the first time that a mutation in GTS1 has been implicated in early germination, growth and development in Arabidopsis thaliana. The molecular mechanism by which GIGANTUS1 regulates plant growth development is still unknown. As a member of WD40 protein family, GTS1 is expected to play central roles in different biological processes including cell division and cytokinesis, flowering, floral development, cytoskeleton dynamics, nuclear export to RNA processing, transcriptional mechanism, and protein-protein interactions . We postulate that GIGANTUS1 might primarily function as a site for protein-protein interaction or mediator of transient interplay among other proteins to regulate different biological processes in plants. The development of protein complexes involves regulatory interactions that are mainly controlled by scaffolding proteins, such as WD40 repeat motifs. These motifs are important features of diverse protein-protein interactions , providing an unbending platform for interactions of proteins with other cellular components and controlling therefore several vital functions of the cell, such as signaling cascades, cellular transport and apoptosis [29–31].
The WD40 domains in GTS1 protein are shown to contain seven or multiples of seven repeats forming a highly stable β-propeller structure (Figure 5A). The 7-fold β-propeller is the most stable β-sheet geometry characterizing the resolved WD40 structures and also used to identify WD40 proteins . However members of this protein family have been also found to contain as high as sixteen repeats . Proteins with less than 7 repeats form an incomplete β-propeller structure and require additional WD-repeats from their neighbors to stabilize themselves, making dimers . There is no apparent folding order for each repeat and the order in which repeats fold might vary among different WD40 proteins, or even within the same protein . These proteins are known be involved in light signaling/photomorphogenesis and flowering , auxin response and cell division , in meristem maintenance , floral development , seed development and flowering , chromatin-based gene silencing and organogenesis , protein turn-over, microtubule dynamics, phospholipid binding and vesicle coating . This justifies the great deal of research interest in the WD40 protein superfamily across plant species. Our data (Figures 2 and 3) suggest that GTS1 belongs to the WD40 protein subfamily regulating auxin response and cell division , meristem maintenance , floral development , seed development and flowering . In this study the Arabidopsis GTS1 ligand-binding domain (main functional domain) lies mainly on the top surface residues, which integrates parts of β-propeller domain (Figure 6B). However, our data revealed that GTS1 WD40 propellers have three distinct surfaces available for interactions: the top region of the propeller, the bottom region (Figure 9), and the circumference [4, 5] (Figures 8 and 9), suggesting the multi-functional properties of GTS1 protein through protein-protein interactions. Indeed, protein-protein and protein-peptide interactions involved the entry site of the central channel of the β-propeller (Figure 6B), where the majority of interaction partners (including small molecules) bind . N- or C-terminal extensions of GTS1 run parallel to the tunnel-like structure, which form the complete 7 WD40 repeat domains, making them accessible for interaction with other partners (Figures 8 and 9). WD40 domains can thus act as large interaction platforms for multiple protein interactions. In comparison to other domains, the proteins containing WD40 motif are components of several interaction pairs [41, 42], and act as scaffolds for larger complex assemblies. This work represents the first time the 3D-structural and molecular features of the GTS1 protein in plants have been examined.
To better understand what area of the GTS1 plays an important role in complex formation for functional conservation across plant species, we identified two fundamentally conserved regions; the first, located on the top rim, constituted by the residues composing blades 4 to 7, including N- and C-terminal arms; and the second large conserved surface is located on the bottom of the propeller and is mainly composed of blades 1, 2 and 3 (Figure 6A). These regions represent potential protein-protein interaction sites .
In general, plant and particularly Arabidopsis-WD repeat proteins are strongly conserved. Most of these proteins are components of basic cellular machinery regulating plant-specific processes. An interesting question arises as to how these proteins evolved into their specific cellular roles. One of the key functional processes of WD proteins is the biogenesis of eukaryotic ribosomes, a highly regulated and dynamic process that begins in the nucleolus with transcription of rRNA precursor (pre-rRNA) and rapidly packaged into the 90S ribonucleoprotein particle containing ribosomal proteins, non-ribosomal proteins, and snoRNA-containing ribonucleoprotein particles (snoRNPs). The 90S pre-RNPs are converted into 43S and 66S ribosome assembly intermediates, which ultimately give rise to mature 40S and 60S ribosomal subunits . It is well known that ribosome biogenesis is driven by a large number of pre-ribosomal factors that associate with and/or dissociate from the pre-ribosomal particles along the maturation pathway. Although there has been much progress to identify ribosome assembly intermediates and their protein and RNA constituents , the information about the architecture of these pre-rRNPs is scarce. It is unclear which proteins are the nearest neighbors within the assembled ribosomes and to what extent neighboring molecules function together.
WD40 protein-protein interaction motifs represent excellent candidates to mediate interactions in the multiprotein subcomplex comprising a neighborhood in assembling ribosomes because of their protein-protein multi-interacting versatility. More than 70 trans-acting factors required for ribosome assembly have been identified , as well as 80 additional assembly factors present in pre-ribosomes . Therefore, such WD40-containing proteins may nucleate the assembly of pre-ribosomes by interacting sequentially or simultaneously with other assembly factors or ribosomal proteins. Among the assembly factors, 17 proteins were found to contain WD40 motifs . Many of the annotated ribosome biogenesis-WD40 repeat proteins were shown to directly interact with, or regulate the levels of other proteins  or to be components of multiprotein subcomplexes. Yeast WD40 protein Ytm1 is a constituent of 66S pre-rRNPs, whose depletion resulted in a deficiency of 60S ribosomal subunits . Its homologue, mammalian WDR12 functions in the maturation of the 60S ribosomal subunit. WDR12 forms a stable complex with a novel member the nucleolar proteins Pes1 and Bop1 (Pe- BoW complex), which are crucial for processing of the 32S precursor ribosomal RNA (rRNA) and cell proliferation . Interestingly, a potential homologous complex of Pes1–Bop1–WDR12 in yeast (Nop7p-Erb1p-Ytm1p) is involved in the control of ribosome biogenesis and S phase entry .
The yeast WD40 repeat protein Mak11 that modulates a p21-activated protein kinase function is an essential factor in nuclear maturation of 60S ribosomal subunits and its depletion led to a cell cycle delay in G1, indicating an early step nucleolar role of Mak11 in ribosome assembly. Another sub-complex, transiently associated with late, nuclear pre-60S precursors, is composed of four proteins and contains Ipi3 as a WD40 repeat member .
In this study, a new interacting counterpart, the Arabidopsis Nop16 protein was identified as a potential ribosome biogenesis factor in plants, which could be implicated in formation of the 60S ribosomal precursor. This process may be regulated by an interaction with GTS1 (Figure 8). This interaction was studied using a docking analysis that showed a stable interaction between GTS1 and Nop16, involving the N-terminal tail and the 4th blade of the first partner, and a cleft formed in the second ribosomal factor by the N-terminal α-helix and the neighboring secondary elements. Another ribosomal protein (L19e protein) was found to interact with GTS1. L19e protein is implicated in the structural stability of ribosome. The interacting area between GTS1 and L19e is very close to that of Nop16 and GTS1 interacting area (Figure 8). This suggests that the interacting mechanism of regulating the biogenesis nucleolar factor Nop16 and the structural ribosome factor L19e may be competitive. Therefore both steps in the 60S ribosomal subunit formation, structural maturation and stabilization may be separate in the time and/or different cellular compartments.
In support of our data are two other examples of WD40 repeat ribosome biogenesis factors, Rrb1 and Sqt1, which interact directly with ribosomal proteins for 60S ribosomal subunit assembly. Rrb1 interacts with the ribosomal protein Rpl3 in the nucleus and regulates its levels , and Sqt1 interacts with Rpl10 in the cytoplasm . Both proteins have a role in the association of the corresponding ribosomal protein with the nascent 60S ribosomal subunits and might regulate the levels of the corresponding ribosomal protein. Other WD40 repeat proteins have been implicated in the formation and stabilization of the small ribosomal subunit 40S. Yeast RACK1 regulates the translation initiation by recruiting PKC to the ribosome [55, 56]. Four RACK1 orthologs identified in Arabidopsis thaliana may have a similar activity . These interactions could provide a mechanism to regulate translation activities of ribosome populations programmed with specific mRNAs .
The present study provides substantial evidence for the role of GIGANTUS1 in controlling seed germination, faster growth and biomass accumulation in plants. The gene is mainly expressed in meristemic regions and is therefore important in cell division. It is postulated to regulate growth development through diverse protein-protein interactions, including those involved in scaffolding and dynamic multi-subunit complexes such as the ribosomal protein biogenesis, stability and activity. Given its rich interaction surfaces, GTS1 functions probably as an adaptor in many different protein complexes or protein-DNA complexes in very diverse cellular processes that needs further research investigation. Our modeling data suggests that GTS1 mediates molecular recognition events mainly through the smaller top surface of domain, which comprises three residues forming a transient complex with other peptides. It would be interesting to further investigate GIGANTUS1 knockout genes in agronomically important crops with the aim of improving crop yield and biomass accumulation for sustainable plant-based biofuel production.
This work was supported by NSF-REU DBI # 1263163 grant (PI, Benedetto Piccoli) and Rutgers-University start-up funds to SOK.
- Gibson TJ: Cell regulation: determined to signal discrete cooperation. Trends Biochem Sci. 2009, 34: 471-482. 10.1016/j.tibs.2009.06.007.View ArticlePubMedGoogle Scholar
- Hurtley S: Spatial cell biology. Location, location, location. Introduction. Science. 2009, 326: 1205-10.1126/science.326.5957.1205.View ArticlePubMedGoogle Scholar
- Lambright DG, Sondek J, Bohm A, Skiba NP, Hamm HE, Sigler PB: The 2.0 A crystal structure of a heterotrimeric G protein. Nature. 1996, 379: 311-319. 10.1038/379311a0.View ArticlePubMedGoogle Scholar
- Stirnimann CU, Petsalaki E, Russell RB, Muller CW: WD40 proteins propel cellular networks. Trends Biochem Sci. 2010, 35: 565-574. 10.1016/j.tibs.2010.04.003.View ArticlePubMedGoogle Scholar
- Russell RB, Sasieni PD, Sternberg MJ: Supersites within superfolds. Binding site similarity in the absence of homology. J Mol Biol. 1998, 282: 903-918. 10.1006/jmbi.1998.2043.View ArticlePubMedGoogle Scholar
- Wilson DK, Cerna D, Chew E: The 1.1-angstrom structure of the spindle checkpoint protein Bub3p reveals functional regions. J Biol Chem. 2005, 280: 13944-13951. 10.1074/jbc.M412919200.View ArticlePubMedGoogle Scholar
- Janda L, Tichý P, Spízek J, Petrícek M: A deduced Thermomonospora curvata protein containing serine/threonine protein kinase and WD-repeat domains. J Bacteriol. 1996, 178: 1487-1489.PubMed CentralPubMedGoogle Scholar
- Grigorieva G, Shestakov S: Transformation in the cyanobacterium Synechocystis sp. 6803. FEMS Microbiol Lett. 1982, 13: 367-370. 10.1111/j.1574-6968.1982.tb08289.x.View ArticleGoogle Scholar
- Mishra AK, Puranik S, Prasad M: Structure and regulatory networks of WD40 protein in plants. J Plant Biochem Biotechnol. 2012, 21: S32-S39. 10.1007/s13562-012-0134-1.View ArticleGoogle Scholar
- Gachomo EW, Jimenez-Lopez JC, Smith SR, Cooksey AB, Oghoghomeh OM, Johnson N, Baba-Moussa L, Kotchoni SO: The cell morphogenesis ANGUSTIFOLIA (AN) gene, a plant homolog of CtBP/BARS, is involved in abiotic and biotic stress response in higher plants. BMC Plant Biol. 2013, 13: 1-11. 10.1186/1471-2229-13-1.View ArticleGoogle Scholar
- Kotchoni SO, Larrimore K, Mukherjee M, Kempinski C, Barth C: Alterations in the endogenous ascorbic acid content affect flowering time in Arabidopsis thaliana. Plant Physiol. 2009, 149: 803-815.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: The European molecular biology open software suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
- Wernersson R, Pedersen AG: RevTrans: multiple alignment of coding DNA from aligned amino acid sequences. Nucleic Acids Res. 2003, 31: 3537-3539. 10.1093/nar/gkg609.PubMed CentralView ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: Testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Burnham KP, Anderson DR: Model selection and multimodel inference, apractical information-theoretic approach. 2nd edition. New York:Springer; 2002.Google Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Antonelli A, Nylander JAA, Persson C, Sanmartín I: Tracing the impact of the Andean uplift on neotropical plant evolution. Proc Natl Acad Sci USA. 2009, 106: 9749-9754. 10.1073/pnas.0811421106.PubMed CentralView ArticlePubMedGoogle Scholar
- Nylander JA, Ronquist F, Huelsenbeck JP, Nieves-ALdrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.View ArticlePubMedGoogle Scholar
- Wu S, Zhang Y: Recognizing protein substructure similarity using segmental threading. Structure. 2010, 18: 858-867. 10.1016/j.str.2010.04.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9: 40-10.1186/1471-2105-9-40.View ArticleGoogle Scholar
- Kozakov D, Hall DR, Beglov D, Brenke R, Comeau SR, Shen Y, Li K, Zheng J, Vakili P, Paschalidis IC, Vajda S: Achieving reliability and high accuracy in automated protein docking: ClusPro, PIPER, SDU, and stability analysis in CAPRI rounds 13-19. Proteins. 2010, 78: 3124-3130. 10.1002/prot.22835.PubMed CentralView ArticlePubMedGoogle Scholar
- Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR: Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136: 2621-2632. 10.1104/pp.104.046367.PubMed CentralView ArticlePubMedGoogle Scholar
- Baulcombe DC, Saunders GS, Bevan MW, Mayo MA, Harrison BD: Expression of biologically active viral satellite RNA from the nuclear genome of transformed plants. Nature. 1986, 321: 446-449. 10.1038/321446a0.View ArticleGoogle Scholar
- Kotchoni SO, Kuhns C, Ditzer A, Kirch H-H, Bartels D: Over-expression of different aldehyde dehydrogenase genes in Arabidopsis thaliana confers tolerance to abiotic stress and protects plants against lipid peroxidation and oxidative stress. Plant Cell Environ. 2006, 29: 1033-1048. 10.1111/j.1365-3040.2005.01458.x.View ArticlePubMedGoogle Scholar
- Clough SJ, Bent AF: Floral dip: A simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998, 16: 735-743. 10.1046/j.1365-313x.1998.00343.x.View ArticlePubMedGoogle Scholar
- Ouyang H, Huang X, Lu Z, Yao J: Genomic survey, expression profile and co-expression network analysis of OsWD40 family in rice. BMC Genomics. 2012, 13: 100-10.1186/1471-2164-13-100.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu C, Min J: Structure and function of WD40 domain proteins. Protein Cell. 2011, 2: 202-214. 10.1007/s13238-011-1018-1.View ArticlePubMedGoogle Scholar
- Neer EJ, Schmidt CJ, Nambudripad R, Smith TF: The ancient regulatory-protein family of WD-repeat proteins. Nature. 1994, 371: 297-300. 10.1038/371297a0.View ArticlePubMedGoogle Scholar
- Hisbergues M, Gaitatzes CG, Joset F, Bedu S, Smith TF: A noncanonical WD-repeat protein from the cyanobacterium Synechocystis PCC6803: structural and functional study. Protein Sci. 2001, 10: 293-300. 10.1110/ps.22701.PubMed CentralView ArticlePubMedGoogle Scholar
- Valeyev NV, Downing AK, Sondek J, Deane C: Electrostatic and functional analysis of the seven-bladed WD beta-propellers. Evol Bioinform Online. 2008, 4: 203-216.PubMed CentralPubMedGoogle Scholar
- Saeki M, Irie Y, Ni L, Yoshida M, Itsuki Y, Kamisaki Y: Monad, a WD40 repeat protein, promotes apoptosis induced by TNF-alpha. Biochem Biophys Res Commun. 2006, 342: 568-572. 10.1016/j.bbrc.2006.02.009.View ArticlePubMedGoogle Scholar
- Thornton C, Tang KC, Phamluong K, Luong K, Vagts A, Nikanjam D, Yaka R, Ron D: Spatial and temporal regulation of RACK1 function and N-methyl-D-aspartate receptor activity through WD40 motif-mediated dimerization. J Biol Chem. 2004, 279: 31357-31364. 10.1074/jbc.M402316200.View ArticlePubMedGoogle Scholar
- Garcia-Higuera I, Gaitatzes C, Smith TF, Neer EJ: Folding a WD repeat propeller. Role of highly conserved aspartic acid residuesin the G protein beta subunit and Sec13. J Biol Chem. 1998, 273: 9041-9049. 10.1074/jbc.273.15.9041.View ArticlePubMedGoogle Scholar
- Lau OS, Deng XW: The photomorphogenic repressors COP1 and DET1: 20 years later. Trends Plant Sci. 2012, 17: 584-593. 10.1016/j.tplants.2012.05.004.View ArticlePubMedGoogle Scholar
- Varaud E, Brioudes F, Szécsi J, Leroux J, Brown S, Perrot-Rechenmann C, Bendahmane M: AUXIN RESPONSE FACTOR8 regulates Arabidopsis petal growth by interacting with the bHLH transcription factor BIGPETALp. Plant Cell. 2011, 23: 973-983. 10.1105/tpc.110.081653.PubMed CentralView ArticlePubMedGoogle Scholar
- Meng X, Muszynski MG, Danilevskaya ON: The FT-like ZCN8 gene functions as a floral activator and is involved in photoperiod sensitivity in maize. Plant Cell. 2011, 23: 942-960. 10.1105/tpc.110.081406.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, He Z, Lu G, Lee SC, Alonso J, Ecker JR, Luan S: A WD40 domain cyclophilin interacts with histone H3 and functions in gene repression and organogenesis in Arabidopsis. Plant Cell. 2007, 19: 2403-2416. 10.1105/tpc.107.053579.PubMed CentralView ArticlePubMedGoogle Scholar
- Jain M, Nijhawan A, Arora R, Agarwal P, Ray S, Sharma P, Kapoor S, Tyagi AK, Khurana JP: F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol. 2007, 143: 1467-1483. 10.1104/pp.106.091900.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu H, Braun P, Yildirim MA, Lemmens I, Venkatesan K, Sahalie J, Hirozane-Kishikawa T, Gebreab F, Li N, Simonis N, Hao T, Rual JF, Dricot A, Vazquez A, Murray RR, Simon C, Tardivo L, Tam S, Svrzikapa N, Fan C, de Smet AS, Motyl A, Hudson ME, Park J, Xin X, Cusick ME, Moore T, Boone C, Snyder M, Roth FP, Barabási AL, Tavernier J, Hill DE, Vidal M: High-quality binary protein interaction map of the yeast interactome network. Science. 2008, 322: 104-110. 10.1126/science.1158684.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins SR, Kemmeren P, Zhao XC, Greenblatt JF, Spencer F, Holstege FC, Weissman JS, Krogan NJ: Toward a comprehensive atlas of the physical interactome of Saccharomyces cerevisiae. Mol Cell Proteomics. 2007, 6: 439-450.View ArticlePubMedGoogle Scholar
- Miele AE, Watson PJ, Evans PR, Traub LM, Owen DJ: Two distinct interaction motifs in amphiphysin bind two independent sites on the clathrin terminal domain beta-propeller. Nat Struct Mol Biol. 2004, 11: 242-248. 10.1038/nsmb736.View ArticlePubMedGoogle Scholar
- Miles TD, Jakovljevic J, Horsey EW, Harnpicharnchai P, Tang L, Woolford JL: Ytm1, Nop7, and Erb1 form a complex necessary for maturation of yeast 66S preribosomes. Mol Cell Biol. 2005, 25: 10419-10432. 10.1128/MCB.25.23.10419-10432.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Fatica A, Tollervey D: Making ribosomes. Curr Opin Cell Biol. 2002, 14: 313-318. 10.1016/S0955-0674(02)00336-8.View ArticlePubMedGoogle Scholar
- Kressler D, Linder P, de La Cruz J: Protein trans-acting factors involved in ribosome biogenesis in Saccharomyces cerevisiae. Mol Cell Biol. 1999, 19: 7897-7912.PubMed CentralView ArticlePubMedGoogle Scholar
- Grandi P, Rybin V, Bassler J, Petfalski E, Strauss D, Marzioch M, Schäfer T, Kuster B, Tschochner H, Tollervey D, Gavin AC, Hurt E: 90S pre-ribosomes include the 35S pre-rRNA, the U3 snoRNP, and 40S subunit processing factors but predominantly lack 60S synthesis factors. Mol Cell. 2002, 10: 105-115. 10.1016/S1097-2765(02)00579-8.View ArticlePubMedGoogle Scholar
- Fromont-Racine M, Senger B, Saveanu C, Fasiolo F: Ribosome assembly in eukaryotes. Gene. 2003, 313: 17-42.View ArticlePubMedGoogle Scholar
- Saveanu C, Rousselle JC, Lenormand P, Namane A, Jacquier A, Fromont-Racine M: The p21-activated protein kinase inhibitor Skb15 and its budding yeast homologue are 60S ribosome assembly factors. Mol Cell Biol. 2007, 27: 2897-2909. 10.1128/MCB.00064-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Harnpicharnchai P, Jakovljevic J, Horsey E, Miles T, Roman J, Rout M, Meagher D, Imai B, Guo Y, Brame CJ, Shabanowitz J, Hunt DF, Woolford JL: Composition and functional characterization of yeast 66S ribosome assembly intermediates. Mol Cell. 2001, 8: 505-515. 10.1016/S1097-2765(01)00344-6.View ArticlePubMedGoogle Scholar
- Holzel M, Rohrmoser M, Schlee M, Grimm T, Harasim T, Malamoussi A, Gruber-Eber A, Kremmer E, Hiddemann W, Bornkamm GW, Eick D: Mammalian WDR12 is a novel member of the Pes1-Bop1 complex and is required for ribosome biogenesis and cell proliferation. J Cell Biol. 2005, 170: 367-378. 10.1083/jcb.200501141.PubMed CentralView ArticlePubMedGoogle Scholar
- Galani K, Nissan TA, Petfalski E, Tollervey D, Hurt E: Rea1, a dynein-related nuclear AAA-ATPase, is involved in late rRNA processing and nuclear export of 60 S subunits. J Biol Chem. 2004, 279: 55411-55418. 10.1074/jbc.M406876200.View ArticlePubMedGoogle Scholar
- Iouk TL, Aitchison JD, Maguire S, Wozniak RW: Rrb1p, a yeast nuclear WD-repeat protein involved in the regulation of ribosome biosynthesis. Mol Cell Biol. 2001, 21: 1260-1271. 10.1128/MCB.21.4.1260-1271.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Hofer A, Bussiere C, Johnson AW: Mutational analysis of the ribosomal protein Rpl10 from yeast. J Biol Chem. 2007, 282: 32630-32639. 10.1074/jbc.M705057200.View ArticlePubMedGoogle Scholar
- Ceci M, Gaviraghi C, Gorrini C, Sala LA, Offenhauser N, Marchisio PC, Biffo S: Release of eIF6 (p27BBP) from the 60S subunit allows 80S ribosome assembly. Nature. 2003, 426: 579-584. 10.1038/nature02160.View ArticlePubMedGoogle Scholar
- Nilsson J, Sengupta J, Frank J, Nissen P: Regulation of eukaryotic translation by the RACK1 protein: a platform for signalling molecules on the ribosome. EMBO Rep. 2004, 5: 1137-1141. 10.1038/sj.embor.7400291.PubMed CentralView ArticlePubMedGoogle Scholar
- Ullah H, Scappini EL, Moon AF, Williams LV, Armstrong DL, Pedersen LC: Structure of a signal transduction regulator, RACK1, from Arabidopsis thaliana. Protein Sci. 2008, 17: 1771-1780. 10.1110/ps.035121.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Baum S, Bittins M, Frey S, Seedorf M: Asc1p, a WD40-domain containing adaptor protein, is required for the interaction of the RNA-binding protein Scp160p with polysomes. Biochem J. 2004, 380: 823-830. 10.1042/BJ20031962.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.