Physiologically-based pharmacokinetic modeling of benzene in humans: a Bayesian approach

Human biomonitoring of low-level benzene exposures

Historically, benzene has been widely used in a large variety of applications. Occupational exposure limits (OELs) were set for benzene as it was found to be acutely toxic, causing central nervous system depression at high exposures. OELs were lowered when it was discovered that chronic exposure to benzene could cause haematotoxicity. After confirmation that benzene is a human carcinogen causing acute myeloid leukaemia and possibly other blood malignancies, OEL were further lowered. The industrial application of benzene as solvent is almost completely discontinued but it is still used as feedstock for the production of other materials, such as styrene. Occupational exposure to benzene may also occur since it is present in crude oil, natural gas condensate and a variety of petroleum products and because benzene can be formed in combustion of organic material. In the past few years, lower OELs for benzene in the range of 0.05-0.25 ppm have been proposed or were already established to protect workers from benzene-induced cancer. The skin is an important potential route of exposure and relatively more important at lower OELs. Consequently, human biomonitoring - which integrates all exposure routes - is routinely applied to control overall exposure to benzene. Several potential biomarkers have been proposed and investigated. For compliance check of the current low OELs, urinary S-phenylmercapturic acid (S-PMA), urinary benzene and blood benzene are feasible biomarkers. S-PMA appears to be the most promising biomarker but proper validation of biomarker levels corresponding to airborne benzene concentrations below 0.25 ppm are needed.

A Scoping Review of Technologies and Their Applicability for Exposome-Based Risk Assessment in the Oil and Gas Industry

INTRODUCTION

Oil and gas workers have been shown to be at increased risk of chronic diseases including cancer, asthma, chronic obstructive pulmonary disease, and hearing loss, among others. Technological advances may be used to assess the external (e.g. personal sensors, smartphone apps and online platforms, exposure models) and internal exposome (e.g. physiologically based kinetic modeling (PBK), biomonitoring, omics), offering numerous possibilities for chronic disease prevention strategies and risk management measures. The objective of this study was to review the literature on these technologies, by focusing on: (i) evaluating their applicability for exposome research in the oil and gas industry, and (ii) identifying key challenges that may hamper the successful application of such technologies in the oil and gas industry.

METHOD

A scoping review was conducted by identifying peer-reviewed literature with searches in MEDLINE/PubMed and SciVerse Scopus. Two assessors trained on the search strategy screened retrieved articles on title and abstract. The inclusion criteria used for this review were: application of the aforementioned technologies at a workplace in the oil and gas industry or, application of these technologies for an exposure relevant to the oil and gas industry but in another occupational sector, English language and publication period 2005-end of 2019.

RESULTS

In total, 72 articles were included in this scoping review with most articles focused on omics and bioinformatics (N = 22), followed by biomonitoring and biomarkers (N = 20), external exposure modeling (N = 11), PBK modeling (N = 10), and personal sensors (N = 9). Several studies were identified in the oil and gas industry on the application of PBK models and biomarkers, mainly focusing on workers exposed to benzene. The application of personal sensors, new types of exposure models, and omics technology are still in their infancy with respect to the oil and gas industry. Nevertheless, applications of these technologies in other occupational sectors showed the potential for application in this sector.

DISCUSSION AND CONCLUSION

New exposome technologies offer great promise for personal monitoring of workers in the oil and gas industry, but more applied research is needed in collaboration with the industry. Current challenges hindering a successful application of such technologies include (i) the technological readiness of sensors, (ii) the availability of data, (iii) the absence of standardized and validated methods, and (iv) the need for new study designs to study the development of disease during working life.

Using quantitative modeling tools to assess pharmacokinetic bias in epidemiological studies showing associations between biomarkers and health outcomes at low exposures

Biomarkers of exposure can be measured at lower and lower levels due to advances in analytical chemistry. Using these sensitive methods, some epidemiology studies report associations between biomarkers and health outcomes at biomarker levels much below those associated with effects in animal studies. While some of these low exposure associations may arise from increased sensitivity of humans compared with animals or from species-specific responses, toxicology studies with drugs, commodity chemicals and consumer products have not generally indicated significantly greater sensitivity of humans compared with test animals for most health outcomes. In some cases, these associations may be indicative of pharmacokinetic (PK) bias, i.e., a situation where a confounding factor or the health outcome itself alters pharmacokinetic processes affecting biomarker levels. Quantitative assessment of PK bias combines PK modeling and statistical methods describing outcomes across large numbers of individuals in simulated populations. Here, we first provide background on the types of PK models that can be used for assessing biomarker levels in human population and then outline a process for considering PK bias in studies intended to assess associations between biomarkers and health outcomes at low levels of exposure. After providing this background, we work through published examples where these PK methods have been applied with several chemicals/chemical classes - polychlorinated biphenyls (PCBs), perfluoroalkyl substances (PFAS), polybrominated biphenyl ethers (PBDE) and phthalates - to assess the possibility of PK bias. Studies of the health effects of low levels of exposure will be improved by developing some confidence that PK bias did not play significant roles in the observed associations.

The shape of low-concentration dose-response functions for benzene: implications for human health risk assessment

Are dose-response relationships for benzene and health effects such as myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) supra-linear, with disproportionately high risks at low concentrations, e.g. below 1 ppm? To investigate this hypothesis, we apply recent mode of action (MoA) and mechanistic information and modern data science techniques to quantify air benzene-urinary metabolite relationships in a previously studied data set for Tianjin, China factory workers. We find that physiologically based pharmacokinetics (PBPK) models and data for Tianjin workers show approximately linear production of benzene metabolites for air benzene (AB) concentrations below about 15 ppm, with modest sublinearity at low concentrations (e.g. below 5 ppm). Analysis of the Tianjin worker data using partial dependence plots reveals that production of metabolites increases disproportionately with increases in air benzene (AB) concentrations above 10 ppm, exhibiting steep sublinearity (J shape) before becoming saturated. As a consequence, estimated cumulative exposure is not an adequate basis for predicting risk. Risk assessments must consider the variability of exposure concentrations around estimated exposure concentrations to avoid over-estimating risks at low concentrations. The same average concentration for a specified duration is disproportionately risky if it has higher variance. Conversely, if chronic inflammation via activation of inflammasomes is a critical event for induction of MDS and other health effects, then sufficiently low concentrations of benzene are predicted not to cause increased risks of inflammasome-mediated diseases, no matter how long the duration of exposure. Thus, we find no evidence that the dose-response relationship is supra-linear at low doses; instead sublinear or zero excess risk at low concentrations is more consistent with the data. A combination of physiologically based pharmacokinetic (PBPK) modeling, Bayesian network (BN) analysis and inference, and partial dependence plots appears a promising and practical approach for applying current data science methods to advance benzene risk assessment.

Challenges Associated With Applying Physiologically Based Pharmacokinetic Modeling for Public Health Decision-Making

The development and application of physiologically based pharmacokinetic (PBPK) models in chemical toxicology have grown steadily since their emergence in the 1980s. However, critical evaluation of PBPK models to support public health decision-making across federal agencies has thus far occurred for only a few environmental chemicals. In order to encourage decision-makers to embrace the critical role of PBPK modeling in risk assessment, several important challenges require immediate attention from the modeling community. The objective of this contemporary review is to highlight 3 of these challenges, including: (1) difficulties in recruiting peer reviewers with appropriate modeling expertise and experience; (2) lack of confidence in PBPK models for which no tissue/plasma concentration data exist for model evaluation; and (3) lack of transferability across modeling platforms. Several recommendations for addressing these 3 issues are provided to initiate dialog among members of the PBPK modeling community, as these issues must be overcome for the field of PBPK modeling to advance and for PBPK models to be more routinely applied in support of public health decision-making.

Characterization of the human kinetic adjustment factor for the health risk assessment of environmental contaminants

A default uncertainty factor of 3.16 (√10) is applied to account for interindividual variability in toxicokinetics when performing non-cancer risk assessments. Using relevant human data for specific chemicals, as WHO/IPCS suggests, it is possible to evaluate, and replace when appropriate, this default factor by quantifying chemical-specific adjustment factors for interindividual variability in toxicokinetics (also referred to as the human kinetic adjustment factor, HKAF). The HKAF has been determined based on the distributions of pharmacokinetic parameters (e.g., half-life, area under the curve, maximum blood concentration) in relevant populations. This article focuses on the current state of knowledge of the use of physiologically based algorithms and models in characterizing the HKAF for environmental contaminants. The recent modeling efforts on the computation of HKAF as a function of the characteristics of the population, chemical and its mode of action (dose metrics), as well as exposure scenario of relevance to the assessment are reviewed here. The results of these studies, taken together, suggest the HKAF varies as a function of the sensitive subpopulation and dose metrics of interest, exposure conditions considered (route, duration, and intensity), metabolic pathways involved and theoretical model underlying its computation. The HKAF seldom exceeded the default value of 3.16, except in very young children (i.e., <≈ 3 months) and when the parent compound is the toxic moiety. Overall, from a public health perspective, the current state of knowledge generally suggest that the default uncertainty factor is sufficient to account for human variability in non-cancer risk assessments of environmental contaminants.

A calibrated human PBPK model for benzene inhalation with urinary bladder and bone marrow compartments

A physiologically-based pharmacokinetic (PBPK) model of benzene inhalation based on a recent mouse model was adapted to include bone marrow (target organ) and urinary bladder compartments. Empirical data on human liver microsomal protein levels and linked CYP2E1 activities were incorporated into the model, and metabolite-specific conversion rate parameters were estimated by fitting to human biomonitoring data and adjusting for background levels of urinary metabolites. Human studies of benzene levels in blood and breath, and phenol levels in urine were used to validate the rate of human conversion of benzene to benzene oxide, and urinary benzene metabolites from Chinese benzene worker populations provided model validation for rates of human conversion of benzene to muconic acid (MA) and phenylmercapturic acid (PMA), phenol (PH), catechol (CA), hydroquinone (HQ), and benzenetriol (BT). The calibrated human model reveals that while liver microsomal protein and CYP2E1 activities are lower on average in humans compared to mice, the mouse also shows far lower rates of benzene conversion to MA and PMA, and far higher conversion of benzene to BO/PH, and of BO/PH to CA, HQ, and BT. The model also differed substantially from existing human PBPK models with respect to several metabolic rate parameters of importance to interpreting benzene metabolism and health risks in human populations associated with bone marrow doses. The model provides a new methodological paradigm focused on integrating linked human liver metabolism data and calibration using biomonitoring data, thus allowing for model uncertainty analysis and more rigorous validation.

The use of biomonitoring data in exposure and human health risk assessment: benzene case study

Abstract A framework of "Common Criteria" (i.e. a series of questions) has been developed to inform the use and evaluation of biomonitoring data in the context of human exposure and risk assessment. The data-rich chemical benzene was selected for use in a case study to assess whether refinement of the Common Criteria framework was necessary, and to gain additional perspective on approaches for integrating biomonitoring data into a risk-based context. The available data for benzene satisfied most of the Common Criteria and allowed for a risk-based evaluation of the benzene biomonitoring data. In general, biomarker (blood benzene, urinary benzene and urinary S-phenylmercapturic acid) central tendency (i.e. mean, median and geometric mean) concentrations for non-smokers are at or below the predicted blood or urine concentrations that would correspond to exposure at the US Environmental Protection Agency reference concentration (30 µg/m(3)), but greater than blood or urine concentrations relating to the air concentration at the 1 × 10(-5) excess cancer risk (2.9 µg/m(3)). Smokers clearly have higher levels of benzene exposure, and biomarker levels of benzene for non-smokers are generally consistent with ambient air monitoring results. While some biomarkers of benzene are specific indicators of exposure, the interpretation of benzene biomonitoring levels in a health-risk context are complicated by issues associated with short half-lives and gaps in knowledge regarding the relationship between the biomarkers and subsequent toxic effects.

Low-dose metabolism of benzene in humans: science and obfuscation

Benzene is a ubiquitous air pollutant that causes human leukemia and hematotoxic effects. Although the mechanism by which benzene causes toxicity is unclear, metabolism is required. A series of articles by Kim et al. used air and biomonitoring data from workers in Tianjin, China, to investigate the dose-specific metabolism (DSM) of benzene over a wide range of air concentrations (0.03-88.9 p.p.m.). Kim et al. concluded that DSM of benzene is greatest at air concentrations <1 p.p.m. This provocative finding motivated the American Petroleum Institute to fund a study by Price et al. to reanalyze the original data. Although their formal 'reanalysis' reproduced Kim's finding of enhanced DSM at sub-p.p.m. benzene concentrations, Price et al. argued that Kim's methods were inappropriate for assigning benzene exposures to low exposed subjects (based on measurements of urinary benzene) and for adjusting background levels of metabolites (based on median values from the 60 lowest exposed subjects). Price et al. then performed uncertainty analyses under alternative approaches, which led them to conclude that '… the Tianjin data appear to be too uncertain to support any conclusions …' regarding the DSM of benzene. They also argued that the apparent low-dose metabolism of benzene could be explained by 'lung clearance.' In addressing these criticisms, we show that the methods and arguments presented by Price et al. are scientifically unsound and that their results are unreliable.

Application of physiologically based pharmacokinetic models in chemical risk assessment

Post-exposure risk assessment of chemical and environmental stressors is a public health challenge. Linking exposure to health outcomes is a 4-step process: exposure assessment, hazard identification, dose response assessment, and risk characterization. This process is increasingly adopting "in silico" tools such as physiologically based pharmacokinetic (PBPK) models to fine-tune exposure assessments and determine internal doses in target organs/tissues. Many excellent PBPK models have been developed. But most, because of their scientific sophistication, have found limited field application-health assessors rarely use them. Over the years, government agencies, stakeholders/partners, and the scientific community have attempted to use these models or their underlying principles in combination with other practical procedures. During the past two decades, through cooperative agreements and contracts at several research and higher education institutions, ATSDR funded translational research has encouraged the use of various types of models. Such collaborative efforts have led to the development and use of transparent and user-friendly models. The "human PBPK model toolkit" is one such project. While not necessarily state of the art, this toolkit is sufficiently accurate for screening purposes. Highlighted in this paper are some selected examples of environmental and occupational exposure assessments of chemicals and their mixtures.

Translational research to develop a human PBPK models tool kit-volatile organic compounds (VOCs)

Toxicity and exposure evaluations remain the two of the key components of human health assessment. While improvement in exposure assessment relies on a better understanding of human behavior patterns, toxicity assessment still relies to a great extent on animal toxicity testing and human epidemiological studies. Recent advances in computer modeling of the dose-response relationship and distribution of xenobiotics in humans to important target tissues have advanced our abilities to assess toxicity. In particular, physiologically based pharmacokinetic (PBPK) models are among the tools than can enhance toxicity assessment accuracy. Many PBPK models are available to the health assessor, but most are so difficult to use that health assessors rarely use them. To encourage their use these models need to have transparent and user-friendly formats. To this end the Agency for Toxic Substances and Disease Registry (ATSDR) is using translational research to increase PBPK model accessibility, understandability, and use in the site-specific health assessment arena. The agency has initiated development of a human PBPK tool-kit for certain high priority pollutants. The tool kit comprises a series of suitable models. The models are recoded in a single computer simulation language and evaluated for use by health assessors. While not necessarily being state-of-the-art code for each chemical, the models will be sufficiently accurate to use for screening purposes. This article presents a generic, seven-compartment PBPK model for six priority volatile organic compounds (VOCs): benzene (BEN), carbon tetrachloride (CCl(4)), dichloromethane (DCM), perchloroethylene (PCE), trichloroethylene (TCE), and vinyl chloride (VC). Limited comparisons of the generic and original model predictions to published kinetic data were conducted. A goodness of fit was determined by calculating the means of the sum of the squared differences (MSSDs) for simulation vs. experimental kinetic data using the generic and original models. Using simplified solvent exposure assumptions for oral ingestion and inhalation, steady-state blood concentrations of each solvent were simulated for exposures equivalent to the ATSDR Minimal Risk Levels (MRLs). The predicted blood levels were then compared to those reported in the National Health and Nutrition Examination Survey (NHANES). With the notable exception of BEN, simulations of combined oral and inhalation MRLs using our generic VOC model yielded blood concentrations well above those reported for the 95th percentile blood concentrations for the U.S. populations, suggesting no health concerns. When the PBPK tool kit is fully developed, risk assessors will have a readily accessible tool for evaluating human exposure to a variety of environmental pollutants.

Development of a human Physiologically Based Pharmacokinetic (PBPK) Toolkit for environmental pollutants

Physiologically Based Pharmacokinetic (PBPK) models can be used to determine the internal dose and strengthen exposure assessment. Many PBPK models are available, but they are not easily accessible for field use. The Agency for Toxic Substances and Disease Registry (ATSDR) has conducted translational research to develop a human PBPK model toolkit by recoding published PBPK models. This toolkit, when fully developed, will provide a platform that consists of a series of priority PBPK models of environmental pollutants. Presented here is work on recoded PBPK models for volatile organic compounds (VOCs) and metals. Good agreement was generally obtained between the original and the recoded models. This toolkit will be available for ATSDR scientists and public health assessors to perform simulations of exposures from contaminated environmental media at sites of concern and to help interpret biomonitoring data. It can be used as screening tools that can provide useful information for the protection of the public.

The impact of saturable metabolism on exposure-response relations in 2 studies of benzene-induced leukemia

Enzymatic saturation of metabolic pathways is one factor that potentially contributes to the nonlinear exposure-response relations that are frequently reported in occupational epidemiologic studies. The authors propose an approach to explore the contribution of saturable metabolism to previously reported exposure-response relations by integrating predictive models of relevant biomarkers of exposure into the epidemiologic analysis. The approach is demonstrated with 2 studies of leukemia in benzene-exposed workers, one conducted in the Australian petroleum industry (1981-1999) and one conducted in a US rubber hydrochloride production factory in Ohio (1940-1996). The studies differed greatly in their magnitudes and durations of exposure. Substitution of biomarker levels for external estimates of benzene exposure reduced the fold difference of the log relative risk of leukemia per unit of cumulative exposure between the 2 studies by 11%-44%. Nevertheless, a considerable difference in the log relative risk per unit of cumulative exposure remained between the 2 studies, suggesting that exposure misclassification, differences in study design, and potential confounding factors also contributed to the heterogeneity in risk estimates.

The role of exposure reconstruction in occupational human health risk assessment: current methods and a recommended framework

Exposure reconstruction for substances of interest to human health is a process that has been used, with various levels of sophistication, as far back as the 1930s. The importance of robust and high-quality exposure reconstruction has been recognized by many researchers. It has been noted that misclassification of reconstructed exposures is relatively common and can result in potentially significant effects on the conclusions of a human health risk assessment or epidemiology study. In this analysis, a review of the key exposure reconstruction approaches described in over 400 papers in the peer-reviewed literature is presented. These approaches have been critically evaluated and classified according to quantitative, semiquantitative, and qualitative approaches. Our analysis indicates that much can still be done to improve the overall quality and consistency of exposure reconstructions and that a systematic framework would help to standardize the exposure reconstruction process in the future. The seven recommended steps in the exposure reconstruction process include identifying the goals of the reconstruction, organizing and ranking the available data, identifying key data gaps, selecting the best information sources and methodology for the reconstruction, incorporating probabilistic methods into the reconstruction, conducting an uncertainty analysis, and validating the results of the reconstruction. Influential emerging techniques, such as Bayesian data analysis, are highlighted. Important issues that will likely influence the conduct of exposure reconstruction into the future include improving statistical analysis methods, addressing the issue of chemical mixtures, evaluating aggregate exposures, and ensuring transparency with respect to variability and uncertainty in the reconstruction effort.

Bayesian algorithm implementation in a real time exposure assessment model on benzene with calculation of associated cancer risks

The objective of the current study was the development of a reliable modeling platform to calculate in real time the personal exposure and the associated health risk for filling station employees evaluating current environmental parameters (traffic, meteorological and amount of fuel traded) determined by the appropriate sensor network. A set of Artificial Neural Networks (ANNs) was developed to predict benzene exposure pattern for the filling station employees. Furthermore, a Physiology Based Pharmaco-Kinetic (PBPK) risk assessment model was developed in order to calculate the lifetime probability distribution of leukemia to the employees, fed by data obtained by the ANN model. Bayesian algorithm was involved in crucial points of both model sub compartments. The application was evaluated in two filling stations (one urban and one rural). Among several algorithms available for the development of the ANN exposure model, Bayesian regularization provided the best results and seemed to be a promising technique for prediction of the exposure pattern of that occupational population group. On assessing the estimated leukemia risk under the scope of providing a distribution curve based on the exposure levels and the different susceptibility of the population, the Bayesian algorithm was a prerequisite of the Monte Carlo approach, which is integrated in the PBPK-based risk model. In conclusion, the modeling system described herein is capable of exploiting the information collected by the environmental sensors in order to estimate in real time the personal exposure and the resulting health risk for employees of gasoline filling stations.

Reconstructing population exposures to environmental chemicals from biomarkers: challenges and opportunities

A conceptual/computational framework for exposure reconstruction from biomarker data combined with auxiliary exposure-related data is presented, evaluated with example applications, and examined in the context of future needs and opportunities. This framework employs physiologically based toxicokinetic (PBTK) modeling in conjunction with numerical "inversion" techniques. To quantify the value of different types of exposure data "accompanying" biomarker data, a study was conducted focusing on reconstructing exposures to chlorpyrifos, from measurements of its metabolite levels in urine. The study employed biomarker data as well as supporting exposure-related information from the National Human Exposure Assessment Survey (NHEXAS), Maryland, while the MENTOR-3P system (Modeling ENvironment for TOtal Risk with Physiologically based Pharmacokinetic modeling for Populations) was used for PBTK modeling. Recently proposed, simple numerical reconstruction methods were applied in this study, in conjunction with PBTK models. Two types of reconstructions were studied using (a) just the available biomarker and supporting exposure data and (b) synthetic data developed via augmenting available observations. Reconstruction using only available data resulted in a wide range of variation in estimated exposures. Reconstruction using synthetic data facilitated evaluation of numerical inversion methods and characterization of the value of additional information, such as study-specific data that can be collected in conjunction with the biomarker data. Although the NHEXAS data set provides a significant amount of supporting exposure-related information, especially when compared to national studies such as the National Health and Nutrition Examination Survey (NHANES), this information is still not adequate for detailed reconstruction of exposures under several conditions, as demonstrated here. The analysis presented here provides a starting point for introducing improved designs for future biomonitoring studies, from the perspective of exposure reconstruction; identifies specific limitations in existing exposure reconstruction methods that can be applied to population biomarker data; and suggests potential approaches for addressing exposure reconstruction from such data.

Whole body physiologically-based pharmacokinetic models: their use in clinical drug development

BACKGROUND

Whole-body physiologically-based pharmacokinetic (WB-PBPK) models mathematically describe an organism as a closed circulatory system consisting of compartments that represent the organs important for compound absorption, distribution, metabolism and elimination.

OBJECTIVES

To review the current state of WB-PBPK model use in the clinical phases of drug development.

METHODS

A qualitative description of the WB-PBPK model structure is included along with a review of the varying methods available for input parameterisation. Current and potential WB-PBPK model application in clinical development is discussed.

CONCLUSIONS

This modelling tool is at present used for small and large molecule drug development primarily as a means to scale pharmacokinetics from animals to humans based on physiology. The pharmaceutical industry is active in employing these models to clinical drug development although the applications in use now are narrow in comparison to the potential. Expanded integration of WB-PBPK models into the drug development process will only be achieved with staff training, managerial will, success stories and regulatory agency openness.

Whole-body physiologically based pharmacokinetic models

This review summarizes the most recent developments in and applications of physiologically based pharmacokinetic (PBPK) modeling methodology originating from both the pharmaceutical and environmental toxicology areas. It focuses on works published in the last 5 years, although older seminal papers have also been referenced. After a brief introduction to the field and several essential definitions, the main body of the text is structured to follow the major steps of a typical PBPK modeling exercise. Various applications of the methodology are briefly described. The major future trends and perspectives are outlined. The main conclusion from the review of the available literature is that PBPK modeling, despite its obvious potential and recent incremental developments, has not taken the place it deserves, especially in pharmaceutical and drug development sciences.