High variability of indinavir and nelfinavir pharmacokinetics in HIV-infected patients with a sustained virological response on highly active antiretroviral therapy

Can we Predict Drug Excretion into Saliva? A Systematic Review and Analysis of Physicochemical Properties

BACKGROUND AND OBJECTIVES

Saliva is a patient-friendly matrix for therapeutic drug monitoring (TDM) but is infrequently used in routine care. This is due to the uncertainty of saliva-based TDM results to inform dosing. This study aimed to retrieve data on saliva-plasma concentration and subsequently determine the physicochemical properties that influence the excretion of drugs into saliva to increase the foundational knowledge underpinning saliva-based TDM.

METHODS

Medline, Web of Science and Embase (1974-2023) were searched for human clinical studies, which determined drug pharmacokinetics in both saliva and plasma. Studies with at least ten subjects and five paired saliva-plasma concentrations per subject were included. For each study, the ratio of the area under the concentration-time curve between saliva and plasma was determined to assess excretion into saliva. Physicochemical properties of each drug (e.g. pKa, lipophilicity, molecular weight, polar surface area, rotatable bonds and fraction of drug unbound to plasma proteins) were obtained from PubChem and Drugbank. Drugs were categorised by their ionisability, after which saliva-to-plasma ratios were predicted with adjustment for protein binding and physiological pH via the Henderson-Hasselbalch equation. Spearman correlation analyses were performed for each drug category to identify factors predicting saliva excretion (α = 5%). Study quality was assessed by the risk of bias in non-randomised studies of interventions tool.

RESULTS

Overall, 42 studies including 40 drugs (anti-psychotics, anti-microbials, immunosuppressants, anti-thrombotic, anti-cancer and cardiac drugs) were included. The median saliva-to-plasma ratios were similar for drugs in the amphoteric (0.59), basic (0.43) and acidic (0.41) groups and lowest for drugs in the neutral group (0.21). Higher excretion of acidic drugs (n = 5) into saliva was associated with lower ionisation and protein binding (correlation between predicted versus observed saliva-to-plasma ratios: R2 = 0.85, p = 0.02). For basic drugs (n = 21), pKa predicted saliva excretion (Spearman correlation coefficient: R = 0.53, p = 0.02). For amphoteric drugs (n = 10), hydrogen bond donor (R = - 0.76, p = 0.01) and polar surface area (R = - 0.69, p = 0.02) were predictors. For neutral drugs (n = 10), protein binding (R = 0.84, p = 0.004), lipophilicity (R = - 0.65, p = 0.04) and hydrogen bond donor count (R = - 0.68, p = 0.03) were predictors. Drugs considered potentially suitable for saliva-based TDM are phenytoin, tacrolimus, voriconazole and lamotrigine. The studies had a low-to-moderate risk of bias.

CONCLUSIONS

Many commonly used drugs are excreted into saliva, which can be partly predicted by a drug's ionisation state, protein binding, lipophilicity, hydrogen bond donor count and polar surface area. The contribution of drug transporters and physiological factors to the excretion needs to be evaluated. Continued research on drugs potentially suitable for saliva-based TDM will aid in adopting this person-centred TDM approach to improve patient outcomes.

An inhibitor of HIV-1 protease modulates constitutive eIF2α dephosphorylation to trigger a specific integrated stress response

Inhibitors of the HIV aspartyl protease [HIV protease inhibitors (HIV-PIs)] are the cornerstone of treatment for HIV. Beyond their well-defined antiretroviral activity, these drugs have additional effects that modulate cell viability and homeostasis. However, little is known about the virus-independent pathways engaged by these molecules. Here we show that the HIV-PI Nelfinavir decreases translation rates and promotes a transcriptional program characteristic of the integrated stress response (ISR). Mice treated with Nelfinavir display hallmarks of this stress response in the liver, including α subunit of translation initiation factor 2 (eIF2α) phosphorylation, activating transcription factor-4 (ATF4) induction, and increased expression of known downstream targets. Mechanistically, Nelfinavir-mediated ISR bypassed direct activation of the eIF2α stress kinases and instead relied on the inhibition of the constitutive eIF2α dephosphorylation and down-regulation of the phophatase cofactor CReP (Constitutive Repressor of eIF2α Phosphorylation; also known as PPP1R15B). These findings demonstrate that the modulation of eIF2α-specific phosphatase cofactor activity can be a rheostat of cellular homeostasis that initiates a functional ISR and suggest that the HIV-PIs could be repositioned as therapeutics in human diseases to modulate translation rates and stress responses.

Differences in the Prediction of Area Under the Curve for a Protease Inhibitor Using Trough Versus Peak Concentration: Assessment Using Published Pharmacokinetic Data for Indinavir

In the present day antiretroviral therapy, Ctrough is a key tool for efficacy assessment. The present work explored the feasibility of using Ctrough or Cmax in the area under the concentration-time curve (AUC) prediction of indinavir. A simple unweighted linear regression model was developed to describe the relationship between Cmax versus AUC (r = 0.8101, P < 0.001) and Ctrough versus AUC (r = 0.8127, P < 0.001) for indinavir. The regression lines were used to predict the AUC values from literature Cmax or Ctrough data of indinavir in HIV and healthy subjects. The fold difference, defined as the quotient of the observed and predicted AUC values, was evaluated along with statistical comparison, including root mean square error (RMSE) prediction for the 2 models. The correlation between Cmax versus AUC and Ctrough versus AUC was established. Majority of the predicted values for Cmax versus AUC were within 0.75- to 1.5-fold differences. However, the Ctrough versus AUC model showed larger variability with approximately one-third of the predictions within 0.75- to 1.5-fold differences. The r value and %RMSE for observed versus predicted AUC for Ctrough (r = 0.5925, n = 65, P < 0.001, and RMSE: 67%) were inferior to the Cmax (r = 0.8773, n = 86, P < 0.001, and RMSE: 46%). In conclusion, Cmax versus AUC and Ctrough versus AUC relationships were established for indinavir showing the utility of a single concentration time point for therapeutic drug monitoring purpose. The Cmax model for indinavir may be more relevant for AUC prediction as determined by the statistical criteria.

Pharmacokinetic modelling of efavirenz, atazanavir, lamivudine and tenofovir in the female genital tract of HIV-infected pre-menopausal women

BACKGROUND AND OBJECTIVES

A previously published study of antiretroviral pharmacokinetics in the female genital tract of HIV-infected women demonstrated differing degrees of female genital tract penetration among antiretrovirals. These blood plasma (BP) and cervicovaginal fluid (CVF) data were co-modelled for four antiretrovirals with varying CVF exposures.

METHODS

Six paired BP and CVF samples were collected over 24 h, and antiretroviral concentrations determined using validated liquid chromatography (LC) with UV detection or LC-mass spectrometry analytical methods. For each antiretroviral, a BP model was fit using Bayesian estimation (ADAPT5), followed by addition of a CVF model. The final model was chosen based on graphical and statistical output, and then non-linear mixed-effects modelling using S-ADAPT was performed. Population mean parameters and their variability are reported. Model-predicated area under the concentration-time curve during the dosing interval (AUC(τ)) and exposure ratios of CVF AUC(τ):BP AUC(τ) were calculated for each drug.

RESULTS

The base model uses first-order absorption with a lag time, a two-compartment model, and a series of transit compartments that transfer the drug from BP to CVF. Protein-unbound drug transfers into CVF for efavirenz and atazanavir; total drug transfers for lamivudine and tenofovir. CVF follows a one-compartment model for efavirenz and atazanavir, and a two-compartment model for lamivudine and tenofovir. As expected, inter-individual variability was high. Model-predicted CVF AUC(τ):BP AUC(τ) ratios are consistent with published results.

CONCLUSIONS

This is the first pharmacokinetic modelling of antiretroviral disposition in BP and CVF. These models will be further refined with tissue data, and used in clinical trials simulations to inform future studies of HIV pre-exposure prophylaxis in women.

Clinical pharmacokinetics of antiretroviral drugs in older persons

INTRODUCTION

Combination antiretroviral therapy has enabled HIV-infected persons to reach older ages in high numbers. Hepatic and renal changes that normally occur with advancing age occur earlier and with higher incidence in HIV-infected individuals. A limited number of prospective controlled studies have demonstrated small reductions (17 to 41%) in lopinavir, atazanavir and lamivudine clearance in older versus younger adults. A much larger number of retrospective studies in adults (age range ∼ 20 to 60 years), including all antiretroviral drugs, have evaluated age as a covariate for pharmacokinetics. Most studies did not detect substantial associations between drug exposures and age.

AREAS COVERED

This review summarizes antiretroviral drug pharmacokinetics in older persons. The authors review articles from PubMed (search terms: elderly, antiretroviral, pharmacokinetics) in addition to the bibliographies of those selected.

EXPERT OPINION

The evidence to date does not support major pharmacokinetic changes in adults between ∼ 20 and 60 years of age. However, additional prospective, well-controlled studies are needed in more persons > 60 years, including those with frailty and comorbidities, with assessment of unbound drug clearance, and incorporation of adherence, pharmacogenetics and concomitant medications. Until then, guidelines for drug-drug interactions and dosing in renal and hepatic impairment should be followed in older HIV-infected individuals.

[Evidence-based therapeutic drug monitoring for indinavir]

The HIV protease inhibitor indinavir presents a wide inter-individual variability related to an intense hepatic metabolism. Published studies were analyzed to establish whether there is evidence that therapeutic drug monitoring of indinavir could improve patient care. It was reported that indinavir virological efficacy in HIV-infected patients with wild-type virus was significantly associated with trough concentrations > 100-150 ng/mL. Concerning the exposure-toxicity relationship, the risk of occurrence of nephrotoxicity was more frequently associated with trough concentrations > 500-1 000 ng/mL. Studies with concentration-controlled indinavir therapy suggest that therapeutic drug monitoring allows to achieve safe and effective concentrations, therefore, the level of evidence of the interest of indinavir therapeutic drug monitoring is highly recommended when indinavir is not associated with ritonavir and recommended when ritonavir is combined with ritonavir.

Population pharmacokinetic-pharmacogenetic study of nevirapine in HIV-infected Cambodian patients

The aims of this ANRS12154 open-label, single-center, multiple-dose pharmacokinetic study were to characterize nevirapine pharmacokinetics in a Cambodian population of HIV-infected patients and to identify environmental and genetic factors of variability, focusing on the CYP2B6, CYP3A5, and ABCB1 (MDR1) genes. A total of 170 Cambodian HIV-infected patients were included. Nevirapine trough concentrations were measured after 18 and 36 months of starting antiretroviral treatment and in samples drawn during a dosing interval in a subset of 10 patients. All data were analyzed by nonlinear mixed-effects modeling. The effect of covariates was investigated using the population pharmacokinetic model. Patients carrying homozygous loss-of-function alleles CYP3A5 6986A>G, CYP2B6 516G>T, CYP2B6 1459C>T, and ABCB1 3435C>T represent 42.4%, 9.2%, 0%, and 18% of the population, respectively. The median nevirapine trough concentrations did not differ after 18 and 36 months of treatment (5,705 ng/ml [range, ≤50 to 13,871] and 5,709 ng/ml [range, ≤50 to 15,422], respectively). Interpatient and intrapatient variabilities of nevirapine apparent clearance were 28% and 17%, respectively. CYP2B6 516G>T and creatinine clearance were found to significantly affect nevirapine apparent clearance. The estimated nevirapine apparent clearances were 2.95 liters/h, 2.62 liters/h, and 1.86 liters/h for CYP2B6 516GG, CYP2B6 516GT, and CYP2B6 516TT genotypes, respectively. The impact of creatinine clearance was small. This study demonstrates that 95% of the patients had sustained nevirapine exposure well above the 3,000-ng/ml threshold. Nevirapine clearance was shown to be affected by CYP2B6 516G>T genetic polymorphism and creatinine clearance, although this explained only part of the interpatient variability, which remains low compared to that for other antiretroviral drugs.

Nelfinavir, an HIV-1 protease inhibitor, induces oxidative stress-mediated, caspase-independent apoptosis in Leishmania amastigotes

Background

Visceral leishmaniasis has now emerged as an important opportunistic disease in patients coinfected with human immunodeficiency virus type-1 (HIV-1). Although the effectiveness of HIV-1 protease inhibitors, such as nelfinavir, in antiretroviral therapies is well documented, little is known of the impact of these drugs on Leishmania in coinfected individuals.

Methodology and Principal Findings

Here, we show that nelfinavir generates oxidative stress in the parasite, leading to altered physiological parameters such as an increase in the sub-G1 DNA content, nuclear DNA fragmentation and loss of mitochondrial potential, which are all characteristics of apoptosis. Pretreatment of axenic amastigotes with the caspase inhibitor z-VAD-fmk did not inhibit the increase in sub-G1 DNA content in nelfinavir-treated parasites, suggesting therefore that this antiviral agent does not kill Leishmania amastigotes in a caspase-dependent manner. Furthermore, we observed that the mitochondrial resident protein endonuclease G is involved. We also demonstrate that parasites overexpressing GSH1 (the rate limiting enzyme of glutathione biosynthesis) were more resistant to nelfinavir when compared to untransfected controls.

Conclusions and Significance

These data suggest that nelfinavir induces oxidative stress in Leishmania amastigotes, culminating in caspase-independent apoptosis, in which DNA is degraded by endonuclease G. This study provides a rationale for future, long-term design of new therapeutic strategies to test nelfinavir as a potential antileishmanial agent as well as for possible future use in Leishmania/HIV-1 coinfections.

Visceral leishmaniasis is the most severe form of disease caused by the parasite Leishmania. It is a major concern in South America, Africa, India and the Middle East. Additionally, it has now emerged as an important opportunistic disease in patients coinfected with HIV-1. This is due, in part, to the increasing overlap between urban centers and rural areas endemic for Leishmania. Although more efficient combinatorial antiviral drug regimens for treating HIV-1 infection have been developed, the impact of such therapies on HIV-1/Leishmania coinfection is yet to be explored. In this study, we investigated the effect of nelfinavir, a well-characterized anti-HIV-1 drug, on Leishmania. Treating the parasite with nelfinavir activates events that are hallmarks of programmed cell death (also called apoptosis). Among these are oxidative stress, changes in DNA replication and fragmentation, and release of mitochondrial enzymes. Furthermore, these events occur without the participation of caspases, which are classically linked to apoptosis; however, this atypical apoptosis requires the translocation of endonuclease G from mitochondria to the cytoplasm. These findings provide insights for the design of new anti-parasitic therapies, particularly in the case of Leishmania/HIV-1 coinfections.

Influence of pharmacogenetics on indinavir disposition and short-term response in HIV patients initiating HAART

AIMS

To assess the relationship between genetic polymorphisms and indinavir pharmacokinetic variability and to study the link between concentrations and short-term response or metabolic safety.

METHODS

Forty protease inhibitor-naive patients initiating highly active antiretroviral therapy (HAART) including indinavir/ritonavir and enrolled in the COPHAR 2-ANRS 111 trial were studied. At week 2, four blood samples were taken before and up to 6 h following drug intake. A population pharmacokinetic analysis was performed using the stochastic approximation expectation maximization (SAEM) algorithm implemented in MONOLIX software. The area under the concentration-time curve (AUC) and maximum (C(max)) and trough concentrations (C(trough)) of indinavir were derived from the population model and tested for their correlation with short-term viral response and safety measurements, while for ritonavir, these same three parameters were tested for their correlation with short-term biochemical safety

RESULTS

A one-compartment model with first-order absorption and elimination best described both indinavir and ritonavir concentrations. For indinavir, the estimated clearance and volume of distribution were 22.2 L/h and 97.3 L, respectively. The eight patients with the *1B/*1B genotype for the CYP3A4 gene showed a 70% decrease in absorption compared to those with the *1A/*1B or *1A/*1A genotypes (0.5 vs. 2.1, P = 0.04, likelihood ratio test by permutation). The indinavir AUC and C(trough) were positively correlated with the decrease in human immunodeficiency virus RNA between week 0 and week 2 (r = 0.4, P = 0.03 and r = -0.4, P = 0.03, respectively). Patients with the *1B/*1B genotype also had a significantly lower indinavir C(max) (median 3.6, range 2.1-5.2 ng/mL) than those with the *1A/*1B or *1A/*1A genotypes (median 4.4, range 2.2-8.3 ng/mL) (P = 0.04) and a lower increase in triglycerides during the first 4 weeks of treatment (median 0.1, range -0.7 to 1.4 vs. median 0.6, range -0.5 to 1.7 mmol/L, respectively; P = 0.02). For ritonavir, the estimated clearance and volume of distribution were 8.3 L/h and 60.7 L, respectively, and concentrations were not found to be correlated to biochemical safety. Indinavir and ritonavir absorption rate constants were found to be correlated, as well as their apparent volumes of distribution and clearances, indicating correlated bioavailability of the two drugs.

CONCLUSION

The CYP3A4*1B polymorphism was found to influence the pharmacokinetics of indinavir and, to some extent, the biochemical safety of indinavir.

Comparison of model-based tests and selection strategies to detect genetic polymorphisms influencing pharmacokinetic parameters

We evaluate by simulation three model-based methods to test the influence of a single nucleotide polymorphism on a pharmacokinetic parameter of a drug: analysis of variance (ANOVA) on the empirical Bayes estimates of the individual parameters, likelihood ratio test between models with and without genetic covariate, and Wald tests on the parameters of the model with covariate. Analyses are performed using the FO and FOCE method implemented in the NONMEM software. We compare several approaches for model selection based on tests and global criteria. We illustrate the results with pharmacokinetic data on indinavir from HIV-positive patients included in COPHAR 2-ANRS 111 to study the gene effect prospectively. Only the tests based on the EBE obtain an empirical type I error close to the expected 5%. The approximation made with the FO algorithm results in a significant inflation of the type I error of the LRT and Wald tests.

Predictive values of the human immunodeficiency virus phenotype and genotype and of amprenavir and lopinavir inhibitory quotients in heavily pretreated patients on a ritonavir-boosted dual-protease-inhibitor regimen

The inhibitory quotient (IQ) of human immunodeficiency virus (HIV) protease inhibitors (PIs), which is the ratio of drug concentration to viral susceptibility, is considered to be predictive of the virological response. We used several approaches to calculate the IQs of amprenavir and lopinavir in a subset of heavily pretreated patients participating in the French National Agency for AIDS Research (ANRS) 104 trial and then compared their potentials for predicting changes in the plasma HIV RNA level. Thirty-seven patients were randomly assigned to receive either amprenavir (600 mg twice a day [BID]) or lopinavir (400 mg BID) plus ritonavir (100 or 200 mg BID) for 2 weeks before combining the two PIs. The 90% inhibitory concentration (IC(90)) was measured using a recombinant assay without or with additional human serum (IC(90+serum)). Total and unbound PI concentrations in plasma were measured. Univariate linear regression was used to estimate the relation between the change in viral load and the IC(90) or IQ values. The amprenavir phenotypic IQ values were very similar when measured with the standard and protein binding-adjusted IC(90)s. No relationship was found between the viral load decline and the lopinavir IQ. During combination therapy, the amprenavir and lopinavir genotypic IQ values were predictive of the viral response at week 6 (P = 0.03). The number of protease mutations (< 5 or > or = 5) was related to the virological response throughout the study. These findings suggest that the combined genotypic IQ and the number of protease mutations are the best predictors of virological response. High amprenavir and lopinavir concentrations in these patients might explain why plasma concentrations and the phenotypic IQ have poor predictive value.

Population pharmacokinetic analysis of lamivudine, stavudine and zidovudine in controlled HIV-infected patients on HAART

OBJECTIVE

This work aimed at building a population pharmacokinetic (PK) model for lamivudine (LMV), stavudine (STV) and zidovudine (ZDV), estimating their inter and intraindividual PK variability and investigating the influence of different covariates.

METHODS

Population PK of LMV, STV and ZDV was separately evaluated from plasma concentrations obtained in 54, 39 and 27 HIV1-infected patients, respectively, enrolled in the COPHAR1-ANRS102 trial. The primary objective of this trial was to study the pharmacokinetics of indinavir (IDV) and nelfinavir (NFV) in treated patients with a sustained virological response. Concentrations of nucleoside analogs (NA) were measured in plasma as a secondary objective. A one-compartment model with first-order elimination was used, with zero-order absorption for LMV and first-order absorption for STV and ZDV.

RESULTS

Mean parameters [interpatient variability in coefficient of variation (CV%)] of LMV, STV and ZDV were: oral volume of distribution (V/F) 145 l (52%), 24 l (81%) and 248 l (80%), oral clearance (Cl/F) 32 l/h, 16 l/h (74%) and 124 l/h (51%), respectively. For LMV, absorption duration (Ta) was 1.46 h (64%). For STV and ZDV, ka was 0.46 h(-1) and 2.9 h(-1), respectively. We found a systematic effect of combination with NFV vs. IDV. We found that intrapatient variability was greater than interpatient variability (except for STV) and greater than 55% for the three drugs.

CONCLUSION

This trial enabled the estimation of the population PK parameters of three NA in patients with a sustained virological response, and the median curves could be used as references for concentration-controlled strategies. We observed, as for the protease inhibitors, a great variability of PK parameters.

Association between pharmacy medication refill-based adherence rates and cd4 count and viral-load responses: A retrospective analysis in treatment-experienced adults with HIV

BACKGROUND

Current guidelines and most contemporary statements in the literature indicate that, like other medical conditions, HIV infection requires exceptionally high adherence to highly active antiretroviral therapy (HAART) for successful treatment.

OBJECTIVE

This study was conducted to determine the association between pharmacy medication refill rates-a surrogate marker for adherence to HAART- and CD4-count/viral-load responses in patients with HIV METHODS: This retrospective study was conducted at the HIV Clinic, Veterans Affairs Medical Center, Durham, North Carolina. Male and female patients aged >/=18 years with a history of HIV who attended clinic appointments on 3 consecutive clinic days were enrolled. Pharmacy medication refill-based adherence rates over the 6 months before the study were determined by examining electronic pharmacy records. The most recent viral load and the change (Delta) in CD4 count over the past year-surrogate measures of outcome-were also collected from each patient's electronic medical record and compared with refill adherence rates. The incidence of AIDS-related events and past antiretroviral experience were also compared with the DeltaCD4 count and adherence rates.

RESULTS

Data from 58 patients were included in the study. Thirty-nine patients were black men; the mean age was 51.5 years. There was a nonsignificant correlation between 6-month pharmacy medication refillbased adherence rates and viral loads (r = 0.10). The relationship between DeltaCD4 count and adherence was complex. With adherence rates >70%, the DeltaCD4 count ranged from +414 to -238, with no indication that increasing adherence led to a greater CD4 count increase. The DeltaCD4 count progressively declined with adherence rates </=70%. Patients' past antiretroviral experience or incidence of AIDS-related events did not significantly affect the distribution of DeltaCD4 counts or adherence rates.

CONCLUSION

Based on our results, in patients with pharmacy medication refill-based adherence rates >70%, there was no significant correlation between adherence rates and DeltaCD4 counts or viral-load responses.

The SAEM algorithm for group comparison tests in longitudinal data analysis based on non-linear mixed-effects model

Non-linear mixed-effects models (NLMEMs) are used to improve information gathering from longitudinal studies and are applied to treatment evaluation in disease-evolution studies, such as human immunodeficiency virus (HIV) infection. The estimation of parameters and the statistical tests are critical issues in NLMEMs since the likelihood and the Fisher information matrix have no closed form. An alternative method to numerical integrations, in which convergence is slow, and to methods based on linearization, in which asymptotic convergence has not been proved, is the Stochastic Approximation Expectation-Maximization (SAEM) algorithm. For the Wald test and the likelihood ratio test, we propose estimating the Fisher information matrix by stochastic approximation and the likelihood by importance sampling. We evaluate these SAEM-based tests in a simulation study in the context of HIV viral load decrease after initiation of an antiretroviral treatment. The results from this simulation illustrate the theoretical convergence properties of SAEM. We also propose a method based on the SAEM algorithm to compute the minimum sample size required to perform a Wald test of a given power for a covariate effect in NLMEMs. Lastly, we illustrate these tests on the evaluation of the effect of ritonavir on the indinavir pharmacokinetics in HIV patients and compare the results with those obtained using the adaptative Gaussian quadrature method implemented in the SAS procedure NLMIXED.