Assessment of a novel multi-array normalization method based on spike-in control probes suitable for microRNA datasets with global decreases in expression.

BACKGROUND: High-quality expression data are required to investigate the biological effects of microRNAs (miRNAs). The goal of this study was, first, to assess the quality of miRNA expression data based on microarray technologies and, second, to consolidate it by applying a novel normalization method. Indeed, because of significant differences in platform designs, miRNA raw data cannot be normalized blindly with standard methods developed for gene expression. This fundamental observation motivated the development of a novel multi-array normalization method based on controllable assumptions, which uses the spike-in control probes to adjust the measured intensities across arrays. RESULTS: Raw expression data were obtained with the Exiqon dual-channel miRCURY LNA™ platform in the "common reference design" and processed as "pseudo-single-channel". They were used to apply several quality metrics based on the coefficient of variation and to test the novel spike-in controls based normalization method. Most of the considerations presented here could be applied to raw data obtained with other platforms. To assess the normalization method, it was compared with 13 other available approaches from both data quality and biological outcome perspectives. The results showed that the novel multi-array normalization method reduced the data variability in the most consistent way. Further, the reliability of the obtained differential expression values was confirmed based on a quantitative reverse transcription-polymerase chain reaction experiment performed for a subset of miRNAs. The results reported here support the applicability of the novel normalization method, in particular to datasets that display global decreases in miRNA expression similarly to the cigarette smoke-exposed mouse lung dataset considered in this study. CONCLUSIONS: Quality metrics to assess between-array variability were used to confirm that the novel spike-in controls based normalization method provided high-quality miRNA expression data suitable for reliable downstream analysis. The multi-array miRNA raw data normalization method was implemented in an R software package called ExiMiR and deposited in the Bioconductor repository.

Construction of a computable cell proliferation network focused on non-diseased lung cells.

BACKGROUND: Critical to advancing the systems-level evaluation of complex biological processes is the development of comprehensive networks and computational methods to apply to the analysis of systems biology data (transcriptomics, proteomics/phosphoproteomics, metabolomics, etc.). Ideally, these networks will be specifically designed to capture the normal, non-diseased biology of the tissue or cell types under investigation, and can be used with experimentally generated systems biology data to assess the biological impact of perturbations like xenobiotics and other cellular stresses. Lung cell proliferation is a key biological process to capture in such a network model, given the pivotal role that proliferation plays in lung diseases including cancer, chronic obstructive pulmonary disease (COPD), and fibrosis. Unfortunately, no such network has been available prior to this work. RESULTS: To further a systems-level assessment of the biological impact of perturbations on non-diseased mammalian lung cells, we constructed a lung-focused network for cell proliferation. The network encompasses diverse biological areas that lead to the regulation of normal lung cell proliferation (Cell Cycle, Growth Factors, Cell Interaction, Intra- and Extracellular Signaling, and Epigenetics), and contains a total of 848 nodes (biological entities) and 1597 edges (relationships between biological entities). The network was verified using four published gene expression profiling data sets associated with measured cell proliferation endpoints in lung and lung-related cell types. Predicted changes in the activity of core machinery involved in cell cycle regulation (RB1, CDKN1A, and MYC/MYCN) are statistically supported across multiple data sets, underscoring the general applicability of this approach for a network-wide biological impact assessment using systems biology data. CONCLUSIONS: To the best of our knowledge, this lung-focused Cell Proliferation Network provides the most comprehensive connectivity map in existence of the molecular mechanisms regulating cell proliferation in the lung. The network is based on fully referenced causal relationships obtained from extensive evaluation of the literature. The computable structure of the network enables its application to the qualitative and quantitative evaluation of cell proliferation using systems biology data sets. The network is available for public use.

Chronic exposure to cigarette smoke condensate in vitro induces epithelial to mesenchymal transition-like changes in human bronchial epithelial cells, BEAS-2B.

Cigarette smoke causes lung tumorigenesis; however, the mechanisms underlying transformation are unknown. We investigated if tobacco compounds induce DNA promoter hypermethylation in BEAS-2B cells treated with low doses of cigarette smoke condensate (CSC) for one month. Transcriptional profiles and anchorage-independent growth were explored using Affymetrix microarray and soft agar assay, respectively. To investigate if tobacco compounds induce hypermethylation, CSC/dimethyl sulfoxide (DMSO)-treated cells were further treated with 5-Aza-2'-deoxycytidine (5AzaC) and trychostatin A (TSA). This treatment was followed by transcriptional profiling. CSC-exposed cells acquired a fibroblast-like shape with enhanced anchorage-independent growth. Silencing of epithelial cadherin, the hallmark of epithelial to mesenchymal transition (EMT), was observed upon exposure to CSC. Changes in the expression of genes involved in epidermal development, intercellular junction formation, and cytoskeleton formation were identified. Gene expression profiles from 5AzaC- and TSA-treated cells revealed 130 genes possibly methylated due to chronic CSC exposure. Our results suggest that E-cadherin may also be silenced by hypermethylation in an in vitro model of chronic exposure to low doses of CSC. This study demonstrates evidence for a tobacco compound induced EMT-like process in vitro and provides insight into possible mechanisms of gene silencing occurring during this treatment.

cis-Effect of DnaJ on DnaK in ternary complexes with chimeric DnaK/DnaJ-binding peptides.

Chimeric peptides, comprising a DnaK-binding sequence of L-amino acid residues (motif k) and an exclusive DnaJ-binding sequence of D-amino acid residues (motif j) connected through a 22-residue linker, were examined as minisubstrates for the DnaK chaperone system. The DnaJ-stimulated ATPase activity of DnaK was three times higher in the presence of the chimeric peptides pjk or pkj than in the simultaneous presence of the corresponding single-motif peptides ala-p5 (k motif) plus D-p5 (j motif). Apparently, pjk and pkj mimic unfolded proteins by forming ternary (ATP x DnaK) x peptide x DnaJ complexes which favor cis-interaction of DnaJ with DnaK. Consistent with this interpretation, the specific stimulatory effect of the chimeric peptides was abolished by either single-motif peptide in excess.

Mechanism of the targeting action of DnaJ in the DnaK molecular chaperone system.

In the DnaK (Hsp70) molecular chaperone system of Escherichia coli, the substrate polypeptide is fed into the chaperone cycle by association with the fast-binding, ATP-liganded form of the DnaK. The substrate binding properties of DnaK are controlled by its two cochaperones DnaJ (Hsp40) and GrpE. DnaJ stimulates the hydrolysis of DnaK-bound ATP, and GrpE accelerates ADP/ATP exchange. DnaJ has been described as targeting the substrate to DnaK, a concept that has remained rather obscure. Based on binding experiments with peptides and polypeptides we propose here a novel mechanism for the targeting action of DnaJ: ATP.DnaK and DnaJ with its substrate-binding domain bind to different segments of one and the same polypeptide chain forming (ATP.DnaK)m.substrate.DnaJn complexes; in these ternary complexes efficient cis-interaction of the J-domain of DnaJ with DnaK is favored by their propinquity and triggers the hydrolysis of DnaK-bound ATP, converting DnaK to its ADP-liganded high affinity state and thus locking it onto the substrate polypeptide.

D-Peptides as inhibitors of the DnaK/DnaJ/GrpE chaperone system.

DnaK, a Hsp70 homolog of Escherichia coli, together with its co-chaperones DnaJ and GrpE protects denatured proteins from aggregation and promotes their refolding by an ATP-consuming mechanism. DnaJ not only stimulates the gamma-phosphate cleavage of DnaK-bound ATP but also binds polypeptide substrates on its own. Unfolded polypeptides, such as denatured luciferase, thus form ternary complexes with DnaJ and DnaK. A previous study has shown that d-peptides compete with l-peptides for the same binding site in DnaJ but do not bind to DnaK (Feifel, B., Schönfeld, H.-J., and Christen, P. (1998) J. Biol. Chem. 273, 11999-12002). Here we report that d-peptides efficiently inhibit the refolding of denatured luciferase by the DnaK/DnaJ/GrpE chaperone system (EC50 = 1-2 microM). The inhibition of the chaperone action is due to the binding of d-peptide to DnaJ (Kd = 1-2 microM), which seems to preclude DnaJ from forming ternary (ATP.DnaK)m.substrate.DnaJn complexes. Apparently, simultaneous binding of DnaJ and DnaK to one and the same target polypeptide is essential for effective chaperone action.

Interdomain communication in the molecular chaperone DnaK.

DnaK, a heat-shock protein 70 (Hsp70) homologue in Escherichia coli, possesses a single tryptophan residue in its ATPase domain. Changes in the intrinsic fluorescence of DnaK offer a simple means not only to follow the binding of ATP and of ADP plus the co-chaperone GrpE to the ATPase domain, but also to investigate the kinetics of peptide binding to the substrate-binding domain of ATP.DnaK and GrpE-liganded ADP.DnaK. Addition of ATP or of ADP plus GrpE to nucleotide-free DnaK resulted in a similar decrease in intrinsic fluorescence, indicating similar open conformations of the ATPase domain under these two conditions. Binding of peptide increased the intrinsic fluorescence of both ATP.DnaK and ADP.DnaK.GrpE and rendered their spectra similar to the spectrum of ADP.DnaK with closed conformation of the ATPase domain. These results, together with the differential kinetics of peptide binding to ADP.DnaK on the one hand, and to ATP.DnaK or ADP.DnaK.GrpE on the other, suggest that ligands for either domain, i.e. ATP or ADP plus GrpE for the ATPase domain and peptides for the substrate-binding domain, shift the conformational equilibrium of both domains of DnaK towards the open and closed forms, respectively, in a concerted and parallel manner.