1 Introduction
Ion channels are transmembrane protein structures that selectively allow given ion species to travel across the cell membrane. Zhou et al. [1] demonstrate that the channel protein transiently stabilize three K+ states, two within the selectivity filter and one within the water basket towards the intracellular side of the selectivity filter. Experimental evidence indicates that the selectivity filter is devoid of water molecules other than single water molecule between K ions [2]. By contrast, the basket is water molecule rich. We propose a Generalized Langevin Equation (GLE) that combines Markovian and non-Markovian processes at different time scales to describe the dynamics of ionic flow in K ion channels. This approach affords the development of an algorithm that describes oscillatory ionic diffusion along the selectivity filter. The resulting oscillatory functional behavior, with exponential decay, is obtained at the weak non-Markovian limit with two distinct time scales. They correspond to the processes of ionic diffusion and ionic drift. This oscillatory motion can be viewed as a time reversible stochastic process where the oscillatory time scale is very small compared to the diffusional time scale. The non-Markovian behavior, having a memory kernel associated with the frictional forces is due to the shear viscosity of water in the basket, considered as glue-like compared with ordinary water, due to the change in length scale [3]. Indeed, it is found in the recent experiments that the viscosity of water increases when the inner lining of nano-cavity is hydrophilic as well as hydrophobic [4]. However, the viscosity increases dramatically in case of hydrophilic surface. As a result of this spatial arrangements of water molecules, at the selectivity filter, ion movement may be considered a time reversible Markovian while at the water basket, it would behave in a non-Markovian regime.
Recent developments concerning quantum processes [5] demonstrate that there exists time reversible Markov process corresponding to Schrödinger's equation (at least in the non-relativistic domain) known as stochastic or Nelson's mechanics. Here, the diffusion constant has been shown to be related to Planck's constant. This framework has been extended [6–9] to address the dissipative forces i.e. velocity dependent forces in terms of quantum dynamics. Following Nelson's framework, a non-linear Schrödinger–Langevin equation has been derived for dissipative forces where the diffusion constant has been assumed to be related to Planck's constant. Recently, Bernroider and Roy [10] estimated the Lagrangian action at the level of a single ion channel using dimensional arguments. By considering empirically established conduction properties, the action is shown to be of the order of Planck's constant. This finding lead us to construct a non-linear Schrödinger–Langevin equation corresponding to this admixture of Markovian and non-Markovian processes for K ion channel dynamics. Moreover, the existence of a time reversible Markovian process associated with K ions in the selectivity filter at a relatively shorter time scale sheds new light on the applicability of quantum theory in the biological domain, at least in case of K ion channel.
The present Note will discuss: results on oscillatory ion flow in Section 2, stochastic quantization for dissipative forces in Section 3, the Schrödinger–Langevin equation approach to ion channel permeability in Section 4 and finally, conclusions and implications in Section 5.
2 Oscillatory ionic sequences
The generalized Langevin equation [11] for interacting ions is given by
(1) |
The Wiener process W is the 3-dimensional Gaussian process of which first moment is zero vector and second moment is a diagonal matrix whose element is minimum time between two Wiener processes. The memory kernel describes the generalized viscosity. Brownian motion triggers a Wiener process that distinguishes the SDE (Stochastic Differential Equation) from the ODE (Ordinary Differential Equation), so the term associated to the random force is referred to as the diffusion term.
We have integrated Eq. (1) numerically using the Euler scheme. In Eq. (1) the frictional force depends on the previous velocities through the integral over the kernel , which is quantified by the fluctuation dissipation theorem as
Here, we attempt to construct a classical description of the invariant measures of 2D Navier Stokes equation including the Stochastic effects. The existence of an invariant measure may be used to represent the asymptotic behavior of the system. If this invariant measure is unique, there is the possibility that the law of the process solution will converge to it. Therefore, when this holds true, this unique invariant measure describes the equilibrium to which the system tends. In fact, a unique invariant measure has been constructed by Li and Wagner [12] and the convergence has been shown when the 2D Navier–Stokes equations are perturbed by a white noise (not degenerate in space) with no limitations on phase space modes. In general, without constraints on the Reynolds number, the deterministic Navier–Stokes equations have many stationary solutions. However, while there is no information about the long time behavior of such deterministic processes, they can be viewed as invariant measures for the Navier–Stokes equations without the noise. Based on their findings, our result indicates that, given a sufficiently distributed random perturbation noise, just one invariant measure exists. The noise effect results in the generation of mixed system dynamics, allowing a unique asymptotic behavior.
At low flow rates, the diffusion term introduces the fluctuation into the ensemble averaged stress tensor, which appears as unwanted noise which severely limits our ability to calculate low flow rate viscosity, where the signal to noise ratio becomes very small. However, it can be reduced by variance reduction method. In the case of simple shear flow, the velocity field is time dependent, along with the dependence on shear rate and the fluctuating viscosity. At the inception of shear flow, the system is initially at equilibrium and the stress tensor vanishes. For time , a constant shear rate is applied and the stresses grow until they reach their steady state values, where the elongation rate is indeed time dependent. Note that, as we make no approximations here, in regards to the relative strength of the solvent molecules in comparison to the ions, and consequently no simplification can be made. The crux here is our assumption that the memory kernel can be written as
(2) |
Eq. (2) contains both Markovian and non-Markovian contributions that allows a continuous change from Markovian to non-Markovian dynamics and enables identification of both the terms. The non-Markovian process has two time scales whose contributions are dominated by the parameters respectively. It is clear from the form of kernel that in the limit of weak non Markovian process (),
(3) |
(4) |
Now we shall consider the existence of glue-like viscous properties [3,4] which are quite different from the water dynamics outside the channel. Doyle et al. [2] justified the existence of one water molecule only between the two K ions in the selectivity filter studying the strong electron density peak in the experimental map. Recently, Saparov and Pohl [13] claimed that the two ions in selectivity filter does not leave enough room to accommodate a vaporized water molecule (within the framework of water-vapor oscillation hypothesis). In fact, it is now known that the distance between two such water molecules exceeds the filter length [13]. The selectivity filter is about 12 Å long, whereas the remainder of the pore is wider and has a hydrophobic lining [2], contains water molecules. Thus, given that the selectivity filter is almost devoid of water and that the water in the basket is glue-like, we propose that such change in length scale at the interface gives rise to shear viscosity and hence to non-Markovian processes with memory. Further more, recent experiments suggest [4] that the viscosity of water becomes even more glue-like in hydrophilic nanoscale constraints, and that a hydrophobic surface resulting in a milder constraint [14]. This gives us the possibility to add a simple viscous term to our simulation. The most challenging aspect of simulation of ion channels has been the implementation of particle coupling and boundary conditions. In Brownian Dynamic (BD) simulations it is particularly difficult to maintain the correct particle concentrations and behavior that occur under extreme conditions in channel simulations. Here we use a simple algorithm in order to solve, the GLE.
We can use the second order Runge–Kutta method by discretizing the above GLE and approximating the potential as . We need n Gaussian random numbers to be picked at each step for the algorithm. Integrations of noise can be simulated by linear combinations of 3 normal Gaussian random numbers. The matrix elements of the coefficients are evaluated by the auto and cross-correlations of the noise integrations. We have used a set of algorithms that allow a comparison between the time scales of diffusion and shift performed along a selectivity filter in a model K-channel. The model describes two time scales for atomic ionic movement along the channel's length. This results in an oscillatory flow in the selectivity filter where the velocity autocorrelation exhibit an exponential decay. The Brownian simulation for GLE has been performed by subdividing the charge of the particle and increasing the number of particles by the same factor. We have obtained the velocity distributions corresponding to different limits of the processes and is explicit that oscillations dominate at the limit of weak non-Markovian process.
The time evolution of the kernel makes it clear that a phase transition of the shear viscosity within the channels may give an important clue as regards the rise of memory effects and the evolution of flip flop motion within ionic channels. It is apparent from the above analysis that using different time scales results oscillatory behavior within selectivity filter. Since this can occur at the weak non-Markovian limit (defined above), this oscillatory behavior can be described in terms of a time-reversible stochastic process. The diffusion coefficient for stochastic process is related to the action relevant to the scale of the occurrence of this process. It is important to estimate the action relevant at the scale of ionic diffusion.
2.1 Appropriate action at the level of ion channel
It is well known from Dimensional analysis that the mechanical action can be written as
3 Stochastic quantization and dissipative forces
Stochastic mechanics provides a deeper understanding to treat velocity dependent forces in terms of quantum theory. Let us consider that the position of a particle executes a stochastic process and the quantization rule is given by
(5) |
(6) |
(7) |
(8) |
Using the following complex function
(9) |
(10) |
(11) |
The function is related to the density . Using this approach Nelson proved that there exists a time reversible Markov process corresponding to Schrödinger's equation by assuming the diffusion constant as related to Planck's constant. This framework has been extended to incorporate the velocity dependent forces and hence for dissipative forces as mentioned be various authors [6–9]. From a physical point of view, this means that if the dynamic system interacts with its chaotic (thermal) environment, it is possible to extend the above framework as introduced by Nelson. In classical paradigm, the dynamical behavior of an open system can be described by Langevin equation
(12) |
The quantum mechanical behavior of an open system can be derived by the quantum mechanical version of Langevin equation. Since there is no well defined procedure to construct a Lagrangian or Hamiltonian for any open system in general, it will be more attractive to use stochastic quantization as proposed by Nelson. Here, there is no need to start from Lagrangian or Hamiltonian but one can begin with Newtonian equation in a more generalized sense. The first basic assumption of stochastic quantization is the following: The quantum mechanical behavior of the coordinate variable for open dynamical system can be represented by an n-dimensional diffusion process described by the following stochastic differential equation.
(13) |
(14) |
The second basic assumption is that the Langevin equation can be expressed in terms of mean velocity and mean acceleration. After some algebraic calculations, one gets quantum mechanical wave equation for the open system:
(15) |
Now we shall consider the GLE and the quantum mechanical wave equation for the velocity dependent forces considered in GLE.
4 Potassium channel and Schrödinger–Langevin equation
The Generalized Langevin Equation (GLE) used for the description of K ion channel is written as
(16) |
In the weak non-Markovian limit i.e.
(17) |
Following the rule of stochastic quantization as described in the previous section, we get
(18) |
Let us discuss the different non-linear terms due to various physical situations. The first nonlinear term is associated with Markov process taking delta function in the memory kernel. The velocity correlation is given by
The first part of the second non-linear term contains a characteristic time scale . Recent experimental observations [4] showed that water exhibits very different properties when it is confined to nano channels. Thus the velocity may be taken as almost constant throughout the integration and this part of the 2nd term can be written as
(19) |
As the second part remains unchanged, the quantum mechanical wave equation can be written as
(20) |
It seems evident, from channels characteristics, that the weak non-Markovian approximation may be valid for the selectivity filter. In this approximation and hence the contribution from the second non-linear term will be negligible leading to oscillatory ionic sequence.
5 Possible implications
The oscillatory ionic dynamics in Kion channels is proposed to occur at the limit of the weak non-Markovian approximation associated with a time reversible Markov process, at the selectivity filter. This reversible stochastic process belongs to a different time scale to that governing diffusion across the rest of the channel, which is determined by the glue-like properties of water at the water basket. The framework of stochastic mechanics provides a model for such dissipative force in terms of quantum theory. That channel ionic permeation can be associated with non-linear Schrödinger equation which addresses the issue of decoherence and time scale considerations. At this point it is worth mentioning that the Nelson process (the time reversible Markov process associated with Schrödinger equation) can, indeed, be considered an intermediate between quantum and classical time reversible processes, at least in the realm of ionic channel permeation.