A Generic Numbering Scheme for TMEM16 Scramblases

The TMEM16 family of calcium-activated phospholipid scramblases (CaPLSs) and chloride channels (CaCCs) performs diverse physiological functions that include regulation of blood coagulation and apoptotic signaling, through a shared ten-transmembrane-helix (TM) architecture organized around a hydrophilic lipid-translocating groove. Mechanistic studies of TMEM16 family members have been hampered by the absence of a unified positional reference framework that would permit direct comparison of structurally equivalent residues across paralogs with different sequence numbering systems. Here we introduce a generic numbering scheme for TMEM16 scramblases (GNS-TMEM16), modeled on the Ballesteros & Weinstein system established for class A G protein-coupled receptors. A reference alignment (TMEM16-RA) was constructed from twelve human and mouse TMEM16 scramblases (TMEM16C/D/E/F/G/J) using structure-based ClustalW alignment of the ten TM helices. From this alignment, a TM-specific reference residue (TsRR) was identified for each helix by hierarchical application of three criteria: (1) 100% conservation in the core TMEM16-RA; (2) conservation in an augmented reference alignment (TMEM16-ARA) incorporating a group of phylogenetically more distant homologs composed of nhTMEM16, afTMEM16, TMEM16K, TMEM16A, and TMEM16B; and (3) structural and functional considerations, including helix-perturbing character, groove localization, conserved motif membership, and central TM position. The resulting ten TsRRs are Y1.50, W2.50, R3.50, E4.50, F5.50, P6.50, E7.50, D8.50, W9.50, and E10.50, and are illustrated in mTMEM16F. Each residue is assigned the identifier N.m(k), where N is the TM number, m is the position relative to the TsRR (for which m = 50), and k is the absolute sequence number. Loop residues receive dual identifiers referenced to the TsRRs of both flanking helices. Application of the GNS-TMEM16 is illustrated with the comparisons of the groove-opening measurements using pairwise distances between residues identified by their N.m indices to be corresponding across mTMEM16F, afTMEM16, and nhTMEM16. The results bring to light the advantages of corresponding residues identification in different TMEM16 proteins and show that the mammalian scramblase undergoes substantially larger separation at the extracellular groove entrance than either fungal homolog. Comparison of mutagenesis data guided by N.m correspondence shows at the conserved (E3.55,R6.26) salt-bridge locus, Ala substitution reduces activity more than 100-fold in nhTMEM16 but less than 2-fold in afTMEM16, illustrating that the GNS identifies structural equivalence of position without implying functional equivalence of the residue, which is a distinct advantage of GNS in providing mechanistic interpretation across paralogs. Also described is a protocol for extending the GNS-TMEM16 to uncharacterized protein sequences, including AlphaFold-predicted models, using structural superposition to mTMEM16F. Thus, the presented GNS-TMEM16 provides a stable positional reference for the integration and comparative analysis of structural, computational, and functional data across the TMEM16 family, utilizing a construction strategy applicable to yet other polytopic membrane protein families sharing a common transmembrane fold.