Naturally occurring mutations at the acetylcholine receptor binding site independently alter ACh binding and channel gating

By defining functional defects in a congenital myasthenic syndrome (CMS), we show that two mutant residues, located in a binding site region of the acetylcholine receptor (AChR) epsilon subunit, exert opposite effects on ACh binding and suppress channel gating. Single channel kinetic analysis reveals that the first mutation, epsilon N182Y, increases ACh affinity for receptors in the resting closed state, which promotes sequential occupancy of the binding sites and discloses rate constants for ACh occupancy of the nonmutant alphadelta site. Studies of the analogous mutation in the delta subunit, deltaN187Y, disclose rate constants for ACh occupancy of the nonmutant alpha epsilon site. The second CMS mutation, epsilon D175N, reduces ACh affinity for receptors in the resting closed state; occupancy of the mutant site still promotes gating because a large difference in affinity is maintained between closed and open states. epsilon D175N impairs overall gating, however, through an effect independent of ACh occupancy. When mapped on a structural model of the AChR binding site, epsilon N182Y localizes to the interface with the alpha subunit, and epsilon D175 to the entrance of the ACh binding cavity. Both epsilon N182Y and epsilon D175 show state specificity in affecting closed relative to desensitized state affinities, suggesting that the protein chain harboring epsilon N182 and epsilon D175 rearranges in the course of receptor desensitization. The overall results show that key residues at the ACh binding site differentially stabilize the agonist bound to closed, open and desensitized states, and provide a set point for gating of the channel.     

Clinicopathological‐genetic features of congenital myasthenic syndrome from a Chinese neuromuscular centre

Congenital myasthenic syndrome (CMS) encompasses a heterogeneous group of inherited disorders affecting nerve transmission across the neuromuscular junction. The aim of this study was to characterize the clinical, physiological, pathohistological and genetic features of nine unrelated Chinese patients with CMS from a single neuromuscular centre. A total of nine patients aged from neonates to 34 years were enrolled who exhibited initial symptoms. Physical examinations revealed that all patients exhibited muscle weakness. Muscle biopsies demonstrated multiple myopathological changes, including increased fibre size variation, myofibrillar network disarray, necrosis, myofiber grouping, regeneration, fibre atrophy and angular fibres. Genetic testing revealed six different mutated genes, including AGRN (2/9), CHRNE (1/9), GFPT1 (1/9), GMPPB (1/9), PLEC (3/9) and SCN4A (1/9). In addition, patients exhibited differential responses to pharmacological treatment. Prompt utilization of genetic testing will identify novel variants and expand our understanding of the phenotype of this rare syndrome. Our findings contribute to the clinical, pathohistological and genetic spectrum of congenital myasthenic syndrome in China.     

Calcium regulation of single ryanodine receptor channel gating analyzed using HMM/MCMC statistical methods

Type-II ryanodine receptor channels (RYRs) play a fundamental role in intracellular Ca(2+) dynamics in heart. The processes of activation, inactivation, and regulation of these channels have been the subject of intensive research and the focus of recent debates. Typically, approaches to understand these processes involve statistical analysis of single RYRs, involving signal restoration, model estimation, and selection. These tasks are usually performed by following rather phenomenological criteria that turn models into self-fulfilling prophecies. Here, a thorough statistical treatment is applied by modeling single RYRs using aggregated hidden Markov models. Inferences are made using Bayesian statistics and stochastic search methods known as Markov chain Monte Carlo. These methods allow extension of the temporal resolution of the analysis far beyond the limits of previous approaches and provide a direct measure of the uncertainties associated with every estimation step, together with a direct assessment of why and where a particular model fails. Analyses of single RYRs at several Ca(2+) concentrations are made by considering 16 models, some of them previously reported in the literature. Results clearly show that single RYRs have Ca(2+)-dependent gating modes. Moreover, our results demonstrate that single RYRs responding to a sudden change in Ca(2+) display adaptation kinetics. Interestingly, best ranked models predict microscopic reversibility when monovalent cations are used as the main permeating species. Finally, the extended bandwidth revealed the existence of novel fast buzz-mode at low Ca(2+) concentrations.     

The role of loop 5 in acetylcholine receptor channel gating

Nicotinic acetylcholine receptor channel (AChR) gating is an organized sequence of molecular motions that couples a change in the affinity for ligands at the two transmitter binding sites with a change in the ionic conductance of the pore. Loop 5 (L5) is a nine-residue segment (mouse alpha-subunit 92-100) that links the beta4 and beta5 strands of the extracellular domain and that (in the alpha-subunit) contains binding segment A. Based on the structure of the acetylcholine binding protein, we speculate that in AChRs L5 projects from the transmitter binding site toward the membrane along a subunit interface. We used single-channel kinetics to quantify the effects of mutations to alphaD97 and other L5 residues with respect to agonist binding (to both open and closed AChRs), channel gating (for both unliganded and fully-liganded AChRs), and desensitization. Most alphaD97 mutations increase gating (up to 168-fold) but have little or no effect on ligand binding or desensitization. Rate-equilibrium free energy relationship analysis indicates that alphaD97 moves early in the gating reaction, in synchrony with the movement of the transmitter binding site (Phi = 0.93, which implies an open-like character at the transition state). alphaD97 mutations in the two alpha-subunits have unequal energetic consequences for gating, but their contributions are independent. We conclude that the key, underlying functional consequence of alphaD97 perturbations is to increase the unliganded gating equilibrium constant. L5 emerges as an important and early link in the AChR gating reaction which, in the absence of agonist, serves to increase the relative stability of the closed conformation of the protein.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations: From tension sensing to channel opening

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

CFTR gating I: Characterization of the ATP-dependent gating of a phosphorylation-independent CFTR channel (DeltaR-CFTR)

The CFTR chloride channel is activated by phosphorylation of serine residues in the regulatory (R) domain and then gated by ATP binding and hydrolysis at the nucleotide binding domains (NBDs). Studies of the ATP-dependent gating process in excised inside-out patches are very often hampered by channel rundown partly caused by membrane-associated phosphatases. Since the severed DeltaR-CFTR, whose R domain is completely removed, can bypass the phosphorylation-dependent regulation, this mutant channel might be a useful tool to explore the gating mechanisms of CFTR. To this end, we investigated the regulation and gating of the DeltaR-CFTR expressed in Chinese hamster ovary cells. In the cell-attached mode, basal DeltaR-CFTR currents were always obtained in the absence of cAMP agonists. Application of cAMP agonists or PMA, a PKC activator, failed to affect the activity, indicating that the activity of DeltaR-CFTR channels is indeed phosphorylation independent. Consistent with this conclusion, in excised inside-out patches, application of the catalytic subunit of PKA did not affect ATP-induced currents. Similarities of ATP-dependent gating between wild type and DeltaR-CFTR make this phosphorylation-independent mutant a useful system to explore more extensively the gating mechanisms of CFTR. Using the DeltaR-CFTR construct, we studied the inhibitory effect of ADP on CFTR gating. The Ki for ADP increases as the [ATP] is increased, suggesting a competitive mechanism of inhibition. Single channel kinetic analysis reveals a new closed state in the presence of ADP, consistent with a kinetic mechanism by which ADP binds at the same site as ATP for channel opening. Moreover, we found that the open time of the channel is shortened by as much as 54% in the presence of ADP. This unexpected result suggests another ADP binding site that modulates channel closing.     

Modulation of proton-induced current fluctuations in the human nicotinic acetylcholine receptor channel

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel that switches upon activation from a closed state to a full conducting state. We found that the mutation δ S268K, located at 12′ position of the second transmembrane domain of the δ subunit of the human nAChR generates a long-lived intermediate conducting state, from which openings to a wild-type like conductance level occur on a submillisecond time scale. Aiming to understand the interplay between structural changes near the 12′ position and channel gating, we investigated the influence of various parameters: different ligands (acetylcholine, choline and epibatidine), ligand concentrations, transmembrane voltages and both fetal and adult nAChRs. Since sojourns in the high conductance state are not fully resolved in time, spectral noise analysis was used as a complement to dwell time analysis to determine the gating rate constants. Open channel current fluctuations are described by a two-state Markov model. The characteristic time of the process is markedly influenced by the ligand and the receptor type, whereas the frequency of openings to the high conductance state increases with membrane hyperpolarization. Conductance changes are discussed with regard to reversible transfer reaction of single protons at the lysine 12′ side chain.     

G-protein-gated TRP-like cationic channel activated by muscarinic receptors: effect of potential on single-channel gating

There is little information about the mechanisms by which G-protein-coupled receptors gate ion channels although many ionotropic receptors are well studied. We have investigated gating of the muscarinic cationic channel, which mediates the excitatory effect of acetylcholine in smooth muscles, and proposed a scheme consisting of four pairs of closed and open states. Channel kinetics appeared to be the same in cell-attached or outside-out patches whether the channel was activated by carbachol application or by intracellular dialysis with GTPgammaS. Since in the latter case G-proteins are permanently active, it is concluded that the cationic channel is the major determinant of its own gating, similarly to the K(ACh) channel (Ivanova-Nikolova, T.T., and G.E. Breitwieser. 1997. J. Gen. Physiol. 109:245-253). Analysis of adjacent-state dwell times revealed connections between the states that showed features conserved among many other ligand-gated ion channels (e.g., nAChR, BK(Ca) channel). Open probability (P(O)) of the cationic channel was increased by membrane depolarization consistent with the prominent U-shaped I-V relationship of the muscarinic whole-cell current at negative potentials. Membrane potential affected transitions within each closed-open state pair but had little effect on transitions between pairs; thus, the latter are likely to be caused by interactions of the channel with its ligands, e.g., Ca(2+) and Galphao-GTP. Channel activity was highly heterogeneous, as was evident from the prominent cycling behavior when P(O) was measured over 5-s intervals. This was related to the variable frequency of openings (as in the K(ACh) channel) and, especially, to the number of long openings between consecutive long shuttings. Analysis of the underlying Markov chain in terms of probabilities allowed us to evaluate the contribution of each open state to the integral current (from shortest to longest open state: 0.1, 3, 24, and 73%) as P(O) increased 525-fold in three stages.     

Independence and cooperativity in rearrangements of a potassium channel voltage sensor revealed by single subunit fluorescence

Voltage-gated potassium channels are composed of four subunits. Voltage-dependent activation of these channels consists of a depolarization-triggered series of charge-carrying steps that occur in each subunit. These major charge-carrying steps are followed by cooperative step(s) that lead to channel opening. Unlike the late cooperative steps, the major charge-carrying steps have been proposed to occur independently in each of the channel subunits. In this paper, we examine this further. We showed earlier that the two major charge-carrying steps are associated with two sequential outward transmembrane movements of the charged S4 segment. We now use voltage clamp fluorometry to monitor these S4 movements in individual subunits of heterotetrameric channels. In this way, we estimate the influence of one subunit's S4 movement on another's when the energetics of their transmembrane movements differ. Our results show that the first S4 movement occurs independently in each subunit, while the second occurs cooperatively. At least part of the cooperativity appears to be intrinsic to the second S4 charge-carrying rearrangement. Such cooperativity in gating of voltage-dependent channels has great physiological relevance since it can affect both action potential threshold and rate of propagation.     

Toward elucidating the heat activation mechanism of the TRPV1 channel gating by molecular dynamics simulation

As a key cellular sensor, the TRPV1 cation channel undergoes a gating transition from a closed state to an open state in response to various physical and chemical stimuli including noxious heat. Despite years of study, the heat activation mechanism of TRPV1 gating remains enigmatic at the molecular level. Toward elucidating the structural and energetic basis of TRPV1 gating, we have performed molecular dynamics (MD) simulations (with cumulative simulation time of 3 μs), starting from the high‐resolution closed and open structures of TRPV1 solved by cryo‐electron microscopy. In the closed‐state simulations at 30°C, we observed a stably closed channel constricted at the lower gate (near residue I679), while the upper gate (near residues G643 and M644) is dynamic and undergoes flickery opening/closing. In the open‐state simulations at 60°C, we found higher conformational variation consistent with a large entropy increase required for the heat activation, and both the lower and upper gates are dynamic with transient opening/closing. Through ensemble‐based structural analyses of the closed state versus the open state, we revealed pronounced closed‐to‐open conformational changes involving the membrane proximal domain (MPD) linker, the outer pore, and the TRP helix, which are accompanied by breaking/forming of a network of closed/open‐state specific hydrogen bonds. By comparing the closed‐state simulations at 30°C and 60°C, we observed heat‐activated conformational changes in the MPD linker, the outer pore, and the TRP helix that resemble the closed‐to‐open conformational changes, along with partial formation of the open‐state specific hydrogen bonds. Some of the residues involved in the above key hydrogen bonds were validated by previous mutational studies. Taken together, our MD simulations have offered rich structural and dynamic details beyond the static structures of TRPV1, and promising targets for future mutagenesis and functional studies of the TRPV1 channel. Proteins 2016; 84:1938–1949. © 2016 Wiley Periodicals, Inc.     

Regulation of K+ flow by a ring of negative charges in the outer pore of BKCa channels. Part II: Neutralization of aspartate 292 reduces long channel openings and gating current slow component

Neutralization of the aspartate near the selectivity filter in the GYGD pore sequence (D292N) of the voltage- and Ca(2+)-activated K+ channel (MaxiK, BKCa) does not prevent conduction like the corresponding mutation in Shaker channel, but profoundly affects major biophysical properties of the channel (Haug, T., D. Sigg, S. Ciani, L. Toro, E. Stefani, and R. Olcese. 2004. J. Gen. Physiol. 124:173-184). Upon depolarizations, the D292N mutant elicited mostly gating current, followed by small or no ionic current, at voltages where the wild-type hSlo channel displayed robust ionic current. In fact, while the voltage dependence of the gating current was not significantly affected by the mutation, the overall activation curve was shifted by approximately 20 mV toward more depolarized potentials. Several lines of evidence suggest that the mutation prevents population of certain open states that in the wild type lead to high open probability. The activation curves of WT and D292N can both be fitted to the sum of two Boltzmann distributions with identical slope factors and half activation potentials, just by changing their relative amplitudes. The steeper and more negative component of the activation curve was drastically reduced by the D292N mutation (from 0.65 to 0.30), suggesting that the population of open states that occurs early in the activation pathway is reduced. Furthermore, the slow component of the gating current, which has been suggested to reflect transitions from closed to open states, was greatly reduced in D292N channels. The D292N mutation also affected the limiting open probability: at 0 mV, the limiting open probability dropped from approximately 0.5 for the wild-type channel to 0.06 in D292N (in 1 mM [Ca2+]i). In addition to these effects on gating charge and open probability, as already described in Part I, the D292N mutation introduces a approximately 40% reduction of outward single channel conductance, as well as a strong outward rectification.     

Voltage and Ca2+ activation of single large-conductance Ca2+-activated K+ channels described by a two-tiered allosteric gating mechanism

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.     

Allosteric regulation of BK channel gating by Ca(2+) and Mg(2+) through a nonselective, low affinity divalent cation site.

The ability of membrane voltage to activate high conductance, calcium-activated (BK-type) K(+) channels is enhanced by cytosolic calcium (Ca(2+)). Activation is sensitive to a range of [Ca(2+)] that spans over four orders of magnitude. Here, we examine the activation of BK channels resulting from expression of cloned mouse Slo1 alpha subunits at [Ca(2+)] and [Mg(2+)] up to 100 mM. The half-activation voltage (V(0.5)) is steeply dependent on [Ca(2+)] in the micromolar range, but shows a tendency towards saturation over the range of 60-300 microM Ca(2+). As [Ca(2+)] is increased to millimolar levels, the V(0.5) is strongly shifted again to more negative potentials. When channels are activated by 300 microM Ca(2+), further addition of either mM Ca(2+) or mM Mg(2+) produces similar negative shifts in steady-state activation. Millimolar Mg(2+) also produces shifts of similar magnitude in the complete absence of Ca(2+). The ability of millimolar concentrations of divalent cations to shift activation is primarily correlated with a slowing of BK current deactivation. At voltages where millimolar elevations in [Ca(2+)] increase activation rates, addition of 10 mM Mg(2+) to 0 Ca(2+) produces little effect on activation time course, while markedly slowing deactivation. This suggests that Mg(2+) does not participate in Ca(2+)-dependent steps that influence current activation rate. We conclude that millimolar Mg(2+) and Ca(2+) concentrations interact with low affinity, relatively nonselective divalent cation binding sites that are distinct from higher affinity, Ca(2+)-selective binding sites that increase current activation rates. A symmetrical model with four independent higher affinity Ca(2+) binding steps, four voltage sensors, and four independent lower affinity Ca(2+)/Mg(2+) binding steps describes well the behavior of G-V curves over a range of Ca(2+) and Mg(2+). The ability of a broad range of [Ca(2+)] to produce shifts in activation of Slo1 conductance can, therefore, be accounted for by multiple types of divalent cation binding sites.     

Activation of the endothelin receptor inhibits the G protein-coupled inwardly rectifying potassium channel by a phospholipase A2-mediated mechanism

To develop a malleable system to model the well-described, physiological interactions between Gq/11 - coupled receptor and Gi/o-coupled receptor signaling, we coexpressed the endothelin A receptor, the mu-opioid receptor, and the G protein-coupled inwardly rectifying potassium channel (Kir 3) heteromultimers in Xenopus laevis oocytes. Activation of the Gi/o-coupled mu-opioid receptor strongly increased Kir 3 channel current, whereas activation of the Gq/11-coupled endothelin A receptor inhibited the Kir 3 response evoked by mu-opioid receptor activation. The magnitude of the inhibition of Kir 3 was channel subtype specific; heteromultimers composed of Kir 3.1 and Kir 3.2 or Kir 3.1 and Kir 3.4 were significantly more sensitive to the effects of endothelin-1 than heteromultimers composed of Kir 3.1 and Kir 3.5. The difference in sensitivity of the heteromultimers suggests that the endothelin-induced inhibition of the opioid- activated current was caused by an effect at the channel rather than at the opioid receptor. The endothelin-1-mediated inhibition was mimicked by arachidonic acid and blocked by the phospholipase A2 inhibitor arachidonoyl trifluoromethyl ketone. Consistent with a possible phospholipase A2-mediated mechanism, the endothelin-1 effect was blocked by calcium chelation with BAPTA-AM and was not affected by kinase inhibition by either staurosporine or genistein. The data suggest the hypothesis that Gq/11-coupled receptor activation may interfere with Gi/o-coupled receptor signaling by the activation of phospholipase A2 and subsequent inhibition of effector function by a direct effect of an eicosanoid on the channel.