Region-dependent dynamics of cAMP response element-binding protein phosphorylation in the basal ganglia

The cAMP response element-binding protein (CREB) is an activity-dependent transcription factor that is involved in neural plasticity. The kinetics of CREB phosphorylation have been suggested to be important for gene activation, with sustained phosphorylation being associated with downstream gene expression. If so, the duration of CREB phosphorylation might serve as an indicator for time-sensitive plastic changes in neurons. To screen for regions potentially involved in dopamine-mediated plasticity in the basal ganglia, we used organotypic slice cultures to study the patterns of dopamine- and calcium-mediated CREB phosphorylation in the major subdivisions of the striatum. Different durations of CREB phosphorylation were evoked in the dorsal and ventral striatum by activation of dopamine D1-class receptors. The same D1 stimulus elicited (i) transient phosphorylation (≤15 min) in the matrix of the dorsal striatum; (ii) sustained phosphorylation (≤2 hr) in limbic-related structures including striosomes, the nucleus accumbens, the fundus striati, and the bed nucleus of the stria terminalis; and (iii) prolonged phosphorylation (up to 4 hr or more) in cellular islands in the olfactory tubercle. Elevation of Ca2+ influx by stimulation of L-type Ca2+ channels, NMDA, or KCl induced strong CREB phosphorylation in the dorsal striatum but not in the olfactory tubercle. These findings differentiate the response of CREB to dopamine and calcium signals in different striatal regions and suggest that dopamine-mediated CREB phosphorylation is persistent in limbic-related regions of the neonatal basal ganglia. The downstream effects activated by persistent CREB phosphorylation may include time-sensitive neuroplasticity modulated by dopamine.

Calcium electrogenesis in distal apical dendrites of layer 5 pyramidal cells at a critical frequency of back-propagating action potentials

Action potentials in juvenile and adult rat layer-5 neocortical pyramidal neurons can be initiated at both axonal and distal sites of the apical dendrite. However, little is known about the interaction between these two initiation sites. Here, we report that layer 5 pyramidal neurons are very sensitive to a critical frequency of back-propagating action potentials varying between 60 and 200 Hz in different neurons. Bursts of four to five back-propagating action potentials above the critical frequency elicited large regenerative potentials in the distal dendritic initiation zone. The critical frequency had a very narrow range (10–20 Hz), and the dendritic regenerative activity led to further depolarization at the soma. The dendritic frequency sensitivity was suppressed by blockers of voltage-gated calcium channels, and also by synaptically mediated inhibition. Calcium-fluorescence imaging revealed that the site of largest transient increase in intracellular calcium above the critical frequency was located 400–700 μm from the soma at the site for initiation of calcium action potentials. Thus, the distal dendritic initiation zone can interact with the axonal initiation zone, even when inputs to the neuron are restricted to regions close to the soma, if the output of the neuron exceeds a critical frequency.

Construction of a high-affinity receptor site for dihydropyridine agonists and antagonists by single amino acid substitutions in a non-l-type Ca2+ channel

The activity of l-type Ca2+ channels is increased by dihydropyridine (DHP) agonists and inhibited by DHP antagonists, which are widely used in the therapy of cardiovascular disease. These drugs bind to the pore-forming α1 subunits of l-type Ca2+ channels. To define the minimal requirements for DHP binding and action, we constructed a high-affinity DHP receptor site by substituting a total of nine amino acid residues from DHP-sensitive l-type α1 subunits into the S5 and S6 transmembrane segments of domain III and the S6 transmembrane segment of domain IV of the DHP-insensitive P/Q-type α1A subunit. The resulting chimeric α1A/DHPS subunit bound DHP antagonists with high affinity in radioligand binding assays and was inhibited by DHP antagonists with high affinity in voltage clamp experiments. Substitution of these nine amino acid residues yielded 86% of the binding energy of the l-type α1C subunit and 92% of the binding energy of the l-type α1S subunit for the high-affinity DHP antagonist PN200–110. The activity of chimeric Ca2+ channels containing α1A/DHPS was increased 3.5 ± 0.7-fold by the DHP agonist (−)Bay K8644. The effect of this agonist was stereoselective as in l-type Ca2+ channels since (+) Bay K8644 inhibited the activity of α1A/DHPS. The results show conclusively that DHP agonists and antagonists bind to a single receptor site at which they have opposite effects on Ca2+ channel activity. This site contains essential components from both domains III and IV, consistent with a domain interface model for binding and allosteric modulation of Ca2+ channel activity by DHPs.

Allosteric modulation of Ca2+ channels by G proteins, voltage-dependent facilitation, protein kinase C, and Cavβ subunits

N-type and P/Q-type Ca2+ channels are inhibited by neurotransmitters acting through G protein-coupled receptors in a membrane-delimited pathway involving Gβγ subunits. Inhibition is caused by a shift from an easily activated “willing” (W) state to a more-difficult-to-activate “reluctant” (R) state. This inhibition can be reversed by strong depolarization, resulting in prepulse facilitation, or by protein kinase C (PKC) phosphorylation. Comparison of regulation of N-type Ca2+ channels containing Cav2.2a α1 subunits and P/Q-type Ca2+ channels containing Cav2.1 α1 subunits revealed substantial differences. In the absence of G protein modulation, Cav2.1 channels containing Cavβ subunits were tonically in the W state, whereas Cav2.1 channels without β subunits and Cav2.2a channels with β subunits were tonically in the R state. Both Cav2.1 and Cav2.2a channels could be shifted back toward the W state by strong depolarization or PKC phosphorylation. Our results show that the R state and its modulation by prepulse facilitation, PKC phosphorylation, and Cavβ subunits are intrinsic properties of the Ca2+ channel itself in the absence of G protein modulation. A common allosteric model of G protein modulation of Ca2+-channel activity incorporating an intrinsic equilibrium between the W and R states of the α1 subunits and modulation of that equilibrium by G proteins, Cavβ subunits, membrane depolarization, and phosphorylation by PKC accommodates our findings. Such regulation will modulate transmission at synapses that use N-type and P/Q-type Ca2+ channels to initiate neurotransmitter release.

Changes in voltage activation, Cs+ sensitivity, and ion permeability in H5 mutants of the plant K+ channel KAT1.

KAT1 is a voltage-dependent inward rectifying K+ channel cloned from the higher plant Arabidopsis thaliana [Anderson, J. A., Huprikar, S. S., Kochian, L. V., Lucas, W. J. & Gaber, R. F. (1992) Proc. Natl. Acad. Sci. USA 89, 3736-3740]. It is related to the Shaker superfamily of K+ channels characterized by six transmembrane spanning domains (S1-S6) and a putative pore-forming region between S5 and S6 (H5). The 115 region between Pro-247 and Pro-271 in KAT1 contains 14 additional amino acids when compared with Shaker [Aldrich, R. W. (1993) Nature (London) 362, 107-108]. We studied various point mutations introduced into H5 to determine whether voltage-dependent plant and animal K+ channels share similar pore structures. Through heterologous expression in Xenopus oocytes and voltage-clamp analysis combined with phenotypic analysis involving a potassium transport-defective Saccharomyces cerevisiae strain, we investigated the selectivity filter of the mutants and their susceptibility toward inhibition by cesium and calcium ions. With respect to electrophysiological properties, KAT1 mutants segregated into three groups: (i) wild-type-like channels, (ii) channels modified in selectivity and Cs+ or Ca2+ sensitivity, and (iii) a group that was additionally affected in its voltage dependence. Despite the additional 14 amino acids in H5, this motif in KAT1 is also involved in the formation of the ion-conducting pore because amino acid substitutions at Leu-251, Thr-256, Thr-259, and Thr-260 resulted in functional channels with modified ionic selectivity and inhibition. Creation of Ca2+ sensitivity and an increased susceptibility to Cs+ block through mutations within the narrow pore might indicate that both blockers move deeply into the channel. Furthermore, mutations close to the rim of the pore affecting the half-activation potential (U1/2) indicate that amino acids within the pore either interact with the voltage sensor or ion permeation feeds back on gating.

Isoform-specific interaction of the alpha1A subunits of brain Ca2+ channels with the presynaptic proteins syntaxin and SNAP-25.

Presynaptic Ca2+ channels are crucial elements in neuronal excitation-secretion coupling. In addition to mediating Ca2+ entry to initiate transmitter release, they are thought to interact directly with proteins of the synaptic vesicle docking/fusion machinery. Here we report isoform-specific, stoichiometric interaction of the BI and rbA isoforms of the alpha1A subunit of P/Q-type Ca2+ channels with the presynaptic membrane proteins syntaxin and SNAP-25 in vitro and in rat brain membranes. The BI isoform binds to both proteins, while only interaction with SNAP-25 can be detected in vitro for the rbA isoform. The synaptic protein interaction ("synprint") site involves two adjacent segments of the intracellular loop connecting domains II and III between amino acid residues 722 and 1036 of the BI sequence. This interaction is competitively blocked by the corresponding region of the N-type Ca2+ channel, indicating that these two channels bind to overlapping regions of syntaxin and SNAP-25. Our results provide a molecular basis for a physical link between Ca2+ influx into nerve terminals and subsequent exocytosis of neurotransmitters at synapses that have presynaptic Ca2+ channels containing alpha1A subunits.

Receptor stimulation causes slow inhibition of IRK1 inwardly rectifying K+ channels by direct protein kinase A-mediated phosphorylation.

Strongly rectifying IRK-type inwardly rectifying K+ channels are involved in the control of neuronal excitability in the mammalian brain. Whole-cell patch-clamp experiments show that cloned rat IRK1 (Kir 2.1) channels, when heterologously expressed in mammalian COS-7 cells, are inhibited following the activation of coexpressed serotonin (5-hydroxytryptamine) type 1A receptors by receptor agonists. Inhibition is mimicked by internal perfusion with GTP[gamma-S] and elevation of internal cAMP concentrations. Addition of the catalytic subunits of protein kinase A (PKA) to the internal recording solution causes complete inhibition of wild-type IRK1 channels, but not of mutant IRK1(S425N) channels in which a C-terminal PKA phosphorylation site has been removed. Our data suggest that in the nervous system serotonin may negatively control IRK1 channel activity by direct PKA-mediated phosphorylation.

Probing the transmembrane topology of cyclic nucleotide-gated ion channels with a gene fusion approach.

Cyclic nucleotide-gated (CNG) cation channels contain two short sequence motifs--a residual voltage-sensor (S4) and a pore-forming (P) segment--that are reminiscent of similar segments in voltage-activated Shaker-type K+ channels. It has been tacitly assumed that CNG channels and this K+ channel subfamily share a common overall topology, characterized by a hydrophobic domain comprising six membrane-spanning segments. We have systematically investigated the topology of CNG channels from bovine rod photoreceptor and Drosophila melanogaster by a gene fusion approach using the bacterial reporter enzymes alkaline phosphatase and beta-galactosidase, which are active only in the periplasm and only in the cytoplasm, respectively. Enzymatic activity was determined after expression of fusion constructs in Escherichia coli. CNG channels were found to have six membrane-spanning segments, suggesting that CNG and Shaker-type K+ channels, albeit distant relatives within a gene superfamily of ion channels, share a common topology.

Modulação da diferenciação neural de células troco embrionárias por transientes de cálcio intracelulares: papéis dos receptores purinérgicos e de canais de cálcio voltagem-dependentes

Receptores purinérgicos e canais de cálcio voltagem-dependentes estão envolvidos em diversos processos biológicos como na gastrulação, durante o desenvolvimento embrionário, e na diferenciação neural. Quando ativados, canais de cálcio voltagem-dependentes e receptores purinérgicos do tipo P2, ativados por nucleotídeos, desencadeiam transientes de cálcio intracelulares controlando diversos processos biológicos. Neste trabalho, nós estudamos a participação de canais de cálcio voltagem-dependentes e receptores do tipo P2 na geração de transientes de cálcio espontâneos e sua regulação na expressão de fatores de transcrição relacionados com a neurogênese utilizando como modelo células tronco (CTE) induzidas à diferenciação em células tronco neurais (NSC) com ácido retinóico. Descrevemos que CTE indiferenciadas podem ter a proliferação acelerada pela ativação de receptores P2X7, enquanto que a expressão e a atividade desse receptor precisam ser inibidas para o progresso da diferenciação em neuroblasto. Além disso, ao longo da diferenciação neural, por análise em tempo real dos níveis de cálcio intracelular livre identificamos 3 padrões de oscilações espontâneas de cálcio (onda, pico e unique), e mostramos que ondas e picos tiveram a frequência e amplitude aumentadas conforme o andamento da diferenciação. Células tratadas com o inibidor do receptor de inositol 1,4,5-trifosfato (IP3R), Xestospongin C, apresentaram picos mas não ondas, indicando que ondas dependem exclusivamente de cálcio oriundo do retículo endoplasmático pela ativação de IP3R. NSC de telencéfalo de embrião de camundongos transgênicos ou pré-diferenciadas de CTE tratadas com Bz-ATP, o agonista do receptor P2X7, e com 2SUTP, agonista de P2Y2 e P2Y4, aumentaram a frequência e a amplitude das oscilações espontâneas de cálcio do tipo pico. Dados, obtidos por microscopia de luminescência, da expressão em tempo real de gene repórter luciferase fusionado à Mash1 e Ngn2 revelou que a ativação dos receptores P2Y2/P2Y4 aumentou a expressão estável de Mash1 enquanto que ativação do receptor P2X7 levou ao aumento de Ngn2. Além disso, células na presença do quelante de cálcio extracelular (EGTA) ou do depletor dos estoques intracelulares de cálcio do retículo endoplasmático (thapsigargin) apresentaram redução na expressão de Mash1 e Ngn2, indicando que ambos são regulados pela sinalização de cálcio. A investigação dos canais de cálcio voltagem-dependentes demonstrou que o influxo de cálcio gerado por despolarização da membrana de NSC diferenciadas de CTE é decorrente da ativação de canais de cálcio voltagem-dependentes do tipo L. Além disso, esse influxo pode controlar o destino celular por estabilizar expressão de Mash1 e induzir a diferenciação neuronal por fosforilação e translocação do fator de transcrição CREB. Esses dados sugerem que os receptores P2X7, P2Y2, P2Y4 e canais de cálcio voltagem-dependentes do tipo L podem modular as oscilações espontâneas de cálcio durante a diferenciação neural e consequentemente alteram o padrão de expressão de Mash1 e Ngn2 favorecendo a decisão do destino celular neuronal.

Estudo de associação entre disfunção neurocognitiva, estresse oxidativa e polimorfismos em pacientes jovens com Transtornos Bipolar tipo I

O Transtorno Bipolar (TB) tipo I é uma doença caracterizada por episódios de mania e depressão recorrentes com importante prejuízo do funcionamento global e comprometimento das funções cognitivas. Além disso, sabe-se que o número de episódios de humor patológico ao longo da vida pode também influenciar o funcionamento cognitivo destes sujeitos. Neste cenário, ocorreu a necessidade de se investigar marcadores genéticos para disfunção cognitiva no TB com o objetivo de estudar este fenômeno. Dentre os potenciais genes responsáveis por influenciar a cognição destacam-se os polimorfismos funcionais do fator neurotrófico derivado do cérebro (BDNF), da catecol-O-metiltransferase (COMT), da apolipoproteína-E (APOE) e do canal de cálcio de baixa voltagem subunidade 1-C (CACNA1C). Sabe-se, também, que no TB os marcadores de estresse oxidativo estão aumentados durante todas as fases da doença, entretanto, não é claro qual impacto destes na disfunção cognitiva de indivíduos com TB. O objetivo dessa tese foi avaliar o desempenho cognitivo de pacientes jovens com bipolaridade tipo I e sua associação com o genótipo de BDNF, COMT, APOE e CACNA1C e também com os níveis plasmáticos de oxidação da guanosina (8-OHdG) e citosina (5-Mec) durante os episódios de humor, eutimia e em controles. Para investigar essa associação foram incluídos 116 pacientes (79 em episódio de humor patológico e 37 eutímicos) com diagnóstico de TB tipo I (DSMIV-TR); 97 controles saudáveis foram submetidos à avaliação neuropsicológica e coleta de sangue para extração de DNA visando genotipagem para BDNF (rs6265), COMT (rs4680; rs165599), APOE (rs429358 e rs7412), CACNA1C (rs1006737), 8-OhdG e 5-Mec. A análise dos dados obtidos revelou que pacientes portadores do genótipo Met/Met rs4680/rs165599 do COMT apresentam comprometimento cognitivo mais grave (função executiva, fluência verbal, memória e inteligência) comparado ao genótipo Val/Met ou Val/Val durante episódios maníacos ou mistos. Na mesma direção destes resultados, verificou-se que pacientes portadores do alelo Met rs4680 do COMT apresentam comprometimento do reconhecimento de emoções faciais em episódios de mania e depressão. Nenhum efeito do COMT foi observado em controles. O alelo de risco Met do CACNA1C se associou a um pior comprometimento executivo independente dos sintomas maníacos ou depressivos no TB, porém nenhum efeito se observou nos controles. O alelo Met do BDNF rs6265 ou a presença do alelo 4 da APOE não representa um fator que identifique um grupo com desempenho cognitivo diferenciado durante as fases do TB ou em controles. Sujeitos com TB apresentaram níveis mais elevados de 8-OHdG e tais níveis eram diretamente proporcionais ao número de episódios maníacos ao longo da vida, sugerindo um papel dos episódios hiperdopaminérgicos na oxidação das bases de DNA. Concluiu-se que a genotipagem para COMT e CACNA1C em pacientes com TB pode identificar um grupo de pacientes associados a pior disfunção cognitiva durante as fases maníacas e mistas do TB. Tal dado pode ser um indicador do envolvimento do sistema dopaminérgico e dos canais de cálcio de baixa voltagem na fisiopatologia da disfunção cognitiva no TB e deve ser explorado em outros estudos

Expression and cellular localization of the voltage-gated calcium channel α2δ3 in the rodent retina

High-voltage-activated calcium channels are hetero-oligomeric protein complexes that mediate multiple cellular processes, including the influx of extracellular Ca2+, neurotransmitter release, gene transcription, and synaptic plasticity. These channels consist of a primary α1 pore-forming subunit, which is associated with an extracellular α2δ subunit and an intracellular β auxiliary subunit, which alter the gating properties and trafficking of the calcium channel. The cellular localization of the α2δ3 subunit in the mouse and rat retina is unknown. In this study using RT-PCR, a single band at ∼305 bp corresponding to the predicted size of the α2δ3 subunit fragment was found in mouse and rat retina and brain homogenates. Western blotting of rodent retina and brain homogenates showed a single 123-kDa band. Immunohistochemistry with an affinity-purified antibody to the α2δ3 subunit revealed immunoreactive cell bodies in the ganglion cell layer and inner nuclear layer and immunoreactive processes in the inner plexiform layer and the outer plexiform layer. α2δ3 immunoreactivity was localized to multiple cell types, including ganglion, amacrine, and bipolar cells and photoreceptors, but not horizontal cells. The expression of the α2δ3 calcium channel subunit to multiple cell types suggests that this subunit participates widely in Ca-channel-mediated signaling in the retina.

Conotoxins as selective inhibitors of neuronal ion channels, receptors and transporters

Cone snails have evolved a vast array of peptide toxins for prey capture and defence. These peptides are directed against a wide variety of pharmacological targets, making them an invaluable source of ligands for studying the properties of these targets in normal and diseased states. A number of these peptides have shown efficacy in vivo, including inhibitors of calcium channels, the norepinephrine transporter, nicotinic acetylcholine receptors, NMDA receptors and neurotensin receptors, with several having undergone pre-clinical or clinical development for the treatment of pain.

Structures of mu O-conotoxins from Conus marmoreus - Inhibitors of tetrodotoxin (TTX)-sensitive and TTX-resistant sodium channels in mammalian sensory neurons

The muO-conotoxins are an intriguing class of conotoxins targeting various voltage-dependent sodium channels and molluscan calcium channels. In the current study, we have shown MrVIA and MrVIB to be the first known peptidic inhibitors of the transient tetrodotoxin-resistant (TTX-R) Na+ current in rat dorsal root ganglion neurons, in addition to inhibiting tetrodotoxin-sensitive Na+ currents. Human TTX-R sodium channels are a therapeutic target for indications such as pain, highlighting the importance of the muO-conotoxins as potential leads for drug development. Furthermore, we have used NMR spectroscopy to provide the first structural information on this class of conotoxins. MrVIA and MrVIB are hydrophobic peptides that aggregate in aqueous solution but were solubilized in 50% acetonitrile/water. The three-dimensional structure of MrVIB consists of a small beta-sheet and a cystine knot arrangement of the three-disulfide bonds. It contains four backbone loops between successive cysteine residues that are exposed to the solvent to varying degrees. The largest of these, loop 2, is the most disordered part of the molecule, most likely due to flexibility in solution. This disorder is the most striking difference between the structures of MrVIB and the known delta- and omega-conotoxins, which along with the muO-conotoxins are members of the O superfamily. Loop 2 of omega-conotoxins has previously been shown to contain residues critical for binding to voltage-gated calcium channels, and it is interesting to speculate that the flexibility observed in MrVIB may accommodate binding to both sodium and molluscan calcium channels.

Tertiapin-Q blocks recombinant and native large conductance K+ channels in a use-dependent manner

Tertiapin, a short peptide from honey bee venom, has been reported to specifically block the inwardly rectifying K+ (Kir) channels, including G protein-coupled inwardly rectifying potassium channel (GIRK) 1 + GIRK4 heteromultimers and ROMK1 homomultimers. In the present study, the effects of a stable and functionally similar derivative of tertiapin, tertiapin-Q, were examined on recombinant human voltage-dependent Ca2+-activated large conductance K+ channel (BK or MaxiK; alpha-subunit or hSlo1 homomultimers) and mouse inwardly rectifying GIRK1 + GIRK2 (i.e., Kir3.1 and Kir3.2) heteromultimeric K+ channels expressed in Xenopus oocytes and in cultured newborn mouse dorsal root ganglion (DRG) neurons. In two-electrode voltage-clamped oocytes, tertiapin-Q (1-100 nM) inhibited BK-type K+ channels in a use- and concentration-dependent manner. We also confirmed the inhibition of recombinant GIRK1 + GIRK2 heteromultimers by tertiapin-Q, which had no effect on endogenous depolarization- and hyperpolarization-activated currents sensitive to extracellular divalent cations (Ca2+, Mg2+, Zn2+, and Ba2+) in defolliculated oocytes. In voltage-clamped DRG neurons, tertiapin-Q voltage- and use-dependently inhibited outwardly rectifying K+ currents, but Cs+-blocked hyperpolarization-activated inward currents including I-H were insensitive to tertiapin-Q, baclofen, barium, and zinc, suggesting absence of functional GIRK channels in the newborn. Under current-clamp conditions, tertiapin-Q blocked the action potential after hyperpolarization (AHP) and increased action potential duration in DRG neurons. Taken together, these results demonstrate that the blocking actions of tertiapin-Q are not specific to Kir channels and that the blockade of recombinant BK channels and native neuronal AHP currents is use-dependent. Inhibition of specific types of Kir and voltage-dependent Ca2+-activated K+ channels by tertiapin-Q at nanomolar range via different mechanisms may have implications in pain physiology and therapy.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.