FRET Experiments Research Articles

Single Molecule FRET Reveals Lipid Induced Conformational Changes in Cytoplasmic Domain of Kir2.1

Anionic phospholipids are common regulators of ion channel activity, including Kir2.1, an inward rectifying potassium channel crucial for setting and maintaining resting membrane potential in many tissues, which is activated specifically by PIP2 and phosphatidyl glycerol (PG). Crystal structures indicate that a simple upward translocation of the cytoplasmic Kir domain is induced by PIP2 and PG binding. However, recent computational and single-molecule FRET experiments show more complex intra-domain movements and rotation in the Kir domain of prokaryotic KirBac1.1 upon PIP2 binding. Single molecule FRET measurements on eukaryotic human Kir2.1 indicate similar complex motions that are not detected in static crystal structures. Our measurements indicate dynamic conformational ‘breathing’ throughout the Kir domain, suggesting greater flexibility in the intracellular domain of ion channels than may be apparent in structural studies. Specifically, our data indicate a PIP2 dependent shifts in average single molecule FRET (smFRET) distributions, with outward motions seen prior to the slide helix (residue 59) and an inward motion seen at the base (residues 286 and 288) consistent with a twisting motion of the Kir domain that was previously shown to be associated with opening in KirBac1.1. In the absence of PIP2, smFRET signals are very dynamic with relatively large (>20A) motions at specific residues. PIP2 stabilizes the conformational dynamics of the Kir domain, with less transitions between structural states, and reduction in the occurrence of large conformational transitions.

Does Membrane Asymmetry Affect Nanoparticle-Membrane Interactions

Plasma membrane damage is one of the primary mechanisms through which engineered nanomaterials cause cytotoxicity. Despite a number of studies, the role of the structural features of the plasma membrane, in particular membrane phospholipid asymmetry in regulating nanoparticle-membrane interactions has remained obscure. In the current study, the role of membrane phospholipid asymmetry, differences in the lipid composition of exofacial and cytofacial leaflet of the plasma membrane in regulating nanoparticle-membrane interactions was investigated. Vesicles mimicking the exofacial (Vexo) and cytofacial leaflets (Vcyto) of the plasma membrane of erythrocytes were exposed to silica nanoparticles (50 nm), with different surface functional groups, plain, or carboxyl-, amine-, and polyethylene glycol (PEG)-modified and the effect of particles on vesicle integrity was examined. Symmetric and asymmetric vesicles with similar outer leaflet, but different inner leaflet lipid compositions, as well as live erythrocytes, were also used to elucidate the role of membrane asymmetry and individual leaflets in regulating nanoparticle-membrane interaction. Vexo and Vcyto vesicles showed significant difference in their interactions with nanoparticles. Vexo vesicles showed drastic leakage of intravesicular content after exposure to plain and amine-modified particles while none of the particles caused significant leakage in Vcyto vesicles. FRET experiments revealed that nanoparticle-induced disruption by plain and amine-modified particles occurs via pore formation. In agreement, FRET assay and GUV images revealed significant localization of plain and amine-modified particles at the surface of the Vexo vesicles. However, no localization was observed with the Vcyto vesicles. Interestingly, symmetric and asymmetric vesicles were disrupted similarly by plain and amine-modified particles. Importantly, hemolysis was observed after the incubation of erythrocytes with plain and amine-modified particles, consistent with leakage assays using Vexo vesicles. In conclusion, these results suggest that the exofacial leaflet of the plasma membrane is the primary regulator of nanoparticle-induced membrane damage in erythrocytes.

Minimum degree of overlap between IL-9R and IL-2R on human T lymphoma cells: A quantitative CLSM and FRET analysis.

The heterodimeric receptor complex of IL-9 consists of the cytokine-specific α-subunit and the common γc -chain shared with other cytokines, including IL-2, a central regulator of T cell function. We have shown previously the bipartite spatial relationship of IL-9 and IL-2 receptors at the surface of human T lymphoma cells: in addition to common clusters, expression of the two receptor kinds could also be observed in segregated membrane areas. Here we analyzed further the mutual cell surface organization of IL-9 and IL-2 receptors. Complementing Pearson correlation data with co-occurrence analysis of confocal microscopic images revealed that a minimum degree of IL-9R/IL-2R co-localization exists at the cell surface regardless of the overall spatial correlation of the two receptor kinds. Moreover, our FRET experiments demonstrated molecular scale assemblies of the elements of the IL-9/IL-2R system. Binding of IL-9 altered the structure and/or composition of these clusters. It is hypothesized, that by sequestering receptor subunits in common membrane areas, the overlapping domains of IL-9R and IL-2R provide a platform enabling both the formation of the appropriate receptor complex as well as subunit sharing between related cytokines. © 2018 International Society for Advancement of Cytometry.

Association of the IP3R to STIM1 provides a reduced intraluminal calcium microenvironment, resulting in enhanced store-operated calcium entry

The involvement of inositol trisphosphate receptor (IP3R) in modulating store-operated calcium entry (SOCE) was established many years ago. Nevertheless, the molecular mechanism responsible for this observation has not been elucidated to this date. In the present study we show that IP3R associates to STIM1 upon depletion of the endoplasmic reticulum (ER) by activation of the inositol trisphosphate signaling cascade via G-protein coupled receptors. IP3R-STIM1 association results in enhanced STIM1 puncta formation and larger Orai-mediated whole-cell currents as well as increased calcium influx. Depleting the ER with a calcium ATPase inhibitor (thapsigargin, TG) does not induce IP3R-STIM1 association, indicating that this association requires an active IP3R. The IP3R-STIM1 association is only observed after IP3R activation, as evidenced by FRET experiments and co-immunoprecipitation assays. ER intraluminal calcium measurements using Mag-Fluo-4 showed enhanced calcium depletion when IP3R is overexpressed. A STIM1-GCaMP fusion protein indicates that STIM1 detects lower calcium concentrations near its EF-hand domain when IP3R is overexpressed when compared with the fluorescence reported by a GCaMP homogenously distributed in the ER lumen (ER-GCaMP). All these data together strongly suggest that activation of inositol trisphosphate signaling cascade induces the formation of the IP3R-STIM1 complex. The activated IP3R provides a reduced intraluminal calcium microenvironment near STIM1, resulting in enhanced activation of Orai currents and SOCE.

The run-on oligomer filament enzyme mechanism of SgrAI: Part 1. Assembly kinetics of the run-on oligomer filament

Filament or run-on oligomer formation by metabolic enzymes is now recognized as a widespread phenomenon having potentially unique enzyme regulatory properties and biological roles, and its dysfunction is implicated in human diseases such as cancer, diabetes, and developmental disorders. SgrAI is a bacterial allosteric type II restriction endonuclease that binds to invading phage DNA, may protect the host DNA from off-target cleavage activity, and forms run-on oligomeric filaments with enhanced DNA-cleavage activity and altered DNA sequence specificity. However, the mechanisms of SgrAI filament growth, cooperativity in filament formation, sequestration of enzyme activity, and advantages over other filament mechanisms remain unknown. In this first of a two-part series, we developed methods and models to derive association and dissociation rate constants of DNA-bound SgrAI in run-on oligomers and addressed the specific questions of cooperativity and filament growth mechanisms. We show that the derived rate constants are consistent with the run-on oligomer sizes determined by EM analysis and are most consistent with a noncooperative growth mode of the run-on oligomer. These models and methods are extended in the accompanying article to include the full DNA-cleavage pathway and address specific questions related to the run-on oligomer mechanism including the sequestration of DNA-cleavage activity and trapping of products.

Response to Comment on "Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water".

Best et al claim that we provide no convincing basis to assert that a discrepancy remains between FRET and SAXS results on the dimensions of disordered proteins under physiological conditions. We maintain that a clear discrepancy is apparent in our and other recent publications, including results shown in the Best et al comment. A plausible origin is fluorophore interactions in FRET experiments.

Dermatan sulfate epimerase 1 and dermatan 4-O-sulfotransferase 1 form complexes that generate long epimerized 4-O-sulfated blocks

During the biosynthesis of chondroitin/dermatan sulfate (CS/DS), a variable fraction of glucuronic acid is converted to iduronic acid through the activities of two epimerases, dermatan sulfate epimerases 1 (DS-epi1) and 2 (DS-epi2). Previous in vitro studies indicated that without association with other enzymes, DS-epi1 activity produces structures that have only a few adjacent iduronic acid units. In vivo, concomitant with epimerization, dermatan 4-O-sulfotransferase 1 (D4ST1) sulfates the GalNAc adjacent to iduronic acid. This sulfation facilitates DS-epi1 activity and enables the formation of long blocks of sulfated iduronic acid-containing domains, which can be major components of CS/DS. In this report, we used recombinant enzymes to confirm the concerted action of DS-epi1 and D4ST1. Confocal microscopy revealed that these two enzymes colocalize to the Golgi, and FRET experiments indicated that they physically interact. Furthermore, FRET, immunoprecipitation, and cross-linking experiments also revealed that DS-epi1, DS-epi2, and D4ST1 form homomers and are all part of a hetero-oligomeric complex where D4ST1 directly interacts with DS-epi1, but not with DS-epi2. The cooperation of DS-epi1 with D4ST1 may therefore explain the processive mode of the formation of iduronic acid blocks. In conclusion, the iduronic acid-forming enzymes operate in complexes, similar to other enzymes active in glycosaminoglycan biosynthesis. This knowledge shed light on regulatory mechanisms controlling the biosynthesis of the structurally diverse CS/DS molecule.

A folded viral noncoding RNA blocks host cell exoribonucleases through a conformationally dynamic RNA structure

Folded RNA elements that block processive 5' → 3' cellular exoribonucleases (xrRNAs) to produce biologically active viral noncoding RNAs have been discovered in flaviviruses, potentially revealing a new mode of RNA maturation. However, whether this RNA structure-dependent mechanism exists elsewhere and, if so, whether a singular RNA fold is required, have been unclear. Here we demonstrate the existence of authentic RNA structure-dependent xrRNAs in dianthoviruses, plant-infecting viruses unrelated to animal-infecting flaviviruses. These xrRNAs have no sequence similarity to known xrRNAs; thus, we used a combination of biochemistry and virology to characterize their sequence requirements and mechanism of stopping exoribonucleases. By solving the structure of a dianthovirus xrRNA by X-ray crystallography, we reveal a complex fold that is very different from that of the flavivirus xrRNAs. However, both versions of xrRNAs contain a unique topological feature, a pseudoknot that creates a protective ring around the 5' end of the RNA structure; this may be a defining structural feature of xrRNAs. Single-molecule FRET experiments reveal that the dianthovirus xrRNAs undergo conformational changes and can use "codegradational remodeling," exploiting the exoribonucleases' degradation-linked helicase activity to help form their resistant structure; such a mechanism has not previously been reported. Convergent evolution has created RNA structure-dependent exoribonuclease resistance in different contexts, which establishes it as a general RNA maturation mechanism and defines xrRNAs as an authentic functional class of RNAs.

Single-Molecule Imaging Reveals Conformational Manipulation of Holliday Junction DNA by the Junction Processing Protein RuvA.

Interactions between DNA and motor proteins regulate nearly all biological functions of DNA such as gene expression, DNA replication and repair, and transcription. During the late stages of homologous recombination (HR), the Escherichia coli recombination machinery, RuvABC, resolves the four-way DNA motifs called Holliday junctions (HJs) that are formed during exchange of nucleotide sequences between two homologous duplex DNA. Although the formation of the RuvA-HJ complex is known to be the first critical step in the RuvABC pathway, the mechanism for the binding interaction between RuvA and HJ has remained elusive. Here, using single-molecule fluorescence resonance energy transfer (smFRET) and ensemble analyses, we show that RuvA stably binds to the HJ, halting its conformational dynamics. Our FRET experiments in different ionic environments created by Mg2+ and Na+ ions suggest that RuvA binds to the HJ via electrostatic interaction. Further, while recent studies have indicated that the HR process can be modulated for therapeutic applications by selective targeting of the HJ by chemotherapeutic drugs, we investigated the effect of drug-modified HJ on binding. Using cisplatin as a proof-of-concept drug, we show that RuvA binds to the cisplatin-modified HJ as efficiently as to the unmodified HJ, demonstrating that RuvA accommodates for the cisplatin-introduced charges and/or topological changes on the HJ.

Alpha-synuclein mitochondrial interaction leads to irreversible translocation and complex I impairment

α-synuclein is involved in both familial and sporadic Parkinson's disease. Although its interaction with mitochondria has been well documented, several aspects remains unknown or under debate such as the specific sub-mitochondrial localization or the dynamics of the interaction. It has been suggested that α-synuclein could only interact with ER-associated mitochondria. The vast use of model systems and experimental conditions makes difficult to compare results and extract definitive conclusions. Here we tackle this by analyzing, in a simplified system, the interaction between purified α-synuclein and isolated rat brain mitochondria. This work shows that wild type α-synuclein interacts with isolated mitochondria and translocates into the mitochondrial matrix. This interaction and the irreversibility of α-synuclein translocation depend on incubation time and α-synuclein concentration. FRET experiments show that α-synuclein localizes close to components of the TOM complex suggesting a passive transport of α-synuclein through the outer membrane. In addition, α-synuclein binding alters mitochondrial function at the level of Complex I leading to a decrease in ATP synthesis and an increase of ROS production.

Incorporation of bridged nucleic acids into CRISPR RNAs improves Cas9 endonuclease specificity

Off-target DNA cleavage is a paramount concern when applying CRISPR-Cas9 gene-editing technology to functional genetics and human therapeutic applications. Here, we show that incorporation of next-generation bridged nucleic acids (2′,4′-BNANC[N-Me]) as well as locked nucleic acids (LNA) at specific locations in CRISPR-RNAs (crRNAs) broadly reduces off-target DNA cleavage by Cas9 in vitro and in cells by several orders of magnitude. Using single-molecule FRET experiments we show that BNANC incorporation slows Cas9 kinetics and improves specificity by inducing a highly dynamic crRNA–DNA duplex for off-target sequences, which shortens dwell time in the cleavage-competent, “zipped” conformation. In addition to describing a robust technique for improving the precision of CRISPR/Cas9-based gene editing, this study illuminates an application of synthetic nucleic acids.

Epoxide hydrolase 1 (EPHX1) hydrolyzes epoxyeicosanoids and impairs cardiac recovery after ischemia

Stimuli such as inflammation or hypoxia induce cytochrome P450 epoxygenase-mediated production of arachidonic acid-derived epoxyeicosatrienoic acids (EETs). EETs have cardioprotective, vasodilatory, angiogenic, anti-inflammatory, and analgesic effects, which are diminished by EET hydrolysis yielding biologically less active dihydroxyeicosatrienoic acids (DHETs). Previous in vitro assays have suggested that epoxide hydrolase 2 (EPHX2) is responsible for nearly all EET hydrolysis. EPHX1, which exhibits slow EET hydrolysis in vitro, is thought to contribute only marginally to EET hydrolysis. Using Ephx1-/-, Ephx2-/-, and Ephx1-/-Ephx2-/- mice, we show here that EPHX1 significantly contributes to EET hydrolysis in vivo Disruption of Ephx1 and/or Ephx2 genes did not induce compensatory changes in expression of other Ephx genes or CYP2 family epoxygenases. Plasma levels of 8,9-, 11,12-, and 14,15-DHET were reduced by 38, 44, and 67% in Ephx2-/- mice compared with wildtype (WT) mice, respectively; however, plasma from Ephx1-/-Ephx2-/- mice exhibited significantly greater reduction (100, 99, and 96%) of those respective DHETs. Kinetic assays and FRET experiments indicated that EPHX1 is a slow EET scavenger, but hydrolyzes EETs in a coupled reaction with cytochrome P450 to limit basal EET levels. Moreover, we also found that EPHX1 activities are biologically relevant, as Ephx1-/-Ephx2-/- hearts had significantly better postischemic functional recovery (71%) than both WT (31%) and Ephx2-/- (51%) hearts. These findings indicate that Ephx1-/-Ephx2-/- mice are a valuable model for assessing EET-mediated effects, uncover a new paradigm for EET metabolism, and suggest that dual EPHX1 and EPHX2 inhibition may represent a therapeutic approach to manage human pathologies such as myocardial infarction.

Protein distance and dynamics characteristics inferred from single-molecule FRET experiment

A novel method for characterizing proteins allows for determining the end-to-end distance, Flory exponent and radius of gyration with FRET data only, in a fraction of the time.

Integrative Molecular Modelling of Biomolecules Guided by FRET Experiments

Integrative Molecular Modelling of Biomolecules Guided by FRET Experiments

Sequence-Encoded Heterogeneity of Interactions Decouples Different Measures of Protein Sizes and Reconciles the Discrepant Inferences from Saxs versus Fret Experiments

Quantifying the conformational properties of folded proteins as a function of denaturant concentration is important for understanding the mechanism of protein folding. Additionally, quantifying the conformational properties of intrinsically disordered proteins (IDPs) can yield insights regarding the functions of IDPs. Radius of gyration, RG, which quantifies the average size of a protein, and end-to-end distance, REE, are two commonly used conformational descriptors as these values can be determined using small-angle x-ray scattering (SAXS) and single molecule Förster resonance energy transfer (smFRET), respectively. However, for several proteins, the RG profile determined directly from SAXS as a function of denaturant concentration differs from the RG profile inferred using REE values determined from smFRET and a given homopolymer model. Here, we show that this discrepancy is not due to any inherent flaws of either measurement, but instead arises as a result of sequence-specific interactions that exist within proteins under low denaturant conditions. Sequence-specific interactions induce heterogeneity within the conformational ensemble of a protein, which in turn decouples globally averaged (e.g., RG) and length-scale-specific (e.g., REE) conformational descriptors. The extent of decoupling of conformational descriptors, and thus the appropriateness of using a homopolymer model to convert between RG and REE, is sequence specific and is often unknown a priori. Here, we investigate the relationship between sequence-complexity and deviations from homopolymer models. Additionally, we show that by combining results from experimental measurements that yield distinct conformational descriptors with molecular simulations we can provide complete conformational descriptions of heterogeneous ensembles.

Identification of Wiskott-Aldrich syndrome protein (WASP) binding sites on the branched actin filament nucleator Arp2/3 complex

Arp2/3 complex nucleates branched actin filaments important for cellular motility and endocytosis. WASP family proteins are Arp2/3 complex activators that play multiple roles in branching nucleation, but little is known about the structural bases of these WASP functions, owing to an incomplete understanding of how WASP binds Arp2/3 complex. Recent data show WASP binds two sites, and biochemical and structural studies led to models in which the WASP C segment engages the barbed ends of the Arp3 and Arp2 subunits while the WASP A segment binds the back side of the complex on Arp3. However, electron microscopy reconstructions showed density for WASP inconsistent with these models on the opposite (front) side of Arp2/3 complex. Here we use chemical cross-linking and mass spectrometry (XL-MS) along with computational docking and structure-based mutational analysis to map the two WASP binding sites on the complex. Our data corroborate the barbed end and back side binding models and show one WASP binding site on Arp3, on the back side of the complex, and a second site on the bottom of the complex, spanning Arp2 and ARPC1. The XL-MS-identified cross-links rule out the front side binding model and show that the A segment of WASP binds along the bottom side of the ARPC1 subunit, instead of at the Arp2/ARPC1 interface, as suggested by FRET experiments. The identified binding sites support the Arp3 tail release model to explain WASP-mediated activating conformational changes in Arp2/3 complex and provide insight into the roles of WASP in branching nucleation.

Slow domain reconfiguration causes power-law kinetics in a two-state enzyme

Protein dynamics are typically captured well by rate equations that predict exponential decays for two-state reactions. Here, we describe a remarkable exception. The electron-transfer enzyme quiescin sulfhydryl oxidase (QSOX), a natural fusion of two functionally distinct domains, switches between open- and closed-domain arrangements with apparent power-law kinetics. Using single-molecule FRET experiments on time scales from nanoseconds to milliseconds, we show that the unusual open-close kinetics results from slow sampling of an ensemble of disordered domain orientations. While substrate accelerates the kinetics, thus suggesting a substrate-induced switch to an alternative free energy landscape of the enzyme, the power-law behavior is also preserved upon electron load. Our results show that the slow sampling of open conformers is caused by a variety of interdomain interactions that imply a rugged free energy landscape, thus providing a generic mechanism for dynamic disorder in multidomain enzymes.

Accepting from the best donor; analysis of long-lifetime donor fluorescent protein pairings to optimise dynamic FLIM-based FRET experiments.

FRET biosensors have proven very useful tools for studying the activation of specific signalling pathways in living cells. Most biosensors designed to date have been predicated on fluorescent protein pairs that were identified by, and for use in, intensity based measurements, however fluorescence lifetime provides a more reliable measurement of FRET. Both the technology and fluorescent proteins available for FRET have moved on dramatically in the last decade. Lifetime imaging systems have become increasingly accessible and user-friendly, and there is an entire field of biology dedicated to refining and adapting different characteristics of existing and novel fluorescent proteins. This growing pool of fluorescent proteins includes the long-lifetime green and cyan fluorescent proteins Clover and mTurquoise2, the red variant mRuby2, and the dark acceptor sREACh. Here, we have tested these donors and acceptors in appropriate combinations against the standard or recommended norms (EGFP and mTFP as donors, mCherry and either Ypet or Venus as acceptors) to determine if they could provide more reliable, reproducible and quantifiable FLIM-FRET data to improve on the dynamic range compared to other donors and breadth of application of biosensor technologies. These tests were performed for comparison on both a wide-field, frequency domain system and a multiphoton, TCSPC time domain FLIM system. Clover proved to be an excellent donor with extended dynamic range in combination with mCherry on both platforms, while mRuby2 showed a high degree of variability and poor FRET efficiencies in all cases. mTFP-Venus was the most consistent cyan-yellow pair between the two FLIM methodologies, but mTurquoise2 has better dynamic range and transfers energy consistently over time to the dark acceptor sRCh. Combination of mTFP-sRCh with Clover-mCherry would allow the simultaneous use of two FLIM-FRET biosensors within one sample by eliminating the crosstalk between the yellow acceptor and green donor emissions.

Human β3-Adrenoreceptor is Resistant to Agonist-Induced Desensitization in Renal Epithelial Cells

Background/Aims: We recently showed that the β₃-adrenoreceptor (β₃AR) is expressed in mouse kidney collecting ducts (CD) cells along with the type-2 vasopressin receptor (AVPR2). Interestingly, a single injection of a β₃AR selective agonist promotes a potent antidiuretic effect in mice. Before considering the feasibility of chronic β₃AR agonism to induce antidiuresis in vivo, we aimed to evaluate in vitro the signaling and desensitization profiles of human β₃AR. Methods: Human β₃AR desensitization was compared with that of human AVPR2 in cultured renal cells. Video imaging and FRET experiments were performed to dissect β₃AR signaling under acute and chronic stimulation. Plasma membrane localization of β₃AR, AVPR2 and AQP2 after agonist stimulation was studied by confocal microscopy. Receptors degradation was evaluated by Western blotting. Results: In renal cells acute stimulation with the selective β₃AR agonist mirabegron, induced a dose-dependent increase in cAMP. Interestingly, chronic exposure to mirabegron promoted a significant increase of intracellular cAMP up to 12 hours. In addition, a slow and slight agonist-induced internalization and a delayed downregulation of β₃AR was observed under chronic stimulation. Furthermore, chronic exposure to mirabegron promoted apical expression of AQP2 also up to 12 hours. Conversely, long-term stimulation of AVPR2 with dDAVP showed short-lasting receptor signaling, rapid internalization and downregulation and apical AQP2 expression for no longer than 3 h. Conclusions: Overall, we conclude that β₃AR is less prone than AVPR2 to agonist-induced desensitization in renal collecting duct epithelial cells, showing sustained cAMP production, preserved membrane localization and delayed degradation after 12 hours agonist exposure. These results may be important for the potential use of chronic pharmacological stimulation of β₃AR to promote antidiuresis overcoming in vivo renal concentrating defects caused by inactivating mutations of the AVPR2.

TMEM16A and TMEM16B channel proteins generate Ca2+-activated Cl− current and regulate melatonin secretion in rat pineal glands

Pinealocytes regulate circadian rhythm by synthesizing and secreting melatonin. These cells generate action potentials; however, the contribution of specific ion channels to melatonin secretion from pinealocytes remains unclear. In this study, the involvement and molecular identity of Ca2+-activated Cl- (ClCa) channels in the regulation of melatonin secretion were examined in rat pineal glands. Treatment with the ClCa channel blockers, niflumic acid or T16Ainh-A01, significantly reduced melatonin secretion in pineal glands. After pineal K+ currents were totally blocked under whole-cell patch clamp conditions, depolarization and subsequent repolarization induced a slowly activating outward current and a substantial inward tail current, respectively. Both of these current changes were dependent on intracellular Ca2+ concentration and inhibited by niflumic acid and T16Ainh-A01. Quantitative real-time PCR, Western blotting, and immunocytochemical analyses revealed that TMEM16A and TMEM16B were highly expressed in pineal glands. siRNA knockdown of TMEM16A and/or TMEM16B showed that both channels contribute to ClCa currents in pinealocytes. Conversely, co-expression of TMEM16A and TMEM16B channels or the expression of this tandem channel in HEK293 cells mimicked the electrophysiological characteristics of ClCa currents in pinealocytes. Moreover, bimolecular fluorescence complementation, FRET, and co-immunoprecipitation experiments suggested that TMEM16A and TMEM16B can form heteromeric channels, as well as homomeric channels. In conclusion, pineal ClCa channels are composed of TMEM16A and TMEM16B subunits, and these fluxes regulate melatonin secretion in pineal glands.