Next Article in Journal
The Changes in the p53 Protein across the Animal Kingdom Point to Its Involvement in Longevity
Next Article in Special Issue
Modulation of the Cardiac Myocyte Action Potential by the Magnesium-Sensitive TRPM6 and TRPM7-like Current
Previous Article in Journal
Time-Dependent Unitary Transformation Method in the Strong-Field-Ionization Regime with the Kramers-Henneberger Picture
Previous Article in Special Issue
TRPM8 Channels: Advances in Structural Studies and Pharmacological Modulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block

1
Department of Physiology, Fukuoka University School of Medicine, Fukuoka 814-0180, Japan
2
Biomedical Information Engineering Lab, Division of Information Systems, the University of Aizu, Aizu-Wakamatsu 965-8580, Japan
*
Authors to whom correspondence should be addressed.
Both authors contributed equally.
Int. J. Mol. Sci. 2021, 22(16), 8513; https://doi.org/10.3390/ijms22168513
Submission received: 25 June 2021 / Revised: 26 July 2021 / Accepted: 4 August 2021 / Published: 7 August 2021
(This article belongs to the Special Issue TRPM Channels)

Abstract

:
In the heart, TRPM4 is most abundantly distributed in the conduction system. Previously, a single mutation, ‘E7K’, was identified in its distal N-terminus to cause conduction disorder because of enhanced cell-surface expression. It remains, however, unclear how this expression increase leads to conduction failure rather than abnormally enhanced cardiac excitability. To address this issue theoretically, we mathematically formulated the gating kinetics of the E7K-mutant TRPM4 channel by a combined use of voltage jump analysis and ionomycin-perforated cell-attached recording technique and incorporated the resultant rate constants of opening and closing into a human Purkinje fiber single-cell action potential (AP) model (Trovato model) to perform 1D-cable simulations. The results from TRPM4 expressing HEK293 cells showed that as compared with the wild-type, the open state is much preferred in the E7K mutant with increased voltage-and Ca2+-sensitivities. These theoretical predictions were confirmed by power spectrum and single channel analyses of expressed wild-type and E7K-mutant TRPM4 channels. In our modified Trovato model, the facilitated opening of the E7K mutant channel markedly prolonged AP duration with concomitant depolarizing shifts of the resting membrane potential in a manner dependent on the channel density (or maximal activity). This was, however, little evident in the wild-type TRPM4 channel. Moreover, 1D-cable simulations with the modified Trovato model revealed that increasing the density of E7K (but not of wild-type) TRPM4 channels progressively reduced AP conduction velocity eventually culminating in complete conduction block. These results clearly suggest the brady-arrhythmogenicity of the E7K mutant channel which likely results from its pathologically enhanced activity.

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 Ca2+ elevation upon transmembrane and intracellular Ca2+ 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 Ca2+-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 Ca2+ influxes via decreasing the Ca2+ driving force. For example, in Jurkat T lymphocytes, the activation of TRPM4 channels was found to attenuate store-operated Ca2+ 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 Ca2+ 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 Ca2+ 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 PIP2 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 Ca2+ 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 Ca2+ concentration ([Ca2+]i) to favor a longer sojourn in the open state and this property renders this mutant channel contributive to conduction failure.

2. Results

2.1. Gating Analysis Reveals That E7K Mutation Facilitates TRPM4 Channel Opening

With normal recording variants of patch clamp technique, the activity of expressed TRPM4 declined quickly because of rapid Ca2+-dependent desensitization and rundown [13,16]. To circumvent this, we developed a new recording technique dubbed ‘ionomycin-perforated cell-attached (Iono-C/A)’ recording [13]. By using this recording procedure, it was possible to stably record TRPM4-mediated currents up to about 10 min upon repeated exposure to different concentrations of Ca2+ in the bath.
We applied this procedure to evaluate the gating behavior of the TRPM4 channel carrying the ‘E7K’ mutation by means of voltage-jump protocols (Figure 1). In response to step depolarizations, TRPM4-mediated currents (induced by Ca2+) showed time-dependent increases in the amplitude. On the contrary, repolarizations back to more negative potentials deactivated the currents (see the tail currents just after the end of depolarizing pulses in Figure 1A). Figure 1B,C, respectively, show the averaged values of variables for steady-state activation [open probability (Po)] and activation/deactivation time courses (time constant (τ) after normalization with respect to three different concentrations of Ca2+ in the bath ([Ca2+]o; note that these are later converted to [Ca2+]i values; for details, see Figure 1 legend). Similar to the wild-type TRPM4 channel [13,17], the relationships of Po and τ against the membrane potential (Vm) for E7K show clear voltage dependencies, and these are modified by Ca2+. Boltzmann fitting of the Po-Vm relationships for the E7K mutant demonstrate that increasing [Ca2+]o shifts the voltage-dependent activation curve toward more negative Vm. Furthermore, the values of the half-activation voltage (V0.5), which is an indication for the dynamic Vm range of a channel being open, are 47.5, 48.8 and 3.7 mV for [Ca2+]o of 0.3,1 and 5 mM, respectively (these values are 54.7, 80.4 and 38.9 mV for the wild-type TRPM4 [13]). This suggests that the voltage dependency of the E7K mutant may be more susceptible to Ca2+ than the wild-type TRPM4 channel [13].
To more unequivocally elucidate the altered gating by the ‘E7K’ mutation, we calculated the rate constants of opening (α) and closing (β) from the voltage and Ca2+ dependencies of Po-Vm and τ-Vm relationships (Figure 1B,C), based on the two-state C-O transition scheme (Figure 2). From these, the final mathematical expressions for α and β were obtained as the functions of Vm and [Ca2+]i:
α ( V , [ C a ] ) = 4.767 · [ C a ] 0.42534 · exp [ ( 0.012121 + 0.000000034609 · [ C a ] ) · V ]
( in s 1 ; [ Ca ] in μ M , V in mV )
β ( V , [ C a ] ) = 0.37975 · exp [ ( 3.046 + 0.07528 · [ C a ] ) + ( 0.0022098 0.0010278 · [ C a ] ) · V + ( 0.000011 0.0000118 · [ C a ] ) · V 2 ]
( in s 1 ; [ Ca ] in μ M , V in mV )
Reconstructed relationships of Po, α and β versus Vm and [Ca2+]i indicate that the E7K mutation renders the channel more readily open at more negative membrane potentials and at lower [Ca2+]i by enhancing the voltage and Ca2+ dependencies of α and decreasing those of β. In other words, this mutation accelerates the closed-to-open state (C-O) transition and simultaneously decelerates the reverse (O-C) transition of TRPM4 channels via facilitated voltage-dependent gating and increased Ca2+-sensitivity (Figure 3).
To confirm this prediction by patch clamp experiments, we next performed the power spectrum analysis of whole-cell TRPM4 currents induced by direct Ca2+ infusion via a sharp electrode [13] and the dwell-time analysis of single TRPM4 channel activities. As demonstrated and summarized in Figure 4, the open life-time of the E7K-mutant was significantly longer than that of the wild-type in single-channel analysis, which appears to match up with the estimation from the power spectrum analysis. In addition, the relationship between the intracellularly perfused [Ca2+]i and the functional density of the resultant whole-cell E7K-TRPM4 current was significantly leftward-shifted with most robust activation near the resting [Ca2+]i level (0.28 μM), as compared with that of the wild-type TRPM4 current (Figure 5). These results are consistent with the theoretical prediction of enhanced Ca2+ susceptibility by the E7K mutation.

2.2. Numerical Simulations Show That Preferred Opening by E7K Mutation Causes Conduction Block

Intuitively, enhanced voltage-and Ca2+-dependency, in particular in their low physiological ranges, may imply more excitatory influences of the E7K mutation on both AP generation and conduction of the heart. In fact, computer simulation of single-cell AP shows that regardless of numerical models employed, the facilitated opening of the E7K-mutant TRPM4 channel much more readily induces AP prolongation as its cell-membrane density increases, compared with the wild-type TRPM4 channel (Figure 6B and Figure S1). Furthermore, when the prolongation is prominent, small oscillatory depolarizations which resemble early after-depolarizations (EADs) appear on the repolarizing phase of AP (Figure 6B).
However, when a 1D-cable model was employed to compute AP conduction, different consequences of the mutation were revealed. Figure 6C and Figure 7A recapitulate the influence of increased maximal channel activity on the conduction velocity (CV) along a 3 cm-long cable. Because TRPM4 is reported to be expressed most abundantly in the conduction system, we employed the most updated human single-cell AP model for Purkinje fiber [15]. The results of simulation clearly indicate that while CV is only slightly affected by 1- to 5-fold increases in wild-type TRPM4 channel density, the same degrees of the density increases cause progressive reductions in CV for the E7K mutant eventually culminating in complete conduction block. The extent of CV reduction appears strengthened by a faster pacing (2 Hz) but little dependent on the intracellular conductivity. The reduction in CV in the E7K mutant channel is proportionate to the decrease in the maximal slope of propagated AP upstroke (dV/dtmax) and appears correlated with the depolarizing shift of the resting membrane potential (RMP) (Figure 6C,D and Figure 7A,B). Interestingly, introducing heterogeneity in the channel density along the cable (see also in [18]) created more complex patterns of conduction block such as 2:1 or 3:1 conductions which were split into non-conductive and full- or sub-conductive states (Figure 7C). These results strongly suggest that the impact of facilitated opening by the E7K mutation would manifest differentially at the single-cell and structurally higher levels.

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 Ca2+ desensitization and rundown [13,16]. The results of gating analysis indicated that both voltage and Ca2+ 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 [Ca2+]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/dtmax (Figure 6D and Figure 7B), the measure for the magnitude of voltage-dependent Na channel (Nav) 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 Nav current flow originating from the excited region. Moreover, the decreases in dV/dtmax 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 Nav availability just before AP generation, and in fact, the extent of the observed RMP shift reasonably accounts for the decrease in Nav 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 Nav 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 Ca2+-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 Ca2+ dependencies as inter-dependent processes. Accordingly, the extent of wild-type TRPM4 channel activation in our model becomes much weaker near RMP and resting [Ca2+]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 NSCCa 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 Ca2+-free, high K+ external solution, and then sequentially exposed to various extracellular concentrations ([Ca2+]o) (0.3, 1 and 5 mM). To null the resting membrane potential, ionomycin and Ca2+ 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 MgCl2, 1.8 CaCl2, 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 MgCl2, 0, 0.3, 1 or 5 CaCl2, 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. Ca2+ 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 (INab) in the original model was replaced by TRPM4 channel conductance [defined as the permeability (PNa, PK) 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 IKr and Ito, respectively. The values of PNa and PK 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 Ca2+ dynamics are illustrated in Figure S4, where except for a newly added TRPM4 current and concomitant reduction in INab, 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

All experimental data in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 are expressed as the mean ± s.e.m. Statistical analysis was performed by the Student t-test.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22168513/s1.

Author Contributions

Conceptualization, Y.H. and R.I.; Methodology, Y.H. and R.I.; Software, R.I., Y.H., Q.L., Y.S. and X.Z.; Validation, Y.H. and R.I.; Formal Analysis, Y.H., Q.L. and R.I.; Investigation, Y.H., Q.L. and R.I.; Resources, Y.H., X.Z. and R.I.; Data Curation, Y.H. and R.I.; Writing—Original Draft Preparation, Y.H. and R.I.; Writing—Review and Editing, T.F. and R.I.; Visualization, Y.H. and R.I.; Supervision, R.I.; Project Administration, R.I.; Funding Acquisition, Y.H. and R.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by JSPS grants-in-aid to R.I. (No.15H04678; No.20K07269) and Y.H. (No.20K16126).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data described in this paper are available on request.

Acknowledgments

We cordially thank Mizue Akaboshi (Department of Physiology, Fukuoka University School of Medicine) for her excellent secretarial assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nilius, B.; Szallasi, A. Transient receptor potential channels as drug targets: From the science of basic research to the art of medicine. Pharmacol. Rev. 2014, 66, 676–814. [Google Scholar] [CrossRef]
  2. Launay, P.; Cheng, H.; Srivatsan, S.; Penner, R.; Fleig, A.; Kinet, J.-P. TRPM4 regulates calcium oscillations after T cell activation. Science 2004, 306, 1374–1377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Feng, J.; Zong, P.; Yan, J.; Yue, Z.; Li, X.; Smith, C.; Ai, X.; Yue, L. Upregulation of transient receptor potential melastatin 4 (TRPM4) in ventricular fibroblasts from heart failure patients. Pflügers Arch.-Eur. J. Physiol. 2021, 473, 521–531. [Google Scholar] [CrossRef]
  4. Hof, T.; Sallé, L.; Coulbault, L.; Richer, R.; Alexandre, J.; Rouet, R.; Manrique, A.; Guinamard, R. TRPM4 non-selective cation channels influence action potentials in rabbit Purkinje fibres. J. Physiol. 2016, 594, 295–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hof, T.; Simard, C.; Rouet, R.; Sallé, L.; Guinamard, R. Implication of the TRPM4 nonselective cation channel in mammalian sinus rhythm. Heart Rhythm 2013, 10, 1683–1689. [Google Scholar] [CrossRef]
  6. Mathar, I.; Kecskes, M.; Van der Mieren, G.; Jacobs, G.; Camacho Londoño, J.E.; Uhl, S.; Flockerzi, V.; Voets, T.; Freichel, M.; Nilius, B. Increased β-adrenergic inotropy in ventricular myocardium from Trpm4−/− mice. Circ. Res. 2014, 114, 283–294. [Google Scholar] [CrossRef] [Green Version]
  7. Pironet, A.; Syam, N.; Vandewiele, F.; Van den Haute, C.; Kerselaers, S.; Pinto, S.; Vande Velde, G.; Gijsbers, R.; Vennekens, R. AAV9-mediated overexpression of TRPM4 increases the incidence of stress-induced ventricular arrhythmias in mice. Front. Physiol. 2019, 10, 802. [Google Scholar] [CrossRef]
  8. Gueffier, M.; Zintz, J.; Lambert, K.; Finan, A.; Aimond, F.; Chakouri, N.; Hédon, C.; Granier, M.; Launay, P.; Thireau, J. The TRPM4 channel is functionally important for the beneficial cardiac remodeling induced by endurance training. J. Muscle Res. Cell Motil. 2017, 38, 3–16. [Google Scholar] [CrossRef] [PubMed]
  9. Kecskés, M.; Jacobs, G.; Kerselaers, S.; Syam, N.; Menigoz, A.; Vangheluwe, P.; Freichel, M.; Flockerzi, V.; Voets, T.; Vennekens, R. The Ca2+-activated cation channel TRPM4 is a negative regulator of angiotensin II-induced cardiac hypertrophy. Basic Res. Cardiol. 2015, 110, 1–14. [Google Scholar] [CrossRef] [Green Version]
  10. Vennekens, R. Recent insights on the role of TRP channels in cardiac muscle. Curr. Opin. Physiol. 2018, 1, 172–184. [Google Scholar] [CrossRef]
  11. Kruse, M.; Schulze-Bahr, E.; Corfield, V.; Beckmann, A.; Stallmeyer, B.; Kurtbay, G.; Ohmert, I.; Schulze-Bahr, E.; Brink, P.; Pongs, O. Impaired endocytosis of the ion channel TRPM4 is associated with human progressive familial heart block type I. J. Clin. Investig. 2009, 119, 2737–2744. [Google Scholar] [CrossRef]
  12. Rezazadeh, S.; Duff, H.J. Genetic determinants of hereditary bradyarrhythmias: A contemporary review of a diverse group of disorders. Can. J. Cardiol. 2017, 33, 758–767. [Google Scholar] [CrossRef]
  13. Hu, Y.; Duan, Y.; Takeuchi, A.; Hai-Kurahara, L.; Ichikawa, J.; Hiraishi, K.; Numata, T.; Ohara, H.; Iribe, G.; Nakaya, M. Uncovering the arrhythmogenic potential of TRPM4 activation in atrial-derived HL-1 cells using novel recording and numerical approaches. Cardiovasc. Res. 2017, 113, 1243–1255. [Google Scholar] [CrossRef] [Green Version]
  14. Hu, Y.; Li, Q.; Kurahara, L.-H.; Shioi, N.; Hiraishi, K.; Fujita, T.; Zhu, X.; Inoue, R. An Arrhythmic Mutation E7K Facilitates TRPM4 Channel Activation via Enhanced PIP2 Interaction. Cells 2021, 10, 983. [Google Scholar] [CrossRef] [PubMed]
  15. Trovato, C.; Passini, E.; Nagy, N.; Varró, A.; Abi-Gerges, N.; Severi, S.; Rodriguez, B. Human Purkinje in silico model enables mechanistic investigations into automaticity and pro-arrhythmic abnormalities. J. Mol. Cell. Cardiol. 2020, 142, 24–38. [Google Scholar] [CrossRef]
  16. Nilius, B.; Mahieu, F.; Prenen, J.; Janssens, A.; Owsianik, G.; Vennekens, R.; Voets, T. The Ca2+-activated cation channel TRPM4 is regulated by phosphatidylinositol 4, 5-biphosphate. EMBO J. 2006, 25, 467–478. [Google Scholar] [CrossRef] [Green Version]
  17. Nilius, B.; Talavera, K.; Owsianik, G.; Prenen, J.; Droogmans, G.; Voets, T. Gating of TRP channels: A voltage connection? J. Physiol. 2005, 567, 35–44. [Google Scholar] [CrossRef] [PubMed]
  18. Gaur, N.; Hof, T.; Haissaguerre, M.; Vigmond, E.J. Propagation failure by TRPM4 overexpression. Biophys. J. 2019, 116, 469–476. [Google Scholar] [CrossRef] [Green Version]
  19. Weidmann, S. The effect of the cardiac membrane potential on the rapid availability of the sodium-carrying system. J. Physiol. 1955, 127, 213–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kléber, A.G. The shape of the electrical action-potential upstroke: A new aspect from optical measurements on the surface of the heart. Circ. Res. 2005, 97, 204–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Li, P.; Rudy, Y. A model of canine purkinje cell electrophysiology and Ca2+ cycling: Rate dependence, triggered activity, and comparison to ventricular myocytes. Circ. Res. 2011, 109, 71–79. [Google Scholar] [CrossRef] [Green Version]
  22. Prescott, S.A.; De Koninck, Y.; Sejnowski, T.J. Biophysical basis for three distinct dynamical mechanisms of action potential initiation. PLoS Comput. Biol. 2008, 4, e1000198. [Google Scholar] [CrossRef] [Green Version]
  23. Skelly, D.A.; Squiers, G.T.; McLellan, M.A.; Bolisetty, M.T.; Robson, P.; Rosenthal, N.A.; Pinto, A.R. Single-cell transcriptional profiling reveals cellular diversity and intercommunication in the mouse heart. Cell Rep. 2018, 22, 600–610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Inoue, R.; Ito, Y. Intracellular ATP slows time-dependent decline of muscarinic cation current in guinea pig ileal smooth muscle. Am. J. Physiol.-Cell Physiol. 2000, 279, C1307–C1318. [Google Scholar] [CrossRef]
  25. A Computational Human Cardiac Purkinje Electrophysiological Model. Available online: https://models.physiomeproject.org/e/5f0/ModTrovato2020.cellml/view (accessed on 10 May 2021).
  26. Available online: http://cor.physiol.ox.ac.uk/ (accessed on 29 May 2021).
  27. Hu, Y.; Kaschitza, D.R.; Essers, M.; Arullampalam, P.; Fujita, T.; Abriel, H.; Inoue, R. Pathological activation of CaMKII induces arrhythmogenicity through TRPM4 overactivation. Pflügers Arch.-Eur. J. Physiol. 2021, 473, 507–519. [Google Scholar] [CrossRef] [PubMed]
  28. Takeuchi, A.; Kim, B.; Matsuoka, S. The mitochondrial Na+-Ca2+ exchanger, NCLX, regulates automaticity of HL-1 cardiomyocytes. Sci. Rep. 2013, 3, 1–11. [Google Scholar] [CrossRef] [Green Version]
  29. Luo, C.-H.; Rudy, Y. A dynamic model of the cardiac ventricular action potential. II. Afterdepolarizations, triggered activity, and potentiation. Circ. Res. 1994, 74, 1097–1113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Chaste. Available online: http://www.cs.ox.ac.uk/chaste/ (accessed on 10 May 2021).
  31. Callewaert, G.; Carmeliet, E.; Vereecke, J. Single cardiac Purkinje cells: General electrophysiology and voltage-clamp analysis of the pace-maker current. J. Physiol. 1984, 349, 643–661. [Google Scholar] [CrossRef] [PubMed]
  32. Clayton, R.; Bernus, O.; Cherry, E.; Dierckx, H.; Fenton, F.H.; Mirabella, L.; Panfilov, A.V.; Sachse, F.B.; Seemann, G.; Zhang, H. Models of cardiac tissue electrophysiology: Progress, challenges and open questions. Prog. Biophys. Mol. Biol. 2011, 104, 22–48. [Google Scholar] [CrossRef]
  33. Stinstra, J.G.; Hopenfeld, B.; MacLeod, R.S. On the passive cardiac conductivity. Ann. Biomed. Eng. 2005, 33, 1743–1751. [Google Scholar] [CrossRef]
Figure 1. Evaluation of E7K-TRPM4 channel gating by voltage-jump experiments. (A) actual traces of E7K-mutant TRPM4 channel currents responding to a set of voltage pulses from−100 to 60 mV (left row) and 80 to 200 mV (right row) by 20 mV steps. In each row, top panels show voltage pulses, while middle and bottom panels show corresponding TRPM4-mediated currents at 0.3 and 5 mM Ca2+ in the bath. The TRPM4-mediated currents were recorded by Ion-C/A recording where 5 μM ionomycin-containing Ca2+-free solution was applied onto the cell immediately after giga-seal formation to permeabilize the cell membrane, which was then switched to solutions with various concentrations of extracellular Ca2+ ([Ca2+]o) to activate the current. All these agents were applied through a fast solution device ‘Y-tube’. It should be noted that the level of a rise in intracellular Ca2+ concentration ([Ca2+]i) corresponding to a given [Ca2+]o, which was separately determined by Ca2+ imaging experiments, was almost stable and did not involve Ca2+ release from internal stores [13]. (B,C) the relationship of open probability (Po) versus membrane potential (B) and that of time constant (τ) of activation/deactivation time courses during voltage pulses (C). Different symbols and bars in the figure denote mean ±S.E.M. from five independent experiments. Smooth curves in A are drawn according to the best fits of data points (normalized to the maximum) by the Boltzmann equation: 1/[1 + exp(−(Vm − V0.5)/s), where Vm, V0.5 and s denote membrane potential, half-activation voltage and slope factor, respectively. The values of these parameters are 47.542, 48.825 and 3.697 mV, and 56.582, 67.327 and 52.844 mV for 0.3, 1 and 5 mM [Ca2+]o, respectively.
Figure 1. Evaluation of E7K-TRPM4 channel gating by voltage-jump experiments. (A) actual traces of E7K-mutant TRPM4 channel currents responding to a set of voltage pulses from−100 to 60 mV (left row) and 80 to 200 mV (right row) by 20 mV steps. In each row, top panels show voltage pulses, while middle and bottom panels show corresponding TRPM4-mediated currents at 0.3 and 5 mM Ca2+ in the bath. The TRPM4-mediated currents were recorded by Ion-C/A recording where 5 μM ionomycin-containing Ca2+-free solution was applied onto the cell immediately after giga-seal formation to permeabilize the cell membrane, which was then switched to solutions with various concentrations of extracellular Ca2+ ([Ca2+]o) to activate the current. All these agents were applied through a fast solution device ‘Y-tube’. It should be noted that the level of a rise in intracellular Ca2+ concentration ([Ca2+]i) corresponding to a given [Ca2+]o, which was separately determined by Ca2+ imaging experiments, was almost stable and did not involve Ca2+ release from internal stores [13]. (B,C) the relationship of open probability (Po) versus membrane potential (B) and that of time constant (τ) of activation/deactivation time courses during voltage pulses (C). Different symbols and bars in the figure denote mean ±S.E.M. from five independent experiments. Smooth curves in A are drawn according to the best fits of data points (normalized to the maximum) by the Boltzmann equation: 1/[1 + exp(−(Vm − V0.5)/s), where Vm, V0.5 and s denote membrane potential, half-activation voltage and slope factor, respectively. The values of these parameters are 47.542, 48.825 and 3.697 mV, and 56.582, 67.327 and 52.844 mV for 0.3, 1 and 5 mM [Ca2+]o, respectively.
Ijms 22 08513 g001
Figure 2. Mathematical formulation of voltage-and Ca2+-dependent gating of E7K-TRPM4 channels. (A) and (B); The voltage-dependencies of α and β were empirically fitted by the exponentials with the first and second-order Vm polynominals, respectively. (CE); the coefficients for each polynominal (C0, C1; K0–K2) obtained in A and B were further empirically fitted with respect to intracellular Ca2+ concentration ([Ca2+]i) by linear or exponential fittings. [Ca2+]i was determined by separately preformed [Ca2+]i measurements with the fura-2 digital fluorescence imaging [13].
Figure 2. Mathematical formulation of voltage-and Ca2+-dependent gating of E7K-TRPM4 channels. (A) and (B); The voltage-dependencies of α and β were empirically fitted by the exponentials with the first and second-order Vm polynominals, respectively. (CE); the coefficients for each polynominal (C0, C1; K0–K2) obtained in A and B were further empirically fitted with respect to intracellular Ca2+ concentration ([Ca2+]i) by linear or exponential fittings. [Ca2+]i was determined by separately preformed [Ca2+]i measurements with the fura-2 digital fluorescence imaging [13].
Ijms 22 08513 g002
Figure 3. Voltage- and Ca2+-dependent profiles of Po, α and β. Po, α and β were recalculated over wide ranges of Vm and [Ca2+]i by using the mathematical expressions of α and β shown above. For wild-type TRPM4, similar mathematical expressions for α and β which had been obtained from wild-type TRPM4 expressing HEK293 cells [13] were used to calculate Po, α and β. A and B; surface plots of Po, α and β (from left to right) for wild-type TRPM4 channel against two variables Vm and [Ca2+]i (A) and those for E7K-mutant TRPM4 channel (B). Note that Po contour is less steep for E7K mutant with increased voltage-dependence of α and prominently decreased voltage and Ca2+ dependency of β which stay low over the whole Vm and [Ca2+]i ranges. These results indicate a preferred sojourn in the open state of E7K channel due to accelerated C-O transition and greatly suppressed O-C transition.
Figure 3. Voltage- and Ca2+-dependent profiles of Po, α and β. Po, α and β were recalculated over wide ranges of Vm and [Ca2+]i by using the mathematical expressions of α and β shown above. For wild-type TRPM4, similar mathematical expressions for α and β which had been obtained from wild-type TRPM4 expressing HEK293 cells [13] were used to calculate Po, α and β. A and B; surface plots of Po, α and β (from left to right) for wild-type TRPM4 channel against two variables Vm and [Ca2+]i (A) and those for E7K-mutant TRPM4 channel (B). Note that Po contour is less steep for E7K mutant with increased voltage-dependence of α and prominently decreased voltage and Ca2+ dependency of β which stay low over the whole Vm and [Ca2+]i ranges. These results indicate a preferred sojourn in the open state of E7K channel due to accelerated C-O transition and greatly suppressed O-C transition.
Ijms 22 08513 g003
Figure 4. Noise spectrum and single channel analyses of E7K channel gating. (A) representative spectral density plots of wild-type (a) and E7K-TRPM4 (b) whole-currents at −60 mV at 0.52 μM [Ca2+]i. Each plot is fitted by the equation: s(f) = s(0)/[1 + (f/fc)2] (see the Methods). The value of time constant (τ) is calculated according to the relationship: τ = 1/2πfc. The τ value would provide a reasonable estimate for the open life-time, provided that the rate of closing is much larger than that of opening, i.e., near the resting membrane potential (see Figure 3). To minimize the desensitization/rundown of currents, sharp patch electrodes were used to infuse desired [Ca2+]i into the cell. (B) summary of spectrum analysis such as shown in A. Open and filled circles and bars indicate the means ± S.E.M.s of τ at 0.28, 0.52 and 1.17 μM [Ca2+]i (n = 4–8). Each condition was tested independently. Only the data at the same [Ca2+]i was compared by unpaired t-test. *, **: p < 0.05 and 0.005. (C) representative traces of single channel activities at −60 mV (upper) and their respective dwell-time histograms (lower) for wild-type (a) and E7K- TRPM4 (b) channels. Their dwell-time histograms (shown in dots) are fitted by a sum of two exponentials with fast and slow time constants of τ1 and τ2. (D) averaged values of τ1 and τ2 with S.E.M bars for wild-type and E7K-mutant TRPM4 channels. n = 7. p values are the results of unpaired t-tests.
Figure 4. Noise spectrum and single channel analyses of E7K channel gating. (A) representative spectral density plots of wild-type (a) and E7K-TRPM4 (b) whole-currents at −60 mV at 0.52 μM [Ca2+]i. Each plot is fitted by the equation: s(f) = s(0)/[1 + (f/fc)2] (see the Methods). The value of time constant (τ) is calculated according to the relationship: τ = 1/2πfc. The τ value would provide a reasonable estimate for the open life-time, provided that the rate of closing is much larger than that of opening, i.e., near the resting membrane potential (see Figure 3). To minimize the desensitization/rundown of currents, sharp patch electrodes were used to infuse desired [Ca2+]i into the cell. (B) summary of spectrum analysis such as shown in A. Open and filled circles and bars indicate the means ± S.E.M.s of τ at 0.28, 0.52 and 1.17 μM [Ca2+]i (n = 4–8). Each condition was tested independently. Only the data at the same [Ca2+]i was compared by unpaired t-test. *, **: p < 0.05 and 0.005. (C) representative traces of single channel activities at −60 mV (upper) and their respective dwell-time histograms (lower) for wild-type (a) and E7K- TRPM4 (b) channels. Their dwell-time histograms (shown in dots) are fitted by a sum of two exponentials with fast and slow time constants of τ1 and τ2. (D) averaged values of τ1 and τ2 with S.E.M bars for wild-type and E7K-mutant TRPM4 channels. n = 7. p values are the results of unpaired t-tests.
Ijms 22 08513 g004
Figure 5. Ca2+-dependent activation curves for wild-type and E7K-TRPM4 currents. Internal solutions containing three different [Ca2+]i were introduced via sharp thin patch electrodes into the cell by breaking the giga-sealed membrane at a holding potential of −60 mV. (A) representative traces of 0.52 μM [Ca2+]i -induced currents from HEK293 cells expressing either wild-type (upper) or E7K-mutant (lower) TRPM4 channels. The arrows indicate the timepoints of ‘break-in’. Vertical positive deflections in the figure show the membrane currents evoked by ramp voltages. (B) the relationships between the TRPM4 current amplitude (at −60 mV) and intracellularly perfused [Ca2+]i. Open and filled circles represent wild-type and E7K-mutant TRPM4 channels. n = 4–12. Each condition was tested independently. Only the data at the same [Ca2+]i was compared by unpaired t-test. *, **: p < 0.05 and 0.005.
Figure 5. Ca2+-dependent activation curves for wild-type and E7K-TRPM4 currents. Internal solutions containing three different [Ca2+]i were introduced via sharp thin patch electrodes into the cell by breaking the giga-sealed membrane at a holding potential of −60 mV. (A) representative traces of 0.52 μM [Ca2+]i -induced currents from HEK293 cells expressing either wild-type (upper) or E7K-mutant (lower) TRPM4 channels. The arrows indicate the timepoints of ‘break-in’. Vertical positive deflections in the figure show the membrane currents evoked by ramp voltages. (B) the relationships between the TRPM4 current amplitude (at −60 mV) and intracellularly perfused [Ca2+]i. Open and filled circles represent wild-type and E7K-mutant TRPM4 channels. n = 4–12. Each condition was tested independently. Only the data at the same [Ca2+]i was compared by unpaired t-test. *, **: p < 0.05 and 0.005.
Ijms 22 08513 g005
Figure 6. E7K mutation causes conduction block. (A) modified human Purkinje fiber Trovato2020 model (for details, see the Methods) was employed to simulate action potential (AP) and conduction velocity (CV). A; adaption of TRPM4 gating kinetics to the original Trovato model. (B) single-cell AP simulations by modified Trovato model incorporated with the gating kinetics of either wild-type (left) or E7K-mutant (right). The results are overlaid for each. (C) relationships between conduction velocity (CV) and channel density (or maximal activity) with intracellular conductivity of 1.7 mS/cm at 1 or 2 Hz pacing. The channel density is presented as a multiple (‘a’ ×) of control (i.e., 7.02 × 10−8 Litre/farad/ms) for wild-type and E7K-TRPM4 channels. The horizontal bar on the top panel in C denotes a 3 cm-long cable (stimulating point of 0–0.05 cm is shown as a large filled circle on the left). To avoid boundary effects, CV was calculated between 0.2 and 2.8 cm. (D) relationships of maximal upstroke slope of propagated AP (dV/dtmax; measured at 2.8 cm) or resting membrane potential (RMP) against TRPM4 channel density which is shown as the fold increase (‘a’ ×) as in (C).
Figure 6. E7K mutation causes conduction block. (A) modified human Purkinje fiber Trovato2020 model (for details, see the Methods) was employed to simulate action potential (AP) and conduction velocity (CV). A; adaption of TRPM4 gating kinetics to the original Trovato model. (B) single-cell AP simulations by modified Trovato model incorporated with the gating kinetics of either wild-type (left) or E7K-mutant (right). The results are overlaid for each. (C) relationships between conduction velocity (CV) and channel density (or maximal activity) with intracellular conductivity of 1.7 mS/cm at 1 or 2 Hz pacing. The channel density is presented as a multiple (‘a’ ×) of control (i.e., 7.02 × 10−8 Litre/farad/ms) for wild-type and E7K-TRPM4 channels. The horizontal bar on the top panel in C denotes a 3 cm-long cable (stimulating point of 0–0.05 cm is shown as a large filled circle on the left). To avoid boundary effects, CV was calculated between 0.2 and 2.8 cm. (D) relationships of maximal upstroke slope of propagated AP (dV/dtmax; measured at 2.8 cm) or resting membrane potential (RMP) against TRPM4 channel density which is shown as the fold increase (‘a’ ×) as in (C).
Ijms 22 08513 g006
Figure 7. Heterogeneous expression of E7K causes partial conduction block. (A,B) relationships of CV (A), dV/dtmax (B) or RMP (B) against the density of wild-type or E7K mutant TRPM4 channels with the intracellular conductivity of 8.5 mS/cm. Definitions of the ordinate and abscissa are the same as in Figure 6C,D. (C) CV-channel density relationship for E7K-mutant TRPM4 channel under spatial heterogenous conditions consisting of two different densities at the middle (‘b’ ×) and both ends (‘a’ ×). For better visualization, only the cases for 1× to 3× with intracellular conductivity of 8.5 mS/cm and 2 Hz-pacing are shown. (D) spatiotemporal profile (surface plot) of partial conduction block (2:1) observed for the case of open squares in C in which the densities in the middle and ends are set to 4× and 3×, respectively. The last 25 out of 50 stimulations are shown.
Figure 7. Heterogeneous expression of E7K causes partial conduction block. (A,B) relationships of CV (A), dV/dtmax (B) or RMP (B) against the density of wild-type or E7K mutant TRPM4 channels with the intracellular conductivity of 8.5 mS/cm. Definitions of the ordinate and abscissa are the same as in Figure 6C,D. (C) CV-channel density relationship for E7K-mutant TRPM4 channel under spatial heterogenous conditions consisting of two different densities at the middle (‘b’ ×) and both ends (‘a’ ×). For better visualization, only the cases for 1× to 3× with intracellular conductivity of 8.5 mS/cm and 2 Hz-pacing are shown. (D) spatiotemporal profile (surface plot) of partial conduction block (2:1) observed for the case of open squares in C in which the densities in the middle and ends are set to 4× and 3×, respectively. The last 25 out of 50 stimulations are shown.
Ijms 22 08513 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hu, Y.; Li, Q.; Shen, Y.; Fujita, T.; Zhu, X.; Inoue, R. Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block. Int. J. Mol. Sci. 2021, 22, 8513. https://doi.org/10.3390/ijms22168513

AMA Style

Hu Y, Li Q, Shen Y, Fujita T, Zhu X, Inoue R. Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block. International Journal of Molecular Sciences. 2021; 22(16):8513. https://doi.org/10.3390/ijms22168513

Chicago/Turabian Style

Hu, Yaopeng, Qin Li, Yanghua Shen, Takayuki Fujita, Xin Zhu, and Ryuji Inoue. 2021. "Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block" International Journal of Molecular Sciences 22, no. 16: 8513. https://doi.org/10.3390/ijms22168513

APA Style

Hu, Y., Li, Q., Shen, Y., Fujita, T., Zhu, X., & Inoue, R. (2021). Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block. International Journal of Molecular Sciences, 22(16), 8513. https://doi.org/10.3390/ijms22168513

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop