A generic whole body physiologically based pharmacokinetic model for therapeutic proteins in PK-Sim.

Proteins are an increasingly important class of drugs used as therapeutic as well as diagnostic agents. A generic physiologically based pharmacokinetic (PBPK) model was developed in order to represent at whole body level the fundamental mechanisms driving the distribution and clearance of large molecules like therapeutic proteins. The model was built as an extension of the PK-Sim model for small molecules incorporating (i) the two-pore formalism for drug extravasation from blood plasma to interstitial space, (ii) lymph flow, (iii) endosomal clearance and (iv) protection from endosomal clearance by neonatal Fc receptor (FcRn) mediated recycling as especially relevant for antibodies. For model development and evaluation, PK data was used for compounds with a wide range of solute radii. The model supports the integration of knowledge gained during all development phases of therapeutic proteins, enables translation from pre-clinical species to human and allows predictions of tissue concentration profiles which are of relevance for the analysis of on-target pharmacodynamic effects as well as off-target toxicity. The current implementation of the model replaces the generic protein PBPK model available in PK-Sim since version 4.2 and becomes part of the Open Systems Pharmacology Suite.

Contribution of rivaroxaban to the international normalized ratio when switching to warfarin for anticoagulation as determined by simulation studies.

AIM: This study evaluated the influence of rivaroxaban 20 mg once daily on international normalized ratio (INR) during the co-administration period when switching from rivaroxaban to warfarin. METHODS: We developed a calibrated coagulation model that was qualified with phase I clinical data. Prothrombin time and INR values were simulated by use of phospholipid concentrations that matched Neoplastin Plus® and Innovin® reagents. To simulate the combined effects of rivaroxaban and warfarin on INR during switching, warfarin initiation was simulated by adjusting the magnitude of the warfarin effect to reach the desired target INRs over the course of 21 days. The warfarin effect values (obtained every 6 h) and the desired rivaroxaban plasma concentrations were used. Nomograms were generated from rivaroxaban induced increases in INR. RESULTS: The simulation had good prediction quality. Rivaroxaban induced increases in the total INR from the warfarin attributed INR were seen, which increased with rivaroxaban plasma concentration. When the warfarin only INR was 2.0-3.0, the INR contribution of rivaroxaban with Neoplastin Plus® was 0.5-1.2, decreasing to 0.3-0.6 with Innovin® at median trough rivaroxaban plasma concentrations (38 µg l(-1) ). CONCLUSIONS: The data indicate that measuring warfarin induced changes in INR are best performed at trough rivaroxaban concentrations (24 h after rivaroxaban dosing) during the co-administration period when switching from rivaroxaban to warfarin. Furthermore, Innovin® is preferable to Neoplastin Plus® because of its substantially lower sensitivity to rivaroxaban, thereby reducing the influence of rivaroxaban on the measured INR.

Computational investigation of potential dosing schedules for a switch of medication from warfarin to rivaroxaban-an oral, direct Factor Xa inhibitor.

The long-lasting anticoagulant effect of vitamin K antagonists can be problematic in cases of adverse drug reactions or when patients are switched to another anticoagulant therapy. The objective of this study was to examine in silico the anticoagulant effect of rivaroxaban, an oral, direct Factor Xa inhibitor, combined with the residual effect of discontinued warfarin. Our simulations were based on the recommended anticoagulant dosing regimen for stroke prevention in patients with atrial fibrillation. The effects of the combination of discontinued warfarin plus rivaroxaban were simulated using an extended version of a previously validated blood coagulation computer model. A strong synergistic effect of the two distinct mechanisms of action was observed in the first 2-3 days after warfarin discontinuation; thereafter, the effect was close to additive. Nomograms for the introduction of rivaroxaban therapy after warfarin discontinuation were derived for Caucasian and Japanese patients using safety and efficacy criteria described previously, together with the coagulation model. The findings of our study provide a mechanistic pharmacologic rationale for dosing schedules during the therapy switch from warfarin to rivaroxaban and support the switching strategies as outlined in the Summary of Product Characteristics and Prescribing Information for rivaroxaban.

Development of a paediatric population-based model of the pharmacokinetics of rivaroxaban.

BACKGROUND: Venous thromboembolism has been increasingly recognised as a clinical problem in the paediatric population. Guideline recommendations for antithrombotic therapy in paediatric patients are based mainly on extrapolation from adult clinical trial data, owing to the limited number of clinical trials in paediatric populations. The oral, direct Factor Xa inhibitor rivaroxaban has been approved in adult patients for several thromboembolic disorders, and its well-defined pharmacokinetic and pharmacodynamic characteristics and efficacy and safety profiles in adults warrant further investigation of this agent in the paediatric population. OBJECTIVE: The objective of this study was to develop and qualify a physiologically based pharmacokinetic (PBPK) model for rivaroxaban doses of 10 and 20 mg in adults and to scale this model to the paediatric population (0-18 years) to inform the dosing regimen for a clinical study of rivaroxaban in paediatric patients. METHODS: Experimental data sets from phase I studies supported the development and qualification of an adult PBPK model. This adult PBPK model was then scaled to the paediatric population by including anthropometric and physiological information, age-dependent clearance and age-dependent protein binding. The pharmacokinetic properties of rivaroxaban in virtual populations of children were simulated for two body weight-related dosing regimens equivalent to 10 and 20 mg once daily in adults. The quality of the model was judged by means of a visual predictive check. Subsequently, paediatric simulations of the area under the plasma concentration-time curve (AUC), maximum (peak) plasma drug concentration (C max) and concentration in plasma after 24 h (C 24h) were compared with the adult reference simulations. RESULTS: Simulations for AUC, C max and C 24h throughout the investigated age range largely overlapped with values obtained for the corresponding dose in the adult reference simulation for both body weight-related dosing regimens. However, pharmacokinetic values in infants and preschool children (body weight <40 kg) were lower than the 90 % confidence interval threshold of the adult reference model and, therefore, indicated that doses in these groups may need to be increased to achieve the same plasma levels as in adults. For children with body weight between 40 and 70 kg, simulated plasma pharmacokinetic parameters (C max, C 24h and AUC) overlapped with the values obtained in the corresponding adult reference simulation, indicating that body weight-related exposure was similar between these children and adults. In adolescents of >70 kg body weight, the simulated 90 % prediction interval values of AUC and C 24h were much higher than the 90 % confidence interval of the adult reference population, owing to the weight-based simulation approach, but for these patients rivaroxaban would be administered at adult fixed doses of 10 and 20 mg. CONCLUSION: The paediatric PBPK model developed here allowed an exploratory analysis of the pharmacokinetics of rivaroxaban in children to inform the dosing regimen for a clinical study in paediatric patients.

Evaluation of the efficacy and safety of rivaroxaban using a computer model for blood coagulation.

Rivaroxaban is an oral, direct Factor Xa inhibitor approved in the European Union and several other countries for the prevention of venous thromboembolism in adult patients undergoing elective hip or knee replacement surgery and is in advanced clinical development for the treatment of thromboembolic disorders. Its mechanism of action is antithrombin independent and differs from that of other anticoagulants, such as warfarin (a vitamin K antagonist), enoxaparin (an indirect thrombin/Factor Xa inhibitor) and dabigatran (a direct thrombin inhibitor). A blood coagulation computer model has been developed, based on several published models and preclinical and clinical data. Unlike previous models, the current model takes into account both the intrinsic and extrinsic pathways of the coagulation cascade, and possesses some unique features, including a blood flow component and a portfolio of drug action mechanisms. This study aimed to use the model to compare the mechanism of action of rivaroxaban with that of warfarin, and to evaluate the efficacy and safety of different rivaroxaban doses with other anticoagulants included in the model. Rather than reproducing known standard clinical measurements, such as the prothrombin time and activated partial thromboplastin time clotting tests, the anticoagulant benchmarking was based on a simulation of physiologically plausible clotting scenarios. Compared with warfarin, rivaroxaban showed a favourable sensitivity for tissue factor concentration inducing clotting, and a steep concentration-effect relationship, rapidly flattening towards higher inhibitor concentrations, both suggesting a broad therapeutic window. The predicted dosing window is highly accordant with the final dose recommendation based upon extensive clinical studies.

A computational systems biology software platform for multiscale modeling and simulation: integrating whole-body physiology, disease biology, and molecular reaction networks.

Today, in silico studies and trial simulations already complement experimental approaches in pharmaceutical R&D and have become indispensable tools for decision making and communication with regulatory agencies. While biology is multiscale by nature, project work, and software tools usually focus on isolated aspects of drug action, such as pharmacokinetics at the organism scale or pharmacodynamic interaction on the molecular level. We present a modeling and simulation software platform consisting of PK-Sim(®) and MoBi(®) capable of building and simulating models that integrate across biological scales. A prototypical multiscale model for the progression of a pancreatic tumor and its response to pharmacotherapy is constructed and virtual patients are treated with a prodrug activated by hepatic metabolization. Tumor growth is driven by signal transduction leading to cell cycle transition and proliferation. Free tumor concentrations of the active metabolite inhibit Raf kinase in the signaling cascade and thereby cell cycle progression. In a virtual clinical study, the individual therapeutic outcome of the chemotherapeutic intervention is simulated for a large population with heterogeneous genomic background. Thereby, the platform allows efficient model building and integration of biological knowledge and prior data from all biological scales. Experimental in vitro model systems can be linked with observations in animal experiments and clinical trials. The interplay between patients, diseases, and drugs and topics with high clinical relevance such as the role of pharmacogenomics, drug-drug, or drug-metabolite interactions can be addressed using this mechanistic, insight driven multiscale modeling approach.