PBPK modeling of the cis- and trans-permethrin isomers and their major urinary metabolites in rats

Toxicokinetic model of the pyrethroid pesticide lambda-cyhalothrin, main exposure route and dose reconstruction predictions in agricultural workers

A toxicokinetic model of the pyrethroid insecticide lambda-cyhalothrin (LCT) was developed to relate absorbed doses to urinary cis-3-(2-chloro-3,3,3-trifluoroprop-1-en-1-yl)-2,2-dimethylcyclopropanecarboxylic acid (CFMP) metabolite levels used as a biomarker of exposure. The model then served to reconstruct absorbed doses in agricultural workers and their probability of exceeding the EFSA Acceptable occupational Exposure Level (AOEL). The toxicokinetic model was able to reproduce the temporal profiles of CFMP in the urine of operators spraying pesticides using the optimized model parameters (adjusted to human volunteer data). Modeling also showed that simulation of an inadvertent oral exposure mainly was the exposure scenario giving the best fit to CFMP urinary time-course data in applicators. With the dermal model parameters optimized from data in volunteers, simulation of a dermal exposure in applicators did not allow to reproduce the observed peak excretions and urinary metabolite levels; extremely high applied dermal doses would be required but still simulated dermal penetration rate would remain too slow. Simulation of an inhalation exposure allowed to reproduce the observed time-courses, but with unrealistic air concentrations. For applicators with the highest urinary concentrations, there was a probability of exceeding the AOEL at some points during the biomonitoring period [>50% probability of exceeding for 27% of 24-h samples]; for non-applicator workers the probability of exceeding the AOEL value was very low [corresponding value of 5%]. Furthermore, the median [95% CI] estimates of 10 000 Monte Carlo simulations led to a biological reference value corresponding to the AOEL of 116 [113-119] ng/kg bw/d and 7.5 [7.3-7.7] μg/L. Overall, 7% of applicators and 1% of workers performing weeding and strawberry picking had a probability of exceeding this biological reference value. As a next step, it would be interesting to apply these methods to multiple exposure to various contaminants.

Development and application of a physiologically-based pharmacokinetic model for ractopamine in goats

Physiologically Based Pharmacokinetic (PBPK) models can provide forecasts of the drug residues within the organism. Ractopamine (RAC) is a typical β-agonist. In this study, we developed a PBPK model for RAC in goats. The goal was to predict the distribution of the drug after multiple oral administrations. The preliminary PBPK model for RAC in goats performed well in predicting the drug's distribution in most tissues. In our sensitivity analysis, we found that the parameter of Qclu (Blood Flow Volume through Lungs) had the greatest impact on the RAC concentrations in plasma, liver, and kidney and was the most sensitive parameter. Furthermore, our study aimed to assess the withdrawal time (WT) of RAC in different tissues after RAC long-term exposure in goats. We found that the WT of RAC in the kidney was the longest, lasting for 13  days. Overall, the insights gained from this study have important implications for optimizing drug administration in goats and ensuring appropriate withdrawal times to prevent any potential risks.

PBPK modeling to support risk assessment of pyrethroid exposure in French pregnant women

BACKGROUND

Pyrethroids are widely used pesticides and are suspected to affect children's neurodevelopment. The characterization of pyrethroid exposure during critical windows of development, such as fetal development and prenatal life, is essential to ensure a better understanding of pyrethroids potential effects within the concept of Developmental Origins of Health and Disease.

OBJECTIVE

The aim of this study was to estimate maternal exposure of French pregnant women from biomonitoring data and simulate maternal and fetal internal concentrations of 3 pyrethroids (permethrin, cypermethrin and deltamethrin) using a multi-substance pregnancy-PBPK (physiologically based pharmacokinetics) model. The estimated maternal exposures were compared to newly proposed toxicological reference values (TRV) children specific also called draft child-specific reference value to assess pyrethroid exposure risk during pregnancy i.e. during the in utero exposure period.

METHODS

A pregnancy-PBPK model was developed based on an existing adult pyrethroids model. The maternal exposure to each parent compound of pregnant women of the Elfe (French Longitudinal Study since Childhood) cohort was estimated by reverse dosimetry based on urinary biomonitoring data. To identify permethrin and cypermethrin contribution to their common urinary biomarkers of exposure, an exposure ratio based on biomarkers in hair was tested. Finally, exposure estimates were compared to current and draft child-specific reference values derived from rodent prenatal and postnatal exposure studies.

RESULTS

The main contributor to maternal pyrethroid diet intake is cis-permethrin. In blood, total internal concentrations main contributor is deltamethrin. In brain, the major contributors to internal pyrethroid exposure are deltamethrin for fetuses and cis-permethrin for mothers. Risk is identified only for permethrin when referring to the draft child-specific reference value. 2.5% of the population exceeded permethrin draft child-specific reference value.

CONCLUSIONS

A new reverse dosimetry approach using PBPK model combined with human biomonitoring data in urine and hair was proposed to estimate Elfe pregnant population exposure to a pyrethroids mixture with common metabolites.

Maternal and fetal tissue distribution of α-cypermethrin and permethrin in pregnant CD-1 mice

Pyrethroid insecticides are widely used throughout agriculture and household products. Recent studies suggest that prenatal exposure to these insecticides may adversely affect fetal development; however, little is known about the distribution of these chemicals in pregnant animals. The present study aimed to address this gap in knowledge by investigating the distribution of two commonly used pyrethroid insecticides, permethrin and α-cypermethrin, in maternal and fetal tissues of pregnant CD-1 mice. Dams were dosed from gestational days 6 to 16 via oral gavage with permethrin (1.5, 15, and 50 mg/kg), α-cypermethrin (0.3, 3, and 10 mg/kg), or corn oil vehicle. Pyrethroid levels were determined in gestational day 16 tissues collected 90 min after the final dose was administered. Across maternal tissues, levels of both pyrethroids were the highest in maternal ovaries, followed by liver and brain, respectively. In addition, levels of both pyrethroids in maternal tissues and placenta were significantly higher than those in the fetal body and amniotic fluid, suggesting that these compounds may exhibit low transfer across the mouse placenta. While additional toxicokinetic studies are needed to verify the time course of pyrethroids in the fetal compartment, these findings support investigation into indirect modes of action relevant to the effects of pyrethroids on mammalian fetal development.

Investigation of steatosis profiles induced by pesticides using liver organ-on-chip model and omics analysis

Several studies have reported a correlation between pesticides exposure and metabolic disorders. Dichlorodiphenyltrichloroethane (DDT) and permethrin (PMT), two pesticides highly prevalent in the environment, have been associated to dysregulation of liver lipids and glucose metabolisms and non-alcoholic fatty liver disease (NAFLD). However, the effects of DDT/PMT mixtures and mechanisms mediating their action remain unclear. Here, we used multi-omic to investigate the liver damage induced by DDT, PMT and their mixture in rat liver organ-on-chip. Organ-on-chip allow the reproduction of in vivo-like micro-environment. Two concentrations, 15 and 150 μM, were used to expose the hepatocytes for 24 h under perfusion. The transcriptome and metabolome analysis suggested a dose-dependent effect for all conditions, with a profile close to control for pesticides low-doses. The comparison between control and high-doses detected 266/24, 256/24 and 1349/30 genes/metabolites differentially expressed for DDT150, PMT150 and Mix150 (DDT150/PMT150). Transcriptome modulation reflected liver inflammation, steatosis, necrosis, PPAR signaling and fatty acid metabolism. The metabolome analysis highlighted common signature of three treatments including lipid and carbohydrates production, and a decrease in amino acids and krebs cycle intermediates. Our study illustrates the potential of organ-on-chip coupled to multi-omics for toxicological studies and provides new tools for chemical risk assessment.

Development and Application of a Life-Stage Physiologically Based Pharmacokinetic (PBPK) Model to the Assessment of Internal Dose of Pyrethroids in Humans

To address concerns around age-related sensitivity to pyrethroids, a life-stage physiologically based pharmacokinetic (PBPK) model, supported by in vitro to in vivo extrapolation (IVIVE) was developed. The model was used to predict age-dependent changes in target tissue exposure of 8 pyrethroids; deltamethrin (DLM), cis-permethrin (CPM), trans-permethrin, esfenvalerate, cyphenothrin, cyhalothrin, cyfluthrin, and bifenthrin. A single model structure was used based on previous work in the rat. Intrinsic clearance (CLint) of each individual cytochrome P450 or carboxylesterase (CES) enzyme that are active for a given pyrethroid were measured in vitro, then biologically scaled to obtain in vivo age-specific total hepatic CLint. These IVIVE results indicate that, except for bifenthrin, CES enzymes are largely responsible for human hepatic metabolism (>50% contribution). Given the high efficiency and rapid maturation of CESs, clearance of the pyrethroids is very efficient across ages, leading to a blood flow-limited metabolism. Together with age-specific physiological parameters, in particular liver blood flow, the efficient metabolic clearance of pyrethroids across ages results in comparable to or even lower internal exposure in the target tissue (brain) in children than that in adults in response to the same level of exposure to a given pyrethroid (Cmax ratio in brain between 1- and 25-year old = 0.69, 0.93, and 0.94 for DLM, bifenthrin, and CPM, respectively). Our study demonstrated that a life-stage PBPK modeling approach, coupled with IVIVE, provides a robust framework for evaluating age-related differences in pharmacokinetics and internal target tissue exposure in humans for the pyrethroid class of chemicals.

Synthetic pyrethroids common metabolite 3-phenoxybenzoic acid induces caspase-3 and Bcl-2 mediated apoptosis in human hepatocyte cells

Synthetic pyrethroids are a group of insecticides frequently used in public health and agriculture, and 3-PBA is a common metabolite of them. Although the liver is the primary organ responsible for metabolizing many compounds including pesticides, to the authors' knowledge there have been no studies on the direct hepatotoxic effects of 3-PBA. Therefore, this study aimed to investigate the possible hepatotoxic effects of 3-PBA on a Human Hepatoma Cell Line (HepG2) and the underlying apoptotic mechanisms. Firstly, an LC₅₀ of 1041.242 µM was calculated for 3-PBA by using the WST-1 test with concentrations ranging between 1 µM and 10 mM. Following that, the HepG2 cells in the experimental group were exposed to 3 different concentrations of 3-PBA (1/5 LC₅₀, 1/10 LC₅₀ and 1/20 LC₅₀) for 24 hours. The apoptotic mechanism was evaluated by using flow cytometry, and immunofluorescence assays for Caspase 3 and Bcl-2. In the flow cytometry assay, the total number of apoptotic cells increased in a dose dependent manner (p < 0.05). In the immunofluorescence assay, the Caspase 3 protein showed strong immunoreactivity in the experimental groups, while the reaction to the Bcl-2 protein was minimal. These results demonstrated that 3-PBA has a significant hepatotoxic effect on HepG2 cells and induces apoptosis via the regulation of Caspase-3 and Bcl-2. Furthermore, our results could further the understanding of the fundamental molecular mechanisms of 3-PBA hepatotoxicity. More studies are needed to determine the effects of long-term exposure to 3-PBA and also the molecular mechanisms underlying hepatotoxicity.

Structural Aspects of Potential Endocrine-Disrupting Activity of Stereoisomers for a Common Pesticide Permethrin against Androgen Receptor

Simple Summary

Human exposure to synthetic or naturally occurring endocrine-disrupting compounds (EDCs) contaminating the environment is associated with disruption in endocrine signaling and homeostatic imbalance of hormones. Pyrethroids constitute an important class of extensively used insecticides reported to have endocrine-disrupting activity. Permethrin is one of the most commonly used pyrethroids and exists in isomeric forms. The aim of this study was to investigate and compare the potential endocrine-disrupting activity of permethrin isomers against the androgen receptor (AR). Structural binding studies showed that all permethrin isomer compounds have the potential to compete with native ligand binding in the AR ligand binding pocket. In conclusion, the results of this study suggest that human exposure to commercially produced isomeric forms of permethrin could potentially interfere with the AR function, which may lead to male reproductive dysfunction.

Abstract

Endocrine-disrupting chemicals (EDCs) are a serious global public health and environmental concern. Pyrethroids are insecticide chemicals that are extensively used for crop protection and household purposes but have been identified as EDCs. On account of their ubiquitous environmental presence, human exposure occurs via food, dermal, or inhalation routes and is associated with health problems, including reproductive dysfunction. Permethrin is the most commonly used pyrethroid, and with two chiral centers in its structure, it has four stereoisomeric forms (two enantiomer pairs), i.e., permethrin (1R,3R)-cis, permethrin (1R,3S)-trans, permethrin (1S,3S)-cis, and permethrin (1S,3R)-trans. The current study was performed for predicting the potential endocrine-disrupting activity of the aforementioned four stereoisomers of permethrin against the androgen receptor (AR). The structural binding characterization and binding energy estimations in the AR binding pocket were done using induced fit docking. The structural binding data indicated that all stereoisomers were placed stably in the AR binding pocket and that the estimated binding energy values were comparable to the AR native ligand, except for permethrin (1S,3S)-cis. Furthermore, the commonality in the amino acid interactions to that of the AR native ligand and the binding energy values suggested the potential AR-disrupting activity of all the stereoisomers; however, stereoselective differences were not observed. Taken together, the results suggest that human exposure to permethrin, either as a racemate mixture or in individual stereoisomer form, could potentially interfere with AR function, which may lead to male reproductive dysfunction.

Estimating human exposure to pyrethroids' mixtures from biomonitoring data using physiologically based pharmacokinetic modeling

Human biomonitoring data provide evidence to exposure of environmental chemicals. Physiologically based pharmacokinetic (PBPK) modelling together with an adequate exposure scenario allows to transpose measured concentrations of chemicals or their metabolites into exposure levels, as daily intakes. In France, high levels of urinary pyrethroids metabolites have been measured in populations. Our work aims at estimating the exposure of the French ENNS cohort to mixtures of four pyrethroids (deltamethrin, permethrin, cypermethrin, and cyfluthrin) from the urinary concentrations of five pyrethroids' metabolites commonly measured in biomonitoring studies. We developed a modelling approach based on a global toxicokinetic model that accounts for the cumulative exposure to pyrethroids as some of the metabolites can be shared by several parent compounds and for human inter-individual variability in metabolism. The median of the individual daily intakes was estimated to 8.1 ng/kg bw/day for permethrin, 17.7 ng/kg bw/day for cypermethrin, 20.4 ng/kg bw/day for cyfluthrin and 34.3 ng/kg bw/day for deltamethrin leading to similar weights for the pair permethrin and cypermethrin (36%), cyfluthrin (31%) and deltamethrin (33%) to the cumulative exposure. Accounting for human variability enabled to explain some of the variations in the metabolites' levels within the cohort. The cumulative exposure was then weighted by their toxicities towards three neurotoxic effects to calculate margins of exposure (MOE). Low MOE values were always associated with high measured concentrations of metabolites in urine and the lowest MOEs were observed for the autonomic division. No risks associated with reconstructed mixtures of pyrethroids were expected for the ENNS cohort. Our approach is an asset to analyse the biomarkers of exposure to pyrethroids simultaneously and could be easily adapted to any local or national specificities in pyrethroids' exposure or populations.

Evaluation of Placental Transfer and Tissue Distribution of cis- and Trans-Permethrin in Pregnant Rats and Fetuses Using a Physiological-Based Pharmacokinetic Model

Biomonitoring studies have highlighted the exposure of pregnant women to pyrethroids based on the measurement of their metabolites in urine. Pyrethroids can cross the placental barrier and be distributed in the fetus as some pyrethroids were also measured in the meconium of newborns. Prenatal exposure to pyrethroids is suspected to alter the neurodevelopment of children, and animal studies have shown that early life exposure to permethrin, one of the most commonly used pyrethroid in household applications, can alter the brain development. This study aimed to characterize the fetal permethrin exposure throughout gestation in rats. We developed a pregnancy physiologically based pharmacokinetic (pPBPK) model that describes the maternal and fetal kinetics of the cis- and trans- isomers of permethrin during the whole gestation period. Pregnant Sprague-Dawley rats were exposed daily to permethrin (50 mg/kg) by oral route from the start of gestation to day 20. Permethrin isomers were quantified in the feces, kidney, mammary gland, fat, and placenta in dams and in both maternal and fetal blood, brain, and liver. Cis- and trans-permethrin were quantified in fetal blood and tissues, with higher concentrations for the cis-isomer. The pPBPK model was fitted to the toxicokinetic maternal and fetal data in a Bayesian framework. Several parameters were adjusted, such as hepatic clearances, partition coefficients, and intestinal absorption. Our work allowed to estimate the prenatal exposure to permethrin in rats, especially in the fetal brain, and to quantitatively estimate the placental transfer. These transfers could be extrapolated to humans and be incorporated in a human pPBPK model to estimate the fetal exposure to permethrin from biomonitoring data.

Differential lymphatic versus portal vein uptake of the synthetic pyrethroids deltamethrin and cis-permethrin in rats

The objective of this study was to obtain data on pathways of absorption of the synthetic pyrethroids deltamethrin (DLM) and cis-permethrin (CPM) following oral administration to rats. Adult male Sprague-Dawley rats with cannulated mesenteric lymph ducts and hepatic portal veins were given single doses of either 5 mg/kg DLM or 60 mg/kg CPM via the duodenum and lymph and portal blood samples collected for up to 300 min. The pyrethroid dosing vehicles (5 mL/kg body weight) were either corn oil or glycerol formal. Levels of DLM and CPM in lymph and portal blood samples were determined by high-performance liquid chromatography-mass spectrometry-mass spectrometry. Over the time period studied, levels of both DLM and CPM following administration in either corn oil or glycerol formal were greater in lymph than in portal blood. Lymphatic uptake of both DLM and CPM was enhanced following dosing in glycerol formal than in corn oil. The results of this study suggest that after oral administration to rats, these two pyrethroids are predominantly absorbed via the lymphatic system rather than via portal blood. The data obtained in this study thus support a recently developed physiologically-based pharmacokinetic (PBPK) model to evaluate age-related differences in pyrethroid pharmacokinetics in the rat, where it was assumed that absorption of pyrethroids was predominantly via lymphatic uptake.

Aggregate and cumulative chronic risk assessment for pyrethroids in the French adult population

Pyrethroids are commonly used as insecticides in households, in agriculture or in veterinary and medicinal products. This study aimed to assess cumulative aggregate exposure to cyfluthrin, cypermethrin, deltamethrin and permethrin in adults in France and the associated health risk, and to identify major contributions of exposure sources and routes. External chronic exposures were estimated from dietary and several environmental sources for the oral, inhalation and dermal routes. Internal concentrations of five associated metabolites were simulated with a physiologically-based pharmacokinetic model. The predicted urinary concentrations were in same order of magnitude as those of the French ENNS biomonitoring survey. Dietary exposure, especially from cereals and animal products, was the major source of exposure. For the 1% of adults most highly exposed, dermal exposure to permethrin through medicinal and veterinary products was an important source of exposure. Considering alterations of motor, sensory and autonomic division, all individual margins of exposure were higher than 100, suggesting that no neurotoxic risk associated with the cumulative aggregate exposure to these four pyrethroids is expected for the French adult population.

Toxicokinetics of cis- and trans-Permethrin: Influence of Isomer, Maturation, and Sex

Permethrin exposure of children and adults is widespread in many populations, but knowledge of its relative toxicokinetics (TK) and health risks in immature age groups is lacking. Studies were conducted in rats to determine the influence of immaturity and sex (on plasma and target organ dosimetry of each of the insecticide’s 2 isomers, cis- and trans-permethrin [CIS and TRANS]). Postnatal day 15, 21, and 90 (adult), Sprague Dawley rats were orally administered a graduated series of doses of CIS and TRANS in corn oil. Serial sacrifices were conducted over 24 h to obtain plasma, brain, liver, skeletal muscle, and fat profiles of CIS and TRANS. Levels of TRANS decreased relatively rapidly, despite administration of relatively high doses. Concentrations of each isomer in plasma, brain, and other tissues monitored were inversely proportional to the animals’ age. The youngest pups exhibited 4-fold higher plasma and brain area under the curves than did adults. Little difference was observed in the TK of CIS or TRANS between adult male and female rats, other than higher initial plasma and liver CIS levels in females. Elevated exposure of the immature brain appears to be instrumental in increased susceptibility to the acute neurotoxicity of high-dose permethrin (Cantalamessa [1993]), but it remains to be established whether age-dependent TK is relevant to long-term, low-level risks.

Determination of maternal and foetal distribution of cis- and trans-permethrin isomers and their metabolites in pregnant rats by liquid chromatography tandem mass spectrometry (LC-MS/MS)

We developed a method to quantify cis-permethrin and trans-permethrin and their metabolites in several biological matrices in pregnant rats and foetuses using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The objective was to quantify cis-permethrin and trans-permethrin in faeces, kidney, mammary gland, fat and placenta in mothers and in both maternal and foetal blood, brain and liver. The metabolites cis-3-(2,2-dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (cis-DCCA), trans-3-(2,2-dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (trans-DCCA) and 3-phenoxybenzoic acid (3-PBA) were measured in blood, liver and urine. Sample preparation was performed by liquid-liquid extraction. A purification step was not carried out except for the more complex biological samples (fat, mammary glands and faeces). Validation parameters including specificity, linearity, matrix effect, limits of quantification (LOQs), accuracy and precision were evaluated. The recoveries of target compounds ranged from 47 to 136%. LOQs were in the range 4 to 80 ng/mL for permethrin isomers and 4 to 800 ng/mL for their respective metabolites. Intra- and inter-batch precision and accuracy in matrix were better than 15%. The validated method was applied in a preliminary toxicokinetic study in pregnant rats with oral dosing of 50 mg/kg permethrin. In pregnant rats, permethrin isomers and their metabolites were quantified in all requested matrices except maternal liver and blood for trans-permethrin and cis-DCCA respectively. In foetuses, cis- and trans-permethrin were also quantified, demonstrating that the method is suitable for the analysis of foetal distribution of permethrin in toxicokinetic studies.

Estimating the cumulative human exposures to pyrethroids by combined multi-route PBPK models: Application to the French population

Human biomarkers of exposure to pyrethroid insecticides are usually urinary concentrations of metabolites that can be specific to a pyrethroid or common to several compounds. We developed a global toxicokinetic model that links the external exposure to four widely-used pyrethroids and their isomers (deltamethrin and cis and trans isomers of permethrin, cypermethrin, and cyfluthrin) to the urinary concentrations of metabolites (cis- and trans-DCCA, 3-PBA, F-PBA and DBCA). This global model includes physiologically based pharmacokinetic models for each parent compound and one-compartment models for the metabolites. Existing in vivo, in vitro and in silico data were used for model calibration, and human toxicokinetic data for model evaluation. Overall, the global model reproduced the data accurately as about 90% of predictions were inside the 3-fold error interval. A sensitivity analysis showed that the most influent parameter for each urinary metabolite concentration was the fraction of parent compound that is transformed into that metabolite. The global model was then tested with realistic exposures for the French population: the predictions were consistent with biomonitoring data. The global model is a tool that will improve the interpretation of biomonitoring data for pyrethroids.

Dose reconstruction in workers exposed to two major pyrethroid pesticides and determination of biological reference values using a toxicokinetic model

A toxicokinetic model has been optimized to describe the time profiles of common biomarkers of exposure to permethrin and cypermethrin: trans-3-(2,2-dichlorovinyl)-2,2-dimethylcyclopropane-1-carboxylic acids (trans-DCCA) and 3-phenoxybenzoic acid (3-PBA). The model then served to reproduce urinary time courses in exposed agricultural workers and predict corresponding absorbed doses. It allowed for the prediction of the main routes of exposure in workers during the study period. Modeling showed that simulating exposure mostly by the oral route, during the 3-day biomonitoring period, provided best-fits to the urinary time courses of most workers. This is compatible with an inadvertent oral exposure during work. According to best-fit scenarios, absorbed doses in workers reconstructed with the model reached a maximum of 2.4 µg/kg bw/day and were below the absorbed dose limits associated with an exposure to the reference dose values established by the U.S. Environmental Protection Agency (0.06 and 0.25 mg/kg bw/day for cypermethrin and permethrin, respectively) and the Acceptable Operator Exposure Level set by the European Commission (0.06 mg/kg bw/day for cypermethrin). Modeling was further used to derive biological reference values for cypermethrin and permethrin exposure. Respective values of 7 and 29 nmol/kg bw/day of trans-DCCA, and 3 and 13 nmol/kg bw/day of 3-PBA were obtained. None of the workers presented values above these biological reference values.

Utilization of data below the analytical limit of quantitation in pharmacokinetic analysis and modeling: promoting interdisciplinary debate

Traditionally, bioanalytical laboratories do not report actual concentrations for samples with results below the LOQ (BLQ) in pharmacokinetic studies. BLQ values are outside the method calibration range established during validation and no data are available to support the reliability of these values. However, ignoring BLQ data can contribute to bias and imprecision in model-based pharmacokinetic analyses. From this perspective, routine use of BLQ data would be advantageous. We would like to initiate an interdisciplinary debate on this important topic by summarizing the current concepts and use of BLQ data by regulators, pharmacometricians and bioanalysts. Through introducing the limit of detection and evaluating its variability, BLQ data could be released and utilized appropriately for pharmacokinetic research.