Development and evaluation of next-generation cardiotoxicity assay based on embryonic stem cell-derived cardiomyocytes.

In accordance with requirements of the ICH S7B safety pharmacology guidelines, numerous next-generation cardiotoxicity studies using human stem cell-derived cardiomyocytes (CMs) are being conducted globally. Although several stem cell-derived CMs are being developed for commercialization, there is insufficient research to verify if these CMs can replace animal experiments. In this study, in vitro high-efficiency CMs derived from human embryonic stem cells (hESC-CMs) were compared with Sprague-Dawley rats as in vivo experimental animals, and primary cultured in vitro rat-CMs for cardiotoxicity tests. In vivo rats were administrated with two consecutive injections of 100 mg/kg isoproterenol, 15 mg/kg doxorubicin, or 100 mg/kg nifedipine, while in vitro rat-CMs and hESC-CMs were treated with 5 μM isoproterenol, 5 μM doxorubicin, and 50 μM nifedipine. We have verified the equivalence of hESC-CMs assessments over various molecular biological markers, morphological analysis. Also, we have identified the advantages of hESC-CMs, which can distinguish between species variability, over electrophysiological analysis of ion channels against cardiac damage. Our findings demonstrate the possibility and advantage of high-efficiency hESC-CMs as next-generation cardiotoxicity assessment. [BMB Reports 2020; 53(8): 437-441].

INTRODUCTION

Cardiac toxicity is a major cause of failure in drug develop-ment (1). Also, drug-induced arrhythmia is the most common cause of restrictions or withdrawal of new drugs in the market. For these reasons, risk assessment for cardiovascular adverse effects is a crucial safety pharmacological element in new drug development. In the Regulatory Guides (2, 3), hERG (hKv11 potassium channel) assays and in vivo ECG assays are recom-mended to predict risk of arrhythmia associated with ventri-cular repolarization (QT prolongation) in humans.

Automated patch clamping is the standard method for high efficiency cardiac safety screening, since most clinical cases of drug-induced cardiac toxicity are associated with hERG inhibi-tion. However, this approach has limitations because securing sufficient primary human myocardial cells and long-term cultivation are difficult. Additionally, QT intervals are the result of interdependent ion channel activity, and inhibitory potential in the hERG does not always lead to arrhythmia for the heart (4). To address this problem, Purkinje fibers of guinea pig or rabbit heartshave been used for patch clamping at later stages of the drug screening process to study the action potential (AP) and ECG changes. However, interspecies variation can interfere with reliable connections to clinical relevance.

As an alternative to primaryhuman myocardial cells and animal experiments, human stem cell-derived cardiomyocyte model systems have been developed as physiologically relevant models for cardiac toxicity testing (5-9). Because of their ability to form synchronously beating cell populations, stem cell-derived cardiomyocytes are very attractive in vitro models that can be used to reliably reproduce drug-related cardiovascular phenotypes (10). Of particular significance is the ability to use stem cell-derived cardiomyocytes to screen for potential toxicities of chemicals that can cause cardiac arrhythmias as they express human ion channels similar to primary human cardiomyocytes and electrophysiological features similar to in vivo electrocardiography (11). Even in Comprehensive in vitro Proarrhythmia Assay (CiPA), developed as a next-generation cardiotoxicity model, cellular electrophysiological functions have been importantly evaluated upon drug effects (12). Although many laboratories, pharmaceutical companies, and research institutes have initiated human stem cell-derived cardiomyocytes research for cardiac toxicity assessments including CiPA (12-15), there is no preclinical study to evaluate if the patch clamping analysis of hERG channels using human embryonic stem cell-derived cardiomyocytes (hESC-CMs) predicts drug-toxic effects compared with standardization techniques such as in vivo animal ECG and hERG assays, which use non-cardiac cell lines overexpressed with the human ion channels (heterologous expression systems) or the disaggregated cardiomyocytes.

In this study, we compared and reported drug-induced cardiac toxicity to characteristic drugs using our high purity hESC-CMs reported in a previous study. For this purpose, we verified various gene panel expression changes and electro-physiological changes of hESC-CMs against animal cell and tissue data. Our study demonstrates that the hESC-CMs are an effective alternative testing method, which can replace the animal cell-based hERG assay or in vivo animal ECG model.

RESULTS

Characterization of hESC-CMs

In this study, hESC-CMs were generated using a small molecules-based differentiation method as previously reported (16-18) (Fig. 1A). Spontaneous beating was observed around differen-tiation days 10-15. Immunostaining for hESC-CMs was performed on day 30 by staining for cardiac specific markers, cardiac troponin T (cTnT) and sarcomeric α-actinin-2 (Fig. 1B) and approximately 80% of purity was confirmed by flow cytometry (Fig. 1C). On day 30, besides, the beating was recorded in video (Supplementary Video 1) and in microelectrode array (MEA) (Supplementary Fig. 1). Developed hESC-CMs were verified to meet the CiPA criteria (12).

Fig. 1

Overview of the hESC-CM differentiation. (A) The hESC-CM differentiation protocol using small-molecules. (B) At differentiation day 30, the cells were stained with cTnT and α-actinin-2. Scale bars: 100 µm. (C) And at day 30, the purity was monitored by the flow cytometry and showed as upper 80%.

Morphological changes in three models

With isoproterenol, doxorubicin, and nifedipine, morphological changes were observed (Fig. 2A). In isoproterenol and doxoru-bicin, rat hearts exhibited damaged cardiac myocyte architectures with cell shrinkage, nuclear pyknosis, karyolysis, densed cyto-plasm, and vacuolation. Coinciding with the histological damages, the cells with isoproterenol and doxorubicin became irregular and exhibited shrinkage and detachment of the cells from cell culture substratum, indicating apoptosis. Additionally, deterioration of α-actinin-2 arrangement was observed in all three models with isoproterenol or doxorubicin (Fig. 2B). Apoptotic indices were also calculated in rat hearts (Fig. 2C), rat-CMs (Fig. 2D), and the hESC-CMs (Fig. 2E) using the TUNEL assay (Supplementary Fig. 2). In the TUNEL results, isoproterenol and doxo-rubicin showed significant apoptotic cell increases in all models.

Fig. 2

Morphological changes in three models with cardiotoxicants. (A) The rat hearts showed damaged cardiac myocyte architectures (arrows), and characteristics of apoptotic cell death showed in the cells (arrowheads). Scale bars: 50 µm. (B) All of samples were stained with α-actinin-2 (green) and DAPI (blue). The arrows and arrow-heads indicate α-actinin-2 deterioration. Scale bars: 100 µm. (C-E) The apoptosis indices of rats (C), rat-CMs (D), and hESC-CMs (E) were assessed with the TUNEL. Data shown are mean ± SD (***P < 0.001; n = 5).

Relative changes of mRNA expression related to cytotoxicity

To assess the general cytotoxic response to drug-induced toxicity, the transcript levels of oxidative stress, cardiotoxicity, inflammation, and apoptosis-related genes were monitored (Fig. 3A). The oxidative stress indicators, iNOS, Hmox1, Sod1, and Gpx1 levels increased in all models with all drugs. Addi-tionally, the cardiotoxicity-related gene levels, β-MHC, ANP, BNP, and Spp1, and inflammation-related genes, IL-1β, IL-6, TLR2, TLR4, and TNF-α, were also increased in all groups. Also, indicators of apoptosis, the ratio of Bax/Bcl-2, and Casp9 and p53, increased significantly in all the groups.

Fig. 3

mRNA expression and electrophysiology in three models with cardiotoxicants. (A) The relative expression of iNOS, Hmox1, Sod1, Gpx1, β-MHC, ANP, BNP, Spp1, IL-1β, IL-6, TLR2, TLR4, TNF-α, ratio of Bax/Bcl-2, Casp9, p53 were monitored. Data shown are mean ± SEM (*, P < 0.05; **, P < 0.01; ***, P < 0.001; n = 5). (B-D) The heat-maps in which the entire set of genes were clus-tered hierarchically to visualize general expression patterns with isoproterenol (B), doxorubicin (C), or nifedipine (D). The profiles are represented by models (columns) and genes (rows). The color scale at the top represents relative expression level, wherein yellow and blue colors indicate upregulated and unaltered expression, respectively (See Supplementary Table 1). (E) The ECG and patch clamp were recorded (see Table 1).

Generally, nifedipine showed weak responses compared to the other drugs in all models. These results show the general cytotoxic response of hESC-CMs to cardiotoxicants is equivalent to rats and rat-CMs. Additionally, overall data was visualized as heat-maps to show general expression patterns of the three models by drugs (Fig. 3B-3D). Additionally, hESC-CMs showed relatively strong correlation with in vivo model with isopro-terenol (Fig. 3B), while rats and rat-CMs showed similarity with the doxorubicin (Fig. 3C) or the nifedipine (Fig. 3D). The gene order is presented in Supplementary Table 1.

Changes in electrophysiological profile

For comparison of quantitative analysis of the cardiac function, beat rate, corrected QT (QTc) intervals, and AP durations (APD) during drug treatment were evaluated respectively. The QTc and APD represent electrical depolarization and repolarization of the ventricles, prolongation of which is a marker of a potential risk factor for sudden death. The QTc was calculated using Bazett’s formula (QTcB), Fridericia’s formula (QTcF), and Sagie’s formula (QTlc) in all groups (19, 20). As shown in Fig. 3E and Table 1, rats and rat-CMs showed significant increases in the beat rate and prolongation of repolarization time with isoproterenol, a β-adrenoceptor agonist. hESC-CMs also showed a significant increase in the beat rates; however, APD of hESC-CMs showed the opposite response to the rat hearts and the rat-CMs. Importantly, this different response is associated with the expression of species-specific β-adrenoceptor subtypes. Next, there was no change in the beat rate with doxorubicin, an anticancer drug, which causes QTc prolonga-tion, but theduration of repolarization was prolonged in the rats. Finally, when treated with the L-type calcium blocker, nifedipine, which causes bradycardia and hypotension, significant reductions in the beat rate and repolarization times were observed in all groups however, no change in amplitude was observed in any group. These results indicate that hESC-CMs are equally applicable as laboratory animal rats and rat cardiomyocytes for conventional cardiac pharmacokinetic studies. Also, given the species different response, it shows the possibility of an optimal tool that can be more suitable for humans.

Table 1

Electrophysiological changes in three models with cardiotoxicants

	Ctrl	ISO	DOX	NIF
ECG results from rats
HR (bpm)	287±2	452±17***	271±16	224±22***
QT (msec)	87.76±12.92	91.84±0.52	130.61±10.43***	65.31±7.78**
QTcB (msec)	191.87±28.24	252.11±1.41***	277.52±22.16***	126.06±15.01***
QTcF (msec)	147.83±21.76	180.06±1.01*	215.87±17.24***	101.25±12.05***
QRlc (msec)	209.54±12.92	225.40±0.51*	250.50±10.43***	177.98±7.78***
Patch clamp results from rat-CMs
HR (bpm)	117±9	142±13**	114±10	38±10***
APD50 (msec)	64.6±5.2	78.4±6.4**	65.1±4.9	47.2±4.5***
Patch clamp results from hESC-CMs
HR (bpm)	83±12	107±9**	87±7	64±8*
APD50 (msec)	295.0±30.9	217.0±22.1**	291.0±21.4	157.0±38.0***

Data shown are mean ± SD (*P < 0.05; **P < 0.01; ***P < 0.001; n = 5). See Fig. 3E.

DISCUSSION

Globally, a variety of drug screening platforms are being studied, and human stem cells have been evaluated as the best alterna-tive to replace conventional animal testing. In particular, the ability to detect toxicity and efficacy of human species-specific responses, which can be missed in an animal model, is a major advantage of stem cells. We have reported a technique for high-efficient differentiation of hESC-CMs (21). In this study, we evaluated the technology of next-generation cardiotoxicity using hESC-CMs.

In this study, hESC-CMs showed very homogenous characteristics over a long period of culture, with the expression of cardiac-specific markers, likewise in the previous studies. The cultured hESC-CMs in sheet-shape additionally exhibited synchronized, uniform myocardial contraction, and stable characteristics in the electrophysiological analysis.

The rats, rat-CMs, and hESC-CMs were treated with reference drugs, isoproterenol, doxorubicin and nifedipine, which cause cardiac damage, respectively. For these cardiotoxicants, all three models showed generally similar responses in morphological and mRNA changes. While isoproterenol and doxorubicin caused significant damages to the tissues and cells, nifedipine showed weak responses, suggesting that nifedipine as a commercialized drug could be safe at concentrations below the lethal dose 50. These results demonstrate hESC-CMs can assess cardiotoxicity at comparable levels to conventional animal experiments or primary cultured cells. To prove further correlations between hESC-CMs with conventional models, heat-maps were assessed. Particularly, hESC-CMs showed strong correlation with in vivo model in isoproterenol than doxorubicin or nifedipine. Correlations between hESC-CMs and in vivo model with doxorubicin or nifedipine, however, seemed weaker than correlations between the rat-CMs and rats.

Different responses between the models were observed, which could be attributed to metabolism and cell composition differences in species. Susceptibilities for oxidative stress, particularly, were somewhat different between the samples originating from other species. For Sod1 and Gpx1 with nifedipine, rats and rat-CMs showed significant increases unlike hESC-CMs. Conversely, for iNOS with doxorubicin, hESC-CMs showed more dramatic change rather than the rat samples. These phenomena could be elaborated with facts that Sod1 and Gpx1 are directly impacted by ROS while iNOS is indirectly influenced via NF-κb. It could explain the different responding quickness for oxidative stress (22). Also, oxidative stress induces lipid peroxidation (23), producing products including malondialdehyde (24). Rats differ from humans within the membrane suscepti-bilities for these products (25, 26). Some differences between the rat samples and the hESC-CMs could be explained with facts about work-timing differences of enzymes in oxidative stress pathways and species differences of susceptibilities for oxidative stress products. Additionally, we observed critical differences by drugs in the electrophysiological results.

To verify the electrophysiological function of the hESC-CMs, we compared the three models. In the heart rate, increases by isoproterenol and decreases by nifedipine were significantly observed in all the models. Unlike morphological and mRNA analyses, doxorubicin did not impact electrophysiologically. Additionally, we analyzed cQT and APD, the main arrhythmia indicators, which correspond to myocardial depolarization and repolarization. cQT and APD reduction by nifedipine was observed in all the models. Importantly, hESC-CM APD was reduced with isoproterenol, unlike rat cQT and APD increase, suggested to be from species differences in β-adrenoceptor distribution. Ruzsnavszky, et al. (27) reported, in large mammals, β1-adrenoceptors are distributed more dominantly than β2-adrenoceptors and cAMP-dependent activities of Ca2+ and K+ currents are different from those of small mammals. These results suggest that the hESC-CMs are capable of assessing cardiotoxicity closer to humans compared to conventional models.

Although further assessment of human-specific toxicants should be accompanied, based on our developed the hESC-CMs, our team will continue to develop the next-generation cardiotoxicity assay through continuous verification.

CONCLUSION

The developed hESC-CMs were verified to have superiority than the conventional models in the toxicological reactions molecularly, morphologically, and electrophysiologically, while the hESC-CM can distinguish between species variabilities. These results demonstrate the possibility and advantage of the next-generation cardiotoxicity assessment using high-efficiency hESC-CMs.

MATERIALS AND METHODS

See supplementary information for Material and Methods.