members of the tmem16 gene family are transmembrane proteins classified into two functional categories : ca - activated ion channels ( caputo et al . , 2008 ; schroeder et al . , 2008 ; yang et al . , 2008 ) and phospholipid scramblases ( suzuki et al . , 2010 , 2013 ; yang et al . , 2012 ; malvezzi et al . , tmem16a and tmem16b in this family form ca - activated cl channels ( caccs ) critical for various physiological processes , including regulation of smooth muscle contraction , transepithelial fluid secretion , and neuronal excitation , etc . in contrast , tmem16c through tmem16j may function as lipid scramblases ( suzuki et al . , 2013 ) , the defect of which causes a variety of diseases , including gnathodiaphyseal dysplasia of bones ( katoh and katoh , 2004 ; tsutsumi et al . , 2004 ; marconi et al . , 2013 ) , muscular dystrophy ( bolduc et al . , 2010 ; mahjneh et al . , 2010 ; hicks et al . , 2011 ) , and a blood coagulation problem known as scott syndrome ( suzuki et al . , 2010 ; lhermusier et al . , 2011 ; bevers and williamson , 2016 ) , to name a few . activation of tmem16a or tmem16b generates mostly cl - selective currents , whereas activation of tmem16f induces ionic currents contributed by anions ( grubb et al . , 2013 ; shimizu et al . , 2013 ) , and possibly also by cations ( yang et al . , 2012 ; adomaviciene et al . , 2013 ; yu et al . , 2015a ) . the functional dichotomy and different ion selectivity associated with these two types of tmem16 members are intriguing research subjects ( kunzelmann et al . , 2014 ; pedemonte and galietta , 2014 ; picollo et al . , 2015 ) . biochemical evidence suggests that a functional tmem16 molecule is composed of two subunits ( fallah et al . , 2011 ; sheridan et al . , 2011 the recently solved crystal structure of a tmem16 family member from the fungus nectria haematococca ( nhtmem16 ) confirmed that tmem16 molecules adopt a homodimeric architecture ( brunner et al . , 2014 ) , similar to that observed in members of the clc cl channel / transporter family ( dutzler et al . 1 a ) revealed that each subunit of the tmem16 molecule binds two ca ions via several chelating acidic residues ( brunner et al . , 2014 ) . mutations of the corresponding acidic residues in tmem16a greatly alter the apparent ca affinity in activating the channel : the half - activation concentration ( k1/2 ) of the wild - type tmem16a channel is in the submicromolar to micromolar range , whereas those of the mutants can increase by two to three orders of magnitude ! thus , these acidic residues , including the first two identified glutamate residues , e698 and e701 ( yu et al . , 2012 ) , along with the three additional acidic residues e650 , e730 , and d734 ( tien et al . , 2014 ) , participate in forming the ca - binding sites . constructing a covalently linked tandem dimer of the tmem16a channel . ( a ) structure of the nhtmem16 molecule , showing the positions of corresponding amino acid residues , the mutations of which alter the ca affinity ( e698 and e701 ) and current rectification ( k584 ) of the tmem16a channel . the nhtmem16 residues corresponding to k584 , e698 , and e701 of tmem16a are n378 , e503 , and e506 , respectively . ( b ) schematic diagrams and the nomenclature of various homodimeric and heterodimeric tmem16a constructs . red stars represent the e698c mutation . as tmem16 proteins adopt a dimeric architecture similar to that of the clc chloride channels / transporters , which are known to consist of two independent ion transport pathways , it is tempting to ask whether the tmem16a channel also consists of two pores , and if so , whether the two pores function independently . these questions were unanswered because the nhtmem16 structure did not provide definitive information for identifying the channel pore or pores . we constructed covalently linked tandem dimers of tmem16a and introduced mutations that alter ca affinity to only one subunit . by doing so , we were able to create tmem16a channels containing two subunits with different ca - binding affinities . our results demonstrate that activation of the two subunits appears to be mostly independent , although cooperativity between activation of the two subunits can not be ruled out . by introducing a pore mutation in only one subunit , we further demonstrate that the current generated by activating the different subunits may come through distinct pores . the wild - type mouse tmem16a ( mtmem16a ) cdna ( ncbi reference sequence : nm_001242349.1 ) , which is the  a  alternative splice variant ( ferrera et al . , 2009 ) , was a courtesy from l.y . jan ( university of california , san francisco , san francisco , ca ) . in comparison with the  a , c  splice variant used in the companion paper by lim et al . , the sequence referred in this paper is four residues shorter . for example , residues k584 and e698 referred in this paper correspond to k588 and e702 , respectively , in lim et al . the monomeric and tandem dimeric wild - type and mutant cdna constructs were subcloned in pegfp - n3 or pires2 expression vectors ( takara bio inc . ) . mutagenesis was achieved by using the quikchange ii site - directed mutagenesis kit ( agilent technologies ) . to make the tmem16a tandem dimer , the cdna of tmem16a without stop codon was subcloned into the 5 end of another tmem16a cdna already in the pires2 construct . the linker between the two subunits was thus encoded by the nucleotide sequence of the 5 untranslated region of the second subunit , leading to an insertion of 33 amino acids of the following sequence linking two subunits : tggpetiaheaqdagtpnsgdatgvvdgereat . for simplicity , the symbols wt and mt will be used to represent the wild - type and the affinity - altered mutant subunit , respectively . furthermore , wtwt represents the wild - type tmem16a channel expressed from the monomeric cdna construct , whereas the channel expressed from the monomeric e698c mutant cdna will be called e698ce698c ( the two subunits naturally assemble together ) . the tandem dimeric channels generated from linking two wt subunits , one e698c and one wt subunit , and two e698c subunits will be called wt - wt , e698c - wt , and e698c - e698c , respectively ( fig . the same principle of nomenclature is used for the e701c constructs . the k584q mutant represents the channel in which both wt subunits contain the k584q mutation . the symbols e698c(k584q)-wt or e698c - wt(k584q ) will be used when the mutation k584q is introduced specifically to the e698c or the wt subunit of the e698c - wt heterodimer . the aforementioned cdna constructs produced channels with a gfp attached to the c terminus of the channel proteins ( in pegfp - n3 construct ) or channels without a gfp tag but with a separate gfp protein expressed in the same cell . transfections of cdnas to human embryonic kidney ( hek ) 293 cells were performed using the lipofectamine transfection method ( ni et al . , 2014 ; yu et al . cells expressing channels were identified by the green fluorescence in an inverted microscope ( dm irb ; leica ) equipped with a fluorescent light source and a gfp filter ( chroma technology corp . ) . excised inside - out patch recordings from gfp - positive hek293 cells were conducted 2472 h after transfections . the extracellular ( pipette ) solution contained 140 mm nacl , 10 mm hepes , and 0.1 mm egta and had a ph of 7.4 ; this solution was defined to be the 0-ca solution . the intracellular solutions contained various concentrations of egta ( [ egta ] ) and [ ca ] , with free [ ca ] ranging from 0.1 m to 20 mm ( ni et al . , 2014 ; yu et al . , 2014 ) . to construct [ ca]-dependent activation curves , all intracellular solutions contained 10 mm hepes and had a final [ cl ] of 140 mm . [ egta ] of 0.1 mm was used for solutions with free [ ca ] between 5 m and 20 mm . for solutions containing 220 mm [ ca ] , the [ nacl ] required for making final 140 mm [ cl ] was reduced according to the required [ cacl2 ] . for example , for the solution containing 20 mm [ ca ] , 100 mm nacl , and 20.1 mm cacl2 were added in the solution in addition to 0.1 mm egta . solutions containing free [ ca ] < 5 m required higher ca - buffered capacity , and thus 1 mm [ egta ] was used in these solutions . the required amounts of total ca for generating specific free [ ca ] were calculated using the maxchelator program with constants from the nist database # 46 v8 available on the internet ( http://maxchelator.stanford.edu/caegta-nist.htm ) . all values of [ ca ] referred in this paper are free [ ca ] , whereas total [ ca ] added to the solutions are stated as [ cacl2 ] . the working solutions for constructing current - voltage relationship ( or i - v curve ) were made from four stock solutions , 1 nacl and 1 nmdg - so4 solutions , which contained 10 mm hepes and 0.1 mm egta with or without 0.12 mm [ cacl2 ] , ph 7.4 . the two 1 nacl solutions contained 140 mm nacl and thus were equivalent to the solutions with 0 and 20 m [ ca ] used in the dose  response experiments ( 140 mm nacl , 10 mm hepes , and 0.1 mm egta , with or without 0.12 mm [ cacl2 ] , ph 7.4 ) . the 1 nmdg - so4 solutions contained 140 mm nmdg buffered with sulfuric acid ( h2so4 ) to ph 7.4 . three working solutions , 1 , 0.6 , and 0.3 nacl solutions , were made by mixing the 1 nacl solution with the 1 nmdg - so4 solution in 1:0 , 0.6:0.4 , and 0.3:0.7 volume ratios , respectively , resulting in solutions containing 140 mm , 84 mm , and 42 mm [ nacl ] , respectively . in separate control experiments , nmdg and so4 between 1 nacl and 1 nmdg - so4 solutions calculated using the clampfit program ( molecular devices ) was only 2 mv . the differences of junction potentials between any two of the three intracellular solutions were thus negligible ( the largest junction potential was 1.6 mv between 1 and 0.3 nacl solutions ) and were not corrected in i - v curve experiments . recording electrodes were made from borosilicate glass capillaries ( world precision instruments ) using the pp830 electrode puller ( narishige ) . when filled with the extracellular pipette solution , the resistance of the electrodes was between 1.5 and 3.0 m. all experiments were performed using the axopatch 200b amplifier and the digidata1440 analogue - digital signal - converting board controlled by the pclamp10 software ( molecular devices ) . various solutions were delivered to the intracellular side of the excised inside - out patch using the sf-77 solution exchanger ( warner instruments ) . to construct dose  response curves , we used a three - pulse protocol described previously ( ni et al . , 2014 ) , in which the cytoplasmic side of the patch was sequentially exposed to saturating [ ca ] ( 20 m for the wild - type channel and 20 mm for the channels containing e698c or e701c mutations ) , a test [ ca ] , and then again the saturating [ ca ] ( fig . , the current amplitude in 0-ca solution was subtracted from that obtained in the presence of ca . the current activated by the test [ ca ] was then normalized to the mean of the two maximum currents activated by the saturating [ ca ] ( or the maximum [ ca ] in the case of e701c ) , and the normalized values were used for generating dose  response curves . to construct i - v curves , two voltage protocols were used in various intracellular [ nacl ] with or without 20 m [ ca ] ( figs . 6 and 7 ) . in one , a ramp voltage change ( duration 1.6 s ) from 80 to 80 mv was applied . the other voltage protocol used step voltage jumps from a holding voltage of 0 mv to a test voltage from 80 to 80 mv in a 20-mv step . instantaneous current within 10 ms of voltage jump to the test voltage was measured . for each patch , the recordings in 0-ca solutions were subtracted from that obtained in [ ca ] ( 20 m or 20 mm ) , and the resulting leak - subtracted i - v curves were used to determine the reversal potential ( fig . to average i - v curves from different membrane patches , the digitized data points from the ramp protocol were normalized to the last digitized point obtained at 80 mv ( fig . 7 , b and c ) . for the current obtained at different step voltages , the instantaneous currents measured at different voltages were normalized to that measured at 80 mv . these two voltage protocols gave nearly identical i - v curves as shown in fig . 7 b. experimental data were analyzed using clampfit 10.4 software ( molecular devices ) and origin pro 8.0 software ( originlab , co. ) . response results obtained in [ ca ] up to 10 mm , after being normalized to the current obtained at 20 mm [ ca ] , were fitted to the hill equation to evaluate the apparent ca affinity : inorm = 1/(1 + ( k1/2/[ca ] ) ) , where inorm is the normalized current and k1/2 and h are the fitted half - effective concentration and the hill coefficient , respectively . to evaluate the permeability ratio pcl / pna , reversal potentials ( erev ) of the i - v curves were determined in the presence of various intracellular solutions containing 1 ( 140 mm ) , 0.6 , and 0.3 [ nacl ] , whereas the extracellular solution contained 1 [ nacl ] . the relationship between erev and pcl / pna follows the goldman - hodgkin - katz equation : the degree of the i - v curve rectification was evaluated by a rectification index ( ri ) , which is defined as the ratio of the absolute current amplitudes at 80 and 80 mv in symmetrical 140 mm [ cl ] . thus , a linear i - v curve will have an ri close to 1 but ri < 1 for outwardly rectified i - v curves . fret experiments were performed according to the previously described spectrofret approach ( bykova et al . , the experiments were conducted with an inverted fluorescence microscope ( ix-81 ; olympus ) equipped with a 100-w mercury lamp . the intensity and the duration of fluorescent light exposure were controlled by a computer - driven neutral density filter and a shutter , using metamorph software ( molecular devices / invitrogen ) . a ccd camera ( roper 128b ) was used , and metamorph user - designed journals initiated automatic collection of fluorescence cell images and their spectral images obtained by feeding the cell images through a spectrograph ( acton spectrapro 2150i ) . for obtaining the spectra of the fret donor ( cerulean , abbreviated as  c  ) and fret acceptor ( eyfp , abbreviated as  y  ) , the excitation light wavelengths were controlled by two filter cubes ( chroma technology corp . ) filtering the excitation light for c ( cube i : excitation filter d436/20 and dichroic mirror 455dclp ) and y ( cube ii : excitation filter hq500/20 and dichroic mirror q515lp ) . background light spectra from the blank region in the same image were subtracted from the spectra of the fluorophores . fluorescence analyses were conducted as described previously ( bykova et al . , 2006 ; yu et al . , 2015b ) . theoretically , the efficiency of the energy transfer ( e ) depends on the distance , r , between the two fluorophores according to the frster equation , e = ( r0)/[(r0 ) + r ] , where r0 is a spectroscopically determined distance when e = 0.5 ( or 50% of fret efficiency ) . in practice , because the fret donor  containing subunit and the fret acceptor  containing subunit are coexpressed , the channel molecules can contain two fret donors , two fret acceptors , or one donor and one acceptor . for cells expressing many more acceptor - containing subunits but a very small number of donor - containing subunits , the apparent fret efficiency ( eapp ) will be very small because of the presence of many molecules containing two fret acceptors , whereas a cell with a high donor / acceptor ratio will give a value of eapp approaching the true fret efficiency . therefore , the fret efficiency depends on the ratio of the expressed donor and acceptor molecules in the cells . we calculated the fret efficiency from individual cells by the approach briefly described below ( for more detailed description , see bykova et al . from cells expressing only the fret acceptor ( y)tagged subunit , ratio a0 is determined as the ratio between the fluorescence excited by the light passing cube i ( f436 ) and the fluorescence excited by the cube ii  filtered light ( f500 ) . namely , ratio a0 = f436/f500 . from cells cotransfected with donor ( c)-containing subunits and acceptor ( y)-containing subunits , we determined another ratio , ratio a. in this case , the spectrum from cells expressing c only was first subtracted from the c excitation spectrum obtained from cells coexpressing c and y , resulting in a y - only spectrum , which includes two components : the y emission caused by fret ( f436 ) and the emission caused by direct excitation of y ( f436 ) . the fret ratio ( fr ) was defined as the ratio of a and a0 : fr = ratio a / ratio a0 = ( f436 + f436)/f436 = 1 + f436/f436 , and the apparent fret efficiency ( eapp ) is calculated according to eapp = ( y/c )  ( fr  1 ) , where c and y are molar extinction coefficients of c and y , respectively . finally , eapp obtained from each cell was plotted against the c / y ratio as shown in fig . the true fret efficiency is the asymptotic value of eapp when the c / y ratio approaches infinity . s1 shows dose  response curves of the wild - type tmem16a channel ( wtwt ) from 0.1 m to 10 mm [ ca ] . table s1 shows reversal potentials in millivolts ( mean  sem ) for the wild - type tmem16a channel and various channel constructs harboring the k584q mutation . the wild - type mouse tmem16a ( mtmem16a ) cdna ( ncbi reference sequence : nm_001242349.1 ) , which is the  a  alternative splice variant ( ferrera et al . , 2009 ) , was a courtesy from l.y . jan ( university of california , san francisco , san francisco , ca ) . in comparison with the  a , c  splice variant used in the companion paper by lim et al . , the sequence referred in this paper is four residues shorter . for example , residues k584 and e698 referred in this paper correspond to k588 and e702 , respectively , in lim et al . the monomeric and tandem dimeric wild - type and mutant cdna constructs were subcloned in pegfp - n3 or pires2 expression vectors ( takara bio inc . ) . mutagenesis was achieved by using the quikchange ii site - directed mutagenesis kit ( agilent technologies ) . to make the tmem16a tandem dimer , the cdna of tmem16a without stop codon was subcloned into the 5 end of another tmem16a cdna already in the pires2 construct . the linker between the two subunits was thus encoded by the nucleotide sequence of the 5 untranslated region of the second subunit , leading to an insertion of 33 amino acids of the following sequence linking two subunits : tggpetiaheaqdagtpnsgdatgvvdgereat . for simplicity , the symbols wt and mt will be used to represent the wild - type and the affinity - altered mutant subunit , respectively . furthermore , wtwt represents the wild - type tmem16a channel expressed from the monomeric cdna construct , whereas the channel expressed from the monomeric e698c mutant cdna will be called e698ce698c ( the two subunits naturally assemble together ) . the tandem dimeric channels generated from linking two wt subunits , one e698c and one wt subunit , and two e698c subunits will be called wt - wt , e698c - wt , and e698c - e698c , respectively ( fig . the same principle of nomenclature is used for the e701c constructs . the k584q mutant represents the channel in which both wt subunits contain the k584q mutation . the symbols e698c(k584q)-wt or e698c - wt(k584q ) will be used when the mutation k584q is introduced specifically to the e698c or the wt subunit of the e698c - wt heterodimer . the aforementioned cdna constructs produced channels with a gfp attached to the c terminus of the channel proteins ( in pegfp - n3 construct ) or channels without a gfp tag but with a separate gfp protein expressed in the same cell . transfections of cdnas to human embryonic kidney ( hek ) 293 cells were performed using the lipofectamine transfection method ( ni et al . , 2014 ; yu et al . cells expressing channels were identified by the green fluorescence in an inverted microscope ( dm irb ; leica ) equipped with a fluorescent light source and a gfp filter ( chroma technology corp . ) . excised inside - out patch recordings from gfp - positive hek293 cells were conducted 2472 h after transfections . the extracellular ( pipette ) solution contained 140 mm nacl , 10 mm hepes , and 0.1 mm egta and had a ph of 7.4 ; this solution was defined to be the 0-ca solution . the intracellular solutions contained various concentrations of egta ( [ egta ] ) and [ ca ] , with free [ ca ] ranging from 0.1 m to 20 mm ( ni et al . , 2014 ; yu et al . , 2014 ) . to construct [ ca]-dependent activation curves , all intracellular solutions contained 10 mm hepes and had a final [ cl ] of 140 mm . [ egta ] of 0.1 mm was used for solutions with free [ ca ] between 5 m and 20 mm . for solutions containing 220 mm [ ca ] , the [ nacl ] required for making final 140 mm [ cl ] was reduced according to the required [ cacl2 ] . for example , for the solution containing 20 mm [ ca ] , 100 mm nacl , and 20.1 mm cacl2 were added in the solution in addition to 0.1 mm egta . solutions containing free [ ca ] < 5 m required higher ca - buffered capacity , and thus 1 mm [ egta ] was used in these solutions . the required amounts of total ca for generating specific free [ ca ] were calculated using the maxchelator program with constants from the nist database # 46 v8 available on the internet ( http://maxchelator.stanford.edu/caegta-nist.htm ) . all values of [ ca ] referred in this paper are free [ ca ] , whereas total [ ca ] added to the solutions are stated as [ cacl2 ] . the working solutions for constructing current - voltage relationship ( or i - v curve ) were made from four stock solutions , 1 nacl and 1 nmdg - so4 solutions , which contained 10 mm hepes and 0.1 mm egta with or without 0.12 mm [ cacl2 ] , ph 7.4 . the two 1 nacl solutions contained 140 mm nacl and thus were equivalent to the solutions with 0 and 20 m [ ca ] used in the dose  response experiments ( 140 mm nacl , 10 mm hepes , and 0.1 mm egta , with or without 0.12 mm [ cacl2 ] , ph 7.4 ) . the 1 nmdg - so4 solutions contained 140 mm nmdg buffered with sulfuric acid ( h2so4 ) to ph 7.4 . three working solutions , 1 , 0.6 , and 0.3 nacl solutions , were made by mixing the 1 nacl solution with the 1 nmdg - so4 solution in 1:0 , 0.6:0.4 , and 0.3:0.7 volume ratios , respectively , resulting in solutions containing 140 mm , 84 mm , and 42 mm [ nacl ] , respectively . in separate control experiments , nmdg and so4 the junction potential between 1 nacl and 1 nmdg - so4 solutions calculated using the clampfit program ( molecular devices ) was only 2 mv . the differences of junction potentials between any two of the three intracellular solutions were thus negligible ( the largest junction potential was 1.6 mv between 1 and 0.3 nacl solutions ) and were not corrected in i - v curve experiments . recording electrodes were made from borosilicate glass capillaries ( world precision instruments ) using the pp830 electrode puller ( narishige ) . when filled with the extracellular pipette solution , the resistance of the electrodes was between 1.5 and 3.0 m. all experiments were performed using the axopatch 200b amplifier and the digidata1440 analogue - digital signal - converting board controlled by the pclamp10 software ( molecular devices ) . various solutions were delivered to the intracellular side of the excised inside - out patch using the sf-77 solution exchanger ( warner instruments ) . to construct dose  response curves , we used a three - pulse protocol described previously ( ni et al . , 2014 ) , in which the cytoplasmic side of the patch was sequentially exposed to saturating [ ca ] ( 20 m for the wild - type channel and 20 mm for the channels containing e698c or e701c mutations ) , a test [ ca ] , and then again the saturating [ ca ] ( fig . , the current amplitude in 0-ca solution was subtracted from that obtained in the presence of ca . the current activated by the test [ ca ] was then normalized to the mean of the two maximum currents activated by the saturating [ ca ] ( or the maximum [ ca ] in the case of e701c ) , and the normalized values were used for generating dose  response curves . to construct i - v curves , two voltage protocols were used in various intracellular [ nacl ] with or without 20 m [ ca ] ( figs . 6 and 7 ) . in one , a ramp voltage change ( duration 1.6 s ) from 80 to 80 mv was applied . the other voltage protocol used step voltage jumps from a holding voltage of 0 mv to a test voltage from 80 to 80 mv in a 20-mv step . instantaneous current within 10 ms of voltage jump to the test voltage was measured . for each patch , the recordings in 0-ca solutions were subtracted from that obtained in [ ca ] ( 20 m or 20 mm ) , and the resulting leak - subtracted i - v curves were used to determine the reversal potential ( fig . to average i - v curves from different membrane patches , the digitized data points from the ramp protocol were normalized to the last digitized point obtained at 80 mv ( fig . 7 , b and c ) . for the current obtained at different step voltages , the instantaneous currents measured at different voltages were normalized to that measured at 80 mv . these two voltage protocols gave nearly identical i - v curves as shown in fig . 7 b. experimental data were analyzed using clampfit 10.4 software ( molecular devices ) and origin pro 8.0 software ( originlab , co. ) . response results obtained in [ ca ] up to 10 mm , after being normalized to the current obtained at 20 mm [ ca ] , were fitted to the hill equation to evaluate the apparent ca affinity : inorm = 1/(1 + ( k1/2/[ca ] ) ) , where inorm is the normalized current and k1/2 and h are the fitted half - effective concentration and the hill coefficient , respectively . to evaluate the permeability ratio pcl / pna , reversal potentials ( erev ) of the i - v curves were determined in the presence of various intracellular solutions containing 1 ( 140 mm ) , 0.6 , and 0.3 [ nacl ] , whereas the extracellular solution contained 1 [ nacl ] . the relationship between erev and pcl / pna follows the goldman - hodgkin - katz equation : the degree of the i - v curve rectification was evaluated by a rectification index ( ri ) , which is defined as the ratio of the absolute current amplitudes at 80 and 80 mv in symmetrical 140 mm [ cl ] . thus , a linear i - v curve will have an ri close to 1 but ri < 1 for outwardly rectified i - v curves . fret experiments were performed according to the previously described spectrofret approach ( bykova et al . , 2006 ; yu et al . , 2015b ) . in brief , the experiments were conducted with an inverted fluorescence microscope ( ix-81 ; olympus ) equipped with a 100-w mercury lamp . the intensity and the duration of fluorescent light exposure were controlled by a computer - driven neutral density filter and a shutter , using metamorph software ( molecular devices / invitrogen ) . a ccd camera ( roper 128b ) was used , and metamorph user - designed journals initiated automatic collection of fluorescence cell images and their spectral images obtained by feeding the cell images through a spectrograph ( acton spectrapro 2150i ) . for obtaining the spectra of the fret donor ( cerulean , abbreviated as c  ) and fret acceptor ( eyfp , abbreviated as  y  ) , the excitation light wavelengths were controlled by two filter cubes ( chroma technology corp . ) filtering the excitation light for c ( cube i : excitation filter d436/20 and dichroic mirror 455dclp ) and y ( cube ii : excitation filter hq500/20 and dichroic mirror q515lp ) . background light spectra from the blank region in the same image were subtracted from the spectra of the fluorophores . fluorescence analyses were conducted as described previously ( bykova et al . , 2006 ; yu et al . , 2015b ) theoretically , the efficiency of the energy transfer ( e ) depends on the distance , r , between the two fluorophores according to the frster equation , e = ( r0)/[(r0 ) + r ] , where r0 is a spectroscopically determined distance when e = 0.5 ( or 50% of fret efficiency ) . in practice , because the fret donor  containing subunit and the fret acceptor  containing subunit are coexpressed , the channel molecules can contain two fret donors , two fret acceptors , or one donor and one acceptor . for cells expressing many more acceptor - containing subunits but a very small number of donor - containing subunits , the apparent fret efficiency ( eapp ) will be very small because of the presence of many molecules containing two fret acceptors , whereas a cell with a high donor / acceptor ratio will give a value of eapp approaching the true fret efficiency . therefore , the fret efficiency depends on the ratio of the expressed donor and acceptor molecules in the cells . we calculated the fret efficiency from individual cells by the approach briefly described below ( for more detailed description , see bykova et al . , 2006 ; yu et al . ratio a0 is determined as the ratio between the fluorescence excited by the light passing cube i ( f436 ) and the fluorescence excited by the cube ii  filtered light ( f500 ) . namely , ratio a0 = f436/f500 . from cells cotransfected with donor ( c)-containing subunits and acceptor ( y)-containing subunits , we determined another ratio , ratio a. in this case , the spectrum from cells expressing c only was first subtracted from the c excitation spectrum obtained from cells coexpressing c and y , resulting in a y - only spectrum , which includes two components : the y emission caused by fret ( f436 ) and the emission caused by direct excitation of y ( f436 ) . the fret ratio ( fr ) was defined as the ratio of a and a0 : fr = ratio a / ratio a0 = ( f436 + f436)/f436 = 1 + f436/f436 , and the apparent fret efficiency ( eapp ) is calculated according to eapp = ( y/c )  ( fr  1 ) , where c and y are molar extinction coefficients of c and y , respectively . finally , eapp obtained from each cell was plotted against the c / y ratio as shown in fig . the true fret efficiency is the asymptotic value of eapp when the c / y ratio approaches infinity . fig . s1 shows dose  response curves of the wild - type tmem16a channel ( wtwt ) from 0.1 m to 10 mm [ ca ] . fig . table s1 shows reversal potentials in millivolts ( mean  sem ) for the wild - type tmem16a channel and various channel constructs harboring the k584q mutation . in the absence of mg , the apparent half - effective [ ca ] ( k1/2 ) for activating wild - type tmem16a channels is in the submicromolar range ( ni et al . , 2014 ) . 2 a illustrates that the wtwt channel can be significantly activated by 0.247 m [ ca ] ( fig . 2 a , left ) , whereas the current activated by 0.562 m was close to that activated by 20 m ( fig . mutations of several acidic residues coordinating ca binding changed the apparent ca affinity by two to three orders of magnitude ( yu et al . , 2012 ; . 2 b , in which 5 m [ ca ] activated very little current in e698ce698c , whereas 500 m [ ca ] was not a saturating concentration . these results corroborate the idea that e698 of tmem16a participates in forming the ca - binding site . channel currents were induced by a test [ ca ] pulse ( the middle pulse , indicated by red horizontal lines ) and two saturating [ ca ] pulses ( indicated by blue horizontal lines ) before and after the test pulse for inducing maximal current . because tmem16a is a dimeric channel , the contribution from each subunit activation to the current induction can be evaluated by using a covalently linked tandem dimeric channel in which the mutation that reduces the apparent ca affinity is introduced into one or both subunits . as shown in fig . 2 c , when two e698c - containing subunits were linked together ( namely the e698c - e698c tandem dimer ) , the ca sensitivity of this concatemer was very similar to that of the channel expressed from the monomeric cdna construct containing the e698c mutation ( namely , e698ce698c ) . in contrast , when e698c was introduced to only one subunit ( e698c - wt ) , the heterodimeric mutant channel was significantly activated by 5 m [ ca ] , whereas 500 m [ ca ] was still not a saturating concentration ( fig . the heterodimeric mutant channel with a reverse arrangement of the two subunits ( namely , wt - e698c ) also showed similar results ( not depicted ) . comprehensive evaluation of the ca sensitivity in the two subunits was obtained by constructing dose  response curves . the e698c mutants expressed either from the monomeric cdna construct ( e698ce698c ) or from the dimeric cdna construct ( e698c - e698c ) exhibited a ca - dependent activation curve easily fit to a single hill equation , with k1/2 in the hundreds of micromolar range . these results are consistent with those reported previously by others ( yu et al . , 2012 ; for the e698c - wt heterodimer , the dose - dependent activation can not be easily fit by a single hill equation . rather , the dose response activation appeared to occur in two different [ ca ] ranges , one in the range that activates wtwt and the other in the range similar to that of activating e698ce698c or e698c - e698c . indeed , the dose  response curve resembles a weighted sum of the wtwt and e698ce698c dose  response curves ( fig . this observation suggests that the two subunits are activated fairly independently ; activation of the two subunits can be separated if the two subunits have significantly different ca - binding affinities . to confirm this conclusion , we repeated similar experiments with a different point mutation , e701c ( fig . 3 b ) , which also altered the apparent ca affinity of the channel ( yu et al . , 2012 ; brunner et al . , 2014 ; tien et al . , 2014 ) . in fact , the e701c mutation led to an even greater decrease in ca affinity : the decrease in the apparent ca affinity was so drastic that the highest [ ca ] used in constructing the dose  response curve ( 10 mm ) was not even a saturating concentration as compared with an even higher [ ca ] ( 20 mm ) . nonetheless , the dose  response relationship obtained from the e701c - wt heterodimer also appears to be a weighted sum of the dose  response curves of the wtwt channel and the e701ce701c mutant ( fig . ( a ) comparison of the dose  response curves of wtwt , e698ce698c , e698c - e698c , and e698c - wt obtained at 40 mv ( top ) and 40 mv ( bottom ) . values of k1/2 ( m ) and h for wtwt , e698ce698c , and e698c - e698c are shown . the curves for the tandem heterodimeric channel ( e698c - wt ) in two regions were scaled from the dose  response curves of the wtwt channel ( i = a  iwtwt , black ) and the e698ce698c channel ( i = a + ( 1  a )  ie698ce698c , red ) , where the values of a ( 0.39 and 0.44 in the top and bottom panel , respectively ) were determined by averaging the values of inorm at [ ca ] = 3.13 m and [ ca ] = 5 m . ( b ) comparison of the dose  response curves of wtwt , e701ce701c , and e701c - wt obtained at 40 mv ( top ) and 40 mv ( bottom ) . values of the k1/2 and h for e701ce701c are not depicted because of the lack of current saturation in the dose the curves for e701c - wt were scaled from those of the wtwt channel ( i = b  iwtwt , black ) and the e701ce701c channel ( i = b + ( 1  b )  ie701ce701c , red ) . the values of b ( 0.542 and 0.566 in the top and bottom panel , respectively ) were obtained by averaging the values of inorm at [ ca ] = 3.13 m and [ ca ] = 5 m . all data points in a and b are depicted as mean  sem . the fitted values of k1/2 and hill coefficient are shown next to the fitted curves . it has been reported that tandem linkage of channel subunits may not ensure the stoichiometry of expressed channels ( mccormack et al . , 1992 ) . to exclude the possibility that the aforementioned data resulted from two separate populations of wtwt and e698ce698c channels because of the assembly of the same subunit ( wt or mt ) from two e698c - wt heterodimeric molecules , we conducted a fret analysis with cerulean ( c ) or eyfp ( y ) tagging on the c terminus of the channel protein . fret is a highly sensitive measurement of distance between the fret donor ( here c ) and fret acceptor ( y ; bykova et al . , 2006 ) . 4 shows that coexpressing wt - c and wt - y subunits resulted in a positive fret signal at high c / y ratios ( black symbols ) , indicating that the c and y tags from different subunits are close enough to allow for fret . this result is consistent with a previous fret study in tmem16a channels ( sheridan et al . , 2011 ) . in comparison , coexpression of e698c - wt - c and e698c - wt - y ( fig . these results indicate that the wt subunit in the heterodimer did not assemble with the wt subunit of another heterodimer to form a wtwt channel , and likewise , the two e698c subunits from two heterodimers do not form an e698ce698c mutant . comparison of the apparent fret efficiencies ( eapp ) between c / y - tagged monomeric wt constructs ( black ) or c / y - tagged tandem e698c - wt dimeric constructs ( red ) . the aforementioned results thus suggest that the current of the e698c - wt channel activated by [ ca ] < 20 m is caused by the activation of the high - affinity wt subunit . response curve of the e698c - wt channel should be similar to that of the wtwt or the wt - wt channel in low [ ca ] where activation of the e698c - containing subunit is minimal . 5 shows a comparison of the dose  response curves of the wt - wt channel and the e698c - wt channel in [ ca ] < 5 m . fitting the data points to the hill equation gave k1/2 = 655 nm and hill coefficient h = 1.8 for the e698c - wt heterodimer and k1/2 = 525 nm and h = 2.7 for the wt - wt channel , indicating that impairing the affinity in one subunit does not drastically affect the other subunit . however , the overall apparent ca sensitivity of the wt subunit of the heterodimer may be slightly diminished compared with that of the wt - wt channel because the heterodimer appears to have a larger k1/2 and a smaller hill coefficient . response curves of the wt - wt linked homodimer and the e698c - wt linked heterodimer at [ ca ] < 5 m . all current was normalized to that obtained with 5 m [ ca ] . the fitted values of k1/2 and hill coefficient are shown next to the fitted curves . the results shown above demonstrate independent activation of the two subunits of the tmem16a channel . however , whether activation of separate subunits is associated with different pore openings is still not clear . to address this question , we initially planned to generate , in a particular subunit , the k584q mutation , which was reported to alter the pcl / pna permeability ratio of the pore ( yang et al . unfortunately , we were unable to observe the change of pcl / pna caused by the k584q mutation ( fig . the reversal potentials measured in different intracellular [ nacl ] ( 140 mm vs. 84 mm or 140 mm vs. 42 mm ) were all similar among various channel constructs whether the k584q mutation was introduced in only one subunit or in both subunits ( fig . the failure of altering pcl / pna by the k584q mutation excluded our initial plan of using the discrimination of cl versus na as a pore property to address the pore architecture of tmem16a . comparison of na and cl permeation properties through the pores of the wtwt channel ( abbreviated as wt ) and the k584q mutant . ( a ) i - v curves of the wild - type channel ( left ) and the k584q mutant ( right ) in three [ nacl]i conditions : 140 mm ( black traces ) , 84 mm ( red traces ) , and 42 mm ( blue traces ) . the reversal potentials ( erev ) in each recording are better shown in the inset of each panel . ( b ) averaged values of erev of wtwt ( black ) , k584q mutant ( red ) , e698c - wt(k584q ) ( green ) , and e698c(k584q)-wt ( blue ) in 140 , 84 , and 42 mm [ nacl]i . values of erev of various channels in the same ionic conditions ( table s1 ) were not significantly different from one another ( p > 0.05 , one - way anova ) . in constructing i - v relationships , we noticed that with 20 m [ ca ] the i - v curve of the k584q mutant ( namely , k584q mutation in both wt subunit ) is outwardly rectified ( fig . the outward rectification of this mutant is not caused by gating change at different voltages because the i - v relationships obtained from a ramp protocol and from measuring the instantaneous current immediately after step voltage jumps are nearly identical . furthermore , the rectification of the i - v curve at a much higher [ ca ] ( 20 mm ) is nearly the same as that obtained in 20 m ( fig . 7 b ) . in comparison , the i - v curve of the wtwt channel is linear ( fig . 7 thus , if activation of individual subunits is associated with distinct pores , we expect that the e698c - wt heterodimer containing the k584q mutant in the e698c subunit , namely , e698c(k584q)-wt ( abbreviated as mt(k584q)-wt in fig . 7 c ) , will have the same i - v curve rectification as the e698c - wt heterodimer when only the wt subunit is activated by low [ ca ] ( such as 20 m ) . in contrast , when the k584q mutation was made in the wt subunit of the e698c - wt heterodimer , namely e698c - wt(k584q ) ( also abbreviated as mt - wt(k584q ) ) , the rectification of the i - v curve was significantly more rectified than that of the e698c - wt channel ( fig . pore properties of the wild - type tmem16a channel , k584q mutant , and the mt - wt(k584q ) and mt(k584q)-wt heterodimeric channels , where mt refers to the e698c subunit . wt(k584q ) and mt(k584q ) , respectively , represent the wt and e698c subunits harboring k584q mutation . ( a ) recordings of the k584q mutant in 20 m [ ca ] ( left ) in response to a voltage step family from 80 to 80 mv ( right ) . ( b ) i - v relationships of the k584q mutant obtained from a ramp voltage protocol ( 80 to 80 mv ) in 20 mm [ ca ] ( black ) and 20 m [ ca ] ( red ) and from a voltage step protocol like that shown in a ( red squares , n = 5 ) . for the two ramp i - v curves , the width of each curve represents the sem of the averaged current of each digitized data point in the i - v curve . ri = 0.45  0.01 ( n = 7 ) and 0.48  0.01 ( n = 5 ) for the i - v curves in 20 m [ ca ] and 20 mm [ ca ] , respectively . ( c ) comparison of the i - v curve rectifications between the wild - type channel and k584q mutant ( left ) and between e698c - wt ( abbreviated as mt - wt ; blue ) and mt(k584q)-wt ( red , middle ) or mt - wt(k584q ) ( red , right ) . the width of each i - v curve represents the sem of the averaged values . ris for wt , k584q , mt - wt , mt(k584q)-wt , and mt - wt(k584q ) were 1.00  0.02 , 0.45  0.01 , 0.78  0.06 , 0.78  0.05 , and 0.47  0.04 , respectively . members of the tmem16 family are known to consist of two identical subunits ( fallah et al . , 2011 ; sheridan et al . , 2011 ; tien et al . , 2013 ; brunner et al . , 2014 ) . the present study examines whether the two subunits of a tmem16a cacc can be independently activated by ca and whether activation of different subunits is associated with the opening of different pores . by linking two tmem16a subunits together , we are able to introduce the mutations that alter the ca - binding affinity to only one subunit . fret experiments showed that coexpressing e698c - wt - c and e698c - wt - y did not produce detectable fret signals between the fret donor and acceptor ( fig . 4 ) , rendering the possibility less likely that wtwt and mtmt channels are formed from assembling two heterodimeric molecules together . rather , it appears that expressing heterodimeric constructs produces a homogeneous population of channels containing one wt subunit and one mt subunit . interestingly , the ca dose  response curves in both the e698c - wt and e701c - wt heterodimers can be described by the weighted sum of two ca activation curves , one with a k1/2 similar to that of the wtwt channel in the submicromolar range and the other with a k1/2 in the hundreds of micromolar range like that of the mutant channel ( fig . 3 ) . the biphasic activation by ca can also be observed in wild - type tmem16a ( fig . s1 however , the biphasic activation in the wild - type channel ( namely , wtwt or wt - wt ) is very different from that of the heterodimeric channel in several aspects . first , the fractions of the high - affinity versus the low - affinity components of the dose  response curves are different between the wild - type channel and the heterodimer . in the wtwt channel , the low - affinity component of the dose  response curve ( namely , the fraction of the current activated by [ ca ] > 20 m ) is 2530% of that of the total current induced by 20 mm [ ca ] ( fig . the fraction of the low - affinity component is much larger , at least 50% , possibly up to 70% , of the maximal current obtained at [ ca ] = 20 mm ( fig . second , the k1/2 of the low - affinity component of the wild - type channel was 2 mm , whereas in the heterodimer the k1/2 appears to depend on the mutation made in the affinity - altered subunit because scaling the dose response curve of the mutant appears to describe the low - affinity component well ( fig . 3 ) . finally , the biphasic activation of the wild - type channel can never explain the different ca activation between the mutant homodimer ( such as e698ce698c or e698c - e698c ) and the heterodimer ( such as e698c - wt ) . for example , with 5 m [ ca ] , the current activated in e698ce698c ( or e698c - e698c ) is very small ( fig . 2 , b and c ) . however , by replacing one e698c subunit with the wt subunit ( the e698c - wt heterodimer ) , the current activated by the same [ ca ] is significantly increased ( fig . 2 d ) . in comparison , 5 m [ ca ] already activates at least 70% of the maximal current in the wtwt channel ( fig . thus , the most straightforward interpretation for the biphasic dose  response curve of the heterodimer is an independent activation of the two subunits of the heterodimeric channel instead of the biphasic activation of the wt subunit . the independent activation of the two tmem16a subunits does not necessarily mean that different subunit activation leads to the opening of distinct pores . it is arguable that only one pore could exist in a channel and that activation of one subunit contributes to a partial opening of the pore . to explore whether the activation of individual subunits is associated with the opening of separate pores , we first attempted to take advantage of the altered discrimination of cl from na by k584q mutation and link this mutation to individual subunits . the ion permeation through tmem16 proteins has been thought not to be particularly selective in discriminating anions from cations . for example , the current associated with the activation of tmem16f is contributed significantly by na ( yang et al . , 2012 ; adomaviciene et al . , 2013 ; yu et al . , 2015a ) , with a pna / pcl ratio estimated to be as high as 6.8 in one previous study ( yang et al . , however , other studies still consider tmem16f as more anion selective ( grubb et al . , 2013 ; shimizu et al . , 2013 ) , and a pna / pcl ratio of 0.3 had been estimated by grubb et al . compared with tmem16f , tmem16a is much more anion selective , but even so , a pcl / pna ratio of only 7 was reported ( yang et al . , 2012 ) . the mechanism underlying the different selectivity for anions versus cations between tmem16a and tmem16f was thought to be caused by a positively charged lysine residue , k584 , in tmem16a versus a corresponding uncharged glutamine residue , q559 , in tmem16f ; mutation of k584 in tmem16a into a glutamine residue ( k584q ) was shown to reduce the pcl / pna ratio from 7.0 to 4.1 , whereas mutating the corresponding residue from q to k in tmem16f ( q559k ) reduced the pna / pcl of tmem16f from 6.8 to 2.2 ( yang et al . , 2012 ) . these results thus suggested the possibility that k584 residue could be located in the tmem16a channel pore to discriminate anions from cations . 6 do not support the conclusion that k584 is critical for discriminating cl from na in the tmem16a channel pore . in fact , our results showed that the i - v curves of the wild - type channel and the k584q mutant have indistinguishable reversal potentials in various intracellular [ nacl ] ( fig . furthermore , the calculated pcl / pna ratio for the wild - type channel from our experiments is much higher than that reported previously . the reversal potentials of the wtwt channel in our measurements were 31.5  0.3 mv , 13.0  0.4 mv , and 0.5  0.5 mv in 0.3 , 0.6 , and 1 intracellular [ nacl ] , respectively ( fig . 0.6 solutions based on the goldman - hodgkin - katz equation would be 21.8 mv and 9.6 mv , 28.2 mv and 12.1 mv , and 29.8 mv and 12.7 mv when pcl / pna ratios are 7 , 30 , and 100 , respectively ( t = 295k ) . our results thus suggest that the permeability of na through the tmem16a pore is negligible in comparison with the cl permeability . studies from others have also suggested that tmem16a has a nearly ideal anionic selectivity ( yang et al . the reasons for the qualitative and quantitative differences between our results and those reported in yang et al . ( 2012 ) are not clear to us , but we suspect that differences in experimental methods , for example , using nmdg - so4 versus sucrose to replace intracellular nacl ( lim et al . , 2016 ) , and/or analyses ( such as subtraction of the endogenous leak current ; fig . although we are unable to confirm that residue k584 helps in discriminating anions versus cations for tmem16a , our experiments do suggest that k584 is a pore residue because of the prominent rectification of the i - v curve of the k584q mutant in symmetrical 140 mm [ nacl ] . the i - v curve of the wild - type tmem16a pore in a saturating [ ca ] is linear with a unity ri . in contrast , the i - v curve of the k584q mutant is outwardly rectified ( ri = 0.45 ; fig . this outward rectification in k584q is not caused by the voltage - dependent change of gating properties because the same rectification is observed in the i - v relationship constructed from instantaneous currents upon jumping membrane potential to various voltages ( fig . furthermore , with 20 mm [ ca ] , the i - v curve of the k584q mutant remains outwardly rectified ( fig . it is possible that k584 is located at the intracellular pore entrance , and therefore the positive charge on the sidechain helps recruit cl to the intracellular pore entrance , similar to the electrostatic control of the conductance of clc-0 by a lysine residue ( middleton et al . if there are separate pores associated with different subunits and if the i - v curve rectification caused by the k584q mutation is indeed a pore property , we expect that introducing k584q mutation to the wt subunit or to the e698c - containing subunit in the e698c - wt heterodimer would produce different degrees of i - v curve rectification when only the high - affinity wt subunit is activated . 7 c indeed shows the fulfillment of this prediction , indicating that activation of different subunits is associated with the opening of distinct pores . thus , it appears that the tmem16a channel may also contain multiple ( most likely two ) pores , and ca activation of an individual subunit opens the pore associated with that activated subunit . double - barrel  feature of clc channels ( miller , 1982 ; miller and white , 1984 ; ludewig et al . , 1996 ; middleton et al . , 1996 ) , except that the opening of the clc channel pores is apparently voltage dependent , whereas the tmem16a pores are opened by ligand binding . although the results presented in this study suggest that activation of different subunits of tmem16a opens distinct pores , the exact mechanism of tmem16a activation is still a subject of much speculation . the fact that the wild - type tmem16a channel has a dose  response curve with a hill coefficient h > 1 indicates cooperative binding of multiple ca in the channel activation process . response curve of the e698c - wt channel in the range of [ ca ] < 5 m still has a hill coefficient substantially greater than 1 ( fig . 5 ) indicates multiple ca - binding sites within one subunit , consistent with the high - resolution crystal structure of the nhtmem16 molecule . nonetheless , it appears that the cooperativity of channel opening may not all come from ca binding within one subunit because the dose  response curve of the e698c - wt in the low [ ca ] range ( < 5 m ) has a lower apparent affinity and a smaller hill coefficient than that of the wt - wt channel . this observation suggests that some cooperativity in the channel activation resulting from cross - talk of the two subunits can not be ruled out . finally , the activation of the wild - type channel by [ ca ] in the mm range ( fig . s1 ) may suggest low - affinity ca - binding sites not seen in the nhtmem16 crystal structure . in conclusion , by constructing covalent tandem dimers for tmem16a channels and introducing mutations to specific subunits , we altered the ca - binding affinity and induced the rectification of the channel current associated with individual subunits . our results suggest that the two subunits with different apparent ca affinities can be independently activated and that the activation of individual subunits is associated with the opening of separate pores . double - barrel  structure found in members of the clc channel / transporter family .
tmem16 lipid scramblases are formed from two identical subunits , but the pore architecture of tmem16 channels remains unclear . 
 jeng et al . 
 use tandem dimers of tmem16a subunits with different mutations and show that each subunit comprises an independent pore and contains more than one ca2 + -binding site .