Intra- and intersubunit cooperativity in activation of BK channels by Ca2+

The activation of BK channels by Ca(2+) is highly cooperative, with small changes in intracellular Ca(2+) concentration having large effects on open probability (Po). Here we examine the mechanism of cooperative activation of BK channels by Ca(2+). Each of the four subunits of BK channels has a large intracellular COOH terminus with two different high-affinity Ca(2+) sensors: an RCK1 sensor (D362/D367) located on the RCK1 (regulator of conductance of K(+)) domain and a Ca-bowl sensor located on or after the RCK2 domain. To determine interactions among these Ca(2+) sensors, we examine channels with eight different configurations of functional high-affinity Ca(2+) sensors on the four subunits. We find that the RCK1 sensor and Ca bowl contribute about equally to Ca(2+) activation of the channel when there is only one high-affinity Ca(2+) sensor per subunit. We also find that an RCK1 sensor and a Ca bowl on the same subunit are much more effective in increasing Po than when they are on different subunits, indicating positive intrasubunit cooperativity. If it is assumed that BK channels have a gating ring similar to MthK channels with alternating RCK1 and RCK2 domains and that the Ca(2+) sensors act at the flexible (rather than fixed) interfaces between RCK domains, then a comparison of the distribution of Ca(2+) sensors with the observed responses suggest that the interface between RCK1 and RCK2 domains on the same subunit is flexible. On this basis, intrasubunit cooperativity arises because two high-affinity Ca(2+) sensors acting across a flexible interface are more effective in opening the channel than when acting at separate interfaces. An allosteric model incorporating intrasubunit cooperativity nested within intersubunit cooperativity could approximate the Po vs. Ca(2+) response for eight possible subunit configurations of the high-affinity Ca(2+) sensors as well as for three additional configurations from a previous study.     

Unraveling Subunit Cooperativity in Homotetrameric HCN2 Channels

In a multimeric receptor protein, the binding of a ligand can modulate the binding of a succeeding ligand. This phenomenon, called cooperativity, is caused by the interaction of the receptor subunits. By using a complex Markovian model and a set of parameters determined previously, we analyzed how the successive binding of four ligands leads to a complex cooperative interaction of the subunits in homotetrameric HCN2 pacemaker channels. The individual steps in the model were characterized by Gibbs free energies for the equilibria and activation energies, specifying the affinity of the binding sites and the transition rates, respectively. Moreover, cooperative free energies were calculated for each binding step in both the closed and the open channel. We show that the cooperativity sequence positive-negative-positive determined for the binding affinity is generated by the combined effect of very different cooperativity sequences determined for the binding and unbinding rates, which are negative-negative-positive and no-negative-no, respectively. It is concluded that in the ligand-induced activation of HCN2 channels, the sequence of cooperativity based on the binding affinity is caused by two even qualitatively different sequences of cooperativity that are based on the rates of ligand binding and unbinding.     

Gating and ionic currents reveal how the BKCa channel's Ca2+ sensitivity is enhanced by its beta1 subunit

Large-conductance Ca(2+)-activated K(+) channels (BK(Ca) channels) are regulated by the tissue-specific expression of auxiliary beta subunits. Beta1 is predominantly expressed in smooth muscle, where it greatly enhances the BK(Ca) channel's Ca(2+) sensitivity, an effect that is required for proper regulation of smooth muscle tone. Here, using gating current recordings, macroscopic ionic current recordings, and unitary ionic current recordings at very low open probabilities, we have investigated the mechanism that underlies this effect. Our results may be summarized as follows. The beta1 subunit has little or no effect on the equilibrium constant of the conformational change by which the BK(Ca) channel opens, and it does not affect the gating charge on the channel's voltage sensors, but it does stabilize voltage sensor activation, both when the channel is open and when it is closed, such that voltage sensor activation occurs at more negative voltages with beta1 present. Furthermore, beta1 stabilizes the active voltage sensor more when the channel is closed than when it is open, and this reduces the factor D by which voltage sensor activation promotes opening by approximately 24% (16.8-->12.8). The effects of beta1 on voltage sensing enhance the BK(Ca) channel's Ca(2+) sensitivity by decreasing at most voltages the work that Ca(2+) binding must do to open the channel. In addition, however, in order to fully account for the increase in efficacy and apparent Ca(2+) affinity brought about by beta1 at negative voltages, our studies suggest that beta1 also decreases the true Ca(2+) affinity of the closed channel, increasing its Ca(2+) dissociation constant from approximately 3.7 microM to between 4.7 and 7.1 microM, depending on how many binding sites are affected.     

Rapid activation of the cardiac ryanodine receptor by submillisecond calcium stimuli

The local control concept of excitation-contraction coupling in the heart postulates that the activity of the sarcoplasmic reticulum ryanodine receptor channels (RyR) is controlled by Ca(2+) entry through adjoining sarcolemmal single dihydropyridine receptor channels (DHPRs). One unverified premise of this hypothesis is that the RyR must be fast enough to track the brief (<0.5 ms) Ca(2+) elevations accompanying single DHPR channel openings. To define the kinetic limits of effective trigger Ca(2+) signals, we recorded activity of single cardiac RyRs in lipid bilayers during rapid and transient increases in Ca(2+) generated by flash photolysis of DM-nitrophen. Application of such Ca(2+) spikes (amplitude approximately 10-30 microM, duration approximately 0.1-0.4 ms) resulted in activation of the RyRs with a probability that increased steeply (apparent Hill slope approximately 2.5) with spike amplitude. The time constants of RyR activation were 0.07-0.27 ms, decreasing with spike amplitude. To fit the rising portion of the open probability, a single exponential function had to be raised to a power n approximately 3. We show that these data could be adequately described with a gating scheme incorporating four sequential Ca(2+)-sensitive closed states between the resting and the first open states. These results provide evidence that brief Ca(2+) triggers are adequate to activate the RyR, and support the possibility that RyR channels are governed by single DHPR openings. They also provide evidence for the assumption that RyR activation requires binding of multiple Ca(2+) ions in accordance with the tetrameric organization of the channel protein.     

Variable Subunit Contact and Cooperativity of Hemoglobins

Tertiary structures of proteins are conserved better than their primary structures during evolution. Quaternary structures or subunit organizations, however, are not always conserved. A typical case is found in hemoglobin family. Although human, Scapharca, and Urechis have tetrameric hemoglobins, their subunit contacts are completely different from each other. We report here that only one or two amino acid replacements are enough to create a new contact between subunits. Such a small number of chance replacements is expected during the evolution of hemoglobins. This result explains why different modes of subunit interaction evolved in animal hemoglobins. In contrast, certain interactions between subunits are necessary for cooperative oxygen binding. Cooperative oxygen binding is observed often in dimeric and tetrameric hemoglobins. Conformational change of a subunit induced by the first oxygen binding to the heme group is transmitted through the subunit contacts and increases the affinity of the second oxygen. The tetrameric hemoglobins from humans and Scapharca have cooperativity in spite of their different modes of subunit contact, but the one from Urechis does not. The relationship between cooperativity and the mode of subunit contacts is not clear. We compared the atomic interactions at the subunit contact surface of cooperative and non-cooperative tetrameric hemoglobins. We show that heme-contact modules M3–M6 play a key role in the subunit contacts responsible for cooperativity. A module was defined as a contiguous peptide segment having compact conformation and its average length is about 15 amino acid residues. We show that the cooperative hemoglobins have interactins involving at least two pairs of modules among the four heme-contact modules at subunit contact. Received: 12 January 2001 / Accepted: 3 April 2001     

Calmodulin oxidation and methionine to glutamine substitutions reveal methionine residues critical for functional interaction with ryanodine receptor-1

Calmodulin (CaM) binds to the skeletal muscle ryanodine receptor Ca(2+) release channel (RyR1) with high affinity, and it may act as a Ca(2+)-sensing subunit of the channel. Apo-CaM increases RyR1 channel activity, but Ca(2+)-CaM is inhibitory. Here we examine the functional effects of CaM oxidation on RyR1 regulation by both apo-CaM and Ca(2+)-CaM, as assessed via determinations of [(3)H]ryanodine and [(35)S]CaM binding to skeletal muscle sarcoplasmic reticulum vesicles. Oxidation of all nine CaM Met residues abolished functional interactions of CaM with RyR1. Incomplete CaM oxidation, affecting 5-8 Met residues, increased the CaM concentration required to modulate RyR1, having a greater effect on the apo-CaM species. Mutating individual CaM Met residues to Gln demonstrated that Met-109 was required for apo-CaM activation of RyR1 but not for Ca(2+)-CaM inhibition of the channel. Furthermore, substitution of Gln for Met-124 increased the apo- and Ca(2+)-CaM concentrations required to regulate RyR1. These results thus identify Met residues critical for the productive association of CaM with RyR1 channels and suggest that oxidation of CaM may contribute to altered regulation of sarcoplasmic reticulum Ca(2+) release during oxidative stress.     

A calmodulin binding domain of RyR increases activation of spontaneous Ca2+ sparks in frog skeletal muscle

The calmodulin C lobe binding region (residues 3614-3643) on the sarcoplasmic reticulum Ca2+ release channel (RyR1) is thought to be a region of contact between subunits within RyR1 homotetramer Ca2+ release channels. To determine whether the 3614-3643 region is a regulatory site/interaction domain within RyR in muscle fibers, we have investigated the effect of a synthetic peptide corresponding to this region (R3614-3643) on Ca2+ sparks in frog skeletal muscle fibers. R3614-3643 (0.2-3.0 microM) promoted the occurrence of Ca2+ sparks in a highly cooperative dose-dependent manner, with a half-maximal activation at 0.47 microM and a maximal increase in frequency of approximately 5-fold. A peptide with a single amino acid substitution within R3614-3643 (L3624D) retained the ability to bind Ca(2+)-free calmodulin but did not increase Ca2+ spark frequency, suggesting that R3614-3643 does not modulate Ca2+ sparks by removal of endogenous calmodulin. Our data support a model in which the calmodulin binding domain of RyR1 modulates channel activity by at least two mechanisms: direct binding of calmodulin as well as interactions with other regions of RyR.     

Kinetics of cardiac RyR channel gating studied at high temporal resolution

Measurements of ryanodine receptor (RyR) activity during dynamic changes of calcium concentration have suggested that RyR has at least four calcium binding sites, and that activation transpires as an increase in the activity within the high open probability H-mode. Binding of several Ca2+ ions within the H-mode should manifest itself in the steady-state RyR activity by the presence of multiple closed times. However, previously only two closed times were detected in the H-mode of RyR activity. Here we recorded steady-state activity of single cardiac RyRs with high temporal resolution and compared it to data simulated under the same conditions using our previously published model of RyR gating. At a 10 kHz resolution, the closed time histograms of both experimental and simulated data had three exponential components. The closed times of simulated data were not significantly different from those obtained experimentally. After filtering at 2 kHz, only two exponential closed time components with time constants not significantly different from those previously published could be detected in both experimental and simulated records. The conformity of the steady-state experimental data to the model derived from the dynamic data provides further support for the idea that RyRs need binding of multiple Ca2+ ions to open.     

An S6 mutation in BK channels reveals beta1 subunit effects on intrinsic and voltage-dependent gating

Large conductance, Ca(2+)- and voltage-activated K(+) (BK) channels are exquisitely regulated to suit their diverse roles in a large variety of physiological processes. BK channels are composed of pore-forming alpha subunits and a family of tissue-specific accessory beta subunits. The smooth muscle-specific beta1 subunit has an essential role in regulating smooth muscle contraction and modulates BK channel steady-state open probability and gating kinetics. Effects of beta1 on channel's gating energetics are not completely understood. One of the difficulties is that it has not yet been possible to measure the effects of beta1 on channel's intrinsic closed-to-open transition (in the absence of voltage sensor activation and Ca(2+) binding) due to the very low open probability in the presence of beta1. In this study, we used a mutation of the alpha subunit (F315Y) that increases channel openings by greater than four orders of magnitude to directly compare channels' intrinsic open probabilities in the presence and absence of the beta1 subunit. Effects of beta1 on steady-state open probabilities of both wild-type alpha and the F315Y mutation were analyzed using the dual allosteric HA model. We found that mouse beta1 has two major effects on channel's gating energetics. beta1 reduces the intrinsic closed-to-open equilibrium that underlies the inhibition of BK channel opening seen in submicromolar Ca(2+). Further, P(O) measurements at limiting slope allow us to infer that beta1 shifts open channel voltage sensor activation to negative membrane potentials, which contributes to enhanced channel opening seen at micromolar Ca(2+) concentrations. Using the F315Y alpha subunit with deletion mutants of beta1, we also demonstrate that the small N- and C-terminal intracellular domains of beta1 play important roles in altering channel's intrinsic opening and voltage sensor activation. In summary, these results demonstrate that beta1 has distinct effects on BK channel intrinsic gating and voltage sensor activation that can be functionally uncoupled by mutations in the intracellular domains.     

Molecular dynamics simulations of NAD+-induced domain closure in horse liver alcohol dehydrogenase

Horse liver alcohol dehydrogenase is a homodimer, the protomer having a coenzyme-binding domain and a catalytic domain. Using all available x-ray structures and 50 ns of molecular dynamics simulations, we investigated the mechanism of NAD+-induced domain closure. When the well-known loop at the domain interface was modeled to its conformation in the closed structure, the NAD+-induced domain closure from the open structure could be simulated with remarkable accuracy. Native interactions in the closed structure between Arg369, Arg47, His51, Ala317, Phe319, and NAD+ were seen to form at different stages during domain closure. Removal of the Arg369 side-chain charge resulted in the loss of the tendency to close, verifying that specific interactions do help drive the domains closed. Further simulations and a careful analysis of x-ray structures suggest that the loop prevents domain closure in the absence of NAD+, and a cooperative mechanism operates between the subunits for domain closure. This cooperative mechanism explains the role of the loop as a block to closure because in the absence of NAD+ it would prevent the occurrence of an unliganded closed subunit when the other subunit closes on NAD+. Simulations that started with one subunit open and one closed supported this.     

The role of calsequestrin, triadin, and junctin in conferring cardiac ryanodine receptor responsiveness to luminal calcium

The level of Ca inside the sarcoplasmic reticulum (SR) is an important determinant of functional activity of the Ca release channel/ryanodine receptor (RyR) in cardiac muscle. However, the molecular basis of RyR regulation by luminal Ca remains largely unknown. In the present study, we investigated the potential role of the cardiac SR luminal auxiliary proteins calsequestrin (CSQ), triadin 1, and junctin in forming the luminal calcium sensor for the cardiac RyR. Recordings of single RyR channels incorporated into lipid bilayers, from either SR vesicle or purified RyR preparations, were performed in the presence of MgATP using Cs+ as the charge carrier. Raising luminal [Ca] from 20 microM to 5 mM increased the open channel probability (Po) of native RyRs in SR vesicles, but not of purified RyRs. Adding CSQ to the luminal side of the purified channels produced no significant changes in Po, nor did it restore the ability of RyRs to respond to luminal Ca. When triadin 1 and junctin were added to the luminal side of purified channels, RyR Po increased significantly; however, the channels still remained unresponsive to changes in luminal [Ca]. In RyRs reassociated with triadin 1 and junctin, adding luminal CSQ produced a significant decrease in activity. After reassociation with all three proteins, RyRs responded to rises of luminal [Ca] by increasing their Po. These results suggest that a complex of CSQ, triadin 1, and junctin confer RyR luminal Ca sensitivity. CSQ apparently serves as a luminal Ca sensor that inhibits the channel at low luminal [Ca], whereas triadin 1 and/or junctin may be required to mediate interactions of CSQ with RyR.     

Shaker potassium channel gating. III: Evaluation of kinetic models for activation

Predictions of different classes of gating models involving identical conformational changes in each of four subunits were compared to the gating behavior of Shaker potassium channels without N-type inactivation. Each model was tested to see if it could simulate the voltage dependence of the steady state open probability, and the kinetics of the single-channel currents, macroscopic ionic currents and macroscopic gating currents using a single set of parameters. Activation schemes based upon four identical single-step activation processes were found to be incompatible with the experimental results, as were those involving a concerted, opening transition. A model where the opening of the channel requires two conformational changes in each of the four subunits can adequately account for the steady state and kinetic behavior of the channel. In this model, the gating in each subunit is independent except for a stabilization of the open state when all four subunits are activated, and an unstable closed conformation that the channel enters after opening. A small amount of negative cooperativity between the subunits must be added to account quantitatively for the dependence of the activation time course on holding voltage.     

Binding of a gating modifier toxin induces intersubunit cooperativity early in the Shaker K channel's activation pathway

Potassium currents from voltage-gated Shaker K channels activate with a sigmoid rise. The degree of sigmoidicity in channel opening kinetics confirms that each subunit of the homotetrameric Shaker channel undergoes more than one conformational change before the channel opens. We have examined effects of two externally applied gating modifiers that reduce the sigmoidicity of channel opening. A toxin from gastropod mucus, 6-bromo-2-mercaptotryptamine (BrMT), and divalent zinc are both found to slow the same conformational changes early in Shaker's activation pathway. Sigmoidicity measurements suggest that zinc slows a conformational change independently in each channel subunit. Analysis of activation in BrMT reveals cooperativity among subunits during these same early steps. A lack of competition with either agitoxin or tetraethylammonium indicates that BrMT binds channel subunits outside of the external pore region in an allosterically cooperative fashion. Simulations including negatively cooperative BrMT binding account for its ability to induce gating cooperativity during activation. We conclude that cooperativity among K channel subunits can be greatly altered by experimental conditions.     

Ablation of skeletal muscle triadin impairs FKBP12/RyR1 channel interactions essential for maintaining resting cytoplasmic Ca2+

Previously, we have shown that lack of expression of triadins in skeletal muscle cells results in significant increase of myoplasmic resting free Ca(2+) ([Ca(2+)](rest)), suggesting a role for triadins in modulating global intracellular Ca(2+) homeostasis. To understand this mechanism, we study here how triadin alters [Ca(2+)](rest), Ca(2+) release, and Ca(2+) entry pathways using a combination of Ca(2+) microelectrodes, channels reconstituted in bilayer lipid membranes (BLM), Ca(2+), and Mn(2+) imaging analyses of myotubes and RyR1 channels obtained from triadin-null mice. Unlike WT cells, triadin-null myotubes had chronically elevated [Ca(2+)](rest) that was sensitive to inhibition with ryanodine, suggesting that triadin-null cells have increased basal RyR1 activity. Consistently, BLM studies indicate that, unlike WT-RyR1, triadin-null channels more frequently display atypical gating behavior with multiple and stable subconductance states. Accordingly, pulldown analysis and fluorescent FKBP12 binding studies in triadin-null muscles revealed a significant impairment of the FKBP12/RyR1 interaction. Mn(2+) quench rates under resting conditions indicate that triadin-null cells also have higher Ca(2+) entry rates and lower sarcoplasmic reticulum Ca(2+) load than WT cells. Overexpression of FKBP12.6 reverted the null phenotype, reducing resting Ca(2+) entry, recovering sarcoplasmic reticulum Ca(2+) content levels, and restoring near normal [Ca(2+)](rest). Exogenous FKBP12.6 also reduced the RyR1 channel P(o) but did not rescue subconductance behavior. In contrast, FKBP12 neither reduced P(o) nor recovered multiple subconductance gating. These data suggest that elevated [Ca(2+)](rest) in triadin-null myotubes is primarily driven by dysregulated RyR1 channel activity that results in part from impaired FKBP12/RyR1 functional interactions and a secondary increased Ca(2+) entry at rest.     

Excitation-contraction coupling gain and cooperativity of the cardiac ryanodine receptor: a modeling approach

During calcium-induced calcium-release, the ryanodine receptor (RyR) opens and releases large amounts of calcium from the sarcoplasmic reticulum into the cytoplasm of the myocyte. Recent experiments have suggested that cooperativity between the four monomers comprising the RyR plays an important role in the dynamics of the overall receptor. Furthermore, this cooperativity can be affected by the binding of FK506 binding protein, and hence, modulated by adrenergic stimulation through the phosphorylating action of protein kinase A. This has important implications for heart failure, where it has been hypothesized that RyR hyperphosphorylation, resulting in a loss of cooperativity, can lead to a persistent leak and a reduced sarcoplasmic-reticula content. In this study, we construct a theoretical model that examines the cooperativity via the assumption of an allosteric interaction between the four subunits. We find that the level of cooperativity, regulated by the binding of FK506 binding-protein, can have a dramatic effect on the excitation-contraction coupling gain and that this gain exhibits a clear maximum. These findings are compared to currently available data from different species and allows for an evaluation of the aforementioned heart-failure scenario.     

Markovian models of low and high activity levels of cardiac ryanodine receptors

The modal gating behavior of single sheep cardiac sarcoplasmic reticulum (SR) Ca2+-release/ryanodine receptor (RyR) channels was assessed. We find that the gating of RyR channels spontaneously shifts between high (H) and low (L) levels of activity and inactive periods where no channel openings are detected (I). Moreover, we find that there is evidence for multiple gating modes within H activity, which we term H1 and H2 mode. Our results demonstrate that the underlying mechanisms regulating gating are similar in native and purified channels. Dwell-time distributions of L activity were best fitted by three open and five closed significant exponential components whereas dwell-time distributions of H1 activity were best fitted by two to three open and four closed significant exponential components. Increases in cytosolic [Ca2+] cause an increase in open probability (Po) within L activity and an increase in the probability of occurrence of H activity. Open lifetime distributions within L activity were Ca2+ independent whereas open lifetime distributions within H activity were Ca2+ dependent. This study is the first attempt to estimate RyR single-channel kinetic parameters from sequences of idealized dwell-times and to develop kinetic models of RyR gating using the criterion of maximum likelihood. We propose distinct kinetic schemes for L, H1, and H2 activity that describe the major features of sheep cardiac RyR channel gating at these levels of activity.     

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.     

Stoichiometry of altered hERG1 channel gating by small molecule activators

Voltage-gated K⁺ channels are tetramers formed by coassembly of four identical or highly related subunits. All four subunits contribute to formation of the selectivity filter, the narrowest region of the channel pore which determines K⁺ selective conductance. In some K⁺ channels, the selectivity filter can undergo a conformational change to reduce K⁺ flux by a mechanism called C-type inactivation. In human ether-a-go-go–related gene 1 (hERG1) K⁺ channels, C-type inactivation is allosterically inhibited by ICA-105574, a substituted benzamide. PD-118057, a 2-(phenylamino) benzoic acid, alters selectivity filter gating to enhance open probability of channels. Both compounds bind to a hydrophobic pocket located between adjacent hERG1 subunits. Accordingly, a homotetrameric channel contains four identical activator binding sites. Here we determine the number of binding sites required for maximal drug effect and determine the role of subunit interactions in the modulation of hERG1 gating by these compounds. Concatenated tetramers were constructed to contain a variable number (zero to four) of wild-type and mutant hERG1 subunits, either L646E to inhibit PD-118057 binding or F557L to inhibit ICA-105574 binding. Enhancement of hERG1 channel current magnitude by PD-118057 and attenuated inactivation by ICA-105574 were mediated by cooperative subunit interactions. Maximal effects of the both compounds required the presence of all four binding sites. Understanding how hERG1 agonists allosterically modify channel gating may facilitate mechanism-based drug design of novel agents for treatment of long QT syndrome.