Development of a quantitative model of pregnane X receptor (PXR) mediated xenobiotic metabolizing enzyme induction.

The pregnane X receptor plays an integral role in the regulation of hepatic metabolism. It has been shown to regulate CYP3A4, which is the most abundant cytochrome P450 in the human liver. With its large and flexible ligand-binding domain, PXR can be activated by an enormous range of relatively small, hydrophobic, exogenous compounds. Upon activation, PXR partners with the retinoid X receptor (RXR) to form a heterodimer. The newly formed heterodimer binds to an appropriate DNA response element, causing increased transcription. This leads to an induction in the level of CYP3A4. These mechanistic steps are included into a biologically-based mathematical model. The quantitative model predicts fold level inductions of CYP3A4 mRNA and protein in response to PXR activation. Model parameter values have been taken from literature when appropriate. Unknown parameter values are estimated by optimizing the model results to published in vivo and in vitro data sets. A sensitivity analysis is performed to evaluate the model structure and identify future data needs which would be critical to revising the model.

Development of a quantitative model incorporating key events in a hepatotoxic mode of action to predict tumor incidence.

Biologically based dose-response (BBDR) modeling of environmental pollutants can be utilized to inform the mode of action (MOA) by which compounds elicit adverse health effects. Chemicals that produce tumors are typically labeled as either genotoxic or nongenotoxic. Though both the genotoxic and the nongenotoxic MOA may be operative as a function of dose, it is important to note that the label informs but does not define a MOA. One commonly proposed MOA for nongenotoxic carcinogens is characterized by the key events cytotoxicity and regenerative proliferation. The increased division rate associated with such proliferation can cause an increase in the probability of mutations, which may result in tumor formation. We included these steps in a generalized computational pharmacodynamic (PD) model incorporating cytotoxicity as a MOA for three carcinogens (chloroform, CHCl(3); carbon tetrachloride, CCL(4); and N,N-dimethylformamide, DMF). For each compound, the BBDR model is composed of a chemical-specific physiologically based pharmacokinetic model linked to a PD model of cytotoxicity and cellular proliferation. The rate of proliferation is then linked to a clonal growth model to predict tumor incidences. Comparisons of the BBDR simulations and parameterizations across chemicals suggested that significant variation among the models for the three chemicals arises in a few parameters expected to be chemical specific (such as metabolism and cellular injury rate constants). Optimization of model parameters to tumor data for CCL(4) and DMF resulted in similar estimates for all parameters related to cytotoxicity and tumor incidences. However, optimization of the CHCl(3) data resulted in a higher estimate for one parameter (BD) related to death of initiated cells. This implies that additional steps beyond cytotoxicity leading to induced cellular proliferation can be quantitatively different among chemicals that share cytotoxicity as a hypothesized carcinogenic MOA.

Employing a Mechanistic Model for the MAPK Pathway to Examine the Impact of Cellular all or None Behavior on Overall Tissue Response.

The mitogen activated protein kinase (MAPK) cascade is a three-tiered phosphorylation cascade that is ubiquitously expressed among eukaryotic cells. Its primary function is to propagate signals from cell surface receptors to various cytosolic and nuclear targets. Recent studies have demonstrated that the MAPK cascade exhibits an all-or-none response to graded stimuli. This study quantitatively investigates MAPK activation in Xenopus oocytes using both empirical and biologically-based mechanistic models. Empirical models can represent overall tissue MAPK activation in the oocytes. However, these models lack description of key biological processes and therefore give no insight into whether the cellular response occurs in a graded or all-or-none fashion. To examine the propagation of cellular MAPK all-or-none activation to overall tissue response, mechanistic models in conjunction with Monte Carlo simulations are employed. An adequate description of the dose response relationship of MAPK activation in Xenopus oocytes is achieved. Furthermore, application of these mechanistic models revealed that the initial receptor-ligand binding rate contributes to the cells' ability to exhibit an all-or-none MAPK activation response, while downstream activation parameters contribute more to the magnitude of activation. These mechanistic models enable us to identify key biological events which quantitatively impact the shape of the dose response curve, especially at low environmentally relevant doses.