Mutiscale Computational Analysis of the Effect on Heart Rate of a HCN4 Gene Double Mutation: from the Single Channel to the Clinical Phenotype

This work aims to assess: I) the effects of the 1479V/A485E HCN4 channel double mutation (DM), II) the role of cellular coupling (p) and III) the role of cellular heterogeneity (σ) on systems of increasing complexity: ionic channel, single cell, 1D fibre and 2D tissue of the human SAN. The Fabbri et al. model was used to describe the human SAN cell and to build the 1D and 2D models, which are constituted of 100 and 2500 (50×50) cells, respectively. σ = 0.05, 0.1, 0.1873, 0.3, 0.4 and p = 10, 100, 1000, 10000 and ∞ MΩ. m were simulated, in order to test their effect on the heart rate (HR). The reduction of I f current due to the heterozygous condition (gj(WT0.5 + DM) = 45% of gf(WT)) leads to an increase of the cycle length of the simulated action potential of a single cell (924 vs 814 ms; + 14%). This WT0.5+DM bradycardic effect is confirmed also by the 1D model (802 vs 690 ms; +14%) as well as by the 2D one (908 vs 794 ms; +14%). These results were obtained forσ = 0.1873 and p = 100 MΩ· m (50 MΩ. m in 1D). Other combinations of σ and p can provide changes in HR greater than 50%, highlighting the importance of these two parameters in the establishment of a physiologic pacing in the human SAN.