1. Global optimization of the Michaelis–Menten parameters using physiologically-based pharmacokinetic (PBPK) modeling and chloroform vapor uptake data in F344 rats
- Author
-
Jane Ellen Simmons, Yusupha M Sey, David N Williams, Marina V. Evans, and Christopher R. Eklund
- Subjects
Male ,Physiologically based pharmacokinetic modelling ,Health, Toxicology and Mutagenesis ,F344 rats ,010501 environmental sciences ,Closed chamber ,Kidney ,Toxicology ,Models, Biological ,01 natural sciences ,Michaelis–Menten kinetics ,Article ,03 medical and health sciences ,chemistry.chemical_compound ,0302 clinical medicine ,Pharmacokinetics ,Administration, Inhalation ,Animals ,Computer Simulation ,0105 earth and related environmental sciences ,Chromatography ,Chloroform ,Inhalation ,Chemistry ,Muscles ,fungi ,food and beverages ,Rats, Inbred F344 ,Adipose Tissue ,Liver ,030228 respiratory system - Abstract
PBPK models are well-established frameworks used to describe absorption, distribution, metabolism, and excretion of xenobiotics. To quantify metabolism, a PBPK model for a volatile compound can be calibrated with 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. 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 64bit, 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. The final values estimated for V(max) (maximum metabolic rate) and K(m) (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 V(max) and K(m) 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 K(m). The importance of an accurate estimation of metabolism using vapor uptake data is discussed. 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.
- Published
- 2020
- Full Text
- View/download PDF