# Humans vary, so cardiac models should account for that too!

## Introduction

### Mechanistic Modeling and Simulation Approach in the Process of Drug Development in Light of the Recent Changes to FDA/EMA/PMDA Regulations

The mathematical modeling and simulation (M&S) approach has held its place in the drug development process since the very beginning. Having moved from academic curiosity to industrial practice the approach is used to both analyze the data and understand the physical mechanisms involved. Model-informed drug development (MIDD) is here to stay and it also has been recently indicated as one of the major areas of scientific priority by the regulators ( Huang et al., 2013 ; Zineh et al., 2017 ). Initially, the outcome of pharmacometric analyses has impacted the decision making process of drug development ( Miller et al., 2005 ; Lee et al., 2011 ). In recent years regulatory support for mechanistic physiological modeling has helped to bridge the gap between preclinical and clinical observations with respect to understanding biological systems ( Rowland et al., 2015 ; Friedrich, 2016 ).

Special interest has been given to physiologically based pharmacokinetic (PBPK) modeling defined by WHO as “ quantitative descriptions of the absorption, distribution, metabolism, and excretion (ADME) of chemicals in biota based on interrelationships among key physiological, biochemical and physicochemical determinants of these processes” 1 . A physiologically based pharmacokinetic approach is based on a combination of the physiology, environment, and drug specific information. These parameters are further utilized in the mechanistic models describing the pharmacokinetics (PK) and/or pharmacodynamics (PD) of a drug(s) of interest ( Rostami-Hodjegan, 2012 ). This is not a new concept and it is suggested that the roots of PBPK models originate from the work of Teorell, published in 1937 ( Teorell, 1937 ). Recent scientific advances and the development of models utilizing PBPK scaffolding are invaluable in the situation where clinical trials are extremely challenging or impossible. This includes specialized populations and situations of special interest such as pregnant woman and pediatric applications ( Lu et al., 2012 ; Abduljalil et al., 2014 ). PBPK models are utilized for various applications throughout a drug’s life cycle. Results of their simulations can be used to support the planning of specialized clinical studies, support dosing recommendations and the labeling of products ( Zhao et al., 2011 ; Jones et al., 2015 ; Shepard et al., 2015 ; Wagner et al., 2015 ). The simulation results are used in lieu of conducting clinical studies or provide information that otherwise would have been missing in some specific situations ( Jamei, 2016 ).

The use of PBPK modeling was included in the guidance documents for industry provided by the U. S. Food and Drug Administration (FDA), the European Medicines Agency (EMA) and the Ministry of Health Labor and Welfare (MHLW) of Japan. The latter, namely the draft of the drug interaction guideline for drug development and labeling recommendations published by the Japanese Ministry of Health, Labor and Welfare in 2014, suggested the PBPK application in the assessment of drug-drug interaction (DDI; Saito et al., 2014 ). Recent guidelines published by the FDA and EMA list several points in the drug development process where PBPK modeling may be applicable in order to support decisions in the premarketing, as well as at postmarketing, stage 2 , 3 . PBPK analyses are currently widely accepted and used not only as a research tool but also to support drug registration applications including investigational new drugs (INDs), new drug applications (NDAs), biologics license applications (BLAs), or abbreviated new drug applications (ANDAs; Sato et al., 2017 ). As of December 2016 there were 36 approved drugs, which provided PBPK M&S results, in the U. S. new drug application labeling procedure. In addition to PBPK, the Advisory Committee mentioned mechanistic safety modeling, particularly, risk prediction/assessment as a promising MIDD area 4 .

It is highly likely that the next application for PBPK models is in the area of precision dosing and personalized medicine ( Hartmanshenn et al., 2016 ). This is achievable due to the specific structure of the physiologically-based models, where system description is clearly separated from the drug and external parameters (i. e., dosing schema). The biological parameters are described by large collections of anatomical and physiological data derived from literature or existing databases. For the assessment of inter- and to some degree intra-individual variability, virtual individuals, and virtual populations are randomly created ( Jamei et al., 2009a ).

### In Vitro—in Vivo Extrapolation (IVIVE) as an Approach

In vitro—in vivo extrapolation (IVIVE) as a phrase covers all techniques utilized for the prediction of human pharmacokinetics and pharmacodynamics based on the ADME information ADME in addition to drug activity and toxicity. It is crucial to mention that, in principle, the data comes from the in vitro studies, where various models and techniques are utilized. The main challenge is to translate the data, often heterogeneous, from the level of a “ Petri dish” to the complex system of the human body. Therefore, what is needed are in vitro methods mimicking basic phenomena occurring in a human body at the cellular or subcellular level, and models describing the human body as a biological system. The latter describing models of human organs, tissues, or the whole body include the incorporation of the above mentioned PBPK systems. There are other approaches currently implemented, which combine and merge traditional compartmental models and various systems biology models, to describe biochemical and physiological phenomena ( Sorger et al., 2011 ). It’s worth adding that the addition of available human in vivo data enriches the model and can improve the degree of predictability ( Tsamandouras et al., 2015 ). Practical utilization of the IVIVE concept, with the use of PBPK models and population data, covers various populations from healthy individuals up to the special populations (e. g., diseased—renal insufficiency, cirrhosis etc.) and pediatric populations ( Sager et al., 2015 ).

Interestingly, conceptually similar approaches were independently introduced in the pharmacokinetics and pharmacodynamics area ( Visser et al., 2014 ). The biophysically detailed models describing cell electrophysiology can be used for the in vitro—in vivo extrapolation, as proposed in the CiPA (Comprehensive in vitro Proarrhythmia Assay) initiative, as described below ( O’Hara et al., 2011 ; Gintant et al., 2016 ). There is although one difference between the IVIVE/PBPK approach and the proposed IVIVE/safety approach. The latter, namely safety in vitro—in vivo extrapolation, has so far been done with the use of models describing cell electrophysiology without accounting for variability.

### CiPA Initiative—What Does It Change, Where We Are, and What Is Potentially Lacking

There are also special populations that should be considered when assessing drug triggered cardiac effects. First of all, the cardiac action potential is affected in patients with cardiac channelopathies associated with genetic mutations. One of the main congenital phenotypes is long QT syndrome (LQTS) which is prevalent in 1: 3, 000–1: 5, 000 in the general population ( Goldenberg and Moss, 2008 ) manifesting with a QTc length of above 460 ms ( Abriel and Zaklyazminskaya, 2013 ). The most frequent mutations that are responsible for congenital LQTS were found in genes that code for proteins in the potassium hERG channel (KCNQ1 and KCNH2) and the Nav1. 5 sodium channel (SCN5A; Bohnen et al., 2016 ). The LQTS patients are said to be particularly vulnerable to drug-related arrhythmias ( Goldenberg et al., 2008 ). Also, comorbidities were found to contribute to QTc-prolongation. The analysis conducted by Vandael et al. (2017) revealed strong evidence for ischemic cardiomyopathy, hypertension, arrhythmia, and thyroid disturbances to be risk factors contributing to QTc interval prolongation.

### Variability in Cardiac Models—Stochastic vs. Deterministic Approach

Mathematical models of cardiac cell electrophysiology have proved their value and now hold an established position in research and drug development ( Amanfu and Saucerman, 2011 ; Davies et al., 2016 ). It all began with the work of Hodgkin and Huxley where they modeled the cell membrane as a capacitor with batteries and resistors and described its electrophysiological behavior ( Hodgkin and Huxley, 1952 ). The following ordinary differential equation was established:

$d V d t = I i o n − I s t i m C m ,$

where V is voltage, t is time, I ion is the sum of all transmembrane ionic currents, I stim is the externally applied stimulus current, and C m is cell capacitance of the membrane per unit surface area ( Ten Tusscher and Panfilov, 2006 ). Beginning from the relatively simple Noble’s model ( Noble, 1962 ), cardiac electrophysiology models have evolved tremendously, and now include a detailed description on cardiac ion channels, pumps and transporters, as well as intracellular calcium handling ( Noble, 2007 ; Fink et al., 2011 ). The models proved successful for studying cardiac physiology and understanding pathological changes (e. g., arrhythmias) associated with diseases over the last decades and currently they are used in the safety assessment of drugs ( Mirams et al., 2012 ; Roberts et al., 2012 ; Davies et al., 2016 ; Gintant et al., 2016 ). To be a vital element of the safety-related decisions made through the various stages of drug development and in the clinic, the models have to ascertain credibility, i. e., properly reproduce the electrophysiology of a population and individual patients. The question arises asking if it is feasible to use a single traditional cardiac model with input parameters being the mean values, averaged across many subjects, and generating an output as a single value, thus presenting behavior a of “ representative” cell. Until recently there wasn’t much interest in the variability in the field of cardiac electrophysiology modeling. Several approaches have now been proposed, and implemented, to introduce variability into cardiac models and account for inter- and intra-patients differences. All of the approaches stem from the belief that the “ average patient” does not exist and a traditional model cannot accurately explain the observed differences between patients.

### Stochastic Approach

Physiological variability can be investigated and modeled by constructing populations of experimentally calibrated models. This approach has been introduced a while ago and is further developed by researchers from various organizations including, but not limited to, the Department of Pharmacology and Systems Therapeutics, Mount Sinai School of Medicine, New York ( Sobie, 2009 ; Sarkar et al., 2012 ), the Department of Computer Sciences at Oxford University ( Britton et al., 2013 ; Sánchez et al., 2014 ; Muszkiewicz et al., 2016 ; Pueyo et al., 2016 ; Zhou et al., 2016 ), and other ( Marder and Taylor, 2011 ). Variability in a model’s behavior is accomplished by the multiplication of its parameter values by sampled scaling factors, which results in an ensemble of possible outputs instead of an average one. Once the baseline model (appropriate for the research question) is selected, the parameters of variation have to be chosen, and the ranges for parameter sampling need to be defined. Both depend on the study aim and the target population (e. g., healthy, diseased) to be investigated.

There are multiple approaches available for generating the population of models. For example, parameters and their ranges generate high-dimensional spaces from which their values are sampled using different sampling algorithms (e. g., the Latin Hypercube, Monte Carlo method) to generate an ensemble of variant models ( Drovandi et al., 2016 ). The candidate models are simulated to mimic, for example, experimental settings regarding bath solution composition, voltage protocol, and pacing rate. The pool of models is pruned according to the boundaries for a range of permitted model outputs, defined by minimal and maximal values observed in the experiment. Such calibration yields the experimentally-calibrated populations of models presenting electrophysiological behavior that reproduce the results of experiments which can be further analyzed for the parameters underlying the observed response variability ( Britton et al., 2013 ). The model could also be employed to assess a range of responses across the population under certain conditions, e. g., diseased, or reflecting drug application ( Muszkiewicz et al., 2016 ).

The other approach is based on parameter sensitivity analysis techniques which can be applied to generate quantitative predictions based on considering behaviors within a population of models ( Romero et al., 2009 ). In the single parameter scanning one can increase or decrease the parameter of interest, run simulation and save the simulation results. Such procedure can be performed for multiple parameters, with the use of various techniques and the outcome is expressed in a quantitative manner. Sobie and colleagues have also utilized more complex procedures and varied all parameters at once ( Sarkar and Sobie, 2010 , 2011 ). Assuming that the endpoint of interest is dichotomous in nature (i. e., EAD occurence) Morotti and Grandi proposed technique based on the multivariate logistic regression analysis, allowing for sensitivity analysis and investigating factors influencing the endpoint occurrence at the same time ( Morotti and Grandi, 2016 ).

The methodology applied to experimentally-calibrated populations of models can be regarded as the refinement and extension of the sensitivity analysis method and populations of models constructed without a calibration step. In the latter two, values of single or multiple parameters are varied in a predefined range, while allowable values of model outputs are not constrained with experimental data thus may not be representative for certain subject groups. All of these models offer valuable insights into sources of variability and the pathological background of different heart conditions. Moreover, they allow the assessment of not only the average drug effect on cardiac electrophysiology, but also allow the screening of drug effects across a certain population. They, however, cannot be employed to evaluate the risk of an individual patient who is to be treated with a certain drug. This is because most variability sources considered by the model (e. g., ion channels conductance, gating kinetics, densities, resting membrane potential, membrane potential, upstroke velocity) are factors whose values cannot be determined for the particular patient. Another problem is the requirement to define the cut off threshold for non-physiological simulation results, since such a decision is always arbitrary.

### Deterministic Approach

The other approach is based on the observation and description of the biological parameters characterizing the system of interest, in this case the human organism. Data describing demographic, genetic, anatomical, and physiological factors are collected and analyzed to describe the distribution of the parameter in a population and allow constructing its virtual counterpart instead of an average patient. Also some of those factors can be correlated with the parameters of a cardiomyocyte model, e. g., age and the cardiomyocyte volume and area ( Polak et al., 2012 ). The main data sources remain the published scientific reports and clinical databases ( ICRP, 2002 ; Valentin, 2003 ). This approach allows for differentiation between different populations; both healthy (e. g., considering ethnicity) and diseased (e. g., obese, renally impaired, diabetic), or of special interest (e. g., pregnant women, pediatric population; Jamei et al., 2009b ). Equations describing the distribution of the system parameters for the PBPK model are derived from the distributions of data based on real populations and patients. Additionally, such an approach takes into account the dynamic changes observed with the parameter of choice (e. g., circadian changes of heart rate). An application includes the ontogeny description and age dependent variation of the chosen parameters ( Salem et al., 2013 ). What is also of importance is the option to define the inter-correlation between various parameters. A virtual population can be generated from values and formulae describing demographic, anatomical, and physiological variables using a correlated Monte Carlo approach which protects from the non-physiological combinations of the model parameters (e. g., kidney size and liver size). This allows the prediction of variability before the clinical study phase, in contrast to a statistical approach (e. g., population PK analysis), which requires prior clinical data to characterize variability. Additionally it allows for the clear separation of information on the system (i. e., human body) from that of the drug (e. g., physicochemical characteristics), and the environment (e. g., dose, concomitant drugs).

It is obvious that the model prediction will depend on the quality of the gathered data, correctness of the analysis, and appropriateness of the conclusions drawn from the distribution of the parameters. Therefore, the criteria of inclusion and exclusion for the reports and papers (in addition to certain values) have to be defined prior to the commencement of the data collection process. The same applies to the statistical analysis, and all methods and tools.

There are multiple biological parameters influencing the ECG and its behavior. From a biological perspective, factors influencing cardiac electrophysiology can be classified into one of three key groups: (i) demography; (ii) anatomy, and/or physiology; and (iii) genetics. Examples of the most important parameters are listed below:

(i) Demography—age, gender

(ii) Anatomy and/or physiology—plasma ions concentrations (K +, Ca 2+), cardiomyocyte size (volume, area), cell electric capacitance, heart wall thickness, heterogeneity of cells across heart wall, heart rate, sex hormones

(iii) Genetics—common polymorphisms and mutations at the level of ion channels/pumps/exchangers

For some of the above mentioned parameters statistical models describing their distribution in the population were developed based on the available data. This includes relationship between age, human left ventricle cardiomyocyte volume, and electric capacitance ( Polak et al., 2012 ), left ventricular heart wall thickness ( Fijorek et al., 2014a ), or plasma ions concentration ( Fijorek et al., 2014b ). The latter, namely plasma ions concentration, together with the heart rate follow the circadian variability ( Massin et al., 2000 ; Sennels et al., 2012 ). Such time of the day dependent variation influences the ECG and the parameters including QT ( Karjalainen et al., 1994 ). Therefore, models describing the diurnal variability of such crucial parameters are also being developed and utilized ( Fijorek et al., 2013a , b ). They are also included in the virtual population generators, which build the virtual individuals being exposed to the drug in the in silico conditions ( Mishra et al., 2014 ; Wiśniowska and Polak, 2016 ). The main limitation of this method is the limited amount of data available for the model building and information relating to the inter-correlation between parameters.

It is undisputable that sources of variability identified and already incorporated into in silico models cannot explain the whole observable inter- and intra-individual variability. The question is if it is already sufficient to provide reliable predictions.

### Drugs-Specific Variability

In the model-based drug safety assessment, apart from system-dependent variability, there is also another source of variability (or rather uncertainty) introduced into the model. The effects caused by exposure to a drug are modeled as concentration-dependent changes of ion currents known to be affected by the drug. The drug-channel interaction is quantified experimentally by constructing a Hill plot characterized by IC50 and the Hill coefficient value. There is no standardized protocol for drug-triggered ion channel block measurement, however multiple cell models and voltage protocols are accepted. This results in significant discrepancies between IC50 values reported for a given compound (amiodarone: 0. 015–38. 3 microM; cisapride: 0. 0051–1. 6 microM; dofetilide 0. 003–25 microM; E4031: 0. 001–15. 8 microM; moxifloxacin: 0. 93–398. 1 microM; propafenone: 0. 085–123 microM; to give just a few examples) 6 . Additional, uncertainty introduced into in silico models comes from uncertainty in the Hill coefficients of reported concentration-inhibition curves, especially when ion channel blocking potency for a compound is estimated by high-throughput screening methods ( Elkins et al., 2013 ). Moreover, even the well-controlled experiments, carefully conducted in the same experimental settings, generate different results. This can be due to the measurement errors, intrinsic- and extrinsic-variability between samples. The subjective selection of IC50 values, which are taken as inputs into in silico models, may lead to misleading results. Drug-specific uncertainty is an undesirable type of variability, however infeasible to eliminate. The use of averaged values has been proposed to minimize drug-specific bias ( Elkins et al., 2013 ). The question remains how, if possible, average results produced with different methods and experimental settings.

The situation is even more complex when the aim is to model a response to a drug in a realistic population where other drug-related factors like concomitant medicines, compliance, food and alcohol consumption have to be accounted for since all these factors can substantially influence the response.

## Conclusions

Drug pro-arrhythmic potency is a function of the intrinsic characteristics of the chemical structure and external parameters associated with the drug. The latter, namely external parameters, includes system-dependent and environment-dependent items. Therefore, to properly predict and assess the potential risk all significant elements should be considered. Translation of the in vitro data to an in vivo situation, e. g., to optimize clinical trials, requires such an approach. There are at least two ways to describe the population specific variability, as described above, and these elements should be accounted for in the in silico based cardiac safety assessment. On the other hand, the choice of an assessment method depends on the actual stage of the drug development cycle. In the early stages, including discovery phase such complexity is not necessary, but as the compound advances along the development process, the assessment should be more detailed and comprehensive.

It is worth noting that the diseased population would probably require even more parameters, and the variability would be larger, as compared against healthy individuals. This is because of the interrelation of parameters directly (e. g., cell volume and capacitance), and indirectly (e. g., thyroid-related diseases and fluid-electrolyte balance disruption) influence the heart cells electrophysiology. Examples of the above-mentioned parameters and their analysis clude hypertrophic cardiomyopathy (HCM) of various character ( Polak and Fijorek, 2012 ). It is worth noting that stochastic approach, based on the virtually simulated population of models can be used to identify ionic mechanisms driving electrophysiological abnormalities not only in HCM, but also atrial fibrillation, and other diseases, accounting for of the disease specific variability ( Liberos et al., 2016 ; Passini et al., 2016 ). Modeling and simulation approach can be also utilized to analyze the influence of genetic modifications at the level of ionic channels on the cardiac myocytes electrophysiology ( Glinka and Polak, 2013 ). This includes the recently published example of the population of models optimized to recapitulate clinical long QT phenotypes ( Mann et al., 2016 ). The detailed discussion of this element, namely the disease related variability is out of the current manuscript scope. There is however clear gap in all current approaches, namely comprehensive, or to be precise—as comprehensive as possible—parametrization of the disease of choice, including all known parameters. This would allowed to simulate population of virtual individuals closely mimicking those met in the clinical settings. The models usage could be than extended from the drugs’ safety assessment to the drugs therapy optimization.

## Author Contributions

BW, ZT, and SP equally contributed to the conception of paper, drafting the manuscript and approved its final version.

## Conflict of Interest Statement

SP is a Simcyp (part of Certara) employee.

The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

Authors would like to thank Dr. Ruth Clayton for professional proof-reading and editorial support.

## References

Abduljalil, K., Jamei, M., Rostami-Hodjegan, A., and Johnson, T. N. (2014). Changes in individual drug-independent system parameters during virtual paediatric pharmacokinetic trials: introducing time-varying physiology into a paediatric PBPK model. AAPS J. 16, 568–576. doi: 10. 1208/s12248-014-9592-9

Abriel, H., and Zaklyazminskaya, E. V. (2013). Cardiac channelopathies: genetic and molecular mechanisms. Gene 517, 1–11. doi: 10. 1016/j. gene. 2012. 12. 061

Amanfu, R. K., and Saucerman, J. J. (2011). Cardiac models in drug discovery and development: a review. Crit. Rev. Biomed. Eng. 39, 379–395. doi: 10. 1615/CritRevBiomedEng. v39. i5. 30

Blinova, K., Stohlman, J., Vicente, J., Chan, D., Johannesen, L., Hortigon-Vinagre, M. P., et al. (2017). Comprehensive translational assessment of human-induced pluripotent stem cell derived cardiomyocytes for evaluating drug-induced arrhythmias. Toxicol. Sci. 155, 234–247. doi: 10. 1093/toxsci/kfw200

Bohnen, M. S., Peng, G., Robey, S. H., Terrenoire, C., Iyer, V., Sampson, K. J., et al. (2016). Molecular pathophysiology of congenital long QT syndrome. Physiol. Rev. 97, 89–134. doi: 10. 1152/physrev. 00008. 2016

Britton, O. J., Rodriguez, B., Bueno-Orovio, A., Van, A. K., Lu, H. R., Towart, R., et al. (2013). Experimentally calibrated population of models predicts and explains intersubject variability in cardiac cellular electrophysiology. Proc. Natl. Acad. Sci. U. S. A 110, E2098–E2105. doi: 10. 1073/pnas. 1304382110

Colatsky, T., Fermini, B., Gintant, G., Pierson, J. B., Sager, P., Sekino, Y., et al. (2016). The comprehensive in vitro proarrhythmia assay (CiPA) initiative – update on progress. J. Pharmacol. Toxicol. Methods 81, 15–20. doi: 10. 1016/j. vascn. 2016. 06. 002

Davies, M. R., Wang, K., Mirams, G. R., Caruso, A., Noble, D., Walz, A., et al. (2016). Recent developments in using mechanistic cardiac modellingmodeling for drug safety evaluation. Drug Discov. Today 21, 924–938. doi: 10. 1016/j. drudis. 2016. 02. 003

Drovandi, C. C., Cusimano, N., Psaltis, S., Lawson, B. A., Pettitt, A. N., Burrage, P., et al. (2016). Sampling methods for exploring between-subject variability in cardiac electrophysiology experiments. J. R. Soc. Interface 13: 20160214. doi: 10. 1098/rsif. 2016. 0214

Elkins, R. C., Davies, M. R., Brough, S. J., Gavaghan, D. J., Cui, Y., Abi-Gerges, N., et al. (2013). Variability in high-throughput ion-channel screening data and consequences for cardiac safety assessment. J. Pharmacol. Toxicol. Methods 68, 112–122. doi: 10. 1016/j. vascn. 2013. 04. 007

Fanoe, S., Jensen, G. B., Sjøgren, P., Korsgaard, M. P., and Grunnet, M. (2009). Oxycodone is associated with dose-dependent QTc prolongation in patients and low-affinity inhibiting of hERG activity in vitro . Br. J. Clin. Pharmacol. 67, 172–179. doi: 10. 1111/j. 1365-2125. 2008. 03327. x

Fijorek, K., Patel, N., Klima, Ł., Stolarz-Skrzypek, K., Kawecka-Jaszcz, K., and Polak, S. (2013a). Age and gender dependent heart rate circadian model development and performance verification on the proarrhythmic drug case study. Theor. Biol. Med. Model. 10: 7. doi: 10. 1186/1742-4682-10-7

Fijorek, K., Püsküllüoglu, M., and Polak, S. (2013b). Circadian models of serum potassium, sodium, and calcium concentrations in healthy individuals and their application to cardiac electrophysiology simulations at individual level. Comput. Math. Methods Med. 2013: 429037. doi: 10. 1155/2013/429037

Fijorek, K., Püsküllüoglu, M., Tomaszewski, R., Tomaszewska, D., Glinka, A., and Polak, S. (2014a). Serum potassium, sodium and calcium levels in healthy individuals – literature review and data analysis. Folia. Med. Cracov. 54, 53–70.

Fijorek, K., Tanner, F. C., Stähli, B. E., Gielerak, G., Krzesinski, P., Uzieblo-Zyczkowska, B., et al. (2014b). Model of the distribution of diastolic left ventricular posterior wall thickess in healthy adults and its impact on the behavior of a string of virtual cardiomyocytes. J. Cardiovasc. Transl. Res. 7, 507–517. doi: 10. 1007/s12265-014-9558-4

Fink, M., Niederer, S. A., Cherry, E. M., Fenton, F. H., Koivumki, J. T., Seemann, G., et al. (2011). Cardiac cell modellingmodeling: observations from the heart of the cardiac physiome project. Prog. Biophys. Mol. Biol. 104, 2–21. doi: 10. 1016/j. pbiomolbio. 2010. 03. 002

Friedrich, C. M. (2016). A model qualification method for mechanistic physiological QSP models to support model-informed drug development. CPT Pharmacometrics Syst. Pharmacol. 5, 43–53. doi: 10. 1002/psp4. 12056

Gintant, G., Sager, P. T., and Stockbridge, N. (2016). Evolution of strategies to improve preclinical cardiac safety testing. Nat. Rev. Drug Discov. 15, 457–471. doi: 10. 1038/nrd. 2015. 34

Glinka, A., and Polak, S. (2013). An analysis of cardiomyocyte electrophysiology in the presence of hERG gene mutations. Bi. Algorithms Med. Syst. 9, 135–140. doi: 10. 1515/bams-2013-0017

Goldenberg, I., and Moss, A. J. (2008). Long QT syndrome. J. Am. Coll. Cardiol. 51, 2291–2300. doi: 10. 1016/j. jacc. 2008. 02. 068

Goldenberg, I., Zareba, W., and Moss, A. J. (2008). Long QT syndrome. Curr. Probl. Cardiol. 33, 629–694. doi: 10. 1016/j. cpcardiol. 2008. 07. 002

Hancox, J. C., Hasnain, M., Vieweg, W. V., Gysel, M., Methot, M., and Baranchuk, A. (2014). Erythromycin, QTc interval prolongation, and torsade de pointes: case reports, major risk factors and illness severity. Ther. Adv. Infect. Dis. 2, 47–59. doi: 10. 1177/2049936114527744

Hartmanshenn, C., Megerle Scherholz, M., and Androulakis, I. P. (2016). Physiologically-based pharmacokinetic models: approaches for enabling personalized medicine. J. Pharmacokinet. Pharmacodyn. 43, 481–504. doi: 10. 1007/s10928-016-9492-y

Hasnain, M., Vieweg, W. V. R., Howland, R. H., Kogut, C., Breden Crouse, E. L., Jayanthi, N., et al. (2014). Quetiapine, QTc interval prolongation, and torsade de pointes: a review of case reports. Ther. Adv. Psychopharmacol. 4, 130–138. doi: 10. 1177/2045125313510194

Hodgkin, A. L., and Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544. doi: 10. 1113/jphysiol. 1952. sp004764

Huang, S.-M., Abernethy, D. R., Wang, Y., Zhao, P., and Zineh, I. (2013). The utility of modeling and simulation in drug development and regulatory review. J. Pharm. Sci. 102, 2912–2923. doi: 10. 1002/jps. 23570

Hulhoven, R., Rosillon, D., Bridson, W. E., Meeus, M. A., Salas, E., and Stockis, A. (2008). Effect of levetiracetam on cardiac repolarization in healthy subjects: a single-dose, randomized, placebo- and active-controlled, four-way crossover study. Clin. Ther. 30, 260–270. doi: 10. 1016/j. clinthera. 2008. 02. 002

ICRP (2002). Basic anatomical and physiological data for use in radiological protection reference values. ICRP Publication 89. Ann. ICRP 32 .

Jamei, M. (2016). Recent advances in development and application of physiologically-based pharmacokinetic (PBPK) models: a transition from academic curiosity to regulatory acceptance. Curr. Pharmacol. Rep. 2, 161–169. doi: 10. 1007/s40495-016-0059-9

Jamei, M., Dickinson, G. L., and Rostami-Hodjegan, A. (2009a). A framework for assessing inter-individual variability in pharmacokinetics using virtual human populations and integrating general knowledge of physical chemistry, biology, anatomy, physiology and genetics: a tale of ‘ bottom-up’ vs. ‘ top-down’ recognition of covariates. Drug Metab. Pharmacokinet. 24, 53–75. doi: 10. 2133/dmpk. 24. 53

Jamei, M., Marciniak, S., Feng, K., Barnett, A., Tucker, G., and Rostami-Hodjegan, A. (2009b). The Simcyp population-based ADME simulator. Expert Opin. Drug Metab. Toxicol. 5, 211–223. doi: 10. 1517/17425250802691074

Jerling, M., and Abdallah, H. (2005). Effect of renal impairment on multiple-dose pharmacokinetics of extended-release ranolazine. Clin. Pharmacol. Ther. 78, 288–297. doi: 10. 1016/j. clpt. 2005. 05. 004

Johnstone, R. H., Chang, E. T. Y., Bardenet, R., de Boer, T. P., Gavaghan, D. J., Pathmanathan, P., et al. (2016). Uncertainty and variability in models of the cardiac action potential: can we build trustworthy models? J. Mol. Cell. Cardiol. 96, 49–62. doi: 10. 1016/j. yjmcc. 2015. 11. 018

Jones, H. M., Chen, Y., Gibson, C., Heimbach, T., Parrott, N., Peters, S. A., et al. (2015). Physiologically based pharmacokinetic modeling in drug discovery and development: a pharmaceutical industry perspective. Clin. Pharmacol. Ther. 97, 247–262. doi: 10. 1002/cpt. 37

Karjalainen, J., Viitasalo, M., Mänttäri, M., and Manninen, V. (1994). Relation between QT intervals and heart rates from 40 to 120 beats/min in rest electrocardiograms of men and a simple method to adjust QT interval values. J. Am. Coll. Cardiol. 23, 1547–1553. doi: 10. 1016/0735-1097(94)90654-8

Krantz, M. J., Kutinsky, I. B., Robertson, A. D., and Mehler, P. S. (2003). Dose-related effects of methadone on QT prolongation in a series of patients with torsade de pointes. Pharmacotherapy 23, 802–805. doi: 10. 1592/phco. 23. 6. 802. 32186

Lee, J. Y., Garnett, C. E., Gobburu, J. V., Bhattaram, V. A., Brar, S., Earp, J. C., et al. (2011). Impact of pharmacometric analyses on new drug approval and labelling decisions: a review of 198 submissions between 2000 and 2008. Clin. Pharmacokinet. 50, 627–635. doi: 10. 2165/11593210-000000000-00000

Li, Z., Dutta, S., Sheng, J., Tran, P. N., Wu, W., Chang, K., et al. (2017). Improving the in silico assessment of proarrhythmia risk by combining herg (Human ether-à-go-go-related gene) channel-drug binding kinetics and multichannel pharmacology. Circ. Arrhythm. Electrophysiol. 10: e004628. doi: 10. 1161/CIRCEP. 116. 004628

Liberos, A., Bueno-Orovio, A., Rodrigo, M., Ravens, U., Hernandez-Romero, I., Fernandez-Aviles, F., et al. (2016). Balance between sodium and calcium currents underlying chronic atrial fibrillation termination: an in silico intersubject variability study. Heart Rhythm. 13, 2358–2365. doi: 10. 1016/j. hrthm. 2016. 08. 028

Lu, G., Abduljalil, K., Jamei, M., Johnson, T. N., Soltani, H., and Rostami-Hodjegan, A. (2012). Physiologically-based pharmacokinetic (PBPK) models for assessing the kinetics of xenobiotics during pregnancy: achievements and shortcomings. Curr. Drug Metab. 13, 695–720. doi: 10. 2174/138920012800840374

Macfarlane, P. W., McLaughlin, S. C., Devine, B., and Yang, T. F. (1994). Effects of age, sex, and race on ECG interval measurements. J. Electrocardiol. 27(Suppl.), 14–19. doi: 10. 1016/S0022-0736(94)80039-1

Mann, S. A., Imtiaz, M., Winbo, A., Rydberg, A., Perry, M. D., Couderc, J. P., et al. (2016). Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes. J. Mol. Cell Cardiol . 100, 25–34. doi: 10. 1016/j. yjmcc. 2016. 09. 011

Marder, E., and Taylor, A. L. (2011). Multiple models to capture the variability in biological neurons and networks. Nat. Neurosci. 14, 133–138. doi: 10. 1038/nn. 2735

Massin, M. M., Maeyns, K., Withofs, N., Ravet, F., and Gerard, P. (2000). Circadian rhythm of heart rate and heart rate variability. Arch. Dis. Child 83, 179–182. doi: 10. 1136/adc. 83. 2. 179

Miller, R., Ewy, W., Corrigan, B. W., Ouellet, D., Hermann, D., Kowalski, K. G., et al. (2005). How modeling and simulation have enhanced decision making in new drug development. J. Pharmacokinet. Pharmacodyn. 32, 185–197. doi: 10. 1007/s10928-005-0074-7

Mirams, G. R., Davies, M. R., Cui, Y., Kohl, P., and Noble, D. (2012). Application of cardiac electrophysiology simulations to pro-arrhythmic safety testing. Br. J. Pharmacol. 167, 932–945. doi: 10. 1111/j. 1476-5381. 2012. 02020. x

Mirams, G. R., Pathmanathan, P., Gray, R. A., Challenor, P., and Clayton, R. H. (2016). Uncertainty and variability in computational and mathematical models of cardiac physiology. J. Physiol. 594, 6833–6847. doi: 10. 1113/JP271671

Mishra, H., Polak, S., Jamei, M., and Rostami-Hodjegan, A. (2014). Interaction between domperidone and ketoconazole: toward prediction of consequent QTc prolongation using purely in vitro information. CPT Pharmacometrics Syst. Pharmacol. 3: e130. doi: 10. 1038/psp. 2014. 26

Molnar, J., Zhang, F., Weiss, J., Ehlert, F. A., and Rosenthal, J. E. (1996). Diurnal pattern of QTc interval: how long is prolonged? Possible relation to circadian triggers of cardiovascular events. J. Am. Coll. Cardiol. 27, 76–83. doi: 10. 1016/0735-1097(95)00426-2

Morotti, S., and Grandi, E. (2016). Logistic regression analysis of populations of electrophysiological models to assess proarrythmic risk. MethodsX 4, 25–34. doi: 10. 1016/j. mex. 2016. 12. 002

Muszkiewicz, A., Britton, O. J., Gemmell, P., Passini, E., Sánchez, C., Zhou, X., et al. (2016). Variability in cardiac electrophysiology: using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm. Prog. Biophys. Mol. Biol. 120, 115–127. doi: 10. 1016/j. pbiomolbio. 2015. 12. 002

Noble, D. (1962). A modification of the Hodgkin–Huxley equations applicable to Purkinje fibre action and pace-maker potentials. J. Physiol. 160, 317–352. doi: 10. 1113/jphysiol. 1962. sp006849

Noble, D. (2007). From the Hodgkin-Huxley axon to the virtual heart. J. Physiol. 1, 15–22. doi: 10. 1113/jphysiol. 2006. 119370

O’Hara, T., Virag, L., Varró, A., and Rudy, Y. (2011). Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS ONE Comput. Biol. 7: e1002061. doi: 10. 1371/journal. pcbi. 1002061

Passini, E., Mincholé, A., Coppini, R., Cerbai, E., Rodriguez, B., Severi, S., et al. (2016). Mechanisms of pro-arrhythmic abnormalities in ventricular repolarisation and anti-arrhythmic therapies in human hypertrophic cardiomyopathy. J. Mol. Cell. Cardiol. 96, 72–81. doi: 10. 1016/j. yjmcc. 2015. 09. 003

Polak, S., and Fijorek, K. (2012). Inter-individual variability in the pre-clinical drug cardiotoxic safety assessment—analysis of the age–cardiomyocytes electric capacitance dependence. J. Cardiovasc. Transl. Res. 5, 321–332. doi: 10. 1007/s12265-012-9357-8

Polak, S., Fijorek, K., Glinka, A., Wisniowska, B., and Mendyk, A. (2012). Virtual population generator for human cardiomyocytes parameters. in silico drug cardiotoxicity assessment. Toxicol Mech. Methods 22, 31–40. doi: 10. 3109/15376516. 2011. 585477

Pueyo, E., Dangerfield, C. E., Britton, O. J., Virág, L., Kistamás, K., Szentandrássy, N., et al. (2016). Experimentally-based computational investigation into beat-to-beat variability in ventricular repolarization and its response to ionic current inhibition. PLoS ONE 11: e0151461. doi: 10. 1371/journal. pone. 0151461

Roberts, B. N., Yang, P. C., Behrens, S. B., Moreno, J. D., and Clancy, C. E. (2012). Computational approaches to understand cardiac electrophysiology and arrhythmias. Am. J. Physiol. Heart Circ. Physiol. 303, H766–H783. doi: 10. 1152/ajpheart. 01081. 2011

Romero, L., Pueyo, E., Fink, M., and Rodríguez, B. (2009). Impact of ionic current variability on human ventricular cellular electrophysiology. Am. J. Physiol. Heart Circ. Physiol. 297, H1436–H1445. doi: 10. 1152/ajpheart. 00263. 2009

Rostami-Hodjegan, A. (2012). physiologically based pharmacokinetics joined with in vitro-in vivo extrapolation of adme: a marriage under the arch of systems pharmacology. Clin. Pharmacol. Ther. 92, 50–61. doi: 10. 1038/clpt. 2012. 65

Rowland, M., Lesko, L. J., and Rostami-Hodjegan, A. (2015). Physiologically based pharmacokinetics is impacting drug development and regulatory decision making. CPT Pharmacometrics Syst. Pharmacol. 4, 313–315. doi: 10. 1002/psp4. 52

Sager, J. E., Yu, J., Ragueneau-Majlessi, I., and Isoherranen, N. (2015). Physiologically based pharmacokinetic (PBPK) modeling and simulation approaches: a systematic review of published models, applications, and model verification. Drug Metab. Dispos. 43, 1823–1837. doi: 10. 1124/dmd. 115. 065920

Saito, Y., Maekawa, K., and Ohno, Y. (2014). Novel Draft guideline for drug interaction studies in the drug development and labeling recommendations. Regul. Sci. Med. Prod. 4, 249–255. doi: 10. 14982/rsmp. 4. 249

Salem, F., Johnson, T. N., Barter, Z. E., Leeder, J. S., and Rostami-Hodjegan, A. (2013). Age related changes in fractional elimination pathways for drugs: assessing the impact of variable ontogeny on metabolic drug-drug interactions. J. Clin. Pharmacol. 53, 857–865. doi: 10. 1002/jcph. 100

Sánchez, C., Bueno-Orovio, A., Wettwer, E., Loose, S., Simon, J., Ravens, U., et al. (2014). Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation. PLoS ONE 9: e105897. doi: 10. 1371/journal. pone. 0105897

Sarkar, A. X., and Sobie, E. A. (2010). Regression analysis for constraining free parameters in electrophysiological models of cardiac cells. PLoS Comput. Biol. 6: e1000914. doi: 10. 1371/journal. pcbi. 1000914

Sarkar, A. X., and Sobie, E. A. (2011). Quantification of repolarization reserve to understand interpatient variability in the response to proarrhythmic drugs: a computational analysis. Heart Rhythm 8, 1749–1755. doi: 10. 1016/j. hrthm. 2011. 05. 023

Sarkar, A. X., Christini, D. J., and Sobie, E. A. (2012). Exploiting mathematical models to illuminate electrophysiological variability between individuals. J. Physiol. 590, 2555–2567. doi: 10. 1113/jphysiol. 2011. 223313

Sato, M., Ochiai, Y., Kijima, S., Nagai, N., Ando, Y., Shikano, M., et al. (2017). Quantitative modeling and simulation in pmda: a japanese regulatory perspective. CPT Pharmacometrics Syst. Pharmacol. 6, 413–415. doi: 10. 1002/psp4. 12203

Sennels, H. P., Jørgensen, H. L., Goetze, J. P., and Fahrenkrug, J. (2012). Rhythmic 24-hour variations of frequently used clinical biochemical parameters in healthy young males–the Bispebjerg study of diurnal variations. Scand. J. Clin. Lab. Invest. 72, 287–295. doi: 10. 3109/00365513. 2012. 662281

Shepard, T., Scott, G., Cole, S., Nordmark, A., and Bouzom, F. (2015). Physiologically based models in regulatory submissions: output from the ABPI/MHRA forum on physiologically based modeling and simulation. CPT Pharmacometrics Syst. Pharmacol. 4, 221–225. doi: 10. 1002/psp4. 30

Sobie, E. A. (2009). Parameter sensitivity analysis in electrophysiological models using multivariable regression. Biophys. J. 96, 1264–1274. doi: 10. 1016/j. bpj. 2008. 10. 056

Sorger, P. K., Allerheiligen, S. R., Abernethy, D. R., Altman, R. B., and Brouwer, K. L. (2011). Quantitative and systems pharmacology in the post-genomic era: new approaches to discovering drugs and understanding therapeutic mechanisms,” in An NIH White Paper by the QSP Workshop Group (Bethesda MD: NIH), 1–48.

Ten Tusscher, K. H., and Panfilov, A. V. (2006). Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions. Phys. Med. Biol. 51, 6141–6156. doi: 10. 1088/0031-9155/51/23/014

Teorell, T. (1937). Studies on the diffusion effect upon ionic distribution: II. Experiments on ionic accumulation. J. Gen. Physiol. 21, 107–122. doi: 10. 1085/jgp. 21. 1. 107

Tsamandouras, N., Rostami-Hodjegan, A., and Aarons, L. (2015). Combining the ‘ bottom up’ and ‘ top down’ approaches in pharmacokinetic modellingmodeling: fitting PBPK models to observed clinical data. Br. J. Clin. Pharmacol. 79, 48–55. doi: 10. 1111/bcp. 12234

Valentin, J. (2003). Basic Anatomical and Physiological Data for use in Radiological Protection: Reference Values . Oxford: Published for the International Commission on Radiological Protection by Pergamon.

Vandael, E., Vandenberk, B., Vandenberghe, J., Willems, R., and Foulon, V. (2017). Risk factors for QTc-prolongation: systematic review of the evidence. Int. J. Clin. Pharm. 39, 16–25. doi: 10. 1007/s11096-016-0414-2

Vieweg, W. V., Hasnain, M., Hancox, J. C., Baranchuk, A., Digby, G. C., Kogut, C., et al. (2013a). Risperidone, QTc interval prolongation, and torsade de pointes: a systematic review of case reports. Psychopharmacology 228, 515–524. doi: 10. 1007/s00213-013-3192-8

Vieweg, W. V., Hasnain, M., Howland, R. H., Clausen, T., Koneru, J. N., Kogut, C., et al. (2013b). Methadone, QTc interval prolongation and torsade de pointes: case reports offer the best understanding of this problem. Ther. Adv. Psychopharmacol. 3, 219–232. doi: 10. 1177/2045125312469982

Visser, S. A. G., de Alwis, D. P., Kerbusch, T., Stone, J. A., and Allerheiligen, S. R. B. (2014). Implementation of quantitative and systems pharmacology in large pharma. CPT Pharmacometrics Syst. Pharmacol. 3: e142. doi: 10. 1038/psp. 2014. 40

Wagner, C., Zhao, P., Pan, Y., Hsu, V., Grillo, J., Huang, S. M., et al. (2015). Application of physiologically based pharmacokinetic (PBPK) modeling to support dose selection: report of an FDA public workshop on PBPK. CPT Pharmacometrics Syst. Pharmacol. 4, 226–230. doi: 10. 1002/psp4. 33

Wiśniowska, B., and Polak, S. (2016). Virtual clinical trial towards polytherapy safety assessment – combination of PBPK/PD based modeling and simulation approach with DDIs involving terfenadine as an example. J. Pharm. Sci. 105, 3415–3424. doi: 10. 1016/j. xphs. 2016. 08. 002

Wiśniowska, B., Tylutki, Z., Wyszogrodzka, G., and Polak, S. (2016). Drug-drug interactions and QT prolongation as a commonly assessed cardiac effect – comprehensive overview of clinical trials. BMC Pharmacol. Toxicol. 17: 12. doi: 10. 1186/s40360-016-0053-1

Wysowski, D. K., Corken, A., Gallo-Torres, H., Talarico, L., and Rodriguez, E. M. (2001). Postmarketing reports of QT prolongation and ventricular arrhythmia in association with cisapride and food and drug administration regulatory actions. Am. J. Gastroenterol. 96, 1698–1703. doi: 10. 1111/j. 1572-0241. 2001. 03927. x

Zhao, P., Zhang, L., Grillo, J. A., Liu, Q., Bullock, J. M., Moon, Y. J., et al. (2011). Applications of physiologically based pharmacokinetic (PBPK) modeling and simulation during regulatory review. Clin. Pharmacol. Ther. 89, 259–267. doi: 10. 1038/clpt. 2010. 298

Zhou, X., Bueno-Orovio, A., Orini, M., Hanson, B., Hayward, M., Taggart, P., et al. (2016). In vivo and in silico investigation into mechanisms of frequency dependence of repolarization alternans in human ventricular cardiomyocytes. Circ. Res. 118, 266–278. doi: 10. 1161/CIRCRESAHA. 115. 307836

Zineh, I., Abernethy, D., Hop, C., Bello, A., McClellan, M. B., Daniel, G. W., et al. (2017). Improving the tools of clinical pharmacology: goals for 2017 and beyond. Clin. Pharmacol. Ther. 101, 22–24. doi: 10. 1002/cpt. 530