The effects of six antipsychotic agents on QTc--an attempt to mimic clinical trial through simulation including variability in the population

Investigating the relevance of CYP2J2 inhibition for drugs known to cause intermediate to high risk torsades de pointes

Cardiac cytochrome P450 2J2 (CYP2J2) metabolizes endogenous polyunsaturated fatty acid, arachidonic acid (AA), to bioactive regioisomeric epoxyeicosatrienoic acid (EET) metabolites. This endogenous metabolic pathway has been postulated to play a homeostatic role in cardiac electrophysiology. However, it is unknown if drugs that cause intermediate to high risk torsades de pointes (TdP) exhibit inhibitory effects against CYP2J2 metabolism of AA to EETs. In this study, we demonstrated that 11 out of 16 drugs screened with intermediate to high risk of TdP as defined by the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative are concurrently reversible inhibitors of CYP2J2 metabolism of AA, with unbound inhibitory constant (K_i,AA,u) values ranging widely from 0.132 to 19.9 µM. To understand the physiological relevancy of K_i,AA,u, the in vivo unbound drug concentration within human heart tissue (C_u,heart) was calculated via experimental determination of in vitro unbound partition coefficient (K_puu) for 10 CYP2J2 inhibitors using AC16 human ventricular cardiomyocytes as well as literature-derived values of fraction unbound in plasma (f_u,p) and plasma drug concentrations in clinical scenarios leading to TdP. Notably, all CYP2J2 inhibitors screened belonging to the high TdP risk category, namely vandetanib and bepridil, exhibited highest K_puu values of 18.2 ± 1.39 and 7.48 ± 1.16 respectively although no clear relationship between C_u,heart and risk of TdP could eventually be determined. R values based on basic models of reversible inhibition as per FDA guidelines were calculated using unbound plasma drug concentrations (C_u,plasma) and adapted using C_u,heart which suggested that 4 out of 10 CYP2J2 inhibitors with intermediate to high risk of TdP demonstrate greatest potential for clinically relevant in vivo cardiac drug-AA interactions. Our results shed novel insights on the relevance of CYP2J2 inhibition in drugs with risk of TdP. Further studies ascertaining the role of CYP2J2 metabolism of AA in cardiac electrophysiology, characterizing inherent cardiac ion channel activities of drugs with risk of TdP as well as in vivo evidence of drug-AA interactions will be required prior to determining if CYP2J2 inhibition could be an alternative mechanism contributing to drug-induced TdP.

Exploring the Potential of Generative Adversarial Networks for Synthesizing Radiological Images of the Spine to be Used in In Silico Trials

In silico trials recently emerged as a disruptive technology, which may reduce the costs related to the development and marketing approval of novel medical technologies, as well as shortening their time-to-market. In these trials, virtual patients are recruited from a large database and their response to the therapy, such as the implantation of a medical device, is simulated by means of numerical models. In this work, we propose the use of generative adversarial networks to produce synthetic radiological images to be used in in silico trials. The generative models produced credible synthetic sagittal X-rays of the lumbar spine based on a simple sketch, and were able to generate sagittal radiological images of the trunk using coronal projections as inputs, and vice versa. Although numerous inaccuracies in the anatomical details may still allow distinguishing synthetic and real images in the majority of cases, the present work showed that generative models are a feasible solution for creating synthetic imaging data to be used in in silico trials of novel medical devices.

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.

Quantitative systems toxicology

The overarching goal of modern drug development is to optimize therapeutic benefits while minimizing adverse effects. However, inadequate efficacy and safety concerns remain to be the major causes of drug attrition in clinical development. For the past 80 years, toxicity testing has consisted of evaluating the adverse effects of drugs in animals to predict human health risks. The U.S. Environmental Protection Agency recognized the need to develop innovative toxicity testing strategies and asked the National Research Council to develop a long-range vision and strategy for toxicity testing in the 21st century. The vision aims to reduce the use of animals and drug development costs through the integration of computational modeling and in vitro experimental methods that evaluates the perturbation of toxicity-related pathways. Towards this vision, collaborative quantitative systems pharmacology and toxicology modeling endeavors (QSP/QST) have been initiated amongst numerous organizations worldwide. In this article, we discuss how quantitative structure-activity relationship (QSAR), network-based, and pharmacokinetic/pharmacodynamic modeling approaches can be integrated into the framework of QST models. Additionally, we review the application of QST models to predict cardiotoxicity and hepatotoxicity of drugs throughout their development. Cell and organ specific QST models are likely to become an essential component of modern toxicity testing, and provides a solid foundation towards determining individualized therapeutic windows to improve patient safety.

A four-compartment PBPK heart model accounting for cardiac metabolism - model development and application

In the field of cardiac drug efficacy and safety assessment, information on drug concentration in heart tissue is desirable. Because measuring drug concentrations in human cardiac tissue is challenging in healthy volunteers, mathematical models are used to cope with such limitations. With a goal of predicting drug concentration in cardiac tissue, we have developed a whole-body PBPK model consisting of seventeen perfusion-limited compartments. The proposed PBPK heart model consisted of four compartments: the epicardium, midmyocardium, endocardium, and pericardial fluid, and accounted for cardiac metabolism using CYP450. The model was written in R. The plasma:tissues partition coefficients (Kp) were calculated in Simcyp Simulator. The model was fitted to the concentrations of amitriptyline in plasma and the heart. The estimated parameters were as follows: 0.80 for the absorption rate [h-1], 52.6 for Kprest, 0.01 for the blood flow through the pericardial fluid [L/h], and 0.78 for the P-parameter describing the diffusion between the pericardial fluid and epicardium [L/h]. The total cardiac clearance of amitriptyline was calculated as 0.316 L/h. Although the model needs further improvement, the results support its feasibility, and it is a first attempt to provide an active drug concentration in various locations within heart tissue using a PBPK approach.

Early assessment of proarrhythmic risk of drugs using the in vitro data and single-cell-based in silico models: proof of concept

BACKGROUND AND PURPOSE

To determine the predictive performance of in silico models using drug-specific preclinical cardiac electrophysiology data to investigate drug-induced arrhythmia risk (e.g. Torsade de pointes (TdP)) in virtual human subjects.

EXPERIMENTAL APPROACH

To assess drug proarrhythmic risk, we used a set of in vitro electrophysiological measurements describing ion channel inhibition triggered by the investigated drugs. The Cardiac Safety Simulator version 2.0 (CSS; Simcyp, Sheffield, UK) platform was used to simulate human left ventricular cardiac myocyte action potential models.

RESULTS

This study shows the impact of drug concentration changes on particular ionic currents by using available experimental data. The simulation results display safety threshold according to drug concentration threshold and log (threshold concentration/ effective therapeutic plasma concentration (ETPC)).

CONCLUSION AND IMPLICATIONS

We reproduced the underlying biophysical characteristics of cardiac cells resulted in effects of drugs associated with cardiac arrhythmias (action potential duration (APD) and QT prolongation and TdP) which were observed in published 3D simulations, yet with much less computational burden.

Virtual Clinical Trial Toward Polytherapy Safety Assessment: Combination of Physiologically Based Pharmacokinetic/Pharmacodynamic-Based Modeling and Simulation Approach With Drug-Drug Interactions Involving Terfenadine as an Example

A Quantitative Systems Pharmacology approach was utilized to predict the cardiac consequences of drug-drug interaction (DDI) at the population level. The Simcyp in vitro-in vivo correlation and physiologically based pharmacokinetic platform was used to predict the pharmacokinetic profile of terfenadine following co-administration of the drug. Electrophysiological effects were simulated using the Cardiac Safety Simulator. The modulation of ion channel activity was dependent on the inhibitory potential of drugs on the main cardiac ion channels and a simulated free heart tissue concentration. ten Tusscher's human ventricular cardiomyocyte model was used to simulate the pseudo-ECG traces and further predict the pharmacodynamic consequences of DDI. Consistent with clinical observations, predicted plasma concentration profiles of terfenadine show considerable intra-subject variability with recorded C_max values below 5 ng/mL for most virtual subjects. The pharmacokinetic and pharmacodynamic effects of inhibitors were predicted with reasonable accuracy. In all cases, a combination of the physiologically based pharmacokinetic and physiology-based pharmacodynamic models was able to differentiate between the terfenadine alone and terfenadine + inhibitor scenario. The range of QT prolongation was comparable in the clinical and virtual studies. The results indicate that mechanistic in vitro-in vivo correlation can be applied to predict the clinical effects of DDI even without comprehensive knowledge on all mechanisms contributing to the interaction.

Systems pharmacology to predict drug safety in drug development

Ensuring that drugs are safe and effective is a very high priority for drug development and the US Food and Drug Administration review process. This is especially true today because of faster approval times and smaller clinical trials, especially in oncology and rare diseases. In light of these trends, systems pharmacology is seen as an essential strategy to understand and predict adverse drug events during drug development by analyzing interactions between drugs and multiple targets rather than the traditional "one-drug-one-target" approach. This commentary offers an overview of the current trends and challenges of using systems pharmacology to reduce the risks of unintended adverse events.

Top-down, Bottom-up and Middle-out Strategies for Drug Cardiac Safety Assessment via Modeling and Simulations

Cardiac safety is an issue causing early terminations at various stages of drug development. Efforts are put into the elimination of false negatives as well as false positives resulting from the current testing paradigm. In silico approaches offer mathematical system and data description from the ion current, through cardiomyocytes level, up to incorporation of inter-individual variability at the population level. The article aims to review three main modelling and simulation approaches, i.e. "top-down" which refers to models built on the observed data, "bottom-up", which stands for a mechanistic description of human physiology, and "middle-out" which combines both strategies. Modelling and simulation is a well-established tool in the assessment of drug proarrhythmic potency with an impact on research and development as well as on regulatory decisions, and it is certainly here to stay. What is more, the shift to systems biology and physiology-based models makes the cardiac effect more predictable.

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.

Computational investigations of hERG channel blockers: New insights and current predictive models

Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.

Modeling and Simulation Approaches for Cardiovascular Function and Their Role in Safety Assessment

Systems pharmacology modeling and pharmacokinetic-pharmacodynamic (PK/PD) analysis of drug-induced effects on cardiovascular (CV) function plays a crucial role in understanding the safety risk of new drugs. The aim of this review is to outline the current modeling and simulation (M&S) approaches to describe and translate drug-induced CV effects, with an emphasis on how this impacts drug safety assessment. Current limitations are highlighted and recommendations are made for future effort in this vital area of drug research.

Interaction Between Domperidone and Ketoconazole: Toward Prediction of Consequent QTc Prolongation Using Purely In Vitro Information

We aimed to investigate the application of combined mechanistic pharmacokinetic (PK) and pharmacodynamic (PD) modeling and simulation in predicting the domperidone (DOM) triggered pseudo-electrocardiogram modification in the presence of a CYP3A inhibitor, ketoconazole (KETO), using in vitro-in vivo extrapolation. In vitro metabolic and inhibitory data were incorporated into physiologically based pharmacokinetic (PBPK) models within Simcyp to simulate time course of plasma DOM and KETO concentrations when administered alone or in combination with KETO (DOM+KETO). Simulated DOM concentrations in plasma were used to predict changes in gender-specific QTcF (Fridericia correction) intervals within the Cardiac Safety Simulator platform taking into consideration DOM, KETO, and DOM+KETO triggered inhibition of multiple ionic currents in population. Combination of in vitro-in vivo extrapolation, PBPK, and systems pharmacology of electric currents in the heart was able to predict the direction and magnitude of PK and PD changes under coadministration of the two drugs although some disparities were detected.

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.