SUMMARY: Highly parallel genomic platforms like microarrays often present researchers with long lists of differentially expressed genes but contain little or no information on how these genes are regulated. rHVDM is a novel R package which uses gene expression time course data to predict the activity and targets of a transcription factor. In the first step, rHVDM uses a small number of known targets to derive the activity profile of a given transcription factor. Then, in a subsequent step, this activity profile is used to predict other putative targets of that transcription factor. A dynamic and mechanistic model of gene expression is at the heart of the technique. Measurement error is taken into account during the process, which allows an objective assessment of the robustness of fit and, therefore, the quality of the predictions. The package relies on efficient algorithms and vectorization to accomplish potentially time consuming tasks including optimization and differential equation integration. We demonstrate the efficiency and accuracy of rHVDM by examining the activity of the tumour-suppressing transcription factor, p53. AVAILABILITY: The version of the package presented here (1.8.1) is freely available from the Bioconductor Web site (http://bioconductor.org/packages/2.3/bioc/html/rHVDM.html).

BACKGROUND: Gene expression microarray data is notoriously subject to high signal variability. Moreover, unavoidable variation in the concentration of transcripts applied to microarrays may result in poor scaling of the summarized data which can hamper analytical interpretations. This is especially relevant in a systems biology context, where systematic biases in the signals of particular genes can have severe effects on subsequent analyses. Conventionally it would be necessary to replace the mismatched arrays, but individual time points cannot be rerun and inserted because of experimental variability. It would therefore be necessary to repeat the whole time series experiment, which is both impractical and expensive. RESULTS: We explain how scaling mismatches occur in data summarized by the popular MAS5 (GCOS; Affymetrix) algorithm, and propose a simple recursive algorithm to correct them. Its principle is to identify a set of constant genes and to use this set to rescale the microarray signals. We study the properties of the algorithm using artificially generated data and apply it to experimental data. We show that the set of constant genes it generates can be used to rescale data from other experiments, provided that the underlying system is similar to the original. We also demonstrate, using a simple example, that the method can successfully correct existing imbalances in the data. CONCLUSION: The set of constant genes obtained for a given experiment can be applied to other experiments, provided the systems studied are sufficiently similar. This type of rescaling is especially relevant in systems biology applications using microarray data.