|Inter-individual variability in the pre-clinical drug cardiotoxic safety assessment--analysis of the age-cardiomyocytes electric capacitance dependence.|
|Jump to Full Text|
|PMID: 22411323 Owner: NLM Status: MEDLINE|
|Electrical phenomena located within the plasma membrane of the mammalian cardiac cells are connected with the cells' main physiological functions--signals processing and contractility. They were extensively studied and described mathematically in so-called Hodgkin-Huxley paradigm. One of the physiological parameters, namely cell electric capacitance, has not been analyzed in-depth. The aim of the study was to validate the mechanistic model describing the capacitive properties of cells, based on a collected experimental dataset which describes the electric capacitance of human ventricular myocytes. The gathered data was further utilized for developing an empirical correlation between a healthy individual's age and cardiomyocyte electric capacitance.|
|Sebastian Polak; Kamil Fijorek|
Related Documents :
|8043593 - Effects of modulations of the energetic metabolism on the mortality of cultured cells.
18193173 - Enhancing glutathione synthesis can decrease zinc-mediated toxicity.
17271993 - A tool for the quantitative spatial analysis of mammary gland epithelium.
|Type: Journal Article; Research Support, Non-U.S. Gov't; Validation Studies Date: 2012-03-13|
|Title: Journal of cardiovascular translational research Volume: 5 ISSN: 1937-5395 ISO Abbreviation: J Cardiovasc Transl Res Publication Date: 2012 Jun|
|Created Date: 2012-05-11 Completed Date: 2012-09-10 Revised Date: 2013-06-26|
Medline Journal Info:
|Nlm Unique ID: 101468585 Medline TA: J Cardiovasc Transl Res Country: United States|
|Languages: eng Pagination: 321-32 Citation Subset: IM|
|Unit of Pharmacoepidemiology and Pharmacoeconomics, Faculty of Pharmacy, Medical College, Jagiellonian University, Medyczna 9 Str, 30-688, Cracow, Poland. email@example.com|
|APA/MLA Format Download EndNote Download BibTex|
Aged, 80 and over
Cell Membrane / drug effects*, pathology
Drug Evaluation, Preclinical / methods*
Heart Diseases / chemically induced*, pathology
Myocytes, Cardiac / drug effects*, pathology
Reproducibility of Results
Journal ID (nlm-ta): J Cardiovasc Transl Res
Journal ID (iso-abbrev): J Cardiovasc Transl Res
Publisher: Springer US, Boston
© The Author(s) 2012
Received Day: 7 Month: 12 Year: 2011
Accepted Day: 29 Month: 2 Year: 2012
Electronic publication date: Day: 13 Month: 3 Year: 2012
pmc-release publication date: Day: 13 Month: 3 Year: 2012
Print publication date: Month: 6 Year: 2012
Volume: 5 Issue: 3
First Page: 321 Last Page: 332
PubMed Id: 22411323
Publisher Id: 9357
|Inter-individual Variability in the Pre-clinical Drug Cardiotoxic Safety Assessment—Analysis of the Age–Cardiomyocytes Electric Capacitance Dependence|
Address: +48-12-6205517 +48-12-6205820 firstname.lastname@example.org
1Unit of Pharmacoepidemiology and Pharmacoeconomics, Faculty of Pharmacy, Medical College, Jagiellonian University, Medyczna 9 Str., 30-688 Cracow, Poland
2Department of Statistics, Faculty of Management, Cracow University of Economics, Rakowicka 27 Str., 31-510 Cracow, Poland
Electrical phenomena located within the plasma membrane of the mammalian cardiac cells are connected with the cell’s main physiological functions—signals processing and contractility. Both of these elements are well-known from an anatomical, physiological, and functional point of view even at the molecular and submolecular level. The first area, namely electrical activity has been studied and described in the form of an equation derived by Hodgkin and Huxley [Eq. 1] which is described as the H-H paradigm . This paradigm is extensively used as a base for multiple models describing heart cells’ electric activity [2–4].
|Vm||Membrane voltage [mV]|
|Iion||Sum of ionic currents [pA/pF]|
|Istim||Stimulation potential [pA]|
|Cm||Cell electric capacitance [pF]|
Models based on the Hodgkin and Huxley paradigm are used in a wide range of areas, from heart diseases and disorders’ simulation to drug cardiotoxicity assessment [5, 6]. The ionic activity (Iion) is subdivided into separate ionic currents described with the use of various mathematical techniques [7, 8]. Being the main triggering points for arrhythmias (including drug-dependent ones) they were widely investigated [9–12]. A large amount of time and effort has been put into studying the ionic activity and to a lesser degree, the stimulation potential. However, cell electric capacitance has not been studied as intensely.
The cell membrane, which consists of a double lipid layer shows a great resistance to the flow of ionic currents. As an effect of the ion’s separation, the electric potential across the membrane is formed. Maintaining this status is possible through the action of the ionic pumps, including the Na+/K+ exchanger. The Na+/K+-ATPase for every three sodium ions pumped out from the cell, transports only two potassium ions into the cell. The electromagnetic field across the membrane causes the attraction of the charged particles. Anions from the extracellular fluid are attracted to the outer side of the membrane. At the same time cations from cytosol accumulate inside. The consequences of such phenomena are demonstrated in other biophysical measures characterizing biological cells, including conductance and capacitance. Capacitance (Cm) defines the amount of electric charge that can be stored in a capacitor at a given voltage. In the biological membranes there are two spaces, which could be recognized as plates, formed by the cytosol and extracellular fluid. These are electricity conductors, separated by an insulator (the phospholipid bilayer). In such a situation, the membrane behaves electrically as a capacitor. Capacitance measurements can be used to study the changes in cell surface area. This is because capacitance is proportional to the cell surface area and therefore depends on the shape of the cell affecting the area . The equation describing this relationship is as follows :
|Cm||Cell electric capacitance [pF]|
|A||Cell surface area [μm2]|
|εr||Dielectric constant of the cell membrane [dimensionless]|
|ε0||Electric constant (vacuum permittivity) [F/m]|
|d||Thickness of the membrane [nm]|
The model presented above can be treated as a mechanistic representation of the observed phenomenon.
In a recent study described elsewhere, it was hypothesized that the electric capacitance of the cardiac cells depends on the cell surface area. Furthermore, according to the study, the cell surface area is subsequently dependent on the age of the individuals . Realistic parameter values are needed in order to mimic the inter-individual variability in the in silico environment and also to simulate cell membrane electric activity. It is especially important when dealing with excitable cells like cardiomyocytes. Their mathematical expressions have to have support in the physiology. This is a crucial factor for the realistic simulations of drugs and pathophysiological-condition-dependent alterations. The above-mentioned previously developed and described model combines the individuals’ age as an independent variable and the cardiomyocyte volume and area as the dependent variables. Based on that correlation, every member of the randomly built virtual population is described by a set of parameters including membrane capacitance calculated with the use of Eq. 2. This equation expresses a mechanistic model whose performance should be further verified. Verification can be done by comparing the results of the simulations with the real values, derived from human cardiac cells, of the cardiomyocyte electric capacitance.
The main aim of the study was to verify the predictive quality of the previously established correlation, combining healthy individuals’ age and cardiomyocyte volume/area/electric capacitance, with the use of the collected and analyzed data describing the electric capacitance of human ventricular myocytes.
The same data, selected based on a strictly prepared procedure, was then used to achieve an additional study aim—a derivation of the novel empirical models, correlating: the healthy individual’s age and cardiomyocyte electric capacitance and the ages and cardiomyocyte electric capacitance of individuals suffering from various heart diseases.
The third aim includes the assessment of the influence of the above-mentioned approaches on the simulated action potential duration.
In order to collect data describing human cardiomyocyte electric capacitance values, database searches were performed. The Scopus, Medline, and Google Scholar databases were searched. There was no time limit for the search query. The key phrases, included in either the article title, keywords, or abstract, were: ‘cardiomyocyte’, ‘cardiomyocyte parameters’, ‘cardiomyocyte capacitance’, ‘cardiac myocyte’, ‘cardiac cell’, ‘ventricular myocyte’, ‘capacitance’, ‘human’, ‘volunteers’, ‘healthy’, and ‘control’. The reports that were found were scrutinized for additional references not covered by the electronic search. Native data was extracted in various ways, depending on the article structure and data presentation method. If not stated otherwise, all data values are presented as algebraic means with the corresponding standard error of the mean (SE). The collected data can be subdivided into two groups—patient data and cell data. The first group covers: sex (male, female), cells’ location (left ventricle, right ventricle; epicardium, endocardium), disease occurrence (healthy, diseased), number of individuals in the study, average, minimum, and maximum age. The latter group includes: measurement temperature and number of probes. Information describing age is presented in years [years] and electric capacitance is described in picofarads [pF].
The collected data was further divided according to the above described parameters. The crucial factor is health status, as the previously established age–cardiomyocyte volume/area–electric capacitance correlation was developed based on healthy individuals’ data. Considering this, it was assumed that the performance verification of the model would be done with the use of identically pre-processed data (healthy individuals only), although a full set of collected data is presented.
The data used during the new empirical models development was carefully analyzed. Inclusion criteria were as follows: clearly defined health status (either health or disease), information regarding average age, and electric capacitance. There are some contradictory reports, although according to the majority of the previously published papers, it was assumed that the sex, cell location, and measurement temperature do not significantly influence the electric capacitance [16, 17]. Two studies describing age-capacitance correlation were based on the biological material derived from young individuals [18, 19]. Both of them describe non-healthy children diagnosed with the tetralogy of Fallot, which is connected with right ventricular hypertrophy [20, 21]. As we were not able to find separate reports describing healthy young individuals, they were intentionally included in the healthy population. This being due to the fact that this was the only available information.
Non-healthy individuals were characterized by various diseases, but the common feature was heart insufficiency. This can furthermore be connected with heart cells' hypertrophy. According to that they were separately analyzed with the use of similar techniques, and the results were compared with those derived from healthy individuals. For the modeling stage, only individuals with an evidently defined health status were included to assure model consistency. No disease-specific analysis was performed, as the data amount for various diseases was not sufficient and the clinical condition was not always clearly defined in the available sources.
All simulations were done either with the use of previously published models and tools or models derived as a part of this study. The tool which enables the generation of virtual populations, with the use of various approaches, was implemented as an easy-to-use Excel file. This is freely available and can be found in the supplementary materials attached to the current article. All simulations were run in ten replications. These were based on the number of individuals involved in the different studies and also on the characteristics regarding the age and gender of the individuals. Such a procedure was introduced to avoid randomization bias. The parameters describing the ages of the groups of individuals (mean, SD), which were derived directly from the analyzed studies, were then used to fully mimic the described populations. The Weibull model was chosen as the default for the description of the age distribution in a healthy population .
Additionally, the percentage of simulation results’ fold errors, which fell below an absolute value of 2, was presented.
The bias of the prediction model was assessed with the use of the geometric mean fold error measure (MFE). This is presented in Eq. 4. Such a measure was chosen due to its characteristics. These allow for the equal weighting of the under and over estimates. A mean fold error ≤2 is generally considered successful, indicating that most of the predicted values fall within the twofold range of the observed ones (between 0.5 and 2.0) [23, 24].
The root mean squared error (RMSE), as presented in Eq. 5, was used to assess the accuracy of the predictions.
It was assumed that the cardiomyocyte capacitance (Cm) distribution has non-negative support (capacitance cannot take negative values). It was also assumed that the capacitance distribution is skewed to the right. The most popular distribution showing this kind of behavior is log-normal which was tested as a default model with the assumption that the natural logarithm of the capacitance can be modeled as a normal random variable with a constant variance and age-dependent mean (log-linear model) as presented in Eq. 6:
The second tested model was a double-log-linear regression model. The model assumes that the natural logarithm of the capacitance can be modeled as a normal random variable. This normal random variable has a constant variance and mean that depends on the natural logarithm of the age (Eq. 7):
Both models parameters were estimated using ordinary least squares. The lm function, available in the R system for statistical computing, was used .
The simulation studies were done with the use of the tenTusscher model of the human left ventricle myocyte . The transmembrane potential was chosen as an investigated endpoint. M cell variation (midmyocardium) was used. Two parameters were modified to apply for physiological changes—the volume of cardiomyocyte (Vc) and electric capacitance (Cm). Volume of the cardiomyocytes was simulated with the previously published tool .
Moxifloxacin, the most commonly used positive control in thorough QTc (TQTc) studies, was chosen as a model drug to simulate the electric capacitance influence on the drug-triggered APD modification . The in vitro data, describing the concentration-dependent ionic currents inhibition, was gathered from the available literature. For the hERG-transfected HEK cells, at a physiological temperature and potassium ions’ concentration similar to the physiological (4 mM) environment which mimics human cardiomyocytes, the reported IC50 value was 35.7 μM . Two different, reported as the maximal in the healthy volunteers clinical trials concentrations, were tested—M1 = 5.57 μM and M2 = 11.21 μM [26, 28]. Assuming that the concentration–inhibition relation is described by the Hill equation (n = 1), the predicted IKr potassium current inhibition was calculated as being 13% and 24%, respectively.
In all cases, the simulation time was set to 500 ms.
R environment was used for the statistical computing. ToxComp system was used to simulate human ventricular cardiomyocyte action potential .
Table 1 presents the full set of data collected during the literature review process. Information is grouped according to the health status. There were 24 studies found combining age and the cardiac myocyte electric capacitance in humans. Among them nine studies fulfilled the inclusion criteria for healthy individuals and were included at the model building stage (Table 2). One study done with the use of biological material derived from the healthy individuals was removed during the analysis as there was no age clearly defined .
The cumulative results of the simulations (sim) are presented in Fig. 1. The results were created with the use of the model described by Eq. 2, in combination with the age–cardiomyocyte volume/area model and then compared with the observed data. The graph presents the average values derived from ten independent simulations linked with the corresponding observed values for all studies separately. Detailed results with error measures are presented in Tables 3 and 4.
The detailed results of simulations mimicking single studies are presented separately and attached as supplementary materials.
Empirical representation of the direct correlation between age and human left ventricle cardiomyocyte electric capacitance was based on the log-linear and double-log-linear models. The estimates of the model’s parameters are presented below [Tables 5 and 6].
Both the estimate of the intercept and the slope are significant (p value < 0.05). It can be inferred that a 1-year increase in age is associated, on average, with an increase of about 0.0115 and 0.022 in log cardiomyocyte capacitance for healthy and diseased individuals, respectively. The residual standard error (the variation of the individual points about the regression line) is 0.3323 and 0.3642. For the models representing healthy and heart-diseased individuals’ the results were R2 = 0.4 and 0.56. This means that 40% and 56% of the variability observed in log cardiomyocyte capacitance can be explained by the estimated regression equation for both populations.
Both the estimate of the intercept and the slope are significant (p value < 0.05). It can be inferred that a 1-year increase in log-age is associated, on average, with an increase of about 0.2522 and 0.3702 in log cardiomyocyte capacitance for healthy and diseased individuals, respectively. The residual standard error (the variation of individual points about the regression line) is 0.3116 and 0.3582. R2 = 0.47 and 0.58, or 47% and 58% of the variability observed in log cardiomyocyte capacitance can be explained by the estimated regression equation for both populations.
The plot below [Fig. 2] shows the relationship between cardiomyocyte capacitance and age, along with a regression curve and 95% pointwise prediction bands for the log-linear model. It is worth noticing that the prediction intervals are asymmetric (right skewed), and they do not allow for negative values.
The plot below [Fig. 3] shows the relationship between cardiomyocyte capacitance and age along with a regression curve and 95% pointwise prediction bands for the double-log-linear model.
A simulation study was performed to examine the influence of the demographic (age) and physiological parameters. These were parameters were the cardiomyocyte volume and the capacitance on the cardiac myocyte transmembrane potential. The remaining parameters were assumed to be constant and had default values (i.e., sarcoplasmic reticulum volume). The APD90 (action potential at 90% of repolarization) was chosen as a measure.
As presented in Eq. 2, the mechanistic model describes the phenomenon responsible for the electric charge storage at the phospholipid border. This border is lying between intra- and extracellular environments. Membrane capacitor efficiency is proportional to the cell surface area and therefore depends on the shape of the cell affecting the area, as well as the membrane thickness. Assuming that the membrane thickness and the dielectric constant of the cell membrane have constant values, then the cell area becomes a pivotal parameter regulating electric capacitance. The combination of the above-mentioned model with empirical correlation combining age and cardiomyocyte area into one semi-mechanistic system can be effectively used in practice . Such an approach has some practical advantages giving the opportunity to have a look inside the process; although its accuracy has to be validated. To perform the validation process, experimental data describing the electric capacitance of human cardiomyocytes was collected. To simulate the described populations, a previously developed virtual population generator was used. The parameters describing the age distribution, namely the mean value and the standard deviation of the mean, were utilized. For seven, from the 11 simulated studies, the fold difference between the observed and predicted data fell below 2. This is generally recognized as the acceptance level. As the studies were simulated in ten replicates, additional accuracy measure was expressed as a percentage number of virtual studies with a fold error below 2. Eight, from the 11 studies, have a minimum of 50% of the simulations which fulfilled such requirements. Such results allow for the making of a statement regarding the high accuracy of the tested tool. It is worth noticing that all the described observations are valid for the populations of healthy individuals which were taken for analysis after applying strict including/excluding criteria. During drug virtual clinical trials, further practical use of the described and validated system should be limited to healthy volunteers. At the same time, it can be clearly seen that the mechanistic model does not differentiate between cells from various locations. Such results may be an indirect indication that such differentiation could be justified.
In nine, from the 11, studies mechanistic model, tested in current study, under predicted the cardiomyocyte electric capacitance. It suggests bias, which can be an effect of the previously mentioned assumption, namely the cylindrical shape of the ventricular cardiac myocytes. Membrane invaginations (transverse tubules) increase the total area of the cell when comparing with a cylinder, and thus increase the electric capacitance. Assuming that ca. 30% of the membrane is present in the t-tubules, a proportional increase in the capacitance can be expected , which is supported by the current study’s findings.
A separate leg of the study was focused on the development of an empirical model which directly correlates age and human cardiomyocyte electric capacitance. Separate models describing healthy volunteers and individuals suffering from heart diseases were developed. In both cases, based on the observations, it was assumed that the capacitance distribution is skewed to the right and has non-negative support. It was also assumed that log-normal and double-log-normal distribution were chosen as the default models.
The double-log-linear model has been presented for information’s sake, as it fits the observed data better than the other tested models; although the log-linear model seems to be more physiological and is recommended by the authors. Due to the model structure, electric capacitance equals to zero (Cm = 0) is assumed for age zero (newborns). This has no support in the physiology.
The main limitation of the proposed models, describing healthy individuals result from the assumptions made at the data analysis stage. The lack of available information in the literature regarding the electric capacitance in non-adults, created the necessity to include data from young non-healthy individuals. The simulation study results show that there was up to a fivefold under prediction, when using the mechanistic model for young individuals. It clearly points out the likely bias of the subsequently developed empirical model. Information describing electric capacitance derived from the young, unaltered, ventricular myocytes would be necessary on order to improve the model’s predictive quality.
Separate models, describing the age and cardiomyocytes electric capacitance for the population with heart diseases, were developed. In all cases, the data was derived from patients suffering from various heart diseases. Regardless of the type of pathology, most of the conditions affecting the heart are connected with cardiac cell hypertrophy. Even when considering the loss of the transverse tubules observed in such conditions, it was expected that heart cell hypertrophy could significantly influence the electric capacitive capabilities of the myocyte, and thus the whole electric behavior. As the elderly population, highly affected by heart disorders, is strongly exposed to drugs, the possibility to virtually test the chemical’s influence on the heart action with the use of reliable data and models is highly beneficial.
The developed models can be considered as being complementary to the previously developed mechanistic model. The latter one is a part of the computational system for the drug cardiotoxicity assessment, based on the mathematical model for human ventricular cardiomyocyte . Both the empirical and mechanistic models represent two approaches, described as top–down and bottom–up, respectively. Both of them have their advantages and disadvantages. As presented in the current study, the direct correlation between age and cell electric capacitance fits better to the experimental data. Therefore, this offers better interpolation capabilities. This is obvious, considering that the same data was used during their development. At the same time they do not offer direct insight into the mechanism lying behind the described phenomena, which is a major advantage of the mechanistic model. The expected increase in the cardiomyocyte area and the following electric capacitance, after implementing a different volume calculation method, can improve the mechanistic model’s predictability. Taking the above-mentioned arguments into account, the authors suggest that the mechanistic model (part of the ToxComp system) should be utilized.
The simulations results presented in Tables 7 and 8 prove that the combination of the physiological parameters (in this case—the cardiomyocyte volume to electric capacitance ratio) ought to be analyzed rather than just single variables. The larger the above-mentioned ratio the higher the APD90 value and the total action potential duration. The combination of such a high-risk element with a concomitantly given drug or drug combination can result in a health risk increase (arrhythmia). It demonstrates the need for an inter-individual variability assessment during the pre-clinical, in-silico-based, drug safety studies.
The latter one, namely the in silico drugs toxicity assessment, offers direct clinical translation of the proposed model. The previously mentioned ToxComp system offers an in silico realized platform for the in vitro–in vivo extrapolation of the cardiotoxic effect with an inter-individual variability analysis. The pharmaceutical industry is viably interested in early cardiac safety assessment, the mechanism of which lies in the acquired long QT syndrome. The wide clinical evidence regarding drug-triggered heart rate disruption has resulted in compulsory cardiac risk assessment during the drug development process. This is strictly required by drug agencies and described in the available guidelines [32, 33]. The best known model of the general population, which is the final target of any drug which is under development, is a cohort included in a specific clinical trial (TQTc study—thorough QTc study) where the heart rate corrected change of the QT interval length is the final point. There are, although, some limitations of different natures—meritorious (the clinical trial is still a model!) and most importantly, financial. The latter one is essential at the early stage of drug development. For such reasons surrogates are recommended. They include in vivo animal studies (ventricular repolarization disruption assessment) and the in vitro evaluation of the drug effects on the ionic current. This is done through a native or heterologously expressed IKr channel protein, such as those encoded by hERG. In both cases an essential question arises—what is the translation between in vitro/animal in vivo results and the human situation? In that regard, the in vitro–in vivo extrapolation approach, with the use of the mathematical models and simulated virtual population, may be a useful solution. Such methods were successfully applied in the pharmacokinetic drug characterization .
Research was funded by the Polish National Centre for Research and Development LIDER (project number 187/L-1/09).
This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
|1..||Hodgkin AL,Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerveJournal of PhysiologyYear: 195211750054412991237|
|2..||Bernus O,Wilders R,Zemlin CW,Verschelde H,Panfilov AV. A computationally efficient electrophysiological model of human ventricular cellsAmerican Journal of Physiology—Heart and Circulation PhysiologyYear: 2002282H2296H2308|
|3..||Luo C,Rudy Y. A model of the ventricular cardiac action potential, depolarization, repolarization, and their interactionCirculation ResearchYear: 199168150115261709839|
|4..||Tusscher KHWJ,Noble D,Noble PJ,Panfilov AV. A model for human ventricular tissueAmerican Journal of Physiology—Heart and Circulation PhysiologyYear: 20042861573158910.1152/ajpheart.00794.2003|
|5..||Moreno JD,Clancy CE. Using computational modeling to predict arrhythmogenesis and antiarrhythmic therapyDrug Discovery Today: Disease ModelsYear: 20096718410.1016/j.ddmod.2010.03.00120652086|
|6..||Obiol-Pardo C,Gomis-Tena J,Sanz F,Saiz J,Pastor MA. Multiscale simulation system for the prediction of drug-induced cardiotoxicityJournal of Chemical Information and ModelingYear: 20115148349210.1021/ci100423z21250697|
|7..||Fink M,Giles WR,Noble D. Contributions of inwardly rectifying K+ currents to repolarization assessed using mathematical models of human ventricular myocytesPhilosophical Transactions of the Royal Society AYear: 20063641207122210.1098/rsta.2006.1765|
|8..||Fink M,Noble D. Markov models for ion channels: Versatility versus identifiability and speedPhilosophical Transactions of the Royal Society AYear: 20093672161217910.1098/rsta.2008.0301|
|9..||Lenkowski PW,Ko SH,Anderson JD,Brown ML,Patel MK. Block of human NaV15 sodium channels by novel alpha-hydroxyphenylamide analogues of phenytoin EuropeanJournal of Pharmaceutical SciencesYear: 20042163564410.1016/j.ejps.2004.01.004|
|10..||Lloyd J,Schmidt JB,Rovnyak G,Ahmad S,Atwal KS,Bisaha SN,Doweyko LM,Lodge NJ. Design and synthesis of 4-substituted benzamides as potent selective and orally bioavailable I(Ks) blockersJournal of Medicinal ChemistryYear: 2001443764376710.1021/jm015505u11689063|
|11..||Polak S,Wisniowska B,Brandys J. Collation assessment and analysis of literature in vitro data on hERG receptor blocking potency for subsequent modeling of drugs cardiotoxic propertiesJournal of Applied ToxicologyYear: 20092918320610.1002/jat.139518988205|
|12..||Wuest M,Christ T,Hiller N,Braeter M,Ravens U. Effects of three metabolites of propiverine on voltage-dependent L-type calcium currents in human atrial myocytesEuropean Journal of PharmacologyYear: 2008598949710.1016/j.ejphar.2008.08.03518796302|
|13..||Lollikea K,Lindaub M. Membrane capacitance techniques to monitor granule exocytosis in neutrophilsJournal of Immunological MethodsYear: 199923211112010.1016/S0022-1759(99)00169-610618513|
|14..||Molleman A. Patch clamping: An introductory guide to patch clamp electrophysiologyYear: 2005New YorkWiley|
|15..||Polak S,Fijorek K,Glinka A,Wisniowska B,Mendyk A. Virtual population generator for human cardiomyocytes parameters in silico drug cardiotoxicity assessmentToxicology Mechanisms and MethodsYear: 201122314010.3109/15376516.2011.58547722150010|
|16..||Nabauer M,Beuckelmann DJ,Uberfuhr P,Steinbeck G. Regional differences in current density and rate-dependent properties of the transient outward current in subepicardial and subendocardial myocytes of human left ventricleCirculationYear: 1996931681778616924|
|17..||Verkerk A,Veldkamp M,Ginneken A,Bouman L. Biphasic response of action potential duration to metabolic inhibition in rabbit and human ventricular myocytes: Role of transient outward current and ATP-regulated potassium currentJournal of Molecular and Cellular CardiologyYear: 1996282443245610.1006/jmcc.1996.02379004161|
|18..||Pelzmann B,Schaffer P,Bernhart E,Lang P,Machler H,Rigler B,Koidl B. L-type calcium current in human ventricular myocytes at a physiological temperature from children with tetralogy of FallotCardiovascular ResearchYear: 19983842443210.1016/S0008-6363(98)00002-99709403|
|19..||Schaffer P,Pelzmann B,Bernhart E,Lang P,Machler H,Rigler B,Koidl B. Repolarizing currents in ventricular myocytes from young patients with tetralogy of FallotCardiovascular ResearchYear: 19994333234310.1016/S0008-6363(99)00118-210536663|
|20..||Apitz C,Webb GD,Redington AN. Tetralogy of fallotLancetYear: 20093741462147110.1016/S0140-6736(09)60657-719683809|
|21..||Farah MC,Castro CR,Moreira VM,Binotto MA,Guerra VC,Riso AA,Marcial MB,Lopes AA,Mathias W Jr,Aiello VD. The impact of preexisting myocardial remodeling on ventricular function early after tetralogy of Fallot repairJournal of the American Society of EchocardiographyYear: 20102391291810.1016/j.echo.2010.06.00820650609|
|22..||Jentsch-Ullrich K,Koenigsmann M,Mohren M,Franke A. Lymphocyte subsets’ reference ranges in an age- and gender-balanced population of 100 healthy adults—A monocentric German studyClinical ImmunologyYear: 200511619219710.1016/j.clim.2005.03.02015993366|
|23..||Obach RS,Baxter JG,Liston TE,Silber BM,Jones BC,MacIntyre F,Rance DJ,Wastall P. The prediction of human pharmacokinetic parameters from preclinical and in vitro metabolism dataJournal of Pharmacology and Experimental TherapeuticsYear: 199728346589336307|
|24..||Parrott N,Paquereau N,Coassolo P,Lave T. An evaluation of the utility of physiologically based models of pharmacokinetics in early drug discoveryJournal of Pharmaceutical SciencesYear: 2005942327234310.1002/jps.2041916136543|
|25..||R Development Core Team R: A language and environment for statistical computing R foundation for statistical computing 2008 Vienna Austria http://wwwR-projectorg/.|
|26..||Bloomfield DM,Kost JT,Ghosh K,Hreniuk D,Hickey LA,Guitierrez MJ,Gottesdiener K,Wagner JA. The effect of moxifloxacin on QTc and implications for the design of thorough QT studiesClinical Pharmacology and TherapeuticsYear: 20088447548010.1038/clpt.2008.3319238652|
|27..||Chen X,Cass JD,Bradley JA,Dahm CM,Sun Z,Kadyszewski E,Engwall MJ,Zhou J. QT prolongation and proarrhythmia by moxifloxacin: Concordance of preclinical models in relation to clinical outcomeBritish Journal of PharmacologyYear: 200514679279910.1038/sj.bjp.070638916158069|
|28..||MacGowan AP. Moxifloxacin (Bay 12-8039): A new methoxy quinolone antibacterialExpert Opinion in Investigative DrugsYear: 1999818119910.1517/135437184.108.40.206|
|30..||Stengl M,Carmeliet E,Mubagwa K,Flameng W. Modulation of transient outward current by extracellular protons and Cd2+ in rat and human ventricular myocytesThe Journal of PhysiologyYear: 199851182783610.1111/j.1469-7793.1998.827bg.x9714863|
|31..||Page E. Quantitative ultrastructural analysis in cardiac membrane physiologyAmerican Journal of Physiology—Cell PhysiologyYear: 1978235C147C158|
|32..||ICH Harmonised Tripartite Guideline. The non-clinical evaluation of the potential for delayed ventricular repolarization (Qt Interval Prolongation) by human pharmaceuticals. S7B (http://www.ich.org/products/guidelines/safety/article/safety-guidelines.html).|
|33..||ICH Harmonised Tripartite Guideline. The clinical evaluation of Qt/Qtc interval prolongation and proarrhythmic potential for non-antiarrhythmic drugs. E14 (http://www.ich.org/products/guidelines/efficacy/article/efficacy-guidelines.html).|
|34..||Rostami-Hodjegan A,Tucker GT. Simulation and prediction of in vivo drug metabolism in human populations from in vitro dataNature Reviews: Drug DiscoveryYear: 2007614014810.1038/nrd217317268485|
|35..||Amos GJ,Wettwer E,Metzger F,Li Q,Himmel HM,Ravens U. Differences between outward currents of human atrial and subepicardial ventricular myocytesThe Journal of PhysiologyYear: 199649131509011620|
|36..||Piacentino V,Weber CR,Chen X,Weisser-Thomas J,Margulies KB,Bers DM,Houser SR. Cellular basis of abnormal calcium transients of failing human ventricular myocytesCirculation ResearchYear: 20039265165810.1161/01.RES.0000062469.83985.9B12600875|
|37..||Sipido KR,Stankovicova T,Flameng W,Vanhaecke J,Verdonck F. Frequency dependence of Ca2+ release from the sarcoplasmic reticulum in human ventricular myocytes from end-stage heart failureCardiovascular ResearchYear: 19983747848810.1016/S0008-6363(97)00280-09614502|
|38..||Wettwer E,Amos GJ,Posival H,Ravens U. Transient outward current in human ventricular myocytes of subepicardial and subendocardial originCirculation ResearchYear: 1994754734828062421|
|39..||Verkerk AO,Veldkamp MW,Jonge N,Wilders R,Ginneken AC. Injury current modulates afterdepolarizations in single human ventricular cellsCardiovascular ResearchYear: 20004712413210.1016/S0008-6363(00)00064-X10869538|
|40..||Koumi S,Backer CL,Arentzen CE,Sato R. beta-Adrenergic modulation of the inwardly rectifying potassium channel in isolated human ventricular myocytes alteration in channel response to beta-adrenergic stimulation in failing human heartsThe Journal of Clinical InvestigationYear: 1995962870288110.1172/JCI1183588675658|
|41..||Jost N,Virág L,Opincariu M,Szécsi J,Varró A,Papp JG. Delayed rectifier potassium current in undiseased human ventricular myocytesCardiovascular ResearchYear: 19984050851510.1016/S0008-6363(98)00204-110070491|
|42..||Jost N,Virág L,Bitay M,Takács J,Lengyel C,Biliczki P,Nagy Z,Bogáts G,Lathrop DA,Papp JG,Varró A. Restricting excessive cardiac action potential and QT prolongation: A vital role for IKs in human ventricular muscleCirculationYear: 20051121392139910.1161/CIRCULATIONAHA.105.55011116129791|
|43..||Konarzewska H,Peeters GA,Sanguinetti MC. Repolarizing K+ currents in nonfailing human hearts similarities between right septal subendocardial and left subepicardial ventricular myocytesCirculationYear: 199592117911877648663|
|44..||Li GR,Feng J,Yue L,Carrier M,Nattel S. Evidence for two components of delayed rectifier K+ current in human ventricular myocytesCirculation ResearchYear: 1996786896968635226|
|45..||Virág L,Jost N,Opincariu M,Szolnoky J,Szécsi J,Bogáts G,Szenohradszky P,Varró A,Papp JG. The slow component of the delayed rectifier potassium current in undiseased human ventricular myocytesCardiovascular ResearchYear: 20014979079710.1016/S0008-6363(00)00306-011230978|
|46..||Koumi S,Backer CL,Arentzen CE. Characterization of inwardly rectifying K+ channel in human cardiac myocytes alterations in channel behavior in myocytes isolated from patients with idiopathic dilated cardiomyopathyCirculationYear: 1995921641747600647|
|47..||Maltsev VA,Sabbah HN,Higgins RS,Silverman N,Lesch M,Undrovinas AI. Novel ultraslow inactivating sodium current in human ventricular cardiomyocytesCirculationYear: 199898254525529843461|
|48..||Hoppe UC,Jansen E,Sudkamp M,Beuckelmann DJ. Hyperpolarization-activated inward current in ventricular myocytes from normal and failing human heartsCirculationYear: 199897613|
Literature data analysis results—full set of data
|Source||Patient-dependent data||Cell-dependent data|
|Location/health status||Number (all)||Number (male)||Number (female)||Average age [years]||SE [years]||Age min [years]||Age max [years]||Average Cm [pF]||SE [pF]||Temp [°C]||No.|
|||Ventricle||Left||Healthy + disease||6||5||1||56||6.3||350.4||29.0||37||28|
|||Ventricle||Left||Healthy + disease||52||49||2.4||114.8||17.9||22||30|
|||Ventricle||Left||Healthy + disease||13||45.8||4.8||245.0||17.0||22||57|
|||Ventricle||Left||Healthy + disease||47||223.7||7.3||22||133|
Dataset selected for the modeling step (healthy individuals)—after applying the inclusion/exclusion criteria
|Source||Patient-dependent data||Cell-dependent data|
|Location/health status||Number (all)||Number male||Number (female)||Average age [years]||SE [years]||Age min [years]||Age max [years]||Average Cm [pF]||SE [pF]||Temp [°C]||No.|
aIntentionally included in the final dataset; please see Materials and methods for further explanation
Dataset selected for the modeling step (heart disease individuals)—after applying the inclusion/exclusion criteria
|Source||Patient-dependent data||Cell-dependent data|
|Location/health status||Number (all)||Number (male)||Number (female)||Average age [years]||SE [years]||Age min [years]||Age max [years]||Cm [pF]||SE [pF]||Temp [°C]||No.|
Simulation results obtained with the use of Eq. 2 in combination with the age–cardiomyocyte volume/area model compared with the observed data
|Age [years]||SE [years]||Cm [pF]||SE [pF]||Fold error (% ≤2 fold)||MFE||RMSE|
|||av sim data||41.48||4.47||103.52||22.72||−1.99 (70)||1.93||93.2|
|||av sim data||44.76||3.96||100.20||20.87||−2.17 (50)||2.13||111.2|
|||av sim data||57.64||2.00||149.26||18.38||1.55 (100)||1.55||55.14|
|||av sim data||45.13||2.17||110.41||14.12||−1.64 (100)||1.63||69.39|
| epi||av sim data||39.25||2.57||81.41||13.36||−2.35 (20)||2.31||104.54|
| endo||av sim data||40.43||2.84||83.31||12.88||−2.12 (30)||2.10||91.31|
|||av sim data||64.73||2.61||205.13||28.22||−1.76 (80)||1.74||152.05|
|||av sim data||43.25||3.32||101.87||18.50||−1.84 (70)||1.81||81.85|
| epi||av sim data||75.10||1.22||295.86||37.93||1.10 (100)||1.29||81.39|
| endo||av sim data||75.35||1.04||331.53||43.89||2.16 (40)||2.12||188.26|
|||av sim data||38.33||1.46||66.83||5.97||−1.89 (70)||1.89||59.49|
|||av sim data||3.44||0.51||18.54||2.11||−5.08 (0)||5.04||74.11|
|||av sim data||3.27||0.55||18.10||2.28||−4.50 (0)||4.47||62.33|
Point and covariate matrix estimates of the log-linear model combining age and cardiomyocyte electric capacitance for the healthy and diseased population
|Population||Parameter||Point estimate||Covariance matrix estimate (×1,000)|
Point and covariate matrix estimates of the double-log-linear model combining age and cardiomyocyte electric capacitance for the healthy and diseased population
|Population||Parameter||Point estimate||Covariance matrix estimate (×1,000)|
Relationship between cardiomyocyte volume, electric capacitance, and APD90 for the log-linear model
|APD90 no drug||286||284||285||285||286||288||289||290||291||292||293||295||296|
Relationship between cardiomyocyte volume, electric capacitance, and APD90 for the double-log-linear model
|APD90 no drug||283.6||284.4||284.4||285.2||286.0||287.9||288.9||289.8||290.8||291.7||293.6||295.5||295.5|
Keywords: Keywords Inter-individual variability, Cardiomyocyte electric capacitance, Human cardiomyocyte, Virtual population.
Previous Document: Embryonic Stem Cells for Severe Heart Failure: Why and How?
Next Document: Contained versus uncontained lesions in juvenile elbow osteochondritis dissecans.