- Research article
- Open Access
Protein–protein interaction of Rv0148 with Htdy and its predicted role towards drug resistance in Mycobacterium tuberculosis
BMC Microbiology volume 20, Article number: 93 (2020)
Mycobacterium tuberculosis resides inside host macrophages during infection and adapts to resilient stresses generated by the host immune system. As a response, M. tuberculosis codes for short-chain dehydrogenases/reductases (SDRs). These SDRs are nicotinamide adenine dinucleotide-reliant oxidoreductases involved in cell homeostasis. The precise function of oxidoreductases in bacteria especially M. tuberculosis were not fully explored. This study aimed to know the detail functional role of one of the oxidoreductase Rv0148 in M. tuberculosis.
In silico analysis revealed that Rv0148 interacts with Htdy (Rv3389) and the protein interactions were confirmed using far western blot. Gene knockout mutant of Rv0148 in M. tuberculosis was constructed by specialized transduction. Macrophage cell line infection with this knockout mutant showed increased expression of pro-inflammatory cytokines. This knockout mutant is sensitive to oxidative, nitrogen, redox and electron transport inhibitor stress agents. Drug susceptibility testing of the deletion mutant showed resistance to first-line drugs such as streptomycin and ethambutol and second-line aminoglycosides such as amikacin and kanamycin. Based on interactorme analysis for Rv0148 using STRING database, we identified 220 most probable interacting partners for Htdy protein. In the Rv0148 knockout mutants, high expression of htdy was observed and we hypothesize that this would have perturbed the interactome thus resulting in drug resistance. Finally, we propose that Rv0148 and Htdy are functionally interconnected and involved in drug resistance and cell homeostasis of M. tuberculosis.
Our study suggests that Rv0148 plays a significant role in various functional aspects such as intermediatory metabolism, stress, homeostasis and also in drug resistance.
Tuberculosis (TB) remains one of the world’s leading cause of death every year among infectious diseases . In 2018, WHO reported that 10.4 million new cases and 1.8 million deaths occur worldwide each year . With the increase in multidrug and extensive drug-resistance TB, there is a need to develop newer diagnostics and vaccines . It is well known that Mycobacterium tuberculosis spreads through aerosol, survives in oxygen and nutrition depleted environments and persists for long periods in the host cells [4,5,6]. The change in the host metabolism is thereby continuously monitored by M. tuberculosis to perform its active replication and persistence [7, 8]. M. tuberculosis during infection is exposed to various anti-bacterial agents secreted by macrophages. In addition, macrophages produce reactive oxygen species (ROS) and reactive nitrogen species (RNS) [9, 10]. The ROS interact directly and the resulting super oxides are converted into various oxidants like hypochlorite (HClO), peroxides (H2O2), peroxynitrite (ONOO−) and hydroxy radicals, which damage the bacterial cells . The bacterial cells resist stress using in-built mechanism. M. tuberculosis uses various defence mechanisms to damage the cell. OxyR and soxR are two prokaryotic regulators against peroxides and superoxides. Due to the absence of oxyR in M. tuberculosis, peroxidase stress is managed by alkylhydroperoxidase reductase (AhpC) by detoxifying the peroxides . In addition, other peroxidases such as KatG, superoxide dismutase, peroxiredoxin (AhpE) and thioredoxin reductase (Tpx) help in reducing peroxidase activity and also in controlling the oxidative and nitrosative stress conditions [13,14,15]. Though the components involved in the oxidative stress of M. tuberculosis have been identified, the functional importance remains unclear. Gene interaction and knockout studies thoroughly predict the functional interconnection between the genes and their role in the metabolism of bacteria.
In this study, for the first time, we attempted to predict the functional role of one of the hypothetical oxidoreductase Rv0148 of M. tuberculosis. Rv0148 is a hypothetical protein/putative short-chain dehydrogenase/reductase (SDR) and exhibits 100% similarity to M. bovis and M. africanum and 87% to M. avium . Rv0148 possesses the conserved SDR domain. As aminoglycosides bind to SDR sites, Rv0148 might neutralize the overexpression of aminoglycosides . Earlier studies reported that over expression of M. tuberculosis Rv0148 from multidrug-resistance isolates in Escherichia coli showed two- to three-folds of higher shift in MIC . Previous studies from our lab identified that PknI, one of the 11 serine-threonine protein kinases (STPKs), interacts with two proteins Rv2159c and Rv0148 . Rv2159c was characterized using gene knockdown studies, and its interaction with PknI increased its peroxidase activity several folds in the mutant strain .
As an extension to the previous study we have chosen Rv0148. It was characterized through bioinformatics tools using sequence and protein interaction analyses. In sequence analysis using Pfam database, we found Rv0148 carrying well-conserved nicotinamide adenine dinucleotide (NAD) domain, whereas other homologues possessed MaoC domain along with SDR. Furthermore, by in silico approach we identified that Rv0148 mainly interacts with hydroxyl acyl thioester dehydratase Htdy, a protein interaction which was confirmed by far western blotting (WB) and pull-down assay. We have constructed the gene knock out mutant of Rv0148 (Δ0148) by specialized transduction to understand the functional role of the gene. In vivo studies confirmed that this mutant induces pro-inflammatory cytokines and is susceptible to oxidative and nitrogen stress compounds. Δ0148 confers drug resistance to streptomycin, ethambutol, amikacin and kanamycin. This might be due to modified functional network in the absence of Rv0148. Overall, the study suggested that Rv0148, though a non-essential gene, is functionally involved in drug resistance and is interconnected with other drug-resistance genes. The absence of Rv0148 displays its prominent role in host immunity.
Sequence homologues of Rv0148
To identify the sequence homologues of Rv0148 in other organisms, blast search was performed against the non-reductant database in NCBI; 104 sequence homologues were selected and the sequences were reconfirmed based on short dehydrogenase domain in these sequences. Most of these sequences had the short dehydrogenase with a few sequences having both short dehydrogenase and MaoC dehydratase domain. A HMM profile based analysis (https://www.ebi.ac.uk/Tools/hmmer/search/hmmsearch) using the adh short domain profile from Pfam database identified 1270 significant sequences from Swissprot database (Table S1). Further, we searched the MaoC dehydratase sequence against mycobacterial genome and identified Htdy protein. Based on the analysis, we identified Rv0148 has only SDR. We searched for Rv0148 in STRING database and identified that Rv0148 interacts with Htdy protein (Fig. 1a). And this formed the basis for the subsequent analysis to predict the interaction of Rv0148 with Hdty using bioinformatics tools and then to validate the interaction using experimental methods.
Homology modelling of Rv0148
Because the structure of Rv0148 was not determined by experimental methods, we attempted to predict its 3D structure using the homology modelling approach. From the blast search results against PDB for a suitable template in complex with NAD, it was observed that Rv0148 showed the best alignment with human 17-beta-hydroxysteroid dehydrogenase type 4 in complex with NAD (PDB code 1ZBQ) and was chosen as the template. The sequence coverage between Rv0148 and 1ZBQ was 95% with a sequence identity of 52%, E-value of 2e-96) and gap of 2%. Both Rv0148 and 1ZBQ have short-chain dehydrogenase domain (PF00106). The 3D structure of Rv0148 was modelled (Fig. 2a) using the model-default python script in MODELLER software and the sequence alignment used for predicting the model structure is shown in Fig. 1b. The predicted model was validated using PROCHECK, profile-3D, and ProSA-web. Stereochemical properties of the predicted model were evaluated using PROCHECK to analyse the Ramachandran plot (Fig. 2b).
The plot revealed that the occurrence of 92.3% amino acids was in the favoured region, 7.3% in the generously allowed region and 0.4% in the disallowed region. Z-score was − 7.89 for the model structure determined using ProSA-web, a program used for validating protein structure, and was located within the space determined by X-ray crystallography (Fig. 2c). Thus, the overall results from structural superimposition, Ramachandran plot, verify-3D, and ProSA-web reveal that the predicted model is satisfactory and thus considered for further analysis. Profile-3D of the Rv1048 model was 90.36, thus suggesting that 90.36% of the residues have an average 3D–1D score ≥ 2 (Fig. 2d).
Docking of NAD within the active site of Rv0148
The docking of NAD co-factor was performed using GOLD software. However, a primary docking of the bound NAD to the template structure 1ZBQ was re-docked in its binding site to evaluate the docking process. The RMSD of the resulting docking of NAD to 1ZBQ was 0.86 Å, suggesting that the docking program is able to reproduce a structure very close to the experimental structure. In order to discover the putative binding pose of NAD co-factor in Rv0148, the binding site for Rv0148 was mapped based on 1ZBQ (Fig. 1b). The best-predicted pose of NAD within the active site of Rv0148 had a score of 95.44 Å (Fig. 1c).
Rv0148 and Htdy interaction analysis
The crystallographically determined structure of Htdy (PDB code 3KHP) and the predicted model structure of Rv0148 were used for docking using the ClusPro server. The structure from the largest cluster and with the lowest energy was selected and examined for the interaction between the two complexes (Fig. 3a). The residues interacting between the two proteins Rv0148 and Htdy from the docked complex were then mapped to identify the active site (Fig. 3b).
Cloning, expression and purification of Rv0148 and Htdy
The coding region of the two genes was amplified by PCR and the product matches the size of Rv0148 with 858 bp and htdy with 969 bp. The E. coli strain (BL21) expressed pGBRJ3 Htdy as His-tagged protein and pGBRJ2 0148 as GST-tagged protein (Fig. S1). The protein was purified. The mass of Htdy is 31 kDa and Rv0148 GST is 55 kDa. These correspond to the expected mass of the above proteins through western blot.
Interaction of Rv0148 and Htdy
The pull-down assay with recombinant His-Htdy (PGBRJ3) lysate using Ni-NTA affinity chromatography confirmed that the Htdy was able to pull down GST-tagged Rv0148 (PGBRJ2). Simultaneously, the GST-tagged lysate was able to pull down His tag recombinant Htdy (Fig. 4). The pull-down assay and far western blot confirmed that Rv0148 interacts with Htdy protein.
Generation of ∆0148 knockout mutant in M. tuberculosis
We constructed knockout mutant of M. tuberculosis for not expressing Rv0148 by phage-facilitated allelic exchange . The coding region of the gene was replaced by the hygromycin cassette and four fragments were generated with 3.6, 1.6, 0.837 and 0.605 bp product sizes (Fig. S1). The 50 kb phAE159 shuttle phasmid with AES retains 6.642 kb (Fig. S2). The generated high-titre transducing phages were used for specialized transduction and the knockout mutant was generated successfully and named as Δ0148 (Fig. 5). We confirmed the knockout mutant using sequencing, PCR for hygromycin cassette and also through real-time PCR (RT-PCR) using gene-specific primers (Table S3).
Colony morphology and growth kinetics
Gene disruption shows no distinct effect in the cell development of the knockout strain compared to H37Rv. To elucidate if the knockout of Rv0148 confers any change in vitro, we compared the growth kinetics of H37Rv, ∆0148 and CΔ0148 in 7H9-enriched media. All the strains showed similar growth and there was no distinct difference observed in in vitro survival of ∆0148.
Cytokines analysis using ELISA
The cell-free supernatants after post-infection were analysed for the expression of pro-inflammatory cytokines IL-1β, IL-6 and TNF-α and anti-inflammatory cytokines IL-10, IL-12p40, IL-12p70 and IL-4. Enhanced expression of IL-1β, IL-6 and TNF-α cytokines was observed in ∆0148 infected with THP-1 macrophages (Fig. 6a–c). Similarly, reduced level of expression of IL-10, IL-12p40, IL-12p70 and IL-4 was observed in ∆0148 compared to wild-type H37Rv (data not shown). The expression was analysed using appropriate statistical methods.
Effect on inhibition of the electron transport chain in ∆0148
Rv0148 is predicted to be involved in the intermediatory metabolism of M. tuberculosis. The ETC plays an important role in the production of proton motive force (PMF) in the metabolism of bacteria. ETC consists of five complexes I, II, III, IV and V, and inhibition of these complexes could accommodate changes in the growth of bacteria. Rotenone is an inhibitor of complex I whereas oxaloacetate inhibits complex II. It was observed that the exposure of Δ0148 to 140 μM rotenone and 2 mM oxaloacetate significantly decreases the growth (Fig. 7a and b). The result suggests that ∆0148 strain has a distinct role in the ETC. The significance was reported using statistical tools.
Effect of oxidative and nitrosative stress in ∆0148
In diseased patients, mycobacterium resides in the alveolar macrophages in the initial stage of infection and secretes various ROS and RNS. We investigated whether these stresses affect the metabolism of ∆0148 strain when exposed to 5 mM hydrogen peroxide (H2O2), 5 mM cumene hydroperoxide, 2 mM thiol-reductant DTT and 1% sulphoxide (DMSO). ∆0148 is sensitive to oxidative stress components and to nitrogen derivatives and there is decrease in growth of Δ0148 compared to wild-type H37Rv (Fig. 8a–d). It was concluded that the growth of ∆0148 declined on exposure to ETC, oxidative and nitrosative stress agents. The result was considered using standard statistical methods.
Sensitivity and resistance to first-line and second-line drugs
BACTEC MGIT was performed with first-line and second-line drugs and shows that ∆0148 is resistant to streptomycin, ethambutol, amikacin, and kanamycin. H37Rv was maintained as a control for all these drug sensitivity assays (Table 1).
The STITCH database was used to identify M. tuberculosis proteins that interact with three anti-TB drugs streptomycin, kanamycin and amikacin. For each of the proteins interacting to the respective drugs identified from the STITCH database, we collected the interactomes with high confidence > 0.7 from STRING database to ensure a level of confidence for the interactions. Simultaneously, we collected the highest confidence interactomes for Rv0148, htdy and EIS genes from STRING (Table 2). Finally, we produced an interactome of 220 proteins that connect Rv0148, Htdy and EIS with the proteins that either interact directly or indirectly with the three drugs. The M. tuberculosis PPI (protein–protein interaction) network incorporating all the identified proteins is shown in Fig. 9. Comparing and analysing the intersection of proteins that connect the three drugs kanamycin, amikacin and streptomycin and the three proteins Rv0148, Htdy and EIS resulted in 57 proteins (Fig. S3). Analysing the 220 genes, except for streptomycin, we identified genes connected to kanamycin and amikacin as an intersection between interacting genes of Rv0148, htdy and EIS.
Differential expression of htdy
The interactome analysis predicted that htdy is also involved in drug resistance. Expression of htdy in the Δ0148 strain of M. tuberculosis was analysed using qPCR as mentioned earlier. The htdy expression in Δ0148 increased to 1.25-fold whereas the wild-type strain showed 0.7-fold (Fig. 10 and Table S4). The fold expression was calculated using SDS 2.4 software and relative quantity (RQ value) method.
In this study, we report the cloning, overexpression, purification, gene knockout and interacting partners of Rv0148, a hypothetical protein that belongs to the oxidoreductase family of M. tuberculosis. Through bio-informatics approaches, we concluded that Rv0148 is a NAD oxidoreductase protein with SDR domain. Rv0148 is predicted as interacting with a protein that has an MaoC dehydratase domain and the subsequent blast analysis led to the identification of Htdy protein. Further, the string analysis of Rv0148 displayed 10 interacting partners such as Rv3538, Rv0149, mbtB, fas, fadB, aacA2, Rv0147, Rv2228c, muT1 and htdy thus confirming our analysis. The predicted Htdy as an interacting partner of Rv0148 was further confirmed through far WB and pull-down assays.
We generated gene knockout mutant of Rv0148 in M. tuberculosis and show that Rv0148 is a non-essential gene as confirmed in the previous report . ∆0148 grew similarly as the wild type under constant shaking in enriched 7H9 media. This indicates clearly that the loss of Rv0148 does not lead to any change in the growth of M. tuberculosis. There were no correlative morphological changes in cell structure owing to the deletion of Rv0148.
We infected the macrophage cell lines with mutants for cytokine profiling. Increased expression of pro-inflammatory cytokines predicts that Δ0148 might be involved in host immune responses [23, 24]. The lower level expression of anti-inflammatory cytokines IL-10, IL-4, IL-12p40 and IL-12p70 indicates that Δ0148 has no role in inflammation .
In aerobic organisms, respiratory chain possesses inner membrane complexes associated with energy production. Complex I known as NADH quinone oxidoreductase transfers electrons from NADH to quinone associated with proton transfer across the membrane. Complex II is involved in tricarboxylic acid cycle without transfer of protons. Complex III transfers electrons from quinone to cytochrome c. Complexes III and IV receive these electrons and transfer them to oxygen. Complex V using PMF synthesizes ATP molecules . Rv0148 is a NAD oxidoreductase involved in the function of intermediary metabolism and respiration. The exposure of Δ0148 to components inhibiting ETC complexes I and II  revealed that ∆0148 responds to changes in the electron transport system, thus predicting its role in cell homeostasis. M. tuberculosis adapts to various stress conditions to survive as a successful pathogen inside the host . Δ0148 interestingly showed decreased growth in hydrogen peroxide, CHP, DTT and DMSO, which indicates its sensitivity to oxidative and nitrogen stress. Thus, in response to stress, the differential growth of ∆0148 demonstrates that the gene plays an important role in the protective efficacy against oxidative and nitrogen stress .
Contradictory to previous findings, Rv0148 showed increased expression in MDR isolates . In the present study, Δ0148 showed resistance to first-line streptomycin and ethambutol and second-line amikacin and kanamycin drugs and the result was correlated with its interaction towards drug-resistance genes. Interactome analysis in the study suggests that Rv0148 interacts with kanamycin resistance genes Rv3538, mutT1, fas, Rv0130, Rv3389 and htdx. The EIS gene is reported earlier to be involved in causing drug resistance of M. tuberculosis . Our analysis reveals that Rv0148 interacts with eis genes accA2, fas, kasA, fabD, KasB and accD3. Rv0148 does not display any direct interaction with amikacin and ethambutol genes; however, the interactomes of EIS and amikacin share common interacting genes such as Rv2417c, Rv2414c, Rv3055, Rv1312, Rv1571, Rv3605c, Rv3839, Rv0428c, Rv0495c and Rv2242. As these proteins are connected to each other indirectly on a larger interactome space, we assume that drug resistance in Δ0148 strain could be due to some disruption of the interactome network.
The major findings from the current study prove that htdy is also involved in drug resistance and the increased expression of htdy in ∆0148 predicts that Rv0148 altered the expression of interacting partners in M. tuberculosis. Interactome analysis using STRING database predicts that Rv0148 has a strong interaction with kanamycin, eis and amikacin resistance genes. Further, interactome analysis reveals that Rv0148 and Htdy are connected to EIS and the interacting proteins of the drugs analysed in this study. Thus, we assumed that the absence of Rv0148 affects the interactome that might bring a major expression difference in these proteins resulting in drug resistance towards respective drugs. An additional transcriptome profiling of these strains will further help in identifying a clear mechanism towards drug resistance. Further, studies like double knockout mutant of Rv0148 and Rv3389 will provide the detailed functional role of these genes.
To summarize, Rv0148 is a non-essential gene belonging to the oxidoreductase family in M. tuberculosis. In silico approach and protein interaction studies reveal that Rv0148 interacts with Htdy, which is a protein of the hydroxyl acyl family. Htdy in the presence of Rv0148 is quite stable and fold expression of the gene is optimal. In the absence of Rv0148, increased fold expression of htdy proves that these two genes are functionally interconnected and furthermore it is experimentally proved. In addition to these findings, Rv0148 not only interacts with htdy, but also has a strong interactome network of 220 genes in M. tuberculosis with a high confidence of 0.7. The absence of Rv0148 showed resistance to first-line drugs streptomycin and ethambutol and second-line drugs amikacin and kanamycin. Increased expression of pro-inflammatory cytokines suggests that Δ0148 might play a role in targeting the immune response. Our study suggests that Rv0148 plays a significant role in various functional aspects such as intermediatory metabolism, stress and homeostasis.
Template selection and homology modelling
A suitable template for Rv0148 was identified by searching in the PDB database using the default settings of BLAST program in February, 2017. The sequence of Rv0148 was searched against Pfam database to further identify available PDB structures having the same domain. The final template (PDB Id: 1ZBQ) was chosen based on sequences having the highest sequence identity and domain coverage against Rv0148 sequence. The template 1ZBQ selected for modelling Rv0148, is the experimentally determined crystal structure of human 17-beta-hydroxysteroid dehydrogenase in complex with cofactor NAD. The template and the target sequence were aligned using the malign program in Modeller software . Homology model of Rv0148 was constructed using Modeller and the quality of the model was further analysed using PROCHECK , ProSA-web  and verify 3D . The model was also validated using root mean square deviation (RMSD).
Docking was performed using NAD cofactor as the ligand with the predicted model of Rv0148 to understand its mode of interaction within the binding site. GOLD v5.2 (genetic optimization for ligand docking) is a widely used protein–ligand docking program . GOLD uses a genetic algorithm methodology for protein–ligand docking that allows full ligand and partial protein flexibility. The homology model of Rv0148 protein was used as the receptor after the addition of hydrogen atoms before docking. The NAD ligand was corrected with Auto Edit Ligand option of GOLD program. Default genetic algorithm (GA) settings that ensure 100% search efficiency were used for docking. A total of 20 poses of the docked molecule were generated. Early termination of the number of GA runs was allowed when the RMSDs of the top three GA solutions were within 1.5 Å. The best pose of the docked ligand was selected based on Goldscore.
The interactions between Rv0148 and Htdy structures were docked using ClusPro . The ClusPro algorithm uses fast fourier transform correlation, making it feasible to generate billions of docked conformations by simple scoring functions. ClusPro filters the docked confirmations with near-native structures, ranks them based on their clustering properties and finally refines a limited number of structures by energy minimization. The server outputs the top predicted complexes based on energy and cluster size.
Bacterial strains and growth conditions
E. coli strains DH5α and BL21 (DE3) (Invitrogen, USA) were used for cloning and expression of recombinant proteins. E. coli cells were grown in Luria Bertani (LB) medium with constant shaking at 37 °C. M. smegmatis (mc2155) was grown in Middlebrook 7H9 broth (Difco, USA) with 0.2% glycerol and 0.05% Tween-80. The media was supplemented with 50 μg/mL ampicillin for E. coli whenever needed. M. tuberculosis was also grown in 7H9 broth with antibiotics carbenicillin 50 μg mL, cyclohexamide 10 μg mL, kanamycin 25 μg mL and hygromycin 50 μg mL, wherever required, along with OADC and Tween-80. Primers used in the study are listed in Table S2. The plasmids and strains constructed in the present study are mentioned in Table 3.
Cloning and expression of Rv0148 and htdy
M. tuberculosis genomic DNA was isolated using CTAB-NaCl standard method . Coding region of Rv0148 from M. tuberculosis was amplified by polymerase chain reaction (PCR) using gene-specific oligonucleotides. The PCR product of htdy was further cloned into pRSET-B harbouring T7 promoter and Rv0148 in pGEX 4 T-1 vector harbouring tac promoter with Bam HI at 5′ and Xho I at 3′ end as restriction sites. The resultant clone was confirmed through restriction digestion and DNA sequencing. For recombinant protein expression, pGBRJ3 and pGBRJ2 clones were transformed into E. coli BL21 (DE3) cells and protein expression was induced with 0.5 mM isopropyl β-d-1 thiogalactopyranoside (Invitrogen, USA). Induced cells were grown for 4 h at 37 °C and harvested by centrifugation at 8000 rpm for 10 min. The cell pellet was re-suspended in ice-cold lysis buffer G (50 mM Tris pH 8, 300 mM NaCl, 20% glycerol,0.5 mM phenyl methyl sulphonyl fluoride and 10 mM imidazole) along with protease inhibitor cocktail Sigma fast (SIGMA, USA) and sonicated on ice using ultrasonic homogenizer (5 cycles: 30 s pulse and 1 min interval). Finally, the lysate was centrifuged at 10,000 rpm for 10 min at 4 °C to separate soluble and insoluble proteins.
Purification of Rv0148 and Htdy
The soluble lysate proteins were analysed on 12% SDS gel for confirmation of recombinant proteins and protein lysis. Lysate was purified using affinity column (Bio-RAD, CA, USA) and 6X histidine-tagged protein was bound to NI-NTA resin (Invitrogen, USA) and incubated at 12 °C for 2 h. The column was washed with 10 column volumes of wash buffer containing 40 mM imidazole and was eluted with a buffer containing 200 mM imidazole. Similarly, GST-tagged protein was bound to glutathione sepharose resin (GST, Takara, Japan) and eluted with 20 mM reduced glutathione. The eluted fractions were then dialysed and analysed by SDS-PAGE and WB.
Pull-down assay was performed using poly-prep chromatography column (Bio-Rad, USA) with NI-NTA resin. Recombinant pGBRJ3 His lysate was incubated with NI-NTA resin for 2 h at 12 °C with gentle agitation. The column was washed with buffer containing 40 mM imidazole and incubated overnight with pGBRJ2 GST lysate from E. coli BL21 (DE3) cells at 12 °C with gentle agitation at 10 rpm in an end-to-end shaker. Unbound proteins were removed by washing with 10 column volumes of buffer. The interacting proteins were eluted with buffer containing 200 mM imidazole.
Far western blot
To recognize the Rv0148 interacting proteins, the fraction eluted in pull-down assay was separated on a 12% SDS-PAGE and transferred to Immobilon-P (PVDF) membrane. The membrane was denatured by incubation with varying concentrations (6, 3, 1 and 0.1 M) of guanidine-HCl in AC buffer containing 5 M NaCl, 1 M Tris pH 7.5, 0.5 M EDTA, 10% Tween-20, 2% skimmed milk powder, 10% glycerol and 1 M dithiothreitol (DTT) for 30 min at room temperature followed by 5 min wash with phosphate buffer saline Tween 20 . Then, the blot was completely renatured in AC buffer without guanidine-HCl and left overnight at 8 °C. The membrane was then blocked with 5% skimmed milk and incubated with bait protein for 7–8 h at 8 °C and further with 1:6000 anti-His monoclonal antibody and developed using enhanced chemiluminescent kit (Thermo Scientific, USA).
Construction of allelic exchange substrate
Gene disruption was carried out using Rv0148 left homology sequence (LHS) and right homology sequence (RHS) flanking regions carrying van911 as restriction site targeted for deletion and substitution. The gene was cloned onto either side of the selectable antibiotic cassette [γδ (sacB-hyg)γδ] to generate allelic exchange substrate for homologous recombination. Van911 restriction enzyme identifies a discontinuous palindrome interjected by five bases of sequence (CCAN_NNN^NTGG) to achieve a four-fragment ligation.
Construction of specialized transducing phages
Shuttle phasmid phAE159 proliferates at 30 °C and does not propagate in the infected mycobacterial cell at 37 °C. Shuttle phasmid phAE159 accommodates up to 10 kb of recombinant DNA into mycobacterial cells. The mycobacteriophage carrying AES disrupts a precise gene in mycobacteria through deletion–substitution. The AES with hygromycin resistance marker differs from the phasmid, and the use of in vitro lambda packaging extract (Epicentre, USA) helps in increasing the efficiency of transformation and also in the selection of phasmid DNA. Then the phasmid DNA was electroporated into M. smegmatis to produce specialized transducing phages (STPs) at permissive temperature of 30 °C. These phages were further propagated using MP buffer (1 M Tris pH 8, 5 M NaCl, 1 M MgCl2, 1 M CaCl2) to obtain high-titre allele-specific STPs.
Specialized transduction was performed at 37 °C for the augmentation of STPs. In this temperature transfer of DNA to the host (mycobacterium) occurs with homologous recombination. The transduced M. tuberculosis cells were selected on hygromycin plate (150 μg/mL) after transduction. A minimum of five colonies appeared on the plate and was confirmed by PCR and Sanger sequencing using hyg as a forward primer and RHS as a reverse primer to confirm homologous recombination. For complementation of Rv0148, the gene was cloned into pMV261 vector harbouring hsp60 promoter. The construct was then electroporated into the ∆0148 strain and colonies were screened on the 7H10 plate with antibiotic selection (kanamycin 20 mg mL). Complementation was confirmed by PCR using Rv0148 primers and the complemented strain labelled as CΔ0148.
In vitro growth kinetics of ∆0148
To check whether the disruption of Rv0148 fetched any changes in the in vitro growth of H37Rv, ∆0148 and CΔ0148 strains, log-phase cultures were grown in Middlebrook 7H9-OADC-Tween media and incubated at 37 °C with 180 rpm shaking. Aliquots of 100 μl culture were taken at different time points and OD values were measured at 600 nm using spectrophotometer (Eppendorf) on days 0, 3, 6, 9, 12 and 17. Survival of wild type, ∆0148 and CΔ0148 strains was analysed by plating the serial diluted broth cultures in 7H10 OADC agar without Tween-80. The plates were incubated for 3–4 weeks and the colony forming units were calculated.
The colony morphology of the mutant was determined by comparing with complementary strain and wild type. The strains were grown in liquid 7H9 media to mid-log phase and 10 μl of culture was inoculated to 7H10 agar plate supplemented with hygromycin. The plates were incubated at 37 °C for 3–4 weeks for colony formation.
RNA extraction and qPCR
Bacterial culture 100 mL was grown in 7H9 media to mid-log phase and the growth was arrested in ice and centrifuged at 4 °C. The pellet was suspended in trizol (Invitrogen, USA). Cells were disrupted using 0.1 mm zirconium beads in a mini bead beater. The total RNA was extracted using the RNeasy purification kit (Qiagen, Germany) according to the manufacturer’s protocol and quantified using ND-1000 Nanodrop spectrophotometer (Nanodrop Technologies). It was stored at − 80 °C. To determine the relative mRNA through qPCR, the first strand cDNA was synthesized from 1 μg of RNA using Maxima first strand cDNA synthesis kit (Thermo Scientific, USA). Quantitative qPCR primers for the Rv0148, htdy and 16S rRNA genes were designed (Table S2). qPCR was carried out using SYBR Green Master Mix (Thermo Scientific, USA) according to manufacturer’s instructions using the Applied Biosystems 7300 real-time PCR system. Control reactions for each sample were carried out in the absence of reverse transcriptase to look over for DNA contamination. The amplification conditions for all reactions were 1 cycle of 50 °C for 2 min, 95 °C for 10 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. qPCR data analysis was carried out by relative quantification of Rv0148, htdy and CΔ0148 gene expression using comparative calculated threshold (CT) method. For each qPCR run, the CT cycle was normalized to the internal control 16S rRNA gene amplified from the corresponding sample. Statistical analysis was carried out using SDS 2.4 software and relative quantity (RQ value) method. The data presented are averages of three independent experiments.
Infection of knockout strain in THP-1 cell lines
Intracellular viability of the knockout strain was determined using macrophage cell line infection studies. THP-1 cells were grown in RPMI media with 10% fetal bovine serum. Cells were further grown to reach 1 × 106 cells, seeded onto 24-well plates and differentiated into macrophages using 50 mM phorbol 12-myristate 13-acetate. Plates were incubated at 37 °C in the presence of 5% CO2 incubator for 2 days, washed with RPMI medium and left overnight. Macrophages were then infected in triplicate with H37Rv, ∆0148 and CΔ0148 with multiplicity of infection of 1:10. Phagocytosis was allowed to take place for 4 h, after which the infected macrophages were incubated with RPMI containing streptomycin to eliminate extracellular bacteria. Further, the cell-free culture supernatants were collected on days 0, 1, 3, 5 and 7 post-infection and stored at − 80 °C for cytokine analysis.
The supernatants derived from THP-1 cell lines infected with H37Rv, ∆0148 and CΔ0148 at different intervals were evaluated for expression of cytokines IL-4, IL-6, IL-10, IL-1β, TNF, IL-2, INF-γ, IL-12 (p40) and IL-12 (p70) by sandwich ELISA using BD ELISA kit (BD pharmingen, USA) according to manufacturer’s protocol.
Role of knockout mutant during oxidative stress
To understand the role of oxidative and nitrogen stress in H37Rv, ∆0148 and CΔ0148, the strains were grown to mid-log phase in 7H9-OADC-T medium and exposed to various stress components such as electron transport chain (ETC) inhibitors, peroxides, sulphoxides, thiol stress with slight modification to the procedure by  and OD values were measured at 600 nm (OD600) for 4, 24 and 48 h time points. The graphs were plotted and significance was reported considering P < 0.01 and P < 0.001.
MGIT and drug susceptibility testing
The H37Rv and knockout strains were grown on LJ solid media, 0.5 McFarland culture suspension was used for MGIT inoculation and drug susceptibility testing (DST) was performed for all the culture-positive MGIT tubes after 4 days. The growth control (GC) was prepared by diluting the cultures in the ratio of 1:100. The test cultures were prepared using saline in 1:5 dilution; 800 μL of supplement was added to tubes according to manufacturer's protocol; 100µL streptomycin 1 μg mL, isoniazid 0.1 μg mL, rifampicin 1 μg mL and ethambutol 5 μg mL (SIRE kit, BD). Kanamycin 2.5 μg mL, ofloxacin 2 μg mL and amikacin 1 μg mL was added to the respective tubes; and finally 500 μL of culture was added to all the tubes and loaded in BACTEC MGIT 960. Tubes were recapped, mixed well and fixed into the BACTEC MGIT 960 using the antimicrobial susceptibility testing set entry protocol. BACTEC MGIT 960 instrument monitored susceptibility test until the growth unit of the GC reaches 400 .
Interactome analysis to predict drug-resistance partners
To predict the interacting genes involved in drug resistance, STITCH database was used and further analysed using a bioinformatics tool STRING.
Expression of htdy in the presence and absence of Rv0148
H37Rv, mutant and complimentary strains were grown in 7H9 media to log phase according to the growth conditions mentioned earlier and RNA was extracted. Individual and co-expression of htdy were analysed through qPCR in mutant and complementary strains and H37Rv. The Ct and R values were further considered for analysing the folds of gene expression.
The data represented all the three independent experiments carried out using replicates. Mean, standard error mean (SEM) and Two-way ANOVA was used to analyse the data obtained from in vitro kinetics, ELISA and oxidative stress assays (**P < 0.01 was said to be significant and ***P < 0.001 was said to be highly significant). Real time PCR data was analysed using SDS software and relative quantity method.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
- M. tuberculosis :
Reactive nitrogen species
Reactive oxygen species
Minimum inhibitory concentration
Serine-threonine protein kinases
Nicotinamide adenine dinucleotide
Real-time Polymerase chain reaction
Enzyme-linked immunosorbent assay
Mycobacteria growth indicator tube
Kyu HH, Maddison ER, Henry NJ, Mumford JE, Mumford JE, Shields C. The global burden of tuberculosis: results from the global burden of disease study. Lancet Infect Dis. 2017;23:261-84.
WHO. Global Tuberculosis Report: World Health Organisation; 2018. https://apps.who.int/iris/handle/10665/274453.
Boshoff HI, Barry CE 3rd. Tuberculosis - metabolism and respiration in the absence of growth. Nat Rev Microbiol. 2005;3(1):70–80.
Vergne I, Chua J, Singh SB, Deretic V. Cell biology of mycobacterium tuberculosis phagosome. Annu Rev Cell Dev Biol. 2004;20:367–94.
Uribe-Querol E, Rosales C. Control of phagocytosis by microbial pathogens. Front Immunol. 2017;8:1368.
Smith I. Mycobacterium tuberculosis pathogenesis and molecular determinants of virulence. Clin Microbiol Rev. 2003;16(3):463–96.
Kumar A, Farhana A, Guidry L, Saini V, Hondalus M, Steyn AJ. Redox homeostasis in mycobacteria: the key to tuberculosis control? Expert Rev Mol Med. 2011;13:e39.
Trivedi A, Singh N, Bhat SA, Gupta P, Kumar A. Redox biology of tuberculosis pathogenesis. Adv Microb Physiol. 2012;60:263–324.
Voskuil MI, Bartek IL, Visconti K, Schoolnik GK. The response of mycobacterium tuberculosis to reactive oxygen and nitrogen species. Front Microbiol. 2011;2:105.
MacMicking JD, North RJ, LaCourse R, Mudgett JS, Shah SK, Nathan CF. Identification of nitric oxide synthase as a protective locus against tuberculosis. Proc Natl Acad Sci U S A. 1997;94(10):5243–8.
Nambi S, Long JE, Mishra BB, Baker R, Murphy KC, Olive AJ, Nguyen HP, Shaffer SA. The Oxidative Stress Network of Mycobacterium tuberculosis Reveals Coordination between Radical Detoxification Systems. Cell Host Microbe. 2015;17(6):829–37.
Deretic V, Philipp W, Dhandayuthapani S, Mudd MH, Curcic R, Garbe T, Heym B, Via LE, Cole ST. Mycobacterium tuberculosis is a natural mutant with an inactivated oxidative-stress regulatory gene: implications for sensitivity to isoniazid. Mol Microbiol. 1995;17(5):889–900.
Edwards KM, Cynamon MH, Voladri RK, Hager CC, DeStefano MS, Tham KT, Lakey DL, Bochan MR, Kernodle DS. Iron-cofactored superoxide dismutase inhibits host responses to Mycobacterium tuberculosis. Am J Respir Crit Care Med. 2001;164(12):2213–9.
Jaeger T, Budde H, Flohe L, Menge U, Singh M, Trujillo M, Radi R. Multiple thioredoxin-mediated routes to detoxify hydroperoxides in Mycobacterium tuberculosis. Arch Biochem Biophys. 2004;423(1):182–91.
Sherman DR, Mdluli K, Hickey MJ, Arain TM, Morris SL, Barry CE 3rd, Stover CK. Compensatory ahpC gene expression in isoniazid-resistant Mycobacterium tuberculosis. Science. 1996;272(5268):1641–3.
Sharma D, Kumar B, Lata M, Joshi B, Venkatesan K, Shukla S, Bisht D. Comparative proteomic analysis of aminoglycosides resistant and susceptible Mycobacterium tuberculosis clinical isolates for exploring potential drug targets. PLoS One. 2015;10(10):e0139414.
Sharma D, Singh R, Deo N, Bisht D. Interactome analysis of Rv0148 to predict potential targets and their pathways linked to aminoglycosides drug resistance: an insilico approach. Microb Pathog. 2018;121:179–83.
Sharma D, Lata M, Faheem M, Khan AU, Joshi B, Venkatesan K, Shukla S, Bisht D: Cloning Expression and Correlation of Rv0148 to Amikacin & Kanamycin Resistance. Curr Proteomics 2015;12:96-100.
Venkatesan A, Hassan S, Palaniyandi K, Narayanan S. In silico and experimental validation of protein-protein interactions between PknI and Rv2159c from Mycobacterium tuberculosis. J Mol Graph Model. 2015;62:283–93.
Venkatesan A, Palaniyandi K, Sharma D, Bisht D, Narayanan S. Functional characterization of PknI-Rv2159c interaction in redox homeostasis of Mycobacterium tuberculosis. Front Microbiol. 2016;7:1654.
Jain P, Hsu T, Arai M, Biermann K, Thaler DS, Nguyen A, Gonzalez PA, Tufariello JM, Kriakov J, Chen B, et al. Specialized transduction designed for precise high-throughput unmarked deletions in Mycobacterium tuberculosis. mBio. 2014;5(3):e01245–14.
Sassetti CM, Boyd DH, Rubin EJ. Genes required for mycobacterial growth defined by high density mutagenesis. Mol Microbiol. 2003;48(1):77–84.
Russell-Goldman E, Xu J, Wang X, Chan J, Tufariello JM. A Mycobacterium tuberculosis Rpf double-knockout strain exhibits profound defects in reactivation from chronic tuberculosis and innate immunity phenotypes. Infect Immun. 2008;76(9):4269–81.
Cooper AM, Mayer-Barber KD, Sher A. Role of innate cytokines in mycobacterial infection. Mucosal Immunol. 2011;4(3):252–60.
Couper KN, Blount DG, Riley EM. IL-10: the master regulator of immunity to infection. J Immunol. 2008;180(9):5771–7.
Ana MP, Melo TMB, Teixeira M. New insights into type II NAD(P)H:Quinone Oxidoreductases; 2004.
Bhat SA, Iqbal IK, Kumar A. Imaging the NADH:NAD(+) homeostasis for understanding the metabolic response of Mycobacterium to physiologically relevant stresses. Front Cell Infect Microbiol. 2016;6:145.
Yano T, Li LS, Weinstein E, Teh JS, Rubin H. Steady-state kinetics and inhibitory action of antitubercular phenothiazines on mycobacterium tuberculosis type-II NADH-menaquinone oxidoreductase (NDH-2). J Biol Chem. 2006;281(17):11456–63.
Bulatovic VM, Wengenack NL, Uhl JR, Hall L, Roberts GD, Cockerill FR 3rd, Rusnak F. Oxidative stress increases susceptibility of Mycobacterium tuberculosis to isoniazid. Antimicrob Agents Chemother. 2002;46(9):2765–71.
Gikalo MB, Nosova EY, Krylova LY, Moroz AM. The role of eis mutations in the development of kanamycin resistance in Mycobacterium tuberculosis isolates from the Moscow region. J Antimicrob Chemother. 2012;67(9):2107–9.
Hatherley R, Brown DK, Glenister M, Tastan Bishop O. PRIMO: an interactive homology modeling pipeline. PLoS One. 2016;11(11):e0166698.
Laskowski RA, Jablonska J, Pravda L, Varekova RS, Thornton JM. PDBsum: structural summaries of PDB entries. Protein Sci. 2018;27(1):129–34.
Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35(Web Server issue):W407–10.
Cloete R, Kapp E, Joubert J, Christoffels A, Malan SF. Molecular modelling and simulation studies of the Mycobacterium tuberculosis multidrug efflux pump protein Rv1258c. PLoS One. 2018;13(11):e0207605.
Kozakov D, Hall DR, Xia B, Porter KA, Padhorny D, Yueh C, Beglov D, Vajda S. The ClusPro web server for protein-protein docking. Nat Protoc. 2017;12(2):255–78.
Baess I. Isolation and purification of deoxyribonucleic acid from mycobacteria. Acta Pathol Microbiol Scand B: Microbiol Immunol. 1974;82(6):780–4.
Wu Y, Li Q, Chen XZ. Detecting protein-protein interactions by far western blotting. Nat Protoc. 2007;2(12):3278–84.
NIRT: NIRT standard operating protocols. 2010.
This work was supported by the ICMR-NIRT intramural research programme. GB awarded with PhD fellowship from the Department of Science and Technology INSPIRE programme. None of the funding bodies were involved in the design of the study, the collection, analysis, or interpretation of data, or in writing and communication of the manuscript.
Ethics approval and consent to participate
No clinical samples were used in this study as this was an entirely in-vitro study using bacterial strains ethical review is not required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 3: Figure S1. Restriction enzyme digestion of Rv0148 &Rv3389 and cloning. (A) Confirmation of recombinant clones in pGEX4T-1 M: 1 kb marker, Lane 1 to 4; pGEX4T-1 as control with 4.9 kb size, Lane 5to 8 recombinant clones, Lane 8 showing insert size 858 bp (B) Confirmation of recombinant clones using restriction digestion in pRSET-B M; marker, Lane 1; Control pRSET-B correlating 2.8 kb, Lane 2 to 6 showing recombinant clones with insert size 969 bp.
Additional file 4: Figure S2. Construction of knockout. (A) allelic exchange substrate construction M: 1 kb Marker, Lane 1: SacB digested with van911 showing four fragments 3.6 kb, 1.6 kb, 979 bp, 567 bp, Lane 2& 3: 0148 AES construct showing four fragments 3.6 kb, 1.6 kb, 837 RA, 605LA. (B) Packaging of AES with phAE159. M: 1 kb ladder, Lane 1: phAE159 digested with pac-I showing 50 kb and insert 3.8 kb, Lane 2& 3 clones without insert, Lane 4& 5 clone digested with pac-I showing phAE159 50 kb and 6.6 kb AES (C) Confirmation of knockout using PCR M: 1 kb ladder, Lane 1: Rv DNA amplified with right arm, Lane 2,4,5: knockout DNA amplified with hyg Forward primer and right arm reverse primer, Lane 6: Rv DNA not showed amplification with hyg & reverse primer.
About this article
Cite this article
Bhargavi, G., Hassan, S., Balaji, S. et al. Protein–protein interaction of Rv0148 with Htdy and its predicted role towards drug resistance in Mycobacterium tuberculosis. BMC Microbiol 20, 93 (2020). https://doi.org/10.1186/s12866-020-01763-1
- Mycobacterium tuberculosis
- Drug resistance
- Stress response