Optimal control for colistin dosage selection.

Optimization of antibiotic administration helps minimizing cases of bacterial resistance. Dosages are often selected by trial and error using a pharmacokinetic (PK) model. However, this is limited to the range of tested dosages, restraining possible treatment choices, especially for the loading doses. Colistin is a last-resort antibiotic with a narrow therapeutic window; therefore, its administration should avoid subtherapeutic or toxic concentrations. This study formulates an optimal control problem for dosage selection of colistin based on a PK model, minimizing deviations of colistin concentration to a target value and allowing a specific dosage optimization for a given individual. An adjoint model was used to provide the sensitivity of concentration deviations to dose changes. A three-compartment PK model was adopted. The standard deviation between colistin plasma concentrations and a target set at 2 mg/L was minimized for some chosen treatments and sample patients. Significantly lower deviations from the target concentration are obtained for shorter administration intervals (e.g. every 8 h) compared to longer ones (e.g. every 24 h). For patients with normal or altered renal function, the optimal loading dose regimen should be divided into two or more administrations to attain the target concentration quickly, with a high first loading dose followed by much lower ones. This regimen is not easily obtained by trial and error, highlighting advantages of the method. The present method is a refined optimization of antibiotic dosage for the treatment of infections. Results for colistin suggest significant improvement in treatment avoiding subtherapeutic or toxic concentrations.

Introduction

The treatment of some infections is a real challenge in the clinical area, since certain bacteria have become resistant to most of available antibiotics, and there are only few new drugs being developed against these bacteria [1]. Therefore, the proper use of the remaining antibiotic options is required to avoid these cases, choosing the most suitable one and the optimal regimen of the treatment [2].

The minimum inhibitory concentration (MIC) is a relevant pharmacodynamic (PD) parameter [3, 4]. The consideration of MIC with other pharmacokinetic (PK) parameters, such as area under the curve (AUC), peak plasma concentration (Cmax) and time period of plasma concentration above the MIC, helps in choosing the proper use of antibiotics [3-5]. However, some antibiotics have a narrow therapeutic window. In these cases, the plasma concentration should be controlled and cannot be several times higher than the MIC [6].

Modern models and simulations integrating PK and PD are being employed in order to predict plasma concentrations and suggest a suitable treatment scheme [7-9]. Usually, a range of doses is tested and the best dose is selected considering the fittest pharmacokinetic profile produced [10, 11]. Such trial and error procedures have undoubtedly contributed to the appropriate use of antibiotics. However, the dose selection is limited to the range of tested doses, which is often based on an initial loading dose followed by a number of equal maintenance doses. While this could lead to satisfactory effects for most drugs, there are some molecules that may require a more refined treatment optimization, especially in cases where there is a narrow therapeutic window and a target concentration should be attained.

Colistin is an old drug which was commercially available in the 1950s but rapidly replaced by other drugs due to its nephrotoxicity [12, 13]. Nonetheless, cases of resistance to these other drugs emerged and colistin has been shown to be a viable therapeutic option against pathogens which are very difficult to treat, such as Gram-negative bacteria [14, 15]. For this reason, interest in colistin was renewed in the 2000s, and it has been used as a last resort antibiotic [12]. Intravenous (IV) administration of colistin is often used in cases of nosocomial pneumonia (ventilator-associated or not) due to multidrug-resistant (MDR) or extensively drug-resistant (XDR) pathogens, such as Acinetobacter baumannii, Pseudomonas aeruginosa and Klebsiella pneumoniae [16, 17], since after IV administration, the molecule is able to cross the blood-air barrier and reach the lung tissues [18, 19]. Colistin is administered as an inactive prodrug, colistimethate sodium (CMS), to mitigate its toxicity [12]. Once administered, CMS is hydrolyzed in colistin, most probably in the bladder of the patients, and approximately 70% of the CMS dose is rapidly eliminated in the urine [20, 21]. Therefore, the dose should be adjusted for patients presenting renal problems [22, 23].

Although it has been extensively used in the recent decades, the optimal dosing regimen of colistin remains to be determined [21, 24]. A key aspect is that plasma concentrations of colistin below 2 mg/L are ineffective [25, 26], and concentrations above 2.42 mg/L lead to toxicity issues [27]. Therefore, colistin administration would greatly benefit from an optimization targeting the safe plasma concentration. Target steady-state concentrations may be obtained using standard methods for linear PK modes [28], but the optimization should also involve the loading dose, in order that the target concentration be rapidly attained. Furthermore, the availability of PK models may allow such optimization to be performed for a specific patient, taking into account how PK parameters vary in the population.

The aim of this study is to use a PK model based optimization to assist the selection of the optimal therapeutic scheme of the colistin. Instead of a selection of doses by trial and error, which invariably limits the possible choices, we formulate an optimization problem that leads to an optimal control formulation. Control theory has been used to optimize treatment in problems with available PK/PD models, such as HIV infections [29, 30], cancer therapy [31-35], or to optimize strategies for the prevention of infectious diseases such as cholera [36], Zika virus infection [37] and COVID-19 [38]. For colistin administration, the optimization problem minimizes the deviations between the plasma concentration of colistin and its target value. Additionally, the optimal regimen was individualized according to the creatinine clearance of the patients. While this work focuses on the specific case of colistin, the formulation is general and allows optimization of administration of other drugs whenever a PK model is available.

Methodology

Definitions

We consider the three-compartment PK model from [22], in a simpler form for patients that are neither on hemodyalisis, nor on continuous renal replacement treatment. For these conditions, the model is given by

The model variables are the mass of CMS in the central compartment (CMSc), the mass of CMS in the peripheral compartment (CMSp) and the mass of colistin in a single colistin compartment (Colistin).V1, V2 and V3 denote volumes of distribution of the CMSc, CMSp and Colistin compartments, respectively, and CLD1 is the distributional clearance between the central and peripheral compartments for CMS. The administration of CMS in the central compartment is given by R in the above model. The variables in the model refer to the total (bound plus unbound) CMS and colistin. Accordingly, in what follows we will deal solely with total concentrations.

CLT and CLT are, respectively, the total intrinsic clearances for CMS and colistin. The total CMS clearance (CLT) was modeled as a function of creatinine clearance (CrCL) and two effects, a renal dependent component (CLR) and nonrenal clearance (CLNR), as

Similarly, the colistin clearance (CLT) was modeled as a sum of a renally dependent component (CLRC) and a nonrenally dependent component (CLNR), as

A schematic pictogram representing the PK model is shown in Fig. 1. The parameter values used in the simulations are detailed in Table 1.

Fig. 1

Structural simultaneous PK modeling:, and volumes of distribution of the central colistimethate sodium , peripheral CMS and colistin compartments, respectively, distributional clearance between the central and peripheral compartments for CMS, administration of CMS in the central compartment, CrCL creatinine clearance, CLR renal dependent component, CLNR nonrenal clearance, total intrinsic clearance for colistin

Table 1

PK parameters and covariate values from previous study (22) after IV infusion (R) of CMS in the central compartment used in simulations

Category	Parameter	(UNITS)	Value
CMS	V1	(L)	11.9
	V2	(L)	18.7
	CLD1	(L/h)	7.98
	CLRslope	(L/h/CrCL)	0.0613
	CLNRCMS	(L/h)	1.9
Colistin	V3	(L)	45.1
	CLTC	(L/h)	2.72
	CLRCSL_POP	(L/h/CrCL)	0.0147
	CLNRC	(L/h)	2.19
	Weight	(Kg)	60

V1, V2, and V3 volumes of distribution of the central colistimethate sodium (CMS), peripheral CMS and colistin compartments, respectively, distributional clearance between the central and peripheral compartments for CMS, CLR renal dependent component for CMS, CLNR nonrenal clearance of CMS, total intrinsic clearance for colistin; CLRC renally dependent component for colistin, CLNR nonrenally dependent component for colistin

It is convenient to write the system in matrix form. We define a state as a column vector,the system matrix asand the input matrix assuch that the system becomes simplywith a scalar input given by , i.e. the infusion rate of CMS.

The quantity of primary interest in the system is the concentration of colistin in its compartment, given by . This output, defined y, may be extracted by defining an output matrix C given bysuch that

The system in Eq. (9). is solved subject to an initial condition, given by . Here, the drug administration is assumed to occur after the initial time , hence the initial state is given by . Once the solution is obtained, the colistin concentration is obtained with Eq. (11).

The definition of the system in matrix form is convenient, as will be shown next, and also provides generality. The optimization problem to be defined may be adapted in straigthforward manner to other linear compartmental PK models once they are cast in matrix form by the definitions of system, input and output matrices , and , respectively.

Optimization problem

We define a target concentration , which may be constant or time-dependent depending on the application. The present optimization problem is set so as to minimize deviations between the model output and the target concentration. For a target concentration that is constant in time, a maintenance dose may be selected by standard methods for linear PK models [28]. On the other hand, obtaining optimal loading doses is not trivial, as it is important to quickly attain the desired concentration for critically ill patients without an overshoot that would possibly be toxic. Selection of does that optimally attain a target concentration that varies with time (i.e. a desired curve of concentration versus time) is also difficult, even for linear models. It is with such purposes that we define the optimization below.

We define the functional as

is related to the mean square value of the deviations between the colistin concentration and its target value. The quantity is thus a standard deviation between the colistin concentration and the target, a metric to evaluate if concentrations are close to the target value. We would thus like to obtain the optimal drug administration , such that is minimized; a minimisation of implies a minimisation of the standard deviation γ to the target concentration. Notice that the specification of a target concentration is arbitrary. For instance, one may specify a desired concentration curve as a function of time, considered to be ideal for clinical application. The procedure will then determine a drug administration that minimizes deviations to such target curve .

The present problem is a constrained optimization, where one wishes to minimize the functional with respect to u such that the constraints given by Eqs. (9) and (11) are satisfied. This may be solved using Lagrange multipliers or adjoint variables, as in optimal control works dealing with cancer therapy [32, 35] or vaccination strategies [37, 38]. The procedure is based on the application of Pontryagin’s Maximum Principle [39], which we describe in what follows. We define an augmented functionalwhere is an adjoint variable, introduced to ensure that the constraint is satisfied, simplifying the solution of the problem. Superscript denotes a matrix transpose. Notice that if satisfies the PK model of Eq. (9) the last integral is equal to zero, and the augmented functional is equal to . Therefore, minimization of will lead to minimisation of , and, in turn, of the standard deviation γ.

The unconstrained minimisation of the augmented functional is equivalent to the constrained minimisation of . To show this, consider that once is minimized its variation with respect to should be zero. Denoting such variation by , defined bywe obtainwhere denotes a small arbitrary variation introduced in . Setting for any implies thatwhich means that is solution to the model in Eq. (9), and the constraint is satisfied.

We now consider the variation of with respect to the state , given asAn integration by parts can be applied to the time derivative of , leading toFor arbitrary , this variation can be set to zero by the definition of an adjoint system. We first consider that the initial condition for is fixed, which implies . We also specify that the adjoint variable is a solution of the adjoint system given byor, taking transposes of all terms,subject to the terminal condition . The adjoint system is solved backwards in time, starting from the terminal condition in order to obtain .

Finally, the variation of with respect to the input provides the sensitivity to drug application. It is given byand thus the changes in input leading to the highest decrease in are

Optimization algorithm

The sensitivities are used in an optimal-control algorithm to find the drug administration that minimizes . We start with zero drug administration. This first input is referred to as . We then set an iterative process, where the variation of with respect to the dose is calculated at each iteration, and an updated guess of the optimal dose is obtained with Eq. (22) in iterations . The steps for iteration are:

Solve the direct PK system given by Eq. (9), with . This provides the direct solution of iteration , , and the corresponding output is obtained from Eq. (11);

Solve the adjoint PK system given by Eq. (20). This provide the adjoint solution of iteration , ;

The change in drug dosage leading to the optimal decrease in the functional , , is given by Eq. (22), and the dose for iteration is updated as , with given as a small constant to ensure small dose updates.

The steps above are shown in schematic form in Fig. 2. Successive iterations are carried out until vanishes. When this happens the variations of with respect to , and are simultaneously zero, and we have thus reached a minimum.

Fig. 2

Schematic representation of the optimization process

Constraints in drug administration

The procedure of the previous section leads to a continuous profile for drug administration, , which is optimal to attain the target concentration but is impractical for application, as it requires a continuous infusion with variable rate. In clinical practice, each infusion is applied with a constant rate over a period of time, which for colistin could be one hour [18, 40]. This amounts to a restriction to , which should be constant during an infusion and zero otherwise. An infusion rate satisfying such constraint may be obtained from the optimal as

where the overbar denotes a time average within an infusion period. Hence, changes to the infusion rate are constant in each specified infusion period, with a value equal to the time average of the optimal change within each period, and are zero for times between infusions. This ensures that the resulting optimization of drug administration satisfy the practical constraint of constant infusion rate within a fixed period. In this study, we have used a fixed infusion period of one hour, and tested different intervals between infusions. Periods of 0.5 to 1 h have been recommended in the literature [26]. If necessary, changing the infusion time in the optimization is straightforward. Some tests were carried out with a 2 h infusion period, leading to minimal differences to the reported optimal doses for 1 h infusions.

Extensions

The above procedure allows extensions for other optimization scenarios of interest, which are described in what follows.

Inter-individual variability

The optimal control problem was here formulated for a deterministic model, with , and matrices held fixed by considering median values of pharmacokinetic parameters. Nonetheless, the optimization may be extended by considering a population analysis. A straightforward approach is the use of a Monte Carlo method to generate a population of patients with varying pharmacokinetic parameters (clearances and volumes of distribution), which in the framework leads to separate , and matrices for each patient i. In this case, an appropriate objective functional iswhich sums deviations of the concentration to the target value for each patient. With this choice, the sensitivity of to changes in the dose is simply the sum of the sensitivity to each patient. Simulations of the direct and adjoint PK systems, considering matrices , and , provide such sensitivities and allow updating the dose in the optimization loop. The optimized dose so obtained minimizes the objection in Eq. (24), thus leading to the dose that leads to the lowest deviations of concentration from its target value once averaged among the considered individuals.

In this case, an interesting possibility is the inclusion of toxicity considerations. The objective functional of Eq. (24) may be supplemented with penalization terms related to patients that reach toxic concentrations, similar to Pierillas et al. [41]. Such penalization terms appear as costs introduced to the objective functional, as in Kirschner et al. [30] and Joshi [29]. The introduction of penalization terms to the objective functional and subsequent derivation of the optimal control problem should reduce toxicity cases. Notice however that the use of a target concentration of 2 mg/L in the present framework already limits toxicity, as the optimization naturally converges to this concentration avoiding values significantly above the target.

Nonlinear models

The optimal control problem was formulated for a linear PK model due to its widespread use for most drugs [28], and, in particular, to colistin [22], which is the subject of this work. This provides a general framework that may be applied for any linear PK model, once written in the form of Eq. (9). Optimal control can also be applied to nonlinear models, such as in the works of Kirschner et al. [30] and Joshi [29] on optimal control of HIV chemotherapy or Sharp et al. [35] on acute myeloid leukaemia. The framework is the same, with the definition of an augmented functional . The adjoint equation is no longer Eq. (20), and needs to be derived for the specific nonlinear model at hand. Once derived and solved, the adjoint solution provides the sensitivity to drug administration.

Concentration at a different compartment

In this work we have considered a target colistin plasma concentration of 2 mg/L, as previous works [25, 26]. For models including an ELF compartment, such as Boisson et al. [18], it is possible to specify a target concentration in the ELF by adopting the appropriate output matrix . More generally, the selection of output compartment and target concentration in a given PK model is arbitrary, and allows dose optimization for any desired concentration profile.

Numerical methods

The procedure outlined in this section was implemented in a MATLAB script. Simulations of the direct and adjoint linear systems were carried out using the “lsim” routine. Time was discretized in intervals of 0.1 h, and reductions of this value did not lead to significant changes in the results. We have checked the implementation of the PK model by simulating cases from [22], with results consistent to the published ones. The model parameters are reported in Table 1. The small constant to adjust the drug administration was chosen to ensure small changes in at each iteration. The process was stopped when relative changes to were lower than .

Results

Optimal dosage for a patient with normal renal function

As defined in the previous section, lower values of γ indicate a lower standard deviation between the total (bound plus unbound) colistin concentration and the target value. Therefore, the lower the γ value, the closer the colistin plasma concentration will be to the target, which was taken here as  = 2 mg/L [25, 26].

Figure 2 shows the results of dosage optimization, considering patients with normal renal function (creatinine clearance = 80 mL/min), with different infusion regimens. The lowest determined value of γ was 0.16 mg/L (Fig. 3b) for a dosage regimen with an infusion every 8 h. For infusions every 12 and 24 h, the values of γ were higher: 0.19 mg/L (every 12 h) and 0.36 mg/L (every 24 h), respectively (Fig. 3d and f). An administration every 6 h for a patient with normal renal function leads to a value of γ equal to 0.15 mg/L (data not shown), a marginal improvement compared to the value determined for a regimen every 8 h.

Fig. 3

Optimal profile for a patient with normal renal function. Dose optimization for patients with normal renal function (creatinine clearance = 80 mL/min) after an administration every 8 h (a), 12 h (c) and 24 h (e). CMS in the central (dotted line) and peripheral (light gray line) compartments and colistin plasma concentrations (dark gray line) over time for a regimen with administration every 8 h (b) 12 h (d) and 24 h (f). γ: Standard deviation from target concentration 2 mg/L; CMS: colistimethate sodium; MIU: million international units

The optimized IV administration of CMS in different dosing regimens (Fig. 3a, c and e) resulted in profiles of CMS and colistin plasma concentrations as a function of the time shown in Fig. 3b, d and f.

For a treatment every 8 h, the optimal control is a loading dose of approximately 8 million international units (MIU) of CMS (equivalent to about 640 mg of CMS and 264 mg of colistin base activity (CBA)), followed by approximately 0.5 MIU of CMS after 8 h and a maintenance dose, starting after 16 h, of approximately 2.4 MIU of CMS, totalizing 7.2 MIU daily of CMS (Fig. 3a).

For a regimen with administration every 12 h, the optimized loading doses were around 8 MIU of CMS followed by 2.5 MIU after 12 h and the maintenance dose was around 3.6 MIU of CMS (7.2 MIU daily of CMS) (Fig. 3c). For a regimen with infusion every 24 h, the optimal loading dose was around 9 MIU and the maintenance dose was around 6.8 MIU (Fig. 3e). In the regimens with more frequent infusions, the first doses obtained from the optimization procedure form a loading profile, with a first dose that quickly leads to the target value, followed by a second, lower dose that maintains the concentration around the target

Optimal dosage regimen of colistin for patients with different renal conditions

As the PK model is a function of creatinine clearance, separate optimizations should be carried out depending on the renal function. Thus, the optimal dose varies (Fig. 4a, c and e) with the renal function to maintain the colistin plasma concentrations close to the target value of 2 mg/mL (Fig. 4b, d and f). For a patient with creatinine clearance equal to 20 mL/min, the doses need to be reduced, and the slower absorption of colistin leads to more complex optimal loading dose profiles; for an infusion every 8 h, the optimal loading profile of CMS in the first 24 h were about 4, 0, 0.125 and 1 MIU at times 0, 8, 16 and 24 h, respectively. In this case the optimal loading profile involves four initial administrations, so as to attain quickly the target concentration with a high initial dose (Fig. 4a and b), and subsequently maintain the concentration with zero or low doses. The optimal maintenance dose was around 0.8 MIU of CMS every 8 h (Fig. 4a). For infusion every 12 h, the optimal loading profile of CMS is of 4.5, 0 and 1.4 MIU at times 0, 12 and 24 h, respectively (Fig. 4c). For infusions occurring every 24 h, the optimal loading profile is of 4.4 and 2 MIU at times 0 and 24 h, respectively (Fig. 4e). The optimal maintenance doses were 1.2 and 2.4 MIU for infusions every 12 and 24 h, respectively (Figs. 4c and e). The optimal values of γ for administrations every 8, 12 or 24 h were, respectively, 0.18, 0.19 and 0.23 mg/L for a patient with creatinine clearance equal to 20 mL/min (Figs. 4b, d and f).

Fig. 4

Optimal profile for patients with different renal conditions. Dose optimization for patients with creatinine clearance equal to 20 mL/min (black bar) 80 mL/min (white bar) and 150 mL/min (gray bar) after an administration every 8 h (a), 12 h (c) and 24 h (e). Colistin plasma concentrations over time for patients with creatinine clearance equal to 20 mL/min (dotted line) 80 mL/min (light gray line) and 150 mL/min (black line) for a regimen with administrations every 8 h (b) 12 h (d) and 24 h (f). CrCL: creatinine clearance; γ: Standard deviation from target concentration 2 mg/L; CMS: colistimethate sodium; MIU: million international units

For patients presenting creatinine clearance equal to 150 mL/min, the optimal dose is higher as shown in the Figs. 4a, c and e. With infusions every 8 h, the first two doses obtained in the optimization are about 12 and 3.5 MIU, followed by a maintenance dose of 5 MIU. Infusions every 12 h lead to initial doses of 13 and 6.5 MIU and a maintenance dose of 7.5 MIU, and for infusions every 24 h the optimization leads to a single loading dose of about 15.5 MIU followed by a maintenance dose of 13 MIU. The values of γ for administrations every 8, 12 or 24 h were 0.18, 0.25 and 0.54 mg/L, respectively (Figs. 4b, d and f).

Discussion

This study used control theory to optimize drug administration. Colistin was chosen for the analysis since: (i) its ideal dosage regimen is not fully established yet; (ii) it is often used in the clinical practice and is frequently the only effective option for the treatment of certain infections; (iii) it possesses a narrow therapeutic window and should be administered carefully in order to avoid toxicity issues.

The target plasma concentration of colistin (2 mg/L) was selected since this concentration seems to be effective against serious infectious diseases by Gram-negative bacteria [25, 26]. Moreover, a previous study showed that plasma colistin concentrations greater than 2.42 mg/L were more related to nephrotoxicity [27]. Therefore, this target seems to be a satisfactory balance between effectiveness and toxicity. The optimization procedure leads to doses that minimize the standard deviation γ between colistin plasma concentration and its target value.

In the optimal control results, the IV administration every 8 h led to a lower γ value (0.16 mg/L) than for every 12 or 24 h (0.19 and 0.36 mg/L, respectively). The difference is less substantial between administrations every 8 and 12 h, but the values are considerably lower than what is obtained with an administration every 24 h, suggesting that a shorter interval between the infusions could favor the success of the treatment with colistin, as a more frequent administration leads to lower fluctuations of the concentrations, keeping the values close to the target plasma concentration (Fig. 3b, d and f). Additionally, a shorter interval between administrations decreases the possibility of developing resistance to colistin [42, 43].

Previous studies showed that in cases of lung infections caused by bacteria presenting high values of MIC for colistin, the target plasma concentration should be greater than 2 mg/L [25, 44, 45]. As the PK model used in this work is linear, optimal doses for other target concentrations may be obtained directly by multiplying the doses reported here by a factor. For instance, a target concentration of 3 mg/L is obtained by multiplying all doses by 1.5, leading to a daily maintenance dose of about 10.5 MIU of CMS. The maintenance doses (around 7 MIU of CMS daily) suggested by this present model for a patient with a normal renal function are higher than the doses frequently used in the clinical practice (3 to 6 MIU of CMS per day) [21, 46]. However, this suggested dose is lower than the dose recently recommended by the International Consensus Guidelines for the Optimal Use of the Polymyxins for a patient presenting a creatinine clearance of 80 mL/min (10.3 MIU of CMS) [26]. This could be explained by the fact that the latter consensus was based on works [22, 25] that developed a dosing algorithm that did not use the median values of obtained data. Instead of this, the authors performed dosing adjustments such that more than 80% of the patients achieve concentrations equal or higher than the target. We have used the same model [22], but we have used the median values, which means that 50% of the patients would be above and the other 50% would be below the target concentration.

The recommended loading dose could vary between 6 to 9 MIU of CMS with the subsequent first maintenance dose starting after 8, 12 or 24 h [10, 26]. Interestingly, the loading dose obtained as optimal control was different from those often used by physicians. The optimal loading regimen is divided into 2 or more administrations, especially for shorter intervals between the administrations. A high loading dose is initially recommended followed by a dose much lower than the first one. This regimen is interesting for the treatment with colistin, since it allows one to attain more rapidly the target plasma concentration, avoiding plasma concentrations much higher than the target. Such loading regimen deserves further investigation.

One of the advantages of the present method is that it allows the selection of doses for specific patients. As most of CMS is excreted in the urine, it is expected that a reduction in the renal function increases the plasma concentration of CMS and, consequently, that of colistin. This requires a dose reduction, as demonstrated in previous works [22, 25], and also obtained in this study. Interestingly, for a patient presenting creatinine clearance equal to 20 mL/min the γ values after an administration every 8, 12 and 24 h were comparable (0.18, 0.19 and 0.23 mg/L) (Figs. 4b, d and f). This can be explained by the fact that this patient eliminates CMS more slowly, which results in a lower fluctuation of the plasma concentrations over time, keeping these concentrations closer to the target values, even after an administration every 24 h. On the other hand, for a patient presenting high creatinine clearance, such as 150 mL/min, the γ values were considerably lower for a shorter interval between the infusions (0.18, 0.25 and 0.54 mg/L) after an administration every 8, 12 and 24 h, respectively, Figs. 4b, d and f). Such augmented renal clearance is often observed for patients with traumatic injury [2]. The faster elimination of CMS suggests a treatment with shorter intervals between infusions, which, when optimized, leads to lower deviations of colistin plasma concentration to its target value. For the larger interval of 24 h (Fig. 4f) the higher clearance leads to a sharp decrease of colistin concentration between infusions, with a significant period with concentrations lower than 2 mg/L.

Here the focus was on the effect of creatinine clearance on the optimal doses, but it is straightforward to perform optimizations varying other parameters of the PK model. For instance, the model in [22] is a function of patient weight. In the present results we have considered a patient weighing 60 kg. Separate optimizations may be performed to select appropriate doses for other chosen patient weights. Moreover, the model also includes parameters for patients undergoing hemodialysis or continuous renal replacement therapy, and inclusion of these in the model enables optimizations for such patients as well.

Conclusion

The present work presented a new method to optimize dose selection based on a PK model, using optimal control. Here, the framework is used to guide the choice of the optimal colistin dose after an IV administration of CMS. This method can be applied for new molecules and other clinically relevant drugs, with clearer benefits for cases with a narrow therapeutic window such as colistin, as the optimization minimizes the deviations between a model result (here, the colistin plasma concentration) and a target value that may be specified based on the available knowledge about the drug at hand. The optimized doses show differences to currently used values, especially for a loading regimen which often uses multiple administrations with varying doses, leading to a faster attainment of the target concentration without toxicity issues. Moreover, the quantification of the optimal standard deviation γ allows an evaluation of potential benefits of lower administration intervals. Further studies with preclinical assays and clinical trials in humans would be relevant in order to confirm these results.