Abstract
An in vitro drug-induced cardiotoxicity assay is a critical step in drug discovery for clinical use. The use of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) is promising for this purpose. However, single hiPSC-CMs are limited in their ability to mimic native cardiac tissue structurally and functionally, and the generation of artificial cardiac tissue using hiPSC-CMs is an ongoing challenging. We therefore developed a new method of constructing three-dimensional (3D) artificial tissues in a short time by coating extracellular matrix (ECM) components on cell surfaces. We hypothesized that 3D cardiac tissues derived from hiPSC-CMs (3D-hiPSC-CT) could be used for an in vitro drug-induced cardiotoxicity assay. 3D-hiPSC-CT were generated by fibronectin and gelatin nanofilm coated single hiPSC-CMs. Histologically, 3D-hiPSC-CT exhibited a sarcomere structure in the myocytes and ECM proteins, such as fibronectin, collagen type I/III, and laminin. The administration of cytotoxic doxorubicin at 5.0 μM induced the release of lactate dehydrogenase, while that at 2.0 μM reduced the cell viability. E-4031, human ether-a-go-go related gene (hERG)-type potassium channel blocker, and isoproterenol induced significant changes both in the Ca transient parameters and contractile parameters in a dose-dependent manner. The 3D-hiPSC-CT exhibited doxorubicin-sensitive cytotoxicity and hERG channel blocker/isoproterenol-sensitive electrical activity in vitro, indicating its usefulness for drug-induced cardiotoxicity assays or drug screening systems for drug discovery.
Introduction
S
The use of stem cell technology is a promising approach for drug discovery, since human induced pluripotent stem cells (hiPSCs) enable us to overcome those species differences. Furthermore, the pluripotency of stem cells may have the potential for high-throughput screening for both safety and therapeutic efficacy of drugs in vitro.17–19 Indeed, hiPSCs established from patients with specific genetic abnormalities, such as long QT syndrome20–22 or congenital chondropathy, 23 were reported to be useful in screening candidate drugs to test their therapeutic efficacy.
Recently, a screening system for cardiotoxicity using cardiomyocytes derived from hiPSCs (hiPSC-CM) was explored. Although single-cell or monolayer cultures are widely used for such screening systems, some studies have suggested that cardiomyocytes cultured in a three-dimensional (3D) environment play a crucial role in pharmacologically induced responses.24–26 Various approaches to construct 3D artificial tissues have been examined to prove the safety of medications. Most of these studies have used specialized molds or scaffolds.26–29 Although specialized scaffolds make it possible to form 3D cardiac structures, those technologies are difficult to reproduce as a consequence of their complexity. The artificial scaffolds require consideration of the inhibition of cell–cell interaction or cell proliferation. Furthermore, in addition to the 3D structure, an extracellular matrix (ECM) is crucial to the function of the heart. 30 We therefore speculate that it may be necessary to develop appropriate cardiac tissues showing microstructures or pathophysiology similar to the native myocardium to evaluate cardiac toxicity accurately.
Our laboratory previously reported a new method of constructing 3D artificial tissues in a short time by coating ECM components onto cell surfaces. 31 This method enabled the construction of 3D structures without the inhibition of cell–cell interactions and culture with ECM components. In this study, we hypothesized that an in vitro drug-induced cardiotoxicity assay using 3D cardiac tissues derived from hiPSCs (3D-hiPSC-CT) may be useful as a new standardized cardiotoxicity test. We examined the mechanical and electrophysiological responses of 3D-hiPSC-CT against known cardiotoxic drugs that had different effects on cardiomyocytes, including doxorubicin to evaluate cytotoxicity, E-4031, hERG-type potassium channel blocker to evaluate the arrhythmic effect causing QT prolongation or torsades de pointes, and isoproterenol, a nonselective beta-adrenergic agonist to examine positive inotropy.
Materials and Methods
Cardiac differentiation of hiPSCs
hiPSCs (253G1; Riken, Ibaraki, Japan) were cultured and maintained in Primate ES Cell medium (ReproCELL, Kanagawa, Japan) with 4 ng/mL human basic fibroblast growth factor (bFGF; Wako, Osaka, Japan) on mouse embryonic fibroblast cells (ReproCELL). Cardiac differentiation was performed as follows: hiPSCs were dissociated using dissociation solution (ReproCELL) and transferred to an ultra-low attachment culture dish (Corning, MA) in mTeSR1 (Stemcell Technologies, Canada) with Y-27632 (Wako). After the formation of embryoid bodies, the culture medium was exchanged for a differentiation medium that contained StemPro34 (ThermoFisher Scientific, MA), 2 mM
Construction of 3D-hiPSC-CT
3D-hiPSC-CT were constructed as previously described. 32 In brief, hiPSC-CMs were dissociated with 0.05% trypsin-EDTA (ThermoFisher Scientific). The cells were suspended in 0.2 mg/mL fibronectin (Sigma-Aldrich) or gelatin (Wako) in phosphate-buffered saline (PBS; Nacalai Tesque) into a 24-mm transwell with 3.0 μm pore size (Corning) and incubated for 1 min at 500 rpm using a shaking incubator (SI-300; AS ONE, Osaka, Japan). The cells in the insert were then centrifuged for 1 min at 1000 rpm. After each coating step, the cells were washed with PBS in the same way. These coating steps were repeated nine times (fibronectin: five times and gelatin: four times), and fibronectin and gelatin nanofilms (FG) were coated over a single-cell surface. To construct the multilayer tissue (3D-hiPSC-CT), FG-coated cells were seeded at 5 × 105 cells into a 6.5-mm transwell with 0.4 μm pore size (Corning), and the culture medium was added to a 6-well plate outside of inserts or FG-coated cells were seeded onto a 96-well plate. The 3D-hiPSC-CT were incubated for 3–5 days, and the culture media were changed to DMEM containing 10% FBS every other day.
Quantitative real-time polymerase chain reaction
Total RNA was isolated from hiPSCs and hiPSC-CMs using a PureLink RNA Mini Kit (ThermoFisher Scientific). Quantitative real-time polymerase chain reaction (qPCR) was performed using a SYBR Green and ViiA 7™ Real-Time PCR System (ThermoFisher Scientific). All data were normalized using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as control and evaluated using the delta-delta CT method. Primer sequences are listed in Table 1.
αMHC, α myosin heavy chain; cTnT, cardiac troponin T.
Histology and immunohistolabeling
hiPSC-CMs and 3D-hiPSC-CT were fixed with 4% paraformaldehyde. hiPSC-CMs and 3D-hiPSC-CT were labeled by primary antibodies such as anti-cardiac troponin T (cTnT, 1:200 dilution; Abcam, Cambridge, UK), anti-sarcomeric alpha actinin (α-actinin, 1:400; Sigma), anti-vimentin (1:100; Dako, Glostrup, Denmark), anti-connexin43 (1:100; Abcam), anti-MLC2a (1:100; Synaptic Systems GmbH), anti-MLC2v (1:200; Proteintech, Rosemont) anti-fibronectin (1:200; Abcam), anti-laminin (1:30; Sigma), anti-collagen type I (1:100; Abcam), or anti-collagen type III (1:100; Abcam), followed by secondary antibodies such as AlexaFluor488 or AlexaFluor555 conjugated goat or donkey anti-mouse or anti-rabbit (ThermoFisher Scientific). Nuclei were counterstained with Hoechst33342 (Dojindo, Kumamoto, Japan) and assessed using confocal microscope (FLUOVIEW FV10i; Olympus, Tokyo, Japan).
Field potential measurements
The 3D-hiPSC-CT in the 6.5-mm transwell was replaced on MED probe (MED-P545A; Alpha MED Scientific, Osaka, Japan) and incubated until cell attachment. Extracellular field potentials were recorded using a multielectrode array (MEA) system (MED64; Alpha MED Scientific) for 30 min at 37°C.
Lactate dehydrogenase release assay
A lactate dehydrogenase (LDH) release assay was performed using an LDH Cytotoxicity Detection Kit (Takara Bio, Shiga, Japan). Doxorubicin (Wako) at 0.1, 1, 2, 5, or 10 μM was added for 24 h. After adding doxorubicin for 24 h, the supernatants were then collected, and a mixture of diaphorase/NAD+ was added to each well. After incubation in the dark room for 30 min at room temperature, the absorbance was recorded on a microplate reader (DS Pharma Biomedical, Osaka, Japan) at a wavelength of 490 nm with a reference wavelength of 600 nm. The experiments were performed in triplicate. The LDH release was determined as the percentage of LDH release compared with that of the vehicle control.
Cell viability assay
The cell viability was assessed using a Cell Counting Kit-8 (CCK8; Dojindo). Doxorubicin at 0.1, 1, 2, 5, or 10 μM was added for 24 h. After adding doxorubicin for 24 h, the media were then changed to fresh media, and 10 μL CCK8 solution was added to each well for incubation for 2 h at 37°C. Absorbance was recorded on a microplate reader at a wavelength of 450 nm with a reference wavelength of 630 nm. The experiments were performed in triplicate. Cell viability was determined as the percentage of surviving cells compared with that of the vehicle control.
Calcium transient measurement
Calcium transient was assessed by FDSS/μCELL (Hamamatsu Photonics, Shizuoka, Japan). The 3D-hiPSC-CT were incubated in a loading buffer containing 5 μM Cal-520™ AA (AAT Bioquest, CA) in PBS for 2 h at 37°C. E-4031 (Calbiochem Merck Millipore, Darmstadt, Germany) at 10, 30, or 100 nM or isoproterenol (Calbiochem Merck Millipore) at 10, 100, or 1000 nM was then added to the culture medium, and the fluorescence intensity of Cal-520 was monitored and recorded until 3 h after adding the drugs.
Cell motion analyses
Contractile property was assessed by a Cell Motion Imaging System (SI8000; SONY, Tokyo, Japan). Videos of the 3D-hiPSC-CT were recorded at a frame rate of 150 fps, a resolution of 1024 × 1024 pixels, and a depth of 8 bits. Subsequently, 10, 30, or 100 nM of E-4031 or 10, 100, or 1000 nM of isoproterenol was added to the culture medium, and the cell motion was monitored and recorded until 1 h after adding the drugs.
Statistical analyses
Data are expressed as the mean ± standard deviation. Statistical significance was determined by Student's t-test (two-tail) for comparisons between two groups. *p < 0.05 was considered statistically significant and ** indicates p < 0.01.
Results
Cardiomyogenic differentiation of hiPSCs
hiPSCs were cardiomyogenically differentiated by the well-established culture protocol (Fig. 1A). After differentiation, hiPSC-CMs showed significantly lower expressions of undifferentiated markers, such as Oct3/4, Nanog, and Lin28, and significantly higher expressions of cardiac markers, such as αMHC, cTnT, and Nkx2.5, than undifferentiated hiPSCs on qPCR (Fig. 1B). In addition, hiPSC-CMs expressed cardiac-specific structural proteins, such as cTnT, α-actinin, MLC2a, and MLC2v, on immunohistolabeling (Fig. 1C). hiPSC-CM included cardiac troponin-positive cells (57% ± 4.0%) and vimentin-positive cells (43% ± 4.0%), which are markers of fibroblast and smooth muscle cells, respectively (Fig. 1D).

Characterization of hiPSC-CMs.
Structural and electromechanical characters of 3D-hiPSC-CT
The structure of 3D-hiPSC-CT was characterized by immunohistolabeling for structural proteins of the myocytes and ECM proteins.
As a result, 3D-hiPSC-CT displayed thick multilayers, with a maximum thickness of 50 μm. Cardiac-specific structural proteins, such as cTnT, α-actinin, and connexin43, were clearly expressed, exhibiting a typical sarcomere structure and gap junctions in the majority of cardiac myocytes in the 3D-hiPSC-CT (Fig. 2A, B). In addition, the contractile protein isoform MLC2v was predominantly expressed in cardiac myocytes (Fig. 2C). Furthermore, cardiac-specific matrix proteins, such as fibronectin, laminin, and collagen, were homogeneously expressed in the extracellular space of 3D-hiPSC-CT (Fig. 2D–G).

Characterization of 3D-hiPSC-CT. Immunostaining of 3D-hiPSC-CT for
The electrical potential of the entire 3D-hiPSC-CT was assessed by the MEA system (Fig. 2H), showing spontaneous, synchronous, regular, and continuous beating of the tissue, indicating electrical linkage throughout the artificial tissue (Fig. 2I).
Drug-induced cytotoxicity of doxorubicin for 3D-hiPSC-CT
Doxorubicin-induced cytotoxicity of 3D-hiPSC-CT was assessed by an LDH assay and CCK8 assay. As a result, doxorubicin induced release of LDH in 3D-hiPSC-CT in a dose-dependent manner. Doxorubicin at 5 μM (145% ± 48%, p < 0.01 vs. vehicle control) and 10 μM (228% ± 97%, p < 0.01 vs. vehicle control) induced a significant increase in LDH release (Fig. 3A). Doxorubicin also reduced the cell viability in 3D-hiPSC-CT in a dose-dependent manner, and the concentrations that reduced cell viability significantly were 2 μM (88.4% ± 13%, p < 0.01 vs. vehicle control), 5 μM (75.3% ± 9.0%, p < 0.01 vs. vehicle control), and 10 μM (58.3% ± 10%, p < 0.01 vs. vehicle control) (Fig. 3B).

Drug-induced cytotoxicity of 3D-hiPSC-CT. The normalized change in
Ca2+ handling of 3D-hiPSC-CT in response to E-4031 and isoproterenol
Ca2+ transient of 3D-hiPSC-CT in response to E-4031 or isoproterenol was assessed by FDSS/μCELL (Fig. 4A).

Ca2+ Handling of 3D-hiPSC-CT by adding E-4031 for 10 min.
E-4031 at 30 and 100 nM for 10 min induced a significant decrease in the Ca2+ peak ratio (22.4 ± 3.4/min, p < 0.01, 4.8 ± 0.3/min, p < 0.01 vs. vehicle control [38.3 ± 2.6/min]) and a significant increase in the Ca2+ peak interval (2728.3 ± 416.6 ms, p < 0.01, 12,563.1 ± 805.9 ms, p < 0.01 vs. vehicle control [1570.2 ± 105.5 ms]) in 3D-hiPSC-CT (Fig. 4B–D). In addition, E-4031 at 10, 30, and 100 nM induced a significant decrease in the downstroke slope (11.2 ± 2.1, p < 0.01, 34.1 ± 0.7, p < 0.01, 1.0 ± 0.5, p < 0.01 vs. vehicle control [16.2 ± 1.5]) and a significant increase in the 90% peak width duration (PWD90) (867.3 ± 37.0 ms, p < 0.01, 1944.3 ± 300.2 ms, p < 0.01, 11,866.6 ± 942.0 ms, p < 0.01 vs. vehicle control [741.7 ± 13.5 ms]) (Fig. 4E, F).
In contrast, isoproterenol at 10, 100, and 1000 nM for 5 min induced a significant increase in the Ca2+ peak ratio (55.0 ± 4.4/min, p < 0.01, 54.5 ± 1.9/min, p < 0.01, 51.6 ± 4.6/min, p < 0.01 vs. vehicle control [35.1 ± 0.9/min]) and in the upstroke slope (47.6 ± 6.2, p < 0.05, 46.7 ± 4.4, p < 0.05, 66.1 ± 10.4, p < 0.01 vs. vehicle control [38.6 ± 3.4]), and a significant decrease in the Ca2+ peak interval (1097.0 ± 85.1 ms, p < 0.01, 1102.8 ± 40.1 ms, p < 0.01, 1168.6 ± 106.5 ms, p < 0.01 vs. vehicle control [1710.9 ± 42.4 ms]) and in the PWD90 (671.9 ± 28.7 ms, p < 0.01, 647.0 ± 21.8 ms, p < 0.01, 660.4 ± 16.4 ms, p < 0.01 vs. vehicle control [805.3 ± 14.6 ms]) (Fig. 5A–E).

Ca2+ Handling of 3D-hiPSC-CT by adding isoproterenol for 5 min.
Contractile property of the 3D-iPSC-CT in response to E-4031 and isoproterenol
Change in contractile property of the 3D-iPSC-CT in response to E-4031 and isoproterenol was assessed by the Cell Motion Imaging System (Fig. 6A).

Contractile property of 3D-hiPSC-CT in response to E-4031 for 10 min.
E-4031 at 30 and 100 nM for 10 min induced a significant decrease in the beating rate (32.7 ± 1.2/min, p < 0.01, 29.6 ± 1.7/min, p < 0.01 vs. vehicle control [37.4 ± 1.5/min]) and a significant increase in the contraction–relaxation duration (CRD) (0.81 ± 0.03 s, p < 0.01, 0.98 ± 0.04 s, p < 0.01 vs. vehicle control [0.71 ± 0.01 s]). In addition, the contraction–relaxation velocities were significantly decreased at 10, 30, and 100 nM of E-4031 for 10 min (contraction velocity; 9.5 ± 1.3 μm/s, p < 0.01, 9.8 ± 2.0 μm/s, p < 0.05, 9.8 ± 1.5 μm/s, p < 0.05 vs. vehicle control [14.7 ± 2.0 μm/s], relaxation velocity; 6.5 ± 1.3 μm/s, p < 0.05, 5.8 ± 1.4 μm/s, p < 0.01, 6.4 ± 0.6 μm/s, p < 0.01 vs. vehicle control [9.8 ± 1.1 μm/s]) (Fig. 6B–E).
In contrast, isoproterenol at 10, 100, and 1000 nM for 10 min induced a significant increase in the beating rate (49.7 ± 1.5/min, p < 0.01, 58.3 ± 1.3/min, p < 0.01, 59.2 ± 1.37/min, p < 0.01 vs. vehicle control [35.0 ± 3.2/min]) and a significant decrease in CRD (0.54 ± 0.008 s, p < 0.01, 0.52 ± 0.01 s, p < 0.01, 0.52 ± 0.02 s, p < 0.01 vs. vehicle control [0.61 ± 0.01 s]). In addition, contraction–relaxation velocities were significantly increased at 10 and 100 nM (contraction velocity; 28.6 ± 1.0 μm/s, p < 0.01, 25.4 ± 1.8 μm/s, p < 0.01 vs. vehicle control [18.2 ± 2.4 μm/s], relaxation velocity; 15.7 ± 2.4 μm/s, p < 0.05, 16.5 ± 2.8 μm/s, p < 0.05 vs. vehicle control [11.5 ± 1.6 μm/s]) (Fig. 7A–D).

Contractile property of 3D-hiPSC-CT in response to isoproterenol for 10 min. Beating rate
Discussion
In this study, we found that the 3D-hiPSC-CT generated from single hiPSC-derived cardiomyocytes coated with fibronectin and gelatin displayed multilayered cells exhibiting cardiac-specific structural proteins and matrix proteins in the extracellular space. The 3D-hiPSC-CT exhibited a cytotoxic response against doxorubicin. Furthermore, 3D-hiPSC-CT exhibited a decrease in the Ca2+ peak ratio and the downstroke slope and an increase in the Ca2+ peak interval and PWD90 in response to E-4031, and an increase in the Ca2+ peak ratio and the upstroke slope and a decrease in the Ca2+ peak interval and PWD90 in response to isoproterenol. The 3D-hiPSC-CT also exhibited a decrease in the beating rate and contraction–relaxation velocities in response to E-4031 and an increase in the beating rate and contraction–relaxation velocities in response to isoproterenol. Given these findings, 3D-hiPSC-CT appear to have potential utility in the evaluation of drug-induced toxicity and pharmacological action in vitro.
The 3D-hiPSC-CT displayed completely synchronized mechanical and electrophysiological properties. Links between the layered cardiac myocytes by the connexin43 gap junction, which was clearly shown immunohistologically, would be a major positive contributor to enhanced properties as cardiac tissue for 3D-hiPSC-CT. Some claim that layered artificial cardiac tissue induces ischemic changes in the centered cardiac myocytes, producing mechanical and electrophysiological heterogeneity without a microvascular network. However, the 3D-hiPSC-CT with a 50 μm thickness did not show heterogeneity in the electrical propagation, suggesting that addition of a microvascular network would not be necessarily required to preserve the functionality of ≤50-μm-thick 3D-hiPSC-CT.
The minimum effective concentration of doxorubicin for 24-h cytotoxicity of 3D-hiPSC-CT was 5.0 μM in the LDH assay and 2.0 μM in the CCK8 assay in this study. Since the reported minimum effective concentration of doxorubicin in single or monolayer hiPSC-derived cardiomyocytes is not substantially different from that obtained in this study, the 3D-hiPSC-CT may carry a comparative tolerance against doxorubicin to single cardiomyocytes or the monolayered cardiomyocytes.33–36 Since the cytotoxic effects of doxorubicin directly target DNA degradation and early apoptosis, these results indicate the homogeneous delivery of doxorubicin into the cardiomyocytes of entire 3D-hiPSC-CT, suggesting the utility of 3D-hiPSC-CT for drug screening tests.
Our study showed that the minimum effective concentration of E-4031 on the Ca2+ transient and contraction–relaxation performance of 3D-hiPSC-CT was 10 nM for 10 min, findings comparable to the reported concentration of E-4031 (10 nM) in the examination using single dissociated cells. 37 Some articles have reported that the minimum effective concentration of E-4031 ranges from 3 to 100 nM and these reported data suggest that the effects of E-4031 on the Ca2+ transient and contraction–relaxation performance in 3D-hiPSC-CT might be comparative to that of standard assays, such as hERG, patch clamp, and MEA.38–44 Our results showed greater drug sensitivity than animal studies, indicating an adequate, at least not inferior, effect of E-4031 on the 3D-hiPSC-CT compared to the standard methods. 45
In addition, the Ca2+ transient levels in 3D-hiPSC-CT changed with the administration of isoproterenol at 10 nM for 5 min. Furthermore, the contraction–relaxation performance of 3D-hiPSC-CT increased with the administration of isoproterenol at 10 nM for 10 min. The effective concentration of isoproterenol has been reported to range from 10 to 1000 nM. Our results therefore achieved equivalent or superior detecting potential for isoproterenol compared with the standard methods or adult human heart slices.46–48
Since most of articles monitored drug responses for 10–60 min after adding E-4031 or isoproterenol, monitoring time in this study is enough. However, it is necessary to consider the monitoring time for evaluating chronic drug responses.
We should also evaluate the appropriateness of this screening system for cardiotoxicity. To validate this system, we should compare the in vitro evaluation using 3D cardiac tissue with an in vivo screening system using a small animal or Langendorff experiments, or compare ours with other in vitro screening systems.
Fibronectin-/gelatin-coated cells were used to generate the 3D-hiPSC-CT, which achieved a uniform artificial cardiac tissue, contributing to the standardization of the protocol of 3D-hiPSC-CT production, as reported by our laboratory for the production of skin, liver, or blood vessel models of 3D artificial tissue.49,50 Although 3D-hiPSC-CT with a 50-μm thickness was used in this study, layered 3D-hiPSC-CT with a developed vascular network might more accurately mimic the microstructure and functionality of human cardiac tissue. Our protocol makes it easy to construct thick tissues without artificial scaffolds by increasing the number of coated cells to be seeded. Moreover, the development of 3D printing techniques to fabricate 3D tissues automatically is now in progress. 51 This protocol therefore has the potential to be scaled and applied to various experimental settings. This study may be limited by the use of a single cell line of iPSCs that were cardiomyogenically differentiated in a single protocol. hiPSC-CM included about 60% cardiomyocytes in this study. Although the differentiation method needs to be improved, our laboratory previously reported that engineered cardiac tissues containing 50–70% cardiomyocytes exhibited stable structures. 52 The optimum ratio of cardiomyocytes in the 3D models should be considered. Drugs act on cardiomyocytes directly or indirectly through surrounding cells. To more closely mimic the drug response in the in vivo environment, it might be useful to construct 3D-hiPSC-CT mixed with fibroblast or endothelial cells and verify their optimum ratio. Furthermore, 3D models generated using various numbers of cells with different ratios may be used to evaluate the effects of drugs not only on cardiotoxicity but also on angiogenesis or fibrosis.
In conclusion, the 3D-hiPSC-CT detected doxorubicin-sensitive cytotoxicity and hERG channel blocker/isoproterenol-sensitive electrical/mechanical activities in vitro, indicating its usefulness for drug-induced cardiotoxicity assays or in drug screening systems for drug discovery.
Footnotes
Acknowledgments
We acknowledge Tsuyoshi Ishikawa and Nobu Miyakawa for expert technical assistance. This work was supported by the SENTAN-JST Program, Research Project for Practical Application of Regenerative Medicine and Health Labour Sciences Research Grant.
Disclosure Statement
No competing financial interests exist.
