RPEM: Randomized Monte Carlo parametric expectation maximization algorithm.

Inspired from quantum Monte Carlo, by sampling discrete and continuous variables at the same time using the Metropolis-Hastings algorithm, we present a novel, fast, and accurate high performance Monte Carlo Parametric Expectation Maximization (MCPEM) algorithm. We named it Randomized Parametric Expectation Maximization (RPEM). We compared RPEM with NONMEM's Importance Sampling Method (IMP), Monolix's Stochastic Approximation Expectation Maximization (SAEM), and Certara's Quasi-Random Parametric Expectation Maximization (QRPEM) for a realistic two-compartment voriconazole model with ordinary differential equations using simulated data. We show that RPEM is as fast and as accurate as the algorithms IMP, QRPEM, and SAEM for the voriconazole model in reconstructing the population parameters, for the normal and log-normal cases.

pyDarwin: A Machine Learning Enhanced Automated Nonlinear Mixed-Effect Model Selection Toolbox.

pyDarwin is an open-source Python package for nonlinear mixed-effect model selection. pyDarwin combines machine-learning algorithms and NONMEM to perform a global search for the optimal model in a user-defined model search space. Compared with traditional stepwise search, pyDarwin provides an efficient platform for conducting an objective, robust, less labor-intensive model selection process without compromising model interpretability. In this tutorial, we will begin by introducing the essential components and concepts within the package. Subsequently, we will provide an overview of the pyDarwin modeling workflow and the necessary files needed for model selection. To illustrate the entire process, we will conclude with an example utilizing quetiapine clinical data.

Respiratory syncytial virus-A dynamics and the effects of lumicitabine, a nucleoside viral replication inhibitor, in experimentally infected humans.

Background: Respiratory syncytial virus (RSV) causes high morbidity, with mortality rates approaching or exceeding that of influenza in adult and infant patient populations, respectively. Lumicitabine (ALS-008176 or JNJ-64041575) is an oral nucleoside analogue prodrug in clinical development to treat RSV infections. This prodrug converts to plasma-circulating ALS-8112, and then to the 5'-active nucleoside triphosphate (NTP) form within host cells. We conducted an RSV-A challenge study in healthy adults to evaluate lumicitabine's activity during an active RSV infection. Objectives: To develop a semi-mechanistic mathematical model describing RSV kinetics, and the pharmacokinetics (PK) and pharmacodynamics (PD) of lumicitabine during treatment. Methods: Nasopharyngeal viral load and concentrations of ALS-8112 and ALS-8144 (uridine metabolite) were measured frequently over the study duration. Population viral kinetic and PK/PD models were developed using NONMEM. The RSV life-cycle was described using a target-cell-limited model that included a physiological delay. Results: The estimated clearances of ALS-8112 and ALS-8144 were 54.2 and 115 L/h/70 kg, respectively. A semi-physiological model was linked to predict ALS-8112 conversion to active intracellular NTP. Extensive and rapid RSV reduction occurred after lumicitabine treatment (EC50 = 1.79 µM), with >99% viral inhibition at 2 h after loading dose. Simulated NTP exposures and time to EC50 attainment suggested that rapid therapeutic effects and reduced dosing frequency are achievable in adult and paediatric patients. Conclusions: The semi-mechanistic model characterizes RSV kinetics and the antiviral effectiveness of lumicitabine in an adult challenge population. This model is applicable to guide dose selection in adult and paediatric patients.

A Regression Approach to Visual Predictive Checks for Population Pharmacometric Models.

A visual predictive check (VPC) is a common diagnostic procedure for population pharmacometric models. Typically, VPCs are generated by specifying intervals, or "bins", of an independent variable (e.g., time). However, bin specification is not always straightforward and the choice of bins may affect the appearance, and possibly conclusions, of VPCs. The objective of this work was to demonstrate how regression techniques can be used to derive VPCs and prediction-corrected VPCs (pcVPCs) for population pharmacometric models. This alternative approach negates the need for empirical bin selection. The proposed method utilizes local and additive quantile regression. Implementation is straightforward and computationally acceptable. This work provides support for deriving VPCs and pcVPCs via regression techniques.

Interdisciplinary pharmacometrics linking oseltamivir pharmacology, influenza epidemiology and health economics to inform antiviral use in pandemics.

AIMS: A modular interdisciplinary platform was developed to investigate the economic impact of oseltamivir treatment by dosage regimen under simulated influenza pandemic scenarios. METHODS: The pharmacology module consisted of a pharmacokinetic distribution of oseltamivir carboxylate daily area under the concentration-time curve at steady state (simulated for 75 mg and 150 mg twice daily regimens for 5 days) and a pharmacodynamic distribution of viral shedding duration obtained from phase II influenza inoculation data. The epidemiological module comprised a susceptible, exposed, infected, recovered (SEIR) model to which drug effect on the basic reproductive number (R0 ), a measure of transmissibility, was linked by reduction of viral shedding duration. The number of infected patients per population of 100 000 susceptible individuals was simulated for a series of pandemic scenarios, varying oseltamivir dose, R0 (1.9 vs. 2.7), and drug uptake (25%, 50%, and 80%). The number of infected patients for each scenario was entered into the health economics module, a decision analytic model populated with branch probabilities, disease utility, costs of hospitalized patients developing complications, and case-fatality rates. Change in quality-adjusted life years was determined relative to base case. RESULTS: Oseltamivir 75 mg relative to no treatment reduced the median number of infected patients, increased change in quality-adjusted life years by deaths averted, and was cost-saving under all scenarios; 150 mg relative to 75 mg was not cost effective in low transmissibility scenarios but was cost saving in high transmissibility scenarios. CONCLUSION: This methodological study demonstrates proof of concept that the disciplines of pharmacology, disease epidemiology and health economics can be linked in a single quantitative framework.

Absorption, distribution, metabolism and excretion (ADME) of the ALK inhibitor alectinib: results from an absolute bioavailability and mass balance study in healthy subjects.

1. Alectinib is a highly selective, central nervous system-active small molecule anaplastic lymphoma kinase inhibitor. 2. The absolute bioavailability, metabolism, excretion and pharmacokinetics of alectinib were studied in a two-period single-sequence crossover study. A 50 µg radiolabelled intravenous microdose of alectinib was co-administered with a single 600 mg oral dose of alectinib in the first period, and a single 600 mg/67 µCi oral dose of radiolabelled alectinib was administered in the second period to six healthy male subjects. 3. The absolute bioavailability of alectinib was moderate at 36.9%. Geometric mean clearance was 34.5 L/h, volume of distribution was 475 L and the hepatic extraction ratio was low (0.14). 4. Near-complete recovery of administered radioactivity was achieved within 168 h post-dose (98.2%) with excretion predominantly in faeces (97.8%) and negligible excretion in urine (0.456%). Alectinib and its major active metabolite, M4, were the main components in plasma, accounting for 76% of total plasma radioactivity. In faeces, 84% of dose was excreted as unchanged alectinib with metabolites M4, M1a/b and M6 contributing to 5.8%, 7.2% and 0.2% of dose, respectively. 5. This novel study design characterised the full absorption, distribution, metabolism and excretion properties in each subject, providing insight into alectinib absorption and disposition in humans.

Comparative benefit of malaria chemoprophylaxis modelled in United Kingdom travellers.

BACKGROUND: Chemoprophylaxis against falciparum malaria is recommended for travellers from non-endemic countries to malarious destinations, but debate continues on benefit, especially with regard to mefloquine. Quantification of benefit for travellers from the United Kingdom (UK) was modelled to assist clinical and public health decision making. METHODS: The model was constructed utilising: World Tourism Organization data showing total number of arrivals from the UK in countries with moderate or high malaria risk; data from a retrospective UK Clinical Practice Research Datalink (CPRD) drug utilisation study; additional information on chemoprophylaxis, case fatality and tolerability were derived from the travel medicine literature. Chemoprophylaxis with the following agents was considered: atovaquone-proguanil (AP), chloroquine with and without proguanil (C ± P), doxycycline (Dx), mefloquine (Mq). The model was validated for the most recent year with temporally matched datasets for UK travel destinations and imported malaria (2007) against UK Health Protection Agency data on imported malaria. RESULTS: The median (mean) duration of chemoprophylaxis for each agent in weeks (CPRD) was: AP 3.3 (3.5), C ± P 9 (12.1), Dx 8 (10.3), Mq 9 (12.3): the maximum duration of use of all regimens was 52 weeks. The model correctly predicted falciparum malaria deaths and gave a robust estimate of total cases--model: 5 deaths from 1118 cases; UK Health Protection Agency: 5 deaths from 1153 cases. The number needed to take chemoprophylaxis (NNP) to prevent a case of malaria considered against the 'background' reported incidence in non-users of chemoprophylaxis deemed in need of chemoprophylaxis was: C ± P 272, Dx 269, Mq 260, AP 252; the NNP to prevent a UK traveller malaria death was: C ± P 62613, Dx 61923, Mq 59973, AP 58059; increasing the 'background' rate by 50% yielded NNPs of: C ± P 176, Dx 175, Mq 171, AP 168. The impact of substituting atovaquone-proguanil for all mefloquine usage resulted in a 2.3% decrease in estimated infections. The number of travellers experiencing moderate adverse events (AE) or those requiring medical attention or drug withdrawal per case prevented is as follows: C ± P 170, Mq 146, Dx 114, AP 103. CONCLUSIONS: The model correctly predicted the number of malaria deaths, providing a robust and reliable estimate of the number of imported malaria cases in the UK, and giving a measure of benefit derived from chemoprophylaxis use against the likely adverse events generated. Overall numbers needed to prevent a malaria infection are comparable among the four options and are sensitive to changes in the background infection rates. Only a limited impact on the number of infections can be expected if Mq is substituted by AP.

Safety and pharmacokinetics of daclizumab in pediatric renal transplant recipients.

This study examined the safety and pharmacokinetics/pharmacodynamics of daclizumab in combination with mycophenolate mofetil (or azathioprine), corticosteroids, and cyclosporine or tacrolimus, in 61 pediatric renal allograft recipients in three age groups: less than or equal to five yr (n = 18), 6-12 yr (n = 18), and 13-17 yr (n = 25). The dosing regimen was daclizumab 1.0 mg/kg before transplantation, followed by four biweekly doses. The pharmacokinetics of daclizumab were described using NONMEM software. Median (range) estimated trough daclizumab levels achieved on day 56 (before dose 5) were 3.88 microg/mL (2.48-8.78), 4.54 microg/mL (1.79-18.7), and 4.94 microg/mL (0.05-10.6) in the less than or equal to five yr (n = 15), 6-12 yr (n = 17), and 13-17 yr (n = 22) age groups, respectively. Steady-state median (range) daclizumab exposures were 2040 mg x h/mL (1585-3778), 2757 mg x h/mL (1873-3494) and 3297 mg x h/mL (1705-6453), respectively. Saturation of the IL-2R occurred rapidly and was maintained for greater than or equal to three months after transplantation. Daclizumab was generally well-tolerated with no acute allergic or anaphylactic reactions, deaths or malignancies during the study. The proportion of patients who developed acute rejection at six and 12 months was 8.5% and 16.7%, respectively. This study shows that adding daclizumab at 1 mg/kg to standard immunosuppressive therapy provides safe and effective IL-2R blockade.

Population pharmacokinetics and exposure-response relationship of enfuvirtide in treatment-experienced human immunodeficiency virus type 1-infected patients.

OBJECTIVE: Our objective was to characterize population pharmacokinetics of enfuvirtide, 90 mg twice daily injected subcutaneously, in treatment-experienced human immunodeficiency virus type 1 (HIV-1)-infected patients, as well as the relationship between exposure and antiviral effect. METHODS: Plasma concentrations of enfuvirtide and HIV-1 ribonucleic acid were obtained from 628 patients in 2 phase III studies. NONMEM software was used for population pharmacokinetic analysis and to assess the effects of age, gender, body weight, anti-gp41 antibodies, and concomitant drugs. Enfuvirtide exposure (area under the plasma concentration-12-hour time curve or steady-state trough concentration) was calculated from individual parameter estimates derived from the model. The decline in HIV-1 ribonucleic acid from baseline at week 2 or 24 was regressed against estimates of enfuvirtide exposure by a maximum effect model. The exposure-response relationship was examined in functional monotherapy (phenotypic sensitivity score of 0) and combination therapy (phenotypic sensitivity score > or = 1). RESULTS: Enfuvirtide population pharmacokinetics was well described by a 1-compartment model with first-order absorption and elimination. Body weight and female gender were identified as affecting apparent clearance but not efficacy and safety. Concomitant medications had no significant effect on enfuvirtide pharmacokinetics. Antiviral response to enfuvirtide was independent of drug exposure, suggesting that the approved 90-mg twice-daily dose was in the plateau portion of the dose-response curve. For functional monotherapy (phenotypic sensitivity score of 0), approximately 66% of estimated maximal effect was achieved at week 2 and 73% at week 24, and for combination therapy, more than 92% was achieved at both weeks 2 and 24. CONCLUSIONS: Body weight and gender affected enfuvirtide clearance, but changes in exposure did not affect efficacy or safety. Efficacy reached a plateau at the 90-mg twice-daily dosage in the exposure-response curve.

Pharmacokinetics, pharmacodynamics and drug interaction potential of enfuvirtide.

Enfuvirtide, the first fusion inhibitor approved for the treatment of HIV-1 infection, is a synthetic peptide that binds to HIV-1 glycoprotein 41, blocking the fusion of viral and cellular membranes. When administered subcutaneously at the recommended dose of 90 mg twice daily with optimised background antiretroviral therapy, enfuvirtide significantly reduces plasma HIV-1 RNA levels up to 48 weeks compared with optimised background therapy alone. Enfuvirtide exhibits a small volume of distribution (5.48 L), low systemic clearance (1.4 L/h) and high plasma protein binding (92%). Less than 17% of enfuvirtide is converted to a minimally active deaminated form of the parent drug. Both enfuvirtide and its metabolite are primarily eliminated via catabolism to amino acid residues. Following subcutaneous administration, enfuvirtide is almost completely absorbed, and exposure increases almost linearly with dose over the range 45-180 mg. When administered at the recommended dose in adults, subcutaneous absorption is slow and protracted, resulting in relatively flat steady-state plasma concentration-time profiles. Bioavailability is high (84.3%) and the elimination half-life (3.8 hours) supports twice-daily administration. Comparable absorption was observed from three different anatomical injection sites. The pharmacokinetic-pharmacodynamic relationship indicates that the recommended dose, in combination with other active antiretrovirals, is optimal. Enfuvirtide clearance is influenced to a small extent by sex and bodyweight but this does not necessitate dosage adjustment. In vitro and in vivo studies indicate that enfuvirtide has a low potential to interact with concomitantly administered drugs. Enfuvirtide did not influence concentrations of drugs metabolised by cytochrome P450 (CYP) 3A4, CYP2D6 or N-acetyltransferase, and had only minimal effects on those metabolised by CYP1A2, CYP2E1 or CYP2C19. Coadministration of ritonavir, ritonavir-boosted saquinavir or rifampicin (rifampin) did not result in clinically significant changes in enfuvirtide pharmacokinetics. In HIV-1-infected paediatric patients, subcutaneous dosages based on bodyweight (2 mg/kg twice daily) produce pharmacokinetics broadly similar to those observed in adults administered 90 mg twice daily.

Pharmacokinetics of enfuvirtide in a patient with impaired renal function.

Enfuvirtide is the first member of the fourth class of antiretroviral drugs to become available. We present a case report of a human immunodeficiency virus (HIV)-positive patient with mild renal impairment who developed severe renal impairment secondary to tenofovir therapy while receiving enfuvirtide. Because the patient's pharmacokinetic profile was not significantly altered, compared with that of HIV-infected patients with normal renal function, no dose modifications were required.

Lack of enzyme-inducing effect of rifampicin on the pharmacokinetics of enfuvirtide.

The primary objective was to determine whether rifampicin influences the pharmacokinetics of enfuvirtide in HIV-1-infected patients. In a single-center, open-label, one-sequence crossover, clinical pharmacology study, 12 HIV-1-infected adults received enfuvirtide (90 mg, twice daily) on days 1 to 3 and days 11 to 13 (morning dose only on days 3 and 13) and rifampicin (600 mg, once daily) from days 4 to 13. Plasma concentrations were measured for enfuvirtide and its metabolite (days 3 and 13) and rifampicin (day 13 only). The ratios of least squares means (LSM) and 90% confidence intervals for enfuvirtide and enfuvirtide metabolite pharmacokinetic parameters (AUC12h, Cmax, Ctrough) were estimated in the presence and absence of rifampicin. Treatments were compared using an analysis of variance for natural log-transformed variables, with factors patient and treatment. Efficacy and safety were also monitored. Steady-state rifampicin had no appreciable effect on any of the pharmacokinetic parameters assessed for either enfuvirtide or its metabolite. The ratio of LSM for AUC12h, Cmax, and Ctrough for enfuvirtide was 97.5%, 103%, and 84.9%, respectively, and 108%, 112%, and 92.9%, for the enfuvirtide metabolite. Rifampicin did not affect the t1/2 of enfuvirtide or its metabolite. There were no unexpected effects of rifampicin on the short-term antiviral effect or safety of the administered antiretroviral treatment. The pharmacokinetics of enfuvirtide are not induced by a 10-day pretreatment with rifampicin.

Pharmacokinetics of daclizumab and mycophenolate mofetil with cyclosporine and steroids in renal transplantation.

Daclizumab and mycophenolate mofetil (MMF) decrease the incidence of acute allograft rejection. This double-blind, randomized, placebo-controlled trial was performed primarily to assess the pharmacokinetics of MMF in an immunosuppressive regimen incorporating daclizumab. At five centers, 75 renal transplant recipients were randomized 2:1 to receive either daclizumab 1 mg/kg or placebo pre-transplantation and every other week, for a total of five doses. All patients received cyclosporine, steroids, and MMF. Levels of mycophenolic acid (MPA), its glucuronide metabolite, and daclizumab were measured after dosing on days 28 and 56. Safety parameters evaluated included: adverse events, laboratory abnormalities, infections, patient/graft survival, incidence of lymphoproliferative disorders, and incidence of acute rejection at 12 months. The concomitant administration of daclizumab and MMF had no effect on the pharmacokinetics of MPA: AUC(0-8) values (microg h/mL +/- SD) on day 28 were 30.1 +/- 13.3 for daclizumab-treat patients vs. 31.1 +/- 12.4 for placebo and on day 56, 37.7 +/- 18.2 for daclizumab-treated patients vs. 35.7 +/- 14.0 for placebo. Adverse events were similar between the two groups. Acute rejection at 12 months occurred in 14% of patients receiving daclizumab and 20% of patients receiving placebo. The coadministration of daclizumab did not result in a pharmacokinetic interaction with MPA, the active metabolite of MMF.

Pharmacokinetics of plasma enfuvirtide after subcutaneous administration to patients with human immunodeficiency virus: Inverse Gaussian density absorption and 2-compartment disposition.

OBJECTIVE: Enfuvirtide (T-20) is the first of a novel class of human immunodeficiency virus (HIV) drugs that block gp41-mediated viral fusion to host cells. The objectives of this study were to develop a structural pharmacokinetic model that would adequately characterize the absorption and disposition of enfuvirtide pharmacokinetics after both intravenous and subcutaneous administration and to evaluate the dose proportionality of enfuvirtide pharmacokinetic parameters at a subcutaneous dose higher than that currently used in phase III studies. METHODS: Twelve patients with HIV infection received 4 single doses of enfuvirtide separated by a 1-week washout period in an open-label, randomized, 4-way crossover fashion. The doses studied were 90 mg (intravenous) and 45 mg, 90 mg, and 180 mg (subcutaneous). Serial blood samples were collected up to 48 hours after each dose. Plasma enfuvirtide concentrations were measured with use of a validated liquid chromatography-tandem mass spectrometry method. RESULTS: Enfuvirtide plasma concentration-time data after subcutaneous administration were well described by an inverse Gaussian density function-input model linked to a 2-compartment open distribution model with first-order elimination from the central compartment. The model-derived mean pharmacokinetic parameters (+/-SD) were volume of distribution of the central compartment (3.8 +/- 0.8 L), volume of distribution of the peripheral compartment (1.7 +/- 0.6 L), total clearance (1.44 +/- 0.30 L/h), intercompartmental distribution (2.3 +/- 1.1 L/h), bioavailability (89% +/- 11%), and mean absorption time (7.26 hours, 8.65 hours, and 9.79 hours for the 45-mg, 90-mg, and 180-mg dose groups, respectively). The terminal half-life increased from 3.46 to 4.35 hours for the subcutaneous dose range from 45 to 180 mg. CONCLUSIONS: An inverse Gaussian density function-input model linked to a 2-compartment open distribution model with first-order elimination from the central compartment was appropriate to describe complex absorption and disposition kinetics of enfuvirtide plasma concentration-time data after subcutaneous administration to patients with HIV infection. Enfuvirtide was nearly completely absorbed from subcutaneous depot, and pharmacokinetic parameters were linear up to a dose of 180 mg in this study.