Effect of Interactions on Molecular Fluxes and Fluctuations in the Transport Across Membrane Channels

Transport of molecules across membrane channels is investigated theoretically using exactly solvable one-dimensional discrete-state stochastic models. An interaction between molecules and membrane pores is modeled via a set of binding sites with different energies. It is shown that the interaction potential strongly influences the particle currents as well as fluctuations in the number of translocated molecules. For small concentration gradients the attractive sites lead to largest currents and fluctuations, while the repulsive interactions yield the largest fluxes and dispersions for large concentration gradients. Interaction energies that lead to maximal currents and maximal fluctuations are the same only for locally symmetric potentials, while they differ for the locally asymmetric potentials. The conditions for the most optimal translocation transport with maximal current and minimal dispersion are discussed. It is argued that in this case the interaction strength is independent of local symmetry of the potential of mean forces. In addition, the effect of the global asymmetry of the interaction potential is investigated, and it is shown that it also strongly affects the particle translocation dynamics. These phenomena can be explained by analyzing the details of the particle entering and leaving the binding sites in the channel.Comment: submitted to J. Chem. Phy

Sensing of Fluctuating Nanoscale Magnetic Fields Using NV Centres in Diamond

New magnetometry techniques based on Nitrogen-Vacancy (NV) defects in diamond allow for the imaging of static (DC) and oscillatory (AC) nanoscopic magnetic systems. However, these techniques require accurate knowledge and control of the sample dynamics, and are thus limited in their ability to image fields arising from rapidly fluctuating (FC) environments. We show here that FC fields place restrictions on the DC field sensitivity of an NV qubit magnetometer, and that by probing the dephasing rate of the qubit in a magnetic FC environment, we are able to measure fluctuation rates and RMS field strengths that would be otherwise inaccessible with the use of DC and AC magnetometry techniques. FC sensitivities are shown to be comparable to those of AC fields, whilst requiring no additional experimental overheads or control over the sample.Comment: 5 pages, 4 figure

New derivation for the equations of motion for particles in electromagnetism

We present equations of motion for charged particles using balanced equations, and without introducing explicitly divergent quantities. This derivation contains as particular cases some well known equations of motion, as the Lorentz-Dirac equations. An study of our main equations in terms of order of the interaction with the external field conduces us to the Landau-Lifshitz equations. We find that the analysis in second order show a special behavior. We give an explicit presentation up to third order of our main equations, and expressions for the calculation of general orders.Comment: 11 pages, 2 figures. Minor changes. Closer to published versio

Steady state existence of passive vector fields under the Kraichnan model

The steady state existence problem for Kraichnan advected passive vector models is considered for isotropic and anisotropic initial values in arbitrary dimension. The model includes the magnetohydrodynamic (MHD) equations, linear pressure model (LPM) and linearized Navier-Stokes (LNS) equations. In addition to reproducing the previously known results for the MHD and linear pressure model, we obtain the values of the Kraichnan model roughness parameter $\xi$ for which the LNS steady state exists.Comment: Improved text & figures, added references & other correction

Monte Carlo simulation for statistical mechanics model of ion channel cooperativity in cell membranes

Voltage-gated ion channels are key molecules for the generation and propagation of electrical signals in excitable cell membranes. The voltage-dependent switching of these channels between conducting and nonconducting states is a major factor in controlling the transmembrane voltage. In this study, a statistical mechanics model of these molecules has been discussed on the basis of a two-dimensional spin model. A new Hamiltonian and a new Monte Carlo simulation algorithm are introduced to simulate such a model. It was shown that the results well match the experimental data obtained from batrachotoxin-modified sodium channels in the squid giant axon using the cut-open axon technique.Comment: Paper has been revise

Kinetic models of ion transport through a nanopore

Kinetic equations for the stationary state distribution function of ions moving through narrow pores are solved for a number of one-dimensional models of single ion transport. Ions move through pores of length $L$, under the action of a constant external field and of a concentration gradient. The interaction of single ions with the confining pore surface and with water molecules inside the pore are modelled by a Fokker-Planck term in the kinetic equation, or by uncorrelated collisions with thermalizing centres distributed along the pore. The temporary binding of ions to polar residues lining the pore is modelled by stopping traps or energy barriers. Analytic expressions for the stationary ion current through the pore are derived for several versions of the model, as functions of key physical parameters. In all cases, saturation of the current at high fields is predicted. Such simple models, for which results are analytic, may prove useful in the study of the current/voltage relations of ion channels through membranes

Non-Markovian Stochastic Resonance: three state model of ion channel gating

Stochastic Resonance in single voltage-dependent ion channels is investigated within a three state non-Markovian modeling of the ion channel conformational dynamics. In contrast to a two-state description one assumes the presence of an additional closed state for the ion channel which mimics the manifold of voltage-independent closed subconformations (inactivated ``state''). The conformational transition into the open state occurs through a domain of voltage-dependent closed subconformations (closed ``state''). At distinct variance with a standard two-state or also three-state Markovian approach, the inactivated state is characterized by a broad, non-exponential probability distribution of corresponding residence times. The linear response to a periodic voltage signal is determined for arbitrary distributions of the channel's recovery times. Analytical results are obtained for the spectral amplification of the applied signal and the corresponding signal-to-noise ratio. Alternatively, these results are also derived by use of a corresponding two-state non-Markovian theory which is based on driven integral renewal equations [I. Goychuk and P. Hanggi, Phys. Rev. E 69, 021104 (2004)]. The non-Markovian features of stochastic resonance are studied for a power law distribution of the residence time-intervals in the inactivated state which exhibits a large variance. A comparison with the case of bi-exponentially distributed residence times possessing the same mean value, i.e. a simplest non-Markovian two-state description, is also presented

Toolbox for analyzing finite two-state trajectories

In many experiments, the aim is to deduce an underlying multi-substate on-off kinetic scheme (KS) from the statistical properties of a two-state trajectory. However, the mapping of a KS into a two-state trajectory leads to the loss of information about the KS, and so, in many cases, more than one KS can be associated with the data. We recently showed that the optimal way to solve this problem is to use canonical forms of reduced dimensions (RD). RD forms are on-off networks with connections only between substates of different states, where the connections can have non-exponential waiting time probability density functions (WT-PDFs). In theory, only a single RD form can be associated with the data. To utilize RD forms in the analysis of the data, a RD form should be associated with the data. Here, we give a toolbox for building a RD form from a finite two-state trajectory. The methods in the toolbox are based on known statistical methods in data analysis, combined with statistical methods and numerical algorithms designed specifically for the current problem. Our toolbox is self-contained - it builds a mechanism based only on the information it extracts from the data, and its implementation on the data is fast (analyzing a 10^6 cycle trajectory from a thirty-parameter mechanism takes a couple of hours on a PC with a 2.66 GHz processor). The toolbox is automated and is freely available for academic research upon electronic request

Ion-channel-like behavior in lipid bilayer membranes at the melting transition

It is well known that at the gel-liquid phase transition temperature a lipid bilayer membrane exhibits an increased ion permeability. We analyze the quantized currents in which the increased permeability presents itself. The open time histogram shows a "-3/2" power law which implies an open-closed transition rate that decreases like $k(t) \propto t^{-1}$ as time evolves. We propose a "pore freezing" model to explain the observations. We discuss how this model also leads to the $1/f^{\alpha}$ noise that is commonly observed in currents across biological and artificial membranes.Comment: 5 pages, 4 figure