seventy patients of ages 7 months to 11 years were included , of whom 59 ( 84% ) were < 5 years of age . during the study ,
of these , one died within 15  min of admission , and two died within the first 24  h after admission .
one patient of age 2.5 years had severe neurologic sequelae 28 days after discharge , comprising spastic hemiparesis , blindness , and hearing impairment . at admission
this child had a count of 880,205 parasites/l and plasma plasmodium falciparum histidine - rich protein-2 concentration of 1,875  ng / ml , both indicating an extremely heavy parasite burden . at the 3-month follow - up ,
demographic , clinical , and laboratory characteristics of the study population are described in table 1 .
severe prostration , severe acidosis , convulsions , and severe anemia were the most common severity criteria .
eleven patients ( 16% ) presented with decompensated or compensated shock . in addition , three patients ( 4.3% ) had blood culture  proven septicemia ( unspeciated gram - negative rods , klebsiella pneumoniae , or staphylococcus aureus ) ; none of these presented with shock .
pretreatment with an oral antimalarial was reported with respect to 31 patients ( 6 with quinine , 3 with amodiaquine , 13 with sulfadoxine - pyrimethamine , 8 with artemether - lumefantrine , and 1 with amodiaquine followed by artemether - lumefantrine ) .
quinine within the 24  h preceding admission , the median ( range ) total dose being 16.1  mg / kg ( 10.153.7  mg / kg ) .
hypoglycemia was corrected with a 10% dextrose bolus at the time of admission and in those who developed hypoglycemia after admission ( eight patients ) .
peripheral blood asexual parasite counts after 24  h were negative in 12/66 ( 18% ) patients ( 24-h slide was not available for four patients , three of them as a result of death ) . the geometric mean ( 95% confidence interval ( ci ) ) parasite count after 24 h was 1,128 ( 5372,368 ) parasites/l in the rest of the study population ( n = 54 ) .
the overall geometric mean ( 95% ci ) fractional reduction in parasite counts at 24 h was 96% ( 9498% ) .
a total of 274 ars and dha postdose samples , randomly distributed over the first 12  h of the study , were included in the model .
a one - compartment disposition model for both ars and dha was adequate to describe the observed plasma concentration  time data ( table 2 ) .
all combinations of two - compartment disposition models displayed significant model misspecification despite significantly lower objective function values .
zero - order distribution / absorption from the injection site(s ) to the central compartment provided the best description of the data , but too few samples were collected during the absorption phase and the absorption rate was therefore fixed for an accurate estimation .
interindividual variability could be estimated reliably for ars and dha clearance and ars volume of distribution , showing correlation between ars clearance and volume .
when body weight was used as a fixed allometric function on all elimination clearance ( power of 0.75 ) and apparent volume of distribution ( power of 1 ) parameters , a significantly better model fit was observed ( objective function value = 14.6 ) .
dha clearance increased 10.2% per unit ( g / dl ) of decrease of hemoglobin .
interindividual variability in clearance decreased from 63.6% coefficient of variation to 55.3% coefficient of variation when this covariate was included , suggesting that it accounts for a limited but significant part of the observed variability .
the final model described the observed data well , with adequate goodness - of - fit diagnostics ( figure 1 ) ; eta - shrinkage : elimination clearance for ars 10.6% , central volume of distribution for ars 12.2% , elimination clearance for dha 6.81 ; epsilon - shrinkage : 49.4 and 22.5% for ars and dha , respectively ) .
a prediction - corrected visual predictive check of the final model resulted in no model misspecification and good simulation properties ( figure 2 ) .
the numerical predictive check ( n = 2,000 ) for ars resulted in 4.35% ( 95% ci 1.459.42 ) and 5.07% ( 95% ci 1.459.42 ) of observations above and below the 90% prediction interval , respectively ; for dha , these values were 2.97% ( 95% ci 1.799.52 ) and 5.95% ( 95% ci 1.798.92 ) , respectively .
ars and dha exposures after the standard 2.4  mg / kg dose were simulated at each body weight level ( 1,000 simulations each at 1-kg intervals from 6 to 25  kg ) using a uniform distribution of hemoglobin concentrations within the observed range ( 2.7213.6  g / dl ) to account for the observed covariate relationship ( figure 3a , b ) .
children with body weights between 6 and 10  kg showed a mean reduction of 20.4% ( p < 0.0001 ) in dha exposure as compared with children with body weights between 21 and 25  kg ( median ( 25th to 75th percentile ) exposure : 3,380 ( 2,1305,470 ) ng  h / ml in the 610  kg patients , 3,780 ( 2,4306,060 ) ng  h / ml in the 1115  kg patients , 4,100 ( 2,5706,590 ) ng  h / ml in the 1620  kg patients , and 4,240 ( 2,7006,840 ) ng  h / ml in the 2125  kg patients ) .
this suggests that smaller children need higher doses of ars to attain ars exposures similar to those in children with higher body weights .
using the final model , we evaluated various dosing regimens based on body weight bands .
the proposed regimen ( table 3 ) resulted in similar exposures in all weight bands after the first dose of i.m .
ars ( figure 3c , d ) . the same simulations were performed assuming a uniform distribution of body weights at different levels of hemoglobin ; these resulted in lower dha exposures with decreasing hemoglobin levels ( figure 4 ) .
there was no statistical difference between survivors and nonsurvivors with respect to total exposure ( ars p = 0.8060 , dha p = 0.4828 ) or maximum concentration ( ars p = 0.7655 , dha p = 0.6865 ) after the first dose .
similarly , an exposure  response relationship could not be established using nonlinear mixed - effects modeling ( nonmem ) in a time - to - event approach . however
, this might be a consequence of the relatively low number of deaths and the high proportion of pretreatments with different antimalarial drugs , doses , and administration routes .
seventy patients of ages 7 months to 11 years were included , of whom 59 ( 84% ) were < 5 years of age . during the study ,
of these , one died within 15  min of admission , and two died within the first 24  h after admission .
one patient of age 2.5 years had severe neurologic sequelae 28 days after discharge , comprising spastic hemiparesis , blindness , and hearing impairment . at admission
this child had a count of 880,205 parasites/l and plasma plasmodium falciparum histidine - rich protein-2 concentration of 1,875  ng / ml , both indicating an extremely heavy parasite burden . at the 3-month follow - up ,
demographic , clinical , and laboratory characteristics of the study population are described in table 1 .
severe prostration , severe acidosis , convulsions , and severe anemia were the most common severity criteria .
eleven patients ( 16% ) presented with decompensated or compensated shock . in addition , three patients ( 4.3% ) had blood culture  proven septicemia ( unspeciated gram - negative rods , klebsiella pneumoniae , or staphylococcus aureus ) ; none of these presented with shock .
pretreatment with an oral antimalarial was reported with respect to 31 patients ( 6 with quinine , 3 with amodiaquine , 13 with sulfadoxine - pyrimethamine , 8 with artemether - lumefantrine , and 1 with amodiaquine followed by artemether - lumefantrine ) .
quinine within the 24  h preceding admission , the median ( range ) total dose being 16.1  mg / kg ( 10.153.7  mg / kg ) .
hypoglycemia was corrected with a 10% dextrose bolus at the time of admission and in those who developed hypoglycemia after admission ( eight patients ) .
peripheral blood asexual parasite counts after 24  h were negative in 12/66 ( 18% ) patients ( 24-h slide was not available for four patients , three of them as a result of death ) . the geometric mean ( 95% confidence interval ( ci ) ) parasite count after 24 h was 1,128 ( 5372,368 ) parasites/l in the rest of the study population ( n = 54 ) .
the overall geometric mean ( 95% ci ) fractional reduction in parasite counts at 24 h was 96% ( 9498% ) .
a total of 274 ars and dha postdose samples , randomly distributed over the first 12  h of the study , were included in the model .
a one - compartment disposition model for both ars and dha was adequate to describe the observed plasma concentration  time data ( table 2 ) .
all combinations of two - compartment disposition models displayed significant model misspecification despite significantly lower objective function values .
zero - order distribution / absorption from the injection site(s ) to the central compartment provided the best description of the data , but too few samples were collected during the absorption phase and the absorption rate was therefore fixed for an accurate estimation .
interindividual variability could be estimated reliably for ars and dha clearance and ars volume of distribution , showing correlation between ars clearance and volume .
when body weight was used as a fixed allometric function on all elimination clearance ( power of 0.75 ) and apparent volume of distribution ( power of 1 ) parameters , a significantly better model fit was observed ( objective function value = 14.6 ) .
dha clearance increased 10.2% per unit ( g / dl ) of decrease of hemoglobin .
interindividual variability in clearance decreased from 63.6% coefficient of variation to 55.3% coefficient of variation when this covariate was included , suggesting that it accounts for a limited but significant part of the observed variability .
the final model described the observed data well , with adequate goodness - of - fit diagnostics ( figure 1 ) ; eta - shrinkage : elimination clearance for ars 10.6% , central volume of distribution for ars 12.2% , elimination clearance for dha 6.81 ; epsilon - shrinkage : 49.4 and 22.5% for ars and dha , respectively ) .
a prediction - corrected visual predictive check of the final model resulted in no model misspecification and good simulation properties ( figure 2 ) .
the numerical predictive check ( n = 2,000 ) for ars resulted in 4.35% ( 95% ci 1.459.42 ) and 5.07% ( 95% ci 1.459.42 ) of observations above and below the 90% prediction interval , respectively ; for dha , these values were 2.97% ( 95% ci 1.799.52 ) and 5.95% ( 95% ci 1.798.92 ) , respectively .
ars and dha exposures after the standard 2.4  mg / kg dose were simulated at each body weight level ( 1,000 simulations each at 1-kg intervals from 6 to 25  kg ) using a uniform distribution of hemoglobin concentrations within the observed range ( 2.7213.6  g / dl ) to account for the observed covariate relationship ( figure 3a , b ) .
children with body weights between 6 and 10  kg showed a mean reduction of 20.4% ( p < 0.0001 ) in dha exposure as compared with children with body weights between 21 and 25  kg ( median ( 25th to 75th percentile ) exposure : 3,380 ( 2,1305,470 ) ng  h / ml in the 610  kg patients , 3,780 ( 2,4306,060 ) ng  h / ml in the 1115  kg patients , 4,100 ( 2,5706,590 ) ng  h / ml in the 1620  kg patients , and 4,240 ( 2,7006,840 ) ng  h / ml in the 2125  kg patients ) .
this suggests that smaller children need higher doses of ars to attain ars exposures similar to those in children with higher body weights .
using the final model , we evaluated various dosing regimens based on body weight bands . the proposed regimen ( table 3 ) resulted in similar exposures in all weight bands after the first dose of i.m . ars ( figure 3c , d ) .
the same simulations were performed assuming a uniform distribution of body weights at different levels of hemoglobin ; these resulted in lower dha exposures with decreasing hemoglobin levels ( figure 4 ) .
there was no statistical difference between survivors and nonsurvivors with respect to total exposure ( ars p = 0.8060 , dha p = 0.4828 ) or maximum concentration ( ars p = 0.7655 , dha p = 0.6865 ) after the first dose .
similarly , an exposure  response relationship could not be established using nonlinear mixed - effects modeling ( nonmem ) in a time - to - event approach . however
, this might be a consequence of the relatively low number of deaths and the high proportion of pretreatments with different antimalarial drugs , doses , and administration routes .
optimal treatment strategies depend on detailed knowledge of the pharmacokinetic properties of drugs in the target population in which the drug is used .
age , disease status , and severity are all factors that may affect drug absorption , distribution , metabolism , and elimination .
the importance of pharmacokinetics in determining the therapeutic response is illustrated by the study of artemether , the first parenteral artemisinin derivative that was compared with quinine for effectiveness in the treatment of severe malaria in large clinical trials . in a meta - analysis of randomized trials in severe malaria ,
artemether significantly reduced mortality in southeast asian adults but did not do so in african children .
subsequent pharmacokinetic studies showed that the oil - based artemether was released slowly and erratically from the i.m .
injection site , and that this likely counterbalanced its pharmacodynamic advantages relative to quinine . despite having pharmacodynamic properties similar to those of artemether ,
ars is superior to artemether as a treatment because of its more favorable pharmacokinetic properties .
dosing regimens for children are often derived from studies in adults , and this practice has led to substantial underdosing of several antimalarials in children .
a pharmacokinetic study of sulfadoxine - pyrimethamine in african children with uncomplicated falciparum malaria showed that , with the usual dose of 25/1.25  mg / kg , the area under the concentration 
time curve in children 25 years of age was half that in adults . this may have caused failure of antimalarial treatment in small children , thereby contributing to the spread of resistance .
it has also been shown that piperaquine exposure levels are lower in small children if a standard body weight - based dose regimen is followed .
only limited data are available on the pharmacokinetics of parenteral , rectal , and oral ars in children .
ars was converted rapidly into dha , the ars elimination half - life being ~26  min .
this is in agreement with ars half - life values reported in other pharmacokinetic studies of i.m .
this is because the elimination rate of ars is limited by the rate of absorption from the i.m .
injection site ( i.e. ,  flip - flop  pharmacokinetics ) . the volume of distribution values , maximum concentration values , and area under the concentration  time curves of ars and dha are also comparable with the findings from a previous small and densely sampled pharmacokinetic study of i.m .
time curve of dha in this study was considerably lower than that reported in adult patients and healthy volunteers after i.v .
administration , this is still far greater than the in vitro  defined dha concentration value of 2.28  ng / ml that is required for 99% inhibition ( ic ) .
injection ( ~90% ) , fast absorption , and comparable estimates for ars and dha exposure support the use of i.m .
our study had only a relatively small number of patients , thereby limiting its power to detect covariate relationships . a parsimonious approach ( p < 0.001 for a covariate to be retained in the model ) was applied in order to avoid artifactual covariate relationships .
this has also been reported in other population pharmacokinetic studies of oral and rectal ars in pediatric and mixed adult  pediatric populations . in general
, physiological processes do not scale linearly with body weight ; consequently children with lower body weights will have higher body weight  normalized elimination clearance values .
dosing simulations using an exact dose of 2.4  mg / kg resulted in lower ars and dha exposures in small children as compared with those with body weights of 25  kg , suggesting the need for higher dosing in small children .
in addition , hemoglobin concentration was also a significant covariate , resulting in lower dha exposure in more anemic children .
severe anemia is a common presenting feature of severe malaria in young children . in severe malarial anemia there is usually substantial intravascular hemolysis , and
this is also the cause of the degradation of ars in hemolyzed plasma samples ex vivo .
hemoglobin concentration was found to be associated with reduced exposure levels of ars and dha , independent of body weight .
this supports the idea that young children with severe malaria , who are generally more anemic than older children , should receive adjusted higher doses .
underdosing in young children with severe malaria may have immediate adverse consequences with respect to outcomes .
a recent study carried out in cambodia has shown that oral ars at a dose of 6  mg / kg / day for 7 days resulted in a reduction in neutrophil counts ( to < 1.0    10/l ) and short - term neutropenia in 19% of the patients .
because oral ars has a bioavailability of ~80% , this corresponds to a total parenteral dose of 33.6  mg / kg , which would be reached only if parenteral treatment were continued for 14 doses . in the aquamat trial , the median ( interquartile range )
number of doses of parenteral treatment in surviving children was 3 ( 24 ) doses .
ars is currently the first - choice treatment for severe malaria , and a product produced according to good manufacturing practices , which is who prequalified and available . to facilitate implementation of an optimized dosing regimen in the treatment of severe malaria in african children , we defined a simplified weight - band  based dosing regimen based on the current population pharmacokinetics model , taking into account accuracy and practicality issues .
we considered 0.5  ml to be the minimum volume of prepared ars solution that can be measured accurately and administered with commonly available types of syringes . because the weight bands are smaller for children with body weights < 14  kg , we propose an incremental ars dose increase of 5 instead of 10  mg and a dilution of 10  mg / ml for i.m .
administration in this group , similar to the dosage regimen in the current study . in children with body weights of < 14  kg ,
doses should be split and administered in both thighs if injection volumes exceed 2  ml . in children of body weights 14  kg ,
larger injection volumes can be avoided by using a dilution of 20  mg / ml . binning of weight bands was also chosen in accordance with the currently available vial , which contains a dose of 60  mg , demarcating the upper limits of weight bands at 25  kg .
the dosing recommendations we propose do not extend beyond the weight ranges of the children included in this study .
no children with body weight < 6.5  kg were included ; therefore more studies of ars population pharmacokinetics are needed to evaluate dosing in these very small children and to support the current dosing recommendations .
more extensive sampling in the first 15  min and during the first 12  h after the dose could give more information than the current study provided . in conclusion ,
ars and dha exposures were lower in small children after i.m . administration of ars in severe malaria , warranting dose adaptation in this group .
this pharmacokinetic assessment of ars was part of the aquamat trial ( registration number isrctn 50258054 ) , a large multinational trial for which the results have been published elsewhere .
this substudy was conducted at teule hospital in muheza , tanzania , from may 2009 to july 2010 . except for the additional blood sampling ,
ethical approval was obtained from the tanzania medical research coordinating committee and the oxford tropical research ethics committee .
children 14 years with a clinical diagnosis of severe malaria confirmed by plasmodium lactate dehydrogenase ( pldh)-based rapid diagnostic test ( optimal , diamed , cressier , switzerland ) were recruited , and written informed consent was obtained from the respective parents or caregivers .
severe malaria was defined as the presence of at least one of the following : coma ( glasgow coma score 10 or blantyre coma score 2 in preverbal children ) , convulsions ( duration > 30  min or 2 episodes in the 24  h preceding admission ) , respiratory distress ( nasal alar flaring , costal indrawing / recession or use of accessory muscles , severe tachypnea ) or acidotic breathing (  deep  breathing ) , shock ( capillary refill time 3  sec and/or temperature gradient and/or systolic blood pressure < 70  mm hg ) , severe symptomatic anemia ( < 5  g / dl with respiratory distress ) , hypoglycemia ( <3
mmol / l ) , hemoglobinuria , severe jaundice , or , in older children , a convincing history of anuria or oliguria . patients who had received full treatment with parenteral quinine or a parenteral artemisinin derivative > 24  h before admission were excluded .
physical examination was carried out at admission , and a venous blood sample was taken for peripheral blood parasite count , quantitative assessment of plasma plasmodium falciparum histidine - rich protein-2 ( a measure of total body parasite burden ) , hiv serology ( sd bio - line hiv 1/2 3.0 ; standard diagnostics , kyonggi - do , korea / determine hiv-1/2 , abbott laboratories , il ) , blood culture , liver function tests ( aspartate aminotransferase , alanine transaminase , -glutamyltransferase , total bilirubin , creatinine , and urea , by reflotron , roche diagnostics , basel , switzerland ) , hematocrit , biochemistry , and acid  base parameters ( ec8 + cartridge for i - stat handheld blood analyzer , abbott laboratories , abbott park , il ) . hematocrit was reported from the i - stat reading or , when not available , measured by hemocue ( hemocue ab , ngelholm , sweden ) ( n = 5 ) . a neurologic examination was conducted at discharge , and repeated at day 28 in children who had not made a full neurologic recovery at discharge .
ars ( guilin pharmaceutical factory , guangxi , china ) was administered as an i.m .
injection ( 2.4  mg / kg ) shortly after admission , again at 12  h and 24  h , and daily thereafter .
the contents of each 60  mg vial of ars were dissolved in 1  ml 5% sodium bicarbonate ( provided with the drug ) and further diluted with 5  ml 5% dextrose ( final concentration of 10  mg / ml ) before administration as a deep i.m . injection into the anterolateral thigh .
dosing was based on the measured body weight of the patient , and injection volumes of > 23  ml were split and divided into two injections , one in each thigh .
when the patient was well enough to take oral medication , but after a minimum of 24  h ( two doses of i.m .
ars ) , a full 3-day course of oral artemether - lumefantrine ( co - artem ; novartis , basel , switzerland ) was given to complete the treatment .
vital signs and glucose were monitored at least every 6  h and also at any sign of deterioration in clinical condition .
a majority of the patients ( i.e. , other than those who were able to be orally fed ) received an infusion with dextrose 5% .
blood transfusion ( 20  ml / kg ) was given to children with hemoglobin concentrations of < 5  g / dl .
blood samples ( 1.5  ml ) were drawn from an indwelling catheter into prechilled fluoride oxalate tubes for ars and dha quantification before the first dose ( at baseline ) .
four subsequent samples were taken from each patient at preset random time points within the following time windows : 01 , 14 , 412 , and 1224  h after the first dose .
randomization of sampling times was done by computer - generated randomization ( stata version 12 ( stata , college station , tx ) . immediately after blood collection
, the blood samples for drug measurements were centrifuged at 4 c at 2,000 g for 7  min .
plasma samples ( 0.5  ml ) were stored at 80 c and shipped on dry ice to the moru department of clinical pharmacology , bangkok , thailand , for drug quantification .
ars drug content and quality were checked in vials taken randomly from the purchase lots ( see supplementary data online ) . drug analysis .
ars and dha plasma concentrations were measured using liquid chromatography  tandem mass - spectrometry .
quality control samples at low , middle , and high concentrations were analyzed in triplicate within each analytical batch to ensure accuracy and precision during the analysis .
the lower limit of quantification was set at 1.2  ng / ml for ars and 2.0  ng / ml for dha .
venous plasma concentrations were transformed into molar units and modeled as natural logarithms , using nonmem v.7 ( icon development solutions , ellicott city , md ) .
ars and dha were modeled simultaneously , using a drug  metabolite model with complete in vivo conversion of ars into dha ( for details , see supplementary data online ) . the first - order conditional estimation method with interaction
model selection was based on the objective function values computed by nonmem , goodness - of - fit graphical analysis , and physiological plausibility .
a p value of 0.05 was used in the forward step and a p value of 0.001 in the backward step to compensate for the relatively small population studied .
simulation - based diagnostics ( visual and numerical predictive checks ) and bootstrap diagnostics were used to evaluate the performance of the final model .
monte carlo simulations using the final model with the observed variability were performed for different body weights to obtain representative population estimates of the exposure levels during the first day of dosing ( area under the concentration 
time curve from time point 0 to 12  h ) after prospective dose regimens .
no drug exposure target is defined for parenteral ars ; therefore , for the purpose of arriving at a practical parenteral dosing regimen , different body weight  bins  were evaluated to ensure similar target exposures in all weight bands in agreement with the exposure in children of body weight 25  kg .
the same simulations were used to evaluate the effect of other significant covariates on drug exposure .
peripheral blood smears were taken at admission and after 24  h. reduction in parasite load over 24  h , survival , and severe neurologic sequelae were evaluated as part of the pharmacodynamic analysis .
the effects of ars and dha exposures on outcomes were investigated using a time - to - event analysis in nonmem , with predicted drug concentrations being used to modulate the hazard function in a traditional maximum effect ( emax ) relationship .
 parenteral artesunate ( ars ) is the drug of choice for the treatment of severe malaria . 
 pharmacokinetics data on intramuscular ars are limited with respect to the main treatment group that carries the highest mortality , namely , critically ill children with severe malaria . a population pharmacokinetic study of ars and dihydroartemisinin ( dha ) 
 was conducted from sparse sampling in 70 tanzanian children of ages 6 months to 11 years . 
 all the children had been admitted with severe falciparum malaria and were treated with intramuscular ars ( 2.4  mg / kg at 0 , 12 , and 24  h ) . 
 venous plasma concentration  
 time profiles were characterized using nonlinear mixed - effects modeling ( nonmem ) . 
 a one - compartment disposition model accurately described first - dose population pharmacokinetics of ars and dha . 
 body weight significantly affected clearance and apparent volume of distribution ( p < 0.001 ) , resulting in lower ars and dha exposure levels in smaller children . 
 an adapted dosing regimen including a practical dosing table per weight band is proposed for young children based on the pharmacokinetic model . 
Results
Clinical details
Population pharmacokinetic-pharmacodynamic analysis
Discussion
Methods
Author Contributions
Study Highlights
Supplementary Material