Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, February 2005, p. 749-759, Vol. 49, No. 2
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.2.749-759.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Bayer HealthCare AG, Pharma Research, Wuppertal, Germany,1 Genedata, Basel, Switzerland2
Received 4 June 2004/ Returned for modification 16 August 2004/ Accepted 3 October 2004
|
|
|---|
|
|
|---|
However, as antibiotics are supposed to inhibit bacterial growth, they also change the expression of a wide range of genes not directly linked to the target's function, which in many cases obscures the primary antibiotic effect. Therefore, in order to deduce characteristic signatures for antibiotics, comparative global expression analyses need to be performed on the basis of a large collection of diverse mRNA profiles that represent different cellular stress states. A so-called reference compendium of expression profiles makes it possible to classify the MOA of antibiotics and to pinpoint the typically very few compound- or MOA-specific marker genes. Recently, Hutter et al. (30) used a reference compendium based on the gram-positive model bacterium Bacillus subtilis to classify the MOAs of antibacterial compounds by using support vector machine algorithms. Complementary to this work, Fischer et al. (16) proposed an integrated genome expression data analysis strategy to identify and characterize antibiotic-responsive regulatory networks.
The present mRNA reference compendium approaches to antibiotic MOA characterization suffer from the limitation that only mechanisms that are employed by already existing antibiotics can be predicted. The expansion of a reference compendium with expression profiles of bacterial mutants underexpressing potential target proteins could therefore be a helpful means to overcome this hurdle in identifying novel antibacterial MOAs. In this study, we generated a compendium of genome-wide expression profiles, representing the response of B. subtilis to treatment with 14 chemically diverse antibiotics. The data representing chemically induced mRNA profiles were supplemented by conditional mutant profiles, representing the effect of underexpressing four distinct essential genes or operons encoding antibacterial target proteins. Here, we focus on the target areas of aminoacyl-tRNA synthesis (AATS) and fatty acid biosynthesis (FAB). In mutant ILES() and mutant FABF(), genes were downregulated which encode two classical antibacterial targets: the isoleucyl-tRNA synthetase and the type II ß-ketoacyl-acyl carrier protein synthetase. Similarly, in mutants PHEST() and ACCDA() the amounts of transcripts were reduced; these transcripts code for two emerging targets: the phenylalanyl-tRNA synthetase (PheST) and the carboxyltransferase subunits (AccDA) of the acetyl coenzyme A (acetyl-CoA) carboxylase (ACC).
First, we show that downregulation of target genes causes transcriptional responses similar to those generated by respective antibiotic treatments as demonstrated with the well-investigated cerulenin (known to target FabF in B. subtilis) and mupirocin (known to inhibit IleS) (29, 34, 40). Second, we demonstrate that a reference compendium including underexpression mutant profiles enables automated MOA classification, by using mRNA profile classification algorithms. Third, we show that prediction of the molecular target of novel antibacterials is possible with such a combined chemical and mutant data compendium. By applying our analysis strategy, we are able to independently support the recently published PheST inhibition mechanism of a novel phenyl-thiazolylurea (PTU) derivative (5). Finally, we find that one of our mutant profiles matches the mRNA profile induced by the natural product moiramide B (MOIB) (42). Based on our MOA analysis techniques, the mRNA profiles triggered by this compound suggest that the antibiotic activity of MOIB is most likely due to inhibition of ACC, representing a completely novel antibacterial mechanism.
Our MOA analysis strategy consists of three major levels. First, we automatically classify compounds and underexpression effects based on a reference compendium of expression profiles, using an algorithm based on a nonparametric statistical test. Second, we analyze the expression patterns of major regulons and biosynthetic pathways by investigating the statistical overrepresentation of key pathways. Third, we focus on the detailed expression activity of individual pathways to identify the actual molecular target of the compounds under investigation.
|
|
|---|
Four genes or operons (ileS, fabF, pheST, and accDA) were each put under the control of a xylose-inducible promoter, resulting in strains called B. subtilis ILES(), FABF(), PHEST(), and ACCDA(), respectively. The first three strains were generated as follows. The 5' ends of the genes were amplified from genomic B. subtilis DNA (Table 1) and cloned into the vector pJH101 (15) via the BamHI restriction site. B. subtilis was transformed with the resulting constructs and selected on Luria-Bertani (LB) medium with 0.25% (wt/vol) xylose and 3 µg of chloramphenicol/ml, the marker for pJH101 integration (single recombination). While the genes of interest were all located downstream of the xylose-dependent promoter, 5' fragments of the respective genes remained under the control of the wild-type promoters.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotides used in this study
|
Growth experiments and RNA isolation. Cultures of B. subtilis 168 were grown in Belitzky minimal medium (58) to an optical density at 600 nm (OD600) of 0.4 at 37°C. After removal of an aliquot (time zero, untreated sample), the remaining culture was subjected to compound stress by addition of antibiotics (Table 2) and incubated for a further 10 and 40 min to obtain aliquots of treated cells (time points 1 and 2, respectively). Empirically, we found that the applied compound concentrations should be as high as possible but should not cause more than 25% growth reduction compared to the control at time point 2. Cell harvesting and RNA isolation were performed as described recently (16).
|
View this table: [in a new window] |
TABLE 2. Compounds and their concentrations used for transcriptome analysis
|
cDNA labeling and microarray experiments. Generation of fluorescence-labeled cDNA, hybridization with B. subtilis whole-genome arrays (Eurogentec), and image analysis with the Axon GenePix 400A confocal laser scanner (Axon Instruments) were performed as described previously (16).
Expression data analysis. The quality of the expression profiling experiments was checked as outlined previously (16) by using the Expressionist 4.1 software (Genedata, Basel, Switzerland). The quality-checked signals of all features were used for a global microarray normalization, with the scaling factor being calculated with the requirement that the sum of all logarithmized values in the Cy3 channel equals the sum of all logarithmized values in the Cy5 channel. For each condition, duplicates or triplicates of microarray experiments were performed, whereas for each gene the average fold factor was calculated on the basis of logarithms of the signal ratios.
For an estimate of the inherent biological variability of our B. subtilis cultures, we treated three independently grown B. subtilis cultures with cerulenin (16). We then calculated for this representative experiment set the standard deviation of the fold factor as 0.164. This means that a fold factor of 1.4 (or 1/1.4) already represents 3 standard deviations of the Gaussian distribution function (equivalent to a P value of >99.73%).
Automated MOA classification. The averaged expression fold factors for the 40-min compound exposure experiments together with the mutant data were used for MOA classification. The a priori assignments of the reference compounds and underexpression effects to categorical MOA classes were performed according to the literature and present biological knowledge (see Results). A consistency check of the precategorized compound-mRNA profiles was performed (cross-validation) by removing individual compound profiles or pairs of compound- and mutant-based profiles (see Results) one after another from the reference set and reclassifying the profile(s) based on the remaining (N-1 or N-2) reference profiles.
Cross-validation as well as classification of novel profiles was performed using the Wilcoxon classifier implemented in Expressionist (Genedata). In brief, this statistical method needs a "query expression profile" as input (e.g., representing an expression profile triggered by a novel compound). In a first step, correlation coefficients of this profile are calculated with regard to all reference experiments in the reference compendium database. Ranking the reference experiments by their corresponding correlation coefficients leads to a sorted list reflecting their global similarity to the query gene expression profile. In an ideal case, one would expect that all the top experiments belong to one MOA class, with the other MOA classes being distributed in the lower parts of the rank list. In order to define a quantitative measure for the rank separation of the gene expression correlations and their associated a priori MOA classes, a Wilcoxon test was employed. Here, the Wilcoxon test was used to check for the separation of two classes, which are the correlation coefficients belonging to a given MOA class compared to the correlations of all other MOA classes. A one-sided Wilcoxon test probes how well these two groups separate in terms of their ranking distribution. In this way, all MOA categories represented in our reference compendium were investigated. Applying a one-sided Wilcoxon test enables us to calculate P values for all query experiments and all represented MOAs. The so-called "affinity values" (A) are then calculated as the inverse error function of the Wilcoxon test's P value, A = erf1(pWilcoxon), reflecting the MOA assignment affinity.
Gene expression correlation analyses and promoter prediction. Correlation analyses based on gene expression data were performed using Expressionist 4.1 (Genedata). Operon and promoter prediction was performed using Phylosopher 4.0 (Genedata). The algorithms have been described previously (16).
Differential MIC determination. Microdilution MICs were determined against B. subtilis strains in 96-well microtiter plates in LB medium containing a 0 or 0.25% xylose concentration and serial dilutions (twofold) of antibiotics. The precultures of the conditional mutants ILES(), FABF(), and PHEST() additionally contained 3 µg of chloramphenicol/ml. Various starting inocula were used in the range of 104 to 106 CFU ml1 derived from precultures grown for 3 h in LB medium with 0 or 0.25% (wt/vol) xylose. The MIC was the lowest concentration of drug that yielded no visible growth after incubation for 18 to 24 h at 37°C. Endpoints were determined by measuring the OD600 with the microtiter plate reader EL312e (Bio-Tec Instruments).
qRT-PCR.
For selected genes, measurements of relative transcript amounts were performed by quantitative reverse transcription-PCR (qRT-PCR) with the QuantiTect SYBR Green RT-PCR kit (Qiagen) according to the manufacturer's instructions. The same RNA samples were used for cDNA synthesis and subsequent quantitative PCR as for microarray experiments. Specific primer pairs were selected for the genes accA, fabF, pheS, pheT, ileS, fabHB, and yurP (Table 1). In addition, a piece of the 16S rRNA was quantified, in order to normalize the obtained expression values of selected genes to the amount of 16S rRNA in each RNA sample (Table 1). The reactions were performed on an ABI Prism 7700 Sequence Detection system (PE Applied Biosystems), during which the fluorescence signal from SYBR Green intercalation was monitored to quantify the double-stranded DNA product formed after each PCR cycle. The difference (n = fold) in the initial concentration of each transcript (normalized to 16S rRNA) with respect to the wild type was calculated according to the comparative Ct method according to the equation 2
Ct (ABI user bulletin).
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Relative expression levels of genes in conditional mutants
|
|
View this table: [in a new window] |
TABLE 4. Differential sensitivity of conditional mutants derived from B. subtilis 168
|
|
View this table: [in a new window] |
TABLE 5. MOA classification of antibiotics and underexpression effects in conditional mutants based on whole-genome expression profilesa
|
Antibiotic-induced regulon and pathway activation pattern analysis. The transcriptome data used for MOA classification were analyzed with a focus on the deregulation of major regulons and biosynthetic pathways. The genes representing 16 regulatory units and 12 (generalized) pathways or groups of pathways from B. subtilis were compiled based on the whole-genome annotation from Phylosopher (Genedata) and on the functional classifications derived from the COG database at the National Center for Biotechnology Information (60) (see http://www.genedata.com/publications/supplements/). Additionally, we compiled and incorporated information from selected publications. In this work, the following regulons were considered: relA-dependent stringent response (13); fapR-dependent FAB and phospholipid biosynthesis genes (16, 52); lexA-dependent SOS response (41); the hrcA- and ctsR-dependent class I and III heat shock regulons (53); and alternative RNA polymerase sigma factor-dependent regulons sigB (21, 46, 47), sigD (39), sigH (8), sigF (1, 14, 25, 45, 54), sigE (31, 32, 35, 45, 49, 59, 64), sigG (31, 45), sigK (18, 31, 44), sigW (10, 11, 48, 62), and sigX (9, 11, 28, 48). The not-yet-completely characterized sigM regulon (9, 11, 27, 61) was extended by additional candidate genes by a combined gene expression-genome sequence analysis. For this, we focused on the gene ypuA, a gene of unknown function that was among the genes most strongly induced by cell wall biosynthesis inhibitors. We ran correlation analyses against the whole reference compendium transcriptome data set. Prediction of operons composed of correlated genes, which are adjacent on the B. subtilis chromosome, was done using Phylosopher (Genedata). The 400-bp upstream regions of the 20 genes and operons, which represented the most highly correlated transcription units with regard to the ypuA profile, were searched for sigma factor binding sites similar to the one predicted upstream of ypuA (5'-TGAAAC-N17-19-GTCT-3'), allowing for maximally two mismatches. Indeed, all top-20 genes and operons possess such binding sites, among these the sigma factor gene sigM itself as well as six additional genes and operons, for which sigM-dependent expression has been previously published (2) (see http://www.genedata.com/publications/supplements/).
In order to get a comprehensive overview of the induced and repressed regulons and pathways, we defined cutoff values for considering genes induced (expression values of
1.8) or repressed (expression values of
1/1.8). When at least 20% of genes belonging to a regulon or pathway were induced or repressed according to the definition mentioned above, the respective regulon or pathway was defined to be "deregulated." Figure 1A gives a comprehensive overview of the B. subtilis global response to a wide spectrum of antibiotics with different MOAs and a comparison of these effects to the corresponding mutant profiles. The detailed view of the pathway activation patterns revealed that FAB and AATS inhibition indeed exhibited striking similarities in their expression profiles. The majority of genes induced or repressed by mupirocin and PTU belong to the relA-dependent stringent response network. This result is consistent with recent findings for the response of B. subtilis to norvaline, another AATS inhibitor (13). The most strongly induced gene cluster in our IleS-PheST inhibition and downregulation experiments was the yur cluster, including the putative operon yurPONML (Table 3). These genes are predicted to encode a sugar ABC transport system and sugar-metabolizing enzymes. Among the repressed genes, not only the stringently controlled genes could be found (especially the ones encoding the translation machinery), but also many genes involved in amino acid and nucleotide biosynthesis. Such a broad inhibition profile across different major biosynthetic pathways together with the stringent response could not be observed with antibiotics employing MOAs other than FAB and AATS. This explains the Wilcoxon classifier-based categorization of the FAB and AATS inhibitors into the same MOA category.
![]() View larger version (47K): [in a new window] |
FIG.1. Antibiotic-induced regulon and pathway activation pattern analysis of B. subtilis. The percentages of transcriptionally induced or repressed genes per biosynthetic pathway or regulon are represented in color codes (see color legend at the bottom of the figure). A deregulation of at least 20% of the genes of a pathway-regulon is required for a compound or underexpression effect to be indicated in the figure. Compound and mutant abbreviations can be found in Tables 2 and 3. The pathway and regulon assignments of the genes can be found online at http://www.genedata.com/publications/supplements/. For clarity, the most striking pathway activation patterns have been framed (black). (A) The compound-triggered mRNA profiles all represent long-time exposure experiments (t, 40-min treatments). The compounds and mutants are sorted according to their MOA categories as known from the literature and as concluded from our transcriptome studies (abbreviations for MOA classes are in Table 5). (B) Time-dependent antibiotic-induced regulon and pathway activation pattern analysis of B. subtilis. The numbers after the compound abbreviations symbolize compound treatment times (1 = 10 min; 2 = 40 min). The compounds and mutants are sorted according to their target area assignment (see text). # indicates that, while a significant proportion of the aromatic biosynthesis genes are repressed under condition PTU1 (10-min treatment with the PheST inhibitor), the phenylalanine biosynthesis operon pheAB is slightly induced, which is not the case in any other measured stress condition.
|
AATS target area analysis: target identification for compound PTU. While initially focusing on the global cellular responses, we then tried to utilize the large-scale expression data to identify the actual molecular target of a given antibiotic. For this purpose we also considered the 10-min data in addition to the 40-min compound treatment experiments, so that time-dependent changes in gene expression could be observed (Fig. 1B). Here, we wanted to exemplify this approach for the recently characterized MOA of PTU as a proof-of-concept case study.
Strikingly, the isoleucyl-tRNA synthetase inhibition by mupirocin led to induction of the target gene ileS (approximately four- to sixfold) but not of any other tRNA synthetase gene. On the other hand, PTU induced only pheS and pheT among all 26 tRNA synthetase genes (approximately twofold). No other antibiotic effect, including FAB inhibition, led to induction of any tRNA synthetase genes (exceptions: the threonyl-tRNA synthetase genes thrZ and thrS). In addition, the isoleucine-leucine biosynthesis genes were induced by mupirocin already after 10 min but not after 40 min (Fig. 1B). Consistent with these results, the ileS underexpression mutant showed slight induction of this amino acid biosynthesis pathway. Significantly fewer genes of this pathway were induced by PheST and FAB inhibition; the other tested antibiotics, independently of their underlying MOA, did not induce those genes at all. Correspondingly, the phenylalanine biosynthesis operon pheAB was slightly induced by PTU after 10 min (
1.6- to 1.9-fold) but not after 40 min of treatment. In summary, these results suggested that PheST is the molecular target of PTU. Indeed, PTU is a substance that has been identified in enzymatic high-throughput screening. By using complementary experimental approaches, only recently could it be independently shown that PTU kills bacteria by inhibition of PheST (5).
FAB target area analysis: target identification for MOIB. Next, we sought to identify the molecular target of the so-far-uncharacterized natural product MOIB. For this, we applied an approach similar to that used for PTU but also included a detailed pathway analysis. First, we mapped the mRNA expression fold factors onto the FAB-phospholipid biosynthesis pathway. This analysis revealed the sections of the pathway that were transcriptionally deregulated (Fig. 2). The fabF downregulation by a factor of 10 in mutant FABF() led to strong induction of fabHB (>10-fold; Table 3 and Fig. 2) and induction of fabHA, plsC, fabD, fabG, and fabI (>1.8-fold). All these genes are known to be transcriptionally controlled by the recently identified regulator FapR (16, 52). Exposure of B. subtilis to the FabF-targeting cerulenin also led to induction of these genes (Fig. 2). By contrast, the FapR-regulated genes are not induced by AATS inhibitors and other tested antibiotics (Fig. 1). They are partially repressed by non-FAB inhibitors such as, e.g., PTU, gentamicin, and chloramphenicol. Remarkably, MOIB represses all measurable FapR-regulated genes already after 10 min (Fig. 2). It is known that reduction of the cellular malonyl-CoA concentration leads to repression of FapR-regulated genes (52). The findings that the global transcriptional response of MOIB resembled more strongly a profile characteristic of FAB inhibitors than one of AATS inhibitors (Fig. 1A) and that the FapR-regulated genes were all repressed instead of being induced suggested that the B. subtilis malonyl-CoA pool might be reduced by this antibiotic. As the key enzyme responsible for malonyl-CoA synthesis is the ACC, we hypothesized that ACC inhibition might be the reason for the observed characteristic mRNA profiles.
![]() View larger version (47K): [in a new window] |
FIG. 2. Transcriptional activation pattern of the B. subtilis fatty acid and phospholipid biosynthesis pathway for antibiotic treatment and the modified gene expression due to target underexpression. (A) Mutant with downregulated fabF; (B) treatment of B. subtilis 168 with cerulenin for 10 min; (C) mutant with downregulated accDA; (D) treatment of B. subtilis 168 with MOIB for 10 min (Phylosopher; Genedata). FabF repression and cerulenin treatment mainly induce genes encoding enzymes downstream of the metabolite malonyl-CoA. In contrast to this finding, treatment with MOIB and accDA repression mainly repress those genes. This pathway analysis indicated that cerulenin targets FabF, while MOIB possibly inhibits the ACC. Upregulated genes are colored in red, while downregulated genes are represented by green boxes; see color scale at left. Nonresponsive genes are colored in gray (metabolite names: R, acetyl, isovaleryl, 2-methylbuturyl, or isobutyryl).
|
10-fold downregulation of these two carboxyltransferase subunits of ACC, all FapR-dependent genes were strongly repressed (Fig. 2), as in the case of the MOIB treatment experiments. This reasoning suggested that MOIB represents an antibacterial agent acting via inhibition of ACC. In the meantime, the ACC has been independently confirmed to be the target of MOIB, by using complementary experimental techniques (17). |
|
|---|
4,100 B. subtilis genes. Thus, the large data volumes generated by microarray technologies require systematic data analysis techniques to characterize the MOA of antibiotics. Previously, it was shown that simple expression clustering or standard correlation analyses performed only poorly for compound MOA analyses (6, 23, 30, 50, 57). More appropriate algorithms that have been suggested for MOA classification purposes are support vector machines, k-nearest neighbor analysis, sparse linear discriminant analysis, and Fisher linear discriminant analysis (65). The results presented here are based on the Wilcoxon classifier, as implemented in the Expressionist system (Genedata). A cross-validation analysis showed that this classifier worked very well with the data produced in this study. A comparison with predictions based on other classification methods gave an overall consistent picture (data to be published elsewhere). The application of classification tools is the first step in the interpretation of complex transcriptome data in drug discovery. As each expression value is coupled to individual genes embedded in known functional and regulatory networks, the data can be alternatively interpreted from the physiological point of view. Therefore, we also investigated the transcriptional activity of the major regulatory units and biosynthetic pathways of B. subtilis. This helped to biologically interpret the computer-based MOA classification results. For instance, the stringent response regulon represents a global regulatory unit whose deregulation was indicative of AATS as well as FAB inhibition, leading to the classification of FAB and AATS inhibitors into one category (55).
In this study, we investigated the value of conditional mutants for the detection of completely novel MOAs. Applying classification algorithms together with regulon and pathway analyses, we were able to show that the directed downregulation of genes indeed led to transcriptional responses being equivalent to chemical inhibition of the corresponding gene products. We demonstrated this by measuring the mRNA profiles of mutants whose target genes were downregulated by xylose depletion. For instance, the FABF() mutant showed an induction of the genes regulated by the fatty acid-phospholipid biosynthesis regulator FapR, while accDA downregulation led to repression of those genes. The ileS underexpression led to an induction of isoleucine biosynthesis genes, which is consistent with the view that the cells sense starvation of amino acid biosynthesis in the case of inhibition of the respective AATS. Only the PHEST() mutant did not reveal deregulation of genes directly linked to the repressed target operon. One reason for this observation might be the very strong repression of the pheST operon (
50-fold), which already led to a strong global stress response including activation of sporulation genes (Fig. 1B). In addition, all four mutants showed secondary responses equivalent to the effects of the respective antibiotics, including the stringent response as well as induction of sigM-sigX-dependent genes in the case of fabF and accDA downregulation.
Briefly put, the primary effects affecting the target area and the broad secondary responses appear to be superposed in the expression profiles of the underexpression mutants. In contrast, the antibiotic treatment experiments revealed that the primary effects on the target area were often seen already after 10 min but were no longer detectable after 40 min (e.g., for PTU and MOIB), while the global secondary responses were more pronounced after 40 min. However, the data analysis techniques used in this work were still able to detect the underlying MOAs of the compounds, as was demonstrated by classifying the late-stage stress response profiles. The difference between the continuous underexpression of essential target genes and the inhibition of cell growth by antibiotic treatment for a defined time frame might explain the observed concurrence or chronology of primary and secondary responses. Nevertheless, we could show that the global characteristics of underexpression mutant profiles are equivalent to antibiotic treatment profiles, enabling the classification of antibiotics employing novel MOAs, as was shown for PTU and MOIB.
Comprehensive cross-validation studies showed that automated MOA classification of the 40-min mRNA profiles consistently resulted in superior classification rates compared to the cross-validation of the 10-min mRNA profiles (data not shown). Apparently, the higher number of deregulated genes after 40 min than at the earlier time points provides a more reliable basis for MOA classification purposes. In fact, the global MOA classification approach was an essential prerequisite to identify, for instance, the FAB pathway as the primary target area of MOIB. Once the broad target area had been pinpointed, early time point experiments were essential for studying the primary antibiotic effects and predicting the molecular targets such as ACC in the case of MOIB.
Transcriptome analyses as well as differential sensitivity tests with underexpression systems always represent indirect methods to prove inhibition effects of antibiotics. Therefore, it is difficult to decide whether other gene products, the expression of which is strictly coupled to the predicted target (here, e.g., accDA), might be the actual molecular target of the investigated antibiotic (here, e.g., MOIB). In the context of MOIB, it is worthwhile to mention that underexpression of the accBC operon leads to repression of the two FapR-dependent operons fabHAF and fabHB (52). The AccBC subunits complete the ACC complex, consisting of the carboxyltransferase subunits AccDA, the biotin carboxylase AccC, and the biotin-carboxyl carrier protein AccB. We therefore cannot discriminate which part of the ACC complex might be targeted by MOIB. In addition, blockage of the biotin synthesis or biotinylation of ACC might cause effects on the transcriptome similar to the ones caused by accDA downregulation. Although this issue has not been addressed in this study, it would be straightforward to apply our proposed analysis strategy to conditional mutants regulating the upstream biotin pathway. Nevertheless, our reference compendium approach together with the differential sensitivity test with the ACCDA() mutant gave consistent results, which narrow down the target area of MOIB to the start of the FAB pathway, with ACC as the most likely target candidate. To our knowledge, no antibiotic that targets ACC has been described yet. This means that our suggested approach with an extended mutant reference compendium can indeed be used to predict novel mechanisms of antibiotics, as demonstrated by MOIB's ACC inhibition.
In conclusion, we show that a combined statistical data analysis of a compendium of diverse global mRNA profiles, together with a thorough pathway and regulon analysis, can pinpoint the molecular target of antibiotics. We systematically analyzed the transcriptional activity of 28 major B. subtilis regulons and groups of biosynthetic pathways. Consideration of more regulatory units will most likely improve the MOA prediction accuracy and further our understanding of the response of bacteria to antibiotic stress. Obviously, addressing novel MOAs is very important, as many promising target areas are not exploited by the currently available antibiotics and can therefore not be used for a reference compendium approach. Thus, the inclusion of conditional mutants in antibiotic-based reference compendia significantly widens existing opportunities for these technologies. We were able to show that, whenever reference antibiotics are not available, conditional mutants can be used to mimic chemical inhibition of the corresponding gene products. The confirmation of the MOA of PTU as a PheST inhibitor and the prediction of the ACC as the target of MOIB validate the usefulness of such mutants for MOA predictions of uncharacterized bioactive compounds.
|
|
|---|
F of Bacillus subtilis by using random-sequence oligonucleotides. J. Bacteriol. 183:3623-3630.
W regulon: a comparative analysis of promoter consensus search, run-off transcription/macroarray analysis (ROMA), and transcriptional profiling approaches. J. Mol. Biol. 316:443-457.[CrossRef][Medline]
W and
M regulons. Mol. Microbiol. 45:1267-1276.[CrossRef][Medline]
B regulons. Mol. Microbiol. 49:1477-1491.[CrossRef][Medline]
D-dependent gene expression. J. Bacteriol. 182:3055-3062.
X and
W. J. Bacteriol. 183:1921-1927.
W regulon. Mol. Microbiol. 41:59-71.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»