Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases

Calcium-activated chloride channel TMEM16A opens via pi-helical transition in transmembrane segment 4

TMEM16A is a Ca2+-activated Cl- channel that has crucial roles in various physiological and pathological processes. However, the structure of the open state of the channel and the mechanism of Ca2+-induced pore opening have remained elusive. Using extensive molecular dynamics simulations, protein structure prediction, and patch-clamp electrophysiology, we demonstrate that TMEM16A opens a hydrated Cl--conductive pore via a pi-helical transition in transmembrane segment 4 (TM4). We also describe a coupling mechanism that links pi-helical transition and pore opening to the Ca2+-induced conformational changes in TMEM16A. Furthermore, we designed a pi-helix-stabilizing mutation (I551P) that facilitates TMEM16A activation, revealing atomistic details of the ion-conduction mechanism. Finally, AlphaFold2 structure predictions revealed the importance of the pi helix in TM4 to structure-function relations in TMEM16 and the related OSCA/TMEM63 family, further highlighting the relevance of dynamic pi helices for gating in various ion channels.

Simulation-based survey of TMEM16 family reveals that robust lipid scrambling requires an open groove

Biological membranes are complex and dynamic structures with different populations of lipids in their inner and outer leaflets. The Ca2+-activated TMEM16 family of membrane proteins plays an important role in collapsing this asymmetric lipid distribution by spontaneously, and bidirectionally, scrambling phospholipids between the two leaflets, which can initiate signaling and alter the physical properties of the membrane. While evidence shows that lipid scrambling can occur via an open hydrophilic pathway ("groove") that spans the membrane, it remains unclear if all family members facilitate lipid movement in this manner. Here we present a comprehensive computational study of lipid scrambling by all TMEM16 members with experimentally solved structures. We performed coarse-grained molecular dynamics (MD) simulations of 27 structures from five different family members solved under activating and non-activating conditions, and we captured over 700 scrambling events in aggregate. This enabled us to directly compare scrambling rates, mechanisms, and protein-lipid interactions for fungal and mammalian TMEM16s, in both open (Ca2+-bound) and closed (Ca2+-free) conformations with statistical rigor. We show that all TMEM16 structures thin the membrane and that the majority of scrambling (>90%) occurs at the groove only when TM4 and TM6 have sufficiently separated. Surprisingly, we also observed 60 scrambling events that occurred outside the canonical groove, over 90% of which took place at the dimer-dimer interface in mammalian TMEM16s. This new site suggests an alternative mechanism for lipid scrambling in the absence of an open groove.

Wide-ranging cellular functions of ion channels and lipid scramblases in the structurally related TMC, TMEM16 and TMEM63 families

Calcium (Ca2+)-activated ion channels and lipid scramblases in the transmembrane protein 16 (TMEM16) family are structurally related to mechanosensitive ion channels in the TMEM63 and transmembrane channel-like (TMC) families. Members of this structurally related superfamily share similarities in gating transitions and serve a wide range of physiological functions, which is evident from their disease associations. The TMEM16, TMEM63 and TMC families include members with important functions in the cell membrane and/or intracellular organelles such as the endoplasmic reticulum, membrane contact sites, endosomes and lysosomes. Moreover, some members of the TMEM16 family and the TMC family perform dual functions of ion channel and lipid scramblase, leading to intriguing physiological implications. In addition to their physiological functions such as mediating phosphatidylserine exposure and facilitation of extracellular vesicle generation and cell fusion, scramblases are involved in the entry and replication of enveloped viruses. Comparisons of structurally diverse scramblases may uncover features in the lipid-scrambling mechanisms that are likely shared by scramblases.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces phosphatidylserine externalization, a key step in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements underlying groove opening and how lipids reorganize outside the closed groove remain unknown. Here we directly visualize how lipids associate at the closed groove of Ca2+-bound fungal nhTMEM16 in nanodiscs using cryo-EM. Functional experiments pinpoint lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryo-EM structure determination.

Structure and function of the human apoptotic scramblase Xkr4

Phosphatidylserine externalization on the surface of dying cells is a key signal for their recognition and clearance by macrophages and is mediated by members of the X-Kell related (Xkr) protein family. Defective Xkr-mediated scrambling impairs clearance, leading to inflammation. It was proposed that activation of the Xkr4 apoptotic scramblase requires caspase cleavage, followed by dimerization and ligand binding. Here, using a combination of biochemical approaches we show that purified monomeric, full-length human Xkr4 (hXkr4) scrambles lipids. CryoEM imaging shows that hXkr4 adopts a novel conformation, where three conserved acidic residues create an electronegative surface embedded in the membrane. Molecular dynamics simulations show this conformation induces membrane thinning, which could promote scrambling. Thinning is ablated or reduced in conditions where scrambling is abolished or reduced. Our work provides insights into the molecular mechanisms of hXkr4 scrambling and suggests the ability to thin membranes might be a general property of active scramblases.

In or out of the groove? Mechanisms of lipid scrambling by TMEM16 proteins

Phospholipid scramblases mediate the rapid movement of lipids between membrane leaflets, a key step in establishing and maintaining membrane homeostasis of the membranes of all eukaryotic cells and their organelles. Thus, impairment of lipid scrambling can lead to a variety of pathologies. How scramblases catalyzed the transbilayer movement of lipids remains poorly understood. Despite the availability of direct structural information on three unrelated families of scramblases, the TMEM16s, the Xkrs, and ATG-9, a unifying mechanism has failed to emerge thus far. Among these, the most extensively studied and best understood are the Ca2+ activated TMEM16s, which comprise ion channels and/or scramblases. Early work supported the view that these proteins provided a hydrophilic, membrane-exposed groove through which the lipid headgroups could permeate. However, structural, and functional experiments have since challenged this mechanism, leading to the proposal that the TMEM16s distort and thin the membrane near the groove to facilitate lipid scrambling. Here, we review our understanding of the structural and mechanistic underpinnings of lipid scrambling by the TMEM16s and discuss how the different proposals account for the various experimental observations.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces the externalization of phosphatidylserine, a key molecule in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove, and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements involved in groove opening and of how lipids reorganize outside the closed groove remain unknown. Using cryogenic electron microscopy, we directly visualize how lipids associate at the closed groove of Ca2+-bound nhTMEM16 in nanodiscs. Functional experiments pinpoint the lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryoEM structure determination.

Recent developments in molecular modeling tools and applications related to pharmaceutical and biomedical research

In modern pharmaceutical and biomedical research, molecular modeling represents a useful tool to explore processes and their mechanistic bases at the molecular level. Integrating experimental and virtual analysis is a fruitful approach to study ligand-receptor interaction in chemical, biochemical and biological environments. In these fields, molecular docking and molecular dynamics are considered privileged techniques for modeling (bio)macromolecules and related complexes. This review aims to present the current landscape of molecular modeling in pharmaceutical and biomedical research by examining selected representative applications published in the last years and highlighting current topics and trends of this field. Thus, a systematic compilation of all published literature has not been attempted herein. After a brief overview of the main theoretical and computational tools used to investigate mechanisms at molecular level, recent applications of molecular modeling in drug discovery, ligand binding and for studying protein conformation and function will be discussed. Furthermore, specific sections will be devoted to the application of molecular modeling for unravelling enantioselective mechanisms underlying the enantioseparation of chiral compounds of pharmaceutical and biomedical interest as well as for studying new forms of noncovalent interactivity identified in biochemical and biological environments. The general aim of this review is to provide the reader with a modern overview of the topic, highlighting advancements and outlooks as well as drawbacks and pitfalls still affecting the applicability of theoretical and computational methods in the field of pharmaceutical and biomedical research.

Structure and Function of Calcium-Activated Chloride Channels and Phospholipid Scramblases in the TMEM16 Family

The transmembrane protein 16 (TMEM16) family consists of Ca2+-activated chloride channels and phospholipid scramblases. Ten mammalian TMEM16 proteins, TMEM16A-K (with no TMEM16I), and several non-mammalian TMEM16 proteins, such as afTMEM16 and nhTMEM16, have been discovered. All known TMEM16 proteins are homodimeric proteins containing two subunits. Each subunit consists of ten transmembrane helices with Ca2+-binding sites and a single ion-permeation/phospholipid transport pathway. The ion-permeation pathway and the phospholipid transport pathway of TMEM16 proteins have a wide intracellular vestibule, a narrow neck, and a smaller extracellular vestibule. Interestingly, the lining wall of the ion-permeation/phospholipid transport pathway may be formed, at least partially, by membrane phospholipids, though the degree of pore-wall forming by phospholipids likely varies among TMEM16 proteins. Thus, the biophysical properties and activation mechanisms of TMEM16 proteins could differ from each other accordingly. Here we review the current understanding of the structure and function of TMEM16 molecules.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces the externalization of phosphatidylserine, a key molecule in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove, and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements involved in groove opening and of how lipids reorganize outside the closed groove remain unknown. Using cryogenic electron microscopy, we directly visualize how lipids associate at the closed groove of Ca2+-bound nhTMEM16 in nanodiscs. Functional experiments pinpoint the lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryoEM structure determination.

Domain- and state-specific shape of the electric field tunes voltage sensing in voltage-gated sodium channels

The ability to sense transmembrane voltage underlies most physiological roles of voltage-gated sodium (Nav) channels. Whereas the key role of their voltage-sensing domains (VSDs) in channel activation is well established, the molecular underpinnings of voltage coupling remain incompletely understood. Voltage-dependent energetics of the activation process can be described in terms of the gating charge that is defined by coupling of charged residues to the external electric field. The shape of the electric field within VSDs is therefore crucial for the activation of voltage-gated ion channels. Here, we employed molecular dynamics simulations of cardiac Nav1.5 and bacterial NavAb, together with our recently developed tool g_elpot, to gain insights into the voltage-sensing mechanisms of Nav channels via high-resolution quantification of VSD electrostatics. In contrast to earlier low-resolution studies, we found that the electric field within VSDs of Nav channels has a complex isoform- and domain-specific shape, which prominently depends on the activation state of a VSD. Different VSDs vary not only in the length of the region where the electric field is focused but also differ in their overall electrostatics, with possible implications in the diverse ion selectivity of their gating pores. Due to state-dependent field reshaping, not only translocated basic but also relatively immobile acidic residues contribute significantly to the gating charge. In the case of NavAb, we found that the transition between structurally resolved activated and resting states results in a gating charge of 8e, which is noticeably lower than experimental estimates. Based on the analysis of VSD electrostatics in the two activation states, we propose that the VSD likely adopts a deeper resting state upon hyperpolarization. In conclusion, our results provide an atomic-level description of the gating charge, demonstrate diversity in VSD electrostatics, and reveal the importance of electric-field reshaping for voltage sensing in Nav channels.

Activation of TMEM16F by inner gate charged mutations and possible lipid/ion permeation mechanisms

Transmembrane protein 16F (TMEM16F) is a ubiquitously expressed Ca2+-activated phospholipid scramblase that also functions as a largely non-selective ion channel. Though recent structural studies have revealed the closed and intermediate conformations of mammalian TMEM16F (mTMEM16F), the open and conductive state remains elusive. Instead, it has been proposed that an open hydrophilic pathway may not be required for lipid scrambling. We previously identified an inner activation gate, consisting of F518, Y563, and I612, and showed that charged mutations of the inner gate residues led to constitutively active mTMEM16F scrambling. Herein, atomistic simulations show that lysine substitution of F518 and Y563 can indeed lead to spontaneous opening of the permeation pore in the Ca2+-bound state of mTMEM16F. Dilation of the pore exposes hydrophilic patches in the upper pore region, greatly increases the pore hydration level, and enables lipid scrambling. The putative open state of mTMEM16F resembles the active state of fungal scramblases and is a meta-stable state for the wild-type protein in the Ca2+-bound state. Therefore, mTMEM16F may be capable of supporting the canonical in-groove scrambling mechanism in addition to the out-of-groove one. Further analysis reveals that the in-groove phospholipid and ion transduction pathways of mTMEM16F overlap from the intracellular side up to the inner gate but diverge from each other with different exits to the extracellular side of membrane.

Evolutionary history of metazoan TMEM16 family

Most of Transmembrane protein 16 (TMEM16) proteins function as either a Ca²⁺-activated Cl^- channel (CaCC) or phospholipid scramblase (CaPLSase) and play diverse physiological roles. It is well conserved in eukaryotes; however, the origin and evolution of different subfamilies in Metazoa are not yet understood. To uncover the evolutionary history of the TMEM16 family, we analyzed 398 proteins from 74 invertebrate species using evolutionary genomics. We found that the TMEM16C-F and J subfamilies are vertebrate-specific, but the TMEM16A/B, G, H, and K subfamilies are ancient and present in many, but not all metazoan species. The most ancient subfamilies in Metazoa, TMEM16L and M, are only maintained in limited species. TMEM16N and O are Cnidaria- and Ecdysozoa-specific subfamilies, respectively, and Ctenophora, Xenacoelomorpha, and Rotifera contain species-specific proteins. We also identified TMEM16 genes that are closely linked together in the genome, suggesting that they have been generated via recent gene duplication. The anoctamin domain structures of invertebrate-specific TMEM16 proteins predicted by AlphaFold2 contain conserved Ca²⁺-binding motifs and permeation pathways with either narrow or wide inner gates. The inner gate distance of TMEM16 protein may have frequently switched during metazoan evolution, and thus determined the function of the protein as either CaCC or CaPLSase. These results demonstrate that TMEM16 family has evolved by gene gain and loss in metazoans, and the genes have been generally under purifying selection to maintain protein structures and physiological functions.

The permeation of potassium ions through the lipid scrambling path of the membrane protein nhTMEM16

The TMEM16 family of transmembrane proteins includes Ca2+-activated phospholipid scramblases (PLS) that can also function as non-selective ion channels. Extensive structural and functional studies have established that a membrane-exposed hydrophilic groove in TMEM16 PLS can serve as a translocation pathway for lipids. However, it is still unclear how the TMEM16 PLS conduct ions. A “protein-delimited pore” model suggests that ions are translocated through a narrow opening of the groove region, which is not sufficiently wide to allow lipid movement, whereas a “proteolipidic pore” model envisions ions and lipids translocating through an open conformation of the groove. We investigated the dynamic path of potassium ion (K+) translocation that occurs when an open groove state of nhTMEM16 is obtained from long atomistic molecular dynamics (MD) simulations, and calculated the free energy profile of the ion movement through the groove with umbrella sampling methodology. The free energy profile identifies effects of specific interactions along the K+ permeation path. The same calculations were performed to investigate ion permeation through a groove closed to lipid permeation in the nhTMEM16 L302A mutant which exhibits a stable conformation of the groove that does not permit lipid scrambling. Our results identify structural and energy parameters that enable K+ permeation, and suggest that the presence of lipids in the nhTMEM16 groove observed in the simulations during scrambling or in/out diffusion, affect the efficiency of K+ permeation to various extents.

TMEM16 scramblases thin the membrane to enable lipid scrambling

TMEM16 scramblases dissipate the plasma membrane lipid asymmetry to activate multiple eukaryotic cellular pathways. Scrambling was proposed to occur with lipid headgroups moving between leaflets through a membrane-spanning hydrophilic groove. Direct information on lipid-groove interactions is lacking. We report the 2.3 Å resolution cryogenic electron microscopy structure of the nanodisc-reconstituted Ca2+-bound afTMEM16 scramblase showing how rearrangement of individual lipids at the open pathway results in pronounced membrane thinning. Only the groove's intracellular vestibule contacts lipids, and mutagenesis suggests scrambling does not require specific protein-lipid interactions with the extracellular vestibule. We find scrambling can occur outside a closed groove in thinner membranes and is inhibited in thicker membranes, despite an open pathway. Our results show afTMEM16 thins the membrane to enable scrambling and that an open hydrophilic pathway is not a structural requirement to allow rapid transbilayer movement of lipids. This mechanism could be extended to other scramblases lacking a hydrophilic groove.

Polymodal Control of TMEM16x Channels and Scramblases

The TMEM16A/anoctamin-1 calcium-activated chloride channel (CaCC) contributes to a range of vital functions, such as the control of vascular tone and epithelial ion transport. The channel is a founding member of a family of 10 proteins (TMEM16x) with varied functions; some members (i.e., TMEM16A and TMEM16B) serve as CaCCs, while others are lipid scramblases, combine channel and scramblase function, or perform additional cellular roles. TMEM16x proteins are typically activated by agonist-induced Ca2+ release evoked by Gq-protein-coupled receptor (GqPCR) activation; thus, TMEM16x proteins link Ca2+-signalling with cell electrical activity and/or lipid transport. Recent studies demonstrate that a range of other cellular factors—including plasmalemmal lipids, pH, hypoxia, ATP and auxiliary proteins—also control the activity of the TMEM16A channel and its paralogues, suggesting that the TMEM16x proteins are effectively polymodal sensors of cellular homeostasis. Here, we review the molecular pathophysiology, structural biology, and mechanisms of regulation of TMEM16x proteins by multiple cellular factors.

Gating and Regulatory Mechanisms of TMEM16 Ion Channels and Scramblases

The transmembrane protein 16 (TMEM16) family consists of Ca2+-activated ion channels and Ca2+-activated phospholipid scramblases (CaPLSases) that passively flip-flop phospholipids between the two leaflets of the membrane bilayer. Owing to their diverse functions, TMEM16 proteins have been implicated in various human diseases, including asthma, cancer, bleeding disorders, muscular dystrophy, arthritis, epilepsy, dystonia, ataxia, and viral infection. To understand TMEM16 proteins in health and disease, it is critical to decipher their molecular mechanisms of activation gating and regulation. Structural, biophysical, and computational characterizations over the past decade have greatly advanced the molecular understanding of TMEM16 proteins. In this review, we summarize major structural features of the TMEM16 proteins with a focus on regulatory mechanisms and gating.

g_elpot: A Tool for Quantifying Biomolecular Electrostatics from Molecular Dynamics Trajectories

Electrostatic forces drive a wide variety of biomolecular processes by defining the energetics of the interaction between biomolecules and charged substances. Molecular dynamics (MD) simulations provide trajectories that contain ensembles of structural configurations sampled by biomolecules and their environment. Although this information can be used for high-resolution characterization of biomolecular electrostatics, it has not yet been possible to calculate electrostatic potentials from MD trajectories in a way allowing for quantitative connection to energetics. Here, we present g_elpot, a GROMACS-based tool that utilizes the smooth particle mesh Ewald method to quantify the electrostatics of biomolecules by calculating potential within water molecules that are explicitly present in biomolecular MD simulations. g_elpot can extract the global distribution of the electrostatic potential from MD trajectories and measure its time course in functionally important regions of a biomolecule. To demonstrate that g_elpot can be used to gain biophysical insights into various biomolecular processes, we applied the tool to MD trajectories of the P2X3 receptor, TMEM16 lipid scramblases, the secondary-active transporter Glt_Ph, and DNA complexed with cationic polymers. Our results indicate that g_elpot is well suited for quantifying electrostatics in biomolecular systems to provide a deeper understanding of its role in biomolecular processes.