Document Detail

In silico identification and characterization of the ion transport specificity for P-type ATPases in the Mycobacterium tuberculosis complex.
Jump to Full Text
MedLine Citation:
PMID:  23031689     Owner:  NLM     Status:  MEDLINE    
BACKGROUND: P-type ATPases hydrolyze ATP and release energy that is used in the transport of ions against electrochemical gradients across plasma membranes, making these proteins essential for cell viability. Currently, the distribution and function of these ion transporters in mycobacteria are poorly understood.
RESULTS: In this study, probabilistic profiles were constructed based on hidden Markov models to identify and classify P-type ATPases in the Mycobacterium tuberculosis complex (MTBC) according to the type of ion transported across the plasma membrane. Topology, hydrophobicity profiles and conserved motifs were analyzed to correlate amino acid sequences of P-type ATPases and ion transport specificity. Twelve candidate P-type ATPases annotated in the M. tuberculosis H37Rv proteome were identified in all members of the MTBC, and probabilistic profiles classified them into one of the following three groups: heavy metal cation transporters, alkaline and alkaline earth metal cation transporters, and the beta subunit of a prokaryotic potassium pump. Interestingly, counterparts of the non-catalytic beta subunits of Hydrogen/Potassium and Sodium/Potassium P-type ATPases were not found.
CONCLUSIONS: The high content of heavy metal transporters found in the MTBC suggests that they could play an important role in the ability of M. tuberculosis to survive inside macrophages, where tubercle bacilli face high levels of toxic metals. Finally, the results obtained in this work provide a starting point for experimental studies that may elucidate the ion specificity of the MTBC P-type ATPases and their role in mycobacterial infections.
Lorena Novoa-Aponte; Andrés León-Torres; Miyer Patiño-Ruiz; Jenifer Cuesta-Bernal; Luz-Mary Salazar; David Landsman; Leonardo Mariño-Ramírez; Carlos-Yesid Soto
Related Documents :
24649699 - Gas separation using porous cement membrane.
24731639 - Rheological and structural properties of differently acidified and renneted milk gels.
9591359 - A long-lifetime ru(ii) metal-ligand complex as a membrane probe.
6631959 - Structure and thermotropic properties of 1,3-dipalmitoyl-glycero-2-phosphocholine.
6824339 - A kinetic and spectral study of the alkaline transitions of chloroperoxidase.
17253739 - Calorimetric analysis of thermal phase transitions in functionalized microgels.
16216879 - A unified view of the role of electrostatic interactions in modulating the gating of cy...
1432979 - Uterine surface temperature changes caused by electrosurgical endometrial coagulation.
25275859 - Perfluoroalkane force field for lipid membrane environments.
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Extramural; Research Support, Non-U.S. Gov't     Date:  2012-10-03
Journal Detail:
Title:  BMC structural biology     Volume:  12     ISSN:  1472-6807     ISO Abbreviation:  BMC Struct. Biol.     Publication Date:  2012  
Date Detail:
Created Date:  2013-02-15     Completed Date:  2013-07-01     Revised Date:  2013-07-11    
Medline Journal Info:
Nlm Unique ID:  101088689     Medline TA:  BMC Struct Biol     Country:  England    
Other Details:
Languages:  eng     Pagination:  25     Citation Subset:  IM    
Chemistry Department, Faculty of Sciences, Universidad Nacional de Colombia, Bogotá, Colombia, Carrera 30 # 45–03, Ciudad Universitaria, Bogotá, Colombia.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Adenosine Triphosphatases / chemistry,  metabolism*
Amino Acid Sequence
Computer Simulation
Conserved Sequence
Hydrophobic and Hydrophilic Interactions
Ion Transport
Ions / metabolism
Markov Chains
Models, Molecular
Molecular Sequence Data
Mycobacterium tuberculosis / chemistry,  enzymology*,  physiology
Sequence Alignment
Substrate Specificity
Tuberculosis / microbiology*
Grant Support
Reg. No./Substance:
0/Ions; EC 3.6.1.-/Adenosine Triphosphatases

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): BMC Struct Biol
Journal ID (iso-abbrev): BMC Struct. Biol
ISSN: 1472-6807
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Novoa-Aponte et al.; licensee BioMed Central Ltd.
Received Day: 25 Month: 6 Year: 2012
Accepted Day: 27 Month: 9 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 3 Month: 10 Year: 2012
Volume: 12First Page: 25 Last Page: 25
PubMed Id: 23031689
ID: 3573892
Publisher Id: 1472-6807-12-25
DOI: 10.1186/1472-6807-12-25

In silico identification and characterization of the ion transport specificity for P-type ATPases in the Mycobacterium tuberculosis complex
Lorena Novoa-Aponte1 Email:
Andrés León-Torres1 Email:
Miyer Patiño-Ruiz1 Email:
Jenifer Cuesta-Bernal1 Email:
Luz-Mary Salazar1 Email:
David Landsman2 Email:
Leonardo Mariño-Ramírez23 Email:
Carlos-Yesid Soto1 Email:
1Chemistry Department, Faculty of Sciences, Universidad Nacional de Colombia, Bogotá, Colombia, Carrera 30 # 45–03, Ciudad Universitaria, Bogotá, Colombia
2Computational Biology Branch, NCBI, NLM, NIH, Bethesda, USA
3PanAmerican Bioinformatics Institute, Santa Marta, Magdalena, Colombia


Tuberculosis (TB) is one of the most important challenges in public health maintenance throughout the world. According to the World Health Organization (WHO), 8.5-9.2 million new TB cases were estimated to have occurred in 2010 [1], and 1.2–1.5 million deaths were caused by species of the Mycobacterium tuberculosis complex (MTBC) that includes M. tuberculosis, M. bovis, M. bovis BCG (vaccine strain), M. africanum, M. microti, M. canettii, and M. pinnipedii, which produces TB in humans and some animal hosts [2,3]. Part of the infected population will develop active TB, whereas the majority of cases (approximately 90%) progress to a non-infectious disease or latent TB, where mycobacteria survive in a dormant state inside immune cells [4]. Individuals with latent TB may be asymptomatic during prolonged periods of time; however, TB can reactivate when the host immune response diminishes due to malnutrition, steroid use, and HIV co-infection [5].

The emergence of multidrug and extensively drug-resistant tuberculosis strains (MDR-TB and XDR-TB) and the lack of drugs against latent TB have become serious problems for TB control. Therefore, the identification of new therapeutic targets useful in the development of novel drugs and vaccines against latent TB is essential. New anti-TB drugs, such as diarylquinolines (TMC207) and benzothiazines (BTZ043) target essential membrane proteins that affect mycobacterial viability [6]. Thus, antimicrobials designed against proteins of plasma membrane are ideal because they avoid problems related to membrane permeability.

Ion transport in bacteria is carried out by enzymatic systems that belong to either the P-type ATPase, ATP binding cassettes (ABC transporters) and metallic ion/H+-antiporter systems [7]. In general, ATPases help maintain the ion gradients responsible for cell volume control and transport of nutrients across the cell membrane [8-11]. ATPases hydrolyze ATP, releasing energy that is used in the transport of ions against electrochemical gradients in plasma membranes. The enzymatic mechanisms of P-type ATPases were initially described in eukaryotic cells [12-14]. These enzymes have five different functional and structural domains: three of these domains are cytoplasmic (A, actuator; N, nucleotide binding and P, phosphorylation), and the other two are embedded in the membrane (T, transport and S, class specific support domain) [15]. P-type ATPases have the following two conformational states: E1, which binds ion substrates on one side of the cell membrane, inducing their auto-phosphorylation and generating a new conformational state; and E2, which has a lower affinity for substrates and therefore releases them to the other side of the cell membrane, promoting ion transport and finally recovery of the E1 state [15,16].

Despite P-type ATPases share the same catalytic mechanism based on conformational changes in their five structural domains, their regulation and substrate affinities are different [15,16]. P-type ATPases are phylogenetically classified into five subfamilies (PI-PV), and within these subfamilies are 10 different subtypes that are categorized based on the transported substrate [17].

In this study, probabilistic profiles were constructed to compare and classify all P-type ATPases of the MTBC based on their structural features and ion transport. Twelve possible P-type ATPases were detected in the proteome of M. tuberculosis H37Rv and from other members of the MTBC. The high number of heavy metal transporters discovered in the MTBC suggests an important role for P-type ATPases in M. tuberculosis survival within macrophages.

Construction of hidden Markov models (HMM)

To obtain a representative group of each phylogenetic subfamily, a set of 128 well-characterized P-type ATPases with evidence of existence at the protein level were retrieved from Uniprot (Swiss-Prot section) [18]. Each group of sequences was aligned using the Praline tool ( with the BLOSUM62 matrix and Phobius transmembrane structure predictor, which was developed to improve the multiple alignments of membrane protein sequences [19]. The HMM package of programs [20], available in the Mobyle Pasteur portal (, was used to find these types of pumps in the MTBC proteomes. The Hmmbuild tool with default settings and the multiple sequence alignments was used for HMM building. The default parameters of the Hmmer tools use an ad hoc position-based sequence weighting algorithm that makes the models appropriate for the identification of distant members of the P-type ATPases family. Consensus sequences were generated with the Hmmemit tool.

Search and classification of MTBC P-type ATPases

To date, the following 10 MTBC genomes have been completely sequenced and assembled (NCBI): M. africanum GM041182 (NC_015758), M. bovis AF2122/97 (NC_002945), M. bovis BCG str. Pasteur 1173P2 (NC_008769), M. bovis BCG str. Tokyo 172 (NC_012207), M. canettii CIPT 140010059 (NC_015848) and five M. tuberculosis strains, H37Rv (NC_000962), H37Ra (NC_009525), F11 (NC_009565), CDC1551 (NC_002755) and KZN1435 (NC_012943); these proteomes were obtained from Uniprot ( Because strains of M. microti and M. pinnipedii had not been sequenced, they were not included in this study. HMM and Hmmsearch tool were used to find P-type ATPases in the MTBC proteomes.

Topology prediction

The topology derived from predictions of transmembrane segments (TMS) was made with the following six programs: TopPred (, DAS (, TMpred (, TMHMM 2.0 (, HMMTOP ( and Phobius ( All programs were used with default settings, except TopPred, which allows the user to provide information about the type of organism. TMDET and PPM servers were used in cases where the predictions based on amino acid sequences did not yield reliable results. These tools must be fed with PDB files to identify TMS based on tertiary protein structure and to generate 3D models with TMS located into hypothetical lipid bilayer planes. Tertiary structure models of CtpH and CtpI were made with the I-TASSER tool [21] based on the threading strategy and were validated with the Whatif package of programs (

Hydrophobicity profile construction

The amino acid sequence of the M. tuberculosis H37Rv P-type ATPases and consensus sequences generated by Hmmemit tool were analyzed by TMHMM 2.0. Hydrophobicity profiles of consensus sequences for previously characterized P-type ATPases were used as comparison patterns.

Conserved motif analysis

Amino acids sequences of the identified proteins were manually analyzed to determine the nine conserved motifs typical of the P-type ATPase family, as previously described by Thever et al. [22] and others [16,23].

Results and discussion
An unusually high number of cation transporter P-type ATPases are present in the Mycobacterium tuberculosis complex

Multiple alignments of 128 reported P-type ATPase protein sequences from a representative group of eukaryotic and prokaryotic cells (obtained from a curated database and confirmed at the protein level) allowed the identification of highly conserved regions within the family and the classification of these members according to ion transport. These alignments were used as the starting point for the construction of HMM profiles that represented groups of ion transporters, which were also used to generate a consensus sequence for each group. The designed HMM were then used to identify P-type ATPases in proteomes of the different MTBC species, whose genome sequences have been reported. The proposed classification for the studied sequences was based on the obtained scores using the HMM and the Hmmsearch tool.

Because the P-type ATPases are transmembrane proteins, typical alignments using BLOSUM and PAM matrices were not adequate. In this study, the Praline tool was used because it considers the differences in evolutionary tendencies of transmembranal and non-transmembranal regions. Praline applies an adequate matrix for each protein region based on the previous prediction of TMS (made with the Phobius algorithm in this case). As P-type ATPases contain ion-binding motifs within a transmembrane domain, the strategy that we employed produced more reliable alignments. Twelve hypothetical proteins with a high probability of being P-type ATPase transporters of metallic cations were identified in each of the MTBC proteomes. All proteins identified in the MTBC share at least 98% identity with their orthologs in the scanned proteomes; however, M. canettii and M. africanum displayed slight differences in non-conserved regions. The pumps identified in the M. tuberculosis H37Rv proteome, CtpA, CtpB, CtpC, CtpD, CtpE, CtpF, CtpG, CtpH, CtpI, CtpJ, CtpV and KdpB, were in agreement with the automatic annotation of probable cation transporter P-type ATPases in the M. tuberculosis H37Rv genome [24] and were taken as references to facilitate further analysis of results.

A broad diversity of cations are potentially transported by Mycobacterium tuberculosis P-type ATPases

Figure 1 shows the Hmmsearch scores obtained when the 16 HMM were used to find P-type ATPases in the M. tuberculosis H37Rv proteome. Those scores show the similarity between the identified sequences in the H37Rv proteome and the grouped sequences used for each HMM construction. The patterns thus obtained allowed classification of the 12 P-type ATPases into the following three groups: transporters of heavy metal (HM) cations, transporters of alkaline and alkaline earth metal (AEM) cations, and a group composed only of the β subunit of the prokaryotic K+ transporter, KdpB. The highest similarity with the conventional P-type ATPases used in this study corresponded to KdpB, CtpF and CtpV, whereas CtpE, CtpH and CtpI had the lowest scores.

The HM group is composed of CtpA, CtpB, CtpC, CtpD, CtpG, CtpJ and CtpV, which may transport Cu2+, Cu+, Co2+, Ag+, Hg2+, Cd2+, Pb2+ and Zn2+. Most of the proteins analyzed in this work correspond to the HM group (60%) suggesting that the active transport of heavy metal cations is relevant for the tubercle bacilli persistence, as it has been hypothesized for other prokaryotes and some unicellular eukaryotes [23]. Interestingly, evidence of toxic concentrations of intracellular heavy metals has been recently described in macrophages infected with M. tuberculosis[25].

Alternatively, CtpE, CtpF, CtpH and CtpI belong to the AEM group and may be involved in Na+, K+, Ca2+, H+ and Mg2+ transport. It is possible that some of these proteins correspond to Na+/K+ or H+/K+ ATPase pumps, but they only appeared to contain the α subunits in their tridimensional structure. Their non-catalytic-β-subunit counterpart, which has been correlated with regulatory processes and assembly of P-type ATPases into the cell membrane of eukaryotic cells [10], was not found in any of the MTBC proteomes.

KdpB is very different protein from the other MTBC ATPases. It shares 63% identity with the KdpB of E. coli (P03960), which corresponds to the β subunit of a K+ transporting multimeric ATPase [22,26]. Despite the observation that KdpB has the characteristic DKTGTLT phosphorylation motif of this type of pump, it does not have a typical ion binding motif [26]. In conclusion, the proposed classification for P-type ATPases of the MTBC provides an initial approximation of their functional characterization. It is noteworthy that the constructed HMM in this work became a useful tool for the identification of P-type ATPases in other biological systems.

Three different topologies can be adopted by P-type ATPases in Mycobacterium tuberculosis

The six different algorithms used in the hydrophobicity analysis showed that all of the M. tuberculosis H37Rv P-type ATPases have an α-helix type TMS containing at least 17 amino acid residues, and the agreement of the topology results obtained using these prediction tools increases the confidence in the predicted TMS. Consensus transmembrane regions for each P-type ATPase were obtained if at least four of the algorithms gave similar results. Because three of the tools (TMHMM 2.0, HMMTOP and Phobius) incorporate size and composition restrictions in TMS, the results obtained using these algorithms are significant.

The strategy used for topology analysis showed the following three different types of topology within the M. tuberculosis H37Rv P-type ATPases: type I, which corresponds to HM pumps with eight TMS (A, B and from 1 to 6), type II, which corresponds to AEM pumps with 10 TMS (from 1 to 10), and type III, which corresponds to KdpB with 7 TMS (Figure 2). It was observed that all P-type ATPases have two cytoplasmic loops (small and large) that include the phosphorylation and ATP binding sites. The small cytoplasmic loop is located between TMS2 and TMS3, whereas the largest cytoplasmic loop is situated between TMS4 and TMS5 (Figure 2) [22]. Additionally, as was expected, the N- and C-termini of these proteins are located in the cytosolic side, except in the case of KdpB, in which the N- terminus is located intracellularly and the C- terminus is outside the cytoplasm.

Disagreement in the results was observed in CtpH and CtpI analysis using the different prediction tools. CtpI showed some hydrophobic amino acid sequences that could be considered part of TMS, but they do not fulfill all of the characteristics associated with TMS; meanwhile, CtpH did not show the expected hydrophobicity pattern for P-type ATPases (Figure 3). Therefore, at first it was difficult to determine whether CtpH and CtpI were in fact typical of P-type ATPases. Recent studies have classified CtpH and CtpI proteins of M. bovis as FUPA 24 (TC No. 3.A.3.24), i.e., “Functionally uncharacterized P-type ATPase family 24”, and classify them as transporter proteins in the TCDB-Transporter Classification Database ( FUPA 24 proteins are homologous to P-type ATPases but have an unusually large N-terminal segment that makes them twofold longer than typical P-type ATPases [23] with two TMS. In addition, functional motifs of P-type ATPases are located within the C-terminal region of FUPA 24 [23].

By contrast, topology analysis in this study showed that M. tuberculosis H37Rv CtpH and CtpI contain more than the expected TMS (from 3 to 12) for FUPA 24 proteins. To overcome the discrepancy between TCDB and the six topology tools used in this work, the TMDET and PPM servers were used. These tools allow for the identification of TMS based on the tertiary structure of the protein; thus, not only can the correct size and composition of α-helixes be guaranteed, but also their adequate disposition and organization in the plasma membrane can be determined. Ten TMS were found in CtpH and CtpI using both the TMDET and PPM servers. These segments are located with high probability in a hypothetical lipid bilayer as shown in Figure 4. This result strongly suggests that CtpH and CtpI have a type II topology, as has been determined for other members of the AEM group.

Hydrophobicity profiles suggest insights into ion transport specificity

The “topological clustering” based on the hydrophobicity profiles obtained with the TMHMM 2.0 algorithm can be used as a predictive tool to determine substrate specificity for HM pumps [27]. Here, this strategy was applied to analyze HM, AEM and KdpB groups of M. tuberculosis H37Rv P-type ATPases. The hydrophobicity profiles from each M. tuberculosis P-type ATPase was compared with the obtained profiles from previously characterized P-type ATPases. As shown in Figure 3, the hydrophobicity profiles of CtpA, CtpB and CtpV were similar to the consensus Cu+ P-type ATPases; the CtpC profile was similar to the different Zn2+ transporters Zn2+, Cd2+/Zn2+ and Pb2+/Cd2+/Zn2+/Hg2+. In addition, CtpD, CtpG and CtpJ had regions similar to the Cd2+/Zn2+/Co2+ transporter. In the case of the AEM group, CtpF possesses TMS that were similar to the hydrophobic profiles obtained for consensus of Na+, Mg2+ and Na+/K+ ATPases. Additionally, CtpE was similar to Na+ or H+ transporters; however, its hydrophobicity profile is closer to that of Na+ P-type ATPase. Considering the C-terminus region of CtpH (819 amino acid residues), this pump had a hydrophobicity profile very similar to the Ca2+ PMCA consensus. Finally, KdpB had the exact same hydrophobic profile as that which was determined for E. coli KdpB.

In general, there are the following two types of hydrophobicity profiles: the first comprises the AEM and the second that corresponds to HM. Interestingly, the relative position of TMS along amino acid sequences is highly conserved between both types of hydrophobicity profiles. The profile for AEM is characterized by a single TMS in the C-terminal region and two in the N-terminal region, whereas profiles for HM P-type ATPases have three highly hydrophobic regions, one of them in the N-terminus, one in the middle of the amino acid sequence, and the third in the C-terminal end. Although the hydrophobicity profile of consensus sequences is similar among enzymes in each group, there are significant differences between profiles that allow their differentiation according to the transported ion.

P-type ATPases of the Mycobacterium tuberculosis complex possess the characteristic motifs involved in cation transport

The common catalytic mechanism of P-type ATPases is partially supported by conserved core sequences that are useful in the recognition of these pumps in proteomes [15]. To find the conserved motifs for P-type ATPases along the protein sequences, manual analysis was carried out because the servers commonly used for analysis of conserved motifs do not identify subtle variations within particular motifs of P-type ATPases. Nine conserved motifs typical for P-type ATPases were found in the MTBC pumps (Figure 5). We observed that each motif sequence was identical within orthologs of the MTBC. From N-terminus to the C-terminus, the motifs were as follows: motifs 1 [(PVA)G(DE)] and 2 [P(AS)D], related to conformational changes, together with motif 3 [TGE(SA)], associated with phosphatase activity, were located between TMS 2 and 3 in the actuator domain; the sequence of motif 4, which determines ion specificity, can be [PEG(LM)] or [(CSA)PCA(LV)], which was located at the end of the fourth TMS; the phosphorylation site [DKTGTLT] was found in the P domain; the remaining motifs, motif 6 [(KI)GA(PVA)(EDA)], an ATP binding facilitator, motifs 7 [(DV)(ASPI)(VP)(KAR)] and 8 [(MLV)I(TS)GD], involved in phosphorylation catalysts, and motif 9 [(VTC)AM(TV)GDG(VSAT)ND(AV)(PAL)A(LI)(RKA)(QMAD)A(DNT)(VI)G(VI)(AG)(MV)], the hinge motif, which provides the flexibility necessary to achieve conformational changes during the pumping process, could be found between the fourth and fifth TMS [22]. Motifs 8 and 9 contain amino acid residues [TGDN and GDGXND] responsible for the coordination of Mg2+ as a cofactor of the enzymes [15,16,22,23].

It was observed that characteristic motifs of P-type ATPases exhibit slight variations compared with previously reported sequences for eukaryotic organisms [22]. The most conserved motif (motif 5) was almost the same within the 12 P-type ATPases, except CtpH and KdpB, which contained the conservative substitutions [DKTGTL(TS)] and [DKTGT(LI)T], respectively. The sequence of motif 4 [PEGL] has been previously associated with binding of AEM cations [22], and it has been found in CtpE, CtpF, and CtpI. Variation in the fourth position of this motif [PEGM] in CtpH has not been reported previously. CtpA, CtpB, CtpC, CtpD, CtpJ, CtpG and CtpV have the motif 4 sequence [(CSA)PCA(LV)] characteristic of HM transporters. Alternatively, KdpB, which corresponds to a β subunit of a multimeric P-type ATPase, did not have motif 4. This finding can be explained by the observation that this subunit mediates phosphorylation/dephosphorylation and energy transduction during K+ transport. In E. coli, the coordination function lies mainly in the KdpA subunit [26].

Cu+, Zn2+ and Co2+ pumping is mediated by P-type ATPases in the M. tuberculosis complex

The PIB phylogenetic group of P-type ATPases is composed of HM pumps, and according to ion specificity, the PIB group is subdivided into five subgroups (1 to 5). Alignments of consensus sequences from characterized PIB ATPases with the HM P-type ATPases of MTBC were used for locating additional transmembrane motifs of this type of transporter [28-30], as shown in Figure 6. CtpA, CtpB and CtpV have the characteristic motifs of PIB-1 group members (Table 1), corroborating their ion transport specificity as predicted by hydrophobicity profiles, and suggesting that Cu+ is transported preferentially to Cu2+ by these pumps.

This result is expected, because the special reducing conditions inside mycobacteria [27] are similar to the intra-phagosomal environment in macrophages where the bacilli reside. Recently, it was reported that the M. tuberculosis H37Rv CtpV is a Cu+ exporter P-type ATPase, reinforcing our predictions [31]. CtpD and CtpJ exhibit the [SCP] and [HEGT] motifs of the PIB-4 group in the TMS6 and TMS8, respectively. Some reports have indicated that amino acids involved in Co2+ transport are still not well understood; for this reason, the absence of residue N in TMS7 cannot rule out the possibility that CtpD and CtpJ are Co2+ transporters. CtpC only exhibits the [CPC(X)4S] motif in TMS6, which allows its classification as PIB-2. In fact, it was recently reported as a Zn2+ P-type ATPase in M. tuberculosis H37Rv [25]; deficiency in this pump produced zinc accumulation within the mycobacterial cytoplasm, resulting in impaired intracellular growth of tubercle bacilli [25]. Finally, CtpG does not have any of the additional motifs present in TMS6, TMS7 and TMS8; however, it possesses the [WI(YE)(RG)] sequence just before TMS6, between positions 406 and 409, and the [LS] motif located on TMS7, which is associated with Zn2+ P-type ATPases [32,33].

All of the results from this work provide predictive evidence for experimental studies to establish the ion specificity of MTBC P-type ATPases and their role in mycobacterial infection. It has been observed that some P-type ATPases of tubercle bacilli are over-expressed under conditions that mycobacteria face during infection [25,34,35]. For example, the expression level of CtpF and CtpC increase when M. tuberculosis is exposed to isoxyl, tetrahydrolipstatine and SRI#9190 antimicrobial compounds, indicating that these pumps might contribute to intrinsic resistance of mycobacteria to antimicrobial drugs [36]. Two additional studies indicate that hypoxia induces upregulation of the AEM transporter CtpF [37,38]; this observation is in agreement with those of another study that reported that CtpF is part of the regulon of the DosRS system, a relevant regulator of latency in M. tuberculosis[39]. In addition, upregulation of CtpF is observed in vitro when M. tuberculosis is incubated in the presence of S-nitrogluthatione GSNO, ethanol, H2O2 and nitric oxide [40]. Moreover, CtpA, CtpC, CtpG, CtpV and CtpF are also induced when M. tuberculosis is phagocytized by macrophages [25,34,35].


Mycobacteria are unicellular organisms that respond to environmental stimuli, and the transport of substances across the plasma membrane could play a fundamental role in their adaptability. Computational analysis shows that each MTBC species has a consistent aggrupation of the 12 P-type ATPases involved in ion transport. In this context, M. tuberculosis strains H37Ra and H37Rv share identical sequences for P-type ATPases, facilitating subsequent genetic studies using the attenuated strain H37Ra. The large number of HM P-type ATPases expressed by the MTBC strongly suggests that they could be essential for the bacteria to counteract the increased level of HM accumulated by macrophages after infection with tubercle bacilli. Thus, compensatory ion transport strategies could be used by mycobacteria to survive in host cells.

The different bioinformatics approaches used in this work to analyze the P-type ATPases identified in the MTBC are in agreement with the initial classification from the HMM search. The results obtained show that M. tuberculosis has the following three groups of P-type ATPases: HM transporters (CtpA, CtpB, CtpC, CtpD, CtpG, CtpJ and CtpV), AEM transporters (CtpE, CtpF, CtpH, and CtpI) and the KdpB protein, which corresponds to the β subunit of a multimeric K+ ATPase transporter exclusive to prokaryotes. Hydrophobicity analysis identified α-helix type TMS grouped into the following three topological types: type I (HM group), type II (AEM group) and type III (KdpB group). Interestingly, we report a possible mis-annotation for CtpH and CtpI in the TCDB Database, where they are classified as FUPA 24 type with two TMS, unlike the ten TMS identified for these unusually large transporters in this work. Finally, a counterpart of non-catalytic β subunits of Na+/K+ or H+/K+ ATPases does not exist within the MTBC proteomes.

Competing interest

There are no conflicts of interest to declare.

Authors' contributions

LN-A, AL-T, MP-R, JC-B, L-M S, DL, LM-R and CY-S wrote the manuscript, validated the tools and carried out the data analysis and interpretation. LN-A, AL-T, MP-R and JC-B contributed to the methodological design, supervised its development and critically revised the manuscript's content. All authors read and approved the final version of the manuscript.


This work was supported by the Dirección de Investigación Bogotá (DIB)-Universidad Nacional de Colombia, grants 12351, 12176, 11875, 11882 and Colciencias grant 13802. This research was supported in part by the Intramural Research Program of the National Institutes of Health, National Library of Medicine and National Center for Biotechnology Information.


This work was supported by the Dirección de Investigación Bogotá (DIB)-Universidad Nacional de Colombia, grants 12351, 12176, 11875, 11882 and Colciencias grant 13802. This research was supported in part by the Intramural Research Program of the National Institutes of Health, National Library of Medicine and National Center for Biotechnology Information. The funding sources had no involvement in the study design; collection, analysis and interpretation of data; the writing of the manuscript or the decision to submit the article for publication.

WHOGlobal tuberculosis control: WHO report 2011Year: 2011Switzerland: Publications of the World Health Organization246
McEvoy CR,et al. The role of IS6110 in the evolution of Mycobacterium tuberculosisTuberculosis (Edinb)Year: 200787539340410.1016/
Issa R,et al. Detection and discrimination of Mycobacterium tuberculosis complexDiagn Microbiol Infect DisYear: 201272162710.1016/j.diagmicrobio.2011.09.02122078904
Knechel NA,Tuberculosis: pathophysiology, clinical features, and diagnosisCrit Care NurseYear: 20092923443 quiz 44. 10.4037/ccn200996819339446
Tufariello JM,Chan J,Flynn JL,Latent tuberculosis: mechanisms of host and bacillus that contribute to persistent infectionLancet Infect DisYear: 2003395789010.1016/S1473-3099(03)00741-212954564
Koul A,et al. The challenge of new drug discovery for tuberculosisNatureYear: 201146973314839010.1038/nature0965721270886
Nagata T,et al. Comparative molecular biological analysis of membrane transport genes in organismsPlant Mol BiolYear: 20086665658510.1007/s11103-007-9287-z18293089
Axelsen KB,Palmgren MG,Evolution of substrate specificities in the P-type ATPase superfamilyJ Mol EvolYear: 19984618410110.1007/PL000062869419228
Pedersen PL,Transport ATPases into the year 2008: a brief overview related to types, structures, functions and roles in health and diseaseJ Bioenerg BiomembrYear: 2007395–63495518175209
Rocafull MA,et al. Isolation and cloning of the K+−independent, ouabain-insensitive Na+−ATPaseBiochim Biophys ActaYear: 201118086168470010.1016/j.bbamem.2011.02.01021334305
Pinoni SAL,A A,Na+ ATPase activities in chela muscle of the euryhaline crab Neohelice granulata: differential response to environmental salinityJ Exp Mar Bio EcolYear: 20093721–29197
Jorgensen PL,Purification and characterization of (Na+, K+)-ATPase. V. Conformational changes in the enzyme Transitions between the Na-form and the K-form studied with tryptic digestion as a toolBiochim Biophys ActaYear: 1975401339941510.1016/0005-2736(75)90239-4126698
Jorgensen PL,Hakansson KO,Karlish SJ,Structure and mechanism of Na, K-ATPase: functional sites and their interactionsAnnu Rev PhysiolYear: 2003658174910.1146/annurev.physiol.65.092101.14255812524462
Xu C,et al. A structural model for the catalytic cycle of Ca(2+)-ATPaseJ Mol BiolYear: 200231612011110.1006/jmbi.2001.533011829513
Palmgren MG,Nissen P,P-type ATPasesAnnu Rev BiophysYear: 2011402436610.1146/annurev.biophys.093008.13133121351879
Kuhlbrandt W,Biology, structure and mechanism of P-type ATPasesNat Rev Mol Cell BiolYear: 2004542829510.1038/nrm135415071553
Bublitz M,et al. In and out of the cation pumps: P-type ATPase structure revisitedCurr Opin Struct BiolYear: 2010204431910.1016/
Uniprot cited 2011; Available from:
Pirovano W,Feenstra KA,Heringa J,PRALINETM: a strategy for improved multiple alignment of transmembrane proteinsBioinformaticsYear: 2008244492710.1093/bioinformatics/btm63618174178
Eddy SR,Profile hidden Markov modelsBioinformaticsYear: 19981497556310.1093/bioinformatics/14.9.7559918945
Zhang Y,I-TASSER server for protein 3D structure predictionBMC BioinformaYear: 200894010.1186/1471-2105-9-40
Thever MD,Saier MH Jr,Bioinformatic characterization of P-type ATPases encoded within the fully sequenced genomes of 26 eukaryotesJ Membr BiolYear: 200922931153010.1007/s00232-009-9176-219548020
Chan H,et al. The P-type ATPase superfamilyJ Mol Microbiol BiotechnolYear: 2010191–2510420962537
Cole ST,Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequenceNatureYear: 199839353710.1038/311599634230
Botella H,et al. Mycobacterial p(1)-type ATPases mediate resistance to zinc poisoning in human macrophagesCell Host MicrobeYear: 20111032485910.1016/j.chom.2011.08.00621925112
Bramkamp M,Altendorf K,Greie JC,Common patterns and unique features of P-type ATPases: a comparative view on the KdpFABC complex from Escherichia coli (Review)Mol Membr BiolYear: 2007245–63758617710642
Lewinson O,Lee AT,Rees DC,A P-type ATPase importer that discriminates between essential and toxic transition metalsProc Natl Acad Sci U S AYear: 20091061246778210.1073/pnas.090066610619264958
Argüello JM,Identification of Ion-Selectivity Determinants in Heavy-Metal Transport P1B-type ATPasesJ Membr BiolYear: 20031959310810.1007/s00232-003-2048-214692449
Argüello JM,Eren E,González-Guerrero M,The structure and function of heavy metal transport P1B ATPasesBiometalsYear: 20072023324810.1007/s10534-006-9055-617219055
Argüello JM,Gonzalez-Guerrero M,Raimunda D,Bacterial transition metal P(1B)-ATPases: transport mechanism and roles in virulenceBiochemistryYear: 201150469940910.1021/bi201418k21999638
Ward SK,et al. CtpV: a putative copper exporter required for full virulence of Mycobacterium tuberculosisMol MicrobiolYear: 20107751096111010.1111/j.1365-2958.2010.07273.x20624225
Lewinson O,Lee AT,Rees DC,A P-type ATPase importer that discriminates between essential and toxic transition metalsProc Natl Acad SciYear: 2009106124677468210.1073/pnas.090066610619264958
Futai M,Wada Y,Kaplan JH,Handbook of ATPases: Biochemistry, Cell Biology, Pathophysiology. 1 edYear: 2004Weinheim: WILEY-VCH Verlag GmbH & Co. KGaA
Schnappinger D,et al. Transcriptional Adaptation of Mycobacterium tuberculosis within Macrophages: Insights into the Phagosomal EnvironmentJ Exp MedYear: 2003198569370410.1084/jem.2003084612953091
Kumar M,et al. Identification of Mycobacterium tuberculosis genes preferentially expressed during human infectionMicrob PathogYear: 201150131810.1016/j.micpath.2010.10.00321035536
Waddell SJ,et al. The use of microarray analysis to determine the gene expression profiles of Mycobacterium tuberculosis in response to anti-bacterial compoundsTuberculosis (Edinb)Year: 2004843–42637415207496
Bacon J,et al. The influence of reduced oxygen availability on pathogenicity and gene expression in Mycobacterium tuberculosisTuberculosis (Edinb)Year: 2004843–42051715207490
Sherman DR,et al. Regulation of the Mycobacterium tuberculosis hypoxic response gene encoding alpha -crystallinProc Natl Acad Sci U S AYear: 200198137534910.1073/pnas.12117249811416222
Kendall SL,et al. The Mycobacterium tuberculosis dosRS two-component system is induced by multiple stressesTuberculosis (Edinb)Year: 2004843–42475515207494
Voskuil MI,et al. Inhibition of respiration by nitric oxide induces a Mycobacterium tuberculosis dormancy programJ Exp MedYear: 200319857051310.1084/jem.2003020512953092


[Figure ID: F1]
Figure 1 

Classification of P-type ATPases of M. tuberculosis H37Rv according to ion specificity. Amino acid sequences of M. tuberculosis H37Rv P-type ATPases were compared with the HMM profiles built from characterized P-type ATPases reported by Uniprot. Grouping is based on the scores obtained from Hmmsearch. HM group (blue bars) includes the P-type ATPase transporters of heavy metals (CtpA, CtpB, CtpC, CtpD, CtpG, CtpJ, CtpV). The EAM group (red bars) includes transporters of alkaline and alkaline earth metals (CtpE, CtpF, CtpH, CtpI). The KdpB group (purple bars) includes the β subunit of K+ P-type ATPase. The scores obtained from the Hmmsearch tool are shown in the bottom of each plot.

[Figure ID: F2]
Figure 2 

Predicted topology for P-type ATPases of MTBC. Type I topology (HM group ATPases) (a), type II topology (EAM group ATPases) (b), type III topology (KdpB group ATPases) (c), Ion binding motifs (number 4) are indicated in the corresponding transmembrane domain.

[Figure ID: F3]
Figure 3 

Hydrophobicity profiles of M. tuberculosis H37Rv P-type ATPases. Hydrophobicity profiles allow grouping of P-type ATPases with seven different ion specificities.

[Figure ID: F4]
Figure 4 

The CtpH and CtpI mycobacterial P-type ATPases (type II topology) exhibit ten TMS in their C-terminal half (different color helixes). These tertiary structure models were generated with the I-TASSER tool and were modified with the PPM server [89] to include the possible location of the lipid bilayer. Dummies (green color) correspond to the location of the carbonyl groups in the bilayer. (a) CtpH model generated with the PPM server, (b) CtpI model built based on the TMDET results.

[Figure ID: F5]
Figure 5 

Conserved motifs for the P-type ATPases identified in M. tuberculosis H37Rv. The conserved motifs were compared with the motifs characteristic of the P-type ATPase superfamily. Identical residues, conserved substitutions, semi-conserved substitutions and unrelated residues are indicated in red, blue, green and black, respectively.

[Figure ID: F6]
Figure 6 

Alignments of consensus sequences from characterized PIB ATPases, with the HM P-type ATPases from M. tuberculosis H37Rv. (a) PIB-1 motifs, (b) PIB-2 motifs and (c) PIB-4 motifs

[TableWrap ID: T1] Table 1 

Additional conserved motifs observed in heavy metal transporter P-type ATPases

Phylogenetic subgroup Cation specificity TMS6 TMS7 TMS8 MTBC pump
CtpC *
CtpD *
CtpJ *
PIB-5 Undetermined TPCP QX4GX3SX3M PX5QEX2DX5N -

The additional motifs are located on TMS6, TMS7 and TMS8. *, pumps without a complete set of motifs (see discussion).

Article Categories:
  • Research Article

Keywords: Tuberculosis, Mycobacterium tuberculosis complex, P-type ATPases, Ion transport, Conserved motifs.

Previous Document:  Trigeminal neuralgia post-styloidectomy in Eagle syndrome: a case report.
Next Document:  An adjuvanted respiratory syncytial virus fusion protein induces protection in aged BALB/c mice.