Joint Population Pharmacokinetic Model Research Articles

A joint population pharmacokinetic model to assess the high variability of whole-blood and intracellular tacrolimus in early adult renal transplant recipients

Tacrolimus (TAC) has high pharmacokinetic (PK) variability during the early transplantation period. The relationships between whole-blood and intracellular TAC concentrations and clinical outcomes remain controversial. This study identifies the factors affecting the PK variability of TAC and characterizes the relationships between whole-blood and intracellular TAC concentrations. Data regarding whole-blood TAC concentrations of 1,787 samples from 215 renal transplant recipients (<90 days postoperative) across two centers and intracellular TAC concentrations (648 samples) digitized from previous studies were analyzed using nonlinear mixed-effects modeling. The effects of potential covariates were screened, and the distribution of whole-blood to intracellular TAC concentration ratios (RWB:IC) was estimated. The final model was evaluated using bootstrap, goodness of fit, and prediction-corrected visual predictive checks. The optimal dosing regimens and target ranges for each type of immune cell subsets were determined using Monte Carlo simulations. A two-compartment model adequately described the data, and the estimated mean TAC CL/F was 23.6 L·h−1 (relative standard error: 11.5 %). The hematocrit level, CYP3A5*3 carrier status, co-administration with Wuzhi capsules, and tapering prednisolone dose may contribute to the high variability of TAC PK variability during the early post-transplant period. The estimated RWB:IC of all TAC concentrations in peripheral blood mononuclear cells (PBMCs) was 4940, and inter-center variability of PBMCs was observed. The simulated TAC target range in PBMCs was 20.2–85.9 pg·million cells−1. Inter-center variability in intracellular concentrations should be taken into account in further analyses. TAC dosage adjustments can be guided based on PK/PD variability and simulated intracellular concentrations.

Joint population pharmacokinetic modeling of venlafaxine and O-desmethyl venlafaxine in healthy volunteers and patients to evaluate the impact of morbidity and concomitant medication.

Introduction: Venlafaxine (VEN) is a widely used dual selective serotonin/noradrenaline reuptake inhibitor indicated for depression and anxiety. It undergoes first-pass metabolism to its active metabolite, O-desmethyl venlafaxine (ODV). The aim of the present study was to develop a joint population pharmacokinetic (PPK) model to characterize their pharmacokinetic characters simultaneously. Methods: Plasma concentrations with demographic and clinical data were derived from a bioequivalence study in 24 healthy subjects and a naturalistic TDM setting containing 127 psychiatric patients. A parent-metabolite PPK modeling was performed with NONMEM software using a non-linear mixed effect modeling approach. Goodness of fit plots and normalized prediction distribution error method were used for model validation. Results and conclusion: Concentrations of VEN and ODV were well described with a one-compartment model incorporating first-pass metabolism. The first-pass metabolism was modeled as a first-order conversion. The morbid state and concomitant amisulpride were identified as two significant covariates affecting the clearance of VEN and ODV, which may account for some of the variations in exposure. This model may contribute to the precision medication in clinical practice and may inspire other drugs with pre-system metabolism.

Development of a Population Pharmacokinetic Model for the Diroximel Fumarate Metabolites Monomethyl Fumarate and 2-Hydroxyethyl Succinimide Following Oral Administration of Diroximel Fumarate in Healthy Participants and Patients with Multiple Sclerosis.

IntroductionDiroximel fumarate (DRF) is a next-generation oral fumarate that is indicated in the USA for relapsing forms of multiple sclerosis (MS). A joint population pharmacokinetic model was developed for the major active metabolite (monomethyl fumarate, MMF) and the major inactive metabolite (2-hydroxyethyl succinimide, HES) of DRF.MethodsMMF and HES data were included from 341 healthy volunteers and 48 patients with MS across 11 phase I and III studies in which DRF was administered as single or multiple doses. Population modeling was performed with NONMEM version 7.3 with the first-order conditional estimation method.ResultsEstimated MMF clearance (CLMMF), volume of distribution, and absorption rate constant (Ka) were 13.5 L/h, 30.4 L, and 5.04 h−1, respectively. CLMMF and HES clearance (CLHES) increased with increasing body weight. CLHES decreased with decreasing renal function. CLMMF and CLHES were 28% and 12% lower in patients with MS than in healthy volunteers, respectively. Ka was reduced in the presence of low-, medium-, and high-fat meals by 37%, 51%, and 67%, respectively, for MMF; and by 34%, 49%, and 62%, respectively, for HES.ConclusionsAge, sex, race, and baseline liver function parameters such as total bilirubin, albumin, and aspartate aminotransferase were not considered to be significant predictors of MMF or HES disposition.Supplementary InformationThe online version contains supplementary material available at 10.1007/s40120-021-00316-6.

Application of Population Pharmacokinetic Analysis to Characterize CYP2C19 Mediated Metabolic Mechanism of Voriconazole and Support Dose Optimization.

Purpose: The aims of this study were to establish a joint population pharmacokinetic model for voriconazole and its N-oxide metabolite in immunocompromised patients, to determine the extent to which the CYP2C19 genetic polymorphisms influenced the pharmacokinetic parameters, and to evaluate and optimize the dosing regimens using a simulating approach. Methods: A population pharmacokinetic analysis was conducted using the Phoenix NLME software based on 427 plasma concentrations from 78 patients receiving multiple oral doses of voriconazole (200 mg twice daily). The final model was assessed by goodness of fit plots, non-parametric bootstrap method, and visual predictive check. Monte Carlo simulations were carried out to evaluate and optimize the dosing regimens. Results: A one-compartment model with first-order absorption and mixed linear and concentration-dependent-nonlinear elimination fitted well to concentration-time profile of voriconazole, while one-compartment model with first-order elimination well described the disposition of voriconazole N-oxide. Covariate analysis indicated that voriconazole pharmacokinetics was substantially influenced by the CYP2C19 genetic variations. Simulations showed that the recommended maintenance dose regimen would lead to subtherapeutic levels in patients with different CYP2C19 genotypes, and elevated daily doses of voriconazole might be required to attain the therapeutic range. Conclusions: The joint population pharmacokinetic model successfully characterized the pharmacokinetics of voriconazole and its N-oxide metabolite in immunocompromised patients. The proposed maintenance dose regimens could provide a rationale for dosage individualization to improve clinical outcomes and minimize drug-related toxicities.

Population Pharmacokinetics of Phosphocreatine and Its Metabolite Creatine in Children With Myocarditis.

Objective: This study aimed to develop a parent-metabolite joint population pharmacokinetic model to characterize the pharmacokinetic (PK) profile for phosphocreatine (PCr) and its metabolite creatine (Cr) in children with myocarditis and to use this model to study the PK profile of different dosing schemes. Methods: One hundred pediatric patients with myocarditis were enrolled. Blood samples were collected at baseline and approximately 30, 40 or 50, 75, and 180 min after a single dose of phosphocreatine sodium. Plasma PCr and Cr concentrations were determined using an HPLC-MS/MS method. A nonlinear mixed effect model approach was used to build the population pharmacokinetic model. After validation, the model was used for simulations to evaluate the PK profile of different dosing schemes. Results: A total of 997 plasma concentrations (498 for PCr and 499 for Cr) were included in the analysis. A four-compartment chain model (central and peripheral compartments for both PCr and Cr) with the first-order elimination adequately characterized the in vivo process of PCr and Cr. Allometric scaling based on bodyweight was applied to the PK parameters. The covariate analysis identified that the glomerular filtration rate (GFR) was strongly associated with Cr clearance. Bootstrapping and visual predictive checks suggested that a robust and reliable pharmacokinetic model was developed. The simulation results showed that PCr had no accumulation in vivo. With the infusion of PCr, the concentration of Cr increased rapidly. Conclusion: A joint population pharmacokinetic model for PCr and Cr in pediatric patients with myocarditis was successfully developed for the first time.

Quantification of In Vivo Metabolic Activity of CYP2D6 Genotypes and Alleles Through Population Pharmacokinetic Analysis of Vortioxetine.

Assignment of CYP2D6 phenotype from genotype data can be challenging and despite efforts to standardize translation, there is currently no universally accepted method. To facilitate standardization, there remains a need to precisely quantify the in vivo function of different CYP2D6 genotypes. Vortioxetine is metabolized to its major metabolite, Lu AA34443, primarily via CYP2D6. The aim of this study was to quantify the in vivo CYP2D6 activity of different CYP2D6 alleles and genotypes through population pharmacokinetic (PopPK) modeling of vortioxetine and Lu AA34443. Plasma concentration data of vortioxetine and Lu AA34443 from 1,140 subjects originating from 29 clinical pharmacology studies were pooled for the analysis. A joint PopPK model described the pharmacokinetics of vortioxetine and Lu AA34443 simultaneously and provided estimates of the CYP2D6‐mediated metabolism for each subject. Subjects normally classified as CYP2D6 intermediate metabolizers (IMs) showed different levels of CYP2D6 activity with carriers of one fully functional allele and one null function allele having 77% higher CYP2D6 activity compared with carriers of two decreased function alleles (P < 0.0001). The decreased function alleles were associated with different levels of reduction of CYP2D6 activity. Fixing the activity of fully functional alleles to 1.0, the relative activities of CYP2D6*9, CYP2D6*10, CYP2D6*17, and CYP2D6*41 were 0.22, 0.37, 0.17, and 0.21, respectively. The activity of CYP2D6*10 was shown to be significantly greater than that of CYP2D6*17 (P = 0.01) and CYP2D6*41 (P = 0.02). These results warrant further discussion of current CYP2D6 genotype‐phenotype classification systems particularly regarding decreased function alleles and the IM phenotype.

Amount of Cycloserine Emanating from Terizidone Metabolism and Relationship with Hepatic Function in Patients with Drug-Resistant Tuberculosis

Background and objectivesThe dosing of cycloserine and terizidone is the same, as both drugs are considered equivalent or used interchangeably. Nevertheless, it is not certain from the literature that these drugs are interchangeable. Therefore, the amount of cycloserine resulting from the metabolism of terizidone and the relationship with hepatic function were determined.MethodsThis prospective clinical study involved 39 patients with drug-resistant tuberculosis admitted for an intensive phase of treatment. Cycloserine pharmacokinetic parameters for individual patients, like area under the curve (AUC), clearance (CLm/F), peak concentration (Cmax) and trough concentration (Cmin), were calculated from a previously validated joint population pharmacokinetic model of terizidone and cycloserine. Correlation and regression analyses were performed for pharmacokinetic parameters and unconjugated bilirubin (UB), conjugated bilirubin (CB), albumin, the ratio of aspartate transaminase to alanine aminotransferase (AST/ALT), or binding affinity of UB to albumin (Kaf), using R statistical software version 3.5.3.ResultsThirty-eight patients took a daily dose of 750 mg terizidone, while one took 500 mg. The amount of cycloserine [median (range)] that emanated from terizidone metabolism was 51.6 (0.64–374) mg. Cmax (R2 = 22%, p = 0.003) and Cmin (R2 = 10.6%, p = 0.044) were significantly associated with increased CB concentration. Cmax was significantly associated with increased Kaf (R2 = 10.1%, p = 0.048), while high CLm/F was significantly associated with decreased AST/ALT (R2 = 21%, p = 0.003).ConclusionsCycloserine is not interchangeable with terizidone, as amounts are lower than expected. Cycloserine may be a predisposing factor to the development of hyperbilirubinaemia, as CLm/F is affected by hepatic function.

A population pharmacokinetic model for simvastatin and its metabolites in children and adolescents

PurposePoor adherence to dietary/behaviour modifications as interventions for hypercholesterolemia in paediatric patients often necessitates the initiation of statin therapy. The aim of this study was to develop a joint population pharmacokinetic model for simvastatin and four metabolites in children and adolescents to investigate sources of variability in simvastatin acid exposure in this patient population, in addition to SLCO1B1 genotype status.MethodsPlasma concentrations of simvastatin and its four metabolites, demographic and polymorphism data for OATP1B1 and CYP3A5 were analysed utilising a population pharmacokinetic modelling approach from an existing single oral dose (10 mg < 17 years and 20 mg ≥ 18 years) pharmacokinetic dataset of 32 children and adolescents.ResultsThe population PK model included a one compartment disposition model for simvastatin with irregular oral absorption described by two parallel absorption processes each consisting of sequential zero and first-order processes. The data for each metabolite were described by a one-compartment disposition model with the formation and elimination apparent parameters estimated. The model confirmed the statistically significant effect of c.521T>C (rs4149056) on the pharmacokinetics of the active metabolite simvastatin acid in children/adolescents, consistent with adult data. In addition, age was identified as a covariate affecting elimination clearances of 6-hydroxymethyl simvastatin acid and 3, 5 dihydrodiol simvastatin metabolites.ConclusionThe model developed describes the pharmacokinetics of simvastatin and its metabolites in children/adolescents capturing the effects of both c.521T>C and age on variability in exposure in this patient population. This joint simvastatin metabolite model is envisaged to facilitate optimisation of simvastatin dosing in children/adolescents.

Mining Small Routine Clinical Data: A Population Pharmacokinetic Model and Optimal Sampling Times of Capecitabine and its Metabolites.

The present study was performed to demonstrate that small amounts of routine clinical data allow to generate valuable knowledge. Concretely, the aims of this research were to build a joint population pharmacokinetic model for capecitabine and three of its metabolites (5-DFUR, 5-FU and 5-FUH2) and to determine optimal sampling times for therapeutic drug monitoring. We used data of 7 treatment cycles of capecitabine in patients with metastatic colorectal cancer. The population pharmacokinetic model was built as a multicompartmental model using NONMEM and was internally validated by visual predictive check. Optimal sampling times were estimated using PFIM 4.0 following D-optimality criterion. The final model was a multicompartmental model which represented the sequential transformations from capecitabine to its metabolites 5-DFUR, 5-FU and 5-FUH2 and was correctly validated. The optimal sampling times were 0.546, 0.892, 1.562, 4.736 and 8 hours after the administration of the drug. For its correct implementation in clinical practice, the values were rounded to 0.5, 1, 1.5, 5 and 8 hours after the administration of the drug. Capecitabine, 5-DFUR, 5-FU and 5-FUH2 can be correctly described by the joint multicompartmental model presented in this work. The aforementioned times are optimal to maximize the information of samples. Useful knowledge can be obtained for clinical practice from small databases.

Population pharmacokinetic model of irinotecan and its metabolites in patients with metastatic colorectal cancer.

Irinotecan (CPT-11) is a drug used against a wide range of tumor types. The individualized dosing of CPT-11 is essential to ensure optimal pharmacotherapy in cancer patients, given the wide interindividual pharmacokinetic variability of this drug and its active metabolite SN-38. Moreover, the reabsorption from SN-38-G to SN-38, by enterohepatic recirculation, is critical due to its influence in the treatment tolerance. The aim of this research was to build a joint population pharmacokinetic model for CPT-11 and its metabolites (SN-38, and its glucuronide, SN-38-G) that enabled an individualized posology adjustment. We used data of 53 treatment cycles of FOLFIRINOX scheme corresponding to 20 patients with metastatic colorectal cancer. In order to build the population pharmacokinetic model, we implemented parametric and non-parametric methods using the Pmetrics library package for R. We also built multivariate regression models to predict the area under the curve and the maximum concentration using basal covariates. The final model was a multicompartmental model which represented the transformations from CPT-11 to its active metabolite SN-38 and from SN-38 to inactive SN-38-G. Besides, the model also represented the extensive elimination of SN-38-G and the reconversion of the remaining SN-38-G to SN-38 by enterohepatic recirculation. We carried out internal validation with 1000 simulations. The regression models predicted the PK parameters with R squared adjusted up to 0.9499. CPT-11, SN-38, and SN-38-G can be correctly described by the multicompartmental model presented in this work. As far as we know, it is the first time that a joint model for CPT-11, SN-38, and SN-38-G that includes the process of reconversion from SN-38-G to SN-38 is characterized.

Development of a joint population pharmacokinetic model of ezetimibe and its conjugated metabolite

Ezetimibe (EZE) is an extensively used antihyperlipidemic drug with an important cholesterol lowering activity. It undergoes extensive first-pass metabolism to form its active glucuronide metabolite (EZEG). Both drugs exhibit complex pharmacokinetic profiles attributed mainly to repetitive enterohepatic kinetics. The aim of the present study was the investigation of EZE and EZEG pharmacokinetics (PK), through the development of a joint population pharmacokinetic model able to characterize their kinetic processes and enterohepatic recirculation simultaneously. Concentration-time data derived from a bioequivalence study in 28 healthy subjects were used for the analysis. Population PK modeling was performed on the obtained data using nonlinear mixed effect modeling approach, where different methodologies were applied for the description of the complex metabolism and recirculation processes of the two compounds. EZE and EZEG concentrations were best described by a population PK model incorporating first-pass metabolism and an enterohepatic recirculation loop, accounting for the recycling process of the two moieties. This is the first joint population pharmacokinetic model describing the kinetics of both EZE and EZEG.

Model-based approaches for ivabradine development in paediatric population, part II: PK and PK/PD assessment.

The objectives of this work were first to describe the pharmacokinetic (PK) of ivabradine and its active metabolite in a paediatric patient population after repeated oral administration of ivabradine using a population PK approach, and secondly to assess whether the blood/plasma ratio and the pharmacokinetic/pharmacodynamic (PK/PD) relationship are preserved in the paediatric population in comparison to adult. PK data for 70 patients were obtained after blood sampling using dried blood spot and one plasma sample in order to assess the relationship between blood and plasma concentration. In order to describe ivabradine and its metabolite blood concentrations in children, a joint population PK model was developed taking into account weight & age effects on PK parameters. Plasma PK exposure parameters were calculated in children using plasma PK profiles. In order to assess the PK/PD relationship in children, an adult PK/PD model was used. The relationship between blood and plasma concentrations was described using linear mixed effect models. Two and one-compartment models best described parent and metabolite dispositions. Weight effects were fixed to the allometric values of ¾ on clearance (CL) and 1 on volume. A maturation function was added on metabolite formation clearance (CL PM ) reflecting enzyme maturation. Plasma exposure comparison indicated that higher dose/kg were necessary to achieve a similar exposure between younger and older children. No differences between age classes were observed in terms of range of exposure at the maintenance dose. The PK/PD relationship in adult patients is conserved in children.

Identification of the Effect of Multiple Polymorphisms on the Pharmacokinetics of Simvastatin and Simvastatin Acid Using a Population-Modeling Approach

The aim of this work was to develop a joint population pharmacokinetic model for simvastatin (SV) and its active metabolite, simvastatin acid (SVA), that incorporates the effects of multiple genetic polymorphisms and clinical/demographic characteristics. SV/SVA plasma concentrations, demographic/clinical data, and genotypes for 18 genetic variants were collected from 74 individuals (three clinical trials) and analyzed using a nonlinear mixed-effects modeling approach. The structural model that best described the data included a two- and a one-compartment disposition model for SV and SVA, respectively. Age, weight, Japanese ethnicity, and seven genetic polymorphisms-rs4149056 (SLCO1B1), rs776746 (CYP3A5), rs12422149 (SLCO2B1), rs2231142 (ABCG2), rs4148162 (ABCG2), rs4253728 (PPARA), and rs35599367 (CYP3A4)-were identified as significantly affecting model parameters. The developed model was used to assess combinations of these covariates, highlighting specific risk factors associated with altered SV/SVA pharmacokinetics, and consequently myopathy cases that cannot be solely attributed to the rs4149056 CC genotype.

Abstract 2213: Influence of genetic polymorphism in UGT1A1, UGT1A7 and UGT1A9 on irinotecan, SN-38 and SN-38G pharmacokinetics.

Abstract Background and objective: SN-38 is the active metabolite of irinotecan (CPT-11) and mainly responsible of hematological and intestinal toxicity after CPT-11 treatments. SN-38 is inactivated into SN-38G by uridine diphosphate glucuronosyltransferase (UGT) enzymes. Some polymorphism of these enzymes leads lower clearance (CL) of SN-38 and therefore bigger exposition and bigger risk of hematological and/or intestinal toxicity (1). Therefore, the goal of this study has been to evaluate the influence of genetic polymorphism in UGT1A1, UGT1A7 and UGT1A9 on the population pharmacokinetics of CPT-11 and its metabolites, SN-38 and SN-38G on Caucassian cancer population. Methods: Plasma concentrations of CPT-11, SN-38 and SN-38G from 72 patients were pooled to develop a joint population pharmacokinetic model using NONMEM 7. The effect of age, sex, body surface area, total bilirubin, co-medication, tumor type, and UGT1A1, UGT1A7 and UGT1A9 genotypes on the model parameters was evaluated. Results: The typical values (between-subject variability; %) of the irinotecan, SN-38 and SN-38G clearances were 42,9 L/h (56,4%), 1340 L/h (76,8%) and 188 L/h (70,1%), respectively. Clearance in patients with some allele (heterozygous or homozygous) associated with low enzymatic activity of UGT1A1*28 was a 35,7% lower (CI 95%: 6,7 - 64,7; p=0,003) than in wild-type patients. Additionally, linkage dissequilibrium indicates a strong interaction between functional UGT1A1, UGT1A7 and UGT1A9 polymorphisms related to the alteration of the UGT enzyme activity. Conclusion: The impact of UGT1A1*28 genotype on the systemic exposure of SN-38 justifies its routine determination in patients receiving CPT-11 and subsequent dose adjustment. Citation Format: Belen Valenzuela, Mario González-Sales, Vanessa Escudero, Elena Navarro, Carlos Perez-Ruixo, Joseba Rebollo, Ramon González-Manzano, Juan José Pérez-Ruixo. Influence of genetic polymorphism in UGT1A1, UGT1A7 and UGT1A9 on irinotecan, SN-38 and SN-38G pharmacokinetics. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2213. doi:10.1158/1538-7445.AM2013-2213

Development of a Complex Parent-Metabolite Joint Population Pharmacokinetic Model

This study aimed to develop a joint population pharmacokinetic model for an antipsychotic agent in development (S33138) and its active metabolite (S35424) produced by reversible metabolism. Because such a model leads to identifiability problems and numerical difficulties, the model building was performed using the FOCE-I and the Stochastic Approximation Expectation Maximization (SAEM) estimation algorithms in NONMEM and MONOLIX, respectively. Four different structural models were compared based on Bayesian information criteria. Models were first written as ordinary differential equations systems and then in closed form (CF) to facilitate further analyses. The impact of polymorphisms on genes coding for the CYP2C19 and CYP2D6 enzymes, respectively involved in the parent drug and the metabolite elimination were investigated using permutation Wald test. The parent drug and metabolite plasma concentrations of 101 patients were analyzed on two occasions after 4 and 8 weeks of treatment at 1, 3, 6, and 24 h following daily oral administration. All configurations led to a two compartment model with back-transformation of the metabolite into the parent drug and a first-pass effect. The elimination clearance of the metabolite through other processes than back-transformation was decreased by 35% [9-53%] in CYP2D6 poor metabolizer. Permutation tests were performed to ensure the robustness of the analysis, using SAEM and CF. In conclusion, we developed a complex joint pharmacokinetic model adequately predicting the impact of CYP2D6 polymorphisms on the parent drug and its metabolite concentrations through the back-transformation mechanism.

Joint Population Pharmacokinetic Analysis of Zidovudine, Lamivudine, and Their Active Intracellular Metabolites in HIV Patients

The population pharmacokinetic parameters of zidovudine (AZT), lamivudine (3TC), and their active intracellular metabolites in 75 naïve HIV-infected patients receiving an oral combination of AZT and 3TC twice daily as part of their multitherapy treatment in the COPHAR2-ANRS 111 trial are described. Four blood samples per patient were taken after 2 weeks of treatment to measure drug concentrations at steady state. Plasma AZT and 3TC concentrations were measured in 73 patients, and among those, 62 patients had measurable intracellular AZT-TP and 3TC-TP concentrations. For each drug, a joint population pharmacokinetic model was developed and we investigated the influence of different covariates. We then studied correlations between the mean plasma and intracellular concentrations of each drug. A one-compartment model with first-order absorption and elimination best described the plasma AZT concentration, with an additional compartment for intracellular AZT-TP. A similar model but with zero-order absorption was found to adequately described concentrations of 3TC and its metabolite 3TC-TP. The half-lives of AZT and 3TC were 0.81 h (94.8%) and 2.97 h (39.2%), respectively, whereas the intracellular half-lives of AZT-TP and 3TC-TP were 10.73 h (69%) and 21.16 h (44%), respectively. We found particularly a gender effect on the apparent bioavailability of AZT, as well as on the mean plasma and intracellular concentrations of AZT, which were significantly higher in females than in males. Relationships between mean plasma drug and intracellular metabolite concentrations were also highlighted both for AZT and for 3TC. Simulation with the model of plasma and intracellular concentrations for once- versus twice-daily regimens suggested that a daily dosing regimen with double doses could be appropriate.