Engineering topics

Shankaran, Harish

Publications and source records attributed to Shankaran, Harish.

7 records

DOE OSTI2018

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

High-content imaging of embryonic 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 so far 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. Here, we have developed and parallelized a workflow to calculate benchmark dose (BMD) values from dynamic, high-content zebrafish behavior data that scales for high-throughput chemical screening. To appropriately capture the dynamics of zebrafish response from light to dark stimulus, we summarized time-dependent data using both area under the curve and 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 BMD 10 values of 1,060 ToxCast chemicals for 24 zebrafish endpoints, including behavior, mortality and morphology. The BMD 10 values provided better precision and separation than LELs for clustering chemicals since they were derived from models that best-fit their concentration-response curves. Clustering analysis revealed behavioral signatures as the most sensitive endpoints. High concordance in chemical activity was observed between ToxCast assays and zebrafish behavioral endpoints, however ToxPi analysis indicated that rankings based on in vitro data was not a reliable predictor of in vivo rankings for lower potency chemicals. The methods developed in this work will enable the use of high-content zebrafish behavioral screening data for BMD analysis in toxicological hazard assessment.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

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, JNK (SAPK), the PI3K/AKT, as well as various transcriptional regulators such as STAT3. When co-expressed, first three members of HER family (HER1-3) can form homo- and hetero-dimers. Based on the considerable evidence which suggest that every receptor dimer activates intracellular signaling pathways differentially, we hypothesized that the HER dimerization pattern is a better predictor of downstream signaling than the total receptor activation levels. We validated our hypothesis using a combination of model-based analysis to quantify the HER dimerization patterns and multi-factorial experiments where HER dimerization patterns and signaling crosstalk were rationally perturbed. We have measured the activation of HER1-3 receptors and of the sentinel signaling proteins ERK, AKT, p38, 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. Our analysis using multiple ways of clustering the activation data has confirmed that the HER receptor dimer is a better predictor of the signaling through p38, ERK and AKT pathways than the total HER receptor expression and activation levels. Targeted inhibition studies to identify the causal effects allowed us to obtain a consensus regulatory interaction model, which revealed that STAT3 occupies a central role in the crosstalk between the studied pathways.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

ERK Oscillation-Dependent Gene Expression Patterns and Deregulation by Stress-Response

Studies were undertaken to determine whether 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 non-oscillating 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 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 non-oscillating cells.

carcinogenesis↗

DOE OSTI2013

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. - Highlights: • We develop a Fentanyl PBPK model for relating external dose to internal levels. • We calibrate the model to oral and inhalation exposures using > 50 human datasets. • Model predictions are in good agreement with the available pharmacokinetic data. • The model can be used for extrapolating across routes, doses and exposure durations. • We illustrate how the model can be used for developing Provisional Advisory Levels.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

Model-based Analysis of HER Activation in Cells Co-Expressing EGFR, HER2 and HER3.

The HER/ErbB family of receptor tyrosine kinases drive 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, 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 epithelial cells lines with known HER expression levels. Simulations over the range of expression levels seen in various cell lines indicate that: i) EGFR phosphorylation is driven by HER1/1 and HER1/2 dimers, and not HER1/3 dimers, ii) HER1/2 and HER2/3 dimers both contribute significantly to HER2 activation with the EGFR expression level determining the relative importance of these species, and iii) the HER2/3 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.

ErbB↗

DOE OSTI2013

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 both space 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), 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.

von Neubeck, Claere↗

DOE OSTI2012

Modeling Dynamic Regulatory Processes in Stroke.

The ability to examine in silico the behavior of biological systems can greatly accelerate the pace of discovery in disease pathologies, such as stroke, where in vivo experimentation is lengthy and costly. In this paper we describe an approach to in silico examination of blood genomic responses 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) relating regulators and functional clusters from the data. These ODEs were used to develop dynamic models that simulate the expression of regulated functional clusters using system dynamics as the modeling paradigm. The dynamic model has the considerable advantage of only requiring an initial starting state, and does not require measurement of regulatory influences at each time point in order to make accurate predictions. The manipulation of input model parameters, such as changing the magnitude of gene expression, made it possible to assess the behavior of the networks through time under varying conditions. We report that an optimized dynamic model can provide accurate predictions of overall system behavior under several different preconditioning paradigms.

dynamic modeling↗