Effect Of Site Selective Ion Channel Blocking on Action Potential Krishnendu Pal1,2,+ and Gautam Gangopadhyay2,∗ 1Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India. 2S N Bose National Centre for Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata-700106, India. (Dated: October 2, 2018) In this work we have theoretically investigated how the action potential generation and its as- sociated intrinsic properties are affected in presence of ion channel blockers by adapting Gillepie’s stochastic simulation technique on a very basic neuron of Hodgkin-Huxley type. With a simple extension of the Hodgkin-Huxley Markov model we have mainly investigated three types of drug blocking mechanism such as (i) only sodium channel blocking, (ii) only potassium channel blocking and (iii) dual type blocking and showed that the major experimental and physiological observations such as ionic currents, spiking frequency trends, change in action potential shape and duration, altered gating dynamics etc due to the presence of ion channel blockers can be well reproduced. Our results show that the nature of action potential termination process in presence of sodium and potassium channel blockers are distinct and physiologically very different from each other. We have found that although the sodium and potassium ionic currents have interdependent relationship over a course of an action potential but sodium and potassium channel blockers have distinct signatures on ionic currents. In presence of only sodium channel blockers the frequency of action potential generation falls off exponentially with increasing drug affinity, whereas in contrast, for only potas- sium channel blockers initially an enhanced spiking activity of action potential is found followed by a gradual decrease of the spiking frequency as the drug affinity increases. In case of dual type blockers with equal sodium and potassium channel binding affinity, the spiking frequency passes through maxima and minima due to the competition between channel number fluctuation and over- all sodium and potassium conductances. We have found that sodium channel blockers shorten the duration of action potential while the potassium channel blockers delay it which are of great phys- iological and pharmacological importance. We have also shown how the ion channel blockers alter the gating dynamics and such altered gating itself modulates the ion channel blocking which opens the possibility of finding fundamental informations regarding probabilistic and dynamical features. Some experimental results of ion channel blocking in diverse systems have been validated through our site selected binding scheme. Many other types of blocking mechanisms such as closed state blocking, inactive state blocking etc can also be explored using our method with desired level of structural and functional details. Keywords: Ion Channels; Sodium Channel Blockers; Potassium Channel Blockers; Action po- tential, Local Anesthetics, Gillespie simulation. I. INTRODUCTION responsible for generating action potentials which are the basic requirement of cell to cell communication and sig- arXiv:1810.00701v1 [physics.bio-ph] 1 Oct 2018 nal propagation in nervous system[2]. Among the var- Ion channels are typically very complex transmem- ious cation, anion and neutral ion channels in general brane proteins[1] exhibiting a high degree of both struc- sodium and potassium channels are mainly responsible tural and functional diversity[2]. They have very distinct for action potential generation and its termination in electrical potential dependent gating mechanism where most of the neurons. There exists certain compounds the protein structures can adopt several conformational or molecules, typically called as ion channel blockers states such as closed, open and inactive states in which which selectively bind to specific protein conformations states they can either transport ions across the mem- of the ion channels and regulate their gating mechanism brane creating pores or restrict the ion permeability by by blocking the passage of ions across the membrane. closing the pore when required[3]. The ion channels are The high degree of specificity of these channel blockers on certain channels make them a valuable tool to treat numerous neural disorders[4–6]. Channel blockers of dif- ∗Electronic address: +[email protected],*[email protected] ferent types such as cationic blockers, anionic blockers, 2 amino acids, and other chemicals regulate the functional small, the effect of internal noise become more and more properties of the channel or prevent them to respond nor- important[18, 21–23]. Channel number fluctuations can mally. The naturally occurring sodium channel selective itself cause spontaneous spiking activity[23] without any blocker TTX[7] was known since 1964, followed by Sax- stimulus and thereby the stochastic version of Hodgkin- itoxin (STX), Neosaxitoxin (NSTX). Local anesthetics Huxley model is proposed[18, 20, 23]. Therefore, the ef- such as Lidocain, Phenytoin, Amiloride, Bupivacaine and fect of channel noise can not be ignored also in the study Tetracaine etc are clinically used[8, 9] for sodium channel of drug binding. Although innumerable number of exper- blocking. On the other hand potassium channel blockers, iments have been performed to study the channel block- such as 4-aminopyridine and 3, 4-diaminopyridine etc are ing phenomena, the theoretical investigation and com- used as local anaesthetics[10] and Tetraethylammonium prehensive understanding of the state specific blocking (TEA) is used only for experimental purpose[6, 11]. A still requires detailed investigation from which the corre- sodium blocker that blocks the open pore of the channel sponding physiological consequences of drug blocking can are called open state blocker[8, 12, 13]. Also there exists be deduced that too from the level of action potential de- closed state blocker, inactive state stabilizing blockers scription itself. From the vast number of studies already etc[8, 12, 13]. Potassium channel blockers either bind made it is very difficult to comprehend the mechanistic with the selectivity filter or at the central cavity(open link between the site specific blocking and the biophysical state) of the channel[14]. From the discovery of TTX consequences of that. So a microscopic reverse analysis toxicity[7] the present knowledge of ion channel block- is necessary from theoretical perspective to validate the ers have been a very vast, complicated and rigorous known observations. journey which was mostly developed for pharmacologi- Inspired by the work done by Schmid, Goychuk and cal interests[4–6, 9, 11]. H¨anggi[24] on the effect of ion channel blockers using Langevin type of stochastic description, we present here Recently Markov modelling and computer simula- a simple yet detailed and realistic approach for chan- tions prove to be promising techniques which provide nel blocking kinetics using the standard Markovian squid new insights into the fundamental principle and mech- axon model of Hodgkin-Huxley[17] which represents a anism about how these drugs alter the normal biologi- basic neuron cell type. We have incorporated new drug cal process[5, 15]. Single channel Markov models with blocked states in the original Hodgkin-Huxley model to discrete protein conformational states can simulate state theoretically extend it to a drug binding model. We then specific channel properties and their alterations by mu- have suitably adapted Gillespie[25] stochastic simulation tations, disease or drug binding[15]. As Markov mod- algorithm to study the effect of ion channel blockers on els can be developed both at the level of single chan- action potential, ionic current, spiking frequency, action nel activity and at the macroscopic state, they pro- potential duration and gating dynamics in moderately vide an implicit relationship between the single channel stochastic limit of channel noise. We have proposed recording and the macroscopic current[16]. Now the fa- three types of drug blocking mechanisms, e.g, sodium mous neuronal action potential model of Hodgkin and channel only blocking, potassium channel only blocking Huxley[17] is based on an assumption that the ion per- and a dual type blocking scheme with comparative bind- meation process can be approximated as both continu- ing affinity to these two channels. We have shown here ous and deterministic[17] as it considers infinite number through a few experimental evidences[26–38] that our ap- of ion channels inside a neuron. However, the perme- proach successfully correlates to the known physiological ation process existing within active membrane is now effects of ion channel blocking. known to be neither continuous nor deterministic. Ac- tive membrane consists of finite number of ion channels The Layout of the paper is as follows. In Section (II) which undergo random fluctuations between open and we have discussed the kinetic scheme of drug blocking and closed states[2] which scale inversely proportional to the the suitable modification of the simulation algorithm for number of channels present in a particular patch of a the present purpose. Effect of channel noise on the action neuron. Recent works also reveal that fluctuation in the potential has been discussed also under this section. In states of these channels are physiologically important in Section (III) the sodium and potassium channel blocking small neuronal structures[18–21]. When the number of effects on the trains of action potential have been dis- ion