Global optimization of the Michaelis-Menten parameters using physiologically-based pharmacokinetic (PBPK) modeling and chloroform vapor uptake data in F344 rats.

Objective: To quantify metabolism, a physiologically based pharmacokinetic (PBPK) model for a volatile compound can be calibrated with the closed chamber (i.e. vapor uptake) inhalation data. Here, we introduce global optimization as a novel component of the predictive process and use it to illustrate a procedure for metabolic parameter estimation.Materials and methods: Male F344 rats were exposed in vapor uptake chambers to initial concentrations of 100, 500, 1000, and 3000 ppm chloroform. Chamber time-course data from these experiments, in combination with optimization using a chemical-specific PBPK model, were used to estimate Michaelis-Menten metabolic constants. Matlab® simulation software was used to integrate the mass balance equations and to perform the global optimizations using MEIGO (MEtaheuristics for systems biology and bIoinformatics Global Optimization - Version 64 bit, R2016A), a toolbox written for Matlab®. The cost function used the chamber time-course data and least squares to minimize the difference between data and simulation values.Results and discussion: The final values estimated for Vmax (maximum metabolic rate) and Km (affinity constant) were 1.2 mg/h and a range between 0.0005 and 0.6 mg/L, respectively. Also, cost function plots were used to analyze the dose-dependent capacity to estimate Vmax and Km within the experimental range used. Sensitivity analysis was used to assess identifiability for both parameters and show these kinetic data may not be sufficient to identify Km.Conclusion: In summary, this work should help toxicologists interested in optimization techniques understand the overall process employed when calibrating metabolic parameters in a PBPK model with inhalation data.

Development and application of a human PBPK model for bromodichloromethane to investigate the impacts of multi-route exposure.

As a result of its presence in water as a volatile disinfection byproduct, bromodichloromethane (BDCM), which is mutagenic, poses a potential health risk from exposure via oral, dermal and inhalation routes. We developed a refined human physiologically based pharmacokinetic (PBPK) model for BDCM (including new chemical-specific human parameters) to evaluate the impact of BDCM exposure during showering and bathing on important measures of internal dose compared with oral exposure. The refined model adequately predicted data from the published literature for oral, dermal and bathing/showering exposures. A liter equivalency approach (L-eq) was used to estimate BDCM concentration in a liter of water consumed by the oral route that would be required to produce the same internal dose of BDCM resulting from a 20-min bath or a 10-min shower in water containing 10 µg l(-1) BDCM. The oral liter equivalent concentrations for the bathing scenario were 605, 803 and 5 µg l(-1) BDCM for maximum venous blood concentration (Cmax), the area under the curve (AUCv) and the amount metabolized in the liver per hour (MBDCM), respectively. For a 10-min showering exposure, the oral L-eq concentrations were 282, 312 and 2.1 µg l(-1) for Cmax, AUC and MBDCM, respectively. These results demonstrate large contributions of dermal and inhalation exposure routes to the internal dose of parent chemical reaching the systemic circulation, which could be transformed to mutagenic metabolites in extrahepatic target tissues. Thus, consideration of the contribution of multiple routes of exposure when evaluating risks from water-borne BDCM is needed, and this refined human model will facilitate improved assessment of internal doses from real-world exposures. Published 2015. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

The impact of variation in scaling factors on the estimation of internal dose metrics: a case study using bromodichloromethane (BDCM).

A rate for hepatic metabolism (Vmax) determined in vitro must be scaled for in vivo use in a physiologically based pharmacokinetic (PBPK) model. This requires the use of scaling factors such as mg of microsomal protein per gram of liver (MPPGL) and liver mass (FVL). Variation in MPPGL and FVL impacts variation in Vmax, and hence PBPK model-derived estimates of internal dose used in dose response analysis. The impacts of adult human variation in MPPGL and FVL on estimates of internal dose were assessed using a human PBPK model for bromodichloromethane (BDCM), a water disinfection byproduct, for multiple internal dose metrics for two exposure scenarios (single 0.25 liter drink of water or 10 min shower) under plausible (5 µg/L) and high level (20 µg/L) water concentrations. For both concentrations, all internal dose metrics were changed less than 5% for the showering scenario (combined inhalation and dermal exposure). In contrast, a 27-fold variation in area under the curve (AUC) for BDCM in venous blood was observed at both oral exposure concentrations, whereas total amount of BDCM metabolized in liver was relatively unchanged. This analysis demonstrates that variability in the scaling factors used for in vitro to in vivo extrapolation (IVIVE) for metabolic rate parameters can have a significant route-dependent impact on estimates of internal dose under environmentally relevant exposure scenarios. This indicates the need to evaluate both uncertainty and variability for scaling factors used for IVIVE.

Comparison of in vivo derived and scaled in vitro metabolic rate constants for several volatile organic compounds (VOCs).

Metabolic rate parameters estimation using in vitro data is necessary due to numbers of chemicals for which data are needed, trend towards minimizing laboratory animal use, and limited opportunity to collect data in human subjects. We evaluated how well metabolic rate parameters derived from in vitro data predict overall in vivo metabolism for a set of environmental chemicals for which well validated and established methods exist. We compared values of VmaxC derived from in vivo vapor uptake studies with estimates of VmaxC scaled up from in vitro hepatic microsomal metabolism studies for VOCs for which data were available in male F344 rats. For 6 of 7 VOCs, differences between the in vivo and scaled up in vitro VmaxC estimates were less than 2.6-fold. For bromodichloromethane (BDCM), the in vivo derived VmaxC was approximately 4.4-fold higher than the in vitro derived and scaled up VmaxC. The more rapid rate of BDCM metabolism estimated based in vivo studies suggests other factors such as extrahepatic metabolism, binding or other non-specific losses making a significant contribution to overall clearance. Systematic and reliable utilization of scaled up in vitro biotransformation rate parameters in PBPK models will require development of methods to predict cases in which extrahepatic metabolism and binding as well as other factors are likely to be significant contributors.

Modeling the toxicokinetics of inhaled toluene in rats: influence of physical activity and feeding status.

Toluene is found in petroleum-based fuels and used as a solvent in consumer products and industrial applications. The critical effects following inhalation exposure involve the brain and nervous system in both humans and experimental animals, whether exposure duration is acute or chronic. The goals of this physiologically based pharmacokinetic (PBPK) model development effort were twofold: (1) to evaluate and explain the influence of feeding status and activity level on toluene pharmacokinetics utilizing our own data from toluene-exposed Long Evans (LE) rats, and (2) to evaluate the ability of the model to simulate data from the published literature and explain differing toluene kinetics. Compartments in the model were lung, slowly and rapidly perfused tissue groups, fat, liver, gut, and brain; tissue transport was blood-flow limited and metabolism occurred in the liver. Chemical-specific parameters and initial organ volumes and blood flow rates were obtained from the literature. Sensitivity analysis revealed that the single most influential parameter for our experimental conditions was alveolar ventilation; other moderately influential parameters (depending upon concentration) included cardiac output, rate of metabolism, and blood flow to fat. Based on both literature review and sensitivity analysis, other parameters (e.g., partition coefficients and metabolic rate parameters) were either well defined (multiple consistent experimental results with low variability) or relatively noninfluential (e.g. organ volumes). Rats that were weight-maintained compared to free-fed rats in our studies could be modeled with a single set of parameters because feeding status did not have a significant impact on toluene pharmacokinetics. Heart rate (HR) measurements in rats performing a lever-pressing task indicated that the HR increased in proportion to task intensity. For rats acclimated to eating in the lab during the day, both sedentary rats and rats performing the lever-pressing task required different alveolar ventilation rates to successfully predict the data. Model evaluation using data from diverse sources together with statistical evaluation of the resulting fits revealed that the model appropriately predicted blood and brain toluene concentrations with some minor exceptions. These results (1) emphasize the importance of experimental conditions and physiological status in explaining differing kinetic data, and (2) demonstrate the need to consider simulation conditions when estimating internal dose metrics for toxicity studies in which kinetic data were not collected.

A General Physiological and Toxicokinetic (GPAT) Model for Simulating Complex Toluene Exposure Scenarios in Humans.

Realistic simulation of environmental exposure scenarios requires dynamic methods in which exposures and human activities vary continuously as a function of time. Simulation of such complex scenarios is, with conventional physiologically based methods, a complex and programming-intensive task. The goal of the present effort was to simplify this task by combining a commercially available general whole-body human physiological model (QCP2004) with a slightly extended physiologically based toxicokinetic (PBTK) model from the literature. The QCP2004 model is a differential equation-based model similar to PBTK models except that normal organ function is simulated and the body organs are appropriately interlinked. Here QCP2004 provided estimates of physiological parameters required by the PBTK model. These were updated as the model was iteratively executed appropriate to the varying activity of the human subject. The combined general physiological model and the PBTK model was called a general physiological and toxicokinetic (GPAT) model. The GPAT model was tested and (within the constraints of available toluene exposure experiments in the literature) found to predict toluene blood concentrations, even in dynamic situations. A model of the structure used in the present work is capable of expansion as new knowledge is developed and greater detail is desired. Similarly, multiple toxicant PBTK models can be developed and incorporated for applications to mixtures risk assessment. Additionally, toxicant effects on organ systems can be achieved by altering organ function during a simulation as a function of the internal dose of toxicants. By cumulatively adding detail to the model as new physiological and chemical-specific information becomes available, the model can become a repository of knowledge for increasingly sophisticated risk-assessment applications.

Applications of dosimetry modeling to assessment of neurotoxic risk.

Risk assessment procedures can be improved through better understanding and use of tissue dose information and linking tissue dose level to adverse outcomes. For volatile organic compounds, such as toluene and trichloroethylene (TCE), blood and brain concentrations can be estimated with physiologically based pharmacokinetic (PBPK) models. Acute changes in the function of the nervous system can be linked to the concentration of test compounds in the blood or brain at the time of neurological assessment. This set of information enables application to a number of risk assessment situations. For example, we have used this approach to recommend duration adjustments for acute exposure guideline levels (AEGLs) for TCE such that the exposure limits for each exposure duration yield identical tissue concentrations at the end of the exposure period. We have also used information on tissue concentration at the time of assessment to compare sensitivity across species, adjusting for species-specific pharmacokinetic differences. Finally this approach has enabled us to compare the relative sensitivity of different compounds on a tissue dose basis, leading to expression of acute solvent effects as ethanol-dose equivalents for purposes of estimating cost-benefit relationships of various environmental control options.

Developing an exposure-dose-response model for the acute neurotoxicity of organic solvents: overview and progress on in vitro models and dosimetry.

We are developing an exposure-dose-response (EDR) model for volatile organic compounds (VOCs) to predict acute effects of VOCs on nervous system function from exposure data (concentration and duration of inhalation). This model contains both toxicokinetic and toxicodynamic components. One advantage of the EDR model will be its ability to relate in vitro effects of solvents on cellular ion channels (putative targets) to in vivo effects, using a combination of physiologically-based toxicokinetic (PBTK) modeling (to estimate VOC concentrations in the blood and brain) and in vitro studies to clarify the mode of action of the VOCs. Recent work in vitro has focused on quantifying the inhibitory effects of toluene, trichloroethylene (TCE) and perchloroethylene (PERC) on ion channel currents. All three VOCs inhibit current through voltage-sensitive calcium channels (VSCCs) in pheochromocytoma cells; PERC blocked calcium currents and altered the current-voltage relationship at lower concentrations than did toluene or TCE. Recombinant nicotinic acetylcholine receptors (nAChRs), expressed in Xenopus oocytes, were also inhibited by PERC and toluene in a concentration-dependent manner. PERC inhibited α7 receptors more than α4ß2 receptors in recombinant human and rat nAChRs. However, human and rat α7 receptors were equally sensitive to PERC and TOL. These in vitro studies will be used to identify an appropriate neuronal receptor system to serve as an index of acute effects of VOCs in vivo. The PBTK model incorporates physiological input parameters derived from radiotelemetered heart rate data from rats performing operant tests of cognitive and motor functions. These studies should improve predictions of target organ concentrations of inhaled VOCs in subjects actively performing behavioral tests over a range of physical activity levels.

Duration adjustment of acute exposure guideline level values for trichloroethylene using a physiologically-based pharmacokinetic model.

Acute Exposure Guideline Level (AEGL) recommendations are developed for 10-minute, 30-minute, 1-hour, 4-hours, and 8-hours exposure durations and are designated for three levels of severity: AEGL-1 represents concentrations above which acute exposures may cause noticeable discomfort including irritation; AEGL-2 represents concentrations above which acute exposure may cause irreversible health effects or impaired ability to escape; and AEGL-3 represents concentrations above which exposure may cause life-threatening health effects or death. The default procedure for setting AEGL values across durations when applicable data are unavailable involves estimation based on Haber's rule, which has an underlying assumption that cumulative exposure is the determinant of toxicity. For acute exposure to trichloroethylene (TCE), however, experimental data indicate that momentary tissue concentration, and not the cumulative amount of exposure, is important. We employed an alternative approach to duration adjustments in which a physiologically-based pharmacokinetic (PBPK) model was used to predict the arterial blood concentrations [TCE(a)] associated with adverse outcomes appropriate for AEGL-1, -2, or -3-level effects. The PBPK model was then used to estimate the atmospheric concentration that produces equivalent [TCE(a)] at each of the AEGL-specific exposure durations. This approach yielded [TCE(a)] values of 4.89 mg/l for AEGL-1, 18.7 mg/l for AEGL-2, and 310 mg/l for AEGL-3. Duration adjustments based on equivalent target tissue doses should provide similar degrees of toxicity protection at different exposure durations.