Skip to main content

A population pharmacokinetic model is beneficial in quantifying hair concentrations of ritonavir-boosted atazanavir: a study of HIV-infected Zimbabwean adolescents



Adolescents experience higher levels of non-adherence to HIV treatment. Drug concentration in hair promises to be reliable for assessing exposure to antiretroviral (ARV) drugs. Pharmacokinetic modelling can explore utility of drug in hair. We aimed at developing and validating a pharmacokinetic model based on atazanavir/ritonavir (ATV/r) in hair and identify factors associated with variabilities in hair accumulation.


We based the study on secondary data analysis whereby data from a previous study on Zimbabwean adolescents which collected hair samples at enrolment and 3 months follow-up was used in model development. We performed model development in NONMEM (version 7.3) ADVAN 13.


There is 16% / 18% of the respective ATV/r in hair as a ratio of steady-state trough plasma concentrations. At follow-up, we estimated an increase of 30% /42% of respective ATV/r in hair. We associated a unit increase in adherence score with 2% increase in hair concentration both ATV/r. Thinner participants had 54% higher while overweight had 21% lower atazanavir in hair compared to normal weight participants. Adolescents receiving care from fellow siblings had atazanavir in hair at least 54% less compared to other forms of care.


The determinants of increased ATV/r concentrations in hair found in our analysis are monitoring at follow up event, body mass index, and caregiver status. Measuring drug concentration in hair is feasibly accomplished and could be more accurate for monitoring ARV drugs exposure.

Peer Review reports


About 36.9 million people were living with Human Immunodeficiency Virus (HIV) worldwide in 2017. Of these, approximately 3 million were children and adolescents under 20 years of age [1]. Zimbabwe has a prevalence of 13.3%, with 1.3 million people living with HIV including 77,000 children and adolescents [2]. Poor adherence to treatment leads to sub-optimal drug exposure limiting treatment efficacy [3, 4]. It has been estimated that between 20 to 50% of adolescents experience adherence-related antiretroviral (ARV) drug treatment failure [5,6,7,8,9,10,11].

It is desirable to have a routine assessment of adherence and exposure to ARV drugs available for use by healthcare providers. The methods used for assessing adherence and exposure to ARV drugs include self-reported missed doses, monitoring pharmacy refills and conducting pill counts, use of electronic monitoring devices, measuring ARV concentration in plasma or hair [12,13,14]. Quantifying ARV drugs in hair provides information of both steady-state pharmacokinetics and long term adherence and has shown to predict well the relationship between drug exposure and treatment outcomes when compared to other approaches [15,16,17,18,19,20,21].

Some suggest that hair uptake most external substances or their metabolites from the systemic circulation through the hair bulb blood supply by passive diffusion from blood into growing hair cells at the base of the follicle and then bound in the hair shaft [22,23,24,25]. Once the drug accumulates into the growing hair, we can detect it long after elimination from the systematic circulation, unlike in conventional biological samples such as blood and urine [26,27,28,29,30]. The scalp hair fibre grows at an average rate of 0.5 to 1.5 cm per month [31]. Thus the amount of drug in hair is constantly increasing until the next hair cut or when all the drug is removed from the systematic circulation.

In Zimbabwe and other resources limited settings, pharmacokinetic (PK) modelling applied focused primarily on systemic exposure to ARV drugs. It based the models used on data generated by quantifying drugs mostly from single time-point plasma samples [32,33,34,35]. Hypothetically hair PK parameters can provide additional information about the patient’s drug exposure overtime, hence the need to determine and apply hair PK for prediction of drug amount in the hair in relation to exposure in plasma.

Using single time-point plasma samples is considered unreliable when assessing drug exposure in populations at risk of non-adherence. Measuring ARV concentrations in hair promises to be more accurate and feasibly accomplished. However, there are very few studies which analysed drug concentration measured in hair using PK approaches. The aim of the work reported in this paper was to develop and validate a population PK model for ATV/r concentrations in hair and explore factors associated with increased or reduced concentrations assuming a direct relationship between ratio in plasma and hair.


Source of data

The study used secondary data from a previously published study conducted in Zimbabwe [15]. They collected the data between January 2015 and May 2016. It comprised 50 adolescents aged between 10 to 18 years and were on ATV/r (300/100 mg) based 2nd line HIV treatment for at least 6 months. They enrolled these participants at a public health hospital in Harare, Zimbabwe, and randomized to either adherence intervention or standard of care arms. They excluded participants if they were on anti-TB treatment, did not prefer home be followed-ups, had viral load < 1000 copies/ml within the previous 2 months, or were on ATV/r as 1st line treatment. The goal of the primary study was to test the impact of a home-based modified directly administered adherence intervention on virologic outcome. They collected questionnaire data, blood samples and hair samples cut closest to the scalp at baseline and at 90 days follow-up. ATV/r in hair were measured using liquid chromatography/mass spectrometry/ mass spectrometry (LC/MS/MS) and the assay range for ATV/r was 0.05–20.0 /0.01–4.0 ng/mg hair, respectively, and a correlation co-efficient of 0.99 for both. Additional details about hair sample preparation and analysis are described in the primary study [15].

Pharmacokinetic modelling

We developed a population PK model to describe ATV/r concentrations in hair. We fixed parameters describing the PK of atazanavir and ritonavir in plasma based on previously published estimates got from studies conducted in almost similar settings [36, 37]. Some of the inter-individual (IIV) and inter-occasion (IOV) variability parameters were estimated while others were fixed in order to improve model fit or stability. We did model development using the first-order conditional estimation method with interaction (FOCE-I) in NONMEM (version 7.3) ADVAN 13 [38]. We schematically presented the structural population PK model applied to both atazanavir and ritonavir concentrations using Fig. 1.

Fig. 1

Schematic representation of the structural population PK model used to predict atazanavir and ritonavir concentrations measured in hair

The model describes the concentration of drug in hair at steady-state trough plasma concentrations in the body 24 h after a dose. Given that there were no plasma concentration data, we used a simplified model which estimate ratios of concentrations between hair and plasma. We describe the rate of change of amount of drug between compartments in Fig. 1 using the differential equations:

$$ \frac{dA_1}{dt}=-{k}_{12}.{A}_1, $$
$$ \frac{dA_2}{dt}={k}_{12}.{A}_1-{k}_{20}.{A}_2, $$
$$ \frac{dA_3}{dt}={F}_{rac}\left({k}_{20}\right).{A}_2, $$
$$ C=\frac{A_3}{V_h}, $$

Equation 1 describe drug absorption into the plasma circulation i.e. central compartment at a rate (k12) proportional to the dose amount (A1); Equation 2 describe input from Equation 1 and total elimination of the drug at rate (\( {k}_{20}=\frac{CL}{V_c} \), where CL and Vc represents clearance and apparent volume of distribution of the bioavailable drug proportional to the amount in the central compartment (A2); Equation 3 describe a ratio (Frac) of hair concentration relative to steady-state plasma trough concentration. Equation 4 predicts the drug concentration C in hair using the ratio of amount of drug in the hair (A3) and apparent volume of distribution (Vh) of the bioavailable drug in the hair.

We tested all covariates in Table 1 during covariate analysis. We selected the optimal covariates relationships through clinical and prior assessment of statistical significance testing using the Stepwise Covariate Model building (SCM) method as implemented in Perl-speaks-NONMEM (PsN). We tested relations on the Frac parameter only.

Table 1 Summary statistics describing data variables of the original study

Model evaluation

We used the change in objective function value (∆OFV) provided from NONMEM model output at 5% level of significance (i.e. ΔOFV > 3.83, Chi-square 1-degree of freedom) in forward selection process and then at 1% level of significance (i.e. ΔOFV > 6.64, Chi-square 1-degree of freedom) in the backward deletion process, to make discriminations between hierarchical models. We performed bootstrap analysis and 90% confidence intervals on the final covariate models by re-sampling 1000 times in PsN as part of model evaluation numerically. We used graphical assessment of the standard goodness-of-fit plots [39]. Both proportional and additive or combined error models were tested and discriminated by means of change in objective function value (∆OFV).


Study participant details

We used 50 participant data in the analysis. The mean (standard deviation) age in years of the study participants was 15.8 (1.8). Fifty-four percent were female. The majority (89%) of these adolescents were attending secondary school, while others were still primary school. Ten percent of the study participants were under the care of fellow siblings below the age of 19 years while others were under the care of parents or relatives. Eighty-two percent of the adolescents were on tenofovir, lamivudine and ATV/r. ATV/r in hair was measured at baseline and 90 days follow-up for every participant and for both drugs, out of the 100 hair samples collected for ATV/r, only 82% / 88%, respectively, were considered for pharmacokinetic modelling, while the remainder 18% / 12% were below the limit of quantification, The median length of hair samples was 1 cm (range 0.5 cm to 1.5 cm). The mean (standard deviation) weight of hair samples was 2.0 g (0.15 g). The mean (standard deviation) drug concentration for atazanavir / ritonavir was 2.5 ng/mg hair (2 ng/mg hair)/0.5 ng/mg hair (0.4 ng/mg hair) respectively. The mean (standard deviation) adherence level of a visual analogue scale of 0 to 100% was 84.1% (18.1%). Further we present details on characteristics of the study population in Table 1.

Population pharmacokinetic modelling

Atazanavir model

We fixed the parameters describing the steady-state population pharmacokinetics of atazanavir in the plasma to 0.44 l per hour for k12, 10 l per hour for CL/F and 63.4 l for Vc [36, 37]. We included body weight as a covariate on CL/F and Vc through allometric scaling, fixing the exponents to 0.75 and 1 for CL/F and Vc respectively [40]. We initially estimated Vh for atazanavir in hair but later on fixed it to 1 because the model estimates were close to 1 and also to stabilise the final model results. We estimated that ATV concentration measured in hair is approximately16% the amount of atazanavir plasma trough concentration after adjusting for covariate effects. Covariate model results show that participants had ATV concentration 30% less at enrolment than that at follow-up event (p-value < 0.0001). A unit increase in self-reported adherence score increased ATV concentration by 2% (p-value = 0.0004). Thinner participants had 54% higher ATV concentration, while overweight participants had 21% lower ATV concentration compared to participants with normal body mass index-for-age (p-value = 0.0165). Participants receiving care from a parent and uncle or aunt had atazanavir in hair 53 and 12% higher respectively, while those receiving care from fellow siblings had atazanavir in hair 54% lower compared to participants receiving care from grandparents (p-value = 0.0406). Based on the change in the OFV value, the most significant covariate was the follow-up occasion (ΔOFV = 47.8, d.f = 1); followed by adherence score (ΔOFV = 18.4, d.f = 1); Body Mass Index-for-age (ΔOFV = 7.5, d.f = 2); and guardian status (ΔOFV = 14.0, d.f = 2), respectively. We present the detailed results in Table 2 and Table 3.

Table 2 Effect of covariate inclusion on the OFV for the atazanavir in hair model
Table 3 Final model parameters describing joint fixed plasma and hair pharmacokinetics of atazanavir

Ritonavir model

We fixed the parameters describing the steady-state population pharmacokinetics of ritonavir in the plasma to 2.31 l per hour for k12, 12.8 l per hour for CL/F and 105 l for Vc [36, 37]. We initially estimated Vh for ritonavir in hair but later on fixed it to 1 because the model estimates were close to 1 and also to stabilise the final model results. We estimated that ritonavir concentration measured in hair is approximately 18% the amount of ritonavir plasma trough concentration after adjusting for covariate effects. Covariate model results show that participants had ritonavir fraction 42% less at enrolment than that at follow-up event (p-value = 0.0003). A unit increase in self-reported adherence score increased ritonavir concentrations by 2% (p-value = 0.0245). Based on the change in the OFV value, the most significant covariate was the follow-up occasion (ΔOFV = 14.8, d.f = 1) and then followed by adherence score (ΔOFV = 5.3, d.f = 1). We present the detailed results in Tables 4 and 5.

Table 4 Effect of covariate inclusion on the OFV for the ritonavir in hair model
Table 5 Final model parameters describing joint fixed plasma and hair pharmacokinetics of ritonavir

Model diagnostics

There was no huge variation between all the estimated final model parameters and those got using 1000 samples bootstrap. All the estimated final model were falling within the 90% confidence intervals. Figure 2 presents the basic goodness-of-fit plots showing the population model predictions versus observations and the residual error plots for both atazanavir and ritonavir final models. The results show low bias, and fairly good precision showing fairly acceptable predictive performance.

Fig. 2

Basic goodness-of-fit plots for the final model for atazanavir 300 mg (a) and ritonavir 100 mg (b). Upper left panel: The observations are plotted versus the population predictions. Upper right panel: The observations are plotted against the individual predictions. Lower left panel: The individually weighted residuals are plotted versus the individual predictions. Lower right panel: The conditional weighted residuals are shown versus time (in hours). The open black circles represents observed data. The bold-dashed line is a locally weighted scatter-plot smoother (LOESS), while the solid line is identity or zero


This is a breakthrough study to perform joint pharmacokinetic modelling of plasma and hair drug concentrations, determine the relationship between exposure of the drug in hair and that to plasma. Several studies have used drug concentration as a tool for measuring antiretroviral drug exposure in situations where non-adherence to treatment maybe a challenge. However the choice of the multivariate statistical models involving hair concentration as the outcome variable in these studies lacked the dose component which plays a critical role when optimising the relationship between drug exposure and treatment outcomes [5, 15, 17,18,19,20]. A non-linear mixed effect PK model has an advantage that includes the dose component. The main purpose of the current model is basically to inform future study design that involve measuring drug concentrations in hair. Later in the discussion, we will present some limitations and recommendations that can improve the power of this method.

Some of the plasma pharmacokinetic parameter we reported while fixing to constant values for the drugs varied from those reported in studies conducted in almost similar settings [36, 37], this could be as a result that some of these values were adjusted by the median body weight observed in our study using allometric scaling. While it is novel to use transit compartment as applied in one of these studies versus the conventional approach (Tlag) to cater for delay in absorption, the disease severity experienced was different in our case due to lower age groups in the previously published study [37]. There is a possibility that certain enzymes and or transporters played a role which physiologically were not possible to incorporate in our model. Hence these two studies were used as a better reference in our study because they were comparable to our study population geographically but had no influence in selection of all parameters describing plasma pharmacokinetics.

The amount of atazanavir concentration determined in hair is approximately 16% of steady-state plasma trough. We estimated an almost similar ratio of 18% for ritonavir as part of model testing. In our conceptual framework, we are interested in finding covariates that affect drug exposure to improve the dosing strategies. We used the SCM in identifying covariates associated with variation in hair drug exposure. The drug that accumulates in hair comes from plasma, therefore one of the major assumption is that the covariates found to have an association with accumulation of drug in hair in our results are a function of an altered plasma PK profile.

ATV/r concentrations increased on follow-up occasion irrespective of study arm which could result from design biases the original study could not eliminate. By being involved in a study, participants are more aware that they are under investigation, hence they adhere more to treatment, increasing hair drug concentrations in both arms. The primary study randomized participants to study arms without blinding, so this could have introduced the biases. High body mass index (BMI)-for-age decreased atazanavir concentrations in hair, while low body mass index-for-age increased atazanavir accumulation in hair. These findings concur with earlier studies in adults [41, 42]. They associated low BMI with high plasma drug concentrations, often resulting in supra-therapeutic drug concentrations, with subsequent drug toxicity, side effects and defaulting treatment. Furthermore, the same literature associates high BMI with low plasma drug concentrations, often resulting in sub-therapeutic drug concentrations and subsequent treatment failure and drug resistance. To the best of our knowledge, our study is the first to show this association in adolescents and using drug concentrations in hair.

We associated receiving care from siblings with lower drug concentrations in hair. Based on current knowledge, this is the 1st study to prove this association using hair samples. Results from a different previous study showed higher self-reported adherence in children and youth who stayed with their parents and grandparents, than those who stayed with siblings [43]. Siblings of HIV-infected children and youth are often immature themselves and still coming to terms with burdens associated with child-headed families, orphan-hood and poverty. The needs of HIV-infected children come in as an extra burden that the siblings may not manage the pressure that comes with the burden, leading to missed doses and hospital visits, and subsequent treatment failure in the HIV-infected adolescents.

A major limitation of the modelling approach applied in the article is we had a small sample size (n = 50) and that we got only two times-points of drug concentration data from each participant. An additional number of participants coupled with having several or segmental measurement of drug concentration from the hair and additionally measuring drug concentration in plasma will improve the power of the modelling framework that used in this paper. Using prior estimates on the plasma PK model could have led to an underestimation of steady state trough concentration because of unavailability of adherence data, however including prior estimates in the form of both fixed and random effects on the plasma PK model could have reduced the bias. Most of the study participants (82%) were on a uniform drug combination, however the unavailability of data on non-HIV/AIDS linked co-therapies limited investigations of the drug-drug interactions which is also very critical to test during PK analysis. Also, apart from the low sample size, the unavailability of data about how the participants cosmetically treated their hair before it was sampled for the study can possible explain why some of the model parameters were reported with notable very high residual standard errors, hence limiting the accuracy of the model estimates.


We have showed some work which can complement the efforts being taken by other scientists to establish the use of measuring drug concentration in hair at HIV/AIDS points of care. Most important determinants of increased concentrations in hair were monitoring at follow up event, BMI-for-age and caregiver. Measuring ARV concentrations in hair promises to be more accurate and feasibly accomplished. It is crucial to perform follow-up work which involves establishing the relationship between hair drug concentrations and a measure of treatment response such as viral loads. Comparing the predictive accuracy for exposure-response models when exposure of interest is plasma or hair drug concentrations is necessary to perform.

Availability of data and materials

The primary study did not publish the data to the public, however the data can be available upon request and approval from the principal investigators of the primary study.







Nonlinear mixed effect modelling


Human immunodeficiency virus


Acquired immunodeficiency syndrome


Body mass index




Inter-individual variability


Inter-occasion variability


First-order conditional estimation method with interaction




Stepwise covariate model building


Change in objective function value


Nanogram per milligram

k 12 :

Absorption rate constant


Apparent drug clearance

V c :

Apparent volume of distribution in the central compartments

V h :

Apparent volume of distribution in the hair compartments

F rac :

Amount of drug cleared into the hair as a proportion of the amount of drug in plasma at steady-state troughs

FACTOR_ Frac :

effect of covariate on Frac

ɛADD :

Additive error term


Proportional error term


Residual error


Standard error


Probability value


Confidence interval


  1. 1.

    Global and Regional Trends. UNICEF DATA. Available from: // [cited 2017 Dec 28].

  2. 2.

    ZWE_2018_countryreport.pdf. Available from: [cited 2019 Jan 21].

  3. 3.

    Paterson DL, Swindells S, Mohr J, et al. Adherence to protease inhibitor therapy and outcomes in patients with HIV infection. Ann Intern Med. 2002;136(3):21–30.

    Google Scholar 

  4. 4.

    Cohen MS, Chen YQ, McCauley M, et al. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365(6):493–505.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Chawana TD, Katzenstein D, Nathoo K, Ngara B, CFB N. Evaluating an enhanced adherence intervention among HIV positive adolescents failing atazanavir/ritonavir-based second line antiretroviral treatment at a public health clinic. J AIDS HIV Res. 2017;9(1):17–30.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Makadzange AT, Higgins-Biddle M, Chimukangara B, et al. Clinical, Virologic, Immunologic Outcomes and Emerging HIV Drug Resistance Patterns in Children and Adolescents in Public ART Care in Zimbabwe. Plos One. 2015;10(12) Available from: [cited 2017 June 13].

  7. 7.

    Zeleke A. Prevalence of antiretroviral treatment failure and associated factors in HIV infected children on antiretroviral therapy at Gondar University hospital, retrospective cohort study. Int J Med Med Sci. 2016;8(11):125–32.

    Article  Google Scholar 

  8. 8.

    Yassin S, Gebretekle GB. Magnitude and predictors of antiretroviral treatment failure among HIV-infected children in Fiche and Kuyu hospitals, Oromia region, Ethiopia: a retrospective cohort study. Pharmacol Res Perspect. 2017;5(1) Available from: [cited 2017 June 13].

  9. 9.

    Adejumo OA, Malee KM, Ryscavage P, Hunter SJ, Taiwo BO. Contemporary issues on the epidemiology and antiretroviral adherence of HIV-infected adolescents in sub-Saharan Africa: a narrative review. J Int AIDS Soc. 2015;18(1) Available from: [cited 2017 June 13].

  10. 10.

    Nglazi MD, Kranzer K, Holele P, et al. Treatment outcomes in HIV-infected adolescents attending a community-based antiretroviral therapy clinic in South Africa. BMC Infect Dis. 2012;12:21.

    PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Okawa S, Chirwa M, Ishikawa N, et al. Longitudinal adherence to antiretroviral drugs for preventing mother-to-child transmission of HIV in Zambia. BMC Pregnancy Childbirth. 2015;15:258.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Bonner K, Mezochow A, Roberts T, Ford N, Cohn J. Viral load monitoring as a tool to reinforce adherence: a systematic review. J Acquir Immune Defic Syndr 1999. 2013;64(1):74–8.

    Article  Google Scholar 

  13. 13.

    Müller AD, Jaspan HB, Myer L, et al. Standard measures are inadequate to monitor pediatric adherence in a resource-limited setting. AIDS Behav. 2011;15(2):422–31.

    PubMed  Article  Google Scholar 

  14. 14.

    Burack G, Gaur S, Marone R, Petrova A. Adherence to antiretroviral therapy in pediatric patients with human immunodeficiency virus (HIV-1). J Pediatr Nurs. 2010;25(6):500–4.

    PubMed  Article  Google Scholar 

  15. 15.

    Chawana TD, Gandhi M, Nathoo K, et al. Defining a cut-off for atazanavir in hair samples associated with virological failure among adolescents failing second-line antiretroviral treatment. J Acquir Immune Defic Syndr. 2017;76(1):55–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Beumer JH, Bosman IJ, Maes RA. Hair as a biological specimen for therapeutic drug monitoring. Int J Clin Pract. 2001;55(6):353–7.

    CAS  PubMed  Google Scholar 

  17. 17.

    Gandhi M, Greenblatt RM. Hair it is: the long and short of monitoring antiretroviral treatment. Ann Intern Med. 2002;137(8):696–7.

    PubMed  Article  Google Scholar 

  18. 18.

    Gandhi M, Ameli N, Bacchetti P, et al. Protease inhibitor levels in hair samples strongly predict Virologic responses to HIV treatment. AIDS Lond Engl. 2009;23(4):471–8.

    CAS  Article  Google Scholar 

  19. 19.

    Gandhi M, Ameli N, Bacchetti P, et al. Atazanavir concentration in hair is the strongest predictor of outcomes on antiretroviral therapy. Clin Infect Dis Off Publ Infect Dis Soc Am. 2011;52(10):1267–75.

    CAS  Article  Google Scholar 

  20. 20.

    Hickey MD, Salmen CR, Tessler RA, et al. Antiretroviral concentrations in small hair samples as a feasible marker of adherence in rural Kenya. J Acquir Immune Defic Syndr 1999. 2014;66(3):311–5.

    CAS  Article  Google Scholar 

  21. 21.

    Duval X, Peytavin G, Breton G, et al. Hair versus plasma concentrations as indicator of indinavir exposure in HIV-1-infected patients treated with indinavir/ritonavir combination. AIDS Lond Engl. 2007;21(1):106–8.

    CAS  Article  Google Scholar 

  22. 22.

    Kidwell DA, Blank DL. Comments on the paper by W.a. Baumgartner and V.a. Hill: sample preparation techniques. Forensic Sci Int. 1993;63(1):137–43.

    CAS  Article  Google Scholar 

  23. 23.

    DuPont RL, Baumgartner WA. Drug testing by urine and hair analysis: complementary features and scientific issues. Forensic Sci Int. 1995;70(1):63–76.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Baumgartner WA, Hill VA. Comments on the paper by David L. Blank and David a. Kidwell: external contamination of hair by cocaine: an issue in forensic interpretation. Forensic Sci Int. 1993;63(1):157–60.

    CAS  Article  Google Scholar 

  25. 25.

    Al-Delaimy WK. Hair as a biomarker for exposure to tobacco smoke. Tob Control. 2002;11(3):176–82.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Khajuria H, Nayak BP. Detection and accumulation of morphine in hair using GC–MS. Egypt J Forensic Sci. 2016;6(4):337–41.

    Article  Google Scholar 

  27. 27.

    Kronstrand R, Grundin R, Jonsson J. Incidence of opiates, amphetamines, and cocaine in hair and blood in fatal cases of heroin overdose. Forensic Sci Int. 1998;92(1):29–38.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Minoli M, Angeli I, Ravelli A, Gigli F, Lodi F. Detection and quantification of 11-nor-Δ9-tetrahydrocannabinol-9-carboxylic acid in hair by GC/MS/MS in negative chemical ionization mode (NCI) with a simple and rapid liquid/liquid extraction. Forensic Sci Int. 2012;218(1–3):49–52.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Lee D, Milman G, Barnes AJ, Goodwin RS, Hirvonen J, Huestis MA. Oral fluid cannabinoids in chronic, daily Cannabis smokers during sustained, monitored abstinence. Clin Chem. 2011;57(8):1127–36.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Schaffer M, Hill V, Cairns T. Hair analysis for cocaine: the requirement for effective wash procedures and effects of drug concentration and hair porosity in contamination and decontamination. J Anal Toxicol. 2005;29(5):319–26.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Drugs-in-Hair-FAQ.pdf. Available from: [cited 2017 Jul 4].

  32. 32.

    Dhoro M, Zvada S, Ngara B, et al. CYP2B6*6, CYP2B6*18, Body weight and sex are predictors of efavirenz pharmacokinetics and treatment response: population pharmacokinetic modeling in an HIV/AIDS and TB cohort in Zimbabwe. BMC Pharmacol Toxicol. 2015;16:4.

    PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Nyakutira C, Röshammar D, Chigutsa E, et al. High prevalence of the CYP2B6 516G→T(*6) variant and effect on the population pharmacokinetics of efavirenz in HIV/AIDS outpatients in Zimbabwe. Eur J Clin Pharmacol. 2008;64(4):357–65.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Mukonzo JK, Röshammar D, Waako P, et al. A novel polymorphism in ABCB1 gene, CYP2B6*6 and sex predict single-dose efavirenz population pharmacokinetics in Ugandans. Br J Clin Pharmacol. 2009;68(5):690–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Nemaura T, Nhachi C, Masimirembwa C. Impact of gender, weight and CYP2B6 genotype on efavirenz exposure in patients on HIV/AIDS and TB treatment: implications for individualising therapy. Afr J Pharm Pharmacol. 2012;6(29):2188–93.

    CAS  Google Scholar 

  36. 36.

    Foissac F, Blanche S, Dollfus C, et al. Population pharmacokinetics of atazanavir/ritonavir in HIV-1-infected children and adolescents. Br J Clin Pharmacol. 2011;72(6):940–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Zhang C, Denti P, Decloedt EH, Ren Y, Karlsson MO, McIlleron H. Model-based evaluation of the pharmacokinetic differences between adults and children for lopinavir and ritonavir in combination with rifampicin. Br J Clin Pharmacol. 2013;76(5):741–51.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    NONMEM® - ICON plc. Available from: [cited 2017 Jun 20].

  39. 39.

    Keizer RJ, Karlsson MO, Hooker A. Modeling and simulation workbench for NONMEM: tutorial on Pirana, PsN, and Xpose. CPT Pharmacomet Syst Pharmacol. 2013;2(6):e50.

    CAS  Article  Google Scholar 

  40. 40.

    Anderson BJ, Holford NHG. Mechanism-based concepts of size and maturity in pharmacokinetics. Annu Rev Pharmacol Toxicol. 2008;48:303–32.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Lowenthal ED, Bakeera-Kitaka S, Marukutira T, Chapman J, Goldrath K, Ferrand RA. Perinatally acquired HIV infection in adolescents from sub-Saharan Africa: a review of emerging challenges. Lancet Infect Dis. 2014;14(7):627–39.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Agwu AL, Fairlie L. Antiretroviral treatment, management challenges and outcomes in perinatally HIV-infected adolescents. J Int AIDS Soc. 2013;16:18579.

    PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Gichane MW, Sullivan KA, Shayo AM, et al. Caregiver role in HIV medication adherence among HIV-infected orphans in Tanzania. AIDS Care. 2018;30(6):701–5.

    PubMed  Article  Google Scholar 

Download references


Authors acknowledge the entire research team of the original study which include researchers, research nurses, the trained community workers, participants, and a team from the University of California San Francisco (United States of America) which worked on the hair assays for their effort in producing the data which we used for secondary analysis. We also acknowledge the University of Cape Town (South Africa) for providing a computing facility with the NONMEM software.


The corresponding author obtained doctorate scholarship from Letten Foundation, Norway. This research is also commissioned by the National Institute for Health Research, using Official Development Assistance (ODA) funding 16/136/33. The views expressed in this publication are those of the authors only and the funders did not play any role in the design of the study; collection, analysis, interpretation of data, and in writing the manuscript.

Author information




B.N, S.P.Z, and S. R model development and validation and statistical interpretation of model results. T.D.C and C. N provided the data. T.D.C, B.S.P and C. N clinical and pharmacological interpretation and review of model results. B. N, S.P.Z and T.D.C wrote the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Bernard Ngara.

Ethics declarations

Ethics approval and consent to participate

The study got approval to use secondary data from the institutional and national ethical review committees [Joint Research Ethics Committee (JREC/101/18) and Medical Research Council of Zimbabwe (MRCZ/A/2301) respectively].

Consent for publication


Competing interests

The authors declare that there is no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ngara, B., Zvada, S., Chawana, T.D. et al. A population pharmacokinetic model is beneficial in quantifying hair concentrations of ritonavir-boosted atazanavir: a study of HIV-infected Zimbabwean adolescents. BMC Pharmacol Toxicol 21, 58 (2020).

Download citation


  • Pharmacokinetic modelling
  • Adolescents
  • Adherence
  • Hair