Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS.     

Phylogeny, sequence conservation, and functional complementation of the SBDS protein family

The Shwachman-Bodian-Diamond syndrome (SBDS) protein family occurs widely in nature, although its function has not been determined. Comprehensive database searches revealed SBDS homologues from 159 species, including examples from all sequenced archaeal and eukaryotic genomes and all eukaryotic kingdoms. Sequence alignment with ClustalX and MUSCLE algorithms led to the identification of conserved residues that occurred predominantly in the amino-terminal FYSH domain where they appeared to contribute to protein folding or stability. Only SBDS residue Gly91 was invariant in all species. Four distantly related protists were found to have two divergent SBDS genes in their genomes. In each case, phylogenetic analyses and the identification of shared sequence features suggested that one gene was derived from lateral gene transfer. We also identified a shared C-terminal zinc finger domain fusion in flowering plants and chromalveolates that may shed light on the function of the protein family and the evolutionary histories of these kingdoms. To assess the extent of SBDS functional conservation, we carried out complementation studies of SBDS homologues and interspecies chimeras in Saccharomyces cerevisiae. We determined that the FYSH domain was widely interchangeable among eukaryotes, while domain 2 imparted species specificity to protein function. Domain 3 was largely dispensable for function in our yeast complementation assay. Overall, the phylogeny of SBDS was shared with a group of proteins that were markedly enriched for RNA metabolism and/or ribosome-associated functions. These findings link Shwachman-Diamond syndrome to other bone marrow failure syndromes with defects in nucleolus-associated processes, including Diamond-Blackfan anemia, cartilage-hair hypoplasia, and dyskeratosis congenita.     

Structural and mutational analysis of the SBDS protein family. Insight into the leukemia-associated Shwachman-Diamond Syndrome

Shwachman-Diamond Syndrome (SDS) is an autosomal recessive disorder characterized by bone marrow failure with significant predisposition to the development of poor prognosis myelodysplasia and leukemia, exocrine pancreatic failure and metaphyseal chondrodysplasia. Although the SBDS gene mutated in this disorder is highly conserved in Archaea and all eukaryotes, the function is unknown. To interpret the molecular consequences of SDS-associated mutations, we have solved the crystal structure of the Archaeoglobus fulgidus SBDS protein orthologue at a resolution of 1.9 angstroms, revealing a three domain architecture. The N-terminal (FYSH) domain is the most frequent target for disease mutations and contains a novel mixed alpha/beta-fold identical to the single domain yeast protein Yhr087wp that is implicated in RNA metabolism. The central domain consists of a three-helical bundle, whereas the C-terminal domain has a ferredoxin-like fold. By genetic complementation analysis of the essential Saccharomyces cerevisiae SBDS orthologue YLR022C, we demonstrate an essential role in vivo for the FYSH domain and the central three-helical bundle. We further show that the common SDS-related K62X truncation is non-functional. Most SDS-related missense mutations that alter surface epitopes do not impair YLR022C function, but mutations affecting residues buried in the hydrophobic core of the FYSH domain severely impair or abrogate complementation. These data are consistent with absence of homozygosity for the common K62X truncation mutation in individuals with SDS, indicating that the SDS disease phenotype is a consequence of expression of hypomorphic SBDS alleles and that complete loss of SBDS function is likely to be lethal.     

The crystal structure of Haloferax volcanii proliferating cell nuclear antigen reveals unique surface charge characteristics due to halophilic adaptation

Background

Defects in the human Shwachman-Bodian-Diamond syndrome (SBDS) protein-coding gene lead to the autosomal recessive disorder characterised by bone marrow dysfunction, exocrine pancreatic insufficiency and skeletal abnormalities. This protein is highly conserved in eukaryotes and archaea but is not found in bacteria. Although genomic and biophysical studies have suggested involvement of this protein in RNA metabolism and in ribosome biogenesis, its interacting partners remain largely unknown.

Results

We determined the crystal structure of the SBDS orthologue from Methanothermobacter thermautotrophicus (mthSBDS). This structure shows that SBDS proteins are highly flexible, with the N-terminal FYSH domain and the C-terminal ferredoxin-like domain capable of undergoing substantial rotational adjustments with respect to the central domain. Affinity chromatography identified several proteins from the large ribosomal subunit as possible interacting partners of mthSBDS. Moreover, SELEX (Systematic Evolution of Ligands by EXponential enrichment) experiments, combined with electrophoretic mobility shift assays (EMSA) suggest that mthSBDS does not interact with RNA molecules in a sequence specific manner.

Conclusion

It is suggested that functional interactions of SBDS proteins with their partners could be facilitated by rotational adjustments of the N-terminal and the C-terminal domains with respect to the central domain. Examination of the SBDS protein structure and domain movements together with its possible interaction with large ribosomal subunit proteins suggest that these proteins could participate in ribosome function.     

Crystal structure of KD93, a novel protein expressed in human hematopoietic stem/progenitor cells

The crystal structure of a novel hypothetical protein, KD93, expressed in human hematopoietic stem/progenitor cells, was determined at 1.9A resolution using the multiple-wavelength anomalous dispersion (MAD) method. The protein KD93, which is encoded by the open reading frame HSPC031, is a NIP7 homologue and belongs to the UPF0113 family. The structural and functional information for the group of homologues has not yet been determined. Crystallographic analysis revealed that the overall fold of KD93 consists of two interlinked alpha/beta domains. Structure-based homology analysis with DALI revealed that the C domain of KD93 matches the PUA domain of some RNA modification enzymes, especially that of archaeosine tRNA-ribosyltransferase (ArcTGT), which suggests that its possible molecular function is related to RNA binding. The difference between the RNA binding regions of KD93 and ArcTGT in amino acid constitution and surface electrostatic potential indicate that they may have different RNA binding modes. The N domain of KD93 is a unique structure with no obvious similarity to other proteins with known three-dimensional structures. The high-resolution structure of KD93 provides a first view of a member of the family of hypothetical proteins. And the structure provides a framework to deduce and assay the molecular function of other proteins of the UPF0113 family.     

Conserved structural elements in glutathione transferase homologues encoded in the genome of Escherichia coli

Multiple sequence alignments of the eight glutathione (GSH) transferase homologues encoded in the genome of Escherichia coli were used to define a consensus sequence for the proteins. The consensus sequence was analyzed in the context of the three-dimensional structure of the gst gene product (EGST) obtained from two different crystal forms of the enzyme. The enzyme consists of two domains. The N-terminal region (domain I) has a thioredoxin-like alpha/beta-fold, while the C-terminal domain (domain II) is all alpha-helical. The majority of the consensus residues (12/17) reside in the N-terminal domain. Fifteen of the 17 residues are involved in hydrophobic core interactions, turns, or electrostatic interactions between the two domains. The results suggest that all of the homologues retain a well-defined group of structural elements both in and between the N-terminal alpha/beta domain and the C-terminal domain. The conservation of two key residues for the recognition motif for the gamma-glutamyl-portion of GSH indicates that the homologues may interact with GSH or GSH analogues such as glutathionylspermidine or alpha-amino acids. The genome context of two of the homologues forms the basis for a hypothesis that the b2989 and yibF gene products are involved in glutathionylspermidine and selenium biochemistry, respectively.     

Voltage-gated H+ channels associated with human phagocyte superoxide-generating NADPH oxidases: sequence comparisons,structural predictions,and phylogenetic analyses

The N-terminal domain of the human phagocyte flavocytochrome b558 NADPH oxidase, gp91phox, is believed to be a heme-containing voltage-gated H+ channel. The authors have conducted structural, sequence and phylogenetic analyses of the putative transmembrane channel/heme-binding domains of all homologous proteins in the NCBI GenBank database as of May 2001, as well as of the full-length proteins. Fifty-six homologues were identified, including 26 from animals, 19 from plants, seven from yeast, one from a slime mould and three from bacteria. Six well-defined sub-families were revealed by phylogenetic tree construction, two consisting of animal proteins, two of plant proteins, and one each of yeast and bacterial homologues, with the slime mould protein clustering loosely with one of the animal clusters. Signature sequences for the entire family as well as for the sub-families were determined. Most proteins have six putative TMSs, four of which may comprise the heme-binding H+ channel. The hydrophobic and amphipathic characteristics of each of the putative alpha-helical transmembrane segments were defined, and conserved residues that may be involved in heme binding, channel formation, and/or conformational changes were identified. The analyses lead to the suggestion that the oxidase domain became associated with the channel/heme-binding domain to form a single polypeptide chain early in evolutionary history, before eukaryotes diverged from prokaryotes, and that genetic transmission to present day organisms occurred primarily by vertical descent.     

Voltage-gated H + channels associated with human phagocyte superoxide-generating NADPH oxidases: Sequence comparisons,structural predictions,and phylogenetic analyses

The N-terminal domain of the human phagocyte flavocytochrome b 558 NADPH oxidase, gp91 phox, is believed to be a heme-containing voltage-gated H + channel. The authors have conducted structural, sequence and phylogenetic analyses of the putative transmembrane channel/heme-binding domains of all homologous proteins in the NCBI GenBank database as of May 2001, as well as of the full-length proteins. Fifty-six homologues were identified, including 26 from animals, 19 from plants, seven from yeast, one from a slime mould and three from bacteria. Six well-defined sub-families were revealed by phylogenetic tree construction, two consisting of animal proteins, two of plant proteins, and one each of yeast and bacterial homologues, with the slime mould protein clustering loosely with one of the animal clusters. Signature sequences for the entire family as well as for the sub-families were determined. Most proteins have six putative TMSs, four of which may comprise the heme-binding H + channel. The hydrophobic and amphipathic characteristics of each of the putative &#102 -helical transmembrane segments were defined, and conserved residues that may be involved in heme binding, channel formation, and/or conformational changes were identified. The analyses lead to the suggestion that the oxidase domain became associated with the channel/heme-binding domain to form a single polypeptide chain early in evolutionary history, before eukaryotes diverged from prokaryotes, and that genetic transmission to present day organisms occurred primarily by vertical descent.     

Identification and characterization of a novel endoplasmic reticulum (ER) DnaJ homologue, which stimulates ATPase activity of BiP in vitro and is induced by ER stress

The activity of Hsp70 proteins is regulated by accessory proteins, which include members of the DnaJ-like protein family. Characterized by the presence of a highly conserved 70-amino acid J domain, DnaJ homologues activate the ATPase activity of Hsp70 proteins and stabilize their interaction with unfolded substrates. DnaJ homologues have been identified in most organelles where they are involved in nearly all aspects of protein synthesis and folding. Within the endoplasmic reticulum (ER), DnaJ homologues have also been shown to assist in the translocation, secretion, retro-translocation, and ER-associated degradation (ERAD) of secretory pathway proteins. By using bioinformatic methods, we identified a novel mammalian DnaJ homologue, ERdj4. It is the first ER-localized type II DnaJ homologue to be reported. The signal sequence of ERdj4 remains uncleaved and serves as a membrane anchor, orienting its J domain into the ER lumen. ERdj4 co-localized with GRP94 in the ER and associated with BiP in vivo when they were co-expressed in COS-1 cells. In vitro experiments demonstrated that the J domain of ERdj4 stimulated the ATPase activity of BiP in a concentration-dependent manner. However, mutation of the hallmark tripeptide HPD (His --> Gln) in the J domain totally abolished this activation. ERdj4 mRNA expression was detected in all human tissues examined but showed the highest level of the expression in the liver, kidney, and placenta. We found that ERdj4 was highly induced at both the mRNA and protein level in response to ER stress, indicating that this protein might be involved in either protein folding or ER-associated degradation.     

The SF3b155 N-terminal domain is a scaffold important for splicing

Recruitment of U2 snRNP to the branch point sequence of introns is a necessary step in pre-mRNA splicing. In the current model, U2AF65, bound at the polypyrimidine tract of the intron, recruits U2 snRNP to the branch point sequence by interacting with the U2 snRNP protein SF3b155. We demonstrate that the N-terminal domain of SF3b155 contains multiple U2AF65 binding sites that are distinct from the binding site for the U2 snRNP protein p14, mapped to amino acids 396-424 of SF3b155. The N-terminal domain of SF3b155 appears to adopt a primarily unfolded structure but is functional to inhibit splicing in vitro. RNA binding studies show that the N-terminal domain of SF3b155 binds RNA nonspecifically and that the sites for U2AF65 binding and RNA binding are overlapping (or the same) within SF3b155. We propose that the N-terminal domain of SF3b155 adopts a primarily unfolded structure that functions as a scaffold to facilitate SF3b155's multiple protein-protein and protein-RNA interactions. The multiple U2AF65 binding sites on SF3b155 suggest a model in which multiple U2AF65 molecules bound to the intron could enhance U2 snRNP recruitment to the branch point sequence.     

The Drosophila homologue of the 64 kDa subunit of cleavage stimulation factor interacts with the 77 kDa subunit encoded by the suppressor of forked gene

During mRNA 3′ end formation, cleavage stimulation factor (CstF) binds to a GU-rich sequence downstream from the polyadenylation site and helps to stabilise the binding of cleavage-polyadenylation specificity factor (CPSF) to the upstream polyadenylation sequence (AAUAAA). The 64 kDa subunit of CstF (CstF-64) contains an RNA binding domain and is responsible for the RNA binding activity of CstF. It interacts with CstF-77, which in turn interacts with CPSF. The Drosophila suppressor of forked gene encodes a homologue of CstF-77, and mutations in it affect mRNA 3′ end formation in vivo. A Drosophila homologue for CstF-64 has now been isolated, both through homology with the human protein and through protein–protein interaction in yeast with the suppressor of forked gene product. Alignment of CstF-64 homologues shows that the proteins have a conserved N-terminal 200 amino acids, the first half of which is the RNA binding domain with the second half likely to contain the CstF-77 interaction domain; a central region variable in length and rich in glycine, proline and glutamine residues and containing an unusual degenerate repeat motif; and then a conserved C-terminal 50 amino acids. In Drosophila, the CstF-64 gene has a single 63 bp intron, is transcribed throughout development and probably corresponds to l(3)91Cd.     

The mammalian homologue of Prp16p is overexpressed in a cell line tolerant to Leflunomide, a new immunoregulatory drug effective against rheumatoid arthritis.

Prp2p, Prp16p, Prp22p, and Prp43p are members of the DEAH-box family of ATP-dependent putative RNA helicases required for pre-mRNA splicing in Saccharomyces cerevisiae. Recently, mammalian homologues of Prp43p and Prp22p have been described, supporting the idea that splicing in yeast and man is phylogenetically conserved. In this study, we show that a murine cell line resistant to the novel immunoregulatory drug Leflunomide (Arava) overexpresses a 135-kDa protein that is a putative DEAH-box RNA helicase. We have cloned the human counterpart of this protein and show that it shares pronounced sequence homology with Prp16p. Apart from its N-terminal domain, which is rich in RS, RD, and RE dipeptides, this human homologue of Prp16p (designated hPrp16p) is 41% identical to Prp16p. Significantly, homology is not only observed within the phylogenetically conserved helicase domain, but also in Prp16p-specific sequences. Immunofluorescence microscopy studies demonstrated that hPrp16p co-localizes with snRNPs in subnuclear structures referred to as speckles. Antibodies specific for hPrp16p inhibited pre-mRNA splicing in vitro prior to the second step. Thus, like its yeast counterpart, hPrp16p also appears to be required for the second catalytic step of splicing. Taken together, our data indicate that the human 135-kDa protein identified here is the structural and functional homologue of the yeast putative RNA helicase, Prp16p.     

Characterization of the RNA-binding domains in the replicase proteins of tomato bushy stunt virus

Tomato bushy stunt virus (TBSV), a tombusvirus with a nonsegmented, plus-stranded RNA genome, codes for two essential replicase proteins. The sequence of one of the replicase proteins, namely p33, overlaps with the N-terminal domain of p92, which contains the signature motifs of RNA-dependent RNA polymerases (RdRps) in its nonoverlapping C-terminal portion. In this work, we demonstrate that both replicase proteins bind to RNA in vitro based on gel mobility shift and surface plasmon resonance measurements. We also show evidence that the binding of p33 to single-stranded RNA (ssRNA) is stronger than binding to double-stranded RNA (dsRNA), ssDNA, or dsDNA in vitro. Competition experiments with ssRNA revealed that p33 binds to a TBSV-derived sequence with higher affinity than to other nonviral ssRNA sequences. Additional studies revealed that p33 could bind to RNA in a cooperative manner. Using deletion derivatives of the Escherichia coli-expressed recombinant proteins in gel mobility shift and Northwestern assays, we demonstrate that p33 and the overlapping domain of p92, based on its sequence identity with p33, contain an arginine- and proline-rich RNA-binding motif (termed RPR, which has the sequence RPRRRP). This motif is highly conserved among tombusviruses and related carmoviruses, and it is similar to the arginine-rich motif present in the Tat trans-activator protein of human immunodeficiency virus type 1. We also find that the nonoverlapping C-terminal domain of p92 contains additional RNA-binding regions. Interestingly, the location of one of the RNA-binding domains in p92 is similar to the RNA-binding domain of the NS5B RdRp protein of hepatitis C virus.     

Structural similarity between the DnaA-binding proteins HobA (HP1230) from Helicobacter pylori and DiaA from Escherichia coli

In prokaryotes, DNA replication is initiated by the binding of DnaA to the oriC region of the chromosome to load the primosome machinery and start a new replication round. Several proteins control these events in Escherichia coli to ensure that replication is precisely timed during the cell cycle. Here, we report the crystal structure of HobA (HP1230) at 1.7 A, a recently discovered protein that specifically interacts with DnaA protein from Helicobacter pylori (HpDnaA). We found that the closest structural homologue of HobA is a sugar isomerase (SIS) domain containing protein, the phosphoheptose isomerase from Pseudomonas aeruginosa. Remarkably, SIS proteins share strong sequence homology with DiaA from E. coli; yet, HobA and DiaA share no sequence homology. Thus, by solving the structure of HobA, we unexpectedly discovered that HobA is a H. pylori structural homologue of DiaA. By comparing the structure of HobA to a homology model of DiaA, we identified conserved, surface-accessible residues that could be involved in protein-protein interaction. Finally, we show that HobA specifically interacts with the N-terminal part of HpDnaA. The structural homology between DiaA and HobA strongly supports their involvement in the replication process and these proteins could define a new structural family of replication regulators in bacteria.     

Solution structure of the GUCT domain from human RNA helicase II/Gu beta reveals the RRM fold, but implausible RNA interactions

Human RNA helicase II/Gu alpha (RH-II/Gu alpha) and RNA helicase II/Gu beta (RH-II/Gu beta) are paralogues that share the same domain structure, consisting of the DEAD box helicase domain (DEAD), the helicase conserved C-terminal domain (helicase_C), and the GUCT domain. The N-terminal regions of the RH-II/Gu proteins, including the DEAD domain and the helicase_C domain, unwind double-stranded RNAs. The C-terminal tail of RH-II/Gu alpha, which follows the GUCT domain, folds a single RNA strand, while that of RH-II/Gu beta does not, and the GUCT domain is not essential for either the RNA helicase or foldase activity. Thus, little is known about the GUCT domain. In this study, we have determined the solution structure of the RH-II/Gu beta GUCT domain. Structural calculations using NOE-based distance restraints and residual dipolar coupling-based angular restraints yielded a well-defined structure with beta-alpha-alpha-beta-beta-alpha-beta topology in the region for K585-A659, while the Pfam HMM algorithm defined the GUCT domain as G571-E666. This structure-based domain boundary revealed false positives in the sequence homologue search using the HMM definition. A structural homology search revealed that the GUCT domain has the RRM fold, which is typically found in RNA-interacting proteins. However, it lacks the surface-exposed aromatic residues and basic residues on the beta-sheet that are important for the RNA recognition in the canonical RRM domains. In addition, the overall surface of the GUCT domain is fairly acidic, and thus the GUCT domain is unlikely to interact with RNA molecules. Instead, it may interact with proteins via its hydrophobic surface around the surface-exposed tryptophan.     

Schistosoma mansoni: SmLIMPETin, a member of a novel family of invertebrate-only regulatory proteins

Eukaryotic LIM domain proteins contain zinc finger forming motifs rich in cysteine and histidine that enable them to interact with other proteins. A cDNA clone isolated from an adult schistosome cDNA library revealed a sequence that coded for a novel class of proteins bearing 6 LIM domains and an N-terminal PET domain, SmLIMPETin. Phylogeny reconstruction of SmLIMPETin and comparison of its sequence to invertebrate homologues and to the vertebrate four-and-a-half LIM domains protein family (FHLs), uncovered a novel LIM domain protein family, the invertebrate LIM and PET domain protein family (LIMPETin). Northern blots, RT-PCR and Western blot showed that SmLIMPETin gene was less expressed in sexually mature adult females compared to sexually immature adult females and sexually mature and immature adult males, and not expressed in schistosomula.     

Altered intracellular localization and mobility of SBDS protein upon mutation in Shwachman-Diamond syndrome

Shwachman-Diamond Syndrome (SDS) is a rare inherited disease caused by mutations in the SBDS gene. Hematopoietic defects, exocrine pancreas dysfunction and short stature are the most prominent clinical features. To gain understanding of the molecular properties of the ubiquitously expressed SBDS protein, we examined its intracellular localization and mobility by live cell imaging techniques. We observed that SBDS full-length protein was localized in both the nucleus and cytoplasm, whereas patient-related truncated SBDS protein isoforms localize predominantly to the nucleus. Also the nucleo-cytoplasmic trafficking of these patient-related SBDS proteins was disturbed. Further studies with a series of SBDS mutant proteins revealed that three distinct motifs determine the intracellular mobility of SBDS protein. A sumoylation motif in the C-terminal domain, that is lacking in patient SBDS proteins, was found to play a pivotal role in intracellular motility. Our structure-function analyses provide new insight into localization and motility of the SBDS protein, and show that patient-related mutant proteins are altered in their molecular properties, which may contribute to the clinical features observed in SDS patients.