-
On the selectivity of KcsA potassium channel: asymptotic analysis and computation
Authors:
Zilong Song,
Xiulei Cao,
Tzyy-Leng Horng,
Huaxiong Huang
Abstract:
Potassium (K$^+$) channels regulate the flux of K$^+$ ions through cell membranes and plays significant roles in many physiological functions. This work studies the KcsA potassium channel, including the selectivity and current-voltage (IV) relations. A modified Poisson-Nernst-Planck system is employed, which include the size effect by Bikerman model and solvation energy by Born model. The selectiv…
▽ More
Potassium (K$^+$) channels regulate the flux of K$^+$ ions through cell membranes and plays significant roles in many physiological functions. This work studies the KcsA potassium channel, including the selectivity and current-voltage (IV) relations. A modified Poisson-Nernst-Planck system is employed, which include the size effect by Bikerman model and solvation energy by Born model. The selectivity of KcsA for various ions (K$^+$, Na$^+$, Cl$^-$, Ca$^{2+}$ and Ba$^{2+}$) is studied analytically, and the profiles of concentrations and electric potential are provided. The selectivity is mainly influenced by permanent negative charges in filter of channel and the ion sizes. K$^+$ is always selected compared with Na$^+$ (or Cl$^-$), as smaller ion size of Na$^+$ causes larger solvation energy. There is a transition for selectivity among K$^+$ and divalent ions (Ca$^{2+}$ and Ba$^{2+}$), when negative charge in filter exceeds a critical value determined by ion size. This explains why divalent ions can block the KcsA channel. The profiles and IV relations are studied by analytical, numerical and hybrid methods, and are cross-validated. The results show the selectivity of the channel and also the saturation of IV curve. A simple strategy is given to compute IV relations analytically, as first approximation. The numerical method deals with general structure or parameters, but the limitations and difficulties of pure numerical simulation are also pointed out. The hybrid method provides IV relations most effectively for comparison. The reason for saturation of IV relation is illustrated, and the IV curve shows agreement with the profile and scale of experimental results.
△ Less
Submitted 12 February, 2019;
originally announced February 2019.
-
Continuum gating current models computed with consistent interactions
Authors:
Tzyy-Leng Horng,
Robert S. Eisenberg,
Chun Liu,
Francisco Bezanilla
Abstract:
The action potential signal of nerve and muscle is produced by voltage sensitive channels that include a specialized device to sense voltage. Gating currents of the voltage sensor are now known to depend on the back-and-forth movements of positively charged arginines through the hydrophobic plug of a voltage sensor domain. Transient movements of these permanently charged arginines, caused by the c…
▽ More
The action potential signal of nerve and muscle is produced by voltage sensitive channels that include a specialized device to sense voltage. Gating currents of the voltage sensor are now known to depend on the back-and-forth movements of positively charged arginines through the hydrophobic plug of a voltage sensor domain. Transient movements of these permanently charged arginines, caused by the change of transmembrane potential, further drag the S4 segment and induce opening/closing of ion conduction pore by moving the S4-S5 linker. The ion conduction pore is a separate device from the voltage sensor, linked (in an unknown way) by the mechanical motion and electric field changes of the S4-S5 linker. This moving permanent charge induces capacitive current flow everywhere. Everything interacts with everything else in the voltage sensor so everything must interact with everything else in its mathematical model, as everything does in the whole protein. A PNP-steric model of arginines and a mechanical model for the S4 segment are combined using energy variational methods in which all movements of charge and mass satisfy conservation laws of current and mass. The resulting 1D continuum model is used to compute gating currents under a wide range of conditions, corresponding to experimental situations. Chemical-reaction-type models based on ordinary differential equations cannot capture such interactions with one set of parameters. Indeed, they may inadvertently violate conservation of current. Conservation of current is particularly important since small violations (<0.01%) quickly (<< 10-6 seconds) produce forces that destroy molecules. Our model reproduces signature properties of gating current: (1) equality of on and off charge in gating current (2) saturating voltage dependence in QV curve and (3) many (but not all) details of the shape of gating current as a function of voltage.
△ Less
Submitted 9 July, 2017;
originally announced July 2017.
-
Dynamics of slow light and light storage in a Doppler-broadened electromagnetically-induced-transparency medium: A numerical approach
Authors:
Shih-Wei Su,
Yi-Hsin Chen,
Shih-Chuan Gou,
Tzyy-Leng Horng,
Ite A. Yu
Abstract:
We present a numerical scheme to study the dynamics of slow light and light storage in an electromagneticallyinduced- transparency (EIT) medium at finite temperatures. Allowing for the motional coupling, we derive a set of coupled Schrödinger equations describing a boosted closed three-level EIT system according to the principle of Galilean relativity. The dynamics of a uniformly moving EIT medium…
▽ More
We present a numerical scheme to study the dynamics of slow light and light storage in an electromagneticallyinduced- transparency (EIT) medium at finite temperatures. Allowing for the motional coupling, we derive a set of coupled Schrödinger equations describing a boosted closed three-level EIT system according to the principle of Galilean relativity. The dynamics of a uniformly moving EIT medium can thus be determined by numerically integrating the coupled Schrödinger equations for atoms plus one ancillary Maxwell-Schrödinger equation for the probe pulse. The central idea of this work rests on the assumption that the loss of ground-state coherence at finite temperatures can be ascribed to the incoherent superposition of density matrices representing the EIT systems with various velocities. Close agreements are demonstrated in comparing the numerical results with the experimental data for both slow light and light storage. In particular, the distinct characters featuring the decay of ground-state coherence can be well verified for slow light and light storage. This warrants that the current scheme can be applied to determine the decaying profile of the ground-state coherence as well as the temperature of the EIT medium.
△ Less
Submitted 31 May, 2011;
originally announced May 2011.