Calcium-dependent electrostatic control of anion access to the pore of the calcium-activated chloride channel TMEM16A
Abstract
TMEM16A is a ligand-gated anion channel that is activated by intracellular Ca2+. This channel comprises two independent pores and closely apposed Ca2+ binding sites that are contained within each subunit of a homodimeric protein. Previously we characterized the influence of positively charged pore-lining residues on anion conduction (Paulino et al., 2017a). Here, we demonstrate the electrostatic control of permeation by the bound calcium ions in mouse TMEM16A using electrophysiology and Poisson-Boltzmann calculations. The currents of constitutively active mutants lose their outward rectification as a function of Ca2+ concentration due to the alleviation of energy barriers for anion conduction. This phenomenon originates from Coulombic interactions between the bound Ca2+ and permeating anions and thus demonstrates that an electrostatic gate imposed by the vacant binding site present in the sterically open pore, is released by Ca2+ binding to enable an otherwise sub-conductive pore to conduct with full capacity.
https://doi.org/10.7554/eLife.39122.001Introduction
The calcium-activated chloride channel TMEM16A is part of a large family of membrane proteins that encompasses ion channels and lipid scramblases with a common conserved molecular architecture (Brunner et al., 2016; Caputo et al., 2008; Picollo et al., 2015; Schroeder et al., 2008; Terashima et al., 2013; Whitlock and Hartzell, 2017; Yang et al., 2008). With respect to their fold, the TMEM16 family is related to mechanosensitive channels of the OSCA and TMC families (Ballesteros et al., 2018; Pan et al., 2018; Zhang et al., 2018). TMEM16A is widely expressed and contributes to important physiological processes including the transport of chloride across epithelia and the control of electrical signal transduction in smooth muscle and certain neurons (Huang et al., 2012; Oh and Jung, 2016; Scudieri et al., 2012). Recent investigations have defined the structural basis for ion permeation and gating in TMEM16A and revealed features that distinguish this ion channel from other homologues working as lipid scramblases (Brunner et al., 2014; Dang et al., 2017; Paulino et al., 2017a; Paulino et al., 2017b). TMEM16A harbors two ion conduction pores, each contained within a single subunit of a homodimeric protein (Figure 1A). Both pores function independently and are activated by the binding of two Ca2+ ions to a site embedded within the membrane-inserted part of each subunit close to the ion conduction path (Jeng et al., 2016; Lim et al., 2016). In the open conformation, anions access the narrow neck of an hourglass-shaped pore via water-filled vestibules from the extra- and intracellular sides and permeate through the constricted part presumably after shedding most of their coordinating water molecules (Betto et al., 2014; Dang et al., 2017; Ni et al., 2014; Paulino et al., 2017a; Qu and Hartzell, 2000). During this process, the energetic penalty for dehydration is surmounted by positive charges placed on both sides of the neck (Paulino et al., 2017b). In the apo conformation, the Ca2+ binding site is accessible from the cytoplasm (Paulino et al., 2017a). The binding of Ca2+ favors conformational changes in the pore-lining helix α6, which provides polar and acidic sidechains that coordinate the bound cations and thus directly couples ligand binding to pore opening (Paulino et al., 2017a; Peters et al., 2018). The subsequent closure of the aqueous access pathway to the ligand binding site buries the bound calcium ions within the transmembrane electric field, which mechanistically accounts for the observed voltage dependence of activation by modulating the binding affinity of Ca2+ (Arreola et al., 1996; Brunner et al., 2014; Xiao et al., 2011). The movement of the helix relays a conformational change towards the neck to release a steric barrier, which acts as a gate in the closed conformation of the channel (Paulino et al., 2017a). Additionally, the bound Ca2+ ions change the charge distribution at the wide intracellular vestibule thereby removing a second, electrostatic barrier that specifically hinders the access of anions in the ligand-free state (Paulino et al., 2017a). Here we have used electrophysiology and Poisson-Boltzmann calculations to characterize the electrostatic control of anion permeation by the bound Ca2+ ions. Our results reveal a strong and favorable Coulombic interaction of permeant anions with the positively charged Ca2+ at the intracellular pore entry. These long-range interactions underlie an electrostatic gating mechanism that acts through space and that is synergistic with the opening of a steric gate, allowing an otherwise sub-conductive pore to conduct with full capacity.
Results
Ca2+ binding to the transmembrane site alleviates energy barriers for anion conduction
We have previously used a simple rate model to characterize ion conduction in TMEM16A (Läuger, 1973; Paulino et al., 2017b). In this model, a central small energy barrier originating from the diffusion of an anion across the narrow neck is sandwiched between two larger barriers resulting from the desolvation of the anion upon entering the constricted part of the channel from the aqueous vestibules located on either side of the pore (Figure 1A, Figure 1—figure supplement 1A,B). The diffusion path does not contain deep energy wells and the model does not account for saturation of the pore, which is generally consistent with the high KM for chloride conduction. This model allowed us to phenomenologically interpret the effect of mutations of positively charged residues on current-voltage (I-V) relationships. Whereas single mutations in the wide vestibules do not exert a recognizable effect on conduction, mutations at the border of the neck region result in a pronounced rectification of currents whose shape depends on the location of the altered residues (Paulino et al., 2017b). The rectification is a consequence of both pore occupancy and the rate of barrier crossing at the applied potential. Here we use a similar analysis to investigate the influence of the bound calcium ions on anion conduction by characterizing the properties of the mutant G644P. In this mutant, the replacement of a flexible glycine at a hinge in α6 (Gly 644, Figure 1B) with a rigid proline increases the potency of Ca2+ and concomitantly results in constitutive activity (Paulino et al., 2017a). The basal current of this mutant is highly outwardly rectifying but it progressively loses its rectification at increasing Ca2+ concentrations until it becomes pseudo-linear in saturating conditions (Figure 1C, Figure 1—figure supplement 1C). This differs from the instantaneous currents of WT, which are linear in the entire Ca2+ concentration range (Figure 1—figure supplement 2). When analyzed with our described minimal model of ion permeation, this effect seems to originate from the alleviation of local energy barriers for anion conduction at the intracellular entrance and the center of the pore as a function of Ca2+ concentration (Figure 1D). This is consistent with a change in long-range interactions affecting permeant anions entering and traversing the narrow neck region upon binding of the positively charged Ca2+. The lowering of the energy barriers can be described by a binding isotherm with EC50 values of 60 nM and 170 nM, and Hill coefficients of 2.1 and 1.7 at 80 mV for the inner and the central barriers respectively (Figure 1E). As the rate constants are not strongly coupled, the similarity of their respective EC50’s and their consistency with data from steady-state current responses provide independent evidence for a saturable effect associated with the binding of Ca2+ to the transmembrane site identified in the structure.
Besides G644P, we have characterized the mutant Q649A that also causes constitutive activity. Located about one helix turn below Gly 644, the residue is not directly involved in Ca2+ coordination (Figure 1B). Like G644P, this mutant displays a left-shifted EC50 and basal current that is equally outwardly rectifying (Figure 1—figure supplement 3A). Since the current magnitude of this mutant is larger and the kinetics of irreversible rundown is slower compared to G644P, we were able to examine the influence of the membrane potential on Ca2+ binding by recording instantaneous currents with pre-pulses at different voltages on the same patch (Figure 1F and Figure 1—figure supplement 3B—D). At the membrane potentials tested, we observed a similar effect of Ca2+ on the energetics of anion conduction as in G644P (Figure 1C–F) and as with steady-state concentration-responses, the potency of Ca2+ is increased upon depolarization (Figure 1F–G, Figure 1—figure supplement 3B–D). Therefore, these results provide further evidence that Ca2+ alleviates local energy barriers for anion conduction at the intracellular entrance and the middle of the pore by binding to the transmembrane site.
The effect of the competitive antagonist Mg2+, and the trivalent cation Gd3+ on conduction
To confirm that the charge of Ca2+, rather than its ability to promote gating-associated conformational changes, accounts for the observed effect, we characterized the I-V relations of G644P in the presence of the divalent cation Mg2+ (Figure 2). Although Mg2+ is incapable of activating WT even at high concentrations (Figure 2—figure supplement 1A), it acts as a low affinity competitive antagonist that lowers the potency of Ca2+, indicating that it might occupy the Ca2+ binding site (Ni et al., 2014). This is confirmed in the mutant G644P, for which we observed a concentration-dependent reduction of the outward rectification in the presence of Mg2+ (Figure 2A). Mg2+ thus appears to ‘activate’ G644P, by primarily increasing its conductance but not its open probability (Figure 2—figure supplement 1B and C). Although we cannot exclude a non-specific effect of Mg2+ due to the screening of surface charges, this is unlikely, as alanine mutations of positive charges in the wide vestibule have previously been shown to exert little effect on conduction (Paulino et al., 2017b). Interestingly, although Mg2+ exerts a qualitatively similar effect on G644P as Ca2+, the current remains rectifying even at saturating concentrations (Figure 2B). Moreover, in contrast to Ca2+, the activation by Mg2+ proceeds with a Hill coefficient close to unity and a maximum conduction rate at the intracellular entrance (σβ) that saturates at an intermediate value, suggesting that the binding site is likely occupied by a single Mg2+ ion (Figure 2B and C). Taken together, our data show three discrete levels of conductance that correspond to the possible occupancies of the binding site with zero, one and two divalent cations bound (Figure 2C), consistent with a scenario where sequential Ca2+ binding lowers the energy barriers for anion conduction in a stepwise manner, as predicted by a simple Coulombic interaction. The model also describes the effect upon addition of the trivalent cation Gd3+, which strongly increases the conductance at saturating concentrations even beyond the level induced by Ca2+. This is likely a consequence of the increased positive charge density in the binding site and is manifested in the observed inward rectification of currents (Figure 2—figure supplement 2).
Ca2+ increases anion conductance by an electrostatic mechanism
Our results suggest that Ca2+ and other di- and trivalent cations gate anion conduction in both constitutively active mutants by an electrostatic mechanism that involves neutralization of the negative charges in the vacant binding site. To gain further insight into this process, we investigated the effect of mutations of residues at the Ca2+ binding site introduced in the G644P background and initially focused on E654Q located on α6. Whereas on a WT background, E654Q exhibits the most severe phenotype among binding site mutants and does not show any activity even at high Ca2+ concentrations (Brunner et al., 2014; Lim et al., 2016; Tien et al., 2014), we observed basal and outwardly-rectifying currents in the mutant G644P/E654Q (PQ) whose steady-state current increases in response to Ca2+ addition (Figure 3—figure supplement 1A). We examined the instantaneous I-V relations at increasing Ca2+ concentrations and found a moderate concentration-dependent reduction of rectification that eventually saturates (Figure 3A). The much lower conduction rate of the PQ mutant compared to G644P at saturating Ca2+ concentrations and a Hill coefficient close to one both suggest that in this mutant the binding of a single calcium ion to the transmembrane site promotes anion permeation, although less efficiently than the binding of two calcium ions in G644P, similar to the effect of Mg2+ (Figures 2 and 3A).
We extended our investigation on Glu 654 and mutated this residue to an arginine (PR) to add a further positive charge to the binding site. Similar to PQ, this mutant likely binds only one calcium ion as the Hill coefficient is close to unity and the I-V relation retains its strong rectification even at saturating Ca2+ concentrations (Figure 3B, Figure 3—figure supplement 1B). However, consistent with the hypothesis that a more positive potential at the binding site lowers the energy barrier for anion permeation, the rate of chloride conduction at the intracellular entrance at saturating Ca2+ concentrations is significantly higher than that observed for PQ (Figure 3B). Next, we removed an additional negative charge in the binding site of the PR mutant by exchanging Glu 702 with Gln (PR2Q), which, similar to the equivalent mutation on the WT background (Brunner et al., 2014), reduces the potency of Ca2+ (Figure 3C, Figure 3—figure supplement 1C). As expected, the rate of chloride conduction at the intracellular pore entrance in the triple mutant PR2Q is increased compared to the double mutant PR both in the absence of Ca2+ and at saturating Ca2+ concentrations (Figure 3C). Together these observations consolidate the notion of the dependence of anion conduction on the electrostatic potential around the binding site.
For electrostatic interactions, we expect the contribution of individual charges in the binding site to conduction to be additive. This hypothesis was tested by analyzing the effect of the removal of a negative charge in the binding site mutant E705Q in the PQ (PQ5Q) and PR (PR5Q) backgrounds. At maximum Ca2+ occupancy, E705Q increases the rate of anion conduction to the same extent in both PQ5Q and PR5Q irrespectively of its background (Figure 3D,E), which emphasizes the orthogonal contributions of individual mutations (Figure 3F). The same mutations also exhibit a similar degree of additivity in the Ca2+-free state (Figure 3G), further supporting an electrostatic mechanism for the functional interaction between the Ca2+ binding site and the anion conduction path.
Electrostatic interactions are evident in calculations
To further strengthen our proposal of an electrostatic interaction between the bound Ca2+ and permeating anions, we investigated the effect of the binding site mutations and the altered Ca2+ occupancy on the electrostatic profile of the pore using Poisson-Boltzmann calculations on the Ca2+-bound mouse TMEM16A structure (Figure 4A). Consistent with the quantal effect of Ca2+ on anion conduction, removal of a single calcium ion from the binding site partially decreases the positive potential observed in the fully bound structure (Figure 4B). We also found that stepwise removal of charges in the binding site exerts progressively larger effects according to the number of charges neutralized (Figure 4C), as observed in our experiments. Although we cannot exclude the possibility that the detailed geometric arrangement of the Ca2+ binding site and nearby helices might be influenced by the introduced mutations, local changes in backbone and sidechain conformations are unlikely to drastically affect the calculated electrostatic potential given the long-range nature of Coulombic interactions. Because several mutants described above seem to bind only one calcium ion (Figure 3 and Figure 3—figure supplement 1), it is possible to estimate from our experimental data the effective dielectric constant of the intracellular vestibule that propagates the electric field originating from the binding site to the pore. When the above mutants are analyzed collectively, we observed an inverse and linear relationship between the energy barrier for chloride conduction at the intracellular pore entrance and the valence of the binding site (Figure 4D). Assuming a Coulombic potential and using the distance measured from the Ca2+-bound structure, we estimated an effective dielectric constant of around 160 in the vacant and 90 for the single Ca2+-bound state. Although potentially inaccurate and physically unreasonable due to experimental limitations, these values indicate bulk solution-like aqueous properties that are consistent with the large water-filled intracellular vestibule observed in the Ca2+-bound structure (Figure 4D). This also suggests that the vestibule is functionally shielded from the membrane.
Bidirectional electrostatic interaction between the binding site and the intracellular pore entrance
In previous experiments, we have demonstrated that the bound Ca2+ ions influence ion permeation by affecting the electrostatics at the narrow neck region via long-range Coulombic interactions. Reciprocally, changes in the charge distribution at the neck should affect Ca2+ binding. To test this hypothesis, we have investigated the effect of mutations of two residues, Lys 588 and Lys 645 (Figure 5A) located at the boundary between the intracellular vestibule and the neck, which lower the energy barrier for chloride conduction by contributing to the positive electrostatics of the pore (Paulino et al., 2017b). As expected, we found an increase in the potency of Ca2+ in neutralizing mutations of the respective residues to serine (Figure 5B,C), which is further enhanced upon the reversal of the charge by mutations to glutamate (Figure 5D,E). Assuming that only the binding affinity of Ca2+ is affected, we observed a linear relationship between the change in binding energy and the valence at the intracellular pore entrance (Figure 5F), with a stronger sensitivity for Lys 645, which is in closer proximity to the Ca2+ binding site (Figure 5A) and thus might experience less solvent screening. In summary, the described effects further confirm the long-range interactions between the Ca2+ binding site and the narrow neck region, which underlie the electrostatic control of ion permeation.
Discussion
Our study demonstrates the strong electrostatic influence of bound Ca2+ ions on anion conduction in the Ca2+-activated Cl- channel TMEM16A. A direct effect of the ligand on the permeation properties of the open state is unique among ligand-gated ion channels and it is a consequence of the high positive charge density of Ca2+ and the location of its binding site in the immediate vicinity of the permeation path. The described effect is synergistic with a conformational change upon Ca2+ binding that opens a steric gate in the narrow part of the pore (Paulino et al., 2017a). In the absence of Ca2+, even in a sterically released conformation, which in this study was stabilized by mutations, an electrostatic gate imposed by the vacant binding site would impede the access of chloride ions to the intracellular vestibule (Figure 6). The binding of Ca2+ releases this second gate by neutralizing the negative charges, enabling the open channel to conduct with full capacity (Figure 6). Our findings thus provide novel mechanistic insight into the regulation of anion permeation in Ca2+-activated Cl- channels of the TMEM16 family and demonstrate that these channels are dually gated by both steric and electrostatic mechanisms. Moreover, they also provide a conceivable explanation for how the oppositely charged substrates, in this case permeating anions and Ca2+ ions, access their respective site of action via the same intracellular vestibule (Figure 6). The observation that the conductance of the open channel can be electrostatically regulated by Ca2+ in a stepwise manner also offers a unique opportunity to examine the relative contribution of open states with partial Ca2+ occupancy during activation. The currents in such partially occupied states would be rectifying and thus be readily identifiable in the current-voltage relationships. In WT, the instantaneous I-V plots are all approximately linear and their shape is only mildly sensitive to Ca2+ in the entire activation range (Figure 1—figure supplement 2). Since these plots reflect the weighted sum of the I-V relations of the open states with different Ca2+ occupancy, our results indicate that, during activation, the major conducting state in WT is the open state with two Ca2+ ions bound, whereas the contribution of open states with submaximal Ca2+ occupancy is minimal (Figure 1—figure supplement 2 and Figure 6—figure supplement 1A,B). A similar but less pronounced distribution of states is also evident in the constitutively active mutant G644P (Figure 6—figure supplement 1C), where the apo and the singly occupied channel make a larger contribution to conduction at low Ca2+ concentrations. These observations thus provide direct evidence for the high cooperativity of TMEM16A activation.
Aside from their impact on membrane biophysics, our findings are particularly relevant for the development of allosteric activators of TMEM16A that may be beneficial for the treatment of obstructive airway diseases such as Cystic Fibrosis (Li et al., 2017; Mall and Galietta, 2015). Although we anticipate that small molecules have the potential to increase the low open probability of the channel at low Ca2+ concentration, this might not suffice to promote robust Cl- efflux required to compensate for the loss of the chloride channel CFTR unless at the same time the electrostatic barrier at the intracellular pore entry is released, which would be the case for compounds with a positive net-charge or ones that at the same time increase the potency of Ca2+.
Materials and methods
Molecular biology and cell culture
Request a detailed protocolHEK293T cells (ATCC CRL-1573) were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Sigma-Aldrich) supplemented with 10 U/ml penicillin, 0.1 mg/ml streptomycin (Sigma-Aldrich), 2 mM L-glutamine (Sigma-Aldrich), and 10% FBS (Sigma-Aldrich) in a humidified atmosphere containing 5% CO2 at 37 °C. HEK293T cells were transfected with 3 μg DNA per 6 cm Petri dish using the calcium phosphate co-precipitation method, and were used within 24–96 hr after transfection. Mutants were generated with a modified QuikChange method (Zheng et al., 2004) using the wild-type mouse TMEM16A(ac) as the template. Double and triple mutants were generated sequentially. All constructs were verified by sequencing.
Electrophysiology
Request a detailed protocolRecordings were performed on inside-out patches excised from cells expressing the construct of interest. Recording pipettes were pulled from borosilicate glass capillaries (O.D. 1.5, I.D. 0.86, Sutter Instrument) and were fire-polished with a microforge (Narishige) before use. Pipette resistance was typically 3–8 MΩ when standard recording solutions were used. Seal resistance was typically 4 GΩ or higher. Voltage-clamp recordings were performed using Axopatch 200B and Digidata 1550 (Molecular devices). Analogue signals were filtered through the in-built 4-pole Bessel filter at 5 kHz and were digitized at 10–20 kHz. Data acquisition was performed using Clampex 10.6 (Molecular devices). Solution exchange was achieved using a double-barreled theta glass pipette mounted on an ultra-high speed piezo-driven stepper (Siskiyou). Liquid junction potential was found to be consistently negligible given the ionic composition of the solutions and was therefore not corrected. All experiments were performed at 20 °C.
Recordings were performed under symmetrical ionic conditions. Stock solution with Ca2+-EGTA contained 150 mM NaCl, 5.99 mM Ca(OH)2, 5 mM EGTA and 10 mM HEPES at pH 7.40. Stock solution with EGTA contained 150 mM NaCl, 5 mM EGTA and 10 mM HEPES at pH 7.40. Free Ca2+ concentrations were adjusted by mixing the stock solutions at the required ratios, which were calculated using the WEBMAXC program (http://web.stanford.edu/~cpatton/webmaxcS.htm). Recording pipettes were filled with stock solution with Ca2+-EGTA, which has a free Ca2+ concentration of 1 mM.
For Mg2+ experiments, stock solution with Mg2+-EGTA contained 75 mM MgCl2, 23 mM (NMDG)2SO4, 5 mM EGTA and 10 mM HEPES at pH 7.40. Free Mg2+ concentrations were adjusted by mixing the stock solutions with Mg2+-EGTA and EGTA at the required ratios. For Gd3+ experiments, the final Gd3+ concentrations were reconstituted from a 1 M GdCl3 stock in a solution containing 150 mM NaCl, 10 mM HEPES at pH 7.40 with no added Ca2+ and metal chelators. In all cases, solutions were prepared in ultrapure molecular biology grade water having a resistivity of 15–18 MΩ cm at 25 °C (ELGA or Millipore).
Concentration-response relations were constructed from steady-state current responses recorded using a rundown-correction protocol as described previously (Lim et al., 2016; Paulino et al., 2017a). For instantaneous I-V relations, current responses were measured at the time points (within milliseconds) where the capacitive current at 0 mV, the membrane potential at which the current reverses, has decayed. The sequence of the applied voltage steps was from −100 to +120 mV to minimize the effect of rundown on the smaller inward current (Figure 1—figure supplements 1–3, Figure 2—figure supplement 1). In a typical recording, the magnitude of the steady-state current at +80 mV would have decayed irreversibly by 10–15% by the time the test pulse at +120 mV is applied. To correct for this effect, the I-V plots were normalized to the current amplitude of the pre-pulses (Paulino et al., 2017b) before their normalization (I/I+120). When used judiciously, this procedure allows one to recover the true I-V relation and to correct for current fluctuation. Before and after each recording at the Ca2+ concentration of interest, we also recorded in zero Ca2+ to ensure that the quality of the patch did not deteriorate over time. Background subtraction was performed for constructs that do not display visible basal activity but was not possible for the mutants that do. In all cases, leaky patches were discarded.
A similar procedure was used for experiments with Mg2+ and Gd3+. In contrast to the divalent cations Ca2+ and Mg2+, Gd3+ binding appears to be biphasic consisting of a reversible phase and an irreversible component that both affect the I-V relation of the instantaneous current (Figure 2—figure supplement 2A). In order to titrate the second reversible site, we first incubated the patch with a submaximal concentration of Gd3+ (0.4 mM) to saturate the irreversible binding site. We then recorded at the test concentration of Gd3+ to obtain the corresponding instantaneous I-V relations (Figure 2—figure supplement 2B,C).
Data analysis
Request a detailed protocolConcentration-response data were fitted to the Hill equation
where is the normalized current, is the concentration at which equals 0.5 and is the Hill coefficient. For the case of one binding site, h was set to equal 1. The voltage dependence of was fitted to
where is the valence of Ca2+, is the electrical distance, is the membrane potential, , and have their usual meanings, and is the value of when = 0.
I-V data were fitted to a minimal permeation model that accounts for the most fundamental biophysical behavior of mouse TMEM16A as described previously (Läuger, 1973; Paulino et al., 2017b),
where is the current, is the number of barriers, and are the intracellular and extracellular concentrations of the charge carrier, is the valence of Cl- and , , and are defined as above. is a proportionality factor where is the value of when = 0 and is a proportionality coefficient that has a dimension of volume. and are respectively the rate of barrier crossing at the middle and the innermost barriers relative to that at the outermost barrier (),
where is the rate of barrier crossing at the middle barrier and at the innermost barrier in the absence of voltage. Thus, for a linear voltage drop as is the case in our model, and are intrinsically voltage-independent and represent the relative rates when = 0. appears to be voltage-dependent in Figure 1G because is a function of Ca2+ binding whose EC50 exhibits voltage dependence.
In this model, current rectification is governed by the two parameters and . Values of and that alter the symmetry of the energy profile (as exemplified in Figure 1—figure supplement 1A ) results in asymmetric I-V relations depending on the directionality of symmetry mismatch (Paulino et al., 2017b). We have previously shown that the model describes the behavior of WT best when the number of barriers equals 3 (Paulino et al., 2017b), which was used as a fixed parameter. This leaves , and the amplitude factor as the only free parameters to be fitted. In this study, we focused on because its best-fit value is generally better defined and it increases monotonically with the asymmetry of the I-V relation (rectification index (RI), I-100/+120), which allows a straightforward comparison with the I-V data.
The Ca2+ dependence of and was fitted to
where the subscript indicates or . We observed for mutants G644P/E654Q, G644P/E654R, G644P/E654Q/E705Q and G644P/E654R/E705Q a low affinity decay of at high Ca2+ concentrations (Figure 3A, B, D, E). Since this decay might be related to processes not intrinsic to the transmembrane binding site (Lim et al., 2016), only the high affinity phase was analyzed. The best-fit values of at zero and saturating Ca2+ concentrations were used to calculate , the difference between the activation energy at the innermost barrier relative to that of the outermost, using
The relationship between and the valence of the binding site () at different discrete Ca2+ occupancies (0, 1 or 2) was fitted to the Coulombic potential
where is the Avogadro’s number, is the electronic charge, is the permittivity of vacuum, is the relative permittivity of the medium, and is the distance between the binding site and the location of the innermost barrier. was considered 0 when the five acidic residues (Glu 654, 702, 705, 734 and Asp 738) in the binding site are neutralized. was estimated by setting to the measured distance between the center of the two calcium ions in the binding site and the sidechain oxygen atom of Ser 592 (= 13.6 Å) in the structure of mouse TMEM16A in the Ca2+-bound state.
An analogous expression was used to account for the changes in the energy of Ca2+ binding as a function of valence of the intracellular pore entrance,
where is the observed change in the free energy of Ca2+ binding when only affinity is affected, is the valence at the intracellular pore entrance and is the valence of Ca2+. was considered 0 for the WT protein. was estimated by setting to the measured distance between the center of the two calcium ions in the binding site and the sidechain nitrogen atoms of Lys 588 or 645 (= 11.9 Å or 10.6 Å) in the structure of mouse TMEM16A in the Ca2+-bound state. was calculated using
where is the of the background (WT or E702Q) on which the mutation of interest was constructed. We have constructed the K588E and K645E mutations on the E702Q background to improve the expression of the constructs. is equivalent to when the affinities of all binding reactions are affected equally without affecting conformational transitions, which is a general property of linked equilibria.
To estimate the contribution of open states with various Ca2+ occupancy to the ensemble current, the I-V data for WT were fitted to a weighted sum of I-V curves corresponding to zero, one and two calcium ions bound at the binding site,
where , and are in the form of Equation 2 with parameters from G644P at zero Ca2+, G644P/E654Q at saturating Ca2+ and G644P at saturating Ca2+ respectively. Since WT does not display noticeable basal activity, the contribution of current corresponding to zero occupancy was considered negligible and the term was omitted.
Data analysis was performed using Clampfit 10.6 (Molecular devices), Excel (Microsoft), and Prism 5 and 6 (GraphPad). Numerical calculations were performed using the NumPy package (http://www.numpy.org/). Data are presented as mean ± s.e.m. Fitted parameters are plotted and are reported as best-fit value ± 95% confidence interval.
Poisson-Boltzmann calculations
Request a detailed protocolThe electrostatic potential along the pore, identified using Hole (Smart et al., 1996), was calculated by solving the linearized Poisson-Boltzmann equation in CHARMM (Im et al., 1998; Brooks et al., 1983) on a 240 Å x 240 Å x 260 Å grid (1 Å grid spacing) followed by focusing on a 160 Å x 160 Å x 160 Å grid (0.5 Å grid spacing). Partial protein charges were derived from the CHARMM36 all-hydrogen atom force field (MacKerell et al., 1998). Hydrogen positions were generated in CHARMM. In silico mutagenesis was carried out in CHARMM. The protein was assigned a dielectric constant () of 2. The membrane was represented as a 35-Å-thick slab (= 2). A 5-Å-thick slab was included on each side of the membrane to account for the headgroup region (= 30). The bulk solvent on either side of the membrane and the solvent-filled conduit were represented as an aqueous medium (= 80) containing 150 mM mobile monovalent ions. Electrostatic calculations were carried out in the absence of an applied voltage. Although the structure of the protein used for calculations might not display a fully activated conformation, conformational changes are assumed to be small and would mainly affect the narrow neck region for which the dielectric constant is uncertain (Paulino et al., 2017a). The region most relevant for this study concerns the wide intracellular vestibule, which was assigned bulk-like dielectric properties. In calculations of mutants, positively charged residues were introduced as protonated histidines instead of arginines for modeling purposes.
Data availability
All data generated or analysed during this study are included in the manuscript.
References
-
Activation of calcium-dependent chloride channels in rat parotid acinar cellsThe Journal of General Physiology 108:35–47.https://doi.org/10.1085/jgp.108.1.35
-
Interactions between permeation and gating in the TMEM16B/anoctamin2 calcium-activated chloride channelThe Journal of General Physiology 143:703–718.https://doi.org/10.1085/jgp.201411182
-
CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry 4:187–217.https://doi.org/10.1002/jcc.540040211
-
Structural basis for phospholipid scrambling in the TMEM16 familyCurrent Opinion in Structural Biology 39:61–70.https://doi.org/10.1016/j.sbi.2016.05.020
-
Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equationComputer Physics Communications 111:59–75.https://doi.org/10.1016/S0010-4655(98)00016-2
-
Independent activation of distinct pores in dimeric TMEM16A channelsThe Journal of General Physiology 148:393–404.https://doi.org/10.1085/jgp.201611651
-
Ion transport through pores: a rate-theory analysisBiochimica et Biophysica Acta (BBA) - Biomembranes 311:423–441.https://doi.org/10.1016/0005-2736(73)90323-4
-
Bypassing CFTR dysfunction in cystic fibrosis with alternative pathways for anion transportCurrent Opinion in Pharmacology 34:91–97.https://doi.org/10.1016/j.coph.2017.10.002
-
Independent activation of ion conduction pores in the double-barreled calcium-activated chloride channel TMEM16AThe Journal of General Physiology 148:375–392.https://doi.org/10.1085/jgp.201611650
-
All-atom empirical potential for molecular modeling and dynamics studies of proteinsThe Journal of Physical Chemistry B 102:3586–3616.https://doi.org/10.1021/jp973084f
-
Targeting ion channels in cystic fibrosisJournal of Cystic Fibrosis 14:561–570.https://doi.org/10.1016/j.jcf.2015.06.002
-
Cellular functions of TMEM16/anoctaminPflügers Archiv - European Journal of Physiology 468:443–453.https://doi.org/10.1007/s00424-016-1790-0
-
TMEM16 proteins: unknown structure and confusing functionsJournal of Molecular Biology 427:94–105.https://doi.org/10.1016/j.jmb.2014.09.028
-
Anion permeation in Ca(2+)-activated Cl(-) channelsThe Journal of General Physiology 116:825–844.https://doi.org/10.1085/jgp.116.6.825
-
The anoctamin family: TMEM16A and TMEM16B as calcium-activated chloride channelsExperimental Physiology 97:177–183.https://doi.org/10.1113/expphysiol.2011.058198
-
HOLE: a program for the analysis of the pore dimensions of ion channel structural modelsJournal of Molecular Graphics 14:354–360.https://doi.org/10.1016/S0263-7855(97)00009-X
-
Structure of the mechanosensitive OSCA channelsNature Structural & Molecular Biology 25:850–858.https://doi.org/10.1038/s41594-018-0117-6
-
An efficient one-step site-directed and site-saturation mutagenesis protocolNucleic Acids Research 32:e115.https://doi.org/10.1093/nar/gnh110
Article and author information
Author details
Funding
H2020 European Research Council (ERC no 339116 AnoBest)
- Raimund Dutzler
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank all members of the Dutzler lab for help at various stages of the project.
Copyright
© 2018, Lam and Dutzler
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,020
- views
-
- 430
- downloads
-
- 35
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.
-
- Structural Biology and Molecular Biophysics
Osmotic stress and chloride regulate the autophosphorylation and activity of the WNK1 and WNK3 kinase domains. The kinase domain of unphosphorylated WNK1 (uWNK1) is an asymmetric dimer possessing water molecules conserved in multiple uWNK1 crystal structures. Conserved waters are present in two networks, referred to here as conserved water networks 1 and 2 (CWN1 and CWN2). Here, we show that PEG400 applied to crystals of dimeric uWNK1 induces de-dimerization. Both the WNK1 the water networks and the chloride-binding site are disrupted by PEG400. CWN1 is surrounded by a cluster of pan-WNK-conserved charged residues. Here, we mutagenized these charges in WNK3, a highly active WNK isoform kinase domain, and WNK1, the isoform best studied crystallographically. Mutation of E314 in the Activation Loop of WNK3 (WNK3/E314Q and WNK3/E314A, and the homologous WNK1/E388A) enhanced the rate of autophosphorylation, and reduced chloride sensitivity. Other WNK3 mutants reduced the rate of autophosphorylation activity coupled with greater chloride sensitivity than wild-type. The water and chloride regulation thus appear linked. The lower activity of some mutants may reflect effects on catalysis. Crystallography showed that activating mutants introduced conformational changes in similar parts of the structure to those induced by PEG400. WNK activating mutations and crystallography support a role for CWN1 in WNK inhibition consistent with water functioning as an allosteric ligand.