|Home | About | Journals | Submit | Contact Us | Français|
Nucleotide biosynthesis pathways have been reported to be essential in some protozoan pathogens. Hence, we evaluated the essentiality of one enzyme in the pyrimidine biosyn-thetic pathway, dihydroorotate dehydrogenase (DHODH) from the eukaryotic parasite Trypanosoma brucei through gene knockdown studies. RNAi knockdown of DHODH expression in bloodstream-form T. brucei did not inhibit growth in normal medium, but profoundly retarded growth in pyrimidine-depleted media or in the presence of the known pyrimidine uptake antagonist 5-fluoruracil (5-FU). These results have significant implications for the development of therapeutics to combat T. brucei infection. Specifically, a combination therapy including a T. brucei-specific DHODH inhibitor plus 5-FU may prove to be an effective therapeutic strategy. We also show that this trypanosomal enzyme is inhibited by known inhibitors of bacterial Class 1A DHODH, in distinction to the sensitivity of DHODH from human and other higher eukaryotes. This selectivity is supported by the crystal structure of the T. brucei enzyme, which is reported here at a resolution of 1.95 Å. Additional research, guided by the crystal structure described herein, is needed to identify potent inhibitors of T. brucei DHODH.
Pyrimidine biosynthesis is a vital biological process, achieved both by de novo synthesis and by pyrimidine salvage pathways. However, many parasitic protozoa are deficient in pyrimidine salvage (Chaudhary and Roos, 2005) and thus enzymes of the pyrimidine synthesis pathway have been considered potential drug targets. The fourth enzyme of this pathway is dihydroorotate dehydrogenase (DHODH), which oxidizes dihydroorotate to orotate, the first aromatic intermediate in this pathway, with consequent transfer of two electrons to the protein cofactor FMN (Figure 1). Recent experiments to validate DHODH as a drug target in Trypanosoma cruzi showed that insect-stage epimastigotes are dependent on de novo synthesis of pyrimidines (Annoura et al., 2005). Specifically, gene interruption of the three gene loci encoding cytosolic DHODH in T. cruzi led to cell non-viability; addition of pyrimidine nucleosides to the medium did not rescue the DHODH knockdown cells. We therefore undertook a more extensive investigation of the essentiality of the closely homologous gene in T. brucei, the causative agent of African Sleeping Sickness.
DHODHs are categorized into two classes, Class 1 and Class 2, based upon their sequence, localization in the cell, and the electron acceptor used to reoxidize FMN in a second half-reaction (Nielsen, Andersen and Jensen, 1996; Nielsen, Rowland, Larsen and Jensen, 1996). Both classes belong to Pfam sequence family PF01180. Class 1 DHODHs are primarily bacterial. They are cytosolic and further divided into 1A (E.C. 184.108.40.206) and 1B (E.C. 220.127.116.11) depending upon their electron acceptor and oligomeric state. Class 2 DHODHs, including the human and plasmodial enzymes, are monomeric and contain mitochondrial targeting sequences located in the N-terminal domain (Rowland et al., 2000). Furthermore, they are found in the inner membrane of mitochondria (Knecht et al., 1996) or associated with the cytoplasmic membranes in the case of E. coli (Karibian, 1978). Specific inhibitors are known for both classes (Palfey et al., 2001; Baldwin et al., 2005).
We present genomic, biochemical, and structural evidence that TbDHODH is a Class 1A DHODH, and in particular that it can be inhibited by Class 1A specific inhibitors. The structural analysis is based on a 1.95Å resolution crystal structure of the T. brucei enzyme co-crystallized with the reaction product orotate. We also present kinetic, inhibition, biophysical, and RNAi knockdown data for the T. brucei DHODH, which has not previously been characterized. Furthermore we show that RNAi knockdown of T. brucei DHODH is synergistic with the effects of 5-fluorouracil (5-FU), an inhibitor of pyrimidine uptake in T. brucei (de Koning and Jarvis, 1998). This finding provides a rationale to develop DHODH inhibitors that could be used in combination with 5-FU or another inhibitor of pyrmidine uptake to yield a potent regimen. We are encouraged that the structure and inhibition studies offer insight into the stereochemical and functional details of the active site of the trypanosomal enzyme, and may facilitate the development of lead compounds for anti-trypanosomal drug discovery. Because our results on the essentiality of this enzyme are divergent from those reported for the closely related T. cruzi parasite, they also underscore the importance of validating anti-protozoan drug targets separately in each target organism (Wang,1997).
Only one DHODH gene was found in the T. brucei genome (GeneDB Tb927.5.3830). This gene was initially identified as a structural genomics target (TargetDB Tbru015978) by the Structural Genomics of Pathogenic Protozoa consortium (SGPP). At the time this work was initiated no trypanosomal DHODH structures were available, although a preliminary characterization of Leishmania major DHODH had been presented (Feliciano et al., 2006) and preliminary crystallization of the T. cruzi enzyme had been reported (Inaoka et al., 2005). Thus we undertook de novo crystallographic structure determination of the T. brucei DHODH, and also undertook biochemical characterization of this enzyme to confirm the assignment of the trypanosomal enzymes to Class 1A.
The crystal structure of TbDHODH was determined to 1.95Å resolution from a single crystal in space group C2221. Consistent with previously determined structures from this sequence family, the overall fold of the TbDHODH monomer is an αβ barrel with the prosthetic FMN group located near the C-terminus (Figure 2a).
The asymmetric unit of the crystal structure contains four monomers that are arranged as two distinct dimers. The monomers of each dimer are related by two-fold symmetry and the dimer overlays well with the L. lactis 1A dimer. The area of contact at the dimer interface of is 2300 Å2 as analyzed by PDBsum (Laskowski et al., 2005). The interface includes 19 inter-monomer hydrogen bonds and a salt bridge between Glu208 of monomer A and Lys297 of monomer B. This salt bridge is also observed in the L. lactis structure. Lys297 is strictly conserved among Class 1A DHODHs; residue 208 is conserved as either glutamate or aspartate (Fig. S2). Thus, this salt bridge is a characteristic feature of the dimer interface in Class 1A DHODHs.
The oligomeric state of the protein was also examined by analytical ultracentrifugation. TbD-HODH was observed by both protein absorption at 280 nm and FMN absorption at 454 nm. Throughout centrifugation experiments, a single ideal species model consistently yielded an average mass indicative of a dimeric enzyme, while a monomer/dimer model reliably produced sub-micromolar dissociation constants for dimerization (Fig. S1). This biophysical result supports the crystallographic observation of TbDHODH as a homodimer.
The TbDHODH structure contains two cis-peptides, Glu56-Pro57 and Cys193-Val194. These cis-peptides are also observed in structures of the L. lactis 1A, L. lactis 1B, and human Class 2 homologs. The two residues following the first cis-peptide, Arg58 and Tyr59, are in close proximity to the active site and are highly conserved among Class 1A DHODH sequences. No other exceptional stereochemistry is observed in the refined TbDHODH model.
The FMN prosthetic group is located near the dimer interface in a pocket formed predominantly by strands β4, β6–7, β11, and β17. FMN forms hydrogen bonds with residues Ala20, Ser45, Lys44, Cys249, Val194, and Lys165 (Fig. S3); and hydrophobic interactions with residues Ala19, Gly21, Met70, Tyr59, Asn68, and with orotate. The isoalloxazine ring of FMN and the aromatic ring of orotate are parallel to each other as in other DHODH-orotate structures. The dimethylbenzene moiety of the isoalloxazine ring system lacks close contacts to any residue and is solvent accessible. This differs from the FMN binding site in the human structure (25% sequence identity), where the Class 2-specific N-terminal domain blocks access to the active site. There are similar hydrogen bonding interactions to the FMN in both the TbDHODH and human structures and these residues are conserved between the two structures except for Lys44 and Cys249. The sulfhydryl of Cys249 forms a hydrogen bond to O3 of the ribityl group of FMN. The equivalent interaction is absent in the human structure, where there is only one hydrogen bond to O3. In a sequence-based structure alignment, Lys44 is a glycine in the human structure but there is a nearby lysine (hLys100), which forms hydrogen bonds to O4 of the isoalloxazine ring and plays the same role in binding to FMN.
Of the five residues in TbDHODH that form hydrophobic interactions to FMN, two are different in the human protein. These differences may be of potential use in designing an inhibitor targeting the FMN binding site with differential specificity for the human and trypanosomal enzymes.
Class 1 and Class 2 DHODHs differ at the active site base responsible for the abstraction of a proton from the C5 of orotate. In Class 1A DHODHs the base is a cysteine; in Class 2 enzymes, including human DHODH, it is a serine (Nørager et al., 2003). In TbDHODH, Sγ of the active site base Cys131 is 3.65 Å away from C5 of orotate. This distance is slightly longer than the 3.3 Å observed in the L. lactis DHODH:orotate structure. In the present crystal structure at pH 4 the cysteine sulfur is most likely protonated. Residues 126–143, spanning the region to either side of the active site Cys, are strongly conserved among all Class 1A DHODHs (Fig. S2). These residues can be described as a pyrimidine binding motif: EZNLSCPNZPUKZQXXYD, where U is a polar uncharged residue, X is a non-specific residue and Z is a hydrophobic residue.
A structure-based sequence alignment of the T. brucei, L. lactis 1B, and human DHODHs, representing Classes 1A, 1B and 2, respectively, shows that the Class 1 sequences have an insertion of 15 residues (residues 202–216) relative to the human Class 2 sequence (Figure 3). These residues form a strand-loop-strand-loop (β14-loop-β15-loop) structural unit. Residues 195–203 preceding the insertion are conserved among Class 1A DHODHs and form a NSXGNGLXI motif (1A motif 1) in close proximity to the active site. Asn195 is one of four conserved Asn residues that interact with orotate.
The orotate binding site is flanked by a loop consisting of residues 132–138 (Fig. S4). This loop is highly conserved in Class 1A DHODHs (Fig. S2). In crystal structures of human DHODH the corresponding loop is poorly ordered, and loop flexibility is hypothesized to allow the active site to open so that dihydroorotate can enter and orotate can exit (Liu et al., 2000). In contrast to the human structure, the loop in the present TbDHODH structure is in a relatively closed conformation and is well defined by electron density. Nevertheless, the orotate and FMN binding sites are more accessible in TbDHODH than in the human structure due to the absence of the N-terminal domain characteristic of Class 2 DHODH (Figure 2). This difference in accessibility may contribute to the differential sensitivity of the Class 1A and Class 2 enzymes to orotate analogues such as hydroxybenzoate.
Gene expression knockdown of the T. brucei DHODH was achieved by over-expressing a stem-loop RNA of the target gene in bloodstream form T. brucei. Messenger RNA levels were decreased by 95% at 72 hours following induction of the RNA interference construct by the addition of tetracycline (Figure 4a). Knockdown of T. brucei DHODH protein appeared to be similarly substantial as observed on Western blots. This was monitored over seven days, with protein levels detected at <5% of levels observed without RNAi induction (Figure 4b). When the T. brucei cells were grown in standard medium (HMI-9), growth rates for cells when RNAi was turned on (TCN+) or turned off (TCN−) were comparable (Figure 5a). The growth rates of these cultures were virtually the same as the control cells (’single-marker’) containing only the T7 RNA polymerase and Tet repressor genes maintained on a neomycin marker. The surprising observation of normal cell growth led us to investigate if the parasites were surviving by scavenging pyrimidines from the medium. Note that HMI-9 medium includes 16 µM thymidine. When the additional thymidine was omitted from the medium, we again saw no change in the growth phenotype upon induction of RNA interference against the DHODH (data not shown). Since the standard medium was made up with 10% fetal bovine serum, we repeated the experiment using commercially available dialyzed serum, which was depleted of small molecules such as pyrimidines. Under these conditions we observed reduced growth upon induction of the RNAi construct. By day 4 the growth rate was significantly reduced and was arrested by day 6 (Figure 5b). Most of the cells expressing the RNAi construct in the pyrimidine-depleted medium took on a bloated, balloon-shaped morphology by day 5. There was a small degree of growth suppression in the TCN− cells compared to the single-marker controls, which may have been due to leaky expression of the RNAi construct.
We therefore proceeded to test whether RNAi knockdown of DHODH would impair growth in standard medium if an inhibitor of pyrimidine uptake were also present. 5-fluorouracil (5-FU) was tested since it has been shown to be a potent inhibitor of pyrimidine uptake by T. brucei (Gudin et al., 2006). The inhibitory concentration of 5-FU causing 50% growth inhibition (IC50) was shifted 16-fold lower in parasites undergoing knockdown of DHODH (IC50 = 0.64 µM) compared to parasites not undergoing knockdown (IC50 = 10.4 µM) using standard medium (Figure 6). In follow up to this finding, we tested the uptake of 3H-5-FU into the DNA of T. brucei cells and observed that uptake was more than 2-fold higher in the DHODH-knockdown parasites than in controls (Figure 7).
TbDHODH activity was examined in the presence of several standard DHODH terminal electron acceptors (Figure 8) and was found to catalyze fumarate most efficiently, typical of Family 1A enzymes (Andersen et al., 1994). The enzyme catalyzes the reduction of menadione (vitamin K3) at ~20% the rate of fumarate catalysis, but does not appear to perform electron transfer to menaquinone (vitamin K2), to the ubiquinone analogues Q0, Q1, Q2, or QD, or to NAD+, further characteristics of a Family 1A DHODH enzyme (Nara et al., 2000; Nielsen, Andersen and Jensen, 1996; Marcinkeviciene et al., 2000). Catalysis of both DHO and fumarate displays saturation kinetics with Km values of 14 µM for DHO and 80 µM for fumarate and a Vmax of 8.5 s−1 (Figure 9). The pH dependence of enzyme activity was assayed (Figure 9 inset) and peak activity was found to occur at pH 7.8, indicative of two ionizable groups with pKa values of 6.3 and 9.2 being involved in substrate binding and catalysis. Analysis of TbDHODH activity across a range of DHO and fumarate concentrations produces a double-reciprocal plot with parallel lines (Figure 10), indicative of a ping-pong mechanism and consistent with results found for other Family 1A enzymes (Marcinkeviciene et al., 2000).
Inhibition of TbDHODH by benzoate pyrimidine analogues was initially performed with 1 mM inhibitor (data not shown). Of the benzoate pyrimidine analogues tested, only 3,4- and 3,5-dihydroxybenzoate displayed inhibition at 1 mM concentration. Subsequent analysis of TbD-HODH inhibition by 3,4- and 3,5-dihydroxybenzoate across various inhibitor and DHO concentrations (Figure 11) revealed competitive inhibition by 3,4-dihydroxybenzoate with a Ki of ~58 µM (determined by IC50 re-plot, not shown) and weaker inhibition by 3,5-dihydroxybenzoate, which appeared non-competitive and yielded a Ki of 200 µM. These values are comparable to that reported for other Class 1A DHODHs (Feliciano et al., 2006; Takashima et al., 2002).
The Trypanosoma brucei genome contains all of the enzymes required for de novo pyrimidine biosynthesis, indicating that T. brucei is capable of providing its own pyrimidine pools (Berriman, 2004). The fourth enzyme in the pathway, DHODH, has been considered a potential drug target in parasitic protozoa including P. falciparum and T. cruzi (McRobert and McConkey, 2002; Annoura et al., 2005; Baldwin et al., 2005). In particular, the Class 1A DHODH from T. cruzi has been shown by gene disruption experiments to be necessary for growth of epimastigotes (Annoura et al., 2005). We performed functional, enzymatic, and structural studies of this protein in the related parasite T. brucei, the causative agent of African sleeping sickness, in an attempt to confirm similar essentiality.
Functional studies using RNA interference to block TbDHODH revealed that bloodstream T. brucei are able to grow normally without this enzyme, apparently by scavenging pyrimidines from the medium (Figure 5). This ability to scavenge pyrimidines is not entirely unexpected, as it has been shown that T. brucei has surface transporters capable of importing uracil and cytosine (de Koning and Jarvis, 1998; Papageorgiou et al., 2005; Gudin et al., 2006). In addition, the T. brucei genome contains at least some of the enzymes associated with pyrimidine salvage, including cytidine deam-inase (GeneDB Tb09.160.1680), uracil phosphoribosyltransferase (GeneDB Tb927.4.3320), and thymidine kinase (GeneDB Tb10.70.7270).
Thus, for T. brucei, a chemotherapeutic strategy to deprive the cells of pyrimidine pools would appear to require both inhibition of pyrimidine biosynthesis and of pyrimidine uptake/salvage. We tested the feasibility of such a two-pronged attack by repeating the DHODH RNAi knockdown experiments in the presence of 5-FU, a compound previously demonstrated to inhibit pyrimidine uptake by T. brucei (Gudin et al., 2006). 5-FU is itself cytotoxic at high concentrations due to inhibition of thymidylate synthase and incorporation into DNA and RNA (de Bono and Twelves, 2001). We observed that cells undergoing this dual treatment were indeed affected by DHODH knockdown, and were 16-fold more sensitive to 5-FU than control cells (Figure 6). This suggests that the cells could not tolerate both inhibition of pyrimidine synthesis and inhibition of pyrimidine uptake.
Although 5-FU inhibits pyrimidine uptake in T. brucei, it is itself transported into the cells with relatively high efficiency (Gudin et al., 2006). A possible cellular response to the suppression of pyrimidine biosynthesis, in this case through interference with DHODH, is to up-regulate the import pathway. Thus a follow-up experiment was performed to determine the relative uptake of radiolabelled 5-FU into the DNA of the T. brucei cells in the presence or absence of DHODH suppression. The experiment demonstrated >2x higher 5-FU incorporation into DHODH knockdown cells than into controls (Figure 7). It seems likely that the increased incorporation of 5-FU into DNA could partly account for the enhanced sensitivity to this drug. We speculate that the increased incorporation of 5-FU into the DNA of DHODH knockdown cells results from one or more of the following factors: 1) the depletion of intracellular pyrimidine pools favors the incorporation of the 5-FU, 2) possible up-regulation of enzymes that convert 5-FU to 5-fluorouracil monophosphate (FUMP), 5-fluorouracil diphosphate (FUDP), or 5-fluorouracil triphosphate (FUTP), or 3) possible increased uptake of 5-FU due to up-regulated pyrimidine transporter activity at the cell surface. Future experiments will be necessary to better establish if any of these mechanisms are responsible for our observations.
Whichever of the above mechanisms may be in play, the apparent synergy of DHODH knockdown plus 5-FU treatment has significant implications for pursuing the pyrimidine biosynthesis pathway as a drug target in T. brucei. One can imagine a combination treatment including 5-FU and a new compound that specifically inhibits the T. brucei DHODH (or possibly another enzyme in the pyrimidine biosynthesis pathway). A new DHODH inhibitor would thus be a “chemosensitizer” for 5-FU and could dramatically enhance the activity of 5-FU for treating T. brucei infection. 5-FU is a well established clinical drug that has been in use since the 1950’s. It has known toxicities, but certainly has many advantages over existing human African trypanosomiasis drugs (e.g. melarso-prol, suramin, pentamidine, etc.) that are extremely toxic or require parenteral administration. It would be important to develop a partner drug that inhibits the T. brucei DHODH selectively. Since the human DHODH is a Class 2 enzyme and thus functionally and structurally distinct from the Class 1 T. brucei enzyme, we expect that finding selective inhibitors is quite feasible. The data presented here provide a strong rationale for further pursuing research into the T. brucei DHODH in order to develop potent inhibitors as drug leads.
The previous focus of inhibitor studies has largely been directed to Class 2 DHODHs. The P. falciparum enzyme (PfDHODH) is the focus of significant interest for its potential as a new target for the development of anti-malarials (Baldwin et al., 2005), and the human enzyme is targeted for treatment of cancer, gastric ulcers, and rheumatoid arthritis (Palfey et al., 2001). Structural studies of PfDHODH inhibitors based upon brequinar analogues (Boa et al., 2005), biphenylamides (Baldwin et al., 2005), A771726 (Heikkilä et al., 2006), atovaquone (Hansen et al., 2004), and lefluno-mide (Hurt et al., 2006) show that all these analogues bind similarly, attacking the ubiquinone binding site formed jointly by the N- and C-terminal domains (Figure 2c). Unfortunately, given the involvement of the Class 2-specific N-terminal domain, this class of inhibitors is unlikely to be relevant to Class 1 enzymes.
However, inhibitors specific to Class 1A enzymes are also known, including hydroxybenzoates. These are orotate analogs that have been shown to specifically inhibit the L. lactis Class 1A DHODH with minimal effect on the L. lactis Class 1B enzyme or the human and E. coli Class 2 enzymes (Palfey et al., 2001). Although the stereochemistry of the protein residues forming the orotate binding site Class 1A and Class 2 enzymes are nearly identical, the redox mechanism differs in the two classes (Fagan et al., 2007). Fagan et al. (2007) have attributed this to subtle dynamic features of the proteins. On the basis of the present structure, we suggest that these features include differential access to the orotate and FMN sites from the "front" side of the protein (Figure 2) and a difference in the relative mobility observed crystallographically for residues 132–138 in the pyr-binding motif (Figure 3).
Wolfe et al. (2007) further investigated the structural details and kinetics of hydroxybenzoate inhibition of the L. lactis Class 1A DHODH. The inhibitor was observed crystallographically only in one monomer of the homodimer, consistent with kinetic data showing cooperativity between the two monomers of the dimer. The structural basis for this cooperativity is not evident (Wolfe et al., 2007). Given the high structural conservation between the L. lactis structure and the present T. brucei structure, it seemed plausible that this class of inhibitors was also relevant to trypanosomal homologs. We have now shown through inhibition assays on TbDHODH that hydroxybenzoate inhibitors act indeed not only on bacterial DHODHs but also on this Class 1A DHODH from a eukaryotic parasite. In particular, we found competitive inhibition by 3,4-dihydroxybenzoate (Ki of 58 µM) and non-competitive inhibition by 3,5-dihydroxybenzoate (Ki of 200 µM). These Ki values are higher than those reported for L. lactis but remain possible starting points for anti-parasitic drug design.
T. brucei is a representative of a larger class of pathogenic protozoan parasites responsible for a substantial burden of human disease. In evaluating the importance of TbDHODH as a possible drug target, it is therefore of interest whether drugs acting on TbDHODH would also be effective against related pathogens. Class 1A DHODHs are primarily bacterial, although they have been acquired by some yeasts (Saccharomyces, Kluyveromyces) apparently through lateral gene transfer from Lactococcus (Hall et al., 2005). The appearance of Class 1A DHODH also in trypanosomatids may similarly be due to lateral transfer, but the point of introduction into the clade of protozoan lineages is not known.
The sequences of T. brucei Class 1A, L. lactis Class 1B, and P. falciparum Class 2 DHODH were submitted to the GeneDB omni-BLAST server (Hertz-Fowler et al., 2004) to find homologs in other protozoa. No species was found to contain more than one DHODH class. No protozoan Class 1B homologs were detected. The Plasmodium genomes contain only Class 2 sequences. The sequences from T. brucei, T.b. gambiense, T. cruzi, T. congolense, T. vivax, L. major, and L. infantum have >75% sequence identity to the genome strain TbDHODH. The residues involved in binding FMN and orotate are conserved, suggesting that drug design based on the TbDHODH structure is relevant to other pathogenic trypanosomatids.
We have shown that RNAi interference with dihydroorotate dehydrogenase in bloodstream form Trypanosoma brucei inhibits growth when access to pyrimidines in the environment is also blocked. This calls into question the general essentiality of the pyrimidine biosynthesis pathway in T. brucei, in distinction to a previous report of essentiality in T. cruzi. However, we have also shown that RNAi suppression of DHODH is synergistic with the effect of the complementary growth-blocking agent 5-FU, which inhibits pyrimidine uptake. This demonstrates that interference with T. brucei DHODH can potentiate the activity of the established clinical drug 5-FU against bloodstream form T. brucei.
The crystal structure and enzymatic properties of the T. brucei enzyme show that its active site is more similar to those of other Class 1A DHODHs than to the Class 2 DHODHs of higher eukaryotes, including human DHODH. Encouragingly, it is inhibited by hydroxybenzoates, in common with other Class 1A DHODHs but in distinction to the Class 2 human DHODH.
Thus T. brucei DHODH, and by extension other enzymes on the pyrimidine biosynthetic pathway, remains a possible target for the development of a two-pronged therapeutic approach attacking both pyrimidine biosynthesis and pyrimidine uptake. Furthermore, we propose that the molecular details provided by the present TbDHODH structure may be useful in designing selective inhibitors and potential lead compounds targeting other protozoan dihydroorotate dehydrogenases.
The gene encoding TbDHODH was PCR amplified from genomic DNA of the strain T. brucei TREU927 GUTat 10.1 and cloned into a vector AVA421, which is based on E. coli expression vector pET14b (Alexandrov et al., 2004). The AVA421 vector contains a cleavable N-terminal His-tag. The protein was purified using a Ni-NTA column and the bound protein was cleaved by protease 3C overnight at 4°C. The released protein was further purified by gel filtration on a HiLoad Superdex 200.
The DHODH gene was amplified from T. brucei brucei 427 strain genomic DNA by polymerase chain reaction using the primers 5’- CGCAAAGGGAAAAGTCATATGAGTTTGAAAGTTAAT-3’ and 5’-CATCAACACCCATGTATCTCGAGATCCATCGTCTTCACGC-3’, which introduced NdeI and XhoI restriction sites. The product was subcloned, then ligated into pET22b (EMD Biosciences) providing a C-terminus His tag. The resulting TbDHODH-pET22b plasmid was then transformed into BL21-DE3 E. coli cells, and grown in LB-Miller medium to and OD600nm of 0.8, then supplemented with 100 µM FMN and induced with 200 µM IPTG. A cell pellet was collected 3 hours after induction and resuspended in lysis buffer (50 mM sodium phosphate, pH 8.0, 10 mM imidazole, 5 mM 2-mercaptoethanol, 10% glycerol, 1 mg/ml lysozyme, 0.1 mg/ml DNAse) and a protease inhibitor cocktail was added with a final concentration of 200 µM phenylmethylsulfonyl fluoride, 1 µg/ml leupeptin, 2 µg/ml antipain, 10 µg/ml benzamidine, 1 µg/ml pepstatin and 1 µg/ml chymostatin. The lysis mixture was then incubated for 2 hours on ice with constant stirring, sonicated and centrifuged at 20,000 × g for 30 minutes at 4°C. The supernatant from centrifugation was loaded onto a 10 ml column of nickel agarose resin (Qiagen) equilibrated with Wash Buffer (50 mM sodium phosphate, pH 8.0, 300 mM sodium chloride, 20 mM imidazole, 5 mM 2-mercaptoethanol, 10% glycerol). The column was washed with ten volumes of Wash Buffer, then His6-tagged protein was eluted with Elution Buffer (50 mM sodium phosphate, pH 8.0, 300 mM sodium chloride, 300 mM imidazole, 5 mM 2-mercaptoethanol, 10% glycerol). The eluted protein was pooled and concentrated to less than 2 ml with a Centriprep YM-10 centrifugal filtration device (Amicon Millipore) and applied to a HiPrep 26/10 Desalting Column (Amersham Biosciences) equilibrated with 50 mM sodium phosphate, pH 8.0, 150 mM sodium chloride, and 10% glycerol. Protein containing fractions were pooled and concentrated as above.
Purified TbDHODH was shipped on dry ice for High Throughput Crystallization Screening at the Hauptman-Woodward Medical Research Institute (Luft et al., 2003). The sample was rapidly thawed in a 30°C water bath and orotic acid (8 µl of a 200mM orotic acid stock solution) was added to 392 µl of the protein solution. Crystallization experiments were set up using the microbatch-under-oil technique (Chayen et al., 1992) with a Robbins Scientific Tango liquid handling system. Each of the 1536 experiments contained 200 nl of crystallization cocktail solution combined with 200 nl of protein solution under USP grade mineral oil (Sigma cat #M1180) contained in a low birefringence 1536 well microassay plate (Greiner BioOne cat# 790801). Manual review of plate images after 1 week at 4°C identified 124 crystallization conditions with outcomes considered suitable for optimization trials. Eight initial hits were optimized and crystals grown by the sitting drop method. Crystals used for data collection grew from a drop consisting of 0.4 µl of protein solution (14.1 mg/ml), mixed with 0.4 µl of reservoir solution containing 40% PEG 1000, 0.1M Na citrate, 0.1M KBr, pH 4.0 at room temperature.
Native crystals were flash frozen directly in liquid nitrogen. Data were collected to 1.95Å resolution from a single crystal maintained at 100K, at an X-ray wavelength corresponding to the bromine K edge (λ = 0.92010Å), on ALS beamline 8.2.1. Data were processed and scaled using MOSFLM and the CCP4 program package (Collaborative Computational Project No. 4, 1994) via the automated Wedger and Scaler scripts of Elves (Holton and Alber, 2004). Data processing statistics are shown in Table 1. The contribution of anomalous dispersion from ordered Br atoms to the measured Bragg intensities was insufficient to solve the structure using SAD phasing, although 11 partially-occupied sites were later identified in the refined structure. Therefore the structure of L. lactis DHODH (PDB accession code 1jue, 51% sequence identity) was submitted to the Swiss-model server (Guex and Peitsch, 1997; Schwede et al., 2003) to create a hybrid homology model based upon the structure of L. lactis DHODH but containing the TbDHODH residues. This model was submitted to Molrep (Vagin and Teplyakov, 1997), which placed four copies of the monomeric model into the a.s.u. yielding a correlation coefficient 0.60. The TbDHODH structure was refined using Refmac5 (Murshudov et al., 1997) against all reflections between 21 Å and 1.95 Å resolution, keeping 5% of the reflections for Rf ree. Model building and validation were accomplished using Coot (Emsley and Cowtan, 2004) and Molprobity (Lovell et al., 2003). Electron density in the region of the orotate binding site is shown in Fig. S5. The final model consists of four molecules forming two dimers. The final model for each monomer consists of residues 1–313, one orotate, and one FMN molecule. The overall model also contains 11 bromide atoms, 3 glycerol molecules and 638 water molecules. The final model and crystallographic structure factors have been deposited with the Protein Data Bank (accession code 2b4g).
Enzyme kinetic assays were performed in enzyme assay buffer (50 mM HEPES, pH 8.0, 150 mM NaCl, 10% glycerol) at 25°C. Analysis of pH dependence utilized MES (pH 5.5–6.5), MOPS (6.5–7.5), Tris (7.5–9.0), or CHES (9.0–10.0) buffer salts in place of HEPES. Production of orotate was monitored at 300 nm (ε300 = 2.65 mM−1cm−1) when using fumarate, at 282 nm (ε282 = 6.60 mM−1cm−1) when using vitamin K analogues, and at 296 nm (ε296 = 4.30 mM−1cm−1) when using ubiquinone analogues. Reduction of NAD+ was monitored at 340 nm (ε340 = 6.22 mM−1cm−1). Strong absorption in the far-UV by the inhibitors tested required the use of the traditional electron acceptor dichlorophenolindophenol (DCIP) in inhibition studies where reduction of DCIP was followed at 600 nm (ε600 = 18.8 mM−1cm−1).
DHO-dependent activity was examined in the presence of 500 µM fumarate. Fumarate-dependent activity was assayed at a constant DHO concentration of 250 µM. Data were analyzed using the software package GraphPad Prism and SigmaPlot and fitted to the Michaelis-Menten equation.
Alternative electron acceptors (vitamin K and ubiquinone analogues, NAD+) were assayed at 100 µM electron acceptor (dissolved oxygen was estimated to be ~210 µM) and 250 µM DHO. Inhibitors were initially examined at 1 mM in the presence of 25 µM DHO and 120 µM DCIP, and subsequent assays utilized 120 µM DCIP and various DHO and inhibitor concentrations. All experiments were performed in at least duplicate; combined data are presented, and kinetic parameters are reported as the mean ± standard error of the fit.
Analytical ultracentrifugation was performed on a Beckman XL-I analytical ultracentrifuge and enzyme was observed at 280 nm (protein) and 454 nm (FMN). Enzyme samples of 5–25 µM were centrifuged at 14000 and 20000 rpm at 4°C, equilibrium sedimentation scans were obtained every two hours, and equilibrium was determined to be reached when contiguous scans were indistinguishable. Data were fitted to both a single ideal species model and a monomer-dimer model and one representative result from multiple experiments is presented.
A region of the TbDHODH gene sequence was selected for RNAi using the program RNAit (Redmond et al., 2003). Bases 413–859 of the gene were amplified from T. brucei 927 genomic DNA using the primers 5’-ATACCAATGTGATGGCACAGGTCGGATATGACTTTGA-3’ and 5’-ATACCATAGAGTTGGGTTGAGTCGAGCAAATATGATG-3’. The amplicon was ligated using TA cloning into the vector pGEM-T (Promega, cat# A3600), then excised with the enzyme BstXI. This insert was ligated into the stem-loop RNAi vector, pQuadra3, as described elsewhere (Inoue et al., 2005). The construct was sequenced to verify the identity of the insert, then linearized with NotI in preparation for electroporation.
T. brucei bloodstream-form parasites expressing the T7 RNA polymerase and the Tet repressor under a single selection marker were provided by G. Cross (Rockefeller University) (Wirtz et al., 1999). The cells were grown in HMI-9 with 10% heat-inactivated fetal bovine serum at 37°C and 5% CO2 with G418 at 2.5 µg/ml (Hirumi and Hirumi, 1989). Mid-log phase T. brucei (2.5 × 107) were suspended in 500 µl of cytomix (120mM KCl, 0.15 mM CaCl2, 10 mM K2HPO4/KH2PO4 pH 7.6, 25 mM Hepes, 2 mM Na2EDTA, 5 mM MgCl2) containing 10 µg of NotI linearized pQuadra3/DHODH DNA. The mixture was electroporated in a 4mm gap cuvette with 1.6 kV and 24Ω resistance. Cells were re-suspended in HMI-9 medium in the presence of 2.5 µg/ml phleomycin and 2.5 µg/ml G418. Individual clones were selected for subsequent RNAi studies. The expression of dsRNA was induced by addition of 1 µg/ml tetracycline to the cultures diluted to 1 × 105 cells/ml. Cultures were passed at a 1:10–1:20 dilution daily and cell concentrations were monitored using a ATPLite Luminescence ATP detection Assay System (PerkinElmer, cat #6016941).
Separate growth experiments were performed using HMI-9 medium that omitted thymidine (normally included at 16 µM). There are no other pyrimidine bases in the medium; it is supplemented with 100 µM hypoxanthine as a source of purines. The medium was supplemented with 10% dialyzed fetal bovine serum (Hyclone) to complete the pyrimidine depleted medium.
RNA was isolated from induced and non-induced cultures after 72 h of growth using a RNEasy kit (Qiagen). 16 µg of total RNA was electrophoresed on a formaldehyde gel and blotted using standard procedures. The RNA membrane was then analyzed with a DNA probe amplified from genomic DNA using the primers 5’-CTCGGACATGAGTTTTCCAACC-3’ and 5’-GGCTTCCTACTGTAATCGTG 3’ that correspond to bases 22–281 of the DHODH open reading frame. The blot was stripped and re-probed with DNA of the T. brucei β-tubulin gene for standardization.
Sensitivity to 5-FU was tested against trypomastigotes of the single marker parent line and the single marker cells containing the pQuadra3/DHODH RNAi vector in both uninduced cells and cells induced for 48 hours with 1 µg/ml tetracycline. Cells were added to 96-well plates in triplicate with an initial inoculum of 2 × 104 trypomastigotes per well. The compound stock solution of 5-fluorouracil (Sigma, F6627) was prepared using sterile water at an initial concentration of 20 mM. This was diluted in HMI-9 to a concentration of 120 µM and added in serial dilutions for a final volume of 200ul/well. Parasite growth was quantified at 48 hours by the addition of Alamar Blue (Alamar Biosciences, DAL1100) (Raz et al., 1997).
Incorporation of 5-FU into T. brucei genomic DNA was measured for the bloodstream forms of the single marker parent cell line and for the bloodstream forms containing the pQuadra3/DHODH RNAi construct. Measurements were taken for both uninduced cultures and cultures that had been induced for 72 hours with 1 µg/ml tetracycline. The trypomastigotes were all grown in standard HMI-9 containing thymidine and hypoxanthine, as well as 2.5 µg/ml neomycin. In addition, the cells containing the pQuadra3/DHODH RNAi vector were also grown in the presence of 2.5 µg/ml phleomycin. After the 72 hour induction, cells were transferred to a 96-well plate and pulsed with 5 µCi of [6-3H]-5-FU (Moravek Biochemicals, MT-686) for 24 hours. Cell density quantification, in order to standardize results, was done with 30 ul of sample using the ATPLite Luminescence ATP detection Assay System (PerkinElmer, cat #6016941) as described above. The incorporation of 3H-5-FU was measured by binding the DNA from the remaining 175 µl onto glass fiber filters with a cell harvester (Inotech Biosystems International, Inc., Rockville, MD), and the radioactivity incorporated into the parasites was counted on a Chameleon 425-104 multilabel plate scintillation counter (Hidex Oy, Turku, Finland).
Purified TbDHODH was sent to Affinity BioReagents (Golden, CO) for antibody production in rabbits as previously described (www.bioreagents.com/index.cfm/fuseaction/content.aod). Antibody specificity was confirmed against purified TbDHODH and T. brucei cells lysates using pre-immune and post-challenge serum and found to be specific. The anti-TbDHODH antibodies were used in Western blots on lysates of T. brucei bloodstream form cells collected before and at daily time intervals after induction of RNAi. Standard Western blot procedures were followed using goat-horseradish peroxidase-conjugated secondary antibody and ECL detection (Amersham Biosciences). Relative TbDHODH protein levels were determined by densitometry. The numerical values were normalized to stained bands from the transferred SDS-PAGE gel (not shown).
We thank Jonathan Caruthers and Meg Holmes for aid in data collection. We are also reliant on the efforts of other members of the SGPP consortium, including Eric Phizicky, Erin Quartley, Angela Lauricella, Lori Schoenfeld, Jennifer Ross, Frank Zucker, and Thomas Earnest. We thank Dr. HP de Koning for helpful discussions about T. brucei transporter activities. This work was supported by the Protein Structure Initiative, award NIGMS GM64655; NIH award AI067921; funds from the Howard Hughes Medical Institute (HHMI) to WGJH; and by the National Institutes of Health Grants AI053680 to MAP. Portions of this work were carried out at the Advanced Light Source, which is supported by the Director, Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.