# The Conserved Candida albicans CA3427 Gene Product Defines a New Family of Proteins Exhibiting the Generic Periplasmic Binding Protein Structural Fold

Sébastien Santini, Jean-Michel Claverie, Nicolas Mouz, Tristan Rousselle, Caroline Maza, Vincent Monchois, Chantal Abergel

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Sébastien Santini, Jean-Michel Claverie, Nicolas Mouz, Tristan Rousselle, Caroline Maza, et al. The Conserved Candida albicans CA3427 Gene Product Defines a New Family of Proteins Exhibiting the Generic Periplasmic Binding Protein Structural Fold. PLoS ONE, 2011, 6 (4), pp.e18528. 10.1371/journal.pone.0018528. hal-05370973

HAL Id: hal-05370973

https://amu.hal.science/hal-05370973v1

Submitted on 18 Nov 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The Conserved *Candida albicans* CA3427 Gene Product Defines a New Family of Proteins Exhibiting the Generic Periplasmic Binding Protein Structural Fold

Sébastien Santini¹, Jean-Michel Claverie¹, Nicolas Mouz², Tristan Rousselle², Caroline Maza¹, Vincent Monchois², Chantal Abergel¹*

¹Information Génomique et Structurale (CNRS UPR2589), Aix-Marseille Université, Mediterranean Institute of Microbiology, Parc Scientifique de Luminy, Marseille, France,

²PX'Therapeutics - 7 Parvis Louis Néel - BP50, Grenoble, France

## Abstract

Nosocomial diseases due to *Candida albicans* infections are in constant rise in hospitals, where they cause serious complications to already fragile intensive care patients. Antifungal drug resistance is fast becoming a serious issue due to the emergence of strains resistant to currently available antifungal agents. Thus the urgency to identify new potential protein targets, the function and structure of which may guide the development of new antifungal drugs. In this context, we initiated a comparative genomics study in search of promising protein coding genes among the most conserved ones in reference fungal genomes. The CA3427 gene was selected on the basis of its presence among pathogenic fungi contrasting with its absence in the non pathogenic *Saccharomyces cerevisiae*. We report the crystal 3D-structure of the *Candida albicans* CA3427 protein at 2.1 Å resolution. The combined analysis of its sequence and structure reveals a structural fold originally associated with periplasmic binding proteins. The CA3427 structure highlights a binding site located between the two protein domains, corresponding to a sequence segment conserved among fungi. Two crystal forms of CA3427 were found, suggesting that the presence or absence of a ligand at the proposed binding site might trigger a "Venus flytrap" motion, coupled to the previously described activity of bacterial periplasmic binding proteins. The conserved binding site defines a new subfamily of periplasmic binding proteins also found in many bacteria of the bacteroidetes division, in a choanoflagellate (a free-living unicellular and colonial flagellate eukaryote) and in a placozoan (the closest multicellular relative of animals). A phylogenetic analysis suggests that this gene family originated in bacteria before its horizontal transfer to an ancestral eukaryote prior to the radiation of fungi. It was then lost by the Saccharomycetales which include *Saccharomyces cerevisiae*.

**Citation:** Santini S, Claverie J-M, Mouz N, Rousselle T, Maza C, et al. (2011) The Conserved *Candida albicans* CA3427 Gene Product Defines a New Family of Proteins Exhibiting the Generic Periplasmic Binding Protein Structural Fold. PLoS ONE 6(4): e18528. doi:10.1371/journal.pone.0018528

**Editor:** Jörg D. Hoheisel, Deutsches Krebsforschungszentrum, Germany

**Received** January 24, 2011; **Accepted** March 11, 2011; **Published** April 11, 2011

**Copyright:** © 2011 Santini et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** The PROFUN project was supported by a grant from the French ministry of research and education and by Centre National de la Recherche Scientifique (CNRS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** All authors have declared that no competing interests exist, including TR, VM and NM who are employees of Protein'eXpert. There are no other financial and non-financial competing interests. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials.

* E-mail: chantal.abergel@igs.cnrs-mrs.fr

## Introduction

*Candida* spp are ubiquitous commensal organisms that can cause serious disseminated infections, particularly in immunocompromised and intensive care patients. *Candida* spp. are the fourth leading cause of nosocomial bloodstream infections in the United States, with treatment costs estimated to be more than $2–$4 billion annually [1] and with mortality rates estimated between 38% to 49% [2]. Candidiasis is the most common invasive fungal infection reported in cancer patients (58%–69%) [3], and over the past decade, the incidence of these fungal infections has increased significantly [4]. Although resistance to antifungal drugs remains uncommon on community acquired infections, they are in constant rise in nosocomial infections [5]. Since it has been demonstrated that clinical isolates of the *Candida* species *C. albicans*, *C. glabrata*, *C. tropicalis*, and *C. krusei* have acquired resistance against first-line agents for treatment of invasive candidiasis by mutations in the gene encoding the target enzyme (glucan synthase) [5,6,7], it appears important to anticipate and enlarge

the antifungal drug spectrum by identifying new original targets. In this context, our laboratory led a prospective structural genomics project (PROFUN [8]) in search of new antifungal targets.

The gene selection was based on the comparison of the following fungi genomes: *Candida albicans* (SC5314), *Saccharomyces cerevisiae* (S288C), *Neurospora crassa* (OR74A), *Magnaporthe grisea* (70-15 (Mat1-1)), *Schizosaccharomyces pombe*, *Aspergillus fumigatus* (Af293), *Phanerochaete chrysosporium*, *Cryptococcus neoformans* ((Serotype D) JEC21+B3501), *Cryptococcus neoformans* ((Serotype A) H99). This study aimed to identify virulence-related targets by focusing on genes conserved in pathogenic fungi and absent from the *Saccharomyces cerevisiae* genome. The CA3427 gene belongs to this category and encodes a 299 amino acid-long, 33.7 kDa molecular weight protein of unknown function (UNIPROT: Q59X88).

The comparison of the CA3427 sequence with its database homologs clearly highlights a new functional family conserved (>30% identity over its entire length) across most fungi genomes

<figure><img src="image_3.png" /></figure>
and present in some flavobacteria. It only shares a weak similarity (<20% identity over the full length sequences) with the Pyrimidine precursor biosynthesis THI13 enzyme from *S. cerevisiae*. To gain insights into the function and druggability of CA3427 we determined its crystal structure by the multi-wavelength anomalous dispersion (MAD) method [9]. Interestingly, two crystal forms were obtained which seem to correspond to a large conformational change induced by the binding of a small ligand at a specific site of the protein.

## Methods

### Cloning expression and purification

As part of the larger structural genomics PROFUN project, CA3427 was produced using the protocol previously described for other targets [10]. Briefly, the cDNA was isolated by PCR using sequence primers specific to the CA3427 gene preceded by 5'-CATCACCATCAATTG (Direct primer) and 5'-TCACCATCCAATTG (Reverse Primer) applied to a template of purified genomic DNA from the *Candida albicans* strain NIH 3147 (ATCC number MYA-2876D). Gene cloning was performed using the ligation-independent cloning (LIC) method and our pSF-04 expression vector [10]. The PCR products were directly purified using the NucleoSpin Extract kit (Macherey Nagel). Then, 0.2 pmol of the purified PCR product was treated with T4 DNA polymerase in the presence of 2.5 mM of dCTP for 30 minutes at 22°C before inactivating the enzyme (20 minutes at 75°C). In a parallel procedure, the pSF-04 expression vector was digested with the *Mfel* restriction enzyme to excise the insert bearing the lacZ encoding sequence. pSF-04 was then purified on agarose gel using the NucleoSpin Extract kit (Macherey Nagel) and treated with T4 DNA polymerase in the presence of 2.5 mM of dGTP for 30 minutes at 22°C before inactivating the enzyme (20 minutes at 75°C).

The CA3427 cloning was performed as follows. A hybridization reaction was carried out by mixing 0.01 pmol of pSF-04 and 0.02 pmol of the insert in a reaction volume of 3 µl, followed by a 5 minutes incubation at 22°C and the subsequent addition of 1 µl of 25 mM EDTA. After a second incubation of 5 minutes at 22°C, the resulting product was used to transform *E. coli* DH5α. Transformants were selected on LB plates containing 100 µg/ml ampicillin, and positive colonies were isolated. This cloning procedure allowed the addition of a (His)₆ tag followed by the GHHHQL sequence to the N-terminal of the CA3427 gene product and of a C-terminal QLDGDLLEAA linker to the GFP protein.

An expression screen was then performed using our standard procedure [11]. The GFP reporter was used to quantify (and determine the optimal condition for) the soluble expression of the CA3427 protein through fluorescence measurements [12]. The subsequent removal of the GFP-encoding gene was done by *NolI* digestion followed by the circularization of the plasmid.

The plasmid born CA3427 gene was over-expressed in *E. coli* BL21 in 1L flasks containing TB medium over one night at 17°C after induction with IPTG (500 µM) at OD₆₀₀ nm = 0.5. The selenomethionine-substituted protein was produced using the appropriate protocol to inhibit methionine synthesis in the presence of selenomethionine and M9 minimal medium [13].

After centrifugation, the pellet was resuspended in buffer A (50 mM NaH₂PO₄, 300 mM NaCl pH 8.0) with 5% glycerol and 0.1% Triton X-100 then sonicated and centrifuged again.

The cleared lysate was applied to a 5 ml HiTrap Chelating Column (GE Healthcare) charged with Ni²⁺ and equilibrated with buffer A. The column was washed with 10 column volumes of buffer A, 10 column volumes of buffer A containing 25 mM

Imidazole and 5 column volumes of buffer A containing 50 mM Imidazole at a flow rate of 1 ml.min⁻¹. Elution was performed with a linear gradient over 7 column volumes from 50 mM to 500 mM Imidazole. The fractions corresponding to the elution of CA3427 with 150–200 mM Imidazole were run on a desalting column (Fast Desalting Column HR 10/10, Pharmacia) and we controlled the recombinant protein sequence by mass spectroscopy and N-terminal Edman sequencing. After purification, the fractions contained at least 98% pure protein in 10 mM Tris buffer at pH 7. Upon isoelectric focusing chromatography, the recombinant CA3427 protein showed a PI≈5. The analysis by dynamic light scattering of the purified recombinant CA3427 protein indicated a monodisperse solution with a gyration radius of ≈2.5 nm, compatible with a monomer.

### Crystallization

The *C. albicans* CA3427 recombinant protein was concentrated to 18.5 g/L in 10 mM MOPS buffer at pH 7.5 using a centrifugal filter device (Ultrafree Biomax 10K, Millipore, Bedford MA, USA). The screening for crystallization conditions was performed using a standard strategy [10].

The best crystals were obtained using the hanging drop vapor diffusion method with a 1 ml reservoir. Crystallization droplets were made of 0.5 µl of protein mixed with 0.5 µl of the reservoir solution made of 21% PEG8000, 0.2 M Calcium Acetate, 0.1 M Tris, 30% Glycerol at a pH of 7.0 (structure 1) and 13% PEG8000, 0.2 M Calcium Acetate, 0.1 M Tris, 10% Glycerol at a pH of 7.0 (structure 2). Crystals appeared within a few days.

To explore the CA3427 specificity, we performed co-crystallization experiments with a variety of ligands at a concentration of 1 mM (pyridoxal phosphate, histidine, lysine, arginine, Glutamine, Leucine, Isoleucine). None of them resulted in a liganded structure with extra electron density in the CA3427 binding site.

### Data collection

Crystals of the CA3427 protein were mounted in a Hampton Research 0.2 mm³ loop, flash frozen to 100K in a cold nitrogen gas stream and subjected to X-rays. The two datasets were collected on a MarCCD (165 mm) camera at the European Synchrotron Radiation Facility (ESRF) on the BM30A-FIP beamline.

The first *C. albicans* CA3427 structure (PDB: 2X7P) was determined using the MAD method based on a two-wavelength data set (Table 1) obtained with a selenomethionine-substituted protein crystal. The crystals belong to the orthorhombic space group P2₁2₁2₁, with unit cell parameters *a* = 42.588 Å, *b* = 66.849 Å, *c* = 113.990 Å, α = β = γ = 90°.

The second dataset (PDB: 2X7Q) was collected at a wavelength of 0.975627 Å. The crystals belong to the P2₁2₁2₁ space group with unit cell parameters *a* = 41.411, *b* = 65.724, *c* = 128.203, α = β = γ = 90°.

### Structures determination and refinement

The diffraction data were indexed with MOSFLM [14] and scaled with the SCALA [15] software from the CCP4 suite [16].

Phase determination was performed by using the SOLVE program [17] on two wavelengths corresponding to the peak (0.979774 Å) and the inflexion point (0.979958 Å) in the 43.437 to 2.341 Å resolution range. A single solution was found with a mean figure of merit of 0.4 for all the data between 35 and 2.5 Å. The phases obtained were improved by using autoSHARP [18]. The electron-density map was used to construct the main chain of the molecules by using COOT [19]. Refinement was performed using the Phenix software [20] including manual rebuilding and rigid body refinement followed by several cycles of positional refinement.

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1.** X-ray data collection (ESRF) and refinement statistics.</caption><thead><tr><th></th><th>2X7P</th><th>2X7Q</th></tr></thead><tbody><tr><td>Beam line</td><td>BM30A</td><td>ID29</td></tr><tr><td>Method</td><td>MAD</td><td>Molecular Replacement</td></tr><tr><td>Wavelength (Å)</td><td>0.979774</td><td>0.975627</td></tr><tr><td>Space group</td><td>P 21 21 21</td><td>P 21 21 21</td></tr><tr><td>Unit cell (Å)</td><td>a=42.588 b=66.849 c=113.990 α=90 β=90 γ=90</td><td>a=41.411 b=65.724 c=128.203 α=90 β=90 γ=90</td></tr><tr><td>Resolution range (Å)</td><td>43.337 to 2.341</td><td>46.127 to 2.0</td></tr><tr><td>(highest resolution shell)</td><td>(2.62 to 2.51)</td><td>(2.11 to 2.0)</td></tr><tr><td>Observations</td><td>61514 (1615)</td><td>147779 (16688)</td></tr><tr><td>Unique reflections</td><td>10802 (534)</td><td>17895 (2314)</td></tr><tr><td>Multiplicity<sup>1</sup></td><td>5.5 (2.4)</td><td>8.0 (6.8)</td></tr><tr><td>Completeness<sup>1</sup></td><td>98.8 (89.6)</td><td>76.7 (72.6)</td></tr><tr><td>&lt;I/σI&gt;<sup>1,2</sup></td><td>7.6 (3.9)</td><td>6.5 (2.7)</td></tr><tr><td>R<sub>sym</sub> (%)<sup>1,3</sup></td><td>6 (16.5)</td><td>9.6 (26.5)</td></tr><tr><td><strong>Refinement</strong></td><td></td><td></td></tr><tr><td>R<sub>cryst</sub> (%)<sup>4</sup></td><td>20.9</td><td>21.8</td></tr><tr><td>R<sub>free</sub> (%)</td><td>25.6</td><td>23.6</td></tr><tr><td>Δ<sub>bond</sub> (Å)</td><td>0.001</td><td>0.002</td></tr><tr><td>Δ<sub>angle</sub> (°)</td><td>0.394</td><td>0.478</td></tr><tr><td>N° Protein atoms</td><td>2463</td><td>2466</td></tr><tr><td>N° water</td><td>140</td><td>220</td></tr><tr><td>N° Heterogen atoms</td><td>46</td><td>8</td></tr><tr><td><strong>Average B factor (Å<sup>2</sup>)</strong></td><td></td><td></td></tr><tr><td>All atoms</td><td>22.9</td><td>14.3</td></tr><tr><td>Protein main chain</td><td>22.2</td><td>13.7</td></tr><tr><td>Water</td><td>26</td><td>18.1</td></tr><tr><td>Ligand</td><td>28</td><td>18.8</td></tr><tr><td><strong>Ramachandran plot (%)</strong></td><td></td><td></td></tr><tr><td>Most favored</td><td>261</td><td>258</td></tr><tr><td>Allowed</td><td>18</td><td>21</td></tr><tr><td>Generously allowed</td><td>0</td><td>0</td></tr><tr><td>Disallowed regions</td><td>1</td><td>1</td></tr></tbody></table>

For the second crystal form, we used molecular replacement on the CaspR server [21] and the MAD-solved three-dimensional structure of CA3427 as template. The structure was refined using COOT and iterative steps of manual rebuilding and positional refinement using Phenix. PROCHECK [22] was used to assess the quality of the structures. All statistics are presented in Table 1.

The atomic coordinates and structure factors for the crystal structures of the CA3427 protein from *Candida albicans* are available in the RCSB Protein Data Bank under PDB id 2X7P and 2X7Q.

## Phylogenetic analysis

The evolutionary relationship of CA3427 with its homologs was assessed as follows. We searched for orthologous sequences against the 82 available reference fungi genomes [23]. All BLAST [24]

searches were performed on the servers hosting the corresponding fungi genomes with default parameters: BROAD Institute Fungal Genome Initiative [25], Department of Energy Joint Genome Institute [26], National Center for Biotechnology Information [27], Resources for Fungal Comparative Genomics [28] and Fungal Genome Research website [29]. Only 60 species showed an unambiguous homolog which were used to build the phylogenetic tree and compare it with a reference tree. To collect a larger panel of homologous sequences, the CA3427 sequence was used as a seed for BLAST search against the Ref-Seq database (NCBI). Sequences of bacterial origin, all belonging to the bacteroidetes/Flavobacteria clade, were readily identified as best matching the CA3427 protein sequence (E value<10<sup>-20</sup>), as well the two additional unexpected homologs of eukaryotic origin, one from the choanoflagellate *Monosiga brevicollis* and the other one

<figure><img src="image_5.png" /></figure>

<footnote n="1">¹values in parentheses are for the highest resolution shell.</footnote>

<footnote n="2">²<i>l</i> is the mean signal to noise ratio, where <i>l</i> is the integrated intensity of a measured reflection and σ is the estimated error in the measurement.</footnote>

<footnote n="3">³R<sub>sym</sub> = Σ<sub>h</sub> Σ<sub>i</sub> |I<sub>h,i</sub> - <i>I</i>|/ Σ<sub>h</sub> Σ<sub>i</sub> |I<sub>h,i</sub>|, where <i>I</i> is the integrated intensity of reflection <i>h</i> having <i>i</i> observations and <i>I</i> is the mean recorded intensity of reflection <i>h</i> over multiple recording.</footnote>

<footnote n="4">⁴R<sub>cryst</sub> = Σ<sub>h</sub> ||F<sub>o</sub> - |F<sub>c</sub>|/ Σ<sub>h</sub> |F<sub>o</sub>|, where F<sub>o</sub> are observed and F<sub>c</sub> calculated structure factor amplitudes. R<sub>free</sub> is calculated from a randomly chosen 9.9% of reflections.</footnote>

<footnote>doi:10.1371/journal.pone.0018528.t001</footnote>
from the placozoan *Trichoplax adhaerens*. All those sequences turned out to respect the residue conservation previously identified for the fungal CA3427 homologs strengthening the definition of a new periplasmic binding protein (PBP)-like subfamily. We then selected 15 non redundant representatives of the fungi that were retained for a detailed phylogenetic analysis together with 9 bacterial sequences as well as the choanoflagellate and the placozoan sequences. To identify a suitable outgroup of PBP-related sequences for rooting purpose, we selected more divergent, yet highly significant (E value<10<sup>-6</sup>), bacterial homolog sequences from a cyanobacterium, a firmicute, a beta and an alpha proteobacterium, none of them sharing the new subfamily signature. To optimize the multiple alignment, proteins were truncated at the domain boundaries of the CA3427 sequence. This dataset was used to study the evolutionary relationship within this new PBP-like subfamily using the phylogeny.fr web server [30]. Details of the parameter used for the computation are provided in the figures legends.

## Results and Discussion

### Overall structure of CA3427

CA3427 is an α/β protein with two domains organized into a C-clamp shape (Fig. 1A–B). Domain I, encompassing residues 1 to 81 and 190 to 299, is composed of a 5 stranded β-sheet (β₂ β₁ β₃ β₁₀ β₄) with β₁₀ anti-parallel to the others, surrounded by 10 helices (α₁ to α₄ and α₈ to α₁₃). The smaller domain II, encompassing residues 88 to 183, is also arranged in a 5 stranded β-sheet (β₇ β₆ β₈ β₅ β₉) with β₉ anti-parallel to the others, surrounded by 3 α-helices (α₅ α₆ α₇). The two domains delimit a large groove and are linked by a hinge region (residues 82–87 and 184–189).

### Comparison of the two CA3427 structures

The two crystal forms correspond to distinct conformations of the CA3427 protein with a root mean square deviation (RMSD) of

1.55 Å based on α-carbons superimposition of the overall structures. Most of the RMSD value results from a change in the relative position of the two domains rather than from local rearrangements (Fig. 1C, see Movie S1 for an animated view). This is well demonstrated by measuring the distances between three α-carbons delimiting the groove: E40, S140 and D236. The distances separating E40 from S140 varies from 12.44 Å to 8.41 Å between the two structures, and from 20.93 Å to 18.33 Å for the distance between D236 and S140.

To determine the motion best describing the transition between the two crystal structures, we performed a normal mode analysis on the El Némo server [31] using the two structures. The C-terminal tag was truncated in order to avoid irrelevant motions. Normal modes were computed on one structure and for each mode, we computed the RMSD of each model fitted onto the alternative structure. The first 6 modes corresponding to self rotations and translations applied to the whole system, were not taken into account further. The lowest RMSD value (1.029) was found to correspond to normal mode number 7, exhibiting a small torsion and a closure of the two domains. It clearly corresponds to a clamp motion, also known as a “Venus flytrap” motion, folding the two domains onto each other using the flexibility of the hinge region [32].

### Analysis of the CA3427 putative binding site

The two CA3427 structures exhibit extra electron density within the groove between the two domains, suggesting the localization of a ligand binding site. One region of extra density is common to both structure, and can be interpreted as a glycerol molecule (present in the crystallization medium). The other one, only showing in the open conformation, was interpreted as acetates, a PEG fragment (also present in the crystallization medium) and a carbon dioxide molecule. Although these molecules present in the crystallization medium are probably not the functional CA3427 ligands, we used them to identify the putative binding site consisting of the residues less than 5 Å apart from the co-

<figure><img src="image_6.png" /></figure>

**Figure 1. CA3427 structures.** A) 2x7p is represented with alpha helices in red and beta-sheet in yellow. Secondary structures are numbered along the protein sequence. Carbon dioxide, glycerol and PEG molecules are in ball and sticks with green carbon and red oxygen atoms. 2 Water molecules are represented as blue spheres. B) Surface representation of the open conformation, 2x7p. The strictly conserved residues in the groove are colored in red and residues with conserved properties are marked in yellow. Carbon dioxide, glycerol and carboxyl molecules are in ball and sticks with green carbon and red oxygen atoms. Two water molecules are represented as blue spheres. C) The two conformations (2x7q: cyan, 2x7p: red) are superposed on domain I (bottom). The venus flytrap motion is illustrated by a black arrow (see Movie S1 for an animated view). doi:10.1371/journal.pone.0018528.g001

<figure><img src="image_7.png" /></figure>
crystallized molecules (Fig. 2A, Tables 2 and 3). Except for a small
cluster of strictly conserved polar residues (Glu 11, His 12, Glu 164
and Thr 167) located at the PEG/Glycerol interface, this putative
binding site is mainly hydrophobic (Fig. 2A). If we take into
account the water molecules (2048, 2049) filling the remaining
space of the groove between the two domains, the hypothetical
binding site can be extended to include 4 supplementary
conserved hydrophobic residues (G114, V119, L273, L279). This
putative binding site could accommodate extended hydrophobic
molecules such as long acyl chains (>C18) or carotenoids (Fig. 2B).
We noticed that the position of the conserved histidine (H11),

glutamate (E164) and threonine (T167) residues are not consistent
with the usual geometry of catalytic triads in hydrolases. The
precise function of the CA3427 protein thus remains to be
determined.

The differences between the two CA3427 crystal forms were
also analyzed in greater detail by comparing each domain
separately. Superimposition of the domain I (RMSD = 0.53 Å)
revealed only one major side chain reorganization within the
predicted binding site. When the PEG/acetates molecules are
present in the structure, the E11-O<sup>e2</sup> forms a hydrogen bond with
Y237-OH (distance 2.7 Å). Upon pointing outwards from the

<figure><img src="image_8.png" /><figcaption>**Figure 2. CA3427 putative binding site.** A) Stereo view of the binding site with non polar residues and tyrosine in yellow, polar residues in light blue, acidic residues in red and proline in green sticks. Carbon dioxide, glycerol and PEG molecules are in ball and sticks with green carbon and red oxygen atoms. Two water molecules are represented as blue spheres. (B) Surface representation of the CA3427 structure with a modeled C18 acyl chain (blue ball and stick) fitted in the Fo-Fc electron density map (green) computed on the open conformation structure (1.5σ). The conserved residues in the predicted binding site are colored in red.
doi:10.1371/journal.pone.0018528.g002</figcaption></figure>

<figure><img src="image_9.png" /></figure>
<table><caption>**Table 2.** Putative ligand-protein interactions.</caption><thead><tr><th>LI</th><th>LN</th><th>AL</th><th>RI</th><th>RN</th><th>D</th><th>LI</th><th>LN</th><th>AL</th><th>RI</th><th>RN</th><th>D</th></tr></thead><tbody><tr><td>309</td><td>CA</td><td>CA</td><td>-</td><td>97</td><td>ASP</td><td><strong>3.29</strong></td><td>315</td><td>ACY</td><td>O</td><td>-</td><td>44</td><td>ARG <strong>2.94</strong></td></tr><tr><td>310</td><td>CL</td><td>CL</td><td>-</td><td>71</td><td>GLY</td><td><strong>2.86</strong></td><td>315</td><td>ACY</td><td>OXT</td><td>-</td><td>37</td><td>LYS <strong>3.48</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>73</td><td>GLU</td><td><strong>3.84</strong></td><td></td><td></td><td></td><td>-</td><td>38</td><td>VAL <strong>3.83</strong></td></tr><tr><td>311</td><td>GOL</td><td>C3</td><td>-</td><td>154</td><td>GLY</td><td><strong>3.24</strong></td><td></td><td></td><td></td><td>-</td><td>40</td><td>GLU <strong>2.79</strong></td></tr><tr><td>311</td><td>GOL</td><td>O3</td><td>-</td><td>154</td><td>GLY</td><td><strong>2.83</strong></td><td></td><td></td><td></td><td>-</td><td>44</td><td>ARG <strong>2.83</strong></td></tr><tr><td>312</td><td>GOL</td><td>C1</td><td>-</td><td>13</td><td>PHE</td><td><strong>3.36</strong></td><td>316</td><td>ACY</td><td>C</td><td>-</td><td>60</td><td>LEU <strong>3.94</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>3.57</strong></td><td></td><td></td><td></td><td>-</td><td>62</td><td>GLU <strong>3.92</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>166</td><td>PHE</td><td><strong>3.89</strong></td><td></td><td></td><td></td><td>-</td><td>116</td><td>GLY <strong>3.48</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>245</td><td>TRP</td><td><strong>3.52</strong></td><td>316</td><td>ACY</td><td>CH3</td><td>-</td><td>63</td><td>ALA <strong>3.77</strong></td></tr><tr><td>312</td><td>GOL</td><td>C2</td><td>-</td><td>87</td><td>LEU</td><td><strong>3.63</strong></td><td>316</td><td>ACY</td><td>O</td><td>-</td><td>116</td><td>GLY <strong>3.16</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>3.37</strong></td><td>316</td><td>ACY</td><td>OXT</td><td>-</td><td>60</td><td>LEU <strong>3.65</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>166</td><td>PHE</td><td><strong>3.9</strong></td><td></td><td></td><td></td><td>-</td><td>62</td><td>GLU <strong>2.78</strong></td></tr><tr><td>312</td><td>GOL</td><td>C3</td><td>-</td><td>85</td><td>PRO</td><td><strong>3.34</strong></td><td></td><td></td><td></td><td>-</td><td>116</td><td>GLY <strong>3.29</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>86</td><td>LEU</td><td><strong>3.97</strong></td><td>317</td><td>ACY</td><td>C</td><td>-</td><td>252</td><td>ASN <strong>3.82</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU</td><td><strong>3.35</strong></td><td>317</td><td>ACY</td><td>O</td><td>-</td><td>252</td><td>ASN <strong>2.93</strong></td></tr><tr><td>312</td><td>GOL</td><td>O1</td><td>-</td><td>12</td><td>HIS</td><td><strong>3.64</strong></td><td></td><td></td><td></td><td>-</td><td>254</td><td>ARG <strong>3.6</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>13</td><td>PHE</td><td><strong>3.73</strong></td><td>317</td><td>ACY</td><td>OXT</td><td>-</td><td>259</td><td>PRO <strong>3.75</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>2.69</strong></td><td>318</td><td>ACY</td><td>C</td><td>-</td><td>52</td><td>ASN <strong>3.35</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>166</td><td>PHE</td><td><strong>3.18</strong></td><td>318</td><td>ACY</td><td>C</td><td>-</td><td>197</td><td>ASP <strong>2.98</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>245</td><td>TRP</td><td><strong>3.17</strong></td><td>318</td><td>ACY</td><td>CH3</td><td>-</td><td>197</td><td>ASP <strong>3.92</strong></td></tr><tr><td>312</td><td>GOL</td><td>O2</td><td>-</td><td>86</td><td>LEU</td><td><strong>3.31</strong></td><td>318</td><td>ACY</td><td>O</td><td>-</td><td>52</td><td>ASN <strong>2.66</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU</td><td><strong>2.84</strong></td><td></td><td></td><td></td><td>-</td><td>197</td><td>ASP <strong>2.57</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>2.71</strong></td><td>318</td><td>ACY</td><td>OXT</td><td>-</td><td>50</td><td>ASN <strong>3.43</strong></td></tr><tr><td>312</td><td>GOL</td><td>O3</td><td>-</td><td>85</td><td>PRO</td><td><strong>3.87</strong></td><td></td><td></td><td></td><td>-</td><td>52</td><td>ASN <strong>3.3</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>248</td><td>THR</td><td><strong>2.74</strong></td><td></td><td></td><td></td><td>-</td><td>197</td><td>ASP <strong>3.31</strong></td></tr><tr><td>313</td><td>AE3</td><td>C1</td><td>-</td><td>42</td><td>SER</td><td><strong>3.53</strong></td><td>2592</td><td>CO2</td><td>C</td><td>-</td><td>88</td><td>TRP <strong>3.61</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>60</td><td>LEU</td><td><strong>3.63</strong></td><td></td><td></td><td></td><td>-</td><td>163</td><td>TRP <strong>3.77</strong></td></tr><tr><td>313</td><td>AE3</td><td>C2</td><td>-</td><td>9</td><td>ILE</td><td><strong>3.71</strong></td><td>2592</td><td>CO2</td><td>O1</td><td>-</td><td>88</td><td>TRP <strong>3.99</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>41</td><td>GLY</td><td><strong>3.74</strong></td><td></td><td></td><td></td><td>-</td><td>163</td><td>TRP <strong>3.76</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>42</td><td>SER</td><td><strong>3.41</strong></td><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU <strong>3.31</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>60</td><td>LEU</td><td><strong>3.83</strong></td><td>2592</td><td>CO2</td><td>O2</td><td>-</td><td>88</td><td>TRP <strong>3.63</strong></td></tr><tr><td>313</td><td>AE3</td><td>C3</td><td>-</td><td>9</td><td>ILE</td><td><strong>3.49</strong></td><td></td><td></td><td></td><td>-</td><td>163</td><td>TRP <strong>3.88</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>190</td><td>TRP</td><td><strong>3.53</strong></td><td></td><td></td><td></td><td>-</td><td></td><td></td></tr><tr><td>313</td><td>AE3</td><td>C4</td><td>-</td><td>11</td><td>GLU</td><td><strong>3.7</strong></td><td>1307</td><td>GOL</td><td>C1</td><td>-</td><td>13</td><td>PHE <strong>3.85</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>190</td><td>TRP</td><td><strong>3.74</strong></td><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU <strong>3.97</strong></td></tr><tr><td>313</td><td>AE3</td><td>C5</td><td>-</td><td>190</td><td>TRP</td><td><strong>3.37</strong></td><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU <strong>3.63</strong></td></tr><tr><td>313</td><td>AE3</td><td>C6</td><td>-</td><td>12</td><td>HIS</td><td><strong>3.85</strong></td><td></td><td></td><td></td><td>-</td><td>245</td><td>TRP <strong>3.95</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>13</td><td>PHE</td><td><strong>3.88</strong></td><td>1307</td><td>GOL</td><td>C2</td><td>-</td><td>85</td><td>PRO <strong>3.49</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>3.72</strong></td><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU <strong>3.95</strong></td></tr><tr><td>313</td><td>AE3</td><td>O2</td><td>-</td><td>60</td><td>LEU</td><td><strong>3.98</strong></td><td></td><td></td><td></td><td>-</td><td>248</td><td>THR <strong>3.75</strong></td></tr><tr><td>313</td><td>AE3</td><td>O4</td><td>-</td><td>11</td><td>GLU</td><td><strong>3.83</strong></td><td>1307</td><td>GOL</td><td>C3</td><td>-</td><td>85</td><td>PRO <strong>3.73</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU</td><td><strong>3.08</strong></td><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU <strong>3.44</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>167</td><td>THR</td><td><strong>3.44</strong></td><td></td><td></td><td></td><td>-</td><td>248</td><td>THR <strong>3.64</strong></td></tr><tr><td>314</td><td>ACY</td><td>C</td><td>-</td><td>83</td><td>LYS</td><td><strong>3.92</strong></td><td>1307</td><td>GOL</td><td>O1</td><td>-</td><td>86</td><td>LEU <strong>3.95</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>264</td><td>THR</td><td><strong>3.92</strong></td><td></td><td></td><td></td><td>-</td><td>87</td><td>LEU <strong>2.87</strong></td></tr><tr><td>314</td><td>ACY</td><td>CH3</td><td>-</td><td>85</td><td>PRO</td><td><strong>3.79</strong></td><td></td><td></td><td></td><td>-</td><td>164</td><td>GLU <strong>3.11</strong></td></tr><tr><td>314</td><td>ACY</td><td>O</td><td>-</td><td>186</td><td>PRO</td><td><strong>3.67</strong></td><td>1307</td><td>GOL</td><td>O2</td><td>-</td><td>248</td><td>THR <strong>2.89</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>187</td><td>TRP</td><td><strong>3.65</strong></td><td></td><td></td><td></td><td>-</td><td>249</td><td>VAL <strong>3.86</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>188</td><td>SER</td><td><strong>3.79</strong></td><td>1307</td><td>GOL</td><td>O3</td><td>-</td><td>85</td><td>PRO <strong>3.25</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>264</td><td>THR</td><td><strong>3.23</strong></td><td></td><td></td><td></td><td>-</td><td>184</td><td>TYR <strong>3.67</strong></td></tr></tbody></table>

<table><caption>**Table 2.** Cont.</caption><thead><tr><th>LI</th><th>LN</th><th>AL</th><th>RI</th><th>RN</th><th>D</th><th>LI</th><th>LN</th><th>AL</th><th>RI</th><th>RN</th><th>D</th></tr></thead><tbody><tr><td>314</td><td>ACY</td><td>OXT</td><td>-</td><td>83</td><td>LYS <strong>3.44</strong></td><td></td><td></td><td></td><td></td><td>-</td><td>248</td><td>THR <strong>3.76</strong></td></tr><tr><td>315</td><td>ACY</td><td>C</td><td>-</td><td>37</td><td>LYS <strong>3.58</strong></td><td>1308</td><td>CA</td><td>CA</td><td>-</td><td>52</td><td>ASN <strong>2.9</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>40</td><td>GLU <strong>3.93</strong></td><td></td><td></td><td></td><td></td><td>-</td><td>197</td><td>ASP <strong>3.05</strong></td></tr><tr><td></td><td></td><td></td><td>-</td><td>44</td><td>ARG <strong>3.33</strong></td><td>1309</td><td>CA</td><td>CA</td><td>-</td><td>71</td><td>GLY <strong>2.91</strong></td></tr><tr><td>315</td><td>ACY</td><td>CH3</td><td>-</td><td>37</td><td>LYS <strong>3.29</strong></td><td></td><td></td><td></td><td></td><td>-</td><td>73</td><td>GLU <strong>3</strong></td></tr></tbody></table>

List of the shortest interatomic distances (<4 Å) between each of the ligands and the putative binding site amino acids. The distances were computed using NCONT (CCP4). Distance values are in bold for 2X7P, for 2X7Q otherwise. LI: ligand residue number, LN: ligand residue name, AL: ligand atom name, RI: protein residue number, RN: protein residue name, D: shortest interatomic distance between the ligand and the protein residue.
doi:10.1371/journal.pone.0018528.t002

binding pocket, it clears the space needed to accommodate the ligand. In the second structure, the above distance becomes 4.5 Å and the E11 side chain now points towards the inner part of the cavity. The domain II superimposition of the two structures (RMSD = 0.51 Å) again revealed a single noticeable difference within the binding site. The R112 side chain points towards the inside of the groove in the unliganded structure and outwards in the liganded one. Finally, in the closed conformation (i.e. without PEG/acetates molecules), D236-Oδ¹ (Domain I) and K170-N²⁵ (Domain II) are linked by a salt bridge (distance 3 Å) that is disrupted in the presence of ligand (distance 4.8 Å), thus opening the “Venus flytrap”.

## CA3427 exhibits a PBP fold

In search for hints about the biochemical function of the CA3427 protein, we compared the newly determined structures against those in the Protein Data Bank [33] using Dali [34,35] and VAST [36,37] through their online servers. The best matching structural homologs all correspond to Periplasmic Binding Proteins (PBP) with RMSD between 2.9 and 4 Å and very low sequence similarity (lower than 16% identical residues) with CA3427 (Table 4). All these matching

<table><caption>**Table 3.** Average B factor values for the interpreted ligands.</caption><thead><tr><th>PDB</th><th>LIGAND</th><th>B factor*</th><th>Occupancy</th></tr></thead><tbody><tr><td rowspan="5">2X7P</td><td>CA 309</td><td>57.1</td><td>1</td></tr><tr><td>CL 310</td><td>55.9</td><td>1</td></tr><tr><td>GLYCEROL 311</td><td>32.4</td><td>1</td></tr><tr><td>GLYCEROL 312</td><td>20.7</td><td>1</td></tr><tr><td>PEG 313</td><td>28.1</td><td>1</td></tr><tr><td rowspan="8">2X7Q</td><td>ACETYL 314</td><td>28.3</td><td>1</td></tr><tr><td>ACETYL 315</td><td>30.9</td><td>1</td></tr><tr><td>ACETYL 316</td><td>24.9</td><td>1</td></tr><tr><td>ACETYL 317</td><td>24.3</td><td>1</td></tr><tr><td>ACETYL 318</td><td>24.6</td><td>1</td></tr><tr><td>CO2 2592</td><td>22.8</td><td>1</td></tr><tr><td>GLYCEROL 307</td><td>14.5</td><td>1</td></tr><tr><td>CA 308</td><td>20.5</td><td>1</td></tr><tr><td>CA 309</td><td>24.2</td><td>1</td></tr></tbody></table>

*computed using CCP4.
doi:10.1371/journal.pone.0018528.t003

<figure><img src="image_10.png" /></figure>
<table><caption>**Table 4.** Closest structural homologs of 2X7P and 2X7Q (Dali server).</caption><thead><tr><th>PDB ID</th><th>Z-score</th><th>RMSD (Å)</th><th>LALI</th><th>NRES</th><th>ID (%)</th><th>PDB DESCRIPTION</th></tr></thead><tbody><tr><td>2X26</td><td>25.4/24.7</td><td>3.3/3.6</td><td>283/280</td><td>308</td><td>12/12</td><td>Periplasmic aliphatic sulphonates-binding</td></tr><tr><td>3IX1</td><td>24.8/26</td><td>3.1/2.8</td><td>276/279</td><td>301</td><td>16/15</td><td>N-formyl-4-amino-5-aminomethyl-2-methylpyrimidine</td></tr><tr><td>3E4R</td><td>23.3/25.5</td><td>3/3</td><td>274/277</td><td>291</td><td>14/13</td><td>Nitrate transport protein</td></tr><tr><td>2I48</td><td>21.9/23.1</td><td>3.3/3.4</td><td>286/286</td><td>399</td><td>10/10</td><td>Bicarbonate transporter</td></tr></tbody></table>

Left and right values correspond to 2X7P and 2X7Q, respectively. RMSD (root mean square deviation): average distance between the Cα backbone atoms of the superimposed proteins. LALI: total number of aligned residues. NRES: length of the homologous protein sequence. ID: percentage of identical residues within the optimal alignment.

doi:10.1371/journal.pone.0018528.t004

proteins are of bacterial origin and members of the class II PBP-like fold family: two similar intertwined domains of 3 layers (α/β/α) each. The β part is a duplication of mixed beta-sheet of 5 strands ordered as 2-1-3-5-4 with strand 5 antiparallel to the others. This PBP structural module, associated with a large variety of functions [38], is found in prokaryotic and eukaryotic protein families, as well as in the soluble part of the eukaryotic ionotropic glutamate

receptors [39]. Due to their functional versatility, PBP have been considered promising protein engineering targets for biotechnology and drug delivery applications [38,39,40].

The multiple sequence alignment of fungal proteins homologous to CA3427 highlights the conservation of the residues delimiting the groove (Fig. 3) further supporting its functional relevance. Moreover, the known binding sites of other PBP proteins are

<figure><img src="image_11.png" /><figcaption>**Figure 3.** Multiple alignment of 29 selected sequences used for phylogeny. All fungi sequence names are in black except *C. albicans* (in red) *D. hanenii* and *Y. lipolitica* (in green). Other eukaryote sequence names are in magenta, bacterial sequence names are in light blue for Flavobacteria/Bacteroidetes and dark blue for the last 4 bacteria (outgroup). All the sequences are aligned together in one shot but similarity at each position is computed separately for the outgroup and the other sequences except for position 41 and 163 where strict conservation (white letters with red background) is observed. Red or orange letters in the alignment represent similarity (computed using a Risler matrix with a 0.9 threshold in ESPrift) within each group. The secondary structure of CA3427 computed with DSSP, is also represented above the alignment. Purple triangles represent the binding site defined in the CA3427 structure. This alignment shows clearly the conservation of the binding site among different eukaryotes and Flavobacteria/Bacteroidetes but not in the outgroup.
doi:10.1371/journal.pone.0018528.g003</figcaption></figure>

<figure><img src="image_12.png" /></figure>
always located at the interface between the two domains as seen in the CA3427 structure. The corresponding conserved residues (Fig. 2A) thus define a new subfamily of PBP-related proteins, likely sharing a common ligand and a related function.

It hence appears that CA3427 was wrongly annotated (UNIPROT: Q59X88) as involved in the synthesis of thiamine, on the basis of its weak homology with the *Saccharomyces cerevisiae*

ThI13 protein (UNIPROT: Q07748) (<20% identity over the full sequences length, blastp E value 0.29), since none of the binding site residues are conserved between the two sequences. In fact, the THI13 orthologue in *Candida albicans* corresponds to the Q5A3Y5 protein with which it shares 75% identity over 338 residues (E value: $10^{-156}$). The function of the new PBP-like family defined by CA3427 remains to be determined.

<figure><img src="image_13.png" /><figcaption>**Figure 4. Phylogeny.** A) **Phylogeny of the CA3427 homologs in Fungi.** The reference list of fungal species is from [23]. Multiple strains of the same species have been removed for clarity. The CA3427 sequence (in red in all trees) is from C. albicans strain SC5314. This unrooted tree was computed on the Phylogeny.fr web server [30], using the default option of the "advanced mode" w/o Gblocks. The final alignment includes 209 ungapped positions. Branch support estimates are indicated in red, and branches have been collapsed for values <50%. CA3427 homologs were found in all species except for saccharomycetales (in green), and cluster according to [23]. No species from the WGD clade (e.g. *Saccharomyces cerevisiae*) appear to possess a CA3427-like protein. B) **Evolutionary relationships between the bacterial PBPs and their eukaryotic CA3427-like homolog.** The phylogenetic analysis includes representative sequences from Cytophaga-Flavobacteria (light blue), Fungi (black and green) (as in Fig. 2), other eukaryotes (in magenta), and more remote bacterial sequences defining and outgroup (dark blue). This tree was computed on Phylogeny.fr web server [30], using the default option of the "advanced mode" without Gblocks. Branch support estimates are indicated in red, and branches have been collapsed for values <50%. The topology of this tree is consistent with the hypothesis that the original CA3427-like PBP was transferred into the eukaryotic gene pool from a cytophagia/flavobacteria into an ancestral opisthokont.
doi:10.1371/journal.pone.0018528.g004</figcaption></figure>

<figure><img src="image_14.png" /></figure>
Interestingly, this new protein family is also well represented in flavobacteria and other members of the Cytophaga-Flexibacter-Bacteroides phylum of eubacteria to the exclusion of other prokaryotes. None of these proteins were functionally characterized. The binding site of CA3427 is strongly hydrophobic and can accommodate acyl chains much longer than C18 (Fig. 2B). It could thus participate in the detection, transport and/or processing of high molecular weight lipids (or carotenoids) in flavobacteria and fungi. Members of the CA3427 family are also present in two ancestral eukaryotes, the choanoflagellate *Monosiga brevicollis* and the placozoan *Trichoplax adhaerens*. These sequences share more than 30% identical residues with the CA3427 protein.

In order to investigate on the evolutionary origin of the C3427 protein family, we performed phylogenetic reconstructions. Figure 4A shows that the CA3427 phylogeny precisely follows the reference fungi classification [23]. A single inconsistency is the absence of a CA3427 homolog in *Saccharomyces cerevisiae*. In fact, all known species from the WGD (Whole Genome Duplication) saccharomycetales clade lack a CA3427 homolog, strongly suggesting that the loss of this gene coincided with the separation of the WGD clade from the CTG clade (the species that translate CTG as Serine instead of leucine) of saccharomycetales (e.g. Debaryomyces). In a more comprehensive phylogenetic reconstruction we included all members of the CA3427 family and, as an outgroup, representatives of bacterial PBP sequences not exhibiting the CA3427 binding site. The resulting tree (Figure 4B) strongly suggests that the eukaryotic and flavobacteria members of the CA3427 family originated from a common ancestral gene.

The presence of CA3427 homologs in the Bacteroides phylum of eubacteria strongly suggest that the PBP-like CA3427 protein has a very ancient bacterial origin. The divergence between the mainstream PBPs and the CA3427-like subfamily probably occurred early on the branch leading to the Bacteroidetes, after its separation from the branch leading to the other main groups of eubacteria (i.e. Proteobacteria, Planctomyces, Firmicutes).

Finally, the surprising presence of a CA3427 homolog in the genome of the choanoflagellate *Monosiga brevicollis* as well as in the genome of placozoan *Trichoplax adhaerens*, the most basal invertebrate form, supports a scenario of horizontal transfer by which all eukaryotic CA3427 homologs originated from the above Bacteroidetes ancestor. Interestingly, Bacteroidetes species such as Algoriphagus are commonly found in association with modern choanoflagellates [41], thus providing opportunities for gene exchanges. The shared presence of CA3427 homologs in most fungi, the only sequenced choanoflagellate (closest unicellular relative of animals) and the only known placozoan (the closest multicellular relative of animals) strongly suggests that the transfer of the bacterial gene to an ancestral eukaryote occurred at the very basis of the opisthokont lineage, before the radiation of fungi [42] (Fig. 5). This gene was later lost in the branch leading to modern animals. Unexpectedly, the CA3427-like family of PBP thus provides a new marker to probe the early scheme of eukaryotic evolution [42].

We solved the 3-D structure of the Candida albicans CA3427 gene product, selected on the basis of its conservation among pathogenic fungi, and thus a potential target for new antifungal drugs. The structure of the protein unambiguously revealed a PBP fold, despite a low level of sequence similarity with previously known members of this family. In addition, the 3D structures allowed the precise delineation of a binding site, defined by highly conserved residues in the vicinity of co-crystallized ligands. The conformational change of the CA3427 protein upon ligand binding illustrates the venus fly trap motion already documented in other PBP structures [38, 39, 40, 43]. A phylogenetic analysis of

<figure><img src="image_15.png" /><figcaption>**Figure 5. Parsimonious evolutionary scenario for the presence of CA3427-homologs in modern eukaryotes.** A horizontal transfer (HGT) is proposed to have occurred from a bacteroidetes to an ancestral unicellular heterotrophic eukaryote prior to the divergence of the main branches leading to fungi and animalia. The branching of the various phyla is adapted from [42].</figcaption></figure>

doi:10.1371/journal.pone.0018528.g005

the CA3427 protein family indicates that it originated in Bacteroidetes before being transferred to an ancestral eukaryote prior to the divergence between the fungi and animal lineages. The intriguing (albeit remote) possibility that the acquisition of this gene might be linked to the evolution towards multicellularity is a strong incentive for further functional studies. Furthermore, the conservation of this family of proteins in all pathogenic fungi coupled to its absence in animals makes it a good target for the design of new drugs against candidiasis and other diseases caused by fungi.

## Supporting Information

**Movie S1 Animated gif for the Morphing of the CA3427 venus flytrap motion.** The two structures were submitted to the Morph Server [44], the pictures were generated by pymol and concatenated to an animated gif with the ImageMagick convert function. The movie illustrates the venus flytrap motion of 2x7q leading to the 2x7p conformation. The strictly conserved residues in the groove are colored in red and residues with conserved properties are marked in yellow. Ligands are modeled as they appear in the opened conformation (2x7p). Carbon dioxide, glycerol and carboxyl molecules are in ball and sticks representation with green carbon and red oxygen atoms. Two water molecules are represented as blue spheres.

(GIF)

## Acknowledgments

We acknowledge the use of the PACA-Bioinfo (IBISA) bioinformatic platform, the IMM-IFR88 proteomic platform and the support of the experimental team of the IGS laboratory. We thank Dr. Pascal Hingamp for reading the manuscript.

## Author Contributions

Conceived and designed the experiments: CA JMC VM NM TR. Performed the experiments: SS CM CA. Analyzed the data: SS JMC CA. Contributed reagents/materials/analysis tools: VM NM TR SS JMC CA. Wrote the paper: SS JMC CA.

<figure><img src="image_16.png" /></figure>
References

1. Wilson LS, Reyes CM, Stolpman M, Speckman J, Allen K, et al. (2002) The direct cost and incidence of systemic fungal infections. Value Health 5: 26-34.

2. Gudlaugsson O, Gillespie S, Lee K, Vande Berg J, Hu J, et al. (2003) Attributable mortality of nosocomial candidemia, revisited. Clin Infect Dis 37: 1172-1177.

3. Bodey G, Bueltmann B, Duguid W, Gibbs D, Hanak H, et al. (1992) Fungal infections in cancer patients: an international autopsy survey. Eur J Clin Microbiol Infect Dis 11: 99-109.

4. Groll AH, Shah PM, Mentzel C, Schneider M, Just-Nuebling G, et al. (1996) Trends in the postmortem epidemiology of invasive fungal infections at a university hospital. J Infect 33: 23-32.

5. Pfaller MA, Castanheira M, Diekema DJ, Messer SA, Moet GJ, et al. (2010) Comparison of European Committee on Antimicrobial Susceptibility Testing (EUCAST) and Etest methods with the CLSI broth microdilution method for echinocandin susceptibility testing of Candida species. J Clin Microbiol 48: 1592-1599.

6. Pfaller MA, Diekema DJ (2007) Epidemiology of invasive candidiasis: a persistent public health problem. Clin Microbiol Rev 20: 133-163.

7. Vanden Bossche H, Dromer F, Improvisi I, Lozano-Chiu M, Rex JH, et al. (1998) Antifungal drug resistance in pathogenic fungi. Med Mycol 36 Suppl 1: 119-128.

8. The PROFUN project. Available: http://www.igs.cnrs-mrs.fr/Bigs/. Accessed 2011 March 22.

9. Hendrickson WA (1991) Determination of macromolecular structures from anomalous diffraction of synchrotron radiation. Science 254: 51-58.

10. Santini S, Monchois V, Mouz N, Sigoillot C, Rousselle T, et al. (2008) Structural characterization of CA1462, the Candida albicans thiamine pyrophosphokinase. BMC Struct Biol 8: 33.

11. Abergel C, Coutard B, Byrne D, Chenivesse S, Claude JB, et al. (2003) Structural genomics of highly conserved microbial genes of unknown function in search of new antibacterial targets. J Struct Funct Genomics 4: 141-157.

12. Waldo GS, Standish BM, Berendzen J, Terwilliger TC (1999) Rapid protein-folding assay using green fluorescent protein. Nat Biotechnol 17: 691-695.

13. Doubie S (1997) Preparation of selenomethionyl proteins for phase determination. Methods Enzymol 276: 523-530.

14. Leslie AGW (1992) Recent changes to the MOSFLM package for processing film and image plate data. Protein Crystallography 26.

15. Kabsch W (1988) Evaluation of Single-Crystal X-ray Diffraction Data from a Position-Sensitive Detector. J Appl Cryst 21: 916-924.

16. Collaborative Computational Project N (1994) The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr 50: 760-763.

17. Terwilliger TC, Berendzen J (1999) Automated MAD and MIR structure solution. Acta Crystallogr D Biol Crystallogr 55: 849-861.

18. Vonrhein C, Blanc E, Roversi P, Bricogne G (2007) Automated structure solution with autoSHARP. Methods Mol Biol 364: 215-230.

19. Emsley P, Cowtan K (2004) Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 60: 2126-2132.

20. Adams PD, Grosse-Kunstleve RW, Hung LW, Ioerger TR, McCoy AJ, et al. (2002) PHENIX: building new software for automated crystallographic structure determination. Acta Crystallogr D Biol Crystallogr 58: 1948-1954.

21. Claude JB, Suhre K, Notredame C, Claverie JM, Abergel C (2004) CaspR: a web server for automated molecular replacement using homology modelling. Nucleic Acids Res 32: W606-609.

22. Laskowski RA, MacArthur MW, Moss DS, M TJ (1993) PROCHECK - a program to check the stereochemical quality of protein structures. J App Cryst 26: 283-291.

23. Wang H, Xu Z, Gao L, Hao B (2009) A fungal phylogeny based on 82 complete genomes using the composition vector method. BMC Evol Biol 9: 195.

24. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403-410.

25. BROAD-FGI. Broad Institute Fungal Genome Initiative. Available: http://www.broad.mit.edu/annotation/fgi/. Accessed 2011 March 22.

26. JGI. Department of Energy Joint Genome Institute. Available: http://genome.jgi-psf.org/. Accessed 2011 March 22.

27. NCBI. National Center for Biotechnology Information. Available: ftp://ftp.ncbi.nih.gov/genomes/. Accessed 2011 March 22.

28. RFCG. Resources for Fungal Comparative Genomics. Available: http://fungal.genome.duke.edu/.

29. FGR. Fungal Genome Research. Available: http://fungalgenomes.org/. Accessed 2011 March 22.

30. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, et al. (2008) Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res 36: W465-469.

31. Suhre K, Sanejouand YH (2004) ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement. Nucleic Acids Res 32: W610-614.

32. Mao B, Pear MR, McCammon JA, Quiocho FA (1982) Hinge-bending in L-arabinose-binding protein. The "Venus's-flytrap" model. J Biol Chem 257: 1131-1133.

33. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, et al. (2000) The Protein Data Bank. Nucleic Acids Res 28: 235-242.

34. Holm L (1998) Unification of protein families. Curr Opin Struct Biol 8: 372-379.

35. Holm L, Kaariainen S, Rosenstrom P, Schenkel A (2008) Searching protein structure databases with DaliLite v.3. Bioinformatics 24: 2780-2781.

36. Gibrat JF, Madej T, Bryant SH (1996) Surprising similarities in structure comparison. Curr Opin Struct Biol 6: 377-385.

37. Madej T, Gibrat JF, Bryant SH (1995) Threading a database of protein cores. Proteins 23: 356-369.

38. Felder CB, Graul RC, Lee AY, Merkle HP, Sadee W (1999) The Venus flytrap of periplasmic binding proteins: an ancient protein module present in multiple drug receptors. AAPS PharmSci 1: E2.

39. Chen PE, Wyllie DJ (2006) Pharmacological insights obtained from structure-function studies of ionotropic glutamate receptors. Br J Pharmacol 147: 839-853.

40. Oh BH, Kang CH, De Bondt H, Kim SH, Nikaido K, et al. (1994) The bacterial periplasmic histidine-binding protein. structure/function analysis of the ligand-binding site and comparison with related proteins. J Biol Chem 269: 4135-4143.

41. Alegado RA, Ferriera S, Nusbaum C, Young SK, Zeng Q, et al. Complete Genome Sequence of Algoriphagus sp. PR1, Bacterial Prey of a Colony-forming Chaoanoflagellate. J Bacteriol.

42. Cavalier-Smith T (2004) Only six kingdoms of life. Proc Biol Sci 271: 1251-1262.

43. Quiocho FA, Ledvina PS (1996) Atomic structure and specificity of bacterial periplasmic receptors for active transport and chemotaxis: variation of common themes. Mol Microbiol 20: 17-25.

44. Krebs WG, Gerstein M (2000) The morph server: a standardized system for analyzing and visualizing macromolecular motions in a database framework. Nucleic Acids Res 28: 1665-1675.

<figure><img src="image_17.png" /></figure>