Characterization of the first beta-class carbonic anhydrase from an arthropod (Drosophila melanogaster) and phylogenetic analysis of beta-class carbonic anhydrases in invertebrates
© Syrjänen et al; licensee BioMed Central Ltd. 2010
Received: 22 February 2010
Accepted: 26 July 2010
Published: 26 July 2010
The β-carbonic anhydrase (CA, EC 22.214.171.124) enzymes have been reported in a variety of organisms, but their existence in animals has been unclear. The purpose of the present study was to perform extensive sequence analysis to show that the β-CAs are present in invertebrates and to clone and characterize a member of this enzyme family from a representative model organism of the animal kingdom, e.g., Drosophila melanogaster.
The novel β-CA gene, here named DmBCA, was identified from FlyBase, and its orthologs were searched and reconstructed from sequence databases, confirming the presence of β-CA sequences in 55 metazoan species. The corresponding recombinant enzyme was produced in Sf9 insect cells, purified, kinetically characterized, and its inhibition was investigated with a series of simple, inorganic anions. Holoenzyme molecular mass was defined by dynamic light scattering analysis and gel filtration, and the results suggested that the holoenzyme is a dimer. Double immunostaining confirmed predictions based on sequence analysis and localized DmBCA protein to mitochondria. The enzyme showed high CO2 hydratase activity, with a kcat of 9.5 × 105 s-1 and a kcat/KM of 1.1 × 108 M-1s-1. DmBCA was appreciably inhibited by the clinically-used sulfonamide acetazolamide, with an inhibition constant of 49 nM. It was moderately inhibited by halides, pseudohalides, hydrogen sulfide, bisulfite and sulfate (KI values of 0.67 - 1.36 mM) and more potently by sulfamide (KI of 0.15 mM). Bicarbonate, nitrate, nitrite and phenylarsonic/boronic acids were much weaker inhibitors (KIs of 26.9 - 43.7 mM).
The Drosophila β-CA represents a highly active mitochondrial enzyme that is a potential model enzyme for anti-parasitic drug development.
Carbonic anhydrases (CAs, EC 126.96.36.199) catalyze the reversible hydration of carbon dioxide according to the following reaction: CO2 + H2O ↔ HCO3- + H+ . CAs are zinc-containing metalloenzymes, except for the ζ form, which uses cadmium [2, 3] as alternative metal cofactor. Additionally γ-CAs contain iron in vivo, at least in anaerobic Archaea [4, 5]. The reaction catalyzed by CAs is crucial in the regulation of acid-base balance in organisms. In addition, CAs participate in many other physiological processes such as bone resorption in vertebrates, gluconeogenesis, production of body fluids, and transport of CO2 and HCO3- to name but a few .
To date, five different classes of CAs have been identified: α, β, γ, δ and ζ . A previously known ε-CA class  has been reclassified as a new type of β-CA based on its crystallographic structure , which shows a fold nearly identical to those of the archaeal cab-type  and plant-type  β-CAs. In ζ-CA, the geometry of the active site is nearly identical to that of β-CAs, and there is also some similarity in the fold, which has led to the suggestion that even ζ-CAs might represent a distantly diverged subtype of β-CAs .
The ζ-CAs are found only in diatoms, and the δ-CAs in diatoms and other marine phytoplankton, whereas the three major classes (α, β and γ) are widely distributed. γ-CA orthologs are present in Archaea, Bacteria and plants [6, 11] but have been secondarily lost in animals and fungi  (M. Tolvanen, unpublished observation). The α class is missing from Archaea (M. Tolvanen, unpublished observation) but is nearly ubiquitously present in Bacteria and Eukarya, with the exception of Fungi, in which only filamentous ascomycetes have α-CAs . In plants and animals, α-CAs exist as numerous isoforms. For example, 16 and 15 α-CAs have been described in non-primate mammals and primates including humans, respectively [1, 13]. β-CAs appear to be the class with the widest distribution as they have been described in multiple lineages throughout the tree of life, including most species belonging to the Archaea and Bacteria domains and probably all species of plants and fungi among Eukarya ( and M. Tolvanen, unpublished data). β-CAs have also been characterized in a number of human pathogens such as fungi/yeasts (e.g., Candida albicans, Candida glabrata, Cryptococcus neoformans and Saccharomyces cerevisiae) [15–20] and bacteria (Helicobacter pylori, Mycobacterium tuberculosis, Haemophilus influenzae, Brucella suis etc.) [21–25], and their inhibition profiles by various agents such as sulfonamides, anions, carboxylates and boronic acids have been explored [15–24, 26–30]. Given that β-CA genes appear to be missing in invertebrates, novel antimicrobial compounds based on the inhibition of β-CAs from pathogenic organisms may soon become available.
Although β-CAs catalyze the same reaction as α-CAs and other CA forms, important structural differences between these classes exist. First, instead of functioning as obligate monomers like most α-forms or trimers like γ-forms, β-CAs are found in many oligomerization states. Crystal structures of dimeric, tetrameric and octameric β-CAs have been reported [9, 10, 31]. In the active site of β-CAs, the zinc atom is coordinated by one histidine and two cysteine residues instead of three histidine residues present in α-, γ- and δ-CAs . Despite this difference, β-class CAs basically share the same molecular mechanism for reversible hydration of carbon dioxide as α-class CAs . β-CAs possess a highly conserved dyad comprising an aspartate and an arginine residue that seem to be crucial for the catalytic mechanism since mutation of these residues severely reduces the catalytic activity of the enzyme . The aspartate makes a hydrogen bond with the Zn(II) coordinated water molecule, activating it for nucleophilic attack of the CO2 molecule .
β-CAs have been reported in many photosynthetic organisms, including plants and algae . In plants, β-class CAs are strongly expressed in both roots and green tissues and are located in chloroplasts, cytoplasm and mitochondria with isozyme-specific patterns [13, 34]. The cytoplasmic and chloroplastic CAs are suggested to be crucial in CO2 accumulation and carbon fixation. The green algae Chlamydomonas reinhardtii also contains β-class CAs that are localized to mitochondria . In fungi like Cryptococcus neoformans and Candida albicans, β-CAs have an important role in CO2 sensing and, consequently, in the pathogenesis of these species . In addition, many fungal β-CAs have been shown to be mitochondrial . These and many other similar findings confirm that β-CAs are physiologically important enzymes with variable localization and function like α-CAs in vertebrates.
The presence of β-CAs in the animal kingdom has been controversial or ignored due to the paucity and poor quality of the available sequences [6, 14]. Here we show, however, that β-CAs are widespread among invertebrates. The aim of this study was to express, purify and characterize a β-CA enzyme from fruit fly (D. melanogaster), a commonly used model organism in biological sciences. The recombinant enzyme was produced in Sf9 insect cells using the baculovirus/insect cell expression system. Along with the characterization of the first arthropod β-CA, this study also describes its inhibition profile with inorganic anions. These results might open new strategies for developing novel anti-parasitic drugs against common diseases like schistosomiasis and malaria.
We found β-CAs in all complete non-chordate animal genomes and in almost all invertebrates with at least 20,000 EST sequences in the NCBI database, plus in some with fewer sequences. We confirmed the existence of β-CA sequences in Placozoa, Cnidaria, Platyhelminthes, Nematoda, Arthropoda, and Annelida, and even in many classes of Deuterostomia, namely in Hemichordata, Echinodermata, and Xenoturbellida. The only major taxon of Protostomia with poor evidence for β-CA is Mollusca, in which we found only one EST to match 55 residues in other β-CAs. Of special interest were numerous pathogenic helminth species with complete β-CA sequences or substantial fragments, namely, filaria-causing Brugia malayi, mouse whipworm Trichuris muris, dog hookworm Ancylostoma caninum, and the flukes Schistosoma mansoni and Schistosoma sinensis.
In the case of chordates, the existence of a functional β-CA is currently unclear. The genomes and NCBI EST sequence collections (as of 1st Oct, 2009) of Ciona intestinalis and Ciona savignyi (tunicates) lack β-CA, whereas we found two recognizable but incomplete β-CA sequences in the genome of the cephalochordate Branchiostoma floridae (of one locus from both haplotypes in the genome). The encoded proteins seem to lack more than 60 residues at the N-terminus, including the active site. In addition, there are two even less complete partial EST transcripts, GenBank BW824885 and BW803919. The latter contains an unrelated sequence in place of the active-site-containing exon. We can only conclude that it remains an open question whether β-CA in B. floridae is a pseudogene or an incompletely sequenced active gene.
Identification of sequences and species in the sequence alignment and phylogenetic tree
Common name and classification
gi or accession numbers
Yellow fever mosquito, insects
157110803 + 77891004
Southern house mosquito, insects
Malaria mosquito, insects
Fruit fly, insects
Fruit fly, insects
Red flour beetle, insects
Pea aphid, insects
Water flea, crustaceans
FE417346 + FE409868
Medical leech, Annelid worms
EY481200 + EY505051 + EY490477
Acorn worm, hemichordate
Caenorhabditis elegans BCA1
Pristionchus pacificus BCA1
FG098717 + GeneWise
Dog hookworm, nematodes
FC551456 + FC550353
Caenorhabditis elegans BCA2
Pristionchus pacificus BCA2
FG106379 + GeneWise
Sea louse, crustaceans
Sea louse, crustaceans
Sea anemone, Cnidaria
Sea anemone, Cnidaria
DB424979 + DB440523
Sea urchin, Echinodermata
139313180 + 139245724
Purple sea urchin, Echinodermata
No β-CA sequences exist in vertebrate genomes. A false positive is found in the X. tropicalis genome, scaffold_1719, which we attribute to bacterial contamination, since the β-CA sequence and every other gene in this scaffold are highly similar (>80% identical) with known genes of Pseudomonas-related bacteria.
All of the complete sequences are classified as β-CA domains of type "BCA_CladeB" in the NCBI Conserved Domain Database. All plant β-CAs are also classified as BCA_CladeB domains, so animal β-CAs can well be said to be "plant-type β-CAs". There are also hundreds of bacterial β-CA sequences of type "BCA_CladeB", and some of them cluster closest to animal sequences, and some closest to plant sequences in phylogenetic trees (data not shown).
Expression of β-CA in Sf9 insect cells
Subcellular localization of DmBCA
When the collection of 38 complete animal β-CA sequences was analyzed for subcellular targeting with TargetP, 22 sequences were predicted to be mitochondrial, with low reliability, and most of the remaining sequences were predicted to be cytoplasmic, again with low reliability. The results from WolF PSORT, Secretome 2.0 and MitoProt II v. 1.101 also supported the hypothesis that it is a mitochondrial enzyme, with a possibility for non-classical secretion (without a signal peptide) suggested by Secretome 2.0. A mitochondrial targeting signal sequence is also predicted in the N-terminus (prediction lengths varying from 14 in MitoProt to 49 in TargetP). Based on these findings, a DmBCA-GFP construct was designed to study the subcellular localization of the protein in Sf9 cells.
Catalytic activity and inhibition of DmBCA
Kinetic parameters for the CO2 hydration reaction catalyzed by various CAs
k cat (s -1 )
k cat /K m (M -1 s -1 )
K I (nM) c
HCA I a
2.0 × 105
5.0 × 107
HCA II a
1.4 × 106
1.5 × 108
3.9 × 105
4.3 × 107
8.0 × 105
9.7 × 107
9.4 × 105
9.8 × 107
3.8 × 105
4.8 × 107
9.5 × 105
1.1 × 108
Inhibition constants of anionic inhibitors against various CAs
K I [mM] #
N 3 -
HCO 3 -
NO 3 -
NO 2 -
HSO 3 -
SO 4 2-
ClO 4 -
H 2 NSO 2 NH 2
H 2 NSO 3 H §
Ph-AsO 3 H 2 §
Dynamic light scattering and gel filtration analysis
The hydrodynamic diameters of the proteins were measured by dynamic light scattering (DLS). At a temperature range from 4°C to 37°C the hydrodynamic diameter of DmBCA was 6.3 ± 0.8 nm and the diameter was found to slightly increase during the heating. According to the hydrodynamic diameter the average molecular weight of DmBCA was estimated to be 48.9 kDa, the lowest estimate being 35.4 kDa, and the highest 64.9 kDa (n = 30). When temperature was raised above 43°C the hydrodynamic diameter of DmBCA started to clearly increase and DmBCA had a transition state at temperature range 46-49°C, where large protein aggregates appeared. This might be associated with protein denaturation.
At the temperature range from 4°C to 37°C the hydrodynamic diameter of HCA II was 3.9 ± 0.6 nm leading to estimated average molecular weight of 16.1 kDa, the lowest estimate being 10.8 kDa, and the highest 22.8 kDa (n = 12). When HCA II was heated from 4°C to 43°C there was a slight increase in hydrodynamic diameter. However, no transition states were observed.
The estimated molecular weights determined by DLS are in agreement with the values obtained by analytical gel filtration which gave a molecular weight estimate for DmBCA of 38.1 ± 3.0 kDa and for HCA II of 23.4 ± 4.0 kDa. The molecular weight of HCA II monomer is 29 kDa. Therefore, gel filtration analysis appears to underestimate the molecular size of CAs. Overall, DLS and gel filtration analyses suggest dimeric state for DmBCA whereas HCA II appears predominantly monomeric in solution.
In the present study, we identified and characterized a novel β-CA enzyme (DmBCA) from an arthropod, D. melanogaster. Although β-CAs have been previously reported in Archaea and Bacteria domains, in addition to plants and fungi in Eukarya, our results suggest the widespread occurrence of at least a single-copy β-CA gene among animal species distinct from chordates. The loss of β-CA gene in the chordate lineage may have occurred either in the last common ancestor of all chordates or in the last common ancestor of tunicates and craniates. Whether cephalochordates have a functional β-CA gene remains an unresolved question.
Using bioinformatic tools, we discovered a single β-CA gene in most invertebrates with reasonable sequence coverage. The exception is nematodes, which seem to have two β-CA genes in their genomes. A very recent article  reported the cloning and characterization of C. elegans β-CAs and the authors found that one of the two isoforms, BCA-1, has no activity and does not work in complementation assay. This might, however, be due to incorrect sequence: the sequence they show for C. elegans BCA-1 contains the coding sequence of the preceding MTP18 gene fused to the β-CA reading frame. This sequence is a result of an incorrect gene prediction which has no support at the mRNA level and which remained in WormBase and UniProt until May 2008, subsequently corrected from our initiative. The previously fused UniProt entry Q8T3C8 now represents only MTP18, and a new entry BCA1_CAEEL contains the more plausible BCA-1 sequence. Since C. elegans BCA-1 has all of the active-site residues and is also well conserved in other nematodes, we think it is probable that both isozymes in nematodes would be functional β-CAs.
Our larger phylogenetic analysis (M. Tolvanen, unpublished data) and NCBI domain classification of animal β-CAs clearly show them to be "plant-type" β-CAs. The "plant-type" group also includes hundreds of bacterial β-CAs, and this group contains plant/bacterial and animal/bacterial subgroups. The presence of such polyphyletic subgroups may indicate horizontal gene transfer, and further investigations are underway to study this phenomenon.
Fasseas et al.  state that the 3D structure of their C. elegans β-CA models would resemble archaeal "cab-type" β-CAs, but we cannot agree with their conclusion. In our eyes, the models are nearly identical to the pea BCA (panel a in their Figure 2), and in our hands SwissModel http://swissmodel.expasy.org/ in fact chooses the pea β-CA structure as a template for both insect and nematode β-CAs.
In the recent paper by Fasseas et al. , the enzyme activity of C. elegans Y116A8C.28 was reported to be typical for β-CAs, with a kcat and kcat/Km of 2.77 × 104 s-1 and 6,383 × 105M-1s-1, respectively. These values are significantly lower than the kcat of 9.5 × 105 s-1, and kcat/Km of 1.1 × 108 M-1s-1 that we report here for DmBCA. The enzymatic activity of DmBCA is one of the highest reported so far for a β-CA, suggesting an important physiological role for this enzyme.
Our studies with different anionic CA inhibitors did not reveal specific inhibitors of DmBCA, which is in fact normal for such simple inorganic ions . The inhibition profile of DmBCA was unique for this enzyme, however, and differed significantly from other β-CAs studied previously.
Indeed, fluoride, chloride, cyanate, cyanide and sulfamide showed submillimolar inhibitory activity against DmBCA, with inhibition constants in the range of 150 μM - 970 μM. Another groups of anions, including bromide, iodide, thiocyanate, azide, hydrogen sulfide, bisulfite and sulfate showed inhibition constants close to 1 mM (KIs in the range of 1.01 - 1.36 mM), thus making them weak DmBCA inhibitors. Other anions such as bicarbonate, nitrate and nitrite, as well as the non-anionic species phenylboronic acid and phenylarsonic acid, were much less effective inhibitors, with inhibition constants of 22.4 - 43.7 mM. It is interesting to note that bicarbonate, a CA substrate, is a weak DmBCA inhibitor, but it appreciably inhibits the fungal enzymes Nce103 (C. Albicans), ScCA and CgCA, with KIs in the submillimolar range (of 86 μM - 0.78 mM). Even though the low apparent KI value for bicarbonate could be a reflection of allosteric inactivation of the enzyme ([39, 40]), the value is still in the same range with the other anion inhibitors, suggesting that their mechanism of inhibition might be quite similar. The present inhibition data may suggest a different physiological role for DmBCA in the insect compared to the corresponding β-CAs in fungal or yeast species studied previously. One should also note that DmBCA has a completely different inhibition profile not only compared to other β-CAs but also compared to the highly investigated α-CA of human origin, HCA II (Table 3).
Animal β-CAs constitute a separate subgroup in the β-CA family according to our bioinformatic analysis. Because the main target, β-CA, is present in many parasites and disease carrying vectors but absent in humans, this discovery clearly carries the significant potential for the design of parasite-specific β-CA inhibitors. Such drugs would presumably combine high specificity with a low incidence of side-effects in humans. These drugs could, for example, provide novel opportunities to restrict malaria propagation and to treat patients suffering from helminth infections.
Predictions of subcellular localization placed most, but not all, animal β-CAs in the mitochondria. Our hypothesis is that all of them will be mitochondrial based on three main reasons. First, we have direct experimental evidence of mitochondrial localization of DmBCA. Second, all of the TargetP predictions were either mitochondrial or cytoplasmic, both with low-reliability. Third, the neural network used in TargetP is trained with human and Arabidopsis mitochondrial proteins; thus, it is perhaps not perfectly adjusted to detect the signals in invertebrate mitochondrial transit peptides. Because of this mitochondrial localization of animal β-CAs, we speculate that the β-CA gene in chordates might have been silenced and ultimately lost after the emergence of mitochondrial α-CAs (CA V) to substitute for this function. The presence of an active CA in mitochondria has been linked to maintaining fundamental metabolic functions such as gluconeogenesis, ureagenesis and lipogenesis [1, 41]. Although our experiments were not focused on the role of the DmBCA enzyme, its mitochondrial localization and high enzymatic activity suggest that knockdown studies in D. melanogaster are warranted to further investigate the physiological function of β-CAs in animals.
The expression of DmBCA has been studied at mRNA level in two publicly available DNA microarray datasets. In FlyAtlas (http://flyatlas.org/atlas.cgi?name=CG11967-RA, ) expression levels in 17 adult and 8 larval D. melanogaster tissues are given. The highest upregulated values in adult are found in spermatheca (female), fat body, and heart. In larval tissues (third instar) downregulated or average expression levels are seen for all tissues. A time course study published in FlyBase (http://flybase.org/reports/FBgn0037646.html, under Microarray data, Personal communications to FlyBase, Gaurav et al. 2008), shows upregulated expression of DmBCA in early larval stages and late in metamorphosis, but downregulated levels in embryonal, late larval, early metamorphosis, and adult stages. The specific temporal and tissue patterns of expression suggest relevant functional and/or physiological roles for DmBCA.
In conclusion, fruit fly (Drosophila melanogaster) β-CA (DmBCA) is an active mitochondrial enzyme for the physiological reaction catalyzed by CAs, the hydration of CO2 to bicarbonate and protons. It is inhibited by various inorganic anions, boronic/arsonic acids or sulfonamides. Mammals do not possess β-class CAs, but these enzymes are widespread throughout the phylogenetic tree, making them exciting new targets for parasitic drug development. Indeed, β-CAs are found in many pathogenic organisms and pathogen vectors of the animal kingdom, including the malaria mosquito Anopheles, the yellow fever mosquito Aedes, the filariasis vector Culex, Ancylostoma hookworms, Brugia filarial nematodes, the blood fluke Schistosoma mansoni and the liver fluke Schistosoma sinensis. Because animal-derived β-CAs probably have a different predicted structure compared to the β-CAs investigated so far in detail from Archaea, plants, algae and fungi, as well as the α-CAs, specific inhibitors against these enzymes could be designed with minimal effect on human CAs and normal bacterial flora.
Taking advantage of the known pea β-CA (UniProt P17067) as an initial query, and subsequently the found invertebrate β-CAs (D. melanogaster, AAF54311; C. elegans CAJ43916), the animal CA sequences were retrieved from NCBI protein databases using Blast  http://blast.ncbi.nlm.nih.gov/Blast.cgi and from complete genomes at the UCSC Genome Bioinformatics Site http://genome.ucsc.edu using BLAT search algorithms . Hits were taken through iterated cycles of multiple sequence alignment (ClustalW; ), evaluation and revision. For revision, sequences with poorly matching or missing regions were subjected to gene model generation with GeneWise (; http://www.ebi.ac.uk/Tools/Wise2/), taking the genomic sequences from the UCSC site. EST and mRNA sequence data from NCBI were used to confirm gene models, sometimes to bridge gaps or fill ends in the genomic sequences, and to discover and assemble β-CAs from less than genome-wide sequenced organisms.
Phylogenetic trees were constructed from the multiple sequence alignments using MEGA 4 . Preliminary Neighbour-Joining trees  were constructed with all sequences, and a representative set of 26 sequences was then selected for a final tree, eliminating excessive insect and nematode sequences and severely incomplete sequences. The final tree was inferred using the Minimum evolution method  from a multiple sequence alignment containing columns corresponding to positions 31 to 255 (of 255) of DmBCA. The bootstrap consensus tree inferred from 1,000 replicates  was taken to represent the evolutionary history of the analyzed sequences. The tree was drawn to scale, with branch lengths proportional to the evolutionary distances used to infer the phylogenetic tree and rooted using the Cnidarian sequences as outgroup. The evolutionary distances were computed using the Poisson correction method and are presented in the units of the number of amino acid substitutions per site. The minimum evolution tree was searched using the Close-Neighbor-Interchange algorithm  at a search level of 3. The Neighbor-Joining algorithm  was used to generate the initial tree. All positions containing alignment gaps and missing data were eliminated only in pairwise sequence comparisons (Pairwise deletion option).
Predictions of subcellular localization were made using TargetP v1.1 (; http://www.cbs.dtu.dk/services/TargetP/), Secretome 2.0 http://www.cbs.dtu.dk/services/SecretomeP/, WolF PSORT, http://wolfpsort.org/ and MitoProt II v. 1.101 http://ihg2.helmholtz-muenchen.de/ihg/mitoprot.html. Conserved domain types were identified at the NCBI Conserved Domains Database http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml.
Construction of recombinant baculoviruses
Two constructs were engineered to study functional aspects of D. melanogaster β-CA (DmBCA). The first one contained a GFP (Green Fluorescent Protein) sequence fused C-terminally to the DmBCA cDNA for detection of the recombinant protein with confocal microscopy. The second construct contained a C-terminal histidine tag for protein purification. Both constructs contained full length β-CA gene, including the proposed N-terminal signal sequence. Bgl II and Xho I restriction sites and a thrombin cleavage site (for tag removal) were introduced into both constructs.
Total RNA extracted from D. melanogaster S2 cells (TRIzol® reagent and protocol, Invitrogen) was precipitated using sodium acetate at a final concentration of 100 mM and 100% ethanol. The solution was centrifuged at 13,000×g for 15 min at +4°C. The RNA sample was washed once with 70% ethanol and recentrifuged in the same manner. The sample was evaporated at room temperature and then suspended in sterile water. Precipitated RNA was transcribed into cDNA using a First strand cDNA Synthesis Kit #K1612 (Fermentas) according to the manufacturer's instructions.
The DmBCA gene was identified and amplified from cDNA by PCR using Phusion™ Hot Start High Fidelity DNA Polymerase (Finnzymes, Espoo, Finland). Sequence-specific primers were ordered from Sigma-Aldrich (St. Louis, MO). The forward primer was 5'-ATGGAGCGTATTTTGAGGGGAATC-3' (F1), and the reverse primer was 3'-CTACGAATAGAATCTTCTGACCTC-5' (R1). PCR was performed in a PTC 2000 thermal cycler (MJ Research, Waltham, MA), and the program consisted of a single 98°C denaturation step for 30 s, followed by 33 cycles of denaturation at 98°C for 10 s, annealing at 53°C for 15 s and extension at 72°C for 25 s. A final extension was performed at 72°C for 5 min. The PCR product band was separated from the gel and dissolved using Illustra™ GFX PCR DNA and GEL Band Purification Kit (GE Healthcare Life Sciences, Buckinghamshire, UK).
To engineer the DmBCA-GFP construct, the sequences of DmBCA and GFP were first amplified separately using sequence specific primers. The templates used for DmBCA and GFP were DmBCA gene product obtained from cDNA and pEGFP-N1™ Vector (Clontech), respectively. The constructs were then combined using PCR reactions whose R2 and F3 primer sequences partly overlapped (bolded in the primer sequences, see below). This allowed the ends of the PCR products to recognize each other. The forward primer used for DmBCA amplification was 5'- GGCCAGATCTATGGAGCGTATTTTGAGGGGA-3' (F2), and the reverse primer was 5'-CACGGAACCACGGGGCACCAG CGAATAGAATCTTCTGACCTC-3' (R2). The bolded sequence was designed to recognize the thrombin site and part of the GFP PCR product, while the other half recognized the DmBCA PCR product.
The forward primer used for GFP amplification was 5'-TCGCTGGTGCCCCGTGGTTCCGTG AGCAAGGGCGAGGAGCTG-3' (F3), and the reverse primer was 5'-CCGCTCGAGTTACTTGTACAGCTCGTCCAT-3' (R3). The bolded sequence in the F3 primer was designed to recognize the thrombin site and part of the DmBCA PCR product, while the other half recognized the GFP PCR product. The PCR program was as follows: 98°C for 2 min; then 33 cycles of 98°C for 10 s, 55°C for 15 s, and 72°C for 30 s; and finally 72°C for 7 min.
Both PCR products were run on agarose gel, purified and used as templates in the next PCR reaction. The PCR program was as follows: 98°C for 2 min; then 33 cycles of 98°C for 10 s, 55°C for 15 s and 72°C for 40 s; and finally 72°C for 8 min. The forward primer used in this PCR was (F2), and the reverse primer was (R3). The extra sequence overlaps at the 3' end of DmBCA and at the 5' end of the GFP allowed these PCR products to anneal to each other.
The PCR product was run on an agarose gel, and the obtained band was purified. pFastBac™ 1 plasmid (Invitrogen) and the PCR product were digested at +37°C overnight with BamHI and XhoI restriction enzymes (New England Biolabs). The digested plasmid and DmBCA-GFP construct were purified and then ligated overnight at +4°C using T4 DNA ligase (New England Biolabs). The ligated product was transformed into TOP10 bacteria (Invitrogen). Overnight cultures (8 ml) were made from these colonies, and plasmids were purified using a QIAprep Spin Miniprep Kit™ (Qiagen, Hilden, Germany). Sequencing was performed to verify the validity of the DmBCA-GFP construct. The construction of baculoviral genomes encoding the recombinant proteins has been described previously .
For recombinant protein production, the DmBCA gene construct with a C-terminal polyhistidine tag of six histidines was constructed and cloned into the pFastBac1™ vector. The forward primer used in the initial amplification of the DmBCA gene was (F2), and the reverse primer was 5'- GCCCTCGAGTTAATGGTGGTGATGGTGGTG GGAACCACGGGGCACCAGCGAATAGAATCTTCTGACCTC -3' (R4). The latter primer contains nucleotide repeats to create the polyhistidine tag (bolded in the primer sequence). The PCR program was as follows: 98°C for 60 s; then 35 cycles of 98°C for 10 s, 66°C for 15 s, and 72°C for 60 s; and finally 72°C for 5 min. Otherwise, the construct was made essentially in the same way as the DmBCA-GFP construct described above.
Study of the subcellular localization of DmBCA
A total of 100 μl of Sf9 insect cells (2 million cells/ml) were infected with 10 μl of baculovirus stock. The cells were kept at +27°C in incubator for three days in Lab-Tek™ Chamber Slide™ System™ plates (Nunc). The medium was then removed, and the cells were incubated in 600 μl of medium containing 100 nM Mitotracker Red CMXros™ (Invitrogen) for 20 minutes at +27°C. The cells were washed three times with 600 μl of medium and kept at +27°C in an incubator for two hours. The cells were then washed with PBS, fixed with 4% paraformaldehyde for five minutes and washed again with PBS. The cells were mounted in VectaShield® Mounting Medium (Vector Laboratories), covered with cover slips and analyzed using a confocal scanning laser microscope (Perkin Elmer-Cetus/Wallac UltraView LCI system™) with two different wavelengths: 488 nm for GFP detection and 579 nm for MitoTracker™. Image acquisition was performed with an Andor iXon™ DV885 EMCCD camera and the Andor iQ™ software (Andor).
Production and Purification of Recombinant Insect β-CA
The Sf9 insect cells were grown in HyQ SFX-Insect serum-free cell culture medium (HyClone, Logan, UT) in an orbital shaker at 27°C (125 rpm) for three days after infection. Although much of the DmBCA was associated with the cell pellet, protein purification was performed after centrifugation (5000 × g, 20°C, 8 min) from the supernatant and yielded highly pure DmBCA protein for characterization. Purification was performed using the Probond™ Purification System (Invitrogen) under native binding conditions, with wash and elution buffers made according to the manufacturer's instructions. The purification procedure per 500 ml of insect cell medium was as follows: 1 liter of native binding buffer and 25 ml of the nickel-chelating resin were added to the medium, and the His-tagged protein was then allowed to bind to the resin on a magnetic stirrer at 25°C for 3 h. The resin was washed with 100 + 30 ml of washing buffer (Invitrogen). The protein was then eluted with elution buffer (50 mM NaH2PO4, 500 mM NaCl, 250 mM imidazole, pH 8.0).
The purified DmBCA recombinant protein was transferred to a buffer of 50 mM Tris-HCl, pH 7.5, using an Amicon Ultracel™ - 10 k centrifugal filter device (Millipore) according to the manufacturer's instructions. To remove the His tag, the recombinant protein was treated with 60 μl of resin-coupled thrombin (Thrombin CleanCleave KIT™, Sigma) per 1 mg of protein with gentle shaking at 25°C for 1 h, according to the manufacturer's instructions. Protein concentration was determined using the DC Protein Assay™ (Bio-Rad) with three different dilutions. Purified recombinant DmBCA proteins were analyzed using 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) under reducing conditions. The gels were stained using the Colloidal Blue Staining Kit™ (Invitrogen).
CA activity measurements
An Applied Photophysics stopped-flow instrument was used to assay the CA-catalyzed CO2 hydration activity. Phenol red (at a concentration of 0.2 mM) was used as an indicator, working at the absorbance maximum of 557 nm, with 10 - 20 mM HEPES (pH 7.5) or Tris(pH 8.3) as buffers and 20 mM Na2SO4 or 20 mM NaClO4 (for maintaining constant ionic strength), following the initial rates of the CA-catalyzed CO2 hydration reaction for a period of 10 - 100 s. The CO2 concentrations ranged from 1.7 to 17 mM for the determination of kinetic parameters and inhibition constants. For each inhibitor at least six traces of the initial 5-10% of the reaction were used to determine the initial velocity. The uncatalyzed rates were determined in the same manner and subtracted from the total observed rates. Stock solutions of inhibitor (100 mM) were prepared in distilled-deionized water, and dilutions up to 0.01 μM were made thereafter with distilled-deionized water. Inhibitor and enzyme solutions were preincubated together for 15 min at room temperature prior to the assay to allow for the formation of the E-I complex. The inhibition constants were obtained by non-linear least-squares methods using PRISM 3, whereas the kinetic parameters for the uninhibited enzymes were obtained from Lineweaver-Burk plots, each representing the mean of at least three different determinations.
Kinetic measurements have been performed also with m-cresol purple (as indicator) - bicine (as buffer) (data not shown), and the results were the same (±5-10% of the reported values, which is the error range of this method) both for the kinetic parameters of CO2 hydration and for the inhibition constants of anionic inhibitors investigated here. Thus, the standard method reported in this paper is reliable for the investigation of β-CAs (in addition to the α-CAs) as reported by this group for several enzymes, such as the three β-class enzymes from Mycobacterium tuberculosis, Helicobacter pylori and Brucella suis, and the fungal class enzymes (for example ).
Dynamic light scattering analysis
The hydrodynamic diameters of the proteins were determined by dynamic light scattering (DLS) using Zetasizer ZS (Malvern Instruments Ltd., Worcestershire, United Kingdom). A 100 μl sample of DmBCA (200 μg/ml) in elution buffer (50 mM NaH2PO4, 500 mM NaCl, 250 mM imidazole, pH 8.0) was analyzed. Human CA II (HCA II) (180 μg/ml) in 0.1 M Tris, 0.4 M NaN3, 1 mM benzamidine, 20% glycerol, pH 7.0 was analyzed to support the results of DmBCA analysis. The small molecules such as imidazole and glycerol of the elution buffers seemed to dominate in the DLS analysis, since most of the light scattering was from particles having diameter smaller than one nm. Therefore, proteins were exchanged to 50 mM Na2HPO4 pH 7.0 containing 100 mM NaCl using protein desalting spin columns (Pierce). DLS analysis was then performed for 100 μl sample by using temperature scanning mode where the temperature was raised from 4°C to 50°C at 3°C intervals. Sample was let to equilibrate to each measurement temperature two minutes before data acquisition. For DmBCA three parallel measurements were carried out at each temperature and for HCA II only one measurement was performed at each temperature. The molecular weight of the protein was estimated from hydrodynamic diameter using globular protein standard curve provided by the manufacturer.
Analytical gel filtration
The molecular size of the protein in solution was measured by size-exclusion chromatography using Superdex200 10/300GL column (GE Healthcare) connected to ÄKTA ™purifier-100 equipped with UV-900 monitor (GE Healthcare). The analysis was done using 50 mM Na2HPO4, 650 mM NaCl (pH 7.0) as mobile phase. 20-30 μg protein was injected per run. All the analyses were done with flow rate 0.3 ml/min, and absorbencies at 280 nm and 205 nm were used to locate the protein peaks in the chromatograms. Molecular weight calibration curve was prepared by analyzing gel filtration standard protein mixture containing thyroglobulin (670kDa), γ-globulin (158kDa), ovalbumin (44kDa), myoglobin (17kDa) and vitamin B12 (1,35kDa) (Bio-Rad).
We thank Professor Hans Spelbrink and Dr. Peter Martinsson for help with fluorescence confocal microscopy and Aulikki Lehmus and Jukka Lehtonen for skilful technical assistance. We thank Joanna Zmurko for assistance in gel filtration analysis. This research was financed in part by grants from the 6th Framework Programme of the European Union (DeZnIT project) and the Competitive Research Funding of the Tampere University Hospital (9L071).
- Sly WS, Hu PY: Human carbonic anhydrases and carbonic anhydrase deficiencies. Annu Rev Biochem. 1995, 64: 375-401. 10.1146/annurev.bi.64.070195.002111.PubMedView ArticleGoogle Scholar
- Lane TW, Saito MA, George GN, Pickering IJ, Prince RC, Morel FM: Biochemistry: a cadmium enzyme from a marine diatom. Nature. 2005, 435 (7038): 42-10.1038/435042a.PubMedView ArticleGoogle Scholar
- Xu Y, Feng L, Jeffrey PD, Shi Y, Morel FM: Structure and metal exchange in the cadmium carbonic anhydrase of marine diatoms. Nature. 2008, 452 (7183): 56-61. 10.1038/nature06636.PubMedView ArticleGoogle Scholar
- Tripp BC, Bell CB, Cruz F, Krebs C, Ferry JG: A role for iron in an ancient carbonic anhydrase. J Biol Chem. 2004, 279 (8): 6683-6687. 10.1074/jbc.M311648200.PubMedView ArticleGoogle Scholar
- Macauley SR, Zimmerman SA, Apolinario EE, Evilia C, Hou YM, Ferry JG, Sowers KR: The archetype gamma-class carbonic anhydrase (Cam) contains iron when synthesized in vivo. Biochemistry. 2009, 48 (5): 817-819. 10.1021/bi802246s.PubMedView ArticleGoogle Scholar
- Zimmerman SA, Ferry JG: The beta and gamma classes of carbonic anhydrase. Curr Pharm Des. 2008, 14 (7): 716-721. 10.2174/138161208783877929.PubMedView ArticleGoogle Scholar
- So AK, Espie GS, Williams EB, Shively JM, Heinhorst S, Cannon GC: A novel evolutionary lineage of carbonic anhydrase (epsilon class) is a component of the carboxysome shell. J Bacteriol. 2004, 186 (3): 623-630. 10.1128/JB.186.3.623-630.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Sawaya MR, Cannon GC, Heinhorst S, Tanaka S, Williams EB, Yeates TO, Kerfeld CA: The structure of beta-carbonic anhydrase from the carboxysomal shell reveals a distinct subclass with one active site for the price of two. J Biol Chem. 2006, 281 (11): 7546-7555. 10.1074/jbc.M510464200.PubMedView ArticleGoogle Scholar
- Strop P, Smith KS, Iverson TM, Ferry JG, Rees DC: Crystal structure of the "cab"-type beta class carbonic anhydrase from the archaeon Methanobacterium thermoautotrophicum. J Biol Chem. 2001, 276 (13): 10299-10305. 10.1074/jbc.M009182200.PubMedView ArticleGoogle Scholar
- Kimber MS, Pai EF: The active site architecture of Pisum sativum beta-carbonic anhydrase is a mirror image of that of alpha-carbonic anhydrases. EMBO J. 2000, 19 (7): 1407-1418. 10.1093/emboj/19.7.1407.PubMedPubMed CentralView ArticleGoogle Scholar
- Parisi G, Perales M, Fornasari MS, Colaneri A, Gonzalez-Schain N, Gomez-Casati D, Zimmermann S, Brennicke A, Araya A, Ferry JG, Echave J, Zabaleta E: Gamma carbonic anhydrases in plant mitochondria. Plant Mol Biol. 2004, 55 (2): 193-207. 10.1007/s11103-004-0149-7.PubMedView ArticleGoogle Scholar
- Elleuche S, Poggeler S: Evolution of carbonic anhydrases in fungi. Curr Genet. 2009, 55 (2): 211-222. 10.1007/s00294-009-0238-x.PubMedView ArticleGoogle Scholar
- Fabre N, Reiter IM, Becuwe-Linka N, Genty B, Rumeau D: Characterization and expression analysis of genes encoding alpha and beta carbonic anhydrases in Arabidopsis. Plant Cell Environ. 2007, 30 (5): 617-629. 10.1111/j.1365-3040.2007.01651.x.PubMedView ArticleGoogle Scholar
- Hewett-Emmett D: Evolution and distribution of the carbonic anhydrase gene families. EXS. 2000, 90 (90): 29-76.PubMedGoogle Scholar
- Klengel T, Liang WJ, Chaloupka J, Ruoff C, Schroppel K, Naglik JR, Eckert SE, Mogensen EG, Haynes K, Tuite MF, Levin LR, Buck J, Muhlschlegel FA: Fungal adenylyl cyclase integrates CO2 sensing with cAMP signaling and virulence. Curr Biol. 2005, 15 (22): 2021-2026. 10.1016/j.cub.2005.10.040.PubMedPubMed CentralView ArticleGoogle Scholar
- Bahn YS, Cox GM, Perfect JR, Heitman J: Carbonic anhydrase and CO2 sensing during Cryptococcus neoformans growth, differentiation, and virulence. Curr Biol. 2005, 15 (22): 2013-2020. 10.1016/j.cub.2005.09.047.PubMedView ArticleGoogle Scholar
- Innocenti A, Muhlschlegel FA, Hall RA, Steegborn C, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors: inhibition of the beta-class enzymes from the fungal pathogens Candida albicans and Cryptococcus neoformans with simple anions. Bioorg Med Chem Lett. 2008, 18 (18): 5066-5070. 10.1016/j.bmcl.2008.07.122.PubMedView ArticleGoogle Scholar
- Isik S, Kockar F, Arslan O, Guler OO, Innocenti A, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of the beta-class enzyme from the yeast Saccharomyces cerevisiae with anions. Bioorg Med Chem Lett. 2008, 18 (24): 6327-6331. 10.1016/j.bmcl.2008.10.100.PubMedView ArticleGoogle Scholar
- Zimmerman SA, Ferry JG, Supuran CT: Inhibition of the archaeal beta-class (Cab) and gamma-class (Cam) carbonic anhydrases. Curr Top Med Chem. 2007, 7 (9): 901-908. 10.2174/156802607780636753.PubMedView ArticleGoogle Scholar
- Isik S, Kockar F, Aydin M, Arslan O, Guler OO, Innocenti A, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors: inhibition of the beta-class enzyme from the yeast Saccharomyces cerevisiae with sulfonamides and sulfamates. Bioorg Med Chem. 2009, 17 (3): 1158-1163. 10.1016/j.bmc.2008.12.035.PubMedView ArticleGoogle Scholar
- Nishimori I, Minakuchi T, Kohsaki T, Onishi S, Takeuchi H, Vullo D, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors: the beta-carbonic anhydrase from Helicobacter pylori is a new target for sulfonamide and sulfamate inhibitors. Bioorg Med Chem Lett. 2007, 17 (13): 3585-3594. 10.1016/j.bmcl.2007.04.063.PubMedView ArticleGoogle Scholar
- Minakuchi T, Nishimori I, Vullo D, Scozzafava A, Supuran CT: Molecular cloning, characterization, and inhibition studies of the Rv1284 beta-carbonic anhydrase from Mycobacterium tuberculosis with sulfonamides and a sulfamate. J Med Chem. 2009, 52 (8): 2226-2232. 10.1021/jm9000488.PubMedView ArticleGoogle Scholar
- Nishimori I, Minakuchi T, Vullo D, Scozzafava A, Innocenti A, Supuran CT: Carbonic anhydrase inhibitors. Cloning, characterization, and inhibition studies of a new beta-carbonic anhydrase from Mycobacterium tuberculosis. J Med Chem. 2009, 52 (9): 3116-3120. 10.1021/jm9003126.PubMedView ArticleGoogle Scholar
- Cronk JD, Rowlett RS, Zhang KY, Tu C, Endrizzi JA, Lee J, Gareiss PC, Preiss JR: Identification of a novel noncatalytic bicarbonate binding site in eubacterial beta-carbonic anhydrase. Biochemistry. 2006, 45 (14): 4351-4361. 10.1021/bi052272q.PubMedView ArticleGoogle Scholar
- Joseph P, Turtaut F, Ouahrani-Bettache S, Montero JL, Nishimori I, Minakuchi T, Vullo D, Scozzafava A, Kohler S, Winum JY, Supuran CT: Cloning, characterization, and inhibition studies of a beta-carbonic anhydrase from Brucella suis. J Med Chem. 2010, 53 (5): 2277-2285. 10.1021/jm901855h.PubMedView ArticleGoogle Scholar
- Innocenti A, Hall RA, Schlicker C, Muhlschlegel FA, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of the beta-class enzymes from the fungal pathogens Candida albicans and Cryptococcus neoformans with aliphatic and aromatic carboxylates. Bioorg Med Chem. 2009, 17 (7): 2654-2657. 10.1016/j.bmc.2009.02.058.PubMedView ArticleGoogle Scholar
- Innocenti A, Hall RA, Schlicker C, Scozzafava A, Steegborn C, Muhlschlegel FA, Supuran CT: Carbonic anhydrase inhibitors. Inhibition and homology modeling studies of the fungal beta-carbonic anhydrase from Candida albicans with sulfonamides. Bioorg Med Chem. 2009, 17 (13): 4503-4509. 10.1016/j.bmc.2009.05.002.PubMedView ArticleGoogle Scholar
- Innocenti A, Winum JY, Hall RA, Muhlschlegel FA, Scozzafava A, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of the fungal beta-carbonic anhydrases from Candida albicans and Cryptococcus neoformans with boronic acids. Bioorg Med Chem Lett. 2009, 19 (10): 2642-2645. 10.1016/j.bmcl.2009.03.147.PubMedView ArticleGoogle Scholar
- Schlicker C, Hall RA, Vullo D, Middelhaufe S, Gertz M, Supuran CT, Muhlschlegel FA, Steegborn C: Structure and inhibition of the CO2-sensing carbonic anhydrase Can2 from the pathogenic fungus Cryptococcus neoformans. J Mol Biol. 2009, 385 (4): 1207-1220. 10.1016/j.jmb.2008.11.037.PubMedView ArticleGoogle Scholar
- Innocenti A, Leewattanapasuk W, Muhlschlegel FA, Mastrolorenzo A, Supuran CT: Carbonic anhydrase inhibitors. Inhibition of the beta-class enzyme from the pathogenic yeast Candida glabrata with anions. Bioorg Med Chem Lett. 2009, 19 (16): 4802-4805. 10.1016/j.bmcl.2009.06.048.PubMedView ArticleGoogle Scholar
- Smith KS, Cosper NJ, Stalhandske C, Scott RA, Ferry JG: Structural and kinetic characterization of an archaeal beta-class carbonic anhydrase. J Bacteriol. 2000, 182 (23): 6605-6613. 10.1128/JB.182.23.6605-6613.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Cox EH, McLendon GL, Morel FM, Lane TW, Prince RC, Pickering IJ, George GN: The active site structure of Thalassiosira weissflogii carbonic anhydrase 1. Biochemistry. 2000, 39 (40): 12128-12130. 10.1021/bi001416s.PubMedView ArticleGoogle Scholar
- Smith KS, Ingram-Smith C, Ferry JG: Roles of the conserved aspartate and arginine in the catalytic mechanism of an archaeal beta-class carbonic anhydrase. J Bacteriol. 2002, 184 (15): 4240-4245. 10.1128/JB.184.15.4240-4245.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Tetu SG, Tanz SK, Vella N, Burnell JN, Ludwig M: The Flaveria bidentis beta-carbonic anhydrase gene family encodes cytosolic and chloroplastic isoforms demonstrating distinct organ-specific expression patterns. Plant Physiol. 2007, 144 (3): 1316-1327. 10.1104/pp.107.098152.PubMedPubMed CentralView ArticleGoogle Scholar
- Mitra M, Lato SM, Ynalvez RA, Xiao Y, Moroney JV: Identification of a new chloroplast carbonic anhydrase in Chlamydomonas reinhardtii. Plant Physiol. 2004, 135 (1): 173-182. 10.1104/pp.103.037283.PubMedPubMed CentralView ArticleGoogle Scholar
- Bahn YS, Muhlschlegel FA: CO2 sensing in fungi and beyond. Curr Opin Microbiol. 2006, 9 (6): 572-578. 10.1016/j.mib.2006.09.003.PubMedView ArticleGoogle Scholar
- Fasseas MK, Tsikou D, Flemetakis E, Katinakis P: Molecular and biochemical analysis of the beta class carbonic anhydrases in Caenorhabditis elegans. Mol Biol Rep. 2010, 37 (6): 2941-50. 10.1007/s11033-009-9857-z.PubMedView ArticleGoogle Scholar
- Supuran CT: Carbonic anhydrases: novel therapeutic applications for inhibitors and activators. Nat Rev Drug Discov. 2008, 7 (2): 168-181. 10.1038/nrd2467.PubMedView ArticleGoogle Scholar
- Rowlett RS, Tu C, Lee J, Herman AG, Chapnick DA, Shah SH, Gareiss PC: Allosteric site variants of Haemophilus influenzae beta-carbonic anhydrase. Biochemistry. 2009, 48 (26): 6146-6156. 10.1021/bi900663h.PubMedPubMed CentralView ArticleGoogle Scholar
- Rowlett RS, Hoffmann KM, Failing H, Mysliwiec MM, Samardzic D: Evidence for a bicarbonate "escort" site in Haemophilus influenzae beta-carbonic anhydrase. Biochemistry. 2010, 49 (17): 3640-3647. 10.1021/bi100328j.PubMedPubMed CentralView ArticleGoogle Scholar
- Hazen SA, Waheed A, Sly WS, LaNoue KF, Lynch CJ: Differentiation-dependent expression of CA V and the role of carbonic anhydrase isozymes in pyruvate carboxylation in adipocytes. FASEB J. 1996, 10 (4): 481-490.PubMedGoogle Scholar
- Chintapalli VR, Wang J, Dow JA: Using FlyAtlas to identify better Drosophila melanogaster models of human disease. Nat Genet. 2007, 39 (6): 715-720. 10.1038/ng2049.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664.PubMedPubMed CentralView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Rzhetsky A, Nei M: METREE: a program package for inferring and testing minimum-evolution trees. Comput Appl Biosci. 1994, 10 (4): 409-412.PubMedGoogle Scholar
- Zharkikh A, Li WH: Estimation of confidence in phylogeny: the complete-and-partial bootstrap technique. Mol Phylogenet Evol. 1995, 4 (1): 44-63. 10.1006/mpev.1995.1005.PubMedView ArticleGoogle Scholar
- Nei M: Phylogenetic analysis in molecular evolutionary genetics. Annu Rev Genet. 1996, 30: 371-403. 10.1146/annurev.genet.30.1.371.PubMedView ArticleGoogle Scholar
- Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol. 2000, 300 (4): 1005-1016. 10.1006/jmbi.2000.3903.PubMedView ArticleGoogle Scholar
- Hilvo M, Baranauskiene L, Salzano AM, Scaloni A, Matulis D, Innocenti A, Scozzafava A, Monti SM, Di Fiore A, De Simone G, Lindfors M, Janis J, Valjakka J, Pastorekova S, Pastorek J, Kulomaa MS, Nordlund HR, Supuran CT, Parkkila S: Biochemical characterization of CA IX, one of the most active carbonic anhydrase isozymes. J Biol Chem. 2008, 283 (41): 27799-27809. 10.1074/jbc.M800938200.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.