1. Introduction
TRPM4 is a melastatin subfamily member of the transient receptor potential (TRP) superfamily and acts as a plasmalemmal route selective for Na
+ influx and K
+ efflux. In physiological settings, the TRPM4 channel is directly activated by intracellular Ca
2+ elevation upon transmembrane and intracellular Ca
2+ mobilizations. The distribution of this channel protein is ubiquitous across the whole body including both excitable (neurons, muscles) and non-excitable (secretory glands, blood cells, etc.) tissues. Thus, the TRPM4 channel is regarded as the most likely molecular identification of a broad class of Ca
2+-activated nonselective cation channels [
1]. Because of its predominant Na
+ permeability, major consequences of TRPM4 channel activation are thought to be two-fold, i.e., membrane depolarization and intracellular Na
+ loading. Generally, in non-excitable cells, the depolarizing effect of TRPM4 channel activation reduces non-voltage-gated Ca
2+ influxes via decreasing the Ca
2+ driving force. For example, in Jurkat T lymphocytes, the activation of TRPM4 channels was found to attenuate store-operated Ca
2+ influx, thereby inhibiting interleukin-2 synthesis [
2]. In excitable cells, membrane depolarization resulting from TRPM4 channel activation secondarily elicits action potentials (APs) via activation of voltage-dependent Na
+ and/or Ca
2+ channels to facilitate neurotransmission and cause muscle contraction. In the heart, TRPM4 has been suggested to play pleiotropic roles; while it acutely modulates inotropic, chronotropic and dromotropic properties of the heart in both positive and negative fashions [
3,
4,
5,
6,
7], it also chronically modifies remodeling processes such as physiological and pathological hypertrophies via regulation of non-voltage-gated Ca
2+ influxes [
8,
9].
In the past decades, genetic linkage analyses and subsequent cohort studies have identified dozens of
trpm4 gene mutations associated with conduction disorders such as progressive familial heart block type I, isolated cardiac conduction disorder, atrio-ventricular block, right-bundle branch block, and Brugada syndrome. However, somewhat contra-intuitively, in vitro assays show that many of these mutations exhibited increased rather than decreased TRPM4 channel expression/activity [
10]. For instance, the first identified TRPM4 mutation (E7K) was reported to increase cell-surface expression of TRPM4 channel protein (about twice) due to impaired SUMOylation without noticeable changes in its gating kinetics [
11]. This might cause the depolarizing shift of resting membrane potential which would in turn facilitate the inactivation of voltage-dependent Na
+ channel and decelerate AP propagation [
12], but there is little evidence to validate this possibility. Furthermore, in our numerical model simulations, simple doubling of the maximal wild-type TRPM4 channel activity produced only marginal changes in the resting membrane potential and the shape of AP [
13]. In addition, our recent study disclosed the functional abnormality of E7K mutation that unusually strengthens the interaction of TRPM4 channel activity with endogenous PIP
2 thereby increasing the risk of generating triggered activities [
14]. These facts raise the question of whether additional functional changes may also be involved in the pathogenesis of E7K-associated conduction disorders.
In the present study, to pursue this possibility, we adopted the following experimental approaches. First, we rigorously evaluated the gating kinetics of the ‘E7K’ mutant TRPM4 channel by use of ionomycin-perforated cell-attached (Iono-C/A) recording technique that allowed to stably record desensitization/rundown-prone TRPM4 channel activities [
13]. Kinetic data obtained from these experiments were then mathematically formulated as the rate constants of opening and closing, and they were incorporated into the most updated single-cell AP model reflecting the unique electrophysiology and intracellular Ca
2+ dynamics of human cardiac Purkinje fiber [
15]. Finally, this modified model was used to perform 1D-cable simulations to investigate the arrhythmogenic impact of the E7K mutation. The results indicate that the E7K mutation greatly increases the sensitivity of TRPM4 channels to voltage and intracellular Ca
2+ concentration ([Ca
2+]
i) to favor a longer sojourn in the open state and this property renders this mutant channel contributive to conduction failure.
3. Discussion
In the present study, we theoretically investigated the gating kinetics of a gain-of-function mutant of TRPM4 ‘E7K’ to show its pathophysiological significance in conduction disorder by means of gating analysis and numerical simulation. For this purpose, we adopted the Iono-C/A recording technique that allowed us to stably record TRPM4-mediated currents by minimizing rapid Ca
2+ desensitization and rundown [
13,
16]. The results of gating analysis indicated that both voltage and Ca
2+ dependencies of opening and closing rate constants are markedly affected by the E7K mutation so as to increase the open probability of TRPM4 channels in particular near the resting membrane potential and resting [Ca
2+]
i (
Figure 3). Furthermore, this theoretical prediction was confirmed by the power spectrum and single channel analyses, both of which showed prominent prolongation of open-life times of the E7K mutant channel (
Figure 4 and
Figure 5).
Incorporating the altered gating of the E7K mutant into a single-cell Purkinje fiber model demonstrated that increasing the channel density 1-to 5-fold produces density-dependent AP prolongation which ultimately converges to partially depolarized levels. In contrast, the same degrees of increases only slightly affect the shape and duration of AP as well as the resting membrane potential (RMP) in the wild-type TRPM4 channel (
Figure 6B). These disparate effects on AP and RMP are clearly reflected in AP conduction. While the speed of AP conduction along the 1D-cable stays almost constant in the wild-type, that of the E7K mutant progressively diminishes as its density (maximal activity) increases (
Figure 6C and
Figure 7A). The slowing of conduction (CV reduction) is well proportionate to the reduction in dV/dt
max (
Figure 6D and
Figure 7B), the measure for the magnitude of voltage-dependent Na channel (Na
v) current primarily contributing to AP upstroke [
19]. Although there is considerable complexity due to local loading effects and structural discontinuities at tissue levels [
20], this finding is consistent with the general idea that the velocity of AP propagation is correlated with the magnitude of a local Na
v current flow originating from the excited region. Moreover, the decreases in dV/dt
max and CV also coincide with the depolarizing shifts of membrane potential prior to AP upstroke, i.e., RMP. The level of RMP is crucial to determine Na
v availability just before AP generation, and in fact, the extent of the observed RMP shift reasonably accounts for the decrease in Na
v availability at AP upstroke which is estimated from its voltage-dependent inactivation curve (not shown). In aggregate, these results provide compelling evidence for the previous speculation [
12] that excessive E7K-TRPM4 activities at resting conditions would facilitate Na
v inactivation during diastole thereby slowing the generation and subsequent propagation of AP. However, it should be noted that the excessive activity of the E7K mutant enough large to produce conduction block would result from its greatly facilitated C-O gating (
Figure 5,
Figure 6 and
Figure 7) rather than its moderately increased expression [
11].
The channel density (maximal activity)-dependent CV reduction was reported by a previous simulation study which adopted a different mathematical formulation to describe the wild-type TRPM4 channel gating [
18]. However, there is a major difference between this and our studies in whether overexpression of wild-type TRPM4 can cause conduction block. While the former study demonstrated conduction failure with increased wild-type TRPM4 channel activity, our results show only negligible effects (
Figure 6C and
Figure 7A). There are at least three factors that could account for this discrepancy. First, Gaur et al. formulated voltage- and Ca
2+-dependent gating of the wild-type TRPM4 channel as independent processes [
18], but our model derived more complex formulations for TRPM4 channel gating by treating voltage and Ca
2+ dependencies as inter-dependent processes. Accordingly, the extent of wild-type TRPM4 channel activation in our model becomes much weaker near RMP and resting [Ca
2+]
i (which is, however, disrupted in the E7K mutant), as compared with the previous models which rather exaggerate the gating of TRPM4 or that of its native counterpart NSC
Ca around RMP [
13,
18]. Secondly, the Pan-Rudy model adopted by the previous study [
18,
21] appears more sensitive to increased TRPM4 channel activity than the human Purkinje fiber Trovato2020 model adopted in the present study [
15]. In our early simulations using the Pan-Rudy model, we noticed that increasing TRPM4 channel density destabilizes the model to generate frequent forward/backward spontaneous AP propagations along the cable. This made the exact evaluation of conduction velocity difficult. Thirdly, the maximal activity or density of TRPM4 channels defined in the present study may be low (it is set to be twice as large as that of atrial myocytes; see the methods), although we referred to the fact that expression of this channel is at least a few fold higher in Purkinje fiber than in the human atrium [
11]. However, doubling the maximal density merely shifts the CV density curves for E7K (
Figure 6C and
Figure 7A) to the left with least changes in that of the wild-type (
Figure S2), and this makes the impact of E7K mutation even more prominent. Therefore, even though the assumptions made for our present simulations are not perfectly realistic, the conclusions discussed above will still hold valid.
Introducing heterogeneity in TRPM4 channel density in the cable produces not only monotonic decrease in conduction velocity but also partial blocks that split into non-conductive and full- or sub-conductive states (
Figure 7B). This phenomenon was already reported by the preceding simulation work which observed more complex patterns of conduction failure classified as the first- to third-degree blocks [
18]. Indeed, in our simulations as well, various types of partial blocks are observed when the linear gradient exists in the channel density along the cable or fibroblasts are inter-placed in the cable (
Figure S3). The mechanism (s) underlying these complex phenomena remains entirely unclear but might involve the bifurcative properties of the model adopted and/or mutant TRPM4 channel gating per se [
22]. In real settings, spatial distribution of proteins would not be homogenous even within the same tissues [
23]. Thus, intercellular variations in transcriptional levels could introduce further complexities into conduction disorders in concert with inherited arrhythmogenicity.
In summary, the present study has disclosed a new pathogenic mechanism by which the E7K mutation induces brady-arrhythmogenicity, by use of the Iono-C/A recording-based gating analysis and numerical simulations with the most updated human Purkinje fiber model. The obtained results clearly show that at the single-cell level this mutation induces an excessive TRPM4 channel activity via acceleration and deceleration of its opening and closing transitions, respectively, but that the same alteration of gating simultaneously produces a variety of AP conduction blocks in cable-like multicellular arrangements. Further work will be needed to explore the exact clinical significance of these theoretical observations in future.
4. Materials and Methods
4.1. Cell Culture and Gene Transfection
Human embryonic kidney cells 293 (HEK293) were purchased from ATCC (Manassas, VA, USA) and maintained in Dulbecco’s modified Eagle medium supplemented with 10% fetal bovine albumin and a mixture of penicillin/streptomycin in a 100%-humidified, 5%CO2-gassed incubator at 35–36 °C, and were passaged every 3–4 days (up to 10–15 times). When reaching 70–90% confluency, HEK293 cells were dispersed by short trypsin treatment and gentle pipetting, and re-plated on cover slips pre-coated with poly L-lysin (30 μg/mL) for transfection. About 12 h later, the coverslips were incubated in a special medium (either DMEM or Opti-MEMTM, Gibco Life Technologies, Carlsbad, CA, USA) containing 1 μg pcI-neo vector encoding the CDS of htrpm4b gene or that of E7K mutant, with the aid of transfection agents superfect (Qiagen, Hilden, Germany) or lipofectatmine2000 TM (Invitrogen, Waltham, MA, USA) according to the manufactures’ instructions. Thirty-six to forty-eight hours after the transfection, electrophysiological recordings were carried out at room temperature. The human trpm4b cDNA (Gene ACC. No: AF497623) inserted in the pcDNA4TO-Flag vector was kindly provided by Profs. J.-P. Kinet (Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, MA, USA) and P. Launay (INSERM, Paris, France). For experimental use, this was subcloned into the pcI-neo vector.
4.2. Electrophysiology
The whole-cell and cell-attached (C/A) variants of the patch-clamp technique were applied. Patch electrodes were fabricated from 1.5 mm borosilicate glass capillaries (Sutter Instrument), and joined to the headstage of a low noise, high impedance patch clamp amplifier (EPC10, HEKA Elektronik, Ludwigshafen, Germany). When filled with the internal solution, the input resistance of the electrodes ranged between 4–7 MΩ. An automated multi-channel data acquisition software ‘Patchmaster’ (HEKA, Germany) was used to control a patch amplifier, and >60% of series resistance was electronically offset. For long-term recordings, current and voltage signals were sampled by the Power Lab data acquisition system (AD Instruments, Sydney, Australia), and obtained data were analyzed offline.
The details of ionomycin-perforated cell-attached (Iono-C/A) recording were described in our former study [
13]. Briefly, after ‘giga’ seal formation, the cell was quickly exposed to 5μM ionomycin-containing Ca
2+-free, high K
+ external solution, and then sequentially exposed to various extracellular concentrations ([Ca
2+]
o) (0.3, 1 and 5 mM). To null the resting membrane potential, ionomycin and Ca
2+ were administered in the presence of high K
+ throughout.
Single channel activities were recorded by the Iono-C/A method at 2 kHz digitization after 1 kHz low-pass filtering. The obtained data with no simultaneous multiple openings were re-filtered at 200 Hz (8-pole Bessel; rise time:1.66 ms) and subjected to dwell-time analysis. Dwell-time histograms were constructed after 2 ms binning and fitted by a sum of two exponentials, where events faster than 2 ms were ignored.
Power spectrum analysis was performed for the steady segments of macroscopic TRPM4 currents which were sampled at 1 kHz after 500 Hz low-pass filtering. The obtained data were re-filtered through 8-pole Butterworth at 200 Hz, from which power spectral density was calculated. The relationship between spectral density and frequency was then fitted to the equation: s(f) = s(0)/[1 + (f/fc)2], where s(f), s(0) and fc denote spectral density at a frequency of interest (in Hz), that of 0 Hz, and corner frequency, respectively.
For data analyses and illustrations, commercial software such as Clampfit v.10 (Axon Instruments, Foster City, CA, USA), Excel 2016 (Microsoft Office) and KaleidaGraph v.4 (Hulinks, Tokyo, Japan) were employed.
4.3. Solution
Standard external solution used for patch clamp experiments consisted of (in mM): 140 NaCl, 5 KCl, 1.2 MgCl
2, 1.8 CaCl
2, 10 Hepes, 10 glucose (adjusted to pH 7.4 with Tris base). High-K
+ external solution for the Iono-C/A recording (in mM): 145 KCl, 1.2 MgCl
2, 0, 0.3, 1 or 5 CaCl
2, 10 Hepes, 10 glucose (adjusted to pH 7.4 with Tris base); pipette solution for the Iono-C/A recording: the standard external solution with 1 mM tetraethylammonium (TEA) and 100μM 4,4’-Diisothiocyano-2,2’-stilbenedisulfonic acid (DIDS) to block K and Cl channels. Ca
2+ concentration of pipette solution was calculated by using a self-written program following Fabiato & Fabiato’s algorithm with the enthalpic and ionic strength corrections of association constants [
24]. All test solutions were applied via a hand-made local perfusion system controlled by electrically driven solenoid valves (time for completing solution change: ~1 s).
4.4. Numerical Model Simulation
For single-cell action potential (AP) simulation, the code of a most updated human Purkinje fiber model [
15] was downloaded from the CELLML repository [
25], and run by Cor1.1 or OpenCOR [
26]. The downloaded code was corrected for unit definition inconsistency. In order to incorporate the gating kinetics of TRPM4 channel and reproduce the time course of AP of the original Trovato 2020 model as exactly as possible, 40% of background Na conductance (I
Nab) in the original model was replaced by TRPM4 channel conductance [defined as the permeability (P
Na, P
K) of 7.02 × 10
−8 Litre/Farad/ms in Hodgkin-Huxley-Katz formalism] together with 5% increase and 5% decrease in the maximal conductances of I
Kr and I
to, respectively. The values of P
Na and P
K chosen for TRPM4 are twice as large as those electrophysiologically determined for HL-1 atrial cardiomyocyte model [
13,
27,
28], and this reflects the fact that expression level of TRPM4 protein is at least a few-fold higher in the Purkinje fiber than the atrium [
11]. A very high correlation (0.9999) was obtained by these alterations between the original and modified models (
Figure 6A). Corresponding changes in the other membrane conductance and intracellular Ca
2+ dynamics are illustrated in
Figure S4, where except for a newly added TRPM4 current and concomitant reduction in I
Nab, differences between before and after the modification are only marginal. All models used were stabilized by at least 1000 runs after any modifications. As for the other single-cell AP models, we adopted modified Luo-Rudy 2000 ventricular AP model and Pan-Rudy 2011 Purkinje fiber model [
21,
29]. The 4th-order Runge-Kutta algorithm was used to solve ordinary differential equations in these models. Action potentials and membrane currents were iterated at a 0.001 ms interval. For illustrative purposes, computed results were output in 1 ms resolution.
For 1D-cable simulations, the ‘Chaste’ simulation package [
30] was employed with a 3 cm-long cable discretized at every 0.01 cm (i.e., 301 nodes or 300 inter-nodal spaces; based on the data from [
31]). A multi-core (Intel
® Core™ i7-6800K CPU @ 3.40 GHz × 12, Intel Corporation, Santa Clara, CA, USA), 64-bit parallel computer (FRGBX911/A, Frontier, Tokyo, Japan; OS: Ubuntu 14.04 LTS, Canonical Ltd., London, UK) was used to speed up the simulation. Since there are considerable variations in reported values [
32,
33], intracellular conductivity was set to 1.7 mS/cm or its five times larger value 8.5 mS/cm with the extracellular/intracellular conductivity ratio of 1.0 [
32,
33]. It should be mentioned, however, that increasing this ratio or the value of intracellular conductivity several-fold, or adopting a longer cable (4.8 cm) with a twice larger inter-nodal distance of 0.02 cm, albeit small quantitative differences, did not essentially affect the conclusions from the simulations (not shown). To attain steady AP conduction, simulations were performed for 25 s with time steps of 0.01 or 0.005 ms and 0.01 or 0.05 ms for ordinary and partial differential equations, respectively (solved by forward Euler and finite element methods, respectively), and the results were output at 0.1 ms resolution as a binary file containing all voltage data at every node and timepoint. These data were used to calculate conduction velocity along the cable and maximal upstroke slope of propagated AP and visualize the spatiotemporal profile of AP propagation as surface plots by self-written programs in Python3 (Python Software Foundation, Beaverton, OR, USA) or Matlab2021a (Mathworks, Natick, MA, USA). The results of analysis were illustrated by KaleidaGraph v.4 (Hulinks, Tokyo, Japan) or Excel 2016 (Microsoft, Redmond, WA, USA).
4.5. Statistical Evaluation