CardiacPBPK: A tool for the prediction and visualization of time-concentration profiles of drugs in heart tissue

BACKGROUND AND OBJECTIVE

Prediction of drug concentration in heart tissue is important in terms of drug safety and efficacy. This work presents the Open-Source CardiacPBPK platform for the prediction of the time-concentration profile of drugs, which could potentially reduce the risk of drug development failure due to cardiotoxicity. The objective of the CardiacPBPK development is to accelerate and simplify the in-silico toxicological assessment of new drugs, and to provide supportive material for the research community to use.

METHODS

The CardiacPBPK software provides a modular implementation of the PBPK model of heart tissue. It can be easily accessed via the Internet or installed locally. The graphical user interface and tabular design are easy to configure and use.

RESULTS

CardiacPBPK is a tool designed to predict and visualize the time-concentration profiles of a parent compound, and one metabolite, in venous plasma and heart tissue after oral or intravenous drug administration. CardiacPBPK is built on the R-environment framework and supports shiny application features such as interactive visualization of the results, and web applications interface by default. A shiny application refers to a computer program created with the use of shiny package in R. The application is freely available at https://github.com/jszlek/CardiacPBPK and https://sourceforge.net/projects/cardiacpbpk/. This open-source application runs on all platforms supporting R-environment (Linux, Windows, Mac OS X, Solaris).

CONCLUSIONS

We demonstrate the application of CardiacPBPK by simulating the study of amitriptyline intoxication in the case of CYP2D6 genetic polymorphism.

Better prediction of the local concentration-effect relationship: the role of physiologically based pharmacokinetics and quantitative systems pharmacology and toxicology in the evolution of model-informed drug discovery and development

Model-informed drug discovery and development (MID3) is an umbrella term under which sit several computational approaches: quantitative systems pharmacology (QSP), quantitative systems toxicology (QST) and physiologically based pharmacokinetics (PBPK). QSP models are built using mechanistic knowledge of the pharmacological pathway focusing on the putative mechanism of drug efficacy; whereas QST models focus on safety and toxicity issues and the molecular pathways and networks that drive these adverse effects. These can be mediated through exaggerated on-target or off-target pharmacology, immunogenicity or the physiochemical nature of the compound. PBPK models provide a mechanistic description of individual organs and tissues to allow the prediction of the intra- and extra-cellular concentration of the parent drug and metabolites under different conditions. Information on biophase concentration enables the prediction of a drug effect in different organs and assessment of the potential for drug-drug interactions. Together, these modelling approaches can inform the exposure-response relationship and hence support hypothesis generation and testing, compound selection, hazard identification and risk assessment through to clinical proof of concept (POC) and beyond to the market.

Physiologically based pharmacokinetic-quantitative systems toxicology and safety (PBPK-QSTS) modeling approach applied to predict the variability of amitriptyline pharmacokinetics and cardiac safety in populations and in individuals

The physiologically based pharmacokinetic (PBPK) models allow for predictive assessment of variability in population of interest. One of the future application of PBPK modeling is in the field of precision dosing and personalized medicine. The aim of the study was to develop PBPK model for amitriptyline given orally, predict the variability of cardiac concentrations of amitriptyline and its main metabolite-nortriptyline in populations as well as individuals, and simulate the influence of those xenobiotics in therapeutic and supratherapeutic concentrations on human electrophysiology. The cardiac effect with regard to QT and RR interval lengths was assessed. The Emax model to describe the relationship between amitriptyline concentration and heart rate (RR) length was proposed. The developed PBPK model was used to mimic 29 clinical trials and 19 cases of amitriptyline intoxication. Three clinical trials and 18 cases were simulated with the use of PBPK-QSTS approach, confirming lack of cardiotoxic effect of amitriptyline in therapeutic doses and the increase in heart rate along with potential for arrhythmia development in case of amitriptyline overdose. The results of our study support the validity and feasibility of the PBPK-QSTS modeling development for personalized medicine.

Mechanistic Physiologically Based Pharmacokinetic (PBPK) Model of the Heart Accounting for Inter-Individual Variability: Development and Performance Verification

Modern model-based approaches to cardiac safety and efficacy assessment require accurate drug concentration-effect relationship establishment. Thus, knowledge of the active concentration of drugs in heart tissue is desirable along with inter-subject variability influence estimation. To that end, we developed a mechanistic physiologically based pharmacokinetic model of the heart. The models were described with literature-derived parameters and written in R, v.3.4.0. Five parameters were estimated. The model was fitted to amitriptyline and nortriptyline concentrations after an intravenous infusion of amitriptyline. The cardiac model consisted of 5 compartments representing the pericardial fluid, heart extracellular water, and epicardial intracellular, midmyocardial intracellular, and endocardial intracellular fluids. Drug cardiac metabolism, passive diffusion, active efflux, and uptake were included in the model as mechanisms involved in the drug disposition within the heart. The model accounted for inter-individual variability. The estimates of optimized parameters were within physiological ranges. The model performance was verified by simulating 5 clinical studies of amitriptyline intravenous infusion, and the simulated pharmacokinetic profiles agreed with clinical data. The results support the model feasibility. The proposed structure can be tested with the goal of improving the patient-specific model-based cardiac safety assessment and offers a framework for predicting cardiac concentrations of various xenobiotics.

Humans Vary, So Cardiac Models Should Account for That Too!

The utilization of mathematical modeling and simulation in drug development encompasses multiple mathematical techniques and the location of a drug candidate in the development pipeline. Historically speaking they have been used to analyze experimental data (i.e., Hill equation) and clarify the involved physical and chemical processes (i.e., Fick laws and drug molecule diffusion). In recent years the advanced utilization of mathematical modeling has been an important part of the regulatory review process. Physiologically based pharmacokinetic (PBPK) models identify the need to conduct specific clinical studies, suggest specific study designs and propose appropriate labeling language. Their application allows the evaluation of the influence of intrinsic (e.g., age, gender, genetics, disease) and extrinsic [e.g., dosing schedule, drug-drug interactions (DDIs)] factors, alone or in combinations, on drug exposure and therefore provides accurate population assessment. A similar pathway has been taken for the assessment of drug safety with cardiac safety being one the most advanced examples. Mechanistic mathematical model-informed safety evaluation, with a focus on drug potential for causing arrhythmias, is now discussed as an element of the Comprehensive in vitro Proarrhythmia Assay. One of the pillars of this paradigm is the use of an in silico model of the adult human ventricular cardiomyocyte to integrate in vitro measured data. Existing examples (in vitro—in vivo extrapolation with the use of PBPK models) suggest that deterministic, epidemiological and clinical data based variability models can be merged with the mechanistic models describing human physiology. There are other methods available, based on the stochastic approach and on population of models generated by randomly assigning specific parameter values (ionic current conductance and kinetic) and further pruning. Both approaches are briefly characterized in this manuscript, in parallel with the drug-specific variability.

Thorough QT (TQT) studies: concordance with torsadogenesis and an evolving cardiac safety testing paradigm

Since 2005, when the International Conference on Harmonisation (ICH) E14 guideline was adopted, no drug has been withdrawn because of QTc prolongation or torsade de pointes arrhythmia. There are, however, costs associated with this success. In addition to the time and money invested, thorough QT (TQT) studies have limited the efficiency of the drug development pipeline. In this paper, we discuss the relevance of TQT trials as a tool for proarrhythmic risk prediction as a part of the debate regarding their usefulness.

Drug-drug interactions and QT prolongation as a commonly assessed cardiac effect - comprehensive overview of clinical trials

Background

Proarrhythmia assessment is one of the major concerns for regulatory bodies and pharmaceutical industry. ICH guidelines recommending preclinical tests have been established in attempt to eliminate the risk of drug-induced arrhythmias. However, in the clinic, arrhythmia occurrence is determined not only by the inherent property of a drug to block ion currents and disturb electrophysiological activity of cardiac myocytes, but also by many other factors modifying individual risk of QT prolongation and subsequent proarrhythmia propensity. One of those is drug-drug interactions. Since polypharmacy is a common practice in clinical settings, it can be anticipated that there is a relatively high risk that the patient will receive at least two drugs mutually modifying their proarrhythmic potential and resulting either in triggering the occurrence or mitigating the clinical symptoms. The mechanism can be observed either directly at the pharmacodynamic level by competing for the molecular targets, or indirectly by modifying the physiological parameters, or at the pharmacokinetic level by alteration of the active concentration of the victim drug.

Methods

This publication provides an overview of published clinical studies on pharmacokinetic and/or pharmacodynamic drug-drug interactions in humans and their electrophysiological consequences (QT interval modification). Databases of PubMed and Scopus were searched and combinations of the following keywords were used for Title, Abstract and Keywords fields: interaction, coadministration, combination, DDI and electrocardiographic, QTc interval, ECG. Only human studies were included. Over 4500 publications were retrieved and underwent preliminary assessment to identify papers accordant with the topic of this review. 76 papers reporting results for 96 drug combinations were found and analyzed.

Results

The results show the tremendous variability of drug-drug interaction effects, which makes one aware of complexity of the problem, and suggests the need for assessment of an additional risk factors and careful ECG monitoring before administration of drugs with anticipated QT prolongation.

Conclusions

DDIs can play significant roles in drugs’ cardiac safety, as evidenced by the provided examples. Assessment of the pharmacodynamic effects of the drug interactions is more challenging as compared to the pharmacokinetic due to the significant diversity in the endpoints which should be analyzed specifically for various clinical effects. Nevertheless, PD components of DDIs should be accounted for as PK changes alone do not allow to fully explain the electrophysiological effects in clinic situations.

Electronic supplementary material

The online version of this article (doi:10.1186/s40360-016-0053-1) contains supplementary material, which is available to authorized users.

The effect of increasing amitriptyline doses on cardiomyocytes’ electrophysiology – simulation study

Overdoses of tricyclic antidepressants may lead to arrhythmia. The aim of the study was to simulate the effect of increasing concentrations of amitriptyline (AMI) and its metabolite, nortriptyline, on the action potential of human ventricular cell.

Simulations were performed in Cardiac Safety Simulator platform with the use of the O’Hara-Rudy model. Input data included literature-derived, drug-specific IC50 values for I

The values of simulated endpoints (APD50, APD90, triangulation, and ΔAPD90) increase with drug concentrations. ΔAPD90 was statistically significant for doses up from 1000 mg. EADs were observed after administration of 10,000-mg AMI.

The consequences of various doses of AMI on the single cardiac myocytes were simulated in our study. Repolarization abnormalities were not expected for the therapeutic doses. EADs may be observed for very high doses of AMI.

Plasma vs heart tissue concentration in humans - literature data analysis of drugs distribution

Little is known about the uptake of drugs into the human heart, although it is of great importance nowadays, when science desires to predict tissue level behavior rather than to measure it. Although the drug concentration in cardiac tissue seems a better predictor for physiological and electrophysiological changes than its level in plasma, knowledge of this value is very limited. Tissue to plasma partition coefficients (Kp) come to rescue since they characterize the distribution of a drug among tissues as being one of the input parameters in physiologically based pharmacokinetic (PBPK) models. The article reviews cardiac surgery and forensic medical studies to provide a reference for drug concentrations in human cardiac tissue. Firstly, the focus is on whether a drug penetrates into heart tissue at a therapeutic level; the provided values refer to antibiotics, antifungals and anticancer drugs. Drugs that directly affect cardiomyocyte electrophysiology are another group of interest. Measured levels of amiodarone, digoxin, perhexiline and verapamil in different sites in human cardiac tissue where the compounds might meet ion channels, gives an insight into how these more lipophilic drugs penetrate the heart. Much data are derived from postmortem studies and they provide insight to the cardiac distribution of more than 200 drugs. The analysis depicts potential problems in defining the active concentration location, what may indirectly suggest multiple mechanisms involved in the drug distribution within the heart. Copyright © 2015 John Wiley & Sons, Ltd.

Abnormal olanzapine toxicokinetic profiles--population pharmacokinetic analysis

CONTEXT

Olanzapine is widely used in the treatment of schizophrenia and it is becoming more frequently responsible for overdoses. Standard pharmacokinetic models do not fit to the toxic concentration data.

OBJECTIVE

The aim of present study is to investigate the reasons for an abnormal olanzapine plasma concentration time curve in the range of toxic concentrations. Two hypotheses were verified: entering the enterohepatic cycle, and drug deposition and its desorption from activated charcoal used for gastrointestinal decontamination.

MATERIALS AND METHODS

One-hundred thirty-five plasma concentration data from 21 patients hospitalized for acute olanzapine poisoning were analyzed with the use of the population pharmacokinetic approach. A non-linear mixed-effects modeling approach with Monolix 4.3.1 was employed.

RESULTS

A model assuming gallbladder emptying at irregular intervals was developed. Also, a model that describes desorption of olanzapine from the charcoal surface, in which the dose is divided into two absorbed fractions, was constructed. The analysis has found gastrointestinal decontamination and previous olanzapine treatment, as the significant covariates for toxicokinetic parameters of olanzapine.

CONCLUSION

Our study provides interesting models for investigation of toxic concentration of olanzapine, which may also be used as the basis for further model development for other drugs as well. The investigated population was not large enough to reliably confirm any of the proposed models. It would be well worth continuing this study with more substantial data. Also, any additional information about olanzapine metabolite concentration could be vital.

Model of the distribution of diastolic left ventricular posterior wall thickness in healthy adults and its impact on the behavior of a string of virtual cardiomyocytes

Correlation of the thickness of the left ventricular posterior wall (LVPWd) with various parameters, including age, gender, weight and height, was investigated in this study using regression models. Multicenter derived database comprised over 4,000 healthy individuals. The developed models were further utilized in the in vitro-in vivo (IVIV) translation of the drug cardiac safety data with use of the mathematical model of human cardiomyocytes operating at the virtual healthy population level. LVPWd was assumed to be equivalent to the length of one-dimensional string of virtual cardiomyocyte cells which was presented, as other physiological factors, to be a parameter influencing the simulated pseudo-ECG (pseudoelectrocardiogram), QTcF and ∆QTcF, both native and modified by exemplar drug (disopyramide) after I Kr current disruption. Simulation results support positive correlation between the LVPWd and QTcF/∆QTc. Developed models allow more detailed description of the virtual population and thus inter-individual variability influence on the drug cardiac safety.