Physiologically based pharmacokinetic modelling of human exposure to 2-butoxyethanol

Development, Testing, Parameterization, and Calibration of a Human Physiologically Based Pharmacokinetic Model for the Plasticizer, Hexamoll® Diisononyl-Cyclohexane-1, 2-Dicarboxylate Using In Silico, In Vitro, and Human Biomonitoring Data

A physiologically based pharmacokinetic model for Hexamoll^® diisononyl-cyclohexane-1, 2-dicarboxylate was developed to interpret the biokinetics in humans after single oral doses. The model was parameterized with in vitro and in silico derived parameters and uncertainty and sensitivity analysis was used during the model development process to assess structure, biological plausibility and behavior prior to simulation and analysis of human biological monitoring data. The model provided good simulations of the urinary excretion (Curine) of two metabolites; cyclohexane-1,2-dicarboxylic acid mono hydroxyisononyl ester (OH-MINCH) and cyclohexane-1, 2-dicarboxylic acid mono carboxyisononyl ester (cx-MINCH) from the biotransformation of mono-isononyl-cyclohexane-1, 2-dicarboxylate (MINCH), the monoester metabolite of di-isononyl-cyclohexane-1,2-dicarboxylate. However, good simulations could be obtained, with and without, a lymphatic compartment. Selection of an appropriate model structure was informed by sensitivity analysis which could identify and quantify the contribution to variability in Curine by parameters, such as, the fraction of oral dose that directly entered the lymphatic compartment and therefore by-passed the liver and the fraction of MINCH bio-transformed to cx-MINCH and OH-MINCH. By constraining these parameters within biologically plausible limits the presence of a lymphatic compartment was deemed an important component of model structure. Furthermore, the use of sensitivity analysis is important in the evaluation of uncertainty around in silico derived parameters. By quantifying their impact on model output sufficient confidence in the use of a model should be afforded. This type of approach could expand the use of physiologically based pharmacokinetic models since parameterization with in silico techniques allows for rapid model development. This in turn could assist in reducing the use of animals in toxicological evaluations by enhancing the utility of “read across” techniques.

The application of global sensitivity analysis in the development of a physiologically based pharmacokinetic model for m-xylene and ethanol co-exposure in humans

Global sensitivity analysis (SA) was used during the development phase of a binary chemical physiologically based pharmacokinetic (PBPK) model used for the analysis of m-xylene and ethanol co-exposure in humans. SA was used to identify those parameters which had the most significant impact on variability of venous blood and exhaled m-xylene and urinary excretion of the major metabolite of m-xylene metabolism, 3-methyl hippuric acid. This analysis informed the selection of parameters for estimation/calibration by fitting to measured biological monitoring (BM) data in a Bayesian framework using Markov chain Monte Carlo (MCMC) simulation. Data generated in controlled human studies were shown to be useful for investigating the structure and quantitative outputs of PBPK models as well as the biological plausibility and variability of parameters for which measured values were not available. This approach ensured that a priori knowledge in the form of prior distributions was ascribed only to those parameters that were identified as having the greatest impact on variability. This is an efficient approach which helps reduce computational cost.

Reprint of PopGen: A virtual human population generator

The risk assessment of environmental chemicals and drugs is moving towards a paradigm shift in approach which seeks the full replacement animal testing with high throughput, mechanistic, in vitro systems. This new vision will be reliant on the measurement in vitro, of concentration-dependent responses where prolonged excessive perturbations of specific biochemical pathways are likely to lead to adverse health effects in an intact organism. Such an approach requires a framework, into which disparate data generated using in vitro, in silico and in chemico systems, can be integrated and utilised for quantitative in vitro-to-in vivo extrapolation (QIVIVE), ultimately to the human population level. Physiologically based pharmacokinetic (PBPK) models are ideally suited for this and are obligatory in order to translate in vitro concentration-response relationships to an exposure or dose, route and duration regime in people. In this report we describe PopGen a virtual human population generator which is a user friendly, open access web-based application for the prediction of realistic anatomical, physiological and phase 1 metabolic variation in a wide range of healthy human populations. We demonstrate how PopGen can be used for QIVIVE by providing input to a PBPK model, at an appropriate level of detail, to reconstruct exposure from human biomonitoring data. We discuss how the process of exposure reconstruction from blood biomarkers, in general, is analogous to exposure or dose reconstruction from concentration-response measurements made in proposed in vitro cell based systems which are assumed to be surrogates for target organs.

PopGen: A virtual human population generator

The risk assessment of environmental chemicals and drugs is moving towards a paradigm shift in approach which seeks the full replacement animal testing with high throughput, mechanistic, in vitro systems. This new vision will be reliant on the measurement in vitro, of concentration-dependent responses where prolonged excessive perturbations of specific biochemical pathways are likely to lead to adverse health effects in an intact organism. Such an approach requires a framework, into which disparate data generated using in vitro, in silico and in chemico systems, can be integrated and utilised for quantitative in vitro-to-in vivo extrapolation (QIVIVE), ultimately to the human population level. Physiologically based pharmacokinetic (PBPK) models are ideally suited for this and are obligatory in order to translate in vitro concentration-response relationships to an exposure or dose, route and duration regime in people. In this report we describe PopGen, a virtual human population generator which is a user friendly, open access web-based application for the prediction of realistic anatomical, physiological and phase 1 metabolic variation in a wide range of healthy human populations. We demonstrate how PopGen can be used for QIVIVE by providing input to a PBPK model, at an appropriate level of detail, to reconstruct exposure from human biomonitoring data. We discuss how the process of exposure reconstruction from blood biomarkers, in general, is analogous to exposure or dose reconstruction from concentration-response measurements made in proposed in vitro cell based systems which are assumed to be surrogates for target organs.

Reconstruction of Exposure to m-Xylene from Human Biomonitoring Data Using PBPK Modelling, Bayesian Inference, and Markov Chain Monte Carlo Simulation

There are numerous biomonitoring programs, both recent and ongoing, to evaluate environmental exposure of humans to chemicals. Due to the lack of exposure and kinetic data, the correlation of biomarker levels with exposure concentrations leads to difficulty in utilizing biomonitoring data for biological guidance values. Exposure reconstruction or reverse dosimetry is the retrospective interpretation of external exposure consistent with biomonitoring data. We investigated the integration of physiologically based pharmacokinetic modelling, global sensitivity analysis, Bayesian inference, and Markov chain Monte Carlo simulation to obtain a population estimate of inhalation exposure to m-xylene. We used exhaled breath and venous blood m-xylene and urinary 3-methylhippuric acid measurements from a controlled human volunteer study in order to evaluate the ability of our computational framework to predict known inhalation exposures. We also investigated the importance of model structure and dimensionality with respect to its ability to reconstruct exposure.

A generic, cross-chemical predictive PBTK model with multiple entry routes running as application in MS Excel; design of the model and comparison of predictions with experimental results

AIM

Physiologically based toxicokinetic (PBTK) models are computational tools, which simulate the absorption, distribution, metabolism, and excretion of chemicals. The purpose of this study was to develop a physiologically based pharmacokinetic (PBPK) model with a high level of transparency. The model should be able to predict blood and urine concentrations of environmental chemicals and metabolites, given a certain environmental or occupational exposure scenario.

MODEL

The model refers to a reference human of 70 kg. The partition coefficients of the parent compound and its metabolites (blood:air and tissue:blood partition coefficients of 11 organs) are estimated by means of quantitative structure-property relationship, in which five easily available physicochemical properties of the compound are the independent parameters. The model gives a prediction of the fate of the compound, based on easily available chemical properties; therefore, it can be applied as a generic model applicable to multiple compounds. Three routes of uptake are considered (inhalation, dermal, and/or oral) as well as two built-in exercise levels (at rest and at light work). Dermal uptake is estimated by the use of a dermal diffusion-based module that considers dermal deposition rate and duration of deposition. Moreover, evaporation during skin contact is fully accounted for and related to the volatility of the substance. Saturable metabolism according to Michaelis-Menten kinetics can be modelled in any of 11 organs/tissues or in liver only. Renal tubular resorption is based on a built-in algorithm, dependent on the (log) octanol:water partition coefficient. Enterohepatic circulation is optional at a user-defined rate. The generic PBTK model is available as a spreadsheet application in MS Excel. The differential equations of the model are programmed in Visual Basic. Output is presented as numerical listing over time in tabular form and in graphs. The MS Excel application of the PBTK model is available as freeware.

EXPERIMENTAL

The accuracy of the model prediction is illustrated by simulating experimental observations. Published experimental inhalation and dermal exposure studies on a series of different chemicals (pyrene, N-methyl-pyrrolidone, methyl-tert-butylether, heptane, 2-butoxyethanol, and ethanol) were selected to compare the observed data with the model-simulated data. The examples show that the model-predicted concentrations in blood and/or urine after inhalation and/or transdermal uptake have an accuracy of within an order of magnitude.

CONCLUSIONS

It is advocated that this PBTK model, called IndusChemFate, is suitable for 'first tier assessments' and for early explorations of the fate of chemicals and/or metabolites in the human body. The availability of a simple model with a minimum burden of input information on the parent compound and its metabolites might be a stimulation to apply PBTK modelling more often in the field of biomonitoring and exposure science.

In vivo monitoring of epidermal absorption of hazardous substances by confocal Raman micro-spectroscopy

BACKGROUND

Presently, percutaneous absorption of potentially hazardous chemicals in humans can only be assessed in animal experiments, in vitro, or predicted mathematically. Our aim was to demonstrate the proof-of-principle of a novel quantitative in vivo assay for percutaneous absorption: confocal Raman micro-spectroscopy (CRS). The advantages and limitations of CRS for health risk assessments are discussed.

PATIENTS AND METHODS

2-butoxyethanol, toluene, and pyrene were applied in pure form, diluted in water, or in ethanol on the skin of three healthy volunteers. CRS measurements were done following application for 15 min and 3 hours. The concentrations of the three substances as a function of distance to the skin surface were calculated and further analyzed with regard to mass transport into the stratum corneum (μg/cm(2)) and the flux through the stratum corneum (μg/cm(2)h). The results were compared with the available data from literature.

RESULTS

Considering the preliminary nature of these data, good accordance with data from the literature was observed. In addition, we observed that 2-butoxyethanol penetrates markedly faster when dissolved in water as compared to ethanol. This observation is also in agreement with previous results.

CONCLUSIONS

CRS has the potential to provide fast, accurate and reliable results for advanced studies of in vivo percutaneous absorption kinetics of hazardous substances in human skin. This will require further research with other substances and under differing conditions.

A Workflow for Global Sensitivity Analysis of PBPK Models

Physiologically based pharmacokinetic (PBPK) models have a potentially significant role in the development of a reliable predictive toxicity testing strategy. The structure of PBPK models are ideal frameworks into which disparate in vitro and in vivo data can be integrated and utilized to translate information generated, using alternative to animal measures of toxicity and human biological monitoring data, into plausible corresponding exposures. However, these models invariably include the description of well known non-linear biological processes such as, enzyme saturation and interactions between parameters such as, organ mass and body mass. Therefore, an appropriate sensitivity analysis (SA) technique is required which can quantify the influences associated with individual parameters, interactions between parameters and any non-linear processes. In this report we have defined the elements of a workflow for SA of PBPK models that is computationally feasible, accounts for interactions between parameters, and can be displayed in the form of a bar chart and cumulative sum line (Lowry plot), which we believe is intuitive and appropriate for toxicologists, risk assessors, and regulators.

Pharmacokinetic predictions in children by using the physiologically based pharmacokinetic modelling

Nowadays, 50-90% of drugs used in children have never been actually studied in this population. Consequently, either our children are often exposed to the risk of adverse drug events or to lack of efficacy, or they are unable to benefit from a number of therapeutic advances offered to adults, as no clinical study has been properly performed in children. Actually the main methods used to calculate the dose for a child are based on allometric methods taking into account different categories of age, the body weight and/or the body surface area. Unfortunately, these calculation methods consider the children as small adults, which is not the case. Physiologically based pharmacokinetics is one way to integrate the physiological changes occurring in the childhood and to anticipate their impact on the pharmacokinetic processes: absorption, distribution, metabolism and excretion/elimination. From different examples, the application of this modelling approach is discussed as a possible and valuable method to minimize the ethical and technical difficulties of conducting research in children.

Biomonitoring as a tool in the human health risk characterization of dermal exposure

Dermal exposure is an important factor in risk characterization. In occupational settings it becomes relatively more important because of the continuous reduction in inhalation exposure. In the public health arena, dermal exposure may also form a significant contribution to the total exposure. Dermal exposure, however, is difficult to assess directly because it is determined by a host of factors, which are difficult to quantify. As a consequence, dermal exposure is often estimated by application of models for external exposure. In combination with modeled or measured data for percutaneous penetration, these provide an estimate for the internal exposure that is directly related to the systemic effects. The advantages and drawbacks of EASE (Estimation and Assessment of Substance Exposure) and RISKOFDERM (Risk Assessment of Occupational Dermal Exposure), two models for external exposure that are mentioned in the Technical Guidance Document for the European Union risk assessments performed under the Existing Substances Regulation (EEC/793/93), are discussed. Although new chemicals regulation (REACh, 1907/2006/EC) is now in place in the European Union, the principles applied under the previous legislation do not change and the same models will continue to be used. The results obtained with these models for styrene, 2-butoxyethanol, and 1-methoxy-2-propanol in specific exposure scenarios are compared with an alternative method that uses biomonitoring data to assess dermal exposure. Actual external exposure measurements combined with measured or modeled percutaneous penetration data give acceptable results in risk assessment of dermal exposure, but modeled data of external dermal exposure should only be used if no other data are available. However, if available, biomonitoring should be considered the method of choice to assess (dermal) exposure.

Projection of Exposure and Efficacious Dose Prior to First-in-Human Studies: How Successful Have We Been?

PURPOSE

Preclinical and clinical data for 35 proprietary Bristol-Myers Squibb discovery compounds (years 1997 to 2005) were collected and analyzed. In each case, exposure and efficacy in human subjects were projected at the time of nomination (for development) prior to first-in-human dosing.

MATERIALS AND METHODS

Projections of area under the plasma concentration-time curve (AUC) in humans involved the use of one or more methods: (1) allometric scaling of animal pharmacokinetic data; (2) clearance projection employing in vitro data (liver microsomes and hepatocytes); (3) chimpanzee as an animal model; (4) the species-invariant time method; and (5) the Css-mean residence time or "Css-MRT" method. Whenever possible, prior clinical experience with lead compounds enabled the selection of the most appropriate method(s). Multiple approaches were also available at the time of the human efficacious dose projections: (1) efficacious exposure from animal efficacy models; (2) in vitro potency; and (3) prior experience with clinical leads.

RESULTS

Over the 8 year period described, AUC in humans was projected within 2-fold (20 out of 35 compounds; 57%), greater than 2-fold to 4-fold (11 out of 35 compounds; 32%), and greater than 4-fold (4 out of 35 compounds; 11%) of the observed value. At the time of writing, clinical efficacy data were available for 10 compounds only. In this instance, the efficacious doses were also projected within 2-fold (7 out of 10 compounds; 70%), greater than 2-fold to 4-fold (2 out of 10 compounds; 20%), and greater than 4-fold (1 out of 10 compounds; 10%) of the actual clinical dose.

CONCLUSION

Overall, it was possible to project human exposure and efficacious dose within 4-fold of observed clinical values for about 90% of the compounds.