Estimating danofloxacin withdrawal time in broiler chickens based on physiologically based pharmacokinetics modeling

A brief review on models for birds exposed to chemicals

"A Who's Who of pesticides is therefore of concern to us all. If we are going to live so intimately with these chemicals eating and drinking them, taking them into the very marrow of our bones - we had better know something about their nature and their power."-Rachel Carson, Silent Spring. In her day, Rachel Carson was right: plant protection products (PPP), like all the other chemical substances that humans increasingly release into the environment without further precaution, are among our worst enemies today (Bruhl and Zaller, 2019; Naidu et al., 2021; Tang et al., 2021; Topping et al., 2020). All compartments of the biosphere, air, soil and water, are potential reservoirs within which all species that live there are impaired. Birds are particularly concerned: PPP are recognized as a factor in the decline of their abundance and diversity predominantly in agricultural landscapes. Due to the restrictions on vertebrates testing, in silico-based approaches are an ideal choice alternative given input data are available. This is where the problem lies as we will illustrate in this paper. We performed an extensive literature search covering a long period of time, a wide diversity of bird species, a large range of chemical substances, and as many model types as possible to encompass all our future need to improve environmental risk assessment of chemicals for birds. In the end, we show that poultry species exposed to pesticides are the most studied at the individual level with physiologically based toxicokinetic models. To go beyond, with more species, more chemical types, over several levels of biological organization, we show that observed data are crucially missing (Gilbert, 2011). As a consequence, improving existing models or developing new ones could be like climbing Everest if no additional data can be gathered, especially on chemical effects and toxicodynamic aspects.

Development and Application of a Physiologically Based Pharmacokinetic Model for Diclazuril in Broiler Chickens

Withdrawal periods for diclazuril in broilers have traditionally been determined through regression analysis. However, over the last two decades, the physiologically based pharmacokinetic (PBPK) model has gained prominence as a predictive tool for veterinary drug residues, which offers an alternative method for establishing appropriate withdrawal periods for veterinary drugs. In this current study, a flow-limited PBPK model was developed to predict diclazuril concentrations in broilers following long-duration administration via medicated feed and water. This model consists of nine compartments, including arterial and venous plasma, lung, muscle, skin + fat, kidney, liver, intestine contents, and the rest of the body compartment. Physiological parameters such as tissue weights (V_cxx) and blood flow (Q_cxx) were gathered from published studies, and tissue/plasma partition coefficients (P_xx) were calculated through the area method or parameter optimization. Published diclazuril concentrations were compared to the predicted values, indicating the accuracy and validity of the model. The sensitivity analysis showed that parameters associated with cardiac output, drug absorption, and elimination significantly affected diclazuril concentrations in the muscle. Finally, a Monte Carlo analysis, consisting of 1000 iterations, was conducted to calculate the withdrawal period. Based on the Chinese MRL values, we calculated a withdrawal period of 0 days for both recommended dosing regimens (through mediated water and feed at concentrations of 0.5-1 mg/L and 1 mg/kg, respectively). However, based on the European MRLs, longer periods were determined for the mediated feed dosing route. Our model provides a foundation for scaling other coccidiostats and poultry species.

Quantitative Comparison of Avian and Mammalian Physiologies for Parameterization of Physiologically Based Kinetic Models

Physiologically based kinetic (PBK) models facilitate chemical risk assessment by predicting in vivo exposure while reducing the need for animal testing. PBK models for mammals have seen significant progress, which has yet to be achieved for avian systems. Here, we quantitatively compare physiological, metabolic and anatomical characteristics between birds and mammals, with the aim of facilitating bird PBK model development. For some characteristics, there is considerable complementarity between avian and mammalian species with identical values for the following: blood hemoglobin and hemoglobin concentrations per unit erythrocyte volume together with relative weights of the liver, heart, and lungs. There are also systematic differences for some major characteristics between avian and mammalian species including erythrocyte volume, plasma concentrations of albumin, total protein and triglyceride together with liver cell size and relative weights of the kidney, spleen, and ovary. There are also major differences between characteristics between sexually mature and sexually immature female birds. For example, the relative weights of the ovary and oviduct are greater in sexually mature females compared to immature birds as are the plasma concentrations of triglyceride and vitellogenin. Both these sets of differences reflect the genetic "blue print" inherited from ancestral archosaurs such as the production of large eggs with yolk filled oocytes surrounded by egg white proteins, membranes and a calciferous shell together with adaptions for flight in birds or ancestrally in flightless birds.

Apply a Physiologically Based Pharmacokinetic Model to Promote the Development of Enrofloxacin Granules: Predict Withdrawal Interval and Toxicity Dose

Enrofloxacin (ENR) granules were developed to prevent and control the infections caused by foodborne zoonotic intestinal pathogens in our previous studies. To promote the further development of ENR granules and standardize their usage in pigs, a physiologically based pharmacokinetic (PBPK) model of the ENR granule in pigs was built to determine the withdrawal time (WT) and evaluate the toxicity to pigs. Meanwhile, the population WT was determined by a Monte Carlo analysis to guarantee pork safety. The fitting results of the model showed that the tissue residual concentrations of ENR, ciprofloxacin, and ENR plus ciprofloxacin were all well predicted by the built PBPK model (R² > 0.82). When comparing with the EMA's WT1.4 software method, the final WT (6 d) of the ENR granules in the population of pigs was well predicted. Moreover, by combining the cytotoxicity concentration (225.9 µg/mL) of ENR against pig hepatocytes, the orally safe dosage range (≤130 mg/kg b.w.) of the ENR granules to pigs was calculated based on the validated PBPK model. The well-predicted WTs and a few uses in animals proved that the PBPK model is a potential tool for promoting the judicious use of antimicrobial agents and evaluating the toxicity of the veterinary antimicrobial products.

A comparative pharmacokinetic study of a novel sustained release danofloxacin formulation and the conventional product in rabbits

Sustained release drug formulations are frequently developed to reduce dosage frequency and to improve outcomes of drug therapy. This study evaluates the pharmacokinetic (PK) parameters of a novel injectable danofloxacin (DANO) formulation in comparison with a conventional product in an animal model. A recently synthesized DANO formulation, prepared by incorporation of DANO-loaded mesoporous silica nanoparticles in liposomes and integration of liposomes in chitosan and β-glycerophosphate solution (lipogel) along with the conventional DANO product were injected subcutaneously (SC) in rabbits. Blood samples were collected at specific time points and DANO concentrations in plasma samples were measured. The PK parameters including maximum concentration (Cmax), time to reach Cmax (Tmax), area under the concentration versus time curves (AUC), area under the first moment concentration-time curve (AUMC) and mean residence time (MRT) were studied by non-compartmental analyses. The values of MRT (156.00 ± 20.00 hr), AUC (15.30 ± 3.00 µg mL^-1 per hr) and Tmax (4.70 ± 1.60 hr) for lipogel formulation were higher than those of the conventional product (8.50 ± 3.60 hr, 3.70 ± 2.00 µg mL^-1 per hr and 0.80 ± 0.26 hr, respectively). However, Cmax values for lipogel formulation (0.41 ± 0.15 µg mL^-1) were significantly lower than those of the conventional drug product (0.68 ± 0.09 µg mL^-1). It was concluded that the novel DANO lipogel effectively slowed down the drug absorption and the incorporation of liposomes in hydrogel could be a useful approach to maintain the therapeutic drug level for a longer period; however, more studies are needed in this field.

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.

Danofloxacin pharmacokinetics and tissue residues in Bilgorajska geese

Danofloxacin is a fluoroquinolone developed for veterinary medicine and used in avian species for the treatment of numerous bacterial infections. However, no pharmacokinetic data have been reported in geese. The aim of the study was three-fold: (i) to evaluate the pharmacokinetics of danofloxacin in geese after single oral (PO) and intravenous (IV) administrations; (ii) to define its residue depletion profile in different goose tissues, and (iii) to recreate a multiple-dose simulation in the practical context of large-scale breeding. Twenty-four healthy geese were randomly divided in three groups each composed of eight animals. Group 1 received danofloxacin IV (5 mg/kg) and groups 2 and 3 were treated PO with the same dose. Blood was collected until 24 h (IV; group 1) and 48 h (PO; group 2) after administration. Two animals from group 3 were sacrificed at 6, 10, 24 and 48 h to collect samples of muscle, heart, kidney, liver, and lung. Danofloxacin was quantified in each matrix using a validated high-performance liquid chromatography method with spectrofluorimetric detection and the pharmacokinetic analysis was performed using non-compartmental and compartmental approaches. Danofloxacin showed a moderate elimination half-life (6.61 h), a slow clearance (0.35 mL/g*h) and a large volume of distribution (1.46 mL/g). The peak plasma concentration after PO administration and the time to reach it were 0.96 μg/mL and 1.70 h, respectively. The oral bioavailability was moderate (58%). Higher residue concentration was found in liver and kidney, compared to the other tissues. If the AUC_(0-24) value found in the present study is included in the pharmacokinetic/pharmacodynamic index (AUC_(0-24)/MIC) for the prediction of fluoroquinolones' efficacy, danofloxacin seems to be effective in geese against gram-negative bacteria with a minimum inhibitory concentration (MIC) < 0.076 μg/mL and against S. pneumoniae with a MIC < 0.29 μg/mL after a single PO dose of 5 mg/kg. Liver and kidney showed the highest drug tissue penetration value, with an explorative withdrawal time of 2.6 and 3.8 days, respectively. A practical multiple-dose regimen simulation does not lead to significant plasma drug accumulation.

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.

Physiological parameter values for physiologically based pharmacokinetic models in food-producing animals. Part II: Chicken and turkey

Physiologically based pharmacokinetic (PBPK) models are growing in popularity due to human food safety concerns and for estimating drug residue distribution and estimating withdrawal intervals for veterinary products originating from livestock species. This paper focuses on the physiological and anatomical data, including cardiac output, organ weight, and blood flow values, needed for PBPK modeling applications for avian species commonly consumed in the poultry market. Experimental and field studies from 1940 to 2019 for broiler chickens (1-70 days old, 40 g - 3.2 kg), laying hens (4-15 months old, 1.1-2.0 kg), and turkeys (1 day-14 months old, 60 g -12.7 kg) were searched systematically using PubMed, Google Scholar, ProQuest, and ScienceDirect for data collection in 2019 and 2020. Relevant data were extracted from the literature with mean and standard deviation (SD) being calculated and compiled in tables of relative organ weights (% of body weight) and relative blood flows (% of cardiac output). Trends of organ or tissue weight growth during different life stages were calculated when sufficient data were available. These compiled data sets facilitate future PBPK model development and applications, especially in estimating chemical residue concentrations in edible tissues to calculate food safety withdrawal intervals for poultry.

An open source physiologically based kinetic model for the chicken (Gallus gallus domesticus): Calibration and validation for the prediction residues in tissues and eggs

Xenobiotics from anthropogenic and natural origin enter animal feed and human food as regulated compounds, environmental contaminants or as part of components of the diet. After dietary exposure, a chemical is absorbed and distributed systematically to a range of organs and tissues, metabolised, and excreted. Physiologically based kinetic (PBK) models have been developed to estimate internal concentrations from external doses. In this study, a generic multi-compartment PBK model was developed for chicken. The PBK model was implemented for seven compounds (with log K_ow range -1.37-6.2) to quantitatively link external dose and internal dose for risk assessment of chemicals. Global sensitivity analysis was performed for a hydrophilic and a lipophilic compound to identify the most sensitive parameters in the PBK model. Model predictions were compared to measured data according to dataset-specific exposure scenarios. Globally, 71% of the model predictions were within a 3-fold change of the measured data for chicken and only 7% of the PBK predictions were outside a 10-fold change. While most model input parameters still rely on in vivo experiments, in vitro data were also used as model input to predict internal concentration of the coccidiostat monensin. Future developments of generic PBK models in chicken and other species of relevance to animal health risk assessment are discussed.

Population pharmacokinetics for danofloxacin in the intestinal contents of healthy and infected chickens

Avian pathogenic Escherichia coli could cause localized and systemic infection in the poultry, and danofloxacin is usually used to treat avian colibacillosis through oral administration. To promote prudent use of danofloxacin and reduce the emergence of drug-resistant E. coli strains, it is necessary to understand the population pharmacokinetics (PopPK) of danofloxacin in chicken intestines. In this study, reversed-phase high performance liquid chromatography (HPLC) with fluorescence detection was used to detect the concentrations of danofloxacin in the contents of duodenum, jejunum, and ileum of the healthy and infected chickens after single oral administration (5 mg/kg body weight). Then, the PopPK of danofloxacin in intestines were analyzed using NONMEM software. As a result, a two-compartment PK model best described the time-concentration profile of duodenal, jejunal, and ileal contents. Interestingly, absorption rate (K_a ), distribution volume (V), and clearance (CL) for danofloxacin from duodenal, jejunal to ileal contents were sequentially decreased in the healthy chickens. However, the trend of K_a , V, and CL of danofloxacin was changed dramatically in the intestine of infected chickens. K_a and V of danofloxacin in the jejunum were higher than in the ileum and duodenum. Compared with healthy chickens, K_a and V of danofloxacin in the duodenum decreased significantly, while increased in jejunum, respectively. It has been noted that K_a decreased and V increased in the ileum of infected chickens. Besides, CL in the duodenum, jejunum, and ileum of infected chickens was, respectively, lower than those of healthy chickens. Interestingly, the relative bioavailability (F) of danofloxacin in the ileum was relatively higher in both healthy and infected chickens. In addition, F in the duodenal, jejunal, and ileal contents of infected chickens was respectively higher than healthy chickens. In summary, the PopPK for danofloxacin in infected chicken intestines was quite different from healthy chickens. The absorption, distribution, and clearance of danofloxacin in healthy chickens decreased from duodenum to jejunum and to ileum. Moreover, the pharmacokinetic characteristics in the intestine of infected chickens changed significantly, and the pharmacokinetic characteristics in the ileum can be used as a representative of all intestinal segments.

An integrated experimental and physiologically based pharmacokinetic modeling study of penicillin G in heavy sows

Penicillin G is widely used in food-producing animals at extralabel doses and is one of the most frequently identified violative drug residues in animal-derived food products. In this study, the plasma pharmacokinetics and tissue residue depletion of penicillin G in heavy sows after repeated intramuscular administrations at label (6.5 mg/kg) and 5 × label (32.5 mg/kg) doses were determined. Plasma, urine, and environmental samples were tested as potential antemortem markers for penicillin G residues. The collected new data and other available data from the literature were used to develop a population physiologically based pharmacokinetic (PBPK) model for penicillin G in heavy sows. The results showed that antemortem testing of urine provided potential correlation with tissue residue levels. Based on the United States Department of Agriculture Food Safety and Inspection Service action limit of 25 ng/g, the model estimated a withdrawal interval of 38 days for penicillin G in heavy sows after 3 repeated intramuscular injections at 5 × label dose. This study improves our understanding of penicillin G pharmacokinetics and tissue residue depletion in heavy sows and provides a tool to predict proper withdrawal intervals after extralabel use of penicillin G in heavy sows, thereby helping safety assessment of sow-derived meat products.

Assessing Global Human Exposure to T-2 Toxin via Poultry Meat Consumption Using a Lifetime Physiologically Based Pharmacokinetic Model

Residue depletion of T-2 toxin in chickens after oral gavage at 2.0 mg/kg twice daily for 2 days was determined in this study. A flow-limited physiologically based pharmacokinetic (PBPK) model was developed for lifetime exposure assessment in chickens. The model was calibrated with data from the residue depletion study and then validated with independent data. A local sensitivity analysis was performed, and 16 sensitive parameters were subjected to Monte Carlo analysis. The population PBPK model was applied to estimate daily intake values of T-2 toxin in different countries based on reported consumption factors and the guidance value of 0.25 mg/kg in feed for chickens by the European Food Safety Authority (EFSA). The predicted daily intakes in different countries were all lower than the EFSA's total daily intake, suggesting that the EFSA's guidance value has minimal risk. This model provides a foundation for scaling to other mycotoxins and other food animal species.

In Vivo Pharmacokinetic/Pharmacodynamic Profiles of Danofloxacin in Rabbits Infected With Salmonella typhimurium After Oral Administration

Salmonella typhimurium is a highly transmissible pathogen in rabbits that causes significant losses. Danofloxacin shows excellent efficacy against S. typhimurium infections. However, there are few reports of the pharmacokinetic/pharmacodynamic (PK/PD) modeling of danofloxacin against this pathogen. The aim of this study was to evaluate the in vivo PK/PD relationship of danofloxacin in rabbits infected with S. typhimurium. We used the reduction of bacterial burden in the blood, liver, spleen, and lung as the target PD endpoints, and determined the PK/PD indexes that best correlated with the efficacy and its corresponding magnitude. Danofloxacin was administrated orally to experimentally S. typhimurium-infected rabbits once daily for three successive days. The concentrations of danofloxacin in the serum and the bacterial burden in the blood, liver, spleen, and lung were determined. The PK/PD relationships of danofloxacin against S. typhimurium were evaluated using a Sigmoid E_max model. The results showed that the area under the concentration-time curve from 0 to 24 h/minimum inhibitory concentration (AUC_{24 h}/MIC) ratio correlated well with the in vivo antibacterial effectiveness in different organs, with an r² of 0.8971, 0.9186, 0.9581, and 0.8708 in the blood, liver, spleen, and lung, respectively. The AUC_{24 h}/MIC ratios for the bactericidal effect (3 × Log₁₀ colony forming units/mL reductions) were 121.30, 354.28, 216.64, and 228.66 in the blood, liver, spleen, and lung, respectively, indicating that the in vivo effectiveness of danofloxacin against S. typhimurium using bacterial reduction in different organs as PD endpoints was not identical. This study illustrated that the selection of the target organ for bacterial reduction determination had little effect on best PK/PD parameter determination, but is critical for parameter magnitude calculation in antimicrobial PK/PD modeling, and furthermore, has an impact on the rational dosage optimization process.