Systematic review of physiologically based kinetic lactation models for transfer of xenobiotic compounds to milk.

Lactational elimination has been described mathematically for nearly 50 years. Over 40 published articles, containing >50 physiologically based kinetic (PBK) lactation models were included in the systematic review. These PBK models described the lactational elimination of xenobiotic compounds in humans, rats, mice, and dairy cows and goats. A total of 78 compounds have been modelled, ranging from industrial chemicals, pesticides, to pain medication, antibiotics, and caffeine. Few models included several species or compounds, and models were thus generally not translational or generic. Three dairy cow models mechanistically described the intramammary disposition of pharmaceuticals after intramammary administration, including volume changes caused by milking, while empirically describing the remaining pharmacokinetics. The remaining models were semi- or whole body PBK models, describing long-term exposure of environmental pollutants, or short-term exposure of pharmaceuticals. The absolute majority described the disposition to the mammary gland or milk with perfusion limited compartments, but permeability limited models were available as well. With long-term exposure, models often included changes in milk volume and/or consumption by the offspring, and changes in body weight of offspring. Periodic emptying of the mammary gland, as with feeding or milking, was sparsely applied. Rodent models used similar physiological parameters, while values of physiological parameters applied in human models could range widely. When milk composition was included in the models, it most often included the fat content. The review gives an extensive overview of the applied functions and modelling strategies of PBK lactation models.

Addressing systemic problems with exposure assessments to protect the public's health.

BACKGROUND: Understanding, characterizing, and quantifying human exposures to environmental chemicals is critical to protect public health. Exposure assessments are key to determining risks to the general population and for specific subpopulations given that exposures differ between groups. Exposure data are also important for understanding where interventions, including public policies, should be targeted and the extent to which interventions have been successful. In this review, we aim to show how inadequacies in exposure assessments conducted by polluting industries or regulatory agencies have led to downplaying or disregarding exposure concerns raised by communities; that underestimates of exposure can lead regulatory agencies to conclude that unacceptable risks are, instead, acceptable, allowing pollutants to go unregulated; and that researchers, risk assessors, and policy makers need to better understand the issues that have affected exposure assessments and how appropriate use of exposure data can contribute to health-protective decisions. METHODS: We describe current approaches used by regulatory agencies to estimate human exposures to environmental chemicals, including approaches to address limitations in exposure data. We then illustrate how some exposure assessments have been used to reach flawed conclusions about environmental chemicals and make recommendations for improvements. RESULTS: Exposure data are important for communities, public health advocates, scientists, policy makers, and other groups to understand the extent of environmental exposures in diverse populations. We identify four areas where exposure assessments need to be improved due to systemic sources of error or uncertainty in exposure assessments and illustrate these areas with examples. These include: (1) an inability of regulatory agencies to keep pace with the increasing number of chemicals registered for use or assess their exposures, as well as complications added by use of 'confidential business information' which reduce available exposure data; (2) the failure to keep assessments up-to-date; (3) how inadequate assumptions about human behaviors and co-exposures contribute to underestimates of exposure; and (4) that insufficient models of toxicokinetics similarly affect exposure estimates. CONCLUSION: We identified key issues that impact capacity to conduct scientifically robust exposure assessments. These issues must be addressed with scientific or policy approaches to improve estimates of exposure and protect public health.

A Novel Multispecies Toxicokinetic Modeling Approach in Support of Chemical Risk Assessment.

Standardized laboratory tests with a limited number of model species are a key component of chemical risk assessments. These surrogate species cannot represent the entire diversity of native species, but there are practical and ethical objections against testing chemicals in a large variety of species. In previous research, we have developed a multispecies toxicokinetic model to extrapolate chemical bioconcentration across species by combining single-species physiologically based toxicokinetic (PBTK) models. This "top-down" approach was limited, however, by the availability of fully parameterized single-species models. Here, we present a "bottom-up" multispecies PBTK model based on available data from 69 freshwater fishes found in Canada. Monte Carlo-like simulations were performed using statistical distributions of model parameters derived from these data to predict steady-state bioconcentration factors (BCFs) for a set of well-studied chemicals. The distributions of predicted BCFs for 1,4-dichlorobenzene and dichlorodiphenyltrichloroethane largely overlapped those of empirical data, although a tendency existed toward overestimation of measured values. When expressed as means, predicted BCFs for 26 of 34 chemicals (82%) deviated by less than 10-fold from measured data, indicating an accuracy similar to that of previously published single-species models. This new model potentially enables more environmentally relevant predictions of bioconcentration in support of chemical risk assessments.

Tissue-Specific Accumulation, Biotransformation, and Physiologically Based Toxicokinetic Modeling of Benzotriazole Ultraviolet Stabilizers in Zebrafish (Danio rerio).

Benzotriazole ultraviolet stabilizers (BUVSs) are high-production-volume chemicals with ubiquitous occurrence in the aquatic environment. However, little is known about their bioconcentration and biotransformation, and physiologically based toxicokinetic (PBTK) models for BUVSs are lacking. This study selected six BUVSs for which experiments were performed with zebrafish (Danio rerio) exposed to two different levels (0.5 and 10 µg·L-1). Higher kinetic bioconcentration factors (BCFs) were observed at the lower exposure level with environmental relevance, with BCF of 3.33 × 103 L·kg-1 for 2-(2-hydroxy-3,5-di-tert-butylphenyl)-5-chlorobenzotriazole (UV-327). This phenomenon was interpreted by a nonlinear adsorption mechanism, where binding with specific protein sites contributes to bioconcentration. Muscle exhibited the lowest accumulation, in which depuration half-life of UV-327 was 19.5 d. In kidney, muscle, ovary, gill, and skin, logBCF increased with increase in log KOW of the BUVSs until log KOW was ca. 6.5, above which logBCF decreased. However, the trend was not observed in the liver and intestine. Six biotransformation products were identified and mainly accumulated in the liver and intestine. Considering the nonlinear adsorption mechanism in the PBTK model, the prediction accuracy of the model was improved, highlighting the binding of xenobiotics with specific protein sites in assessing the bioconcentration of chemicals for their risk assessment.

Evaluation of toxicokinetics of nonylphenol in the adult female Sprague-Dawley rats using a physiologically based toxicokinetic model.

The physiologically based toxicokinetic (PBTK) model was firstly developed in female rats to quantitatively evaluate toxicokinetics of nonylphenol (NP). Changes in NP serum concentrations over time of single oral NP administration experiments in Sprague-Dawley rats and literature data were collected to establish and calibrate the PBTK model in the SimBiology framework. The calibrated model predicted the serum and tissue NP concentrations of repeat oral NP administration for model evaluation. NP concentrations in serum and tissues (liver, kidneys, adipose, brain, uterus and ovaries) were quantified using ultra-high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The model output of the time course data (values are within the standard deviation defined for each data point) indicated proximity of predictions to reality. The coefficients of determination r2 were all greater than 0.9, and the root mean squared error (RMSE) were within 0.177-2.027, which means the model predicted and observed serum NP concentrations were in excellent agreement. The results indicated that the model could contribute to a simplification of the future exposure risk assessments of NP in a more realistic scenario and provided a better understanding of the disposition process of NP in human.

[Physiologically based toxicokinetic model for nickel].

OBJECTIVE: To estimate the metabolic parameters in different tissues and organs, build the physiologically based pharmacokinetic( PBPK) model of rat and occupational population, and predict the toxic dynamic characteristics exposure to nickel. METHODS: The partition coefficients in different tissues and organs were estimated using vector datas of nickel by the optimization and statistics files of acslx software. The PBTK model of occupational population exposure to nickel was built according to the metabolic parameters by acslx software. RESULTS: The evaluated partition coefficient of nickel were kidney blood( 0. 668), lung blood( 0. 102), spleen blood( 0. 037), liver blood( 0. 028), heart blood( 0. 022), and brain blood( 0. 006). The constructed successful PBPK model of occupational population exposed to 0. 1 mg/m~3 nickel for 8 hours showed that the nickel concentration is higher in kidney reached at 3. 328 µg/kg, followed by the spleen( 0. 185 µg/kg), liver( 0. 140 µg/kg) and heart( 0. 110 µg/kg). The content of nickel is lower in the brain( 0. 030 µg/kg). The kidneys is the major metabolic organs for nickel. CONCLUSION: The PBPK model can be used to convert the nickel levels from external exposure to internal exposure for each organ and to evaluate the time-dose relationship exposure to nickel in both rat and occupational population studies.

A physiologically based toxicokinetic model of P-glycoprotein transporter-mediated placenta perfusion of dexamethasone in the pregnant rat.

The present dosage of Dexamethasone (DEX) administered to pregnant women may pose a risk of toxicity to their unborn offspring. We aimed to develop a maternal-fetal physiologically based toxicokinetic (PBTK) model for DEX in pregnant rats, with a specific focus on the role of the P-glycoprotein (P-gp) transporter in placenta perfusion, and finally facilitate the optimization of clinical DEX dosage. We conducted animal experiments to determine DEX concentrations in various rat tissues, and constructed the PBTK model using MATLAB software. Sensitivity analysis was performed to assess input parameters and the model stability, with fold error (FE) values serving as evaluation indices. Our results indicate the successful construction of the PBTK model, with the fitting key parameters such as the absorption rate constant (Ka), intrinsic hepatic clearance (CLh,int) and intrinsic P-gp clearance (CLint,P-gp). The median concentration of DEX in maternal plasma, fetal plasma, fetal lung, and fetal brain were determined, which allowed us to fit the tissue-to-plasma partition coefficients for the fetal lung (Kp,lung,f) and fetal brain (Kp,brain,f). After making adjustments, all calculated FE values were found to be less than 2, demonstrating the acceptability and accuracy of our model's predictions. Our model integrated external literature data and internal animal experimentation to comprehensively evaluate the maternal-fetal PK characteristics of DEX. These findings provide valuable support for the optimization of clinical DEX dosing.

Bayesian analysis of physiologically based toxicokinetic (PBTK) modeling for pentachlorophenol exposure in pregnant women.

Pentachlorophenol (PCP) is a persistent organic compound that is widely present in the environment. The estimation of internal exposure levels for a given external exposure using toxicokinetic models is key to the human health risk assessment of PCP. The present study developed a physiologically based multicompartmental pharmacokinetic (PBTK) model to describe and predict the behavior of pentachlorophenol (PCP) in an organism. The model consists of stomach, intestines, adipose tissue, kidneys and fast- and poorly perfused tissues that are interconnected via blood circulation. We constructed a PBTK model of PCP in rats and extrapolated it to human dietary PCP exposure. The toxicokinetic data of PCP in human tissues and excreta were obtained from the published literature. Based on the collected PCP dietary survey and internal exposure data of pregnant women in Shanghai, Bayesian statistical analysis was performed for the model using Markov chain Monte Carlo (MCMC) simulation. The posterior distributions of the sensitive parameters were estimated, and the model was parameter optimized and validated using the pregnant women's test dataset. The results showed that the root mean square error (RMSE) improved 37.3% compared to the original model, and a systematic literature search revealed that the optimized model achieved acceptable prediction results for other datasets in China. A PCP metabolism model based on the exposure characteristics of pregnant women in China was constructed in the present study. The model provides a theoretical basis for the study of PCP toxicity and risk assessment.

Human risk assessment through development and application of a physiologically based toxicokinetic model for 4-tert-octylphenol.

4-tert-octylphenol (4-tert-OP) is an ecologically hazardous substance, and exposure to it in the environment has been consistently reported in the past. Despite the hazards and widespread exposure to 4-tert-OP, tools for scientific assessment of 4-tert-OP exposure risk level in humans are lacking. The main purpose of this study was to develop a physiologically-based-toxicokinetic (PBTK) model for 4-tert-OP and to perform quantitative risk assessment of 4-tert-OP in various population groups using the established model. Based on the results of toxicokinetic experiments on male rats, the PBTK model for 4-tert-OP was established and verified, and this was converted to a model for humans through interspecies extrapolation. Based on the previously reported no-observed-adverse-effect-levels for rats, it was possible to estimate the 4-tert-OP reference dose in humans through reverse dosimetry using the model. Biomonitoring data derived from various population groups were applied to the human PBTK model to calculate external exposures and margin of safety for 4-tert-OP for each population group. The PBTK model established in this study adequately explained the toxicokinetic experimental values at acceptable levels and was able to quantitatively predict the 4-tert-OP exposure level in the testes related to male reproductive toxicity. In addition, the degree of external exposure to 4-tert-OP could be scientifically estimated based on biomonitoring values derived from various biological matrices. The reference doses for systemic and reproductive toxicity caused by 4-tert-OP in male humans were calculated to be 0.16 and 1.12 mg/kg/day, respectively. The mean external exposure to 4-tert-OP in each population group estimated based on plasma and urine biomonitoring data was 0.04-66.24 mg/kg/day, showing very large exposure diversity between groups. Exposure risks to 4-tert-OP in populations ranged from safe to risky, suggesting the need for continued monitoring and risk management of 4-tert-OP worldwide. This study provides valuable scientific insight regarding the 4-tert-OP human risk assessment.

Recommendations for the reference concentration of cadmium exposure based on a physiologically based toxicokinetic model integrated with a human respiratory tract model.

Cadmium (Cd) poses a significant threat to human health. However, chronic toxicity parameters for inhalation exposure are lacking, especially for noncritical systemic toxic effects. A physiologically based toxicokinetic (PBTK) model can be used to extrapolate toxicity parameters across various exposure routes. We combined a PBTK model with a human respiratory tract (HRT) model, which is applicable to the general population and capable of simulating the deposition and clearance processes of various airborne Cd compounds in the respiratory tract. Monte Carlo analysis was used to simulate the distribution of sensitive parameters to reflect individual variability. Validation based on datasets from general and occupational populations showed that the improved model had acceptable or better predictive performance, outperforming the original model with a 14.45 % decrease in the root mean square error (RMSE). Using this PBTK-HRT model, we extrapolated toxicity parameters from oral exposure to inhalation exposure for four systemic toxic effects with doseresponse relationships but no known inhalation toxicity parameters, and ultimately recommended reference concentrations (RfCs) for four diseases (chronic kidney disease: 0.01 µg/m3, osteoporosis: 0.01 µg/m3, stroke: 0.04 µg/m3, diabetes mellitus: 0.13 µg/m3), contributing to a comprehensive assessment of the health risks of Cd inhalation exposure. ENVIRONMENTAL IMPLICATION: Cadmium (Cd), a heavy metal, can cause lung cancer, chronic kidney disease, and osteoporosis and pose a significant threat to human health. We combined a physiologically based toxicokinetic (PBTK) model with a human respiratory tract (HRT) model to achieve better predictive performance and wider applicability; this model was subsequently employed for route-to-route extrapolation of toxicity parameters. Additionally, for the first time, we focused on multiple subchronic and chronic systemic toxic effects in addition to critical effects and derived their reference concentrations (RfCs), which can be used to assess the health risk of Cd inhalation exposure more comprehensively and accurately.

Health risk assessment of cadmium exposure by integration of an in silico physiologically based toxicokinetic model and in vitro tests.

Cadmium (Cd) is a common environmental pollutant that can damage multiple organs, including the kidney. To prevent renal effects, international authorities have set health-based guidance values of Cd from epidemiological studies. To explore the health risk of Cd exposure and whether human equivalent doses (HEDs) derived from in vitro tests match the current guidance values, we integrated renal tubular epithelial cell-based assays with a physiologically based toxicokinetic model combined with the Monte Carlo method. For females, the HEDs (µg/kg/week) derived from KE2 (DNA damage), KE3 (cell cycle arrest), and KE4 (apoptosis) were 0.20 (2.5th-97.5th percentiles: 0.09-0.48), 0.52 (0.24-1.26), and 2.73 (1.27-6.57), respectively; for males the respective HEDs were 0.23 (0.10-0.49), 0.60 (0.27-1.30), and 3.11 (1.39-6.78). Among them, HEDKE4 (female) was close to the tolerable weekly intake (2.5 µg/kg/week) set by the European Food Safety Authority. The margin of exposure (MOE) derived from HEDKE4 (female) indicated that risks of renal toxicity for populations living in cadmium-contaminated regions should be of concern. This study provided a new approach methodology (NAM) for environmental chemical risk assessment using in silico and in vitro methods.

Physiologically-based toxicokinetic model for the prediction of perchlorate distribution and its application.

Perchlorate is a stable and readily transportable thyroid hormone disruptor, and prevalent exposure to perchlorate through food and drinking water has raised public concern about its health effects. The physiologically based toxicokinetic (PBTK) model as a dose prediction method is effective to predict the toxicant exposure dose of an organism and helps quantitatively assess the dose-dependent relationship with toxic effects. The current study aimed to establish a multi-compartment PBTK model based on updated time-course datasets of single oral exposure to perchlorate in rats. With adjustment of the kinetic parameters, the model fitted well the toxicokinetic characteristics of perchlorate in urine, blood, and thyroid from our experiments and the literature, and the coefficient of determination (R2) between the fitting values and the experimental data in regression analysis was greater than 0.91, indicating the robustness of the current model. The results of sensitivity analysis and daily repeated exposure simulations together confirmed its effective renal clearance. According to the distribution characteristic of perchlorate, a correlation study of internal and external exposure was conducted using urinary perchlorate as a bioassay indicator. The developed multi-compartment model for perchlorate updates important toxicokinetic data and kinetic parameters, providing analytical and modeling tools for deriving total exposure levels in the short term.

Inter-individual exposure variability interpretation through reflection of biological age algorithm in physiologically based toxicokinetic model: Application to human risk assessment of di-isobutyl-phthalate.

Age-related changes and interindividual variability in the degree of exposure to hazardous substances in the environment are pertinent factors to be considered in human risk assessment. Existing risk assessments remain in a one-size-fits-all approach, often without due consideration of inter-individual toxicokinetic variability factors, such as age. The purpose of this study was to advance from the existing risk assessment of hazardous substances based on toxicokinetics to a precise human risk assessment by additionally considering the effects of physiologic and metabolic fluctuations and interindividual variability in age. Qualitative age-associated physiologic and metabolic changes in humans, obtained through a meta-analysis, were quantitatively modeled to produce the final biological age algorithm (BAA). The developed BAAs (for males) were extended and applied to the reported testicular reproductive toxicity-focused di-isobutyl-phthalate (DiBP)-mono-isobutyl-phthalate (MiBP) physiologically based toxicokinetic (PBTK) model in males. The advanced PBTK model combined with the BAA was applied to the human risk assessment based on MiBP biomonitoring data. As a result, the specialized DiBP external exposure values for each age could be estimated. Additionally, by applying the Monte Carlo simulation, the distribution of internal exposure diversity among individuals according to the same external exposure dose could be estimated. The contributions of physiologic and metabolic factors to the age-dependent toxicokinetic changes were approximately 93.41-99.99 and 0.01-6.59%, respectively. In addition, the relative contribution of metabolic factors was major in infants and continued to decrease as age increased (up to about age 30 years). This study provides a step-by-step platform that can be widely applied to overcome the limitations of existing toxicokinetic models that still require interindividual pharmacokinetic variability explanations. This will be important for the rationalization and explanation of inter-individual variability in the pharmacokinetics of many substances.

Development of physiologically based toxicokinetic models for 3-monochloropropane-1,2-diol and glycidol.

3-Monochloropropane-1,2-diol (3-MCPD), glycidol, together with their fatty acid esters are commonly presented in various food and have shown carcinogenicity in various laboratory animals. Public health risk assessment of 3-MPCD and glycidol exposure relies on quantitative tools that represent their in vivo toxicokinetics. In order to better understand the absorption, distribution, metabolism, and excretion profiles of 3-MCPD and glycidol in male rats, a physiologically based pharmacokinetic (PBTK) model was developed. The model's predictive power was evaluated by comparing in silico simulations to in vivo time course data obtained from experimental studies. Results indicate that our PBTK model successfully captured the toxicokinetics of both free chemicals in key organs, and their metabolites in accessible biological fluids. With the validated PBTK model, we then gave an animal-free example on how to extrapolate the toxicological knowledge acquired from a single gavage to a realistic dietary intake scenario. Three biomarkers, free compound in serum, urinary metabolite DHPMA, and glycidol-hemoglobin adduct (diHOPrVal) were selected for in silico simulation following constant dietary intakes, and their internal levels were correlated with proposed external daily exposure via reverse dosimetry approaches. Taken together, our model provides a computational approach for extrapolating animal toxicokinetic experiments to biomonitoring measurement and risk assessment.

Enhanced prediction of internal concentrations of phenolic endocrine disrupting chemicals and their metabolites in fish by a physiologically based toxicokinetic incorporating metabolism (PBTK-MT) model.

Bisphenol A (BPA), 4-nonylphenol (4-NP), and triclosan (TCS) are phenolic endocrine disrupting chemicals (EDCs), which are widely detected in aquatic environments and further bioaccumulated and metabolized in fish. Physiologically based toxicokinetic (PBTK) models have been used to describe the absorption, distribution, metabolism, and excretion (ADME) of parent compounds in fish, whereas the metabolites are less explored. In this study, a PBTK incorporating metabolism (PBTK-MT) model for BPA, 4-NP, and TCS was established to enhance the performance of the traditional PBTK model. The PBTK-MT model comprised 16 compartments, showing great accuracy in predicting the internal concentrations of three compounds and their glucuronidated and sulfated conjugates in fish. The impact of typical hepatic metabolism on the PBTK-MT model was successfully resolved by optimizing the mechanism for deriving the partition coefficients between the blood and liver. The PBTK-MT model exhibited a potential data gap-filling capacity for unknown parameters through a backward extrapolation approach of parameters. Model sensitivity analysis suggested that only five parameters were sensitive in at least two PBTK-MT models, while most parameters were insensitive. The PBTK-MT model will contribute to a well understanding of the environmental behavior and risks of pollutants in aquatic biota.

Pregnancy-specific physiologically-based toxicokinetic models for bisphenol A and bisphenol S.

Predictions from physiologically based toxicokinetic (PBTK) models can help inform human health risk assessment for potentially toxic chemicals in the environment. Bisphenol S (BPS) is the second most abundant bisphenol detected in humans in the United States, after bisphenol A (BPA). We have recently demonstrated that BPS, much like BPA, can cross the placental barrier and disrupt placental function. Differences in physicochemical properties, toxicokinetics, and exposure outcomes between BPA and other bisphenols prevent direct extrapolation of existing BPA PBTK models to BPS. The current study aimed to develop pregnancy-specific PBTK (p-PBTK) models for BPA and BPS, using a common p-PBTK model structure. Novel paired maternal and fetal pregnancy data sets for total, unconjugated, and conjugated BPA and BPS plasma concentrations from three independent studies in pregnant sheep were used for model calibration. The nine-compartment (maternal blood, liver, kidney, fat, placenta and rest of body, and fetal liver, blood and rest of body) models simulated maternal and fetal experimental data for both BPA and BPS within one standard deviation for the majority of the experimental data points, highlighting the robustness of both models. Simulations were run to examine fetal exposure following daily maternal exposure to BPA or BPS at their tolerable daily intake dose over a two-week period. These predictive simulations show fetal accumulation of both bisphenols over time. Interestingly, the steady-state approximation following this dosing strategy achieved a fetal concentration of unconjugated BPA to levels observed in cord blood from human biomonitoring studies. These models advance our understanding of bisphenolic compound toxicokinetics during pregnancy and may be used as a quantitative comparison tool in future p-PBTK models for related chemicals.

Bayesian toxicokinetic modeling of cadmium exposure in Chinese population.

Cadmium (Cd) is a toxic heavy metal widely present in the environment. Estimating its internal levels for a given external exposure using toxicokinetic (TK) models is key to the human health risk assessment of Cd. In this study, existing Cd TK models were adapted to develop a one-compartment TK model and a multi-compartment physiologically based toxicokinetic (PBTK) model by estimating the characteristics of Cd kinetics based on Cd exposure data from 814 Chinese residents. Both models not only considered the effect of gender difference on Cd kinetics, but also described the model parameters in terms of distributions to reflect individual variability. For both models, the posterior distributions of sensitive parameters were estimated using the Markov chain-Monte Carlo method (MCMC) and the approximate Bayesian computation-MCMC algorithm (ABC-MCMC). Validation with the test dataset showed 1.4-22.5% improvement in the root mean square error (RMSE) over the original models. After a systematic literature search, the optimized models showed acceptable prediction on other Chinese datasets. The study provides a method for parameter optimization of TK models under different exposure environment, and the validated models can serve as new quantitative assessment tools for the risk assessment of Cd in the Chinese population.

Can the inhalation exposure of a specific worker in a cross-ventilated factory be evaluated by time- and spatial-averaged contaminant concentration?

Industry implies economic growth; however, outdoor and indoor air pollution generated by industrial activities represents a widespread problem for the environment and human beings. In terms of human health, indoor air quality assessment has become essential in a society where people spend most of their time in indoor dwellings, as in the case of industry workers. Because indoor air quality is strongly affected by the outdoor environment, especially under natural ventilation conditions (e.g., cross-ventilation), a comprehensive analysis that includes outdoor atmospheric-urban environment is needed to reproduce realistic scenarios. In this context, computational fluid dynamics (CFD) is a useful tool. To perform a precise analysis of the inhalation exposure of factory workers to potential gas-phase contaminants in the working environment (i.e., inhaled dose of contaminants and potential effects), the human body and respiratory tract need to be integrated in the analysis. Therefore, in this study, we performed an integrated occupational inhalation exposure/toxicology assessment in a factory building that applies a computer simulated person (CSP), a virtual human respiratory tract and integrated physiologically-based toxicokinetic (PBTK) model to predict tissue dosimetry distribution. Outdoor airflow variation was transported into the enclosure through an hourly change in wind pressure coefficient to calculate transient ventilation rate and indoor contaminant concentration between 08:00 and 17:00 h. Thereafter, the time-averaged contaminant concentration calculated at the nares of the human body was employed in a steady state calculation of airflow and contaminant distribution inside the virtual respiratory tract. Subsequently, we predicted adsorbed contaminant in the first layer of tissue of the human airways; highest adsorption took place in the nasal cavity. Finally, based on the results of the comprehensive coupled numerical analysis performed using the CFD-CSP-PBTK model, we quantitatively discussed differences between the inhalation exposure concentration and representative contaminant concentration in the factory space (e.g., time and volume-averaged concentration).

Computational Analysis of Deposition and Translocation of Inhaled Nicotine and Acrolein in the Human Body with E-cigarette Puffing Topographies.

Recently, toxicants such as formaldehyde and acrolein were detected in electronic cigarette (EC) aerosols. It is imperative to conduct research and provide sufficient quantitative evidence to address the associated potential health risks. However, it is still a lack of informative data, i.e., high-resolution local dosimetry of inhaled aerosols in lung airways and other systemic regions, due to the limited imaging resolutions, restricted operational flexibilities, and invasive nature of experimental and clinical studies. In this study, an experimentally validated multiscale numerical model, i.e., Computational Fluid-Particle Dynamics (CFPD) model combined with a Physiologically Based Toxicokinetic (PBTK) model is developed to predict the systemic translocation of nicotine and acrolein in the human body after the deposition in the respiratory system. In-silico parametric analysis is performed for puff topography influence on the deposition and translocation of nicotine and acrolein in human respiratory systems and the systemic region. Results indicate that the puff volume and holding time can contribute to the variations of the nicotine and acrolein plasma concentration due to enhanced aerosol deposition in the lung. The change in the holding time has resulted in significant difference in the chemical translocation which was neglected in a large group of experimental studies. The capability of simulating multiple puffs of the new CFPD-PBTK model paves the way to a valuable computational simulation tool for assessing the chronic health effects of inhaled EC toxicants.

A physiologically based toxicokinetic model for inhaled ethylene and ethylene oxide in mouse, rat, and human.

Ethylene (ET) is the largest volume organic chemical. Mammals metabolize the olefin to ethylene oxide (EO), another important industrial chemical. The epoxide alkylates macromolecules and has mutagenic and carcinogenic properties. In order to estimate the EO burden in mice, rats, and humans resulting from inhalation exposure to gaseous ET or EO, a physiological toxicokinetic model was developed. It consists of the compartments lung, richly perfused tissues, kidneys, muscle, fat, arterial blood, venous blood, and liver containing the sub-compartment endoplasmic reticulum. Modeled ET metabolism is mediated by hepatic cytochrome P450 2E1, EO metabolism by hepatic microsomal epoxide hydrolase or cytosolic glutathione S-transferase in various tissues. EO is also spontaneously hydrolyzed or conjugated with glutathione. The model was validated on experimental data collected in mice, rats, and humans. Modeled were uptake by inhalation, wash-in-wash-out effect in the upper respiratory airways, distribution into tissues and organs, elimination via exhalation and metabolism, and formation of 2-hydroxyethyl adducts with hemoglobin and DNA. Simulated concentration-time courses of ET or EO in inhaled (gas uptake studies) or exhaled air, and of EO in blood during exposures to ET or EO agreed excellently with measured data. Predicted levels of adducts with DNA and hemoglobin, induced by ET or EO, agreed with reported levels. Exposures to 10000 ppm ET were predicted to induce the same adduct levels as EO exposures to 3.95 (mice), 5.67 (rats), or 0.313 ppm (humans). The model is concluded to be applicable for assessing health risks from inhalation exposure to ET or EO.