mrgsolve is an open‐source R package available on the Comprehensive R Archive Network. It combines R and C++ coding for simulation from hierarchical, ordinary differential equation–based models. Its efficient simulation engine and integration into a parallelizable, R‐based workflow makes mrgsolve a convenient tool both for simple and complex models and thus is ideal for physiologically‐based pharmacokinetic (PBPK) and quantitative systems pharmacology (QSP) model. This tutorial will first introduce the basics of the mrgsolve simulation workflow, including model specification, the introduction of interventions (dosing events) into the simulation, and simulated results postprocessing. An applied simulation example is then presented using a PBPK model for voriconazole, including a model validation step against adult and pediatric data sets. A final simulation example is then presented using a previously published QSP model for mitogen‐activated protein kinase signaling in colorectal cancer, illustrating population simulation of different combination therapies.

IntroductionThe aim of the analysis was to characterize the population pharmacokinetics (PKs) and exposure–response (E–R) for efficacy (fasting plasma glucose, glycated hemoglobin) and safety/tolerability [hypoglycemia, genital infections, urinary tract infection (UTI), and volume depletion] of the sodium glucose cotransporter 2 inhibitor, empagliflozin, in patients with type 2 diabetes mellitus. This study extends the findings of previous analyses which described the PK and pharmacodynamics (PD) using early clinical studies of up to 12 weeks in duration.MethodsPopulation pharmacokinetic and E–R models were developed based on two Phase I, four Phase II, and four Phase III studies.ResultsVariability in empagliflozin exposure was primarily affected by estimated glomerular filtration rate (eGFR) (less than twofold increase in exposure in patients with severe renal impairment). Consistent with its mode of action, the efficacy of empagliflozin was increased with elevated baseline plasma glucose levels and attenuated with decreasing renal function, but was still maintained to nearly half the maximal effect with eGFR as low as 30 mL/min/1.73 m2. All other investigated covariates, including sex, body mass index, race, and age did not alter the PK or efficacy of empagliflozin to a clinically relevant extent. Compared with placebo, empagliflozin administration was associated with an exposure-independent increase in the incidence of genital infections and no significant change in the risk of UTI, hypoglycemia, or volume depletion.ConclusionBased on the results from the PK and E–R analysis, no dose adjustment is required for empagliflozin in the patient population for which the drug is approved.FundingBoehringer Ingelheim.Electronic supplementary materialThe online version of this article (doi:10.1007/s13300-016-0174-y) contains supplementary material, which is available to authorized users.

The goal of this work was to extend a mathematical, multiscale systems model of bone function, remodeling, and health in order to explore hypotheses related to therapeutic modulation of sclerostin and quantitatively describe purported osteocyte activity within bone remodeling events. A pharmacokinetic model with first-order absorption and dual elimination pathways was used to describe the kinetics of romosozumab, a monoclonal antibody (mAb) against sclerostin. To describe total circulating sclerostin, an extended indirect response model of inhibition of offset was developed. These models were subsequently linked to the systems model, with sclerostin signaling changes in resorption and formation through established osteocyte-mediated mechanisms. The model proposes relative contributions of the osteocyte to the RANKL pool, a major player in feedback signaling, and is used to explore hypotheses surrounding attenuation of anabolic activity after multiple doses of sclerostin mAbs, a phenomenon whose mechanism is poorly understood.

To identify linzagolix doses, an oral GnRH receptor antagonist, that effectively lower oestradiol (E2) to relieve endometriosis-related pelvic pain without compromising bone health.Methods: Integrated statistical, pharmacokinetic-pharmacodynamic and systems pharmacology models were developed from Phase 1 and 2 clinical trial data in healthy volunteers and patients, receiving linzagolix 25-200 mg daily or placebo, and analysed simultaneously. The main outcome measures were pelvic pain scores for dysmenorrhoea, nonmenstrual pelvic pain (NMPP), uterine bleeding and lumbar spine bone mineral density (BMD).Results: Linzagolix pharmacokinetics were described by a 2-compartment model with sequential zero/first-order absorption process (CL/F: 0.422 L/h). E2 changes over time were well described as a function of linzagolix 24-hour AUC (AUC 50 : 1.68 Â 10 5 ng h/mL). For a Caucasian reference patient, a change in E2 from 50-20 pg/mL at 24 weeks increased the odds of relief of dysmenorrhoea 1.33-fold and NMPP 1.07-fold (95% CI: 1.22-1.47 and 1.02-1.12, respectively) and decreased bleeding days by 1.55 (95% CI: 1.39-1.72). A previously validated quantitative systems pharmacology BMD model was adjusted to the clinical data. The mean week 24 lumbar spine BMD change from baseline ranged from À0.092% in the 50 mg dose, À1.30% in the 100 mg dose group and À2.67% in the 200 mg dose group.Discussion: The previously-reported E2 target range (20-50 pg/mL) to balance efficacy and safety endpoints was confirmed. Linzagolix once daily doses between 75-125 mg daily were expected to meet endometriosis-associated pain, efficacy, and BMD loss targets in Caucasian patients.

No abstract