The HHpred interactive server for protein homology detection and structure prediction

Drosophila multiplexin (Dmp) modulates motor axon pathfinding accuracy

Multiplexins are multidomain collagens typically composed of an N-terminal thrombospondin-related domain, an interrupted triple helix and a C-terminal endostatin domain. They feature a clear regulatory function in the development of different tissues, which is chiefly conveyed by the endostatin domain. This domain can be found in proteolytically released monomeric and trimeric versions, and their diverse and opposed effects on the migratory behavior of epithelial and endothelial cell types have been demonstrated in cell culture experiments. The only Drosophila multiplexin displays specific features of both vertebrate multiplexins, collagens XV and XVIII. We characterized the Drosophila multiplexin (dmp) gene and found that three main isoforms are expressed from it, one of which is the monomeric endostatin version. Generation of dmp deletion alleles revealed that Dmp plays a role in motor axon pathfinding, as the mutants exhibit ventral bypass defects of the intersegmental nerve b (ISNb) similar to other motor axon guidance mutants. Transgenic overexpression of monomeric endostatin as well as of full-length Dmp, but not trimeric endostatin, were able to rescue these defects. In contrast, trimeric endostatin increased axon pathfinding accuracy in wild type background. We conclude that Dmp plays a modulating role in motor axon pathfinding and may be part of a buffering system that functions to avoid innervation errors.

COMPASS server for homology detection: improved statistical accuracy, speed and functionality

COMPASS is a profile-based method for the detection of remote sequence similarity and the prediction of protein structure. Here we describe a recently improved public web server of COMPASS, http://prodata.swmed.edu/compass. The server features three major developments: (i) improved statistical accuracy; (ii) increased speed from parallel implementation; and (iii) new functional features facilitating structure prediction. These features include visualization tools that allow the user to quickly and effectively analyze specific local structural region predictions suggested by COMPASS alignments. As an application example, we describe the structural, evolutionary and functional analysis of a protein with unknown function that served as a target in the recent CASP8 (Critical Assessment of Techniques for Protein Structure Prediction round 8). URL: http://prodata.swmed.edu/compass

Vaccinia virus G8R protein: a structural ortholog of proliferating cell nuclear antigen (PCNA)

Background

Eukaryotic DNA replication involves the synthesis of both a DNA leading and lagging strand, the latter requiring several additional proteins including flap endonuclease (FEN-1) and proliferating cell nuclear antigen (PCNA) in order to remove RNA primers used in the synthesis of Okazaki fragments. Poxviruses are complex viruses (dsDNA genomes) that infect eukaryotes, but surprisingly little is known about the process of DNA replication. Given our previous results that the vaccinia virus (VACV) G5R protein may be structurally similar to a FEN-1-like protein and a recent finding that poxviruses encode a primase function, we undertook a series of in silico analyses to identify whether VACV also encodes a PCNA-like protein.

Results

An InterProScan of all VACV proteins using the JIPS software package was used to identify any PCNA-like proteins. The VACV G8R protein was identified as the only vaccinia protein that contained a PCNA-like sliding clamp motif. The VACV G8R protein plays a role in poxvirus late transcription and is known to interact with several other poxvirus proteins including itself. The secondary and tertiary structure of the VACV G8R protein was predicted and compared to the secondary and tertiary structure of both human and yeast PCNA proteins, and a high degree of similarity between all three proteins was noted.

Conclusions

The structure of the VACV G8R protein is predicted to closely resemble the eukaryotic PCNA protein; it possesses several other features including a conserved ubiquitylation and SUMOylation site that suggest that, like its counterpart in T4 bacteriophage (gp45), it may function as a sliding clamp ushering transcription factors to RNA polymerase during late transcription.

webPRC: the Profile Comparer for alignment-based searching of public domain databases

Profile–profile methods are well suited to detect remote evolutionary relationships between protein families. Profile Comparer (PRC) is an existing stand-alone program for scoring and aligning hidden Markov models (HMMs), which are based on multiple sequence alignments. Since PRC compares profile HMMs instead of sequences, it can be used to find distant homologues. For this purpose, PRC is used by, for example, the CATH and Pfam-domain databases. As PRC is a profile comparer, it only reports profile HMM alignments and does not produce multiple sequence alignments. We have developed webPRC server, which makes it straightforward to search for distant homologues or similar alignments in a number of domain databases. In addition, it provides the results both as multiple sequence alignments and aligned HMMs. Furthermore, the user can view the domain annotation, evaluate the PRC hits with the Jalview multiple alignment editor and generate logos from the aligned HMMs or the aligned multiple alignments. Thus, this server assists in detecting distant homologues with PRC as well as in evaluating and using the results. The webPRC interface is available at http://www.ibi.vu.nl/programs/prcwww/.

HorA web server to infer homology between proteins using sequence and structural similarity

The biological properties of proteins are often gleaned through comparative analysis of evolutionary relatives. Although protein structure similarity search methods detect more distant homologs than purely sequence-based methods, structural resemblance can result from either homology (common ancestry) or analogy (similarity without common ancestry). While many existing web servers detect structural neighbors, they do not explicitly address the question of homology versus analogy. Here, we present a web server named HorA (Homology or Analogy) that identifies likely homologs for a query protein structure. Unlike other servers, HorA combines sequence information from state-of-the-art profile methods with structure information from spatial similarity measures using an advanced computational technique. HorA aims to identify biologically meaningful connections rather than purely 3D-geometric similarities. The HorA method finds approximately 90% of remote homologs defined in the manually curated database SCOP. HorA will be especially useful for finding remote homologs that might be overlooked by other sequence or structural similarity search servers. The HorA server is available at http://prodata.swmed.edu/horaserver.

Structural insights into chondroitin sulphate A binding Duffy-binding-like domains from Plasmodium falciparum: implications for intervention strategies against placental malaria

BACKGROUND

Placental malaria is typified by selective clustering of Plasmodium falciparum in the intervillous blood spaces of the placenta. Sequestration of malaria parasite in the human placenta is mediated by interactions between chondroitin sulphate A (CSA) on the syncytiotrophoblasts and proteins expressed on the surface of infected human erythrocytes. Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) encoded by the var2CSA gene is believed to be the main parasite ligand for CSA-mediated placental binding.

METHODS

Extensive sequence and structure comparisons of the various CSA-binding and non-binding DBL domains from the var2CSA gene from A4 and 3D7 strains of P. falciparum were performed. Three-dimensional structural models of various DBL domains were built and analysed with a view to assessing conservation of CSA interaction sites across various DBL domains.

RESULTS

Each of the six DBL domains from var2CSA are likely to retain the disulfide linkages evident from previously published DBL domain crystal structures. The number of disulfide linkages between the various DBL domains analysed varies from three to seven, of which two are conserved across all DBL domains. The conserved disulfide linkages are distributed within the respective three sub-domains and only one linkage is shared by sub-domains I and II. Major differences between CSA-binding DBL domains are in the loop regions, which tie the alpha helices together, and in variable length terminal extensions. Intriguingly, a crucial loop from A4 DBL 3X which provides the important Gly and Lys residues that chelate the bound sulphate is missing or significantly altered in all other DBL domains that interact with CSA. Further analysis of the proposed sulphate and predicted CSA-binding site indicates either none or very low level of conservation among the critical interacting residues.

CONCLUSION

Structural comparisons of the three-dimensional structures of CSA-binding DBL domains indicates that the proposed CSA interaction site on A4 DBL 3X is unlikely to be conserved across the other CSA-binding DBL domains from var2CSA. Therefore, the 4 CSA-binding DBL domains encoded by var2CSA are unlikely to have common architectures to their CSA recognition sites. These structural insights have clear implications in using CSA-binding DBL domains for vaccines against placental malaria as it is proposed that the various CSA-binding DBL domains on var2CSA will recognize their CSA ligands differently.

Analysis of CASP8 targets, predictions and assessment methods

Results of the recent Critical Assessment of Techniques for Protein Structure Prediction, CASP8, present several valuable sources of information. First, CASP targets comprise a realistic sample of currently solved protein structures and exemplify the corresponding challenges for predictors. Second, the plethora of predictions by all possible methods provides an unusually rich material for evolutionary analysis of target proteins. Third, CASP results show the current state of the field and highlight specific problems in both predicting and assessing. Finally, these data can serve as grounds to develop and analyze methods for assessing prediction quality. Here we present results of our analysis in these areas. Our objective is not to duplicate CASP assessment, but to use our unique experience as former CASP5 assessors and CASP8 predictors to (i) offer more insights into CASP targets and predictions based on expert analysis, including invaluable analysis prior to target structure release; and (ii) develop an assessment methodology tailored towards current challenges in the field. Specifically, we discuss preparing target structures for assessment, parsing protein domains, balancing evaluations based on domains and on whole chains, dividing targets into categories and developing new evaluation scores. We also present evolutionary analysis of the most interesting and challenging targets.

Database URL: Our results are available as a comprehensive database of targets and predictions at http://prodata.swmed.edu/CASP8.

ORF-C4 from the early branching eukaryote Giardia lamblia displays characteristics of alpha-crystallin small heat-shock proteins

Giardia lamblia is a medically important protozoan parasite with a basal position in the eukaryotic lineage and is an interesting model to explain the evolution of biochemical events in eukaryotic cells. G. lamblia trophozoites undergo significant changes in order to survive outside the intestine of their host by differentiating into infective cysts. In the present study, we characterize the previously identified Orf-C4 (G. lamblia open reading frame C4) gene, which is considered to be specific to G. lamblia. It encodes a 22 kDa protein that assembles into high-molecular-mass complexes during the entire life cycle of the parasite. ORF-C4 localizes to the cytoplasm of trophozoites and cysts, and forms large spherical aggregates when overexpressed. ORF-C4 overexpression and down-regulation do not affect trophozoite viability; however, differentiation into cysts is slightly delayed when the expression of ORF-C4 is down-regulated. In addition, ORF-C4 protein expression is modified under specific stress-inducing conditions. Neither orthologous proteins nor conserved domains are found in databases by conventional sequence analysis of the predicted protein. However, ORF-C4 contains a region which is similar structurally to the alpha-crystallin domain of sHsps (small heat-shock proteins). In the present study, we show the potential role of ORF-C4 as a small chaperone which is involved in the response to stress (including encystation) in G. lamblia.

Multiple Motif Scanning to identify methyltransferases from the yeast proteome

A new program (Multiple Motif Scanning) was developed to scan the Saccharomyces cerevisiae proteome for Class I S-adenosylmethionine-dependent methyltransferases. Conserved Motifs I, Post I, II, and III were identified and expanded in known methyltransferases by primary sequence and secondary structural analysis through hidden Markov model profiling of both a yeast reference database and a reference database of methyltransferases with solved three-dimensional structures. The roles of the conserved amino acids in the four motifs of the methyltransferase structure and function were then analyzed to expand the previously defined motifs. Fisher-based negative log statistical matrix sets were developed from the prevalence of amino acids in the motifs. Multiple Motif Scanning is able to scan the proteome and score different combinations of the top fitting sequences for each motif. In addition, the program takes into account the conserved number of amino acids between the motifs. The output of the program is a ranked list of proteins that can be used to identify new methyltransferases and to reevaluate the assignment of previously identified putative methyltransferases. The Multiple Motif Scanning program can be used to develop a putative list of enzymes for any type of protein that has one or more motifs conserved at variable spacings and is freely available (www.chem.ucla.edu/files/MotifSetup.Zip). Finally hidden Markov model profile clustering analysis was used to subgroup Class I methyltransferases into groups that reflect their methyl-accepting substrate specificity.

Modeling of glutamate GluR6 receptor and its interactions with novel noncompetitive antagonists

The study proposes the first complete model of an ionotropic glutamate receptor (GluR6). The model is in accordance with available experimental data from single-particle electron microscopy images and exhibits correct shape and dimensions and the appropriate symmetry: 2-fold in the N-terminal domain (NTD), ligand-binding domain (LBD), and external part of the transmembrane region, whereas it is 4-fold deeper in the channel. The methodology applied for GluR6 receptor model building was validated in the docking procedure of competitive and uncompetitive antagonists. The constructed model was used to study molecular interactions of novel noncompetitive GluR6 antagonists with their molecular target. A new binding site in the GluR6 receptor transduction domain has been identified. It is situated between two subunits in the receptor dimer. The following residues were recognized as crucial for interactions: Arg663A, Arg663B (M3-S2 linker), Ser809B (S2-M4 linker), and Phe553A (S1-M1 linker).

DNA-binding Domain within the Brh2 N Terminus Is the Primary Interaction Site for Association with DNA

The C-terminal region of Brh2 (Brh2(CT)), the BRCA2 homolog in Ustilago maydis, is highly conserved and aligns with the DSS1/DNA-binding domain (DBD) of mammalian BRCA2, while the N-terminal region (Brh2(NT)) is poorly conserved and has no obvious functional domain except for the single Rad51-interacting BRC element. Paradoxically, Brh2(NT), but not Brh2(CT), complements the DNA repair and recombination deficiency of the brh2 mutant. We show here that Brh2(NT) exhibits an unexpected DNA binding activity with properties similar to that of the full-length protein. Deletion mapping localized the region responsible for the DNA binding activity to a stretch of residues between the BRC element and the canonical DBD. A heterologous DNA-binding domain from the large subunit of replication protein A substituted for the endogenous binding region within Brh2(NT) in supporting DNA repair. Rad51-promoted strand invasion was stimulated by Brh2(NT), but required the presence of the BRC element. The findings suggest a model in which Brh2(NT) serves as the principal site for association with DNA, while the Brh2(CT) provides a means for regulation.

Evolution of DNA polymerases: an inactivated polymerase-exonuclease module in Pol epsilon and a chimeric origin of eukaryotic polymerases from two classes of archaeal ancestors

Background

Evolution of DNA polymerases, the key enzymes of DNA replication and repair, is central to any reconstruction of the history of cellular life. However, the details of the evolutionary relationships between DNA polymerases of archaea and eukaryotes remain unresolved.

Results

We performed a comparative analysis of archaeal, eukaryotic, and bacterial B-family DNA polymerases, which are the main replicative polymerases in archaea and eukaryotes, combined with an analysis of domain architectures. Surprisingly, we found that eukaryotic Polymerase ε consists of two tandem exonuclease-polymerase modules, the active N-terminal module and a C-terminal module in which both enzymatic domains are inactivated. The two modules are only distantly related to each other, an observation that suggests the possibility that Pol ε evolved as a result of insertion and subsequent inactivation of a distinct polymerase, possibly, of bacterial descent, upstream of the C-terminal Zn-fingers, rather than by tandem duplication. The presence of an inactivated exonuclease-polymerase module in Pol ε parallels a similar inactivation of both enzymatic domains in a distinct family of archaeal B-family polymerases. The results of phylogenetic analysis indicate that eukaryotic B-family polymerases, most likely, originate from two distantly related archaeal B-family polymerases, one form giving rise to Pol ε, and the other one to the common ancestor of Pol α, Pol δ, and Pol ζ. The C-terminal Zn-fingers that are present in all eukaryotic B-family polymerases, unexpectedly, are homologous to the Zn-finger of archaeal D-family DNA polymerases that are otherwise unrelated to the B family. The Zn-finger of Polε shows a markedly greater similarity to the counterpart in archaeal PolD than the Zn-fingers of other eukaryotic B-family polymerases.

Conclusion

Evolution of eukaryotic DNA polymerases seems to have involved previously unnoticed complex events. We hypothesize that the archaeal ancestor of eukaryotes encoded three DNA polymerases, namely, two distinct B-family polymerases and a D-family polymerase all of which contributed to the evolution of the eukaryotic replication machinery. The Zn-finger might have been acquired from PolD by the B-family form that gave rise to Pol ε prior to or in the course of eukaryogenesis, and subsequently, was captured by the ancestor of the other B-family eukaryotic polymerases. The inactivated polymerase-exonuclease module of Pol ε might have evolved by fusion with a distinct polymerase, rather than by duplication of the active module of Pol ε, and is likely to play an important role in the assembly of eukaryotic replication and repair complexes.

Reviewers

This article was reviewed by Patrick Forterre, Arcady Mushegian, and Chris Ponting. For the full reviews, please go to the Reviewers' Reports section.

Verification of a topology model of PorT as an integral outer-membrane protein in Porphyromonas gingivalis

PorT is a membrane-associated protein shown to be essential for the maturation and secretion of a class of cysteine proteinases, the gingipains, from the periodontal pathogen Porphyromonas gingivalis. It was previously reported that PorT is located on the periplasmic surface of the inner membrane to function as a chaperone for the maturing proteinases. Our modelling suggested it to be an integral outer-membrane protein with eight anti-parallel, membrane-traversing beta-strands. In this report, the outer-membrane localization model was confirmed by the structural and functional tolerance of PorT to hexahistidine (6xHis) tag insertions at selected locations within the protein using site-directed mutagenesis. Interestingly, those PorT mutations adversely affecting gingipain secretion enhanced expression of the porT gene but at the same time suppressed the transcription of the gingipain rgpB gene. Further, PorT mutants deficient in gingipain activities produced significantly more di- and triaminopeptidase activities. PorT homologues have been found in restricted members of the Bacteroidetes phylum where there is potential for PorT to participate in the maturation and secretion of proteins with characteristic C-terminal domains (CTDs). Knowledge of the cellular localization of PorT will enable analysis of the role of this protein in a new secretory pathway for the export of gingipains and other CTD-class proteins.

Protein structure homology modeling using SWISS-MODEL workspace

Homology modeling aims to build three-dimensional protein structure models using experimentally determined structures of related family members as templates. SWISS-MODEL workspace is an integrated Web-based modeling expert system. For a given target protein, a library of experimental protein structures is searched to identify suitable templates. On the basis of a sequence alignment between the target protein and the template structure, a three-dimensional model for the target protein is generated. Model quality assessment tools are used to estimate the reliability of the resulting models. Homology modeling is currently the most accurate computational method to generate reliable structural models and is routinely used in many biological applications. Typically, the computational effort for a modeling project is less than 2 h. However, this does not include the time required for visualization and interpretation of the model, which may vary depending on personal experience working with protein structures.

SimShiftDB; local conformational restraints derived from chemical shift similarity searches on a large synthetic database

We present SimShiftDB, a new program to extract conformational data from protein chemical shifts using structural alignments. The alignments are obtained in searches of a large database containing 13,000 structures and corresponding back-calculated chemical shifts. SimShiftDB makes use of chemical shift data to provide accurate results even in the case of low sequence similarity, and with even coverage of the conformational search space. We compare SimShiftDB to HHSearch, a state-of-the-art sequence-based search tool, and to TALOS, the current standard tool for the task. We show that for a significant fraction of the predicted similarities, SimShiftDB outperforms the other two methods. Particularly, the high coverage afforded by the larger database often allows predictions to be made for residues not involved in canonical secondary structure, where TALOS predictions are both less frequent and more error prone. Thus SimShiftDB can be seen as a complement to currently available methods.

Evidence for novel tomato seed allergens: IgE-reactive legumin and vicilin proteins identified by multidimensional protein fractionation-mass spectrometry and in silico epitope modeling

Tomato fruit and seed allergens were detected by IgE-immunoblotting using sera from 18 adult tomato-sensitized patients selected based on a positive history skin prick test (SPT) and specific Immunglobulin (Ig) E-levels. Isolated tomato seed total protein showed high SPT activity comparable or even higher than tomato fruit protein. For the molecular characterization of tomato seed allergens, a multidimensional protein fractionation strategy and LC-MS/MS was used. Two legumin- and vicilin-proteins were purified and showed strong IgE-reactivity in immunoblots. Individual patient sera exhibited varying IgE-sensitivity against the purified proteins. In silico structural modeling indicates high homology between epitopes of known walnut allergens and the detected IgE-crossreactive tomato proteins.

Type VI secretion apparatus and phage tail-associated protein complexes share a common evolutionary origin

Protein secretion is a common property of pathogenic microbes. Gram-negative bacterial pathogens use at least 6 distinct extracellular protein secretion systems to export proteins through their multilayered cell envelope and in some cases into host cells. Among the most widespread is the newly recognized Type VI secretion system (T6SS) which is composed of 15-20 proteins whose biochemical functions are not well understood. Using crystallographic, biochemical, and bioinformatic analyses, we identified 3 T6SS components, which are homologous to bacteriophage tail proteins. These include the tail tube protein; the membrane-penetrating needle, situated at the distal end of the tube; and another protein associated with the needle and tube. We propose that T6SS is a multicomponent structure whose extracellular part resembles both structurally and functionally a bacteriophage tail, an efficient machine that translocates proteins and DNA across lipid membranes into cells.

Conservation of helical bundle structure between the exocyst subunits

BACKGROUND

The exocyst is a large hetero-octomeric protein complex required for regulating the targeting and fusion of secretory vesicles to the plasma membrane in eukaryotic cells. Although the sequence identity between the eight different exocyst subunits is less than 10%, structures of domains of four of the subunits revealed a similar helical bundle topology. Characterization of several of these subunits has been hindered by lack of soluble protein for biochemical and structural studies.

METHODOLOGY/PRINCIPAL FINDINGS

Using advanced hidden Markov models combined with secondary structure predictions, we detect significant sequence similarity between each of the exocyst subunits, indicating that they all contain helical bundle structures. We corroborate these remote homology predictions by identifying and purifying a predicted domain of yeast Sec10p, a previously insoluble exocyst subunit. This domain is soluble and folded with approximately 60% alpha-helicity, in agreement with our predictions, and capable of interacting with several known Sec10p binding partners.

CONCLUSIONS/SIGNIFICANCE

Although all eight of the exocyst subunits had been suggested to be composed of similar helical bundles, this has now been validated by our hidden Markov model structure predictions. In addition, these predictions identified protein domains within the exocyst subunits, resulting in creation and characterization of a soluble, folded domain of Sec10p.

Ab initio protein modelling reveals novel human MIT domains

Database searches can fail to detect all truly homologous sequences, particularly when dealing with short, highly sequence diverse protein families. Here, using microtubule interacting and transport (MIT) domains as an example, we have applied an approach of profile-profile matching followed by ab initio structure modelling to the detection of true homologues in the borderline significant zone of database searches. Novel MIT domains were confidently identified in USP54, containing an apparently inactive ubiquitin carboxyl-terminal hydrolase domain, a katanin-like ATPase KATNAL1, and an uncharacterized protein containing a VPS9 domain. As a proof of principle, we have confirmed the novel MIT annotation for USP54 by in vitro profiling of binding to CHMP proteins.

Comparative genomic analyses of nickel, cobalt and vitamin B12 utilization

BACKGROUND

Nickel (Ni) and cobalt (Co) are trace elements required for a variety of biological processes. Ni is directly coordinated by proteins, whereas Co is mainly used as a component of vitamin B12. Although a number of Ni and Co-dependent enzymes have been characterized, systematic evolutionary analyses of utilization of these metals are limited.

RESULTS

We carried out comparative genomic analyses to examine occurrence and evolutionary dynamics of the use of Ni and Co at the level of (i) transport systems, and (ii) metalloproteomes. Our data show that both metals are widely used in bacteria and archaea. Cbi/NikMNQO is the most common prokaryotic Ni/Co transporter, while Ni-dependent urease and Ni-Fe hydrogenase, and B12-dependent methionine synthase (MetH), ribonucleotide reductase and methylmalonyl-CoA mutase are the most widespread metalloproteins for Ni and Co, respectively. Occurrence of other metalloenzymes showed a mosaic distribution and a new B12-dependent protein family was predicted. Deltaproteobacteria and Methanosarcina generally have larger Ni- and Co-dependent proteomes. On the other hand, utilization of these two metals is limited in eukaryotes, and very few of these organisms utilize both of them. The Ni-utilizing eukaryotes are mostly fungi (except saccharomycotina) and plants, whereas most B12-utilizing organisms are animals. The NiCoT transporter family is the most widespread eukaryotic Ni transporter, and eukaryotic urease and MetH are the most common Ni- and B12-dependent enzymes, respectively. Finally, investigation of environmental and other conditions and identity of organisms that show dependence on Ni or Co revealed that host-associated organisms (particularly obligate intracellular parasites and endosymbionts) have a tendency for loss of Ni/Co utilization.

CONCLUSION

Our data provide information on the evolutionary dynamics of Ni and Co utilization and highlight widespread use of these metals in the three domains of life, yet only a limited number of user proteins.

Organ-specific attenuation of murine hepatitis virus strain A59 by replacement of catalytic residues in the putative viral cyclic phosphodiesterase ns2

The Murine hepatitis virus (MHV) strain A59 ns2 protein is a 30-kDa nonstructural protein that is expressed from a subgenomic mRNA in the cytoplasm of virus-infected cells. Its homologs are also encoded in other closely related group 2a coronaviruses and more distantly related toroviruses. Together, these proteins comprise a subset of a large superfamily of 2H phosphoesterase proteins that are distinguished by a pair of conserved His-x-Thr/Ser motifs encompassing catalytically important residues. We have used a vaccinia virus-based reverse genetic system to produce recombinant viruses encoding ns2 proteins with single-amino-acid substitutions in, or adjacent to, these conserved motifs, namely, inf-ns2 H46A, inf-ns2 S48A, inf-ns2-S120A, and inf-ns2-H126R. All of the mutant viruses replicate in mouse 17 clone 1 fibroblast cells and mouse embryonic cells to the same extent as the parental wild-type recombinant virus, inf-MHV-A59. However, compared to inf-MHV-A59, the inf-ns2 H46A and inf-ns2-H126R mutants are highly attenuated for replication in mouse liver following intrahepatic inoculation. Interestingly, none of the mutant viruses were attenuated for replication in mouse brain following intracranial inoculation. These results show that the ns2 protein of MHV-A59 has an important role in virus pathogenicity and that a substitution of the histidine residues of the MHV-A59 ns2 His-x-Thr/Ser motifs is critical for virus virulence in the liver but not in the brain. This novel phenotype suggests a strategy to investigate the function of the MHV-A59 ns2 protein involving the search for organ-specific proteins or RNAs that react differentially to wild-type and mutant ns2 proteins.

HPC2 and ubinuclein define a novel family of histone chaperones conserved throughout eukaryotes

Using sensitive protein sequence profile analyses we investigated the evolution of some histone chaperones and showed that Hir3p and Hpc2p have a much wider phyletic pattern than was previously known.

While histone chaperones have been intensely studied, the roles of components of the Hir–Asf1 histone chaperone complex such as Hir3p and Hpc2p are poorly understood. Using sensitive protein sequence profile analyses we investigated the evolution of these proteins and showed that Hir3p and Hpc2p have a much wider phyletic pattern than was previously known. We established the animal histone-deacetylase-complex-interacting proteins, CAIN/CABIN, to be orthologs of Hir3p. They contain a conserved core of around 30 TPR-like bi-helical repeats that are likely to form a super-helical scaffold. We identified a conserved domain, the HUN domain, in all Hpc2p homologs, including animal ubinuclein/yemanuclein and the recently discovered vertebrate cell-cycle regulator FLJ25778. The HUN domain has a characteristic pattern of conserved acidic residues based on which we predict that it is a previously unrecognized histone-tail-binding chaperone. By analyzing various high-throughput data sets, such as RNAi knock-downs, genetic and protein interaction maps and cell-cycle-specific gene expression data, we present evidence that Hpc2p homologs might be deployed in specific processes of chromatin dynamics relating to cell-cycle progression in vertebrates and schizogony in Plasmodium. Beyond the conserved HUN domain these proteins show extensive divergence patterns in different eukaryotic lineages. Hence, we propose that Hpc2p homologs are probably involved in recruitment of the ancient conserved histone-loading Hir–Asf1 complex to different lineage-specific chromatin reorganization processes.

Protein domain boundary predictions: a structural biology perspective

One of the important fields to apply computational tools for domain boundaries prediction is structural biology. They can be used to design protein constructs that must be expressed in a stable and functional form and must produce diffraction-quality crystals. However, prediction of protein domain boundaries on the basis of amino acid sequences is still very problematical. In present study the performance of several computational approaches are compared. It is observed that the statistical significance of most of the predictions is rather poor. Nevertheless, when the right number of domains is correctly predicted, domain boundaries are predicted within very few residues from their real location. It can be concluded that prediction methods cannot be used yet as routine tools in structural biology, though some of them are rather promising.

The type VI secretion system: translocation of effectors and effector-domains

A number of prominent Gram-negative bacteria use the type VI secretion system (T6SS) to transport proteins across the bacterial envelope. Rapid progress is being made in elucidating the structural components of the T6SS apparatus, and a few effectors have been reported to pass through it. However, this is not the complete story: a family of T6SS proteins, the VgrGs, share structural features with the cell-puncturing device of the T4 bacteriophage, and may be used in a similar fashion by bacteria to puncture host cell membranes and insert the T6SS apparatus into the host cytosol. Interestingly, a number of VgrGs contain C-terminal extensions with effector-domains. Thus, the T6SS may translocate soluble effectors, as well as VgrG effector-domains.

The human Ago2 MC region does not contain an eIF4E-like mRNA cap binding motif

Background

Argonaute (Ago) proteins interact with small regulatory RNAs to mediate gene regulatory pathways. A recent report by Kiriakidou et al. [1] describes an MC sequence region identified in Ago2 that displays similarity to the cap-binding motif in translation initiation factor 4E (eIF4E). In a cap-bound eIF4E structure, two important aromatic residues of the motif stack on either side of a 7-methylguanosine 5'-triphosphate (m⁷Gppp) base. The corresponding Ago2 aromatic residues (F450 and F505) were hypothesized to perform the same cap-binding function. However, the detected similarity between the MC sequence and the eIF4E cap-binding motif was questionable.

Results

A number of sequence-based and structure-based bioinformatics methods reveal the reported similarity between the Ago2 MC sequence region and the eIF4E cap-binding motif to be spurious. Alternatively, the MC sequence region is confidently assigned to the N-terminus of the Ago piwi module, within the mid domain of experimentally determined prokaryotic Ago structures. Confident mapping of the Ago2 MC sequence region to the piwi mid domain results in a homology-based structure model that positions the identified aromatic residues over 20 Å apart, with one of the aromatic side chains (F450) contributing instead to the hydrophobic core of the domain.

Conclusion

Correct functional prediction based on weak sequence similarity requires substantial evolutionary and structural support. The evolutionary context of the Ago mid domain suggested by multiple sequence alignment is limited to a conserved hydrophobicity profile required for the fold and a motif following the MC region that binds guide RNA. Mapping of the MC sequence to the mid domain structure reveals Ago2 aromatics that are incompatible with eIF4E-like mRNA cap-binding, yet display some limited local structure similarities that cause the chance sequence match to eIF4E.

Reviewers

This article was reviewed by Arcady Mushegian, Chris Ponting, and Igor Jouline (nominated by Igor Zhulin).

Conformational flexibility of metazoan fatty acid synthase enables catalysis

The metazoan cytosolic fatty acid synthase (FAS) contains all of the enzymes required for de novo fatty acid biosynthesis covalently linked around two reaction chambers. Although the three-dimensional architecture of FAS has been mostly defined, it is unclear how reaction intermediates can transfer between distant catalytic domains. Using single-particle EM, we have identified a near continuum of conformations consistent with a remarkable flexibility of FAS. The distribution of conformations was influenced by the presence of substrates and altered by different catalytic mutations, suggesting a direct correlation between conformation and specific enzymatic activities. We interpreted three-dimensional reconstructions by docking high-resolution structures of individual domains, and they show that the substrate-loading and condensation domains dramatically swing and swivel to access substrates within either reaction chamber. Concomitant rearrangement of the beta-carbon-processing domains synchronizes acyl chain reduction in one chamber with acyl chain elongation in the other.

Phylogenetic analysis of the triterpene cyclase protein family in prokaryotes and eukaryotes suggests bidirectional lateral gene transfer

Functional constraints to modifications in triterpene cyclase amino acid sequences make them good candidates for evolutionary studies on the phylogenetic relatedness of these enzymes in prokaryotes as well as in eukaryotes. In this study, we used a set of identified triterpene cyclases, a group of mainly bacterial squalene cyclases and a group of predominantly eukaryotic oxidosqualene cyclases, as seed sequences to identify 5288 putative triterpene cyclase homologues in publicly available databases. The Cluster Analysis of Sequences software was used to detect groups of sequences with increased pairwise sequence similarity. The sequences fall into two main clusters, a bacterial and a eukaryotic. The conserved, informative regions of a multiple sequence alignment of the family were used to construct a neighbour-joining phylogenetic tree using the AsaturA and maximum likelihood phylogenetic tree using the PhyML software. Both analyses showed that most of the triterpene cyclase sequences were similarly grouped to the accepted taxonomic relationships of the organism the sequences originated from, supporting the idea of vertical transfer of cyclase genes from parent to offspring as the main evolutionary driving force in this protein family. However, a small group of sequences from three bacterial species (Stigmatella, Gemmata and Methylococcus) grouped with an otherwise purely eukaryotic cluster of oxidosqualene cyclases, while a small group of sequences from seven fungal species and a sequence from the fern Adiantum grouped consistently with a cluster of otherwise purely bacterial squalene cyclases. This suggests that lateral gene transfer may have taken place, entailing a transfer of oxidosqualene cyclases from eukaryotes to bacteria and a transfer of squalene cyclase from bacteria to an ancestor of the group of Pezizomycotina fungi.

Non-LTR retrotransposons encode noncanonical RRM domains in their first open reading frame

Non-LTR retrotransposons (NLRs) are a unique class of mobile genetic elements that have significant impact on the evolution of eukaryotic genomes. However, the molecular details and functions of their encoded proteins, in particular of the accessory ORF1p proteins, are poorly understood. Here, we identify noncanonical RNA-recognition-motifs (RRMs) in several phylogenetically unrelated NLR ORF1p proteins. This provides an explanation for their RNA-binding properties and clearly shows that they are not related to the retroviral nucleocapsid protein Gag, despite the frequent presence of CCHC zinc knuckles. In particular, we characterize the ORF1p protein of the human long interspersed nuclear element 1 (LINE-1 or L1). We show that L1ORF1p is a multidomain protein, consisting of a coiled coil (cc), RRM, and C-terminal domain (CTD). Most importantly, we solved the crystal structure of the RRM domain, which is characterized by extended loops stabilized by unique salt bridges. Furthermore, we demonstrate that L1ORF1p trimerizes via its N-terminal cc domain, and we suggest that this property is functionally important for all homologues. The formation of distinct complexes with single-stranded nucleic acids requires the presence of the RRM and CTD domains on the same polypeptide chain as well as their close cooperation. Finally, the phylogenetic analysis of mammalian L1ORF1p shows an ancient origin of the RRM domain and supports a modular evolution of NLRs.

Identification, characterization and essentiality of the unusual peroxin 13 from Trypanosoma brucei

Peroxin 13 (PEX13) is one of the components of a peroxisomal membrane complex involved in import of proteins into the matrix of the organelles and has previously been characterized in a variety of organisms. Trypanosomatids (Trypanosoma, Leishmania), protozoan parasites having peroxisome-like organelles designated glycosomes, possess an unusual PEX13 which shares very low sequence identity with others and lacks some typical PEX13 characteristics. It was identified in the databases through its multiple YGx motifs present in a glycine-rich N-terminal region of low sequence complexity. Like other PEX13s, it contains predicted transmembrane segments and a SH3 domain in its C-terminal half. The localization of T. brucei PEX13 in the glycosomal membrane was confirmed by expression of a fusion construct with Green Fluorescent Protein, and western blot analysis of purified organelles and membranes. The C-terminal half of the protein was shown to interact with the third of three pentapeptide repeats of the previously characterized PEX5, the receptor of glycosomal proteins with a type 1 peroxisome-targeting signal, and with PEX14, another component of the same peroxisomal protein import complex in the membrane. PEX13 is essential for the parasite; depletion by RNA interference results in mislocalization of glycosomal proteins and death of the parasites.

Proteomic tools for the analysis of cytoskeleton proteins

Proteomic tools have become an essential part of the tool kit of the molecular biologist, and provide techniques for detecting homologous sequences, recognizing functional domains, modeling, and analyzing the three-dimensional structure for any given protein sequence. Although a wealth of structural and functional information is available for a large number of members of the various classes of cytoskeletal proteins, many more members remain uncharacterized. These computational tools that are freely and easily accessible to the scientific community provide an excellent starting point to predict the structural and functional properties of such partially or fully uncharacterized protein sequences, and can lead to elegantly designed experiments to probe the hypothesized function. This chapter discusses various proteomic analysis tools with a focus on protein structure and function predictions.

The tumor marker human placental protein 11 is an endoribonuclease

Human PP11 (placental protein 11) was previously described as a serine protease specifically expressed in the syncytiotrophoblast and in numerous tumor tissues. Several PP11-like proteins were annotated in distantly related organisms, such as worms and mammals, suggesting their involvement in evolutionarily conserved processes. Based on sequence similarity, human PP11 was included in a protein family whose characterized members are XendoU, a Xenopus laevis endoribonuclease involved in small nucleolar RNA processing, and Nsp15, an endoribonuclease essential for coronavirus replication. Here we show that the bacterially expressed human PP11 displays RNA binding capability and cleaves single stranded RNA in a Mn(2+)-dependent manner at uridylates, to produce molecules with 2',3'-cyclic phosphate ends. These features, together with structural and mutagenesis analyses, which identified the potential active site residues, reveal striking parallels to the amphibian XendoU and assign a ribonuclease function to PP11. This newly discovered enzymatic activity places PP11-like proteins in a completely new perspective.

PDBalert: automatic, recurrent remote homology tracking and protein structure prediction

BACKGROUND

During the last years, methods for remote homology detection have grown more and more sensitive and reliable. Automatic structure prediction servers relying on these methods can generate useful 3D models even below 20% sequence identity between the protein of interest and the known structure (template). When no homologs can be found in the protein structure database (PDB), the user would need to rerun the same search at regular intervals in order to make timely use of a template once it becomes available.

RESULTS

PDBalert is a web-based automatic system that sends an email alert as soon as a structure with homology to a protein in the user's watch list is released to the PDB database or appears among the sequences on hold. The mail contains links to the search results and to an automatically generated 3D homology model. The sequence search is performed with the same software as used by the very sensitive and reliable remote homology detection server HHpred, which is based on pairwise comparison of Hidden Markov models.

CONCLUSION

PDBalert will accelerate the information flow from the PDB database to all those who can profit from the newly released protein structures for predicting the 3D structure or function of their proteins of interest.

Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily, but exhibits an unusual active site

BACKGROUND

Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds.

RESULTS

R.Hpy188I is a Type II REase with unknown structure. PSI-BLAST searches of the non-redundant protein sequence database reveal only 1 homolog (R.HpyF17I, with nearly identical amino acid sequence and the same DNA sequence specificity). Standard application of state-of-the-art protein fold-recognition methods failed to predict the relationship of R.Hpy188I to proteins with known structure or to other protein families. In order to increase the amount of evolutionary information in the multiple sequence alignment, we have expanded our sequence database searches to include sequences from metagenomics projects. This search resulted in identification of 23 further members of R.Hpy188I family, both from metagenomics and the non-redundant database. Moreover, fold-recognition analysis of the extended R.Hpy188I family revealed its relationship to the GIY-YIG domain and allowed for computational modeling of the R.Hpy188I structure. Analysis of the R.Hpy188I model in the light of sequence conservation among its homologs revealed an unusual variant of the active site, in which the typical Tyr residue of the YIG half-motif had been substituted by a Lys residue. Moreover, some of its homologs have the otherwise invariant Arg residue in a non-homologous position in sequence that nonetheless allows for spatial conservation of the guanidino group potentially involved in phosphate binding.

CONCLUSION

The present study eliminates a significant "white spot" on the structural map of REases. It also provides important insight into sequence-structure-function relationships in the GIY-YIG nuclease superfamily. Our results reveal that in the case of proteins with no or few detectable homologs in the standard "non-redundant" database, it is useful to expand this database by adding the metagenomic sequences, which may provide evolutionary linkage to detect more remote homologs.

Phospholipid scramblases and Tubby-like proteins belong to a new superfamily of membrane tethered transcription factors

Motivation: Phospholipid scramblases (PLSCRs) constitute a family of cytoplasmic membrane-associated proteins that were identified based upon their capacity to mediate a Ca²⁺-dependent bidirectional movement of phospholipids across membrane bilayers, thereby collapsing the normally asymmetric distribution of such lipids in cell membranes. The exact function and mechanism(s) of these proteins nevertheless remains obscure: data from several laboratories now suggest that in addition to their putative role in mediating transbilayer flip/flop of membrane lipids, the PLSCRs may also function to regulate diverse processes including signaling, apoptosis, cell proliferation and transcription. A major impediment to deducing the molecular details underlying the seemingly disparate biology of these proteins is the current absence of any representative molecular structures to provide guidance to the experimental investigation of their function.

Results: Here, we show that the enigmatic PLSCR family of proteins is directly related to another family of cellular proteins with a known structure. The Arabidopsis protein At5g01750 from the DUF567 family was solved by X-ray crystallography and provides the first structural model for this family. This model identifies that the presumed C-terminal transmembrane helix is buried within the core of the PLSCR structure, suggesting that palmitoylation may represent the principal membrane anchorage for these proteins. The fold of the PLSCR family is also shared by Tubby-like proteins. A search of the PDB with the HHpred server suggests a common evolutionary ancestry. Common functional features also suggest that tubby and PLSCR share a functional origin as membrane tethered transcription factors with capacity to modulate phosphoinositide-based signaling.

Contact:agb@sanger.ac.uk

SHOC1, an XPF endonuclease-related protein, is essential for the formation of class I meiotic crossovers

Crossovers (COs) are essential for the completion of meiosis in most species and lead to new allelic combinations in gametes. Two pathways of meiotic crossover formation have been distinguished. Class I COs, which are the major class of CO in budding yeast, mammals, Caenorhabditis elegans, and Arabidopsis, depend on a group of proteins called ZMM and rely on specific DNA structure intermediates that are processed to form COs. We identified a novel gene, SHOC1, involved in meiosis in Arabidopsis. Shoc1 mutants showed a striking reduction in the number of COs produced, a similar phenotype to the previously described Arabidopsis zmm mutants. The early steps of recombination, revealed by DMC1 foci, and completion of synapsis are not affected in shoc1 mutants. Double mutant analysis showed that SHOC1 acts in the same pathway as AtMSH5, a conserved member of the ZMM group. SHOC1 is thus a novel gene required for class I CO formation in Arabidopsis. Sequence similarity studies detected putative SHOC1 homologs in a large range of eukaryotes including human. SHOC1 appears to be related to the XPF endonuclease protein family, which suggests that it is directly involved in the maturation of DNA intermediates that lead to COs.

Unraveling the biochemistry and provenance of pupylation: a prokaryotic analog of ubiquitination

Recently Mycobacterium tuberculosis was shown to possess a novel protein modification, in which a small protein Pup is conjugated to the epsilon-amino groups of lysines in target proteins. Analogous to ubiquitin modification in eukaryotes, this remarkable modification recruits proteins for degradation via archaeal-type proteasomes found in mycobacteria and allied actinobacteria. While a mycobacterial protein named PafA was found to be required for this conjugation reaction, its biochemical mechanism has not been elucidated. Using sensitive sequence profile comparison methods we establish that the PafA family proteins are related to the gamma-glutamyl-cysteine synthetase and glutamine synthetase. Hence, we predict that PafA is the Pup ligase, which catalyzes the ATP-dependent ligation of the terminal gamma-carboxylate of glutamate to lysines, similar to the above enzymes. We further discovered that an ortholog of the eukaryotic PAC2 (e.g. cg2106) is often present in the vicinity of the actinobacterial Pup-proteasome gene neighborhoods and is likely to represent the ancestral proteasomal chaperone. Pup-conjugation is sporadically present outside the actinobacteria in certain lineages, such as verrucomicrobia, nitrospirae, deltaproteobacteria and planctomycetes, and in the latter two lineages it might modify membrane proteins.

REVIEWERS

This article was reviewed by M. Madan Babu and Andrei Osterman.

Mutations in AtPS1 (Arabidopsis thaliana parallel spindle 1) lead to the production of diploid pollen grains

Polyploidy has had a considerable impact on the evolution of many eukaryotes, especially angiosperms. Indeed, most--if not all-angiosperms have experienced at least one round of polyploidy during the course of their evolution, and many important crop plants are current polyploids. The occurrence of 2n gametes (diplogametes) in diploid populations is widely recognised as the major source of polyploid formation. However, limited information is available on the genetic control of diplogamete production. Here, we describe the isolation and characterisation of the first gene, AtPS1 (Arabidopsis thaliana Parallel Spindle 1), implicated in the formation of a high frequency of diplogametes in plants. Atps1 mutants produce diploid male spores, diploid pollen grains, and spontaneous triploid plants in the next generation. Female meiosis is not affected in the mutant. We demonstrated that abnormal spindle orientation at male meiosis II leads to diplogamete formation. Most of the parent's heterozygosity is therefore conserved in the Atps1 diploid gametes, which is a key issue for plant breeding. The AtPS1 protein is conserved throughout the plant kingdom and carries domains suggestive of a regulatory function. The isolation of a gene involved in diplogamete production opens the way for new strategies in plant breeding programmes and progress in evolutionary studies.

Crystal structure of SpoVT, the final modulator of gene expression during spore development in Bacillus subtilis

Endospore formation in Bacillus subtilis is orchestrated by five developmental sigma factors and further modulated by several auxiliary transcription factors. One of these, SpoVT, regulates forespore-specific sigma(G)-dependent genes and plays a key role in the final stages of spore formation. We have determined the crystal structure of the isolated C-terminal domain of SpoVT at 1.5 A by experimental phasing techniques and used this model to solve the structure of the full-length SpoVT at 2.6 A by molecular replacement. SpoVT is a tetramer that shows an overall significant distortion mediated by electrostatic interactions. Two monomers dimerize via the highly charged N-terminal domains to form swapped-hairpin beta-barrels. These asymmetric dimers further tetramerize through the formation of mixed helix bundles between their C-terminal domains, which themselves fold as GAF (cGMP-specific and cGMP-stimulated phosphodiesterases, Anabaena adenylate cyclases, and Escherichia coli FhlA) domains. The combination of a swapped-hairpin beta-barrel with a GAF domain represents a novel domain architecture in transcription factors. The occurrence of SpoVT homologs throughout Bacilli and Clostridia demonstrates the ancestral origin of this factor in sporulation.

Characterization of SNAREs determines the absence of a typical Golgi apparatus in the ancient eukaryote Giardia lamblia

Giardia is a eukaryotic protozoal parasite with unusual characteristics, such as the absence of a morphologically evident Golgi apparatus. Although both constitutive and regulated pathways for protein secretion are evident in Giardia, little is known about the mechanisms involved in vesicular docking and fusion. In higher eukaryotes, soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) of the vesicle-associated membrane protein and syntaxin families play essential roles in these processes. In this work we identified and characterized genes for 17 SNAREs in Giardia to define the minimal set of subcellular organelles present during growth and encystation, in particular the presence or not of a Golgi apparatus. Expression and localization of all Giardia SNAREs demonstrate their presence in distinct subcellular compartments, which may represent the extent of the endomembrane system in eukaryotes. Remarkably, Giardia SNAREs, homologous to Golgi SNAREs from other organisms, do not allow the detection of a typical Golgi apparatus in either proliferating or differentiating trophozoites. However, some features of the Golgi, such as the packaging and sorting function, seem to be performed by the endoplasmic reticulum and/or the nuclear envelope. Moreover, depletion of individual genes demonstrated that several SNAREs are essential for viability, whereas others are dispensable. Thus, Giardia requires a smaller number of SNAREs compared with other eukaryotes to accomplish all of the vesicle trafficking events that are critical for the growth and differentiation of this important human pathogen.

Unassigned MURF1 of kinetoplastids codes for NADH dehydrogenase subunit 2

Background

In a previous study, we conducted a large-scale similarity-free function prediction of mitochondrion-encoded hypothetical proteins, by which the hypothetical gene murf1 (maxicircle unidentified reading frame 1) was assigned as nad2, encoding subunit 2 of NADH dehydrogenase (Complex I of the respiratory chain). This hypothetical gene occurs in the mitochondrial genome of kinetoplastids, a group of unicellular eukaryotes including the causative agents of African sleeping sickness and leishmaniasis. In the present study, we test this assignment by using bioinformatics methods that are highly sensitive in identifying remote homologs and confront the prediction with available biological knowledge.

Results

Comparison of MURF1 profile Hidden Markov Model (HMM) against function-known profile HMMs in Pfam, Panther and TIGR shows that MURF1 is a Complex I protein, but without specifying the exact subunit. Therefore, we constructed profile HMMs for each individual subunit, using all available sequences clustered at various identity thresholds. HMM-HMM comparison of these individual NADH subunits against MURF1 clearly identifies this hypothetical protein as NAD2. Further, we collected the relevant experimental information about kinetoplastids, which provides additional evidence in support of this prediction.

Conclusion

Our in silico analyses provide convincing evidence for MURF1 being a highly divergent member of NAD2.

A channel connecting the mother cell and forespore during bacterial endospore formation

At an early stage during Bacillus subtilis endospore development the bacterium divides asymmetrically to produce two daughter cells. The smaller cell (forespore) differentiates into the endospore, while the larger cell (mother cell) becomes a terminally differentiated cell that nurtures the developing forespore. During development the mother cell engulfs the forespore to produce a protoplast, surrounded by two bilayer membranes, which separate it from the cytoplasm of the mother cell. The activation of sigma(G), which drives late gene expression in the forespore, follows forespore engulfment and requires expression of the spoIIIA locus in the mother cell. One of the spoIIIA-encoded proteins SpoIIIAH is targeted specifically to the membrane surrounding the forespore, through an interaction of its C-terminal extracellular domain with the C-terminal extracellular domain of the forespore membrane protein SpoIIQ. We identified a homologous relationship between the C-terminal domain of SpoIIIAH and the YscJ/FliF protein family, members of which form multimeric rings involved in type III secretion systems and flagella. If SpoIIIAH forms a similar ring structure, it may also form a channel between the mother cell and forespore membranes. To test this hypothesis we developed a compartmentalized biotinylation assay, which we used to show that the C-terminal extracellular domain of SpoIIIAH is accessible to enzymatic modification from the forespore cytoplasm. These and other results lead us to suggest that SpoIIIAH forms part of a channel between the forespore and mother cell that is required for the activation of sigma(G).

An unusual RNA recognition motif acts as a scaffold for multiple proteins in the pre-mRNA retention and splicing complex

The yeast pre-mRNA retention and splicing complex counteracts the escape of unspliced pre-mRNAs from the nucleus and activates splicing of a subset of Mer1p-dependent genes. A homologous complex is present in activated human spliceosomes. In many components of the spliceosome, RNA recognition motifs (RRMs) serve as versatile protein-RNA or protein-protein interaction platforms. Here, we show that in the retention and splicing complex, an atypical RRM of the Snu17p (small nuclear ribonucleoprotein-associated protein 17) subunit acts as a scaffold that organizes the other two constituents, Bud13p (bud site selection 13) and Pml1p (pre-mRNA leakage 1). GST pull-down experiments and size exclusion chromatography revealed that Snu17p constitutes the central platform of the complex, whereas Bud13p and Pml1p do not interact with each other. Fluorimetric structure probing showed the entire Bud13p and the N-terminal third of Pml1p to be natively disordered in isolation. Mutational analysis and tryptophan fluorescence confirmed that a conserved tryptophan-containing motif in the C terminus of Bud13p binds to the core RRM of Snu17p, whereas a different interaction surface encompassing a C-terminal extension of the Snu17p RRM is required to bind an N-terminal peptide of Pml1p. Isothermal titration calorimetry revealed 1:1 interaction stoichiometries, large negative binding entropies, and dissociation constants in the low nanomolar and micromolar ranges for the Snu17p-Bud13p and the Snu17p-Pml1p interactions, respectively. Our results demonstrate that the noncanonical Snu17p RRM concomitantly binds multiple ligand proteins via short, intrinsically unstructured peptide epitopes and thereby acts as a platform that displays functional modules of the ligands, such as a forkhead-associated domain of Pml1p and a conserved polylysine motif of Bud13p.

The Est3 protein associates with yeast telomerase through an OB-fold domain

The Est3 protein is a small regulatory subunit of yeast telomerase which is dispensable for enzyme catalysis but essential for telomere replication in vivo. Using structure prediction combined with in vivo characterization, we show here that Est3 consists of a predicted OB (oligo-saccharide/oligo-nucleotide binding) fold. Mutagenesis of predicted surface residues was used to generate a functional map of one surface of Est3, which identified a site that mediates association with the telomerase complex. Surprisingly, the predicted OB-fold of Est3 is structurally similar to the OB-fold of the mammalian TPP1 protein, despite the fact that Est3 and TPP1, as components of telomerase and a telomere capping complex, respectively, perform functionally distinct tasks at chromosome ends. The analysis performed on Est3 may be instructive in generating comparable missense mutations on the surface of the OB-fold domain of TPP1.

AMIN domains have a predicted role in localization of diverse periplasmic protein complexes

We describe AMIN (Amidase N-terminal domain), a novel protein domain found specifically in bacterial periplasmic proteins. AMIN domains are widely distributed among peptidoglycan hydrolases and transporter protein families. Based on experimental data, contextual information and phyletic profiles, we suggest that AMIN domains mediate the targeting of periplasmic or extracellular proteins to specific regions of the bacterial envelope.

Contact:fgueiros@iq.usp.br

Supplementary information: Supplementary data are available at Bioinformatics online.

Remote homology detection of integral membrane proteins using conserved sequence features

Compared with globular proteins, transmembrane proteins are surrounded by a more intricate environment and, consequently, amino acid composition varies between the different compartments. Existing algorithms for homology detection are generally developed with globular proteins in mind and may not be optimal to detect distant homology between transmembrane proteins. Here, we introduce a new profile-profile based alignment method for remote homology detection of transmembrane proteins in a hidden Markov model framework that takes advantage of the sequence constraints placed by the hydrophobic interior of the membrane. We expect that, for distant membrane protein homologs, even if the sequences have diverged too far to be recognized, the hydrophobicity pattern and the transmembrane topology are better conserved. By using this information in parallel with sequence information, we show that both sensitivity and specificity can be substantially improved for remote homology detection in two independent test sets. In addition, we show that alignment quality can be improved for the most distant homologs in a public dataset of membrane protein structures. Applying the method to the Pfam domain database, we are able to suggest new putative evolutionary relationships for a few relatively uncharacterized protein domain families, of which several are confirmed by other methods. The method is called Searcher for Homology Relationships of Integral Membrane Proteins (SHRIMP) and is available for download at http://www.sbc.su.se/shrimp/.

A novel pathway of intercellular signalling in Bacillus subtilis involves a protein with similarity to a component of type III secretion channels

During spore formation in Bacillus subtilis, sigma(E)-directed gene expression in the mother-cell compartment of the sporangium triggers the activation of sigma(G) in the forespore by a pathway of intercellular signalling that is composed of multiple proteins of unknown function. Here, we confirm that the vegetative protein SpoIIIJ, the forespore protein SpoIIQ and eight membrane proteins (SpoIIIAA through SpoIIIAH) produced in the mother cell under the control of sigma(E) are ordinarily required for intercellular signalling. In contrast, an anti-sigma(G) factor previously implicated in the pathway is shown to be dispensable. We also present evidence suggesting that SpoIIIJ is a membrane protein translocase that facilitates the insertion of SpoIIIAE into the membrane. In addition, we report the isolation of a mutation that partially bypasses the requirement for SpoIIIJ and for SpoIIIAA through SpoIIIAG, but not for SpoIIIAH or SpoIIQ, in the activation of sigma(G). We therefore propose that under certain genetic conditions, SpoIIIAH and SpoIIQ can constitute a minimal pathway for the activation of sigma(G). Finally, based on the similarity of SpoIIIAH to a component of type III secretion systems, we speculate that signalling is mediated by a channel that links the mother cell to the forespore.