1
|
Almog M, Korngreen A. Is realistic neuronal modeling realistic? J Neurophysiol 2016; 116:2180-2209. [PMID: 27535372 DOI: 10.1152/jn.00360.2016] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2016] [Accepted: 08/17/2016] [Indexed: 11/22/2022] Open
Abstract
Scientific models are abstractions that aim to explain natural phenomena. A successful model shows how a complex phenomenon arises from relatively simple principles while preserving major physical or biological rules and predicting novel experiments. A model should not be a facsimile of reality; it is an aid for understanding it. Contrary to this basic premise, with the 21st century has come a surge in computational efforts to model biological processes in great detail. Here we discuss the oxymoronic, realistic modeling of single neurons. This rapidly advancing field is driven by the discovery that some neurons don't merely sum their inputs and fire if the sum exceeds some threshold. Thus researchers have asked what are the computational abilities of single neurons and attempted to give answers using realistic models. We briefly review the state of the art of compartmental modeling highlighting recent progress and intrinsic flaws. We then attempt to address two fundamental questions. Practically, can we realistically model single neurons? Philosophically, should we realistically model single neurons? We use layer 5 neocortical pyramidal neurons as a test case to examine these issues. We subject three publically available models of layer 5 pyramidal neurons to three simple computational challenges. Based on their performance and a partial survey of published models, we conclude that current compartmental models are ad hoc, unrealistic models functioning poorly once they are stretched beyond the specific problems for which they were designed. We then attempt to plot possible paths for generating realistic single neuron models.
Collapse
Affiliation(s)
- Mara Almog
- The Leslie and Susan Gonda Interdisciplinary Brain Research Centre, Bar-Ilan University, Ramat Gan, Israel; and.,The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Alon Korngreen
- The Leslie and Susan Gonda Interdisciplinary Brain Research Centre, Bar-Ilan University, Ramat Gan, Israel; and .,The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| |
Collapse
|
2
|
Moreno JD, Lewis TJ, Clancy CE. Parameterization for In-Silico Modeling of Ion Channel Interactions with Drugs. PLoS One 2016; 11:e0150761. [PMID: 26963710 PMCID: PMC4786197 DOI: 10.1371/journal.pone.0150761] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2015] [Accepted: 02/18/2016] [Indexed: 12/31/2022] Open
Abstract
Since the first Hodgkin and Huxley ion channel model was described in the 1950s, there has been an explosion in mathematical models to describe ion channel function. As experimental data has become richer, models have concomitantly been improved to better represent ion channel kinetic processes, although these improvements have generally resulted in more model complexity and an increase in the number of parameters necessary to populate the models. Models have also been developed to explicitly model drug interactions with ion channels. Recent models of drug-channel interactions account for the discrete kinetics of drug interaction with distinct ion channel state conformations, as it has become clear that such interactions underlie complex emergent kinetics such as use-dependent block. Here, we describe an approach for developing a model for ion channel drug interactions. The method describes the process of extracting rate constants from experimental electrophysiological function data to use as initial conditions for the model parameters. We then describe implementation of a parameter optimization method to refine the model rate constants describing ion channel drug kinetics. The algorithm takes advantage of readily available parallel computing tools to speed up the optimization. Finally, we describe some potential applications of the platform including the potential for gaining fundamental mechanistic insights into ion channel function and applications to in silico drug screening and development.
Collapse
Affiliation(s)
- Jonathan D. Moreno
- Division of Cardiology, Department of Medicine, Barnes-Jewish Hospital, Washington University in St. Louis, St. Louis, MO, United States of America
- * E-mail:
| | - Timothy J. Lewis
- Department of Mathematics, University of California Davis, Davis, CA, United States of America
| | - Colleen E. Clancy
- Department of Pharmacology, University of California Davis, Davis, CA, United States of America
| |
Collapse
|
3
|
Abstract
The ability to record the currents from single ion channels led to the need to extract the underlying kinetic model from such data. This inverse hidden Markov problem is difficult but led to the creation of a software suite called QuB utilizing likelihood optimization. This review presents the software. The software is open source and, in addition to solving kinetic models, has many generic database operations including report generation with publishable graphics, function fitting and scripting for new and repeated processing and AD/DA I/O. The core algorithms allow for constraints such as fixed rates or maintaining detailed balance in the model. All rate constants can be driven by a stimulus and the system can analyze nonstationary data. QuB also can analyze the kinetics of multichannel data where individual events cannot be discriminated, but the fitting algorithms utilize the signal variance as well as the mean to fit models. QuB can be applied to any data appropriately modeled with Markov kinetics and has been utilized to solve ion channels but also the movement of motor proteins, the sleep cycles in mice, and physics processes. [Formula: see text]Special Issue Comment: This is a review about the software QuB that can extract a model from the trajectory. It is connected with the review about treatments when solving single molecules,60 and the reviews about enzymes.61,62
Collapse
Affiliation(s)
- CHRISTOPHER NICOLAI
- Physiology and Biophysics, SUNY Buffalo, 301 Cary Hall, Buffalo, NY 14214, USA
| | - FREDERICK SACHS
- Physiology and Biophysics, SUNY Buffalo, 301 Cary Hall, Buffalo, NY 14214, USA
| |
Collapse
|
4
|
Thomson AS, Rothberg BS. Voltage-dependent inactivation gating at the selectivity filter of the MthK K+ channel. ACTA ACUST UNITED AC 2010; 136:569-79. [PMID: 20937694 PMCID: PMC2964515 DOI: 10.1085/jgp.201010507] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Voltage-dependent K+ channels can undergo a gating process known as C-type inactivation, which involves entry into a nonconducting state through conformational changes near the channel’s selectivity filter. C-type inactivation may involve movements of transmembrane voltage sensor domains, although the mechanisms underlying this form of inactivation may be heterogeneous and are often unclear. Here, we report on a form of voltage-dependent inactivation gating observed in MthK, a prokaryotic K+ channel that lacks a canonical voltage sensor and may thus provide a reduced system to inform on mechanism. In single-channel recordings, we observe that Po decreases with depolarization, with a half-maximal voltage of 96 ± 3 mV. This gating is kinetically distinct from blockade by internal Ca2+ or Ba2+, suggesting that it may arise from an intrinsic inactivation mechanism. Inactivation gating was shifted toward more positive voltages by increasing external [K+] (47 mV per 10-fold increase in [K+]), suggesting that K+ binding at the extracellular side of the channel stabilizes the open-conductive state. The open-conductive state was stabilized by other external cations, and selectivity of the stabilizing site followed the sequence: K+ ≈ Rb+ > Cs+ > Na+ > Li+ ≈ NMG+. Selectivity of the stabilizing site is weaker than that of sites that determine permeability of these ions, suggesting that the site may lie toward the external end of the MthK selectivity filter. We could describe MthK gating over a wide range of positive voltages and external [K+] using kinetic schemes in which the open-conductive state is stabilized by K+ binding to a site that is not deep within the electric field, with the voltage dependence of inactivation arising from both voltage-dependent K+ dissociation and transitions between nonconducting (inactivated) states. These results provide a quantitative working hypothesis for voltage-dependent, K+-sensitive inactivation gating, a property that may be common to other K+ channels.
Collapse
Affiliation(s)
- Andrew S Thomson
- Department of Biochemistry, Temple University School of Medicine, Philadelphia, PA 19140, USA
| | | |
Collapse
|
5
|
Fink M, Noble D. Markov models for ion channels: versatility versus identifiability and speed. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2009; 367:2161-2179. [PMID: 19414451 DOI: 10.1098/rsta.2008.0301] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Markov models (MMs) represent a generalization of Hodgkin-Huxley models. They provide a versatile structure for modelling single channel data, gating currents, state-dependent drug interaction data, exchanger and pump dynamics, etc. This paper uses examples from cardiac electrophysiology to discuss aspects related to parameter estimation. (i) Parameter unidentifiability (found in 9 out of 13 of the considered models) results in an inability to determine the correct layout of a model, contradicting the idea that model structure and parameters provide insights into underlying molecular processes. (ii) The information content of experimental voltage step clamp data is discussed, and a short but sufficient protocol for parameter estimation is presented. (iii) MMs have been associated with high computational cost (owing to their large number of state variables), presenting an obstacle for multicellular whole organ simulations as well as parameter estimation. It is shown that the stiffness of models increases computation time more than the number of states. (iv) Algorithms and software programs are provided for steady-state analysis, analytical solutions for voltage steps and numerical derivation of parameter identifiability. The results provide a new standard for ion channel modelling to further the automation of model development, the validation process and the predictive power of these models.
Collapse
Affiliation(s)
- Martin Fink
- Department of Physiology, Anatomy and Genetics, University of Oxford, Parks Road, Oxford OX1 3PT, UK.
| | | |
Collapse
|
6
|
Abstract
This unit provides detailed descriptions for the steps of patch excision, data acquisition, and data analysis, and elaborates upon the relevant issues discussed in other units from Chapter 3. It includes discussions of the instrumentation for single-channel recording and the key concepts necessary for the interpretation of single-channel data.
Collapse
Affiliation(s)
- M B Jackson
- University of Wisconsin Medical School, Madison, Wisconsin, USA
| |
Collapse
|
7
|
Abstract
We suggest that bubbles are the bistable hydrophobic gates responsible for the on-off transitions of single channel currents. In this view, many types of channels gate by the same physical mechanism-dewetting by capillary evaporation-but different types of channels use different sensors to modulate hydrophobic properties of the channel wall and thereby trigger and control bubbles and gating. Spontaneous emptying of channels has been seen in many simulations. Because of the physics involved, such phase transitions are inherently sensitive, unstable threshold phenomena that are difficult to simulate reproducibly and thus convincingly. We present a thermodynamic analysis of a bubble gate using morphometric density functional theory of classical (not quantum) mechanics. Thermodynamic analysis of phase transitions is generally more reproducible and less sensitive to details than simulations. Anesthetic actions of inert gases-and their interactions with hydrostatic pressure (e.g., nitrogen narcosis)-can be easily understood by actions on bubbles. A general theory of gas anesthesia may involve bubbles in channels. Only experiments can show whether, or when, or which channels actually use bubbles as hydrophobic gates: direct observation of bubbles in channels is needed. Existing experiments show thin gas layers on hydrophobic surfaces in water and suggest that bubbles nearly exist in bulk water.
Collapse
|
8
|
Gurkiewicz M, Korngreen A. A numerical approach to ion channel modelling using whole-cell voltage-clamp recordings and a genetic algorithm. PLoS Comput Biol 2007; 3:e169. [PMID: 17784781 PMCID: PMC1963494 DOI: 10.1371/journal.pcbi.0030169] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2006] [Accepted: 07/16/2007] [Indexed: 11/23/2022] Open
Abstract
The activity of trans-membrane proteins such as ion channels is the essence of neuronal transmission. The currently most accurate method for determining ion channel kinetic mechanisms is single-channel recording and analysis. Yet, the limitations and complexities in interpreting single-channel recordings discourage many physiologists from using them. Here we show that a genetic search algorithm in combination with a gradient descent algorithm can be used to fit whole-cell voltage-clamp data to kinetic models with a high degree of accuracy. Previously, ion channel stimulation traces were analyzed one at a time, the results of these analyses being combined to produce a picture of channel kinetics. Here the entire set of traces from all stimulation protocols are analysed simultaneously. The algorithm was initially tested on simulated current traces produced by several Hodgkin-Huxley–like and Markov chain models of voltage-gated potassium and sodium channels. Currents were also produced by simulating levels of noise expected from actual patch recordings. Finally, the algorithm was used for finding the kinetic parameters of several voltage-gated sodium and potassium channels models by matching its results to data recorded from layer 5 pyramidal neurons of the rat cortex in the nucleated outside-out patch configuration. The minimization scheme gives electrophysiologists a tool for reproducing and simulating voltage-gated ion channel kinetics at the cellular level. Voltage-gated ion channels affect neuronal integration of information. Some neurons express more than ten different types of voltage-gated ion channels, making information processing a highly convoluted process. Kinetic modelling of ion channels is an important method for unravelling the role of each channel type in neuronal function. However, the most commonly used analysis techniques suffer from shortcomings that limit the ability of researchers to rapidly produce physiologically relevant models of voltage-gated ion channels and of neuronal physiology. We show that conjugating a stochastic search algorithm with ionic currents measured using multiple voltage-clamp protocols enables the semi-automatic production of models of voltage-gated ion channels. Once fully automated, this approach may be used for high throughput analysis of voltage-gated currents. This in turn will greatly shorten the time required for building models of neuronal physiology to facilitate our understanding of neuronal behaviour.
Collapse
Affiliation(s)
- Meron Gurkiewicz
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, Israel
- The Leslie and Susan Gonda Multidisciplinary Brain Research Center, Bar-Ilan University, Ramat-Gan, Israel
| | - Alon Korngreen
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, Israel
- The Leslie and Susan Gonda Multidisciplinary Brain Research Center, Bar-Ilan University, Ramat-Gan, Israel
- * To whom correspondence should be addressed. E-mail:
| |
Collapse
|
9
|
Gurkiewicz M, Korngreen A. Recording, analysis, and function of dendritic voltage-gated channels. Pflugers Arch 2006; 453:283-92. [PMID: 16604366 DOI: 10.1007/s00424-006-0076-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2006] [Revised: 03/07/2006] [Accepted: 03/13/2006] [Indexed: 11/25/2022]
Abstract
Ever since the publication of the Hamill et al. [Hamill et al., Pflügers Arch, 391:85-100, 1981] paper and the following increase in popularity of acute brain slice preparations, there has been a large increase in the volume of publications investigating voltage-gated channels in the central nervous system using the patch-clamp technique. In the preceding decade, investigations of voltage-gated channels have moved out of the somatic region into dendrites providing much needed information about dendritic voltage-gated channels. In this study, we review some aspects related to the investigation of voltage-gated ion channels in dendrites: recording, analysis, and function.
Collapse
Affiliation(s)
- Meron Gurkiewicz
- The Mina and Everand Goodman Faculty of Life Sciences and the Leslie and Susan Gonda Multidisciplinary Brain Research Center, Bar-Ilan University, Ramat-Gan, 52900, Israel
| | | |
Collapse
|
10
|
Rothberg BS, Magleby KL. Voltage and Ca2+ activation of single large-conductance Ca2+-activated K+ channels described by a two-tiered allosteric gating mechanism. J Gen Physiol 2000; 116:75-99. [PMID: 10871641 PMCID: PMC2229615 DOI: 10.1085/jgp.116.1.75] [Citation(s) in RCA: 129] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
The voltage- and Ca2+-dependent gating mechanism of large-conductance Ca2+-activated K+ (BK) channels from cultured rat skeletal muscle was studied using single-channel analysis. Channel open probability (Po) increased with depolarization, as determined by limiting slope measurements (11 mV per e-fold change in Po; effective gating charge, q(eff), of 2.3 +/- 0.6 e(o)). Estimates of q(eff) were little changed for intracellular Ca2+ (Ca2+(i)) ranging from 0.0003 to 1,024 microM. Increasing Ca2+(i) from 0.03 to 1,024 microM shifted the voltage for half maximal activation (V(1/2)) 175 mV in the hyperpolarizing direction. V(1/2) was independent of Ca2+(i) for Ca2+(i) < or = 0.03 microM, indicating that the channel can be activated in the absence of Ca2+(i). Open and closed dwell-time distributions for data obtained at different Ca2+(i) and voltage, but at the same Po, were different, indicating that the major action of voltage is not through concentrating Ca2+ at the binding sites. The voltage dependence of Po arose from a decrease in the mean closing rate with depolarization (q(eff) = -0.5 e(o)) and an increase in the mean opening rate (q(eff) = 1.8 e(o)), consistent with voltage-dependent steps in both the activation and deactivation pathways. A 50-state two-tiered model with separate voltage- and Ca2+-dependent steps was consistent with the major features of the voltage and Ca2+ dependence of the single-channel kinetics over wide ranges of Ca2+(i) (approximately 0 through 1,024 microM), voltage (+80 to -80 mV), and Po (10(-4) to 0.96). In the model, the voltage dependence of the gating arises mainly from voltage-dependent transitions between closed (C-C) and open (O-O) states, with less voltage dependence for transitions between open and closed states (C-O), and with no voltage dependence for Ca2+-binding and unbinding. The two-tiered model can serve as a working hypothesis for the Ca2+- and voltage-dependent gating of the BK channel.
Collapse
Affiliation(s)
- Brad S. Rothberg
- Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430
| | - Karl L. Magleby
- Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430
| |
Collapse
|
11
|
Nimigean CM, Magleby KL. Functional coupling of the beta(1) subunit to the large conductance Ca(2+)-activated K(+) channel in the absence of Ca(2+). Increased Ca(2+) sensitivity from a Ca(2+)-independent mechanism. J Gen Physiol 2000; 115:719-36. [PMID: 10828246 PMCID: PMC2232893 DOI: 10.1085/jgp.115.6.719] [Citation(s) in RCA: 76] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
Coexpression of the beta(1) subunit with the alpha subunit (mSlo) of BK channels increases the apparent Ca(2+) sensitivity of the channel. This study investigates whether the mechanism underlying the increased Ca(2+) sensitivity requires Ca(2+), by comparing the gating in 0 Ca(2+)(i) of BK channels composed of alpha subunits to those composed of alpha+beta(1) subunits. The beta(1) subunit increased burst duration approximately 20-fold and the duration of gaps between bursts approximately 3-fold, giving an approximately 10-fold increase in open probability (P(o)) in 0 Ca(2+)(i). The effect of the beta(1) subunit on increasing burst duration was little changed over a wide range of P(o) achieved by varying either Ca(2+)(i) or depolarization. The effect of the beta(1) subunit on increasing the durations of the gaps between bursts in 0 Ca(2+)(i) was preserved over a range of voltage, but was switched off as Ca(2+)(i) was increased into the activation range. The Ca(2+)-independent, beta(1) subunit-induced increase in burst duration accounted for 80% of the leftward shift in the P(o) vs. Ca(2+)(i) curve that reflects the increased Ca(2+) sensitivity induced by the beta(1) subunit. The Ca(2+)-dependent effect of the beta(1) subunit on the gaps between bursts accounted for the remaining 20% of the leftward shift. Our observation that the major effects of the beta(1) subunit are independent of Ca(2+)(i) suggests that the beta(1) subunit mainly alters the energy barriers of Ca(2+)-independent transitions. The changes in gating induced by the beta(1) subunit differ from those induced by depolarization, as increasing P(o) by depolarization or by the beta(1) subunit gave different gating kinetics. The complex gating kinetics for both alpha and alpha+beta(1) channels in 0 Ca(2+)(i) arise from transitions among two to three open and three to five closed states and are inconsistent with Monod-Wyman-Changeux type models, which predict gating among only one open and one closed state in 0 Ca(2+)(i).
Collapse
Affiliation(s)
- C M Nimigean
- Department of Physiology and Biophysics, University of Miami School of Medicine, FL 33101-6430, USA.
| | | |
Collapse
|
12
|
Rothberg BS, Magleby KL. Gating kinetics of single large-conductance Ca2+-activated K+ channels in high Ca2+ suggest a two-tiered allosteric gating mechanism. J Gen Physiol 1999; 114:93-124. [PMID: 10398695 PMCID: PMC2229641 DOI: 10.1085/jgp.114.1.93] [Citation(s) in RCA: 114] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
The Ca2+-dependent gating mechanism of large-conductance calcium-activated K+ (BK) channels from cultured rat skeletal muscle was examined from low (4 microM) to high (1,024 microM) intracellular concentrations of calcium (Ca2+i) using single-channel recording. Open probability (Po) increased with increasing Ca2+i (K0. 5 11.2 +/- 0.3 microM at +30 mV, Hill coefficient of 3.5 +/- 0.3), reaching a maximum of approximately 0.97 for Ca2+i approximately 100 microM. Increasing Ca2+i further to 1,024 microM had little additional effect on either Po or the single-channel kinetics. The channels gated among at least three to four open and four to five closed states at high levels of Ca2+i (>100 microM), compared with three to four open and five to seven closed states at lower Ca2+i. The ability of kinetic schemes to account for the single-channel kinetics was examined with simultaneous maximum likelihood fitting of two-dimensional (2-D) dwell-time distributions obtained from low to high Ca2+i. Kinetic schemes drawn from the 10-state Monod-Wyman-Changeux model could not describe the dwell-time distributions from low to high Ca2+i. Kinetic schemes drawn from Eigen's general model for a ligand-activated tetrameric protein could approximate the dwell-time distributions but not the dependency (correlations) between adjacent intervals at high Ca2+i. However, models drawn from a general 50 state two-tiered scheme, in which there were 25 closed states on the upper tier and 25 open states on the lower tier, could approximate both the dwell-time distributions and the dependency from low to high Ca2+i. In the two-tiered model, the BK channel can open directly from each closed state, and a minimum of five open and five closed states are available for gating at any given Ca2+i. A model that assumed that the apparent Ca2+-binding steps can reach a maximum rate at high Ca2+i could also approximate the gating from low to high Ca2+i. The considered models can serve as working hypotheses for the gating of BK channels.
Collapse
Affiliation(s)
- Brad S. Rothberg
- From the Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430
| | - Karl L. Magleby
- From the Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430
| |
Collapse
|
13
|
Nimigean CM, Magleby KL. The beta subunit increases the Ca2+ sensitivity of large conductance Ca2+-activated potassium channels by retaining the gating in the bursting states. J Gen Physiol 1999; 113:425-40. [PMID: 10051518 PMCID: PMC2222905 DOI: 10.1085/jgp.113.3.425] [Citation(s) in RCA: 114] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Coexpression of the beta subunit (KV,Cabeta) with the alpha subunit of mammalian large conductance Ca2+- activated K+ (BK) channels greatly increases the apparent Ca2+ sensitivity of the channel. Using single-channel analysis to investigate the mechanism for this increase, we found that the beta subunit increased open probability (Po) by increasing burst duration 20-100-fold, while having little effect on the durations of the gaps (closed intervals) between bursts or on the numbers of detected open and closed states entered during gating. The effect of the beta subunit was not equivalent to raising intracellular Ca2+ in the absence of the beta subunit, suggesting that the beta subunit does not act by increasing all the Ca2+ binding rates proportionally. The beta subunit also inhibited transitions to subconductance levels. It is the retention of the BK channel in the bursting states by the beta subunit that increases the apparent Ca2+ sensitivity of the channel. In the presence of the beta subunit, each burst of openings is greatly amplified in duration through increases in both the numbers of openings per burst and in the mean open times. Native BK channels from cultured rat skeletal muscle were found to have bursting kinetics similar to channels expressed from alpha subunits alone.
Collapse
Affiliation(s)
- C M Nimigean
- Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430, USA
| | | |
Collapse
|
14
|
Rothberg BS, Magleby KL. Kinetic structure of large-conductance Ca2+-activated K+ channels suggests that the gating includes transitions through intermediate or secondary states. A mechanism for flickers. J Gen Physiol 1998; 111:751-80. [PMID: 9607935 PMCID: PMC2217154 DOI: 10.1085/jgp.111.6.751] [Citation(s) in RCA: 59] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/1998] [Accepted: 04/13/1998] [Indexed: 11/20/2022] Open
Abstract
Mechanisms for the Ca2+-dependent gating of single large-conductance Ca2+-activated K+ channels from cultured rat skeletal muscle were developed using two-dimensional analysis of single-channel currents recorded with the patch clamp technique. To extract and display the essential kinetic information, the kinetic structure, from the single channel currents, adjacent open and closed intervals were binned as pairs and plotted as two-dimensional dwell-time distributions, and the excesses and deficits of the interval pairs over that expected for independent pairing were plotted as dependency plots. The basic features of the kinetic structure were generally the same among single large-conductance Ca2+-activated K+ channels, but channel-specific differences were readily apparent, suggesting heterogeneities in the gating. Simple gating schemes drawn from the Monod- Wyman-Changeux (MWC) model for allosteric proteins could approximate the basic features of the Ca2+ dependence of the kinetic structure. However, consistent differences between the observed and predicted dependency plots suggested that additional brief lifetime closed states not included in MWC-type models were involved in the gating. Adding these additional brief closed states to the MWC-type models, either beyond the activation pathway (secondary closed states) or within the activation pathway (intermediate closed states), improved the description of the Ca2+ dependence of the kinetic structure. Secondary closed states are consistent with the closing of secondary gates or channel block. Intermediate closed states are consistent with mechanisms in which the channel activates by passing through a series of intermediate conformations between the more stable open and closed states. It is the added secondary or intermediate closed states that give rise to the majority of the brief closings (flickers) in the gating.
Collapse
Affiliation(s)
- B S Rothberg
- Department of Physiology and Biophysics, University of Miami School of Medicine, Miami, Florida 33101-6430, USA
| | | |
Collapse
|
15
|
Pastushenko VP, Schindler H. Level detection in ion channel records via idealization by statistical filtering and likelihood optimization. Philos Trans R Soc Lond B Biol Sci 1997; 352:39-51. [PMID: 9051714 PMCID: PMC1691908 DOI: 10.1098/rstb.1997.0004] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
A parameter-free method is presented for the level detection in ion channel records via recovery of step wise current changes. No assumptions about ion channel mechanism are made. The primary detection of the transitions is made by statistical filtering the data using the Student's t-test. The event currents are calculated as the average value of the current between two adjacent transitions. An optimal ideal trace is found by maximization of a likelihood function. The distribution of event currents recovered from the raw data is then analysed, again by using the Student's t-test, for their grouping into separate statistical ensembles, defining current levels. The method is subjected to rigorous test using simulated data, and is compared with several other methods. It produces the levels of channel current, their noise amplitudes and distributions of dwell times, the desired information for constructing the channel mechanism.
Collapse
Affiliation(s)
- V P Pastushenko
- Institute for Biophysics, Johannes Kepler University of Linz, Linz-Auhof, Austria
| | | |
Collapse
|
16
|
Currás MC, Pallotta BS. Single-channel evidence for glycine and NMDA requirement in NMDA receptor activation. Brain Res 1996; 740:27-40. [PMID: 8973795 DOI: 10.1016/s0006-8993(96)00845-1] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
N-Methyl-D-aspartate (NMDA) receptor dose-response relationships that are based on macroscopic currents suggest that NMDA and a different agonist molecule, glycine, must together activate the channel. Since single-channel recordings have a much higher resolution than whole-cell currents, they provide a highly sensitive test for the absolute requirement of NMDA channel opening for glycine. Rapid application of 10-300 microM NMDA to outside-out patches from cultured cortical neurons evoked substantial single-channel activity in the absence of added glycine. However, in the presence of a high affinity and highly selective glycine-site antagonist, 5,7-dichlorokynurenate (DCK), NMDA failed to evoke any openings on its own. Channel openings could not be produced by saturating concentrations of NMDA (up to 1 mM) but were evoked when glycine was added to the test solution. Glycine alone (up to 100 microM) was similarly ineffective in the continuous presence of D(-)-2-amino-5-phosphonovaleric acid (D-APV), an NMDA-site antagonist. Reversal of antagonist blockade by the appropriate ligand (glycine or NMDA) and the normal appearance and duration of channel openings evoked in the presence of either antagonist ruled out open channel block. These single-channel data confirm the hypothesis that both NMDA and glycine are coagonists of the NMDA receptor. Furthermore, the coagonist requirement increases the potential targets for therapeutic drugs aimed at blocking the pathologies resulting from overactivation of NMDA receptors.
Collapse
Affiliation(s)
- M C Currás
- Department of Pharmacology, University of North Carolina School of Medicine, Chapel Hill 27599-7365, USA.
| | | |
Collapse
|
17
|
Manivannan K, Mathias RT, Gudowska-Nowak E. Description of interacting channel gating using a stochastic Markovian model. Bull Math Biol 1996; 58:141-74. [PMID: 8819758 DOI: 10.1007/bf02458286] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Single-channel recordings from membrane patches frequently exhibit multiple conductance levels. In some preparations, the steady-state probabilities of observing these levels do not follow a binomial distribution. This behavior has been reported in sodium channels, potassium channels, acetylcholine receptor channels and gap junction channels. A non-binomial distribution suggests interaction of the channels or the presence of channels or the presence of channels with different open probabilities. However, the current trace sometimes exhibits single transitions spanning several levels. Since the probability of simultaneous transitions of independent channels is infinitesimally small, such observations strongly suggest a cooperative gating behavior. We present a Markov model to describe the cooperative gating of channels using only the all-points current amplitude histograms for the probability of observing the various conductance levels. We investigate the steady-state (or equilibrium) properties of a system of N channels and provide a scheme to express all the probabilities in terms of just two parameters. The main feature of our model is that lateral interaction of channels gives rise to cooperative gating. Another useful feature is the introduction of the language of graph theory which can potentially provide a different avenue to study ion channel kinetics. We write down explicit expressions for systems of two, three and four channels and provide a procedure to describe the system of N channels.
Collapse
Affiliation(s)
- K Manivannan
- Department of Physiology and Biophysics, State University of New York, Stony Brook 11794, USA
| | | | | |
Collapse
|
18
|
Townsend C, Rosenberg RL. Characterization of a chloride channel reconstituted from cardiac sarcoplasmic reticulum. J Membr Biol 1995; 147:121-36. [PMID: 8568849 DOI: 10.1007/bf00233541] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We have characterized a voltage-sensitive chloride channel from cardiac sarcoplasmic reticulum (SR) following reconstitution of porcine heart SR into planar lipid bilayers. In 250 mM KCl, the channel had a main conductance level of 130 pS and exhibited two substrates of 61 and 154 pS. The channel was very selective for Cl- over K+ or Na+ (PK+/PCl- = 0.012 and PNa+/PCl- approximately 0.040). It was permeable to several anions and displayed the following sequence of anion permeability: SCN- > I- > NO3- approximately Br- > Cl- > F- > HCOO-. Single-channel conductance saturated with increasing Cl- concentrations (Km = 900 mM and gamma max = 488 pS). Channel activity was voltage dependent, with an open probability ranging from approximately 1.0 around 0 mV to approximately 0.5 at +80 mV. From -20 to +80 mV, channel gating was time-independent. However, at voltages below -40 mV the channel entered a long-lasting closed state. Mean open times varied with voltage, from approximately 340 msec at -20 mV to approximately 6 msec at +80 mV, whereas closed times were unaffected. The channel was not Ca(2+)-dependent. Channel activity was blocked by disulfonic stilbenes, arylaminobenzoates, zinc, and cadmium. Single-channel conductance was sensitive to trans pH, ranging from approximately 190 pS at pH 5.5 to approximately 60 pS at pH 9.0. These characteristics are different from those previously described for Cl- channels from skeletal or cardiac muscle SR.
Collapse
Affiliation(s)
- C Townsend
- Department of Pharmacology, University of North Carolina, Chapel Hill 27599-7365, USA
| | | |
Collapse
|
19
|
Kleckner NW, Pallotta BS. Burst kinetics of single NMDA receptor currents in cell-attached patches from rat brain cortical neurons in culture. J Physiol 1995; 486 ( Pt 2):411-26. [PMID: 7473207 PMCID: PMC1156531 DOI: 10.1113/jphysiol.1995.sp020822] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
1. The patch-clamp technique was used to record single-channel currents from cell-attached patches on rat brain cortical neurons in culture. The composition of the open and shut intervals during bursts of openings was studied in N-methyl-D-aspartate (NMDA) receptors exposed to 1 microM NMDA and 10 microM glycine at a membrane potential of -70 mV. 2. Open interval histograms were constructed for openings at each position (first, second, third, etc.) during all bursts. Distributions from openings two to five were fitted with similar (two or three) exponential components. The first opening from all bursts tended to be of shorter duration than the other openings. 3. Bursts were sorted according to the number of openings they contained and duration histograms were obtained for the openings at each position (one to five) during bursts of up to five openings. For bursts containing two or more openings, the distribution of open durations obtained at a given position were similar to each other regardless of the number of openings in the burst. 4. In bursts of two or more openings (up to five), duration histograms from the openings at each position in the burst were fitted with two or three exponential components that were similar for each opening. Bursts consisting of a single opening had a different distribution, having a relatively larger component of short duration. 5. Shut intervals during bursts were described by two exponential components with average time constants (and relative areas) (means +/- S.E.M.) of 0.06 +/- 0.01 ms (0.59 +/- 0.02) and 0.64 +/- 0.02 ms (0.41 +/- 0.02). Their distribution was independent of the numbers of openings in the bursts, their position within the burst, and the types of openings (long or short duration) contained within the burst. 6. These results suggest that each opening in bursts of two or more openings is kinetically similar to every other opening regardless of burst length. Analogously, each shut period during a burst was similar to every other. A kinetic model with three open and four closed intraburst states is shown to be consistent with these results for bursts of two or more openings.
Collapse
Affiliation(s)
- N W Kleckner
- Department of Pharmacology, School of Medicine, University of North Carolina, Chapel Hill 27599-7365, USA
| | | |
Collapse
|
20
|
Donnelly JL, Pallotta BS. Single-channel currents from diethylpyrocarbonate-modified NMDA receptors in cultured rat brain cortical neurons. J Gen Physiol 1995; 105:837-59. [PMID: 7561746 PMCID: PMC2216961 DOI: 10.1085/jgp.105.6.837] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
The role of histidine residues in the function of N-methyl-D-aspartate (NMDA)-activated channels was tested with the histidine-modifying reagent diethylpyrocarbonate (DEP) applied to cells and membrane patches from rat brain cortical neurons in culture. Channels in excised outside-out patches that were treated with 3 mM DEP for 15-30 s (pH 6.5) showed an average 3.4-fold potentiation in steady state open probability when exposed to NMDA and glycine. Analysis of the underlying alterations in channel gating revealed no changes in the numbers of kinetic states: distributions of open intervals were fitted with three exponential components, and four components described the shut intervals, in both control and DEP-modified channels. However, the distribution of shut intervals was obviously different after DEP treatment, consistent with the single-channel current record. After modification, the proportion of long shut states was decreased while the time constants were largely unaffected. Burst kinetics reflected these effects with an increase in the average number of openings/burst from 1.5 (control) to 2.2 (DEP), and a decrease in the average interburst interval from 54.1 to 38.2 ms. These effects were most likely due to histidine modification because other reagents (n-acetylimidazole and 2,4,6-trinitrobenzene 1-sulfonic acid) that are specific for residues other than histidine failed to reproduce the effects of DEP, whereas hydroxylamine could restore channel open probability to control levels. In contrast to these effects on channel gating, DEP had no effect on average single-channel conductance or reversal potential under bi-ionic (Na+:Cs+) conditions. Inhibition by zinc was also unaffected by DEP. We propose a channel gating model in which transitions between single- and multi-opening burst modes give rise to the channel activity observed under steady state conditions. When adjusted to account for the effects of DEP, this model suggests that one or more extracellular histidine residues involved in channel gating are associated with a single kinetic state.
Collapse
Affiliation(s)
- J L Donnelly
- Department of Pharmacology, School of Medicine, University of North Carolina at Chapel Hill 27599-7365, USA
| | | |
Collapse
|
21
|
Nonaka T, Matsuzaki K, Kawahara K, Suzuki K, Hoshino M. Monovalent cation selective channel in the apical membrane of rat inner medullary collecting duct cells in primary culture. BIOCHIMICA ET BIOPHYSICA ACTA 1995; 1233:163-74. [PMID: 7532435 DOI: 10.1016/0005-2736(94)00241-g] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Ion channels in the apical membrane of rat inner medullary collecting duct (IMCD) were investigated by the patch clamp technique. Owing to the histological heterogeneity of IMCD, cells were cultured from the lower half of the inner medulla of Wistar rat kidney. Channel activity was rarely seen in cell attached patch, but membrane excision activated multiple units of 28.2 +/- 0.7 pS cation selective channel. A Na or K selective channel was not found. The 28 pS channel showed membrane voltage dependency, no rectification, almost equal permeability to monovalent cations (Na/K/Li/Cs/Rb/NH4 = 1:1.00:0.82:0.97:1.10:1.71) and no significant permeation to anions or divalent cations. Calcium of the cytoplasmic side from 10(-7) M to 10(-4) M affected the mean number of open channels (nPo) dose-dependently in excised patch (IC50 = 5 x 10(-6) M). 1 mM of ATP, ADP, AMP and gadolinium reversibly suppressed nPo to near zero whereas amiloride, cAMP or cGMP had no effect. Multiple conductance substates were frequently observed. These results suggested that this channel belongs to the nonselective cation channels which has been identified in other epithelia and is not responsible for amiloride sensitive Na transport through IMCD cells.
Collapse
Affiliation(s)
- T Nonaka
- Tenjinbashi Clinic, Tokyo, Japan
| | | | | | | | | |
Collapse
|
22
|
Milone M, Hutchinson DO, Engel AG. Patch-clamp analysis of the properties of acetylcholine receptor channels at the normal human endplate. Muscle Nerve 1994; 17:1364-9. [PMID: 7526207 DOI: 10.1002/mus.880171204] [Citation(s) in RCA: 27] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Normative data were obtained on the kinetic properties of the acetylcholine receptor (AChR) channel at the human motor endplate by patch-clamp analysis. Single channel currents were recorded from 34 endplates of 8 nonweak subjects in the presence of 1 micron acetylcholine (ACh) at 22 +/- 0.5 degrees C. The vast majority of channels opened to a conductance of about 60 pS. The dwell-time distributions of these channels were well described as the sum of two exponential functions. The mean duration of the dominant longer component was 1.9 ms for the open intervals and 3.04 ms for the bursts. At three endplates, a small proportion of the channels had lower conductance and longer open time, resembling immature AChR channels. At 28 endplates it was also possible to obtain an estimate of the rate constant for channel closure (alpha) and approximate estimates for the rate constants of channel opening (beta) and ACh dissociation (k-2). Estimates of k-2 varied by 15% with methods of estimation. This is attributed to errors inherent in estimating the duration of the briefest channel events. The normative data will be useful for evaluating pathologic alterations in the kinetic properties of the AChR channel found in some congenital myasthenic syndromes.
Collapse
Affiliation(s)
- M Milone
- Department of Neurology, Mayo Clinic, Rochester, Minnesota 55905
| | | | | |
Collapse
|