CREST - a large and diverse superfamily of putative transmembrane hydrolases
© Pei et al; licensee BioMed Central Ltd. 2011
Received: 12 April 2011
Accepted: 6 July 2011
Published: 6 July 2011
A number of membrane-spanning proteins possess enzymatic activity and catalyze important reactions involving proteins, lipids or other substrates located within or near lipid bilayers. Alkaline ceramidases are seven-transmembrane proteins that hydrolyze the amide bond in ceramide to form sphingosine. Recently, a group of putative transmembrane receptors called progestin and adipoQ receptors (PAQRs) were found to be distantly related to alkaline ceramidases, raising the possibility that they may also function as membrane enzymes.
Using sensitive similarity search methods, we identified statistically significant sequence similarities among several transmembrane protein families including alkaline ceramidases and PAQRs. They were unified into a large and diverse superfamily of putative membrane-bound hydrolases called CREST (alkaline ceramidase, PAQR receptor, Per1, SID-1 and TMEM8). The CREST superfamily embraces a plethora of cellular functions and biochemical activities, including putative lipid-modifying enzymes such as ceramidases and the Per1 family of putative phospholipases involved in lipid remodeling of GPI-anchored proteins, putative hormone receptors, bacterial hemolysins, the TMEM8 family of putative tumor suppressors, and the SID-1 family of putative double-stranded RNA transporters involved in RNA interference. Extensive similarity searches and clustering analysis also revealed several groups of proteins with unknown function in the CREST superfamily. Members of the CREST superfamily share seven predicted core transmembrane segments with several conserved sequence motifs.
Universal conservation of a set of histidine and aspartate residues across all groups in the CREST superfamily, coupled with independent discoveries of hydrolase activities in alkaline ceramidases and the Per1 family as well as results from previous mutational studies of Per1, suggests that the majority of CREST members are metal-dependent hydrolases.
This article was reviewed by Kira S. Markarova, Igor B. Zhulin and Rob Knight.
Membrane proteins play important roles in various biological processes such as transportation and signal transduction across the membrane. Although a significant fraction of proteins in genomes are predicted to be integral membrane proteins , our understandings of their function, mechanism and evolution remain limited due to experimental difficulties in assaying their activities and obtaining their structures. The low complexity of the hydrophobic membrane-spanning segments also renders transmembrane proteins difficult targets for computational sequence analysis and structure prediction .
Membrane-bound enzymes catalyze a variety of key reactions involving proteins, lipids or other substrates within or near lipid bilayers. A limited number of membrane-bound enzymes with known structures exist, including several intramembrane proteases  and members of the MAPEG (membrane-associated proteins in eicosanoid and glutathione metabolism) family . A family of alkaline ceramidases, which are seven-transmembrane proteins without known structures, catalyze the hydrolysis of the amide bond in ceramide to form sphingosine . Sphingosine, ceramide and their phosphorylated products are bioactive sphingolipid molecules involved in a number of biological processes . Recently, a large group of putative transmembrane hormone receptors called progestin and adipoQ receptors (PAQRs) were found to be distantly related to alkaline ceramidases, raising the possibility that they may also function as membrane enzymes [7, 8]. Interestingly, the bacterial homologs of these receptors have hemolysin activities .
In this study, we employed sensitive sequence comparison methods to search for homologs of alkaline ceramidases and PAQRs. Three additional known protein families were identified: the Per1 family of fatty acid remodeling hydrolases for GPI-anchored proteins [10, 11], the SID-1 family of putative RNA transporters involved in systematic RNA interference [12, 13], and the TMEM8 family of putative tumor suppressors [14, 15]. In addition, we identified distant homologs comprising five distinct groups of proteins with unknown function. We unified these proteins into a superfamily of putative transmembrane hydrolases called CREST (alkaline ceramidase, PAQR receptor, Per1, SID-1 and TMEM8). They share a core structure of seven predicted transmembrane segments and five conserved residues (three histidines, one aspartic acid and one serine). Sequence conservation and mutational studies in the Per1 family  suggest that most members of this superfamily are metal-dependent hydrolases. Our analyses and predictions offer insights into functions and mechanisms for members in this diverse superfamily.
Sequence similarity searches and homology inference for the CREST superfamily
Extensive transitive PSI-BLAST  searches (see Methods for details) uncovered statistically significant sequence similarities between alkaline ceramidases, PAQR receptors, bacterial hemolysins III, the SID-1 family, and the TMEM8 family. For example, a PSI-BLAST search starting with a human alkaline ceramidase ACER3 (NCBI gene identification (gi) number: 296439452, with an e-value inclusion cutoff of 1e-4) found a bacterial hypothetical protein (gi: 114706410 from Fulvimarina pelagi) with an e-value of 6e-7 in the fourth iteration. The PSI-BLAST search using this bacterial protein identified statistically significant similarities to PAQR receptors and various bacterial hemolysins III (e.g. a human PAQR receptor (gi: 38018661) was found in the ninth iteration with an e-value of 1e-6 and a sequence identity of 13%). Similarly, the homologous relationship between alkaline ceramidases and the SID-1 family was supported by statistically significant hits through intermediate protein sequences (e.g. a hypothetical protein from Nitrosomonas eutropha (gi: 114331832) was found by human ACER3 with an e-value of 2e-6 in the seventh iteration, and it in turn found a sea urchin SID-1 protein (gi: 115686293) with an e-value of 4e-06 in the third iteration). A PSI-BLAST search starting with a bacterial hemolysin III (gi: 1708219 from Bacillus cereus) found a TMEM8 protein (gi: 47227992 from Tetraodon nigroviridis) with an e-value of 9e-6 and a sequence identity of 12% at the fourth iteration. Conversely, transitive PSI-BLAST searches starting from TMEM8 proteins also found PAQRs and bacterial hemolysins III with statistically significant e-values (less than 1e-4).
Additional support of homology among these families was obtained by online CSI-BLAST  searches (e-value inclusion cutoff 1e-4, against nr database). For example, a CSI-BLAST search starting from the human alkaline ceramidase ACER3 found PAQR, TMEM8, SID-1 and Per1 family members with statistically significant e-values (less than 1e-4).
Sequence groups of the CREST superfamily
Three types of ceramidases exist in mammalian genomes and are classified according the the pH optima of their catalytic activity . Acid ceramidase resides in the lysosome and its catalytic domain exhibits similarity to members of the Ntn (N-terminal nucleophile)-hydrolases . Neutral ceramidase is a metal-dependent hydrolase with a different structural fold . The three copies of alkaline ceramidases in mammalian genomes are transmembrane proteins located in the ER and Golgi that have been indicated in intracellular sphingolipid regulation [5, 22]. Close homologs of mammalian alkaline ceramidases form a tight group in the CLANS clustering result (Figure 2). Alkaline ceramidases are populated in other eukaryotic organisms, such as insects, fungi and plants. Interestingly, alkaline ceramidases from both fungi  and insects  have been shown to have a negative effect on life span, suggesting that they are potential targets in ageing studies.
Two distinct groups of PAQR receptors
Several members of the PAQR receptor family  have been identified as putative receptors for progestin [26, 27] and the adipose-derived protein adiponectin (also called adipoQ) . Adiponectin and its receptors are associated with a number of pathophysiological conditions including obesity and type 2 diabetes [29, 30]. The CLANS sequence clustering revealed two distinct groups in the large PAQR family (Figure 2) comprising nearly 2,000 eukaryotic and bacterial proteins.
The first PAQR group consists of only eukaryotic proteins and includes the originally identified adiponectin receptors in mammals  and the putative membrane-bound progestin receptors in fish  and their mammalian orthologs  (their subcellular localization and function as progestin receptors have been controversial and under debate ). This group also includes members from fungi , plants  and various protists. An expansion of this group occurs in many eukaryotic species, e.g., the human genome has nine members and the budding yeast Saccharomyces cerevisiae has four members (Izh1-4p, Figure 1), suggesting that these putative receptors can potentially accept a variety of ligands. The four S. cerevisiae members have been shown to play a role in zinc metabolism . One of them (Izh2p) was also shown to be involved in lipid and phosphate metabolism  and to inhibit the expression of a gene involved in iron uptake .
The second PAQR group consists of bacterial and eukaryotic proteins. The bacterial homologs are frequently annotated as 'hemolysin III', suggesting the role of a cytotoxin. A bacterial member from Bacillus cereus was shown to be a pore-forming hemolysin . The same activity was demonstrated for a homolog from the Gram-negative pathogen Vibrio vulnificus . Members of this group are widely distributed in various Gram-positive and Gram-negative bacterial phyla (see Additional file 1). The STRING protein-protein interaction server  revealed that some of the bacterial hemolysin III genes are neighbors to a gene encoding a hypothetical protein that contains a DegV domain (Pfam entry: PF02645) known to bind fatty acids (e.g., gi: 81428574 from Lactobacillus sakei found a DegV domain-containing protein LSA0964 with a score of 0.752). For some other bacterial members, the STRING server reported the gene neighborhood association with a DUF1836 domain-containing protein (e.g. gi: 157151576 from Streptococcus gordonii found a DUF1836 domain-containing hypothetical protein SGO_1308 with a score of 0.877). DUF1836 is classified in the helix-turn-helix (HTH) clan in the Pfam database and could serve as a transcriptional regulator for some of the bacterial hemolysins III. The two human proteins from the second group (PAQR10 and PAQR11) are annotated as 'monocyte to macrophage differentiation factors'  with limited experimental studies . Compared to the first PAQR group, the second group has a more restricted distribution in eukaryotes. This group has no plant members and the only two species from fungi in this group belong to the parasitic Microsporidia phylum.
The Per1 family
The Per1p protein in S. cerevisiae was originally identified as a suppressor of the cdc1 mutant  and was suggested to play a role in manganese homeostasis . Later studies showed that Per1p and its mammalian ortholog PGAP3 are involved in the lipid remodeling of GPI-anchored proteins [10, 11]. This remodeling process involves the detachment of lipid chains from the intermediate GPI-anchored proteins and the reattachment of new, usually saturated lipid chains, thus facilitating the localization of GPI-anchored proteins to specialized microdomains ('lipid rafts'). Both the yeast Per1p and the mammalian PGAP3 were shown to be required for the hydrolysis of lipid moieties from GPI-anchored proteins. Direct evidence of hydrolase activity in Per1p or PGAP3 is lacking. However, their homology relationship to alkaline ceramidase (which catalyzes a similar hydrolysis reaction) and the conservation of similar residues (discussed below) together suggest that Per1p and PGAP3 are hydrolases. The Per1 family (Figure 2) comprises about 150 closely related eukaryotic proteins mainly from fungi, plants and metazoans.
The SID-1 family
The SID-1 protein from Caenorhabditis elegans was characterized as a transmembrane protein that facilitates double-stranded RNA (dsRNA) uptake and thus plays a key role in the systematic RNAi response [12, 13]. Its mammalian homologs were also implicated in the RNAi pathways [42, 43]. SID-1 members have a restricted distribution in eukaryotes, as they are only found in metazoan species and a few lower organisms such as Dictyostelium discoideum and Monosiga brevicollis. Mammalian genomes contain two closely related SID-1 homologs, while copy number varies in insects (e.g., no copies in Drosophila melanogaster but up to three copies in Tribolium castaneum) [44, 45]. SID-1 proteins are usually long with more than 800 amino acids and have four additional predicted membrane segments in addition to the seven core membrane segments (discussed below).
Lineage-specific expansion of the SID-1 family is observed in C. elegans and Caenorhabditis briggsae, with several divergent copies of SID-1-like proteins in both genomes. In C. elegans, only the originally identified SID-1 protein and another SID-1-like protein tag-130 possess the complete set of conserved residues (Figure 1), while the other SID-1-like proteins could have lost the proposed hydrolase activity. The SID-1-like protein tag-130 in C. elegans exhibits higher sequence similarity to SID-1-like proteins in mammals and insects than the originally identified C. elegans SID-1 protein. Experimental studies have shown that unlike SID-1, tag-130 does not contribute to systematic RNAi response in C. elegans  and its function remains to be investigated.
The TMEM8 family
The TMEM8 family proteins share a Pfam domain of unknown function (DUF3522) corresponding to the transmembrane region that is homologous to other CREST members. They are present in metazoans and plants. Vertebrate genomes have three copies of TMEM8 proteins (TMEM8a-c) while most insects only have one copy. TMEM8a and TMEM8b in vertebrates, as well as the TMEM8 family members in other metazoan species and plants, are multi-domain proteins usually with more than 400 amino acids. The transmembrane domains are located at their C-termini. They also possess a divergent EGF domain with six conserved cysteines (detected by HHpred) and an unannotated N-terminal region. TMEM8c proteins in vertebrates are single-domain proteins of about 220 amino acids. The conserved serine and aspartic acid are absent in TMEM8c proteins (Figure 1), possibly resulting in the loss of hydrolase activity. Recently, several studies have indicated that human TMEM8b (also named NGX6 for nasopharyngeal carcinoma associated gene 6) was down-regulated in several cancer cell lines and may be a tumor suppressor [14, 15].
Five groups of proteins with unknown function
Hypothetical group 1 consists of proteins from bacteria and plants. One member in this group is SAG18 (gi: 15223970) from Arabidopsis thaliana, which was found in a screen of senescence-associated genes . The bacterial members of this group are mainly from Proteobacteria (see Additional file 1). Hypothetical group 2 consists of only bacterial proteins. They are mainly from the Alphaproteobacteria class (see Additional file 1). These two hypothetical groups appear to be more closely related to alkaline ceramidases (with a modest bootstrap percentage value of 82% for grouping them with alkaline ceramidases) than to the other known groups according to the phylogenetic tree (Figure 3). Therefore, it would be interesting to experimentally test the ceramidase activity for these hypothetical proteins. Proteins of hypothetical group 3 are from various protists of genera such as Trypanosoma, Toxoplasma and Paramecium, as well as from the Choanoflagellate Monosiga brevicollis. They form a well-supported clade with the TMEM8 group (Figure 3). Hypothetical groups 4 and 5 consist of proteins from green plants (Viridiplantae) and green algae (Chlorophyta) respectively. Hypothetical group 5 and the two PAQR groups form a well-supported clade, whereas the relationship of hypothetical group 4 to the other groups remains unclear (it clusters with the Per1 group with a low bootstrap support value) (Figure 3).
The phylogenetic tree revealed two separate locations of bacterial sequences (blue names in Figure 3). One cluster of bacterial sequences (from hypothetical groups 1 and 2) are clustered with alkaline ceramidases from eukaryotes, and the other cluster (bacterial hemolysins III) is grouped with PAQRs from eukaryotes. Such a distribution suggests that at least two copies of CREST members are present in the last common ancestor of bacteria and eukaryotes. No CREST members were found from archaeal organisms, likely reflecting the difference of membrane lipid composition of archaea from bacteria and eukaryotes . Alkaline ceramidase and Per1 are hydrolases acting on the amide or ester bond between the hydrophobic acyl group and the hydrophilic head group (sphingosine or glycerol). However, ether bond exists between the acyl group and the head group of archaeal membrane lipids, in contrast to the ester or amide bond in bacterial and eukaryotic membrane lipids.
Domain architecture and sequence motifs of the CREST superfamily
Three conserved sequence motifs with semi-invariant residues were identified in all CREST groups (Figure 1). The first motif SxxxH ('x' is any amino acid) is located at the end of the second predicted core transmembrane segment. The second motif with a conserved aspartic acid is located at the beginning of the third predicted core transmembrane segment. The third motif with two conserved histidines (HxxxH) resides at the beginning of the seventh predicted core transmembrane segment. We also observed a conserved position comprising mainly small residues in the sixth predicted core transmembrane segment (see Additional file 2). Topologically, the three motifs are predicted to be located near the same side of the membrane and thus may be proximal to each other and form the putative active site (Figure 4a). The three predicted core transmembrane segments (TM2, TM3 and TM7) harboring the motifs could also be spatially close to each other (shown in Figure 4a). Coincidently, in the seven-transmembrane GPCR structures , transmembrane segments 2, 3 and 7 are also spatially close and interact with each other. We did not detect statistically significant sequence similarities between CREST members and GPCRs. The set of CREST-specific sequence motifs was not found in GPCRs, and GPCRs have a different topology than PAQRs in the plasma membrane (discussed below). Structure determination of CREST members could help elucidate whether CREST members exhibit a similar structural fold as GPCRs and whether they are evolutionarily related to GPCRs.
Cataloguing distantly related protein families into superfamilies helps understanding their evolution and deducing useful information about their structure and function. Such a level of classification is used in protein structure classification databases such as SCOP  and CATH  (the superfamily level), the Pfam protein families database  (the clan level), the MEROPS peptidase database  (the clan level), and in membrane transporter classification  (the superfamily level). Distantly related protein families can result from gene duplication events followed by sequence and functional divergence, or from accelerated evolution in certain evolutionary lineages. Although they can have different cellular functions and/or phylogenetic distributions, distantly related protein families often share similar biochemical activities, exemplified by enzymes performing similar reactions. The identities and placements of active site residues responsible for their common biochemical activities are usually conserved as well. Cross-referencing information of distantly related protein families can shed light on their function and help experimental design.
In this study, we unified five known protein families as well as several groups of proteins with unknown function into a large and diverse superfamily of putative transmembrane hydrolases comprising nearly 3000 sequences in the current sequence database. The statistically significant sequence similarities detected by multiple similarity search methods, together with shared sequence motifs, support the homologous relationships among members of the CREST superfamily. The most striking feature of the CREST superfamily of transmembrane proteins is the large functional capacity revealed by independent experimental studies of its members, with diverse and seemingly unrelated functions of hormone receptors (PAQRs), bacterial hemolysins, dsRNA transporters (SID-1) and putative tumor suppressers (TMEM8). It also remains unclear how the putative hydrolase activity contributes to these functions. The hydrolase prediction could open new directions for the future research of CREST members.
Based on the inferred homologous relationship between PAQRs and alkaline ceramidases, it has been speculated that PAQRs could function as ceramidases and generate sphingosine as a second messenger for downstream signaling . Although no direct evidence of these receptors being hydrolases currently exist, experimental results showed that sphingolipids function as downstream effectors of a PAQR member in yeast , and that a ceramidase inhibitor can antagonize human PAQRs . Furthermore, a latest study showed that adiponectin can induce increased ceramidase activity and this effect is associated with adiponectin receptors . As two of the CREST groups, alkaline ceramidase and Per1, perform similar yet distinct hydrolysis reactions on different substrates, it is also possible that the substrates and reaction products of the other CREST groups differ, which could contribute to their functional diversity. Therefore, the future experimental studies of PAQRs should not be restricted to just the potential effects of ceramidase activity. Other hydrolase activities on various lipid molecules, especially the phospholipase A2 activity (suggested for Per1), and their cellular effects can be tested for PAQR members. The products of the phospholipase A2 activity on a phospholipid substrate are a lysophospholipid molecule and a free fatty acid , both of which are potential precursors of second messengers that could regulate downstream signaling events [59, 60].
Classic seven-transmembrane GPCR receptors transduce signals across the plasma membrane and rely on a separate effector protein (such as adenylate cyclase or phospholipase C) to produce second messengers. Well-studied enzyme-linked receptors, such as receptor tyrosine kinases, receptor guanylate cyclases and histidine kinases, are multi-domain proteins that have soluble domains in the cytosol with enzymatic activities. If PAQRs do prove to be receptors with hydrolase activity, they would represent a unique type of receptors that uses the enzymatic activity in the transmembrane domain to generate second messengers (such as sphingosine resulted from ceramidase activity) or their precursors (such as lysophospholipid molecules resulted from phospholipase A2 activity). Further studies are required to determine if PAQRs possess intrinsic ceramidase or other hydrolase activity. A recent genetic linkage analysis revealed that one single nucleotide polymorphism (rs10920533) in human ADIPOR1 (the gene encoding adiponectin receptor 1) interacts with plasma saturated fatty acids , indirectly supporting its role in lipid modification. Another interesting question is if and how the hydrolase activity can be modulated by the binding of extracellular ligands such as adiponectin. The second group of PAQRs includes bacterial hemolysins. Besides pore forming , these toxins could exert their virulence through the putative hydrolase activity, e.g., by changing membrane lipid composition or generating second messengers to modify signaling events in host cells.
In contrast to GPCR receptors, human adiponectin receptors have an internal N-terminus and an extracellular C-terminus , thus positioning the conserved, putative active site residues near the cytosol side of the plasma membrane (Figure 5). The other human PAQR members were also predicted to have their N-termini residing in the cytosol (predictions made by Phobius ). The C. elegans SID-1, also located in the plasma membrane, adopts a topology similar to adiponectin receptors with its N-terminus in the cytosol and its putative active site near the inner lipid layer of the plasma membrane  (Figure 5). SID-1 was proposed to be a dsRNA transporter as its transfection to Drosophila (lacking an endogenic SID-1 ortholog) cells is sufficient to enable dsRNA uptake [13, 62]. In this case, SID-1 has the pore-forming ability similar to bacterial hemolysins III. In a less likely scenario, SID-1 itself is not a dsRNA transporter, but indirectly causes dsRNA uptake by regulating the activity of a separate dsRNA transporter . The proposed hydrolase activity of SID-1 could be involved in regulating its dsRNA transport ability, whether such a transport process is directly or indirectly mediated by SID-1. It would also be interesting to test the hydrolase activity and its roles in the TMEM8 family of putative tumor suppressors.
The specific roles of the conserved residues in CREST proteins have not been reported, and the catalytic mechanisms of alkaline ceramidases and Per1 are yet to be determined. Mutational studies performed on yeast Per1p  with the proposed phospholipase A2 activity offer valuable information about the relative importance of some conserved residues. Alanine mutation of either the conserved histidine in the first motif ('SxxxH') or the first conserved histidine in the third motif ('HxxxH') abolished Per1p activity, suggesting that these histidines are candidates of active site residues that directly contribute to catalysis. On the other hand, alanine mutation of the conserved serine residue in the first motif (the only serine conserved among all CREST groups), as well as alanine mutations of several other residues that are not conserved among all CREST groups, did not affect the activity of Per1p. Therefore, the proposed phospholipase A2 activity of Per1p is unlikely to have a serine esterase-type mechanism that is utilized in certain groups of soluble phospholipase A2 enzymes with an α/β hydrolase fold . Although the conserved serine in the first motif may not directly involve in catalysis, it could participate in hydrogen-bonding interactions with an active site residue to help maintain its optimal orientation for metal-binding or catalysis. This serine appears to be important for the activity of SID-1, as one of the RNA interference-defective SID-1 mutants has a single mutation of this serine to an isoleucine (see the supplemental Figure S3 in reference ).
Three histidines and one aspartate are conserved in CREST proteins, reminiscent of the active site composition of several well-studied metal-dependent hydrolases such as thermolysin , bovine carboxypeptidase A , D-Ala-D-Ala carboxypeptidase  and the intramembrane S2P-like proteases . These enzymes adopt different structural folds, yet they all utilize histidines and negatively charged residues for their catalytic activity (three residues used for metal binding and one residue used for activating a metal-bound water molecule as the nucleophile to attack the scissile bond). The presence of such a set of conserved residues in the CREST superfamily and the mutational results of yeast Per1p support the hypothesis that most CREST members are metal-dependent hydrolases. A large group of secreted enzymes with phospholipase A2 activity (proposed for Per1p) are also metal-dependent hydrolases with histidines in their active sites .
We inferred distant homologous relationships among five transmembrane protein families (alkaline ceramidase, PAQR, Per1, SID-1 and TMEM8) and several groups of proteins with unknown function. Members of this superfamily share seven predicted core transmembrane segments and a set of conserved histidine and aspartate residues. Such a conservation pattern coupled with experimental evidence suggests that they are putative metal-dependent hydrolases acting on molecules with fatty acid moieties. It remains unclear how such putative hydrolase activity contributes to the diverse and seemingly unrelated functions of hormone receptors, bacterial hemolysins, the SID-1 dsRNA transporters and the TMEM8 putative tumor suppressers. Further experimental investigations could reveal new aspects of lipid biology in the actions of these transmembrane proteins.
Sequence similarity searches
PSI-BLAST  was used to search for homologs of the alkaline ceramidase family starting with the human alkaline ceramidase protein ACER3 (NCBI gene identification (gi) number: 296439452) against the non-redundant (nr) protein database (e-value inclusion cutoff: 1e-4, more stringent than the default value in NCBI PSI-BLAST web server (0.005)). To perform transitive searches, the protein hits found by PSI-BLAST were grouped by BLASTCLUST (with the score coverage threshold (-S, defined as the bit score divided by alignment length) set to 1, length coverage threshold (-L) set to 0.5, and no requirement of length coverage on both sequences (-b F)) and a representative sequence from each group was used to initiate new PSI-BLAST searches. The HHpred web server  was used for profile-against-profile-based similarity searches using several members found by PSI-BLAST against the Pfam database  and the human proteome database with default parameter settings. Newly identified human homologs were again subject to transitive PSI-BLAST searches. Several CREST members were also submitted to the CSI-BLAST server  (database: nr; e-value inclusion cutoff: 1e-4; maximum iteration number: 20) to further verify the findings of PSI-BLAST and HHpred.
Domain architecture analysis
HMMER3  and HHpred were used to detect known Pfam domains (Pfam version: 24.0) in CREST members with default parameter settings. Phobius  was used to predict transmembrane segments and membrane topology. Phobius and SignalP 3.0  were used to predict signal peptides.
Sequence clustering, alignment and phylogenetic reconstruction
For sequence clustering, highly similar sequences were removed at the 95% identity level by CD-HIT . Sequence fragments with less than 100 residues were removed. CLANS  was used to cluster the reduced non-redundant set of domains based on pairwise BLAST p-values. The clustering was run to equilibrium in a 2D representation under default settings. Connections between points (sequences) in the final diagram were set by a BLAST P-value cutoff of 1e-6. For sequence alignment, representative CREST members in each group were manually selected to sample sequences from diverse organisms with emphasis on those that have been experimentally studied and described in this article. The multiple sequence alignment of these representatives was made by PROMALS that uses information of predicted secondary structures and database homologs . The alignment was then improved by manual curation. The MOLPHY package  was used for phylogenetic reconstruction based on this alignment (positions with gap fractions larger than 0.1 were removed). The JTT amino acid substitution model  was used in MOLPHY. The local estimates of bootstrap percentages were obtained by the RELL method  (-R option in the ProtML program of MOLPHY). We also used PhyML (version 3.0)  for phylogenetic reconstruction with default parameter settings (the LG amino acid substitution model , four rate categories of a discrete gamma model to take into account rate variability among sites, and an approximate likelihood-ratio test for branch support ).
Reviewer's report 1
Kira S. Markarova, National Center for Biotechnology Information
This paper describes new superfamily of membrane proteins predicted to be metal-dependent hydrolases. Among this superfamily there are several important eukaryotic families that before were not recognized as being homologous. Generally this is a good paper with all conclusions well-supported and justified. In my opinion, however, the important message of the paper could be fit in a shorter format, like, for example, the "Discovery note" in Biology Direct. For instance, the Results section of the paper describing the support for unification of these families is an "overkill". It would be enough to briefly describe the results of either PSI-BLAST or HHpred and mention that other programs support them. The description of the others work related to these families also could be briefer, since not everything what was shown is related to this work. There are other areas for omitting in family description and discussion, just because that they reiterate the same information that is already present in the text somewhere.
We considered other article types when preparing the manuscript, but decided that the length (less than 1500 words) and figure number (no more than one figure or table) limits of the discovery note format could not accommodate the complete contents of this work. In this revised version, one more figure is added to make a total of five figures. We feel that the new figure showing the phylogenetic tree is necessary since it provided additional information about the relationships among the CREST groups and helped address comments raised by all three reviewers. Researchers can have different opinions on the necessity of detailed description of similarity search results. For example, the third reviewer asked for more support from other similarity search methods. We feel that providing details of the similarity search results makes the descriptions scientifically more rigorous and would be appealing to audiences who are interested in these details. We agree that some of the cited works appeared to be not closely related to this work. However, we still considered them worthwhile since not much is known for many CREST groups and the limited existing literature could potentially be useful for cross reference when new experimental data are available in the future. We checked the manuscript and deleted some texts that were reiterated.
For people who are interested in these proteins it would be helpful to provide a supplementary material with all ~ 3000 identified protein IDs and other related information like organism where it was found and assignment to a family.
There are a lot of bacterial proteins in CREST superfamily. For most of them no functional information is available. Neighborhood analysis can provide some insights (same "guilt by association approach" as the analysis of domain fusions) into their function and point to the potential functionally related gene families. Have you tried to look into it? If not, it is worth checking and reporting.
The CREST superfamily has two nice motifs: SxxxH and HxxxH. Have you tried to apply motif search in order to identify families with the same motifs but not found by sequence similarity approach? Surprisingly this approach still can be useful (as we recently showed for FtsZ superfamily - PMID: 20459678). In any event it would be interesting to know the outcome of this search.
What is your opinion on the origin of CREST family? Could the ancestor be placed in LUCA and if yes, are there any ideas, on its ancestral function. Have you noticed any representatives of CREST family in archaea?
We built a phylogenetic tree using MOLPHY for the eleven groups of representative sequences shown in Figure 1. According to this tree, there are two separate groups of bacterial proteins associated with eukaryotic alkaline ceramidases and PAQR receptors, respectively. We thus proposed in the revised manuscript that there might be at least two copies of CREST members in the last common ancestor of eukaryotes and bacteria. We did not find members of the CREST superfamily in archaea. We proposed in the revised manuscript that the lack of CREST members in archaea is likely related to the difference of archaeal membrane lipid composition from those of the bacteria and eukaryotes. CREST members such as alkaline ceramidase and Per1 are proposed to hydrolyze the amide or ester bond between the hydrophobic acyl group and the hydrophilic head group (sphingosine or glycerol). However, in archaea, such amide and ester bonds are replaced by the ether bonds in membrane lipids. Therefore, the reactions catalyzed by CREST members may not be useful in archaea.
Reviewer's report 2
Igor B. Zhulin, University of Tennessee at Knoxville
The novelty is rather marginal. A relationship between eukaryotic PAQR proteins and bacterial hemolysin-type proteins has been captured in the Pfam HlyIII family (PF03006), which includes more than 2000 transmembrane proteins. Current work does expand this family into a superfamily by essentially combining it with other Pfam families, such as Ceramidase (PF05875), and DUF3522 (PF12036) and also identifies several novel families within this superfamily. This is definitely a finding, but it does not sound too exciting. Authors postulate that "Cataloguing distantly related protein families into superfamilies helps understanding their evolution and deducing useful information about their structure and function". I agree, but this work does not suggest too much about the evolution of this superfamily and useful information about function is rather limited. It appears to me that the most useful information is the identification of conserved residues that suggest the hydrolase activity - something, which is testable experimentally.
More serious concern is that no actual data is provided other than some examples. I understand that there are more than 3000 protein sequences, but there must be some way of presenting it. Figure 1 shows a representative alignment of less than 60 sequences. As for the rest, we have to take authors' word for it. Providing the actual data (as a supplementary file, I suppose) is especially important in this case, because (i) cataloguing is the whole purpose of this work (but "the full catalog" is missing, we've got only "the advertisement with 'best buy' items"), and (ii) in the absence of a model deposited to a public resource (e.g. Pfam) this is the only way for biologists to obtain this information. I think both of these concerns can be easily addressed by placing more emphasis on novelty and potential impact and by providing data.
In the revised manuscript, we provided supplementary data listing the eleven groups of CREST proteins found by transitive PSI-BLAST searches. Besides the linking of various known protein families and hypothetical groups in the CREST superfamily and the identification of the conserved and potential active site residues, we feel that one novel aspect of this work lies in the discovery of a diverse array of cellular functions that this superfamily encompasses. While the alkaline ceramidase and Per1 were proposed to possess hydrolase activities in ceramide or GPI-anchored protein biosynthesis pathways, it is quite interesting how such putative hydrolase activity contributes to the diverse and seemingly unrelated functions of hormone receptors, bacterial hemolysins, the SID1 dsRNA transporters and the TMEM8 putative tumor suppressers. The hydrolase prediction can potentially provide new directions in the future experimental studies of these proteins. We proposed that PAQRs could represent a novel type of receptors that use their transmembrane domains to generate second messengers. We also suggested that the experimental work on these receptors should not be restricted to testing the ceramidase activity, as CREST members can catalyze different reactions exemplified by alkaline ceramidases and Per1. We emphasized these points in the discussion and conclusions of the revised manuscript.
Reviewer's report 3
Rob Knight, University of Colorado at Boulder
In this manuscript, the authors unite several families of putative transmembrane hydrolases into one superfamily using similarity search methods, specifically iterated PSI-BLAST. The superfamily contains several protein families of unknown function, providing a clear hypothesis concerning the function of these new members (although no experimental data are provided for validation). Overall, the paper is competently written and executed, though not especially groundbreaking. One puzzle is that the authors use iterated PSI-BLAST rather than profile-profile techniques such as COMPASS, which might be worth commenting on as the same group has previously argued that profile-profile methods provide substantially better results than sequence-profile matching, especially for detecting similarities in highly diverged protein families. The technique is also similar to the Shotgun algorithm from Patsy Babbitt's lab and implemented in some of our own software, which the authors might consider citing (although the approach is straightforward and has likely been independently invented many times).
For similarity searches, we used HHpred, a profile-profile-based method similar to COMPASS, but with the predicted secondary structure in its scoring to obtain better sensitivity. Transitive PSI-BLAST searches were still used since (1) they helped linking most of the groups with statistical support; (2) they enabled exhaustive search for members of the CREST superfamily in the NCBI nonredundant sequence database (nr). No profile-profile based searches on the nr database currently exist as far as we know, possibly due to the time-consuming nature of such searches. We submitted some CREST members to the COMPASS web server . Using the alkaline ceramidase ACER3 from human as a query against the Pfam 23.0 database, COMPASS did identify Hemolysin III (e-value: 8.05e-13) as the third best hit and Per1 (e-value: 2.94e-06) among the top 20 hits. However, most of the top 20 hits with statistically significant e-values (less than 1e-4), such as Herpes_LMP1 and Macoilin, appear to be false positives as they do not possess the conserved motifs unique to the CREST superfamily. A COMPASS search using yeast Per1p as the query also retrieved many false positives as significant hits. These results reflect possible issues with the profiles of membrane proteins on which COMPASS has not been extensively tested. We thus did not include the COMPASS results in the manuscripthttp://prodata.swmed.edu/compass.
It would be interesting to compare an explicit phylogenetic tree with the clans output shown in Figure 1. Doug Theobald, now at Brandeis, developed a method for converting COMPASS profile-profile significance scores into a phylogenetic tree (PMID 16266719) which it might be interesting to use here.
As the COMPASS results of several CREST members frequently included false positives among the top hits, we chose not to use the method suggested by the reviewer. We instead added the phylogenetic analysis of CREST members by MOLPHY, a frequently used tree reconstruction method based on maximum likelihood. The tree reported by MOLPHY is consistent with the CLANS grouping and provided more information about the positioning of the hypothetical groups relative to the known groups. We added the phylogenetic analysis results in the revised manuscript. We also provided the phylogenetic tree reconstructed by another program PhyML in the supplementary data.
The discussion of the relationships among protein families and conclusions generally seem appropriate, and the manuscript should be a useful contribution to the literature, especially if specific experimental tests could be proposed that would confirm the function of the unknown protein families and if the results could be confirmed with the additional methods suggested above.
In the revised manuscript, we proposed to test the ceramidase activity for hypothetical groups 1 and 2 since they are more closely related to alkaline ceramidases than to the other known groups. Since two of the CREST groups, alkaline ceramidase and Per1, catalyze similar yet distinct reactions on different substrates, we reason that the substrates for the other CREST groups, such as PAQRs, could also differ. We thus proposed in the revised manuscript that experimental studies of PAQRs should not be restricted to just the ceramidase activity. Other potential lipid hydrolase activities should also be considered. Specifically, we proposed in the revised manuscript that the phospholipase A2 activity (previously suggested for Per1) could be tested. The reaction products (a lysophospholipid molecule and a free fatty acid) of the phospholipase A2 activity on a phospholipid substrate can be precursors of second messengers that could function in PAQR signaling.
We would like to thank Lisa Kinch and Jeremy Semeiks for critical reading of the manuscript. This work was supported by Howard Hughes Medical Institute, National Institute of Health (GM094575 to NVG) and the Welch Foundation (I-1505 to NVG).
- Fagerberg L, Jonasson K, von Heijne G, Uhlen M, Berglund L: Prediction of the human membrane proteome. Proteomics. 2010, 10: 1141-1149. 10.1002/pmic.200900258.PubMedView ArticleGoogle Scholar
- Hurwitz N, Pellegrini-Calace M, Jones DT: Towards genome-scale structure prediction for transmembrane proteins. Philos Trans R Soc Lond B Biol Sci. 2006, 361: 465-475. 10.1098/rstb.2005.1804.PubMedPubMed CentralView ArticleGoogle Scholar
- Urban S, Shi Y: Core principles of intramembrane proteolysis: comparison of rhomboid and site-2 family proteases. Curr Opin Struct Biol. 2008, 18: 432-441. 10.1016/j.sbi.2008.03.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Martinez Molina D, Eshaghi S, Nordlund P: Catalysis within the lipid bilayer-structure and mechanism of the MAPEG family of integral membrane proteins. Curr Opin Struct Biol. 2008, 18: 442-449. 10.1016/j.sbi.2008.04.005.PubMedView ArticleGoogle Scholar
- Mao C, Obeid LM: Ceramidases: regulators of cellular responses mediated by ceramide, sphingosine, and sphingosine-1-phosphate. Biochim Biophys Acta. 2008, 1781: 424-434.PubMedPubMed CentralView ArticleGoogle Scholar
- Fyrst H, Saba JD: An update on sphingosine-1-phosphate and other sphingolipid mediators. Nat Chem Biol. 2010, 6: 489-497.PubMedPubMed CentralView ArticleGoogle Scholar
- Villa NY, Kupchak BR, Garitaonandia I, Smith JL, Alonso E, Alford C, Cowart LA, Hannun YA, Lyons TJ: Sphingolipids function as downstream effectors of a fungal PAQR. Mol Pharmacol. 2009, 75: 866-875. 10.1124/mol.108.049809.PubMedPubMed CentralView ArticleGoogle Scholar
- Holland WL, Miller RA, Wang ZV, Sun K, Barth BM, Bui HH, Davis KE, Bikman BT, Halberg N, Rutkowski JM, Wade MR, Tenorio VM, Kuo MS, Brozinick JT, Zhang BB, Birnbaum MJ, Summers SA, Scherer PE: Receptor-mediated activation of ceramidase activity initiates the pleiotropic actions of adiponectin. Nat Med. 2011, 17: 55-63. 10.1038/nm.2277.PubMedPubMed CentralView ArticleGoogle Scholar
- Baida GE, Kuzmin NP: Mechanism of action of hemolysin III from Bacillus cereus. Biochim Biophys Acta. 1996, 1284: 122-124. 10.1016/S0005-2736(96)00168-X.PubMedView ArticleGoogle Scholar
- Fujita M, Umemura M, Yoko-o T, Jigami Y: PER1 is required for GPI-phospholipase A2 activity and involved in lipid remodeling of GPI-anchored proteins. Mol Biol Cell. 2006, 17: 5253-5264. 10.1091/mbc.E06-08-0715.PubMedPubMed CentralView ArticleGoogle Scholar
- Maeda Y, Tashima Y, Houjou T, Fujita M, Yoko-o T, Jigami Y, Taguchi R, Kinoshita T: Fatty acid remodeling of GPI-anchored proteins is required for their raft association. Mol Biol Cell. 2007, 18: 1497-1506. 10.1091/mbc.E06-10-0885.PubMedPubMed CentralView ArticleGoogle Scholar
- Winston WM, Molodowitch C, Hunter CP: Systemic RNAi in C. elegans requires the putative transmembrane protein SID-1. Science. 2002, 295: 2456-2459. 10.1126/science.1068836.PubMedView ArticleGoogle Scholar
- Feinberg EH, Hunter CP: Transport of dsRNA into cells by the transmembrane protein SID-1. Science. 2003, 301: 1545-1547. 10.1126/science.1087117.PubMedView ArticleGoogle Scholar
- Su Z, Wang X, Shen S, Wang L, Li Y, Li N, Li Z: Expression of 2 transcripts of NGX6 gene in colorectal cancer and the correlation with carcinoembryonic antigen. Zhong Nan Da Xue Xue Bao Yi Xue Ban. 2010, 35: 401-408.PubMedGoogle Scholar
- Li J, Tan C, Xiang Q, Zhang X, Ma J, Wang JR, Yang J, Li W, Shen SR, Liang S, Li G: Proteomic detection of changes in protein synthesis induced by NGX6 transfected in human nasopharyngeal carcinoma cells. J Protein Chem. 2001, 20: 265-271. 10.1023/A:1010912311564.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Soding J, Biegert A, Lupas AN: The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33: W244-248. 10.1093/nar/gki408.PubMedPubMed CentralView ArticleGoogle Scholar
- Biegert A, Soding J: Sequence context-specific profiles for homology searching. Proc Natl Acad Sci USA. 2009, 106: 3770-3775. 10.1073/pnas.0810767106.PubMedPubMed CentralView ArticleGoogle Scholar
- Frickey T, Lupas A: CLANS: a Java application for visualizing protein families based on pairwise similarity. Bioinformatics. 2004, 20: 3702-3704. 10.1093/bioinformatics/bth444.PubMedView ArticleGoogle Scholar
- Pei J, Grishin NV: Peptidase family U34 belongs to the superfamily of N-terminal nucleophile hydrolases. Protein Sci. 2003, 12: 1131-1135. 10.1110/ps.0240803.PubMedPubMed CentralView ArticleGoogle Scholar
- Inoue T, Okino N, Kakuta Y, Hijikata A, Okano H, Goda HM, Tani M, Sueyoshi N, Kambayashi K, Matsumura H, Kai Y, Ito M: Mechanistic insights into the hydrolysis and synthesis of ceramide by neutral ceramidase. J Biol Chem. 2009, 284: 9566-9577. 10.1074/jbc.M808232200.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu R, Sun W, Jin J, Obeid LM, Mao C: Role of alkaline ceramidases in the generation of sphingosine and its phosphate in erythrocytes. FASEB J. 2010, 24: 2507-2515. 10.1096/fj.09-153635.PubMedPubMed CentralView ArticleGoogle Scholar
- Aerts AM, Zabrocki P, Francois IE, Carmona-Gutierrez D, Govaert G, Mao C, Smets B, Madeo F, Winderickx J, Cammue BP, Thevissen K: Ydc1p ceramidase triggers organelle fragmentation, apoptosis and accelerated ageing in yeast. Cell Mol Life Sci. 2008, 65: 1933-1942. 10.1007/s00018-008-8129-8.PubMedView ArticleGoogle Scholar
- Yang Q, Gong ZJ, Zhou Y, Yuan JQ, Cheng J, Tian L, Li S, Lin XD, Xu R, Zhu ZR, Mao C: Role of Drosophila alkaline ceramidase (Dacer) in Drosophila development and longevity. Cell Mol Life Sci. 2010, 67: 1477-1490. 10.1007/s00018-010-0260-7.PubMedView ArticleGoogle Scholar
- Tang YT, Hu T, Arterburn M, Boyle B, Bright JM, Emtage PC, Funk WD: PAQR proteins: a novel membrane receptor family defined by an ancient 7-transmembrane pass motif. J Mol Evol. 2005, 61: 372-380. 10.1007/s00239-004-0375-2.PubMedView ArticleGoogle Scholar
- Zhu Y, Rice CD, Pang Y, Pace M, Thomas P: Cloning, expression, and characterization of a membrane progestin receptor and evidence it is an intermediary in meiotic maturation of fish oocytes. Proc Natl Acad Sci USA. 2003, 100: 2231-2236. 10.1073/pnas.0336132100.PubMedPubMed CentralView ArticleGoogle Scholar
- Fernandes MS, Brosens JJ, Gellersen B: Honey, we need to talk about the membrane progestin receptors. Steroids. 2008, 73: 942-952. 10.1016/j.steroids.2007.12.004.PubMedView ArticleGoogle Scholar
- Yamauchi T, Kamon J, Ito Y, Tsuchida A, Yokomizo T, Kita S, Sugiyama T, Miyagishi M, Hara K, Tsunoda M, Murakami K, Ohteki T, Uchida S, Takekawa S, Waki H, Tsuno NH, Shibata Y, Terauchi Y, Froguel P, Tobe K, Koyasu S, Taira K, Kitamura T, Shimizu T, Nagai R, Kadowaki T: Cloning of adiponectin receptors that mediate antidiabetic metabolic effects. Nature. 2003, 423: 762-769. 10.1038/nature01705.PubMedView ArticleGoogle Scholar
- Yamauchi T, Kadowaki T: Physiological and pathophysiological roles of adiponectin and adiponectin receptors in the integrated regulation of metabolic and cardiovascular diseases. Int J Obes (Lond). 2008, 32 (Suppl 7): S13-18.View ArticleGoogle Scholar
- Ouchi N, Parker JL, Lugus JJ, Walsh K: Adipokines in inflammation and metabolic disease. Nat Rev Immunol. 2011, 11: 85-97. 10.1038/nri2921.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhu Y, Bond J, Thomas P: Identification, classification, and partial characterization of genes in humans and other vertebrates homologous to a fish membrane progestin receptor. Proc Natl Acad Sci USA. 2003, 100: 2237-2242. 10.1073/pnas.0436133100.PubMedPubMed CentralView ArticleGoogle Scholar
- Lyons TJ, Villa NY, Regalla LM, Kupchak BR, Vagstad A, Eide DJ: Metalloregulation of yeast membrane steroid receptor homologs. Proc Natl Acad Sci USA. 2004, 101: 5506-5511. 10.1073/pnas.0306324101.PubMedPubMed CentralView ArticleGoogle Scholar
- Hsieh MH, Goodman HM: A novel gene family in Arabidopsis encoding putative heptahelical transmembrane proteins homologous to human adiponectin receptors and progestin receptors. J Exp Bot. 2005, 56: 3137-3147. 10.1093/jxb/eri311.PubMedView ArticleGoogle Scholar
- Karpichev IV, Cornivelli L, Small GM: Multiple regulatory roles of a novel Saccharomyces cerevisiae protein, encoded by YOL002c, in lipid and phosphate metabolism. J Biol Chem. 2002, 277: 19609-19617. 10.1074/jbc.M202045200.PubMedView ArticleGoogle Scholar
- Kupchak BR, Garitaonandia I, Villa NY, Mullen MB, Weaver MG, Regalla LM, Kendall EA, Lyons TJ: Probing the mechanism of FET3 repression by Izh2p overexpression. Biochim Biophys Acta. 2007, 1773: 1124-1132. 10.1016/j.bbamcr.2007.04.003.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen YC, Chang MC, Chuang YC, Jeang CL: Characterization and virulence of hemolysin III from Vibrio vulnificus. Curr Microbiol. 2004, 49: 175-179.PubMedView ArticleGoogle Scholar
- Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, von Mering C: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011, 39: D561-568. 10.1093/nar/gkq973.PubMedPubMed CentralView ArticleGoogle Scholar
- Kall L, Krogh A, Sonnhammer EL: Advantages of combined transmembrane topology and signal peptide prediction--the Phobius web server. Nucleic Acids Res. 2007, 35: W429-432. 10.1093/nar/gkm256.PubMedPubMed CentralView ArticleGoogle Scholar
- Gonez LJ, Naselli G, Banakh I, Niwa H, Harrison LC: Pancreatic expression and mitochondrial localization of the progestin-adipoQ receptor PAQR10. Mol Med. 2008, 14: 697-704.PubMedPubMed CentralView ArticleGoogle Scholar
- Paidhungat M, Garrett S: Cdc1 and the vacuole coordinately regulate Mn2+ homeostasis in the yeast Saccharomyces cerevisiae. Genetics. 1998, 148: 1787-1798.PubMedPubMed CentralGoogle Scholar
- Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, Chen Y, Cheng X, Chua G, Friesen H, Goldberg DS, Haynes J, Humphries C, He G, Hussein S, Ke L, Krogan N, Li Z, Levinson JN, Lu H, Menard P, Munyana C, Parsons AB, Ryan O, Tonikian R, Roberts T, Sdicu AM, Shapiro J, Sheikh B, Suter B, Wong SL, Zhang LV, Zhu H, Burd CG, Munro S, Sander C, Rine J, Greenblatt J, Peter M, Bretscher A, Bell G, Roth FP, Brown GW, Andrews B, Bussey H, Boone C: Global mapping of the yeast genetic interaction network. Science. 2004, 303: 808-813. 10.1126/science.1091317.PubMedView ArticleGoogle Scholar
- Duxbury MS, Ashley SW, Whang EE: RNA interference: a mammalian SID-1 homologue enhances siRNA uptake and gene silencing efficacy in human cells. Biochem Biophys Res Commun. 2005, 331: 459-463. 10.1016/j.bbrc.2005.03.199.PubMedView ArticleGoogle Scholar
- Wolfrum C, Shi S, Jayaprakash KN, Jayaraman M, Wang G, Pandey RK, Rajeev KG, Nakayama T, Charrise K, Ndungo EM, Zimmermann T, Koteliansky V, Manoharan M, Stoffel M: Mechanisms and optimization of in vivo delivery of lipophilic siRNAs. Nat Biotechnol. 2007, 25: 1149-1157. 10.1038/nbt1339.PubMedView ArticleGoogle Scholar
- Tomoyasu Y, Miller SC, Tomita S, Schoppmeier M, Grossmann D, Bucher G: Exploring systemic RNA interference in insects: a genome-wide survey for RNAi genes in Tribolium. Genome Biol. 2008, 9: R10-10.1186/gb-2008-9-1-r10.PubMedPubMed CentralView ArticleGoogle Scholar
- Huvenne H, Smagghe G: Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review. J Insect Physiol. 2010, 56: 227-235. 10.1016/j.jinsphys.2009.10.004.PubMedView ArticleGoogle Scholar
- Miller JD, Arteca RN, Pell EJ: Senescence-associated gene expression during ozone-induced leaf senescence in Arabidopsis. Plant Physiol. 1999, 120: 1015-1024. 10.1104/pp.120.4.1015.PubMedPubMed CentralView ArticleGoogle Scholar
- De Rosa M, Gambacorta A, Gliozzi A: Structure, biosynthesis, and physicochemical properties of archaebacterial lipids. Microbiol Rev. 1986, 50: 70-80.PubMedPubMed CentralGoogle Scholar
- Sun W, Jin J, Xu R, Hu W, Szulc ZM, Bielawski J, Obeid LM, Mao C: Substrate specificity, membrane topology, and activity regulation of human alkaline ceramidase 2 (ACER2). J Biol Chem. 2010, 285: 8995-9007. 10.1074/jbc.M109.069203.PubMedPubMed CentralView ArticleGoogle Scholar
- Rosenbaum DM, Rasmussen SG, Kobilka BK: The structure and function of G-protein-coupled receptors. Nature. 2009, 459: 356-363. 10.1038/nature08144.PubMedPubMed CentralView ArticleGoogle Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol. 1995, 247: 536-540.PubMedGoogle Scholar
- Cuff AL, Sillitoe I, Lewis T, Redfern OC, Garratt R, Thornton J, Orengo CA: The CATH classification revisited--architectures reviewed and new ways to characterize structural divergence in superfamilies. Nucleic Acids Res. 2009, 37: D310-314. 10.1093/nar/gkn877.PubMedPubMed CentralView ArticleGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-222. 10.1093/nar/gkp985.PubMedPubMed CentralView ArticleGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A: MEROPS: the peptidase database. Nucleic Acids Res. 2010, 38: D227-233. 10.1093/nar/gkp971.PubMedPubMed CentralView ArticleGoogle Scholar
- Saier MH, Yen MR, Noto K, Tamang DG, Elkan C: The Transporter Classification Database: recent advances. Nucleic Acids Res. 2009, 37: D274-278. 10.1093/nar/gkn862.PubMedPubMed CentralView ArticleGoogle Scholar
- Sun W, Xu R, Hu W, Jin J, Crellin HA, Bielawski J, Szulc ZM, Thiers BH, Obeid LM, Mao C: Upregulation of the human alkaline ceramidase 1 and acid ceramidase mediates calcium-induced differentiation of epidermal keratinocytes. J Invest Dermatol. 2008, 128: 389-397.PubMedView ArticleGoogle Scholar
- Mayor S, Riezman H: Sorting GPI-anchored proteins. Nat Rev Mol Cell Biol. 2004, 5: 110-120. 10.1038/nrm1309.PubMedView ArticleGoogle Scholar
- Kupchak BR, Garitaonandia I, Villa NY, Smith JL, Lyons TJ: Antagonism of human adiponectin receptors and their membrane progesterone receptor paralogs by TNFalpha and a ceramidase inhibitor. Biochemistry. 2009, 48: 5504-5506. 10.1021/bi9006258.PubMedPubMed CentralView ArticleGoogle Scholar
- Burke JE, Dennis EA: Phospholipase A2 biochemistry. Cardiovasc Drugs Ther. 2009, 23: 49-59. 10.1007/s10557-008-6132-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Funk CD: Prostaglandins and leukotrienes: advances in eicosanoid biology. Science. 2001, 294: 1871-1875. 10.1126/science.294.5548.1871.PubMedView ArticleGoogle Scholar
- Noguchi K, Herr D, Mutoh T, Chun J: Lysophosphatidic acid (LPA) and its receptors. Curr Opin Pharmacol. 2009, 9: 15-23. 10.1016/j.coph.2008.11.010.PubMedView ArticleGoogle Scholar
- Ferguson JF, Phillips CM, Tierney AC, Perez-Martinez P, Defoort C, Helal O, Lairon D, Planells R, Shaw DI, Lovegrove JA, Gjelstad IM, Drevon CA, Blaak EE, Saris WH, Leszczynska-Golabek I, Kiec-Wilk B, Riserus U, Karlstrom B, Miranda JL, Roche HM: Gene-nutrient interactions in the metabolic syndrome: single nucleotide polymorphisms in ADIPOQ and ADIPOR1 interact with plasma saturated fatty acids to modulate insulin resistance. Am J Clin Nutr. 2010, 91: 794-801. 10.3945/ajcn.2009.28255.PubMedView ArticleGoogle Scholar
- Shih JD, Fitzgerald MC, Sutherlin M, Hunter CP: The SID-1 double-stranded RNA transporter is not selective for dsRNA length. RNA. 2009, 15: 384-390. 10.1261/rna.1286409.PubMedPubMed CentralView ArticleGoogle Scholar
- Matthews BW, Colman PM, Jansonius JN, Titani K, Walsh KA, Neurath H: Structure of thermolysin. Nat New Biol. 1972, 238: 41-43. 10.1038/238041a0.PubMedView ArticleGoogle Scholar
- Christianson DW, David PR, Lipscomb WN: Mechanism of carboxypeptidase A: hydration of a ketonic substrate analogue. Proc Natl Acad Sci USA. 1987, 84: 1512-1515. 10.1073/pnas.84.6.1512.PubMedPubMed CentralView ArticleGoogle Scholar
- Dideberg O, Charlier P, Dive G, Joris B, Frere JM, Ghuysen JM: Structure of a Zn2+-containing D-alanyl-D-alanine-cleaving carboxypeptidase at 2.5 A resolution. Nature. 1982, 299: 469-470. 10.1038/299469a0.PubMedView ArticleGoogle Scholar
- Feng L, Yan H, Wu Z, Yan N, Wang Z, Jeffrey PD, Shi Y: Structure of a site-2 protease family intramembrane metalloprotease. Science. 2007, 318: 1608-1612. 10.1126/science.1150755.PubMedView ArticleGoogle Scholar
- Eddy SR: A new generation of homology search tools based on probabilistic inference. Genome Inform. 2009, 23: 205-211.PubMedGoogle Scholar
- Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007, 2: 953-971. 10.1038/nprot.2007.131.PubMedView ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Pei J, Grishin NV: PROMALS: towards accurate multiple sequence alignments of distantly related proteins. Bioinformatics. 2007, 23: 802-808. 10.1093/bioinformatics/btm017.PubMedView ArticleGoogle Scholar
- Adachi J, Hasegawa M: MOLPHY version 2.3, programs for molecular phylogenetics based on maximum likelihood. Computer Science Monographs (The Institute of Statistical Mathematics). 1996, 28: 1-150.Google Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Hasegawa M, Kishino H, Saitou N: On the maximum likelihood method in molecular phylogenetics. J Mol Evol. 1991, 32: 443-445. 10.1007/BF02101285.PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.PubMedView ArticleGoogle Scholar
- Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55: 539-552. 10.1080/10635150600755453.PubMedView ArticleGoogle Scholar
- Le SQ, Gascuel O: An improved general amino acid replacement matrix. Mol Biol Evol. 2008, 25: 1307-1320. 10.1093/molbev/msn067.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.