VetCAST Method for Determination of the Pharmacokinetic-Pharmacodynamic Cut-Off Values of a Long-Acting Formulation of Florfenicol to Support Clinical Breakpoints for Florfenicol Antimicrobial Susceptibility Testing in Cattle

The PK/PD cut-off (PK/PD_CO) value of florfenicol for calf pathogens was determined for long acting formulations (MSD Nuflor^® and a bioequivalent generic product). PK/PD_CO is one of the three MICs considered by VetCAST, a sub-committee of the European Committee on Susceptibility Testing (EUCAST), to establish a Clinical Breakpoint for interpreting Antimicrobial Susceptibility Testing (AST). A population model was built by pooling three pharmacokinetic data sets, obtained from 50 richly sampled calves, receiving one of two formulations (the pioneer product and a generic formulation). A virtual population of 5,000 florfenicol disposition curves was generated by Monte Carlo Simulations (MCS) over the 96 h of the assumed duration of action of the formulations. From this population, the maximum predicted MIC, for which 90% of calves can achieve some a priori selected critical value for two PK/PD indices, AUC/MIC and T>MIC, was established. Numerical values were established for two bacterial species of the bovine respiratory disease (BRD) complex, Pasteurella multocida and Mannheimia haemolytica. It was concluded that the PK/PD_CO of florfenicol for both AUC/MIC and T>MIC was 1 mg/L.

Introduction

Florfenicol is an antimicrobial drug (AMD) used extensively to treat Bovine Respiratory Disease (BRD). It's prudent and rational use should be based on the results of Antimicrobial Susceptibility Testing (AST). Clinical breakpoints (CBP) are the MIC values (units mg/L) used by antimicrobial testing laboratories to report qualitatively the results of AST as Susceptible or not. The Veterinary Committee on Antimicrobial Susceptibility Testing (VetCAST) is a recently established sub-committee of the European Union Committee on Susceptibility Testing (EUCAST). EUCAST is the reference committee for AST in human medicine for the EU and VetCAST operates within the guidelines and structure of EUCAST. The VetCAST remit encompasses all aspects of AST of bacterial pathogens of animal origin and animal bacteria with zoonotic potential. In the VetCAST approach (Toutain et al., 2017), CBPs are determined by taking into account at least an epidemiological cut-off (ECOFF) and a PK/PD cut-off (PK/PD_CO). In addition, a clinical cut-off can also be considered when clinical data are available to link MICs to clinical efficacy (Turnidge and Martinez, 2017). PK/PD_CO is defined as the highest possible MIC for which a given percentage of animals in the target population (e.g., 90%) achieve a pre-defined target value, hereafter named PDT (pharmacodynamic target) according to European Medicines Agency (EMA) terminology (European Medicines Agency, 2015). For Histophilus somni (HS), Pasteurella multocida (PM), and Mannheimia haemolytica (MH), possible florfenicol MICs for the wild populations ranged from 0.12 to 2 mg/L; MIC90 values were 0.25 mg/L (HS), 0.5 mg/L (PM), and 1 mg/L (MH) (de Jong et al., 2014).

The Veterinary Antimicrobial Susceptibility Testing (VAST) sub-committee of the Clinical and Laboratory Standards Institute (CLSI), hereafter named CLSI/VAST, historically approved CBP for florfenicol for BRD treatment; selected values were 2, 4, and 8 mg/L, respectively, for Susceptible (S), Intermediate (I), or Resistant (R) (Clinical Laboratory Standards Institute, 2018). To our knowledge, these CBPs were not accompanied by a CLSI/VAST explanatory document to justify the selected values for BRD, although this is now the case for all new CLSI/VAST CBPs. At the time of ascribing these values to florfenicol, CLSI/VAST did not consider PK/PD relationships in the decision taking methods for establishing the CBP. In veterinary medicine, publicly available clinical data on AMD efficacy are generally scarce or non-existent. For florfenicol, several publications have described the clinical efficacy of the formulations (MSD Nuflor^® and its generics) considered in the present paper. Results of clinical trials for florfenicol were comprehensively analyzed using a mixed-treatment comparison meta-analysis, which combined evidence from published trials and published estimates of comparative efficacy for 12 AMDs registered for use in the USA (O'Connor et al., 2016). It was concluded that florfenicol was efficacious, ranking fourth of the 12 AMDs investigated. VetCAST, having no access to the company files describing the results of these clinical trials, considers, as do others, that the pivotal information required to establish a CBP is embedded in a PK/PD breakpoint (Turnidge and Paterson, 2007). This is because the PK/PD breakpoint is a hybrid value, incorporating all three principal components (microbiological, pharmacological, and clinical) predicting clinical efficacy. Hence, EUCAST relies on such PK/PD breakpoints to establish CBPs. PK/PD breakpoints should be clearly distinguished from a PK/PD_CO, in that the latter is derived from PK data only, without any clinical data input. PK/PD_CO is established solely by exploring a range of possible (not probable) MICs, and the VetCAST methodology involves computing a series of Probability of Target Attainments (PTA) from plasma concentration-time profiles. This is also the procedure adopted by CLSI/VAST under the name PD_CO (Clinical Laboratory Standards Institute, 2018).

VetCAST has re-evaluated the CBP for florfenicol in cattle, in order to provide a proof of concept of its scientific approach, which may differ significantly from that of VAST/CLSI in several respects, including determination of a PK/PD_CO. Being pivotal for the VetCAST approach, a robust estimation of PK/PD_CO requires first the building of a valid population pharmacokinetic (POP PK) model from individual animal data collected from differing sources to quantify typical PK parameters and their between-subject variability (BSV). Simply retrieving, from literature publications, PK parameters estimated by others and aggregating them is not used by VetCAST. Florfenicol and calf pathogens have been selected to illustrate the VetCAST method of meta-analysis (Li et al., 2015). The Non-Linear Mixed Effect Model (NLME) is used to handle unbalanced data (Schoemaker and Cohen, 1996) with one data set having been analyzed using a mono-compartmental model (Sidhu et al., 2014), while more recent data sets have been obtained with a lower limit of quantification (LLOQ) of the analytical technique, thereby providing an extended terminal half-life. This is a very common situation in veterinary medicine, as long-acting (LA) formulations are used extensively (Toutain and Bousquet-Mélou, 2004).

A further aspect of data analysis, specific to VetCAST, is the rationale for selecting the most appropriate PK/PD index, either the time for which plasma concentration remains above the MIC during the dosage interval (f T>MIC) or the ratio of Area Under the plasma concentration-time Curve divided by the MIC (f AUC/MIC), f Cmax/MIC being not considered by EUCAST. The term f indicates that these indices should be computed in terms of plasma free drug concentrations. For florfenicol in cattle, the binding to plasma proteins has been reported in several publications, with very disparate results. At the time of model building, the most recent protein binding data for florfenicol were those published by Foster et al. (2016) who concluded “Florfenicol protein binding was only 5% at the high concentration and was negligible at the low concentrations, representing a fu of essentially 1.0.” However, others have reported values ranging from 10 to ~25% (Adams et al., 1987; Bretzlaff et al., 1987; Lobell et al., 1994; Sidhu et al., 2014). In light of these data heterogeneity, it was decided to ignore the extent of drug binding in making the present computations, as further explained the Discussion.

Florfenicol is often classified as time-dependent in its killing action and, as for chloramphenicol, T>MIC has been reported as the appropriate PK/PD index (Martinez et al., 2013). However, AUC/MIC has also been proposed as the most appropriate index predictive of clinical efficacy, especially for PM and MH (Sidhu et al., 2014). Actually, it has been shown, using a semi-mechanistic in silico model, that AUC/MIC (and not T>MIC) is the most appropriate index, when terminal half-life is relatively long relative to the dosing interval, even for beta-lactam drugs (Nielsen et al., 2011; Kristoffersson et al., 2016). This is the case likewise for florfenicol LA formulations. In the VetCAST project, the best predictive index for florfenicol and its magnitude were investigated from in silico simulations using a semi-mechanistic PK/PD model (Nielsen and Friberg, 2013) to replace the classical in vivo rodent infection model that, for several decades (Craig, 1998; Andes and Craig, 2002), was used to select the best PK/PD index. VetCAST calculates PDT through an in silico dose-fractionation approach (Pelligand et al., 2019).

The aim of the present investigation was to build a population model for florfenicol in cattle, generating by Monte Carlo simulations (MCS) a large number of plasma florfenicol disposition curves (n = 5,000). This virtual in silico meta-population was used to determine the percentages of animals (PTA) for which a series of possible PDT values would be attainable with differing possible MICs (actually 0.25, 0.5, 1, and 2 mg/L).

Materials and Methods

Individual calf PK data from three different sources (A = 10, B = 32, C = 8) were used for the POP PK analysis (see Supplementary Material, giving raw data). Source A consisted of 10 calves from a published study (Sidhu et al., 2014). Source B was a drug company (Norbrook Laboratories Limited); it comprised 16 calves enrolled in a cross-over bioequivalence study (MSD Nuflor^® and Norbrook Norfenicol^® formulations) Norfenicol^® being a FDA and EMA approved generic product (Anonymous, 2018a). The 32 data sets were provided by 16 sets for each product, so that for this analysis each of these calves provided two data sets. The third source comprised data from 8 calves in an unpublished study. All calves were in good health and all received a subcutaneous florfenicol dose of 40 mg/kg. Table 1 gives details for the three sources of individual animal data.

TABLE 1

Table 1. The three data sets considered for florfenicol population pharmacokinetic analysis.

Pharmacokinetic data analyses were carried out using Phoenix^® WinNonlin^® 8.0 (Pharsight Corporation St. Louis, MO, USA). Data sets obtained from the three sources were analyzed using a NLME model. A two-compartmental model was selected, based on the Likelihood Ratio Test (LRT), the Akaike Information Criterion (AIC), and inspection of different diagnostic plots (vide infra). For the LRT test, the critical value of the χ² distribution considered for a given nominal risk of 0.05, and a given number of degrees of freedom, was obtained using the Excel function CHISQ.INV.RT().

The parameterization of the structural two-compartmental model was of the closed form (Equation 1):

$\begin{array}{l}C\left(t\right)=A\times EXP\left(-Alpha\times t\right)+B\times EXP\left(-Beta\times t\right)\\ -(A+B)\times EXP(-Ka\times t)& (1)\end{array}$

where t is the time (h) macroconstants, A and B (μg/ml) are intercepts and Alpha, Beta, and Ka are rate constants (1/h) associated with the phases of plasma concentration-time profile. Parameterization was in terms of macroconstants and rate constants rather than in terms of clearance and volume of distribution for reasons explained in the Discussion. The five fixed parameters (described as vector Thetas) were estimated and reported as typical values (tv) with coefficient of variation as a measure of precision of the estimate. The random component that describes biological variability around the structural fixed parameters i.e., the BSV across individuals was described by an exponential model of the form (Equation 2):

$\begin{array}{ll}{\theta }_{1i}={\theta }_1\times Exp({\eta }_{1i})& (2)\end{array}$

where θ₁ is the typical population value of theta (A, B, Alpha, Beta, or Ka), θ_1i the value of theta in the ith animal, and η_1i (eta) the deviation associated with the ith animal from the corresponding theta population value. This exponential model assumes a log-normal distribution of parameters, i.e., that the distribution of the etas is normal in the log-domain, with a mean of 0 and a variance ω² where:

$\eta \approx N(0,{\omega }^2)$

Each eta distribution associated to each theta with its own variance ${ \omega }_A^2$, ${\omega }_{Alpha}^2 {\omega }_{B }^2$, ${\omega }_{Beta, }^2$or ${\omega }_{Ka}^2$was computed, but covariance terms between etas have been ignored (diagonal matrix) to ensure identifiability of the parameters.

The BSV was reported as coefficient of variation in the original scale with the following equation that converts the variance terms (ω²) to a coefficient of variation (CV%).

$\begin{array}{ll}CV\left(\%\right)=100\times \sqrt{\exp \left({\omega }^2\right)-1}& (3)\end{array}$

The residual variability was modeled with an additive and a multiplicative component. Like other random-effects, the residual error can be dependent on subject-specific covariate of the analytical technique used to generate plasma concentration (Bonate, 2011). Assuming that the residual mainly reflects variability of the analytical technique, we explored, as a part of the quality control of the merged data sets, what might be the precision of the three analytical techniques used to generate the data i.e., included in the error model was the source of the data as a covariate. It was concluded that differences were not sufficiently large to retain this covariate. Therefore, in the final model, a single residual without covariate was used.

The residual error model without covariate was of the form (Equation 4):

$\begin{array}{ll}Y=f\left(\theta ,Time\right)\times \left(1+{\epsilon }_1\right)+{\epsilon }_2\text{'}& (4)\end{array}$

with ε1 the multiplicative error term having a mean of 0 and a variance of σ1

$\epsilon 1\approx N(0,\sigma 1^2)$

and ε2 the common additive error term having a mean of 0 and a variance noted σ2

$\epsilon 2\approx N(0,\sigma 2^2)n$

Sigma1 and Sigma2 were estimated by Phoenix and reported as a CV% for sigma1 and as a STDV for sigma2.

No covariates were included in the final model, as the computed PK/PD_CO is expected to cover all sources of biological variability across animals. However, in a preliminary analysis, two covariates were explored, in order to support the merging of the three data sets (A, B, and C) and the two formulations (Nuflor^® and generic). There was no major influence of these covariates (results not shown) and no specific issue linked to the merging of the data sets.

Parameter estimations, with their associated SE and coefficient of variation as a measure of the precision of the estimate, were based on minimizing an objective function value (OFV), using Laplace engine for the Maximum Likelihood Estimation.

As only 22 florfenicol concentrations were reported as BLQ (comprising 2.6% of the whole data set), BLQ data were discounted in the analysis without the risk of introducing bias in the parameter estimates leading to model mis-specification (Byon et al., 2008). For the two-compartment model, when the BLQ incidence was <5%, it was shown that omission of the BLQ data generally did not inflate the bias in the fixed-effect parameters (Xu et al., 2011).

The shrinkage for the etas was estimated by the equation (Karlsson and Savic, 2007):

$\begin{array}{ll}Etashrinkage=1-\frac{SD(EB{E}_{\eta })}{\omega }& (5)\end{array}$

where ω is the estimated variability for the population and SD is the SD of the individual values of the Empirical Bayesian Estimates (EBE) of η.

Different diagnostic plots were reviewed to determine whether or not a model was adequate. These included PRED (Population Predicted Value based on population parameter estimates) and IPRED (Individual Predicted value based on individual's ETAs) vs. the DV (Dependent variable) (with and without a log scale) Conditional weighted residuals (CWRES) and individual fitting. The overall adequacy of the 2-compartment PK model was established by plotting the Visual Predictive Check (VPC) i.e., a graphical comparison between the observed data and prediction intervals derived from the simulated data.

Secondary parameters were also computed (terminal half-lives for the first and second phases of drug disposition and contribution of the first and second phases to drug absorption).

Monte Carlo simulations (MCS) of the predicted concentration (IPRED) from the model i.e., simulation of concentration without the error term from 0 to 96 h post administration, with a step of 1 h, were used to generate a meta-population of 5,000 calves. These curves were analyzed using the Non-Compartmental tool of Phoenix to compute the areas under the curve and the time above selected MICs from 0 to 96 h, 96 h being the claimed duration of florfenicol activity after a single SC administration of Nuflor^® (Anonymous, 2018b). A PDT of 40% was selected for T>MIC as a default value (Mouton et al., 2012). These metrics were then analyzed with the statistical tool of Phoenix to compute the quantiles of interest (90th) to establish PK/PD_COs.

In human medicine, PK/PD indices and their PDT are established primarily in rodent models over a fractionated-dosing interval of 24 h. For florfenicol, such data are not available. Therefore, in this project, an in silico approach was used as a surrogate for the dose fractionation trial. Briefly, PD parameters for florfenicol were first estimated by modeling killing curves obtained with PM and MH with a semi-mechanistic model described by others (Nielsen and Friberg, 2013). Then, the selected PD model was solved with average plasma concentrations predicted by the population model of the present investigation. This in silico approach established, retrospectively for the main human AMD classes, all indices derived using the animal model (Nielsen et al., 2011). It was concluded that the best index for florfenicol was AUC/MIC. This component of the project is fully described in a companion paper (Pelligand et al., 2019).

Results

Figure 1 displays the 50 curves used in the POP PK analysis, sorted either by sources (n = 3) or by formulations (n = 2). Figures 2–5 are Goodness-of-fit (GOF) plots supporting the 2-comparmental structural model; the exponential model for the random component; and the additive plus multiplicative model for the error sub-model used to analyze the data. To evaluate the adequacy of the developed population model, the VPC plots are presented in Figure 5, which illustrates the 10th, 50th, and 90th percentiles of the simulated distribution compared to the observed values. Typical values of the primary structural parameters of the model (thetas), the secondary parameters (half-life and percentage of the bioavailable dose absorbed during Alpha and Beta phases), their associated Standard Error (SE) and the SD of the residual for the basic model are presented in Table 2.

FIGURE 1

Figure 1. Semi-logarithmic spaghetti plot for 50 calves sorted by sources (left) (Red = A, Gray = B, Blue = C) or by formulation (right) (Black = Nuflor^®, Red = generic). Visual inspection of the plots does not suggest major differences either between the three sources of data or for the two formulations, as seen from the intermingling of the curves.

FIGURE 2

Figure 2. Plot of the dependent variable (DV) i.e., plasma florfenicol concentration (μg/ml) vs. population predicted plasma florfenicol concentrations (PRED) (no random component). The plot illustrates observed vs. fitted values of the model function. Ideally, they should fall close to the line of unity y = x. Arithmetic scale (left) and logarithmic scale (right). For both arithmetic and logarithmic scales, data were evenly distributed about the line of identity, indicating no major bias in the population component of the model.

FIGURE 3

Figure 3. Plot of the dependent variable (DV) i.e., observed plasma florfenicol concentration (μg/ml) vs. individual predicted plasma florfenicol values (IPRED). Individual prediction was obtained by setting random effects to the “post-hoc” or empirical Bayesian estimate of the random effects for the individual from which the DV observation was made. Thus, the plots illustrate observed vs. fitted values of the model function. Ideally they should fall close to the line of unity y = x. Arithmetic scale (left) and logarithmic scale (right). For both the arithmetic and logarithmic scales, data were evenly distributed about the line of identity, indicating no major bias in the population component of the model.

FIGURE 4

Figure 4. Plot of CWRES (conditional weighted residuals) against time after dose (h). Values of CWRES should be approximately N (0, 1) and hence concentrated between y = −2 and y = +2. Inspection of the figure shows that data were evenly distributed about zero (see the trends as given by the blue line) and the red line (with its negative reflection) did not show any fanning, indicating no bias in the structural model.

FIGURE 5

Figure 5. Visual Predictive Check (VPC) obtained with 100 replicates of each animal. The observed quantiles (10, 50, and 90%) were well superimposed on the corresponding predictive check quantiles over the observed data. Theoretically, ~20% of the data should be outside the plotted quantiles. The red lines are 10, 50, and 90% quantiles from the actual observed values. The black lines are the 10, 50, and 90% quantiles from the simulated observations (left). Blue and red shaded (right) areas correspond to the 95% confidence interval of the three predicted quantiles.

TABLE 2

Table 2. Population primary (Thetas) and secondary parameters and random effects (Omega) for florfenicol in calves obtained with a 2-compartment model.

The coefficient of variation of the multiplicative component of the residual was 14%. The BSV for the estimated parameters was ~20–30% but was 57% for ka, suggesting a homogeneous exposure between animals for these formulations. This is consistent with BSV of AUC, as estimated approximately by others when reporting observed AUC mean and SD (Soback et al., 1995; Sidhu et al., 2014).

Using the developed population pharmacokinetic model and estimated parameters, 5,000 curves were generated by Monte Carlo Simulation (simulated IPREDs taking into account Thetas and Omega but not Sigma, the residual error) over 96 h with a step of 1 h; the simulated dosage regimen 40 mg/kg (single subcutaneous administration). The corresponding AUC from 0 to 96 h and the time for which plasma concentrations remained above selected MICs are given in Table 3.

TABLE 3

Table 3. AUC (0–96 h), average concentration (μg/ml), and Time (h) above possible MICs ranging from 0.25 to 2 μg/ml for selected quantiles and corresponding value of the T>MIC in % of 96 h, the claimed duration of action of Nuflor^®.

Data presented in Table 3 indicates that, for a MIC of 1.0 μg/ml, 90% of calves achieved a time above the MIC of at least 38.70 h i.e., a T>MIC of 40.31% of the duration of the assumed florfenicol activity of 96 h. Accepting the claim of the company licensing the pioneer product that the duration of action of Nuflor^® is 96 h (Anonymous, 2018b) and a default PDT value of 40% (Mouton et al., 2012), the florfenicol PK/PD_CO for T>MIC was 1.0 μg/ml, because, for a higher MIC of 2 μg/ml, a T>MIC of 40% was achieved in only 10% of calves. In accepting AUC/MIC as the appropriate index, the average concentration over 96 h achieved by at least 90% of calves was 1.18 μg/ml. Considering the nearest two-fold MIC value, the PK/PD_CO for this index was also 1 μg/ml. This is equivalent to a classical AUC/MIC of 24 h per day in steady-state conditions, as traditionally expressed in human medicine (Toutain et al., 2007). For AUC/MIC values >24 h, the current dosage regimen would not cover 90% of the population; only 10% of calves would be able to achieve an AUC/MIC of 48 h (equivalent to an average plasma concentration of 2 μg/ml over the 96 h interval). An average plasma concentration of 2 μg/ml is equal to the VAST/CLSI CBP.

Discussion

FDA guidance indicates that population PK modeling (Food Drug Administration, 1999) is the only appropriate tool to allow the meta-analysis of data retrieved from different unbalanced designs i.e., study designs in which all individuals do not supply the same amount of information. For the present analysis, the differences in LLOQ of the analytical technique initially prevented direct comparison of the data set obtained by Sidhu et al. (2014) which fitted a one-compartment model (results not shown) with more recent data obtained with a more sensitive analytical technique and best fitted to a 2-exponential model. Population modeling enabled the older, but nevertheless informative data, to be used to generate a single set of parameters (with SE) for florfenicol. This further enabled generation by MCS of a virtual in silico calf population for PK/PD_COs.

Florfenicol disposition in calves has been investigated following administration by the intravenous route (Varma et al., 1986); PK parameters were estimated with a plasma clearance of 2.85 ml/kg/min, a steady-state volume of distribution, Vss, of 0.75 L/kg, and an elimination t1/2 of 2.86 h. Similar results were reported for different types of cattle, including dairy cattle (Soback et al., 1995), dry cows (Bretzlaff et al., 1987), and steers (Lobell et al., 1994), suggesting no major differences in the florfenicol disposition profile in different classes of cattle. Hence, it is likely that the present findings will be representative of and applicable to differing types of cattle. It is also concluded that a single CBP for cattle can be proposed for these LA formulations. In the present analysis, t1/2 values were much longer than after IV administration, with t1/2 of 16 and 67 h for the Alpha and Beta phases, respectively. Just as a CBP depends on a specific dosage regimen (Heil and Johnson, 2016), the computed PK/PD_CO (in the present analysis) is specific for these LA formulations, administered SC as a single dose of 40 mg/kg with an assumed duration of effect of 4 days. It cannot be assumed that the findings apply to any other dosage regimens and/or other formulations and/or other routes of administration. For example, for another LA florfenicol formulation, it has been shown that the mean differences between a SC and intramuscular (IM) administration were as high as 35 and 63%, respectively, for AUC and Cmax, the IM administration route thus leading to higher florfenicol exposure than subcutaneous dosing (Lacroix et al., 2011). This is typical for veterinary medicine, in which many modalities of AMD administration exist, rendering a universal and robust CBP difficult to propose (Toutain et al., 2017).

The very long terminal t1/2 is explained by flip-flop PK, with the Alpha phase corresponding to a first process of relatively slow absorption and the terminal Beta phase corresponding to a very slow absorption process. It is concluded that the respective contributions of the Alpha and Beta phases to the total AUC were ~60% of the bioavailable florfenicol fraction absorbed in the Alpha phase and 40% in the Beta phase. This second phase is not well-characterized in several publications having a rather high analytical method LLOQ. Nevertheless, the population model allows incorporation of all data in calculating the PK/PD_CO.

The flip-flop PK profile of the investigated florfenicol formulation is also the basis for choosing to parameterize the model in terms of macroconstants, rather than in terms of clearance and volume of distribution, as is usually the case in population modeling. Indeed, the aim was to simulate 5,000 curves and, whatever the parameterization, the plasma concentration vs. time curves will be the same. From a mechanistic point of view, however, it is important to explore, when estimating PK/PD_CO, the influence (or not) of the two major covariates involved, namely a possible “formulation” effect (here Nuflor^® vs. generic) and a possible “source” effect (here three sources). For both covariates, it is the relative bioavailability that may differ rather than clearance, which is not determinable when only extravascular data are available. To explore the influence of the covariates in question, the Alpha and Beta slopes are the two parameters to be estimated as primary rather than secondary parameters—hence our parameterization.

The final objective of this population pharmacokinetic analysis was to determine a possible PK/PD_CO for florfenicol, this being the pivotal parameter considered by VetCAST in the decision making process for establishing a CBP. PK/PD_CO provides insight into the overall PK variability across the targeted populations, because of the relationship between drug exposure and efficacy. This relationship is expressed through PK/PD indices (AUC/MIC ratio or T>MIC) which should achieve critical values to predict clinical efficacy. The magnitude of a PK-PD index providing an appropriate level of predicted response is the PDT (European Medicines Agency, 2015).

As explained in Materials and Methods, no dose-fractionation has been conducted in rodents to determine the best PK/PD index predictive of florfenicol efficacy and in this project an in silico approach was used as a surrogate for a dose fractionation trial. It was concluded that the best index was AUC/MIC. This is consistent with the opinion that AUC/MIC is always the most relevant index, when the terminal half-life is long (Nielsen et al., 2011). In addition, it was established that the PDT should be ~24 h per day, indicating that, to achieve an in silico bacteriological eradication, the average florfenicol concentration over the 4 days should be equal to the MIC [see Toutain et al. (2007) for explanation of the relationship between PDT expressed in h vs. as a scaling factor]. This is slightly lower than the bactericidal PDT reported from the killing action of florfenicol against MH and PM from modeling of the time-kill data after 24 h exposure of florfenicol to a constant concentration (Illambas et al., 2013).

From the 5,000 curves generated by MCS, the average plasma concentration was estimated to be 1.2 μg/ml (Table 2) and 1 μg/ml is the PK/PD_CO value for the AUC/MIC index. An identical PK/PD_CO of 1 μg/ml has been derived for florfenicol in pigs, but for a dose of 30 mg/kg (Lei et al., 2018). Florfenicol has been classified as a bacteriostatic drug. It can therefore be argued that T>MIC is also valuable in respect of detection of resistance. However, as quoted by others (Dudley and Ambrose, 2000), the dual aim of achieving a single breakpoint to predict both clinical outcome and avoidance of resistance is likely to fail in many circumstances and constitutes a source of confusion. Nevertheless, the PK/PD_CO for T>MIC was computed, assuming that the current dosage regimen should ensure a T>MIC for ~40% of the duration of treatment in 90% of animals: a critical MIC of 1.0 μg/ml was obtained, a value identical with the critical value for the AUC/MIC index. In the present data analysis a PTA of 90% was used to compute the PK/PD_CO, a quantile that is routinely used for PTA analysis (Turnidge and Paterson, 2007). It should be noted that the quantile 90% is related to the concept of prediction interval and not to the concept of confidence interval. Moreover, the PK/PD_CO as applied by VetCAST is not equivalent to the EUCAST PK/PD breakpoint, as the latter takes account additionally of clinical data (Mouton et al., 2012).

In this investigation florfenicol binding to plasma protein was discounted; as discussed in the Introduction, very disparate figures have been reported in cattle, ranging from no binding to binding of 10 to ~25%. Recently, binding was reported with a wide BSV: values in a cohort of 20 calves ranged from 1.88 to 57.5% in 7 day old calves and from 16.8 to 27.8% in 46 day old calves, both at a florfenicol concentration of 1 mg/L (Mzyk et al., 2018). In addressing these differences, we considered that, for the present study, it was appropriate to consider that plasma protein binding of florfenicol was negligible, as suggested by Foster et al. (2016). This approach will render easier any possible future update of computed PK/PD_CO for the selected PK/PD index (i.e., AUC/MIC), because the extent of binding is simply a scaling factor for this PK/PD index. Equally important, it seems probable that the variability reported by Mzyk et al. (2018) is not simply associated with some technical issue, but rather actually reflects a true BSV. At present, to the best of our knowledge, this variability is not factored into models used, in veterinary medicine, to compute the PTA using MCS. It is the average value which is adopted. If a wide BSV for protein binding was, in due course, confirmed for florfenicol in cattle, it would be necessary not only to scale our results but to re-run the population model to include this source of variability. To summarize, what has been determined in this paper, as a PK/PD_CO, is the simplest hypothesis of no plasma protein binding (and thus no variability for this factor) for florfenicol.

In conclusion, any CBP is both dose- and exposure-dependent. In human medicine most AMDs are administered by the oral route and CBPs have a generic value for oral formulations that are relatively similar in terms of the internal exposure they provide. This is unfortunately not the case for veterinary medicine, where CBPs can also be “formulation-dependent.” The formulations, Nuflor^® from MSD and its generics, evaluated in this study were all administered by the subcutaneous route, these being the most extensively used formulations and route of administration for florfenicol in cattle. However, other florfenicol formulations and other routes of administration are used in cattle, so that VetCAST CBP is not guaranteed to be applicable to other formulations and/or other routes of administration. These issues are discussed in the VetCAST position paper (Toutain et al., 2017). Finally from a pooled raw data analysis, using a NLME and MCS for florfenicol, a PK/PD_CO of 1 mg/L is proposed for the extensively used LA florfenicol formulations investigated.

Data Availability

The florfenicol raw data are available in the Supplemental File.

Author Contributions

PL, AR, and PS generated raw data. LP retrieved and validated raw data. P-LT performed the modeling analysis and drafted the paper. All co-authors critically reviewed several drafts of the manuscript.

Funding

This work was partly supported by the Direction Générale de l'Alimentation (DGAL) of the French Ministry of Agriculture and Food. DGAL has no role in data collection, interpretation and the decision to submit this work for publication. A part of this study was supported by a generous grant from the Department for the Environment, Food and Rural Affairs (United Kingdom). PS was sponsored by a UKIERI fellowship.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The VetCAST steering committee is acknowledged for supporting this publication.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01310/full#supplementary-material

References

Adams, P. E., Varma, K. J., Powers, T. E., and Lamendola, J. F. (1987). Tissue concentrations and pharmacokinetics of florfenicol in male veal calves given repeated doses. Am. J. Vet. Res. 48, 1725–1732.

PubMed Abstract | Google Scholar

Andes, D., and Craig, W. A. (2002). Animal model pharmacokinetics and pharmacodynamics: a critical review. Int. J. Antimicrob. Agents 19, 261–268. doi: 10.1016/S0924-8579(02)00022-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Anonymous (2018a). NORFENICOL- Florfenicol Injection, Solution Norbrook Laboratories Limited. Available online at: https://www.norbrook.com/products/united-states/norfenicol-injectable-solution-florfenicol (accessed August 22, 2018).

Google Scholar

Anonymous (2018b). Nuflor - Frequently Asked Questions. Available online at: http://www.nuflor.com/nuflor_glance/faq.asp (accessed August 25, 2018).

Google Scholar

Bonate, P. L. (2011). Pharmacokinetic-Pharmacodynamic Modeling and Simulation, 2nd Edn. New York, NY: Springer.

Google Scholar

Bretzlaff, K. N., Neff-Davis, C. A., Ott, R. S., Koritz, G. D., Gustafsson, B. K., and Davis, L. E. (1987). Florfenicol in non-lactating dairy cows: pharmacokinetics, binding to plasma proteins, and effects on phagocytosis by blood neutrophils. J. Vet. Pharmacol. Ther. 10, 233–240. doi: 10.1111/j.1365-2885.1987.tb00534.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Byon, W., Fletcher, C. V., and Brundage, R. C. (2008). Impact of censoring data below an arbitrary quantification limit on structural model misspecification. J. Pharmacokinet. Pharmacodyn. 35, 101–116. doi: 10.1007/s10928-007-9078-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Clinical and Laboratory Standards Institute (2018). Performance Standards for Antimicrobial Disk and Dilution Susceptibility Tests for Bacteria Isolated From Animals. 4th Edn. CLSI supplement VET08. Wayne, PA: Clinical and Laboratory Standard Institute.

Google Scholar

Craig, W. A. (1998). Pharmacokinetic/pharmacodynamic parameters: rationale for antibacterial dosing of mice and men. Clin. Infect. Dis. 26, 1–10; quiz 11–12. doi: 10.1086/516284

PubMed Abstract | CrossRef Full Text | Google Scholar

de Jong, A., Thomas, V., Simjee, S., Moyaert, H., El Garch, F., Maher, K., et al. (2014). Antimicrobial susceptibility monitoring of respiratory tract pathogens isolated from diseased cattle and pigs across Europe: the VetPath study. Vet. Microbiol. 172, 202–215. doi: 10.1016/j.vetmic.2014.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudley, M. N., and Ambrose, P. G. (2000). Pharmacodynamics in the study of drug resistance and establishing in vitro susceptibility breakpoints: ready for prime time. Curr. Opin. Microbiol. 3, 515–521. doi: 10.1016/S1369-5274(00)00132-6

PubMed Abstract | CrossRef Full Text | Google Scholar

European Medicines Agency (2015). Guideline on the Use of Pharmacokinetics and Pharmacodynamics in the Development of Antimicrobial Medicinal Products. Available online at: arxiv:/www.ema.europa.eu/docs/en_GB/document_library/Scientific_guideline/2016/07/WC500210982.pdf

Google Scholar

Food Drug Administration (1999). Guidance for Industry; Population Pharmacokinetics.U.S. Department of Health and Human Services Food and Drug Administration ;Center for Drug Evaluation and Research (CDER); Center for Biologics Evaluation and Research (CBER). Available online at: https://www.fda.gov/downloads/drugs/guidances/UCM072137.pdf (accessed August 22, 2018).

Google Scholar

Foster, D. M., Martin, L. G., and Papich, M. G. (2016). Comparison of active drug concentrations in the pulmonary epithelial lining fluid and interstitial fluid of calves injected with enrofloxacin, florfenicol, ceftiofur, or tulathromycin. PLoS ONE 11:e0149100. doi: 10.1371/journal.pone.0149100

CrossRef Full Text | Google Scholar

Heil, E. L., and Johnson, J. K. (2016). Impact of CLSI breakpoint changes on microbiology laboratories and antimicrobial stewardship programs. J. Clin. Microbiol. 54, 840–844. doi: 10.1128/JCM.02424-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Illambas, J., Potter, T., Sidhu, P., Rycroft, A. N., Cheng, Z., and Lees, P. (2013). Pharmacodynamics of florfenicol for calf pneumonia pathogens. Vet. Rec. 172:340. doi: 10.1136/vr.101155

PubMed Abstract | CrossRef Full Text | Google Scholar

Karlsson, M. O., and Savic, R. M. (2007). Diagnosing model diagnostics. Clin. Pharmacol. Ther. 82, 17–20. doi: 10.1038/sj.clpt.6100241

PubMed Abstract | CrossRef Full Text | Google Scholar

Kristoffersson, A. N., David-Pierson, P., Parrott, N. J., Kuhlmann, O., Lave, T., Friberg, L. E., et al. (2016). Simulation-based evaluation of PK/PD Indices for meropenem across patient groups and experimental designs. Pharm. Res. 33, 1115–1125. doi: 10.1007/s11095-016-1856-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lacroix, M. Z., Gayrard, V., Picard-Hagen, N., and Toutain, P. L. (2011). Comparative bioavailability between two routes of administration of florfenicol and flunixin in cattle. Rev. Med. Vet. 162, 321–324.

Google Scholar

Lei, Z., Liu, Q., Yang, S., Yang, B., Khaliq, H., Li, K., et al. (2018). PK-PD Integration Modeling and Cutoff Value of Florfenicol against Streptococcus suis in Pigs. Front. Pharmacol. 9:2. doi: 10.3389/fphar.2018.00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Gehring, R., Lin, Z., and Riviere, J. (2015). A framework for meta-analysis of veterinary drug pharmacokinetic data using mixed effect modeling. J. Pharm. Sci. 104, 1230–1239. doi: 10.1002/jps.24341

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobell, R. D., Varma, K. J., Johnson, J. C., Sams, R. A., Gerken, D. F., and Ashcraft, S. M. (1994). Pharmacokinetics of florfenicol following intravenous and intramuscular doses to cattle. J. Vet. Pharmacol. Ther. 17, 253–258. doi: 10.1111/j.1365-2885.1994.tb00241.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez, M. N., Toutain, P. L., and Turnidge, J. (2013). “The pharmacodynamics of antimicrobial agents,” in Antimicrobial Therapy in Veterinary Medicine, 5th Edn, eds S. Giguère, J. F. Prescott, and P. M. Dowling (Ames, IA: Wiley Blackwell, 79–103.

Google Scholar

Mouton, J. W., Brown, D. F. J., Apfalter, P., Cantón, R., Giske, C. G., Ivanova, M., et al. (2012). The role of pharmacokinetics/pharmacodynamics in setting clinical MIC breakpoints: the EUCAST approach. Clin. Microbiol. Infect. 18, E37–45. doi: 10.1111/j.1469-0691.2011.03752.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mzyk, D. A., Bublitz, C. M., Hobgood, G. D., Martinez, M. N., Smith, G. W., and Baynes, R. E. (2018). Effect of age on the pharmacokinetics and distribution of tulathromycin in interstitial and pulmonary epithelial lining fluid in healthy calves. Am. J. Vet. Res. 79, 1193–1203. doi: 10.2460/ajvr.79.11.1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, E. I., Cars, O., and Friberg, L. E. (2011). Pharmacokinetic/pharmacodynamic (PK/PD) indices of antibiotics predicted by a semimechanistic PKPD model: a step toward model-based dose optimization. Antimicrob Agents Chemother. 55, 4619–4630. doi: 10.1128/AAC.00182-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, E. I., and Friberg, L. E. (2013). Pharmacokinetic-pharmacodynamic modeling of antibacterial drugs. Pharmacol. Rev. 65, 1053–1090. doi: 10.1124/pr.111.005769

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Connor, A. M., Yuan, C., Cullen, J. N., Coetzee, J. F., da Silva, N., and Wang, C. (2016). A mixed treatment meta-analysis of antibiotic treatment options for bovine respiratory disease - an update. Prev. Vet. Med. 132, 130–139. doi: 10.1016/j.prevetmed.2016.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelligand, L., Lees, P., Sidhu, P. K., and Toutain, P-L. (2019). Semi-mechanistic modelling of florfenicol time-kill curves and in silico dose fractionation for calf respiratory pathogens. Front. Microbiol. doi: 10.3389/fmicb.2019.01237. [Epub ahead of print].

CrossRef Full Text | Google Scholar

Schoemaker, R. C., and Cohen, A. F. (1996). Estimating impossible curves using NONMEM. Br. J. Clin. Pharmacol. 42, 283–290. doi: 10.1046/j.1365-2125.1996.04231.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sidhu, P., Rassouli, A., Illambas, J., Potter, T., Pelligand, L., Rycroft, A., et al. (2014). Pharmacokinetic-pharmacodynamic integration and modelling of florfenicol in calves. J. Vet. Pharmacol. Ther. 37, 231–242. doi: 10.1111/jvp.12093

PubMed Abstract | CrossRef Full Text | Google Scholar

Soback, S., Paape, M. J., Filep, R., and Varma, K. J. (1995). Florfenicol pharmacokinetics in lactating cows after intravenous, intramuscular and intramammary administration. J. Vet. Pharmacol. Ther. 18, 413–417. doi: 10.1111/j.1365-2885.1995.tb00618.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Toutain, P.-L., Bousquet-Mélou, A., Damborg, P., Ferran, A. A., Mevius, D., Pelligand, L., et al. (2017). En route towards european clinical breakpoints for veterinary antimicrobial susceptibility testing: a position paper explaining the VetCAST approach. Front. Microbiol. 8:2344. doi: 10.3389/fmicb.2017.02344

PubMed Abstract | CrossRef Full Text | Google Scholar

Toutain, P.-L., Bousquet-Mélou, A., and Martinez, M. (2007). AUC/MIC: a PK/PD index for antibiotics with a time dimension or simply a dimensionless scoring factor? J. Antimicrob. Chemother. 60, 1185–1188. doi: 10.1093/jac/dkm360

PubMed Abstract | CrossRef Full Text | Google Scholar

Toutain, P. L., and Bousquet-Mélou, A. (2004). Plasma terminal half-life. J. Vet. Pharmacol. Ther. 27, 427–439. doi: 10.1111/j.1365-2885.2004.00600.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnidge, J., and Paterson, D. L. (2007). Setting and revising antibacterial susceptibility breakpoints. Clin. Microbiol. Rev. 20, 391–408. doi: 10.1128/CMR.00047-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnidge, J. D., and Martinez, M. N. (2017). Proposed method for estimating clinical cut-off (COCL) values: an attempt to address challenges encountered when setting clinical breakpoints for veterinary antimicrobial agents. Vet. J. 228, 33–37. doi: 10.1016/j.tvjl.2017.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Varma, K. J., Adams, P. E., Powers, T. E., Powers, J. D., and Lamendola, J. F. (1986). Pharmacokinetics of florfenicol in veal calves. J. Vet. Pharmacol. Ther. 9, 412–425. doi: 10.1111/j.1365-2885.1986.tb00062.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X. S., Dunne, A., Kimko, H., Nandy, P., and Vermeulen, A. (2011). Impact of low percentage of data below the quantification limit on parameter estimates of pharmacokinetic models. J. Pharmacokinet. Pharmacodyn. 38, 423–432. doi: 10.1007/s10928-011-9201-9

PubMed Abstract | CrossRef Full Text | Google Scholar