A systems biology approach reveals the dose- and time-dependent effect of primary human airway epithelium tissue culture after exposure to cigarette smoke in vitro.

To establish a relevant in vitro model for systems toxicology-based mechanistic assessment of environmental stressors such as cigarette smoke (CS), we exposed human organotypic bronchial epithelial tissue cultures at the air liquid interface (ALI) to various CS doses. Previously, we compared in vitro gene expression changes with published human airway epithelia in vivo data to assess their similarities. Here, we present a follow-up evaluation of these in vitro transcriptomics data, using complementary computational approaches and an integrated mRNA-microRNA (miRNA) analysis. The main cellular pathways perturbed by CS exposure were related to stress responses (oxidative stress and xenobiotic metabolism), inflammation (inhibition of nuclear factor-κB and the interferon gamma-dependent pathway), and proliferation/differentiation. Within post-exposure periods up to 48 hours, a transient kinetic response was observed at lower CS doses, whereas higher doses resulted in more sustained responses. In conclusion, this systems toxicology approach has potential for product testing according to "21st Century Toxicology".

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.

A modular cell-type focused inflammatory process network model for non-diseased pulmonary tissue.

Exposure to environmental stressors such as cigarette smoke (CS) elicits a variety of biological responses in humans, including the induction of inflammatory responses. These responses are especially pronounced in the lung, where pulmonary cells sit at the interface between the body's internal and external environments. We combined a literature survey with a computational analysis of multiple transcriptomic data sets to construct a computable causal network model (the Inflammatory Process Network (IPN)) of the main pulmonary inflammatory processes. The IPN model predicted decreased epithelial cell barrier defenses and increased mucus hypersecretion in human bronchial epithelial cells, and an attenuated pro-inflammatory (M1) profile in alveolar macrophages following exposure to CS, consistent with prior results. The IPN provides a comprehensive framework of experimentally supported pathways related to CS-induced pulmonary inflammation. The IPN is freely available to the scientific community as a resource with broad applicability to study the pathogenesis of pulmonary disease.

Construction of a computable network model for DNA damage, autophagy, cell death, and senescence.

Towards the development of a systems biology-based risk assessment approach for environmental toxicants, including tobacco products in a systems toxicology setting such as the "21st Century Toxicology", we are building a series of computable biological network models specific to non-diseased pulmonary and cardiovascular cells/tissues which capture the molecular events that can be activated following exposure to environmental toxicants. Here we extend on previous work and report on the construction and evaluation of a mechanistic network model focused on DNA damage response and the four main cellular fates induced by stress: autophagy, apoptosis, necroptosis, and senescence. In total, the network consists of 34 sub-models containing 1052 unique nodes and 1538 unique edges which are supported by 1231 PubMed-referenced literature citations. Causal node-edge relationships are described using the Biological Expression Language (BEL), which allows for the semantic representation of life science relationships in a computable format. The Network is provided in .XGMML format and can be viewed using freely available network visualization software, such as Cytoscape.

Human bronchial epithelial cells exposed in vitro to cigarette smoke at the air-liquid interface resemble bronchial epithelium from human smokers.

Organotypic culture of human primary bronchial epithelial cells is a useful in vitro system to study normal biological processes and lung disease mechanisms, to develop new therapies, and to assess the biological perturbations induced by environmental pollutants. Herein, we investigate whether the perturbations induced by cigarette smoke (CS) and observed in the epithelium of smokers' airways are reproducible in this in vitro system (AIR-100 tissue), which has been shown to recapitulate most of the characteristics of the human bronchial epithelium. Human AIR-100 tissues were exposed to mainstream CS for 7, 14, 21, or 28 min at the air-liquid interface, and we investigated various biological endpoints [e.g., gene expression and microRNA profiles, matrix metalloproteinase 1 (MMP-1) release] at multiple postexposure time points (0.5, 2, 4, 24, 48 h). By performing a Gene Set Enrichment Analysis, we observed a significant enrichment of human smokers' bronchial epithelium gene signatures derived from different public transcriptomics datasets in CS-exposed AIR-100 tissue. Comparison of in vitro microRNA profiles with microRNA data from healthy smokers highlighted various highly translatable microRNAs associated with inflammation or with cell cycle processes that are known to be perturbed by CS in lung tissue. We also found a dose-dependent increase of MMP-1 release by AIR-100 tissue 48 h after CS exposure in agreement with the known effect of CS on this collagenase expression in smokers' tissues. In conclusion, a similar biological perturbation than the one observed in vivo in smokers' airway epithelium could be induced after a single CS exposure of a human organotypic bronchial epithelium-like tissue culture.

Nrf2: friend and foe in preventing cigarette smoking-dependent lung disease.

Chronic exposure to cigarette smoke (CS) generally confronts cellular defense systems with one of the strongest known environmental challenges. In particular, the continuous exposure of tissues of the respiratory tract to abundant concentrations of radicals; volatile compounds of the gas phase, mainly reactive oxygen and nitrogen species; and CS condensate deposits trigger a pleiotropic adaptive response, generally aimed at restoring tissue homeostasis. As documented by numerous studies published over the past decade, a hallmark of this defense system is the activation of the transcription factor NF-E2-related factor 2 (Nrf2), which, consequent to its established role as master regulator of the cellular antioxidant response, has been shown to orchestrate the first line of defense against cell- and tissue-damaging components present in CS. The key to CS-dependent Nrf2 activation is assumed to be based on the long-known phenomenon of a general strong sulfhydryl (-SH) reactivity inherent to CS. This chemical trait is virtually predestined to be sensitized by the major route leading to Nrf2 activation, characterized by its dependence on the interaction of electrophiles with specific cysteine residues inherited by Nrf2's negative cytosolic regulator Keap1 (Kelch-like ECH-associated protein 1). In addition, other pathways involving CS-activated protein kinases implicated in the upstream regulation of Nrf2, such as protein kinase C, represent an alternative/complementary mechanism of CS-induced Nrf2 activation. Because of the outstanding function of the Nrf2-Keap1 axis in defending cells and tissues against oxidant and chemical stress, either directly or indirectly via cross-talking with other defense pathways, changes in the Nrf2 or Keap1 genotype have long been associated with disease development. In terms of the two major smoking-related diseases of the lung, that is, emphysema and lung cancer, a fully functional Nrf2 genotype seems to be necessary, although not sufficient by itself, to protect the smoker from acquiring emphysema. Contrasting with this protective role, however, Nrf2 function may be potentially fatal in smoking-related lung tumorigenesis: as concluded from recent clinical investigations, lung tumor tissues harbor increased mutation or, alternatively, aberrant expression rates in either the KEAP1 or the NRF2 gene, generally resulting in constitutive Nrf2 activation, suggesting that "abuse" of Nrf2 function is an advantageous strategy of the (developing) tumor to protect itself against oxidative stress in general. On the basis of the fundamental significance of the Nrf2 pathway in smoking-dependent disease development, several attempts have been described for dietary and pharmacological intervention, the majority of which are intended to activate Nrf2 aiming at emphysema prevention. The intention of this review is to compile and discuss the various aspects of CS-Nrf2/Keap1 interaction in terms of mechanism, disease development, and chemoprevention.

A computable cellular stress network model for non-diseased pulmonary and cardiovascular tissue.

BACKGROUND: Humans and other organisms are equipped with a set of responses that can prevent damage from exposure to a multitude of endogenous and environmental stressors. If these stress responses are overwhelmed, this can result in pathogenesis of diseases, which is reflected by an increased development of, e.g., pulmonary and cardiac diseases in humans exposed to chronic levels of environmental stress, including inhaled cigarette smoke (CS). Systems biology data sets (e.g., transcriptomics, phosphoproteomics, metabolomics) could enable comprehensive investigation of the biological impact of these stressors. However, detailed mechanistic networks are needed to determine which specific pathways are activated in response to different stressors and to drive the qualitative and eventually quantitative assessment of these data. A current limiting step in this process is the availability of detailed mechanistic networks that can be used as an analytical substrate. RESULTS: We have built a detailed network model that captures the biology underlying the physiological cellular response to endogenous and exogenous stressors in non-diseased mammalian pulmonary and cardiovascular cells. The contents of the network model reflect several diverse areas of signaling, including oxidative stress, hypoxia, shear stress, endoplasmic reticulum stress, and xenobiotic stress, that are elicited in response to common pulmonary and cardiovascular stressors. We then tested the ability of the network model to identify the mechanisms that are activated in response to CS, a broad inducer of cellular stress. Using transcriptomic data from the lungs of mice exposed to CS, the network model identified a robust increase in the oxidative stress response, largely mediated by the anti-oxidant NRF2 pathways, consistent with previous reports on the impact of CS exposure in the mammalian lung. CONCLUSIONS: The results presented here describe the construction of a cellular stress network model and its application towards the analysis of environmental stress using transcriptomic data. The proof-of-principle analysis described here, coupled with the future development of additional network models covering distinct areas of biology, will help to further clarify the integrated biological responses elicited by complex environmental stressors such as CS, in pulmonary and cardiovascular cells.

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.

Endoplasmic reticulum stress induced by aqueous extracts of cigarette smoke in 3T3 cells activates the unfolded-protein-response-dependent PERK pathway of cell survival.

Cigarette smoke (CS) generally places severe stress on cells, as reflected by gene expression profiling and pathway analysis, which, among other effects, also suggested activation of the unfolded protein response pathway triggered by the stressed endoplasmic reticulum (ER stress). Here, we present data indicating that noncytotoxic concentrations of aqueous extracts of CS induce a distinct ER stress response in immortalized nontransformed Swiss 3T3 cells, primarily by activating the PERK pathway of global protein synthesis inhibition. Activation of PERK and PERK-dependent signaling by aqueous extracts of CS was demonstrated by (i) the inhibition of protein synthesis, (ii) the phosphorylation of PERK and its substrate eIF2alpha, (iii) the activation of ATF4, and (iv) the expression of ATF4-dependent target genes chop, gadd34, BiP, and atf3. Within the dose range tested, all effects appeared to be transient in nature, while the periods of recovery from ER stress were clearly concentration dependent. In contrast to these data and to the effects seen with thapsigargin (used as positive control), only minor effects were observed for the activation of xbp-1, a common target of the other two canonical sensors of ER stress, i.e., ATF6 and IRE1. In mechanistic terms, neither the disruption of energy levels nor a contribution of arylating quinones played a major role under the experimental conditions tested. Notably however, the effects of aqueous extracts of CS on the ER could be mimicked in the presence of acrolein at CS-relevant concentrations, indicating that CS interferes with proper ER function, presumably due mainly to changes in cellular redox homeostasis. Since ER stress has been linked to diseases that are also related to CS exposure, these data are relevant in the discussion of a general molecular mechanism of CS-induced disease.

Characterization of Nrf2 activation and heme oxygenase-1 expression in NIH3T3 cells exposed to aqueous extracts of cigarette smoke.

Cigarette smoke (CS) is a complex chemical mixture estimated to be composed of up to 5000 different chemicals, many of which are prooxidant. Here we show that, at least in vitro, the cellular response designed to combat oxidative stress resulting from CS exposure is primarily controlled by the transcription factor Nrf2, a principal inducer of antioxidant and phase II-related genes. The prominent role of Nrf2 in the cellular response to CS is substantiated by the following observations: In NIH3T3 cells exposed to aqueous extracts of CS (i) Nrf2 is strongly stabilized and becomes detectable in nuclear extracts. (ii) Nuclear localization of Nrf2 coincides with increased DNA binding of a putative Nrf2/MafK heterodimer to its cognate cis-regulatory site, i.e., the antioxidant-responsive element (ARE). (iii) Studies on the regulatory elements of the oxidative stress-inducible gene heme oxygenase-1 (hmox1) using various hmox1 promoter/luciferase reporter constructs revealed that the strong CS-dependent expression of this gene is primarily governed by the distal enhancers 1 ("E1") and 2 ("E2"), which both contain three canonical ARE-like stress-responsive elements (StREs). Notably, depletion of Nrf2 levels caused by RNA interference significantly compromised CS-induced hmox1 promoter activation, based on the distinct Nrf2 sensitivity exhibited by E1 and E2. Finally, (iv) siRNA-dependent knock-down of Nrf2 completely abrogated CS-induced expression of phase II-related genes. Taken together, these results confirm the outstanding role of Nrf2 both in sensing (oxidant) stress and in orchestrating an efficient transcriptional response aimed at resolving the stressing conditions.

Growth suppression induced by downregulation of E6-AP expression in human papillomavirus-positive cancer cell lines depends on p53.

The ubiquitin-protein ligase E6-AP is utilized by the E6 oncoprotein of human papillomaviruses (HPVs) associated with cervical cancer to target the tumor suppressor p53 for degradation. Here, we report that downregulation of E6-AP expression by RNA interference results in both the accumulation of p53 and growth suppression of the HPV-positive cervical cancer cell lines HeLa and SiHa. In addition, HeLa cells, in which p53 expression was suppressed by RNA interference, are significantly less sensitive to the downregulation of E6-AP expression with respect to growth suppression than parental HeLa cells. These data indicate that the anti-growth-suppressive properties of E6-AP in HPV-positive cells depend on its ability to induce p53 degradation.

Involvement of the DNA repair protein hHR23 in p53 degradation.

The stability of the tumor suppressor protein p53 is regulated via the ubiquitin-proteasome-dependent proteolytic pathway. Like most substrates of this pathway, p53 is modified by the attachment of polyubiquitin chains prior to proteasome-mediated degradation. However, the mechanism(s) involved in the delivery of polyubiquitylated p53 molecules to the proteasome are currently unclear. Here, we show that the human DNA repair protein hHR23 binds to polyubiquitylated p53 via its carboxyl-terminal ubiquitin-associated (Uba) domain shielding p53 from deubiquitylation in vitro and in vivo. In addition, downregulation of hHR23 expression within cells by RNA interference results in accumulation of p53. Since the Ubl domain of hHR23 has been shown to interact with the 26S proteasome, we propose that hHR23 is intrinsically involved in the delivery of polyubiquitylated p53 molecules to the proteasome. In this model, the Uba domain of hHR23 binds to polyubiquitin chains formed on p53 and protects them from deubiquitylation, while the Ubl domain delivers the polyubiquitylated p53 molecules to the proteasome.

HdmX stimulates Hdm2-mediated ubiquitination and degradation of p53.

The RING finger proteins HdmX and Hdm2 share significant structural and functional similarity. Hdm2 is a member of the RING finger family of ubiquitin-protein ligases E3 and targets the tumor suppressor protein p53 for degradation. Although HdmX also binds to p53, HdmX does not induce p53 degradation. Moreover, HdmX has been reported to interfere with p53 degradation in overexpression experiments. To obtain insight into the mechanism by which HdmX interferes with p53 degradation, we studied the effect of HdmX on the E3 activity of Hdm2 in vitro. Surprisingly, this revealed that HdmX stimulates Hdm2-mediated ubiquitination of p53 and that HdmX facilitates ubiquitination of Hdm2 and vice versa. In addition, down-regulation of HdmX expression within cells results in the accumulation of both p53 and Hdm2. Because HdmX alone does not have appreciable E3 activity, these data indicate that HdmX acts as a stimulator, rather than as an inhibitor, of the E3 activity of Hdm2 and that, at least under certain conditions, HdmX is actively involved in the degradation of both p53 and Hdm2.

siRNA targeting of the viral E6 oncogene efficiently kills human papillomavirus-positive cancer cells.

The targeted inhibition of antiapoptotic factors in tumour cells may provide a rational approach towards the development of novel anticancer therapies. Using human papillomavirus (HPV)-transformed cells as a model system, we investigated if RNA interference (RNAi)-mediated gene silencing can be employed in order to overcome the apoptosis resistance of cancer cells. We found that both vector-borne and synthetic small interfering (si)RNAs, specifically directed against the antiapoptotic HPV E6 oncogene, restored dormant tumour suppressor pathways in HPV-positive cancer cells that are otherwise inactive in the presence of E6. This ultimately resulted in massive apoptotic cell death, selectively in HPV-positive tumour cells. These findings show that RNAi provides a powerful molecular strategy to inactivate intracellular E6 function efficiently. Moreover, they define E6 as a most promising therapeutic target to eliminate HPV-positive tumour cells specifically by RNAi. Thus, by sequence-specific targeting of antiapoptotic genes, siRNAs may be developed into novel therapeutics that can efficiently correct the apoptosis deficiency of cancer cells.