De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Identification of a novel C2 domain factor in ovaries of orange-spotted grouper (Epinephelus coioides)

Follicle consists of an oocyte and a lot of surrounding follicular cells, and significant interactions exist between the oocyte and the somatic cells. In this study, a novel cDNA has been screened from a subtractive cDNA library between tail bud embryos and blastula embryos in the protogynous hermaphrodite orange-spotted grouper (Epinephelus coioides). Its full-length cDNA is 821 bp, and has an ORF of 414 by for encoding a peptide of 137 aa, which shows 38%, 37%, 33%, and 33% homology with 4 putative proteins screened from zebrafish (Danio rerio). Conserved domain search in NCBI reveals a single C2 domain existing in the C2 domain superfamily proteins, and has only 7 beta strands in comparison with 8 beta strands of C2 domains in other C2 domain superfamily proteins. Artificial sex reversal, RT-PCR analysis and Western blot detection demonstrated ovary-specific expression of the C2 domain factor, and therefore the novel gene was designated as E. coioides ovary-specific C2 domain factor, EcOC2 factor. Moreover, predominant expression of EcOC2 factor was further revealed in grouper mature ovary, and its strong immunofluorescence signals were located between granulosa cells and oocyte zona radiata in grouper mature follicles. The data indicate that the novel EcOC2 factor might be a main component that associates between granulosa cells and the oocyte during oocyte maturation, and might play significant roles in regulating oocyte maturation and ovulation. Further studies on its developmental behaviour and physiological functions will elucidate the interactions between oocyte and the surrounding somatic cells and the underlying molecular mechanisms. (C) 2005 Elsevier Inc. All rights reserved.

De novo RNA synthesis and homology modeling of the classical swine fever virus RNA polymerase

Classical swine fever virus (CSFV) non-structural protein 5B (NS5B) encodes an RNA-dependent RNA polymerase (RdRp), a key enzyme which initiates RNA replication by a de novo mechanism without a primer and is a potential target for anti-virus therapy. We expressed the NS5B protein in Escherichia coli. The rGTP can stimulate de novo initiation of RNA synthesis and mutation of the GDD motif to Gly-Asp-Asp (GAA) abolishes the RNA synthesis. To better understand the mechanism of viral RNA synthesis in CSFV, a three-dimensional model was built by homology modeling based on the alignment with several virus RdRps. The model contains 605 residues folded in the characteristic fingers, palm and thumb domains. The fingers domain contains an N-terminal region that plays an important role in conformational change. We propose that the experimentally observed promotion of polymerase efficiency by rGTP is probably due to the conformational changes of the polymerase caused by binding the rGTP. Mutation of the GDD to GAA interferes with the interaction between the residues at the polymerase active site and metal ions, and thus renders the polymerase inactive. (c) 2005 Elsevier B.V. All rights reserved.

A novel C1q-domain-containing protein from Zhikong scallop Chlamys farreri with lipopolysaccharide binding activity

The C1q-domain-containing (C1qDC) proteins are a family of proteins characterized by a globular C1q (gC1q) domain in their C-terminus. They are involved in various processes of vertebrates and supposed to be an important pattern recognition receptor in innate immunity of invertebrates. In this study, a novel member of C1q-domain-containing protein family was identified from Zhikong scallop Chlamys farreri (designated as CfC1qDC) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of CfC1qDC was of 777 bp, consisting of a T-terminal untranslated region (UTR) of 62 bp and a 3' UTR of 178 bp with a polyadenylation signal sequence AATAAA and a poly (A) tail. The CfC1qDC cDNA encoded a polypeptide of 178 amino acids, including a signal peptide and a C1q-domain of 158 amino acids with the theoretical isoelectric point of 5.19 and the predicted molecular weight of 17.2 kDa. The C1q-domain in CfC1qDC exhibited homology with those in sialic acid binding lectin from mollusks and C1qDC proteins from higher vertebrates. The typical 10 beta-strand jelly-roll folding topology structure of C1q-domain and the residues essential for effective packing of the hydrophobic core were well conserved in CfC1qDC. By fluorescent quantitative real-time PCR, mRNA transcripts of CfC1qDC were mainly detected in kidney, mantle, adductor muscle and gill, and also marginally detectable in hemocytes. In the bacterial challenge experiment, after the scallops were challenged by Listonella anguillarum, there was a significant up-regulation in the relative expression level of CfC1qDC and at 6 h post-injection, the mRNA expression reached the maximum level and was 4.55-fold higher than that of control scallops. Similarly, the expression of CfC1qDC mRNA in mixed primary cultures of hemocytes stimulated by lipopolysaccharides (LPS) was up-regulated and reached the maximum level at 6 h post-stimulation, and then dropped back to the original level gradually. In order to investigate its function, the cDNA fragment encoding the mature peptide of CfC1qDC was recombined and expressed in Escherichia coli BL21 (DE3). The recombinant CfC1qDC protein displayed a significantly strong activity to bind LIDS from E. coli, although no obvious antibacterial or agglutinating activity toward Gram-negative bacteria E. coli JM109, L. anguillarum and Gram-positive bacteria Micrococcus luteus was observed. These results suggested that CfC1qDC was absolutely a novel member of the C1qDC protein family and was involved in the recognition of invading microorganisms probably as a pattern recognition molecule in mollusk. (c) 2008 Elsevier Ltd. All rights reserved.

AiC1qDC-1, a novel gC1q-domain-containing protein from bay scallop Argopecten irradians with fungi agglutinating activity

The globular C1q-domain-containing (C1qDC) proteins are a family of versatile pattern recognition receptors via their globular C1q (gC1q) domain to bind various ligands including several PAMPs on pathogens. In this study, a new gC1q-domain-containing protein (AiC1qDC-1) gene was cloned from Argopecten irradians by rapid amplification of cDNA ends (RACE) approaches and expressed sequence tag (EST) analysis. The full-length cDNA of AiC1qDC-1 was composed of 733 bp, encoding a signal peptide of 19 residues and a typical gC1q domain of 137 residues containing all eight invariant amino acids in human C1qDC proteins and seven aromatic residues essential for effective packing of the hydrophobic core of AiC1qDC-1. The gC1q domain of AiC1qDC-1, which possessed the typical 10-stranded beta-sandwich fold with a jelly-roll topology common to all C1q family members, showed high homology not only to those of Cl qDC proteins in mollusk but also to those of C1qDC proteins in human. The AiC1qDC-1 transcripts were mainly detected in the tissue of hepatopancreas and also marginally detectable in adductor, heart, mantle, gill and hemocytes by fluorescent quantitative real-time PCR. In the microbial challenge experiment, there was a significant up-regulation in the relative expression level of AiC1qDC-1 in hepatopancreas and hemocytes of the scallops challenged by fungi Pichia pastoris GS115, Gram-positive bacteria Micrococcus luteus and Gram-negative bacteria Listonella anguillarum. The recombinant AiC1qDC-1 (rAiC1qDC-1) protein displayed no obvious agglutination against M. luteus and L. anguillarum, but it aggregated P. pastoris remarkably. This agglutination could be inhibited by D-mannose and PGN but not by LPS, glucan or D-galactose. These results indicated that AiC1qDC-1 functioned as a pattern recognition receptor in the immune defense of scallops against pathogens and provided clues for illuminating the evolution of the complement classical pathway. (C) 2010 Elsevier Ltd. All rights reserved.

15N NMR relaxation data reveal significant chemical exchange broadening in the alpha-domain of human α-lactalbumin

Human alpha-lactalbumin (alpha-LA), a 123-residue calcium-binding protein, has been studied using (15)N NMR relaxation methods in order to characterize backbone dynamics of the native state at the level of individual residues. Relaxation data were collected at three magnetic field strengths and analyzed using the model-free formalism of Lipari and Szabo. The order parameters derived from this analysis are generally high, indicating a rigid backbone. A total of 46 residues required an exchange contribution to T(2); 43 of these residues are located in the alpha-domain of the protein. The largest exchange contributions are observed in the A-, B-, D-, and C-terminal 3(10)-helices of the alpha-domain; these residues have been shown previously to form a highly stable core in the alpha-LA molten globule. The observed exchange broadening, along with previous hydrogen/deuterium amide exchange data, suggests that this part of the alpha-domain may undergo a local structural transition between the well-ordered native structure and a less-ordered molten-globule-like structure.

Structure of the adenylation domain of an NAD+-dependent DNA ligase

Background: DNA ligases catalyse phosphodiester bond formation between adjacent bases in nicked DNA, thereby sealing the nick. A key step in the catalytic mechanism is the formation of an adenylated DNA intermediate. The adenyl group is derived from either ATP (in eucaryotes and archaea) or NAD+4 (in bacteria). This difference in cofactor specificity suggests that DNA ligase may be a useful antibiotic target.

Results: The crystal structure of the adenylation domain of the NAD+-dependent DNA ligase from Bacillus stearothermophilus has been determined at 2.8 Å resolution. Despite a complete lack of detectable sequence similarity, the fold of the central core of this domain shares homology with the equivalent region of ATP-dependent DNA ligases, providing strong evidence for the location of the NAD+-binding site.

Conclusions: Comparison of the structure of the NAD+4-dependent DNA ligase with that of ATP-dependent ligases and mRNA-capping enzymes demonstrates the manifold utilisation of a conserved nucleotidyltransferase domain within this family of enzymes. Whilst this conserved core domain retains a common mode of nucleotide binding and activation, it is the additional domains at the N terminus and/or the C terminus that provide the alternative specificities and functionalities in the different members of this enzyme superfamily.

Functional analysis of the glycero-manno-heptose 7-phosphate kinase domain from the bifunctional HldE protein, which is involved in ADP-L-glycero-D-manno-heptose biosynthesis

The core oligosaccharide component of the lipopolysaccharide can be subdivided into inner and outer core regions. In Escherichia coli, the inner core consists of two 3-deoxy-d-manno-octulosonic acid and three glycero-manno-heptose residues. The HldE protein participates in the biosynthesis of ADP-glycero-manno-heptose precursors used in the assembly of the inner core. HldE comprises two functional domains: an N-terminal region with homology to the ribokinase superfamily (HldE1 domain) and a C-terminal region with homology to the cytidylyltransferase superfamily (HldE2 domain). We have employed the structure of the E. coli ribokinase as a template to model the HldE1 domain and predict critical amino acids required for enzyme activity. Mutation of these residues renders the protein inactive as determined in vivo by functional complementation analysis. However, these mutations did not affect the secondary or tertiary structure of purified HldE1, as judged by fluorescence spectroscopy and circular dichroism. Furthermore, in vivo coexpression of wild-type, chromosomally encoded HldE and mutant HldE1 proteins with amino acid substitutions in the predicted ATP binding site caused a dominant negative phenotype as revealed by increased bacterial sensitivity to novobiocin. Copurification experiments demonstrated that HldE and HldE1 form a complex in vivo. Gel filtration chromatography resulted in the detection of a dimer as the predominant form of the native HldE1 protein. Altogether, our data support the notions that the HldE functional unit is a dimer and that structural components present in each HldE1 monomer are required for enzymatic activity.

Inactivating mutations in an SH2 domain-encoding gene in X-linked lymphoproliferative syndrome

X-linked lymphoproliferative syndrome (XLP) is an inherited immunodeficiency characterized by increased susceptibility to Epstein-Barr virus (EBV). In affected males, primary EBV infection leads to the uncontrolled proliferation of virus-containing B cells and reactive cytotoxic T cells, often culminating in the development of high-grade lymphoma. The XLP gene has been mapped to chromosome band Xq25 through linkage analysis and the discovery of patients harboring large constitutional genomic deletions. We describe here the presence of small deletions and intragenic mutations that specifically disrupt a gene named DSHP in 6 of 10 unrelated patients with XLP. This gene encodes a predicted protein of 128 amino acids composing a single SH2 domain with extensive homology to the SH2 domain of SHIP, an inositol polyphosphate 5-phosphatase that functions as a negative regulator of lymphocyte activation. DSHP is expressed in transformed T cell lines and is induced following in vitro activation of peripheral blood T lymphocytes. Expression of DSHP is restricted in vivo to lymphoid tissues, and RNA in situ hybridization demonstrates DSHP expression in activated T and B cell regions of reactive lymph nodes and in both T and B cell neoplasms. These observations confirm the identity of DSHP as the gene responsible for XLP, and suggest a role in the regulation of lymphocyte activation and proliferation. Induction of DSHP may sustain the immune response by interfering with SHIP-mediated inhibition of lymphocyte activation, while its inactivation in XLP patients results in a selective immunodeficiency to EBV.

Isolation and analysis of a corprinus cinereus DNA fragment showing homology to a fimbrial cDNA from ustilago violacea

Surface proteinaceous fibrils, termed fimbriae, were first identified on gram negative bacteria in the 1940s. Fungal fimbriae, discovered some 25 years later, are found on members of all fungal classes. In the present study, polyclonal antiserum raised against the fimbrial proteins of U. vio/acea were used in order to identify antigenically related proteins from Coprinus cinereus and Schizophy//um commune. Two polypeptides with molecular masses of 37 and 39 kDa from C. cinereus were observed and confirm earlier results. A single previously unidentified 50 kDa polypeptide in S. commune crossreacted with the antiserum. The 50 kDa protein was found to consist of 3 isoforms with isoelectric points ranging from 5.6 to 5.8. A fimbrial cDNA derived from U. vio/acea was used to identify DNA restriction fragments from C. cinereus and S. commune showing homology to the fimbrial transcript of U. vio/acea. Heterologous hybridization with this cDNA was used in order to screen a C. cinereus genomic DNA library. A single clone, A2-3A, with a 14 kbp insert showed strong homology to the pfim3-1 cDNA. The region of homology, a 700 bp Xba I fragment, was subcloned into pUG19. This plasmid was refered to as pXX8. DNA sequence determinations of pXX8 and adjacent fragments from A2-3A suggested that the cloned DNA was a portion of the rONA repeat encoding the small subunit rRNA. DNA sequence analysis of pfim3-1 yielded an incomplete open reading frame. The predicted amino acid sequence codes for a 206 amino acid, 22 kDa polypeptide which contains a domain similar to a transmembrane domain from rat leukocyte antigen, GDS3. As well, an untranslated 576 nucleotide domain showed 81 % homology to pXX8 and 830/0 homology to the 188 rRNA sequence of Ustilago maydis. This sequence was found adjacent to a region of adenine-thymine base pairs presumed to represent the polyadenylation sequence of the fimbrial transcript. The size and extent of homology is sufficient to account for the hybridization of pfim3-1 to rDNA. It is suggested that this domain represents a completely novel regulatory domain within eukaryotes that may enable the observed rapid regeneration of fimbriae in U. violacea.

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

The nuclear magnetic resonance (NMR) structure of a central segment of the previously annotated severe acute respiratory syndrome (SARS)-unique domain (SUD-M, for "middle of the SARS-unique domain") in SARS coronavirus (SARS-CoV) nonstructural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3 residues 528 to 648, and there is a flexibly extended N-terminal tail with the residues 513 to 527 and a C-terminal flexible tail of residues 649 to 651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527 to 651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly(A) and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In a further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1 ''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows three-dimensional structure homology with several helicases and nucleoside triphosphate-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

NMR structure shows that the SARS-unique domain contains a macrodomain fold

The NMR structure of a central segment of the previously annotated "SARS-unique domain" (SUD-M; "middle of the SARS-unique domain") in the SARS coronavirus (SARS-CoV) non-structural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3-residues 528-648, and there is a flexibly extended N-terminal tail with the residues 513-527 and a C-terminal flexible tail of residues 649-651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527-651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly-A and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows 3D structure homology with several helicases and NTP-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

Hepatitis C virus NS5A protein binds the SH3 domain of the Fyn tyrosine kinase with high affinity: mutagenic analysis of residues within the SH3 domain that contribute to the interaction

Background: The hepatitis C virus (HCV) non-structural 5A protein (NS5A) contains a highly conserved C-terminal polyproline motif with the consensus sequence Pro-X-X- Pro-X-Arg that is able to interact with the Src-homology 3 (SH3) domains of a variety of cellular proteins. Results: To understand this interaction in more detail we have expressed two N-terminally truncated forms of NS5A in E. coli and examined their interactions with the SH3 domain of the Src-family tyrosine kinase, Fyn. Surface plasmon resonance analysis revealed that NS5A binds to the Fyn SH3 domain with what can be considered a high affinity SH3 domain-ligand interaction (629 nM), and this binding did not require the presence of domain I of NS5A (amino acid residues 32-250). Mutagenic analysis of the Fyn SH3 domain demonstrated the requirement for an acidic cluster at the C-terminus of the RT-Src loop of the SH3 domain, as well as several highly conserved residues previously shown to participate in SH3 domain peptide binding. Conclusion: We conclude that the NS5A: Fyn SH3 domain interaction occurs via a canonical SH3 domain binding site and the high affinity of the interaction suggests that NS5A would be able to compete with cognate Fyn ligands within the infected cell.

Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.