Optimal estimation of ion-channel kinetics from macroscopic currents.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang W
PROVIDER: S-EPMC3335051 | biostudies-literature | 2012
REPOSITORIES: biostudies-literature
Publications
Optimal estimation of ion-channel kinetics from macroscopic currents.
Wang Wei W Xiao Feng F Zeng Xuhui X Yao Jing J Yuchi Ming M Ding Jiuping J
PloS one 20120420 4
Markov modeling provides an effective approach for modeling ion channel kinetics. There are several search algorithms for global fitting of macroscopic or single-channel currents across different experimental conditions. Here we present a particle swarm optimization(PSO)-based approach which, when used in combination with golden section search (GSS), can fit macroscopic voltage responses with a high degree of accuracy (errors within 1%) and reasonable amount of calculation time (less than 10 hou ...[more]
Similar Datasets
Project description:A new method is described that accurately estimates kinetic constants, conductance and number of ion channels from macroscopic currents. The method uses both the time course and the strength of correlations between different time points of macroscopic currents and utilizes the property of semiseparability of covariance matrix for computationally efficient estimation of current likelihood and its gradient. The number of calculation steps scales linearly with the number of channel states as opposed to the cubic dependence in a previously described method. Together with the likelihood gradient evaluation, which is almost independent of the number of model parameters, the new approach allows evaluation of kinetic models with very complex topologies. We demonstrate applicability of the method to analysis of synaptic currents by estimating accurately rate constants of a 7-state model used to simulate GABAergic macroscopic currents.
| S-EPMC3248447 | biostudies-literature
Project description:Electrochemical signaling in the brain depends on pentameric ligand-gated ion channels (pLGICs). Recently, crystal structures of prokaryotic pLGIC homologues from Erwinia chrysanthemi (ELIC) and Gloeobacter violaceus (GLIC) in presumed closed and open channel states have been solved, which provide insight into the structural mechanisms underlying channel activation. Although structural studies involving both ELIC and GLIC have become numerous, thorough functional characterizations of these channels are still needed to establish a reliable foundation for comparing kinetic properties. Here, we examined the kinetics of ELIC and GLIC current activation, desensitization, and deactivation and compared them to the GABAA receptor, a prototypic eukaryotic pLGIC. Outside-out patch-clamp recordings were performed with HEK-293T cells expressing ELIC, GLIC, or α1β2γ2L GABAA receptors, and ultra-fast ligand application was used. In response to saturating agonist concentrations, we found both ELIC and GLIC current activation were two to three orders of magnitude slower than GABAA receptor current activation. The prokaryotic channels also had slower current desensitization on a timescale of seconds. ELIC and GLIC current deactivation following 25 s pulses of agonist (cysteamine and pH 4.0 buffer, respectively) were relatively fast with time constants of 24.9 ± 5.1 ms and 1.2 ± 0.2 ms, respectively. Surprisingly, ELIC currents evoked by GABA activated very slowly with a time constant of 1.3 ± 0.3 s and deactivated even slower with a time constant of 4.6 ± 1.2 s. We conclude that the prokaryotic pLGICs undergo similar agonist-mediated gating transitions to open and desensitized states as eukaryotic pLGICs, supporting their use as experimental models. Their uncharacteristic slow activation, slow desensitization and rapid deactivation time courses are likely due to differences in specific structural elements, whose future identification may help uncover mechanisms underlying pLGIC gating transitions.
| S-EPMC3833957 | biostudies-literature
Project description:Transcranial focused ultrasound (US) has been demonstrated to stimulate neurons in animals and humans, but the mechanism of this effect is unknown. It has been hypothesized that US, a mechanical stimulus, may mediate cellular discharge by activating mechanosensitive ion channels embedded within cellular membranes. To test this hypothesis, we expressed potassium and sodium mechanosensitive ion channels (channels of the two-pore-domain potassium family (K2P) including TREK-1, TREK-2, TRAAK; NaV1.5) in the Xenopus oocyte system. Focused US (10 MHz, 0.3-4.9 W/cm(2)) modulated the currents flowing through the ion channels on average by up to 23%, depending on channel and stimulus intensity. The effects were reversible upon repeated stimulation and were abolished when a channel blocker (ranolazine to block NaV1.5, BaCl2 to block K2P channels) was applied to the solution. These data reveal at the single cell level that focused US modulates the activity of specific ion channels to mediate transmembrane currents. These findings open doors to investigations of the effects of US on ion channels expressed in neurons, retinal cells, or cardiac cells, which may lead to important medical applications. The findings may also pave the way to the development of sonogenetics: a non-invasive, US-based analogue of optogenetics.
| S-EPMC4845013 | biostudies-other
Project description:Ion currents associated with channel proteins in the presence of membrane potential are ubiquitous in cellular and organelle membranes. When an ion current occurs through a channel protein, Joule heating should occur. However, this Joule heating seems to have been largely overlooked in biology. Here we show theoretical investigation of Joule heating involving channel proteins in biological processes. We used electrochemical potential to derive the Joule's law for an ion current through an ion transport protein in the presence of membrane potential, and we suggest that heat production and absorption can occur. Simulation of temperature distribution around a single channel protein with the Joule heating revealed that the temperature increase was as small as <10-3 K, although an ensemble of channel proteins was suggested to exhibit a noticeable temperature increase. Thereby, we theoretically investigated the Joule heating of systems containing ensembles of channel proteins. Nerve is known to undergo rapid heat production followed by heat absorption during the action potential, and our simulation of Joule heating for a squid giant axon combined with the Hodgkin-Huxley model successfully reproduced the feature of the heat. Furthermore, we extended the theory of Joule heating to uncoupling protein 1 (UCP1), a solute carrier family transporter, which is important to the non-shivering thermogenesis in brown adipose tissue mitochondria (BATM). Our calculations showed that the Joule heat involving UCP1 was comparable to the literature calorimetry data of BATM. Joule heating of ion transport proteins is likely to be one of important mechanisms of cellular thermogenesis.
| S-EPMC10728626 | biostudies-literature
Project description:A soft repulsion (SR) model of short-range interactions between mobile ions and protein atoms is introduced in the framework of continuum representation of the protein and solvent. The Poisson-Nernst-Plank (PNP) theory of ion transport through biological channels is modified to incorporate this soft wall protein model. Two sets of SR parameters are introduced. The first is parametrized for all essential amino acid residues using all atom molecular dynamic simulations; the second is a truncated Lennard-Jones potential. We have further designed an energy-based algorithm for the determination of the ion accessible volume, which is appropriate for a particular system discretization. The effects of these models of short-range interactions were tested by computing current-voltage characteristics of the ?-hemolysin channel. The introduced SR potentials significantly improve prediction of channel selectivity. In addition, we studied the effect of the choice of some space-dependent diffusion coefficient distributions on the predicted current-voltage properties. We conclude that the diffusion coefficient distributions largely affect total currents and have little effect on rectifications, selectivity, or reversal potential. The PNP-SR algorithm is implemented in a new efficient parallel Poisson, Poisson-Boltzmann, and PNP equation solver, also incorporated in a graphical molecular modeling package HARLEM.
| S-EPMC3059120 | biostudies-literature
Project description:Acid-sensing ion channel-2a (ASIC2a) is an amiloride-blockable proton-gated cation channel, probably contributing to sour-taste detection in rat taste cells. To isolate another subtype of the sour-taste receptor, we screened a rat circumvallate papilla cDNA library and identified ASIC2b, an N-terminal splice variant of ASIC2a. Reverse transcription-PCR analyses confirmed the expression of ASIC2b transcripts in the circumvallate papilla and, moreover, demonstrated its expression in the foliate and fungiform papillae. Immunohistochemical analyses revealed that ASIC2b, as well as ASIC2a, was expressed in a subpopulation of taste cells in the circumvallate, foliate, and fungiform papillae, and some of the cells displayed both ASIC2a and ASIC2b immunoreactivities. Subsequent coimmunoprecipitation studies with circumvallate papillae extracts indicated that ASIC2b associated with ASIC2a to form assemblies and, together with our immunohistochemical findings, strongly suggested that both ASIC2 subunits formed heteromeric channels in taste cells in the circumvallate, foliate, and fungiform papillae. Oocyte electrophysiology demonstrated that the ASIC2a/ASIC2b channel generated maximal inward currents at a pH of < or =2.0, which is in agreement with the in vivo pH sensitivity of rat taste cells, and that the amiloride sensitivity of the heteromer decreased with decreasing pH and was almost completely abolished at a pH of 2.0. These findings provide persuasive explanations for the amiloride insensitivity of acid-induced responses of rat taste cells.
| S-EPMC6742202 | biostudies-literature
Project description:Ion channels play a central role in setting gradients of ion concentration and electrostatic potentials, which in turn regulate sensory systems and other functions. Based on the structure of the open configuration of the Kv1.2 channel and the suggestion that the two ends of the N-terminal inactivating peptide form a bivalent complex that simultaneously blocks the channel pore and binds to the cytoplasmic T1 domain, we propose a six state kinetic model that for the first time reproduces the kinetics of recovery of the Drosophila Shaker over the full range of time scales and hyperpolarization potentials, including tail currents. The model is motivated by a normal mode analysis of the inactivated channel that suggests that a displacement consistent with models of the closed state propagates to the T1 domain via the S1-T1 linker. This motion stretches the bound (inactivating) peptide, hastening the unblocking of the pore. This pulling force is incorporated into the rates of the open to blocked states, capturing the fast recovery phase of the current for repolarization events shorter than 1 ms. If the membrane potential is hyperpolarized, essential dynamics further suggests that the T1 domain returns to a configuration where the peptide is un-stretched and the S1-T1 linker is extended. Coupling this novel hyperpolarized substate to the closed, open and blocked pore states is enough to quantitatively estimate the number of open channels as a function of time and membrane potential. A straightforward prediction of the model is that a slow ramping of the potential leads to very small currents.
| S-EPMC2551740 | biostudies-literature
Project description:Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour.
| S-EPMC5599066 | biostudies-literature
Project description:KEY POINTS:Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long. Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions. The new model predicts the current under traditional square-wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well. The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell-cell variability in current kinetics for the first time. ABSTRACT:Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics - the voltage-dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell-specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell-cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.
| S-EPMC5978315 | biostudies-literature
Project description:The prokaryotic K(+) channel KcsA is activated by intracellular protons and its gating is modulated by transmembrane voltage. Typically, KcsA functions have been studied under steady-state conditions, using macroscopic Rb(+)-flux experiments and single-channel current measurements. These studies have provided limited insights into the gating kinetics of KcsA due to its low open probability, uncertainties in the number of channels in the patch, and a very strong intrinsic kinetic variability. In this work, we have carried out a detailed analysis of KcsA gating under nonstationary conditions by examining the influence of pH and voltage on the activation, deactivation, and slow-inactivation gating events. We find that activation and deactivation gating of KcsA are predominantly modulated by pH without a significant effect of voltage. Activation gating showed sigmoidal pH dependence with a pKa of approximately 4.2 and a Hill coefficient of approximately 2. In the sustained presence of proton, KcsA undergoes a time-dependent decay of conductance. This inactivation process is pH independent but is modulated by voltage and the nature of permeant ion. Recovery from inactivation occurs via deactivation and also appears to be voltage dependent. We further find that inactivation in KcsA is not entirely a property of the open-conducting channel but can also occur from partially "activated" closed states. The time course of onset and recovery of the inactivation process from these pre-open closed states appears to be different from the open-state inactivation, suggesting the presence of multiple inactivated states with diverse kinetic pathways. This information has been analyzed together with a detailed study of KcsA single-channel behavior (in the accompanying paper) in the framework of a kinetic model. Taken together our data constitutes the first quantitative description of KcsA gating.
| S-EPMC2151670 | biostudies-literature