Abstract
Free full text
Simultaneous electromechanical monitoring in engineered heart tissues using a mesoscale framework
Abstract
Engineered heart tissues (EHTs) generated from human induced pluripotent stem cell–derived cardiomyocytes (hiPSC-CMs) represent powerful platforms for human cardiac research, especially in drug testing and disease modeling. Here, we report a flexible, three-dimensional electronic framework that enables real-time, spatiotemporal analysis of electrophysiologic and mechanical signals in EHTs under physiological loading conditions for dynamic, noninvasive, longer-term assessments. These electromechanically monitored EHTs support multisite measurements throughout the tissue under baseline conditions and in response to stimuli. Demonstrations include uses in tracking physiological responses to pharmacologically active agents and in capturing electrophysiological characteristics of reentrant arrhythmias. This platform facilitates precise analysis of signal location and conduction velocity in human cardiomyocyte tissues, as the basis for a broad range of advanced cardiovascular studies.
INTRODUCTION
Human induced pluripotent stem cell–derived cardiomyocytes (hiPSC-CMs) have changed the landscape of cardiovascular research by providing access to a reliable source of human cardiomyocytes that can be reproducibly differentiated from hiPSCs by well-established, reproducible protocols (1, 2). Small-animal models of cardiovascular disease often fall short in recapitulating key aspects of human cardiac disease, partly mediated by considerable differences in heart rate, heart size, and protein isoform expression (3). For example, in the mouse α-myosin heavy chain is the predominant myosin heavy chain isoform in the ventricle, whereas β-myosin heavy chain in humans, is the predominant human ventricular isoform (4, 5). In addition, the high heart rate of small-animal models reflects physiologically critical differences in ion channel expression and Ca2+ handling. Larger-animal models can improve phenotyping but are impractical, cost-prohibitive, and unable to replicate many aspects of human myocardial disorders (6). Given the limitations of animal models in predicting the safety and efficacy of human response to drugs, the FDA Reauthorization Act 2.0 now allows for animal-free alternatives in developing new drugs (7, 8). hiPSC-CMs are well positioned to serve as a valuable model for both mechanistic and translational cardiovascular science including drug development, especially as they have a human genomic context (9). Successful translation requires, however, improved physiological readouts from hiPSC-CMs.
Most hiPSC-CM work has been conducted in a two-dimensional (2D) cell culture monolayer format. These models suffer from nonphysiologic culture conditions, immaturity, and a limited capacity to measure clinically relevant outputs such as contractility and arrhythmia propensity (9–11). Engineered heart tissues (EHTs) are created by embedding hiPSC-CMs with (12) or without fibroblasts or other cells (13) into a hydrogel matrix or decellularized tissue scaffolds (14). Such EHTs promote improved cardiomyocyte alignment and maturity. These systems allow for the application of load using flexible support posts, which is not possible using typical rigid 2D culture platforms (15). In this way, physiologic outputs such as contractile force and contraction velocity can be directly measured (12, 13). EHTs promote hiPSC-CM maturation with improved sarcomere alignment and protein localization and, by some reports, nascent T-tubule formation (15). While approaches for the dynamic monitoring of contractile force have been developed (16, 17), continuous, nonterminal integration with electrophysiologic monitoring remains a key challenge. Arrhythmia modeling in EHTs has largely been limited to terminal experiments using membrane voltage–sensitive dyes requiring toxic myosin inhibition with agents like blebbistatin (18). Multi-electrode arrays (MEAs) can support continuous field potential and impedance measurements, but such technologies are limited to the 2D cell culture format (19). Micropatterning, pacing, and application of mechanical stress can enhance the translational potential of 2D culture hiPSC-CMs (20), and advanced monitoring of spherical cardiac organoid formats has been reported (21). However, the resulting outputs are limited due to the aforementioned shortcomings of 2D culture methods and the lack of cellular alignment in spherical organoids relative to the 3D EHT format. Adaptation of conventional electronics for electrophysiologic interrogation, electrical stimulation, and mechanical monitoring has been hampered by the mechanical mismatch between rigid conventional sensors and soft, mechanically active tissues.
Through advances in materials science and engineering, soft and flexible electronics can now be created with a variety of sensing capabilities, including force, strain, temperature, pH, and electrical potential (22, 23) in both 2D and 3D architectures. In particular, beyond conventional 2D MEAs, which are limited to the bottom contact area of 3D organoids, geometries across the entire 3D surface offer great potential for organoid recording. For example, a shell-shaped MEA based on a bilayer, self-folding polymer for brain organoids provides a larger recording contact area compared to conventional MEAs (24). Another example involves a 3D MEA formed using a mechanical assembly process similar to the one used here, designed to gently wrap the entire spheroid surface with a collection of thin, flexible electronic and optoelectronic elements (23). In addition, self-rolling biosensor arrays have been used to support multiplexed monitoring with high sensitivity and spatiotemporal resolution for cardiac organoids (25). In addition, these platforms allow for electrical stimulation, which has clinical therapeutic relevance in pacemakers and defibrillators (26). This report leverages previously reported 3D systems, referred to as 3D multifunctional mesoscale frameworks (3D-MMFs), to create a 3D electromechanically monitored EHT platform that integrates electrodes and strain sensors distributed onto constructs mechanically tuned to support auxotonic EHT contraction. A range of experiments demonstrates reproducible, noninvasive, spatiotemporally monitoring during which field potential and contraction dynamics are correlated in EHTs over a period of weeks. The direct integration of flexible electronics with EHTs enables nonterminal electrophysiology (EP) studies in these tissues using stimulation techniques similar to those used on human hearts in the clinical EP laboratory. Here, we demonstrate an ability to provoke and sustain reentrant rhythms, highlighting the potential utility of this platform to model clinically relevant cardiac tachyarrhythmias.
RESULT
EHTs integrated with flexible 3D-MMFs
The device introduced here allows for real-time electrophysiologic sensing, pacing, and contractility assessments of EHTs under physiologic loading conditions. Figure 1A depicts the integration of a ~5.1-mm ring-shaped EHT onto a flexible 3D-MMF consisting of a multipost structure that includes six electrodes and six strain sensors distributed in a circular array, with defined physiologic auxotonic loading conditions. This system enables simultaneous acquisition of electrophysiological and mechanical behaviors at multiple, well-defined locations along circular-shaped EHTs, in real time. This 3D-MMF follows from geometrical transformation of a symmetrical 2D precursor with a custom design guided by finite element analysis (fig. S1), using adapted versions of concepts reported elsewhere (fig. S2) (23, 27, 28). Thin, mechanically compliant arc-shaped features define posts that include pairs of extended features and tissue supports to stabilize the mounted cardiac ring during culture. These structures ensure optimal, repeatable positioning of the ring with stable contact to each post and electrode (Fig. 1B) for measurements of local field potentials (Fig. 1C, left). Strain sensors, also integrated with each post, respond to bending induced by tissue contraction (Fig. 1C, right). The electrodes have dimensions of 0.82 × 0.22 mm and consist of electroplated nanocluster platinum (0.5 to 1.0 μm) as a uniform (fig. S4) and low-impedance interface [~4 kilohm at 1 kHz in a phosphate-buffered saline (PBS) solution (fig. S5)] for high-quality electrophysiological recordings (Fig. 1D, top). The overall framework relies on a multilayer of polyimide (PI; 5 μm)/Cr/Au (10/200 nm)/PI (5 μm) (fig. S3).
Mechanical loading is a normal part of myocardial physiology and, in excess, it drives pathophysiology, as in hypertensive cardiomyopathy or, in animal models, transverse aortic constriction to induce a pressure-overloaded state (29). Longitudinal tissue loading and force generation depends on the tissue and loading geometry. The posts have a low bending stiffness (k = 1.32 μN/μm; fig. S8) such that, using the geometric constraints particular to our platform and others, the loading conditions for this 3D-MMF (4.63 × 10−5 mN·mm/cell) are comparable to those of a standard commercial system (3.75 × 10−5 mN·mm/cell; 2-post, Myriamed; details of calculations are in Materials and Methods) and similar to those reported by others (30). Modified engineering designs can deterministically address a range of loading conditions relevant to study requirements.
The lower part of Fig. 1D shows a serpentine trace (width of 3 μm; initial resistance, 1.4 to 1.7 kilohm) as a strain gauge that exhibits a linear change in resistance with applied force [formula, 15.7x − 0.04 (fig. S8)]. This gauge integrates into a Wheatstone bridge with two identical fixed resistors. The output of this circuit, driven by a regulated dropout voltage (5 V), passes through an instrument amplifier (INA322, Texas Instruments) with a gain of 1000 (v/v) to yield a voltage difference (V) between the initial voltage and those induced by deformation of the framework. The signal then passes through a low-pass filter (fig. S6), to allow comprehensive real-time characterization of tissue contraction and relaxation dynamics under defined systolic (Fig. 1F) and diastolic (Fig. 1G) loading conditions.
EHTs form by embedding 750,000 cells in a collagen-based hydrogel matrix using 90% hiPSC-CMs differentiated through a standard ventricular differentiation protocol (2, 31) and 10% human cardiac fibroblasts (30). Tissues are initially matured in ring-shaped polydimethylsiloxane (PDMS) molds, followed by 7 days of remodeling (fig. S8) and then are manually mounted onto the 3D-MMF. These ring-shaped tissues have inner and outer diameters of ~4.4 and ~5.1 mm, respectively. The shape-matched 3D-MMF provides a mechanism for gentle tissue contact with the recording electrodes, secured with two tissue supports on each post of 3D-MMF (Fig. 1, B and E). Figure 1H shows a representative demonstration of time-synchronized field potential (red line) from one electrode and the signal (blue line) from the strain sensor on the same post. These synchronized data allow precise correlation of electrical and mechanical responses during different phases of the cardiac cycle.
Simultaneous spatiotemporal monitoring
The platform allows real-time electrophysiologic and contractility assessment under physiologic loading conditions. The time-synchronized electrical measurement system includes an amplifier (INTAN 128 ch RHS Stim/Recording controller with I/O expander; Intan Technology) capable of simultaneously capturing local tissue field potentials and tissue contractile dynamics from spontaneously contractive cardiac tissue. A separate optical recording setup precisely captures displacements of the EHTs through a microscope (Fig. 2A) for validating the responses of the strain gauges. Figure 2B shows the positioning of the six electrodes (E1 to E6), the six strain sensors (S1 to S6), and the reference electrode. Figure 2 (C and D) displays consistent field potentials and contractile waveforms. The signal to noise can be further improved by signal averaging a series of contractions (shown here for a 1-s interval, with ~30 cycles of activity). The potentials are ~±200 μV, and the corresponding forces have amplitudes of ~150 μN. Information on the conduction propagation and dynamic activity follow from recordings at sampling rates of 20 kHz from each of the six symmetrically positioned posts (Fig. 2, E and H). The time between electrical activation and peak systolic force generation (time to peak) is ~400 ms (Fig. 2, F and I), and the half-relaxation time is ~220 ms. The field potential measurements shown in Fig. 2F demonstrate bidirectional conduction propagation from a focus of automaticity that initiates near E5 and terminates near E2 ~150 ms later (Fig. 2G) with a conduction velocity of ~5.65 cm/s. The time delay of the propagation in mechanical systolic peak initiated by electrical activation is approximately 88 ms (Fig. 2J). Displacement profiles determined optically at D1 to D6 exhibit trends similar to those of the strain signals at corresponding locations of S1 to S6 (Fig. 2K). The results indicate that the maximum displacement during systolic movement is 80 μm, corresponding to 106 μN (Fig. 2L and fig. S8). During this period, the strain sensor captures a force of 105 μN, with an average match rate across all locations of 88 ± 2%. Unlike traditional optical voltage or calcium mapping approaches, which require inhibition of mechanical contraction, the EHT platform allows rapid, noninvasive measurements of both electrophysiologic and mechanical performance simultaneously across a single EHT over extended periods of observation.
Drug responses
A key aspect of hIPSC-CM derived EHTs is their utility for pharma-cological assessment. Using our electromechanically monitored EHT platform, we assessed the electrophysiologic and contractile effects of drugs used clinically to augment or inhibit cardiac rhythm and function. Figure 3 (A and B) shows a representative response to the β-agonist isoproterenol (3 μM), demonstrating an increase in beat rate [beats per minute (BPM)] and in contractile force (32), mirroring human cardiac responses to these agents. Figure 3 (C to E) highlight the dose response of isoproterenol with respect to BPM, QT interval, and contractile force (average value from 30 peaks, data from three tissues). The stepwise increase in beat rate and contractile force with a concomitant reduction in the QT interval are all consistent with known in vivo human cardiac response to isoproterenol (18, 33, 34). We also assessed the effect of sotalol, a commonly used antiarrhythmic drug with both potassium channel and β-adrenergic–blocking activity. As shown in Fig. 3 (F and G), the QT interval significantly prolonged after 10 μM sotalol treatment (30 overlaid plots of field potentials). Figure 3 (H to J) shows the dose response curve for sotalol for BPM, QT, and corrected QT (QTc). Increasing sotalol leads to a stepwise decrease in the beat rate with the expected increase in the QT interval (35, 36). We also examined the effect of omecamtiv mecarbil, a myosin inhibitor studied in the GALATIC-HF trial (37). Omecamtiv is known to prolong systole (38–40), with a response that depends on loading conditions, even in EHTs (32). The average tissue response (Fig. 3K) shows a delayed time to peak with increasing omecamtiv dose (Fig. 3L), consistent with prior observations in hiPSC-CM EHTs (32). Together, these findings demonstrate successful simultaneous electromechanical monitoring of hiPSC-CM EHT drug responses to clinically relevant adrenergic and sarcomeric agents.
Pacing and stimulation and termination of reentrant arrhythmias
The MMF platform also provides an opportunity to study the effects of tissue stimulation, both nonterminally and noninvasively, without the need for contamination-prone field effect electrodes or tissue-penetrating contact electrodes. Each electrode within the 3D-MMF can also serve as a pacing electrode, as shown in Fig. 4A. In this case, electrode 2 stimulates the EHT for 30 s at a frequency of 60 bpm (300 μs, 200 μA, biphasic pulse) while the electrical and mechanical response is monitored at electrode 5 and strain sensor 5. After fast pacing ceases, the slower automatic rhythm resumes as expected. Figure 4B displays a typical period of field potential and conduction propagation in a circular array. EHTs are known to occasionally exhibit spontaneously reentrant rhythms (18), which we also observed (Fig. 4C) with minimal mechanical deflection due to the dyssynchronous contraction. As in the clinical setting, reentrant rhythms are sustainable when the wavelength (refractory period × conduction velocity) is shorter than the path length, or alternatively said, an excitable gap is present between the wavefront front and wave tail.
Experiments that further leverage the pacing and the electrophysiologic and mechanical sensing capabilities of the EHT platform involve nonterminal studies to systematically induce and detect clinically relevant reentrant rhythms, without the need for membrane-sensitive dyes or toxic myosin inhibitors. The standard clinical protocol used on human hearts uses this same programmed electrical stimulation to induce reentrant rhythms (41). Figure 4D shows the induction of reentry in EHTs at a basic cycle length (BCL) of 600 ms (S1) with an S2 beat at 200 ms (Fig. 4G). Also shown in Fig. 4D, the reentrant rhythm terminates with rapid pacing. Figure 4E shows the activation map of the tissue under its spontaneous rhythm induced by an autonomic focus near E4. Conduction velocity, in this case, was ~8 cm/s. During reentry (Fig. 4F), the conduction velocity decreased markedly to ~4 cm/s, consistent with the expected conduction slowing at shorter cycle lengths. Reentry susceptibility increased with 25 μM pinacidil, which shortens the action potential duration. For some tissues, the excitable gap was too small with the path length close to the wavelength, such that a reentrant rhythm could not be induced.
DISCUSSION
Soft, biocompatible forms of electronics, especially those in 3D designs, that support sensors capable of monitoring biological processes have broad relevance in both clinical and research applications. This report substantially extends previously reported 3D-MMF technologies tailored to monitor electrical activity of neuronal organoids (23) and contraction of engineered skeletal muscle rings (42). The advances are in 3D framework designs tailored specifically to EHTs as both mechanical supports with well-defined loading conditions and as mounting sites for strain sensors and electrodes for both simulation and recording. The resulting technology enables continuous, simultaneous noninvasive monitoring of both field potentials and contraction dynamics at multiple locations in a precise, time-synchronized fashion. This EHT platform can noninvasively interrogate soft, mechanically active tissues, not only in circular ring geometries as demonstrated here but in other formats including strips or wedges. The noncompliant nature of conventional rigid electronics in standard 2D formats does not allow for adequate signal transduction in these contexts.
Over 100 years ago, Mines (43) first showed that reentrant rhythms could be developed in cardiac tissue and postulated a theoretical framework defining the requirements for initiation and maintenance of reentry. Garrey (44) also promoted the concept of a critical tissue mass for sustaining fibrillatory rhythms, highlighting the importance of tissue size or path length in reentry. Lewis later established that reentry was the mechanism of atrial fibrillation, which was further extended to ventricular fibrillation (45). Spontaneous reentrant rhythms have been observed in hiPSC-CM EHTs using optical mapping with voltage-sensitive membrane dyes in terminal experiments (18). While there are several alternatives for monitoring native action potentials with microelectrodes, voltage-sensitive membrane dyes, or genetically encoded transmembrane potential reporters, electromechanically monitored EHTs have the advantage of longer-term field potential monitoring, as a valuable proxy for clinically relevant parameters such as the repolarization interval (QTc) and conduction velocity. On the basis of our findings and the improved maturity of the electromechanically monitored EHT format, the EHT platform helps predict safety, arrhythmia susceptibility, response to drugs, and serves as a platform for testing gene-directed therapeutic strategies. Because EHTs are generated from human cells, EHTs can be used to advance personalized treatments like “n of 1” antisense treatment or gene editing (46).
Electromechanically monitored EHTs are also useful for measuring conduction velocities. When comparing the tissue conduction velocity to the native myocardium, the EHTs in this study are comparable to others (~6 cm/s), but do not reach the longitudinal conduction velocities measured in the normal ventricular myocardium (~50 cm/s). The EHTs examined here have a path length of 15.7 mm (diameter of 5 mm), requiring the refractory period to be <~260 ms to sustain reentry. This value is, however, an estimate, due to tissue-to-tissue variability. Also, both conduction velocity and action potential duration decrease with decreasing cycle length, leading to a concomitant shortening of wavelength with an increase in pacing rate. As EHT formation methods improve to enhance cardiomyocyte maturity, intercellular connectivity, and incorporate additional cell types found in the human heart, conduction velocity and, therefore, wavelength will likely increase. Previous work on human heart wedge preparations demonstrated a longitudinal wavelength of ~10 cm at a BCL 250 ms (240 bpm). The addition of 25 μM pinacidil reduces the wavelength to ~5 cm at a BCL = 250 ms (240 bpm) and ~3 cm at BCL = 125 ms (480 bpm) (47). The ability to simulate reentry could be enhanced by increasing the circumferential path lengths of the tissues and devices. For EHTs that approach native myocardium conduction velocities, this path length requirement would be >~10 cm, which is 6.3× the length of our system. Alternatively, regional manipulation of EHT conduction velocity to slow conduction in a single tissue region is a viable alternative. While technically demanding, this approach may provide a feasible solution as larger tissue sizes will necessarily require increased hiPSC-CM input and considerable additional cost.
By comparison to other work examining reentry in an EHT format (18), this format has several key advantages. First, the MMF system eliminates the need for specialized optical equipment to track movements and allowing for parallelization of measurements across multiple tissues. Second, signals are monitored noninvasively without the need to remove the tissue from the incubator environment, thereby facilitating time-dependent and longer-term experiments. Third, experimental induction of reentry can be accomplished using the same techniques applied to patients undergoing EP studies, namely, S1S2 stimulation, in a way that avoids the mechanical manipulation of the fragile EHT tissue required with point electrode stimulation. Fourth, the MMF platform provides a ready mechanism to document arrhythmias without the need for toxic dyes and myosin inhibitors, as the field potential can be monitored from electrodes simultaneously. Together, the MMF platform serves a critical niche for the study of arrhythmia biology in hiPSC-CMs for both basic and translational applications.
3D flexible electronic frameworks coupled with human hiPSC-CM–derived EHTs form the basis of this electromechanically monitored EHT system, enabling immediate spatiotemporal observation and manipulation through noninvasive, continuous electrophysiological and mechanical means. Experiments presented here demonstrate that electromechanically monitored EHTs provide a potential for efficient disease modeling and developing and evaluating candidate pharmacological compounds. These multisensor platforms, as applied to circular EHT rings, offer the capability to measure many sites simultaneously and instantly, providing immediate observation of signal propagation and signal intensity for drug reactions and other forms of stimulation. In particular, this system permits the study of cardiac tissue derived from a human genetic background, avoiding animal experiments. Also, each electrode can be used to deliver site-specific electrical stimulation to the tissue. EP studies on this MMF platform of ring-shaped EHTs demonstrate reentry rhythms, underscoring the ability to model human arrhythmias. This concept extends beyond the monitoring of biological signals under physiological loading conditions in EHTs to other forms of engineered soft tissues. These approaches thus have potential to improve the development of treatment options and drugs by integration with disease models.
MATERIALS AND METHODS
Fabrication of the 2D precursor
Techniques in photolithography and related methods in microfabrication yielded multilayered, flexible 2D precursors to the 3D-MMFs. The process began with spin coating (3000 rpm for 40 s) and baking (180°C for 150 s) of poly(methyl methacrylate) (PMMA) on a glass slide, followed by spin coating (1000 rpm for 40 s) and baking (120°C for 120 s, 150°C for 120 s, and 260°C for 1 hour with vacuum) of PI (PI-2545, HD Microsystems). A maskless aligner (MLA150, Heidelberg Instruments) defined patterns of a negative photoresist (AZ nLOF 2035) spin cast (5000 rpm for 40 s) and baked (110°C for 60 s). Postbaking (110°C for 90s) after photo exposure prepared the resist for development (immersion AZ MIF300 for 35 s) and rinsing with distilled water. Electron beam evaporation defined a thin, uniform layer of chrome (10 nm) and gold (200 nm), for lift off by immersion in a stripper (AZ400T, overnight). Spin coating and baking defined an encapsulating layer of PI on these patterned traces of Cr/Au. Patterning of a thin layer of copper (50 nm) deposited by sputtering on top of the PI involved first forming a layer of positive photoresist (S1805) by spin-casting (4000 rpm for 40 s) and baking (110°C for 60 s). Exposing the resist with a maskless aligner followed by immersion in a developer (AZ 917 for 30 s) exposed regions of the PI for removal by exposure to an oxygen plasma etching (220 mT, 200 W, 60 min). Immersion in a copper etchant diluted with distilled water (ratio of 10:1) for 2 s removed the copper hard mask. Rising in distilled water completed the fabrication.
Assembly of 3D MEA with electronics
Immersion in acetone for 30 min removed the sacrificial underlayer of PMMA, to allow release of the 2D precursor from the glass substrate. Transfer of the precursor onto a water-soluble polyvinyl alcohol (PVA) tape allowed exposure of the backside for patterning bonding sites by deposition through a shadow mask of PI. These sites consisted of Ti (10 nm)/SiO2 (50 nm) formed by electron beam evaporation through a shadow mask. Contact with a prestretched elastomer substrate (Dragon skin 10 medium, Smooth-On) after exposing both the elastomer and the backside of the 2D precursor to an ultraviolet plasma for 4 min activated the surface. Bringing these two elements into contact followed by heating to 70°C for 10 min formed strong covalent linkages at the bonding locations.
Dissolving the PVA tape in distilled water, and release of the prestretched elastomer, led to the formation of a 3D structure via mechanical buckling. A medical-grade silicone adhesive (Kiwk-sil, World Precision Instrument) bonded a well structure formed in PDMS to the elastomer substrate at distance slightly beyond the perimeter of the 3D-MMF. Interconnect contact pads served as an interface to a printed circuit board via anisotropic conductive film cables. A coating of Pt black on the electrodes resulted from electrodeposition controlled by a potentiostat (Autolab PGSTAT128N, Metrohm AG) in a three-electrode configuration (−0.1 V for 60 s) with a Pt wire as a counter electrode and Ag/AgCl as a reference electrode in a Pt solution (3 wt % of chloroplatinic acid and 0.1 wt % of lead acetate). The surface morphology of the Pt black was observed by field-emission scanning electron microscopy (S-3400N, Hitachi, Japan) (fig. S4). The impedance was measured in a PBS (pH 7.4) with a frequency range from 1 Hz to 1 mHz (fig. S5). All samples were sterilized using an autoclave at 124°C for 1 hour. Additional sterilization involved exposure to ultraviolet light for 30 min and to EtOH before mounting the EHTs on the device.
Determination of forces from strain gauge measurements
The output voltage of the Wheatstone bridge V0 follows
where VDD is the applied voltage, R0 is initial resistance of the strain gauge, R1,2 are fixed 1-megohm resistors, and R3 is a potentiometer for controlling bias voltage (fig. S6). As cardiac activation induces bending of a post, and its associated strain gauge, the resistance changes to
Quantitative relationships between the contraction force of the EHT and the change in output voltage follow from measurements of the gauge factors of the strain sensors using direct, nanoindentation (Hysitron Biosoft Indenter, Bruker, USA) tests. Measurement of the indentation force involved 10 repeated linear cycles with a cross-head speed of 20 μm/s and displacement ranging from 0 to 160 μm. Recording changes in the output voltage of the strain sensor established correlations to the deformations of the post. Experiments under uniaxial indentation indicate linear responses with gauge factors of 15.7x × 0.04 (fig. S8).
Culturing hiPSC-CMs, forming EHTs, and transfer to 3D-MMFs
Control hiPSCs GM03348 [healthy male, Corriel, previously reprogrammed (48)] and SCVI2372 (healthy male, Stanford SCVI Biobank) were cultured and differentiated into hiPSC-CMs using methods described elsewhere (49), which are derived from the protocols of Burridge et al. (50) and Buikema et al. (31). EHTs were created by combining 90% hiPSC-CMS with 10% human cardiac fibroblasts (Promocell, C-12375) and embedding in a collagen-based hydrogel per prior methods (30, 51) with the modification of an increase in the total number of cells per EHT to 750,000 from 500,000 (total volume of 180 μl). Tissues are cast on custom fabricated circular PMDS molds. Tissues are allowed to gel for 1 hour at 37 C, after which they were cultured in 10% horse serum (Gibco, 26050088), 1% penicillin/streptomycin, and 0.1% insulin in low-glucose Dulbecco’s minimum essential medium (Sigma-Aldrich, D5546) with media exchanged on a Monday-Wednesday-Friday schedule. For the first 3 days of culture, transforming growth factor–β1 (5 ng/ml; Peprotech 100-21C) was added to improve EHT remodeling and consolidation per prior reports (12, 30). After 1 week of culture, EHTs are washed with calcium- and magnesium-free PBS (Gibco, 14190144) and then treated with 30 μM N-benzyl-p-toluenesulfonamide (TCI, B3082) to inhibit tissue contraction for 5 min. The tissues were then gently lifted of the PDMS posts with the blunt side of forceps and transferred to the previously prepared flexible electronic. The unique layout matches requirements for use with the ring-shaped EHTs reported here, with inward slanting cantilever post design with an upper support structure to stabilize the tissue and ensure firm contact between the tissue and the electrodes. Media is exchanged for the EHT culture media and tissues were maintained for downstream studies.
Quantifying cell loading conditions
The cell loading conditions in length and geometry provided by a commercial tissue support structure with a traditional 2-post format are in fig. S11. The following describes a general approach to characterize auxotonic tissue loading conditions for 3D-MMFs with an arbitrary number of posts.
The following equation defines the auxotonic loading of the tissue
where sFtis is the circumferential force of the tissue normalized to cell/length or equivalently to cross-sectional area, K is the loading constant of the tissue, Ltis is the change in the circumferential length from the resting length, and Ltis is the circumferential resting (diastolic) length of the tissue
where Ftis is the circumferential tissue force, and CD is the cell density (cells per unit length). When K is equal between arbitrary tissue geometries, cross-sectional areas, and lengths, then the auxotonic tissue load is equivalent. Figure S11 shows free body diagrams for the two specific cases of interest, the 2-post commercial device and the 6-post 3D-MMF.
The equation for the 2-post case is summarized as follows
Summing the forces along the y axis, yields the following
Rearranging yields
The equation for the 6-post 3-MMF is summarized as follows
Where Δr is the diastole radius (rd) minus the systole radius (rs). The following relationship therefore holds
Then, combining from above
And therefore K for the 6-post 3D-MMF is the following
The contraction force information from the commercially available 2-post Myriamed system is as follows
The contraction force calculation of the 6-post 3D MMF system is as follows
Therefore, the loading conditions into the 3D MMF (K = 4.63 × 10−5 mN mm cells−1) are similar to the EHT created on the commercial Myriamed system (TM5: 3.75 × 10−5 mN mm cells−1).
Signal acquisition and data processing
A commercial amplifier system (INTAN 128 ch RHS Stim/Recording controller with I/O expander, Intan Technology) and its associated software (Intan RHX) served as the data acquisition platform for the time-synchronized recording of electrophysiological and mechanical responses of the EHT in a 37°C incubator (fig. S14). The RHS2116 amplifier chip allowed collection of electrophysiological data with 16-bit resolution at a sampling rate of 20 kHz by linking to a digital port. The strain gauges yielded changes in output voltage captured through an instrument amplifier (INA322, Texas Instrument) with a gain of 1000× from a quarter Wheatstone bridge circuit. A hardware low-pass filter was connected to analog ports for recording. In addition, real-time display capabilities allowed for immediate monitoring and data visualization. The raw data were processed using a digital band-stop filter constructed in MATLAB ranging from 59.5 to 60.5 Hz and a digital high-pass filter having cutoff frequency of 300 Hz. An inherent peak finding function in MATLAB defined electrophysiological peaks with minimum peak intervals of 500 ms and a threshold of 80% value from minimum to maximum for each peak. These peak indexes served as criteria for overlapping both electrophysical and mechanical signals from EHT in the MATLAB script, allowing for overlap of 30 peaks with minimal artifacts.
Optical video analysis
A custom optical tracking method quantified the deformation fields of EHTs and validated the correlated mechanical contraction signals from the 3D-MMF. The experiments used a camera (Motic SMZ-171 head and Motic Moticam 1080BHM, Motic, Canada) with a resolution of 1920 × 1080 pixels at 25 fps, mounted perpendicular to the sensor. The tracking method consisted of three steps. First, the image sequence was converted into a binary sequence using the Renyi Entropy method, where the pixel values representing the EHT were set to 255, while the background was set to 0. Second, the sectional area of the EHT with respect to angle (Δθ = 1 ± 3°) was tracked in the spatiotemporal domain (52). Last, the deformation fields were quantified by tracking the radial distance of the angular area over time. Movie S2 shows the overall process of the optical method.
Pharmacological responses on EHT
Pharmacologic studies were performed between week ~3 of EHT culture. In some cases, the same tissues were used for pharmacologic testing, in which case >48 hours of drug wash out was allowed. Stock solutions were prepared and diluted to the specified concentration in EHT media. Stock solutions were as follows: 10 mM isoproterenol (Sigma-Aldrich, I5627) in deionized water, 100 mM sotalol (Sigma-Aldrich, PHR2798) in dimethyl sulfoxide (DMSO), and 10 mM omecamtiv (Cayman Chemical, 21074) in DMSO. Dose response curves were obtained by incubating different concentrations of isoproterenol (0, 0.01, 0.1, 1, 3, and 10 μM), sotalol (0, 0.01, 0.1, 1, 10, and 100 μM), and omecamtiv (0, 0.1, 0.5, 1.0, 2.0, and 3.0 μM) in a stepwise fashion. Tissues were equilibrated for 20 min at each concentration in the incubator before measuring the tissue response.
Rapid pacing for EP study
For programmable S1S2 stimulation, a 2-ch oscilloscope/waveform generator (Analog Discovery 2, Digilent) was connected to the input terminal of the Intan Stimulating/Recorder, which was programmed to follow the output of the Analog Discovery 2 waveform generator that was programmed to produce a standard S1S2 stimulation protocol. The stimulation protocol involves eight cycles of S1 stimulation with an adjustable BCL and the S2 stimulation follows the last S1 stimulation in a stepwise manner with an adjustable S2 cycle length. The biphasic stimulus has a phase duration of 300 μs and an amplitude 400 μA. The cycle length of S2 is adjusted while recording the patterns of the six field potentials in real time.
Acknowledgments
D.E.F., I.R.E., and E.M.M. are CZ Biohub Chicago Investigators.Funding: This work is supported by the Querrey-Simpson Institute for Bioelectronics at Northwestern University and funding by NIH HL168758, NIH HL128075, NIH K08HL 163392, NIH HL167813, and Leducq Foundation. Y.P. acknowledges support by National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (RS-2023-00211261) and Korea Institute of Science and Technology (KIST) through 2E32960 and the Basic Science Research Capacity Enhancement Project through Korea Basic Science Institute (National Research Facilities and Equipment Center) grant funded by the Ministry of Education (grant no. 2019R1A6C1010052). S.O. acknowledges support by National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (2021R1C1C2010180). K.S.K and J.K acknowledge support by Ministry of Health & Welfare, Republic of Korea (HI22C0467) and National Research Foundation of Korea (NRF) grant funded by the Korea government (MIST) (RS-2024-00345402). This work made use of the Northwestern University Micro/nano Fabrication facility (NUFAB) of Northwestern University’s Atomic and Nanoscale Characterization Experimental (NUANCE) Center, which has received support from the SHyNE Resource (NSF ECCS-2025633), the INN, and Northwestern’s MRSEC program (NSF DMR-1720139).
Author contributions: D.E.F.: Conceptualization, methodology, writing–review and editing, and writing–original draft; W.-Y.M.: methodology, software, validation, data analysis, engineering investigation, visualization, data curation, writing–review and editing, and writing–original draft; S.O.: conceptualization, methodology, software, data analysis, engineering investigation, writing–review and editing; H.L.: theoretical simulation and software; K.S.K.: methodology, software, validation, data analysis, and writing–review and editing; I.A.C.: methodology and investigation; J.-T.K.: data analysis and investigation; J.-Y.Y.: data analysis and visualization; C.W.H.: methodology; A.R.D.: methodology; S.G.: methodology; B.L.: methodology; Y.H.: methodology; G.C.: visualization; J.Y.: visualization; J.K.: writing–review and editing; Y.P.: conceptualization, methodology, engineering investigation, writing–review and editing; I.R.E.: conceptualization and writing–review and editing; E.M.M.: conceptualization, project administration, supervision, writing–review and editing, and writing–original draft; J.A.R.: conceptualization, project administration, supervision, writing–review and editing, and writing–original draft.
Competing interests: The other authors declare that they have no competing interests.
Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.
Supplementary Materials
The PDF file includes:
Supplementary Text
Figs. S1 to S15
Legends for movies S1 to S4
Legends for data S1 to S5
References
Other Supplementary Material for this manuscript includes the following:
Movies S1 to S4
Data S1 to S5
REFERENCES AND NOTES
Articles from Science Advances are provided here courtesy of American Association for the Advancement of Science
Citations & impact
This article has not been cited yet.
Impact metrics
Alternative metrics
Discover the attention surrounding your research
https://www.altmetric.com/details/167234614
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
IGSR: The International Genome Sample Resource
- (1 citation) IGSR/1000 Genomes - GM03348
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Engineered heart tissue models from hiPSC-derived cardiomyocytes and cardiac ECM for disease modeling and drug testing applications.
Acta Biomater, 92:145-159, 07 May 2019
Cited by: 93 articles | PMID: 31075518
Afterload promotes maturation of human induced pluripotent stem cell derived cardiomyocytes in engineered heart tissues.
J Mol Cell Cardiol, 118:147-158, 28 Mar 2018
Cited by: 103 articles | PMID: 29604261 | PMCID: PMC5940558
Generation of a human iPSC-derived cardiomyocyte/fibroblast engineered heart tissue model.
F1000Res, 12:1224, 01 Jan 2023
Cited by: 0 articles | PMID: 38298530 | PMCID: PMC10828555
Engineered heart tissues and induced pluripotent stem cells: Macro- and microstructures for disease modeling, drug screening, and translational studies.
Adv Drug Deliv Rev, 96:234-244, 30 Sep 2015
Cited by: 78 articles | PMID: 26428619 | PMCID: PMC4698222
Review Free full text in Europe PMC