Identification of Transport-critical Residues in a Folate Transporter from the Folate-Biopterin Transporter (FBT) Family*

The Synechocystis Slr0642 protein and its plastidial Arabidopsis (Arabidopsis thaliana) ortholog At2g32040 belong to the folate-biopterin transporter (FBT) family within the major facilitator superfamily. Both proteins transport folates when expressed in Escherichia coli. Because the structural requirements for transport activity are not known for any FBT protein, we applied mutational analysis to identify residues that are critical to transport and interpreted the results using a comparative structural model based on E. coli lactose permease. Folate transport was assessed via the growth of an E. coli pabA abgT strain, which cannot synthesize or take up folates or p-aminobenzoylglutamate. In total, 47 residues were replaced with Cys or Ala. Mutations at 22 positions abolished folate uptake without affecting Slr0642 expression in membranes, whereas other mutations had no effect. Residues important for function mostly line the predicted central cavity and are concentrated in the core α-helices H1, H4, H7, and H10. The essential residue locations are consistent with a folate-binding site lying roughly equidistant from both faces of the transporter. Arabidopsis has eight FBT proteins besides At2g32040, often lacking conserved critical residues. When six of these proteins were expressed in E. coli or in Leishmania folate or pterin transporter mutants, none showed evidence of folate or pterin transport activity, and only At2g32040 was isolated by functional screening of Arabidopsis cDNA libraries in E. coli. Such negative data could reflect roles in transport of other substrates. These studies provide the first insights into the native structure and catalytic mechanism of FBT family carriers.

Carrier-mediated folate transport into cells and between subcellular compartments is vital in organisms that require folates and is also vital in plants due to the way the folate synthesis pathway and folates are compartmented. Folates are made in mitochondria and exported to the cytosol, plastids, and vacuole (5)(6)(7). This implies the existence of several organellar folate transporters. Because plant cells can import folates, there must also be a transporter in the plasmalemma (8,9).
Taking all organisms together, six types of folate transporter have been cloned so far. Five occur in mammals (10 -12); of these, one (from the mitochondrial carrier family) occurs also in plants (13). The sixth type, the folate-biopterin transporter (FBT) family, also called the BT1 family, occurs in protists, cyanobacteria, and plants (3,14,15). The FBT family belongs to the major facilitator superfamily (MFS) (16). MFS proteins consist of a single polypeptide chain containing 12 transmembrane ␣-helices, with both N and C termini on the cytoplasmic side of the membrane (17). MFS transporters can have uniport, symport, or antiport mechanisms (18). Structures have now been determined for MFS proteins, including the lactose permease LacY, the glycerol-3-phosphate antiporter GlpT, and the multidrug transporter EmrD (19 -21). These structures imply a general architecture for the superfamily, with substantial conservation of secondary and tertiary structure elements, and a common substrate translocation mechanism (17,21).
FBT family members shown experimentally to transport folate include the FT1 and FT5 proteins from the protistan parasite Leishmania (22,23), Slr0642 from the cyanobacterium Synechocystis, and its Arabidopsis ortholog At2g32040, which is located in the plastid envelope (15). Nothing is known for any of them about their structure in the membrane, the basis for sub-strate preference, or the catalytic mechanism, despite the centrality of FBTs in Leishmania to antifolate drug resistance (24).
The absence of information on the structure and mechanism of FBT carriers is due partly to the paucity of satisfactory eukaryotic or prokaryotic expression systems. Although Slr0642 and At2g32040 have both been functionally expressed in Escherichia coli, the experiments that could be done were limited by occasional spontaneous acquisition of folate uptake by the host cells (15). For example, the use of site-directed mutagenesis to probe the transport mechanism was precluded.
Besides At2g32040, the Arabidopsis genome encodes eight other FBT proteins, all phylogenetically more divergent from the cyanobacterial and Leishmania proteins than At2g32040 (15). Their function is unknown. They might a priori be supposed to transport folates, because folates (or pterins) are the only known substrates of the FBT family, and mitochondrial, vacuolar, and plasmalemmal folate carriers remain to be identified.
In this study, we first improved an E. coli expression system for FBT proteins. Using this system along with site-directed mutagenesis led to the definition of transport-critical residues in Slr0642. The expression system also enabled testing of more Arabidopsis FBT proteins for folate transport and functional screening of cDNA libraries. Last, we explored whether several FBT proteins could rescue the folate or biopterin growth requirements of Leishmania folate or biopterin transporter mutants.
Site-directed Mutagenesis of Slr0642-Site-directed mutagenesis was performed using the QuikChange II XL site-directed mutagenesis kit (Stratagene, La Jolla, CA). The pLOI707HE vector containing the slr0642 cDNA plus a C-terminal hexahistidine tag sequence was used as template. Synthetic oligonucleotides used for mutagenesis are listed in supplemental Table  S1. Mutations were confirmed by sequencing.
Arabidopsis FBT Expression in E. coli-Full-length cDNAs encoding At5g25040, At5g25050, At1g79710, At5g10820, At5g54860, and At1g04570 were isolated by PCR using cDNA prepared from Arabidopsis rosette leaves as a template and the primers listed in supplemental Table S1. The SuperScript TM III first-strand synthesis system (Invitrogen) was used for reverse transcription of mRNA. For expression in E. coli, full-length cDNAs were cloned into pLOI707HE as described (15). Chimeric constructs (with FBT sequences truncated just before the first predicted transmembrane ␣-helix and fused to the N-terminal 27 residues of Slr0642) were also generated as described (15), using primers listed in supplemental Table S1. Constructs were made in E. coli strain DH5␣, sequence-verified, and then introduced into the E. coli double mutant ⌬pabA abgT::kan, harboring pACYC-RIL (Stratagene, La Jolla, CA), which contains three rare E. coli tRNA genes and the cat marker. Transformants were grown on LB plates containing 10 g/ml tetracycline, 50 g/ml kanamycin, and 20 g/ml chloramphenicol.
Arabidopsis FBT Expression in Leishmania donovani Mutants-Two versions of At2g32040 and the protein-coding regions of At5g25040, At5g25050, At1g79710, At5g10820, and At5g54860 were cloned into the unique BglII site of pIR1SAT, inserting a Kozak sequence (CCACC) just before the start codon. pIR1SAT is an integrating Leishmania expression vector containing streptothricin acetyltransferase as a selectable marker. Transfection into Leishmania of the linear targeting fragment of the vector and selection for streptothricin resistance yields parasites where the construct has replaced a copy of the ribosomal small subunit gene and is transcribed at high levels from the polI promoter. 6 All sequences were amplified by PCR using primers to add 5Ј and 3Ј BglII, BclI, or BamHI sites and were digested with the respective enzyme before ligating into pIR1SAT. For At1g79710, which has a BglII site at bp 81, the first 69 bp were replaced by a start codon using a primer that ablated the BglII site; this deleted about half the region before the first transmembrane domain. One At2g32040 construct was truncated to remove the plastid targeting sequence by replacing the first 309 bp with a start codon. For the other, the first 357 bp were replaced by the first 26 codons of Leishmania major BT1, in which 12 G or C residues were changed to A or T (to eliminate a GC-rich region), and a BspEI site was introduced to join the fragments. All constructs were sequence-verified. The plasmids were linearized with SwaI and transfected into bt1 Ϫ or MTXA5 mutant L. donovani promastigotes, and transformants were isolated by plating on semisolid M199 medium containing 100 g/ml nourseothricin, using published methods (28,29). A positive control episomal add-back line was generated by the same procedure for the MTXA5 null mutant by introducing the L. major FT1 gene in the pXGNEO vector (B2437) and plating on medium containing 10 g/ml G418 (3,30). Episomal addback lines for the bt1 Ϫ null mutants were as described (22). Growth of wild-type L. donovani and bt1 Ϫ and MTXA5 mutants and transformants under semidefined conditions used folate-free RPMI medium (Invitrogen) supplemented with 30 mM HEPES, pH 7.4, 63 M adenosine, 5 g/ml hemin, and 1% (v/v) heat-inactivated fetal calf serum. The growth of each transport-deficient parental line was compared with that of its episomal add-back and the Arabidopsis FBT transformants.  Screening of Arabidopsis cDNA Libraries-Arabidopsis cDNA expression libraries from 3-day-old hypocotyls (31) were obtained from the Arabidopsis Biological Resource Center (stock numbers CD4-14, CD4-15, and CD4-16). Mass excision was made using the ExAssist Interference-Resistance Helper Phage kit (Stratagene, La Jolla, CA) according to the manufacturer's instructions. Briefly, a portion of the bacteriophage libraries was combined with E. coli XL1-blue MRFЈ cells and ExAssist helper phage at a 1:10:100 ratio, incubated 15 min at 37°C, mixed with LB for a 3-h incubation at 37°C with shaking, heated at 70°C for 20 min, and centrifuged for 10 min at 1,000 ϫ g. A portion of the resulting supernatant (excised phagemids) was incubated with E. coli SOLR cells for 20 min at 37°C, plated onto 10 15-cm LB-ampicillin (100 g/ml) agar plates, and incubated overnight at 37°C. The resulting clones (ϳ10 6 ) were harvested in liquid LB and pooled for large scale plasmid DNA isolation by alkaline lysis. E. coli double mutant (⌬pabA abgT::kan) cells were transformed with 5 g of plasmid DNA and plated on selective medium (M9 minimal medium containing 0.5 mM IPTG, 100 g/ml ampicillin, 50 g/ml kanamycin, 20 g/ml chloramphenicol, 2 g/ml sulfathiazole, and 11 M 5-FTHF) for screening.
Bioinformatics-Because Slr0642 belongs to the MFS superfamily of transporters (16), a homology model of the transporter was built using the known structure of the lactose permease LacY (19) as template. However, the sequence homology between LacY and Slr0642 is only 9.5%, precluding an accurate alignment based on sequence identity alone. Therefore, additional criteria were applied to obtain a more reliable alignment, because the biophysical and chemical properties of related membrane proteins are often better conserved than the sequence identity. First, the secondary structure elements were determined using PSIPRED (32) to obtain a coarse alignment of the ␣-helices (Fig. 1A). Second, hydropathy profiles of LacY and Slr0642 were used to improve the alignment, because they are sensitive measures of the underlying structural similarity even when sequence similarity is virtually absent (33) (Fig. 1B). The orientation of the ␣-helices of Slr0642 was then improved in an iterative way by applying two further principles: (a) that residues that interact with the alkyl chains of the lipid bilayer are hydrophobic, whereas those in protein interactions or in the water phase are hydrophilic; and (b) that residues that interact with the membrane are less conserved than those involved in substrate and protein interactions, because there is less selective pressure on hydrophobic interactions (34). A hydropathy score was assigned to each residue by using the hydrophobicity scale of Kyte and Doolittle (35). A conservation score was assigned by ConSurf (36) by using an alignment of Slr0642 orthologs. For the iterative process, first a structure was built based on the alignment by using Swiss-Model (37). Second, the conservation and hydropathy scores were assigned to residues by using a color scale. Third, the structure was inspected visu-ally with PyMol (38), and the alignment was altered to correct the inconsistencies. In total, 18 rounds of refinement were applied to generate a model that was most consistent with these criteria (see supplemental Fig. S1 and its interactive version, supplemental Fig. S2, A and B). The final pairwise alignment was used to calculate a comparative model of Slr0642 with MODELLER (39) by using the LacY structure (19) as template. The highest ranking model was analyzed with MOLPROBITY (40) to check the , , and C␤ angles of the residues.

RESULTS
Construction and Validation of an Improved E. coli pabA abgT Strain Suitable for Folate Uptake Studies-Because wildtype E. coli cells cannot take up folates, the functional expression of plasmid-encoded folate transporters can be detected in a pabA host strain (unable to synthesize the folate precursor pABA) by monitoring growth on medium containing a folate, 5-FTHF (15). However, this assay suffers interference from the occasional appearance of cells that spontaneously acquire the ability to grow on 5-FHTF (15), presumably due to sporadic manifestation of a latent capacity for uptake of folate or folate breakdown products that give rise to pABA.
The E. coli abgT gene encodes a cryptic transporter of pABAglutamate (a folate breakdown product) that is activated by point mutations in its promoter region (25,41,42), making it a prime suspect to explain the interference problem. We therefore constructed a pabA abgT strain and compared its propensity to grow in liquid minimal medium plus 5-FHTF with that of the parent pabA strain (Fig. 2A). A strain harboring the folate transporter Slr0642 was included for comparison. In triplicate observations, the double mutant cultures showed only slight and uniform growth, most likely sustained by traces of pABA produced by 5-FHTF breakdown in the medium.
In contrast, all of the single mutant cultures began to grow exponentially after various lag times. This result was verified by streaking the double mutant on minimal medium plates containing 5-FHTF; in multiple trials, no colonies appeared. We therefore concluded that AbgT expression is responsible for the interference problem and used the pabA abgT mutant for all further work.

Functionality of a His-tagged Slr0642 Protein and Vector
Choice-Synechocystis Slr0642 was chosen for site-directed mutagenesis because it is more readily expressed in E. coli than the At2g32040 protein (15). A C-terminal His tag was added to allow immunodetection, which did not alter performance; the tagged and The E. coli pabA abgT strain was transformed with pLOI707HE with no insert (V) or containing native slr0642 (WT) or 19 mutants thereof. Two independent clones of each mutant construct were streaked on minimal medium plus 0.5 mM IPTG alone (control) or containing 3.6 M pABA or 11 M 5-FTHF. Clones were streaked once on pABA-supplemented plates and three times in succession on the others. The pABA plates were photographed 1 day after streaking, and others were photographed after 2 days. Red, transport-critical residues; blue, non-critical residues. B, Western blot analysis of Slr0642 protein levels in membrane preparations from the above clones. Membrane proteins (25 g/track) were separated by SDS-PAGE. The position of the 52-kDa molecular size marker is shown. Strictly conserved residues are indicated by a plus sign, conservatively replaced residues by "c", and non-conserved residues by a minus sign. untagged versions of Slr0642 gave identical growth of the pabA abgT mutant on 5-FHTF medium (Fig. 2B). In this and subsequent work, we retained the pLOI707HE vector used previously (15). This is a 13-kb, medium copy plasmid with a tac promoter and a copy of the lacI q repressor gene to ensure tight control of expression by IPTG (43). The high copy, less tightly IPTG-regulated pBluescript plasmid gave poorer results; colonies grew only half as fast as when pLOI707HE was used. Others have noted that high copy plasmids can be unsuitable for membrane protein expression in E. coli (44). Sensitivity of the Growth Assay-To enhance the sensitivity of the growth assay to reductions of transport activity in mutant Slr0642 proteins, it was carried out at a 5-FTHF concentration determined to be subsaturating (11 M). Thus, after a standard 2-day incubation, growth was markedly less at 11 M than at 50 M 5-FTHF (Fig. 2C). When 5-FTHF concentration was lowered in 1-M steps from 11 M, growth ceased to be evident at a concentration of 4 M (Fig. 2C). If the 5-FTHF concentration is limiting for growth, it follows that transporter activity is likewise limiting and that if a concentration drop from 11 to 4 M (ϳ60%) suffices to halt growth, then a similar activity drop would also do so. That a mutant protein retaining as much as 40% of wild type activity would fail to support growth makes the assay a responsive one.

Residue
Mutagenesis Round 1-Seventeen potential active site residues were selected (Table 1) based on their predicted location in transmembrane ␣-helices, their chemistry, and their pattern of conservation among FBT family members shown to transport folates. Polar residues were emphasized because these could in principle face the hydrophilic cavity of the carrier and interact with the folate substrate. Twelve of the residues are conserved or conservatively replaced in all of the folate-transporting FBT proteins and five others are conserved between Slr0642 and At2g32040. Two non-conserved Arg residues (Arg-158 and Arg-301) were included for comparison. All residues were replaced by Cys except for Cys-371, where Ala was used. Cys substitution mutants have been extensively used to investigate the structure and function of transporters (e.g. see Ref. 45).
Of the 19 mutations tested, 12 were unable to support growth of the pabA abgT mutant on 5-FTHF medium, whereas the others had no noticeable effect (Fig. 3A). Western blotting of membrane fractions from the cells harboring the mutant constructs showed that Slr0642 protein was always expressed, albeit at varying levels (Fig. 3B). However, it varied about as much among the constructs that were functional as among those that were not, making it unlikely that the inactive mutant phenotype arose from insufficient protein. Although the protein expression results do not exclude the possibility that the deleterious mutations resulted in incorrect insertion into the membrane, not inhibition of transport per se, aberrant insertion is the less economical explanation and is inconsistent with the high incidence of criticality.
Of the 12 target amino acids conserved in known FBT folate transporters, nine were critical; of the five less conserved residues only three were critical; and of the two non-conserved residues, neither was critical (Table 1). There is thus a general correlation between conservation and importance for transport function. (Table 2) was selected based on (a) a high degree of conservation among cyanobacterial and plant orthologs of Slr0642, (b) a predicted location in the aqueous cavity, and (c) a polar or charged nature (or commonly found to be interacting with folate in other proteins of known structure) (supplemental Fig.  S3). All were changed to Cys. Of the 16 residues, seven were important for transport (supplemental Fig. S4), six of which were predicted to face the central cavity.

Mutagenesis Rounds 2 and 3-A second set of 16 residues
Collectively, the results of rounds 1 and 2 provided good support for the overall structure of the comparative model because a high proportion of the residues that are important for transport line the predicted central cavity. In the model, ␣-helices H1 and H7 are dominant in forming the cavity, followed by H4 and H10, and together these four ␣-helices accounted for two-thirds of the critical residues. The third round of mutagenesis therefore focused mainly on these ␣-helices and sought to confirm their orientation and to better define the extent of the folate binding site within the cavity (Table 3). Of the 12 residues mutated (all to Cys), only three were critical (supplemental Fig. S5).
Prediction and Testing of Folate Transport Activity of Other Arabidopsis FBT Proteins-The above analysis identified seven residues that are not only important to Slr0642 function but also

Residues of Synechocystis Slr0642 mutagenized in round 2
Strictly conserved residues are indicated by a plus sign, conservatively replaced residues by "c", and non-conserved residues by a minus sign.

Conservation pattern in proven folate carriers Critical Slr0642
At2g32040 FT1 FT5

Residues of Synechocystis Slr0642 mutagenized in round 3
Strictly conserved residues are indicated by a plus sign, conservatively replaced residues by "c", and non-conserved residues by a minus sign.

Conservation pattern in proven folate carriers Critical Slr0642
At2g32040 FT1 FT5 Leu-54 strictly conserved in known folate transporters as well as in cyanobacterial and plant orthologs of Slr0642 (Tables 1 and 2). These residues are Trp-82, Lys-85, Asp-93, Asp-145, Phe-267, Asp-348, and Cys-371. The absence of any of the seven might thus suggest lack of folate transport activity. Inspection of the aligned sequences of the eight unstudied Arabidopsis FBT proteins (15) showed that five of them (At5g10820, At5g54860, At2g33280, At1g04570, and At1g64890) lack up to four of the seven residues and that only three (At5g25040 and At5g25050, which are 88% similar, and At1g79710) have the whole set. Folate transport activity was tested for five of these proteins (At5g25050, At1g79710, At5g10820, At5g54860, and At1g04570) by expressing them in the pabA abgT strain and monitoring growth on 5-FTHF medium. Two types of expression construct were used: native sequences and sequences truncated before the first predicted transmembrane ␣-helix and fused to the N-terminal 27 residues of Slr0642. A strategy similar to the latter was used to express At2g32040 in E. coli (15). Neither version of any of the five proteins conferred growth on 5-FTHF, although positive controls with At2g32040 did so (data not shown). Negative results were also obtained from functional complementation tests in Leishmania folate and biopterin transport mutants with At5g25040, At5g25050, At1g79710, At5g10820, and At5g54860 proteins (data not shown). It should be noted, however, that constructs of the known folate transporter At2g32040 gave no complementation in the Leishmania tests, although positive Leishmania FT1 and BT1 controls behaved as expected.
As an indirect way to further test Arabidopsis FBT proteins for folate transport activity, we screened three cDNA expression libraries for functional complementation of the pabA abgT mutant on 5-FHTF medium. The only folate transporter found by screening 10 7 recombinants from each library was At2g32040, which was recovered three times from the same library. All three clones were truncated at residue 43, close to the predicted start of the mature protein.

DISCUSSION
This study was made possible by the development of a robust E. coli expression system for FBT carriers, comprising a pabA  JANUARY 22, 2010 • VOLUME 285 • NUMBER 4 abgT double mutant host strain, and the medium copy, tightly IPTG-regulated pLOI707HE plasmid (43). This expression system could potentially be used for other FBT proteins and has already been used for another class of folate transporter (46). When used with a folate concentration that was subsaturating for growth, this system enabled a sensitive growth assay to detect transport defects in mutant proteins.

Functional Analysis of an FBT Transporter
It should be said that extensive efforts were made to express Slr0642 in E. coli, Lactococcus lactis, or yeast at levels sufficient for transport assays with 3 H-labeled folates. These efforts included recoding the gene, using specialized membrane protein expression strains, and modifying the N terminus to enhance membrane insertion. Protein expression in all systems was detectable immunologically, but none gave enough transport activity to measure on a time scale less than that of growth. We therefore measured growth itself.
In the course of this work, we tested other members of the Arabidopsis FBT family in our E. coli expression system and also in a heterologous eukaryotic system, employing Leishmania folate or pterin transporter mutants where only a low level of active protein expression is required to relieve folate or pterin auxotrophy. Although positive and negative controls were successful, no evidence of activity could be found in either system. These negative findings may have technical explanations, but other explanations cannot be excluded. Perhaps these transporters act on other substrates or function in compartments incompatible with the heterologous systems used. Further studies will be required to address these possibilities.
Our comparative structural model of Slr0642 is necessarily speculative due to the lack of close homology between FBT proteins and other MFS proteins whose three-dimensional structures have been determined. However, criteria beyond homology were carefully combined to guide model development, and the final result is quite consistent with the pattern of residue conservation and with the results of site-directed mutagenesis. Moreover, our mutagenesis data are consistent with the results of a recent study in which 10 conserved polar residues in the Leishmania FT1 transporter were changed to Leu or Val (47). Of the six residues common to both studies, five (Lys-85, Asp-93, Asp-145, Arg-309, and Asp-348 in Slr0642) were found to be essential or important to transport function in both. As summarized in Fig. 4 (and its interactive version, supplemental Fig. S6, A-C), the majority of critical residues line the predicted central cavity and are concentrated in ␣-helices H1, H4, H7, and H10, which form the cavity's core.
Slr0642 most probably operates via an alternate access mechanism in which the transporter has a central substrate binding site and two gates on either side of the central cavity that regulate access to the binding site in an alternating fashion. Since the model is only an approximation of the real structure, it is not possible to carry out substrate docking or molecular dynamics studies to obtain the binding conformation of the substrate. However, the substrate binding site is likely to be in the middle of the central cavity, corresponding to the area for binding of galactopyranosylthiogalactopyranoside in the structure of LacY (19). Several conserved residues that are found in the central cavity are critical for transport: Arg-53, Lys-85, Asp-145, Asp-149, Phe-252, Thr-259, Glu-263, Gln-357, Phe-380, Met-384, and Asn-388. The interactions of folate with the side chains of amino acids were inspected (using the Protein Data Bank) to understand the chemistry of folate binding (supplemental Fig. S3). Generally, the carboxyl groups of the glutamate moiety interact with Lys or Arg, the N3 and amine group of the pterin moiety interact with Asp or Glu, the ring structure of the pterin moiety interacts with the aromatic Phe or Tyr or hydrophobic Val or Ala, and the pABA interacts with the hydrophobic Leu or Ile. These interactions are to be expected based on the chemical properties of the substrate and side chains. On this basis, Asp-145 and/or Asp-149 could be interacting with N3 and the amine group of the pterin moiety, and Arg-53 and/or Lys-85 could have electrostatic interactions with the carboxyl groups of the glutamate group of folate. The aromatic Phe-252 and/or Phe-380 could have aromatic stacking interactions with the pterin ring, and the hydrophobic residue Met-384 could interact with the pABA or the pterin ring. The polar residues Thr-259, Gln-357, and Asn-388 could be interacting with the amino, carboxylate, or keto groups.
Several other residues could be involved in the gating aspects of the mechanism, and they line the interface between the two domains. On the cytoplasmic side are residues Glu-263 and Arg-53, which might form a salt bridge interaction, linking the two domains. Other residues in the interface are Lys-61, Asp-62, and Phe-267, which have no obvious interaction partners. Facing the outside of the cell is residue Asp-93, which is in the interface of the domains and is conserved not just in all FBT proteins (15) but in all other members of the major facilitator superfamily, pointing to a universal catalytic or structural role. In addition, there are Arg-369 and Phe-380, which could be important in gating. Trp-82 and Arg-309 may be closing the gaps between the two domains of the transporter, whereas the roles of Arg-155, Thr-259, Cys-371, Arg-309, and Asp-348 are unclear. Arg-335 is a conserved residue present in a long loop that does not exist in LacY.
Finally, the identification of conserved, transport-critical residues is a step toward prediction of substrate preference for FBT proteins. In this respect, it is encouraging that of the seven residues that are essential in Slr0642 and strictly conserved in known folate transporters, two are missing from the BT1 transporter of Leishmania, which strongly prefers biopterin (22,24).