A physiologically based pharmacokinetics model for florfenicol in crucian carp and oral-to-intramuscular extrapolation

Predicting PFAS fate in fish: Assessing the roles of dietary, respiratory, and dermal uptake in bioaccumulation modeling

An increasing number of per- and polyfluoroalkyl substances (PFAS) exposed to the environment may pose a threat to organisms and human beings. However, there is a lack of simulations comprehensively addressing and comparing the bioaccumulation of PFAS across all three major exposure routes (oral, inhalation, and dermal), especially for dermal uptake. In this study, we proposed a physiologically based kinetic (PBK) model for PFAS, aiming to predict bioaccumulation factors (BAF) in fish by considering these diverse exposure routes. 15 PFAS were used for model validation, and 11 PFAS from Taihu Lake were used for exposure contribution modeling. Approximately 64% of estimations fell within 10-fold model bias from measurements in Taihu Lake, underscoring the potential efficacy of the developed PBK model in predicting BAFs for fish. The dermal route emerges as a contributor to short-chain PFAS exposure. For example, it ranged widely from 46% to 75% (mean) for all modeling short-chain PFAS (C6-C7) in Taihu Lake. It indicated the criticality of considering dermal exposure for PFAS in fish, highlighting a gap in field studies to unravel cutaneous intake mechanisms and contributions. For longer carbon chains of PFAS (C8-C12), dermal exposure accounted for 2%-27% for all species of aquatic organisms. The fish's lipid fraction and water content played a significant role in the contribution of PFAS intake through cutaneous exposure and inhalation. Kow had a significant positive correlation with skin intake rate (p < 0.05) and gill intake rate (p < 0.001), while having a significant negative correlation with skin intake (p < 0.05) and skin intake contribution (p < 0.001). Based on the proposed modeling approach, we have introduced a simulation spreadsheet for projecting PFAS BAFs in fish tissues, hopefully broadening the predictive operational tool for a variety of chemical species.

Development and application of a physiologically based pharmacokinetic model for orbifloxacin in crucian carp (Carassius auratus)

A flow-limited physiologically based pharmacokinetic (PBPK) model consisting of seven compartments was established for orbifloxacin in crucian carp to predict drug concentrations after intravenous or intramuscular injections. Physiological and anatomical parameters, including tissue weights and blood flow through different tissues, were obtained from previous literature. The tissue/plasma partition coefficients for orbifloxacin were calculated using the area method or parameter optimization. In addition, their values were 0.9326, 1.1204, 1.1644, 1.3514, and 2.0057 in the liver, skin, muscle, kidney, and the rest of the body compartment, respectively. Based on the current PBPK model, orbifloxacin concentrations were predicted and compared with those previously reported for further validation. In addition, the mean absolute percentage error (MAPE) values were also calculated, with values ranging from 10.21% in plasma to 42.37% in kidneys, indicating acceptable predictions for all tissues and plasma. A local sensitivity analysis was performed, which showed that the parameters related to elimination and distribution were most influential on orbifloxacin concentrations in muscle. This model was finally used to predict plasma and tissue concentrations after multiple intramuscular dosing. The current PBPK model provided a valuable tool for predicting the tissue residues of orbifloxacin in crucian carp following intramuscular injection.

Application of physiologically based pharmacokinetic models to promote the development of veterinary drugs with high efficacy and safety

Physiologically based pharmacokinetic (PBPK) models have become important tools for the development of novel human drugs. Food-producing animals and pets comprise an important part of human life, and the development of veterinary drugs (VDs) has greatly impacted human health. Owing to increased affordability of and demand for drug development, VD manufacturing companies should have more PBPK models required to reduce drug production costs. So far, little attention has been paid on applying PBPK models for the development of VDs. This review begins with the development processes of VDs; then summarizes case studies of PBPK models in human or VD development; and analyzes the application, potential, and advantages of PBPK in VD development, including candidate screening, formulation optimization, food effects, target-species safety, and dosing optimization. Then, the challenges of applying the PBPK model to VD development are discussed. Finally, future opportunities of PBPK models in designing dosing regimens for intracellular pathogenic infections and for efficient oral absorption of VDs are further forecasted. This review will be relevant to readers who are interested in using a PBPK model to develop new VDs.

Bayesian population physiologically-based pharmacokinetic model for robustness evaluation of withdrawal time in tilapia aquaculture administrated to florfenicol

The antimicrobial residues of aquacultural production is a growing public concern, leading to reexamine the method for establishing robust withdrawal time and ensuring food safety. Our study aims to develop the optimizing population physiologically-based pharmacokinetic (PBPK) model for assessing florfenicol residues in the tilapia tissues, and for evaluating the robustness of the withdrawal time (WT). Fitting with published pharmacokinetic profiles that experimented under temperatures of 22 and 28 °C, a PBPK model was constructed by applying with the Bayesian Markov chain Monte Carol (MCMC) algorithm to estimate WTs under different physiological, environmental and dosing scenarios. Results show that the MCMC algorithm improves the estimates of uncertainty and variability of PBPK-related parameters, and optimizes the simulation of the PBPK model. It is noteworthy that posterior sets generated from temperature-associated datasets to be respectively used for simulating residues under corresponding temperature conditions. Simulating the residues under regulated regimen and overdosing scenarios for Taiwan, the estimated WTs were 12-16 days at 22 °C and 9-12 days at 28 °C, while for the USA, the estimated WTs were 14-18 and 11-14 days, respectively. Comparison with the regulated WT of 15 days, results indicate that the current WT has well robustness and resilience in the environment of higher temperatures. The optimal Bayesian population PBPK model provides effective analysis for determining WTs under scenario-specific conditions. It is a new insight into the increasing body of literature on developing the Bayesian-PBPK model and has practical implications for improving the regulation of food safety.

Development and Application of a Water Temperature Related Physiologically Based Pharmacokinetic Model for Enrofloxacin and Its Metabolite Ciprofloxacin in Rainbow Trout

Enrofloxacin (ENR) has been approved for the treatment of infections in aquaculture, but it may cause tissue residue. This research aimed to develop and validate a water temperature related PBPK model, including both ENR and ciprofloxacin (CIP), in rainbow trout, and to predict further their residue concentrations and the withdrawal periods for ENR at different water temperatures. With the published concentrations data, a flow-limited PBPK model including both ENR and CIP sub-models was developed to predict ENR and CIP concentrations in tissues and plasma/serum after intravenous, oral, or immersion administration. A Monte Carlo simulation including 500 iterations was further incorporated into this model. Based on the model and Monte Carlo analysis, the withdrawal intervals were estimated for different dosage regimens and at different water temperatures, ranging from 80 to 272 degree-days. All of these values were shorter than the labeled withdrawal period (500 degree-days) in fish. This model provided a useful tool for predicting the tissue residues of ENR and CIP in rainbow trout under different dosage regimens and at different water temperatures.

A physiologically based pharmacokinetic model of doxycycline for predicting tissue residues and withdrawal intervals in grass carp (Ctenopharyngodon idella)

The extensive use of doxycycline in aquaculture results in drug residue violations that negatively impact human food safety. This study aimed to develop a physiologically based pharmacokinetic (PBPK) model for doxycycline to predict drug residues and withdrawal times (WTs) in grass carp (Ctenopharyngodon idella) after daily oral administration for 3 days. Physiological parameters including cardiac output and organ weights were measured experimentally. Chemical-specific parameters were obtained from the literature or estimated by fitting to the observed data. The model properly captured the observed kinetic profiles of doxycycline in tissues (i.e., liver, kidney, muscle + skin and gill). The predicted WT in muscle + skin by Monte Carlo analysis based on sensitive parameters identified at 24 h after drug administration was 41 d, which was similar to 43 d calculated using the tolerance limit method. Sensitivity analysis identified two additional sensitive parameters at 6 weeks: intestinal transit rate constant and urinary elimination rate constant. The predicted WT in muscle + skin based on sensitive parameters identified at 6 weeks was 54 d. This model provides a useful tool to estimate tissue residues and withdrawal times for doxycycline in grass carp and also serves a foundation for extrapolation to other fish species and other tetracyclines.

Depletion study and withdrawal period calculation of florfenicol in the crucian carp (Carassius auratus) following multiple intramuscular injections

The previously adopted marker residue for florfenicol (FF) in China was only florfenicol amine (FFA); however, the marker residue has been changed to FF plus FFA since the end of 2017. The previous official withdrawal period determined based on the only concentration of FFA may no longer be suitable. Therefore, the present study aimed to determine the depletion profiles of FF and FFA and further calculate the withdrawal period in the crucian carp (Carassius auratus) based on the new marker residues. Florfenicol was intramuscularly administered at 10 mg/kg bodyweight daily for five consecutive days to crucian carps reared in freshwater at 10°C. After the last dose, plasma and tissue samples were randomly collected from 10 fish at different time points. The FF and FFA concentrations were simultaneously determined by high-performance liquid chromatography (HPLC) with a fluorescence detector and further subjected to noncompartmental analysis. The elimination half-life (h) of FF in different tissues decreased as follows: liver (39.1) > kidney (36.3) > skin plus muscle (34.6) > plasma (31.7), whereas that of FFA decreased as follows: kidney (41.4) > skin plus muscle (39.4) > liver (39.3) > plasma (35.7). Considering a maximum residue limit of 1 μg/g for the total concentration of FF and FFA in the skin plus muscle, a withdrawal period of 6 days was calculated based on the upper limit of the one-sided 95% confidence interval.

Pharmacokinetics of florfenicol in blunt-snout bream (Megalobrama amblycephala) at two water temperatures with single-dose oral administration

The pharmacokinetics and residue elimination of florfenicol (FFC) and its metabolite florfenicol amine (FFA) were studied in healthy blunt-snout bream (Megalobrama amblycephala, 50 ± 10 g). The study was conducted with a single-dose (25 mg/kg) oral administration at a water temperature of 18 or 28°C, while in the residue elimination study, fish were administered at 25 mg/kg daily for three consecutive days by oral gavage to determine the withdrawal period (WDT) at 28°C. The FFC and FFA levels in plasma and tissues (liver, kidneys and muscle) were analysed using high-performance liquid chromatography (HPLC). A no-compartment model was used to analyse the concentration versus time data of M. amblycephala. In the two groups at 18 and 28°C, the maximum plasma concentration (C_max ) of FFC was 5.89 and 6.21 μg/ml, while the time to reach C_max (T_max ) was 5.97 and 2.84 hr, respectively. These suggested that higher temperature absorbed more drug and more quickly at M. amblycephala. And the elimination half-life (T_{1/2 kβ} ) of FFC was calculated as 26.75 and 16.14 hr, while the total body clearance (CL) was 0.09 and 0.15 L kg^-1  hr^-1 , and the areas under the concentration-time curves (AUCs) were 265.87 and 163.31 μg hr/ml, respectively. The difference demonstrated that the elimination rate of FFC in M. amblycephala at 28°C was more quickly than that at 18°C. The results of FFA showed the same trend in tissues of M. amblycephala. After multiple oral doses (25 mg/kg daily for 3 days), the k (eliminate rate constant) of FFA in M. amblycephala muscle was 0.017, the C₀ (initial concentration) was 3.07 mg/kg, and the WDT was 10 days (water temperature 28°C).

Physiologically based toxicokinetics (PBTK) models for pharmaceuticals and personal care products in wild common carp (Cyprinus carpio)

Pharmaceutical and personal care products (PPCPs) are regarded as "pseudo-persistent" pollutants due to their being continuously loaded into the aquatic environment. Physiologically based toxicokinetics (PBTK) models that can quantitatively describe absorption, distribution, metabolism and excretion processes of chemicals in biota are of importance to predict internal exposure (e.g. doses at specific target tissues/organs) from external exposure concentrations. In this study, PBTK models with up to six compartments including brain, liver, kidney, gills, richly perfused tissues and poorly perfused tissues, were developed for predicting internal distribution of 10 PPCPs in wild common carp (Cyprinus carpio). The PBTK predicted concentrations were close to the measured ones, with deviations less than 1 log unit for most of PPCPs. Sensitivity analysis showed that various partition coefficients of the chemicals exerted significant influence on model outputs.

Pharmacokinetics of florfenicol and its metabolite florfenicol amine in crucian carp (Carassius auratus) at three temperatures after one single intramuscular injection

The pharmacokinetic profiles of florfenicol (FF) or florfenicol amine (FFA) in crucian carp were compared at different water temperatures after single intramuscular administration of FF at 10 mg/kg bodyweight. The concentrations of FF and FFA were determined by a high-performance liquid chromatography method, and then, the concentration versus time data were subjected to compartmental analysis using a one-compartment open model. At the water temperatures of 10, 20, and 25°C, the peak concentrations (C_max s) of FF were 2.28, 2.29, and 2.34 μg/ml, respectively, while those of FFA were 0.42, 0.71, and 0.82 μg/ml, respectively. And the absorption half-life (t_1/2ka ) of FF was 0.21, 0.19, and 0.21 hr, while the elimination half-life (t_1/2kel ) was 31.66, 24.77, and 21.48 hr, respectively. For FFA, the formation half-life (t_1/2kf ) was 3.85, 8.97, and 12.43 hr, while the t_1/2kel was 58.34, 30.27, and 21.22 hr, respectively. The results presented here demonstrated that the water temperature had effects on the elimination of both FF and FFA and the formation of FFA. Based on the T > MIC values calculated here, to treat the infections of bacterial with MIC value ≤ 0.5 μg/ml, FF intramuscularly given at 10 mg/kg bodyweight with a 72-hr interval is sufficient at the water temperature of 10°C, while the intervals of 60 and 48 hr were needed at 20 and 25°C, respectively. But to treat bacterial with higher MIC values, more FF or FF at 10 mg/kg BW but with shorter intervals should be intramuscularly given to the infected fish.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

Toxicokinetic models and related tools in environmental risk assessment of chemicals

Environmental risk assessment of chemicals for the protection of ecosystems integrity is a key regulatory and scientific research field which is undergoing constant development in modelling approaches and harmonisation with human risk assessment. This review focuses on state-of-the-art toxicokinetic tools and models that have been applied to terrestrial and aquatic species relevant to environmental risk assessment of chemicals. Both empirical and mechanistic toxicokinetic models are discussed using the results of extensive literature searches together with tools and software for their calibration and an overview of applications in environmental risk assessment. These include simple tools such as one-compartment models, multi-compartment models to physiologically-based toxicokinetic (PBTK) models, mostly available for aquatic species such as fish species and a number of chemical classes including plant protection products, metals, persistent organic pollutants, nanoparticles. Data gaps and further research needs are highlighted.