TMEM16 and OSCA/TMEM63 proteins share a conserved potential to permeate ions and phospholipids

  1. Augustus J Lowry
  2. Pengfei Liang
  3. Mo Song
  4. Yuichun Wan
  5. Zhen-Ming Pei
  6. Huanghe Yang  Is a corresponding author
  7. Yang Zhang  Is a corresponding author
  1. Department of Biochemistry, Duke University School of Medicine, United States
  2. Institute of Molecular Physiology, Shenzhen Bay Laboratory, China
  3. Department of Biology, Duke University, United States
  4. Department of Neurobiology, Duke University School of Medicine, United States

eLife Assessment

This important study advances our understanding of the mechanisms controlling lipid flux and ion permeation in the TMEM16 and OSCA/TMEM63 family channels. The study provides compelling new evidence indicating that side chains along the TM4/6 interface play a key role in gating lipid and ion fluxes in these channels. The authors suggest that the transmembrane channel/scramblase family proteins may have originally functioned as scramblases but lost this capacity over evolution.

https://doi.org/10.7554/eLife.96957.3.sa0

Abstract

The calcium-activated TMEM16 proteins and the mechanosensitive/osmolarity-activated OSCA/TMEM63 proteins belong to the Transmembrane Channel/Scramblase (TCS) superfamily. Within the superfamily, OSCA/TMEM63 proteins, as well as TMEM16A and TMEM16B, are thought to function solely as ion channels. However, most TMEM16 members, including TMEM16F, maintain an additional function as scramblases, rapidly exchanging phospholipids between leaflets of the membrane. Although recent studies have advanced our understanding of TCS structure–function relationships, the molecular determinants of TCS ion and lipid permeation remain unclear. Here, we show that single mutations along the transmembrane helix (TM) 4/6 interface allow non-scrambling TCS members to permeate phospholipids. In particular, this study highlights the key role of TM 4 in controlling TCS ion and lipid permeation and offers novel insights into the evolution of the TCS superfamily, suggesting that, like TMEM16s, the OSCA/TMEM63 family maintains a conserved potential to permeate ions and phospholipids.

Introduction

Ion channels and phospholipid scramblases catalyze the passive flux of ions and phospholipids down their respective chemical gradients. Compared to ion channels, both our understanding of scramblases and the evolutionary origins of phospholipid scrambling are underdeveloped. Discoveries of two scramblase families—the TMEM16 calcium-activated phospholipid scramblases (CaPLSases) and XKR caspase-dependent phospholipid scramblases—have revealed roles for scramblases in processes such as blood coagulation (Suzuki et al., 2010; Yang et al., 2012; Fujii et al., 2015), angiogenesis (Shan et al., 2024), cell death signaling (Kunzelmann et al., 2019; Suzuki et al., 2013a), phagocytosis (Suzuki et al., 2013a), cell migration (Jacobsen et al., 2013; Chang et al., 2017), membrane repair (Griffin et al., 2016; Chandra et al., 2019), microparticle release (Fujii et al., 2015), cell–cell fusion (Whitlock et al., 2018; Zhang et al., 2020; Zhang et al., 2022), and viral infection (Zaitseva et al., 2017; Nanbo et al., 2018). Among identified scramblases, TMEM16 CaPLSases are the most extensively studied at the molecular level (Le et al., 2021). The TMEM16 family (Figure 1—figure supplement 1) was originally classified as calcium-activated chloride channels (CaCCs) based on the first-discovered members, TMEM16A and TMEM16B (Caputo et al., 2008; Yang et al., 2008; Schroeder et al., 2008). Most remaining members, however, exhibit CaPLSase activity, with TMEM16F representing the canonical CaPLSase (Suzuki et al., 2013b; Figure 1a). Uniquely, CaPLSases also function as non-selective ion channels along with their non-specific permeability to phospholipids (Yang et al., 2012; Grubb et al., 2013; Wong et al., 2013; Yu et al., 2015; Di Zanni et al., 2018; Le et al., 2019b). Substrate permeation is facilitated by calcium-induced conformational changes causing the clamshell-like separation of TMs 4 and 6 to facilitate ion and phospholipid permeation (Le et al., 2019a; Lam et al., 2021; Lam and Dutzler, 2021; Figure 1—figure supplement 1b–d). These conformational changes further catalyzed phospholipid permeation by thinning the membrane near the permeation pathway (Falzone et al., 2019; Falzone et al., 2022).

Figure 1 with 2 supplements see all
Lysine mutations along TM 4 enable TMEM16F channel and scramblase activities in the absence of calcium stimulation.

(a) Top: TMEM16F is a calcium-activated phospholipid scramblase. Bottom: TM 4 mutant locations mapped on TMEM16F CaPLSase structure with side chains shown as yellow sticks (PDB 6QPB). (b) Representative images of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged TMEM16F wildtype (WT), I521K, M522K, and T526K (center column). CF 594-conjugated Annexin V (AnV, right column) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left column). (c) Quantification of the percentage of cells with AnV labeling for TMEM16F WT (n = 5), I521K (n = 5), M522K (n = 5), and T526K (n = 5) transfected cells (Figure 1—source data 1). Values were derived from images of biological replicates. Statistical comparisons to T526K were done using unpaired t-tests with Welch’s correction (ns: p > 0.05, ****p < 0.0001). (d) Representative current recordings and (e) current–voltage (IV) relationships (Figure 1—source data 1) of whole-cell patches from TMEM16F KO HEK293T cells expressing eGFP-tagged TMEM16F WT (n = 6), I521K (n = 5), M522K (n = 5), and T526K (n = 7). Currents were elicited by the voltage protocol shown with the pipette solution containing 5 mM EGTA(Ethylene glycol tetraacetic acid). Dotted line denotes zero current. (f) Quantification of current at +160 mV (Figure 1—source data 1). Currents in (e) and (f) were normalized to cell capacitance. Statistical comparisons to T526K were done using unpaired t-tests with Welch’s correction (*p < 0.05, **p <0 .01). (g) Lysine mutations along TM 4 in TMEM16F enable spontaneous phospholipid permeation in the absence of calcium. All error bars represent standard error of the mean (SEM) calculated from replicate images or independent patches.

Figure 1—source data 1

Results from the lysine mutations on TM 4 of TMEM16F.

https://cdn.elifesciences.org/articles/96957/elife-96957-fig1-data1-v1.xlsx

The TMEM16 family together with the OSCA/TMEM63 family and TMC family (Keresztes et al., 2003; Hahn et al., 2009; Medrano-Soto et al., 2018; Ballesteros et al., 2018) collectively form the Transmembrane Channel/Scramblase (TCS) superfamily (Le et al., 2021). Despite the conserved 10-TM architectural core (Jojoa-Cruz et al., 2018; Liu et al., 2018; Zhang et al., 2018; Maity et al., 2019; Jeong et al., 2022), phospholipid permeability has not been experimentally demonstrated in mechanosensitive OSCA/TMEM63 or TMC proteins, raising the intriguing question of how TCS proteins discriminate between ion and ion/phospholipid substrates. Interestingly, studies have shown that the TMEM16A CaCC can be genetically modified to enable phospholipid permeability, either by substituting TMEM16F domains as small as 15 amino acids (Yu et al., 2015), or through single mutations at the hydrophobic gate (Le et al., 2019a). These findings suggest a modest energetic barrier for scramblase activity and lead us to hypothesize that evolutionary relatives of the TMEM16 family may maintain a conserved potential for both ion and phospholipid permeation.

To test this hypothesis, we introduced single mutations along TM 4 of TMEM16F, TMEM16A, OSCA1.2, and TMEM63A at sites near the corresponding hydrophobic gate of TMEM16F (Le et al., 2019a). Mutations along the TM 4/6 interface in TMEM16F and TMEM16A resulted in constitutive phospholipid scramblase activity. Strikingly, equivalent TM 4 mutations in OSCA1.2 and TMEM63A also converted these ion channels into phospholipid scramblases, which were either constitutive or activated by osmotic stimulation. Individual mutations also resulted in gain-of-function (GOF) ion channel activity, suggesting that these mutations may disrupt a conserved activation gate in these evolutionarily related proteins. Together, our findings (1) advance the mechanistic understanding of gating and substrate permeation in the TMEM16 and OSCA/TMEM63 families, (2) underscore a key design principle for phospholipid scramblases, and (3) support our hypothesis that the TCS superfamily maintains a conserved potential for both ion and phospholipid permeation.

Results

TM 4 mutations result in a constitutively active TMEM16F scramblase

Our previous study identified the activation gate of TMEM16F, consisting of hydrophobic residues from TMs 4, 5, and 6, which collectively govern ion and phospholipid permeation (Le et al., 2019a). Mutations to the activation gate result in GOF CaPLSases, some of which are constitutively activated without requiring calcium (Le et al., 2019a). The hydrophobic gate represents the most constricted point along the TM 4/6 interface in TMEM16F and mutations along this interface likely alter substrate permeation and gating. Structural and computational studies on the fungal NhTMEM16 (Kalienkova et al., 2019) and AfTMEM16 (Falzone et al., 2019) orthologs, as well as human TMEM16K (Bushell et al., 2019) further support a clamshell-like gating model (Le et al., 2019a) whereby scrambling activity is promoted by separation of the TM 4/6 interface (Le et al., 2021; Figure 1—figure supplement 1b–d). The conformational transition involves the N-terminal half of TM 4, which bends away from TM 6, and the C-terminal half of TM 6, which collapses onto the calcium-binding sites formed together with TMs 7 and 8 (Le et al., 2021). Given that the N-termini of TM 4 are largely hydrophobic among TMEM16 CaPLSases (Figure 1—figure supplement 1f) and mutations at F518 lead to constitutive TMEM16F activity (Le et al., 2019a), we hypothesized that introducing a charged side chain along the TM 4/6 interface might disrupt the hydrophobic gate and result in TMEM16 scramblases with GOF scramblase activity. To test this hypothesis, we overexpressed eGFP-tagged mutant constructs in TMEM16F knockout (KO) HEK293T cells and used an established scramblase assay (Le et al., 2019b; Le et al., 2019a; Le et al., 2020) that detects phosphatidylserine (PS) exposure using fluorophore-tagged Annexin V (AnV) as a reporter (Figure 1b). In the absence of calcium stimulation, PS is predominantly in the inner leaflet of the plasma membrane of TMEM16F wildtype (WT) expressing cells and is therefore not labeled by extracellular AnV. Similar to TMEM16F F518K (Le et al., 2019a), overexpressing the single lysine mutations TMEM16F I521K and M522K led to spontaneous, global exposure of PS on the plasma membrane without requiring calcium stimulation (Figure 1b, c). In contrast, T526K, which is one helical turn below I521 and M522 (Figure 1a), caused minimal spontaneous PS exposure (Figure 1b, c). We further tested whether I521K, M522K, and T526K have GOF ion channel activity using whole-cell patch clamp. In the absence of calcium, I521K and M522K, but not WT TMEM16F or T526K, showed robust depolarization-activated outward rectifying currents (Figure 1d–f). Our imaging and patch clamp experiments thus demonstrate that I521K and M522K are functionally expressed on the plasma membrane and their GOF phospholipid scramblase activity results in spontaneous, global PS exposure on the cell surface. Conversely, despite apparent plasma membrane localization, T526K exhibited strongly attenuated GOF ion channel and scramblase activities. The TM 4/6 interaction is notably less prominent near T526 (Figure 1a), two helical turns below the hydrophobic gate residue, F518 (Le et al., 2019b; Figure 1a) and may represent a lower limit for this charge-induced effect. On the other hand, I521 and M522 are approximately one helical turn below the hydrophobic gate (Figure 1a) supporting the idea that charged mutations along the TM 4/6 interface promote gate opening and substrate permeation (Figure 1g).

Our previous study also showed that glutamate, but not alanine, substitution at F518 led to spontaneous scramblase activity in TMEM16F (Le et al., 2019a). To assess whether alternative side chains cause spontaneous scrambling below the hydrophobic gate in TMEM16F, we tested I521A and I521E. I521E led to spontaneous, global PS exposure in all transfected cells, whereas I521A failed to cause spontaneous PS exposure (Figure 1—figure supplement 2). This result is consistent with our previous finding, providing additional support that charged TM 4 mutations disrupt the TM 4/6 interface.

I611K on TM 6 results in a constitutively active TMEM16F scramblase

We next tested whether a TM 6 mutation along the TM 4/6 interface could promote scrambling. Directly adjacent to I521 and M522 are a pair of glycine residues on TM 6 that are proposed to function as a hinge during calcium-dependent activation (Paulino et al., 2017). We therefore chose to introduce a lysine mutation one helical turn above this site at I611, a residue adjacent to the pore-facing inner gate residue I612 (Le et al., 2019b; Figure 2a). Comparable to I521K (Figure 1b, c), I521E (Figure 1—figure supplement 2), and M522K (Figure 1b, c), I611K resulted in spontaneous, global PS exposure (Figure 2b, c). Likewise, I611K-expressing cells exhibited robust, depolarization-activated outward rectifying currents in the absence of calcium (Figure 2d–f). Together with our previous report that I612K is a GOF phospholipid scramblase and ion channel (Le et al., 2019a), I611K further supports that charged mutations in TM 6 near the activation gate can also promote substrate permeation in TMEM16F (Figure 2g).

I611K on TM 6 enables TMEM16F channel and scramblase activities in the absence of calcium stimulation.

(a) I611 highlighted on the TMEM16F CaPLSase structure with side chains shown as yellow sticks (PDB 6QPB). (b) Representative images of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged TMEM16F wildtype (WT) and I611K (center column). CF 594-conjugated Annexin V (AnV, right column) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left column). (c) Quantification of the percentage of cells with AnV labeling for TMEM16F WT (n = 5) and I611K (n = 5) transfected cells (Figure 2—source data 1). Values were derived from images of biological replicates. Statistical comparison was done using an unpaired t-test with Welch’s correction (****p < 0.0001). (d) Representative current recordings and (e) current–voltage (IV) relationships (Figure 2—source data 1) of whole-cell patches from TMEM16F KO HEK293T cells expressing eGFP-tagged TMEM16F WT (n = 6) and I611K (n = 5). Currents were elicited by the voltage protocol shown with the pipette solution containing 5 mM EGTA. Dotted line denotes zero current. (f) Quantification of current at +160 mV (Figure 2—source data 1). Currents in (e) and (f) were normalized to cell capacitance. Statistical comparison was done using an unpaired t-test with Welch’s correction (***p < 0.001). (g) A lysine mutation on TM 6 in TMEM16F enables spontaneous phospholipid permeation in the absence of calcium. All error bars represent standard error of the mean (SEM) calculated from replicate images or independent patches.

TM 4 lysine mutations convert TMEM16A into a constitutively active scramblase

TMEM16A is a CaCC (Caputo et al., 2008; Yang et al., 2008; Schroeder et al., 2008) without scramblase activity (Yu et al., 2015; Le et al., 2019a; Figure 3a–c). We previously reported that a single TM 4 lysine mutation, L543K, at the hydrophobic gate of TMEM16A (Figure 3a) allows the CaCC to permeate phospholipids spontaneously, analogous to the F518K mutation in TMEM16F (Le et al., 2019a). Thus, we reasoned that lysine mutations equivalent to TMEM16F I521K and M522K might also result in spontaneous phospholipid permeability. To test this hypothesis, we overexpressed eGFP-tagged TMEM16A I546K and I547K mutants in TMEM16F KO HEK293T cells and assessed their ability to expose PS on the plasma membrane with confocal microscopy. Similar to TMEM16A L543K (Le et al., 2019a), I546K- and I547K-expressing cells exhibit spontaneous, global PS exposure (Figure 3b, c), analogous to the equivalent TMEM16F mutations I521K and M522K (Figure 1—figure supplement 1f; Figure 1b, c). Similarly, whole-cell patch clamp revealed GOF ion channel activity at depolarizing potentials, even in the absence of calcium stimulation (Figure 3d–f). In contrast, E551K, equivalent to T526K in TMEM16F, failed to exhibit PS exposure, but showed modest ion channel activity at depolarizing potentials in the absence of calcium, perhaps reflecting the gate is open sufficiently for ion but not phospholipid permeation. In total, these results suggest that TM 4 mutations in this region destabilize the interface to endow the TMEM16A CaCC with GOF ion channel activity, and, in the case of I546K and I547K, constitutive phospholipid permeability (Figure 3g).

Figure 3 with 1 supplement see all
Lysine mutations along TM 4 enable TMEM16A channel and scramblase activities in the absence of calcium stimulation.

(a) Top: TMEM16A is a calcium-activated chloride channel (CaCC). Bottom: TM 4 mutant locations mapped on TMEM16A CaCC structure with side chains shown as yellow sticks (PDB 5OYG). (b) Representative images of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged TMEM16A wildtype (WT), I546K, I547K, and E551K (center column). CF 594-conjugated Annexin V (AnV, right column) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left column). (c) Quantification of the percentage of cells with AnV labeling for TMEM16A WT (n = 4), I546K (n = 4), I547K (n = 4), and E551K (n = 5) transfected cells (Figure 3—source data 1). Values were derived from images of biological replicates. Statistical comparisons were done using unpaired t-tests with Welch’s correction (*p < 0.05, **p < 0.01). (d) Representative whole-cell current recordings and (e) current–voltage (IV) relationships (Figure 3—source data 1) of whole-cell patches from TMEM16F KO HEK293T cells expressing eGFP-tagged TMEM16A WT (n = 14), I546K (n = 5), I547K (n = 5), and E551K (n = 6). Currents were elicited by the voltage protocol shown with the pipette containing 5 mM EGTA. Dotted line denotes zero current. (f) Quantification of current at +160 mV (Figure 3—source data 1). Currents in (e) and (f) were normalized to cell capacitance with the mean ± SEM calculated from independent patches. Statistical comparisons were done using unpaired t-tests with Welch’s correction (*p < 0.05, ****p < 0.0001). (g) Lysine mutations along TM 4 in TMEM16A enable spontaneous phospholipid permeation in the absence of calcium. All error bars represent standard error of the mean (SEM) calculated from replicate images or independent patches.

Figure 3—source data 1

Results from the lysine mutations on TM 4 of TMEM16A.

https://cdn.elifesciences.org/articles/96957/elife-96957-fig3-data1-v1.xlsx

Ani9, a selective TMEM16A inhibitor, attenuates I546K ion and phospholipid permeation

We then tested whether Ani9, a selective, extracellular TMEM16A inhibitor (Seo et al., 2016), could attenuate phospholipid and ion permeation through I546K. Similar to our previous observation of TMEM16A L543K (Le et al., 2019a), chronic incubation of Ani9 at 10 μM also prevented spontaneous PS exposure for I546K (Figure 3—figure supplement 1d). Interestingly, 10 μM Ani9 reversibly inhibited voltage-elicited TMEM16A I546K currents by over 50% (Figure 3—figure supplement 1a–c). These results demonstrate that although the TM 4 GOF mutations alter TMEM16A substrate permeation, they retain their sensitivity to Ani9 inhibition (Figure 3—figure supplement 1e). This further supports that the TM 4 GOF mutations of TMEM16A mediate the spontaneous, global PS exposure.

L438K on TM 4 converts OSCA1.2 into a constitutively active scramblase

We next applied our approach to a TMEM16 relative from the OSCA/TMEM63 family (Figure 1—figure supplement 1a), which was first discovered in plants as a family of mechanosensitive and osmolarity-activated cation non-selective ion channels (Yuan et al., 2014; Figure 4a). We hypothesized that analogous lysine mutations on TM 4 along the TM 4/6 interface of OSCA/TMEM63 proteins would result in GOF channels that may also permeate phospholipids. Within the family, we chose OSCA1.2 from Arabidopsis thaliana due to previous structural (Jojoa-Cruz et al., 2018) and biophysical (Murthy et al., 2018) characterization demonstrating that the channel is activated directly by membrane tension. Overexpressing eGFP-tagged OSCA1.2 WT in TMEM16F KO HEK293T cells did not induce PS exposure, demonstrating that OSCA1.2 WT lacks spontaneous phospholipid scramblase activity (Figure 4b). We then introduced TM 4 lysine mutations at analogous sites, as identified by structural (Figure 4a) and sequence alignment (Figure 1—figure supplement 1f). Strikingly, a single-point mutation, L438K, causes cells overexpressing the OSCA1.2 mutant to exhibit spontaneous and global PS exposure (Figure 4b, c). This mirrors our results with both TMEM16F (Figure 1) and TMEM16A (Figure 3), suggesting that the L438K mutation allows the OSCA1.2 channel to permeate phospholipids. Next, we used inside-out patch clamp to examine if L438K also enhances OSCA1.2 channel activity. We found that the mutant significantly left-shifts the conductance–voltage (G–V) relationship (Figure 4e, f) and accelerates channel activation kinetics (Figure 4g) compared to WT. Under −50 mmHg of pressure, L438K has a half-maximal voltage (V0.5) of 66.7 ± 3.7 mV, while the WT V0.5 is nearly 108.7 ± 5.6 mV (Figure 4f), underscoring that this mutation also promotes channel gating. Together, these experiments show that, like TMEM16A, a single lysine mutation near the putative gate allows the mechanosensitive and osmolarity-activated OSCA1.2 channel to become spontaneously permeable to phospholipids (Figure 4h).

Figure 4 with 3 supplements see all
Lysine mutations along TM 4 enable OSCA1.2 channel and scramblase activities.

(a) Top: OSCA1.2 is a cation non-selective ion channel gated by membrane tension. Bottom: TM 4 mutant locations mapped onto the TM 4/6 interface of OSCA1.2 (PDB 6MGV) with key residues shown as yellow sticks. (b) Representative images of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged OSCA1.2 wildtype (WT), L438K, or A439K mutants (center column). CF 594-conjugated Annexin V (AnV) (right column) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left column). Asterisk highlights a PS-positive cell for the A439K mutant. (c) Quantification of the percentage of cells with AnV labeling for OSCA1.2 WT (n = 4), L438K (n = 7), and A439K-transfected cells (n = 6) (Figure 4—source data 1). Statistical comparisons were conducted with unpaired t-tests with Welch’s correction (*p < 0.05, ****p < 0.0001). (d) Representative current recordings and (e) normalized conductance–voltage (GV) relationships of inside-out patches from TMEM16F KO HEK293T cells expressing eGFP-tagged OSCA1.2 WT (n = 8), L438K (n = 8), and A439K (n = 6) (Figure 4—source data 1). Currents were elicited by the voltage protocol shown next to the listed pressures. Dotted lines denote zero current. (f) Quantification of half-maximal voltage at −50 mmHg for WT (109 mV), L438K (67 mV), and A439K (63 mV) (Figure 4—source data 1). Statistical comparison was conducted with unpaired t-tests with Welch’s correction (***p < 0.001, ****p < 0.0001). (g) Quantification of activation τon at −50 mmHg and 160 mV for WT (41 ms), L438K (13 ms), and A439K (16 ms) (Figure 4—source data 1). Statistical comparison was conducted with unpaired t-tests with Welch’s correction (***p < 0.001, ****p < 0.0001). (h) A lysine mutation along TM 4 converts the OSCA1.2 channel into a phospholipid scramblase with spontaneous phospholipid permeability. All error bars represent standard error of the mean (SEM) calculated from replicate images or independent patches.

Figure 4—source data 1

Results from the lysine mutations on TM 4 of OSCA1.2.

https://cdn.elifesciences.org/articles/96957/elife-96957-fig4-data1-v1.xlsx

To help visualize how L438K alters the TM 4/6 interface of OSCA1.2, we created a homology model of the mutant using the Swiss-PDB server and inserted it into a phospholipid membrane with the CHARMM-GUI webserver (Jo et al., 2008). After equilibration, we ran 900 ns atomistic molecular dynamics (MD) simulations for both WT and L438K in GROMACS, using a combination of conventional and metadynamics simulations. Remarkably, the L438K trajectory showed considerably increased hydration of the pore region compared to WT (Figure 4—figure supplement 1b). This enhanced hydration is attributed to a modest expansion of the TM 4/6 interface throughout the simulation (Figure 4—figure supplement 1a). Together with our functional results, these MD simulations reinforce our hypothesis that disruption of the TM 4/6 interface of OSCA1.2 promotes GOF ion and phospholipid permeation.

One helical turn above OSCA1.2 L438, at the corresponding site to the activation gate residue F518 in TMEM16F, is an endogenous lysine at position 435 (Figure 4—figure supplement 2a). This positioning places two positively charged residues in close proximity within TM 4 of the OSCA1.2 L438K mutant. Previous studies on model transmembrane helical peptides have demonstrated that hydrophilic side chains within the hydrophobic core of the membrane promote trans-bilayer phospholipid transport and this effect is further enhanced when these hydrophilic residues are stacked within the helix (Nakao and Nakano, 2022). To explore the role of the endogenous K435 in facilitating L438K-mediated spontaneous phospholipid permeability, we assessed PS exposure in the K435L and K435L/L438K mutants. Remarkably, the double mutant K435L/L438K, which shifts the lysine to position 438, resulted in a modest level of spontaneous PS exposure (31%, Figure 4—figure supplement 2b, c). This level is intermediate between the WT (K435) with no PS exposure and the L438K mutant (K435/K438) with 63% spontaneous PS exposure (Figure 4c; Figure 4—figure supplement 2b, c). In contrast, K435L (L435/L438) did not induce significant spontaneous PS exposure (3%, Figure 4—figure supplement 2b, c). These findings, consistent with previous studies on model transmembrane peptides (Nakao and Nakano, 2022), highlight the role of multiple charged resides in TM 4 in promoting phospholipid scrambling. Moreover, our results further support the hypothesis that the instability of the TM 4/6 interface in OSCA1.2 is crucial for controlling phospholipid and ion permeation.

A439K on TM4 converts OSCA1.2 into an osmolarity-sensing scramblase

Interestingly, mutating the neighboring amino acid (Figure 4a), A439, to lysine resulted in minimal spontaneous PS exposure (10%, Figure 4b, c). We thus reasoned that A439K scramblase activity may require additional stimulation to trigger scramblase activity. Given that WT OSCA1.2 ion channel activity can be induced by hypotonic treatment (Figure 4—figure supplement 3), we acutely treated WT- and A439K-expressing cells with a hypotonic solution (120 mOsm/kg) and assessed scramblase activity using time-lapse imaging. Indeed, PS exposure was robustly induced for the A439K mutant (Figure 5b, c) but not WT (Figure 5a, c) in response to hypotonic stimulation. Inside-out patch clamp further demonstrated that A439K enhances OSCA1.2 ion channel activity as evidenced by the accelerated activation kinetics (Figure 4g) and left-shifted G–V relationship (Figure 4f). Our experiments thus indicate that A439K modestly disrupts OSCA1.2 gating and converts the osmolarity-activated ion channel (Figure 4—figure supplement 3) into an osmolarity-sensing phospholipid scramblase (Figure 5d).

OSCA1.2 A439K is an osmolarity-activated scramblase.

Representative images of hypotonic osmolarity stimulation of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged OSCA1.2 (a) wildtype (WT) or (b) the A439K mutant (center columns). CF 594-conjugated Annexin V (AnV) (right columns) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left columns). Each row of representative images corresponds to the indicated time after hypo-osmotic stimulation. (c) Quantification of AnV intensity for OSCA1.2 WT (n = 5) and A439K (n = 5) after hypo-osmotic stimulation (Figure 5—source data 1). Statistical comparison was conducted with an unpaired t-test with Welch’s correction (**p < 0.01). (d) The A439K mutation converts OSCA1.2 to an osmolarity-activated phospholipid scramblase. Error bars represent standard error of the mean (SEM).

Figure 5—source data 1

Results from osmolarity activation of OSCA1.2 A439K scramblase.

https://cdn.elifesciences.org/articles/96957/elife-96957-fig5-data1-v1.xlsx

TM 4 lysine mutations convert TMEM63A into a constitutively active scramblase

Finally, we turned our attention to TMEM63A to further investigate the evolutionary conservation of our observation that TM 4 lysine mutations convert TMEM16 and OSCA members into scramblases. TMEM63s represent the mammalian members of the OSCA/TMEM63 family with three members present in humans (TMEM63A–C). Recent structural and functional characterizations indicate that TMEM63s function as mechanosensitive ion channels gated by high-threshold membrane tension and, in notable contrast to all other structurally resolved TCS superfamily members, they likely function as monomers (Zheng et al., 2023; Qin et al., 2023; Wu et al., 2024). Given their structural homology to TMEM16s and OSCA1.2 (Zheng et al., 2023), we hypothesized that TM 4 mutations in TMEM63A would also result in GOF activity. We again identified residues near the putative gate by structural (Figure 6a) and sequence alignment (Figure 1—figure supplement 1f), selecting W472 (equivalent to F518 in TMEM16F, L543 in TMEM16A, and K435 in OSCA1.2), S475 (equivalent to I521 in TMEM16F, I546 in TMEM16A, and L438 in OSCA1.2), and A476 (equivalent to M522 in TMEM16F, I547 in TMEM16A, and A439 in OSCA1.2). Indeed, overexpressing eGFP-tagged mouse TMEM63A with single lysine mutations at either W472 or S475 led to spontaneous PS exposure (Figure 6b, c), though the AnV staining revealed punctate rather than global patterns of PS exposure induced by their OSCA1.2 counterparts. Notably, A476K failed to show obvious PS exposure (Figure 6b, c), even after 120 mOsm/kg hypotonic treatment (Figure 6—figure supplement 1b, c), perhaps reflecting differences in mechanical pressure threshold between OSCA1.2 and TMEM63A (Zheng et al., 2023). To confirm membrane localization and further probe mutant effects, we exploited TMEM63A ion channel function using the cell-attached patch clamp configuration. As TMEM63A exhibits voltage-dependent activity under high pressure, we compared mutant and WT I–V relationships at −80 mmHg. W472K, S475K, and A476K caused marked reductions in V0.5 from 122.3 ± 3.5 mV for WT to 95.9 ± 6.3, 92.1 ± 8.6, and 96.8 ± 4.2 mV, respectively (Figure 6d–f). By comparison, mock-transfected controls failed to elicit current (Figure 6d, e). Together, these results are consistent with our observations in TMEM16F, TMEM16A, and OSCA1.2, indicating that single lysine mutations along TM 4 of TMEM63A can facilitate both ion and phospholipid permeation (Figure 6g).

Figure 6 with 1 supplement see all
Lysine mutations along TM 4 enable TMEM63A channel and scramblase activities.

(a) Top: TMEM63A is an ion channel gated by high-threshold membrane tension. Bottom: the TM 4/6 interface of HsTMEM63A (PDB 8GRS) with key residues shown as yellow sticks using amino acid numbering corresponding to the mouse ortholog. (b) Representative images of TMEM16F knockout (KO) HEK293T cells expressing eGFP-tagged TMEM63A wildtype (WT), W472K, S475K, and A476K (center column). CF 594-conjugated Annexin V (AnV) (right column) labeled phosphatidylserine (PS) exposing cells. BF denotes bright field images (left column). (c) Quantification of the percentage of cells with AnV labeling for TMEM63A WT (n = 4), W472K (n = 4), S475K (n = 4), and A476K-transfected cells (n = 4) (Figure 6—source data 1). Statistical comparisons were conducted with unpaired t-tests with Welch’s correction (**p < 0.01, ****p < 0.0001). (d) Representative current recordings and (e) normalized conductance–voltage (IV) relationships of cell-attached patches from TMEM16F KO HEK293T cells expressing either eGFP mock-transfected (n = 7) or eGFP-tagged TMEM63A WT (n = 7), W472K (n = 7), S475K (n = 5), or A476K (n = 7) (Figure 6—source data 1). Currents represent the subtraction of voltage alone from currents elicited by the voltage and pressure protocols shown. Dotted line denotes zero current. Note that the mock control was normalized to the mean maximal current elicited from WT-transfected cells. (f) Quantification of half-maximal voltage at −80 mmHg for WT (122 mV), W472K (96 mV), S475K (92 mV), and A476K (97 mV) (Figure 6—source data 1). Statistical comparison was conducted with unpaired t-tests with Welch’s correction (*p < 0.05, **p < 0.01, ***p < 0.001). (g) Lysine mutations along TM 4 in TMEM63A enable spontaneous phospholipid permeability. All error bars represent standard error of the mean (SEM) calculated from replicate images or independent patches.

Figure 6—source data 1

Results from the lysine mutations on TM 4 of TMEM63A.

https://cdn.elifesciences.org/articles/96957/elife-96957-fig6-data1-v1.xlsx

Discussion

Mechanistically, our study improves models of TMEM16 substrate permeation and gating. We identified multiple mutations in TMEM16F and TMEM16A that promote phospholipid permeation and cause commensurate changes in ion channel activities. Our findings complement functional characterizations of TMEM16F gating mutants where gate destabilization is inversely correlated with side chain hydropathy (Le et al., 2019a; Arndt et al., 2022). Lysine is above only arginine on the hydropathy index and thus likely explains why it readily destabilizes the gate. For instance, F518K exhibits spontaneous PS exposure, even when the calcium-binding site is destroyed (Le et al., 2019a), whereas F518H does not exhibit spontaneous PS exposure (Arndt et al., 2022). Interestingly, the recent TMEM16F F518H structure (Figure 1—figure supplement 1e) shows local membrane thinning due in part to unexpected conformational changes in TM 3 (Arndt et al., 2022). Future structural studies are needed to assess whether TMEM16A or OSCA/TMEM63 mutant scramblases also promote membrane thinning and/or conformational rearrangements in TM 3. More broadly, our results highlight an increasingly appreciated design principle of scramblases where polar and charged residues often line a membrane-spanning groove. This observation has been noted for TMEM16 (Bushell et al., 2019; Brunner et al., 2014; Alvadia et al., 2019), Xkr (Straub et al., 2021; Sakuragi et al., 2021), and opsin (Morra et al., 2018) scramblases and should be a key criterion for identifying and characterizing new scramblases.

Our findings also advance our understanding of evolutionary relatives of the TCS superfamily (Le et al., 2021) and help uncover their roles in human diseases. Although OSCA/TMEM63 proteins are not known to scramble phospholipids and we did not detect obvious hypotonicity-induced PS exposure in cells overexpressing WT OSCA1.2 (Figure 5a) or TMEM63A (Figure 1—figure supplement 1a), we show that single mutations in TM 4 of OSCA1.2 and TMEM63A convert these osmolarity-activated and/or mechanosensitive ion channels into phospholipid scramblases, similar to our findings with TMEM16A mutants. We thus speculate that the conserved structural architecture in the transmembrane region endows TCS proteins with a potential to scramble phospholipids, though this capability may have been lost by some members during evolution. This hypothesis is especially intriguing given the recent OSCA open state structures, which detail dramatic conformational rearrangements of TMs 3–6 leading to a phospholipid-lined (or proteolipidic) pore near TMs 4 and 6 (Han et al., 2024). It will be interesting to test whether equivalent mutations can convert transmembrane channel-like (TMC) proteins—the third TCS relative of TMEM16 and OSCA/TMEM63 (Le et al., 2021)—into phospholipid scramblases. TMC1 is best known for its role in auditory sensation, and thus far has mostly been characterized in vivo due to expression difficulties in heterologous systems (Pan et al., 2018). However, recent in vivo characterization of mouse TMC1 M412K, known as the Beethoven mutation, provided an important insight (Ballesteros and Swartz, 2022). The deafness-associated mutation is located in TM 4 (Figure 1—figure supplement 1f; equivalent of M522 in TMEM16F, I547 in TMEM16A, A439 in OSCA1.2, and A476 in TMEM63A) and results in constitutive PS exposure when expressed in the hair cell membranes of both heterozygous and homozygous mice (Ballesteros and Swartz, 2022). This raises the intriguing possibility that the Beethoven mutation may enable TMC1 to spontaneously permeate phospholipids, leading to a loss of membrane homeostasis and ultimately, auditory sensation (Sakuragi and Nagata, 2023). The possibility of converting TMC proteins into phospholipid scramblases should be thoroughly investigated. Additionally, disease-associated mutations in TMEM63 proteins are present along the TM 4/6 interface, such as TMEM63B T481N (Zheng et al., 2023). We speculate that introducing more hydropathic side chains along this interface may lead to spontaneous ion and or phospholipid permeability, perhaps contributing to underlying pathophysiology.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Mus musculus)MmTMEM16FNCBINP_780553.2Contains a three amino
acid N-terminal truncation (MQM)
Gene (M. musculus)MmTMEM16ANCBINP_001229278
Gene (Arabidopsis thaliana)AtOSCA1.2GeneBankAIU34614.1
Gene (M. musculus)MmTMEM63ANCBINP_001404481.1
Commercial assay or kitIn-Fusion Snap AssemblyTakara638947
Recombinant DNA reagentpeGFP-N1Addgene: Vector Database—pEGFP-N1
Recombinant DNA reagentpeGFP-N1_MmTMEM16FLe et al., 2019a
Recombinant DNA reagentpeGFP-N1_MmTMEM16ALe et al., 2019c
Recombinant DNA reagentpeGFP-N1_AtOSCA1_2This paperSubcloned with In-fusion (IF) primers
Recombinant DNA reagentpeGFP-N1_MmTMEM63AThis paperSubcloned with In-fusion (IF) primers
Sequence-based reagenteGFP-N1_678–679_FThis paperPCR primersATGGTGAGCAAGGGCGAGG
Sequence-based reagenteGFP-N1_678–679_RThis paperPCR primersGGTGGCGACCGGTGGATCC
Sequence-based reagentIFpeGFP-N1_OSCA1_2-eGFP_FThis paperPCR primersCCACCGGTCGCCACCATGGC
GACACTTCAGGAT
Sequence-based reagentIFpeGFP-N1_OSCA1_2-eGFP_RThis paperPCR primersGCCCTTGCTCACCATGACTAG
TTTACCACTAAAGGGC
Sequence-based reagentIFpeGFP-N1_MmTMEM63A-eGFP_FThis paperPCR primersCCACCGGTCGCCACCATGACC
AGCTCCCCGTTCC
Sequence-based reagentIFpeGFP-N1_MmTMEM63A-eGFP_RThis paperPCR primersGCCCTTGCTCACCATGGATTC
CTGGTAGGCATAAGC
Sequence-based reagentMmTMEM16F_I521K_FThis paperPCR primersTTCATCATCAAGATG
ATCCTCAACACG
Sequence-based reagentMmTMEM16F_I521K_RThis paperPCR primersGAGGATCATCTTGA
TGATGAAGCTGATG
Sequence-based reagentMmTMEM16F_I521A_FThis paperPCR primersCTTCATCATCGCC
ATGATCCTCAACAC
Sequence-based reagentMmTMEM16F_I521A_RThis paperPCR primersGAGGATCATGGCG
ATGATGAAGCTG
Sequence-based reagentMmTMEM16F_I521E_FThis paperPCR primersTCATCATCGAGATGATCCTCAACACG
Sequence-based reagentMmTMEM16F_I521E_RThis paperPCR primersTGAGGATCATCTCGATGATGAAGCTG
Sequence-based reagentMmTMEM16F_M522K_FThis paperPCR primersATCATCATCAAGATCCTCAACACG
Sequence-based reagentMmTMEM16F_M522K_RThis paperPCR primersTGAGGATCTTGATGATGATGAAGC
Sequence-based reagentMmTMEM16F_T526K_FThis paperPCR primersTCCTCAACAAGATCTACGAGAAGGTGG
Sequence-based reagentMmTMEM16F_T526K_RThis paperPCR primersCTCGTAGATCTTGTTGAGGATCATGATG
Sequence-based reagentMmTMEM16F_I611K_FThis paperPCR primersGCTGACGAAGATCATGGGGGGA
Sequence-based reagentMmTMEM16F_I611K_RThis paperPCR primersCCCCCATGATCTTCGTCAGCTGTG
Sequence-based reagentMmTMEM16A_I546K_FThis paperPCR primersCGTGGTCAAGATTCTGCTGGATGAAG
Sequence-based reagentMmTMEM16A_I546K_RThis paperPCR primersCCAGCAGAATCTTGACCACGAGGT
Sequence-based reagentMmTMEM16A_I547K_FThis paperPCR primersTGGTCATCAAGCTGCTGGATGAAG
Sequence-based reagentMmTMEM16A_I547K_RThis paperPCR primersTCCAGCAGCTTGATGACCACGAG
Sequence-based reagentMmTMEM16A_E551K_FThis paperPCR primersTCTGCTGGATAAGGTTTACGGCTGC
Sequence-based reagentMmTMEM16A_E551K_RThis paperPCR primersGCCGTAAACCTTATCCAGCAGAATGATG
Sequence-based reagentAtOSCA1_2_K435L_FThis paperPCR primersTTGCACTGCTGCTTTTCCTCGC
Sequence-based reagentAtOSCA1_2_K435L_RThis paperPCR primersGAGGAAAAGCAGCAGTGCAATACCC
Sequence-based reagentAtOSCA1_2_K435L_L438K_FThis paperPCR primersATTGCACTGCTGCTTTTCAAGGCC
Sequence-based reagentAtOSCA1_2_K435L_L438K_RThis paperPCR primersCTTGAAAAGCAGCAGTGCAATACCCG
Sequence-based reagentAtOSCA1_2_L438K_FThis paperPCR primersAAGCTTTTCAAGGCC
TTTCTGCCATC
Sequence-based reagentAtOSCA1_2_L438K_RThis paperPCR primersCAGAAAGGCCTTGA
AAAGCTTCAGTGC
Sequence-based reagentAtOSCA1_2_A439K_FThis paperPCR primersGCTTTTCCTCAAGT
TTCTGCCATCC
Sequence-based reagentAtOSCA1_2_A439K_RThis paperPCR primersGGCAGAAACTTGAGGAAAAGCTTCAG
Sequence-based reagentMmTMEM63A_W472K_FThis paperPCR primersCTCCTGCTGAAGTCCTTCTCTGCG
Sequence-based reagentMmTMEM63A_W472K_RThis paperPCR primersAGAGAAGGACTTCAGCAGGAGTGTGG
Sequence-based reagentMmTMEM63A_S475K_FThis paperPCR primersGGTCCTTCAAGGCGCTGCTTCCG
Sequence-based reagentMmTMEM63A_S475K_RThis paperPCR primersAAGCAGCGCCTTGAAGGACCACAGC
Sequence-based reagentMmTMEM63A_A476K_FThis paperPCR primersTCCTTCTCTAAGCT
GCTTCCGTCCA
Sequence-based reagentMmTMEM63A_A476K_RThis paperPCR primersCGGAAGCAGCTT
AGAGAAGGACCAC
Cell line (Homo sapiens)TMEM16F KO HEK293TLe et al., 2019b
Peptide, recombinant proteinAnnexinV-594Biotium#29011
Chemical compound, drugAni9Sigma-AldrichSML1813
Software, algorithmScrambling assayhttps://github.com/yanghuanghe/scrambling_activity, copy archived at Huanghe, 2020

Cloning and mutagenesis

Request a detailed protocol

All constructs used a peGFP-N1 vector backbone. GFP mock controls used the empty peGFP-N1 vector. WT sequence and mutation numbers correspond to NCBI: NP_780553.2 (Mus musculus TMEM16F) with a three amino acid (MQM) N-terminal truncation, NCBI: NP_001229278 (M. musculus TMEM16A), GenBank: AIU34614.1 (A. thaliana OSCA1.2), and NCBI: NP_001404481.1 (M. musculus TMEM63A). AtOSCA1.2 and MmTMEM63A cDNAs were subcloned using In-Fusion Snap Assembly (Takara Bio, Kusatsu, Japan, Cat. #638947). Point mutants were generated by PCR site-directed mutagenesis with primers from IDT DNA Technologies (Coralville, IA). Sequences were confirmed by Sanger sequencing (Azenta, Burlington, MA).

Bioinformatics

Request a detailed protocol

The following sequences were obtained from UniProt and aligned using Clustal Omega: Q8BHY3-2 (MmTMEM16A), Q9NQ90 (HsTMEM16B), Q9BYT9 (HsTMEM16C), Q32M45 (HsTMEM16D), Q75V66 (HsTMEM16E), Q6P9J9 (MmTMEM16F), Q6IWH7 (HsTMEM16G), Q9HCE9 (HsTMEM16H), A1A5B4 (HsTMEM16J), Q9NW15 (HsTMEM16K), C7Z7K1 (NhTMEM16), Q4WA18 (AfTMEM16), Q9XEA1 (AtOSCA1.1), Q5XEZ5 (AtOSC1.2), B5TYT3 (AtOSCA1.3), A0A097NUQ0 (AtOSCA1.4), A0A097NUS0 (AtOSCA1.5), A0A097NUP1 (AtOSCA1.6), A0A097NUP8 (AtOSCA1.7), A0A097NUQ2 (AtOSCA1.8), A0A097NUQ5 (AtOSCA2.1), A0A097NUS5 (AtOSCA2.2), A0A097NUP6 (AtOSCA2.3), A0A097NUQ3 (AtOSCA2.4), A0A097NUQ7 (AtOSCA2.5), Q9C8G5 (AtOSCA3.1), A0A097NUT0 (AtOSCA4.1), Q91YT8 (MmTMEM63A), Q5T3F8 (HsTMEM63B), and Q9P1W3 (HsTMEM63C), Q8R4P5 (MmTMC1), Q8R4P4 (MmTMC2), Q7TQ69 (MmTMC3), Q7TQ65 (MmTMC4), Q32NZ6 (MmTMC5), Q7TN60 (MmTMC6), Q8C428 (MmTMC7), and Q7TN58 (MmTMC8). A subset of the alignment was selected for Figure 1—figure supplement 1. Structural models were obtained from the PDB, aligned, and visualized using Pymol (Schrödinger, New York, NY).

Cell culture

Request a detailed protocol

The HEK293T cell line was authenticated and tested negative of mycoplasma contamination by the Duke Life Science Facility. The TMEM16F KO HEK293T cell line was generated by the Duke Functional Genomics Core and characterized in previous studies (Le et al., 2019b; Le et al., 2019a). All cells were cultured with DMEM (Dulbecco's Modified Eagle Medium) (Gibco/Thermo Fisher Scientific, Waltham, MA, Cat. #11995-065) supplemented with 10% fetal bovine serum (Sigma-Aldrich, St. Louis, MO, Cat. # F2442) and 1% penicillin/streptomycin (Gibco/Thermo Fisher Scientific, Waltham, MA, Cat. #15-140-122) at 37°C in 5% CO2-95% air. Cells tested negative for mycoplasma contamination (Applied Biological Materials, Vancouver, Canada, Cat. #G239).

Transfection

Request a detailed protocol

Plasmids were transiently transfected into TMEM16F KO HEK 293T cells by using X-tremeGENE9 (MilliporeSigma, Burlington, MA, Cat. #XTG9-RO), X-tremeGENE360 (MilliporeSigma, Burlington, MA, Cat. #XTG360-RO), or Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA, Cat. # 11668027). Media was replaced with calcium-free DMEM (Gibco/Thermo Fisher Scientific, Waltham, MA, Cat. #21068-028) 3–4 hr after transfection. The transfected cells were imaged or patched 18–24 or 18–48 hr after transfection, respectively.

Fluorescence imaging of scramblase-mediated PS exposure

Request a detailed protocol

A Zeiss 780 inverted confocal microscope was used to monitor scramblase activity in live cells using the methods described in previous publications (Zhang et al., 2020; Le et al., 2019b; Le et al., 2019a). 18–24 hr after transfection, the cells were incubated in AnV buffer (1:175 dilution of the fluorescently tagged AnV (Biotium, Fremont, CA, Cat. #29011) in Hank’s balanced salt solution) immediately before and throughout the duration of the imaging experiment without a formal incubation period. For Figure 3—figure supplement 1, 10 µM Ani9 (Sigma-Aldrich, St. Louis, MO, Cat. #SML1813) was added into the AnV solution before application. Spontaneous PS-positive cells were readily labeled by fluorescently tagged AnV. Results were quantified as a percentage of PS-positive cells among all GFP-positive cells. For osmolarity activation, 2 mM CaCl2 in ddH2O was added to the AnV buffer at a 2:1 ratio. The final osmolarity was ~120 mOsm/kg as measured by a micro-osmometer (Advanced Instrument, Norwood, MA). Cells overexpressing WT or mutant were treated with low osmolarity AnV buffer, and the scramblase activity was measured by recording fluorescent AnV surface accumulation at 5- (Figure 5) or 7- (Figure 6—figure supplement 1) s intervals. A custom MATLAB code was used to quantify AnV signal and is available at GitHub (https://github.com/yanghuanghe/scrambling_activity, copy archived at Huanghe, 2020; Le et al., 2019a).

Electrophysiology

All electrophysiology recordings were conducted using an Axopatch 200B amplifier with the signal digitally sampled at 10 kHz using an Axon Digidata 1550A (Molecular Devices, Inc, San Jose, CA). All electrophysiology recordings were carried out at room temperature 18–48 hr after transfection. Glass pipettes were pulled from borosilicate capillaries (Sutter Instruments, Novato, CA) and fire-polished using a microforge (Narishige, Amityville, NY). Pipettes had resistances of 2–4 MΩ in the physiological bath solution.

Inside-out patch clamp recordings

Request a detailed protocol

The pipette solution (external) contained 140 mM NaCl, 10 mM HEPES (Hydroxyethylpiperazine Ethane Sulfonic Acid), 2 mM MgCl2, adjusted to pH 7.3 (NaOH), and the bath solution (internal) contained 140 mM NaCl, 10 mM HEPES, 5 mM EGTA, adjusted to pH 7.3 (NaOH). OSCA1.2 WT and mutations were held at constant pressure administrated using a syringe calibrated with a manometer, similar to a previous study (Tsuchiya et al., 2018). In our experience, the constant pressure and variable voltage protocol achieved more consistent measurements. Patches were held at a membrane potential of −60  mV and at the indicated pressure, then stimulated using the indicated voltage protocol, taking advantage of the larger OSCA1.2 currents elicited by depolarizing potentials.

Cell-attached patch clamp recordings

Request a detailed protocol

For TMEM63 recordings, cell-attach mode was used instead of inside-out to avoid breaking the giga-ohm seal under higher applied pressure and voltage. The bath solution contained (in mM): 140 KCl, 10 HEPES, 2 MgCl2, 10 glucose, pH 7.3 adjusted with KOH. The pipette solution contained (in mM): 130 NaCl, 5 KCl, 10 HEPES, 10 TEA-Cl (Tetraethylammonium chloride), 1 CaCl2, 1 MgCl2, pH 7.3 (with NaOH). The mechano-activated current was evoked with a 200-ms pressure pulse at −80 mmHg using a high-speed pressure clamp system (ALA Scientific Instruments, Farmingdale, NY, Cat. #HSPC-1). The membrane potential inside the patch was held at −60 mV. The voltage pulse alone was run first followed by voltage pulse with pressure. The mechanosensitive current was obtained by subtracting the voltage pulse from the voltage pulse with pressure.

Whole-cell TMEM16 patch clamp recordings

Request a detailed protocol

The pipette solution (internal) contained 140 mM CsCl, 10 mM HEPES, 5 mM EGTA, adjusted to pH 7.3 (CsOH), and the bath solution (external) contained 140 mM NaCl, 10 mM HEPES, 5 mM EGTA, adjusted to pH 7.3 (NaOH). Patches were held at a membrane potential of −60  mV, then stimulated using the indicated voltage protocol. For Figure 3—figure supplement 1, 10µM Ani9 (Sigma-Aldrich, St. Louis, MO, Cat. #SML1813) was added to bath solution and applied using a perfusion manifold (ALA Scientific Instruments, Cat. #ALA-VM8).

Whole-cell hypotonic patch clamp recordings

Request a detailed protocol

The pipette solution (internal) contained 140 mM Na gluconate, 10 mM HEPES, 1 mM MgCl2, and 0.2 mM EGTA, adjusted to pH 7.3 (NaOH). Cells seeded on a small section of cover glass were placed in the bath solution (external) containing 140 mM Na gluconate, 10 mM HEPES, 1 mM MgCl2, and 0.2 mM EGTA adjusted to pH 7.3 (NaOH). Patches were held at a membrane potential of −60 mV, then stimulated using the indicated voltage protocol using a 2-s sweep protocol. After five sweeps, ddH2O was slowly and gently added to the bath at a ratio of 2:1 using a hypodermic needle to induce hypotonic cell swelling. The sweep protocol continued for a minimum of 4 min after hypotonic stimulation.

Data analysis for electrophysiology

Request a detailed protocol

All data analysis was performed using Clampfit (Molecular Devices, Inc, San Jose, CA), Excel (Microsoft, Redmond, WA), MATLAB (MathWorks, Portola Valley, CA), and Prism software (GraphPad, Boston, MA). Individual G–V curves were fitted with a Boltzmann function,

(1) G=Gmax1+e-ZFV-V0.5RT

where Gmax denotes the fitted value for maximal conductance, V0.5 denotes the voltage of half-maximal activation of conductance, Z denotes the net charge moved across the membrane during the transition from the closed to the open state, and F denotes the Faraday constant.

Model building and MD simulations

Request a detailed protocol

The initial model of OSCA1.2 for atomistic MD simulations and homology modeling was constructed based on PDB: 6MGV. Missing segments were modeled using the Swiss-PDB server (https://swissmodel.expasy.org/). The OSCA1.2 L438K mutation structure was modeled using 6MGV as the template, and the RMSD between the mutation and the original structure of 6MGV was less than 0.5 Å. To mimic the state of transmembrane proteins in the plasma membrane, the model was inserted into a lipid membrane consisting of POPC (Palmitoyl oleoyl phosphatidyl choline) (exoplasmic leaflet) and 1:1 POPS/POPE (Palmitoyl oleoyl phosphatidylserine/palmitoyl oleoyl phosphatidylethanolamine) (cytoplasmic leaflet) using the CHARMM-GUI web server (Jo et al., 2008). To neutralize charge levels in the simulation system, 150 mM KCl solution was added. The final simulation boxes contained ~600 lipid molecules, with ~290,000 atoms. During the simulation, the CHARMM36m all-atom force field was used. The temperature was maintained at 310 K to maintain the fluidity of lipids. During the atomistic simulation, long-range electrostatic interactions were described by the Particle Mesh Ewald algorithm with a cutoff of 12 Å. Van der Waals interactions were cut off at 12 Å and the MD time step was set at 2 fs. All atomistic MD simulations were conducted using GROMACS-2022.3 software package and enhanced sampling MD simulations were conducted using GROMACS-2022.3 with the COLVARS module. Before the production process, the system performed 5000 steps of energy minimization. Subsequently, all systems were equilibrated using the default equilibration parameters provided by CHARMM-GUI. After equilibration, the velocities of all atoms were randomly generated, and the product simulation system consisted of a 2-stage MD simulation. In step 1, the system used a conventional MD simulation run for 200 ns. In step 2, the system used metadynamics simulations run for ~700 ns. During enhanced sampling MD simulations, the sampling interval of the collective variable was adjusted to adapt the lipid movement trajectory. VMD and PYMOL were used to analyze simulation trajectories and water occupancy.

Materials availability

Request a detailed protocol

All plasmids generated in this study can be requested from the corresponding authors.

Data availability

A custom MATLAB code to quantify AnV signal is available at GitHub, copy archived at Huanghe, 2020. All other data are included in the manuscript and/or Supplementary files. Source data files have been provided for Figures 1–6.

References

Article and author information

Author details

  1. Augustus J Lowry

    Department of Biochemistry, Duke University School of Medicine, Durham, United States
    Contribution
    Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    Contributed equally with
    Pengfei Liang and Yang Zhang
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3524-3008
  2. Pengfei Liang

    Department of Biochemistry, Duke University School of Medicine, Durham, United States
    Contribution
    Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology
    Contributed equally with
    Augustus J Lowry and Yang Zhang
    Competing interests
    No competing interests declared
  3. Mo Song

    Institute of Molecular Physiology, Shenzhen Bay Laboratory, Guangdong, China
    Contribution
    Data curation, Software, Formal analysis, Visualization, Methodology
    Competing interests
    No competing interests declared
  4. Yuichun Wan

    Department of Biochemistry, Duke University School of Medicine, Durham, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  5. Zhen-Ming Pei

    Department of Biology, Duke University, Durham, United States
    Contribution
    Resources
    Competing interests
    No competing interests declared
  6. Huanghe Yang

    1. Department of Biochemistry, Duke University School of Medicine, Durham, United States
    2. Department of Neurobiology, Duke University School of Medicine, Durham, United States
    Contribution
    Conceptualization, Resources, Data curation, Supervision, Funding acquisition, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    huanghe.yang@duke.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9521-9328
  7. Yang Zhang

    1. Department of Biochemistry, Duke University School of Medicine, Durham, United States
    2. Institute of Molecular Physiology, Shenzhen Bay Laboratory, Guangdong, China
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Investigation, Methodology, Writing – review and editing
    Contributed equally with
    Augustus J Lowry and Pengfei Liang
    For correspondence
    zhangyang@szbl.ac.cn
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3625-9965

Funding

National Institute of General Medical Sciences (R21GM146152)

  • Huanghe Yang

National Institute of General Medical Sciences (DP2GM126898)

  • Huanghe Yang

National Institute of General Medical Sciences (R35GM153196)

  • Huanghe Yang

Duke University (Science & Technology SPARK Award)

  • Huanghe Yang

National Science Foundation (DGE 2139754)

  • Augustus J Lowry

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by the National Institutes of Health (DP2GM126898, R21GM146152, R35GM153196, and Duke Science & Technology SPARK Award to HY) and the National Science Foundation Graduate Research Fellowship Program (DGE 2139754 to AJL). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

Version history

  1. Preprint posted:
  2. Sent for peer review:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. Version of Record published:

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.96957. This DOI represents all versions, and will always resolve to the latest one.

Copyright

© 2024, Lowry, Liang et al.

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

  • 1,287
    views
  • 94
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Augustus J Lowry
  2. Pengfei Liang
  3. Mo Song
  4. Yuichun Wan
  5. Zhen-Ming Pei
  6. Huanghe Yang
  7. Yang Zhang
(2024)
TMEM16 and OSCA/TMEM63 proteins share a conserved potential to permeate ions and phospholipids
eLife 13:RP96957.
https://doi.org/10.7554/eLife.96957.3

Share this article

https://doi.org/10.7554/eLife.96957

Further reading

    1. Structural Biology and Molecular Biophysics
    Johannes Elferich, Lingli Kong ... Nikolaus Grigorieff
    Research Advance

    Images taken by transmission electron microscopes are usually affected by lens aberrations and image defocus, among other factors. These distortions can be modeled in reciprocal space using the contrast transfer function (CTF). Accurate estimation and correction of the CTF is essential for restoring the high-resolution signal in cryogenic electron microscopy (cryoEM). Previously, we described the implementation of algorithms for this task in the cisTEM software package (Grant et al., 2018). Here we show that taking sample characteristics, such as thickness and tilt, into account can improve CTF estimation. This is particularly important when imaging cellular samples, where measurement of sample thickness and geometry derived from accurate modeling of the Thon ring pattern helps judging the quality of the sample. This improved CTF estimation has been implemented in CTFFIND5, a new version of the cisTEM program CTFFIND. We evaluated the accuracy of these estimates using images of tilted aquaporin crystals and eukaryotic cells thinned by focused ion beam milling. We estimate that with micrographs of sufficient quality CTFFIND5 can measure sample tilt with an accuracy of 3° and sample thickness with an accuracy of 5 nm.

    1. Structural Biology and Molecular Biophysics
    Sneha Menon, Subinoy Adhikari, Jagannath Mondal
    Research Article

    The mis-folding and aggregation of intrinsically disordered proteins (IDPs) such as α-synuclein (αS) underlie the pathogenesis of various neurodegenerative disorders. However, targeting αS with small molecules faces challenges due to the lack of defined ligand-binding pockets in its disordered structure. Here, we implement a deep artificial neural network-based machine learning approach, which is able to statistically distinguish the fuzzy ensemble of conformational substates of αS in neat water from those in aqueous fasudil (small molecule of interest) solution. In particular, the presence of fasudil in the solvent either modulates pre-existing states of αS or gives rise to new conformational states of αS, akin to an ensemble-expansion mechanism. The ensembles display strong conformation-dependence in residue-wise interaction with the small molecule. A thermodynamic analysis indicates that small-molecule modulates the structural repertoire of αS by tuning protein backbone entropy, however entropy of the water remains unperturbed. Together, this study sheds light on the intricate interplay between small molecules and IDPs, offering insights into entropic modulation and ensemble expansion as key biophysical mechanisms driving potential therapeutics.