OptiMal-PK: an internet-based, user-friendly interface for the mathematical-based design of optimized anti-malarial treatment regimens
Malaria Journal volume 15, Article number: 344 (2016)
The search for highly effective anti-malarial therapies has gathered pace and recent years have seen a number of promising single and combined therapies reach the late stages of development. A key drug development challenge is the need for early assessment of the clinical utility of new drug leads as it is often unclear for developers whether efforts should be focused on efficacy or metabolic stability/exposure or indeed whether the continuation of iterative QSAR (quantitative structure–activity and relationships) cycles of medicinal chemistry and biological testing will translate to improved clinical efficacy. Pharmacokinetic and pharmacodynamic (PK/PD)-based measurements available from in vitro studies can be used for such clinical predictions. However, these predictions often require bespoke mathematical PK/PD modelling expertise and are normally performed after candidate development and, therefore, not during the pre-clinical development phase when such decisions need to be made.
An internet-based tool has been developed using STELLA® software. The tool simulates multiple differential equations that describe anti-malarial PK/PD relationships where the user can easily input PK/PD parameters. The tool utilizes a simple stop-light system to indicate the efficacy of each combination of parameters. This tool, called OptiMal-PK, additionally allows for the investigation of the effect of drug combinations with known or custom compounds.
The results of simulations obtained from OptiMal-PK were compared to a previously published and validated mathematical model on which this tool is based. The tool has also been used to simulate the PK/PD relationship for a number of existing anti-malarial drugs in single or combined treatment. Simulations were predictive of the published clinical parasitological clearance activities for these existing therapies.
OptiMal-PK is designed to be implemented by medicinal chemists and pharmacologists during the pre-clinical anti-malarial drug development phase to explore the impact of different PK/PD parameters upon the predicted clinical activity of any new compound. It can help investigators to identify which pharmacological features of a compound are most important to the clinical performance of a new chemical entity and how partner drugs could potentially improve the activity of existing therapies.
Anti-malarial drug development has entered a new era where the community can boast to be in possession of thousands of anti-malarial active compounds with in vitro growth inhibition IC50 values in the sub-micro molar range [1–3]. The increased activity in anti-malarial drug development within the context of the malaria eradication agenda, has resulted in the recommendation of a series of target candidate profiles (TCPs) and target product profiles (TPPs). The TCPs include; fast-parasite clearance drug profiles (TCP-1), long-duration of drug action profiles (TCP-2), liver-stage (including hypnozoites) and sexual-stage acting drug profiles (TCP-3), and chemoprotection profiles (TCP-4). The ambition is to combine drugs with one or more TCP attributes to either (i) generate a single exposure radical cure and prophylaxis (TPP-1, SERCaP) treatment- or (ii) generate a single exposure chemoprotection treatment (TPP-2, SECA) .
The challenge for the drug discovery community is to prioritize anti-malarial hits and leads for subsequent pre-clinical development. Lead optimization studies are used to enhance the effectiveness of the most promising compounds, however it is not always apparent whether a drug’s progress is directly in line with desired TCPs or TPPs. For example, despite the availability of in vitro ADME data (absorption, disposition, metabolism and excretion) or in vivo DMPK (pharmacokinetics) studies are unable to show whether a drug’s less desirable pharmacokinetic attributes can be counter-balanced by other more desirable features such as potency. Slow parasite killing kinetics (measured as the parasite reduction ratio, PRR) are clearly undesirable but can be overcome if the drug has an extensive therapeutic half-life. Similarly, a rapid-killing drug can be effective even with a short therapeutic half-life. Therefore to be confident of a drug’s ultimate clinical performance, extensive in vivo experimentation or pharmacokinetic–pharmacodynamic (PK/PD) modelling is required.
The major factors influencing the overall clinical performance of a drug include; (i) the drug’s therapeutic window describing the maximum dose tolerated by patients and the minimum dose required to produce a therapeutic effect, (ii) the drug’s PK properties describing its clearance and volume of distribution (ultimately reflected in the elimination half-life) and (iii) the drug’s PD properties, which include intrinsic potency (as determined in static growth inhibition assays) and time-dependent kill dynamics (in the case of malaria the PRR which is a fixed characteristic of the biological target/pathway affected by the anti-malarial drug class).
The complexity of interacting factors that affect the overall performance of a drug make it a challenge during lead optimization and pre-clinical drug development to fully understand when appropriate PK/PD optimization has been (or can ever be) achieved. Pharmacological models have the potential to bridge this knowledge gap and provide quantitative predictions of drug performance. In silico PK/PD models have historically been used to investigate anti-malarial monotherapies [5–10] and more recently combination therapies [11–14]. However, these studies do not offer a platform for non-specialists to explore the dynamic consequences of different pharmacological drug properties in order to support decision-making in the discovery process.
Presented here, is a user-friendly, internet-based platform “OptiMal-PK” that provides performance data that can support pre-clinical anti-malarial drug development initiatives. OptiMal-PK allows the user to input drug PK/PD parameters and utilizes a simple traffic-light display system to indicate the minimum number of days the drug should be administered to completely clear parasites form a patient. OptiMal-PK also allows for the optional addition of a partner drug from either a list of currently available drugs or a second user-defined drug. This allows the users to investigate how drug combinations could help shorten treatment times or alter regimens if needed. This tool serves two major purposes (i) to develop a better understanding of the complex interaction of anti-malarial drug PK and PD properties and (ii) to identify and quantify the key features of any drug under development that drive efficacy and therefore represent a focus for improvement in the drug development cycle.
The OptiMal-PK tool is an open access tool available online  which has been designed to be user friendly. Users can easily alter PK and PD values and generate instantaneous answers for any combination of PK/PD parameters with or without the addition of a partner drug. The model design, validation and examples of the uses of OptiMal-PK are described in the sections that follow.
The software used for the construction of this model is STELLA® (Systems Thinking for Education and Research) which is a package that allows for the dynamic visualization and communication of complex differential equation systems . The model itself is hosted on forio.com, a platform for enabling simulations, data exploration and data analytics.
The model, based on Winter and Hastings  as previously described, connects a PK component describing the drug concentration in the systemic circulation over time and a PD component defining the effect of the drug concentration at any given time on the total parasite count. This allows the user to input the PK/PD parameters of their experimental drug and roughly estimate its predictive clinical performance in patients.
The PK model within OptiMal-PK assumes the drug is given orally, and because it is based on a previous model, it also assumes first-order absorption kinetics and distribution into one PK compartment. The model was hence only defined by two standard differential equations:
where X 1 represents the mass of drug in the gut at time t and is at its maximum when a dose is administered. X2 is the mass of drug in the blood at any given time, it increases as the drug is absorbed from the gut at rate k a and decreases as the drug is eliminated at a rate k e . Blood concentration levels are used throughout the model instead of plasma as the model assumes that drug elimination is driven by total clearance where protein binding does not alter overall exposure.
Setting the dose in the interface will then change the value of X1 accordingly to match the required dosage. Initial value for X2 is set to be 0 by default as this represents initial drug amount in blood circulation as can be seen in Fig. 1.
The drug concentration (C) at any given time (t) is determined by dividing the drug mass in blood (X 2 ) by the volume of distribution (V d ).
This relationship can be simplified into terms for the apparent elimination half-life and the clearance of the drug; these parameters are required as user inputs are converted by the program into the micro-parameters as follows:
The pharmacodynamic component of the model takes into account the parasite multiplication rate per parasite life cycle (PMR = parasite multiplication rate/48 h), the concentration of drug in blood at any specific time, the PRR (the maximum Parasite Reduction Ratio occurring in 48 h [11, 17–19]), C50 of the drug, defined as the concentration of the drug required to achieve half of the maximal parasite reduction rate and the Hill slope, which defines the slope parameter of the drug’s C50. The C50 value is the concentration required to achieve half the maximal PRR which is defined by the user. This value is ideally extracted from in vitro PRR studies , and if not available, the IC50 value derived from traditional static 48 or 72 h assays can be used as a suitable surrogate in OptiMal-PK . Utilizing the IC50 values from static models have resulted in realistic predictions that match clinical results obtained from various drugs as previously shown .
The PD model relates drug concentration C to its effect on parasite viability. The concentration and time-dependent killing function f(C) for each drug is described using the standard Michaelis–Menten equation, i.e.
where V max is the maximal drug-killing rate, n is the slope of the dose response curve, and C 50 is the concentration of drug at which 50 % of the maximal PRR occurs.
The change in the number of parasites P over time t can be found with the standard differential equation.
where (a) is the parasite growth rate determined by the user-defined parasite multiplication rate (PMR). PMR is set by default to ten based on previous evidence , but could be altered by the user to reflect the different PMR values that have been reported in different regions .
The model additionally calculates the minimum parasiticidal concentration (MPC), a term often used to describe the minimum concentration needed to achieve a net decrease in parasite count over time. MPC is directly calculated from the drug concentration (C) that results in a net reduction in parasite load (e.g. rate of parasite kill (f(C)) > PMR, Eq. 8).
The model’s work-flow follows the schematic shown in Fig. 1.
Parameter values for all built in partner drugs supplied in the table (see OptiMal-PK website) were taken from the paper on which OptimMal-PK is based  except for atovaquone where the PK parameters were taken from , the IC50 data from  and the PRR values obtained from clinical data  which matches the in vivo PRR of drugs with similar mode of action .
Stage specificity within OptiMal-PK. A recent paper by Hodel et al.  investigated the accuracy of this methodology by modelling drugs with long and short half-lives, with and without stage specificity. The study found stage-specificity was only important for short half-life drugs with stage-specific killing (e.g. the artemisinins) because, depending on the timing of treatment, parasites might be in highly drug-tolerant stages or in much less tolerant stages. When modelling drugs with very short half-lives and stage-specific killing users should note that their results could vary but in all other instances, the model was shown to be very robust without the addition of stage-specificity.
OptiMal-PK has been designed to be a user-friendly tool for non-specialists allowing the potential activity of a developing drug or new drug combination to be assessed but care should be taken into over-interpretation of any results due to the previous complexities that have been necessarily simplified in the model.
OptiMal-PK is an internet-based program that utilizes STELLA® (Systems Thinking for Education and Research) software. The user-friendly interface allows scientists with limited knowledge of mathematical modelling and/or the PK-PD dynamics of anti-malarial drugs to easily assess the predicted clinical performance of any anti-malarial drug/drug combination and explore the consequence of modifying key pharmacological properties of the drug.
The software works in two simple steps. First the user inputs the parameters describing drug PK, PD and initial parasite load with user-friendly dials, shown in Fig. 2a, and then runs the simulation by pressing the green “run” button (Fig. 2b). The OptiMal-PK simulation calculates the minimum treatment time (days) required to clear an infection and displays the results using a simple traffic light system. A green light indicates an adequate number of treatment days, i.e. the number of days of treatment is sufficient to achieve a cure (defined as reaching a total parasite count of less than one, where the simulation stops and indicates achievement of cure). While a red light indicates that the number of days of treatment administration is insufficient to achieve a cure.
Results and discussion
The OptiMal-PK model was based on a validated mechanistic PK/PD model implemented in R  and previously published by Winter and Hastings  To validate the methodology of the OptiMal-PK model, the drug concentration and parasite number for the first 30 days post-treatment were simulated using both the OptiMal-PK and R models. The results of both simulations for the seven anti-malarial monotherapies previously calibrated and validated in  were compared in house and found to be largely consistent.
Small differences in the absolute values appear to be the result of the methods used by the modelling platforms in integrating differential equations. Stella® uses the Euler’s method to solve differential equations in 1 h time steps while the Winter and Hastings model uses an algebraic solution for the differential equations. However, parasite levels at the selected time points were consistently within ±10 % across the two models.
The assumption of a one compartment model has been used previously to simulate three widely used anti-malarial drug combinations (artemether–lumefantrine, artesunate–mefloquine and DHA–piperaquine) . However it is noted that subsequent pharmacological models have shown that using a two- or three-compartment PK model, or allowing for conversion to an active metabolite are more appropriate for some anti-malarials (see for example [12, 29]). The purpose of the PK model within OptiMal-PK is to inform the development of a drug candidate so the assumption of a one-compartment model is a reasonable and pragmatic first step.
It is important to note that the model makes a large number of other simplifications similar to previous models. For example, OptiMal-PK does not currently include bioavailability, dose dependency, active metabolites and intercessional/intersubjective variability which are known for many drugs. However, the simple nature of the interface will allow the user to probe the effects of these variables and assess their effect. For instance, the effect of bioavailability can be easily explored by reducing or increasing the drug dosage, or by altering the clearance values. Active metabolites can be modelled by adding an extra drug with parameters that match that of the metabolite and whose absorption rate is equal to the conversion rate to the metabolite.
For combination therapy simulations, a model of dominant killing was utilised, where the rate of kill is equal to that of the drug inferring the fastest kill rate at each time point as previously described in .
A case study—Is OZ439 sufficiently effective for delivery as a one dose therapy and how it could it be improved?
OZ439 is a peroxide anti-malarial drug that has been designed to produce a single oral dose cure for malaria . Unlike other peroxides, OZ439 showed excellent activity in Plasmodium berghei-infected mice producing a complete cure after a single 30 mg/kg dose administered orally. The clinical PK data for the drug were described in a later publication  and showed superior PK properties when compared to similar drugs within the same class, with an elimination half-life that exceeds 24 h due partially to a relatively slow clearance. However, the question remains of whether OZ439 could truly deliver a single dose cure (defined as the reduction of total number of parasites to less than one) when administered as a single therapy in clinical studies, although it is accepted that the drug will be eventually deployed as a combination. OptiMal-PK was utilised to predict the clinical performance of OZ439 in future human studies. A number of simulations were executed assuming different dosages and feeding states and then allowed the program to assess the minimum number of days needed to achieve total elimination of parasites in the subjects under each scenario.
The OptiMalPK model calculated minimum number of days needed for treatment using the PK and PD parameters that are presented in Table 1 (derived from [30–32]) In this example, different doses of OZ439 were compared in a capsule and oral dispersion forms on empty (lower exposure) or full (higher exposure) stomachs based on clinical data. The results demonstrate that a single monotherapy dose of 800 or 1200 mg (as a capsule, it is important to note that the disposition characteristics of OZ439 capsules differ significantly from a drug dispersion formulation, see below) would be insufficient for the successful elimination of parasites in all patients and would result in recrudescence within 9 days after the treatment. The model predicts that 1200 mg OZ439 capsule would need at least 2 days of treatment to completely cure a malaria infection, as defined by depletion of parasite total number to one.
However, the simulations show that a single 800 mg dose of the oral dispersion formulation in fed or fasted subjects would lead to a cure following a single administration; administering the drug in the capsule form, not the oral dispersion formulation, seems to fail a 1 day treatment regimen according to the model. Nevertheless, even for the superior oral dispersion formulation, if metabolism is slightly accelerated (decreasing half-life by 20 %) in a proportion of the target population, the single dose for the oral dispersion would be rendered insufficient in either fed or fasted subjects. These data suggest that the single dose regimen is very unforgiving which could have very serious consequences in field deployment in large patient populations if OZ439 were to be deployed alone. Additionally a modest two to threefold decrease in potency (increase in IC50 from x–y; i.e. emergence of moderate resistance) would also make the single dose insufficient for an absolute cure (Fig. 3). OptiMal-PK highlights the sensitivity of single dose OZ439 to emergence of low level parasite resistance and moderate population PK variations. Given the large variability of PK parameters shown in Table 1 (30–115 %) the simulations here show that a single dose of OZ439 dispersion would leave a significant proportion of subjects at risk of incomplete elimination of malaria parasites, treatment failures and parasite recrudescence and increased pressure for resistance development.
Case study 2—effect of resistance on drug activity and mitigating strategies
Chloroquine (CQ) was for many decades a primary chemotherapeutic agent in anti-malarial therapy. Resistance to the 4-aminoquinolines was initially reported in Southeast Asia and South America and subsequently spread to almost every malaria endemic region in the world . CQ was used as an example on how drug resistance would affect the outcome of treatment. The original chloroquine parameters are all based on those previously reported as described earlier . The number of days it takes to achieve a cure was investigated as drug resistance increases (IC50 was elevated in 40 % increments, Fig. 4). Using CQ sensitive IC50 values (e.g. ≤62.7 nM ), the model predicts successful outcome with a 3-day treatment, consistent with published clinical data , and consistent with WHO guidelines for CQ . When parasite CQ sensitivity values are increased to resistant values e.g. generally considered ≥100 nM [37, 38], the model predicts a requirement of 4 or more days of treatment, depending on the fold increase of CQ C50 (Fig. 4), all PK parameters derived from  and PD parameters from . These data are in agreement with clinical trials performed in CQ resistant regions in Afghanistan where 60 % of the population failed a 5-day treatment of CQ .
Parameters that could be altered to reverse the effect of increasing the IC50 of CQ were then investigated. Assuming a 4.6-fold increase of IC50 from 62.7 to 288.4 nM the drug will not be able to deliver a cure within 14 days, altering the dosing regimen from once daily dosing to twice daily brings the number of days needed to achieve a cure to five to 6 days according to the model, while artificially increasing the half-life by 40 % would reduce the predicted number of days needed to achieve a cure down to seven to 8 days. Increasing the PRR by tenfold from 1000 only improves the activity of CQ slightly according to the model. Although these changes are purely an academic exercise it highlights the key parameters that need to be considered in the discovery process where candidate selection could take into account these pharmacological features. Importantly, even if it were possible to achieve radical changes in the drug properties, such as increasing its PRR by tenfold or its half-life by 40 %, these changes are unable to mitigate resistance resulting from a 4.6-fold increase in the CQ IC50, emphasizing the need to avoid resistance development in the first case by appropriate use of forgiving dosage regimens.
Case study 3—studying the potential of long-duration dihydroorotate dehydrogenase (DHODH) inhibitor (DSM265) for treatment of malaria
OptiMal-PK is used here to study the therapeutic potential of a drug, DMS265, with slow elimination kinetics and killing dynamics. Potency and clinical PK properties were used as previously reported  and shown in Table 2 (PRR value was estimated to be 17 based on data shown in ).
The simulations show that using the parameters in Table 2, at least four doses would be needed to achieve complete elimination of parasites in the circulation of patients. However, if metabolism is slowed down by a factor of only 10–20 % (by increasing half-life from 130 to 170 h) the drug would be capable of achieving a cure in one to two doses. These data, therefore, suggest that in back-up programmes aimed at delivering second-generation DHODH inhibitors, a focus on improving the therapeutic half-life of the drug may be more rewarding than a focus on potency alone.
Like OZ439, albeit with a completely different mode of action and different PK/PD properties, DSM 265 seems to be on the cusp of delivering a single dose cure. As with OZ439, the laudable aim of achieving a single dose cure, needs to take into account inherent parasite and patient variability at the population level, as the models suggest that these drugs may end up being deployed in very unforgiving dosage regimens, setting them up for premature failure due to either inadequate clinical performance in some patients or the early development of parasite resistance, something that is readily achieved with DSM265 in laboratory settings under moderate drug pressure.
OptiMal-PK is an open access, internet-based tool designed for the wider anti-malarial drug discovery community. It allows non-specialists to identify which pharmacological features of a compound are most important to the clinical performance of a new chemical entity or new drug combination. OptiMal-PK should therefore serve as an aid to critical decision-making during anti-malarial drug development programmes.
absorption, disposition, metabolism, excretion and toxicity
disposition, metabolism and pharmacokinetics
minimum parasitic concentration
parasite multiplication ratio
parasite reduction ratio
minimum parasitic concentration
- k a :
- k e :
quantitative structure–activity and relationships
- T1/2 :
target candidate profile
target product profile
- Vmax :
maximal kill rate
- Vd :
volume of distribution
Gamo FJ, Sanz LM, Vidal J, de Cozar C, Alvarez E, Lavandera JL, et al. Thousands of chemical starting points for antimalarial lead identification. Nature. 2010;465:305–10.
Guiguemde WA, Shelat AA, Bouck D, Duffy S, Crowther GJ, Davis PH, et al. Chemical genetics of Plasmodium falciparum. Nature. 2010;465:311–5.
Plouffe D, Brinker A, McNamara C, Henson K, Kato N, Kuhen K, et al. In silico activity profiling reveals the mechanism of action of antimalarials discovered in a high-throughput screen. Proc Natl Acad Sci USA. 2008;105:9059–64.
Burrows JN, van Huijsduijnen RH, Mohrle JJ, Oeuvray C, Wells TN. Designing the next generation of medicines for malaria control and eradication. Malar J. 2013;12:187.
Austin DJ, White NJ, Anderson R. The dynamics of drug action on the within-host population growth of infectious agents: melding pharmacokinetics with pathogen population dynamics. J Theor Biol. 1998;194:313–39.
Hoshen MB, Na-Bangchang K, Stein WD, Ginsburg H. Mathematical modelling of the chemotherapy of Plasmodium falciparum malaria with artesunate: postulation of ‘dormancy’, a partial cytostatic effect of the drug, and its implication for treatment regimens. Parasitology. 2000;121:237–46.
Hoshen MB, Stein WD, Ginsburg H. Modelling the chloroquine chemotherapy of falciparum malaria: the value of spacing a split dose. Parasitology. 1998;116:407–16.
Hoshen MB, Stein WD, Ginsburg H. Pharmacokinetic-pharmacodynamic modelling of the anti-malarial activity of mefloquine. Parasitology. 2001;123:337–46.
Saralamba S, Pan-Ngum W, Maude RJ, Leea SJ, Tarning J, Lindegardh N, et al. Intrahost modeling of artemisinin resistance in Plasmodium falciparum. Proc Natl Acad Sci USA. 2011;108:397–402.
Simpson JA, Watkins ER, Price RN, Aarons L, Kyle DE, White NJ. Mefloquine pharmacokinetic–pharmacodynamic models: implications for dosing and resistance. Antimicrob Agents Chemother. 2000;44:3414–24.
Winter K, Hastings IM. Development, evaluation, and application of an in silico model for antimalarial drug treatment and failure. Antimicrob Agents Chemother. 2011;55:3380–92.
Kay K, Hastings IM. Improving pharmacokinetic–pharmacodynamic modeling to investigate anti-infective chemotherapy with application to the current generation of antimalarial drugs. PLoS Comput Biol. 2013;9:e1003151.
Hodel EM, Kay K, Hayes DJ, Terlouw DJ, Hastings IM. Optimizing the programmatic deployment of the anti-malarials artemether–lumefantrine and dihydroartemisinin–piperaquine using pharmacological modelling. Malar J. 2014;13:138.
Zaloumis S, Humberstone A, Charman SA, Price RN, Moehrle J, Gamo-Benito J, et al. Assessing the utility of an anti-malarial pharmacokinetic–pharmacodynamic model for aiding drug clinical development. Malar J. 2012;11:303.
OptiMal-PK. 2016. http://optimalpk.lstmed.ac.uk. Accessed 24 May 2016.
Systems I. 2006. http://www.iseesystems.com. Accessed 24 May 2016.
Edwards G, Krishna S. Pharmacokinetic and pharmacodynamic issues in the treatment of parasitic infections. Eur J Clin Microbiol. 2004;23:233–42.
Pukrittayakamee S, Chantra A, Simpson JA, Vanijanonta S, Clemens R, Looareesuwan S, et al. Therapeutic responses to different antimalarial drugs in vivax malaria. Antimicrob Agents Chemother. 2000;44:1680–5.
White NJ. Assessment of the pharmacodynamic properties of antimalarial drugs in vivo. Antimicrob Agents Chemother. 1997;41:1413–22.
Sanz LM, Crespo B, De-Cozar C, Ding XC, Llergo JL, Burrows JN, et al. P. falciparum in vitro killing rates allow to discriminate between different antimalarial mode-of-action. PLoS One. 2012;7:e30949.
Smilkstein M, Sriwilaijaroen N, Kelly JX, Wilairat P, Riscoe M. Simple and inexpensive fluorescence-based technique for high-throughput antimalarial drug screening. Antimicrob Agents Chemother. 2004;48:1803–6.
Chotivanich K, Udomsangpetch R, Simpson JA, Newton P, Pukrittayakamee S, Looareesuwan S, et al. Parasite multiplication potential and the severity of falciparum malaria. J Infect Dis. 2000;181:1206–9.
Deans AM, Lyke KE, Thera MA, Plowe CV, Kone A, Doumbo OK, et al. Low multiplication rates of African Plasmodium falciparum isolates and lack of association of multiplication rate and red blood cell selectivity with malaria virulence. Am J Trop Med Hyg. 2006;74:554–63.
Hussein Z, Eaves J, Hutchinson DB, Canfield CJ. Population pharmacokinetics of atovaquone in patients with acute malaria caused by Plasmodium falciparum. Clin Pharmacol Ther. 1997;61:518–30.
Lacy MD, Maguire JD, Barcus MJ, Ling J, Bangs MJ, Gramzinski R, et al. Atovaquone/proguanil therapy for Plasmodium falciparum and Plasmodium vivax malaria in Indonesians who lack clinical immunity. Clin Infect Dis. 2002;35:e92–5.
Phillips MA, Lotharius J, Marsh K, White J, Dayan A, White KL, et al. A long-duration dihydroorotate dehydrogenase inhibitor (DSM265) for prevention and treatment of malaria. Sci Transl Med. 2015;7:296ra111.
Hodel EM, Kay K, Hastings IM. Incorporating stage-specific drug action into pharmacological modeling of antimalarial drug treatment. Antimicrob Agents Chemother. 2016;60:2747–56.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. http://www.R-project.org/. Accessed 24 May 2016.
Kay K, Hodel EM, Hastings IM. Altering drug regimens to restore and enhance antimalarial drug effectiveness. Antimicrob Agents Chemother. 2015;59:6419–27.
Charman SA, Arbe-Barnes S, Bathurst IC, Brun R, Campbell M, Charman WN, et al. Synthetic ozonide drug candidate OZ439 offers new hope for a single-dose cure of uncomplicated malaria. Proc Natl Acad Sci USA. 2011;108:4400–5.
Moehrle JJ, Duparc S, Siethoff C, van Giersbergen PL, Craft JC, Arbe-Barnes S, et al. First-in-man safety and pharmacokinetics of synthetic ozonide OZ439 demonstrates an improved exposure profile relative to other peroxide antimalarials. Br J Clin Pharmacol. 2013;75:524–37.
Phyo AP, Jittamala P, Nosten FH, Pukrittayakamee S, Imwong M, White NJ, et al. Antimalarial activity of artefenomel (OZ439), a novel synthetic antimalarial endoperoxide, in patients with Plasmodium falciparum and Plasmodium vivax malaria: an open-label phase 2 trial. Lancet Infect Dis. 2016;16:61–9.
Ridley RG. Medical need, scientific opportunity and the drive for antimalarial drugs. Nature. 2002;415:686–93.
Mayxay M, Barends M, Brockman A, Jaidee A, Nair S, Sudimack D, et al. In vitro antimalarial drug susceptibility and pfcrt mutation among fresh Plasmodium falciparum isolates from the Lao PDR (Laos). Am J Trop Med Hyg. 2007;76:245–50.
Trape J-F. The public health impact of chloroquine resistance in Africa. Am J Trop Med Hyg. 2001;64(1 suppl):12–7.
Lalloo DG, Shingadia D, Pasvol G, Chiodini PL, Whitty CJ, Beeching NJ, et al. UK malaria treatment guidelines. J Infect. 2007;54:111–21.
Brasseur P, Kouamouo J, Moyou-Somo R, Druilhe P. Multi-drug resistant falciparum malaria in Cameroon in 1987–1988. I. Stable figures of prevalence of chloroquine- and quinine-resistant isolates in the original foci. Am J Trop Med Hyg. 1992;46:1–7.
Cremer G, Basco LK, Le Bras J, Camus D, Slomianny C. Plasmodium falciparum: detection of P-glycoprotein in chloroquine-susceptible and chloroquine-resistant clones and isolates. Exp Parasitol. 1995;81:1–8.
White NJ. Antimalarial pharmacokinetics and treatment regimens. Br J Clin Pharmacol. 1992;34:1–10.
Howard N, Durrani N, Sanda S, Beshir K, Hallett R, Rowland M. Clinical trial of extended-dose chloroquine for treatment of resistant falciparum malaria among Afghan refugees in Pakistan. Malar J. 2011;10:171.
GA model design and construction. KK performed model validation. GB conceived the project. GA, KK, SW and GB contributed to the writing of the manuscript. All authors read and approved the final version of the manuscript.
The Authors would like to thank forio (Training Simulations & Predictive Analytics Applications) for hosting OptiMalPK on its platform.
The authors declare that they have no competing interests.
Authors wish to acknowledge support from the Medical Research Council (MRC, MC_PC_12017, MC_PC_13069, MC_PC_14111 and MR/L000644/1).
About this article
Cite this article
Aljayyoussi, G., Kay, K., Ward, S.A. et al. OptiMal-PK: an internet-based, user-friendly interface for the mathematical-based design of optimized anti-malarial treatment regimens. Malar J 15, 344 (2016). https://doi.org/10.1186/s12936-016-1401-8