Development of a physiologically based pharmacokinetic model for flunixin in cattle (Bos taurus)

Application of a generic PBK model for beef cattle: Tissue/fluid correlations of paracetamol and NSAIDs

Nonsteroidal Anti-Inflammatory Drugs (NSAIDs) and paracetamol can be administered off-label to cattle. Since the use of these veterinary medicines in cattle may pose a public health risk after meat consumption, it is important to translate measured concentrations in urine and tissues into concentrations in meat for human consumption. A generic physiologically-based kinetic (PBK) model for cattle can enable this translation. In this work, a beef cattle PBK model was applied to calculate the relationships between concentrations in different bovine tissues and those were compared to measured concentrations in different matrices. Sixty-seven kidney samples, the corresponding urine and meat samples, and available 19 serum samples were analysed. Overall, 70% of the PBK model predictions are within a 2-fold factor and relationships for kidney/meat, urine/meat, and plasma/meat ratios were established. The conversions of measured kidney concentrations into meat concentrations were mostly within a factor two, while those based on plasma and urine were underpredicted. Based on these ratios, plasma and urine could be used as an appropriate surrogate matrix for a fast, simple in vivo sample screening test under field conditions, such as in local farms and slaughterhouses, to predict a maximum residue level exceedance in meat, reducing the number of test samples.

A physiologically based pharmacokinetic/pharmacodynamic model to determine dosage regimens and withdrawal intervals of aditoprim against Streptococcus suis

Introduction: Streptococcus suis (S. suis) is a zoonotic pathogen threatening public health. Aditoprim (ADP), a novel veterinary medicine, exhibits an antibacterial effect against S. suis. In this study, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model was used to determine the dosage regimens of ADP against S. suis and withdrawal intervals. Methods: The PBPK model of ADP injection can predict drug concentrations in plasma, liver, kidney, muscle, and fat. A semi-mechanistic pharmacodynamic (PD) model, including susceptible subpopulation and resistant subpopulation, is successfully developed by a nonlinear mixed-effect model to evaluate antibacterial effects. An integrated PBPK/PD model is conducted to predict the time-course of bacterial count change and resistance development under different ADP dosages. Results: ADP injection, administrated at 20 mg/kg with 12 intervals for 3 consecutive days, can exert an excellent antibacterial effect while avoiding resistance emergence. The withdrawal interval at the recommended dosage regimen is determined as 18 days to ensure food safety. Discussion: This study suggests that the PBPK/PD model can be applied as an effective tool for the antibacterial effect and safety evaluation of novel veterinary drugs.

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

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

The construction and application of physiologically based toxicokinetic models for acrylamide, glycidamide and their biomarkers in rats and humans

Acrylamide (AA), a class 2A probable carcinogen to humans classified by the International Agency for Research on Cancer, has attracted extensive attention worldwide since it was widely used in industrial and domestic water treatment and detected in thermal processing foods. The metabolic adducts of AA and its primary metabolite glycidamide (GA) have been served as biomonitoring markers of AA intake, but the physiologically based toxicokinetics (PBTK) models to estimate internal dosimetry still remain unclear. An updated PBTK model for AA, GA and their metabolic biomarkers in rats and humans was developed and extended with time-course datasets from both literatures and our experiments. With adjustments to the model parameters, linear regression correlation coefficient (R²) between the fitting values and the validation datasets of rats and humans was greater than 0.76. The current model fits well with the experimental datasets of urinary N-acetyl-S-(2-carbamoylethyl)-l-cysteine (AAMA) and (N-(R,S)-acetyl-S-(carbamoyl-2-hydroxyethyl)-l-cysteine) (GAMA) of rats exposed to AA from 0.1 to 50 mg/kg b.w. and humans exposed to AA from 0.0005 to 0.020 mg/kg b.w., indicating the robustness of the current models. Parameters for adduct of AA with N-terminal valine of hemoglobin (AAVal) were extended to humans and validated. Kinetic parameters for rats were assessed and validated based upon fit to the experimental datasets for liver N3-(2-carbamoyl-2-hydroxyethyl)-adenine (N3-GA-Ade) and N7-(2-carbamoyl-2-hydroxyethyl)-guanine (N7-GA-Gua) adducts. Compared with the previous model, the developed model included the correlation between AA intake and its mercapturic acid adducts, AAMA and GAMA, in a larger dose range with new experimental data, and parameters for AAVal, N3-GA-Ade and N7-GA-Gua were improved and verified. The current multi-component PBTK models provide a superior foundation for the estimation of short-term to medium and long-term intake levels of human exposure to AA.

Assessment of pain associated with bovine respiratory disease and its mitigation with flunixin meglumine in cattle with induced bacterial pneumonia

Pleuritic chest pain from bacterial pneumonia is often reported in human medicine. However, studies investigating pain associated with bovine respiratory disease (BRD) are lacking. The objectives of this study were to assess if bacterial pneumonia elicits a pain response in calves with experimentally induced BRD and to determine the analgesic effects of transdermally administered flunixin. Twenty-six calves, 6-7 months of age, with no history of BRD were enrolled into 1 of 3 treatment groups: (1) experimentally induced BRD + transdermal flunixin at 3.3 mg/kg twice, 24 h apart (BRD + FTD); (2) experimentally induced BRD + placebo (BRD + PLBO); and (3) sham induction + placebo (CNTL + PLBO). Calves induced with BRD were inoculated with Mannheimia haemolytica via bronchoalveolar lavage. Outcomes were collected from -48 to 192 hours post-treatment and included serum cortisol; infrared thermography; mechanical nociceptive threshold; substance P; kinematic gait analysis; visual analog scale (VAS); clinical illness score; computerized lung score; average activity and rumination level; prostaglandin E2 metabolite; plasma serum amyloid A and rectal temperature. Outcomes were evaluated using either a generalized logistic mixed model for categorical variables or a generalized linear mixed model for continuous variables. Right front force differed by treatment (P = 0.01). The BRD + PLBO had lower mean force applied to the right front limb (85.5 kg) compared to BRD + FTD (96.5 kg) (P < 0.01). Average VAS differed by a treatment by time interaction (P = 0.01). The VAS scores differed for BRD + PLBO at -48 (3.49 mm) compared to 168 and 192 h (13.49 and 13.64 mm, respectively) (P < 0.01). Activity for BRD + PLBO was higher at -48 h (27 min/h) compared to 48, 72, 120 and 168 h (≤ 22.24 min/h) (P < 0.01). Activity differed by a treatment by time interaction (P = 0.01). Activity for BRD + FTD was higher at -48 and 0 h (28.2 and 28.2 min/h, respectively) compared to 48, 72, 96 and 168 h (≤ 23.7 min/h) (P < 0.01). Results show a combination of reduced activity levels, decreased force on the right front limb, and increased visual analog scale pain scores all support that bacterial pneumonia in cattle is painful. Differences in right front force indicate that flunixin transdermal may attenuate certain pain biomarkers in cattle with BRD. These findings suggest that BRD is painful and analgesic drugs may improve the humane aspects of care for cattle with BRD.

A comprehensive review on non-clinical methods to study transfer of medication into breast milk - A contribution from the ConcePTION project

Breastfeeding plays a major role in the health and wellbeing of mother and infant. However, information on the safety of maternal medication during breastfeeding is lacking for most medications. This leads to discontinuation of either breastfeeding or maternal therapy, although many medications are likely to be safe. Since human lactation studies are costly and challenging, validated non-clinical methods would offer an attractive alternative. This review gives an extensive overview of the non-clinical methods (in vitro, in vivo and in silico) to study the transfer of maternal medication into the human breast milk, and subsequent neonatal systemic exposure. Several in vitro models are available, but model characterization, including quantitative medication transport data across the in vitro blood-milk barrier, remains rather limited. Furthermore, animal in vivo models have been used successfully in the past. However, these models don't always mimic human physiology due to species-specific differences. Several efforts have been made to predict medication transfer into the milk based on physicochemical characteristics. However, the role of transporter proteins and several physiological factors (e.g., variable milk lipid content) are not accounted for by these methods. Physiologically-based pharmacokinetic (PBPK) modelling offers a mechanism-oriented strategy with bio-relevance. Recently, lactation PBPK models have been reported for some medications, showing at least the feasibility and value of PBPK modelling to predict transfer of medication into the human milk. However, reliable data as input for PBPK models is often missing. The iterative development of in vitro, animal in vivo and PBPK modelling methods seems to be a promising approach. Human in vitro models will deliver essential data on the transepithelial transport of medication, whereas the combination of animal in vitro and in vivo methods will deliver information to establish accurate in vitro/in vivo extrapolation (IVIVE) algorithms and mechanistic insights. Such a non-clinical platform will be developed and thoroughly evaluated by the Innovative Medicines Initiative ConcePTION.

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.

Population physiologically based modeling of pirlimycin milk concentrations in dairy cows

Predictions of drug residues in milk are critical in food protection and are a major consideration in the economics of treatment of mastitis in dairy cows. Nonlinear mixed-effects modeling (NLME) has been advocated as a suitable pharmaco-statistical method for the study of drug residues in milk. Recent developments in physiologically based pharmacokinetic (PBPK) modeling of intramammary drugs allow the combination of a mechanistic description of milk pharmacokinetics with NLME methods. The PBPK model was applied to NLME analysis of a data set consisting of milk drug concentrations from 78 healthy cows and 117 with clinical mastitis. Pirlimycin milk pharmacokinetics were adequately described by the model across the range of observed concentrations. Mastitis was characterized by increased variance in milk production volume. Udder residual volume was larger in cows with 1, or 2 or greater diseased mammary glands than in the healthy cows. Low-producing cows had a greater risk of prolonged milk residues. With the exclusion of the low-production cows, the model predicted that healthy cows required a milk discard time 12 h longer than that indicated by the label, and the diseased cows 36 h longer than indicated by the label. More pirlimycin was systemically absorbed in the gram-positive infected compared with the gram-negative infected or healthy cows, suggesting a greater risk of violative meat residues in gram-positive infected cows. Using NLME and PBPK models, we identified factors associated with changes in pirlimycin milk residues that may affect food safety. This model extends the verification of a simple physiologically based framework for the study of intramammary drugs.

A study to assess the correlation between plasma, oral fluid and urine concentrations of flunixin meglumine with the tissue residue depletion profile in finishing-age swine

Background

Flunixin meglumine (FM) was investigated for the effectiveness of plasma, oral fluid, and urine concentrations to predict tissue residue depletion profiles in finishing-age swine, along with the potential for untreated pigs to acquire tissue residues following commingled housing with FM-treated pigs. Twenty pigs were housed in groups of three treated and one untreated control. Treated pigs received one 2.2 mg/kg dose of FM intramuscularly. Before treatment and at 1, 3, 6, 12, 24, 36, and 48 h (h) after treatment, plasma samples were taken. At 1, 4, 8, 12 and 16 days (d) post-treatment, necropsy and collection of plasma, urine, oral fluid, muscle, liver, kidney, and injection site samples took place. Analysis of flunixin concentrations using liquid chromatography/tandem mass spectrometry was done. A published physiologically based pharmacokinetic (PBPK) model for flunixin in cattle was extrapolated to swine to simulate the measured data.

Results

Plasma concentrations of flunixin were the highest at 1 h post-treatment, ranging from 1534 to 7040 ng/mL, and were less than limit of quantification (LOQ) of 5 ng/mL in all samples on Day 4. Flunixin was detected in the liver and kidney only on Day 1, but was not found 4–16 d post-treatment. Flunixin was either not seen or found less than LOQ in the muscle, with the exception of one sample on Day 16 at a level close to LOQ. Flunixin was found in the urine of untreated pigs after commingled housing with FM-treated pigs. The PBPK model adequately correlated plasma, oral fluid and urine concentrations of flunixin with residue depletion profiles in liver, kidney, and muscle of finishing-age pigs, especially within 24 h after dosing.

Conclusions

Results indicate untreated pigs can be exposed to flunixin by shared housing with FM-treated pigs due to environmental contamination. Plasma and urine samples may serve as less invasive and more easily accessible biological matrices to predict tissue residue statuses of flunixin in pigs at earlier time points (≤24 h) by using a PBPK model.

Fate of the nonsteroidal, anti-inflammatory veterinary drug flunixin in agricultural soils and dairy manure

A large percentage of flunixin, a nonsteroidal anti-inflammatory drug widely used for treating livestock, is excreted in intact form and thus potentially available for environmental transport. As the fate of flunixin in the environment is unknown, our objective was to quantify sorption, desorption, and transformation in five agricultural soils and manure using batch equilibrium methods. Concentrations of flunixin and degradation products were determined by high performance liquid chromatography time of flight mass spectrometry. For all studied soils, sorption of flunixin exhibited linear character, with both linear and Freundlich models providing adequate fit. Linear sorption coefficients varied from 8 to 112 L kg^-1. The strongest Pearson correlations with sorption coefficients were for clay content (r = 0.8693), total nitrogen (r = 0.7998), and organic carbon (r = 0.6291). Desorption of the reversibly bound fraction (3-10% of total sorbed flunixin) from all five studied soils exhibited non-hysteretic character suggesting low affinity of this fraction of flunixin to soil. Flunixin degradation in soils was relatively slow, exhibiting half-lives of 39-203 days, thus providing time for off-site transport and environmental contamination. The biological impacts of flunixin at environmentally relevant concentrations must be determined given its environmental behavior and extensive use as a nonsteroidal anti-inflammatory drug in livestock. Graphical abstract.

Physiological parameter values for physiologically based pharmacokinetic models in food-producing animals. Part I: Cattle and swine

Physiologically based pharmacokinetic (PBPK) models for chemicals in food animals are a useful tool in estimating chemical tissue residues and withdrawal intervals. Physiological parameters such as organ weights and blood flows are an important component of a PBPK model. The objective of this study was to compile PBPK‐related physiological parameter data in food animals, including cattle and swine. Comprehensive literature searches were performed in PubMed, Google Scholar, ScienceDirect, and ProQuest. Relevant literature was reviewed and tables of relevant parameters such as relative organ weights (% of body weight) and relative blood flows (% of cardiac output) were compiled for different production classes of cattle and swine. The mean and standard deviation of each parameter were calculated to characterize their variability and uncertainty and to allow investigators to conduct population PBPK analysis via Monte Carlo simulations. Regression equations using weight or age were created for parameters having sufficient data. These compiled data provide a comprehensive physiological parameter database for developing PBPK models of chemicals in cattle and swine to support animal‐derived food safety assessment. This work also provides a basis to compile data in other food animal species, including goats, sheep, chickens, and turkeys.

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.

Integration of Food Animal Residue Avoidance Databank (FARAD) empirical methods for drug withdrawal interval determination with a mechanistic population-based interactive physiologically based pharmacokinetic (iPBPK) modeling platform: example for flunixin meglumine administration

Violative chemical residues in animal-derived food products affect food safety globally and have impact on the trade of international agricultural products. The Food Animal Residue Avoidance Databank program has been developing scientific tools to provide appropriate withdrawal interval (WDI) estimations after extralabel drug use in food animals for the past three decades. One of the tools is physiologically based pharmacokinetic (PBPK) modeling, which is a mechanistic-based approach that can be used to predict tissue residues and WDIs. However, PBPK models are complicated and difficult to use by non-modelers. Therefore, a user-friendly PBPK modeling framework is needed to move this field forward. Flunixin was one of the top five violative drug residues identified in the United States from 2010 to 2016. The objective of this study was to establish a web-based user-friendly framework for the development of new PBPK models for drugs administered to food animals. Specifically, a new PBPK model for both cattle and swine after administration of flunixin meglumine was developed. Population analysis using Monte Carlo simulations was incorporated into the model to predict WDIs following extralabel administration of flunixin meglumine. The population PBPK model was converted to a web-based interactive PBPK (iPBPK) framework to facilitate its application. This iPBPK framework serves as a proof-of-concept for further improvements in the future and it can be applied to develop new models for other drugs in other food animal species, thereby facilitating the application of PBPK modeling in WDI estimation and food safety assessment.

Transport of Potential Manure Hormone and Pharmaceutical Contaminants through Intact Soil Columns

Although adding manure to agricultural soils is a commonly practiced disposal method and a means to enhance soil productivity, potential environmental contamination by any associated chemicals of emerging concern (CECs) such as hormones and pharmaceuticals is not well understood. Our objective was to provide field-relevant predictions of soil transport and attenuation of 19 potential manure CECs using undisturbed soil columns irrigated under unsaturated conditions. The CEC concentrations in leached water were monitored for 13 wk using high performance liquid chromatography-time of flight-mass spectrometry (HPLC-TOF-MS), after which time soil in the cores was removed and sampled for extractable CECs. Compounds quantified in column leachate included all four of the added sulfonamide antibiotics and the nonsteroidal, anti-inflammatory drug flunixin. Only trace amounts of several of the seven hormones, five remaining antibiotics, and two antimicrobials leached from the columns from exogenous soil additions. Soil residues of all 19 compounds were detected, with highest extractable amounts for 17α-hydroxyprogesterone > triclosan (antimicrobial) > flunixin > oxytetracycline. Those CECs with the highest recoveries as calculated by summing leached and extractable amounts were flunixin (14.5%), 17α-hydroxyprogesterone (5.3%), triclosan (4.6%), and sulfadimethoxine (4.8%). Manure management to prevent CEC contamination should consider the potential environmental problems caused by negatively charged compounds with the greatest mobility (flunixin and sulfadimethoxine) and those that have long residence times in soil (triclosan, 17α-hydroxyprogesterone, flunixin, and oxytetracycline). Flunixin is particularly important given its mobility and long residence time in soil.

Predicting mammalian metabolism and toxicity of pesticides in silico

Pesticides must be effective to be commercially viable but they must also be reasonably safe for those who manufacture them, apply them, or consume the food they are used to produce. Animal testing is key to ensuring safety, but it comes late in the agrochemical development process, is expensive, and requires relatively large amounts of material. Surrogate assays used as in vitro models require less material and shift identification of potential mammalian toxicity back to earlier stages in development. Modern in silico methods are cost-effective complements to such in vitro models that make it possible to predict mammalian metabolism, toxicity and exposure for a pesticide, crop residue or other metabolite before it has been synthesized. Their broader use could substantially reduce the amount of time and effort wasted in pesticide development. This contribution reviews the kind of in silico models that are currently available for vetting ideas about what to synthesize and how to focus development efforts; the limitations of those models; and the practical considerations that have slowed development in the area. Detailed discussions are provided of how bacterial mutagenicity, human cytochrome P450 (CYP) metabolism, and bioavailability in humans and rats can be predicted. © 2018 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Mathematical modeling and simulation in animal health. Part III: Using nonlinear mixed-effects to characterize and quantify variability in drug pharmacokinetics

A common feature of human and veterinary pharmacokinetics is the importance of identifying and quantifying the key determinants of between-patient variability in drug disposition and effects. Some of these attributes are already well known to the field of human pharmacology such as bodyweight, age, or sex, while others are more specific to veterinary medicine, such as species, breed, and social behavior. Identification of these attributes has the potential to allow a better and more tailored use of therapeutic drugs both in companion and food-producing animals. Nonlinear mixed effects (NLME) have been purposely designed to characterize the sources of variability in drug disposition and response. The NLME approach can be used to explore the impact of population-associated variables on the relationship between drug administration, systemic exposure, and the levels of drug residues in tissues. The latter, while different from the method used by the US Food and Drug Administration for setting official withdrawal times (WT) can also be beneficial for estimating WT of approved animal drug products when used in an extralabel manner. Finally, NLME can also prove useful to optimize dosing schedules, or to analyze sparse data collected in situations where intensive blood collection is technically challenging, as in small animal species presenting limited blood volume such as poultry and fish.

Estimation of tulathromycin depletion in plasma and milk after subcutaneous injection in lactating goats using a nonlinear mixed-effects pharmacokinetic modeling approach

BACKGROUND

Extra-label use of tulathromycin in lactating goats is common and may cause violative residues in milk. The objective of this study was to develop a nonlinear mixed-effects pharmacokinetic (NLME-PK) model to estimate tulathromycin depletion in plasma and milk of lactating goats. Eight lactating goats received two subcutaneous injections of 2.5 mg/kg tulathromycin 7 days apart; blood and milk samples were analyzed for concentrations of tulathromycin and the common fragment of tulathromycin (i.e., the marker residue CP-60,300), respectively, using liquid chromatography mass spectrometry. Based on these new data and related literature data, a NLME-PK compartmental model with first-order absorption and elimination was used to model plasma concentrations and cumulative excreted amount in milk. Monte Carlo simulations with 100 replicates were performed to predict the time when the upper limit of the 95% confidence interval of milk concentrations was below the tolerance.

RESULTS

All animals were healthy throughout the study with normal appetite and milk production levels, and with mild-moderate injection-site reactions that diminished by the end of the study. The measured data showed that milk concentrations of the marker residue of tulathromycin were below the limit of detection (LOD = 1.8 ng/ml) 39 days after the second injection. A 2-compartment model with milk as an excretory compartment best described tulathromycin plasma and CP-60,300 milk pharmacokinetic data. The model-predicted data correlated with the measured data very well. The NLME-PK model estimated that tulathromycin plasma concentrations were below LOD (1.2 ng/ml) 43 days after a single injection, and 62 days after the second injection with a 95% confidence. These estimated times are much longer than the current meat withdrawal time recommendation of 18 days for tulathromycin in non-lactating cattle.

CONCLUSIONS

The results suggest that twice subcutaneous injections of 2.5 mg/kg tulathromycin are a clinically safe extra-label alternative approach for treating pulmonary infections in lactating goats, but a prolonged withdrawal time of at least 39 days after the second injection should be considered to prevent violative residues in milk and any dairy goat being used for meat should have an extended meat withdrawal time.

Human Food Safety Implications of Variation in Food Animal Drug Metabolism

Violative drug residues in animal-derived foods are a global food safety concern. The use of a fixed main metabolite to parent drug (M/D) ratio determined in healthy animals to establish drug tolerances and withdrawal times in diseased animals results in frequent residue violations in food-producing animals. We created a general physiologically based pharmacokinetic model for representative drugs (ceftiofur, enrofloxacin, flunixin, and sulfamethazine) in cattle and swine based on extensive published literature. Simulation results showed that the M/D ratio was not a fixed value, but a time-dependent range. Disease changed M/D ratios substantially and extended withdrawal times; these effects exhibited drug- and species-specificity. These results challenge the interpretation of violative residues based on the use of the M/D ratio to establish tolerances for metabolized drugs.

Estimation of residue depletion of cyadox and its marker residue in edible tissues of pigs using physiologically based pharmacokinetic modelling

Physiologically based pharmacokinetic (PBPK) models are powerful tools to predict tissue distribution and depletion of veterinary drugs in food animals. However, most models only simulate the pharmacokinetics of the parent drug without considering their metabolites. In this study, a PBPK model was developed to simultaneously describe the depletion in pigs of the food animal antimicrobial agent cyadox (CYA), and its marker residue 1,4-bisdesoxycyadox (BDCYA). The CYA and BDCYA sub-models included blood, liver, kidney, gastrointestinal tract, muscle, fat and other organ compartments. Extent of plasma-protein binding, renal clearance and tissue-plasma partition coefficients of BDCYA were measured experimentally. The model was calibrated with the reported pharmacokinetic and residue depletion data from pigs dosed by oral gavage with CYA for five consecutive days, and then extrapolated to exposure in feed for two months. The model was validated with 14 consecutive day feed administration data. This PBPK model accurately simulated CYA and BDCYA in four edible tissues at 24-120 h after both oral exposure and 2-month feed administration. There was only slight overestimation of CYA in muscle and BDCYA in kidney at earlier time points (6-12 h) when dosed in feed. Monte Carlo analysis revealed excellent agreement between the estimated concentration distributions and observed data. The present model could be used for tissue residue monitoring of CYA and BDCYA in food animals, and provides a foundation for developing PBPK models to predict residue depletion of both parent drugs and their metabolites in food animals.

A physiologically based pharmacokinetic model for polyethylene glycol-coated gold nanoparticles of different sizes in adult mice

Nanoparticles (NPs) are widely used in various fields of nanomedicine. A systematic understanding of NP pharmacokinetics is crucial in their design, applications, and risk assessment. In order to integrate available experimental information and to gain insights into NP pharmacokinetics, a membrane-limited physiologically based pharmacokinetic (PBPK) model for polyethylene glycol-coated gold (Au) NPs (PEG-coated AuNPs) was developed in mice. The model described endocytosis of the NPs in the liver, spleen, kidneys, and lungs and was calibrated using data from mice that were intravenously injected with 0.85 mg/kg 13 nm and 100 nm PEG-coated AuNPs. The model adequately predicted multiple external datasets for PEG-coated AuNPs of similar sizes (13-20 nm; 80-100 nm), indicating reliable predictive capability in suitable size ranges. Simulation results suggest that endocytosis of NPs is time and size dependent, i.e. endocytosis of larger NPs occurs immediately and predominately from the blood, whereas smaller NPs can diffuse through the capillary wall and their endocytosis appears mainly from the tissue with a 10-h delay, which may be the primary mechanism responsible for the reported size-dependent pharmacokinetics of NPs. Several physiological parameters (e.g. liver weight fraction of body weight) were identified to have a high influence on selected key dose metrics, indicating the need for additional interspecies comparison and scaling studies and to conduct pharmacokinetic studies of NPs in species that are more closely related to humans in these parameters. This PBPK model provides useful insights into the size, time, and species dependence of NP pharmacokinetics.

A framework for meta-analysis of veterinary drug pharmacokinetic data using mixed effect modeling

Combining data from available studies is a useful approach to interpret the overwhelming amount of data generated in medical research from multiple studies. Paradoxically, in veterinary medicine, lack of data requires integrating available data to make meaningful population inferences. Nonlinear mixed-effects modeling is a useful tool to apply meta-analysis to diverse pharmacokinetic (PK) studies of veterinary drugs. This review provides a summary of the characteristics of PK data of veterinary drugs and how integration of these data may differ from human PK studies. The limits of meta-analysis include the sophistication of data mining, and generation of misleading results caused by biased or poor quality data. The overriding strength of meta-analysis applied to this field is that robust statistical analysis of the diverse sparse data sets inherent to veterinary medicine applications can be accomplished, thereby allowing population inferences to be made.

Development and application of a multiroute physiologically based pharmacokinetic model for oxytetracycline in dogs and humans

Oxytetracycline (OTC) is a commonly used tetracycline antibiotic in veterinary and human medicine. To establish a quantitative model for predicting OTC plasma and tissue exposure, a permeability-limited multiroute physiologically based pharmacokinetic model was developed in dogs. The model was calibrated with plasma pharmacokinetic data in beagle dogs following single intravenous (5 mg/kg), oral (100 mg/kg), and intramuscular (20 mg/kg) administrations. The model predicted other available dog data well, including drug concentrations in the liver, kidney, and muscle after repeated exposure, and data in the mixed-breed dog. The model was extrapolated to humans and the human model adequately simulated measured plasma OTC concentrations after intravenous (7.14 mg/kg) and oral exposures (6.67 mg/kg). The dog model was applied to predict 24-h OTC area-under-the-curve after three therapeutic treatments. Results were 27.75, 51.76, and 64.17 μg/mL*h in the plasma, and 120.93, 225.64, and 279.67 μg/mL*h in the kidney for oral (100 mg/kg), intravenous (10 mg/kg), and intramuscular (20 mg/kg) administrations, respectively. This model can be used to predict plasma and tissue concentrations to aid in designing optimal therapeutic regimens with OTC in veterinary, and potentially, human medicine; and as a foundation for scaling to other tetracycline antibiotics and to other animal species. © 2014 Wiley Periodicals, Inc. and the American Pharmacists Association J Pharm Sci 104:233-243, 2015.