Current noise around steady states in discrete transport systems

Subject of this paper is the transport noise in discrete systems. The transport systems are given by a number (n) of binding sites separated by energy barriers. These binding sites may be in contact with constant outer reservoirs. The state of the system is characterized by the occupation numbers of particles (current carriers) at these binding sites. The change in time of the occupation numbers is generated by individual “jumps” of particles over the energy barriers, building up the flux matter (for charged particles: the electric current). In the limit n → ∞ continuum processes as e.g. usual diffusion are included in the transport model. The fluctuations in occupation numbers and other quantities linearly coupled to the occupation numbers may be treated with the usual master equation approach. The treatment of the fluctuations in fluxes (current) makes necessary a different theoretical approach which is presented in this paper under the assumption of vanishing interactions between the particles. This approach may be applied to a number of different transport systems in biology and physics (ion transport through porous channels in membranes, carrier mediated ion transport through membranes, jump diffusion e.g. in superionic conductors). As in the master equation approach the calculation of correlations and noise spectra may be reduced to the solution of the macroscopic equations for the occupation numbers. This result may be regarded as a generalization to non-equilibrium current fluctuations of the usual Nyquist theorem relating the current (voltage) noise spectrum in thermal equilibrium to the macroscopic frequency dependent admittance.The validity of the general approach is demonstrated by the calculation of the autocorrelation function and spectrum of current noise for a number of special examples (e.g, pores in membrances, carrier mediated ion transport).     

Current fluctuations in discrete transport systems far from equilibrium. Breakdown of the fluctuation dissipation theorem

A recently developed theoretical approach to transport fluctuations around stable steady states in discrete biological transport systems is used in order to investigate general fluctuation properties at nonequilibrium. An expression for the complex frequency dependent admittance at nonequilibrium is derived by calculation of the linear current response of the transport systems to small disturbances in the applied external voltage. It is shown that the Nyquist or fluctuation dissipation theorem, by which at equilibrium the macroscopic admittance or linear response can be expressed in terms of fluctuation properties of the system, breaks down at nonequilibrium. The spectral density of current fluctuations is decomposed into one term containing the macroscopic admittance and a second term which is bilinear in current. This second term is generated by microscopic disturbances, which cannot be excited by external macroscopic perturbations. At special examples it is demonstrated that this second term is decisive for the occurrence of excess noise e.g. the 1/f(2)-Iorentzian noise generated by the opening and closing of nerve channels in biological membranes.     

Theory of single-file noise

A general theoretical approach to the analysis of electric fluctuations generated by the so-called single-file diffusion through narrow channels is presented. The formalism is a slight extension of an approach to electric fluctuations in discrete transport systems with negligible interactions between the particles recently developed by one of the authors. In the single-file transport mechanism interactions between the particles must be taken into account. Three main results of principal interest are: (a) the electric fluctuations around stationary states (at equilibrium and non-equilibrium) are determined by the time-dependent solutions of the macroscopic single-file transport equations, (b) as a direct consequence of the interactions between the ions in the single-file transport the macroscopic time-dependent current and the autocorrelation function of the microscopic current fluctuations can exhibit damped oscillatory behavior, and the current noise spectrum can show peaking, (c) the number of binding sites for the ions within the pores seems to have a strong influence on the oscillatory behavior: with increasing number of binding sites the damping of the oscillations decreases and the peaking of the spectrum becomes stronger.     

Theory of single-file noise.

A general theoretical approach to the analysis of electric fluctuations generated by the so-called single-file diffusion through narrow channels is presented. The formalism is a slight extension of an approach to electric fluctuations in discrete transport systems with negligible interactions between the particles recently developed by one of the authors. In the single-file transport mechanism interactions between the particles must be taken into account. Three main results of principal interest are: (a) the electric fluctuations around stationary states (at equilibrium and non-equilibrium) are determined by the time-dependent solutions of the macroscopic single-file transport equations, (b) as a direct consquence of the interactions between the ions in the single-file transport the macroscopic time-dependent current and the autocorrelation function of the microscopic current fluctuations can exhibit damped oscillatory behavior, and the current noise spectrum can show peaking, (c) the number of binding sites for the ions within the pores seems to have a strong influence on the oscillatory behavior: with increasing number of binding sites the damping of the oscillations decreases and the peaking of the spectrum becomes stronger.     

Rate theory models for ion transport through rigid pores. I. Time-dependent analysis in the case of vanishing interactions

In this first of a series of papers concerning the theoretical analysis of rate theory models for ion transport through rigid pores, the case of vanishing interactions is investigated. "Rigidity" means that ions crossing membranes through pores see a fixed structure of the pores, not changing in time. A single pore is considered to be a sequence of (n + 1) activation barriers separated by n energy minima. The explicit analytical treatment is restricted to pores with regular internal barrier structure, including the nonequilibrium situation of an applied electric field. In this case the connection with continuum diffusion models is demonstrated by performing in the limit n leads to infinity (n = number of binding sites within the pores) the transition to continuum. Thus, from diffusion equations describing a discrete number of jumps, the corresponding diffusion-like partial differential equations and boundary conditions are generated. For regular pores, from the time dependent solutions of the discrete equations, the corresponding solutions of the continuum equations are explicitly generated. The time-dependent relaxation behaviour of the discrete model is in good agreement with the continuum model if one assumes more than two binding sites in the pores.     

Nonequilibrium voltage fluctuations in biological membranes. I. General framework of charge transport in discrete systems and related voltage noise

This paper continues our work on the theory of nonequilibrium voltage noise generated by electric transport processes in membranes. Introducing the membrane voltage as a further variable, a system of kinetic equations linearized in voltage is derived by which generally the time-dependent behaviour of charge-transport processes under varying voltage can be discussed. Using these equations, the treatment of voltage noise can be based on the usual master equation approach to steady-state fluctuations of scalar quantities. Thus, a general theoretical approach to nonequilibrium voltage noise is presented, completing our approach to current fluctuations which had been developed some years ago. It is explicitly shown that at equilibrium the approach yields agreement with the Nyquist relation, while at nonequilibrium this relation is not valid. A further general property of voltage noise is the reduction of low-frequency noise with increasing number of transport units as a consequence of the interactions via the electric field. In a second paper, the approach will be applied for a number of special transport mechanisms, such as ionic channels, carriers or electrogenic pumps.     

Non-equilibrium voltage noise generated by ion transport through pores

In this paper, we describe a systematic approach to the theoretical analysis of non-equilibrium voltage noise that arises from ions moving through pores in membranes. We assume that an ion must cross one or two barriers in the pore in order to move from one side of the membrane to the other. In our analysis, we consider the following factors: a) surface charge as a variable in the kinetic equations, b) linearization of the kinetic equations, c) master equation approach to fluctuations. To analyze the voltage noise arising from ion movement through a two barrier (i.e., one binding site) pore, we included the effects of ions in the channel's interior on the voltage noise. The current clamp is considered as a white noise generating additional noise in the system. In contrast to what is found for current noise, at low frequencies the voltage noise intensity is reduced by increasing voltage across the membrane. With this approach, we demonstrate explicity for the examples treated that, apart from additional noise generated by the current clamp, the non-equilibrium voltage fluctuations can be related to the current fluctuations by the complex admittance.     

Pushing the limits of luminescence thermometry: probing the temperature of proteins in cells

Proteins are involved in numerous cellular activities such as transport and catalysis. Misfolding during biosynthesis and malfunctioning as a molecular machine may lead to physiological disorders and metabolic problems. Protein folding and mechanical work may be viewed as thermodynamic energetically favorable processes in which stochastic nonequilibrium intermediate states may be present with conditions such as thermal fluctuations. In my opinion, measuring those thermal fluctuations may be a way to access the energy exchange between the protein and the physiological environment and to better understand how those nonequilibrium states may influence the misfolding/folding process and the efficiency of the molecular engine cycle. Here, I discuss luminescence thermometry as a possible way to measure those temperature fluctuations from a single-molecule experimental perspective with its current technical limitations and challenges.

     

Fluctuation of surface charge in membrane pores

Surface charge in track-etched polyethylene terephthalate (PET) membranes with narrow pores has been probed with a fluorescent cationic dye (3,3'-diethyloxacarbocyanine iodide (diO-C2-(3))) using confocal microscopy. Staining of negatively charged PET membranes with diO-C2-(3) is a useful measure of surface charge for the following reasons: 1) the dye inhibits K(+) currents through the pores and reduces their selectivity for cations; 2) it inhibits [3H]-choline+ transport and promotes 36Cl- transport across the membrane in a pH- and ionic-strength-dependent fashion; and 3) staining of pores by diO-C2-(3) is reduced by low pH and by the presence of divalent cations such as Ca2+ and Zn2+. Measurement of the time dependence of cyanine staining of pores shows fluctuations of fluorescence intensity that occur on the same time scale as do fluctuations of ionic current in such pores. These data support our earlier proposal that fluctuations in ionic current across pores in synthetic and biological membranes reflect fluctuations in the surface charge of the pore walls in addition to molecular changes in pore proteins.     

A parallel path model forNecturus proximal tubule

Summary A parallel path model based on the principles of nonequilibrium thermodynamics was developed for theNecturus proximal tubule. The cellular path was represented as a luminal membrane followed by an irreversible active NaCl transport system in the peritubular barrier. The shunt pathway was described as three coarse barriers in series: tight junction, lateral intercellular spaces, and basement membrane with connective tissue. Volume and solute flows were predicted by the model equations as a function of applied electric current. Variations of the model parameters revealed the quantitative importance of the shunt path properties and the relative insensitivity of epithelial transport to changes in most cell parameters. Circulation of electric current and solute within the epithelium were shown to significantly influence the bahavior of the tubule in the presence of an electric field. Values for all transport parameters of the shunt path and epithelium were calculated and compared with available experimental evidence. Volume flow and electric currents predicted by the model compared favorably with experimental observations.     

Attenuation of current and voltage noise signals recorded from epithelia

Transepithelially recorded current and voltage fluctuations are filtered by the impedance of the electrical equivalent parameters of the preparation, in series or in parallel, with the noise source. Fluctuations in voltage and current are assumed to be caused by fluctuations in conductance of the apical membrane. In order to obtain an estimation of the intrinsic noise amplitudes, calculations are presented to correct the transepithelial fluctuations. The influence of different model parameters on the recorded noise spectra is investigated. It is shown that the shape of the transepithelially recorded noise spectra may differ from the intrinsic ones, e.g. “peaking” in the power spectra may be explained by the assumption of a positive (referred to cell inside) e.m.f. at the basolateral membrane or a polarization impedance in series with the epithelium. Furthermore it is demonstrated that the ratio of voltage to current noise power may differ from the squared magnitude of the impedance.     

Nonequilibrium voltage fluctuations in biological membranes. II. Voltage and current noise generated by ion carriers, channels and electrogenic pumps

As applications of the general theoretical framework of charge transport in biological membranes and related voltage and current noise, a number of model calculations are presented for ion carriers, rigid channels, channels with conformational substates and electrogenic pumps. The results are discussed with special reference to the problem of threshold values for sensory transduction processes and their limitations by voltage fluctuations. Furthermore, starting from the special results of model calculations, an attempt is made to determine more general aspects of electric fluctuations generated by charge-transport processes in biological membranes: different frequency dependences of voltage and current noise, and dependence of noise intensities with increasing distance from the equilibrium state.     

Location of a constriction in the lumen of a transmembrane pore by targeted covalent attachment of polymer molecules

Few methods exist for obtaining the internal dimensions of transmembrane pores for which 3-D structures are lacking or for showing that structures determined by crystallography reflect the internal dimensions of pores in lipid bilayers. Several approaches, involving polymer penetration and transport, have revealed limiting diameters for various pores. But, in general, these approaches do not indicate the locations of constrictions in the channel lumen. Here, we combine cysteine mutagenesis and chemical modification with sulfhydryl-reactive polymers to locate the constriction in the lumen of the staphylococcal alpha-hemolysin pore, a model protein of known structure. The rates of reaction of each of four polymeric reagents (MePEG-OPSS) of different masses towards individual single cysteine mutants, comprising a set with cysteines distributed over the length of the lumen of the pore, were determined by macroscopic current recording. The rates for the three larger polymers (1.8, 2.5, and 5.0 kD) were normalized with respect to the rates of reaction with a 1.0-kD polymer for each of the seven positions in the lumen. The rate of reaction of the 5.0-kD polymer dropped dramatically at the centrally located Cys-111 residue and positions distal to Cys-111, whether the reagent was applied from the trans or the cis side of the bilayer. This semi-quantitative analysis sufficed to demonstrate that a constriction is located at the midpoint of the pore lumen, as predicted by the crystal structure, and although the constriction allows a 2.5-kD polymer to pass, transport of a 5.0-kD molecule is greatly restricted. In addition, PEG chains gave greater reductions in pore conductance when covalently attached to the narrower regions of the lumen, permitting further definition of the interior of the pore. The procedures described here should be applicable to other pores and to related structures such as the vestibules of ion channels.     

Internal elastic lamina affects the distribution of macromolecules in the arterial wall: a computational study

The internal elastic lamina (IEL), which separates the arterial intima from the media, affects macromolecular transport across the medial layer. In the present study, we have developed a two-dimensional numerical simulation model to resolve the influence of the IEL on convective-diffusive transport of macromolecules in the media. The model considers interstitial flow in the medial layer that has a complex entrance condition because of the presence of leaky fenestral pores in the IEL. The IEL was modeled as an impermeable barrier to both water and solute except for the fenestral pores that were assumed to be uniformly distributed over the IEL. The media were modeled as a heterogeneous medium composed of an array of smooth muscle cells (SMCs) embedded in a continuous porous medium representing the interstitial proteoglycan and collagen fiber matrix. Results for ATP and low-density lipoprotein (LDL) demonstrate a range of interesting features of molecular transport and uptake in the media that are determined by considering the balance among convection, diffusion, and SMC surface reaction. The ATP concentration distribution depends strongly on the IEL pore structure because ATP fluid-phase transport is dominated by diffusion emanating from the fenestral pores. On the other hand, LDL fluid-phase transport is only weakly dependent on the IEL pore structure because convection spreads LDL laterally over very short distances in the media. In addition, we observe that transport of LDL to SMC surfaces is likely to be limited by the fluid phase (surface concentration less than bulk concentration), whereas ATP transport is limited by reaction on the SMC surface (surface concentration equals bulk concentration).     

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.     

Linear noise approximation method for calculating nonequilibrium fluctuations of the occupation numbers in radiative-collisional average ion models

A novel method is proposed for calculating nonequilibrium fluctuations of the mean occupation numbers of the electron shells in the radiative-collisional average-ion models of multicharged plasma kinetics. For the class of Slater ionic models, equations are derived for the mean occupation numbers of the electron shells and their fluctuations in the Fokker-Planck approximation. To calculate the fluctuations, the Fokker-Planck equation is linearized in the vicinity of the steady-state nonequilibrium solution to the kinetic equations (linear noise approximation). The method proposed allows one to take into account both the nonequilibrium correlations of the occupation-number fluctuations and the thermodynamically equilibrium statistical correlations related to the Coulomb interaction among bound electrons. The relation among the coefficients in the Fokker-Planck equation for the occupation-number fluctuations of the electron shells is discussed based on the fluctuation-dissipative theorem.     

Nonequilibrium response spectroscopy of voltage-sensitive ion channel gating.

We describe a new electrophysiological technique called nonequilibrium response spectroscopy, which involves application of rapidly fluctuating (as high as 14 kHz) large-amplitude voltage clamp waveforms to ion channels. As a consequence of the irreversible (in the sense of Carnot) exchange of energy between the fluctuating field and the channel protein, the gating response is exquisitely sensitive to features of the kinetics that are difficult or impossible to adequately resolve by means of traditional stepped potential protocols. Here we focus on the application of dichotomous (telegraph) noise voltage fluctuations, a broadband Markovian colored noise that fluctuates between two values. Because Markov kinetic models of channel gating can be embedded within higher-dimensional Markov models that take into account the effects of the voltage fluctuations, many features of the response of the channels can be calculated algebraically. This makes dichotomous noise and its generalizations uniquely suitable for model selection and kinetic analysis. Although we describe its application to macroscopic ionic current measurements, the nonequilibrium response method can also be applied to gating and single channel current recording techniques. We show how data from the human cardiac isoform (hH1a) of the Na+ channel expressed in mammalian cells can be acquired and analyzed, and how these data reveal hidden aspects of the molecular kinetics that are not revealed by conventional methods.