Services on Demand
- Cited by SciELO
- Access statistics
- Cited by Google
- Similars in SciELO
- Similars in Google
Print version ISSN 1415-4757
On-line version ISSN 1678-4685
Genet. Mol. Biol. vol.32 no.3 São Paulo 2009
GENOMICS AND BIOINFORMATICS
Michely C. DinizI; Marcília P. CostaI; Ana C.L. PachecoI; Michel T. KamimuraI; Samara C. SilvaI; Laura D.G. CarneiroI; Ana P.L. SousaI; Carlos E.A. SoaresI; Celeste S.F. SouzaII; Diana Magalhães de OliveiraI
INúcleo Tarcísio Pimenta de Pesquisa Genômica e Bioinformática, Universidade Estadual do Ceará, Fortaleza, CE, Brazil
IILaboratório de Imunomodulação e Protozoologia, Fundação Oswaldo Cruz, Manguinhos, RJ, Brazil
Several motile processes are responsible for the movement of proteins into and within the flagellar membrane, but little is known about the process by which specific proteins (either actin-associated or not) are targeted to protozoan flagellar membranes. Actin is a major cytoskeleton protein, while polymerization and depolymerization of parasite actin and actin-interacting proteins (AIPs) during both processes of motility and host cell entry might be key events for successful infection. For a better understanding the eukaryotic flagellar dynamics, we have surveyed genomes, transcriptomes and proteomes of pathogenic Leishmania spp. to identify pertinent genes/proteins and to build in silico models to properly address their putative roles in trypanosomatid virulence. In a search for AIPs involved in flagellar activities, we applied computational biology and proteomic tools to infer from the biological meaning of coronins and Arp2/3, two important elements in phagosome formation after parasite phagocytosis by macrophages. Results presented here provide the first report of Leishmania coronin and Arp2/3 as flagellar proteins that also might be involved in phagosome formation through actin polymerization within the flagellar environment. This is an issue worthy of further in vitro examination that remains now as a direct, positive bioinformatics-derived inference to be presented.
Key words: actin-interacting proteins (AIPs), flagellar proteins, Leishmania, coronin and Arp2/3, phagosome.
Leishmania spp. is a trypanosomatid protozoan responsible for a complex group of diverse clinical forms of classically neglected diseases collectively known as leishmaniases (Peters and Killick-Kendrick, 1987). As an insinuating and persistent pathogen, Leishmania has a specialized organelle for motility, the flagellum, which is essential for parasite migration, invasion and persistence on host tissues. Several motile processes have been reported in the literature that could be responsible for the movement of proteins into and within the flagellar membrane (Kozminski et al., 1993). However, little is known about the process by which specific proteins (either actin-associated or not) are targeted to protozoan flagellar membranes, with few exceptions such as the work by Snapp and Landfear (1999) which targets a motif for a flagellar integral membrane protein in Leishmania enriettii. Actin is a major cytoskeleton protein, while polymerization and depolymerization of parasite actin and actin motor-associated proteins during both processes of motility and host cell entry might be key events for successful infection.
To gain a better understanding of flagellar dynamics in Leishmania, we have surveyed genomes, transcriptomes and proteomes of these pathogens to identify pertinent genes/proteins and to build in silico models to properly address their putative roles in trypanosomatid virulence. We have applied computational tools (hidden Markov models, Viterbi algorithm and comparative modeling) to infer biological meaning through detailed sequence-structural-functional analyses on actin-interacting proteins (AIPs) involved in flagellar activities of Leishmania spp. The actin system evolved to be fail-safe with multiple proteins sharing overlapping but novel roles, whereas these multifaceted roles are likely to provide the versatile scenario involving actin-interacting activities.
Proteins such as coronin and Arp2/3 complex have a common feature of actin-binding activity and might be involved in Leishmania intraflagellar pathways (Broadhead et al., 2006; Costa et al., 2007). The coronin family comprises two groups of evolutionary conserved WD-repeat proteins known to interact with Arp2/3 complex that help regulate the nucleation dynamics of actin filaments (Rybakin and Clemen, 2005). Invasion and differentiation processes in Leishmania directly involve phagocytosis and the formation of an intracellular membrane-bounded organelle called a phagosome, since Leishmania promastigotes must first evade complement-mediated lysis until they are engulfed by a mammalian host macrophage (Puentes et al., 1989; Olivier et al., 2005).
By phagocytosis we mean the process of ingestion per se. A phagosome is the endocytic compartment that contains a non-interfering particle, as opposed to a vacuole which describes a compartment containing a particle, such as a pathogen like Leishmania that diverts normal phagosome maturation. Phagosome maturation refers to the process of intracellular phagosome development after closure of the phagocytic cup (reviewed by Haas, 2007). Phagocytosis has evolved into a highly complex and regulated process in multicellular eukaryotes by which microbes and other particles are taken up into the phagosome (Griffiths and Mayorga, 2007). The precise mechanisms by which coronin and Arp2/3 contribute to phagocytosis are not known (Yan et al., 2005), but it is crucial to investigate elements that will shed more light on these mechanisms.
Although not clearly understood, it has been shown in vitro that, after the phagocytosis of L. donovani by macrophage-like cells J774, parasites are transiently located in phagosomes with poor fusogenic properties towards late endocytic compartments (Desjardins and Descoteaux, 1997). In contrast, after internalization, they are found in compartments that rapidly fuse with late endocytic organelles (Dermine et al., 2000). A possibility is that these distinctive features of the early phagosomes could be linked to different proteins expressed on the parasite surface, which may modify the fusion capacity of the phagosomal membrane (Desjardins and Descoteaux, 1997; Dermine et al., 2000). Furthermore, most of the Leishmania-containing phagosomes have been shown to accumulate F-actin, which is noted around these phagosomes (Holm et al., 2001).
We should recall, then, the early events following phagocytosis of Leishmania, pre-adapted to the encounter within intracellular conditions of mammals. As reported (Courret et al., 2002), young phagosomes containing Leishmania rapidly acquire a competence to fuse with late endosomes/lysosomes. If we take this version, AIPs have a greater chance of being key elements on the phagosome maturation than any of the many surface membrane proteins that have been widely investigated. For that reason, we have focused our studies on AIPs because they are rich in motif-binding activities, besides being excellent models for sequence and structural comparisons, predictions and examinations to be carried out with all kinds of restraints on spatial structure of the amino acid sequence(s) and ligands. Our coronin and Arp2/3 model restraints were derived from known related protein structures (comparative modeling) and from rules of secondary structure packing (combinatorial modeling). Our in silico examinations provide information about AIPs in the specific domain of phagosome function in Leishmania spp. infection. The role of the phagosome is to deliver particulate material to a hydrolytic environment that will lead to its degradation. The maturation process of the phagocytic compartment is linked intimately to both digestive housekeeping processes and innate sensing of molecules associated with infection (Vieira et al., 2002). Therefore, mechanisms and regulation behind this maturation process are strategic for both hosts and pathogens, whereas here we present and discuss results concerning Leishmania AIPs that are also flagellar proteins that might be targeted in phagosome formation.
Biological databases and bioinformatics tools
As previously (Oliveira et al., 2005; Gouveia et al., 2007; Vasconcelos et al., 2007, 2008), we have used publicly available datasets of individual/clusters of gene/protein data on Leishmania spp. (from GeneDB, a core part of Sanger Institute Pathogen Genomics) and from other available eukaryotic organisms (at NCBI and UniProt/Swiss-Prot/trEMBL knowledge DB), including sequences from the genome projects of C. reinhardtii in the U.S. Department of Energy (JGI) and all available data at The Chlamydomonas Flagellar Proteome Project. For databases (DB) searches, accession numbers/identifiers are those used in four main sources: NCBI, GeneDB, PDB and UniProt. As previously described (Oliveira et al., 2005; Gouveia et al., 2007; Vasconcelos et al., 2007, 2008), BLAST and its variants (Altschul et al., 1997) and MUSCLE (Edgar, 2004) were used for sequence similarity searches and comparisons/analysis through pairwise and multiple alignments. As previously detailed (Costa et al., 2007), for pattern recognition tasks of gene/protein predictions, motif finding and refinements for defining core domains of interest, HMM implementations such as HMMER (Eddy, 1998) and a structural prediction adaptation (Yoshizawa et al., 2006) were employed, with the additional improvement of the Viterbi algorithm (VA) (Forney, 1973; Black, 2004) for best possible alignments after recursion (Alexandrov and Gerstein, 2004). As previously (Costa et al., 2007), we assume that the HMMs generated the input sequence and are looking for the highest probability path. For sequence position i = 0, 1 , , L+1; for state l = 0, 1, , n:
Vl(i + 1) = max[Vk(i) akl] el(Xi + 1)
As described in previous studies (Oliveira et al., 2005; Gouveia et al., 2007; Vasconcelos et al. 2008), we took alignments created with FASTA/BLAST as input and computed alignment tables, providing hierarchical and successive correlations between each of the two sets of sequences. FASTA files for amino-acid (aa) sequences of coding regions were downloaded from sources detailed above. All flagellar proteins and AIPs in Leishmania were collected from experimental papers and from the comprehensive GeneDB database. An illustration of the tools used is seen in Figure 1. Briefly, these sequences were used as queries against genomic datasets with PSI-BLAST (BLASTP2.2.10) (Altschul et al., 1997) to determine sequence similarity among all possible sequences, followed by multiple sequence alignments (MSA) performed with MUSCLE (Edgar, 2004) on target entries of main source DBs searched against various collections of protein motifs and families. Results of MSA were used as training datasets for HMM profiles (Eddy, 1998; Yoshizawa et al., 2006) and the SMART (Letunic et al., 2006) motif patterns of cluster motifs. The quality of trained HMM predictions and SMART pattern matching was examined. The HMM profiles for cluster motifs were trained, calibrated, and used for DB searches using HMMER package ver.2.3.1 with default parameters. Gene ontology (GO) terms were assigned, based on top matches to proteins with GO annotations from Swiss-Prot/trEMBL, AMIGO after GeneDB and TargetP access. Functional assignment of these genes/gene products was inferred using an RPS-BLAST search against conserved domain DBs (CDD) (Marchler-Bauer et al., 2005); information was taken into account about subcellular localization (Emanuelsson et al., 2000), sequence and structural features, domains/motifs conservation (von Mering et al., 2005; Letunic et al, 2006), and in vitro characterization (Avidor-Reiss et al., 2004; Tull et al., 2004). For three-dimensional (3D) modeling of AIPs we employed Modeller (Sali and Blundell, 1993; Marti-Renom et al., 2000) and ESyPred3D (Lambert et al., 2002). For visualization of 3D modeled structures, we employed the browser plug-in Chime. For addressing particular sites of structural-functional relationships of AIPs we used Cn3D (NCBI) and STING (Higa et al., 2004). Best models were evaluated by using PROCHECK (Laskowski et al., 1993) and a thorough analysis of consensus AIPs structures was carried out for model dissection, superposition and rms deviation (rmsd) calculations.
Samples. Purified fractions of flagella from Leishmania amazonenesis, strain H21 (MHOM/BR/76/ MA-76), were kindly provided by Fiocruz (Fundação Oswaldo Cruz, Manguinhos, Rio de Janeiro). An amount of 250 µg of flagellar protein extract was normalized for use.
Two-dimensional (2D) gel electrophoresis. We performed proteomic analyses through 2D polyacrylamide gel electrophoresis (2D-PAGE), as previously described (Drummelsmith et al., 2003, Brobey et al., 2006). Briefly, proteins were treated with a lysis buffer with 7 M urea, 2 M thiourea, 2% NP-40, 2% DTT and 2% ampholytes (pH 3 to 10). In the first dimension, a protein separation by isoelectric point was performed applying IPG strips, soaked in the sample solution, at 300V for 5 min, 3,500 V for 5 h through Ettan IPGphor II Isoelectric Focusing System (GE Healthcare). To promote efficient protein transfer from the first to the second dimension, IPG strips were incubated in reducing buffer (1% DTT, 6M urea, 30% glycerol, 2% SDS, 50 mM Tris, pH 8.8) for 15 min, followed by 15 min incubation in alkylation buffer (2.5% iodoacetamide, 6 M urea, 30% glycerol, 2% SDS, 50 mM Tris, pH 8.8). The strips were immediately applied to the second dimension 12.5% polyacrylamide gel. Strips were, then, overlaid on freshly poured 10% Tricine-SDS gels (18 x 16 cm) and sealed with agarose solution (0.5% agarose, plus a few grains of bromophenol blue in a Tris-tricine cathode buffer). Protein standards of molecular mass range 14400 Da to 97000 Da were used (GE Healthcare). Gel electrophoresis was carried out at 250 V for 1 h, and then at 500 V until the dye front reached the bottom of the gel (around 8 h). Gels were silver-stained as previously described (Gromova and Celis, 2006).
After running and scanning the gels, the resulting protein gel was analyzed using ImageMaster 2-D Platinum 6.0® (Amersham Biosciences, Uppsala, Sweden). The spots detected automatically by the software were visually inspected. Spot filtering and editing were performed manually to remove artifacts and to correct for spots that did not split correctly or were not detected by the software's automatic spot detection process. Molecular mass (MW) and isoeletric point (pI) were predicted for the best spots.
To identify putative flagellar proteins of Leishmania participating in the actin system, we have taken into account all actin-binding, interacting or related/regulated proteins as source data to investigate their possible involvement in assembly of F-actin (the principal driving force behind many forms of parasite locomotion, including those regarding or underlying internalization into the host cell). Our underlying hypothesis is that polymerization and depolymerization of parasite actin and actin motor-associated proteins, during both processes of motility and host cell entry, might be key events for successful infection, including the parasite survival within phagosomes. Recently (Vasconcelos et al., 2007, 2008) we reported an actin-polymerization Leishmania protein, profilin, and its partner, formin, as putative flagellar proteins due to their likely involvement in axonemal assembly/disassembly, therefore, flagellar dynamics/remodeling. In addition, we predicted other Leishmania AIPs, coronins and Arp2/3 complex as being flagellar proteins after detailed in silico structural alignments (Costa et al., 2007). Now, we focus on Leishmania coronin and Arp2/3 complex proteins as flagellar elements to phagosome formation in an attempt to directly associate the flagellum to the phagosomes, which are pivotal organelles in the ability of macrophages to perform several of their key functions, such as the handling of apoptotic cells, tissue remodeling, and restriction of the establishment and spread of intracellular pathogens (Méresse et al., 1999), such as Leishmania. Comparative bioinformatics analyses (performed after in silico genomic searches and in vitro proteomic screenings) have led to the identification of flagellar proteins that are involved not only in cytoskeleton activity, signaling, endocytosis, and lytic activity, but also putatively in the phagosome formation. With this study, we are able to confirm previous data on proteins already shown to be involved in actin systems, but not previously reported as Leishmania- or as flagellar-associated.
Flagellar coronin identification in Leishmania
Coronins constitute an evolutionary conserved family of WD40-repeat (WDR) actin-binding proteins (reviewed by de Hostos, 1999; Rybakin and Clemen, 2005), originally described in the social amoeba Dictyostelium discoideum (de Hostos et al., 1991). They can be divided, according to structure and function, into two subfamilies: one of shorter proteins (45-50kDa) implicated in, among other functions, nucleation of F-actin; and the other of longer, around 900 aa length, highly homologous coronins (POD-1 and Crn7) found in C. elegans, Drosophila and D. discoideum (Appleton et al., 2006; Rybakin et al., 2006). Until now, biochemical activities of coronins have been largely unknown, although one report (Humphries et al., 2002) identified how coronin and Arp2/3 complex interact in vitro. The comparative alignment of primary aa sequences of Leishmania coronins showed that they are highly conserved and have similar domain structures (Figure 2a), sharing ~45% aa identity with D. discoideum. Our models for the predicted 3D structure of Leishmania coronins, seen in Figure 2b, also reinforce their close proximity to PDB solved structures for coronin.
The proteomic screening performed on the flagellar fraction of L. amazonensis enabled us to anticipate evidence of coronin, among other proteins, as at least six putative spots, as shown in Figure 3 which provides a fairly good insight on Leishmania flagellar proteins (manuscript in preparation) that are also classical AIPs, such as coronin. Since coronin can be required for an early step in phagosome formation (Yan et al., 2005), which is consistent with its role in actin polymerization and accumulation around phagosomes formed during the ingestion of mycobacteria (Ferrari et al., 1999) and also in preventing phagosome maturation and mycobacterial killing (Ferrari et al., 1999), here we have considered this information as an important element for assigning a possible phagosome role for Leishmania coronin as well. We must recall that Leishmania has differential levels of survival inside phagosomes, as has been elegantly reported by Gueirard et al. (2007), whose work with L. donovani shows degraded parasites in spacious phagosomes, in contrast to morphologically intact parasites in tight compartments within neutrophils. These results suggest that the survival of parasites is linked to their ability to be targeted to tight/non-cidal and non-degradative compartments.
It is of paramount interest to mention that a phagosome, when it normally matures into a phagolysosome, undergoes a transition that correlates with functional changes from an organelle with early endosomal characteristics to a compartment with lysosomal, degradative properties. If we consider that persistent accumulation of coronin around phagosomes containing pathogens is believed to prevent phagosome maturation and mycobacterial killing (Ferrari et al., 1999), it would be conceivable to determine how coronin, actively recruited by these pathogens, may be involved in the modulation of the actin cytoskeleton, thereby influencing intracellular trafficking and survival (Jayachandran et al., 2008). While broad investigations into this are ongoing, the few pieces of information available stem largely from studies on phagosomes containing pathogens such as Leishmania that force the maturation machinery to slow down but will not permanently inhibit phagolysosome formation, with growth dependent upon eventual phagolysosome formation (Haas, 2007).
To what extent these different killing pathways (delayed or not) are used in a given phagosome certainly depends on various cellular and environmental conditions (usually a plethora of macrophage receptors at variable proportions that may influence phagosome fate to variable degrees), including interactions with the flagellum apparatus. We must bear in mind that flagella are dynamic structures that exchange up to 20% of their polypeptides within 3 h without any changes in their overall length (Song and Dentler, 2001; Wiese et al., 2003). Because flagellar microtubules are constantly turning over, the phagosomal membrane remains tightly associated with the parasites but the compartments become shorter (Marshal and Rosembaum, 2001). This reduction in size correlates with the progressive loss of the flagellum, suggesting that the parasite remodeling is accompanied by the removal of the phagosomal membrane.
Proteins involved in actin cytoskeleton dynamics that were recently implicated in phagocytosis in another protozoan, Entamoeba histolytica, are actobindin, coactosin and formin-like proteins (Marion et al., 2005), members of the so-called "core machinery" involved in de novo polymerization of actin filaments, which includes several subunits of the Arp2/3 complex, CAP protein, profilin and ADF/cofilin protein family (Loisel et al., 1999; Marion et al., 2005). ADF/cofilin proteins are essential to activate the turnover of actin filaments in dynamic regions of the cell and have already been shown to promote the remodeling of the actin network beneath the phagocytic cup (Bamburg, 1999; Bierne et al., 2001) and they have also been the subject of our own bioinformatics analyses (Pacheco et al., 2007) and of a recent in vitro identification (Tammana et al., 2008) in Leishmania flagellum.
Here we have chosen to follow the hypothesis that the actin filaments beneath the phagocytic cup could be organized as a bundle-like structure, rather than a dendritic network, which is necessary to form a stiff actin network around the phagocytic cup thereby facilitating the contractile activity that closes the vacuole (Marion et al., 2005). It is tempting to recall the model in which Toxoplasma gondii-containing vacuoles possess deep invaginations supported by host microtubules that deliver host endosomes and lysosomes into the parasitophorous vacuole (Coppens et al., 2006), where the parasites feed on these endocytic organelles and not, as expected, on the Toxoplasma vacuole arising completely separate from the endocytic pathway (Sibley, 2003).
Leishmania Arp2/3 complex genes
Arp2/3 complex has emerged as a central effector of actin assembly (Higgs and Pollard, 2001; Humphries et al., 2002). It is composed of seven evolutionarily conserved subunits: two actin-related proteins (Arp2 and Arp3) and five other subunits (yeast Arc40, Arc35, Arc18, Arc19, and Arc15). Since infection of mammalian hosts with Leishmania protozoa does depend on the ability of these parasites to replicate within macrophage phagolysosomes, it is noteworthy to link the function of Arps, such as the multiprotein Arp2/3 complex, and the phagosome formation, as a key step in parasite persistence within the host. Previously we have catalogued a list of Leishmania Arp2/3 complex genes comprising six putative sequences in L. braziliensis and L. infantum, and five in L. major(Costa et al., 2007), which are highly homologous with well-characterized Arp2/3 proteins in different organisms (as seen in Figure 4a). In Figure 4b we show a panel of four modeled Arp2/3 proteins in Leishmania where all of them possess the WDR, but no other known binding sites. WDR domains dominantly interfere with coronin function, acting as multimolecular scaffolding domains by bringing together interacting proteins on a single surface (Oku et al., 2003; Yan et al., 2005). They are likely to interfere with normal cell function by scavenging interacting proteins away from endogenous, full-length coronin.
If Arp2/3 complex is, in fact, one of those proteins that interact with coronin (Humphries et al., 2002; Yan et al., 2005), other portions (except WDR) of the protein are responsible for coronin recruitment to sites of actin remodeling, which would, then, link coronin and Arp2/3 to the sites of dynamic actin remodeling, with the contradiction that the WDR domain of coronin would act in a dominant-negative manner by preventing Arp2/3 accumulation. Conflicting interpretations (Humphries et al., 2002; Yan et al., 2005) pose a likely, unforeseen alternative: that coronin might play both positive (stimulating) and negative (inhibiting) roles in actin polymerization, depending on whether Arp2/3 functions prevent branching of F-actin or not.
Since our flagellar proteomic analyses on L. amazonensis did demonstrate clear support for the presence of at least four Arp2/3 proteins, as shown in Figure 3, we put forward an interesting point of view regarding the possibility that different sources of Arp2/3 proteins, whether flagellar or from the cytoskeleton, might be responsible in part for such paradox on activities of the coroninArp2/3 interaction.
Some Leishmania flagellar AIPs are homologous to phagosome proteins
To find out which AIPs in Leishmania were homologous to known phagosome proteins, as characterized in other studies (Desjardins et al., 1994; Garin et al., 2001; Marion et al., 2005; Stuart et al., 2007; Haas, 2007), we have searched lists of reported in vitro phagosome assays that gave us a panel ranging from 140-617 proteins (Garin et al., 2001; Gotthardt et al., 2006; Stuart et al., 2007) that are studied due to their direct involvement in phagosome formation. Therefore, we have compared their data to that found in Leishmania genomes, proteomes and particularly in our L. amazonensis flagellar proteome. We then collected a selected list of Leishmania proteins, including coronin and Arp2/3, ascribed as highly homologous to mammalian phagosome proteins (Table 1). We built a chart to emphasize these comparisons (Figure 5), which is a modified view of the full virtual phagosome (Garin et al., 2001) to cope with our predictive analysis of parasite flagellar proteins as homologous to mammalian phagosome counterparts. Therefore, in this study, we attempt to make a small contribution to further comprehensive analysis of Leishmania-containing phagosomes from which most of proteins have yet to be identified.
Many groups are currently looking into these various aspects of phagocyte biology and progress is steady, although no comprehensive phagosome proteome or lipidome has been published for any pathogen-containing vacuole (Haas, 2007). There is an amazing diversity of compartments that intracellular pathogens inhabit, and they appear to use a stunning number of different factors to establish these compartments. The real challenge remains in understanding how phagocytes and their phagosomes co-operate in pathological situations to phagosome biogenesis, as well as how to turn a phagosome from a hospitable into a killing environment for Leishmania.
Phagocytosis is initiated by recognition of a pathogen by a host cell (usually macrophages and neutrophils) receptors that trigger its engulfment into the phagosome. Following their attachment to the macrophage, Leishmania promastigotes are internalized to the relatively benign environment of the endosome, where they begin to differentiate into amastigotes. Unlike amastigotes, promastigotes are vulnerable to degradation by the acidic and hydrolytic environment of the phagolysosome (Olivier et al., 2005). They must therefore retard endosome maturation and phagosome-endosome fusion, a process that has been observed by the absence or delayed arrival of late endosomal markers such as rab7 and LAMP-1 (Scianimanico et al., 1999) and is related to the accumulation of F-actin (Holm et al., 2001). The mechanism is not completely understood (Olivier et al., 2005), but it has been shown that the delay in phagolysosome maturation provides a window during which promastigotes can differentiate into the more resistant amastigotes.
Nascent phagosomes have a composition that resembles that of the plasmalemma and are unable to digest the vacuolar contents. This capability is acquired by remodeling of the phagosomal membrane and contents through a complex series of fusion events with compartments of the endocytic and secretory pathways (Desjardins and Descoteaux, 1997). The components delivered to maturing phagosomes include acid hydrolases, as well as vacuolar-type ATPases that are responsible for acidification of the phagosomal lumen. This maturation process culminates in the formation of the phagolysosome, a highly acidic organelle (pH < 5.5) where degradation occurs.
Here we report the association of two classical AIPs (coronins and Arp2/3) to the flagellar compartment in Leishmania, certainly a reasonable indication to link flagellar proteins to a possible F-actin binding in this protozoan. In vitro, purified coronin is said to bind specifically to F-actin, to bundle actin filaments, and to weakly promote actin assembly (David et al., 1998; Fukui et al., 1999), whereas, in Acanthamoeba, coronin was localized in the cell's periphery (on the leading edge) consistent with that of actin located around the phagocytic cups (Baldo et al., 2005). Understanding the molecular mechanisms of phagosome maturation is critical because a number of pathogens, particularly Leishmania, survive inside host cells through subversion of this process (Rosenberger and Finlay, 2003). Keeping these facts in mind, together with our own results, we, then, believe that flagellar AIPs, such as coronin and Arp2/3/, might play a role in these subversion mechanisms that maintain the parasite viable within the phagosome.
Some reports have shown that Leismania can persist in mammalian neutrophils for up to 48 h, a period that largely exceeds the normal life span of neutrophils, indicating that the phagocytosis/infectious process might alter the ability of neutrophils to initiate their programmed cell death, increasing the lifetime of neutrophils by delaying apoptosis (Aga et al., 2002; Gueirard et al., 2007). It seems that Leishmania may have co-evolved with their mammalian hosts (Table 1 and Figure 5) to take advantage of phagosomes to establish a privileged niche for the transient parasitism of leukocytes and their subsequent invasion of macrophages, since it has been shown that the delay in phagolysosome maturation provides a window during which promastigotes can differentiate into the more resistant amastigotes (Olivier et al., 2005).
One could expect that phagosomes containing a given type of particle that entered cells simultaneously via the same receptor would behave the same, at least in a single cell. Surprisingly, however, phagosomes formed via the same receptors are found in different chemical states even within the same macrophage (Griffiths, 2004), a notion that is very useful to help our speculations on Leishmania AIPs rich content and at the same time homologous to a few phagosome proteins. Each phagosome is an individual entity whose behavior depends on a finite number of stable equilibrium states in its membrane signaling networks (Griffiths, 2004), which would, thus, underline a diversity of available elements to disturb effective phagosome formation.
As stated by Rasmusson and Descoteaux (2004), L. donovani promastigotes inhibit phagolysosome biogenesis in a lipophosphoglycan (LPG)-dependent manner, which correlates with an accumulation of periphagosomal F-actin, forming a physical barrier that prevents L. donovani promastigote-containing phagosomes from interacting with endocytic vacuoles. Such inhibition of phagosome maturation may constitute a strategy to provide an environment propitious to promastigote-to-amastigote differentiation. Interestingly, some molecules, known to be involved in apoptotic signaling, are thought to accumulate in phagosome membrane microdomains, suggesting that this location might be important to their specific function, such as LPG in a disruptive process in the phagosome (Gueirard et al., 2007), which could imply an involvement in the ability of Leishmania to inhibit neutrophil apoptosis and, thus, the posterior macrophage activity that also depends on actin.
Assuming that the driving force for host cell entry involves polymerization of parasite actin and its AIPs, plus the recognized subcellular localization of functional coronins and Arps shown in the phagosome (Bricheux et al., 2000; Asano et al., 2001; Baldo et al., 2005), then we can propose a possible viable role for flagellar coronin and Arp2/3 within phagosome. As Love et al. (1998) have shown, periphagosomal actin is rapidly lost when parasites are internalized. These authors have mentioned a higher percentage of L. major promastigotes surrounded by actin after experimentally adding parasites to macrophages. These findings suggest a clear binding to F-actin and implicit evidence of intense actin-interacting activity during phagosome formation after host cell entry (promastigote phagocytosis).
It seems quite obvious that intracellular parasites living within the harsh environment of phagocytes have developed strategies that allow them to adapt quickly, escape from first-line defense systems, and inhibit several functions of their host cells (Olivier et al., 2005). One of these strategies might be explained through the complex interactions of the actin system (Sturgill-Koszycki et al., 1996; Tardieux et al., 1998; Ullrich et al., 1999; Vieira et al., 2002; Yan et al., 2005). Coronins were first identified as AIPs, but, although they are crucial to the dynamics of actin filaments and cell movement, their actin-binding sites have been difficult to pin down, which reinforces the significance of computational biology and proteomic predictions that focus on coronin direct ligands, such as Arp2/3. Moreover, establishing relationships of subcellular/organellar localizations, such as those we now report for coronin and Arp2/3 on the flagellum of L. amazonensis, will help to clarify unforeseen roles that parasite AIPs might play on phagocytosis and phagosome formation through direct actin polymerization within the flagellar dynamic environment. This is an issue worthy of further in vitro examination that remains now as a direct, positive bioinformatics-derived inference to be presented. Research into phagosome biogenesis has flourished in recent years (Haas, 2007) and it will surely lead to a better understanding of Leishmania pathogenesis and mechanisms of immune invasion that involve the flagellar apparatus.
We wish to thank Fundação Oswaldo Cruz (FIOCRUZ, Brazil) for providing biological samples and GeneDB (The Wellcome Trust Sanger Institute, Pathogen Sequencing Unit) for full access to bioinformatics tools and databases on Leishmania genomes. This work is supported in part by CNPq and FUNCAP through individual grants to DMO and graduate fellowships to MCD, MPC, ACLP and SCS. MCD is a PhD student from RENORBIO (Rede Nordeste de Biotenologia) and this research is part of her thesis work.
Aga E, Katschinski DM, van Zandbergen G, Laufs H, Hansen B, Muller K, Solbach W and Laskay T (2002) Inhibition of the spontaneous apoptosis of neutrophil granulocytes by the intracellular parasite Leishmania major. J Immunol 169:898-905. [ Links ]
Alexandrov V and Gerstein M (2004) Using 3D Hidden Markov Models that explicitly represent spatial coordinates to model and compare protein structures. BMC Bioinform 5:2. [ Links ]
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W and Lipman DJ (1997) Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res 25:3389-402. [ Links ]
Appleton BA, Wu P and Wiesmann C (2006) The crystal structure of murine coronin-1: A regulator of actin cytoskeletal dynamics in lymphocytes. Structure 14:87-96. [ Links ]
Asano S, Mishima M and Nishida E (2001) Coronin forms a stable dimer through its C-terminal coiled coil region: An implicated role in its localization to cell periphery. Genes Cells 6:225-235. [ Links ]
Avidor-Reiss T, Maer AM, Koundakjian E, Polyanovsky A, Keil T, Subramaniam S and Zuker CS (2004) Decoding cilia function: Defining specialized genes required for compartmentalized cilia biogenesis. Cell 117:527-539. [ Links ]
Baldo ET, Moon EK, Kong HH and Chung DI (2005) Acanthamoeba healyi: Molecular cloning and characterization of a coronin homologue an actin-related protein. Exp Parasitol 110:114-22. [ Links ]
Bamburg JR (1999) Proteins of the ADF/cofilin family: Essential regulators of actin dynamics. Annu Rev Cell Dev Biol 15:185-230. [ Links ]
Bierne H, Gouin E, Roux P, Caroni P, Yin HL and Cossart P (2001) A role for cifinin and LIM kinase in Listeria-induced phagocytosis. J Cell Biol 155:101-112. [ Links ]
Bricheux G, Coffe G, Bayle D and Brugerolle G (2000) Characterization, cloning and immunolocalization of a coronin homologue in Trichomonas vaginalis. Eur J Cell Biol 79:413-22. [ Links ]
Broadhead R, Dawe HR, Farr H, Griffiths S, Hart SR, Portman N, Shaw MK, Ginger ML, Gaskell SJ, McKean PG, et al. (2006) Flagellar motility is required for the viability of the bloodstream trypanosome. Nature 440:224-227. [ Links ]
Brobey RKB, Mei FC, Cheng X and Soong L (2006) Comparative two-dimensional gel electrophoresis maps for promastigotes of Leishmania amazonensis and Leishmania major. Braz J Infect Dis 10:1-6. [ Links ]
Coppens I, Dunn JD, Romano JD, Pypaert M, Zhang H, Boothroyd JC and Joiner KA (2006) Toxoplasma gondii sequesters lysosomes from mammalian hosts in the vacuolar space. Cell 125:261-274. [ Links ]
Costa MP, Lucas HR, Maia ARS, Pacheco ACL, Pinheiro DP, Kamimura MT, Araújo-Filho R and Oliveira DM (2007) Flagellar proteins prediction after sequence-structure alignments of Coronin and Arp2/3 complex in Leishmania spp. In: Proceedings of the IEEE International Conference on Bioinformatics and Biomedicine - Workshop on Computational Structural Bioinformatics. IEEE Press, Piscitaway, pp 80-88. [ Links ]
Courret N, Fréhel C, Gouhier N, Pouchelet M, Prina E, Roux P and Antoine JC (2002) Biogenesis of Leishmania-harbouring parasitophorous vacuoles following phagocytosis of the metacyclic promastigote or amastigote stages of the parasites. J Cell Sci 115:2303-2316. [ Links ]
David V, Gouin E, Troys MV, Grogan A, Segal AW, Ampe C and Cossart P (1998) Identification of cofilin, coronin, Rac and capZ in actin tails using a Listeria affinity approach. J Cell Sci 111:2877-2884. [ Links ]
de Hostos EL (1999) The coronin family of actin-associated proteins. Trends Cell Biol 9:345-350. [ Links ]
de Hostos EL, Bradtke B, Lottspeich F, Guggenheim R and Gerisch G (1991) Coronin, an actin binding protein of Dictyostelium discoideum localized to cell surface projections has sequence similarities to G protein beta subunits. EMBO J 10:4097-104. [ Links ]
Dermine JF, Scianimanico S, Privé C, Descoteaux A and Desjardins M (2000) Leishmania promastigotes require lipophosphoglycan to actively modulate the fusion properties of phagosomes at an early step of phagocytosis. Cell Microbiol 2:115-126. [ Links ]
Desjardins M and Descoteaux AJ (1997) Inhibition of phagolysosomal biogenesis by the Leishmania lipophosphoglycan. Exp Med 185:2061-2068. [ Links ]
Desjardins M, Huber LA, Parton RG and Griffiths G (1994) Biogenesis of phagosomes proceeds through a sequential series of interactions with the endocytic apparatus. J Cell Biol 124:677-688. [ Links ]
Drummelsmith J, Brochu V, Girard I, Messier N and Ouellette M (2003) Proteome mapping of the protozoan parasite Leishmania and application to the study of drug targets and resistance mechanisms. Mol Cell Proteomics 2:146-155. [ Links ]
Eddy SR (1998) Profile hidden Markov models. Bioinformatics 14:755-763. [ Links ]
Edgar RC (2004) MUSCLE: A multiple sequence alignment method with reduced time and space complexity. Nucleic Acids Res 32:1792-1797. [ Links ]
Emanuelsson OH, Nielsen S, Brunak and Heijne GV (2000) Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol 300:1005-1016. [ Links ]
Ferrari G, Langen H, Naito M and Pieters J (1999) A coat protein on phagosomes involved in the intracellular survival of mycobacteria. Cell 97:435-447. [ Links ]
Forney Jr DG (1973) The Viterbi algorithm. Proceedings of the IEEE 61 3:268-278. [ Links ]
Fukui Y, Engler S, Inoue S and de Hostos EL (1999) Architectural dynamics and gene replacement of coronin suggest its role in cytokinesis. Cell Motil Cytoskel 42:204-217. [ Links ]
Garin J, Diez R, Kieffer S, Dermine J-F, Duclos S, Gagnon E, Sadoul R, Rondeau C and Desjardins M (2001) The phagosome proteome: Insight into phagosome functions. J Cell Biol 152:165-180. [ Links ]
Gotthardt D, Blancheteau V, Bosserhoff A, Ruppert T, Delorenzi M and Soldati T (2006) Proteomic fingerprinting of phagosome maturation and evidence for the role of a G alpha during uptake. Mol Cell Proteomics 5:2228-2243. [ Links ]
Gouveia JJ, Vasconcelos EJR, Pacheco ACL, Araújo-Filho R, Maia AR, Kamimura MT, Costa MP, Viana DA, Costa RB, Maggioni R, et al. (2007) Intraflagellar transport (IFT) complex in Leishmania spp: In silico genome-wide screening and annotation of gene function. Genet Mol Res 6:675-689. [ Links ]
Gueirard P, Laplante A, Rondeau C, Milon G and Desjardins M (2007) Trafficking of Leishmania donovani promastigotes in non-lytic compartments in neutrophils enables the subsequent transfer of parasites to macrophages. Cell Microbiol 10:100-111. [ Links ]
Griffiths G (2004) On phagosome individuality and membrane signalling networks. Trends Cell Biol 14:343-51. [ Links ]
Griffiths G and Mayorga L (2007) Phagosome proteomes open the way to a better understanding of phagosome function. Genome Biology 8:207. [ Links ]
Gromova I and Celis JE (2006) Protein detection in gels by silver staining: A procedure compatible with mass-spectrometry. In: Celis JE, Carter N, Hunter T, Simons K, Small JV and Shotton D (eds) Cell Biology: A Laboratory Handbook. 3rd edition. v. 4. Elsevier. Academic Press, San Diego, pp 421-429. [ Links ]
Haas A (2007) The Phagosome: Compartment with a license to kill. Traffic 8:311-330. [ Links ]
Higa RH, Togawa RC, Montagner AJ, Palandrani JC, Okimoto IK, Kuser PR, Yamagishi ME, Mancini AL and Neshich G (2004) STING Millennium Suite: Integrated software for extensive analyses of 3d structures of proteins and their complexes. BMC Bioinform 5:107-111. [ Links ]
Higgs HN and Pollard TD (2001) Regulation of actin filament network formation through ARP2/3 complex: Activation by a diverse array of proteins. Annu Rev Biochem 70:649-676. [ Links ]
Holm A, Tejle K, Magnusson KE, Descoteaux A and Rasmusson B (2001) Leishmania donovani lipophosphoglycan causes periphagosomal actin accumulation: Correlation with impaired translocation of PKCalpha and defective phagosome maturation. Cell Microbiol 3:439-447. [ Links ]
Humphries CL, Balcer HI, DAgostino JL, Winsor B, Drubin DG, Barnes G, Andrews BJ and Goode BL (2002) Direct regulation of Arp2/3 complex activity and function by the actin binding protein coronin. J Cell Biol 159:993-1004. [ Links ]
Jayachandran R, Gatfield J, Massner J, Albrecht I, Zanolari B and Pieters J (2008) RNA interference in J774 macrophages reveals a role for coronin 1 in mycobacterial trafficking but not in actin-dependent processes. Mol Biol Cell 19:1241-1251. [ Links ]
Kozminski KG, Johnson KA, Forscher P and Rosenbaum JL (1993) A motility in the eukaryotic flagellum unrelated to flagellar beating. Proc Natl Acad Sci USA 90:5519-5523. [ Links ]
Lambert C, Leonard N, De Bolle X and Depiereux E (2002) ESyPred3D: Prediction of proteins 3D structures. Bioinformatics 18:1250-1256. [ Links ]
Laskowski RA, MacArthur MW, Moss DS and Thornton JM (1993) PROCHECK: A program to check the stereochemical quality of protein structures. J Appl Cryst 26:283-291. [ Links ]
Letunic L, Copley RC, Pils B, Pinkert S, Schultz J and Bork P (2006) SMART 5: Domains in the context of genomes and networks. Nucleic Acids Res 34:257-260. [ Links ]
Loisel TP, Boujemaa R, Pantaloni D and Carlier M-F (1999) Reconstitution of actin-based motility of Listeria and Shigella using pure proteins. Nature 401:613-616. [ Links ]
Love DC, Kane MM and Mosser DM (1998) Leishmania amazonensis: The phagocytosis of amastigotes by macrophages. Exp Parasitol 88:161-171. [ Links ]
Marchler-Bauer A, Anderson JB, Cherukuri PF, DeWeese-Scott C, Geer LY, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, et al. (2005) CDD: A conserved domain database for protein classification. Nucleic Acids Res 33:D192-196. [ Links ]
Marion S, Laurent C and Guillén N (2005) Signalization and cytoskeleton activity through myosin IB during the early steps of phagocytosis in Entamoeba histolytica: A proteomic approach. Cell Microbiol 7:1504-1518. [ Links ]
Marshall WF and Rosenbaum JL (2001) Intraflagellar transport balances continuous turnover of outer doublet microtubules: Implications for flagellar length control. J Cell Biol 155:405-414. [ Links ]
Marti-Renom MA, Stuart A, Fiser A, Sánchez R, Melo F and Sali A (2000) Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct 29:291-325. [ Links ]
Méresse S, Steele-Mortimer O, Moreno E, Desjardins M, Finlay BB and Gorvel JP (1999) Controlling the maturation of pathogen-containing vacuoles: A matter of life or death. Nat Cell Biol 1:E183-E188. [ Links ]
Oku T, Itoh SM, Okano M, Suzuki A, Suzuki K, Nakajin S, Tsuji T, Nauseef WM and Toyoshima S (2003) Two regions responsible for the actin binding of p57, a mammalian coronin family actin-binding protein. Biol Pharm Bull 26:409-416. [ Links ]
Oliveira DM, Gouveia JJS, Diniz NB, Pacheco AC, Vasconcelos EJ, Diniz MC, Viana DA, Ferreira TD, Albuquerque MC, Fortier DC, et al. (2005) Pathogenomics analysis of Leishmania spp: Flagellar gene families of putative virulence factors. OMICS 9:171-191. [ Links ]
Olivier M, Gregory DJ and Forget G (2005) Subversion mechanisms by which Leishmania parasites can escape the host immune response: A signaling point of view. Clin Microbiol Rev 18:293-305. [ Links ]
Pacheco ACL, Araujo FF, Kamimura MT, Medeiros SR, Viana DA, Oliveira FCE, Araújo-Filho R, Costa MP and Oliveira DM (2007) Following the Viterbi path to deduce flagellar actin-interacting proteins of Leishmania spp.: Report on Cofilins and Twinfilins. In: Proceedings of Computational Models For Life Sciences - International Symposium on Computational Models of Life Sciences - AIP Proceedings of International Symposium on Computational Models of Life Sciences. American Institute of Physics, Melville, v. 952, pp 315-324. [ Links ]
Peters W and Killick-Kendrick R (1987) The Leishmaniases in Biology and Medicine v. II: Clinical Aspects and Control. Peters W and Killick-Kendrick R (eds) Academic Press, London, pp 551-941. [ Links ]
Puentes SM, Dwyer DM, Bates PA and Joiner KA (1989) Binding and release of C3 from Leishmania donovani promastigotes during incubation in normal human serum. J Immunol 143:3743-3749. [ Links ]
Rasmusson B and Descoteaux A (2004) Contribution of electron and confocal microscopy in the study of Leishmania-macrophage interactions. Microsc Microanal 10:656-61. [ Links ]
Rosenberger CM and Finlay BB (2003) Phagocyte sabotage: Disruption of macrophage signalling by bacterial pathogens. Nat Rev Mol Cell Biol 4:385-396. [ Links ]
Rybakin V and Clemen CS (2005) Coronin proteins as multifunctional regulators of the cytoskeleton and membrane trafficking. BioEssays 27:625-632. [ Links ]
Rybakin V, Gounko NV, Späte K, Höning S, Majoul IV, Duden R and Noegel AA (2006) Crn7 interacts with AP-1 and is required for the maintenance of Golgi morphology and protein export from the Golgi. Biol Chem 281:31070-31078. [ Links ]
Sali A and Blundell TLJ (1993) Comparative protein modelling by satisfaction of spatial restraints. Mol Biol 234:779-815. [ Links ]
Scianimanico S, Desrosiers M, Dermine JF, Meresse S, Descoteaux A and Desjardins M (1999) Impaired recruitment of the small GTPase rab7 correlates with the inhibition of phagosome maturation by Leishmania donovani promastigotes. Cell Microbiol 1:19-32. [ Links ]
Sibley LD (2003) Toxoplasma gondii: Perfecting an intracellular life style. Traffic 4:581-586. [ Links ]
Snapp EL and Landfear SMJ (1999) Characterization of a targeting motif for a flagellar membrane protein in Leishmania enriettii. Biol Chem 274:29543-29548. [ Links ]
Song L and Dentler WL (2001) Flagellar protein dynamics in Chlamydomonas. J Biol Chem 276:29754-29763. [ Links ]
Stuart LM, Boulais J, Charriere GM, Hennessy EJ, Brunet S, Jutras I, Goyette G, Rondeau C, Letarte S, Huang H, et al. (2007) A systems biology analysis of the Drosophila phagosome. Nature 445:95-101. [ Links ]
Sturgill-Koszycki S, Schaible UE and Russell DG (1996) Mycobacterium-containing phagosomes are accessible to early endosomes and reflect a transitional state in normal phagosome biogenesis. EMBO J 15:6960-696. [ Links ]
Tammana TV, Sahasrabuddhe AA, Mitra K, Bajpai VK and Gupta CM (2008) Actin-depolymerizing factor, ADF/cofilin, is essentially required in assembly of Leishmania flagellum. Mol Microbiol 70:837-852. [ Links ]
Tardieux I, Liu X, Poupel O, Parzy D, Dehoux P and Langsley G (1998) A Plasmodium falciparum novel gene encoding a coronin-like protein which associates with actin filaments. FEBS Lett 441:251-256. [ Links ]
Tull D, Vince JE, Callaghan JM, Naderer T, Spurck T, McFadden GI, Currie G, Ferguson K, Bacic A and McConville MJ (2004) SMP-1 a member of a new family of small myristoylated proteins in kinetoplastid parasites is targeted to the flagellum membrane in Leishmania. Mol Biol Cell 15:4775-4786. [ Links ]
Ullrich HJ, Beatty WL and Russell DG (1999) Direct delivery of procathepsin D to phagosomes: Implications for phagosome biogenesis and parasitism by Mycobacterium. Eur J Cell Biol 78:739-748. [ Links ]
Vasconcelos EJR, Pacheco ACL, Gouveia JJS, Araujo FF, Diniz MC, Kamimura MT, Costa MP, Maggioni R, Araujo-Filho R, Costa RB, et al. (2007) Profilins, formins and katanins as flagellar proteins of Leishmania spp: A genome-based multi-step bioinformatics-driven description. In: Bourbakis NG (ed) Conference Proceedings IEEE 7th International Symposium on Bioinformatics & Bioengineering-2007. IEEE Press, Boston, v. II, pp 880-887. [ Links ]
Vasconcelos EJR, Pacheco ACL, Gouveia JJS, Araujo FF, Diniz MC, Kamimura MT, Costa MP, Araujo-Filho R and Oliveira DM (2008) Actin-interacting proteins in flagellated pathogenic Leishmania spp.: A genome-based bioinformatics report on profilins, formins and katanins. Int J Funct Informat Personal Med 1:234-252. [ Links ]
Vieira OV, Botelho RJ and Grinstein S (2002) Phagosome maturation: Aging gracefully. Biochem J 366:689-704. [ Links ]
von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA and Bork P (2005) STRING: Known and predicted protein-protein associations integrated and transferred across organisms. Nucleic Acids Res 33:D433-D437. [ Links ]
Wiese M, Kuhn D and Grünfelder CG (2003) Protein kinase involved in flagellar-length control. Eukaryot Cell 2:769-777. [ Links ]
Yan M, Collins RF, Grinstein S and Trimble WS (2005) Coronin-1 function is required for phagosome formation. Mol Biol Cell 16:3077-87. [ Links ]
Yoshizawa AC, Kawashima S, Okuda S, Fujita M, Itoh M, Moriya Y, Hattori M and Kanehisa M (2006) Extracting sequence motifs and the phylogenetic features of SNARE-dependent membrane traffic. Traffic 7:1104-1118. [ Links ]
Black PE, "Viterbi algorithm", in Dictionary of Algorithms and Data Structures [online], Paul E. Black, ed., U.S. National Institute of Standards and Technology. 6 July 2004 (accessed 28 June 2009) Available from: http://www.itl.nist.gov/div897/sqg/dads/HTML/viterbiAlgorithm.html. [ Links ]
Send correspondence to:
Diana Magalhães de Oliveira
Núcleo Tarcísio Pimenta de Pesquisa Genômica e Bioinformática, Universidade Estadual do Ceará
Av. Paranjana 1700, Bloco UHV/NUGEN, Campus do Itaperi
60740-000 Fortaleza, CE, Brazil
Received: January 14, 2009; Accepted: May 25, 2009.
Guest Editor: José Carlos Merino Mombach