Immunogenic characterization and protective efficacy of recombinant CsgA, major subunit of curli fibers, against Vibrio parahaemolyticus.

Vibrio parahaemolyticus is one of the major pathogens responsible for vibriosis and zoonotic infections in teleosts, marine invertebrates, and also humans through consumption of contaminated or unprocessed seafood. Emergence of resistance against current accessible antibiotics and spread to the food chain and environment necessitate the development of safe and effective subunit vaccine against this bacterium. Many bacteria including V. parahaemolyticus produce extracellular curli fibrils, heteropolymeric filaments of major and minor subunit, which have been implicated in adhesion, biofilm formation, and virulence. Adhesins are the primary contact points with the host which help in establishing infection and colonization. CsgA, an adhesin, is the major structural component of the curli fiber that forms homopolymers of several hundred units. Due to their exposure on the cell surface, the curli fibers are recognized by the host's immune system, would generate high immune response, and therefore can serve as effective vaccine candidate. In the present study, we describe characterization of the csgA gene, and preparation of recombinant soluble CsgA of V. parahaemolyticus (rVpCsgA), and evaluation of its vaccine potential. Immunization of BALB/c mice with the rVpCsgA mounted a strong immune response. Cellular immune assays such as antibody isotyping, in vitro splenocyte proliferation analysis, and cytokine profiling revealed mixed T-helper cell immune response. The anti-rVpCsgA antiserum was agglutination positive and specifically cross-reacted with the curli CsgA present on the outer membrane of V. parahaemolyticus cells, thus demonstrating its neutralization potential. One hundred percent survival of the immunized mice upon challenge with the lethal dosage of the bacterium established that the rVpCsgA could serve as an effective vaccine against the bacterium. KEY POINTS: • Recombinant histidine-tagged CsgA of V. parahaemolyticus, rVpCsgA, was purified. • The rVpCsgA immunization produced mixed immune response and agglutinating antibodies. • Immunization with the rVpCsgA protected mice against V. parahaemolyticus challenge.

G-protein α-subunit (GPA1) regulates stress, nitrate and phosphate response, flavonoid biosynthesis, fruit/seed development and substantially shares GCR1 regulation in A. thaliana.

Heterotrimeric G-proteins are implicated in several plant processes, but the mechanisms of signal-response coupling and the roles of G-protein coupled receptors in general and GCR1 in particular, remain poorly understood. We isolated a knock-out mutant of the Arabidopsis G-protein α subunit (gpa1-5) and analysed its transcriptome to understand the genomewide role of GPA1 and compared it with that of our similar analysis of a GCR1 mutant (Chakraborty et al. 2015, PLoS ONE 10(2):e0117819). We found 394 GPA1-regulated genes spanning 79 biological processes, including biotic and abiotic stresses, development, flavonoid biosynthesis, transcription factors, transporters and nitrate/phosphate responses. Many of them are either unknown or unclaimed explicitly in other published gpa1 mutant transcriptome analyses. A comparison of all known GPA1-regulated genes (including the above 394) with 350 GCR1-regulated genes revealed 114 common genes. This can be best explained by GCR1-GPA1 coupling, or by convergence of their independent signaling pathways. Though the common genes in our GPA1 and GCR1 mutant datasets constitute only 26% of the GPA1-regulated and 30% of the GCR1-responsive genes, they belong to nearly half of all the processes affected in both the mutants. Thus, GCR1 and GPA1 regulate not only some common genes, but also different genes belonging to the same processes to achieve similar outcomes. Overall, we validate some known and report many hitherto unknown roles of GPA1 in plants, including agronomically important ones such as biotic stress and nutrient response, and also provide compelling genetic evidence to revisit the role of GCR1 in G-protein signalling.

Role of indirect readout mechanism in TATA box binding protein-DNA interaction.

Gene expression generally initiates from recognition of TATA-box binding protein (TBP) to the minor groove of DNA of TATA box sequence where the DNA structure is significantly different from B-DNA. We have carried out molecular dynamics simulation studies of TBP-DNA system to understand how the DNA structure alters for efficient binding. We observed rigid nature of the protein while the DNA of TATA box sequence has an inherent flexibility in terms of bending and minor groove widening. The bending analysis of the free DNA and the TBP bound DNA systems indicate presence of some similar structures. Principal coordinate ordination analysis also indicates some structural features of the protein bound and free DNA are similar. Thus we suggest that the DNA of TATA box sequence regularly oscillates between several alternate structures and the one suitable for TBP binding is induced further by the protein for proper complex formation.

Sorting of LPXTG peptides by archetypal sortase A: role of invariant substrate residues in modulating the enzyme dynamics and conformational signature of a productive substrate.

Transpeptidase sortase catalyzes the covalent anchoring of surface proteins to the cell wall in Gram-positive bacteria. Sortase A (SrtA) of Staphylococcus aureus is a prototype enzyme and considered a bona fide drug target because several substrate proteins are virulence-related and implicated in pathogenesis. Besides, SrtA also works as a versatile tool in protein engineering. Surface proteins destined for cell wall anchoring contain a LPXTG sequence located in their C-terminus which serves as a substrate recognition motif for SrtA. Recent studies have implicated substrate-induced conformational dynamics in SrtA. In the present work, we have explored the roles of invariant Leu and Pro residues of the substrate in modulating the enzyme dynamics with a view to understand the selection process of a catalytically competent substrate. Overall results of molecular dynamics simulations and experiments carried out with noncanonical substrates and site-directed mutagenesis reveal that the kinked conformation due to Pro in LPXTG is obligatory for productive binding but does not per se control the enzyme dynamics. The Leu residue of the substrate appears to play the crucial role of an anchor to the beta6-beta7 loop directing the conformational transition of the enzyme from an "open" to a "closed" state subsequent to which the Pro residue facilitates the consummation of binding through predominant engagement of the loop and catalytic motif residues in hydrophobic interactions. Collectively, our study provides insights about specificity, tolerance, and conformational sorting of substrate by SrtA. These results have important implications in designing newer substrates and inhibitors for this multifaceted enzyme.

Multivariate sequence analysis reveals additional function impacting residues in the SDR superfamily.

The "extended" type of short chain dehydrogenases/reductases (SDR), share a remarkable similarity in their tertiary structures inspite of being highly divergent in their functions and sequences. We have carried out principal component analysis (PCA) on structurally equivalent residue positions of 10 SDR families using information theoretic measures like Jensen-Shannon divergence and average shannon entropy as variables. The results classify residue positions in the SDR fold into six groups, one of which is characterized by low Shannon entropies but high Jensen-Shannon divergence against the reference family SDR1E, suggesting that these positions are responsible for the specific functional identities of individual SDR families, distinguishing them from the reference family SDR1E. Site directed mutagenesis of three residues from this group in the enzyme UDP-Galactose 4-epimerase belonging to SDR1E shows that the mutants promote the formation of NADH containing abortive complexes. Finally, molecular dynamics simulations have been used to suggest a mechanism by which the mutants interfere with the re-oxidation of NADH leading to the formation of abortive complexes.

Analysis of stacking overlap in nucleic acid structures: algorithm and application.

RNA contains different secondary structural motifs like pseudo-helices, hairpin loops, internal loops, etc. in addition to anti-parallel double helices and random coils. The secondary structures are mainly stabilized by base-pairing and stacking interactions between the planar aromatic bases. The hydrogen bonding strength and geometries of base pairs are characterized by six intra-base pair parameters. Similarly, stacking can be represented by six local doublet parameters. These dinucleotide step parameters can describe the quality of stacking between Watson-Crick base pairs very effectively. However, it is quite difficult to understand the stacking pattern for dinucleotides consisting of non canonical base pairs from these parameters. Stacking interaction is a manifestation of the interaction between two aromatic bases or base pairs and thus can be estimated best by the overlap area between the planar aromatic moieties. We have calculated base pair overlap between two consecutive base pairs as the buried van der Waals surface between them. In general, overlap values show normal distribution for the Watson-Crick base pairs in most double helices within a range from 45 to 50 Å(2) irrespective of base sequence. The dinucleotide steps with non-canonical base pairs also are seen to have high overlap value, although their twist and few other parameters are rather unusual. We have analyzed hairpin loops of different length, bulges within double helical structures and pseudo-continuous helices using our algorithm. The overlap area analyses indicate good stacking between few looped out bases especially in GNRA tetraloop, which was difficult to quantitatively characterise from analysis of the base pair or dinucleotide step parameters. This parameter is also seen to be capable to distinguish pseudo-continuous helices from kinked helix junctions.

Multiway Analysis Reveals Hydrophobicity as the Sole Determinant of Dynamic Peptide-Acetonitrile-Water Association Behavior.

Multiway analysis, a class of techniques developed for the purpose of studying multi-dimensional multivariate data, has been applied to study the dynamical structure of the first solvation layer of Ace-Gly-X-Gly-Nme peptides (where X is any amino acid) perturbed with the increase in concentrations of acetonitrile. Separate MD simulations of each peptide were carried out in five different concentrations of acetonitrile. Association of peptide, water, and acetonitrile atoms was quantified in terms of the relative abundance of Delaunay tetrahedra whose vertices could be centered on either the peptide, acetonitrile, or water atoms. A three-way data set comprising nine types of Delaunay tetrahedra in the first dimension, five concentrations of acetonitrile in the second dimension, and 26 different peptides in the third dimension was subjected to two different multiway methods viz., the constrained PARAFAC and the unconstrained Tucker3 analysis. The results unequivocally show that the dynamic peptide-acetonitrile-water association behavior could be solely explained by the hydrophobicity of the central amino acid. The study also demonstrates the utility of multiway analysis for the integration and interpretation of large number of separate MD simulations.

Immunogenic and neutralization efficacy of recombinant perfringolysin O of Clostridium perfringens and its C-terminal receptor-binding domain in a murine model.

Clostridium perfringens is a Gram-positive anaerobe ubiquitously present in different environments, including the gut of humans and animals. C. perfringens have been classified in the seven toxinotypes based on the secreted toxins that cause different diseases in humans and animals. Perfringolysin O (PFO), a cholesterol-dependent pore-forming cytolysin, is one of the potent toxins secreted by almost all C. perfringens isolates. The PFO acts in synergy with α-toxin in the progression of gas gangrene in humans and necrohemorrhagic enteritis in the calves.C. perfringens infections spread very fast, and the animals die within a few hours of the onset of infection. This necessitates the use of vaccines to control clostridial infections. Though the vaccine potential of other toxins has been reported, PFO has remained unexplored. The present study describes the immunogenic and protective potential of native recombinant PFO (WTrPFO). Since the PFO is toxic to the host cells, the non-toxic C-terminal domain of PFO (rPFOC-ter) was also assessed for its immunogenicity and protective efficacy. Immunization of mice with the purified soluble recombinant histidine-tagged WTrPFO and rPFOC-ter, expressed in E. coli, generated robust mixed immune response and T cell memory. Pre-incubation of the WTrPFO with anti-WTrPFO and rPFOC-ter antisera negated its hemolytic activity in mice RBCs, as well as its cytotoxic effect in mice peritoneal macrophages in vitro. Thus, immunization with the WTrPFO and its non-toxic C-terminal domain generated neutralizing antibodies, suggesting their vaccine potential against the PFO. Thus, the non-toxic C-terminal domain of PFO could serve as an alternative to PFO as a vaccine candidate.

Recombinant FimH, a fimbrial tip adhesin of Vibrio parahaemolyticus, elicits mixed T helper cell response and confers protection against Vibrio parahaemolyticus challenge in murine model.

Vibrio parahaemolyticus causes vibriosis in wide range of marine organisms, and is responsible for food borne illnesses in humans through consumption of contaminated uncooked/partially cooked seafood. Continued and widespread antibiotics usage to increase the productivity has led to antibiotics resistance development. This has necessitated the need to develop alternative methods to control its infection. Use of safe and effective vaccines against the virulence factors not only protects from infection, it also minimizes antibiotic usage. The colonization of V. parahaemolyticus in the host and disease development requires several adhesins present on the cell surface, and thereby make them attractive vaccine candidates. V. parahaemolyticus produces extracellular type 1 fimbriae that have been shown to play a role in adhesion, biofilm formation and virulence. FimH is one of the minor components of the type 1 fimbriae occurring on its very tip. Being present on the cell surface, it is highly immunogenic, and can be targeted as a potential vaccine candidate. The present study describes the immunogenic and vaccine potential of recombinant V. parahaemolyticus FimH (rVpFimH) expressed in E. coli. Immunization of BALB/c mice with the rVpFimH elicited a strong mixed immune response, T-cell memory (evidenced by antibody isotyping, cytokine profiling and T-cell proliferation assay), and agglutination positive antibodies. FACS analysis and immunogold labeling showed that the polyclonal anti-rVpFimH antibodies were able to recognize the FimH on V. parahaemolyticus cells. In vivo challenge of the rVpFimH-immunized mice with 2×LD50 dose of live bacteria showed one hundred percent survival. Thus, our findings clearly demonstrate the potential of FimH as an effective vaccine candidate against V. parahaemolyticus.

Genomewide bioinformatic analysis negates any specific role for Dof, GATA and Ag/cTCA motifs in nitrate responsive gene expression in Arabidopsis.

Nitrate response at the plant level is mediated by the transcriptional regulation of several hundreds of genes, but no common cis-acting nitrate-responsive elements (NREs) have been identified so far. Earlier, we bioinformatically ruled out the possibility that the previously published [(a/t)7Ag/cTCA] motif could act as NRE on its own (Das et al., 2007, Mol. Genet. Genomics, 278: 519-525). In the present study, we examined other motifs such as Dof and GATA binding elements in homologous as well as heterologous pairwise combinations in the Arabidopsis genome in silico. None of the above three motifs revealed any unique association with nitrate responsive genes or their subsets in any combination, either within their ORFs or 1 kb flanking sequences on either side. Additionally, twelve new, top-scoring candidate motifs that were generated using different online motif samplers were analyzed in silico using a subset of 21 'early' nitrate responsive genes, but did not reveal any specificity of occurence. These results underscore the need to continue the search for novel candidate NREs, as possible sites of intervention to understand/improve nitrate-responsive gene expression and nitrate use efficiency.

Conformational perturbation of peptides in presence of polar organic solvents.

The critical role played by solvent environment in maintaining the conformational integrity of peptides and proteins is accepted without question. Numerous experiments have suggested that perturbing the solvent environment of peptides and proteins by the addition of polar organic solvents have important consequences for the conformation of these molecules. However, experimental studies of such perturbations often report different kinds of effects depending on the solvent used and/or the sequence/structure of the molecule under study. In this work we report a simulation based comparative study on the effects of adding two common organic solvents viz. Dimethyl sulfoxide (DMSO) and Acetonitrile (MeCN) on the dynamical conformation of a test peptide Ace-Gly-X-Gly-Nme where X is any amino acid. Our studies identify important differences in peptide solvation by these two solvents, which we attempt to correlate with the kinetic stability of the conformation, as well as the identity of the central 'X' residue in the test peptide.

Role of Ser65, His148 and Thr203 in the Organic Solvent-dependent Spectral Shift in Green Fluorescent Protein.

The photophysics of green fluorescent protein (GFP) is remarkable because of its exceptional property of excited state proton transfer (ESPT) and the presence of a functional proton wire. Another interesting property of wild-type GFP is that its absorption and fluorescence excitation spectra are sensitive to the presence of polar organic solvents even at very low concentrations. Here, we use a combination of methodologies including site-specific mutagenesis, absorption spectroscopy, steady-state and time-resolved fluorescence measurements and all-atom molecular dynamics simulations in explicit solvent, to uncover the mechanism behind the unique spectral sensitivity of GFP toward organic solvents. Based on the evidences provided herein, we suggest that organic solvent-induced changes in the proton wire prevent ground state movement of a proton through the wire and thus bring about the spectral changes observed. The present study can not only help to understand the mechanism of proton transfer by further dissecting the intricate steps in GFP photophysics but also encourages to develop GFP-based organic solvent biosensors.

MFS transportome of the human pathogenic yeast Candida albicans.

BACKGROUND: The major facilitator superfamily (MFS) is one of the two largest superfamilies of membrane transporters present ubiquitously in bacteria, archaea, and eukarya and includes members that function as uniporters, symporters or antiporters. We report here the complete transportome of MFS proteins of a human pathogenic yeast Candida albicans. RESULTS: Computational analysis of C. albicans genome enabled us to identify 95 potential MFS proteins which clustered into 17 families using Saier's Transport Commission (TC) system. Among these SP, DHA1, DHA2 and ACS represented major families consisting of 22, 22, 9 and 16 members, respectively. Family designations in C. albicans were validated by subjecting Saccharomyces cerevisiae genome to TC system. Based on the published available genomics/proteomics data, 87 of the putative MFS genes of C. albicans were found to express either at mRNA or protein levels. We checked the expression of the remaining 8 genes by using RT-PCR and observed that they are not expressed under basal growth conditions implying that either these 8 genes are expressed under specific growth conditions or they may be candidates for pseudogenes. CONCLUSION: The in silico characterisation of MFS transporters in Candida albicans genome revealed a large complement of MFS transporters with most of them showing expression. Considering the clinical relevance of C. albicans and role of MFS members in antifungal resistance and nutrient transport, this analysis would pave way for identifying their physiological relevance.

A statistical geometry analysis of simulated water-DMSO and water-MeCN binary mixtures for biomolecular studies.

Water-Dimethylsulfoxide (DMSO) and water-Acetonitrile (MeCN) binary mixtures at various molar ratios ranging from 0 to 1 are studied using Molecular Dynamics (MD) simulations. Hydration properties of water in different regions of MeCN/DMSO are investigated by using the statistical geometry approach. The obtained results reveal that in water-DMSO simulations both water and solvent molecules prefer to be in mixed cluster forms, depending upon the concentration of DMSO. While in case of water-MeCN mixtures, self-association of water and acetonitrile molecules, take place, showing microheterogeneity associated with the water- MeCN binary mixtures. The results highlight the utility of statistical geometric analysis of MD simulation data of binary liquid mixtures for rapid screening of polar organic solvents in non-aqueous enzymology.

Modification of axial fiber contact residues impact sickle hemoglobin polymerization by perturbing a network of coupled interactions.

The identity of intermolecular contact residues in sickle hemoglobin (HbS) fiber is largely known. However, our knowledge about combinatorial effects of two or more contact sites or the mechanistic basis of such effects is rather limited. Lys16, His20, and Glu23 of the alpha-chain occur in intra-double strand axial contacts in the sickle hemoglobin (HbS) fiber. Here we have constructed two novel double mutants, HbS (K16Q/E23Q) and (H20Q/E23Q), with a view to delineate cumulative impact of interactions emanating from the above contact sites. Far-UV and visible region CD spectra of the double mutants were similar to the native HbS indicating the presence of native-like secondary and tertiary structure in the mutants. The quaternary structures in both the mutants were also preserved as judged by the derivative UV spectra of liganded (oxy) and unliganded (deoxy) forms of the double mutants. However, the double mutants displayed interesting polymerization behavior. The polymerization behaviour of the double mutants was found to be non-additive of the individual single mutants. While HbS (H20Q/E23Q) showed inhibitory effect similar to that of HbS (E23Q), the intrinsic inhibitory propensity of the associated single mutants was totally quelled in HbS (K16Q/E23Q) double mutant. Molecular dynamics (MD) simulations studies of the isolated alpha-chains as well as a module of the fiber containing the double and associated single mutants suggested that these contact sites at the axial interface of the fiber impact HbS polymerization through a coupled interaction network. The overall results demonstrate a subtle role of dynamics and electrostatics in the polymer formation and provide insights about interaction-linkage in HbS fiber assembly.

Transcriptome analysis of Arabidopsis GCR1 mutant reveals its roles in stress, hormones, secondary metabolism and phosphate starvation.

The controversy over the existence or the need for G-protein coupled receptors (GPCRs) in plant G-protein signalling has overshadowed a more fundamental quest for the role of AtGCR1, the most studied and often considered the best candidate for GPCR in plants. Our whole transcriptome microarray analysis of the GCR1-knock-out mutant (gcr1-5) in Arabidopsis thaliana revealed 350 differentially expressed genes spanning all chromosomes. Many of them were hitherto unknown in the context of GCR1 or G-protein signalling, such as in phosphate starvation, storage compound and fatty acid biosynthesis, cell fate, etc. We also found some GCR1-responsive genes/processes that are reported to be regulated by heterotrimeric G-proteins, such as biotic and abiotic stress, hormone response and secondary metabolism. Thus, GCR1 could have G-protein-mediated as well as independent roles and regardless of whether it works as a GPCR, further analysis of the organism-wide role of GCR1 has a significance of its own.

Rational design, synthesis, and verification of affinity ligands to a protein surface cleft.

The structure-based design, synthesis, and screening of a glucuronic acid scaffold library of affinity ligands directed toward the catalytic cleft on porcine pancreas alpha-amylase are presented. The design was based on the simulated docking to the enzyme active site of 53 aryl glycosides from the Available Chemicals Directory (ACD) selected by in silico screening. Twenty-three compounds were selected for synthesis and screened in solution for binding toward alpha-amylase using nuclear magnetic resonance techniques. The designed molecules include a handle outside of the binding site to allow their attachment to various surfaces with minimal loss of binding activity. After initial screening in solution, one affinity ligand was selected, immobilized to Sepharose (Amersham Biosciences), and evaluated as a chromatographic probe. A column packed with ligand-coupled Sepharose specifically retained the enzyme, which could be eluted by a known inhibitor.

Genomewide computational analysis of nitrate response elements in rice and Arabidopsis.

Nitrate response element (NRE) was originally reported to be comprised of an Ag/cTCA core sequence motif preceded by a 7-bp AT rich region, based on promoter deletion analyses in nitrate and nitrite reductases from Arabidopsis thaliana and birch. In view of hundreds of new nitrate responsive genes discovered recently, we sought to computationally verify whether the above motif indeed qualifies to be the cis-acting NRE for all the responsive genes. We searched for the specific occurrence of at least two copies of the above motif in and around the nitrate responsive genes and elsewhere in the Arabidopsis and rice (Oryza sativa) genomes, with respect to their positional, orientational and strand-specific bias. This is the first comprehensive analysis of NREs for 625 nitrate responsive genes of Arabidopsis and their rice homologs, representing dicots and monocots, respectively. We report that the above motifs are present almost randomly throughout these genomes and do not reveal any specificity or bias towards nitrate responsive genes. This also seems to be true for smaller subsets of nitrate responsive genes in Arabidopsis, such as the 21 early responsive genes, 261 and 90 genes for root-specific and shoot-specific response, respectively, and 25 housekeeping genes. This necessitates a fresh search for candidate sequences that qualify to be NREs in these and other plants.

The affinity of the FimH fimbrial adhesin is receptor-driven and quasi-independent of Escherichia coli pathotypes.

Type-1 fimbriae are important virulence factors for the establishment of Escherichia coli urinary tract infections. Bacterial adhesion to the high-mannosylated uroplakin Ia glycoprotein receptors of bladder epithelium is mediated by the FimH adhesin. Previous studies have attributed differences in mannose-sensitive adhesion phenotypes between faecal and uropathogenic E. coli to sequence variation in the FimH receptor-binding domain. We find that FimH variants from uropathogenic, faecal and enterohaemorrhagic isolates express the same specificities and affinities for high-mannose structures. The only exceptions are FimHs from O157 strains that carry a mutation (Asn135Lys) in the mannose-binding pocket that abolishes all binding. A high-mannose microarray shows that all substructures are bound by FimH and that the largest oligomannose is not necessarily the best binder. Affinity measurements demonstrate a strong preference towards oligomannosides exposing Manalpha1-3Man at their non-reducing end. Binding is further enhanced by the beta1-4-linkage to GlcNAc, where binding is 100-fold better than that of alpha-d-mannose. Manalpha1-3Manbeta1-4GlcNAc, a major oligosaccharide present in the urine of alpha-mannosidosis patients, thus constitutes a well-defined FimH epitope. Differences in affinities for high-mannose structures are at least 10-fold larger than differences in numbers of adherent bacteria between faecal and uropathogenic strains. Our results imply that the carbohydrate expression profile of targeted host tissues and of natural inhibitors in urine, such as Tamm-Horsfall protein, are stronger determinants of adhesion than FimH variation.

Complete inventory of ABC proteins in human pathogenic yeast, Candida albicans.

The recent completion of the sequencing project of the opportunistic human pathogenic yeast, Candida albicans (http://www.ncbi.nlm.nih.gov/), led us to analyze and classify its ATP-binding cassette (ABC) proteins, which constitute one of the largest superfamilies of proteins. Some of its members are multidrug transporters responsible for the commonly encountered problem of antifungal resistance. TBLASTN searches together with domain analysis identified 81 nucleotide-binding domains, which belong to 51 different putative open reading frames. Considering that each allelic pair represents a single ABC protein of the Candida genome, the total number of putative members of this superfamily is 28. Domain organization, sequence-based analysis and self-organizing map-based clustering led to the classification of Candida ABC proteins into 6 distinct subfamilies. Each subfamily from C. albicans has an equivalent in Saccharomyces cerevisiae suggesting a close evolutionary relationship between the two yeasts. Our searches also led to the identification of a new motif to each subfamily in Candida that could be used to identify sequences from the corresponding subfamily in other organisms. It is hoped that the inventory of Candida ABC transporters thus created will provide new insights into the role of ABC proteins in antifungal resistance as well as help in the functional characterization of the superfamily of these proteins.