key: cord-0966161-k2fc4gxh
authors: Sharma, Aashish; Kumar, Arbind; Rashid, Mudasir; Amnekar, Ramchandra Vijay; Gupta, Sanjay; Kaur, Jagdeep
title: A Phagosomally Expressed Gene, rv0428c, of Mycobacterium tuberculosis Demonstrates Acetyl Transferase Activity and Plays a Protective Role Under Stress Conditions
date: 2022-02-17
journal: Protein J
DOI: 10.1007/s10930-022-10044-x
sha: 08dd15bb805f7fd9e13063bf6e513cd7b1ed8116
doc_id: 966161
cord_uid: k2fc4gxh

Mycobacterium tuberculosis genome is composed of several hypothetical gene products that need to be characterized for understanding the physiology of bacteria. Rv0428c was one of the 11 proteins exclusively identified within the phagosomal compartment of macrophages infected with mycobacteria and marked as hypothetical. The expression of rv0428c gene was upregulated under acidic and nutritive stress conditions in M. tuberculosis H37Ra, which was supported by potential sigma factor binding sites in the region upstream to the rv0428c gene. The bioinformatics analysis predicted it to be a GCN5- acetyl transferase, belonging to the Histone acetyl transferase (HAT) family. The docking analysis predicted formation of hydrogen bonds and hydrophobic interactions between donor acetyl-co-A and histone H3 tail region. rv0428c gene was cloned and expressed in E. coli. The protein was purified to homogeneity and was fairly stable over a wide range of pH 5.0–9.0 and temperature up to 40 °C. The HAT activity of purified Rv0428c was confirmed by in vitro acetylation assay using recombinant H3 histone expressed in bacteria as substrate, which increased in time dependent manner. The results suggested that it is the second confirmed acetyl transferase in M. tuberculosis H37Rv. Furthermore, rv0428c was over expressed in surrogate host M. smegmatis, which led to enhanced growth rate and altered colony morphology. The expression of rv0428c in M. smegmatis promoted the survival of bacteria under acidic and nutritive stress conditions. In conclusion, Rv0428c, a phagosomal acetyl transferase of M. tuberculosis, might be involved in survival under stress conditions. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10930-022-10044-x.

Tuberculosis (TB) is amongst the most infectious diseases that has inflicted mankind. Every third individual is infected with Mycobacterium tuberculosis. According to WHO, an estimated 1.4 million deaths were reported from tuberculosis worldwide with nearly 0.2 million individuals being Human Immunodeficiency Virus (HIV)-positive [1] . There has been no respite from the woes of M. tuberculosis due to its poorly understood biology, survival strategies, life cycle and pathogenesis [2, 3] . The genome of M. tuberculosis was sequenced way back in 1998, identifying approximately 4000 proteins with 40% of these annotated as hypothetical proteins [4] . For better understanding of the physiology and virulence of this bacterium, there is an urgent need for assigning specific roles to these hypothetical proteins by detailed characterization.

M. tuberculosis is an intracellular pathogen with a dynamic proteome, which plays a role in survival of the pathogen under stress conditions encountered within the host macrophages [5] . The immense success of this bacterium as a pathogen is dependent upon it's capability to utilize the host macrophages for survival and proliferation by employing various strategies [6] . These include prevention of the lysosome phagosome fusion and acidification of the phagosome, protection from reactive oxygen radicals and altering the immune response [7] . Intraphagosomal microbes were reported to alter host cell physiology. In some cases it can induce apoptosis, while macrophages containing mycobacterium in phagosomes are known for long term survival. For its intraphagosomal survival, mycobacterium has to deal with various intracellular stress conditions such as hypoxia, acidic, oxidative, nutritive and iron stress [8] .

The 2-dimesional gel electrophoresis and Mass spectrometry analysis of intraphagosomally grown mycobacteria in bone marrow macrophages, identified 11 exclusively present mycobacterial proteins [9] . These proteins needed special attention to pin point their role in intracelluar survival of M. tuberculosis. Six of the total 11 identified phagosomal mycobacterial proteins, i.e., Rv2691, Rv1627c, Rv1191, Rv1130, Rv0489 and Rv0428c, were unique as these were not noticed in in vitro growth of mycobacterium. Out of the above mentioned six proteins, Rv1191, Rv1130 and Rv0428c were reported as conserved hypothetical proteins. As Rv0428c has not been characterized yet, it was selected for detail characterization.

The bioinformatics analysis predicted Rv0428c to be a GCN5-acetyl transferase belonging to the Histone acetyl transferase family (HAT). HAT proteins are involved in acetylation of core histones, which further results in important regulatory effects on chromatin structure, assembly and gene transcription. Till date only one mycobacterium protein, eis (enhanced intracellular protein), has been shown to possess N ε -acetyltransferase activity [10] .

Rv0428c was reported to be specifically present in pathogenic strains like M. tuberculosis, M. bovis and clinical strain CDC1551, but is conspicuously absent in non-pathogenic M. smegmatis strain pointing towards some critical role played by this gene in pathogenesis/virulence. In the present investigation, attempt has been made to characterize Rv0428c by biochemical and biophysical methods. The expression of Rv0428c was studied under normal and stress conditions. The effect of expression of rv0428c on colony morphology and growth pattern was monitored in surrogate host M. smegmatis which lacks rv0428c gene sequence in its genome. The experiments were also carried out to investigate the role of Rv0428c in conferring drug resistance.

The bacterial strains used in the study were procured from IMTECH, Chandigarh. The strains which were used included E. coli DH5α, E. coli BL21 (DE3) and M. tuberculosis H37Ra. The E. coli was grown in 2% Luria-Bertani (LB) broth and Mycobacterium H37Ra was cultured in 0.5% Middlebrook 7H9 (7H9) media supplemented with 1% glycerol and 0.05% Tween-80. 

The promoter analysis of rv0428c gene was done by analyzing the 250 bp nucleotide sequence upstream to the operonic arrangement of the rv0428c. The sigma factor binding sites were identified and marked within this upstream DNA sequence.

Multiple sequence alignment was performed by Clustal Ω, to study the extent of conservation and variation in protein sequences [11] . Espript 3 was used to assign the secondary structure to the alignment file [12] . The genomic organization of rv0428c gene of M. tuberculosis and its orthologs in M. bovis and M. leprae was checked by analyzing the gene sequences upstream and downstream to these genes in their respective genomes. The protein-protein interaction of Rv0428c protein was done by using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database which predicts the possible interaction of the protein with other proteins that uses numerous sources, including experimental data, computational prediction methods and public text collections [13] .

The blast protein hits showed BLASTP analysis of Rv0428c protein was performed against PDB proteins for selection of appropriate templates for generation of 3D structure models. The blast protein hits demonstrated that the identity and query coverage were found to be below 35%. Therefore, ab-initio based approach was utilised for protein modelling using I-TASSER program. Visualization of the 3D model structure of protein was done by using PyMol software [14, 15] .

The analysis of binding pocket of Mtb Rv0428c was carried out by Computed Atlas of Surface Topography of proteins (CASTp) server [16] . This server was endowed with weighted Delaunay triangulation and the alpha complex framework which were used for measuring the shapes of molecules. It unveils inaccessible cavities and solvent accessible surface geometry in protein. The structural pockets and cavities were calculated in relation to area and volume by two approaches; solvent accessibl surface model (Richards' surface) and molecular surface model (Connolly's surface).

The ligand binding site of Rv0428c protein was predicted using COACH program [17] . As predicted, Gly183, Gly186, Ser215, Met137, Ala181, Arg187, Trp190, Thr212 and Val216 from Rv0428c were involved in binding with acetylco-A, whereas, Gly128, Val151, Trp122, Leu125, Ala152 and Arg155 of Rv0428c were found to be involved in binding with histone H3 tail region. A grid was generated in the region of these predicted binding site residues of prepared protein by means of the Autodock tools [18] . For docking procedures, Autodock Vina was used at its default parameters. Through Autodock tools, protein and ligands were generated in pdbqt formats. The gasteiger charge and polar hydrogen were assigned to the receptor and ligands. The grid boxes were set in the range of approximately 30 × 30 × 30 point size, spaced at 0.375 Å in each direction around the cavity for accommodating the ligands.

M. tuberculosis H37Rv chromosomal DNA was a kind gift from Dr. U. D. Gupta, JALMA, Agra. The 18 to 24 bases were selected from the terminal region of gene sequences and analysed by integrated oligoanalyzer tools for optimal designing of primers (Integrated DNA Technologies, Inc., U.S.A., http:// eu. idtdna. com/ site). Optimization of Tm, GC content and hairpin loop structure of primers was carried out [19] . The primers used in the study are specified in Table 1 soluble fraction and was purified to homogeneity by Nickelnitrilotriacetic acid (Ni-NTA) chromatography. The integrity and purity of protein was analyzed on 12% SDS-PAGE.

Biophysical characterization was carried out by performing CD spectroscopy and fluorescent spectroscopy. CD spectroscopy analysis was carried out by measuring Far UV-CD spectra of rRv0428c protein with a JASCO J-815 spectrofluorometer. Far UV-CD spectra were collected over the wavelength range of 195-250 nm at 25 °C. To study the effect of temperature, 200 µg of protein was incubated at temperatures ranging from 30 to 70 °C for 30 min. The spectra were collected over a range of 200-240 nm wavelengths with a band width of 1 nm and response time of two seconds. The CD values were expressed as molar ellipticity ([θ], degcm 2 dmol −1 ). For studying the effect of pH, 200 µg protein was incubated with 50 mM phosphate buffer of varying pH from 5.0 to 11.0 for 30 min and the spectra were collected from 200 to 240 nm wavelengths. The fluorescence spectrum of rRv0428c protein was measured using JASCO J-815 spectrofluorometer using a 10 mm path length quartz cuvette. For determining the intrinsic tryptophan fluorescence, the excitation wavelength used was 295 nm with emission measured in the range of 310-400 nm. The spectra were recorded from 20 to 90 °C by incubating the protein for 30 min.

In vitro acetylation assay for rRv0428c was performed by using bacterially expressed recombinant histone H3 as substrate and acetyl-coA as a donor molecule and mammalian core histones as a positive control. The reaction was set up in acetyl-transferase assay buffer (50 mM Tris-Cl, pH 8, 10% glycerol, 10 mM butyric acid, 0.1 mM EDTA, 1 mM DTT, 1 mM PMSF) with 10 µM acetyl-coA, 10 µg histone H3, with rRv0428c. A time dependent kinetics was studied by incubating the reaction mixture for 1 h, 2 h and 3 h, in a 30 °C water bath followed by addition of 2X SDS-PAGE sample buffers and 10 min boiling for stopping the reaction. The prepared samples were then subjected to 18% SDS-PAGE gel electrophoresis and transferred onto PVDF membranes. The membrane was stained with 0.05% Fast green for validating equal loading of samples followed by western blotting with Anti H3K9ac antibody (1:5000, Millipore).

rv0428c was amplified from the recombinant pET28a plasmid using primers mentioned in Table 1 . Amplified gene product was separated on a 1.2% agarose gel and the fragments were subsequently excised and eluted. Amplified rv0428c and pVV16 plasmid were digested with BamHI restriction enzyme followed by purification of cut fragments. The digested products were ligated using T4 DNA ligase. M. smegmatis mc 2 155 was grown in 7H9 broth supplemented with 1% glycerol, 1% OADC and 0.2% Tween 80 at 37 °C, 180 rpm till the OD 600 reached 0.6-0.8. The cells were harvested by centrifugation at 4000 rpm at 4 °C. Cells were incubated on ice for 2 h. Then, cells were washed three times with 10% glycerol by centrifugation at 4000 rpm for 10 min. Cells were re-suspended in 10% glycerol at 1/500th of original volume. The recombinant plasmid and pVV16 alone (5 ng) were added to the 100 µl electrocompetent M. smegmatis mc 2 155 cells and transferred to 0.2 cm Biorad electroporator cuvette. The cells were then incubated on ice for 10 min. The electroporation was carried out at 25 µF and 2.5 kV for 5 ms using Gene pulser, BioRad, USA followed by addition of 7H9 media immediately. The electroporated cells were then incubated at 37 °C, 180 rpm for 6 h. Finally, the cells were spread onto 7H10 agar plates containing Kan + . M. smegmatis harbouring pVV16-rv0428c or pVV16 alone were inoculated in 10 ml M7H9 media supplemented with 0.1% glycerol, 0.2% Tween 80 and 1% OADC at 180 rpm, 37 °C. Next day, the cultures were sub-cultured in 100 ml flask containing M7H9 media supplemented with 0.1% glycerol, 0.05% Tween 80 and 1% OADC after normalization of A 600 . The cultures were further grown at 37 °C and 180 rpm. The growth pattern was measured by taking absorbance at 600 nm and CFU counting at regular time intervals. To study the differences in colony morphology of M. smegmatis harbouring pVV16-rv0428c or pVV16 alone, both the strains were plated onto M7H10 plates supplemented with 0.1% glycerol, 0.2% Tween 80 and 1% OADC and incubated at 37 °C for 5-6 days. Microphotography of the colonies was carried out.

Conditions In Vitro

The cultures Msmeg-pVV16-rv0428c (test) and Msmeg-pVV16 (control) were grown in M7H9 liquid media as described previously in 10 ml tubes. Next day, sub-culturing of both the cultures was carried out in 20 ml M7H9 media (pH 7.2), M7H9 media (pH 6.0) and M7H9 media (pH 5.0) flasks containing 1% OADC to make the final absorbance 0.5 at 600 nm. The survival was checked by plating suitable dilutions of the cultures onto M7H10 agar plates having kanamycin at 37 °C for 2-3 days.

For induction of nutritive stress conditions, both the cultures Msmeg-pVV16-rv0428c (test) and Msmeg-pVV16 (control) were grown overnight in 10 ml M7H9 tubes at 180 rpm at 37 °C. Further sub-culturing was carried out in 20 ml flasks containing 1X PBS. The final absorbance was normalized to 0.5 at 600 nm. The cultures were then incubated at 180 rpm at 37 °C for 12 h and 24 h followed by spreading of appropriate dilutions onto M7H10 plates containing kanamycin and incubating the plates for 2-3 days at 37 °C.

For oxidative stress administration, the bacterial pellet was re-suspended in M7H9 media containing 1% glycerol, 0.2% Tween 80, 1% OADC and 5 mM H 2 O 2 . The culture was then kept at 37 °C without shaking. The culture without H 2 O 2 served as control. Both the cultures were harvested after 6 h. RNA was isolated from the bacterial cultures and used for relative expression analysis.

Cells were washed thrice with PBS (pH 7.2) followed by two washings with iron deficient media to remove any traces of culture media. The washed cells were further re-suspended in iron deficient media. The control culture was supplemented with 160 µM of FeCl 3 . The cells were grown for 96 h and harvested for RNA isolation.

Resazurin redox indicator test was used for checking the effect of Rv0428c protein on drug susceptibility (Palomino et al., 2002) . The susceptibility of Msmeg-rv0428c and Msmeg-pVV16 was determined against various drugs-streptomycin, chloramphenicol, isoniazid and rifampicin. The bacterial cultures were grown to mid-log phase and diluted for equal dispensation of the cells (4 × 10 5 ) in 48-well plates.

The reason for using mid-log phase bacterial cells was that these exhibit constant growth rate and have an optimal level of expression. The drugs were diluted to working concentration ranges in M7H9 media without the detergent Tween-80. Different concentrations of drugs was added and incubated at 37 °C for 2 h. A working 1:1 dilution of 10X stock resazurin was prepared in 20% tween-80 and 8 µl of it was added per well in 48-well plates. The viable bacteria lead to conversion of resazurin to resurfin which was monitored by the change in colour of the medium from blue to pink. The survival of Msmeg-pVV16-rv0428c and Msmeg-pVV16 was also monitored by counting the CFU/mL after treatment with chloramphenicol.

rv0428c protein was detected exclusively in the phagosomal compartment of macrophages infected with M. tuberculosis, making it an ideal candidate for being involved in playing a regulatory role during stress encounter. We, thereby, proceeded for analyzing the expression of rv0428c in various stress conditions including acidic, oxidative, nutritive and iron stress in-vitro in M. tuberculosis H37Ra. Semi quantitative PCR analysis demonstrated upregulation of rv0428c in acidic and nutritive stress by 5.4 and 3.6 fold, respectively (Fig. 1A, B) . The analysis of upstream regulatory promoter region showed the presence of two putative sigma factors binding sites i.e. sigE (CGA CAT (15) (16) (17) (18) (19) GGTTC) and sigF (GGC GAA (16) (17) (18) (19) (20) SGTTS), which were previously implicated in acidic and nutritive stress conditions [22, 23] (Fig. 1C ).

The gene rv0428c and protein sequence was retrieved from TubercuList database. The size of the gene is 909 bp and it encodes a protein of 302 amino acids. It is predicted to be involved in intermediary metabolism and contains GNAT domain at C-terminal [24] . Rv0428c has orthologs in three other pathogenic mycobacterium species M. bovis, M. leprae and M. canettii. GCN5-related N-acetyltransferase from Kribbella flavida with PDB_ID 4IUS showed maximum identity with Rv0428c protein and was used as a template for performing multiple sequence alignment. Multiple sequence alignment revealed the conserved regions between the Rv0428c protein and its counterparts in other Mycobacterium species (Sup Fig. 1 ). In prokaryotes, the genomic organization of a gene is often considered in speculating on the probable gene function based on its positional counterpart or gene organization. The gene encoding for Rv0428c in M. tuberculosis is flanked by probable polypeptide deformylase (def) and exodeoxyribonuclease III (xthA) in the genomic organization. This genomic organization is also conserved in M. bovis and M. leprae (Fig. 2A) . The Rv0428c protein possessed the VAPTHRRRG sequence similar to the V/I-X-X-X-X-Q/R-X-X-G consensus sequence of GCN5acetyl transferases suggesting that Rv0428c is a probable member of the GCN5-acetyl transferase family. The interaction studies revealed that Rv0428c protein interacted with several proteins including the eis (enhanced intracellular survival) protein of M. tuberculosis which is an N ε -acetyl transferase (Sup Fig. 1 ).

The model of Rv0428c protein was generated by using I-TASSER. A set of 5 models were constructed based on the 10 best templates. The models generated were sorted based on their C-scores, which represent the confidence in the predicted structure on the basis of threading template alignments and convergence parameters involved in the simulation [17] . The range for C-score lied between − 5 and 2 with higher C-score signifying a model with higher confidence. The C-score for Rv0428c model lied in the permissible range and was found to be − 0.40. Depending upon the C-score, TM-score and RMSD values were calculated. The model had a TM-score of 0.66 ± 0.13 and RMSD value of 6.6 ± 4.0 Å. The final PROCHECK program measured the Ramchandran plot statistics, which showed that around 80% of the residues are in the most favoured region. Quality of the generated models was assessed by Verify 3D and QMEAN score 4. Verfy3D result demonstrated that the 80% of the amino acids are in the acceptable range (≥ 0.2 in the 3D/1D profile). QMEAN score is below 0.5 which showed the accuracy of predicted model for further experimental use (Supp. Figure 2) . These statistics confirmed that our predicted model had precise topology. The alignment analysis of Rv0428c protein with the template GCN5-related N-acetyltransferase from Kribbella flavida (PDB_ID: 4IUS) was done (Fig. 2B) . The superimposition of 3D models of Rv0428c protein and template GCN5-related N-acetyltransferase from Kribbella flavida (PDB_ID: 4IUS) revealed overlapping of the α-helices and β-strands (Fig. 2C ).

Acetyl coenzyme-A (AcoA) acted as a donor of acetyl group for conversion of conserved lysine amino acid residues on histone proteins to be acetylated to ε-N-acetyl lysine and hence regulating the gene expression. As Rv0428c is a probable GCN5-acetyl transferase, we docked the acetyl-co-A molecule in the binding pocket of Rv0428c protein. The size of the cavity of Rv042c protein was predicted by using CASTp algorithm. Histone H3 was used as substrate because most of the histone acetylation events take place on H3, for interacting with the cavity as the substrate for GCN5-acetyl transferase is histone H3. The binding pocket size for acetylco-A and histone H3 tail were calculated. The binding pocket volume for acetyl-co-A was observed to be 30.9 Å and area was 94.6 Å, whereas, in case of histone H3 tail region, the volume was 34.9 Å with an area of 88.6 Å ( Table 2 ). The binding energy was found to be − 6.1 and − 5.8 kcal/mol for acetyl-co-A and histone H3 tail, respectively.

Acetyl-co-A has a molecular weight of 809.57 g/mol. The complex of Rv0428c-acetyl-co-A was stabilized by the formation of hydrogen bonds and hydrophobic interactions. The acetyl-co-A was docked against Rv0428c protein (Fig. 3A) .

The H-bonds were formed between Gly185, Asp223 and Gly224 amino acid residues of Rv0428c and acetyl-co-A. The distances between the donor and acceptor molecules was computed on the basis of maximum acceptor (A) and H-bond donor (D) and were found to be 3.07 Å, 2.92 Å and 3.02 Å, respectively. Several hydrophobic interactions were also observed which included Met137, Ala181, Arg192, Asp223, Ala225 and Thr272 (Fig. 3B ). The histone H3 tail was 14 amino acids long with molecular weight of 1.45 kDa. The formation of hydrogen bonds and hydrophobic interactions stabilized the protein-protein interactions between Rv0428c protein and H3 tail region. The H-bonds were formed between Leu183 and Gln188 amino acids of Rv0428c protein and H3 tail. The distance between the donor and acceptor were found to be 2.99 Å and 2.96 Å, respectively. The hydrophobic interactions between Trp182, Ala212, and Arg215 were observed (Fig. 3C, D; Table 3 ). The docking of acetyl-co-A and H3 tail results signified that these molecules interact with our protein of interest, Rv0428c and might help in acetylation of histone proteins further leading to regulation of gene expression.

For expression and purification studies, pET Rv0428c plasmid was transformed in E. coli BL21 DE3 cells. Rv0428c protein was expressed in soluble fraction and was purified by Ni-NTA chromatography to homogeneity as demonstrated by a single band on SDS-PAGE with an approximate molecular mass of 35 kDa (Fig. 4A ). 

Far UV CD spectrum (190-260 nm) of rRv0428c protein was recorded using the J-815 CD spectropolarimeter at room temperature. The spectrum revealed the characteristic negative ellipticity comprising of both α-helix and β-sheets in the secondary structure (Fig. 4B) . The secondary structure of protein was determined by spectra analysis software.

Relative amount of structural element estimated for the rRv0428c protein were 12.6% α-helices, 36% β-sheets, 17.5% turn and 33.9% random coil. This data confirmed that the purified protein is in properly folded state.

The unfolding property of rRv0428c protein was determined by incubating it at temperatures ranging from 30 to 70 °C. The gradual loss of molar ellipticity combined with shift of temperature D Far UV CD spectra of rRV0428c protein after incubating the protein for 1 h with different buffers of pH-5.0-9.0 demonstrated that Rv0428c protein was stable over a wide range of pH E Fluorescence spectra of purified rRv0428c protein at different temperatures (20 to 90 °C) recorded from 310 to 400 nm wavelength demonstrated gradual decrease in the intrinsic fluorescence with the subsequent increase in temperature of minima was used to monitor the state of protein through CD spectroscopic analysis. There was a gradual decrease in negative molar ellipticity of rRv0428c protein with variation in temperature. The Rv0428c protein conformation was stable up to temperatures ≤ 40 °C with subsequent convergence of curves. With increase in temperature beyond 40 °C the protein starting losing its secondary structure (Fig. 4C) . The protein exhibited stability over a wide pH range (5.0-9.0), as was evident from gradual change in the negative ellipticity (Fig. 4D ).

To study the effect of temperature on the tertiary structure of Rv0428c protein, intrinsic tryptophan fluorescence spectroscopy was done. The Rv0428c protein has 11 tryptophan amino acid residues. The maximum fluorescence intensity was observed at 340 nm emission wavelength. There was gradual decrease in the intrinsic fluorescence with the subsequent increase in temperature. The peak maxima shifted from 340 to 344 nm by increase in temperature from 60 to 90 °C indicating red shift (Fig. 4E) . These results pointed towards the fact that intrinsic fluorescence decreased with subsequent increase in temperature.

As bioinformatics analysis predicted that Rv0428c belongs to GNAT family of histone acetyl transferases (HATs), we proceeded for checking the histone acetyl transferase activity of rRv0428c. rRv0428c was purified to homogeneity and in vitro acetylation assay was performed using bacterially purified H3 histone as substrate. The H3 extracted from eukaryotic AGS cell line was used as positive control (lane 8) because histones from eukaryotic system will have endogenous acetylation. The time dependent in vitro acetylation was also performed. With increase in time a sequential increase in the amount of acetylated histone H3 was observed, proving the acetyltransferases activity of rRv0428c (Fig. 5) . The lanes 5-7 demonstrated HAT activity of NC fraction of AGS cell line which was found to be less compared to rRv0428c. Overall, these results signify that rRv0428c protein possessed significant HAT activity, which increased upon subsequent increase in incubation period.

The colony morphology of M. smegmatis harboring rv0428c (Msmeg-pVV16-rv0428c) and pVV16 (Msmeg-pVV16) cultures were compared. Both the cultures were spread at very low cell density onto M7H10-Kan + plates. A significant change in the colony morphology of Msmeg-pVV16-rv0428c was observed as compared to Msmeg-pVV16 (Fig. 6A) . The Mmeg-pVV16 colony was rough with a bulge in the centre, whereas, the Msmeg-pVV16-rv0428c colony was smooth, wet and flattened. The growth pattern of M. smegmatis containing rv0428c was compared with that of the vector pVV16 alone at different time points-24 h, 48 h, 72 h and 96 h by CFU counting. Enhanced growth of the M. smegmatis harbouring rv0428c gene was observed as compared to the pVV16 vector alone. The Msmeg-pVV16-rv0428c in comparison to Msmeg-pVV16 displayed nearly 1.7 fold enhanced growth after 48 h and 2.8 fold increase in growth after 72 h. The difference in growth was 1.4 fold after 96 h (Fig. 6B ). Both the recombinant M. smegmatis culures were also exposed to nutrient stress conditions by incubating the cultures in 1X PBS instead of the M7H9 growth media. There was a 2.5 fold increase in survival of the test culture in comparison to control in the presence of nutrient deprived conditions (Fig. 7B ).

Drug susceptibility of Msmeg-pVV16-rv0428c and Msmeg-pVV16 cultures was checked using resazurin as the indicator of viability. The widely used anti-TB drugs-streptomycin, chloramphenicol, isoniazid and rifampicin were used for DST. The Msmeg-pVV16-rv0428c culture was able to grow in the presence of chloramphenicol with minimum inhibitory concentration (MIC) of 3 μg/ml. This is evident from the colour change from blue to pink in the 48-well plate (Fig. 8A) . The CFU/ml were also monitored which depicted enhanced survival of Msmeg-pVV16-rv0428c in comparison to the Msmeg-pVV16 culture upon incubation with increasing concentration of chloramphenicol (Fig. 8B ). The plates with isoniazid, streptomycin and rifampicin showed no colour change indicating that the M. smegmatis cultures were susceptible to these drugs.

M. tuberculosis is unique in its ability to survive for prolonged period within the harsh environment inside the host, specifically within the phagosomal compartment of the host macrophages by inhibiting the maturation of phagosome [25, 26] . Out of 11 mycobacterial proteins exclusively identified in phagosome, 3 proteins, Rv0428c, Rv1130 and Rv1191 were annotated as hypothetical proteins [9] . Since the identification of hypothetical proteins by Cole et al. in 1998 , attempts were being made to characterize these hypothetical proteins. Rv0428c protein has been predicted as an ideal candidate for playing a protective role in the survival of mycobacteria inside the adverse host cell environment.

The adaptive responses of bacterial species to various environmental stress conditions have been shown to be involving alternative sigma factors [27] . The up regulated expression of rv0428c under acidic and nutrient stress conditions and the potential sigma factor binding sites for sigE and sigF in the nucleotide sequence upstream to the operon suggested their role in the transcriptional regulation under stress conditions. Several mycobacterium genes including PknG, rv3097c and rv1169c, were over expressed under the influence of acidic and nutritive stress, aiding in the survival of M. tuberculosis [22, 23, 28] . Rv0428c of M. tuberculosis was predicted to be a probable acetyl transferase belonging to the GNAT (GCN5related acetyl transferase) family of HATs (Histone Acetyl Transferase). ε-amino lysine acetylation was not limited only to histone modification and regulation of transcription, but was reported to be involved in several cellular processes [29] . This might alter the charge on protein, its conformation and protein stability. Acetylomes or genome-wide identification of lysine-acetylated proteins in bacteria has pointed towards wide diversity of functions for lysine acetylated proteins. The GCN5-related N-acetyltransferase (GNAT) super family is a large group of evolutionarily related acetyl transferases, with multiple paralogs in organisms from all kingdoms of life [24, 30] . Although the functional role of protein acetylation in eukaryotes has long been studied, it was recently discovered that acetylation of proteins is common in bacteria as well [31] [32] [33] [34] . The GNAT family of transferases have been shown to exhibit sequence homology with a class of eukaryotic transcription factors, the first of which was the yeast GCN5 [35] . As Rv0428c has been predicted to be a member of the GNAT family, we proceeded for docking of Rv0428c with eukaryotic histone H3 tail region and acetyl-co-A to establish the interactions between the substrate histone and donor acetyl-co-A with our protein of interest. The docking analysis predicted multiple hydrogen bonds and hydrophobic interactions between the Rv0428c protein-acetyl-co-A complex and Rv0428c-H3 tail region. Previously, GNAT proteins have been implicated in acetylation of lysine in the core histone H3 tail region [36, 37] . These findings provide substance to the fact that Rv0428c is a probable member of the GNAT family, which might be playing an important role in acetylation of proteins.

The thermal unfolding studies of recombinant Rv0428c using CD spectroscopy revealed that the secondary structure of rRv0428c protein was completely stable up to 40 °C. Previously, multiple tryptophan residues have been implicated in providing stability to the M. tuberculosis protein Rv0774c [38] . The Rv0428c protein have 11 tryptophan amino acid residues, however, the conformational changes in the tertiary structure of Rv0428c protein revealed that it was stable up to 60 °C with peak maxima at 340 nm. After 60 °C, red shift is observed leading to shift of peak maxima towards the higher wavelength. An in vitro acetylation assay performed using recombinant histone H3, confirmed the acetyl transferase activity of the protein that increased in time dependent manner. Rv1988, a secretory mycobacterial methyl transferase was localized in the host nucleus and interacted with histone H3 resulting in repression of genes responsible for first line of defence against tuberculosis. The deletion of Rv1988 suppresses the survival of M. tuberculosis inside the host [39] . The protein-protein interaction study revealed that Rv0428c showed significant interactions with Rv2416c (eis). Previously, Rv2416c was designated as enhanced intracellular survival (eis), owing to the fact that it was responsible for enhanced survival of M. smegmatis inside the macrophage cell line [40] . Pink colour indicates live bacteria, while blue indicates dead bacteria B Survival of Msmeg-pVV16-rv0428c and Msmeg-pVV16 was monitored by counting the CFU/mL after treatment with chloramphenicol. Results were expressed in % survival (CFU counts without drug treatment was considered to be the 100% survival). Data are representative of three independent biological replicates and shown as mean ± SD. Statistical analysis was assessed using student's t-test (*p ≤ 0.05 and **p ≤ 0.01) No homolog of rv0428c was found in M. smegmatis making it an ideal host candidate for performing in vitro experiments which mimic the mycobacterium species. The over expression of rv0428c altered the colony morphology besides increasing the growth rate of M. smegmatis. The association of colony morphology with virulence of mycobacterium species is a well established fact. There are previous reports which lend support to the fact that expression of several mycobacterial genes like rv1169c [41] , rv1818c [42] and rv0774c [43] altered the colony morphology and growth rate of M. smegmatis. The advent of drug-resistant strains of M. tuberculosis has led to drug susceptibility testing of individuals complaining of symptoms of tuberculosis. As this makes sure that a particular individual will respond to the prescribed anti-TB drug regimen. The GNATs have been demonstrated to be involved in acetylation of aminoglycoside antibiotics by acting as aminoglycoside modifying enzymes leading to resistance against these antibiotics [44, 45] . revealed that Rv0428c is resistant to chloramphenicol, implicating the role of Rv0428c in conferring drug resistance to Mycobacterium species. The resistance to chloramphenicol could be due to the sequence similarity of Rv0428c with the chloramphenicol acetyl transferase, an enzyme which detoxifies chloramphenicol by acetylation and is responsible for chloramphenicol resistance in bacteria [46] . Further study is required to find out the specific target of this enzyme in M. tuberculosis by making the knockout of the gene followed by acetylome analysis.

In summary, Rv0428c is the active acetyl transferase belonging to GNAT family of HATs identified in M. tuberculosis. It was exclusive to the intraphagosomal compartment of infected macrophages. Rv0428c was stable over a wide range of pH and could retain its tertiary structure upto 60 °C pointing towards its probable role in harsh conditions. It also plays a protective role under nutrient and acidic stress conditions in vitro. The expression of protein resulted in altered colony morphology and enhanced growth of M. smegmatis under various stress conditions. The role played by Rv0428c in survival of M. tuberculosis under the stress conditions makes it a probable candidate for drug targeting.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s10930-022-10044-x.

Population-based resistance of Mycobacterium tuberculosis isolates to pyrazinamide and fluoroquinolones: results from a multicountry surveillance project

Latent tuberculosis: interaction of virulence factors in Mycobacterium tuberculosis

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence

Differential gene expression in response to exposure to antimycobacterial agents and other stress conditions among seven mycobacterium tuberculosis whiB-like genes

Making the most of the host; targeting the autophagy pathway facilitates Staphylococcus aureus intracellular survival in neutrophils

Mycobacterial survival strategies in the phagosome: defense against host stresses

The response of mycobacterium tuberculosis to reactive oxygen and nitrogen species

Proteins unique to intraphagosomally grown Mycobacterium tuberculosis

Mycobacterium tuberculosis Eis protein initiates suppression of host immune responses by acetylation of DUSP16/MKP-7

Clustal Omega, accurate alignment of very large numbers of sequences

Deciphering key features in protein structures with the new ENDscript server

STRING v10: protein-protein interaction networks, integrated over the tree of life

SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information

CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues

I-TASSER: fully automated protein structure prediction in CASP8

AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility

Primer based approach for PCR amplification of high GC content gene: mycobacterium gene as a model

Strategies for optimization of heterologous protein expression in E. coli: roadblocks and reinforcements

Expression and purification of functional epitope of pigment epithelium-derived factor in E. coli with inhibiting effect on endothelial cells

A novel lipase belonging to the hormone-sensitive lipase family induced under starvation to utilize stored triacylglycerol in Mycobacterium tuberculosis

Structure and functions of the GNAT superfamily of acetyltransferases

M tuberculosis PknG manipulates host autophagy flux to promote pathogen intracellular survival

From infection niche to therapeutic target: the intracellular lifestyle of Mycobacterium tuberculosis

Identification and characterization of lipase activity and immunogenicity of lipl from mycobacterium tuberculosis

Evaluation of a nutrient starvation model of Mycobacterium tuberculosis persistence by gene and protein expression profiling

Lysine acetylation targets protein complexes and co-regulates major cellular functions

GCN5-related N-acetyltransferases: a structural overview

Bacterial protein acetylation: the dawning of a new age

The diversity of lysineacetylated proteins in Escherichia coli

Lysine acetylation is a highly abundant and evolutionarily conserved modification in Escherichia Coli

Rv0802c is an acyltransferase that succinylates and acetylates Mycobacterium tuberculosis nucleoid-associated protein HU

Genetic isolation of ADA2: a potential transcriptional adaptor required for function of certain acidic activation domains

Histone H3 specific acetyltransferases are essential for cell cycle progression

Transcription-linked acetylation by Gcn5p of histones H3 and H4 at specific lysines

Rv0774c, an iron stress inducible, extracellular esterase is involved in immune-suppression associated with altered cytokine and TLR2 expression

Mycobacteria modulate host epigenetic machinery by Rv1988 methylation of a non-tail arginine of histone H3

Eis (Enhanced Intracellular Survival) protein of Mycobacterium tuberculosis disturbs the cross regulation of T-cells

PE11, a PE/PPE family protein of Mycobacterium tuberculosis is involved in cell wall remodeling and virulence

Rv1818c-encoded PE_ PGRS protein of Mycobacterium tuberculosis is surface exposed and influences bacterial cell structure

Strategies for optimization of heterologous protein expression in E. coli : roadblocks and reinforcements

Aminoglycosides modified by resistance enzymes display diminished binding to the bacterial ribosomal aminoacyl-tRNA site

Overexpression of the chromosomally encoded aminoglycoside acetyltransferase eis confers kanamycin resistance in Mycobacterium tuberculosis

Recombinant genomes which express chloramphenicol acetyltransferase in mammalian cells

Acknowledgements The authors duly acknowledge Department of Science and Technology (DST) for the financial assistance to Dr Jagdeep Kaur and University Grants Commission (UGC) for providing fellowship to Aashish Sharma.

Author Contributions JK conceived the idea, designed the study and supervised the research work. AS performed most of the experiments and wrote the manuscript. AK performed the Bioinformatics work. Histone acetylation experiments were carried out by MR and RVA under the supervision of SG at ACTREC, Mumbai.Funding This work was supported by the Department of Science and Technology (DST), India.

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

The authors declare that they have no competing interests.Ethical Approval Not applicable.

Not applicable.