High content omic techniques in combination with stable human in vitro cell culture systems have the potential to improve on current pre-clinical safety regimes by providing detailed mechanistic information of altered cellular processes. Here we investigated the added benefit of integrating transcriptomics, proteomics and metabolomics together with pharmacokinetics for drug testing regimes. Cultured human renal epithelial cells (RPTEC/TERT1) were exposed to the nephrotoxin Cyclosporine A (CsA) at therapeutic and supratherapeutic concentrations for 14days. CsA was quantified in supernatants and cellular lysates by LC-MS/MS for kinetic modeling. There was a rapid cellular uptake and accumulation of CsA, with a non-linear relationship between intracellular and applied concentrations. CsA at 15μM induced mitochondrial disturbances and activation of the Nrf2-oxidative-damage and the unfolded protein-response pathways. All three omic streams provided complementary information, especially pertaining to Nrf2 and ATF4 activation. No stress induction was detected with 5μM CsA; however, both concentrations resulted in a maximal secretion of cyclophilin B. The study demonstrates for the first time that CsA-induced stress is not directly linked to its primary pharmacology. In addition we demonstrate the power of integrated omics for the elucidation of signaling cascades brought about by compound induced cell stress.
BackgroundIncorporation of omic data streams for building improved systems biology models has great potential for improving their predictions of biological outcomes. We have recently shown that cyclosporine A (CsA) strongly activates the nuclear factor (erythroid-derived 2)-like 2 pathway (Nrf2) in renal proximal tubular epithelial cells (RPTECs) exposed in vitro. We present here a quantitative calibration of a differential equation model of the Nrf2 pathway with a subset of the omics data we collected.ResultsIn vitro pharmacokinetic data on CsA exchange between cells, culture medium and vial walls, and data on the time course of omics markers in response to CsA exposure were reasonably well fitted with a coupled PK-systems biology model. Posterior statistical distributions of the model parameter values were obtained by Markov chain Monte Carlo sampling in a Bayesian framework. A complex cyclic pattern of ROS production and control emerged at 5 μM CsA repeated exposure. Plateau responses were found at 15 μM exposures. Shortly above those exposure levels, the model predicts a disproportionate increase in cellular ROS quantity which is consistent with an in vitro EC50 of about 40 μM for CsA in RPTECs.ConclusionsThe model proposed can be used to analyze and predict cellular response to oxidative stress, provided sufficient data to set its parameters to cell-specific values. Omics data can be used to that effect in a Bayesian statistical framework which retains prior information about the likely parameter values.
We have integrated in vitro and in silico information to investigate acetaminophen (APAP) and its metabolite N-acetyl-p-benzoquinone imine (NAPQI) toxicity in liver biochip. In previous works, we observed higher cytotoxicity of HepG2/C3a cultivated in biochips when exposed to 1 mM of APAP for 72 h as compared to Petri cultures. We complete our investigation with the present in silico approach to extend the mechanistic interpretation of the intracellular kinetics of the toxicity process. For that purpose, we propose a mathematical model based on the coupling of a drug pharmacokinetic model (PK) with a systemic biology model (SB) describing the reactive oxygen species (ROS) production by NAPQI and the subsequent glutathione (GSH) depletion. The SB model was parameterized using (i) transcriptomic data, (ii) qualitative results of time lapses ROS fluorescent curves for both control and 1-mM APAP-treated experiments, and (iii) additional GSH literature data. The PK model was parameterized (i) using the in vitro kinetic data (at 160 μM, 1 mM, 10 mM), (ii) using the parameters resulting from a physiologically based pharmacokinetic (PBPK) literature model for APAP, and (iii) by literature data describing NAPQI formation. The PK-SB model predicted a ROS increase and GSH depletion due to the NAPQI formation. The transition from a detoxification phase and NAPQI and ROS accumulation was predicted for a NAPQI concentration ranging between 0.025 and 0.25 μM in the cytosol. In parallel, we performed a dose response analysis in biochips that shows a reduction of the final hepatic cell number appeared in agreement with the time and doses associated with the switch of the NAPQI detoxification/accumulation. As a result, we were able to correlate in vitro extracellular APAP exposures with an intracellular in silico ROS accumulation using an integration of a coupled mathematical and experimental liver on chip approach.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.