Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.     

Estimating the dielectric constant of the channel protein and pore

When modelling biological ion channels using Brownian dynamics (BD) or Poisson–Nernst–Planck theory, the force encountered by permeant ions is calculated by solving Poisson’s equation. Two free parameters needed to solve this equation are the dielectric constant of water in the pore and the dielectric constant of the protein forming the channel. Although these values can in theory be deduced by various methods, they do not give a reliable answer when applied to channel-like geometries that contain charged particles. To determine the appropriate values of the dielectric constants, here we solve the inverse problem. Given the structure of the MthK channel, we attempt to determine the values of the protein and pore dielectric constants that minimize the discrepancies between the experimentally-determined current–voltage curve and the curve obtained from BD simulations. Two different methods have been applied to determine these values. First, we use all possible pairs of the pore dielectric constant of water, ranging from 20 to 80 in steps of 10, and the protein dielectric constant of 2–10 in steps of 2, and compare the simulated results with the experimental values. We find that the best agreement is obtained with experiment when a protein dielectric constant of 2 and a pore water dielectric constant of 60 is used. Second, we employ a learning-based stochastic optimization algorithm to pick out the optimum combination of the two dielectric constants. From the algorithm we obtain an optimum value of 2 for the protein dielectric constant and 64 for the pore dielectric constant.     

Integrated allosteric model of voltage gating of HCN channels

Hyperpolarization-activated (pacemaker) channels are dually gated by negative voltage and intracellular cAMP. Kinetics of native cardiac f-channels are not compatible with HH gating, and require closed/open multistate models. We verified that members of the HCN channel family (mHCN1, hHCN2, hHCN4) also have properties not complying with HH gating, such as sigmoidal activation and deactivation, activation deviating from fixed power of an exponential, removal of activation "delay" by preconditioning hyperpolarization. Previous work on native channels has indicated that the shifting action of cAMP on the open probability (Po) curve can be accounted for by an allosteric model, whereby cAMP binds more favorably to open than closed channels. We therefore asked whether not only cAMP-dependent, but also voltage-dependent gating of hyperpolarization-activated channels could be explained by an allosteric model. We hypothesized that HCN channels are tetramers and that each subunit comprises a voltage sensor moving between "reluctant" and "willing" states, whereas voltage sensors are independently gated by voltage, channel closed/open transitions occur allosterically. These hypotheses led to a multistate scheme comprising five open and five closed channel states. We estimated model rate constants by fitting first activation delay curves and single exponential time constant curves, and then individual activation/deactivation traces. By simply using different sets of rate constants, the model accounts for qualitative and quantitative aspects of voltage gating of all three HCN isoforms investigated, and allows an interpretation of the different kinetic properties of different isoforms. For example, faster kinetics of HCN1 relative to HCN2/HCN4 are attributable to higher HCN1 voltage sensors' rates and looser voltage-independent interactions between subunits in closed/open transitions. It also accounts for experimental evidence that reduction of sensors' positive charge leads to negative voltage shifts of Po curve, with little change of curve slope. HCN voltage gating thus involves two processes: voltage sensor gating and allosteric opening/closing.     

A cysteine scan of the inner vestibule of cyclic nucleotide-gated channels reveals architecture and rearrangement of the pore

Cyclic nucleotide-gated (CNG) channels belong to the P-loop-containing family of ion channels that also includes KcsA, MthK, and Shaker channels. In this study, we investigated the structure and rearrangement of the CNGA1 channel pore using cysteine mutations and cysteine-specific modification. We constructed 16 mutant channels, each one containing a cysteine mutation at one of the positions between 384 and 399 in the S6 region of the pore. By measuring currents activated by saturating concentrations of the full agonist cGMP and the partial agonists cIMP and cAMP, we show that mutating S6 residues to cysteine caused both favorable and unfavorable changes in the free energy of channel opening. The time course of cysteine modification with 2-aminoethylmethane thiosulfonate hydrochloride (MTSEA) was complex. For many positions we observed decreases in current activated by cGMP and concomitant increases in current activated by cIMP and cAMP. A model where modification affected both gating and permeation successfully reproduced the complex time course of modification for most of the mutant channels. From the model fits to the time course of modification for each mutant channel, we quantified the following: (a) the bimolecular rate constant of modification in the open state, (b) the change in conductance, and (c) the change in the free energy of channel opening for modification of each cysteine. At many S6 cysteines, modification by MTSEA caused a decrease in conductance and a favorable change in the free energy of channel opening. Our results are interpreted within the structural framework of the known structures of KcsA and MthK. We conclude that: (a) MTSEA modification affects both gating and permeation, (b) the open configuration of the pore of CNGA1 channels is consistent with the structure of MthK, and (c) the modification of S6 residues disrupts the helical packing of the closed channel, making it easier for channels to open.     

A comparative tool for the validity of rate kinetics in ion channels by Onsager reciprocity theorem

The precise form of the rate constant functions of ion channels is very crucial for reproducing the electrophysiological behavior. Therefore, how well they account for experimental data plays an important role in the behavior of the model. In this study, we derive kinetic coefficients of activation and inactivation gates in ion channels by Onsager reciprocity theorem for an ensemble of gating particles, and propose that the obtained kinetic coefficients can be used as a comparative tool for the empirical validity of fitted rate constant functions to experimental data. We also illustrate its applicability based on the activation and inactivation kinetics of T-type calcium channel in thalamic relay neurons. We show that the shape of the steady-state curve by itself seems to be a poor indicator of the functional form of the rate functions, but the time constant curves reflect considerable variation depending on the particular form of the rate functions, and that the kinetic coefficients related to the time constants provide a powerful tool to determine the empirical validity of the fitted rate constants.     

Statistical properties of ion channel records. Part II: estimation from the macroscopic current

Macroscopic ion channel current is the summation of the stochastic records of individual channel currents and therefore relates to their statistical properties. As a consequence of this relationship, it may be possible to derive certain statistical properties of single channel records or even generate some estimates of the records themselves from the macroscopic current when the direct measurement of single channel currents is not applicable. We present a procedure for generating the single channel records of an ion channel from its macroscopic current when the stochastic process of channel gating has the following two properties: (I) the open duration is independent of the time of opening event and has a single exponential probability density function (pdf), (II) all the channels have the same probability to open at time t. The application of this procedure is considered for cases where direct measurement of single channel records is difficult or impossible. First, the probability density function (pdf) of opening events, a statistical property of single channel records, is derived from the normalized macroscopic current and mean channel open duration. Second, it is shown that under the conditions (I) and (II), a non-stationary Markov model can represent the stochastic process of channel gating. Third, the non-stationary Markov model is calibrated using the results of the first step. The non-stationary formulation increases the model ability to generate a variety of different single channel records compared to common stationary Markov models. The model is then used to generate single channel records and to obtain other statistical properties of the records. Experimental single channel records of inactivating BK potassium channels are used to evaluate how accurately this procedure reconstructs measured single channel sweeps.     

Kinetics of 9-aminoacridine block of single Na channels

The kinetics of 9-aminoacridine (9-AA) block of single Na channels in neuroblastoma N1E-115 cells were studied using the gigohm seal, patch clamp technique, under the condition in which the Na current inactivation had been eliminated by treatment with N-bromoacetamide (NBA). Following NBA treatment, the current flowing through individual Na channels was manifested by square-wave open events lasting from several to tens of milliseconds. When 9-AA was applied to the cytoplasmic face of Na channels at concentrations ranging from 30 to 100 microM, it caused repetitive rapid transitions (flickering) between open and blocked states within single openings of Na channels, without affecting the amplitude of the single channel current. The histograms for the duration of blocked states and the histograms for the duration of open states could be fitted with a single-exponential function. The mean open time (tau o) became shorter as the drug concentration was increased, while the mean blocked time (tau b) was concentration independent. The association (blocking) rate constant, kappa, calculated from the slope of the curve relating the reciprocal mean open time to 9-AA concentration, showed little voltage dependence, the rate constant being on the order of 1 X 10(7) M-1s-1. The dissociation (unblocking) rate constant, l, calculated from the mean blocked time, was strongly voltage dependent, the mean rate constant being 214 s-1 at 0 mV and becoming larger as the membrane being hyperpolarized. The voltage dependence suggests that a first-order blocking site is located at least 63% of the way through the membrane field from the cytoplasmic surface. The equilibrium dissociation constant for 9-AA to block the Na channel, defined by the relation of l/kappa, was calculated to be 21 microM at 0 mV. Both tau -1o and tau -1b had a Q10 of 1.3, which suggests that binding reaction was diffusion controlled. The burst time in the presence of 9-AA, which is the sum of open times and blocked times, was longer than the lifetime of open channels in the absence of drug. All of the features of 9-AA block of single Na channels are compatible with the sequential model in which 9-AA molecules block open Na channels, and the blocked channels could not close until 9-AA molecules had left the blocking site in the channels.     

Critical assessment of a proposed model of Shaker

Detailed three-dimensional structures at atomic resolution are essential to understand how voltage-activated K(+) channels function. The X-ray crystallographic structure of the KvAP channel has offered the first view at atomic resolution of the molecular architecture of a voltage-activated K(+) channel. In the crystal, the voltage sensors are bound by monoclonal Fab fragments, which apparently induce a non-native conformation of the tetrameric channel. Thus, despite this significant advance our knowledge of the native conformation of a Kv channel in a membrane remains incomplete. Numerous results from different experimental approaches provide very specific constraints on the structure of K(+) channels in functional conformations. These results can be used to go further in trying to picture the native conformation of voltage-gated K(+) channels. However, the direct translation of all the available information into three-dimensional models is not straightforward and many questions about the structure of voltage-activated K(+) channels are still unanswered. Our aim in this review is to summarize the most important pieces of information currently available and to provide a critical assessment of the model of Shaker recently proposed by Lainé et al.     

Spike latency and jitter of neuronal membrane patches with stochastic Hodgkin-Huxley channels

Ion channel stochasticity can influence the voltage dynamics of neuronal membrane, with stronger effects for smaller patches of membrane because of the correspondingly smaller number of channels. We examine this question with respect to first spike statistics in response to a periodic input of membrane patches including stochastic Hodgkin-Huxley channels, comparing these responses to spontaneous firing. Without noise, firing threshold of the model depends on frequency—a sinusoidal stimulus is subthreshold for low and high frequencies and suprathreshold for intermediate frequencies. When channel noise is added, a stimulus in the lower range of subthreshold frequencies can influence spike output, while high subthreshold frequencies remain subthreshold. Both input frequency and channel noise strength influence spike timing. Specifically, spike latency and jitter have distinct minima as a function of input frequency, showing a resonance like behavior. With either no input, or low frequency subthreshold input, or input in the low or high suprathreshold frequency range, channel noise reduces latency and jitter, with the strongest impact for the lowest input frequencies. In contrast, for an intermediate range of suprathreshold frequencies, where an optimal input gives a minimum latency, the noise effect reverses, and spike latency and jitter increase with channel noise. Thus, a resonant minimum of the spike response as a function of frequency becomes more pronounced with less noise. Spike latency and jitter also depend on the initial phase of the input, resulting in minimal latencies at an optimal phase, and depend on the membrane time constant, with a longer time constant broadening frequency tuning for minimal latency and jitter. Taken together, these results suggest how stochasticity of ion channels may influence spike timing and thus coding for neurons with functionally localized concentrations of channels, such as in “hot spots” of dendrites, spines or axons.     

Habitat change in braided flood plains (Tagliamento, NE-Italy)

1. Relative changes and age distribution of habitats were investigated in the active channel of a bar‐braided and an island‐braided reach of the Tagliamento River (NE‐Italy). Between September 1999 and January 2002, six habitat types were delineated with a differential Global Positioning System on five dates following floods of different magnitude. Overlay maps were employed to calculate age and relative change of habitats. We established exponential decay rates (k‐values) for islands and major aquatic habitats. 2. Relative changes of all aquatic habitats combined were up to 82% between survey dates in the bar‐braided flood plain, with a cumulative rate of 85% over the 2.5‐year period. Relative habitat changes in the island‐braided flood plain were lower with a cumulative change of almost 60% during the study period. In the bar‐braided flood plain significant exponential decay relationships were established for channels, alluvial channels, backwaters, and ponds. 3. Half‐lives were particularly short for backwaters and ponds. In the island‐braided reach, significant relationships existed for channels and alluvial channels. The half‐lives of channels and alluvial channels increased with the presence of vegetated islands. Relative habitat composition within the active corridor remained almost constant, supporting the applicability of the shifting mosaic steady state model to braided floodplain ecosystems. 4. Our results indicate that under natural conditions aquatic floodplain habitats can be highly dynamic over short time‐scales. Even small water level fluctuations (‘flow pulses’) can lead to major habitat changes with important consequences for the fauna and flora.     

The structure of the gramicidin A transmembrane channel

Gramicidin A is a linear polypeptide antibiotic that facilitates the diffusion of monovalent cations across lipid bilayer membranes by forming channels. It has been proposed that the conducting channel is a dimer which is in equilibrium with nonconducting monomers in the membrane. To directly test this model in several independent ways, we have prepared and purified a series of gramicidin C derivatives. All of these derivatives are fully active analogs of gramicidin A, and each derivative has a useful chromophore esterified to the phenolic hydroxyl of tyrosine #11. Simultaneous conductance and fluorescence measurements on planar lipid bi-layer membranes containing dansyl gramicidin C yielded four conclusions: (1) A plot of the logarithm of the membrane conductance versus the logarithm of the membrane fluorescence had a slope of 2.0 ± 0.3, over a concentration range for which nearly all the gramicidin was monomeric. Hence, the active channel is a dimer of the nonconducting species. (2) In a membrane in which nearly all of the gramicidin was dimeric, the number of channels was approximately equal to the number of dimers. Thus, most dimers are active channels and so it should be feasible to carry out spectroscopic studies of the conformation of the transmembrane channel. (3) The association constant for dimerization is more than 1,000-fold larger in a glycerolester membrane with 26 Å-hydrocarbon thickness than in a 47 Å-glycerolester membrane. The dimerization constant in a 48 Å-phosphatidyl choline membrane was 200 times larger than in a 47 Å-glycerolester membrane, showing that it depends on the type of lipid as well as on the thickness of the hydrocarbon core. (4) We were readily able to detect 10^?14 mole cm^?2 of dansyl gramicidin C in a bilayer membrane, which corresponds to 60 fluorescent molecules per square μm. The fluorescent techniques described here should be sufficiently sensitive for fluorescence studies of reconstituted gates and receptors in planar bilayer membranes. An alternative method of determining the number of molecules of gramicidin in the channel is to measure the fraction of hybrid channels present in a mixture of 2 chemically different gramicidins. The single-channel conductance of p-phenylazo-benzene-sulfonyl ester gramicidin C (PABS gramicidin C) was found to be 0.68 that of gramicidin A. In membranes containing a mixture of these 2 gramicidins, a hybrid channel was evident in addition to 2 pure channels. The hybrid channel conductance was 0.82 that of gramicidin A. Fluorescence energy transfer from dansyl gramicidin C to diethylamino-phenylazobenzene-sulfonyl ester gramicidin C (DPBS gramicidin C), provided an independent way to measure the fraction of hybrid channels on liposomes. For both techniques the fraction of hybrid channels was found to be 2ad where a² and d² were the fractions of the 2 kinds of pure channels. This result strongly supports a dimer channel and the hybrid data excludes the possibility of a tetramer channel. The study of hybrid species by conductance and fluorescence techniques should be generally useful in elucidating the subunit structure of oligomeric assemblies in membranes. The various models which have been proposed for the conformation of the gramicidin transmembrane channel are briefly discussed.     

Electrophysiological Properties of Ion Channels in Ascaris suum Tissue Incorporated into Planar Lipid Bilayers

Ion channels are important targets of anthelmintic agents. In this study, we identified 3 types of ion channels in Ascaris suum tissue incorporated into planar lipid bilayers using an electrophysiological technique. The most frequent channel was a large-conductance cation channel (209 pS), which accounted for 64.5% of channels incorporated (n=60). Its open-state probability (P_o) was ~0.3 in the voltage range of −60~+60 mV. A substate was observed at 55% of the main-state. The permeability ratio of Cl⁻ to K⁺ (P_Cl/P_K) was ~0.5 and P_Na/P_K was 0.81 in both states. Another type of cation channel was recorded in 7.5% of channels incorporated (n=7) and discriminated from the large-conductance cation channel by its smaller conductance (55.3 pS). Its Po was low at all voltages tested (~0.1). The third type was an anion channel recorded in 27.9% of channels incorporated (n=26). Its conductance was 39.0 pS and P_Cl/P_K was 8.6±0.8. P_o was ~1.0 at all tested potentials. In summary, we identified 2 types of cation and 1 type of anion channels in Ascaris suum. Gating of these channels did not much vary with voltage and their ionic selectivity is rather low. Their molecular nature, functions, and potentials as anthelmintic drug targets remain to be studied further.     

Entropy effects on the ion-diffusion rate in transmembrane protein channels

We treat the transport of univalent cations through pore-like protein channels in biological membranes analytically, using two models (A + B) for the channel and the ion-channel interaction. A Lennard-Jones-type repulsion between the ions and the pore wall is introduced. We also include Van der Waals- and coulomb-type interactions between polar ligands of the pore-forming protein (e.g., carbonyl groups directed towards the axis of the channel) and the migrating particles. In model A, the polar groups are assumed to occur in pairs of dipoles pointing in opposite directions (as in the gramicidin A channel), while in model B the channel is treated as a pore with a radially isotropic charge distribution. In both models the ion-channel interaction leads to the occurrence of periodic potentials, corresponding to quasi-equilibrium and transition state sites of the ion in the pore. The diffusion rate can be calculated employing rate-theoretical concepts on the basis of microscopic parameters. It is demonstrated that the anomaly (inversion of the normal mass effect) for the transport rates of different ions can be related to differences in the activation entropy. The latter quantity is estimated analytically for both models. As a test, we performed numerical calculations with parameters based on the gramicidin A model. The results are in good agreement with experimental data and data from computer simulations. This shows that simple analytic expressions are well suited for predicting trends in the ionic conductivity of protein channels on the basis of microscopic interactions.     

The gating of single calcium-dependent potassium channels is described by an activation/blockade mechanism

Single calcium dependent potassium channels from cultured rat myoballs have been studied with the patch clamp technique, and current records subjected to statistical analysis. From the dependence of the mean open state probability on the internal calcium concentration, two calcium ions are required to open the channel. The open state and closed state lifetime distributions reveal that the usual activation model is not applicable to these channels. They are consistent with a two step gating mechanism that involves both activation by calcium and blockade by a calcium-sensitive gate.     

Does VDAC insert into membranes in random orientation?

It is widely accepted that voltage-dependent anion-selective channel (VDAC) inserts into planar lipid bilayers in a random orientation. This is in contrast to the well-documented oriented insertion of various channel-forming proteins. Because of the potential importance of this issue, we have examined the orientation of VDAC inserted in membranes. The time constants of the VDAC-current relaxation in response to applied positive and negative voltage pulses were used to characterize the channel orientation.We have found that VDAC channels can be separated into two groups according to differences in the time constant ratio. The difference in time constant ratio between the two main groups of VDAC channels was quantitative, and not qualitative as would be expected for opposite topologies. This finding allows us to hypothesize that both groups of VDAC channels possess a qualitatively similar asymmetry with respect to the localization of voltage-gated domains and, consequently, with respect to its entire molecular structure. The probability of having each type of VDAC channel conformation is predetermined by the protein structure in aqueous solution.A striking resemblance between asymmetry in voltage sensitivity at the single-channel and multi-channel levels was also demonstrated. The first inserted channel seems to direct subsequent insertions of channels with a similar conformation.     

A microscopic electrostatic model for the amphotericin B channel.

A microscopic model of an amphotericin B channel is proposed. The structure of the pores is generated using the atomic coordinates of the molecule in the structure determined experimentally by X-ray diffraction. The net charges of the atoms are determined by Mulliken analysis. With these charges the electrostatic energy profiles are calculated for a monovalent ion passing through the channels formed by different number of antibiotic molecules having different radii. The water inside the channel was considered through a continuum medium using the dielectric constant of the bulk, and the membrane contribution was included using the virtual images of the pore in a dielectric slab of epsilon = 3. The model satisfactorily explains the permeability and selectivity characteristics as well as other observations yet unexplained. The electrostatic profiles obtained reinforce the hypothesis of the existence of channels formed by a variable number of units.     

Anomalous effects of external TEA on permeation and gating of the A-type potassium current in H. aspersa neuronal somata

The model proposed for external TEA block of Shaker K+ channels predicts a proportional relationship between TEA sensitivity and calculated electrical distance derived from measurements of voltage dependence of TEA block. In the present study, we examined this relationship for the A-type K+ current (IA) of Helix aspersa in neuronal somata using the whole-cell patch-clamp technique. External TEA inhibited IA with strong voltage dependence, such that the TEA dissociation constant was increased at depolarized test potentials. The half-inhibition constant (V0.5) for TEA block was approximately 21 mM at 0 mV, and V0.5 increased to approximately 67 mM at 50 mV. The calculated electrical distance for TEA block suggested that TEA traversed 65% of the way into the membrane electrical field. TEA also caused significant shifts in the voltage-dependence of A-type K+ channel gating. For example, at TEA concentrations below that required to fully suppress delayed outward currents, TEA caused depolarizing shifts in the voltage-dependence of A-type channel activation, steady-state inactivation, time for removal of inactivation, and slowed channel activation kinetics. Taken together, these observations suggest that TEA biased the local field potential near voltage-sensing domains of A-type K+ channels, causing the transmembrane electrical field to be relatively hyperpolarized in the presence of TEA. In summary, the calculated electrical distance of TEA block of A-type K+ channels in H. aspersa neurons is unprecedented among other K+ channels. This raises concerns about the conventional interpretation of this value. Furthermore, the voltage-dependent properties of IA are modified by TEA at concentrations previously used to isolate delayed rectifier potassium channels (IKDR) selectively. This lack of specificity has important implications for recent, as well as future studies of IA in H. aspersa and possibly other snail neurons.     

Molecular cloning and functional characterization of a novel delayed rectifier potassium channel from channel catfish (Ictalurus punctatus): expression in taste buds

The gustatory system of channel catfish is widely studied for its sensitivity to amino acids. As a first step in identifying the molecular components that play a role in taste transduction in catfish, we cloned the full-length cDNA for Kv2-catfish, a novel K(+) channel that is expressed in taste buds. The deduced amino acid sequence is 816 residues, and shares a 56-59% sequence identity with Kv2.1 and Kv2.2, the other members of the vertebrate Kv2 subfamily of voltage-gated K(+) channels. The Kv2-catfish RNA was expressed in taste buds, brain, skeletal muscle, kidney, intestine and gills, and its gene is represented as a single copy in the catfish genome. Recombinant channels expressed in XENOPUS: oocytes were selective for K(+), and were inhibited by tetraethylammonium applied to the extracellular side of the membrane during two-electrode voltage clamp analysis with a 50% inhibitory constant of 6.1 mM. The channels showed voltage-dependent activation, and did not inactivate within 200 ms. Functionally, Kv2-catfish is a voltage-gated, delayed rectifier K(+) channel, and its primary structure is the most divergent sequence identified among the vertebrate members of the Kv2 subfamily of K(+) channels, being related equally well to Kv2.1 and Kv2.2.     

Patch-clamp profile of ion channels in resting murine B lymphocytes

Summary Patch-clamp studies of single ion channel currents in freshly isolated murine B lymphocytes are characterized here according to their respective unitary conductances, ion selectivities, regulatory factors, distributions and kinetic behavior. The most prevalent ion channel in murine B lymphocytes is a large conductance (348 pS) nonselective anion channel. This report characterizes additional conductances including: two chloride channels (40 and 128 pS), a calcium-activated potassium channel (93 pS), and an outwardly rectifying potassium channel which displays two distinct conductances (18 and 30 pS). Like the anion channel, both chloride channels exhibit little activity in the cellattached patch configuration. The kinetic behavior of all of these channels is complex, with variable periods of bursting and flickering activity interspersed between prolonged closed/open intervals (dwell times). It is likely that some of these channels play an important role in the signal transduction of B cell activation.     

Comparison of the gramicidin a potassium/sodium permeability and single channel conductance ratio

If the ion concentration is low enough that most channels are unoccupied, then the ‘independence relations’ should be satisfied and the permeability ratio should equal the conductance ratio. It has been previously reported that for the gramicidin A channel these ratios for Na⁺ and K⁺ were not equal at concentrations as low as 10 mM. However, these ratios were not measured at the same applied potential, as is required by the theory. Instead, the conductance ratio was measured at 100 mV and corrected using calculated current-voltage relations. In this report the comparison between permeability and conductance ratios is reexamined using data obtained at the correct potential. There is no significant difference in the ratios at 10 mM when they are measured at the same voltage. This implies that most channels are not occupied by sodium or potassium ions at 10 mM.