Modeling of ion permeation in calcium and sodium channel selectivity filters

On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation

Ion channel selectivity is essential for their function, yet the molecular basis of a channel's ability to select between ions is still rather controversial. In this work, using a combination of molecular dynamics simulations and electrophysiological current measurements we analyze the ability of the NaChBac channel to discriminate between calcium and sodium. Our simulations show that a single calcium ion can access the Selectivity Filter (SF) interacting so strongly with the glutamate ring so as to remain blocked inside. This is consistent with the tiny calcium currents recorded in our patch-clamp experiments. Two reasons explain this scenario. The first is the higher free energy of ion/SF binding of Ca²⁺ with respect to Na⁺. The second is the strong electrostatic repulsion exerted by the resident ion that turns back a second potentially incoming Ca²⁺, preventing the knock-on permeation mechanism. Finally, we analyzed the possibility of the Anomalous Mole Fraction Effect (AMFE), i.e. the ability of micromolar Ca²⁺ concentrations to block Na⁺ currents. Current measurements in Na⁺/Ca²⁺ mixed solutions excluded the AMFE, in agreement with metadynamics simulations showing the ability of a sodium ion to by-pass and partially displace the resident calcium. Our work supports a new scenario for Na⁺/Ca²⁺ selectivity in the bacterial sodium channel, challenging the traditional notion of an exclusion mechanism strictly confining Ca²⁺ ions outside the channel.

Energetics of discrete selectivity bands and mutation-induced transitions in the calcium-sodium ion channels family

We use Brownian dynamics (BD) simulations to study the ionic conduction and valence selectivity of a generic electrostatic model of a biological ion channel as functions of the fixed charge Q(f) at its selectivity filter. We are thus able to reconcile the discrete calcium conduction bands recently revealed in our BD simulations, M0 (Q(f)=1e), M1 (3e), M2 (5e), with a set of sodium conduction bands L0 (0.5e), L1 (1.5e), thereby obtaining a completed pattern of conduction and selectivity bands vs Q(f) for the sodium-calcium channels family. An increase of Q(f) leads to an increase of calcium selectivity: L0 (sodium-selective, nonblocking channel) → M0 (nonselective channel) → L1 (sodium-selective channel with divalent block) → M1 (calcium-selective channel exhibiting the anomalous mole fraction effect). We create a consistent identification scheme where the L0 band is putatively identified with the eukaryotic sodium channel The scheme created is able to account for the experimentally observed mutation-induced transformations between nonselective channels, sodium-selective channels, and calcium-selective channels, which we interpret as transitions between different rows of the identification table. By considering the potential energy changes during permeation, we show explicitly that the multi-ion conduction bands of calcium and sodium channels arise as the result of resonant barrierless conduction. The pattern of periodic conduction bands is explained on the basis of sequential neutralization taking account of self-energy, as Q(f)(z,i)=ze(1/2+i), where i is the order of the band and z is the valence of the ion. Our results confirm the crucial influence of electrostatic interactions on conduction and on the Ca(2+)/Na(+) valence selectivity of calcium and sodium ion channels. The model and results could be also applicable to biomimetic nanopores with charged walls.

Dielectric characterization of costal cartilage chondrocytes

BACKGROUND

Chondrocytes respond to biomechanical and bioelectrochemical stimuli by secreting appropriate extracellular matrix proteins that enable the tissue to withstand the large forces it experiences. Although biomechanical aspects of cartilage are well described, little is known of the bioelectrochemical responses. The focus of this study is to identify bioelectrical characteristics of human costal cartilage cells using dielectric spectroscopy.

METHODS

Dielectric spectroscopy allows non-invasive probing of biological cells. An in house computer program is developed to extract dielectric properties of human costal cartilage cells from raw cell suspension impedance data measured by a microfluidic device. The dielectric properties of chondrocytes are compared with other cell types in order to comparatively assess the electrical nature of chondrocytes.

RESULTS

The results suggest that electrical cell membrane characteristics of chondrocyte cells are close to cardiomyoblast cells, cells known to possess an array of active ion channels. The blocking effect of the non-specific ion channel blocker gadolinium is tested on chondrocytes with a significant reduction in both membrane capacitance and conductance.

CONCLUSIONS

We have utilized a microfluidic chamber to mimic biomechanical events through changes in bioelectrochemistry and described the dielectric properties of chondrocytes to be closer to cells derived from electrically excitably tissues.

GENERAL SIGNIFICANCE

The study describes dielectric characterization of human costal chondrocyte cells using physical tools, where results and methodology can be used to identify potential anomalies in bioelectrochemical responses that may lead to cartilage disorders.

Applied field nonequilibrium molecular dynamics simulations of ion exit from a β-barrel model of the L-type calcium channel

We present results of applied field nonequilibrium molecular dynamics simulations (AF NEMD) of a minimal beta-barrel model channel intended to represent an L-type calcium channel that suggests a possible relationship between glutamate side chain conformational changes and ion flux in calcium channels. The beta-barrel is used to provide a scaffolding for glutamate side chains and a confinement for electrolyte of dimensions similar to the expected channel structure. It was preloaded with ions to explore relative rates of ion exit for different occupancy configurations. Our simulations with an asymmetrical flexible selectivity filter represented by four glutamate side chains (EEEE), one of which differs in initial dihedrals from the other three, indicate a plausible mechanism for the observed anomalous mole fraction effect seen in calcium channels. Apparent rates of electric field-induced exit from channels preloaded with three Na+ ions are much higher than for channels with one Ca2+ followed by two Na+ ions, consistent with the common notion that Ca2+ block of Na+ current is due to competition between the Ca2+ and Na+ ions for the negatively charged (EEEE) locus. In our model, the Ca2+ ion ligates simultaneously to the four negatively charged glutamate side chains and sterically blocks the permeation pathway. Ca2+-relief of Ca2+-block is suggested by a much higher rate of exit for channels preloaded with three Ca2+ ions than for channels with two Ca2+ ions.