- Research article
- Open Access
Predicting transcriptional regulatory interactions with artificial neural networks applied to E. coli multidrug resistance efflux pumps
BMC Microbiology volume 8, Article number: 101 (2008)
Little is known about bacterial transcriptional regulatory networks (TRNs). In Escherichia coli, which is the organism with the largest wet-lab validated TRN, its set of interactions involves only ~50% of the repertoire of transcription factors currently known, and ~25% of its genes. Of those, only a small proportion describes the regulation of processes that are clinically relevant, such as drug resistance mechanisms.
We designed feed-forward (FF) and bi-fan (BF) motif predictors for E. coli using multi-layer perceptron artificial neural networks (ANNs). The motif predictors were trained using a large dataset of gene expression data; the collection of motifs was extracted from the E. coli TRN. Each network motif was mapped to a vector of correlations which were computed using the gene expression profile of the elements in the motif. Thus, by combining network structural information with transcriptome data, FF and BF predictors were able to classify with a high precision of 83% and 96%, respectively, and with a high recall of 86% and 97%, respectively. These results were found when motifs were represented using different types of correlations together, i.e., Pearson, Spearman, Kendall, and partial correlation. We then applied the best predictors to hypothesize new regulations for 16 operons involved with multidrug resistance (MDR) efflux pumps, which are considered as a major bacterial mechanism to fight antimicrobial agents. As a result, the motif predictors assigned new transcription factors for these MDR proteins, turning them into high-quality candidates to be experimentally tested.
The motif predictors presented herein can be used to identify novel regulatory interactions by using microarray data. The presentation of an example motif to predictors will make them categorize whether or not the example motif is a BF, or whether or not it is an FF. This approach is useful to find new "pieces" of the TRN, when inspecting the regulation of a small set of operons. Furthermore, it shows that correlations of expression data can be used to discriminate between elements that are arranged in structural motifs and those in random sets of transcripts.
Unraveling transcriptional regulatory systems is a key step in understanding the regulation of bacterial biological processes, as a whole. In bacteria, transcription factors (TFs) dictate regulation in a great extent because they are directly involved in environmentally-triggered internal and external signals. This is well illustrated by the fact that 3/4 of the transcription factors (TFs) in E. coli have been described as exogenous signal sensing ; i.e., they are promptly activated/repressed by changes in extracellular conditions. Nevertheless, our present knowledge of the E. coli TRN is limited, as only half of the regulators currently known, and only a quarter of the total number of genes, are included in the TRN . Therefore, the use of computational methods can help revealing additional TRNs and contribute to a better understanding of bacterial systems behavior.
As a consequence, several network inference techniques have been introduced to predict regulatory interactions. Initially, motif discovery algorithms and comparative genomic approaches have been widely applied to the prediction of regulatory interactions by identifying binding sites in DNA sequences. Lately, with the increasing use of high-throughput processes for obtaining gene expression data, protein levels, and DNA-protein interactions, network inference techniques started to make use of these data. Identification of causal or dependence relationships in source data rely mostly on static/dynamic Bayesian networks [3–6], graphical Gaussian models [7–9], and relevance networks [10–12].
In general, these graphical models are used to learn causal relationships among genes or proteins, depending on the data available. The results of the analysis is a (un)directed network, in which the links represent regulatory interactions among the entities of the network. Usually, these methods do not use the available information of the organism's TRN as a prior input for the inference procedure. Instead, these methods recover a whole network from scratch, based only on data.
We propose a different view to tackle the problem of inferring regulatory interactions from high-throughput data, which uses pattern identification rather than causal analysis. Over the years, high-quality regulatory interaction data has been collected for E. coli. By using a pattern identification technique, this accumulated knowledge can be easily encoded as learning examples, and then a trained model can infer new hypotheses.
As previous studies have shown , E. coli's TRN is made up of building blocks known as structural network motifs. Network motifs are clusters of transcription factors and target genes that are linked together in a specific pattern that not appear in random networks. In E. coli, the FF and BF arrangement patterns were identified to be these discrete units. These blocks of the network reveal the mechanisms of regulation employed by the bacterium. In the FF, a global TF X regulates a function-specific TF Y, and both cooperatively regulate an operon. In the BF, TFs X and Y regulate a pair of operons involved in the same biological process. Also, several motifs bundle together and distribute in functional modules that were found to control several cellular processes [2, 14].
We present herein a novel approach for inferring regulatory interactions from transcriptome data tailored for E. coli. The method uses the information encoded in the E. coli' s TRN in the form of motifs, and benefits from the fact that experimentally validated E. coli's network is larger than other organisms' TRNs. In this study we show that in the TRN, FF and BF motifs exhibit a joint-expression pattern that differs from a random collection of regulators and targets.
Having this in mind, we designed ANNs classifiers to infer network motifs in E. coli. Those ANNs use the information learned from the collection of FF and BF motifs that appeared in the TRN, as well as from gene expression data, to suggest new motifs. Thus, as FF and BF predictors, these ANNs could be applied for inferring novel regulatory interactions in transcriptome data.
We applied FF and BF predictors to study the regulation of operons that participate in efflux transport systems related to multidrug resistance (MDR) in E. coli. Lately, clinically relevant strains of human pathogens exhibiting an MDR phenotype have been identified, including S. aureus, P. aeruginosa, and E. coli . To expel toxic molecules out of its outer membrane, these strains switch on their MDRs pumps. This mechanism may also be responsible for turning these strains resistant to antibiotics. Despite their potential role in turning bacteria resistant to antibiotics, which is a major health problem, little is known about these transport systems. Only 14 (~12%) of E. coli TFs have been described to modulate MDR pumps, which is a tiny proportion of the TRN . In addition, recent microarray studies [17–25] that analyzed E. coli K-12 grown under anaerobic, aerobic, alkaline, or acid conditions have shown that expression of several MDR pumps depends on a given stress conditions.
We built an expression dataset consisting of 58 Affymetrix chips that contained 95.36% of the bacteria genome and which hybridized with RNA samples from bacteria grown in different stress conditions, such as aerobic knock-out, anaerobic knock-out, and pH changes. Then, motif predictors were applied to investigate the transcriptional control of 16 operons that coded for key transport proteins in MDR pumps. This approach was useful to predict novel TFs that regulate the expression of drug transporters operons and which are involved in responding to those stress conditions.
Results and Discussion
The design of the ANNs as motif predictors comprised building the learning dataset, and training the networks and model assessment, as detailed below.
Assembling the Learning Dataset
The learning dataset is the set of examples that is used by the ANNs to adjust the weight for the classification task. Our learning dataset was built with two classes: true motifs examples that were extracted exclusively from the E. coli TRN available at RegulonDB version 5.0 , and non-motif samples.
Feature vectors were built in three steps. First, we applied Mfinder  to find the motifs (Figure 1A). Mfinder identifies interconnecting patterns between nodes (i.e., TFs and operons) that are more frequent than in random networks. Mfinder identified 199 instances of the FF motif and 926 instances of the BF motif in the TRN. These sets of motifs were labelled as the true motif class (TM).
Next, the examples for the non-motif class were generated (NM). For each true motif extracted, a non-motif example was generated by replacing one or more targets of the true motif, randomly chosen among the rest of the operons (Figure 1B). Thus, to create a non-BF motif we could change either the first, or the second, or both targets of the original sample. For instance, consider the BF motif with Crp and Fnr as common activators of acnA and hlyE. We can modify the targets of these regulators to say, mtr and zwf because it is known that they do not control these operons. Having this in mind, we assembled a balanced learning dataset with 50% of samples from each class.
The final step involved computing the correlations, and assembling the feature vectors. The expression data was used to compute different types of pairwise correlations, as well as partial correlations, for all operons and TFs in the microarray dataset. These correlation values were later used to assemble the feature vectors. Thus, the correlation values between profiles of motif and non-motifs elements mapped to feature vectors. The example in Figure 1C demonstrates how the 1, 2, 3, and 4 BF motifs, as well as the BF motifs 1, 2, 5, and 6 are represented as feature vectors (see legend of Figure 1C for a detailed description). We assembled six different feature vectors: Pearson correlation only (p), Spearman correlation (s), Kendall correlation (k), partial correlation (pc), Spearman/Kendall/Pearson (skp), and a sixth type containing all previous measures (all). Next, principal component analysis (PCA) was applied in order to reduce the dimensionality of the feature vectors and to optimize the classification task.
Using these learning datasets, we trained two types of neural networks, one for predicting FF motifs and one for predicting BF motifs. We trained the predictors using the six feature vectors designs (p, s, k, pc, skp, and all) mentioned earlier in the Assembling Learning Dataset section. All ANNs had the same architecture that included the following: (i) an input layer with a size that depended on the feature vector size after PCA, (ii) a 3-neurons hidden layer, and (iii) an output layer with a single classification neuron (Figure 1D). The over-fitting was avoided using one partition of the data as the validation dataset.
Validating and Assessing ANNs Performance
The FF and BF predictors were trained and evaluated using a 10 × 10-fold cross-validation, which is a widely adopted approach to evaluate classification algorithms.
This approach consists in splitting the learning dataset in 10 folds of equal size (last fold may have a slight different size). Then, eight folds are used for training (adjusting the weights of the ANN), one fold is used for validation and one fold is used for testing. The validation fold is useful for preventing over-fitting. The testing fold is used to compute the classification statistics, i.e. error, precision, recall (sensitivity) and f-measure (see Methods for definitions). In each cross-validation, all different possibilities of folds for training, testing and validation are chosen. The classification statistics for each cross-validation round are taken as the mean of the 10 possible iterations. These results are shown in Figure 2 (see Additional file 1 for the complete data on classification measures for each cross-validation, including confidence intervals). As shown in Figure 2, for the FF ANNs that use only one type of correlation (FFp, FFs, FFk, FFpc), the error rates were similar and around 30%, whereas recalls were nearly 70%. Using hybrid feature vectors, the classification results improved substantially. FFskp and FFall correctly identified 83% of the samples classified as FF motifs (high precision or positive prediction value) and achieved a high sensitivity (recall) of 83% and 86%, respectively. This indicates that hybrid network models had a low number of both false positives and false negatives. The same was observed for BF ANNs. Error and recall rates for single correlation feature vectors, BFp, BFs, BFk and BFpc, were approximately 30% and 80%, respectively. For the mixed correlation networks, BFskp and BFall, the precision also increased to 88% and 96%, respectively, and the recall to 90% and 97%, respectively. Figure 3 depicts the histograms for the output neuron (class neuron) for each network configuration. The neural network assigned to each sample in the training dataset a real value between -1 (motif) and +1 (non-motif). In this study, configurations using hybrid vectors were able to better separate the two classes because the values were clustered around -1 and 1.
Regulation of E. coli Multidrug Resistance Pumps
The emergence of pathogenic bacterial strains that are resistant to available antimicrobial agents is a major issue in public health. Bacteria strains acquire antibiotic resistance through mutation and selection (vertical evolution), or through exchanging of genetic material (horizontal evolution).
To be protected from the action of drugs, bacteria recruit several efflux pumps that may or may not be specific for a single substrate. These pumps are frequently associated with multidrug resistance (also known as MDR pumps). In prokaryotes, five types of MDR pumps [28, 29] have been identified. These are (i) the major facilitator superfamily (MSF), (ii) the ATP-binding cassette superfamily (ABC superfamily), (iii) the small multidrug resistance family (SMR family), (iv) the resistance-nodulation-cell division superfamily (NDR superfamily), and (v) the multidrug and toxic compound extrusion family (MATE) (Table 1). A great interest in how these pumps are regulated has generated substantial research on the regulatory pathways that govern the expression of MDR genes, resulting in the elucidation of transcriptional regulation at both local and global levels. Most efflux pumps known to be under regulatory controls belong to either the MFS or to the RND superfamilies  (Table 1). Both families employ the proton- (or sodium-) motive force to energize drug efflux and to trigger regulatory controls. These regulatory controls act by preventing an excessive production of nonspecific cation transport, and by preventing lost of membrane H+ potential, and even cell death. Unlike the MFS and the RND superfamilies, some drug pumps, including the SMR superfamily, do not have their synthesis controlled at the transcriptional level (Table 1).
Interestingly, there is increasing evidence that besides conferring resistance to drugs, MDR pumps also play other physiological roles that are required for bacteria survival in their ecological niche [29, 31]. In addition, recent studies using microarray have shown that E. coli K-12 that are grown under either anaerobic, aerobic, alkanine, or acid conditions express several MDR pumps [17, 20, 23]. Under oxygen limitations and acid conditions, the expression of the mdtEF multidrug resistance efflux increases. When pH increases to a basic condition, the expression of a different pump, namely the acridine efflux pump (acrAB), increases. These findings indicate their physiological roles in bacteria survival in natural ecosystems.
So far, only 14 (~12%) of E. coli TFs have been described as modulators of MDR pumps, which correspond to a very small proportion of the TRNs. Therefore, a better knowledge about drug transporters and their regulatory networks, along with structural analysis of prominent regulatory proteins, is required for developing drug pump inhibitors. Such pump inhibitors would be helpful to restore bacteria sensitivity to a specific drug or to reduce bacteria ability to colonize and survive in their host.
Investigating MDR Efflux Pumps Regulation
After validation in the learning dataset, we applied the best ANN motif predictors to investigate E. coli MDR efflux pumps, with the purpose to find new regulators for proteins involved in these transportation systems in bacteria. To generate the set of testing motif candidates, we first enumerated a target listing of 16 operons (23 genes), corresponding to all operons coding for transporters and accessory proteins as illustrated in Table 1: acrAB, acrD, acrEF, baeRSmdtABCD, cmr, emrABR, emrErenD, emrKY, macA, macB, mdlAB, mdtEF, mdtK, msbA, tolC, and yojI. The listing contains a comprehensive review of both well-known and poorly characterized MDR proteins identified so far in E. coli. Data was taken from the literature and from organism databases. The next step in our process involved creating motif candidates. For each target operon, large sets of FF and BF candidates were generated making use of the catalog of experimentally characterized inducers and repressors in E. coli. This set of 136 TFs (133 operons) selected from the last version of the TRN has been annotated according to their function and into several categories, from a to g, to facilitates the subsequent analysis of the interactions found by the ANNs (Figure 4, Additional file 5). These categories describe which cellular process(es) a particular regulator is/are controlling and include local and global regulators of MDR efflux pumps (14 instances). TFs directly related to pumping in general (regardless of the type of substract), along with other TFs promoting associated efflux pumps events (categories a to e), account for 25% of the total; the remaining 75% are regulatory proteins required to uptake transport and metabolism. Additionally, TFs were annotated with respect to a number of stress conditions, in which they regulate the expression of their target(s) (Figure 5, Additional file 5). This search revealed that most regulators showed expression dependence under anaerobiosis (32.46%), acid (21.93%) and basic (10.50%) growth conditions, while regulators able to function under aerobiosis (6.14%) and neutral pH (5.26%) are less frequent.
FF candidates were generated by making all 2-combinations from TFs: , i.e., ≅ 9·103 FF testing candidates for each target operon. BF candidates were produced by all 2-combinations from TFs along with operons, except a fixed one: , i.e., ≅ 130·103 BF candidates for every target operon.
In the last step, we applied the ANN predictors FFall and BFall to classify the motif hypotheses. Those are the predictors that outperformed other configurations in the curated dataset and which were trained with hybrid feature vectors. Figure 6 gives an overview of the FFall results and shows the proportion of candidates classified as true FF motifs using a cutoff of -0.9 for the motif class. It was observed that, for half of the operons, the motif class was assigned to < 4% of the inspected candidates. Figure 7 shows the general functioning of BFall. The true positive rates were ≤ 3.07%, for all evaluated MDR operons. This indicates that ANNs are discriminating only a small proportion of the candidates as a motif, which makes it feasible to be analyzed.
To assess the quality of the predictions, we annotated the set of inferred regulatory interactions using the types of evidence defined for the regulators (see Figure 4 for description of evidences). Regulatory interactions were annotated according to the five types of transcription regulators, namely a, b, c, d, and e (Additional file 5). The type a transcription regulator was a local regulator of an MDR efflux pump; type b was a global regulator of an MDR efflux pump, type c was a member of an MDR family of efflux pump regulators, type d was either a local or a global regulator of non-MDR efflux pumps, and type e was a regulator of proteins related to efflux pumps or secretion. Figure 8 presents the profile of the transcriptional interactions recovered by the FFall, for each MDR operon. To construct the chart shown in Figure 8, we considered that any FF motif has two regulatory connections (TF1-target, TF2-target) whereas BF has four, (TF1-target1, TF1-target2, TF2-target1, TF2-target2). Each slice in the charts corresponds to the proportion of a given type of interactions x, multiplied by a weighting factor 1/p(x), where p(x) is the probability of taking a TF of type x in the set of all TFs. This weighted proportion is necessary because the number of examples varies between regulators types. For instance, there are many more candidates involving regulators of type f than regulators of type a. As depicted in the charts, the inner donut shows the percentage for each interaction type, whereas in the outer donut there is the summed a to e and f to g percentages. The proportion of inferred a-e interactions is larger in all cases, especially in operons emrErendD (87%), and in baeRSmdtABCD (80%). Interactions of types a to e are the most promising ones, because they may represent a DNA-protein interaction between regulator and target, since those a to e TFs already have been related to regulation of efflux pumps in a direct way. The same was observed about the relationships inferred with BFall, where a to e slices were larger (Figure 9).
In general, this denotes that designed ANNs are preferably selecting MDR regulators, which is biologically consistent. Therefore, we believe that ANNs may infer very plausible hypotheses to be further inspected. Also, there are a considerable number of f to g interactions recovered by ANNs. Ideally, predictors should identify only a to e interactions, because f to g regulators are not responsible for immediately controlling transcription of MDR operons, apparently. These suggest that such interactions between f to g regulators with the MDR operons were identified because they belong to the same regulatory network of a given stress condition, and their expression profiles correlate likewise real motifs correlate. This also shows that these correlations are far from being the precise representations of molecular binding and they are only an approximate representation for quantifying testable interactions.
In-depth analysis of selected predictions
We selected some promising hypothesis for a more careful analysis. We first analyzed the inferred BF motif MarA and IscR regulating acrEB and mdlAB (classification = 0.99112), among the top ten predictions of acrEB. Currently it is known that upon reaching a high intracellular level, MarA activates its own transcription and several other genes from the Mar regulon, including acrAB (acrEB homolog), by binding to the marbox sequence in the promoter region that is upstream of the RNA polymerase site . The marbox is a 20 basepair (bp) asymmetric consensus sequence found in promoters of the E. coli marA/soxS/rob regulon(s), where each of the three activators, MarA, SoxS and Rob binds as a monomer . The marbox degenerate consensus was recently updated to: A≠GRGCACRWWNNRYYAAA≠GN (N = any base; R = A/G; W = A/T; Y = C/T) . Although AcrEF is homologous to AcrAB, both efflux pumps can function in combination with TolC and acrEF upregulation, which might compensate for an acrAB deletion . The acrEF operon has not been identified yet as a member of the mar regulon and its expression is known to be activated only by Fnr and repressed by AcrS (EnvR) . However, investigation of the promoter of acrEF has revealed a sequence of A A GGCAC AT A AA CA CAAA A A that is in the backward orientation and which extends from position 142 to position 161, upstream of the first acrEF codon (class I promoter ). This sequence matches the mar binding site consensus in precisly 14 positions, as defined in the acrAB promoter region (underlined), and in 15 out of the 17 most important positions (in bold) according to the most recent marbox consensus [32, 33]. Thus, it seems that the acrEF promoter region could be a false negative for the marbox consensus sequence containing only two positions, that are inverted compared with marbox acrAB (in italic), in disagreement with latest consensus . As a result, we suggest that the acrEF promoter region contains a putative mar/sox/rob activator binding site at -142 to -161. As such, the feature predicted here, that MarA regulates transcription of acrEF operon, is feasible for experimental testing. This BF motif also predicts that the regulation of mdlAB is MarA-dependent. The mdlA expression was shown to be enhanced by the MarA homolog Rob, which binds to the mdlA marbox promoter DNA . As mentioned above, the highly overlapping marA/soxS/rob regulon(s) contain promoters whose expression can be activated by the homologous transcriptional activators MarA, SoxS and Rob. Beacuse mdlA is known to be activated by Rob, its expression might also be marA-dependent, but that remains to be shown. At last, IscR (regulator 2) was identified by this hypothesis in combination with mdlA and acrEF. Presently, it is not known if the IscR regulation of the transcription of these targets are done directly or indirectly. In any way, IscR is a member of the Mar/Sox/Rob family of transcriptional regulators and therefore its effect on acrEF and mdlA can be tested in vitro.
Other worth-mentioning case involves the cmr (or mdfA) operon. The top 20 FF predictions (sorted by classification results) identified EmrR as a regulator of cmr, which belongs to the MarR family of transcriptional repressors. EmrR negatively regulates the emrRAB operon by directly binding to its regulatory region . This binding is antagonized by some EmrAB pump inducers (CCC, NA and TCS, Table 1) that interact with EmrR and prevent it from binding to this region [37, 38]. Besides, there is evidence showing that EmrR might play a regulatory role at the transcriptional level of the expression of several E. coli genes . At this point, we suggest that under a given condition of growth and without antibiotics that can antagonize the EmrR function, this regulator acts as a repressor of emrRAB and of potential target genes as well. Interestingly, emrA shows basic and anaerobic or acid and aerobic dependent expression , and thus provide possible evidence that under these conditions some metabolites can act as emrAB inducers that antagonize its repression by EmrR. Regarding the several hypotheses associating emrR and mdfA expressions, we argue that these associations might be representing a repression of EmrR to its targets under acid conditions because the MdfA function is at extreme alkaline pH growth. Both EmrAB and MdfA pumps act as drug/proton antiporters; also MdfA couples proton influx movement to antiport of Na or K ions under extreme alkaline conditions . For example, in an acid/anaerobic growth condition, emrAB can be turned off and this repression is mediated directly by EmrR. In a similar condition (acid), this repressor acts directly on the control of other functionally related pumps, such as mdfA. We found several other FF high-rated hypotheses and one BF motif that suggest this proposition (see Additional files 2 and 3). Finally, the second regulators that were identified in these predictions could regulate mdfA indirectly (positively or negatively), because they have been previously described to have an expression dependence (Additional file 5) on acid (note that the regulator FadR is off when referring to a repressor under a basic condition) neutral (IdnR), stress oxidative (IscR), anaerobiosis (PdhR, Fnr, TorR), and osmotic stress (KdpE).
ANNs presented here successfully recovered structural motifs available in E. coli TRN, reaching high precision and recall rates in this curated dataset. Therefore, the novel TFs predicted to be regulating the expression of MDR pumps proteins are likely candidates to be detected through DNA-protein binding experiments. Enriching the description of FF and BF motifs by using other dependence measures between variables, as mutual information [40, 41], could be a future improvement to increase reliability of the classifier. The Matlab source code and all accompanying data are available for downloading at the project website .
Assembling of the Expression Dataset
The microarray data comprises 58 E. coli Antisense Genome Array chips from Affymetrix with RNA from bacteria grown under different experimental conditions such as aerobic knock-out, anaerobic knock-out, changes of global gene expression during an oxygen shift, and pH changes, was retrieved under the accession code GPL199 from the GEO database . Quantification of probesets and normalization using MAS5 algorithm were performed with the R Bioconductor software package . Probesets annotation was done using the NetAffx tool . In the case of several probesets for the same gene, the mean expression was assigned as the transcript quantity. After computing gene expressions, the last step was to group the genes in operons and compute the expression level of the operon. Considering that an operon is a set of genes that transcribe together, it was reasonable to consider the mRNA levels of the operon as the mean expression of its constitutive genes. The final transcriptome dataset encompasses 784 operons and 2324 individual genes, summing up 4249 genes, or 95.36% of the genome (see Additional file 4).
ANNs Implementation Details
The entire design of the ANN classifiers was done in Matlab® 7.0, using the Neural Network Toolbox v. 4.0.2. Building the learning datasets, including computing of Pearson, Spearman, and Kendall correlations to assemble feature vectors, as well as principal component analysis (PCA), was done using functions of the Statistics Toolbox v. 5.0. The multi-layer perceptron ANNs were designed to have an input layer, a hidden layer with 3 neurons, and an output layer with the classification neuron. The adopted activation function of the perceptrons was the tangent sigmoid. Optimization of the weights was done using the Levenberg-Marquardt algorithm implemented in Matlab (trainlm function). The over-fitting was prevented by using the early stopping criterion. Thus, in each cross-validation, one partition of the data was set apart as validation dataset, and the training stopped when the error in the validation dataset started to increase (defined as the ability of the extrapolation of the network to get worse. Pre-processing of data was performed to scale inputs in the range [-1, 1]. The training parameters learning rate and momentum, were empirically adjusted to 0.5 and 0.1, through minimizing sum squared error (SSE). The size of the the hidden layer was set to 3 based on empirical simulations. The performance of the classification was assessed with the following measures: Error = FP/(TP + FN), Precision = TP/(TP + FP), Recall (sensitivity) = TP/(TP + FN) and F - measure = 2·(Precision·Recall)/(Precision·Recall), where TP is a motif correctly classified, FP is a non-motif classified in the opposite class, and FN is a motif classified in the non-motif group.
Artificial Neural Network(s)
Transcriptional Regulatory Network(s)
Multidrug and Toxic compound Extrusion
Major Facilitator Superfamily
Small Multidrug Resistance
Principal Component Analysis.
Martinez-Antonio A, Janga SC, Salgado H, Collado-Vides J: Internal-sensing machinery directs the activity of the regulatory network in Escherichia coli. Trends Microbiol. 2006, 14 (1): 22-27. 10.1016/j.tim.2005.11.002.
Seshasayee AS, Bertone P, Fraser GM, Luscombe NM: Transcriptional regulatory networks in bacteria: from input signals to output responses. Curr Opin Microbiol. 2006, 9 (5): 511-519. 10.1016/j.mib.2006.08.007.
Dojer N, Gambin A, Mizera A, Wilczynski B, Tiuryn J: Applying dynamic Bayesian networks to perturbed gene expression data. BMC Bioinformatics. 2006, 7: 249-10.1186/1471-2105-7-249.
Nachman I, Regev A, Friedman N: Inferring quantitative models of regulatory networks from expression data. Bioinformatics. 2004, 20 Suppl 1: I248-I256. 10.1093/bioinformatics/bth941.
Yu J, Smith VA, Wang PP, Hartemink AJ, Jarvis ED: Advances to Bayesian network inference for generating causal networks from observational biological data. Bioinformatics. 2004, 20 (18): 3594-3603. 10.1093/bioinformatics/bth448.
Zou M, Conzen SD: A new dynamic Bayesian network (DBN) approach for identifying gene regulatory networks from time course microarray data. Bioinformatics. 2005, 21 (1): 71-79. 10.1093/bioinformatics/bth463.
de la Fuente A, Bing N, Hoeschele I, Mendes P: Discovery of meaningful associations in genomic data using partial correlation coefficients. Bioinformatics. 2004, 20 (18): 3565-3574. 10.1093/bioinformatics/bth445.
Schafer J, Strimmer K: An empirical Bayes approach to inferring large-scale gene association networks. Bioinformatics. 2005, 21 (6): 754-764. 10.1093/bioinformatics/bti062.
Wille A, Zimmermann P, Vranova E, Furholz A, Laule O, Bleuler S, Hennig L, Prelic A, von Rohr P, Thiele L, Zitzler E, Gruissem W, Buhlmann P: Sparse graphical Gaussian modeling of the isoprenoid gene network in Arabidopsis thaliana. Genome Biol. 2004, 5 (11): R92-10.1186/gb-2004-5-11-r92.
Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, Califano A: Reverse engineering of regulatory networks in human B cells. Nature genetics. 2005, 37 (4): 382-390. 10.1038/ng1532.
Butte AJ, Kohane IS: Mutual information relevance networks: functional genomic clustering using pairwise entropy measurements. Pac Symp Biocomput. 2000, 418-429.
Margolin AA, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla Favera R, Califano A: ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics. 2006, 7 Suppl 1: S7-10.1186/1471-2105-7-S1-S7.
Shen-Orr SS, Milo R, Mangan S, Alon U: Network motifs in the transcriptional regulation network of Escherichia coli. Nature genetics. 2002, 31 (1): 64-68. 10.1038/ng881.
Resendis-Antonio O, Freyre-Gonzalez JA, Menchaca-Mendez R, Gutierrez-Rios RM, Martinez-Antonio A, Avila-Sanchez C, Collado-Vides J: Modular analysis of the transcriptional regulatory network of E. coli. Trends Genet. 2005, 21 (1): 16-20. 10.1016/j.tig.2004.11.010.
Tenover FC: Mechanisms of antimicrobial resistance in bacteria. Am J Infect Control. 2006, 34 (5 Suppl 1): S3-10; discussion S64-73. 10.1016/j.ajic.2006.05.219.
Kumar A, Schweizer HP: Bacterial resistance to antibiotics: active efflux and reduced uptake. Adv Drug Deliv Rev. 2005, 57 (10): 1486-1513. 10.1016/j.addr.2005.04.004.
Hayes ET, Wilks JC, Sanfilippo P, Yohannes E, Tate DP, Jones BD, Radmacher MD, BonDurant SS, Slonczewski JL: Oxygen limitation modulates pH regulation of catabolism and hydrogenases, multidrug transporters, and envelope composition in Escherichia coli K-12. BMC Microbiol. 2006, 6: 89-10.1186/1471-2180-6-89.
Kang Y, Weber KD, Qiu Y, Kiley PJ, Blattner FR: Genome-wide expression analysis indicates that FNR of Escherichia coli K-12 regulates a large number of genes of unknown function. J Bacteriol. 2005, 187 (3): 1135-1160. 10.1128/JB.187.3.1135-1160.2005.
Kirkpatrick C, Maurer LM, Oyelakin NE, Yoncheva YN, Maurer R, Slonczewski JL: Acetate and formate stress: opposite responses in the proteome of Escherichia coli. J Bacteriol. 2001, 183 (21): 6466-6477. 10.1128/JB.183.21.6466-6477.2001.
Lewinson O, Padan E, Bibi E: Alkalitolerance: a biological function for a multidrug transporter in pH homeostasis. Proc Natl Acad Sci U S A. 2004, 101 (39): 14073-14078. 10.1073/pnas.0405375101.
Maurer LM, Yohannes E, Bondurant SS, Radmacher M, Slonczewski JL: pH regulates genes for flagellar motility, catabolism, and oxidative stress in Escherichia coli K-12. J Bacteriol. 2005, 187 (1): 304-319. 10.1128/JB.187.1.304-319.2005.
Overton TW, Griffiths L, Patel MD, Hobman JL, Penn CW, Cole JA, Constantinidou C: Microarray analysis of gene regulation by oxygen, nitrate, nitrite, FNR, NarL and NarP during anaerobic growth of Escherichia coli: new insights into microbial physiology. Biochem Soc Trans. 2006, 34 (Pt 1): 104-107.
Salmon K, Hung SP, Mekjian K, Baldi P, Hatfield GW, Gunsalus RP: Global gene expression profiling in Escherichia coli K12. The effects of oxygen availability and FNR. J Biol Chem. 2003, 278 (32): 29837-29855. 10.1074/jbc.M213060200.
Salmon KA, Hung SP, Steffen NR, Krupp R, Baldi P, Hatfield GW, Gunsalus RP: Global gene expression profiling in Escherichia coli K12: effects of oxygen availability and ArcA. J Biol Chem. 2005, 280 (15): 15084-15096. 10.1074/jbc.M414030200.
Yohannes E, Thurber AE, Wilks JC, Tate DP, Slonczewski JL: Polyamine stress at high pH in Escherichia coli K-12. BMC Microbiol. 2005, 5: 59-10.1186/1471-2180-5-59.
Salgado H, Gama-Castro S, Peralta-Gil M, Diaz-Peredo E, Sanchez-Solano F, Santos-Zavaleta A, Martinez-Flores I, Jimenez-Jacinto V, Bonavides-Martinez C, Segura-Salazar J, Martinez-Antonio A, Collado-Vides J: RegulonDB (version 5.0): Escherichia coli K-12 transcriptional regulatory network, operon organization, and growth conditions. Nucleic Acids Res. 2006, 34 (Database issue): D394-7. 10.1093/nar/gkj156.
Kashtan N, Itzkovitz S, Milo R, Alon U: Efficient sampling algorithm for estimating subgraph concentrations and detecting network motifs. Bioinformatics. 2004, 20 (11): 1746-1758. 10.1093/bioinformatics/bth163.
Paulsen IT: Multidrug efflux pumps and resistance: regulation and evolution. Curr Opin Microbiol. 2003, 6 (5): 446-451. 10.1016/j.mib.2003.08.005.
Piddock LJ: Multidrug-resistance efflux pumps - not just for resistance. Nat Rev Microbiol. 2006, 4 (8): 629-636. 10.1038/nrmicro1464.
Grkovic S, Brown MH, Skurray RA: Regulation of bacterial drug export systems. Microbiol Mol Biol Rev. 2002, 66 (4): 671-701, table of contents. 10.1128/MMBR.66.4.671-701.2002.
Krulwich TA, Lewinson O, Padan E, Bibi E: Do physiological roles foster persistence of drug/multidrug-efflux transporters? A case study. Nat Rev Microbiol. 2005, 3 (7): 566-572. 10.1038/nrmicro1181.
Martin RG, Gillette WK, Rhee S, Rosner JL: Structural requirements for marbox function in transcriptional activation of mar/sox/rob regulon promoters in Escherichia coli: sequence, orientation and spatial relationship to the core promoter. Mol Microbiol. 1999, 34 (3): 431-441. 10.1046/j.1365-2958.1999.01599.x.
Martin RG, Rosner JL: Genomics of the marA/soxS/rob regulon of Escherichia coli: identification of directly activated promoters by application of molecular genetics and informatics to microarray data. Mol Microbiol. 2002, 44 (6): 1611-1624. 10.1046/j.1365-2958.2002.02985.x.
Jellen-Ritter AS, Kern WV: Enhanced expression of the multidrug efflux pumps AcrAB and AcrEF associated with insertion element transposition in Escherichia coli mutants Selected with a fluoroquinolone. Antimicrob Agents Chemother. 2001, 45 (5): 1467-1472. 10.1128/AAC.45.5.1467-1472.2001.
Nishino K, Yamaguchi A: Analysis of a complete library of putative drug transporter genes in Escherichia coli. J Bacteriol. 2001, 183 (20): 5803-5812. 10.1128/JB.183.20.5803-5812.2001.
Bennik MH, Pomposiello PJ, Thorne DF, Demple B: Defining a rob regulon in Escherichia coli by using transposon mutagenesis. J Bacteriol. 2000, 182 (13): 3794-3801. 10.1128/JB.182.13.3794-3801.2000.
Xiong A, Gottman A, Park C, Baetens M, Pandza S, Matin A: The EmrR protein represses the Escherichia coli emrRAB multidrug resistance operon by directly binding to its promoter region. Antimicrob Agents Chemother. 2000, 44 (10): 2905-2907. 10.1128/AAC.44.10.2905-2907.2000.
Brooun A, Tomashek JJ, Lewis K: Purification and ligand binding of EmrR, a regulator of a multidrug transporter. J Bacteriol. 1999, 181 (16): 5131-5133.
del Castillo I, Gomez JM, Moreno F: mprA, an Escherichia coli gene that reduces growth-phase-dependent synthesis of microcins B17 and C7 and blocks osmoinduction of proU when cloned on a high-copy-number plasmid. J Bacteriol. 1990, 172 (1): 437-445.
Daub CO, Steuer R, Selbig J, Kloska S: Estimating mutual information using B-spline functions--an improved similarity measure for analysing gene expression data. BMC Bioinformatics. 2004, 5: 118-10.1186/1471-2105-5-118.
Steuer R, Kurths J, Daub CO, Weise J, Selbig J: The mutual information: detecting and evaluating dependencies between variables. Bioinformatics. 2002, 18 Suppl 2: S231-40.
Project website. [http://www.labinfo.lncc.br/publicacao/]
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, 35 (Database issue): D760-5. 10.1093/nar/gkl887.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Cheng J, Sun S, Tracy A, Hubbell E, Morris J, Valmeekam V, Kimbrough A, Cline MS, Liu G, Shigeta R, Kulp D, Siani-Rose MA: NetAffx Gene Ontology Mining Tool: a visual approach for microarray data analysis. Bioinformatics. 2004, 20 (9): 1462-1463. 10.1093/bioinformatics/bth087.
Keseler IM, Collado-Vides J, Gama-Castro S, Ingraham J, Paley S, Paulsen IT, Peralta-Gil M, Karp PD: EcoCyc: a comprehensive database resource for Escherichia coli. Nucleic Acids Res. 2005, 33 (Database issue): D334-7. 10.1093/nar/gki108.
The Universal Protein Resource (UniProt). Nucleic Acids Res. 2007, 35 (Database issue): D193-7.
This work was funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), under grant number 506321/2004-5. We thank many lab colleagues, and Marcia Triunfol for editing the manuscript.
DFTV participated in the design of the study, performed neural networks training and evaluation, contributed to biological analysis, and drafted the manuscript; FFRV conceived the original idea, assembled the learning datasets, and helped to draft the manuscript; MFN designed and wrote the biological analysis; and ATRV managed the project and helped to improve the text. All authors read and approved the final manuscript.
Diogo FT Veiga, Fábio FR Vicente contributed equally to this work.
Electronic supplementary material
Additional file 1: Error, precision, recall, and f-measure rates for the ANN classifiers in the learning dataset. Classification statistics were measured using a 10 × 10-fold cross-validation procedure for the six different feature vector types available: Pearson (p), Spearman (s), Kendall (k), partial correlation (pc), Spearman/Kendall/Pearson (skp), and the last type containing all previous measures (all). Confidence intervals were estimated with α = 0.05 (probability = 95%). (XLS 30 KB)
Additional file 2: Complete list of predicted FF motifs involving interactions of types (a)-(e), which were classified by FFall using a threshold of -0.9. (XLS 672 KB)
Additional file 3: Complete list of predicted BF motifs involving interactions of types (a)-(e), which were classified by BFall using a threshold of -0.9. (XLS 2 MB)
Additional file 4: Microarray dataset used in this study (see Methods for detail on pre-processing routines used to generate this data). (TXT 2 MB)
Functional classification and stress condition-dependent regulatory patterns for transcriptional regulators investigated in the study*
Additional file 5: . Each E. coli regulator was annotated according to its functional characteristics such as whether a gene or a controlling target operon, whether a global or local level of transcription, and whether an activator, a repressor, or both. * As found in expression studies [17–25] and E. coli databases [26, 46, 47] (accessions date: Feb 2007). (DOC 170 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Veiga, D.F., Vicente, F.F., Nicolás, M.F. et al. Predicting transcriptional regulatory interactions with artificial neural networks applied to E. coli multidrug resistance efflux pumps. BMC Microbiol 8, 101 (2008). https://doi.org/10.1186/1471-2180-8-101
- Feature Vector
- Efflux Pump
- Network Motif
- True Motif
- Small Multidrug Resistance