ORIGINAL RESEARCH article

Front. Immunol., 26 April 2019

Sec. Vaccines and Molecular Therapeutics

Volume 10 - 2019 | https://doi.org/10.3389/fimmu.2019.00716

Exploring the Potential Role of Moonlighting Function of the Surface-Associated Proteins From Mycobacterium bovis BCG Moreau and Pasteur by Comparative Proteomic

  • 1. Laboratório de Genômica Funcional e Bioinformática, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil

  • 2. Laboratório Interdisciplinar de Pesquisas Médicas, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz—FIOCRUZ, Rio de Janeiro, Brazil

  • 3. Unidade de Espectrometria de Massas e Proteômica, Instituto de Bioquímica Médica Leopoldo de Meis, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

  • 4. Laboratório de Microbiologia Celular, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil

Abstract

Surface-associated proteins from Mycobacterium bovis BCG Moreau RDJ are important components of the live Brazilian vaccine against tuberculosis. They are important targets during initial BCG vaccine stimulation and modulation of the host's immune response, especially in the bacterial-host interaction. These proteins might also be involved in cellular communication, chemical response to the environment, pathogenesis processes through mobility, colonization, and adherence to the host cell, therefore performing multiple functions. In this study, the proteomic profile of the surface-associated proteins from M. bovis BCG Moreau was compared to the BCG Pasteur reference strain. The methodology used was 2DE gel electrophoresis combined with mass spectrometry techniques (MALDI-TOF/TOF), leading to the identification of 115 proteins. Of these, 24 proteins showed differential expression between the two BCG strains. Furthermore, 27 proteins previously described as displaying moonlighting function were identified, 8 of these proteins showed variation in abundance comparing BCG Moreau to Pasteur and 2 of them presented two different domain hits. Moonlighting proteins are multifunctional proteins in which two or more biological functions are fulfilled by a single polypeptide chain. Therefore, the identification of such proteins with moonlighting predicted functions can contribute to a better understanding of the molecular mechanisms unleashed by live BCG Moreau RDJ vaccine components.

Introduction

Tuberculosis (TB) is one of the 10 major causes of death worldwide. According to the World Health Organization (WHO), TB killed 1.7 million people in 2016 with 10.4 million new cases estimated worldwide ratifying the need for more effective treatment and prevention (1). To date, Bacillus Calmette-Guérin (BCG) is the only widely used prophylactic measure against TB (1). BCG is an attenuated Mycobacterium bovis strain obtained at the beginning of the Twentieth century at the Pasteur Institute, in Lille. It was distributed to more than 34 countries and maintained in culture without adequate standardization for a long time (2). Because of in vitro evolution, slightly different BCG substrains emerged from the parental BCG, such as those with deletions and duplications of genomic regions and/or single nucleotide polymorphism (SNPs) well documented through analysis and genome sequencing (39). These genetic differences among the various BCG strains in use worldwide partially explain the variable efficacy in protection against pulmonary TB in adults (2). In Brazil, the strain used for vaccine production since 1927 is M. bovis BCG Moreau. The genomic comparison between BCG Pasteur, reference strain, and BCG Moreau showed regions of difference (RD), for example the loss of RD2 and RD14 in BCG Pasteur and RD16 in BCG Moreau, leading to unique genomic/proteomic characteristics (5, 10). These particularities justify more detailed proteomic studies in order to elucidate which proteins are effectively expressed by these bacteria.

BCG is a live attenuated vaccine, and the expression of secreted and surface-associated proteins is extremely relevant since these proteins may play a role in the bacteria-host cell interaction at the beginning of infection (11, 12). These proteins can also be released from the surface since they are non-covalently linked to the mycomembrane (13). Many intracellular proteins with a known function in cell metabolism have also been found on the cell surface, such as glutamine synthetase, gamma-glutamyl phosphate reductase, and cysteine desulfurase (14). Different factors may contribute to the switch between functions such as release to the extracellular space, changes in temperature, redox state of the cell, oligomeric state of the protein, direct interactions with a variety of binding partner proteins, or even to changes in the cellular concentration of a ligand/substrate, cofactor or product, bringing to light the importance of surface-associated proteins playing different roles in cell system (15, 16). This switch between functions is a characteristic called moonlighting (15, 17).

Considering this variability on protein expression, localization and function(s) in different strains of BCG, we compared the surface-associated proteome from the Brazilian strain used in TB vaccine production, M. bovis BCG Moreau, to that of BCG Pasteur through 2DE gel electrophoresis combined with mass spectrometry. Complementary to the information already available in the literature, our approach allows a more confident evaluation of expression, abundance, localization and function(s) of proteins between these two BCG strains. The results presented here may lead to the identification of key components of the M. bovis BCG Moreau vaccine strain which can be related to the variability in immunological response observed in vaccinated individuals.

Materials and Methods

Mycobacterial Cultivation

Mycobacterium bovis BCG Pasteur strain 1173P2 was obtained from the Pasteur Institute (Paris, France) and seed-stocks maintained at −80°C. M. bovis BCG Moreau was supplied in Sauton/potato medium by the Ataulpho de Paiva Foundation (FAP), producer of the BCG vaccine in Brazil. Both strains were cultured as surface pellicles for 2 weeks at 37°C in Sauton medium (18).

Surface Fraction Preparation

After the removal of the culture filtrate, the surface associated material was obtained through vigorous manual shaking of the bacterial pellicle with 2 mm glass beads, essentially as described (19). The surface components were recovered in Milli-Q sterile water and centrifuged twice at 2,500 g for 10 min at 4°C. Aliquots of 1 mL were further centrifuged twice at 16,000 g in order to remove any remaining bacteria. The surface-associated protein fraction was obtained using an adaptation of the method described by Wessel and Flugge (20). Briefly, proteins were precipitated with 15% TCA/acetone, the resulting pellets were washed sequentially with 400 μL of 100% cold acetone, 200 μL of diethyl ether, and 200 μL of chloroform. The final pellets were resuspended in 100 μL Isoelectric focusing (IEF) buffer (8M ureum/2% CHAPS).

Bi-Dimensional Electrophoresis

The immobilized pH gradient (IPG) strips and all 2DE reagents were purchased from Bio-Rad, (Hercules, CA, USA). For the first dimension, 500 μg of proteins were diluted to a final volume of 300 μL of rehydration solution (8M urea, 2% CHAPS), 4 mM Tributyl phosphine (TBP), 0.4% ampholytes pH 3-10, trace of bromophenol blue). The samples were applied to IPG strips (17 cm, pH interval of 4–7) by in-gel rehydration and incubated for 1 h at room temperature. All isoelectric focusing was performed on a Protean® IEF cell (Bio-Rad) with a temperature of 20°C and a maximum current of 50 μA/strip. Running conditions: active rehydration (50V) for 11 h; step 1- linear gradient from 1 to 250V over 20 min; step 2 - linear gradient from 250 to 10,000V over 2 h; step 3- constant 10,000V until complete 80,000 V/h. After isoelectric focusing, proteins were reduced in 130 mM dithiothreitol (DTT) and alkylated in 270 mM iodoacetamide, both in equilibration buffer (6M urea, 2% SDS, 375 mM Tris-HCl pH 8.8, 20% glycerol). Second dimension separation was done in 17 cm, 12% SDS-PAGE gels, 1.0 mm thick, using a vertical system (Bio-Rad) in standard Tris/glycine/SDS buffer at 40 mA/gel, 10°C, until the tracking dye left the gel. Proteins were visualized with Comassie Brilliant Blue (CBB) following procedures described elsewhere (21).

Image Analysis

Gel images were documented using a GS-800TM calibrated imaging densitometer (Bio-Rad) and images were analyzed using PDQuestTM software (Bio-Rad). During the alignment of the images, to compensate for subtle differences in sample loading, gel staining, and destaining, the volume of each spot was normalized in relation to the total density of valid spots present in the gel image. Comparison of 2DE maps derived from three independent protein preparations, each one obtained from three independent BCG cultures, was performed. To determine experimental isoelectric point (pI) and molecular mass (Mr) coordinates for each single spot, 2DE gels were calibrated using a select set of reliable identification landmarks distributed throughout the entire gel. The theoretical pI and Mr of proteins identified by mass spectrometry were obtained using the BCG Moreau RDJ genome reference sequence (9; accession number: AM412059.2).

Protein Digestion, Peptide Extraction and MALDI-TOF/TOF Analysis

In-gel digestion of the 2DE SDS-PAGE separated proteins was carried out using the procedure according to Shevchenko (22). Briefly, protein spots were excised, and the gel pieces were washed three times with 50% (v/v) acetonitrile in 25 mM ammonium bicarbonate for 15 min each, dehydrated in acetonitrile, and dried in a vacuum centrifuge. Gel pieces were rehydrated in 15 μL of 50 mM ammonium bicarbonate containing 20 ng of sequencing grade modified trypsin (Promega). After 15 min, 20 μL of 50 mM ammonium bicarbonate was added to keep the gel pieces wet during tryptic digestion (37°C, 16 h). To extract peptides, 20 μL of 0.5% (v/v) trifluoroacetic acid (TFA) in 50% (v/v) acetonitrile were added and samples were sonicated for 30 min. The separated liquid was concentrated under vacuum to an approximate volume of 10 μL. The resulting peptides were extracted, partially dried, and salts were removed using ZipTipC18 columns (Millipore, Bedford, MA) following the manufacturer's instructions. The tryptic peptides were analyzed on a 4700-Proteomics Analyzer MALDI-TOF/TOF (Applied Biosystems, Foster City, CA). All mass spectra were acquired on positive ion reflector mode with 2,000 shots per spot and externally mass calibrated with a peptide mixture. The 10 most intense ion peaks from the peptide mass fingerprinting (or MS run) were further submitted to fragmentation using post source decay (PSD) mode with collision induced dissociation (CID) gas off and 1 keV collision energy.

Data Analyses and Protein Identification

Following MS/MS acquisition, the processed data files (ppw files) from the MALDI-TOF/TOF were analyzed on a Mascot Server license v. 2.2 (23, 24). The mass spectra were searched against the M. bovis BCG str. Moreau protein database (9; accession number: AM412059.2). The parameters used for the search were as follows: peptide and fragment ions mass tolerance was set at 0.5 Da; maximum of one miss cleavage site by trypsin; carbamidomethylation of cysteine residues as fixed modification, whereas oxidation of methionine/tryptophan, acetylation of the N-terminal, pyroglutamic acid, pyroglutamine, and deamidation of asparagine/glutamine were considered as variable modifications. Positive protein hit identification was accepted with at least 1 matched unique peptide. False discovery rate was estimated at <1%. The peptide ion score was considered >15 with a significance threshold of p < 0.05, whereas the protein score was above 20. Finally, a good correlation between the experimental and theoretical molecular mass and pI was also considered for positive identifications. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (25) with the dataset identifier PXD006141.

Statistical Analysis

Differences between spot intensity observed in 2DE gel images of M. bovis BCG strains Moreau and Pasteur were considered statistically significant when *p < 0.05; **p < 0.01; or ***p < 0.001, as determined by Graph Pad Prism 4.0 software (Graph- Pad Software Inc., San Diego, CA, USA). The unpaired Student‘s t-test was used to analyze the significant differences among both strains, using data obtained from at least 3 different sets of independent biological samples.

Bioinformatic Data

Regions of difference RD2 and RD14 both present on BCG Moreau were evaluated in terms of protein identification and predicted function according to BCG Moreau genome while the region of difference RD16 absent from BCG Moreau was evaluated according to M. tuberculosis H37Rv reference genome (accession numbers: AM412059.2 and NC_000962, respectively). For this purpose, TubercuList (http://genolist.pasteur.fr/TubercuList/) was used as the reference site to find the information about predicted function and amino acid sequence from all the proteins in M. tuberculosis H37Rv encompassing RD2, RD14, and RD16 allowing the construction of Table 1.

Table 1

Region of DifferenceORFs completely removed from the genomeOrthologs H37Rv/BCG MoreauProtein identificationPredicted function
RD2 (10,787 bp)Rv1979c–Rv1987Rv1979c/BCGM_1980cPossible conserved permeaseUnknown; Possibly involved in transport of amino acid across the membrane
Rv1980c/mpb64/BCGM_1981cImmunogenic protein mpt64 (antigen mpt64/mpb64)Unknown
Rv1981c/nrdF1/BCGM_1982cRibonucleoside-diphosphate reductase (beta chain)NRDF1 (Ribonucleotide reductase small subunit) (R2F protein)Involved in the DNA replication pathway. Catalyzes the biosynthesis of deoxyribonucleotides from the corresponding ribonucleotides, precursors that are necessary for DNA synthesis (catalytic activity: 2-deoxyribonucleoside diphosphate + oxidized thioredoxin + H2O = ribonucleoside diphosphate + reduced thioredoxin)
Rv1982c/BCGM_1983cConserved hypothetical proteinUnknown
Rv1983/PE_PGRS35/BCGM_1984PE-PGRS family proteinUnknown
Rv1984c/cfp21/BCGM_1985cProbable cutinase precursor cfp21Hydrolyzes cutin
Rv1985c/BCGM_1986cProbable transcriptional regulatory protein (probablyInvolved in transcriptional mechanism
Lysr- family
Rv1986/BCGM_1987Probable conserved integral membrane proteinUnknown; possibly involved in transport of lysine across the membrane
Rv1987/BCGM_1988Possible chitinaseHydrolysis of chitin
RD14 (9,068 bp)Rv1765A–Rv1772Rv1765A/BCGM_1774cPutative transposase (fragment)Possibly required for the transposition of an insertion element
Rv1766/BCGM_1775Conserved hypothetical proteinUnknown
Rv1767/BCGM_1776Conserved hypothetical proteinUnknown
Rv1768/PE_PGRS31/BCGM_1777PE-PGRS family proteinUnknown
Rv1769/BCGM_1778Conserved hypothetical proteinUnknown
Rv1770/BCGM_1779Conserved hypothetical proteinUnknown
Rv1771/BCGM_1780Probable oxidoreductaseUnknown; Probably involved in cellular metabolism
Rv1772/BCGM_1781Hypothetical proteinUnknown
RD16 (7,608 bp)Rv3401–Rv3404cRv3401/*ΔConserved hypothetical proteinUnknown; Probably enzyme involved in cellular metabolism
Rv3402c/*ΔConserved hypothetical proteinUnknown; Thought to be involved in cell process
Rv3403c/*ΔHypothetical proteinUnknown
Rv3404c/*ΔConserved hypothetical proteinUnknown

Proteins encoded by the genes localized in the regions of difference RD2, RD14, and RD16.

means deleted from the genome of BCG Moreau.

Putative signal peptide for protein export were predicted using SignalP 4.1 (2628), LipoP 1.0 (29, 30), TatP 1.0 (31, 32) and SecretomeP 2.0 (33, 34) in order to predict protein localization. Potential transmembrane domains were predicted with TMHMM 2.0 (35, 36). Beta-barrel membrane proteins structural subclass from integral membrane proteins were discriminated using a Hidden Markov Model method with PRED-TMBB (3739).

The protein sequences described with moonlighting function were retrieved from the MoonProt (40, 41) or Multitasking Protein databases (42) as well as references on Table 3 and Table S3. These sequences were used to construct the sequence database used in comparison to the sequences of the proteins identified in this work. The similarity search was performed using BLASTP (BLOSUM62 matrix), assuming the results with bitscore > 50 and E-value < e-10. Domain analyses were accomplished with proteins that obtained hit in the MoonProt and Multitasking databases using CDD (43, 44) and PFAM (45, 46).

Results

Proteins Encoded by the Genes From Regions of Difference RD2, RD14, and RD16

According to the literature, the regions of difference RD2 and RD14 are present in BCG Moreau, whereas RD16 is absent, when compared to BCG Pasteur (5). Table 1 lists the proteins encoded by these RDs. Our analysis on the surface associated proteome map of BCG Moreau (Figure 2 and Table 2) detected Mpb64, encoded in RD2, and revealed the Rv3406 protein, whose regulation is affected by the truncation of gene rv3405c, due to deletion of RD16 (10).

Table 2

Spot no.GeneProtein identificationFold difference M/P**p-value
32*BCGM_3440Alpha-ketoglutarate dependent sulfate ester dioxygenase14.150.0017
33*518.630.0021
39mpb70Secreted immunogenic protein Mpb704.100.1366
49*23.74<0.0001
501585.000.0558
51*75.110.0346
520.380.2983
56*8.990.0031
5819.740.0709
43*mpb64Immunogenic protein Mpt64 (lost in BCG Pasteur due to RD2)9.660.0022
44*21.670.0002
119ahpCAlkyl hydroperoxide reductase C protein1.630.0885
1203.110.0524
121*13.70<0.0001
26*adoKAdenosine kinase7.710.0217
14fadA3Probable beta-ketoacyl CoA thiolase2.560.0538
24*6.260.0104
1451.090.8296
146*2.470.0464
60*clpP2ATP-dependent Clp protease proteolytic subunit 25.570.0302
751.430.1921
42*BCGM_0830cFatty acid binding protein-like protein (UPF0678)5.140.0043
64echA3Probable enoyl-CoA hydratase (crotonase)88.560.4390
118*4.95<0.0001
38*fixBElectron transfer flavoprotein (alpha-subunit)4.120.0361
1331.250.2959
1341.230.5248
67*ppiAIron regulated peptidyl-prolyl cis-trans isomerase A3.810.0267
683.040.1288
691.370.5389
1170.790.3860
9apaAlanine and proline rich secreted protein0.900.7515
147*2.830.0125
148*3.400.0015
46*cfp17. garAGlycogen accumulation regulator GarA1.760.0120
470.470.1294
48*3.390.0160
149*TB39.8. fhaAFHA domain-containing protein3.040.0311
1502.420.1602
1511.040.9438
1521.290.4435
154*BCGM_1880cProbable reductase2.440.0296
1551.460.1402
124*prcB20S Proteasome (beta subunit)1.930.0229
4*dnaKChaperone protein DnaK (Hsp70)1.620.0271
21*fadAPossible acyl-CoA thiolase0.260.0025
107*esxJESAT-6 like protein EsxJ0.510.0081
74TB18.6Conserved hypothetical protein (UPF0098)0.500.2573
87*0.560.0483
890.760.0707
5*groEL2Chaperonin 2, GroEL2 (65 kDa antigen; Hsp65)0.680.0257
66*fbaFructose-biphosphate aldolase0.690.0209
1640.480.3384
122*gpm1Phosphoglycerate mutase 10.710.0310
72*ssbSingle-strand binding protein0.770.0467
730.410.0930

Differential surface-associated proteins between M. bovis BCG Moreau and M. bovis BCG Pasteur.

Differences in the spot intensity of BCG Moreau and Pasteur 2DE were considered by performing unpaired statistical analysis.

*

p < 0.05.

**

Ratio of mean pixel intensity value for the specified protein spot in BCG Moreau (M) vs. Pasteur (P).

Identification of M. bovis BCG Moreau Surface-Associated Proteins From 2DE Gels Using MALDI-TOF-TOF

The first goal of this study was to perform a proteomic analysis of surface-associated proteins of M. bovis BCG Moreau through 2DE gel electrophoresis and MALDI-TOF-TOF and compare it to M. bovis BCG Pasteur. The confirmation of the proteomic profile and the differences in protein expression between the two strains analyzed were done by comparing the surface proteomic maps in the pH 4–7 range through five biological replicates. Figure 1A shows a representative map of the surface-associated proteins of BCG Moreau. According to the statistical analysis and overlapping of the processed gel image (or virtual image), all replicates presented a reproducible profile related to the total number of spots as well as the localization (migration) and intensity (data not shown). A total of 173 protein spots are reported and they ranged in Mr between 19 and 97 kDa, mostly concentrated above the 31 kDa range. Analysis of the 2DE profiles showed that, in some cases, different spots were identified as representing the same peptide sequence (Table S1), possibly due to the occurrence of proteolysis, different protein isoforms, post-translational modifications (PTMs) or the formation of complexes or heterodimers between proteins. These events may cause differences between theoretical and experimental Mr and pI. The occurrence of different proteins identified in the same spot was also observed and can be explained by the limit of resolution in the 2DE technique (47). Among the 173 spots detected, 115 different proteins were identified by mass spectrometry (Table S1), which were classified according to the M. bovis BCG Moreau gene annotation and orthologs in M. tuberculosis H37Rv and M. bovis BCG Pasteur (accession numbers: AM412059.2, NC_000962 and AM408590.1, respectively). Functional classification of all identified proteins per spot, according to the TubercuList databank (http://genolist.pasteur.fr/TubercuList/) revealed that the major group composed by 70 proteins are related to intermediary metabolism/respiration whereas 37 proteins, the second major group, are conserved hypothetical proteins, without associated function (Figure 1B).

Figure 1

Bioinformatics Analyses of the Identified Surface-Associated Proteins

In order to predict the type of secretion process, transmembrane portion unit on the surface-associated proteins and their comparison with the membrane and secreted proteins already described in the literature, a combination of bioinformatic online programs and mycobacterium genome evaluation was used (Table S2). 2DE comparative quantification with non-paired t-test statistical analysis allowed the identification of 31 spots differentially expressed between BCG Moreau and Pasteur that represented 24 distinct proteins (Figure 2 and Table 2)−17 proteins were more abundant in BCG Moreau and 7 proteins in BCG Pasteur (Table 2). All 115 different proteins identified by mass spectrometry were searched against the Moonprot and Multitasking databases and we could predict moonlighting function for 27 of them (Table 3). Moreover, among the 27 moonlighting proteins, 5 of them were more abundant in BCG Moreau whereas 3 protein spots were up-regulated in BCG Pasteur (Table S3). Domain analysis on the moonlighting predicted proteins using CDD and pFam databases allowed the identification of two domain hits for fructose-biphosphate aldolase (Fba) and aldehyde dehydrogenase (AldC), thus, suggesting distinct functions (data not shown). However, the spot corresponding to the Fba protein was found 2-fold more expressed in BCG Pasteur than in BCG Moreau.

Figure 2

Table 3

GeneBCG MoreauPredicted functionMoonlighting functionOrganism/speciesReferences
ppiABCGM_0009Probable iron-regulated peptidyl-prolyl cis-trans isomerase AProinflammatory cytokine-activate endothelial cellsHomo sapiens(48)
Induces apoptosis of gastric epithelial cellsHelicobacter pylori(49)
Ppiases might play a role in the intracellular survival by subverting the host cell defenses, such as oxidative stress as well as by immunomodulationMycobacterium tuberculosis(50)
BCGM_0129Probable serine proteaseChaperone, aids in folding of proteinsEscherichia coli(51)
Escherichia coli(52)
Escherichia coli(53)
BCGM_0152BCGM_0152Probable short-chain type dehydrogenase/reductaseTranscriptional repressorHomo sapiens(54)
gabD1BCGM_0238cProbable succinate-semialdehyde dehydrogenase [nadp+] dependant (ssdh)Tau-crystallin/alpha-enolasePseudemys scripta; Petromyzon marinus(55)
Transcriptional repressorEscherichia coli; Salmonella typhimurium(5659)
fabG4BCGM_0247cProbable 3-oxoacyl-[acyl-carrier protein] reductaseTranscriptional repressorHomo sapiens(54)
dnaKBCGM_0357Probable chaperone proteinPlasminogen bindingBifidobacterium(60)
Neisseria meningitides(61)
fba §BCGM_0369cProbable fructose-biphosphate aldolasePlasminogen bindingCandida albicans(62)
V-ATPase assemblySaccharomyces cerevisiae(6365)
AdhesinStreptococcus pneumoniae; Neisseria meningitides(66, 67)
groEL2BCGM_044760 kDa chaperonin 2, GroEL2ToxinEnterobacter aerogenes(68)
Invasion in a HeLa cell modelLegionella pneumophila(69)
Binds DNAMycobacterium tuberculosis(70)
Regulation of inflammationMycobacterium leprae(71)
Aggravation of atherosclerosisChlamydia pneumoniae(72)
Role in the growth of the bacteriaHelicobacter pylori(73)
Essential for biofilm formationMycobacterium smegmatis(74)
Stimulation of osteoclastogenesisEscherichia coli(75)
LpdBCGM_0470Dihydrolipoamide dehydrogenase LpdProteaseHomo sapiens(76)
gpm1BCGM_0498Probable phosphoglycerate mutase 1Plasminogen binding proteinStreptococcus sp. Oral(77)
echA3BCGM_0644cProbable enoyl-CoA hydrataseRNA and microtubule binding proteinRice seed(78)
tufBCGM_0699Probable elongation factor TU tuf (EF-TU)Attachment to human cells and mucinsLactobacillus johnsonii NCC533 (La1)(79)
Fibronectin, plasminogen, and mucin bindingMycoplasma pneumoniae; Candida albicans; Streptococcus gordonii; Pseudomonas aeruginosa(62, 8082)
Cytoskeleton structure regulationCereal grains(83)
enoBCGM_1046Probable enolaseBinds plasminogen, laminin, fibronectin, mucin, and actinLactobacillus crispatus ST1 and Lactobacillusjohnsonii F133; Streptococcus gordonii; Lactobacillus plantarum; Streptococcus suis; Borrelia burgdorferi; Candida albicans; Staphylococcus aureus; Neisseria meningitidis; Aeromonas hydrophila; Echinococcus granulosus; Trichomonas vaginalis(61, 81, 8492)
It has neurotrophic and neuroprotective effects on rather a broad spectrum of neurons in the central nervous system (promotes cell survival)Rattus norvegicus(93)
Activates homotypic vacuole fusionSaccharomyces cerevisiae(94)
Depending on the localization may be involved in food vacuole formation or may play a role in transcriptionPlasmodium falciparum(95)
Mitochondrial tRNA importSaccharomyces cerevisiae(96)
glyA1BCGM_1119Probable serine hydroxymethiltransferase 1mRNA binds the 5' untranslated region (UTR) of its own mRNAHomo sapiens(97)
fumBCGM_1124cProbable fumarase fum(Fumarate hydratase)Crystallins (alpha-, beta-gama-, delta-, epsilon- and tau- crystallins)Vertebrate lenses(98)
Tumor suppressorHomo sapiens(99)
gapBCGM_1465Probable glyceraldehyde 3 phosphate dehydrogenaseMucin, fibronectin, laminin, plasminogen, and fibrinogen binding proteinMycoplasma genitalium; Candida albicans; Candida albicans; Candida albicans; Streptococcus pneumoniae; enterohemorrhagic and enteropathogenic Escherichia coli(62, 100103)
Single-gene disordersHomo sapiens(104)
Regulatory role in the bundling/unbundling of microtubules; involvements in the assembly of junctional triads from transverse tubules in skeletal muscle cells; protein kinase activity; single- stranded DNA-binding activity regulating transcription; human Nuclear uracil DNA glycocylase activityGroup A Streptococci(105)
Binds uPAR/CD87 receptor on human cellsStreptococcus pyogenes(106)
Gene-specific silencing of translationChordates(107)
Oral biofilmStreptococcus oralis(108)
Transferrin receptor, NAD-ribosylating activityStaphylococci(109)
Role in extracellular polysaccharide synthesisXanthomonas campestris pv. Campestris(110)
pgkBCGM_1466Probable phosphoglycerate kinaseInteracts with plasminogen and its tissue activatorStreptococcus pneumoniae(111)
Acts in tumor angiogenesis as a disulphide reductaseHomo sapiens(112)
canBCGM_1505cProbable aconitate hydratase AcnIron responsive element binding protein—iron uptake and usemammalian tissues (beef liver)(113)
Iron-dependent RNA-binding activity—mRNA binding proteinHomo sapiens; Mycobacterium tuberculosis(114, 115)
Protection to mitochondrial DNA and regulation to mitochondrial DNA maintenanceSaccharomyces cerevisiae(116, 117)
BCGM_1880cBCGM_1880cProbable reductaseProteaseEukaryotes(76)
glnA1BCGM_2221Glutamine synthetaseControl of gene expressionBacteria(118)
Adhesion to hostLactobacillus crispatus ST1(119)
ahpEBCGM_2239cPeroxiredoxin AhpEPlasminogen bindingCandida albicans(62)
Molecular chaperone functionYeast(120)
Cell wall biogenesisCandida albicans(121)
ahpCBCGM_2431Alkyl hydroperoxide reductase c proteinMolecular chaperone functionYeast(120)
Plasminogen bindingCandida albicans(62)
Cell wall biogenesisCandida albicans(121)
Bifunctional enzyme with Glutathione peroxidase and phospholipase A2 activitiesHomo sapiens(122)
ndkABCGM_2449cProbable nucleoside diphosphate kinase ndkARegulate cellular redox state and enhances multiple stress tolerance in transgenic plantsArabidopsis thaliana(123)
aldC #BCGM_2863cProbable aldehyde dehydrogenaseEta-crystallin; Tau-crystallin/alpha-enolaseElephantulus edwardii; Anas platyrhynchos/Pseudemys scripta/Petromyzon marinus;(55, 124)
Transcriptional regulatorSalmonella typhimurium(59)
adhCBCGM_3053Probable NADP-dependent alcohol dehydrogenaseZeta-crystallinHyla japonica; Guinea Pig(125128)
groEL1BCGM_3451c60 kDa chaperonin 1, GroEL1Associated with nucleoidsMycobacterium tuberculosis(70)
ToxinEnterobacter aerogenes(68)
Invasion in a HeLa cell modelLegionella pneumophila(69)
Involved in mycolic acid biosynthesis during biofilm formation in mycobacteriaMycobacterium smegmatis(74)
Regulation of inflammationMycobacterium leprae(71)
Regulates macrophage tumor necrosis factor-alpha and matrix metalloproteinase expressionChlamydia pneumoniae(72)
Role in the growth of the bacteriaHelicobacter pylori(73)
Stimulator of osteoclast formationEscherichia coli(75)
groESBCGM_3452c10 kDa chaperoninActivate human macrophages in synergy with LPSFrancisella tularensis(129)

M. bovis BCG Moreau surface-associated proteins with moonlighting predicted function.

§

Two domain hits: FTBP_aldolase_II (canonical hit) and ICL_KPHMT (moonlighting hit).

#

Two domain hits: ALDH_F1AB_F2_RALDH1 (canonical hit) and putA (moonlighting hit).

Discussion

In this study, the surface-associated proteins from M. bovis BCG Moreau were investigated by 2DE combined with MALDI-TOF-TOF and bioinformatic analysis. Two-dimensional gel electrophoresis maps in the pH range of 4–7 led to the identification of 173 spots (Figure 1A) that could be assigned by mass spectrometry to a total of 115 different proteins (Table S1). This choice of pH range was based on previous analysis carried out in the broader pH range of 3–10, which showed that the majority of protein spots occurred in the pH range of 4–7 (data not shown). These results complement our previous report on the secretome of BCG strains Moreau and Pasteur (12) (Table 1).

The genome of BCG Moreau differs from BCG Pasteur, among others, by the presence of two regions (RD2 and RD14; 5). Several proteins encoded in these two regions have been described as potentially immunogenic (130133), reinforcing the importance of performing studies to address the protein expression in these two strains. In this context, differentially expressed proteins present on the bacterial surface can partially account for some differences observed in individual response to vaccination. Among the proteins identified by mass spectrometry, the Mpb64 protein (encoded in RD2) was observed to be expressed in M. bovis BCG Moreau. The quantification performed by the PDQuest software indicated a very low expression of the Mpb64 protein in BCG Pasteur proteome map, probably due to the background observed in the area where the protein should be located (Figure 2, spots 43 and 44).

Mpb64 is an immunodominant antigen capable of inducing protective immunity by T cell response; however the role of this protein in the pathogenesis of tuberculosis is not known (134, 135). We found Mpb64 expression on the cell surface proteome of BCG Moreau (spots 43 and 44, Figure 1A) with the expected predicted signal peptide and transmembrane domains (Table S2), suggesting that vaccination with BCG Moreau could trigger a better cellular immune response when compared to BCG Pasteur. With respect to RD16, this characteristic deletion found in BCG Moreau (between genes rv3400, and rv3405c) leads to the truncation and functional loss of a TetR transcriptional regulator encoded by rv3405c, resulting in the constitutive expression of the adjacent rv3406 gene (10); the resulting protein (Rv3406) was identified in the surface proteome of BCG Moreau, and is absent in Pasteur (Figure 2).

Considering the surface-associated proteomic profiles from M. bovis BCG Moreau and Pasteur, differences between the experimental and theoretical protein molecular mass and pI were observed for 49% of the identified spots, which may be due to post-translational modifications (PTMs), protein degradation, the presence of isoforms and conformers of the proteins (136, 137). For example, Apa (spots 9, 147, and 148) is already described as glycosylated on diverse threonine residues, which could result in the observation of multiple protein species by 2DE gel (138). In fact, the experimental molecular mass and pI of spots 9, 147, and 148 are in disagreement to the theoretical ones (Table S1). We observed that the ESAT-6 like protein EsxJ and glutamine synthetase (GS) found in spots 107 and 169/170, respectively, were identified as N-formylation-containing proteins (Table S1). Such PTM could contribute to electrophoretic mobility shift in 2DE gels, and therefore, would account for the differences between experimental and theoretical pI. It has been reported that N-formylated peptides may serve as good candidates for a universal vaccine against M. tuberculosis when administered in combination with drugs (139141).

Overall, these PTMs could be important factors for eliciting the immune response. It has been reported that the glycosylated motif found in Apa protein was related to the high capacity of BCG to stimulate an immune response in BCG-immunized guinea pigs (142, 143). Our results indicate that the Apa protein shows higher expression on the cell surface of BCG Moreau compared to Pasteur (Figure 2 and Table 2). Other important proteins identified in the BCG Moreau surface-associated proteome were secreted antigens 85A and B (spots 37, 126, and 127 of Figure 1A). These proteins belong to the antigen Ag-85 complex and are found in association with the mycobacteria cell surface, constituting the major secreted proteins observed on mycobacteria culture filtrate (12, 144, 145). These proteins are strongly immunogenic and can trigger both humoral and cellular immune responses in vivo (144, 146148). Ag85 proteins have mycolyltransferase activity, crucial to maintain the cell wall integrity of mycobacteria (149151). In addition, they can bind to human fibronectin (152) and elastin (153) present in the extracellular matrix having, therefore, moonlighting functions. Recent studies demonstrated that Ag-85 complex was identified as a target for mannose-binding lectin and ficolins (154). These proteins have been extensively studied as potential candidates for new vaccines against tuberculosis (155, 156).

Protein fate in terms of cellular localization is an important aspect that might be explored on surface-associated proteins (Table S2). Non-covalently surface-associated proteins can be secreted through different systems such as via signal peptide (as predicted by Signal P and Lipo P software), non-classical secretory pathway like proteins without an N-terminal signal peptide (according to Secretome P analysis), and twin-arginine translocation pathway (as found by Tat prediction program) in which a twin-arginine consensus motif is located within the signal peptide itself (29, 31, 157159). The non-classical secretory system was predicted for aconitase (spots 172 and 173) and glutamine synthetase (spots 169 and 170). In M. tuberculosis, aconitase is found in the cytosol, cell wall, and cell membrane fractions and it has also been described as a bifunctional protein acting as an enzyme in the presence of iron and RNA-binding in the absence of iron (115, 160). Glutamine synthetase (GS) presented moonlighting features of acyltransferase in M. tuberculosis (161). Previous studies have shown that this enzyme is secreted into the culture medium and plays a crucial role in pathogenicity as well as in bacterial growth (162, 163). In Lactobacillus crispatus, GS is a novel adhesive moonlighting enzyme that associates to the cell surface at an acidic pH (119). Further experiments must be carried out in order to show that aconitase, glutamine synthetase and other surface-associated proteins (Table 3) have moonlighting behavior also in M. bovis BCG Moreau.

Bioinformatic tools have been used to predict protein moonlighting function by primary sequence analysis and hence domain investigations of predicted moonlighting proteins can be pursued to propose novel functions that corroborate with in vitro and in vivo studies (40, 42, 164, 165). Nevertheless, most moonlighting proteins described to date have been identified by chance (166). In general, highly conserved proteins, often metabolic proteins/enzymes (167, 168) or molecular chaperones (169), receptors (170), ribosomal proteins, and transmembrane channels (171), were shown to be moonlighting proteins (172). These findings suggest that the presence of intracellular proteins at “unexpected” locations is not always due to experimental artifacts such as cellular lysis. The methodology used here for obtaining the fraction of surface-associated proteins is well characterized (19), strengthening the fact that the intracellular proteins identified in this study may be performing moonlighting functions. In addition, the appearance of a new function for the same protein can be considered a great advantage for the microorganism since it optimizes the functional repertoire encoded by a compact genome. In this context, bioinformatic investigation using the Moonprot/Multitask programs led to the prediction of moonlighting functions for some of the identified proteins, already described for other organisms/species (Table 3). Besides, Conserved Domains Database (CDD) and pFam were also used to infer the number of domains these proteins could present (data not shown). Most of the identified proteins grouped in the “intermediary metabolism” functional category (Figure 1B), raising the question of why proteins normally encountered intracellularly had been found surface-associated to the mycomembrane. This could be explained by different functions these proteins may have depending on localization. Finally, other features such as the presence of distinct protein surfaces or domains influencing the oligomeric state of the protein, concentration of cellular ligands, substrates and cofactors must also be considered to presume moonlighting function for these kind of proteins (16, 173175).

Interestingly, we found that fructose-biphosphate aldolase (Fba) and the aldehyde dehydrogenase (AldC) contained two domains hit: one representative of a canonical function in mycobacteria and another representing a probable moonlighting function already described in other organisms (Table 3). In mycobacteria, the canonical domain hit for fructose-biphosphate aldolase is FTBP_aldolase_II representative of an enzyme that controls the condensation of dihydroxyacetone phosphate with glyceraldehyde-3-phosphate to yield fructose-1,6-bisphosphate (176). Nevertheless, the moonlighting domain hit is ICL_KPHMT that represents an enzyme superfamily that catalyzes the formation and cleavage of either P-C (proline-cysteine) or C-C (cysteine-cysteine) bonds (data not shown). In other organisms, as described in Table 3, the moonlighting function of fructose- biphosphate aldolase is related to protein binding and cell adherence. (177) confirmed experimentally by in vitro assays that Fba from M. tuberculosis binds to human plasminogen. This generates the proteolytic enzyme plasmin leading to the breakdown of extracellular matrix and basal membrane proteins, contributing to tissue injury in tuberculosis. More recently, de la Paz Santangelo et al. (178) reported that Fba of M. tuberculosis binds to human plasminogen in a dose dependent manner and is important for M. tuberculosis growth. According to our results Fba is less expressed on the cell surface of BCG Moreau than in BCG Pasteur (Table S3). On the other hand, the canonical domain hit for aldehyde dehydrogenase is ALDH_F1AB_F2_RALDH1 that corresponds to NAD+-dependent retinal dehydrogenase 1 also known as aldehyde dehydrogenase family 1 member A1 (ALDH1A1) in humans. It is a cytosolic enzyme that catalyzes the oxidation of retinaldehyde to retinoic acid (RA). RA is the active metabolite of vitamin A and it is required for spermatogenesis and many other biological processes (179). The moonlighting domain hit is PutA, a trifunctional protein in bacteria: transcriptional regulator, proline dehydrogenase and pyrroline-5- carboxylate dehydrogenase (data not shown). Christgen et al. (180) discovered a membrane binding region on the PutA domain from Escherichia coli AldC that explains the PutA functional switch from self-regulating transcriptional repressor to membrane binding domain. Our results indicate that AldC is a surface-associated protein from M. bovis BCG Moreau (spot 171 of Figure 1A) not found in the culture filtrate of BCG Moreau (Table S2). Fba and AldC may play a role in the immunopathology of tuberculosis, but this still needs further investigation.

Altogether, the differences in abundance of surface-associated proteins identified between BCG strains Moreau and Pasteur could have an impact on vaccine efficacy. The finding that some of these proteins have moonlighting functions opens new possibilities for investigating the role of extracellular proteins on the bacterial-host interface.

Statements

Author contributions

TP, MB-P, and LM-L conceived, designed, performed experiments and analyzed results. AG performed bioinformatic analysis. MW and PC performed statistical analysis of results. DK organized, analyzed, processed, and deposited the mass spectrometry data on PRIDE repository. WD conceived and designed experiments. TP, DK, PC, MB-P, and LM-L wrote the paper.

Funding

This work received financial support from the WHO/TDR Special Programme for Research and Training in Tropical Diseases and the following Brazilian agencies: CAPES, CNPq, and PDTIS/FIOCRUZ.

Acknowledgments

We thank Dr. Jonas Perales and Dr. André Teixeira from Oswaldo Cruz Foundation-FIOCRUZ/PDTIS 2DE and Mass Spectrometry platform facilities for the MALDI-TOF-TOF data acquisition and in particular to Dr. André Teixeira for kindly recruit the raw mass spectrometry data. We are also thankful for the Proteomics and Mass Spectrometry Unit platform (UEMP) at the Rio de Janeiro Federal University (UFRJ) run by Dr. Russolina Benedeta Zingali for support in the use of the license Mascot Server for the protein identification searches. We would like to thank the post-doctoral fellow Dr. Marjolly Caruso Brígido and the technician Augusto Vieira Magalhães both members of UEMP for helping us with the PRIDE mass spectrometry repository. We are also in debt to Carolina Zavareze (FAP) that kindly provided the Sauton culture medium and the BCG Moreau vaccine strain.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.00716/full#supplementary-material

    Abbreviations

  • TB

    Tuberculosis

  • WHO

    World Health Organization

  • BCG

    Bacillus Calmette-Guerin

  • SNPs

    single nucleotide polymorphisms

  • RD

    Region of Difference

  • 2DE gel

    Two-dimensional gel electrophoresis

  • PTMs

    Post-translational modifications

  • FH

    Fumarate hydratase

  • TCA

    trichloroacetic acid

  • GS

    glutamine synthetase

  • TDM

    trehalose dimycolate

  • CDD

    Conserved Domains Database

  • Fba

    fructose-biphosphate aldolase

  • AldC

    aldehyde dehydrogenase.

References

  • 1.

    WHO. Tuberculosis. (2018). Available online at: http://www.who.int/mediacentre/factsheets/fs104/en/ (accessed February 26, 2019).

  • 2.

    Benevolo-de-AndradeTCMonteiro-MaiaRCosgroveCCastello-BrancoLR. BCG Moreau Rio de Janeiro: an oral vaccine against tuberculosis–review. Mem Inst Oswaldo Cruz. (2005) 100:45965. 10.1590/S0074-02762005000500002

  • 3.

    MahairasGGSaboPJHickeyMJSinghDCStoverCK. Molecular analysis of genetic differences between Mycobacterium bovis BCG and virulent M. bovis J Bacteriol. (1996) 178:127482. 10.1128/jb.178.5.1274-1282.1996

  • 4.

    BehrMASmallPM. A historical and molecular phylogeny of BCG strains. Vaccine. (1999) 17:91522. 10.1016/S0264-410X(98)00277-1

  • 5.

    BehrMAWilsonMAGillWPSalamonHSchoolnikGKRaneSet al. Comparative genomics of BCG vaccines by whole-genome DNA microarray. Science. (1999) 284:15203. 10.1126/science.284.5419.1520

  • 6.

    BroschRGordonSVGarnierTEiglmeierKFriguiWValentiPet al. Genome plasticity of BCG and impact on vaccine efficacy. Proc Natl Acad Sci USA. (2007) 104:5596601. 10.1073/pnas.0700869104

  • 7.

    LeungASTranVWuZYuXAlexanderDCGaoGFet al. Novel genome polymorphisms in BCG vaccine strains and impact on efficacy. BMC Genomics. (2008) 9:413. 10.1186/1471-2164-9-413

  • 8.

    Garcia PelayoMCUplekarSKeniryAMendoza LopezPGarnierTNunez GarciaJet al. A comprehensive survey of single nucleotide polymorphisms (SNPs) across Mycobacterium bovis strains and M. bovis BCG vaccine strains refines the genealogy and defines a minimal set of SNPs that separate virulent M bovis strains and M bovis BCG strains. Infect Immun. (2009) 77:22308. 10.1128/IAI.01099-08

  • 9.

    GomesLHOttoTDVasconcellosEAFerraoPMMaiaRMMoreiraASet al. Genome sequence of Mycobacterium bovis BCG Moreau, the Brazilian vaccine strain against tuberculosis. J Bacteriol. (2011) 193:56001. 10.1128/JB.05827-11

  • 10.

    GalvaoTCLimaCRGomesLHPaganiTDFerreiraMAGoncalvesASet al. The BCG Moreau RD16 deletion inactivates a repressor reshaping transcription of an adjacent gene. Tuberculosis (Edinb). (2014) 94:2633. 10.1016/j.tube.2013.11.004

  • 11.

    El-EtrSHCirilloJD. Entry mechanisms of mycobacteria. Front Biosci. (2001) 6:D737747. 10.2741/sahar

  • 12.

    Berredo-PinhoMKalumeDECorreaPRGomesLHPereiraMPda SilvaRFet al. Proteomic profile of culture filtrate from the Brazilian vaccine strain Mycobacterium bovis BCG Moreau compared to M. bovis BCG Pasteur BMC Microbiol. (2011) 11:80. 10.1186/1471-2180-11-80

  • 13.

    DaffeMEtienneG. The capsule of Mycobacterium tuberculosis and its implications for pathogenicity. Tuber Lung Dis. (1999) 79:15369. 10.1054/tuld.1998.0200

  • 14.

    WangWJefferyCJ. An analysis of surface proteomics results reveals novel candidates for intracellular/surface moonlighting proteins in bacteria. Mol Biosyst. (2016) 12:142031. 10.1039/C5MB00550G

  • 15.

    JefferyCJ. Moonlighting proteins–an update. Mol Biosyst. (2009) 5:34550. 10.1039/b900658n

  • 16.

    HubertsDHvan der KleiIJ. Moonlighting proteins: an intriguing mode of multitasking. Biochim Biophys Acta. (2010) 1803:5205. 10.1016/j.bbamcr.2010.01.022

  • 17.

    JefferyCJ. Moonlighting proteins. Trends Biochem Sci. (1999) 24:811. 10.1016/S0968-0004(98)01335-8

  • 18.

    EickhoffTC. The current status of BCG immunization against tuberculosis. Annu Rev Med. (1977) 28:41123. 10.1146/annurev.me.28.020177.002211

  • 19.

    Ortalo-MagneADupontMALemassuAAndersenABGounonPDaffeM. Molecular composition of the outermost capsular material of the tubercle bacillus. Microbiology. (1995) 141:160920. 10.1099/13500872-141-7-1609

  • 20.

    WesselDFluggeUI. A method for the quantitative recovery of protein in dilute solution in the presence of detergents and lipids. Anal Biochem. (1984) 138:1413. 10.1016/0003-2697(84)90782-6

  • 21.

    NeuhoffVAroldNTaubeDEhrhardtW. Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. Electrophoresis. (1988) 9:25562. 10.1002/elps.1150090603

  • 22.

    ShevchenkoAShevchenkoA. Evaluation of the efficiency of in-gel digestion of proteins by peptide isotopic labeling and MALDI mass spectrometry. Anal Biochem. (2001) 296:27983. 10.1006/abio.2001.5321

  • 23.

    PerkinsDNPappinDJCreasyDMCottrellJS. Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. (1999) 20:355167. 10.1002/(SICI)1522-2683(19991201)20:18 < 3551::AID-ELPS3551>3.0.CO;2-2

  • 24.

    MatrixScience. Matrix Science. London, UK. (2016). Available online at: http://www.matrixscience.com (accessed June 1, 2018).

  • 25.

    VizcainoJACsordasAdel-ToroNDianesJAGrissJLavidasIet al. 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. (2016) 44:D44756. 10.1093/nar/gkw880

  • 26.

    NielsenHEngelbrechtJBrunakSvon HeijneG. Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. (1997) 10:16. 10.1093/protein/10.1.1

  • 27.

    PetersenTNBrunakSvon HeijneGNielsenH. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. (2011) 8:7856. 10.1038/nmeth.1701

  • 28.

    SignalP4.1. SignalP 4.1 Server. DK. (2016). Available online at: http://www.cbs.dtu.dk/services/SignalP/ (accessed March 14, 2018).

  • 29.

    JunckerASWillenbrockHVon HeijneGBrunakSNielsenHKroghA. Prediction of lipoprotein signal peptides in Gram-negative bacteria. Protein Sci. (2003) 12:165262. 10.1110/ps.0303703

  • 30.

    LipoP1.0. LipoP 1.0 Server. DK. (2015). Available online at: http://www.cbs.dtu.dk/services/LipoP/ (accessed November 25, 2018).

  • 31.

    BendtsenJDNielsenHWiddickDPalmerTBrunakS. Prediction of twin-arginine signal peptides. BMC Bioinformatics. (2005b) 6:167. 10.1186/1471-2105-6-167

  • 32.

    TatP1.0. TatP 1.0 Server. DK. (2016). Available online at: http://www.cbs.dtu.dk/services/TatP/ (accessed September 16, 2018).

  • 33.

    BendtsenJDKiemerLFausbollABrunakS. Non-classical protein secretion in bacteria. BMC Microbiol. (2005a) 5:58. 10.1186/1471-2180-5-58

  • 34.

    SecretomeP2,.0. SecretomeP 2.0 Server. DK. (2015). Available online at: http://www.cbs.dtu.dk/services/SecretomeP/ (accessed December 14, 2018).

  • 35.

    EmanuelssonOBrunakSvon HeijneGNielsenH. Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. (2007) 2:95371. 10.1038/nprot.2007.131

  • 36.

    TMHMM2,.0. TMHMM 2.0 Server. DK. (2015). Available online at: http://www.cbs.dtu.dk/services/TMHMM/ (accessed August 4, 2018).

  • 37.

    BagosPGLiakopoulosTDSpyropoulosICHamodrakasSJ. PRED-TMBB: a web server for predicting the topology of beta-barrel outer membrane proteins. Nucleic Acids Res. (2004b) 32:W40004. 10.1093/nar/gkh417

  • 38.

    BagosPGLiakopoulosTDSpyropoulosICHamodrakasSJ. A hidden markov model method, capable of predicting and discriminating beta-barrel outer membrane proteins. BMC Bioinformatics. (2004a) 5:29. 10.1186/1471-2105-5-29

  • 39.

    PRED-TMBB. PRED-TMBB Server. (2014). Available online at: http://bioinformatics.biol.uoa.gr/PRED-TMBB/input.jsp (accessed November 5, 2018).

  • 40.

    ManiMChenCAmbleeVLiuHMathurTZwickeGet al. MoonProt: a database for proteins that are known to moonlight. Nucleic Acids Res. (2015) 43:D27782. 10.1093/nar/gku954

  • 41.

    MoonProt. MoonProt. (2014). Available online at: http://www.moonlightingproteins.org/ (accessed December 3, 2018).

  • 42.

    HernandezSFerragutGAmelaIPerez-PonsJPinolJMozo-VillariasAet al. MultitaskProtDB: a database of multitasking proteins. Nucleic Acids Res. (2014) 42:D51720. 10.1093/nar/gkt1153

  • 43.

    Marchler-BauerABoYHanLHeJLanczyckiCJLuSet al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. (2017) 45:D2003. 10.1093/nar/gkw1129

  • 44.

    ConservedDomains. Conserved Domains within a Protein or Coding Nucleotide Sequence. (2015). Available online at: https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed December 6, 2018).

  • 45.

    Pfam. Pfam. (2015). Available online at: http://pfam.xfam.org (accessed December 4, 2018).

  • 46.

    FinnRDCoggillPEberhardtRYEddySRMistryJMitchellALet al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. (2016) 44:D27985. 10.1093/nar/gkv1344

  • 47.

    BunaiKYamaneK. Effectiveness and limitation of two-dimensional gel electrophoresis in bacterial membrane protein proteomics and perspectives. J Chromatogr B Analyt Technol Biomed Life Sci. (2005) 815:22736. 10.1016/j.jchromb.2004.08.030

  • 48.

    JinZGLunguAOXieLWangMWongCBerkBC. Cyclophilin A is a proinflammatory cytokine that activates endothelial cells. Arterioscler Thromb Vasc Biol. (2004) 24:118691. 10.1161/01.ATV.0000130664.51010.28

  • 49.

    BasakCPathakSKBhattacharyyaAPathakSBasuJKunduM. The secreted peptidyl prolylcis, trans-isomerase HP0175 of Helicobacter pylori induces apoptosis of gastric epithelial cells in a TLR4- and apoptosis signal- regulating kinase 1-dependent manner. J Immunol. (2005) 174:567280. 10.4049/jimmunol.174.9.5672

  • 50.

    PandeySTripathiDKhubaibMKumarASheikhJASumanlathaGet al. Mycobacterium tuberculosis peptidyl-prolyl isomerases are immunogenic, alter cytokine profile and aid in intracellular survival. Front Cell Infect Microbiol. (2017) 7:38. 10.3389/fcimb.2017.00038

  • 51.

    KrojerTGarrido-FrancoMHuberREhrmannMClausenT. Crystal structure of DegP (HtrA) reveals a new protease-chaperone machine. Nature. (2002) 416:4559. 10.1038/416455a

  • 52.

    SawaJMaletHKrojerTCanellasFEhrmannMClausenT. Molecular adaptation of the DegQ protease to exert protein quality control in the bacterial cell envelope. J Biol Chem. (2011) 286:3068090. 10.1074/jbc.M111.243832

  • 53.

    MaletHCanellasFSawaJYanJThalassinosKEhrmannMet al. Newly folded substrates inside the molecular cage of the HtrA chaperone DegQ. Nat Struct Mol Biol. (2012) 19:1527. 10.1038/nsmb.2210

  • 54.

    MarkovaNGPinkas-SarafovaASimonM. A metabolic enzyme of the short-chain dehydrogenase/reductase superfamily may moonlight in the nucleus as a repressor of promoter activity. J Invest Dermatol. (2006) 126:201931. 10.1038/sj.jid.5700347

  • 55.

    WistowGJLietmanTWilliamsLAStapelSOde JongWWHorwitzJet al. Tau-crystallin/alpha-enolase: one gene encodes both an enzyme and a lens structural protein. J Cell Biol. (1988) 107(6 Pt 2):272936. 10.1083/jcb.107.6.2729

  • 56.

    WoodJM. Genetics of L-proline utilization in Escherichia coli. J Bacteriol. (1981) 146:895901.

  • 57.

    MenzelRRothJ. Enzymatic properties of the purified putA protein from Salmonella typhimurium. J Biol Chem. (1981) 256:97626.

  • 58.

    Ostrovskyde Spicer PO'BrienKMaloyS. Regulation of proline utilization in Salmonella typhimurium: a membrane-associated dehydrogenase binds DNA in vitro. J Bacteriol. (1991) 173:2119. 10.1128/jb.173.1.211-219.1991

  • 59.

    Ostrovskyde Spicer PMaloyS. PutA protein, a membrane-associated flavin dehydrogenase, acts as a redox- dependent transcriptional regulator. Proc Natl Acad Sci USA. (1993) 90:42958. 10.1073/pnas.90.9.4295

  • 60.

    CandelaMBergmannSViciMVitaliBTurroniSEikmannsBJet al. Binding of human plasminogen to Bifidobacterium. J Bacteriol. (2007) 189:592936. 10.1128/JB.00159-07

  • 61.

    KnaustAWeberMVHammerschmidtSBergmannSFroschMKurzaiO. Cytosolic proteins contribute to surface plasminogen recruitment of Neisseria meningitidis. J Bacteriol. (2007) 189:324655. 10.1128/JB.01966-06

  • 62.

    CroweJDSievwrightIKAuldGCMooreNRGowNABoothNA. Candida albicans binds human plasminogen: identification of eight plasminogen-binding proteins. Mol Microbiol. (2003) 47:163751. 10.1046/j.1365-2958.2003.03390.x

  • 63.

    LuMAmmarDIvesHAlbrechtFGluckSL. Physical interaction between aldolase and vacuolar H+- ATPase is essential for the assembly and activity of the proton pump. J Biol Chem. (2007) 282:24495503. 10.1074/jbc.M702598200

  • 64.

    LuMHollidayLSZhangLDunnWAJrGluckSL. Interaction between aldolase and vacuolar H+-ATPase: evidence for direct coupling of glycolysis to the ATP-hydrolyzing proton pump. J Biol Chem. (2001) 276:3040713. 10.1074/jbc.M008768200

  • 65.

    LuMSautinYYHollidayLSGluckSL. The glycolytic enzyme aldolase mediates assembly, expression, and activity of vacuolar H+-ATPase. J Biol Chem. (2004) 279:87329. 10.1074/jbc.M303871200

  • 66.

    BlauKPortnoiMShaganMKaganovichARomSKafkaDet al. Flamingo cadherin: a putative host receptor for Streptococcus pneumoniae. J Infect Dis. (2007) 195:182837. 10.1086/518038

  • 67.

    TunioSAOldfieldNJBerryAAla'AldeenDAWooldridgeKGTurnerDP. The moonlighting protein fructose-1, 6-bisphosphate aldolase of Neisseria meningitidis: surface localization and role in host cell adhesion. Mol Microbiol. (2010) 76:60515. 10.1111/j.1365-2958.2010.07098.x

  • 68.

    YoshidaNOedaKWatanabeEMikamiTFukitaYNishimuraKet al. Protein function. chaperonin turned insect toxin. Nature. (2001) 411:44. 10.1038/35075148

  • 69.

    GardunoRAGardunoEHoffmanPS. Surface-associated hsp60 chaperonin of Legionella pneumophila mediates invasion in a HeLa cell model. Infect Immun. (1998) 66:460210.

  • 70.

    BasuDKhareGSinghSTyagiAKhoslaSMandeSC. A novel nucleoid-associated protein of Mycobacterium tuberculosis is a sequence homolog of GroEL. Nucleic Acids Res. (2009) 37:494454. 10.1093/nar/gkp502

  • 71.

    RhaYHTaubeCHaczkuAJoethamATakedaKDuezCet al. Effect of microbial heat shock proteins on airway inflammation and hyperresponsiveness. J Immunol. (2002) 169:53007. 10.4049/jimmunol.169.9.5300

  • 72.

    KolASukhovaGKLichtmanAHLibbyP. Chlamydial heat shock protein 60 localizes in human atheroma and regulates macrophage tumor necrosis factor-alpha and matrix metalloproteinase expression. Circulation. (1998) 98:3007. 10.1161/01.CIR.98.4.300

  • 73.

    YamaguchiHOsakiTTaguchiHHanawaTYamamotoTFukudaMet al. Growth inhibition of Helicobacter pylori by monoclonal antibody to heat-shock protein 60. Microbiol Immunol. (1997) 41:90916. 10.1111/j.1348-0421.1997.tb01949.x

  • 74.

    OjhaAAnandMBhattAKremerLJacobsWRJrHatfullGF. GroEL1: a dedicated chaperone involved in mycolic acid biosynthesis during biofilm formation in mycobacteria. Cell. (2005) 123:86173. 10.1016/j.cell.2005.09.012

  • 75.

    ReddiKMeghjiSNairSPArnettTRMillerADPreussMet al. The Escherichia coli chaperonin 60 (groEL) is a potent stimulator of osteoclast formation. J Bone Miner Res. (1998) 13:12606. 10.1359/jbmr.1998.13.8.1260

  • 76.

    BabadyNEPangYPElpelegOIsayaG. Cryptic proteolytic activity of dihydrolipoamide dehydrogenase. Proc Natl Acad Sci USA. (2007) 104:615863. 10.1073/pnas.0610618104

  • 77.

    KinnbyBBoothNASvensaterG. Plasminogen binding by oral streptococci from dental plaque and inflammatory lesions. Microbiology. (2008) 154(Pt 3):92431. 10.1099/mic.0.2007/013235-0

  • 78.

    ChuongSDMullenRTMuenchDG. Identification of a rice RNA- and microtubule-binding protein as the multifunctional protein, a peroxisomal enzyme involved in the beta -oxidation of fatty acids. J Biol Chem. (2002) 277:241929. 10.1074/jbc.M109510200

  • 79.

    GranatoDBergonzelliGEPridmoreRDMarvinLRouvetMCorthesy-TheulazIE. Cell surface-associated elongation factor Tu mediates the attachment of Lactobacillus johnsonii NCC533 (La1) to human intestinal cells and mucins. Infect Immun. (2004) 72:21609. 10.1128/IAI.72.4.2160-2169.2004

  • 80.

    DalloSFKannanTRBlaylockMWBasemanJB. Elongation factor Tu and E1 beta subunit of pyruvate dehydrogenase complex act as fibronectin binding proteins in Mycoplasma pneumoniae. Mol Microbiol. (2002) 46:104151. 10.1046/j.1365-2958.2002.03207.x

  • 81.

    KesimerMKilicNMehrotraRThorntonDJSheehanJK. Identification of salivary mucin MUC7 binding proteins from Streptococcus gordonii. BMC Microbiol. (2009) 9:163. 10.1186/1471-2180-9-163

  • 82.

    KunertALosseJGruszinCHuhnMKaendlerKMikkatSet al. Immune evasion of the human pathogen Pseudomonas aeruginosa: elongation factor Tuf is a factor H and plasminogen binding protein. J Immunol. (2007) 179:297988. 10.4049/jimmunol.179.5.2979

  • 83.

    SunYCarneiroNCloreAMMoroGLHabbenJELarkinsBA. Characterization of maize elongationfactor 1A and its relationship to protein quality in the endosperm. Plant Physiol. (1997) 115:11017. 10.1104/pp.115.3.1101

  • 84.

    AntikainenJKuparinenVLahteenmakiKKorhonenTK. Enolases from Gram-positive bacterial pathogens and commensal lactobacilli share functional similarity in virulence-associated traits. FEMS Immunol Med Microbiol. (2007) 51:52634. 10.1111/j.1574-695X.2007.00330.x

  • 85.

    CastaldoCVastanoVSicilianoRACandelaMViciMMuscarielloLet al. Surface displaced alfa-enolase of Lactobacillus plantarum is a fibronectin binding protein. Microb Cell Fact. (2009) 8:14. 10.1186/1475-2859-8-14

  • 86.

    EsgleasMLiYHancockMAHarelJDubreuilJDGottschalkM. Isolation and characterization of alpha-enolase, a novel fibronectin-binding protein from Streptococcus suis. Microbiology. (2008) 154(Pt. 9):2668–79. 10.1099/mic.0.2008/017145-0

  • 87.

    FlodenAMWattJABrissetteCA. Borrelia burgdorferi enolase is a surface-exposed plasminogen binding protein. PLoS ONE. (2011) 6:e27502. 10.1371/journal.pone.0027502

  • 88.

    JongAYChenSHStinsMFKimKSTuanTLHuangSH. Binding of Candida albicans enolase to plasmin(ogen) results in enhanced invasion of human brain microvascular endothelial cells. J Med Microbiol. (2003) 52:61522. 10.1099/jmm.0.05060-0

  • 89.

    CarneiroCRPostolENomizoRReisLFBrentaniRR. Identification of enolase as a laminin-binding protein on the surface of Staphylococcus aureus. Microbes Infect. (2004) 6:6048. 10.1016/j.micinf.2004.02.003

  • 90.

    ShaJErovaTEAlyeaRAWangSOlanoJPPancholiVet al. Surface-expressed enolase contributes to the pathogenesis of clinical isolate SSU of Aeromonas hydrophila. J Bacteriol. (2009) 191:3095107. 10.1128/JB.00005-09

  • 91.

    LorenzattoKRMonteiroKMParedesRPaludoGPda FonsecaMMGalantiNet al. Fructose- bisphosphate aldolase and enolase from Echinococcus granulosus: genes, expression patterns and protein interactions of two potential moonlighting proteins. Gene. (2012) 506:7684. 10.1016/j.gene.2012.06.046

  • 92.

    MundodiVKucknoorASAldereteJF. Immunogenic and plasminogen-binding surface-associated alpha-enolase of Trichomonas vaginalis. Infect Immun. (2008) 76:52331. 10.1128/IAI.01352-07

  • 93.

    HattoriTTakeiNMizunoYKatoKKohsakaS. Neurotrophic and neuroprotective effects of euron-specific enolase on cultured neurons from embryonic rat brain. Neurosci Res. (1995) 21:1918. 10.1016/0168-0102(94)00849-B

  • 94.

    DeckerBLWicknerWT. Enolase activates homotypic vacuole fusion and protein transport to the vacuole in yeast. J Biol Chem. (2006) 281:145238. 10.1074/jbc.M600911200

  • 95.

    BhowmickIPKumarNSharmaSCoppensIJaroriGK. Plasmodium falciparum enolase: stage-specific expression and sub-cellular localization. Malar J. (2009) 8:179. 10.1186/1475-2875-8-179

  • 96.

    EntelisNBrandinaIKamenskiPKrasheninnikovIAMartinRPTarassovI. A glycolytic enzyme, enolase, is recruited as a cofactor of tRNA targeting toward mitochondria in Saccharomyces cerevisiae. Genes Dev. (2006) 20:160920. 10.1101/gad.385706

  • 97.

    LiuXReigBNasrallahIMStoverPJ. Human cytoplasmic serine hydroxymethyltransferase is an mRNA binding protein. Biochemistry. (2000) 39:1152331. 10.1021/bi000665d

  • 98.

    WistowGPiatigorskyJ. Recruitment of enzymes as lens structural proteins. Science. (1987) 236:15546. 10.1126/science.3589669

  • 99.

    TomlinsonIPAlamNARowanAJBarclayEJaegerEEKelsellDet al. Germline mutations in FH predispose to dominantly inherited uterine fibroids, skin leiomyomata and papillary renal cell cancer. Nat Genet. (2002) 30:40610. 10.1038/ng849

  • 100.

    AlvarezRABlaylockMWBasemanJB. Surface localized glyceraldehyde-3-phosphate dehydrogenase of Mycoplasma genitalium binds mucin. Mol Microbiol. (2003) 48:141725. 10.1046/j.1365-2958.2003.03518.x

  • 101.

    GozalboDGil-NavarroIAzorinIRenau-PiquerasJMartinezJPGilML. The cell wall-associated glyceraldehyde-3-phosphate dehydrogenase of Candida albicans is also a fibronectin and laminin binding protein. Infect Immun. (1998) 66:20529.

  • 102.

    BergmannSRohdeMHammerschmidtS. Glyceraldehyde-3-phosphate dehydrogenase of Streptococcus pneumoniae is a surface-displayed plasminogen-binding protein. Infect Immun. (2004) 72:24169. 10.1128/IAI.72.4.2416-2419.2004

  • 103.

    EgeaLAguileraLGimenezRSorollaMAAguilarJBadiaJet al. Role of secreted glyceraldehyde-3- phosphate dehydrogenase in the infection mechanism of enterohemorrhagic and enteropathogenic Escherichia coli: interaction of the extracellular enzyme with human plasminogen and fibrinogen. Int J Biochem Cell Biol. (2007) 39:1190203. 10.1016/j.biocel.2007.03.008

  • 104.

    SriramGMartinezJAMcCabeERLiaoJCDippleKM. Single-gene disorders: what role could moonlighting enzymes play?Am J Hum Genet. (2005) 76:91124. 10.1086/430799

  • 105.

    PancholiVFischettiVA. A major surface protein on group A streptococci is a glyceraldehyde-3-phosphate- dehydrogenase with multiple binding activity. J Exp Med. (1992) 176:41526. 10.1084/jem.176.2.415

  • 106.

    JinHSongYPBoelGKocharJPancholiV. Group A streptococcal surface GAPDH, SDH, recognizes uPAR/CD87 as its receptor on the human pharyngeal cell and mediates bacterial adherence to host cells. J Mol Biol. (2005) 350:2741. 10.1016/j.jmb.2005.04.063

  • 107.

    SampathPMazumderBSeshadriVGerberCAChavatteLKinterMet al. Noncanonical function of glutamyl-prolyl-tRNA synthetase: gene-specific silencing of translation. Cell. (2004) 119:195208. 10.1016/j.cell.2004.09.030

  • 108.

    NagataHIwasakiMMaedaKKuboniwaMHashinoEToeMet al. Identification of the binding domain of Streptococcus oralis glyceraldehyde-3-phosphate dehydrogenase for Porphyromonas gingivalis major fimbriae. Infect Immun. (2009) 77:51308. 10.1128/IAI.00439-09

  • 109.

    ModunBMorrisseyJWilliamsP. The staphylococcal transferrin receptor: a glycolytic enzyme with novel functions. Trends Microbiol. (2000) 8:2317. 10.1016/S0966-842X(00)01728-5

  • 110.

    LuGTXieJRChenLHuJRAnSQSuHZet al. Glyceraldehyde-3-phosphate dehydrogenase of Xanthomonas campestris pv. campestris is required for extracellular polysaccharide production and full virulence. Microbiology. (2009) 155:160212. 10.1099/mic.0.023762-0

  • 111.

    FuldeMBernardo-GarciaNRohdeMNachtigallNFrankRPreissnerKTet al. Pneumococcal phosphoglycerate kinase interacts with plasminogen and its tissue activator. Thromb Haemost. (2014) 111:40116. 10.1160/TH13-05-0421

  • 112.

    LayAJJiangXMKiskerOFlynnEUnderwoodACondronRet al. Phosphoglycerate kinase acts in tumour angiogenesis as a disulphide reductase. Nature. (2000) 408:86973. 10.1038/35048596

  • 113.

    KennedyMCMende-MuellerLBlondinGABeinertH. Purification and characterization of cytosolic aconitase from beef liver and its relationship to the iron-responsive element binding protein. Proc Natl Acad Sci USA. (1992) 89:117304. 10.1073/pnas.89.24.11730

  • 114.

    PhilpottCCKlausnerRDRouaultTA. The bifunctional iron-responsive element binding protein/cytosolic aconitase: the role of active-site residues in ligand binding and regulation. Proc Natl Acad Sci USA. (1994) 91:73215. 10.1073/pnas

  • 115.

    BanerjeeSNandyalaAKRaviprasadPAhmedNHasnainSE. Iron-dependent RNA-binding activity of Mycobacterium tuberculosis aconitase. J Bacteriol. (2007) 189:404652. 10.1128/JB.00026-07

  • 116.

    ChenXJWangXButowRA. Yeast aconitase binds and provides metabolically coupled protection to mitochondrial DNA. Proc Natl Acad Sci USA. (2007) 104:1373843. 10.1073/pnas.0703078104

  • 117.

    ChenXJWangXKaufmanBAButowRA. Aconitase couples metabolic regulation to mitochondrial DNA maintenance. Science. (2005) 307:7147. 10.1126/science.1106391

  • 118.

    CommichauFMStulkeJ. Trigger enzymes: bifunctional proteins active in metabolism and in controlling gene expression. Mol Microbiol. (2008) 67:692702. 10.1111/j.1365-2958.2007.06071.x

  • 119.

    KainulainenVLoimarantaVPekkalaAEdelmanSAntikainenJKylvajaRet al. Glutamine synthetase and glucose-6-phosphate isomerase are adhesive moonlighting proteins of Lactobacillus crispatus released by epithelial cathelicidin LL-37. J Bacteriol. (2012) 194:250919. 10.1128/JB.06704-11

  • 120.

    JangHHLeeKOChiYHJungBGParkSKParkJHet al. Two enzymes in one; two yeast peroxiredoxins display oxidative stress-dependent switching from a peroxidase to a molecular chaperone function. Cell. (2004) 117:62535. 10.1016/j.cell.2004.05.002

  • 121.

    UrbanCXiongXSohnKSchroppelKBrunnerHRuppS. The moonlighting protein Tsa1p is implicated in oxidative stress response and in cell wall biogenesis in Candida albicans. Mol Microbiol. (2005) 57:131841. 10.1111/j.1365-2958.2005.04771.x

  • 122.

    ChenJWDodiaCFeinsteinSIJainMKFisherAB. 1-Cys peroxiredoxin, a bifunctional enzyme with glutathione peroxidase and phospholipase A2 activities. J Biol Chem. (2000) 275:284217. 10.1074/jbc.M005073200

  • 123.

    MoonHLeeBChoiGShinDPrasadDTLeeOet al. NDP kinase 2 interacts with two oxidative stress- activated MAPKs to regulate cellular redox state and enhances multiple stress tolerance in transgenic plants. Proc Natl Acad Sci USA. (2003) 100:35863. 10.1073/pnas.252641899

  • 124.

    BatemanOAPurkissAGvan MontfortRSlingsbyCGrahamCWistowG. Crystal structure of beta-crystallin: adaptation of a class 1 aldehyde dehydrogenase for a new role in the eye lens. Biochemistry. (2003) 42:434956. 10.1021/bi027367w

  • 125.

    FujiiYKimotoHIshikawaKWatanabeKYokotaYNakaiNet al. Taxon-specific zeta -crystallin in Japanese tree frog (Hyla japonica) lens. J Biol Chem. (2001) 276:281349. 10.1074/jbc.M102880200

  • 126.

    RaoPVZiglerJSJr. Zeta-crystallin from guinea pig lens is capable of functioning catalytically as an oxidoreductase. Arch Biochem Biophys. (1991) 284:1815. 10.1016/0003-9861(91)90281-M

  • 127.

    RaoPVKrishnaCMZiglerJSJr. Identification and characterization of the enzymatic activity of zeta- crystallin from guinea pig lens. a novel NADPH:quinone oxidoreductase. J Biol Chem. (1992) 267:96102.

  • 128.

    HuangQLRussellPStoneSHZiglerJSJr. Zeta-crystallin, a novel lens protein from the guinea pig. CurrEye Res. (1987) 6:72532. 10.3109/02713688709034836

  • 129.

    NoahCEMalikMBublitzDCCamenaresDSellatiTJBenachJLet al. GroEL and lipopolysaccharide from Francisella tularensis live vaccine strain synergistically activate human macrophages. Infect Immun. (2010) 78:1797806. 10.1128/IAI.01135-09

  • 130.

    CocklePJGordonSVLalvaniABuddleBMHewinsonRGVordermeierHM. Identification of novel Mycobacterium tuberculosis antigens with potential as diagnostic reagents or subunit vaccine candidates by comparative genomics. Infect Immun. (2002) 70:69967003. 10.1128/IAI.70.12.6996-7003.2002

  • 131.

    ErnstJDLewinsohnDMBeharSBlytheMSchlesingerLSKornfeldHet al. Meeting report: NIH workshop on the tuberculosis immune epitope database. Tuberculosis (Edinb). (2008) 88:36670. 10.1016/j.tube.2007.11.002

  • 132.

    LuMXiaZYBaoL. Enhancement of antimycobacterial Th1-cell responses by a Mycobacterium bovis BCG prime-protein boost vaccination strategy. Cell Immunol. (2013) 285:1117. 10.1016/j.cellimm.2013.10.001

  • 133.

    JiangYLiuHWangXXiaoSLiMLiGet al. Genetic diversity of immune-related antigens in region of difference 2 of Mycobacterium tuberculosis strains. Tuberculosis (Edinb). (2017) 104:17. 10.1016/j.tube.2016.05.002

  • 134.

    AndersenP. Host responses and antigens involved in protective immunity to Mycobacterium tuberculosis. Scand J Immunol. (1997) 45:11531. 10.1046/j.1365-3083.1997.d01-380.x

  • 135.

    MustafaTWikerHGMorkveOSvilandL. Reduced apoptosis and increased inflammatory cytokines in granulomas caused by tuberculous compared to non-tuberculous mycobacteria: role of MPT64 antigen in apoptosis and immune response. Clin Exp Immunol. (2007) 150:10513. 10.1111/j.1365-2249.2007.03476.x

  • 136.

    BackusKMDolanMABarryCSJoeMMcPhiePBoshoffHIet al. The three Mycobacterium tuberculosis antigen 85 isoforms have unique substrates and activities determined by non-active site regions. J Biol Chem. (2014) 289:2504153. 10.1074/jbc.M114.581579

  • 137.

    van ElsCACorbiereVSmitsKvanGaans-van den Brink JAPoelenMCMascartFet al. Toward Understanding the Essence of Post-Translational Modifications for the Mycobacterium tuberculosis immunoproteome. Front Immunol. (2014) 5:361. 10.3389/fimmu.2014.00361

  • 138.

    SmithGTSweredoskiMJHessS. O-linked glycosylation sites profiling in Mycobacterium tuberculosis culture filtrate proteins. J Proteomics. (2014) 97:296306. 10.1016/j.jprot.2013.05.011

  • 139.

    LauvauGPamerEG. CD8 T cell detection of bacterial infection: sniffing for formyl peptides derived from Mycobacterium tuberculosis. J Exp Med. (2001) 193:F3539. 10.1084/jem.193.10.F35

  • 140.

    DoiTYamadaHYajimaTWajjwalkuWHaraTYoshikaiY. H2-M3-restricted CD8+ T cells induced by peptide-pulsed dendritic cells confer protection against Mycobacterium tuberculosis. J Immunol. (2007) 178:380613. 10.4049/jimmunol.178.6.3806

  • 141.

    MirSASharmaS. Immunotherapeutic potential of N-formylated peptides of ESAT-6 and glutamine synthetase in experimental tuberculosis. Int Immunopharmacol. (2014) 18:298303. 10.1016/j.intimp.2013.09.010

  • 142.

    DobosKMKhooKHSwiderekKMBrennanPJBelisleJT. Definition of the full extent of glycosylation of the 45-kilodalton glycoprotein of Mycobacterium tuberculosis. J Bacteriol. (1996) 178:2498506. 10.1128/jb.178.9.2498-2506.1996

  • 143.

    HornCNamaneAPescherPRiviereMRomainFPuzoGet al. Decreased capacity of recombinant 45/47- kDa molecules (Apa) of Mycobacterium tuberculosis to stimulate T lymphocyte responses related to changes in their mannosylation pattern. J Biol Chem. (1999) 274:3202330.

  • 144.

    WikerHGHarboeM. The antigen 85 complex: a major secretion product of Mycobacterium tuberculosis. Microbiol Rev. (1992) 56:64861.

  • 145.

    MalenHSoftelandTWikerHG. Antigen analysis of Mycobacterium tuberculosis H37Rv culture filtrate proteins. Scand J Immunol. (2008) 67:24552. 10.1111/j.1365-3083.2007.02064.x

  • 146.

    WikerHGHarboeMNagaiSPatarroyoMERamirezCCruzN. MPB59, a widely cross-reacting protein of Mycobacterium bovis BCG. Int Arch Allergy Appl Immunol. (1986) 81:30714. 10.1159/000234154

  • 147.

    HarboeMWikerHGDuncanJRGarciaMMDukesTWBrooksBWet al. Protein G-based enzyme-linked immunosorbent assay for anti-MPB70 antibodies in bovine tuberculosis. J Clin Microbiol. (1990) 28:91321.

  • 148.

    BaldwinSLD'SouzaCDOrmeIMLiuMAHuygenKDenisOet al. Immunogenicity and protective efficacy of DNA vaccines encoding secreted and non-secreted forms of Mycobacterium tuberculosis Ag85A. Tuber Lung Dis. (1999) 79:2519. 10.1054/tuld.1998.0196

  • 149.

    BelisleJTVissaVDSievertTTakayamaKBrennanPJBesraGS. Role of the major antigen of Mycobacterium tuberculosis in cell wall biogenesis. Science. (1997) 276:14202. 10.1126/science.276.5317.1420

  • 150.

    KremerLMaughanWNWilsonRADoverLGBesraGS. The M. tuberculosis antigen 85 complex and mycolyltransferase activity. Lett Appl Microbiol. (2002) 34:2337. 10.1046/j.1472-765x.2002.01091.x

  • 151.

    GoinsCMDajnowiczSSmithMDParksJMRonningDR. Mycolyltransferase from Mycobacterium tuberculosis in covalent complex with tetrahydrolipstatin provides insights into antigen 85 catalysis. J Biol Chem. (2018) 293:365162. 10.1074/jbc.RA117.001681

  • 152.

    Abou-ZeidCRatliffTLWikerHGHarboeMBennedsenJRookGet al. Characterization of fibronectin-binding antigens released by Mycobacterium tuberculosis and Mycobacterium bovis BCG. Infect Immun. (1988) 56:304651.

  • 153.

    KuoCJPtakCPHsiehCLAkeyBLChangYF. Elastin, a novel extracellular matrix protein adhering to mycobacterial antigen 85 complex. J Biol Chem. (2013) 288:388696. 10.1074/jbc.M112.415679

  • 154.

    SwierzkoASBartlomiejczykMABrzostekALukasiewiczJMichalskiMDziadekJet al. Mycobacterial antigen 85 complex (Ag85) as a target for ficolins and mannose-binding lectin. Int J Med Microbiol. (2016) 306:21221. 10.1016/j.ijmm.2016.04.004

  • 155.

    PathanAASanderCRFletcherHAPoultonIAlderNCBeveridgeNEet al. Boosting BCG with recombinant modified vaccinia ankara expressing antigen 85A: different boosting intervals and implications for efficacy trials. PLoS ONE. (2007) 2:e1052. 10.1371/journal.pone.0001052

  • 156.

    KarbalaeiZadeh Babaki MSoleimanpourSRezaeeSA. Antigen 85 complex as a powerful Mycobacterium tuberculosis immunogene: biology, immune-pathogenicity, applications in diagnosis, and vaccine design. Microb Pathog. (2017) 112:209. 10.1016/j.micpath.2017.08.040

  • 157.

    BerksBC. A common export pathway for proteins binding complex redox cofactors?Mol Microbiol. (1996) 22:393404. 10.1046/j.1365-2958.1996.00114.x

  • 158.

    BendtsenJDJensenLJBlomNVon HeijneGBrunakS. Feature-based prediction of non-classical and leaderless protein secretion. Protein Eng Des Sel. (2004) 17:34956. 10.1093/protein/gzh037

  • 159.

    PalmerTSargentFBerksBC. The tat protein export pathway. EcoSal Plus. (2010) 4:135. 10.1128/ecosalplus.4.3.2

  • 160.

    MawuenyegaKGForstCVDobosKMBelisleJTChenJBradburyEMet al. Mycobacterium tuberculosis functional network analysis by global subcellular protein profiling. Mol Biol Cell. (2005) 16:396404. 10.1091/mbc.e04-04-0329

  • 161.

    BaghelASTandonRGuptaGKumarASharmaRKAggarwalNet al. Characterization of protein acyltransferase function of recombinant purified GlnA1 from Mycobacterium tuberculosis: a moon lighting property. Microbiol Res. (2011) 166:66272. 10.1016/j.micres.2011.02.001

  • 162.

    HarthGClemensDLHorwitzMA. Glutamine synthetase of Mycobacterium tuberculosis: extracellular release and characterization of its enzymatic activity. Proc Natl Acad Sci USA. (1994) 91:93426. 10.1073/pnas.91.20.9342

  • 163.

    HarthGMaslesa-GalicSTulliusMVHorwitzMA. All four Mycobacterium tuberculosis glnA genes encode glutamine synthetase activities but only GlnA1 is abundantly expressed and essential for bacterial homeostasis. Mol Microbiol. (2005) 58:115772. 10.1111/j.1365-2958.2005.04899.x

  • 164.

    HernandezSFrancoLCalvoAFerragutGHermosoAAmelaIet al. Bioinformatics and moonlighting proteins. Front Bioeng Biotechnol. (2015) 3:90. 10.3389/fbioe.2015.00090

  • 165.

    WongAGehringCIrvingHR. Conserved functional motifs and homology modeling to predict hidden moonlighting functional sites. Front Bioeng Biotechnol. (2015) 3:82. 10.3389/fbioe.2015.00082

  • 166.

    KhanIKKiharaD. Computational characterization of moonlighting proteins. Biochem Soc Trans. (2014) 42:17805. 10.1042/BST20140214

  • 167.

    ChaputMClaesVPortetelleDCludtsICravadorABurnyAet al. The neurotrophic factor neuroleukin is 90% homologous with phosphohexose isomerase. Nature. (1988) 332:4545. 10.1038/332454a0

  • 168.

    XuWSeiterKFeldmanEAhmedTChiaoJW. The differentiation and maturation mediator for human myeloid leukemia cells shares homology with neuroleukin or phosphoglucose isomerase. Blood. (1996) 87:45026.

  • 169.

    SuzukiCKRepMvan DijlJMSudaKGrivellLASchatzG. ATP-dependent proteases that alsochaperone protein biogenesis. Trends Biochem Sci. (1997) 22:11823. 10.1016/S0968-0004(97)01020-7

  • 170.

    SokerSTakashimaSMiaoHQNeufeldGKlagsbrunM. Neuropilin-1 is expressed by endothelial andtumor cells as an isoform-specific receptor for vascular endothelial growth factor. Cell. (1998) 92:73545. 10.1016/S0092-8674(00)81402-6

  • 171.

    StuttsMJCanessaCMOlsenJCHamrickMCohnJARossierBCet al. CFTR as a cAMP-dependent regulator of sodium channels. Science. (1995) 269:84750. 10.1126/science.7543698

  • 172.

    HendersonBMartinA. Bacterial virulence in the moonlight: multitasking bacterial moonlighting proteins are virulence determinants in infectious disease. Infect Immun. (2011) 79:347691. 10.1128/IAI.00179-11

  • 173.

    JefferyCJ. Moonlighting proteins: old proteins learning new tricks. Trends Genet. (2003) 19:4157. 10.1016/S0168-9525(03)00167-7

  • 174.

    JefferyCJ. Molecular mechanisms for multitasking: recent crystal structures of moonlighting proteins. Curr OpinStruct Biol. (2004) 14:6638. 10.1016/j.sbi.2004.10.001

  • 175.

    NobeliIFaviaADThorntonJM. Protein promiscuity and its implications for biotechnology. Nat Biotechnol. (2009) 27:15767. 10.1038/nbt1519

  • 176.

    ShamsFOldfieldNJWooldridgeKGTurnerDP. Fructose-1, 6-bisphosphate aldolase (FBA)-a conserved glycolytic enzyme with virulence functions in bacteria: ‘ill met by moonlight'. Biochem Soc Trans. (2014) 42:17925. 10.1042/BST20140203

  • 177.

    XolalpaWVallecilloAJLaraMMendoza-HernandezGCominiMSpallekRet al. Identification of novel bacterial plasminogen-binding proteins in the human pathogen Mycobacterium tuberculosis. Proteomics. (2007) 7:333241. 10.1002/pmic.200600876

  • 178.

    dela Paz Santangelo MGestPMGuerinMECoinconMPhamHRyanGet al. (2011). Glycolytic and non-glycolytic functions of Mycobacterium tuberculosis fructose-1, 6-bisphosphate aldolase, an essential enzyme produced by replicating and non-replicating bacilli. J Biol Chem.286:4021931. 10.1074/jbc.M111.259440

  • 179.

    ArnoldSLKentTHogarthCASchlattSPrasadBHaenischMet al. Importance of ALDH1A enzymes in determining human testicular retinoic acid concentrations. J Lipid Res. (2015) 56:34257. 10.1194/jlr.M054718

  • 180.

    ChristgenSLZhuWSanyalNBibiBTannerJJBeckerDF. Discovery of the membrane binding domain in trifunctional proline utilization A. Biochemistry. (2017) 56:6292303. 10.1021/acs.biochem.7b01008

Summary

Keywords

Mycobacterium bovis BCG Moreau, bi-dimensional electrophoresis, MALDI-TOF-TOF, surface-associated proteins, moonlighting proteins

Citation

Pagani TD, Guimarães ACR, Waghabi MC, Corrêa PR, Kalume DE, Berrêdo-Pinho M, Degrave WM and Mendonça-Lima L (2019) Exploring the Potential Role of Moonlighting Function of the Surface-Associated Proteins From Mycobacterium bovis BCG Moreau and Pasteur by Comparative Proteomic. Front. Immunol. 10:716. doi: 10.3389/fimmu.2019.00716

Received

11 December 2018

Accepted

18 March 2019

Published

26 April 2019

Volume

10 - 2019

Edited by

Juraj Ivanyi, King's College London, United Kingdom

Reviewed by

Patrick Brennan, Colorado State University, United States; Mahavir Singh, Lionex GmbH, Germany

Updates

Copyright

*Correspondence: Talita Duarte Pagani Leila Mendonça-Lima

†Present Address: Talita Duarte Pagani, Laboratório de Genômica Estrutural, Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

This article was submitted to Vaccines and Molecular Therapeutics, a section of the journal Frontiers in Immunology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics