Integration of pharmacokinetic and NRF2 system biology models to describe reactive oxygen species production and subsequent glutathione depletion in liver microfluidic biochips after flutamide exposure.

We present a systems biology analysis of rat primary hepatocytes response after exposure to 10 μM and 100 μM flutamide in liver microfluidic biochips. We coupled an in vitro pharmacokinetic (PK) model of flutamide to a system biology model of its reactive oxygen species (ROS) production and scavenging by the Nrf2 regulated glutathione production. The PK model was calibrated using data on flutamide kinetics, hydroxyflutamide and glutathione conjugates formation in microfluidic conditions. The parameters of Nrf2-related gene activities and the subsequent glutathione depletion were calibrated using microarray data from our microfluidic experiments and literature information. Following a 10 μM flutamide exposure, the model predicted a recovery time to baseline levels of glutathione (GSH) and ROS in agreement with our experimental observations. At 100 μM, the model predicted that metabolism saturation led to an important accumulation of flutamide in cells, a high ROS production and complete GSH depletion. The high levels of ROS predicted were consistent with the necrotic switch observed by transcriptomics, and the high cell mortality we had experimentally observed. The model predicted a transition between recoverable GSH depletion and deep GSH depletion at about 12.5 μM of flutamide (single perfusion exposure). Our work shows that in vitro biochip experiments can provide supporting information for complex in silico modeling including data from extra cellular and intra cellular levels. We believe that this approach can be an efficient strategy for a global integrated methodology in predictive toxicology.