A physiologically based pharmacokinetic model for the prediction of the depletion of methyl‐3‐quinoxaline‐2‐carboxylic acid,the marker residue of olaquindox,in the edible tissues of pigs

To estimate the consumer exposure to olaquindox (OLA) residues in porcine edible tissues, a physiologically based pharmacokinetic (PBPK) model for methyl‐3‐quinoxaline‐2‐carboxylic acid (MQCA), the marker residue of OLA, was developed in pigs based on the assumptions of the flow‐limited distribution, hepatic metabolism, and renal excretion. The model included separate compartments corresponding to blood, muscle, liver, kidney, adipose, and an extra compartment representing the remaining carcass. Physiological parameters were determined from literatures. Plasma protein binding, partition coefficients, and renal clearance for MQCA were determined in in vitro and in vivo studies. The metabolic conversion of OLA to MQCA was assumed as a simple, one‐step process, and an apparent first‐order rate constant (k) was employed to describe this metabolic process. The PBPK model was optimized and validated with plasma and tissue data from literatures and our study. Sensitivity analysis and Monte Carlo simulation were also implemented to estimate the influence of model parameters on the goodness of fit. When compared with the observed data, the PBPK model underestimated the MQCA level in all compartments at the early time points, whereas gave excellent predictions of MQCA concentration in porcine edible tissues at later time points. The correlation coefficients between the predicted and observed values were over 0.88. The consistency between the model predictions and the real residues of OLA in pigs proved the good applicability of our model in food safety risk assessment.     

Development of a physiologic-based pharmacokinetic model for estimating sulfamethazine concentrations in swine and application to prediction of violative residues in edible tissues

OBJECTIVE: To develop a flow-limited, physiologic-based pharmacokinetic model for use in estimating concentrations of sulfamethazine after IV administration to swine. SAMPLE POPULATION: 4 published studies provided physiologic values for organ weights, blood flows, clearance, and tissue-to-blood partition coefficients, and 3 published studies provided data on plasma and other tissue compartments for model validation. PROCEDURE: For the parent compound, the model included compartments for blood, adipose, muscle, liver, and kidney tissue with an extra compartment representing the remaining carcass. Compartments for the N-acetyl metabolite included the liver and the remaining body. The model was created and optimized by use of computer software. Sensitivity analysis was completed to evaluate the importance of each constant on the whole model. The model was validated and used to estimate a withhold interval after an IV injection at a dose of 50 mg/kg. The withhold interval was compared to the interval estimated by the Food Animal Residue Avoidance Databank (FARAD). RESULTS: Specific tissue correlations for plasma, adipose, muscle, kidney, and liver tissue compartments were 0.93, 0.86, 0.99, 0.94, and 0.98, respectively. The model typically overpredicted concentrations at early time points but had excellent accuracy at later time points. The withhold interval estimated by use of the model was 120 hours, compared with 100 hours estimated by FARAD. CONCLUSIONS AND CLINICAL RELEVANCE: Use of this model enabled accurate prediction of sulfamethazine pharmacokinetics in swine and has applications for food safety and prediction of drug residues in edible tissues.     

Development of a multiroute physiologically based pharmacokinetic model for orbifloxacin in rabbits

To predict the orbifloxacin concentrations in rabbits after multiple routes of administration, a flow‐limited multiroute physiologically based pharmacokinetic (PBPK) model was developed. Three routes of administration (IV, IM, and PO) were incorporated into this model. Physiological parameters including tissue weights and blood flows through different tissues were obtained from the literature. The tissue/plasma partition coefficients (P_Xs) for noneliminating tissues were calculated according to the area method, while the P_Xs for kidney and the rest of the body compartment, together with other parameters for absorption and elimination, were optimized based on the published concentrations. The comparisons between predicted and observed orbifloxacin concentrations proved its validity, and the present model predicted available concentration data well, including those in liver, kidney, muscle, lung, heart, and plasma after oral, intravenous, or intramuscular administration. A local sensitivity analysis was also performed, which showed that the parameters for oral absorption were most influential on the orbifloxacin concentrations. This model was used to predict plasma and tissue concentrations after multiple oral or intramuscular administration. This study demonstrated the feasibility of predicting drug residues in minor species after multiple routes of administration in the extra‐label manner using the PBPK modeling.     

Diffusion‐limited PBPK model for predicting pulmonary pharmacokinetics of florfenicol in pig

For most bacterial lung infections, the concentration of unbound antimicrobial agent in lung interstitial fluid has been thought to be responsible for antimicrobial efficacy. In this study, a diffusion‐limited physiologically based pharmacokinetic (PBPK) model was developed to predict the pulmonary pharmacokinetics of florfenicol (FF) in pigs. The model included separate compartments corresponding to blood, diffusion‐limited lung, flow‐limited muscle, liver, and kidney and an extra compartment representing the remaining carcass. The absorption rate constant and renal and hepatic clearance of FF were determined in vivo. Other parameters were taken from the literature or optimized based on existing pharmacokinetic data. All mathematical operations during the development of the model were performed using acslXtreme version 3.0.2.1 (Aegis Technologies Group, Inc., Huntsville, AL, USA). The model accurately predicted the concentration–time courses of FF in lung interstitial fluid, serum, and plasma following different dosing schedules, except at the dose of 15 mg/kg. When compared with the tissue residue data, the model generally underestimated the FF concentration at the injection site, whereas it gave good predictions of FF concentrations in lung, liver, and kidney at early time points. The model predictions provide a scientific basis for the dosage regimen design of FF.     

A physiologically based pharmacokinetic model linking plasma protein binding interactions with drug disposition

Combination drug therapy increases the chance for an adverse drug reactions due to drug-drug interactions. Altered disposition for sulfamethazine (SMZ) when concurrently administered with flunixin meglumine (FLU) in swine could lead to increased tissue residues. There is a need for a pharmacokinetic modeling technique that can predict the consequences of possible drug interactions. A physiologically based pharmacokinetic model was developed that links plasma protein binding interactions to drug disposition for SMZ and FLU in swine. The model predicted a sustained decrease in total drug and a temporary increase in free drug concentration. An in vivo study confirmed the presence of a drug interaction. Neither the model nor the in vivo study revealed clinically significant changes that alter tissue disposition. This novel linkage approach has use in the prediction of the clinical impact of plasma protein binding interactions. Ultimately it could be used in the design of dosing regimens and in the protection of the food supply through prediction and minimization of tissue residues.     

A physiologically based pharmacokinetic model for oxytetracycline residues in sheep

A physiologically based pharmacokinetic model (PBPK) for oxytetracycline (OTC) residues in sheep was developed using previously published data from a combined serum pharmacokinetic and tissue residue study [Craigmill et al. (2000) J. Vet. Pharmacol. Ther.23, 345]. Physiological parameters for organ weights and tissue blood flows were obtained from the literature. The tissue/serum partition coefficients for OTC were estimated from the serum and tissue residue data obtained at slaughter. The model was developed to include all of the tissues for which residue data were available (serum, kidney, liver, fat, muscle and injection site), and all of the remaining tissues were combined into a slowly perfused compartment with low permeability. Total body clearance of OTC calculated in the previous study was used as the starting value for clearance in the PBPK model, with the kidney being the only eliminating organ. The model was built using ACSL (Advanced Continuous Simulation Language) Graphic Modeler, and the model was fit to the serum and tissue data using the ACSL Math/Optimizer software (AEgis Technologies Group, Inc., Huntsville, AL, USA). A sensitivity analysis was also performed to determine which parameters had the greatest effect on the goodness of fit. Numerous strategies were tested to model the injection site, and a model providing a biexponential absorption of the drug from the injection bolus gave the best fit to the experimental data. The model was validated using the clearance parameters calculated from the traditional pharmacokinetic model for each individual animal in the PBPK model. This simple PBPK model well predicted OTC residues in sheep tissues after intramuscular dosing with a long-acting preparation and may find use for other species and other veterinary drugs.     

Nonlinear mixed‐effects pharmacokinetic modeling of the novel COX‐2 selective inhibitor vitacoxib in dogs

The objective of this study was to develop a nonlinear mixed‐effects model of vitacoxib disposition kinetics in dogs after intravenous (I.V.), oral (P.O.), and subcutaneous (S.C.) dosing. Data were pooled from four consecutive pharmacokinetic studies in which vitacoxib was administered in various dosing regimens to 14 healthy beagle dogs. Plasma concentration versus time data were fitted simultaneously using the stochastic approximation expectation maximization (SAEM) algorithm for nonlinear mixed‐effects as implemented in Monolix version 2018R2. Correlations between random effects and significance of covariates on population parameter estimates were evaluated using multiple samples from the posterior distribution of the random effects. A two‐compartment mamillary model with first‐order elimination and first‐order absorption after P.O. and S.C. administration, best described the available pharmacokinetic data. Final parameter estimates indicate that vitacoxib has a low‐to‐moderate systemic clearance (0.35 L hr^?1 kg^?1) associated with a low global extraction ratio, but a large volume of distribution (6.43 L/kg). The absolute bioavailability after P.O. and S.C. administration was estimated at 10.5% (fasted) and 54.6%, respectively. Food intake was found to increase vitacoxib oral bioavailability by a fivefold, while bodyweight (BW) had a significant impact on systemic clearance, thereby confirming the need for BW adjustment with vitacoxib dosing in dogs.     

Factors that affect drug disposition in food-producing animals during maturation.

Drugs administered to neonatal food-producing animals (cattle, sheep, goats, swine) may exhibit significantly different pharmacokinetic/disposition characteristics than they do in adult animals of the same species. Undesirable consequences such as suboptimum therapeutic concentrations, toxic effects, and violative tissue residues may result if adult dosage regimens are employed in young animals. Using selected drugs as examples, this paper reviews factors that contribute to differences in drug disposition in newborn vs adult animals. Immaturity of mechanisms involved in drug absorption, especially from gastrointestinal and parenteral sites of administration, and of drug distribution to sites such as plasma proteins, adipose tissue, and fluid compartments are considered. The role of developmental changes in drug biotransformation in the liver and other tissues and the maturation of excretory mechanisms, primarily from the kidney, in the increased rate of drug clearance during maturation is described. Pharmacokinetic studies with specific drugs in the target species are an important approach to establishing rational drug use in immature food-producing animals.     

The pharmacokinetics of glycopyrrolate in Standardbred horses

The disposition of plasma glycopyrrolate (GLY) is characterized by a three‐compartment pharmacokinetic model after a 1‐mg bolus intravenous dose to Standardbred horses. The median (range) plasma clearance (Clp), volume of distribution of the central compartment (V₁), volume of distribution at steady‐state (Vss), and area under the plasma concentration–time curve (AUC_0‐inf) were 16.7 (13.6–21.7) mL/min/kg, 0.167 (0.103–0.215) L/kg, 3.69 (0.640–38.73) L/kg, and 2.58 (2.28–2.88) ng*h/mL, respectively. Renal clearance of GLY was characterized by a median (range) of 2.65 (1.92–3.59) mL/min/kg and represented approximately 11.3–24.7% of the total plasma clearance. As a result of these studies, we conclude that the majority of GLY is cleared through hepatic mechanisms because of the limited extent of renal clearance of GLY and absence of plasma esterase activity on GLY metabolism. Although the disposition of GLY after intravenous administration to Standardbred horses was similar to that in Thoroughbred horses, differences in some pharmacokinetic parameter estimates were evident. Such differences could be attributed to breed differences or study conditions. The research could provide valuable data to support regulatory guidelines for GLY in Standardbred horses.     

Carprofen pharmacokinetics in plasma and in control and inflamed canine tissue fluid using in vivo ultrafiltration

Measurement of unbound drug concentrations at their sites of action is necessary for accurate PK/PD modeling. The objective of this study was to determine the unbound concentration of carprofen in canine interstitial fluid (ISF) using in vivo ultrafiltration and to compare pharmacokinetic parameters of free carprofen concentrations between inflamed and control tissue sites. We hypothesized that active concentrations of carprofen would exhibit different dispositions in ISF between inflamed vs. normal tissues. Bilateral ultrafiltration probes were placed subcutaneously in six healthy Beagle dogs 12 h prior to induction of inflammation. Two milliliters of either 2% carrageenan or saline control was injected subcutaneously at each probe site, 12 h prior to intravenous carprofen (4 mg/kg) administration. Plasma and ISF samples were collected at regular intervals for 72 h, and carprofen concentrations were determined using HPLC. Prostaglandin E2 (PGE₂) concentrations were quantified in ISF using ELISA. Unbound carprofen concentrations were higher in ISF compared with predicted unbound plasma drug concentrations. Concentrations were not significantly higher in inflamed ISF compared with control ISF. Compartmental modeling was used to generate pharmacokinetic parameter estimates, which were not significantly different between sites. Terminal half‐life (T½) was longer in the ISF compared with plasma. PGE₂ in ISF decreased following administration of carprofen. In vivo ultrafiltration is a reliable method to determine unbound carprofen in ISF, and that disposition of unbound drug into tissue is much higher than predicted from unbound drug concentration in plasma. However, concentrations and pharmacokinetic parameter estimates are not significantly different in inflamed vs. un‐inflamed tissues.     

Pharmacokinetics and bioavailability after intramuscular injection of the 5‐HT1A serotonin agonist R‐8‐hydroxy‐2‐(di‐n‐propylamino) tetralin (8‐OH‐DPAT) in domestic goats (Capra aegagrus hircus)

To determine the bioavailability and pharmacokinetic properties of the serotonin 5‐HT_1A receptor agonist R‐8‐OH‐DPAT in goats, and 0.1 mg kg^?1 R‐8‐OH‐DPAT hydrobromide was administered intramuscularly (i.m.) and intravenously (i.v.) to six goats in a two‐phase cross‐over design experiment. Venous blood samples were collected from the jugular vein 2, 5, 10, 15, 20, 30, 40 and 60 min following treatment and analysed by liquid chromatography tandem mass spectrometry. Bioavailability and pharmacokinetic parameters were determined by a one‐compartment analysis. Mean bioavailability of R‐8‐OH‐DPAT when injected i.m. was 66%. The mean volume of distribution in the central compartment was 1.47 L kg^?1. The mean plasma body clearance was 0.056 L kg^?1 min^?1. All goats injected i.v. and two of six goats injected i.m. showed signs of serotonin toxicity. In conclusion, R‐8‐OH‐DPAT is well absorbed following i.m. injection and the observed pharmacokinetics suggest that administration via dart is feasible. Administration of R‐8‐OH‐DPAT hydrobromide, at a dosage of 0.1 mg kg^?1, resulted in the observation of clinical signs of serotonin toxicity in the goats. It is suggested that dosages for the clinical use of the compound should be lower in order to achieve the desired clinical effect without causing serotonin toxicity.     

A physiologically based pharmacokinetic model for valnemulin in rats and extrapolation to pigs

A flow-limited, physiologically based pharmacokinetic (PBPK) model for predicting the plasma and tissue concentrations of valnemulin after a single oral administration to rats was developed, and then the data were extrapolated to pigs so as to predict withdrawal interval in edible tissues. Blood/tissue pharmacokinetic data and blood/tissue partition coefficients for valnemulin in rats and pigs were collected experimentally. Absorption, distribution and elimination of the drug were characterized by a set of mass-balance equations. Model simulations were achieved using a commercially available software program. The rat PBPK model better predicted plasma and tissue concentrations. The correlation coefficients of the predicted and experimentally determined values for plasma, liver, kidney, lung and muscle were 0.96, 0.94, 0.96, 0.91 and 0.91, respectively. The rat model parameters were extrapolated to pigs to estimate valnemulin residue withdrawal interval in edible tissues. Correlation (R(2) ) between predicted and observed liver, kidney and muscle were 0.95, 0.97 and 0.99, respectively. Based on liver tissue residue profiles, the pig model estimated a withdrawal interval of 10 h under a multiple oral dosing schedule (5.0 mg/kg, twice daily for 7.5 days). PBPK models, such as this one, provide evidence of the usefulness in interspecies PK data extrapolation over a range of dosing scenarios and can be used to predict withdrawal interval in pigs.     

Mixed effects modeling of the disposition of gentamicin across domestic animal species

An interspecies pharmacokinetic model for gentamicin was developed using the mixed effects modeling approach and serum disposition data obtained from the Food Animal Residue Avoidance Databank (FARAD). Data that met a priori quality criteria was obtained from the database and analysed using the traditional double logarithmic analysis and the mixed effects modeling approach. Body weight, brain weight and fever were the covariates of interest in our study. Population pharmacokinetic models across species were developed and validated with swine data. The parameter volume of distribution was modeled as a function of body weight. The total clearance was initially modeled as a function of body weight. The predictability performance of the model improved dramatically when the parameter brain weight was included in the covariate model for clearance. This was a surprising finding worthy of further study. The covariate fever seemed to influence the magnitude of the volume of distribution, although the scarcity of data pertaining to diseased animals makes this finding uncertain. We conclude that the pharmacokinetic characteristics of drugs such as gentamicin, can be predicted across species using a population pharmacokinetics modeling approach, and that clinical features that affect species in a similar manner can be also explored in this fashion.     

Measurement of carbonic anhydrase isozyme VI (CA‐VI) in swine sera,colostrums, saliva,bile, seminal plasma and tissues

Swine secretory carbonic anhydrase VI (CA‐VI) was purified from swine saliva and an antibody to CA‐VI was generated. A specific and sensitive enzyme‐linked immunosorbent assay (ELISA) has been developed for the measurement of swine CA‐VI. The assay can detect as little as 5 ng/mL of swine CA‐VI. Typical standard curves were determined for a range of CA‐VI solutions (7.8 to 500 ng/mL). The coefficients of variation for these solutions were less than 5%. When 500, 250 or 100 ng/mL of swine CA‐VI was added to swine sera, the recoveries were 102.0%, 109.7% and 100.2%, respectively. The concentrations of CA‐VI in the saliva (26.2 ± 30.4 µg/mL), sera (3.3 ± 4.9 ng/mL), bile (153.0 ± 114.0 ng/mL), seminal plasma (124.0 ± 39.0 ng/mL) and parotid gland (441.3 ± 90.0 µg/g wet tissue), submaxillary gland (88.1 ± 124.4 µg/g wet tissue), sublingual gland (58.6 ± 24.6 µg/g wet tissue) and gallbladder (2.4 ± 1.3 µg/1g wet tissue) were determined by ELISA. The concentration of CA‐VI in colostrum was 163.3 ± 101.4 ng/mL and did not decrease within 10 days following parturition. An immunohistochemical reaction to anti‐CA‐VI antiserum was observed in the columnar epithelial cells lining the gallbladder. These data suggest that secretory CA‐VI plays various roles in pH regulation and the maintenance of ion and fluid balance.     

Tolerability and pharmacokinetic profile of a novel benzene‐poly‐carboxylic acids complex with cis‐diammineplatinum (II) dichloride in dogs with malignant mammary tumours

The pharmacokinetic profile, tolerability and efficacy of benzene‐poly‐carboxylic acids complex with cis‐diammineplatinum (II) dichloride (BP‐C1) were studied in dogs with mammary cancer. A three‐level response surface pathway designed trial was performed on seven dogs. At each level BP‐C1 was administered subcutaneously daily for 7 days followed by a 7‐day rest period in a dose escalating manner. Adverse events according to VCOG‐CTCAE, performance status and tumour progression were recorded. The pharmacokinetic profile followed a two‐compartment model with rapid absorption, short distribution, and a slow elimination phase. The overall elimination half‐life was 125 h. The maximum tolerated dose of BP‐C1 was estimated to be above 0.46 mg kg^?1. A significant reduction in VCOG‐CTCAE toxicity which correlated negatively with increasing dose was found. The dogs' general performance status remained unchanged. No decrease in total tumour burden was found, although temporary tumour reduction was seen in some target tumours.     

A PBPK model for midazolam in four avian species

A physiologically based pharmacokinetic (PBPK) model was developed for midazolam in the chicken and extended to three other species. Physiological parameters included organ weights obtained from 10 birds of each species and blood flows obtained from the literature. Partition coefficients for midazolam in tissues vs. plasma were estimated from drug residue data obtained at slaughter. The avian models include separate compartments for venous plasma, liver, kidney, muscle, fat and all other tissues. An estimate of total body clearance from an earlier in vitro study was used as a starting value in the model, assuming almost complete removal of the parent compound by liver metabolism. The model was optimized for the chicken with plasma and tissue data from a pharmacokinetic study after intravenous midazolam (5 mg/kg) dose. To determine which parameters had the most influence on the goodness of fit, a sensitivity analysis was performed. The optimized chicken model was then modified for the turkey, pheasant and quail. The models were validated with midazolam plasma and tissue residue data in the turkey, pheasant and quail. The PBPK models in the turkey, pheasant and quail provided good predictions of the observed tissue residues in each species, in particular for liver and kidney.     

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

Yang, F., Sun, N., Sun, Y. X., Shan, Q., Zhao, H. Y., Zeng, D. P., Zeng, Z. L. A physiologically based pharmacokinetics model for florfenicol in crucian carp and oral‐to‐intramuscular extrapolation. J. vet. Pharmacol. Therap.  36 , 192–200. In this study, an oral physiologically based pharmacokinetics (PBPK) model was developed for florfenicol in crucian carp (Carassius auratus). Subsequently, oral‐to‐intramuscular extrapolation was performed and the two models were used to predict florfenicol concentrations in the edible tissues of crucian carp. The oral model gave good predictions in most tissues, except for kidney and liver in which the florfenicol concentrations were underestimated at the later time points. In contrast, using the intramuscular model, the concentrations in the kidney were overestimated at the later time points. Both models had the best predictive ability in the main edible tissue, the muscle. The oral model also accurately predicted the florfenicol concentrations in the muscle after multiple doses. The present study demonstrated the feasibility of predicting florfenicol concentrations in the edible tissues of crucian carp using a route‐to‐route extrapolation method.     

Evaluation of enrofloxacin use in koalas (Phascolarctos cinereus) via population pharmacokinetics and Monte Carlo simulation

Clinically normal koalas (n = 6) received a single dose of intravenous enrofloxacin (10 mg/kg). Serial plasma samples were collected over 24 h, and enrofloxacin concentrations were determined via high‐performance liquid chromatography. Population pharmacokinetic modeling was performed in S‐ADAPT. The probability of target attainment (PTA) was predicted via Monte Carlo simulations (MCS) using relevant target values (30–300) based on the unbound area under the curve over 24 h divided by the minimum inhibitory concentration (MIC) (fAUC_0–24/MIC), and published subcutaneous data were incorporated (Griffith et al., 2010). A two‐compartment disposition model with allometrically scaled clearances (exponent: 0.75) and volumes of distribution (exponent: 1.0) adequately described the disposition of enrofloxacin. For 5.4 kg koalas (average weight), point estimates for total clearance (SE%) were 2.58 L/h (15%), central volume of distribution 0.249 L (14%), and peripheral volume 2.77 L (20%). MCS using a target fAUC_0–24/MIC of 40 predicted highest treatable MICs of 0.0625 mg/L for intravenous dosing and 0.0313 mg/L for subcutaneous dosing of 10 mg/kg enrofloxacin every 24 h. Thus, the frequently used dosage of 10 mg/kg enrofloxacin every 24 h subcutaneously may be appropriate against gram‐positive bacteria with MICs ≤ 0.03 mg/L (PTA > 90%), but appears inadequate against gram‐negative bacteria and Chlamydiae in koalas.     

Understanding the pharmacokinetics of tulathromycin: a pulmonary perspective

Tulathromycin is approved in the United States for the treatment of respiratory disease in bovine and swine, infectious bovine keratoconjunctivitis associated with Moraxella bovis, and interdigital necrobacillosis in bovine. This macrolide highly concentrates in lung tissue and persists in the intra‐airway compartment for a long time after a single administration. It also accumulates in inflammatory cells, including neutrophils and macrophages. This article reviews pharmacokinetic information about tulathromycin in different veterinary species with particular emphasis on the respiratory system.     

A possible solution to model nonlinearity in elimination and distributional clearances with α2‐adrenergic receptor agonists: Example of the intravenous detomidine and methadone combination in sedated horses

The alpha(α)₂‐agonist detomidine is used for equine sedation with opioids such as methadone. We retrieved the data from two randomized, crossover studies where detomidine and methadone were given intravenously alone or combined as boli (STUDY 1) (Gozalo‐Marcilla et al., 2017, Veterinary Anaesthesia and Analgesia, 2017, 44 , 1116) or as 2‐hr constant rate infusions (STUDY 2) (Gozalo‐Marcilla et al., 2019, Equine Veterinary Journal, 51 , 530). Plasma drug concentrations were measured with a validated tandem Mass Spectrometry assay. We used nonlinear mixed effect modelling and took pharmacokinetic (PK) data from both studies to fit simultaneously both drugs and explore their nonlinear kinetics. Two significant improvements over the classical mammillary two‐compartment model were identified. First, the inclusion of an effect of detomidine plasma concentration on the elimination clearances (Cls) of both drugs improved the fit of detomidine (Objective Function Value [OFV]: ?160) and methadone (OFV: ?132) submodels. Second, a detomidine concentration‐dependent reduction of distributional Cls of each drug further improved detomidine (OFV: ?60) and methadone (OFV: ?52) submodel fits. Using the PK data from both studies (a) helped exploring hypotheses on the nonlinearity of the elimination and distributional Cls and (b) allowed inclusion of dynamic effects of detomidine plasma concentration in the model which are compatible with the pharmacology of detomidine (vasoconstriction and reduction in cardiac output).