Stochastic models and simulation of ion channel dynamics

The behaviour of ion channels within cardiac and neuronal cells is intrinsically stochastic in nature. When the number of channels is small this stochastic noise is large and can have an impact on the dynamics of the system which is potentially an issue when modelling small neurons and drug block in cardiac cells. While exact methods correctly capture the stochastic dynamics of a system they are computationally expensive, restricting their inclusion into tissue level models and so approximations to exact methods are often used instead. The other issue in modelling ion channel dynamics is that the transition rates are voltage dependent, adding a level of complexity as the channel dynamics are coupled to the membrane potential. By assuming that such transition rates are constant over each time step, it is possible to derive a stochastic differential equation (SDE), in the same manner as for biochemical reaction networks, that describes the stochastic dynamics of ion channels. While such a model is more computationally efficient than exact methods we show that there are analytical problems with the resulting SDE as well as issues in using current numerical schemes to solve such an equation. We therefore make two contributions: develop a different model to describe the stochastic ion channel dynamics that analytically behaves in the correct manner and also discuss numerical methods that preserve the analytical properties of the model.

Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Helix-sheet interconversion in a synthetic pore lining peptide derived from a designed ion channel

We have designed a four-helix protein that is expected to tetramerize in the membrane to form an ion channel with a structurally well defined pore. A synthetic peptide corresponding to the channel lining helix facilitates ion transport across liposomal membranes and largely helical in membranes. Detailed circular dichroism studies of the peptide in methanol, water and methanal-water mixtures reveal that it is helical in methanol, beta-structured in 97.5% water and a combination of these two structures at intermediate compositions of methanol and water. A fluorescence resonance energy transfer study of the peptide shows that the peptide is monomeric in methanol but undergoes extensive anti-parallel aggregation in aqueous solution.

Ion channel gene variants predisposing to severe human ventricular arrhythmias

Congenital long QT syndrome (LQTS) with an estimated prevalence of 1:2000-1:10 000 manifests with prolonged QT interval on electrocardiogram and risk for ventricular arrhythmias and sudden death. Several ion channel genes and hundreds of mutations in these genes have been identified to underlie the disorder. In Finland, four LQTS founder mutations of potassium channel genes account for up to 40-70% of genetic spectrum of LQTS. Acquired LQTS has similar clinical manifestations, but often arises from usage of QT-prolonging medication or electrolyte disturbances. A prolonged QT interval is associated with increased morbidity and mortality not only in clinical LQTS but also in patients with ischemic heart disease and in the general population. The principal aim of this study was to estimate the actual prevalence of LQTS founder mutations in Finland and to calculate their effect on QT interval in the Finnish background population. Using a large population-based sample of over 6000 Finnish individuals from the Health 2000 Survey, we identified LQTS founder mutations KCNQ1 G589D (n=8), KCNQ1 IVS7-2A>G (n=1), KCNH2 L552S (n=2), and KCNH2 R176W (n=16) in 27 study participants. This resulted in a weighted prevalence estimate of 0.4% for LQTS in Finland. Using a linear regression model, the founder mutations resulted in a 22- to 50-ms prolongation of the age-, sex-, and heart rate-adjusted QT interval. Collectively, these data suggest that one of 250 individuals in Finland may be genetically predisposed to ventricular arrhythmias arising from the four LQTS founder mutations. A KCNE1 D85N minor allele with a frequency of 1.4% was associated with a 10-ms prolongation in adjusted QT interval and could thus identify individuals at increased risk of ventricular arrhythmias at the population level. In addition, the previously reported associations of KCNH2 K897T, KCNH2 rs3807375, and NOS1AP rs2880058 with QT interval duration were confirmed in the present study. In a separate study, LQTS founder mutations were identified in a subgroup of acquired LQTS, providing further evidence that congenital LQTS gene mutations may underlie acquired LQTS. Catecholaminergic polymorphic ventricular tachycardia (CPVT) is characterized by exercise-induced ventricular arrhythmias in a structurally normal heart and results from defects in the cardiac Ca2+ signaling proteins, mainly ryanodine receptor type 2 (RyR2). In a patient population of typical CPVT, RyR2 mutations were identifiable in 25% (4/16) of patients, implying that noncoding variants or other genes are involved in CPVT pathogenesis. A 1.1 kb RyR2 exon 3 deletion was identified in two patients independently, suggesting that this region may provide a new target for RyR2-related molecular genetic studies. Two novel RyR2 mutations showing a gain-of-function defect in vitro were identified in three victims of sudden cardiac death. Extended pedigree analyses revealed some surviving mutation carriers with mild structural abnormalities of the heart and resting ventricular arrhythmias suggesting that not all RyR2 mutations lead to a typical CPVT phenotype, underscoring the relevance of tailored risk stratification of a RyR2 mutation carrier.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Calcium-permeable acid-sensing ion channel is a molecular target of the neurotoxic metal ion lead

Acid-sensing ion channels (ASICs) are emerging as fundamental players in the regulation of neural plasticity and in pathological conditions. Here we showed that lead (Pb2+), a well known neurotoxic metal ion, reversibly and concentration-dependently inhib

Coupling between NMDA receptor and acid-sensing ion channel contributes to ischemic neuronal death

Acid-sensing ion channels (ASICs) composed of ASIC1a subunit exhibit a high Ca2+ permeability and play important roles in synaptic plasticity and acid-induced cell death. Here, we show that ischemia enhances ASIC currents through the phosphorylation at Ser478 and Ser479 of ASIC1a, leading to exacerbated ischemic cell death. The phosphorylation is catalyzed by Ca2+/calmodulin-dependent protein kinase II (CaMKII) activity, as a result of activation of NR2B-containing N-methyl-D-aspartate subtype of glutamate receptors (NMDARs) during ischemia. Furthermore, NR2B-specific antagonist, CaMKII inhibitor, or overexpression of mutated form of ASIC1a with Ser478 or Ser479 replaced by alanine (ASICla-S478A, ASIC1a-S479A) in cultured hippocampal neurons prevented ischemia-induced enhancement of ASIC currents, cytoplasmic Ca2+ elevation, as well as neuronal death. Thus, NMDAR-CaMKII cascade is functionally coupled to ASICs and contributes to acidotoxicity during ischemia. Specific blockade of NMDAR/CaMKII-ASIC coupling may reduce neuronal death after ischemia and other pathological conditions involving excessive glutamate release and acidosis.

Upregulation of acid-sensing ion channel ASIC1a in spinal dorsal horn neurons contributes to inflammatory pain hypersensitivity

Development of chronic pain involves alterations in peripheral nociceptors as well as elevated neuronal activity in multiple regions of the CNS. Previous pharmacological and behavioral studies suggest that peripheral acid-sensing ion channels (ASICs) cont

Cell types, network homeostasis, and pathological compensation from a biologically plausible ion channel expression model.

How do neurons develop, control, and maintain their electrical signaling properties in spite of ongoing protein turnover and perturbations to activity? From generic assumptions about the molecular biology underlying channel expression, we derive a simple model and show how it encodes an "activity set point" in single neurons. The model generates diverse self-regulating cell types and relates correlations in conductance expression observed in vivo to underlying channel expression rates. Synaptic as well as intrinsic conductances can be regulated to make a self-assembling central pattern generator network; thus, network-level homeostasis can emerge from cell-autonomous regulation rules. Finally, we demonstrate that the outcome of homeostatic regulation depends on the complement of ion channels expressed in cells: in some cases, loss of specific ion channels can be compensated; in others, the homeostatic mechanism itself causes pathological loss of function.

Study of the ion channel behavior of didodecyldimethylammonium bromide formed bilayer lipid membrane stimulated by PF6-

Bilayer lipid membranes ( BLM) formed from didode-cyldimethylammonium bromide were made on the freshly exposed surface of a glassy carbon (GC) and were demonstrated by the ac impedance spectroscopy. The ion channels of membrane properties induced by PF6- were studied by the cyclic voltammetric methods. Experimental results indicated that the ion channel of BLM was open in the presence of the PF6- due to the interaction of PF6- with the BLM, while it was switched off in the absence of PF6-. Because the ion channel behavior was affected by the concentration of PF6-, a sensor for PF6- can be developed.