A Surrogate Matrix-Based Approach Toward Multiplexed Quantitation of an sGC Stimulator and cGMP in Ocular Tissue and Plasma.

A liquid chromatography-tandem mass spectrometry assay was developed and qualified for the multiplexed quantitation of a small molecule stimulator of soluble guanylate cyclase (sGC) and its target engagement biomarker, 3',5'-cyclic guanosine monophosphate (cGMP), in ocular tissues and plasma from a single surrogate matrix calibration curve. A surrogate matrix approach was used in this assay due to the limited quantities of blank ocular matrices in a discovery research setting. After optimization, the assay showed high accuracy, precision, and recovery as well as parallelism between the surrogate matrix and the biological matrices (rabbit plasma, vitreous, and retina-choroid). This assay provided pharmacokinetic and target engagement data after intravitreal administration of the sGC stimulator. The nitric oxide-sGC-cGMP pathway is a potential target to address glaucoma. Increasing sGC-mediated production of cGMP could improve aqueous humor outflow and ocular blood flow. The sGC stimulator showed dose-dependent exposure in rabbit vitreous, retina-choroid, and plasma. The cGMP exhibited a delayed yet sustained increase in vitreous humor but not retina-choroid. Multiplexed measurement of both pharmacokinetic and target engagement analytes reduced animal usage and provided improved context for interpreting PK and PD relationships.

Investigating nephrotoxicity of polymyxin derivatives by mapping renal distribution using mass spectrometry imaging.

Colistin and polymyxin B are effective treatment options for Gram-negative resistant bacteria but are used as last-line therapy due to their dose-limiting nephrotoxicity. A critical factor in developing safer polymyxin analogues is understanding accumulation of the drugs and their metabolites, which is currently limited due to the lack of effective techniques for analysis of these challenging molecules. Mass spectrometry imaging (MSI) allows direct detection of targets (drugs, metabolites, and endogenous compounds) from tissue sections. The presented study exemplifies the utility of MSI by measuring the distribution of polymyxin B1, colistin, and polymyxin B nonapeptide (PMBN) within dosed rat kidney tissue sections. The label-free MSI analysis revealed that the nephrotoxic compounds (polymyxin B1 and colistin) preferentially accumulated in the renal cortical region. The less nephrotoxic analogue, polymyxin B nonapeptide, was more uniformly distributed throughout the kidney. In addition, metabolites of the dosed compounds were detected by MSI. Kidney homogenates were analyzed using LC/MS/MS to determine total drug exposure and for metabolite identification. To our knowledge, this is the first time such techniques have been utilized to measure the distribution of polymyxin drugs and their metabolites. By simultaneously detecting the distribution of drug and drug metabolites, MSI offers a powerful alternative to tissue homogenization analysis and label or antibody-based imaging.

ERK oscillation-dependent gene expression patterns and deregulation by stress response.

Studies were undertaken to determine whether extracellular signal regulated kinase (ERK) oscillations regulate a unique subset of genes in human keratinocytes and subsequently whether the p38 stress response inhibits ERK oscillations. A DNA microarray identified many genes that were unique to ERK oscillations, and network reconstruction predicted an important role for the mediator complex subunit 1 (MED1) node in mediating ERK oscillation-dependent gene expression. Increased ERK-dependent phosphorylation of MED1 was observed in oscillating cells compared to nonoscillating counterparts as validation. Treatment of keratinocytes with a p38 inhibitor (SB203580) increased ERK oscillation amplitudes and MED1 and phospho-MED1 protein levels. Bromate is a probable human carcinogen that activates p38. Bromate inhibited ERK oscillations in human keratinocytes and JB6 cells and induced an increase in phospho-p38 and a decrease in phospho-MED1 protein levels. Treatment of normal rat kidney cells and primary salivary gland epithelial cells with bromate decreased phospho-MED1 levels in a reversible fashion upon treatment with p38 inhibitors (SB202190; SB203580). Our results indicate that oscillatory behavior in the ERK pathway alters homeostatic gene regulation patterns and that the cellular response to perturbation may manifest differently in oscillating vs nonoscillating cells.

Model-based analysis of HER activation in cells co-expressing EGFR, HER2 and HER3.

The HER/ErbB family of receptor tyrosine kinases drives critical responses in normal physiology and cancer, and the expression levels of the various HER receptors are critical determinants of clinical outcomes. HER activation is driven by the formation of various dimer complexes between members of this receptor family. The HER dimer types can have differential effects on downstream signaling and phenotypic outcomes. We constructed an integrated mathematical model of HER activation, and trafficking to quantitatively link receptor expression levels to dimerization and activation. We parameterized the model with a comprehensive set of HER phosphorylation and abundance data collected in a panel of human mammary epithelial cells expressing varying levels of EGFR/HER1, HER2 and HER3. Although parameter estimation yielded multiple solutions, predictions for dimer phosphorylation were in agreement with each other. We validated the model using experiments where pertuzumab was used to block HER2 dimerization. We used the model to predict HER dimerization and activation patterns in a panel of human mammary epithelial cells lines with known HER expression levels in response to stimulations with ligands EGF and HRG. Simulations over the range of expression levels seen in various cell lines indicate that: i) EGFR phosphorylation is driven by HER1-HER1 and HER1-HER2 dimers, and not HER1-HER3 dimers, ii) HER1-HER2 and HER2-HER3 dimers both contribute significantly to HER2 activation with the EGFR expression level determining the relative importance of these species, and iii) the HER2-HER3 dimer is largely responsible for HER3 activation. The model can be used to predict phosphorylated dimer levels for any given HER expression profile. This information in turn can be used to quantify the potencies of the various HER dimers, and can potentially inform personalized therapeutic approaches.

Physiologically-based pharmacokinetic model for Fentanyl in support of the development of Provisional Advisory Levels.

Provisional Advisory Levels (PALs) are tiered exposure limits for toxic chemicals in air and drinking water that are developed to assist in emergency responses. Physiologically-based pharmacokinetic (PBPK) modeling can support this process by enabling extrapolations across doses, and exposure routes, thereby addressing gaps in the available toxicity data. Here, we describe the development of a PBPK model for Fentanyl - a synthetic opioid used clinically for pain management - to support the establishment of PALs. Starting from an existing model for intravenous Fentanyl, we first optimized distribution and clearance parameters using several additional IV datasets. We then calibrated the model using pharmacokinetic data for various formulations, and determined the absorbed fraction, F, and time taken for the absorbed amount to reach 90% of its final value, t90. For aerosolized pulmonary Fentanyl, F=1 and t90<1 min indicating complete and rapid absorption. The F value ranged from 0.35 to 0.74 for oral and various transmucosal routes. Oral Fentanyl was absorbed the slowest (t90~300 min); the absorption of intranasal Fentanyl was relatively rapid (t90~20-40 min); and the various oral transmucosal routes had intermediate absorption rates (t90~160-300 min). Based on these results, for inhalation exposures, we assumed that all of the Fentanyl inhaled from the air during each breath directly, and instantaneously enters the arterial circulation. We present model predictions of Fentanyl blood concentrations in oral and inhalation scenarios relevant for PAL development, and provide an analytical expression that can be used to extrapolate between oral and inhalation routes for the derivation of PALs.

Modeling dynamic regulatory processes in stroke.

The ability to examine the behavior of biological systems in silico has the potential to greatly accelerate the pace of discovery in diseases, such as stroke, where in vivo analysis is time intensive and costly. In this paper we describe an approach for in silico examination of responses of the blood transcriptome to neuroprotective agents and subsequent stroke through the development of dynamic models of the regulatory processes observed in the experimental gene expression data. First, we identified functional gene clusters from these data. Next, we derived ordinary differential equations (ODEs) from the data relating these functional clusters to each other in terms of their regulatory influence on one another. Dynamic models were developed by coupling these ODEs into a model that simulates the expression of regulated functional clusters. By changing the magnitude of gene expression in the initial input state it was possible to assess the behavior of the networks through time under varying conditions since the dynamic model only requires an initial starting state, and does not require measurement of regulatory influences at each time point in order to make accurate predictions. We discuss the implications of our models on neuroprotection in stroke, explore the limitations of the approach, and report that an optimized dynamic model can provide accurate predictions of overall system behavior under several different neuroprotective paradigms.

Using measures of metabolic flux to align screening and clinical development: Avoiding pitfalls to enable translational studies.

Screening campaigns, especially those aimed at modulating enzyme activity, often rely on measuring substrate→product conversions. Unfortunately, the presence of endogenous substrates and/or products can limit one's ability to measure conversions. As well, coupled detection systems, often used to facilitate optical readouts, are subject to interference. Stable isotope labeled substrates can overcome background contamination and yield a direct readout of enzyme activity. Not only can isotope kinetic assays enable early screening, but they can also be used to follow hit progression in translational (pre)clinical studies. Herein, we consider a case study surrounding lipid biology to exemplify how metabolic flux analyses can connect stages of drug development, caveats are highlighted to ensure reliable data interpretations. For example, when measuring enzyme activity in early biochemical screening it may be enough to quantify the formation of a labeled product. In contrast, cell-based and in vivo studies must account for variable exposure to a labeled substrate (or precursor) which occurs via tracer dilution and/or isotopic exchange. Strategies are discussed to correct for these complications. We believe that measures of metabolic flux can help connect structure-activity relationships with pharmacodynamic mechanisms of action and determine whether mechanistically differentiated biophysical interactions lead to physiologically relevant outcomes. Adoption of this logic may allow research programs to (i) build a critical bridge between primary screening and (pre)clinical development, (ii) elucidate biology in parallel with screening and (iii) suggest a strategy aimed at in vivo biomarker development.

Inhibition of ERK oscillations by ionizing radiation and reactive oxygen species.

The shuttling of activated protein kinases between the cytoplasm and nucleus is an essential feature of normal growth factor signaling cascades. Here we demonstrate that transforming growth factor alpha (TGFα) induces oscillations in extracellular signal regulated kinase (ERK) cytoplasmic-nuclear translocations in human keratinocytes. TGFα-dependent ERK oscillations mediated through the epidermal growth factor receptor (EGFR) are inhibited by low dose X-irradiation (10 cGy) and low concentrations of hydrogen peroxide (0.32-3.26 µM H(2)O(2)) used as a model reactive oxygen species (ROS). A fluorescent indicator dye (H2-DCFDA) was used to measure cellular ROS levels following X-irradiation, 12-O-tetradecanoyl phorbol-13-acetate (TPA) and H(2)O(2). X-irradiation did not generate significant ROS production while 0.32 µM H(2)O(2) and TPA induced significant increases in ROS levels with H(2)O(2) > TPA. TPA alone induced transactivation of the EGFR but did not induce ERK oscillations. TPA as a cotreatment did not inhibit TGFα-stimulated ERK oscillations but qualitatively altered TGFα-dependent ERK oscillation characteristics (amplitude, time-period). Collectively, these observations demonstrate that TGFα-induced ERK oscillations are inhibited by ionizing radiation/ROS and perturbed by epigenetic carcinogen in human keratinocytes.

Rapid and sustained nuclear-cytoplasmic ERK oscillations induced by epidermal growth factor.

Although the ERK pathway has a central role in the response of cells to growth factors, its regulatory structure and dynamics are incompletely understood. To investigate ERK activation in real time, we expressed an ERK-GFP fusion protein in human mammary epithelial cells. On EGF stimulation, we observed sustained oscillations of the ERK-GFP fusion protein between the nucleus and cytoplasm with a periodicity of approximately 15 min. The oscillations were persistent (>45 cycles), independent of cell cycle phase, and were highly dependent on cell density, essentially disappearing at confluency. Oscillations occurred even at ligand doses that elicited very low levels of ERK phosphorylation, and could be detected biochemically in both transfected and nontransfected cells. Mathematical modeling revealed that negative feedback from phosphorylated ERK to the cascade input was necessary to match the robustness of the oscillation characteristics observed over a broad range of ligand concentrations. Our characterization of single-cell ERK dynamics provides a quantitative foundation for understanding the regulatory structure of this signaling cascade.

HER/ErbB receptor interactions and signaling patterns in human mammary epithelial cells.

BACKGROUND: Knowledge about signaling pathways is typically compiled based on data gathered using different cell lines. This approach implicitly assumes that the cell line dependence is not important. However, different cell lines do not always respond to a particular stimulus in the same way, and lack of coherent data collected from closely related cellular systems can be detrimental to the efforts to understand the regulation of biological processes. To address this issue, we created a clone library of human mammary epithelial (HME) cells that expresses different levels of HER2 and HER3 receptors in combination with endogenous EGFR/HER1. Using our clone library, we have quantified the receptor activation patterns and systematically tested the validity of the existing hypotheses about the interaction patterns between HER1-3 receptors. RESULTS: Our study identified HER2 as the dominant dimerization partner for both EGFR and HER3. Contrary to earlier suggestions, we find that lateral interactions with HER2 do not lead to strong transactivation between EGFR and HER3, i.e., EGFR activation and HER3 activation are only weakly linked in HME cells. We also find that observed weak transactivation is uni-directional where stimulation of EGFR leads to HER3 activation whereas HER3 stimulation does not activate the EGFR. Repeating our experiments at lower cell confluency established that cell confluency is not a major factor in the observed interaction patterns. We have also quantified the dependence of the kinetics of Erk and Akt activation on different HER receptors. We found that HER3 signaling makes the strongest contribution to Akt activation and that, stimulation of either EGFR or HER3 leads to significant Erk activation. CONCLUSION: Our study shows that clone cell libraries can be a powerful resource in systems biology research by making it possible to differentiate between various hypotheses in a consistent cellular background. Using our constructed clone library we profiled the cell signaling patterns to establish the role of HER2 in the crosstalk between EGFR and HER3 receptors in HME cells. Our results for HME cells show that the weak linkage between EGFR and HER3 pathways can lead to distinct downstream cellular signaling patterns in response to the ligands of these two receptors.

Time-dependent behavioral data from zebrafish reveals novel signatures of chemical toxicity using point of departure analysis.

High-content imaging of larval zebrafish behavior can be used as a screening approach to rapidly evaluate the relative potential for chemicals to cause toxicity. However, most statistical methods applied to these data transform movement values to incidence-based "hits" and calculate lowest effect levels (LELs), which loses individual fish resolution of behavior and defies hazard ranking due to reliance on applied dose levels. We developed a parallelizable workflow to calculate benchmark dose (BMD) values from dynamic, high-content zebrafish behavior data that scales for high-throughput chemical screening. To capture the zebrafish movement response from light to dark stimulus, we summarized time-dependent data using both area under the curve and the immediate change at the transition point into two novel metrics that characterized abnormal behavior as a function of chemical concentration. The BMD workflow was applied to calculate BMD10 values of 1,060 ToxCast chemicals for 24 zebrafish endpoints, including behavior, mortality and morphology. The BMD10 values provided better precision and separation than LELs for clustering chemicals since they were derived from models that best-fit their concentration-response curves. Analysis of BMD10 values revealed behavioral signatures as the most sensitive endpoints. High concordance in chemical activity was observed between ToxCast in vitro data and zebrafish in vivo behavioral data, however ToxPi analysis indicated that rankings based on in vitro data were not a reliable predictor of in vivo rankings for lower potency chemicals. This analysis method will enable the use of high-content zebrafish behavioral screening data for BMD analysis in toxicological hazard assessment.

A Systems Toxicology Approach for the Prediction of Kidney Toxicity and Its Mechanisms In Vitro.

The failure to predict kidney toxicity of new chemical entities early in the development process before they reach humans remains a critical issue. Here, we used primary human kidney cells and applied a systems biology approach that combines multidimensional datasets and machine learning to identify biomarkers that not only predict nephrotoxic compounds but also provide hints toward their mechanism of toxicity. Gene expression and high-content imaging-derived phenotypical data from 46 diverse kidney toxicants were analyzed using Random Forest machine learning. Imaging features capturing changes in cell morphology and nucleus texture along with mRNA levels of HMOX1 and SQSTM1 were identified as the most powerful predictors of toxicity. These biomarkers were validated by their ability to accurately predict kidney toxicity of four out of six candidate therapeutics that exhibited toxicity only in late stage preclinical/clinical studies. Network analysis of similarities in toxic phenotypes was performed based on live-cell high-content image analysis at seven time points. Using compounds with known mechanism as reference, we could infer potential mechanisms of toxicity of candidate therapeutics. In summary, we report an approach to generate a multidimensional biomarker panel for mechanistic de-risking and prediction of kidney toxicity in in vitro for new therapeutic candidates and chemical entities.

Pharmacokinetic Properties of Humanized IgG1 and IgG4 Antibodies in Preclinical Species: Translational Evaluation.

Assessment of the factors that regulate antibody exposure-response relationships in the relevant animal models is critical for the design of successful translational strategies from discovery to the clinic. Depending on the specific clinical indication, preclinical development paradigms may require that the efficacy or dosing-related attributes for the existing antibody be assessed in various species when cross-reactivity of the lead antibody to the intended species is justified. Additionally, with the success of monoclonal antibodies for management of various human conditions, a parallel interest in therapeutic use of these novel modalities in various veterinary species has followed. The protective role of neonatal Fc receptor (FcRn) in regulation of IgG homeostasis and clearance is now well recognized and the "nonspecific clearance" of antibodies through bone marrow-derived phagocytic and vascular endothelial cells (via lysosomal processes) is modulated by interactions with FcRn receptors. In this study, we have attempted to examine the PK properties of human IgG antibodies in dog and monkey. These studies establish a translational framework for evaluation of IgG antibody PK properties across species.

Quantifying the effects of co-expressing EGFR and HER2 on HER activation and trafficking.

The human epidermal growth factor receptor (HER) system is an intricately regulated system that plays critical roles in development and tumorigenesis. Here, we apply integrated experimentation and modeling to analyze HER receptor activation in a panel of cell lines expressing endogenous levels of EGFR/HER1 and different levels of HER2. A mathematical model that includes the fundamental processes involved in receptor activation and trafficking was used to fit the experimental data, and values of the independent parameters for active receptor dimer formation affinities, trafficking rates and relative phosphorylation levels were estimated. Obtained parameter values quantitatively support the existing ideas on the effect of HER2 on EGFR dynamics, and enable us to predict the abundances of various phosphorylated receptor dimers in the cell lines.

Cell surface receptors for signal transduction and ligand transport: a design principles study.

Receptors constitute the interface of cells to their external environment. These molecules bind specific ligands involved in multiple processes, such as signal transduction and nutrient transport. Although a variety of cell surface receptors undergo endocytosis, the systems-level design principles that govern the evolution of receptor trafficking dynamics are far from fully understood. We have constructed a generalized mathematical model of receptor-ligand binding and internalization to understand how receptor internalization dynamics encodes receptor function and regulation. A given signaling or transport receptor system represents a particular implementation of this module with a specific set of kinetic parameters. Parametric analysis of the response of receptor systems to ligand inputs reveals that receptor systems can be characterized as being: i) avidity-controlled where the response control depends primarily on the extracellular ligand capture efficiency, ii) consumption-controlled where the ability to internalize surface-bound ligand is the primary control parameter, and iii) dual-sensitivity where both the avidity and consumption parameters are important. We show that the transferrin and low-density lipoprotein receptors are avidity-controlled, the vitellogenin receptor is consumption-controlled, and the epidermal growth factor receptor is a dual-sensitivity receptor. Significantly, we show that ligand-induced endocytosis is a mechanism to enhance the accuracy of signaling receptors rather than merely serving to attenuate signaling. Our analysis reveals that the location of a receptor system in the avidity-consumption parameter space can be used to understand both its function and its regulation.

Multiscale Mathematical Model of Drug-Induced Proximal Tubule Injury: Linking Urinary Biomarkers to Epithelial Cell Injury and Renal Dysfunction.

Drug-induced nephrotoxicity is a major cause of acute kidney injury, and thus detecting the potential for nephrotoxicity early in the drug development process is critical. Various urinary biomarkers exhibit different patterns following drug-induced injury, which may provide greater information than traditional biomarkers like serum creatinine. In this study, we developed a multiscale quantitative systems pharmacology model relating drug exposure to proximal tubule (PT) epithelial cell injury and subsequently to expression of multiple urinary biomarkers and organ-level functional changes. We utilized urinary kidney injury molecule-1 (Kim-1), alpha glutathione S-transferase, albumin (αGST), glucose, and urine volume time profiles as well as serum creatinine and histopathology data obtained from rats treated with the nephrotoxicant cisplatin to develop the model. Although the model was developed using single-dose response to cisplatin, the model predicted the serum creatinine response to multidose cisplatin regimens. Further, using only the urinary Kim-1 response to gentamicin (a nephrotoxicant with a distinctly different injury time course than cisplatin), the model detected and predicted mild to moderate PT injury, as confirmed with histopathology, even when serum creatinine was unchanged. Thus, the model is generalizable, and can be used to deconvolute the underlying degree and time course of drug-induced PT injury and renal dysfunction from a small number of urinary biomarkers, and may provide a tool to determine optimal dosing regimens that minimize renal injury.

Systems Pharmacology Model of Gastrointestinal Damage Predicts Species Differences and Optimizes Clinical Dosing Schedules.

Gastrointestinal (GI) adverse events (AEs) are frequently dose limiting for oncology agents, requiring extensive clinical testing of alternative schedules to identify optimal dosing regimens. Here, we develop a translational mathematical model to predict these clinical AEs starting from preclinical GI toxicity data. The model structure incorporates known biology and includes stem cells, daughter cells, and enterocytes. Published data, including cellular numbers and division times, informed the system parameters for humans and rats. The drug-specific parameters were informed with preclinical histopathology data from rats treated with irinotecan. The model fit the rodent irinotecan-induced pathology changes well. The predicted time course of enterocyte loss in patients treated with weekly doses matched observed AE profiles. The model also correctly predicts a lower level of AEs for every 3 weeks (Q3W), as compared to the weekly schedule.

Integrated analysis reveals that STAT3 is central to the crosstalk between HER/ErbB receptor signaling pathways in human mammary epithelial cells.

Human epidermal growth factor receptors (HER, also known as ErbB) drive cellular proliferation, pro-survival and stress responses by activating several downstream kinases, in particular ERK, p38 MAPK, JNK (SAPK), the PI3K/AKT, as well as various transcriptional regulators such as STAT3. When co-expressed, the first three members of HER family (HER1-3) can form homo- and hetero-dimers, and there is considerable evidence suggesting that the receptor dimers differentially activate intracellular signaling pathways. To better understand the interactions in this system, we pursued multi-factorial experiments where HER dimerization patterns and signaling pathways were rationally perturbed. We measured the activation of HER1-3 receptors and of the sentinel signaling proteins ERK, AKT, p38 MAPK, JNK, STAT3 as a function of time in a panel of human mammary epithelial (HME) cells expressing different levels of HER1-3 stimulated with various ligand combinations. We hypothesized that the HER dimerization pattern is a better predictor of downstream signaling than the total receptor activation levels. We validated this hypothesis using a combination of model-based analysis to quantify the HER dimerization patterns, and by clustering the activation data in multiple ways to confirm that the HER receptor dimer is a better predictor of the signaling through p38 MAPK, ERK and AKT pathways than the total HER receptor expression and activation levels. We then pursued combinatorial inhibition studies to identify the causal regulatory interactions between sentinel signaling proteins. Quantitative analysis of the collected data using the modular response analysis (MRA) and its Bayesian Variable Selection Algorithm (BVSA) version allowed us to obtain a consensus regulatory interaction model, which revealed that STAT3 occupies a central role in the crosstalk between the studied pathways in HME cells. Results of the BVSA/MRA and cluster analysis were in agreement with each other.

Integrated experimental and computational approach to understand the effects of heavy ion radiation on skin homeostasis.

The effects of low dose high linear energy transfer (LET) radiation on human health are of concern for space, occupational, and clinical exposures. As epidemiological data for such radiation exposures are scarce for making relevant predictions, we need to understand the mechanism of response especially in normal tissues. Our objective here is to understand the effects of heavy ion radiation on tissue homeostasis in a realistic model system. Towards this end, we exposed an in vitro three dimensional skin equivalent to low fluences of neon (Ne) ions (300 MeV u(-1)), and determined the differentiation profile as a function of time following exposure using immunohistochemistry. We found that Ne ion exposures resulted in transient increases in the tissue regions expressing the differentiation markers keratin 10, and filaggrin, and more subtle time-dependent effects on the number of basal cells in the epidermis. We analyzed the data using a mathematical model of the skin equivalent, to quantify the effect of radiation on cell proliferation and differentiation. The agent-based mathematical model for the epidermal layer treats the epidermis as a collection of heterogeneous cell types with different proliferation-differentiation properties. We obtained model parameters from the literature where available, and calibrated the unknown parameters to match the observed properties in unirradiated skin. We then used the model to rigorously examine alternate hypotheses regarding the effects of high LET radiation on the tissue. Our analysis indicates that Ne ion exposures induce rapid, but transient, changes in cell division, differentiation and proliferation. We have validated the modeling results by histology and quantitative reverse transcription polymerase chain reaction (qRT-PCR). The integrated approach presented here can be used as a general framework to understand the responses of multicellular systems, and can be adapted to other epithelial tissues.

Using imaging methods to interrogate radiation-induced cell signaling.

There is increasing emphasis on the use of systems biology approaches to define radiation-induced responses in cells and tissues. Such approaches frequently rely on global screening using various high throughput 'omics' platforms. Although these methods are ideal for obtaining an unbiased overview of cellular responses, they often cannot reflect the inherent heterogeneity of the system or provide detailed spatial information. Additionally, performing such studies with multiple sampling time points can be prohibitively expensive. Imaging provides a complementary method with high spatial and temporal resolution capable of following the dynamics of signaling processes. In this review, we utilize specific examples to illustrate how imaging approaches have furthered our understanding of radiation-induced cellular signaling. Particular emphasis is placed on protein colocalization, and oscillatory and transient signaling dynamics.