
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

70652
10.1038/s41598-024-70652-8
Article
Mint/X11 PDZ domains from non-bilaterian animals recognize and bind CaV2 calcium channel C-termini in vitro
Harracksingh Alicia N. 1
Singh Anhadvir 1
Mayorova Tatiana D. 12
Bejoy Brian 1
Hornbeck Jillian 3
Elkhatib Wassim 1
McEdwards Gregor 1
Gauberg Julia 1
Taha Abdul 1
Islam Ishrat Maliha 1
Erclik Ted 1
Currie Mark A. 1
Noyes Marcus 3
Senatore Adriano adriano.senatore@utoronto.ca

1
1 https://ror.org/03dbr7087 grid.17063.33 0000 0001 2157 2938 Department of Biology, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, ON L5L 1C6 Canada
2 https://ror.org/01cwqze88 grid.94365.3d 0000 0001 2297 5165 NINDS, National Institutes of Health, Bethesda, MD 20892 USA
3 grid.137628.9 0000 0004 1936 8753 Institute for Systems Genetics, NYU Grossman School of Medicine, 550 1st Ave, New York, NY 10016 USA
16 9 2024
16 9 2024
2024
14 2161527 3 2024
20 8 2024
© Crown 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
PDZ domain mediated interactions with voltage-gated calcium (CaV) channel C-termini play important roles in localizing membrane Ca2+ signaling. The first such interaction was described between the scaffolding protein Mint-1 and CaV2.2 in mammals. In this study, we show through various in silico analyses that Mint is an animal-specific gene with a highly divergent N-terminus but a strongly conserved C-terminus comprised of a phosphotyrosine binding domain, two tandem PDZ domains (PDZ-1 and PDZ-2), and a C-terminal auto-inhibitory element that binds and inhibits PDZ-1. In addition to CaV2 chanels, most genes that interact with Mint are also deeply conserved including amyloid precursor proteins, presenilins, neurexin, and CASK and Veli which form a tripartite complex with Mint in bilaterians. Through yeast and bacterial 2-hybrid experiments, we show that Mint and CaV2 channels from cnidarians and placozoans interact in vitro, and in situ hybridization revealed co-expression in dissociated neurons from the cnidarian Nematostella vectensis. Unexpectedly, the Mint orthologue from the ctenophore Hormiphora californiensis strongly bound the divergent C-terminal ligands of cnidarian and placozoan CaV2 channels, despite neither the ctenophore Mint, nor the placozoan and cnidarian orthologues, binding the ctenophore CaV2 channel C-terminus. Altogether, our analyses suggest that the capacity of Mint to bind CaV2 channels predates bilaterian animals, and that evolutionary changes in CaV2 channel C-terminal sequences resulted in altered binding modalities with Mint.

Subject terms

Protein-protein interaction networks
Molecular evolution
Phylogenetics
Ontario Graduate ScholarshipNSERC Canadian Graduate ScholarshipNSERC DiscoveryRGPIN-2023-05615 RGPIN-2021-03557 Currie Mark A. Senatore Adriano Canadian Foundation for Innovation Grant40684 35297 Currie Mark A. Senatore Adriano NIH GrantR01GM133936 Noyes Marcus NSERC Discovery Accelerator SupplementRGPAS-2021-00002 Senatore Adriano Ontario Early Researcher AwardER17-13-247 Senatore Adriano issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

CaV1 and CaV2 voltage-gated calcium (CaV) channels are high voltage-activated channels that couple membrane excitation with numerous important Ca2+-dependent processes including neurotransmitter release in neurons and contraction in muscle1. Phylogenetic analyses suggest that these two channel types evolved through gene duplication and diversification very early during animal evolution2–4, and accordingly, they share common structural features and biophysical properties that distinguish them from the more phylogenetically distant low voltage-activated CaV3 calcium channels1,5. Despite their similarities, CaV1 and CaV2 channels also exhibit some unique and mostly non-overlapping cellular specializations that are broadly conserved in animals. For example, conserved between mammals and invertebrates including arthropods (i.e., the fly Drosophila melanogaster and the marine isopod Idotea baltica), platyhelminthes (i.e., the flatworm Dugesia japonica), nematodes (the worm Caenorhabditis elegans), and molluscs (the snail Lymnaea stagnalis and the sea hare Aplysia californica), CaV1 calcium channels are the primary drivers of cardiac and body muscle contraction6–13, while CaV2 channels are the primary drivers of neurotransmitter release at fast synapses14–19. In part, these differences in cellular function might be supported by evolved differences in intrinsic ion-conducting properties. For example, while both CaV1 and CaV2 channels are subject to negative feedback regulation by the Ca2+ sensor protein calmodulin, CaV1 channels exhibit an additional and unique form of regulation by the C-terminal lobe of calmodulin that is insensitive to cytoplasmic Ca2+ buffering and accelerates inactivation of the macroscopic current20–22. This feature, combined with the generally lower sensitivity of CaV1 channels to voltage-dependent inactivation23,24, makes them well suited for roles in cells that exhibit prolonged bouts of depolarization and rises in cytoplasmic [Ca2+], as occurs in contracting muscle cells. Mammalian CaV2 channels on the other hand are uniquely subject to negative feedback regulation by G proteins (i.e., Gβγ heterodimers) which become activated by neuromodulatory G protein-coupled receptors25–28, an important neuromodulatory process that can fine tune synaptic transmission and information flow along neural circuits. Gβγ regulation is also apparent for the CaV2 channel from the mollusc L. stagnalis27–29, but not from the early-diverging animal Trichoplax adhaerens30.

With Ca2+ ions generally kept at very low levels in the cytoplasm (i.e., roughly 100 nM)31, CaV channel activation results in a discrete and transient domain of cytoplasmic Ca2+, with a > 1000-fold increase in concentration from baseline near the channel pore that sharply declines over a distance of 100 nm32,33. A key mechanism for proper activation of Ca2+-dependent processes by CaV channels thus involves their close apposition to calcium-responsive machinery, such that target Ca2+-binding sites can become adequately saturated34,35. Interestingly, various protein–protein interactions that distinguish CaV1 and CaV2 channels appear to be broadly conserved, and hence may represent an additional diversification that contributed to the evolution of these two channels types4,36,37. For example, in mouse and C. elegans, CaV1 channels associate with the postsynaptic scaffolding protein Shank, a conserved interaction that permits the calcium channels to activate changes in nuclear gene expression through the transcription factor CREB38–40. Here, mammalian CaV1.3 channels and C. elegans CaV1 have respective C-terminal short linear motifs (SLiMs) of ITTLCOOH and VTTLCOOH that are recognized and bound by the PDZ domain of Shank41, a protein–protein interaction domain named for its shared presence in three proteins: Postsynaptic density protein 95 kDa (PSD95), Drosophila disc large tumor suppressor 1 (Dlg1), and Zonula occludens-1 protein (Zo-1). PDZ domains are commonly found in synaptic scaffolding proteins, where they differentially recognize and bind the extreme C-termini of target proteins via 16 unique ligand specificity classes which are conserved between vertebrates and invertebrates41.

Importantly, like CaV1 channels, CaV2 channels from vertebrates and invertebrates also exhibit conserved interactions with RIM42–44 and Mint (a.k.a. X11)15,45–49, both via PDZ domains that recognize the extreme C-terminus of CaV2 channels at a conserved consensus ligand sequence of DDWCCOOH. Indeed, the first report of a PDZ domain binding the C-terminus of any CaV channel was from Mint46, a neuronal modular adaptor protein better known for its interaction with amyloid beta precursor protein (APP)50–53 and its γ-secretase processing enzyme presenilin54,55, and hence possible involvement in Alzheimer’s disease56. Here, a yeast 2-hybrid assay revealed an interaction between the first of two tandem PDZ domains of the rat Mint-1 protein (i.e., PDZ-1) and the extreme C-terminus of the rat CaV2.2 channel bearing a PDZ ligand sequence of DHWCCOOH46. Subsequently, the in vitro interaction between Mint and CaV2 channels from mammals was corroborated though a series of experiments including yeast 2-hybrid, protein pulldown, and nuclear magnetic resonance (NMR)45,47,48. Furthermore, a direct in vitro interaction between Mint and CaV2 was demonstrated for L. stagnalis via yeast 2-hybrid15, and indirectly for chick (Gallus gallus) where endogenous Mint protein was pulled down from brain lysates using the CaV2.2 channel C-terminus as bait49. Notably, while the physiological function of the RIM-CaV2 channel interaction has been established as important for the proper presynaptic localization of CaV2 channels in mammals43, C. elegans44, D. melanogaster42, and A. californica16, the physiological significance of the Mint-CaV2 channel interaction remains unclear. In mice, double knockout of the two neuronal Mint genes, Mint-1 and Mint-2, did not disrupt the presynaptic localization of CaV2 channels, despite causing significant disruptions in synaptic morphology and function, and lethality in 80% of newborns57. Similarly, genetic ablation of Mint in C. elegans caused significant phenotypic defects, including disrupted vulva development and nervous system function, but did not perturb the presynaptic localization of CaV2 channels18. Lastly, siRNA downregulation of Mint in L. stagnalis did not perturb synaptic transmission in synapses reconstituted in vitro mediated by the CaV2 channel15, altogether suggesting that Mint does not play a role in localizing CaV2 channels to the synapse active zone. Nonetheless, the robust nature of this interaction in vitro and its conservation between vertebrates and invertebrates suggests a yet undiscovered neural function.

Spurred by recent observations that many synaptic genes evolved either prior to animals, or are present in early-diverging animals that lack a nervous system and synapses, we set out to explore the phylogenetic origins of Mint and several of its known interacting proteins (i.e., CASK, Veli, Amyloid-β precursor protein, presenilin, and neurexin), and to directly test the in vitro interaction between Mint and CaV2 channels from non-bilaterian animals of the phyla Cnidaria (jellyfish, sea anemones, corals), Placozoa (flat aneural animals), and Ctenophora (comb jellies). We demonstrate that the capacity for Mint and CaV2 channels to interact in vitro evolved prior to the emergence of bilaterian animals, and that Mint PDZ domains from non-bilaterian animals recognize and bind canonical CaV2 channel C-terminal DDWC-like ligands, as well as the distinct ligand of CTAVCOOH present in placozoan CaV2 channels. Furthermore, we provide evidence for this dual binding capacity of Mint which involves distinct structural elements within its C-terminus, comprised of two tandem PDZ domains (PDZ-1 and PDZ-2) and a C-terminal sequence that binds the PDZ-1 ligand-binding groove to impose autoinhibition of ligand binding48 (hereafter referred to as the PDZ-1 Binding Motif or P1BM).

Results

Phylogenetic analyses point to deep conservation of the Mint PDZ-1/CaV2 calcium channel interaction

CaV2 channels are unique to animals, having emerged alongside CaV1 channels via an ancient gene duplication event2,4. A protein phylogeny of select CaV2 channel sequences produced a branching pattern mostly consistent with the expected species phylogeny rooted on the most divergent orthologues from the ctenophore species Mnemiopsis leidyi and Hormiphora californiensis (Fig. 1A). The three vertebrate CaV2 channel isotypes, P/Q-type (CaV2.1), N-type (CaV2.2), and R-type (CaV2.3) form a strongly supported monophyletic clade reflecting their proposed emergence through gene/genome duplication events that occurred within the phylum Chordata, after tunicates (e.g., Ciona intestinalis and Styela clava) and cephalochordates (Branchiostoma belcheri) diverged. Consistent with previous reports, the well-supported monophyletic clades of CaV2 channels from cnidarians and platyhelminthes reflect independent duplication events that likely gave rise to CaV2a, CaV2b, and CaV2c in cnidarians, and CaV2A and CaV2B in platyhelminthes2,4,36.Fig. 1 The structural elements required for the Mint-CaV2 interaction are conserved in cnidarians and bilaterians. (A) Maximum likelihood protein phylogeny of CaV2 channels from selected animals, with names color-coded according to phylum. Sequences of the extreme C-termini of the channels are shown on the right of the phylogenetic tree, revealing a conserved PDZ ligand motif with a consensus sequence of DEDDWCCOOH. Species for which a direct interaction between CaV2 channels and Mint have been experimentally determined are indicated in red text to the right of select CaV2 channel C-terminal sequences (coIP: co-immunoprecipitation; Y2H: Yeast 2-Hybrid; NMR: Nuclear Magnetic Resonance spectroscopy). (B) Illustration of the CaV2 channel transmembrane topology, depicting four homologous repeat domains (DI-DIV) each with 6 transmembrane alpha helices (in black), as well as the intracellular N- and C-termini and cytoplasmic linkers that separate the four domains (in red). (C) EMBOSS Plotcon alignment similarity plot for single CaV2 channel representatives selected from each phylum from the tree depicted in panel A, revealing extensive divergence in amino acid sequences within the cytoplasmic regions (linkers and N- and C-termini), but strong conservation within transmembrane regions of Domains I-IV. Notably, the conserved DDWC ligand sequence appears as a short linear motif at the extreme end of the highly divergent and structurally disordered C-terminus. (D) Maximum likelihood phylogeny of Mint proteins from select animals. The lengths of the various Mint proteins are indicated with black bars on the left of the phylogenetic tree. The symbols on the black bars denote ubiquitously conserved PTB and tandem PDZ domains (i.e., PDZ-1 and PDZ-2), as well as unique tetratricopeptide (TCP) repeat domains in the placozoan orthologues77, as predicted with InterProScan. Also to the left of the tree is an alignment of the extreme C-terminal sequences of the various Mint proteins, revealing deep conservation of the auto-inhibitory domain (P1BM) that associates with the ligand binding pocket of PDZ-1 causing ligand binding auto-inhibition48. For panels A and D, node support values generated with SH-aLRT, UFBoot2, and aBayes are indicated by symbols according to the provided legend. (E) Illustration of the various domains found within the Mint protein that are involved in interactions with other proteins (MID: Munc-18 Interacting Domain; CID: CASK Interacting Domain; PTB: PhosphoTyrosine-binding Domain; PDZ domains 1 and 2: PDZ1 and PDZ2; P1BM: PDZ-1 auto-inhibitory peptide). The black arrows indicate domain-ligand relationships, with the head pointing towards the ligand. (F) EMBOSS Plotcon alignment similarity plot for single Mint representatives selected from each phylum from the tree depicted in panel B, revealing a highly divergent N-terminus, but strong sequence conservation within the C-terminal PTB domain, the tandem PDZ domains, and the P1BM auto-inhibitory domain.

Alignment of the extreme C-terminal protein sequences of CaV2 channels reveals a conserved amino acid motif with a consensus sequence of DDWCCOOH, shared between cnidarian and bilaterian homologues (Fig. 1A). This PDZ ligand motif can be classified as a SLiM, preceded by a long stretch of highly divergent amino acids that contribute to a mostly disordered CaV2 channel C-terminus (Fig. 1B and C). The most invariable amino acid within this motif is the tryptophan residue (W) in the second last position, a common feature of PDZ ligands being present in 5 of the 16 defined specificity classes41. Nonetheless, the presence of a distal C-terminal cysteine residue indicates that this motif is a non-canonical PDZ domain ligand. Apparent is the independent loss of DDWC-like motifs in several lineages of CaV2 channels. Within chordates, the chick (Gallus gallus) CaV2.1 channel and the C. intestinalis CaV2 channel lost this motif, the former attributed to the loss of a C-terminal exon that encodes the distal portion of the CaV2 channel C-terminus in chick and perhaps all birds4,58. Notably, vertebrate CaV2.3 channels bear cationic lysine/arginine (K/R) residues rather than tryptophan in the second last position, which likely disrupts PDZ domain binding since rat Mint-1 does not bind the human CaV2.3 C-terminus in vitro46, and mutation of the L. stagnalis CaV2 channel DDWC motif to DDKC disrupts its interaction with Mint15. Given the conservation of DDWC-like C-terminal SLiMs in cnidarian and bilaterian CaV2 channels, it appears that platyhelminthes and the nematode species Caenorhabditis elegans and Strongyloides ratti independently lost the DDWCCOOH motif. However, since the CaV2 channel from the early-diverging nematode Trichinella pseudospiralis does bear a DDWC-like SLiM (EDWC), loss of this motif only occurred in a subset of nematode species. Lastly, we note a complete absence of DDWC-like motifs in CaV2 channel C-termini from placozoans and ctenophores, the former lacking a nervous system and synapses59, and the latter comprising the most divergent animal phylum to possess CaV2 channels, and which interestingly, are proposed to have independently evolved a nervous system and synapses60.

Phylogeny of Mint and several key interacting proteins

To better understand the phylogenetic properties of Mint, we mined select genomes and transcriptomes to generate a maximum likelihood phylogeny which was rooted by the ctenophore orthologues (Fig. 1D). Given the absence of CaV2 channels in poriferans (sponges)2,4, we excluded them from our analysis despite identifying bona fide Mint homologues for several species including Amphimedon queenslandica. Like CaV2 channels, most invertebrates included in our analysis appear to possess only single Mint orthologues, with independent duplications apparent for vertebrates (i.e., Mint 1 to 3) and D. melanogaster (Mint A and B). Analysis of the identified Mint protein sequences with InterProScan61 reveals complete conservation of a predicted Phosphotyrosine-Binding (PTB) domain, followed by two tandem PDZ domains (PDZ-1 and PDZ-2) near the C-terminal end of the protein (Fig. 1D and E). Alignment of the extreme C-termini of the analyzed Mint sequences reveals strong sequence conservation (Fig. 1D), notable because this sequence serves an auto-inhibitory function on the Mint protein by binding to the PDZ-1 ligand recognition groove and thus inhibiting ligand binding48. In contrast to the conserved C-terminus, the N-terminus of Mint is quite variable in length (Fig. 1D) and amino acid sequence, evident in a plot of global sequence conservation among single representative Mint orthologues from each phylum included in our analysis (Fig. 1F). This strong sequence divergence is notable because the N-terminus of Mint bears the CASK Interaction Domain (CID), conserved between vertebrates (Mint-1 only), D. melanogaster, and C. elegans62, which contributes to the formation of a tripartite complex with the proteins CASK and Veli63–65. The N-terminus of Mint also harbours the Munc18 Interaction Domain (MID), a protein–protein interaction that associates Mint function with vesicle exocytosis and hormone secretion66. Altogether, the fact that our alignments failed to identify conserved CID and MID signatures suggests these elements are either highly divergent in protein sequence, or are absent in early-diverging, non-bilaterian Mint orthologues. Nonetheless, we note conservation of both CASK and Veli in all metazoans included in our study (Fig. 2A and B), most bearing canonical predicted domain architectures including a Calmodulin Kinase (CaMK) domain at the N-terminus of CASK shown in crystal structures to associate with the Mint CID motif (Protein Data Bank[PDB] accession number 6LNM), as well as Lin-2/Lin-7 (L27) domains, important for the interaction between CASK and Veli64.Fig. 2 Phylogenetic properties of several key proteins that interact with Mint. (A) Maximum likelihood protein phylogeny of CASK homologues from various animals. Protein lengths are indicated with black bars, and domains predicted with InterProScan are illustrated (CaMK: Ca2+/calmodulin-dependent protein kinase domain; L27: Lin-2/Lin-7 domain; PDZ: PSD95/Dlg1/zo-1 domain; SH3: Src Homology 3 domain; GK: guanylate kinase domain). (B) Maximum likelihood protein phylogeny of Veli homologues from various animals. Protein lengths are indicated with black bars, and domains predicted with InterProScan are illustrated (L27: Lin-2/Lin-7 domain; PDZ: PSD95/Dlg1/zo-1 domain). (C) Maximum likelihood protein phylogeny of amyloid β precursor proteins (APP) and amyloid β-like (APPL) homologues from various animals, colour-coded according to phyla. C-terminal sequences are shown on the right, illustrating deep conservation of the YENPTY motif which is bound to by the Mint PTB domain. No ctenophore homologues were identified, and the placozoan and tardigrade sequences lack the canonical YENPTY motif. The lengths of the APP and APPL proteins are indicated with black bars, and protein domains predicted with InterProScan are illustrated (LRR: leucine-rich repeat domain; SH2: Src Homology 2 domain; AgGEx: amyloidogenic glycoprotein, extracellular domain; Kun: pancreatic trypsin inhibitor Kunitz domain; EGFL: EGF-like domain; AgGE2: amyloidogenic glycoprotein, E2 domain; TM: transmembrane helix; AβP: amyloidogenic glycoprotein, amyloid-beta peptide; CTD: beta-amyloid precursor protein C-terminal). (D) Maximum likelihood protein phylogeny of presenilin homologues from various animals, colour-coded according to phyla (legend below). An alignment of extreme C-terminal protein sequences is shown on the right, illustrating deep conservation of the ligand sequence for the Mint PDZ-1 domain. The lengths of the presenilin proteins are indicated with black bars, and transmembrane helices (TMs) predicted with Phobius are illustrated. For all panels, node support values generated with SH-aLRT, UFBoot2, and aBayes are indicated by symbols according to the provided legend.

The deep conservation of the Mint PTB, tandem PDZ, and C-terminal auto-inhibitory PDZ-1 Binding Motif (P1BM), suggests that the various interactions mediated by these domains might also be conserved. The PTB domain, for example, mediates protein–protein interactions with several partners including Amyloid Precursor Protein (APP)50–53, Amyloid Precursor Protein-Like (APPL)67, and the Golgi transport protein Rab668 (Fig. 1E). We therefore sought to survey the phylogeny of some of these proteins, with a primary focus on exploring conserved motifs associated with Mint interactions, and predicted domain architectures. Interestingly, our analysis revealed that APP and APPL homologues are broadly conserved in animals (Fig. 2C), finding homologues in all examined species with the exception of ctenophores, most bearing a YENPTY motif that forms a beta strand which integrates into a binding pocket of the Mint PTB domain50,69. We only found single APP homologues in invertebrates, but a strongly supported monophyletic clade of vertebrate homologues consistent with gene duplications giving rise to APP, AAPL1, and APPL2 (Fig. 2C). As noted, the first PDZ domain of Mint (PDZ-1), recognizes and binds the CaV2 channel DDWC motif15,46,47,70. Mint PDZ-1, as well as PDZ-2, also recognize the C-terminus of the transmembrane protein presenilin48,54, 55 which is part of the γ-secretase complex that regulates target proteins including APP via proteolytic processing71, and PDZ-1 binds the trans-synaptic cell adhesion protein neurexin45,47. Phylogenetically, all examined species possess presenilin homologues with vertebrates uniquely duplicating the gene to give rise to presenilins 1 and 2 (Fig. 2D). Additionally, all of the identified presenilin homologues possess 8 to 10 predicted transmembrane helices, as well as a conserved hydrophobic C-terminus with a consensus sequence of DSLASHQVYICOOH (Fig. 2D), reflecting deep conservation of a putative Mint PDZ-1/PDZ-2 ligand sequence as defined for mammalian Mint 1 and 2 in vitro48,54,55. We also conducted a phylogenetic analysis of synaptic neurexin homologues (referred to as neurexins 1 to 3 in vertebrates and neurexin-1 in invertebrates, collectively referred to as neurexin-1 genes in this study), which belong to a larger group of cell adhesion proteins collectively referred to as the neurexin family72 (Fig. 3). For our analysis, we focused on homologues bearing domain signatures of two LamG2 domains flanked by EGF domains, which are unique to neurexin-1. We also included one of each of the three types of neurexin-1 homologues previously identified from N. vectensis (i.e., neurexins alpha 1, delta 1, and epsilon 1), a species recently shown to possess at least 10 neurexin homologues73. Like the other analyzed genes, the phylogenetic analysis suggests that neurexin-1 duplicated in vertebrates to produce neurexins 1 to 3, and separately in spiralian species within the phyla Mollusca, Brachiopoda, and Platyhelminthes, the chelicerate arthropod Limulus polyphemus, and both placozoans and cnidarians. In molluscs and brachiopods, one of the duplicated neurexin-1 paralogues retained a canonical predicted domain architecture comprised of repeating EGF-LamG2-LamG2 domains and a C-terminal transmembrane helix, while the other incorporated a unique C-terminal WD40 domain and an extended arrangement of LamG2 domains near the N-terminus. This derived isotype from the brachiopod species Lingula anatina also possesses a unique set of 9 predicted tandem fibronectin type 3 domains that are associated with cell adhesion and migration. A similar phenomenon is apparent for the duplicated neurexin-1 paralogues from platyhelminths, with one isotype retaining a more canonical domain arrangement, and the other being considerably shorter in length with less predicted LamG2 and EGF domains. With respect to putative C-terminal PDZ ligands, it is notable that most of the examined neurexin homologues exhibit conserved C-terminal KE[W/Y]YVCOOH motifs (Fig. 3), shown to associate with the PDZ-1 domain of rat Mint-1 via surface plasmon resonance47.Fig. 3 Maximum likelihood protein phylogeny of neurexin-1 homologues from various representative animals. Names of species and genes are colour-coded according to phylum. To the right of the phylogenetic tree is an alignment of the extreme C-terminal sequences of the proteins, illustrating deep conservation of the Mint PDZ-1 ligand motif. The lengths of the neurexin-1 homologues are indicated with black bars, and protein domains predicted with InterProScan are illustrated (LNS: laminin, neurexin, sex-hormone binding globulin domain; EGFL: epidermal growth factor-like domain; FN3: fibronectin type III domain; WD40: tryptophan-aspartic acid 40 repeat domain). Node support values generated with SH-aLRT, UFBoot2, and aBayes are indicated by symbols according to the provided legend.

Like PDZ-1, the Mint PDZ-2 domain mediates several distinct interactions from PDZ-1, including with the copper cofactor to the antioxidant superoxide dismutase-1 (CCS)74 and the Ca2+-regulated transcription factor NK-κB75 (Fig. 1E). Lastly, we note unique tetra-trico-peptide repeat (TCP) protein–protein interaction domains upstream of the PTB domain in the two examined placozoan Mint orthologues, notable because placozoans exhibit an enrichment of this particular domain within their proteomes76, and like PDZ domains, some TPR domains can recognize C-terminal ligands of target proteins77.

Stronger sequence divergence of Mint PDZ-1 versus PDZ-2 suggests placozoans and ctenophores lack the capacity for a Mint-CaV2 interaction

To better understand the conservation between orthologous Mint PDZ domains, we generated protein sequence alignments and predicted secondary structures of select PDZ-1 and PDZ-2 domains (Fig. 4A and B). Both PDZ-1 and PDZ-2 bear a standard arrangement of beta strands (β1 to β6) and alpha helices (α1 and α2), with the exception of a short β5 strand, suggesting both adopt a mostly canonical PDZ domain tertiary structure78. For PDZ-1, we were able to map the positions of six amino acids whose sequence divergence correlates with evolutionary divergence in ligand specificity between different classes of PDZ domains79. These residues, which are clustered within and proximal to the β2 and β3 strands, are highly conserved within the PDZ-1 domains of bilaterian Mint orthologues, with some notable differences apparent for those from the non-bilaterian species T. adhaerens and H. californiensis (Fig. 4A). More variable are five amino acids in the α2 helix that for the rat Mint-1 PDZ-1 domain were shown to exhibit strong spectral shifts in NMR spectrograms upon binding a CaV2.2 channel DDWC ligand48, with the exception of a cationic lysine/arginine residue in the tail end of α2 only absent in the T. adhaerens and H. californiensis orthologues.Fig. 4 The Mint PDZ-1 domain is more divergent in protein sequence compared to PDZ-2. (A) Protein sequence alignment of select Mint PDZ-1 domains revealing deep sequence conservation throughout the secondary structure comprised of 6 β-strands and 3 α-helices. The black chevrons under the alignment indicate the position of the 17 amino acids correlated with ligand specificity41, the red asterisks indicate positions that correlate with the evolution of distinct ligand preferences for PDZ domains in animals79, and the blue asterisks denote amino acid positions that show major spectral shifts upon binding to the CaV2 channel DDWC ligand in NMR studies48. (B) Protein sequence alignment of select Mint PDZ-2 domains. Like panel A, the black chevrons under the alignment indicate the position of the 17 amino acids correlated with ligand specificity41. For panels A and B, the black bars indicate the position of the ligand carboxylate accommodating loop within the PDZ domain ligand-binding groove, the latter comprised of the β2 strand and the α2 helix. The labelled arrows placed under the alignments denote expected locations of β strands, and the labelled bars the locations of alpha helices. (C) Sequence identity matrix (percent identity) of 17 aligned amino acids within the PDZ domain ligand-binding groove shown to correlate with ligand specificity41. Font colors for names of species/genes indicate of phylum according to the legend in Fig. 1A. Percent identity values that fall below 70%, described as the threshold for correlated ligand specificity, are depicted in red font. (D) Similar sequence identity matrix as described in panel C, but for the Mint PDZ-2 domain.

Together, the noted sequences above overlap with 17 amino acid positions that correlate with distinct ligand specificities of PDZ domains determined via high-throughput peptide phage display screening41. Specifically, 17 amino acids clustered within the β2 and β3 strands and the α2 helix (Fig. 4C), which mediate key contacts with C-terminal peptide ligands80, were shown to strongly predict common ligand-binding specificity when conserved above 70% sequence identity41. A sequence identity matrix of these 17 aligned amino acids reveals stronger conservation in Mint PDZ-2 compared to PDZ-1 (Fig. 4C and D). For PDZ-1, notable divergence is evident for rat Mint-3 and the Mint orthologues from T. adhaerens and H. californiensis, with all comparisons yielding alignment scores below 70%. Sequence comparisons of PDZ-1 show an overall reduction in alignment scores compared to PDZ-2, particularly so for rat Mint-3 and orthologues from Strongylocentrotus purpuratus, C. elegans, N. vectensis, and H. californiensis. The non-bilaterian Mint PDZ-1 domains are highly divergent, both from each other and from bilaterian orthologues, producing alignment scores ranging between 64.7% and 17.7%. Notably, a correlation between global sequence identity of entire PDZ domains and ligand specificity was also reported, with a sequence identity cut-off of 50%41. Alignment matrices of entire Mint PDZ-1 and PDZ-2 domains reveal that most comparisons yield sequence identity scores above 50%, with only the PDZ-1 domains from the T. adhaerens and H. californiensis orthologues consistently falling below the 50% cut-off (Figure S1). Altogether, the examined sequence features suggest that Mint PDZ-1 ligand binding is more divergent compared to PDZ-2, particularly so for the Mint orthologues from placozoans and ctenophores which as noted above have CaV2 channels lacking DDWC-like motifs (Fig. 1A).

To date, the structural arrangement of the two tandem PDZ domains of Mint have not been experimentally determined. However, separate NMR structures of the human Mint-1 PDZ-1 and PDZ-2 domains are available (PDB numbers 1U37 and 1Y7N, respectively)48,81. We therefore used AlphaFold282 to predict structures of entire PDZ-1/PDZ-2/P1BM modules of Mint orthologues from R. norvegicus (Mint-1), N. vectensis, T. adhaerens, and H. californiensis (Fig. 5A). All four structures were predicted with high to very high confidence, in particular the H. californiensis orthologue (Figure S2). Structural alignment between the solved Mint NMR structures and the predicted PDZ-1/PDZ-2/P1BM structure of R. norvegicus Mint-1 reveals root mean square deviation (r.m.s.d.) values of 0.87 and 0.67 angstroms for PDZ-1 and PDZ-2, respectively, indicating low deviation in spatial positioning of backbone alpha carbons between the solved and predicted structures. In all predicted structures, the extreme C-terminus of P1BM is found nestled within the ligand-binding groove of PDZ-1 comprised of the β2 strand and the α2 helix, consistent with NMR analyses indicating that the extreme C-terminus of P1BM inserts in a perpendicular orientation into this region of PDZ-1 leading to inhibition of ligand binding (i.e., auto-inhibition)48. Interestingly, within all four structures the P1BM sequence bears a predicted alpha helix positioned between PDZ-1 and PDZ-2, corresponding to highly conserved residues including an arginine/lysine near the center of the alpha helix that is predicted to form a salt bridge with a conserved glutamate/aspartate residue 6 amino acids downstream (Figs. 1D and 5B). The only orthologues lacking this cationic residue are from the ctenophore species M. leidyi and H. californiensis. Nonetheless, the H. californica Mint P1BM bears a polar asparagine residue at this position which is predicted to hydrogen bond with the conserved glutamate downstream. Clearly, future studies to determine whether these predicted structural elements serve to project and/or stabilize the extreme C-terminus of P1BM into the ligand-binding groove of PDZ-1 are warranted.Fig. 5 AlphaFold2-predicted structures of Mint PDZ-1, PDZ-2 and the auto-inhibitory P1BM element are consistent with the separately resolved structures of human Mint-1 PDZ-1 and PDZ-2. (A) Ribbon diagrams of AlphaFold2-predicted structures of the tandem PDZ domains of select Mint proteins, along with the C-terminal P1BM auto-inhibitory domain. Ribbon diagrams of the predicted structures are rainbow colored, with N-terminal elements in blue and the C-terminal in red. The black ribbon structures structurally aligned with the Rattus norvegicus AlphaFold2 structure are NMR structures of the human PDZ-1 and PDZ-2 domains, which were resolved as separate structures in two different studies (PDB accession numbers 1U37 and 1Y7N, respectively). Note that for all examined PDZ-1/PDZ-2/P1BM structures, AlphaFold2 predicts the existence of a novel alpha helix that helps orient the C-terminal auto-inhibitory peptide into the PDZ-1 ligand-binding groove. (B) Depiction of predicted intra-molecular hydrogen bonds between a conserved glutamate near the C-terminus of Mint and amino acids upstream within a novel predicted helix of the P1BM element. (C) Depiction of the R. norvegicus PDZ-1/PDZ-2/P1BM AlphaFold2 structure color-coded according to sequence conservation across phyla using ConSurf83. (D) Primary sequence of the R. norvegicus Mint-1 PDZ1-PDZ2-P1BM elements, color-coded according to the ConSurf sequence conservation analysis. The black chevrons under the alignment indicate the position of the 17 amino acids correlated with ligand specificity41. For panels B and C, purple denotes highly conserved sequences/structural regions, while blue denotes highly variable sequences/regions (legend at the bottom of panel E).

Next, we used the program ConSurf83 to map sequence conservation of aligned PDZ domains (i.e., the single representatives from each phylum in Fig. 4A and B) onto the predicted PDZ-1/PDZ-2/P1BM structure of R. norvegicus Mint-1 (Fig. 5C and D). As expected, PDZ-2 shows stronger conservation compared to PDZ-1 along the entire structure, particularly along the β2 strand and residues in the α2 helix that together comprise the ligand-binding groove. Also deeply conserved is a TMP amino acid motif at the C-terminal end of β6 which demarks the boundary between PDZ-2 and P1BM. In PDZ-1, strong conservation is apparent at the interface between the β1-β2 linker and the β2 strand, a region that harbours the carboxyl accommodating loop for bound C-terminal ligands. In addition, the linkers between β2-β3 and the α1-β4 are considerably conserved within PDZ-1, despite the former being predicted as disordered (Figure S2)84.

Mint and CaV2c channels from the cnidarian sea anemone Nematostella vectensis interact in vitro and are co-expressed in dissociated neurons

We took advantage of the extensive whole animal transcriptome data available for N. vectensis in the NvERTx database85 to explore the expression patterns of Mint, CASK, Veli, and its three CaV2 channel isotypes through embryonic development, along with the transcription factors forkhead box protein L2 (FoxL2) and homeobox protein Otx1 which delineate two distinct populations of neurons in adult N. vectensis86. NvERTx also provides information on the developmental co-expression of genes, by partitioning genes with similar expression profiles into eight defined clusters85. During embryonic development, the CaV2c channel, as well as CASK, OtxC, and FoxL2 exhibit general increases in mRNA expression, most markedly within the first 25 h after fertilization (i.e., the blastula stage; Fig. 6A). This contrasts the mRNA expression of CaV2a, CaV2b, and Mint, which exhibit sharp decreases during the blastula stage, but become stabilized as embryos transition through the gastrula, planula, and juvenile stages. Notably, CaV2a and CaV2b are part of the same co-expression cluster (E-2), suggesting they are part of a common cellular gene regulatory network85. However, both cellular transcriptome analysis and mRNA in situ hybridization indicate these two genes are expressed in different cell types, with CaV2a being restricted to cnidocytes (stinging cells), and CaV2b being expressed more broadly in cells within the mouth and tentacles2,87. Also notable is that Mint, CASK, and Veli, which as noted form a conserved tripartite complex in bilaterians, fall within different co-expression clusters, and that none of the CaV2 channel genes share a cluster with Mint or the neuronal marker genes, but that CaV2c and CASK do fall within a common cluster.Fig. 6 Mint and CaV2c from the cnidarian N. vectensis interact in vitro. (A) Plot of mRNA transcript counts (Log2 of counts per million) of relevant gene transcripts in the developmental transcriptome N. vectensis obtained from the NvERTx database85. Included are the neuronal marker genes FoxL2 and OtxC which delineate two separate populations of neurons in the adult nervous system86. To the right of each line in the plot is an indication of the NvERTx co-expression cluster that each transcript is assigned to. BL = blastula, GA = gastrula, PL = planula, JU = juvenile. (B) Plot of average β-galactosidase activity (plus/minus standard deviation), normalized to a control positive interaction, measured for directed yeast 2-hybrid experiments between the N. vectensis Mint PDZ-1/PDZ-2/P1BM module and the CaV2a to CaV2c C-termini. The dashed red lines and asterisks denote statistical differences for two sample T-tests. (C) Map of the pEGFP-C1 plasmid construct use to express the N. vectensis CaV2c channel C-terminus protein in HEK-293T cells as a HA-tagged GFP fusion protein. (D) Illustration of the in vitro expressed CaV2c protein, depicting inclusion/exclusion of the 9 C-terminal amino acids bearing a putative Mint PDZ-1 domain ligand sequence. (E) Map of the pIRES2-DsRed2 plasmid construct use to express a portion of the Nematostella Mint PDZ1-PDZ2-P1BM protein with a Myc tag at the C-terminus. (F) Illustration of the in vitro expressed Mint protein. The lengths of the recombinant CaV2c and Mint proteins are denoted by a shared scale bar of 100 amino acids. (G) Co-immunoprecipitation experiment showing a robust interaction between the N. vectensis Mint and the CaV2c channel proteins in vitro. This interaction is completely disrupted upon deletion of the CaV2c channel C-terminal FDNDEETWCCOOH sequence. The red chevron denotes non-specific labeling of immunoglobulins by the α-Myc antibodies.

Considering the observed sequence and predicted structural similarity between the N. vectensis Mint PDZ-1 domain with its bilaterian orthologues (Figs. 4, S1, and 5), and the presence of conserved C-terminal DDWC-like motifs among all three cnidarian CaV2c channel isotypes (Fig. 1A), we reasoned that this protein–protein interaction might be conserved in cnidarians. To test this possibility, we conducted directed yeast 2-hybrid experiments to explore whether the C-terminal portion of N. vectensis Mint, containing the PDZ-1/PDZ-2/P1BM elements of the protein (i.e., D1141 to I1350 of NCBI accession number XP_032220836.1; Supplementary File 1), could bind to the 70 most distal amino acids of the three corresponding N. vectensis CaV2a, CaV2b, and CaV2c channel C-termini. Quantification of β-galactosidase activity, normalized to positive control values for known interacting proteins (murine p53 and the SV40 large T-antigen), revealed a strong interaction between the N. vectensis Mint PDZ-1/PDZ-2/P1BM protein and the 70 most distal amino acids of CaV2c, but not CaV2a or CaV2b (Fig. 6B). The interaction between Mint and CaV2c was completely disrupted upon removal of the last four amino acids of the channel C-terminus (i.e., ΔETWCCOOH). We next sought to validate this finding by conducting co-immunoprecipitation experiments with recombinant, epitope-tagged proteins expressed in HEK-293T cells. Through gene synthesis, we cloned the entire C-terminus of the N. vectensis CaV2c channel (i.e., S1632 to C1896 of NCBI accession number XP_032220836.1; Supplementary File 1) into the mammalian expression vector pEGFP-C1, with and without the 9 distal C-terminal amino acids bearing the putative PDZ ligand sequence ETWC (Fig. 6C). Expressed in vitro, this construct produces a chimeric protein with enhanced green fluorescent protein at its N-terminus (Fig. 6D), which we found to enhance in vitro expression of the disordered channel C-terminus, plus a 9 amino acid haemagglutinin (HA) epitope tag used for the pulldown experiments. Similarly, we cloned the C-terminal portion of N. vectensis Mint, bearing the PDZ-1/PDZ-2/P1BM elements (E1076 to I1350 of NCBI accession number XP_032233969.1; Supplementary File 2), into the bicistronic expression vector pIRES2-DsRed2 (Fig. 6E). This construct expresses red fluorescent protein separately from the recombinant Mint protein, the latter engineered to contain a 10 amino acid myc epitope tag at its C-terminus for the pulldown experiments (Fig. 6F). Western blotting of the two recombinant CaV2c channel proteins expressed in HEK-293T cells, using monoclonal anti-HA antibodies, produced bands of appropriate molecular weight for both the wildtype channel C-terminus (i.e., 59.5 kDa), and the C-terminal deletion variant (i.e., ΔPDZ-L; 58.3 kDa; Fig. 6G). Similarly, Western blotting of the recombinant Mint protein with anti-myc antibodies revealed an appropriate band of ~ 30.8 kDa. In three separate pull-down experiments of the Mint protein, we were able co-immunoprecipitate the wildtype CaV2c channel C-terminus, but not the deletion variant lacking 9 C-terminal amino acids (Fig. 6G). Thus, the PTB/PDZ-1/PDZ-2/P1BM module of the N. vectensis Mint can bind the CaV2c channel C-terminus in vitro, but only in the presence of the extreme C-terminus bearing the ETWC motif. Altogether, these analyses indicate that the capacity for Mint to bind the CaV2 channel C-terminus is conserved between cnidarians and bilaterians.

Next, we sought to determine whether the mRNAs of these two genes are co-expressed in dissociated neurons, along with the neuronal marker genes OtxC and FoxL286. Co-labelling of freshly dissociated cells from tentacle bud stage N. vectensis specimens with fluorescently labeled RNAscope probes against CaV2c, FoxL2, and Mint (Fig. 7A to C′) revealed that most cells that express CaV2c also express Mint, and that roughly one fifth of labelled cells express CaV2c, Mint, and FoxL2 together (18.1 ± 1.6%) (Fig. 7G). Separate experiments on cells co-labelled with probes against CaV2c, OtxC, and Mint (Fig. 7D to F′) revealed similar co-expression of CaV2c and Mint (26.8 ± 2.4% of 484 cells), and of CaV2c, Mint, and OtxC together (12.9 ± 1.4%) (Fig. 7H). In both experiments co-expression of all three mRNA targets was observed in cells with morphology characteristic of unipolar neurons (Fig. 7C, D, and E), bipolar neurons (Fig. 7A and F), and tripolar neurons (Fig. 7B), identified as distinct neuron subtypes in transgenic N. vectensis polyps88. Given that FoxL2 and OtxC are expected to label neurons, it appears as though the CaV2c channel is expressed in neurons, together with Mint. Furthermore, Mint appears to have broader cellular expression than the other tested transcripts, since 23.9 ± 4.3% and 26.9 ± 1.9% of all cells in the two separate experiments were labeled only with the Mint probe (Fig. 7G and H). Morphologically, some of Mint-positive cells resemble gland cells (e.g., middle cell on Fig. 7A or upper right cell on Fig. 7E), shown to exhibit high diversity in N. vectensis polyps via single cell RNA-seq86,89.Fig. 7 Mint and CaV2c channels are co-expressed in dissociated N. vectensis neurons. (A to C’) Transmitted light (A to C) and fluorescent (A’ to C’) micrographs of dissociated N. vectensis cells subject to in situ hybridization with RNAscope probes for the CaV2c channel (green), Mint (magenta), and FoxL2 (red). (D to F’) Transmitted light (D to F) and fluorescence (D’ to F’) micrographs of dissociated cells probed for the CaV2c channel (green), Mint (magenta), and OtxC (red). White chevrons denote cells where all three mRNAs could be detected. Scale bars represent 10 µm. (G) Venn diagram showing average single cell counts for co-expression of Mint, CaV2c and FoxL2, showing that most cells expressing CaV2c also express Mint, and that over 1/3 of cells expressing CaV2c also express the neural marker FoxL2. (H) Venn diagram showing average single cell counts for co-expression of Mint, CaV2c and OtxC, showing that most cells expressing CaV2c also express Mint, and under 1/3 of cells expressing CaV2c also express the neural marker OtxC.

To validate these observations, we mined the recently updated single cell transcriptome data available for N. vectensis90. Apparent is that CaV2c mRNA expression is largely restricted to neurons and related secretory cells, as well as mucin-producing secretory gland cells (Fig. 8A and B). Indeed, CaV2c expression is much more prevalent in neuronal metacells comparted to CaV2a and CaV2b, present in 401 out of 3709 neuronal metacells versus only 57 for CaV2a and 56 for CaV2b. Mint also shows prominent expression in neurons (427 neuronal metacells) and secretory cells but is also expressed more broadly in other cell types including cnidocytes, gastrodermis, and ectoderm (Fig. 8C). Consistent with their suggested classification as neuronal marker genes86, FoxL2 and OtxC show enriched expression in neurons (459 and 295 metacells, respectively), but also moderate expression in a few other cell types (e.g., FoxL2 in cnidocytes and retractor muscle, OtxC in gastrodermis and embryonic mesoderm; Fig. 8D and E). Further analysis confirmed co-expression of Mint and CaV2c in 72 metacells, and triple expression of CaV2c, Mint, and either FoxL2 or OtxC in 13 metacells (Fig. 8F and G). Thus, the single cell RNA-Seq data confirms neuronal co-expression of Mint and CaV2c channels. Nonetheless, the low percentage of neuronal metacells co-expressing these two genes, which is considerably lower than what we observed in our in situ experiments (Fig. 7), prompted us to explore whether the RNA-Seq data might underestimate real co-expression values. To do this, we explored the neuronal expression of the two N. vectensis CaVβ ancillary subunits, β1 and β2, which are critical components of CaV1 and CaV2 channel complexes required for proper protein folding, membrane expression, and function of the pore-forming CaV subunit91. Hence, there is an expectation that CaV2c and either β1 or β2 should be highly co-expressed in neurons. Interestingly, 620 neuronal metacells express the β1 subunit (accession number NVE4481), whereas only 18 express β2 (NVE10585), indicating that the former is the predominant neural subunit. An examination of co-expression with CaV2c identified only 76 metacells for β1, and 5 for β2. Hence, of the 401 neuronal metacells that express CaV2c, 18.0% co-express Mint, 19.0% co-express β1, and only 1.2% co-express β2. Based on these numbers, it seems reasonable to infer that the co-expression value for CaV2c and Mint extracted from the RNA-Seq data underestimates co-expression that occurs in neurons.Fig. 8 Mint and CaV2c channels are co-expressed in N. vectensis neuronal metacells. (A) Complete cell state atlas of N. vectensis (from90). (B) Metacell expression of CaV2c mRNA. (C) Metacell expression of Mint mRNA. (D) Metacell expression of FoxL2 mRNA. (E) Metacell expression of OtxC mRNA. (F) Overlay of metacells expressing CaV2c, Mint, and FoxL2 within the neuronal cluster corresponding to the box in panel A. (G) Overlay of metacells expressing CaV2c, Mint, and OtxC within the neuronal cluster corresponding to the box in panel A. For panel A, different colors are used to indicate different metacell clusters. For panels B to E, accession numbers for each gene are listed in brackets on the right of each gene name, and the legends within each panel are used to indicate normalized mRNA expression levels within each metacell. The Venn diagrams in the upper left of panels F and G illustrate the color scheme for co-expression within the different metacells.

Mint PDZ domains from placozoans and ctenophores interact with CaV2 channel C-termini in vitro

Placozoans like T. adhaerens are the most early diverging animals that possess all three types of metazoan CaV channel types, CaV1 to CaV32,4, 92. Previously, the three T. adhaerens CaV channels, TCaV1 to TCaV3, were shown to conduct voltage-gated calcium currents in vitro with properties like those of their bilaterian orthologues20,30, 93. Nonetheless, the absence of a DDWC-like C-terminus in placozoan CaV2 channels (Fig. 1A) suggests that the CaV2-Mint interaction is not conserved in these animals. Alternatively, given the strong sequence divergence of the T. adhaerens Mint PDZ domains (Figs. 2C and S1), the interaction may occur but through unique ligand binding properties of PDZ-1. To explore these possibilities, we sought to emulate the in vitro co-immunoprecipitation experiments done for N. vectensis using the corresponding T. adhaerens proteins. Thus, we similarly cloned the C-terminus of TCaV2 (i.e., R1646 to V2093 of NCBI accession number QOQ37977.1; Supplementary File 1) into the pEGFP-C1 plasmid and tagged it with EGFP and an HA epitope at its N-terminus, and the PTB/PDZ-1/PDZ-2/P1BM region of T. adhaerens Mint (K760 to I1034 of NCBI accession number GHJI01000994.1; Supplementary File 2) into pIRES2-DsRed2 tagged with a C-terminal myc tag (Figure S3A and B). Unfortunately, the recombinant TCaV2 protein proved difficult to express in HEK-293T cells (Figure S3C), perhaps due to its extended length of 448 amino acids compared to only 265 for N. vectensis CaV2c. As a result, an interaction between Mint and TCaV2 could neither be observed nor ruled out via co-immunoprecipitation.

As an alternative approach, we sought to identify putative binding partners of the T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein, as well as PDZ-1 by itself, by generating and screening a custom yeast 2-hybrid library derived from T. adhaerens whole animal mRNA. Searching this library with PDZ-1/PDZ-2/P1BM as bait (i.e., K825 to I1034) identified 12 putative binding partners, all but three bearing a valine residue at the extreme C-terminus (Table 1A). Two of the identified prey proteins, homologues of cathepsin D and a G protein γ subunit, produced 4 and 2 hits respectively, the latter notable in bearing a C-terminal sequence strikingly similar to that of TCaV2 (i.e., SRCTLVCOOH vs. SKCTAVCOOH). Using just Mint PDZ-1 as bait identified a total of 32 putative interactions (Table 1B), 17 of which corresponded to an artificial peptide encoded within the 3′ UTR of a betaine-homocysteine S-methyltransferase 1 homologue. Interestingly, this non-biological peptide has a C-terminal sequence ending with tyrosine and isoleucine (i.e., YICOOH), resembling the Mint P1BM C-terminus of ETPEYICOOH which based on homology and structural prediction is expected to insert into the PDZ-1 ligand-binding groove (Figs. 1D, 5, and S2). Two additional prey sequences pulled from the cDNA library also had C-termini with terminal isoleucine residues, altogether suggesting that the PDZ-1 domain of T. adhaerens Mint has a high affinity for peptides bearing isoleucine in the last position. Although we did not directly test the interaction between the T. adhaerens Mint PDZ-1 domain to the downstream P1BM element, the strong binding of this PDZ-1 to P1BM-like C-terminal peptides suggests that auto-inhibition is an ancient and deeply conserved trait, also suggested by the deep sequence conservation of the P1BM element (Fig. 1D). The second-most prominent identified prey sequences had C-termini ending with WVCOOH (6 hits), suggesting that like the PDZ-1 domain of mammalian Mint-1, the T. adhaerens orthologue has the capacity to bind C-terminal ligands with a tryptophan at position −1, coupled to a valine at position 0. Table 1 Annotation of T. adhaerens proteins found to interact with (A) the Mint PDZ-1/PDZ-2/P1BM and (B) the Mint PDZ-1 proteins in yeast 2-hybrid screens.

Sequence ID	C-terminal Seq.	SmartBlast ID	Search method	E-value	Query Cover	% ID	Species	Accession	Colony count	
(A) PDZ-1/PDZ-2/P1BM	
evg93589	QNRLGFAKAV	Cathepsin D	SmartBlast	1.0E−140	95.0%	53.4%	D. melanogaster	NP_001334713.1	4	
TRINITY_DN22496_c0_g1_i1	FVEKSRCTLV	G protein Gamma 7	SmartBlast	1.0E−17	100.0%	50.0%	D. rerio	NP_001002397.1	2	
evg19742	NSILSLPTTV	H/ACA ribonucleoprotein complex subunit 2	SmartBlast	4.0E−48	71.0%	48.7%	M. musculus	NP_080907.1	1	
evg1173499	IHIYECPSNV	SKI8 subunit of superkiller complex protein	SmartBlast	3.0E−156	97.0%	65.8%	D. rerio	NP_957147.1	1	
evg485730	REQDKRGNLV	4-hydroxyphenylpyruvate dioxygenase	SmartBlast	0.0E+00	98.0%	70.7%	D. rerio	XP_005171786.1	1	
evg1106341	SGGPTIEEVD	Heat shock cognate 71 kDa protein	SmartBlast	0.0E+00	93.0%	88.4%	H. sapiens	NP_006588.1	1	
evg1508575	GKSLRTVLNF	Alcohol dehydrogenase class-3	SmartBlast	9.0E+00	94.0%	74.3%	D. rerio	NP_571924.2	1	
evg264334	RTERPRPRAY	40S ribosomal protein S2	SmartBlast	1.0E−136	84.0%	79.9%	M. musculus	NP_032529.2	1	
(B) PDZ-1	
evg10473	LDINGFIIYI	3′ UTR of: betaine-homocysteine S-methyltransferase 1	SmartBlast	7.0E−176	99.0%	65.9%	D. rerio	NP_001012498.1	17	
evg929466	NKETIDAMEI	60S ribosomal protein L9	SmartBlast	1.0E−89	94.0%	67.0%	D. rerio	XP_005170964.1	1	
evg1428490	LIESRNLAVI	Alt. ORF: Acid sphingomyelinase-like phosphodiesterase 3b	SmartBlast	3.0E−32	78.0%	26.2%	M. musculus	NP_598649.1	1	
Not found	EVLSGSMFWV	Basic proline-rich protein-like	NCBI BLAST	3.0E−10	97.0%	42.4%	C. pelagica	XP_009995607.1	4	
evg1363342	KLFHGWQAWV	Serine-protein kinase ATM	SmartBlast	0.0E+00	100.0%	34.5%	H. sapiens	XP_011541144.1	1	
evg1036905	AQDGVIKVWV	Phosphoinositide 3-kinase regulatory subunit 4	SmartBlast	0.0E+00	99.0%	46.2%	D. rerio	XP_005158356.1	1	
evg1095586	GHLTVLVAMV	Msx2-interacting protein	SmartBlast	3.0E−50	61.0%	46.2%	M. musculus	NP_001334164.1	1	
evg1187940	ILNTTGAFNV	Multidrug and toxin extrusion protein 1	SmartBlast	5.0E−107	87.0%	36.7%	M. musculus	NP_080459.2	1	
evg1106341	SGGPTIEEVD	Heat shock cognate 71 kDa protein	SmartBlast	0.0E+00	93.0%	88.4%	H. sapiens	NP_006588.1	1	
evg108179	PSEISVTLRR	Carboxypeptidase D	SmartBlast	7.0E−138	97.0%	52.6%	H. sapiens	NP_001186704.1	1	
evg1107053	LIKIDRKEEL	Peptidyl-prolyl cis–trans isomerase FKBP2	SmartBlast	1.0E−58	95.0%	64.0%	M. musculus	NP_001159840.1	1	
evg1256402	MCIEINSRLL	Acetyl-CoA acetyltransferase	SmartBlast	9.0E−177	98.0%	59.3%	D. rerio	NP_571445.2	1	
evg1382077	STKLEICSNS	lysine-specific demethylase 5A	SmartBlast	0.0E+00	82.0%	37.3%	H. sapiens	NP_001036068.1	1	

Intrigued by the similarity between the putative PDZ ligand sequence of the identified Gγ protein (SRCTLVCOOH) and the TCaV2 channel (SKCTAVCOOH), we next explored whether T. adhaerens Mint can bind the 70 most distal amino acids of the channel in directed yeast 2-hybrid experiments as we did previously for N. vectensis (Fig. 6B). These experiments revealed a moderate interaction between the T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein and the 70 most distal amino acids of TCaV2 (Fig. 9), which was completely disrupted upon removal of the last four amino acids of the channel C-terminus (i.e., ΔCTAVCOOH). Instead, experiments using just PDZ-1 domain of T. adhaerens Mint failed to show an interaction with the TCaV2 C-terminus. Since our yeast 2-hybrid screen with PDZ-1 on its own identified putative ligands with a tryptophan residue in the −1 position (Table 1), we reasoned that this PDZ domain might share conserved binding properties with mammalian Mint PDZ-1 in this regard. Thus, we expanded these experiments to test in vitro interactions between T. adhaerens Mint and the CaV2c channel from N. vectensis bearing the C-terminal sequence ETWCCOOH. This revealed a moderate but statistically significant interaction (Fig. 9). Interestingly, we did not detect an interaction between the complete T. adhaerens Mint PDZ-1/PDZ-2/P1BM module at the ETWCCOOH ligand, perhaps because the auto-inhibitory sequence of P1BM binds PDZ-1 to restrict its ability to bind the N. vectensis C-terminus (Fig. 9).Fig. 9 Yeast 2-hybrid analyses reveal deep conservation of the Mint-CaV2 interaction. Plot of average β-galactosidase activity (plus/minus standard deviation), normalized to a control positive interaction, measured for various directed yeast 2-hybrid experiments. The dashed red lines and asterisks denote statistical differences for two sample T-tests.

We also tested in vitro interactions between ctenophore Mint and CaV2 from the species H. californiensis (Mint) and M. leidyi (CaV2) which have highly divergent Mint PDZ domains (Figs. 4 and S1) and CaV2 channel C-terminal sequences that are distinct from those of placozoans and cnidarians (Fig. 1A). These experiments failed to detect an interaction between H. californiensis Mint (V344 to L547 of evg127856; Supplementary File 2) and the M. leidyi CaV2 C-terminus (QOY24596.1; Supplementary File 1). Furthermore, neither the T. adhaerens nor the N. vectensis PDZ-1/PDZ-2/P1BM modules bound the M. leidyi CaV2 channel C-terminal motif of SGDICOOH. To gain some insights into the binding properties of ctenophore Mint, we conducted cross-species experiments, finding strong interactions between the H. californiensis Mint PDZ-1/PDZ-2/P1BM and both the ETWCCOOH ligand of N. vectensis CaV2c, and the CTAVCOOH ligand of T. adhaerens CaV2 (Fig. 9). The ability of ctenophore Mint to bind the distinct C-terminal ligands of placozoan and cnidarian CaV2 channels prompted us to test whether N. vectensis Mint, like placozoan and ctenophore Mint, has dual binding capabilities in being able to bind the placozoan CTAVCOOH ligand in addition to its putative cognate ligand of ETWCCOOH. These experiments revealed a moderate interaction between N. vectensis PDZ-1/PDZ-2/P1BM and the TCaV2 channel C-terminus (Fig. 9).

Together, these various experiments reveal that T. adhaerens Mint can bind the TCaV2 channel C-terminus in vitro, albeit in a non-canonical manner as the binding is not mediated by PDZ-1 on its own. Like PDZ-1 in mammals, T. adhaerens Mint PDZ-1 can bind peptides with a tryptophan in the −1 position. Furthermore, the experiments demonstrate that cnidarian, placozoan, and ctenophore Mint PDZ domains can all recognize the placozoan CaV2 ligand of CTAVCOOH, but that none of the Mint PDZ-1/PDZ-2/P1BM modules from these animals bind the ctenophore CaV2 channel C-terminus in vitro.

Bacterial 2-hybrid screening uncovers unique binding properties of the T. adhaerens Mint tandem PDZ domains

Next, we sought to characterize the preferred ligand specificities of Mint PDZ domains from T .adhaerens, N. vectensis, and H. californiensis, as well as the brown rat R. norvegicus, by conducting bacterial 2-hybrid assays94. For each species, the protein-coding sequences of the PDZ-1/PDZ-2/P1BM module or PDZ-1 alone were cloned into the bacterial 2-hybrid vector pGHUC-ω-Erbin. These constructs were then screened for interactions with a library of random hexameric C-terminal peptides fused to a DNA-binding element in vector pB1H2-UV5-Zif268, permitting identification of ligand peptides en masse via Illumina sequencing94. Consistent with the ability of the rat Mint-1 PDZ-1 domain to bind the DDWC and DHWC sequences of CaV2.1 and CaV2.2 channels46,47, bacterial 2-hybrid screening with this domain identified a preferred ligand sequence of E/DWLCOOH, resembling that of CaV2 channels with an invariant tryptophan residue at position −1 preceded by a negatively-charged glutamate (E) or aspartate (D) at position −2 (Fig. 10A). Notable however is the presence of a leucine at position 0 of the consensus ligand sequence, contrasting the cysteine residue conserved in CaV2 channels (Fig. 1A). Screening the entire PDZ-1/PDZ-2/P1BM module of rat Mint-1 revealed a markedly different ligand specificity, with a reduced selection for tryptophan at position −1, consistent with auto-inhibition of PDZ-1 by P1BM, as well as a slight enrichment for hydrophobic residues at positions −2 to −4, and a reduced preference for leucine at position 0 in lieu of the also hydrophobic residue valine.Fig. 10 Bacterial 2-hybrid experiments, structural modelling, and mutagenesis shed light the unique ligand-binding properties of T. adhaerens Mint. (A) Sequence logos for the preferred ligand signatures of the R. norvegicus Mint PDZ-1/PDZ-2/P1BM protein (left), and the PDZ-1 domain by itself (right), determined via bacterial 2-hybrid screening. (B) Sequence logos for the preferred binding signatures of the T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein (left), and the PDZ-1 domain by itself (right), determined via bacterial 2-hybrid screening. (C) Sequence logo for the preferred ligand signature of the N. vectensis Mint PDZ-1/PDZ-2/P1BM protein determined via bacterial 2-hybrid screening. (D) AlphaFold2 model of the R. norvegicus Mint-1 auto-inhibitory peptide of the P1BM element embedded within the PDZ-1 ligand-binding groove predicting. Weak Van der Waals interactions between the C-terminal isoleucine of P1BM and several amino acids within the PDZ-1 β1-β2 loop and β2 stand that would serve to stabilize the interaction are shown. (E) AlphaFold2 model of the T. adhaerens Mint auto-inhibitory peptide-PDZ-1 interaction, revealing unique hydrogen bonds between a glutamate residue two amino acids upstream of the C-terminal isoleucine (E1034) and two lysine residues in the α2 helix pf PDZ-1 (K920 and K904). (F) Plot of average β-galactosidase activity (plus/minus standard deviation), normalized to a control positive interaction measured for various directed yeast 2-hybrid experiments. The dashed red lines and asterisks denote statistical differences for paired T-tests (Legend: wt = wildtype T. adhaerens PDZ-1/PDZ-2/P1BM or PDZ-1; TC2 = TCaV2; NC2c = N. vectensis CaV2c; ΔCT = last four amino acids removed). (G) Solved x-ray crystal structure of the human erbin PDZ domain bound to its preferred peptide ligand ETWVCOOH (PDB number 6Q0N). (H) Modelled structure of the T. adhaerens Mint PDZ-1 domain bound to its recognized ligand ETWCCOOH from the C-terminus of the N. vectensis CaV2c channel. (I) Modelled structure of the T. adhaerens Mint PDZ-1 domain bound to the non-recognized ligand CTAVCOOH from the C-terminus of the T. adhaerens CaV2 channel.

For the T. adhaerens Mint PDZ-1 domain, the screening identified a preferred ligand sequence of GxWVCOOH (Fig. 10B), which is consistent with the observed interaction of this protein with the N. vectensis CaV2c channel C-terminus in our directed yeast 2-hybrid experiments (Fig. 9), and C-terminal sequences of several interacting proteins identified in the yeast 2-hybrid screen (Table 1). Also consistent, the T. adhaerens PDZ-1/PDZ-2/P1BM module produced a ligand specificity of GTxVCOOH, corroborating its ability to bind the TCaV2 channel C-terminus in directed yeast 2-hybrid experiments (Fig. 9), as well as other proteins bearing TxVCOOH or VCOOH C-termini in the yeast 2-hybrid screening (Table 1). Notably, like rat Mint-1, the T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein exhibits a diminished preference for tryptophan in the second last position, consistent with our inability to detect an interaction between this protein and the N. vectensis CaV2c channel C-terminus in the yeast 2-hybrid assays (Fig. 9), and with the occurrence of auto-inhibition of PDZ-1 by P1BM. A notable and unique feature of the T. adhaerens Mint PDZ-1/PDZ-2/P1BM compared to the rat and N. vectensis orthologues is a strong preference for threonine at position −2, and a common preference of PDZ-1/PDZ-2/P1BM and PDZ-1 alone to bind glycine and valine at positions −3 and 0, respectively (Fig. 10A). As such, the most notable difference in ligand selectivity between the two T. adhaerens proteins are a tryptophan at position −1 for PDZ-1, instead of a threonine at position −2 for PDZ-1/PDZ-2/P1BM.

Unfortunately, neither of the two H. californinesis PDZ domain constructs produced viable ligand specificity signatures, nor did the N. vectensis PDZ-1 domain. However, we were able to generate a selectivity signature of WVCOOH for the N. vectensis Mint PDZ-1/PDZ-2/P1BM protein, which contrasts the signatures for rat and T. adhaerens orthologues by bearing a moderate preference for tryptophan at position −1 and corroborates the co-immunoprecipitation and yeast 2-hybrid experiments where we detected robust in vitro interactions of this protein with the ETWC ligand sequence of CaV2c.

Indirectly, our analyses suggest that auto-inhibition of PDZ-1 by the P1BM module is conserved in the T. adhaerens Mint orthologue. To gain more insights into this possibility, we examined the AlphaFold2 structures of the rat and T. adhaerens PDZ-1/PDZ-2/P1BM modules, focusing on the predicted interaction between PDZ-1 and P1BM (Fig. 10D and E). For both proteins, the P1BM element was predicted to insert the C-terminal YICOOH motif into the PDZ-1 ligand-binding groove in a perpendicular orientation compared to the canonical orientation of PDZ domain C-terminal ligands (Fig. 10D and E). For the rat orthologue, the C-terminal isoleucine of P1BM (I841) forms predicted hydrogen bonds with the residues L671, G672, and V673 comprising the highly conserved ligand carboxylate-binding loop in PDZ-1 (Fig. 10D). This predicted arrangement is highly consistent with previous NMR studies which revealed a similar orientation of the C-terminal element of P1BM within the PDZ-1 ligand-binding groove, with hydrogen bonding between I841 and the residues L671 and G67248.

For the T. adhaerens protein, homologous interactions are predicted for the C-terminal residue I1034 of P1BM, which in the model forms hydrogen bonds with the carboxylate-binding loop residues L864, G865, and I866 of PDZ-1 (Fig. 10E). However, there are also some additional hydrogen bonds predicted between Y1033 of P1BM with N863, as well as stronger salt bridges between an upstream glutamate in P1BM (E1032) and a set of lysine residues in the PDZ-1 domain α1 helix (K920 and K924). Interestingly, whereas the basic K924 residue is deeply conserved in the PDZ-1 domains of Mint, the upstream residue K920 appears to be a unique feature of T. adhaerens Mint, with most other orthologues possessing a glutamine in this position (Fig. 4A). We therefore decided to explore whether this unique residue contributes to the unique ligand-binding properties of T. adhaerens Mint in vitro. Specifically, we repeated the directed yeast 2-hybrid experiments to test interactions between the T. adhaerens Mint PDZ-1/PDZ-2/P1BM and PDZ-1 proteins, and the TCaV2 and N. vectensis CaV2c channel C-termini, comparing wildtype Mint variants to ones bearing a K920Q mutation within PDZ-1. Interestingly, this amino acid substitution completely disrupted the interaction between the T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein and the TCaV2 channel C-terminal motif of CTAVCOOH, while only reducing the interaction between PDZ-1 alone and the ETWCCOOH motif of N. vectensis CaV2c (Fig. 10F). Thus, while the unique K920 residue in PDZ-1 of T. adhaerens Mint does not seem critical for PDZ-1 binding to DDWC-like ligands, it plays a major role in the ability of PDZ-1/PDZ-2/P1BM to recognize and bind the TCaV2 CTAVCOOH ligand.

Lastly, to further explore why the T. adhaerens Mint PDZ-1 domain is able bind the N .vectensis CaV2c peptide of ETWCCOOH, but not its cognate CaV2 peptide of CTAVCOOH, we modelled these interactions for comparison with the solved X-ray crystal structure of the erbin PDZ domain bound to its preferred (and similar) ligand peptide ETWVCOOH95 (PDB number 6Q0N). Interestingly, in both the resolved structure of erbin-ETWVCOOH and the modelled structure of Mint PDZ-1-ETWCCOOH (Fig. 10G to I), a lysine residue in the distal end of the PDZ α1 helix (i.e., K87 in erbin and the conserved K924 residue in Mint) hydrogen bonds with the backbone carbonyl oxygen of the peptide tryptophan, W−1, orienting its aromatic side chain over the proximal end of the β2 strand (Fig. 10G and H). Given the deep conservation of an acidic residue in this position of Mint PDZ-1 (Fig. 4A), including in all orthologues from species where an interaction with CaV2 channel C-termini bearing W−1 residues has been confirmed in vitro (Figs. 1A, 10A to C), this putative ligand binding feature appears ancient in nature, perhaps shared with other PDZ domains that fall into the 5 of 16 defined ligand specificity classes that select for W−1 residues41. Furthermore, although distinct from erbin in lacking a histidine residue in α1 helix that hydrogen bonds with the ligand T−2 side chain (Fig. 10G), the unique K920 residue in the central region of the T. adhaerens Mint α1 helix is predicted to form an analogous interaction with T−2 of the ETWCCOOH peptide (Fig. 10H). This modelled interaction is consistent with the reduction in β-galactosidase activity in the directed yeast 2-hybrid experiments upon mutation of K920 to glutamine (Fig. 10F), reflecting reduced peptide binding affinity perhaps caused by disrupted hydrogen bonding with T−2. Instead, the modelled structure of Mint PDZ-1-CTAVCOOH suggests a greater distance between the K924 side chain and the carbonyl oxygen of A−1, precluding hydrogen bonding, and a clash between side chain sulfur atoms from the C−3 residue of the peptide and a methionine residue within the β2 strand of PDZ-1 (M865), altogether consistent with our inability to detect an interaction between these two protein elements in vitro. In sum, our modelling studies indicate that an interaction between T. adhaerens Mint PDZ-1 and the ETWCCOOH peptide is at least possible, providing some putative molecular determinants. However, more direct analyses, such as molecular dynamics and/or structural determination, will be required to confirm the molecular nature of this interaction and the various other interactions between Mint proteins and CaV2 channel C-termini.

Discussion

High-throughput DNA sequencing has allowed us to learn much about the evolutionary origins of genes that give rise to complexity in animals, including those crucial for electrical and synaptic signaling in the nervous system. Lagging however is our understanding of whether such genes exhibit orthologous protein–protein interactions and cellular functions in the earliest animal lineages, or in pre-metazoan organisms, as defined in bilaterians through wet lab research96. One recent pioneering study employed ancetral gene reconstruction to explore the interaction between the third PDZ domain of the postsynaptic scaffolding protein Discs Large Homolgue (DLG) and its binding partner CRIPT, finding it emerged prior to animals in early eukaryotes, but lost binding affinity in select ecdysozoan invertebrates suggesting that even deep ancestral interactions can be subject to dynamic evolution97. A few other studies have similarly employed biochemical techniques to explore the origins of select protein–protein interactions in vitro98,99, or the structure of PDZ-ligand complexes from pre-metazoans100. Indeed, continued efforts in this regard will be important as we seek to integrate knowledge about the co-existence of genes in different organismal genomes, with functional data comparing how these genes function in non-mainstream animal systems in particular early diverging animals. This, combined with explorations into the electrophysiological and genetic properties of non-bilaterian neurons and synapses, has the potential to answer important questions about nervous system evolution.

In this study, we explored the evolutionary origins of the first identified interaction between a PDZ domain and a calcium channel C-terminus, Mint and CaV246, and conducted phylogenetic analyses to explore the origins and domain structure of several genes that associate with Mint in bilaterians. Two caveats that must be considered when exploring the presence and evolution of genes is that sequences can potentially be missed under a given sequence identification strategy, and that future sequencing efforts might uncover homologues in organismal lineages previously reported to lack these genes. Bearing this in mind, our current analysis indicates that Mint/X11 proteins are found broadly in animals, present in all examined animal phyla (Fig. 1D). Notably, we used BLAST to search for Mint homologues outside of animals and were unable to identify sequences bearing a conserved PTB/PDZ-1/PDZ-2/P1BM module. The closest sequence with BLAST homology to Mint was from the holozoan eukaryote Sphaeroforma arctica, bearing two C-terminal PDZ domains but lacking a PTB domain (NCBI accession number XP_014158662.1). Based on these observations, Mint appears to be unique to animals having evolved via fusion of pre-existing protein domains, a rarer form of gene evolution refered to as type III novelty101. Homologues for CASK and Veli, which form a tripartite complex with Mint in vertebrates and bilaterian invertebrates63–65, were also found broadly in animals, most bearing respective CaMK and L27 domains required for complexing with Mint (Fig. 2A and B). However, whether these indeed form complexes in non-bilaterian animals is uncertain given the strong sequence divergence of the Mint N-terminus where the CASK Interaction Domain is found in vertebrates and bilaterian invertebrates (e.g., D. melanogaster and C. elegans62; Fig. 1E). In our sequence alignments, we failed to identify conserved signatures corresponding to the CID in non-bilaterian Mint proteins, comprised of an amphipathic alpha helix, a proline-rich linker, and a core binding region bearing an essential tryptophan-valine (WV) motif62.

Interestingly, we found strong conservation of a YENPTY motif in APP homologues from cnidarians to bilaterians (Fig. 2C). This motif is bound by the PTB domain of Mint50 and is required for membrane internalization and subsequent proteolytic processing of APP into amyloid β (Aβ)102,103, a product which is subsequently secreted and accumulates in the brain of Alzheimer’s patients in the form of plaques104. Placozoan APP homologues on the other hand lack this motif, while ctenophores appear to lack APP homologues altogether. Instead, all examined animals possess presenilin homologues with conserved C-terminal motifs recognized by both PDZ-1 and PDZ-2 of mammalian Mint 1 and 2 proteins48,54, 55 (Fig. 2D). This particular interaction links Mint to the production of Aβ, wherein Mint promotes internalization of presenilin (the core subunit of γ-secretase), as well as activity-dependent colocalization of presenilin and APP which presumably precedes cleavage of the latter into Aβ56,105. With respect to neurexin, we identified homologues in most examined animals bearing conserved C-terminal sequences that associate with the Mint-1 PDZ-1 domain from mammals in vitro45,47 (Fig. 3). Mint has been implicated in the proper pre-synaptic localization of neurexin106, and through its association with presenilin, perhaps also regulates the proteolytic processing of neurexins by this enzyme at synapses107, two important processes for synaptic homeostasis and function. Indeed, future experiments, for example proteomic analysis of Mint-interacting proteins, will be needed to determine whether the various conserved motifs observed in this study mediate interactions with Mint in non-bilaterians and other invertebrates.

This work provides the first experimental evidence that Mint proteins from non-bilaterians can bind CaV2 channel C-termini in vitro. First, we identified a robust interaction between the N. vectensis Mint PDZ-1/PDZ-2/P1BM and the CaV2c channel C-terminus, but not the other two isotypes, CaV2a and CaV2b (Figs. 6B, G and 9). A major difference in the C-terminal ligands of these channels is that CaV2c bears a threonine at position −2 (i.e., ETWCCOOH), while CaV2 and CaV2b have more canonical CaV2 channel motifs bearing −2 residues of aspartate and glutamate respectively (i.e., DDWC and SEWC; Fig. 1A). In this study, we did not delineate whether the interaction between N. vectensis CaV2c and Mint was mediated by PDZ-1 or PDZ-2, since we utilized the full PDZ-1/PDZ-2/P1BM module for all experiments. Our attempts to determine the ligand specificity of PDZ-1 by itself for N. vectensis Mint, which would have provided important insights, were unfortunately unsuccessful. We did however succeed in rendering a bacterial 2-hybrid specificity profile for the N. vectensis PDZ-1/PDZ-2/P1BM module, finding a preference for ligands with tryptophan and valine at positions −1 and 0 respectively (Fig. 10C). Hence, the preference for tryptophan at position −1 is consistent with our ability to detect an interaction with the CaV2c ETWCCOOH ligand, where perhaps, the threonine at position −2 is permissive for binding while the asparate/glutamate residues in CaV2a and CaV2b disrupt binding. This same module bound the T. adhaerens CaV2 ligand sequence of CTAVCOOH, albeit with lower affinity than ETWCCOOH as approximated by β-galactosidase activity (Fig. 9), perhaps attributable to the threonine at position −2 and a preferred valine at position 0.

In addition to the in vitro interaction between N. vectensis Mint and CaV2c, our fluorescence in situ hybridization experiments revealed overlapping mRNA expression of these two genes in neurons expressing the marker genes FoxL2 or OtxC86 (Fig. 7). Furthermore, examination of the available single cell RNASeq data for this species confirmed co-expression of these genes in neurons (Fig. 8). Thus, it is at least possible that these proteins also physically interact in vivo, but future experiments will be required to determine whether this is indeed the case. Expression of CaV2c in N. vectensis neurons is consistent with the key role that CaV2 channels play in bilaterian synapses. However, a homologous role for CaV2 channels in presynaptic exocytosis has not been established for cnidarians or ctenophores36,108, although fast chemical synaptic transmission in cnidarians requires Ca2+ influx109,110. Indeed, while there is general agreement that cnidarian and bilaterian synapses are homologous, ctenophores are proposed to have evolved them independently60. Alternatively, synapses evolved once at the base of animals and were lost in placozoans and poriferans (sponges)111. For considering these alternate hypotheses, resolving the molecular mechanisms for synaptic trasnmission in cnidarians and ctenophores, including roles for CaV2 channels, will be important. Like CaV2 channels, it is unknown whether Mint plays homologous roles in non-bilaterian neurons as it does in bilaterians, where it is a major player in presynaptic function57,112 in part through its selective trafficking of axonal and synaptic proteins105,113, 114. In mice, knockout of the two neuronal Mint genes, Mint-1 and Mint-2, causes lethality in 80% of newborns while the surviving 20% of animals exhibit severe pre-synaptic defects57. Neuronal silencing of the two Mint isotypes in D. melanogaster, X11α and X11β, also causes larval lethality113, but not so in C. elegans114, indicating that while Mint is essential in some animals, it appears to play less essential roles in others.

Prompted by our yeast 2-hybrid screens showing that T. adhaerens Mint PDZ-1/PDZ-2/P1BM binds several proteins with C-terminal signatures similar to the T. adhaerens CaV2 channel sequence of SKCTAVCOOH (Fig. 1A and Table 1), we conducted directed yeast 2-hybrid experiments revealing that the tandem PDZ module could directly bind the CaV2 C-terminal ligand of CTAVCOOH in vitro (Fig. 9). This finding was consistent with our bacterial 2-hybrid experiments revealing a strong preference for threonine and valine at positions −2 and 0 (Fig. 10B). Interestingly, an interaction was not observed for PDZ-1 on its own, making this interaction different from mammals where PDZ-1 mediates binding to CaV2 channels46. However, PDZ-1 from T. adhaerens Mint did bind the C-terminus of N. vectensis CaV2c bearing the ligand sequence ETWCCOOH (Fig. 9). Furthermore, bacterial 2-hybrid screeing of this domain revealed a preference for tryptophan at position −1 similar to rat Mint-1 PDZ1 (Fig. 10A and B), both unexpected findings given that the sequence divergence between T. adhaerens and mammalian Mint PDZ-1 falls well below the cut-off that correlates similar ligand specificities between domains41 (Figs. 4 and S1). Altogether, our findings indicate that the placozoan CaV2 channel and Mint can interact in vitro, though not through PDZ-1, but rather involving the entire PDZ-1/PDZ-2/P1BM module (Fig. 9). This suggests that this non-canonical interaction is being mediated by PDZ-2, or perhaps by a transient locus formed by the association between P1BM and PDZ-1. Consistent with this latter notion is the observation that mutation of K920 to glutamine in PDZ-1, which should have destabilized the PDZ-1-P1BM interaction but minimally impacted PDZ-2 functionality, completely disrupted the PDZ-1/PDZ-2/P1BM interaction with CTAVCOOH, while maintaining PDZ-1 functionality as demonstrated by its retained ability to bind the N. vectensis ETWCCOOH ligand (Fig. 10). Also, the capacity of Mint to interact with CTAVCOOH peptides was also apparent for the cnidarian and ctenophore orthologues, suggesting this mode of binding is conserved. A difference however is that both the cnidarian and ctenophore Mint PDZ-1/PDZ-2/P1BM modules could also bind ETWCCOOH peptides, which was not the case for T. adhaerns Mint. Although we did not determine the locus for binding of ETWCCOOH to cnidarians or ctenophore Mint, our bacterial 2-hybrid data for N. vectensis indicates that PDZ-1/PDZ-2/P1BM has a preference for tryptophan at position −1, consistent with weak auto-inhibition of PDZ-1 by P1BM.

As noted, placozoans are the most early-diverging animals to possess all three types of metazoan CaV channel types, CaV1 to CaV3. Functional characterization of these channels revealed general similarities with mammalian and bilaterian orthologues, including the divergent calmodulin feedback regulation that distinguishes CaV1 and CaV2 channels from each other20,30,93. Unknown however is whether these channels serve similar cellular functions as they do in bilaterians, where for example CaV1 channels are the primary drivers of muscle contraction and CaV2 channels of neurotransmitter release at fast synapses36,108. Furthermore, like cnidarians and bilaterians, the physiological significance of the placozoan CaV2-Mint interaction characterized in this study is unknown. Of note, a recent single cell transcriptome study of four placozoan species, T. adhaerens, the closely-related species Trichoplax species H2, and two more distantly-related species Hoilungia hongkongensis and Cladtertia collaboinventa, revealed conserved expression patterns of CaV channels and Mint, with all three CaV channel types and the single accessory CaVβ subunit of CaV1 and CaV2 channels co-expressed in “peptidergic cells”, a cell type proposed to resemble an evolutionary precursor of neurons and neuroendcorine cells115. In the two placozoan species for which CaV2 expression data is available, Trichoplax sp. 2 and C. collaboinventa, CaV2 expression is mostly restricted to peptidergic cells, mirroring the enriched expression of CaV2 channels in the human brain116 and CaV2c in N. vectensis neurons (Figs. 7 and 8). Instead, both the single cell transcriptome data and immunolabeling with polyclonal antibodies against T. adhaerens CaV120 revealed broader expression including dorsal epithelial cells which are dynamically contractile and require Ca2+ influx to do so117. Although the homology of these cells to muscle is unknown, it is notable that broader expression of CaV1 channels oustide of neurons, including muscle, is similarly observed in humans116. Lastly, the single cell transcriptome data indicate that Mint and CaV2 are co-expressed in peptidergic cells115, indicating that the in vitro interaction identified here at least have the opportunity to interact in vivo, which by extension, provides a putative functional and physiological link between these two proteins.

Our most unexpected finding was that while none of the examined Mint proteins could bind the ctenophore (M. leydyi) CaV2 channel C-terminus of SGDICOOH (Figs. 1A and 9), H. californiensis Mint (PDZ-1/PDZ-2/P1BM) could bind both the N. vectensis ETWCCOOH and T. adhaerenes CTAVCOOH peptide ligands with high affinity (Fig. 9). This is especially notable given the extreme sequence divergence of both PDZ-1 and PDZ-2 of this species (Fig. 4C and D, Figure S1). Indeed, the combined ability of N. vectensis and T. adhaerens PDZ-1/PDZ-2/P1BM to also bind the non-canonical T. adhaerens CaV2 peptide of CTAVCOOH suggests this is a conserved and homologous interaction. Altogether, our various observations permit the following three hypotheses about how the interaction between Mint and CaV2 channels evolved (Fig. 11): (1) CaV2 channels in cnidarians/bilaterians and placozoans independently evolved the capacity to interact with Mint via emergence of distinct PDZ ligand motifs. This would be consistent with a large-scale analysis indicating that roughly one third of PDZ-ligand pairs observed in humans evolved through emergent C-terminal mutations creating nascent ligands recognized by pre-existing PDZ domains118. (2) The placozoan/cnidarian/bilaterian ancestor had a classical type I PDZ ligand of CTAVCOOH (i.e., X–S/T-X-ΦCOOH; where Φ denotes a hydrophobe)119, and subsequent sequence changes occurred in cnidarians/bilaterians produce rare, non-canonical, DDWC-like ligands (i.e., outside of the 16 specificity classes defined by large-scale analysis41). Or (3) the lowest common ancestor of placozoans/cnidarians/bilaterians had a DDWC-like ligand, and placozoans uniquely evolved a type I CTAVCOOH ligand.Fig. 11 Putative phylogeny of Mint and CaV2 channels in metazoans and choanoflagellates. Animal phylogeny under the ctenophore first and single origin of the nervous system hypotheses, illustrating the known presence/absence of Mint and CaV channels in different animal phyla and choanoflagellates. The proposed model for CaV channel evolution is based on previous phylogenetic studies2–4. Putative evolutionary events including gains/losses of genes and synapses are represented by symbols on corresponding branches of the tree, while the summary of presence/absence of genes is illustrated bedside the name of each organismal lineage. Under this scenario, CaV1 and CaV2 channels evolved via gene duplication of a CaV1/2 type channel near the stem metazoan lineage2, along with Mint. This was followed by losses of CaV1/2 channel genes is all animals except poriferans, losses of CaV1 to CaV3 channels in poriferans, and losses of CaV1 and CaV3 channels in ctenophores. The green numbers are used to illustrate three alternative hypotheses for evolution of the Mint-CaV2 interaction as described in the text.

Lastly, it is important to consider possible evolutionary interactions between Mint and RIM, since both are known to bind CaV2 channel DDWCCOOH and related ligands. As noted, while the physiological relevance of the CaV2-Mint interaction is unknown, RIM plays an essential and conserved role in tethering CaV2 channels at the presynaptic active zone in vertebrates and invertebrates16,42–44. Therefore, a question that emerges is whether the selective pressure to retain DDWC-like SLiMs in the otherwise highly divergent C-termini of cnidarian/bilaterian CaV2 channels came primarily from interactions with RIM, or did the CaV2-Mint interaction also contribute? Without understanding whether the Mint-CaV2 interaction serves some functional purpose in vivo, it is not possible to address this question. However, it is notable that in placozoans, an interaction between RIM and CaV2 channels is not expected, since the CaV2 channels have distinct C-terminal motifs that are inconsistent with the known specificity of the RIM PDZ domain. Moreover, the placozoan RIM homologues lack a PDZ domain. Placozoans do possess a second RIM gene which possesses a PDZ domain, referred to as type II RIM4. However, these genes are phylogenetically distinct from canonical type I RIMs that interact with CaV2 channels in synapses, and bear highly divergent PDZ sequences making it highly unlikely that they share common ligand preferences with type I RIMs4. On the other hand, our study revealed that placozoan CaV2 channels and Mint have the capacity to interact in vitro, extending this capability deep into the animal phylogeny, which at least implies that this could be an ancestral interaction with physiological purpose. Lastly, it is noteworthy that none of the examined Mint proteins could bind the ctenophore CaV2 channel C-terminus in vitro. Thus, based on our assays, ctenophores seem to distinguish themselves from other animals by lacking the capacity for Mint and CaV2 to interact. Ctenophores also distinguish themselves with respect to RIM, being the only animals with synapses that lack canonical type I RIM, only possessing the divergent type II RIM4.

Materials and methods

Phylogenetic analyses

Protein sequences used in this study were identified via protein BLAST searches of the NBCI database with CaV2, Mint, CASK, Veli, APP, presenilin, and neurexin orthologues from humans, C. elegans, and D. melanogaster as queries. The exception are proteins sequences from T. adhaerens, M. leidyi, and H. californiensis, which were obtained from mRNA transcriptomes assembled in house120, and H. hongkongensis, obtained from a published transcriptome121. All candidate sequences which produced a BLAST expect value less than 1E−6 were validated as true homologues via reciprocal BLAST of the NCBI database, as well as SmartBLAST122. In addition, InterProScan61 was used to predict and confirm the presence of domains and transmembrane helices expected for the various orthologous proteins (e.g., PDZ domains, SH3 domains, etc.). All extracted and validated sequences used for the presented phylogenetic analyses are provided in separate FASTA files as Supplemental Files (Information) 1 (CaV2), 2 (Mint), 3 (CASK), 4 (Veli), 5 (APP), 6 (presenilin) and 7 (neurexin). For consistency, we sought to identify sequences from a finite set of representative species with high quality gene datasets (i.e., possessing full length protein sequences for all/most extracted genes), spanning some of the major animal groupings. In cases where an identified sequence from a given species was too fragmented and missing key domains, an orthologue from a closely-related species was extracted using the same strategy detailed above.

For phylogenetic inference, homologous protein sequences were first aligned with MUSCLE123, then trimmed with trimAl124 using gap thresholds ranging from 0.6 to 0.95. The program IQ-TREE 2125 was then used to predict Bayesian Information Criterion (BIC) best-fit models and infer maximum-likelihood phylogenetic trees, with 1000 ultrafast bootstrap replicates to generate node support values. The various phylogenetic trees inferred from different trimmed alignments were then analyzed to select a single tree for each protein type having the highest node support values. All parameters used for phylogenetic inference, including the selected gap thresholds (GT), are as follows: CaV2: GT = 0.90, Q.pfam + R7; Mint: GT = 0.90, Q.insect + F + R6; CASK: GT = 0.70, Q.insect + R4; Veli: GT = 0.90, Q.insect + I + G4; APP: GT = 0.80, LG + I + R5; presenilin: GT = 0.60, Q.pfam + R5; and neurexin: GT = 0.90, WAG + F + I + R5. All aligned and trimmed protein sequences are provided in FASTA format within Supplemental Information 8 (CaV2), 9 (Mint), 10 (CASK), 11 (Veli), 12 (APP), 13 (presenilin) and 14 (neurexin).

Protein alignments and structural modeling

All protein alignments were prepared using MUSCLE and visualized with Jalview version 2126. For the alignment similarity plots of CaV2 channels and Mint (Fig. 1A and C, respectively), single homologous/paralogous sequences from representative species from each phylum were selected and aligned with MUSCLE, then analyzed with EMBOSS Plotcon127 using a running windows of 11 amino acids and the EBLOSUM62 comparison matrix. The aligned sequences are provided in FASTA format within Supplemental Information 15 (CaV2) and 16 (Mint). Predictions of MINT protein structure for orthologues from N. Vectensis (NCBI accession number XP_032233970.1, amino acids 1095 to 1306), T. adhaerens (evg1036658, amino acids 823 to 1034; Supplemental Information 2), H. californiensis (evg127856, amino acids 337 to 547; Supplemental Information 2), R. norvegicus (Mint-1, NCBI accession number NP_113967.2, amino acids 630 to 841) comprised of the tandem PDZ-1 and PDZ-2 domains and the P1BM element were generated with Alphafold282 on Compute Canada128. The relaxed structural predictions for each homologue bearing the highest per-residue confidence score (pLDDT) were selected and visualized using UCSF ChimeraX129. Modelling of the T. adhaerens Mint PDZ-1 bound to the ETWCCOOH and CTAVCOOH sequences was achieved by predicting the structure of the PDZ-1/PDZ-2/P1BM sequence with Phyre2130 in the intensive mode (100% of residues modelled at > 90% confidence), positioning the peptide based on a peptide-bound structure of the Erbin PDZ domain (PDB ID: 6Q), followed by alteration of the four C-terminal amino acids of the auto-inhibitory sequence (i.e., PEYICOOH) to that of the two peptides using Coot131, then energy minimization with the Chiron Rapid Energy Minimization Server132. The resulting structures were then analyzed and visualized with UCSF ChimeraX.

Fluorescence in situ hybridization

About one hundred seven-day old N. vectensis polyps at the tentacle bud stage were transferred into a tube with 1:2 diluted artificial sea water (Nematostella medium). To anesthetize/immobilize the polyps, 6.5% MgCl2 w/v was added to the water with animals. After one minute spinning at 500 RCF, Nematostella medium was replaced with Ca and Mg free seawater (CMFSW) with the following components: NaCl (449 mM), KCl (9 mM), Na2SO4 (33 mM), and NaHCO3 (2.5 mM). Polyps were washed in CMFSW three more times with spinning in between each wash. On the final wash, the polyps were macerated by vigorous pipetting until no big tissue chunks were visible. The dissociated cells were then centrifuged for 5 min at 500 RCF and CMFSW was replaced back with Nematostella medium. The pellet was then resuspended in Nematostella medium and transferred onto cover slips where cells were left to adhere to the cover slips for 1–2 h. Nematostella medium was then removed from the cover slips and these were quickly submerged in 4% PFA in phosphate buffered saline (PBS) for one hour. Fixed cells adhered to cover slips were then rinsed in PBS and dehydrated through an ascending methanol concentration series (30%, 60%, and 100% v/v). Samples were then stored at −20 °C in 100% methanol until labelling.

Before labelling, the samples were transferred to ethanol, rinsed one more time with ethanol, and rehydrated in a descending series of ethanol concentrations (100%, 70%, and 50% v/v) in PBS. In situ hybridization was then performed with probes and reagents from Advanced Cell Diagnostics, Inc. (Hayward, USA) according to the manufacturer’s instruction. The cells were then pretreated with Protease III diluted 1:15 in PBS before probe application. Sequences of N. vectensis gene orthologues were retrieved from NvERTx database85 and used to generate Z-probes (Nve-Cav2C, #852861; Nve-MINT, #852881-C3; Nve-Otx-C2, #1006821-C2; and Nve-FoxL2-O1-C2, #1007261), which were designed by ACD using the proprietary RNAscope Probe Design pipeline. Special care was taken during the design of N. vectensis CaV2c so that the probes do not hybridize with CaV2a and CaV2b transcripts. The negative control included application of 3-Plex Negative Control Probe (#320871), which gave no labelling.

Samples were counterstained with DAPI, mounted in ProLong™ Gold antifade reagent (# P36934, Invitrogen, USA), and examined in LSM 800 confocal microscope with a 63X 1.4 NA objective and fluorescence and DIC optics. Labelling with each probe was done at least twice; the results of independently repeated experiments were similar. For each preparation several randomly chosen fields were scanned. Cell counts were then made for each field, and for each cell in the field we evaluated whether it expressed one gene, combination of any two genes, all three genes, or none of the studied genes. A cell was considered labeled if it possessed at least one fluorescent grain. Mean percent and standard deviation of expressing cells was obtained by averaging the data from all fields. To demonstrate the nature of overlapping expression patterns we drew accurate area-proportional Venn diagrams using eulerAPE software133. The fluorescence images were processed with the Fiji image analysis software134.

Co-immunoprecipitation experiments

Co-immunoprecipitation experiments were confirmed in triplicate sets of experiments using recombinant proteins transiently expressed in HEK-293T cells (MilliporeSigma, Canada). Briefly, cDNAs corresponding to the C-terminal portions of the N. vectensis and T. adhaerens Mint PDZ-1/PDZ-2/P1BM protein sequences (as described in the results section) were cloned into the mammalian bicistronic expression vector pIRES2-DsRed2 (Clontech, USA) to express recombinant proteins bearing C-terminal myc epitope tags (Fig. 6K). The N. vectensis sequence was gene synthesized with codon optimization for mammalian cells, whereas the T. adhaerens sequence was PCR amplified from a whole animal cDNA library (all primers and template details for cloning constructs used in these experiments are listed in Supplemental Table (Information) 17). Similarly, the entire C-termini of N. vectensis CaV2c and T. adhaerens CaV2 were PCR amplified from corresponding pIRES2-EGFP plasmid constructs bearing the complete protein coding sequences of the two channels, the former generated via gene synthesis (GenScript), the latter cloned directly from a cDNA library30, and subsequently restriction digest cloned into the vector pEGFP-C1 (Clontech) to express recombinant fusion proteins bearing an N-terminal Enhanced Green Fluorescent Protein (EGFP) followed by an HA epitope tag (Fig. 6J). ΔPDZ ligand variants were created via site directed mutagenesis to introduce premature stop codon upstream of the original stop codon, resulting in the deletion of the 9 most distal amino acids of each CaV2 channel C-terminus. HEK-293T cells were cultured in 25 cm2 vented flasks at 37°C in a humidified 5% CO2 incubator in 6 mL of DMEM supplemented with 10% FBS and 0.5 mL of penicillin/streptomycin (Sigma-Aldrich, Canada; 5000 units/mL penicillin and 5 mg/mL streptomycin). For transfections, cells were grown to 60% confluency then transfected with 3 μg of each relevant DNA construct using the PolyJet™ transfection reagent (SignaGen Laboratories, USA). 6 h later, cells were washed with warm 1 × Phosphate Buffered Saline (PBS) (Gibco, ThermoScientific Scientific, Canada) and placed at 37 °C for 24 to 36 h. Cells were then transferred to a 28 °C 5% CO2 incubator for 2 to 3 days to increase protein expression135, then harvested by first washing with 1 × PBS, then solubilizing in extraction buffer containing 0.5% Triton X-100, 20 mM imidazole at pH 6.8, 100 mM NaCl, 1 mM EDTA, 1 mM DDT, and 1 × Halt™ Protease Inhibitor Cocktail (ThermoFisher Scientific). Lysates were centrifuged at 16,000×g for 15 min and treated with 1 unit/μL Benzonase Nuclease (E1014, Millipore Sigma, Canada) for 30 min on ice. Lysates were then then quantified for total protein content using the Pierce™ BCA Protein Assay Kit (ThermoFisher Scientific) then either snap frozen in liquid nitrogen and stored at −80 °C or immediately used in downstream experiments.

Co-immunoprecipitation experiments were performed using the Pierce™ Crosslink IP Kit (ThermoFisher Scientific) with minor modifications to the manufacturer’s protocol. Briefly, the included Pierce™ Protein A/G Plus agarose beads were incubated with rabbit anti-myc polyclonal antibody (ab9106, Abcam, USA) and cross-linked using 2.5 mM DSS crosslinker. Protein lysates were then pre-cleared of non-specificities using the Pierce™ Control Agarose Resin and incubated with the cross-linked Pierce™ Protein A/G Plus Agarose beads. The beads were then washed with extraction buffer multiple times and immunoprecipitated complexes were boiled off in NuPAGE™ 4 × Sample Loading Buffer (ThermoFisher Scientific, Canada). The co-immunoprecipitated and immunoprecipitated proteins were then characterized by Western Blotting using anti-HA monoclonal antibody (ab18181, Abcam, USA) and anti-myc monoclonal antibody (ab18185, Abcam, USA), respectively, both at a dilution of 1:1000. Input was 1% of the amount of cell lysate used for co-immunoprecipitation and was similarly characterized by Western blotting with anti-myc and anti-HA monoclonal antibodies. Chemiluminescent signals on blots were detected using SuperSignal West Femto Maximum Sensitivity Substrate (Millipore Sigma) imaged an ImageQuant LAS 500 Chemiluminescence CDD Imager (GE Life Sciences, USA). All raw images of Western blots are provided in Supplemental Information 18.

Yeast 2-hybrid and bacterial 2-hybrid experiments

The T. adhaerens yeast 2-hybrid cDNA library was constructed using the Make Your Own “Mate & Plate™” Library System (Takara Bio USA, Mountain View, CA), using whole animal total RNA extracted from approximately 30 animals using a Nucleospin RNA Plus Mini Kit (Macherey–Nagel, Düren, Germany), both according to manufacturer’s instructions. The resulting yeast 2-hybrid library was screened according to the manufacturer’s instructions, using the Mint PDZ-1/PDZ-2, P1BM and PDZ-1 proteins sequences as bait cloned expressed from the bait vector pGBKT7 as fusion proteins N-terminally tagged with the Gal4 DNA binding domain. All primers used for PCR amplification and restriction digest cloning, as well as template details, are provided in Supplemental Information 17. For the directed yeast 2-hybrid screens, the 70 most distal amino acids of T. adhaerens CaV2 and M. leidyi CaV2 were PCR amplified from previously cloned pIRES2-EGFP plasmid constructs containing the full-length protein coding sequences of each channel, then cloned into the prey vector pGADT7 with and without the last four amino acids of the channel C-termini to express fusion proteins with the Gal4 activation domain at the N-terminus. The C-terminus of N. vectensis CaV2c was similarly cloned, using a gene synthesized coding sequence (GenScript, USA) as template for the PCR reaction. All yeast 2-hybrid screens were done using the Matchmaker® Gold Yeast Two-Hybrid System (Takara Bio USA) according to the manufacturer’s instructions, including necessary controls to rule out auto-activation of protein–protein interaction reporter genes by the various Mint insert proteins. Measurements of β-galactosidase reporter gene activity were made with the ThermoScientific Yeast β-Galactosidase Assay Kit according to the manufacturer’s instructions, on flat bottom 96-well plates read using a BioTek Epoch Microplate Spectrophotometer (Agilent Technologies Canada Inc., Canada) to determine cell densities (OD600) and colorimetric enzyme activity (absorbance at 420 nm).

Bacterial 2-hybrid experiments

Bacterial 2-hybrid assays were carried out as previously described94 with slight modifications to the assay. First, we modified the peptide length of the library from 7 amino acids to 6, finding that in general screens of PDZ-peptide specificity rarely enrich any preference at position −6. Reducing the peptide length in this way reduced the complexity of the library, which permitted us to sample potential ligand sequences in greater depth. We used the same library building approach as previously described to minimize bias in the library136. Second, we carried out all selections on agar plates containing 5 mM 3-aminotriazole (3AT), an inhibitor of our HIS3 reporter which minimizes the enrichment of false positives. We used plates rather than liquid culture as we found similar data was produced with significantly less effort. Once grown, the surviving members of the library were pooled, sequenced by Illumina, and motifs generated as previously described94.

Supplementary Information

Supplementary Information 1.

Supplementary Information 2.

Supplementary Information 3.

Supplementary Information 4.

Supplementary Information 5.

Supplementary Information 6.

Supplementary Information 7.

Supplementary Information 8.

Supplementary Information 9.

Supplementary Information 10.

Supplementary Information 11.

Supplementary Information 12.

Supplementary Information 13.

Supplementary Information 14.

Supplementary Information 15.

Supplementary Information 16.

Supplementary Information 17.

Supplementary Information 18.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70652-8.

Acknowledgements

We would like to thank the Light Imaging Facility and Drs. Carolyn Smith and Thomas Reese of the National Institute of Neurological Disorders and Stroke, Drs. Kristen Koenig and Kyle McCulloch for providing the N. vectensis specimens used in this study, Dr. Elise F. Stanley for her feedback during early stages of this study, and Dr. Joseph Ryan for his feedback on the manuscript.

Author contributions

A.Senatore and A.N.S. conceived the project and conducted the phylogenetic and sequence comparison analyses. A.Senatore, G.M., and M.A.C. conducted the in silico structural analyses. T.D.M. conducted the in situ hybridization experiments. A.Senatore conducted the scRNA-Seq analysis. J.H., J.G. and M.N. conducted the bacterial 2-hybrid experiments. A.Singh, B.B., I.M.I., and T.E. conducted the yeast 2-hybrid experiments. A.N.H and A.T. conducted molecular biology and co-immunoprecipitation experiments. A.Senatore and A.N.H. wrote the first draft of the manuscript and all authors contributed to its revision and finalization.

Funding

This research was funded by an NSERC Discovery Grant (RGPIN-2021-03557), an NSERC Discovery Accelerator Supplement (RGPAS-2021-00002), an Ontario Early Researcher Award (ER17-13-247), and a Canadian Foundation for Innovation Grant (35297) to A.Senatore, Ontario Graduate Scholarships to A.Singh, W.E., and G.M., an NSERC Canadian Graduate Scholarship to J.G., an NSERC Discovery Grant (RGPIN-2023-05615) and a Canadian Foundation for Innovation Grant (40684) to M.A.C., a NIH grant (R01GM133936) to M.N., and the Intramural Research Program of the National Institute of Neurological Disorders and Stroke which supported work done by T.D.M.

Data availability

All datasets of protein sequences used for phylogenetic and/or in silico structural analyses are included in the published article and its supplementary information files. All other datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Alicia N. Harracksingh, Anhadvir Singh, Tatiana D. Mayorova and Brian Bejoy.
==== Refs
References

1. Catterall WA Voltage-gated calcium channels Cold Spring Harb. Perspect. Biol. 2011 3 8 a003947 10.1101/cshperspect.a003947 21746798
Catterall, W. A. Voltage-gated calcium channels. Cold Spring Harb. Perspect. Biol. 3(8), a003947 (2011).21746798 10.1101/cshperspect.a003947
2. Moran Y Zakon HH The evolution of the four subunits of voltage-gated calcium channels: Ancient roots, increasing complexity, and multiple losses Genome Biol. Evol. 2014 6 9 2210 2217 10.1093/gbe/evu177 25146647
Moran, Y. & Zakon, H. H. The evolution of the four subunits of voltage-gated calcium channels: Ancient roots, increasing complexity, and multiple losses. Genome Biol. Evol. 6(9), 2210–2217 (2014).25146647 10.1093/gbe/evu177
3. Pozdnyakov I Matantseva O Skarlato S Diversity and evolution of four-domain voltage-gated cation channels of eukaryotes and their ancestral functional determinants Sci. Rep. 2018 8 1 3539 10.1038/s41598-018-21897-7 29476068
Pozdnyakov, I., Matantseva, O. & Skarlato, S. Diversity and evolution of four-domain voltage-gated cation channels of eukaryotes and their ancestral functional determinants. Sci. Rep. 8(1), 3539 (2018).29476068 10.1038/s41598-018-21897-7
4. Piekut T Early metazoan origin and multiple losses of a novel clade of rim presynaptic calcium channel scaffolding protein homologs Genome Biol. Evol. 2020 12 8 1217 1239 10.1093/gbe/evaa097 32413100
Piekut, T. et al. Early metazoan origin and multiple losses of a novel clade of rim presynaptic calcium channel scaffolding protein homologs. Genome Biol. Evol. 12(8), 1217–1239 (2020).32413100 10.1093/gbe/evaa097
5. Perez-Reyes E Molecular physiology of low-voltage-activated t-type calcium channels Physiol. Rev. 2003 83 1 117 161 10.1152/physrev.00018.2002 12506128
Perez-Reyes, E. Molecular physiology of low-voltage-activated t-type calcium channels. Physiol. Rev. 83(1), 117–161 (2003).12506128 10.1152/physrev.00018.2002
6. Senatore A T-type channels become highly permeable to sodium ions using an alternative extracellular turret region (S5-P) Outside the selectivity filter* J. Biol. Chem. 2014 289 17 11952 11969 10.1074/jbc.M114.551473 24596098
Senatore, A. et al. T-type channels become highly permeable to sodium ions using an alternative extracellular turret region (S5-P) Outside the selectivity filter*. J. Biol. Chem. 289(17), 11952–11969 (2014).24596098 10.1074/jbc.M114.551473
7. Hara Y Koganezawa M Yamamoto D The Dmca1D channel mediates Ca2+ inward currents in Drosophila embryonic muscles J. Neurogenet. 2015 29 2–3 117 123 10.3109/01677063.2015.1054991 26004544
Hara, Y., Koganezawa, M. & Yamamoto, D. The Dmca1D channel mediates Ca2+ inward currents in Drosophila embryonic muscles. J. Neurogenet. 29(2–3), 117–123 (2015).26004544 10.3109/01677063.2015.1054991
8. Collet C Excitation-contraction coupling in skeletal muscle fibers from adult domestic honeybee Pflügers Archiv-Eur. J. Physiol. 2009 458 601 612 10.1007/s00424-009-0642-6 19198873
Collet, C. Excitation-contraction coupling in skeletal muscle fibers from adult domestic honeybee. Pflügers Archiv-Eur. J. Physiol. 458, 601–612 (2009).19198873 10.1007/s00424-009-0642-6
9. Eberl DF Genetic and developmental characterization of Dmca1D, a calcium channel α1 subunit gene in Drosophila melanogaster Genetics 1998 148 3 1159 1169 10.1093/genetics/148.3.1159 9539432
Eberl, D. F. et al. Genetic and developmental characterization of Dmca1D, a calcium channel α1 subunit gene in Drosophila melanogaster. Genetics 148(3), 1159–1169 (1998).9539432 10.1093/genetics/148.3.1159
10. Erxleben C Rathmayer W A dihydropyridine-sensitive voltage-dependent calcium channel in the sarcolemmal membrane of crustacean muscle J. Gen. Physiol. 1997 109 3 313 326 10.1085/jgp.109.3.313 9089439
Erxleben, C. & Rathmayer, W. A dihydropyridine-sensitive voltage-dependent calcium channel in the sarcolemmal membrane of crustacean muscle. J. Gen. Physiol. 109(3), 313–326 (1997).9089439 10.1085/jgp.109.3.313
11. Jospin M The L-type voltage-dependent Ca2+ channel EGL-19 controls body wall muscle function in Caenorhabditis elegans J. Cell Biol. 2002 159 2 337 348 10.1083/jcb.200203055 12391025
Jospin, M. et al. The L-type voltage-dependent Ca2+ channel EGL-19 controls body wall muscle function in Caenorhabditis elegans. J. Cell Biol. 159(2), 337–348 (2002).12391025 10.1083/jcb.200203055
12. Di Biase V Franzini-Armstrong C Evolution of skeletal type e–c coupling: A novel means of controlling calcium delivery J. Cell Biol. 2005 171 4 695 704 10.1083/jcb.200503077 16286507
Di Biase, V. & Franzini-Armstrong, C. Evolution of skeletal type e–c coupling: A novel means of controlling calcium delivery. J. Cell Biol. 171(4), 695–704 (2005).16286507 10.1083/jcb.200503077
13. Chan JD Utilizing the planarian voltage-gated ion channel transcriptome to resolve a role for a Ca2+ channel in neuromuscular function and regeneration Biochim. Biophys. Acta (BBA) Mol. Cell Res. 2017 1864 6 1036 1045 10.1016/j.bbamcr.2016.10.010
Chan, J. D. et al. Utilizing the planarian voltage-gated ion channel transcriptome to resolve a role for a Ca2+ channel in neuromuscular function and regeneration. Biochim. Biophys. Acta (BBA) Mol. Cell Res. 1864(6), 1036–1045 (2017).10.1016/j.bbamcr.2016.10.010
14. Spafford JD Expression and modulation of an invertebrate presynaptic calcium channel α1 subunit homolog J. Biol. Chem. 2003 278 23 21178 21187 10.1074/jbc.M302212200 12672808
Spafford, J. D. et al. Expression and modulation of an invertebrate presynaptic calcium channel α1 subunit homolog. J. Biol. Chem. 278(23), 21178–21187 (2003).12672808 10.1074/jbc.M302212200
15. Spafford JD Calcium channel structural determinants of synaptic transmission between identified invertebrate neurons J. Biol. Chem. 2003 278 6 4258 4267 10.1074/jbc.M211076200 12458203
Spafford, J. D. et al. Calcium channel structural determinants of synaptic transmission between identified invertebrate neurons. J. Biol. Chem. 278(6), 4258–4267 (2003).12458203 10.1074/jbc.M211076200
16. Dunn, T. W. et al. The role of specific isoforms of CaV2 and the common C-terminal of CaV2 in calcium channel function in sensory neurons of Aplysia. Sci. Rep. 13(1), 20216 (2023).
17. Edmonds B Contributions of two types of calcium channels to synaptic transmission and plasticity Science 1990 250 4984 1142 1147 10.1126/science.2174573 2174573
Edmonds, B. et al. Contributions of two types of calcium channels to synaptic transmission and plasticity. Science 250(4984), 1142–1147 (1990).2174573 10.1126/science.2174573
18. Saheki Y Bargmann CI Presynaptic CaV2 calcium channel traffic requires CALF-1 and the α2δ subunit UNC-36 Nat. Neurosci. 2009 12 10 1257 1265 10.1038/nn.2383 19718034
Saheki, Y. & Bargmann, C. I. Presynaptic CaV2 calcium channel traffic requires CALF-1 and the α2δ subunit UNC-36. Nat. Neurosci. 12(10), 1257–1265 (2009).19718034 10.1038/nn.2383
19. Kawasaki F Active zone localization of presynaptic calcium channels encoded by the cacophony locus of Drosophila J. Neurosci. 2004 24 1 282 285 10.1523/JNEUROSCI.3553-03.2004 14715960
Kawasaki, F. et al. Active zone localization of presynaptic calcium channels encoded by the cacophony locus of Drosophila. J. Neurosci. 24(1), 282–285 (2004).14715960 10.1523/JNEUROSCI.3553-03.2004
20. Gauberg J Divergent Ca2+/calmodulin feedback regulation of CaV1 and CaV2 voltage-gated calcium channels evolved in the common ancestor of Placozoa and Bilateria J. Biol. Chem. 2022 298 4 101741 10.1016/j.jbc.2022.101741 35182524
Gauberg, J. et al. Divergent Ca2+/calmodulin feedback regulation of CaV1 and CaV2 voltage-gated calcium channels evolved in the common ancestor of Placozoa and Bilateria. J. Biol. Chem. 298(4), 101741 (2022).35182524 10.1016/j.jbc.2022.101741
21. Ben-Johny M Yue DT Calmodulin regulation (calmodulation) of voltage-gated calcium channels J. Gen. Physiol. 2014 143 6 679 692 10.1085/jgp.201311153 24863929
Ben-Johny, M. & Yue, D. T. Calmodulin regulation (calmodulation) of voltage-gated calcium channels. J. Gen. Physiol. 143(6), 679–692 (2014).24863929 10.1085/jgp.201311153
22. Taiakina V The calmodulin-binding, short linear motif, NSCaTE is conserved in L-type channel ancestors of vertebrate Cav1.2 and Cav1.3 channels PLoS ONE 2013 8 4 e61765 10.1371/journal.pone.0061765 23626724
Taiakina, V. et al. The calmodulin-binding, short linear motif, NSCaTE is conserved in L-type channel ancestors of vertebrate Cav1.2 and Cav1.3 channels. PLoS ONE 8(4), e61765 (2013).23626724 10.1371/journal.pone.0061765
23. Zamponi GW An MT Zamponi GW Voltage-Dependent Inactivation of Voltage Gated Calcium Channels. Voltage-Gated Calcium Channels 2005 Springer 194 204
Zamponi, G. W., An, M. T. & Zamponi, G. W. Voltage-Dependent Inactivation of Voltage Gated Calcium Channels. Voltage-Gated Calcium Channels 194–204 (Springer, 2005).
24. Cens T Voltage-and calcium-dependent inactivation in high voltage-gated Ca2+ channels Prog. Biophys. Mol. Biol. 2006 90 1–3 104 117 10.1016/j.pbiomolbio.2005.05.013 16038964
Cens, T. et al. Voltage-and calcium-dependent inactivation in high voltage-gated Ca2+ channels. Prog. Biophys. Mol. Biol. 90(1–3), 104–117 (2006).16038964 10.1016/j.pbiomolbio.2005.05.013
25. Tedford HW Zamponi GW Direct G protein modulation of Cav2 calcium channels Pharmacol. Rev. 2006 58 4 837 862 10.1124/pr.58.4.11 17132857
Tedford, H. W. & Zamponi, G. W. Direct G protein modulation of Cav2 calcium channels. Pharmacol. Rev. 58(4), 837–862 (2006).17132857 10.1124/pr.58.4.11
26. Herlitze S Modulation of Ca2+ channels by G-protein beta gamma subunits Nature 1996 380 6571 258 262 10.1038/380258a0 8637576
Herlitze, S. et al. Modulation of Ca2+ channels by G-protein beta gamma subunits. Nature 380(6571), 258–262 (1996).8637576 10.1038/380258a0
27. Huang X G-proteins modulate invertebrate synaptic calcium channel (LCav2) differently from the classical voltage-dependent regulation of mammalian Cav2.1 and Cav2.2 channels J. Exp. Biol. 2010 213 12 2094 2103 10.1242/jeb.042242 20511524
Huang, X. et al. G-proteins modulate invertebrate synaptic calcium channel (LCav2) differently from the classical voltage-dependent regulation of mammalian Cav2.1 and Cav2.2 channels. J. Exp. Biol. 213(12), 2094–2103 (2010).20511524 10.1242/jeb.042242
28. de Hoog E Lukewich MK Spencer GE Retinoid receptor-based signaling plays a role in voltage-dependent inhibition of invertebrate voltage-gated Ca2+ channels J. Biol. Chem. 2019 294 26 10076 10093 10.1074/jbc.RA118.006444 31048374
de Hoog, E., Lukewich, M. K. & Spencer, G. E. Retinoid receptor-based signaling plays a role in voltage-dependent inhibition of invertebrate voltage-gated Ca2+ channels. J. Biol. Chem. 294(26), 10076–10093 (2019).31048374 10.1074/jbc.RA118.006444
29. Dong N Dopamine-mediated calcium channel regulation in synaptic suppression in L. stagnalis interneurons Channels 2018 12 1 153 173 10.1080/19336950.2018.1457897 29589519
Dong, N. et al. Dopamine-mediated calcium channel regulation in synaptic suppression in L. stagnalis interneurons. Channels 12(1), 153–173 (2018).29589519 10.1080/19336950.2018.1457897
30. Gauberg J Conserved biophysical features of the CaV2 presynaptic Ca2+ channel homologue from the early-diverging animal Trichoplax adhaerens J. Biol. Chem. 2020 295 52 18553 18578 10.1074/jbc.RA120.015725 33097592
Gauberg, J. et al. Conserved biophysical features of the CaV2 presynaptic Ca2+ channel homologue from the early-diverging animal Trichoplax adhaerens. J. Biol. Chem. 295(52), 18553–18578 (2020).33097592 10.1074/jbc.RA120.015725
31. Clapham DE Calcium signaling Cell 2007 131 6 1047 1058 10.1016/j.cell.2007.11.028 18083096
Clapham, D. E. Calcium signaling. Cell 131(6), 1047–1058 (2007).18083096 10.1016/j.cell.2007.11.028
32. Naraghi M Neher E Linearized buffered Ca2+ diffusion in microdomains and its implications for calculation of [Ca2+] at the mouth of a calcium channel J. Neurosci. 1997 17 18 6961 6973 10.1523/JNEUROSCI.17-18-06961.1997 9278532
Naraghi, M. & Neher, E. Linearized buffered Ca2+ diffusion in microdomains and its implications for calculation of [Ca2+] at the mouth of a calcium channel. J. Neurosci. 17(18), 6961–6973 (1997).9278532 10.1523/JNEUROSCI.17-18-06961.1997
33. Parekh AB Ca2+ microdomains near plasma membrane Ca2+ channels: Impact on cell function J. Physiol. 2008 586 13 3043 3054 10.1113/jphysiol.2008.153460 18467365
Parekh, A. B. Ca2+ microdomains near plasma membrane Ca2+ channels: Impact on cell function. J. Physiol. 586(13), 3043–3054 (2008).18467365 10.1113/jphysiol.2008.153460
34. Berridge MJ Bootman MD Roderick HL Calcium signalling: Dynamics, homeostasis and remodelling Nat. Rev. Mol. Cell Biol. 2003 4 7 517 529 10.1038/nrm1155 12838335
Berridge, M. J., Bootman, M. D. & Roderick, H. L. Calcium signalling: Dynamics, homeostasis and remodelling. Nat. Rev. Mol. Cell Biol. 4(7), 517–529 (2003).12838335 10.1038/nrm1155
35. Stanley EF The nanophysiology of fast transmitter release Trends Neurosci. 2016 39 3 183 197 10.1016/j.tins.2016.01.005 26896416
Stanley, E. F. The nanophysiology of fast transmitter release. Trends Neurosci. 39(3), 183–197 (2016).26896416 10.1016/j.tins.2016.01.005
36. Senatore A Spafford JD Voltage-Gated Calcium Channels in Invertebrates Voltage-Gated Calcium Channels 2022 Springer 115 158
Senatore, A. & Spafford, J. D. Voltage-Gated Calcium Channels in Invertebrates. In Voltage-Gated Calcium Channels 115–158 (Springer, 2022).
37. Wang S Cortes CJ Interactions with PDZ proteins diversify voltage-gated calcium channel signaling J. Neurosci. Res. 2021 99 1 332 348 10.1002/jnr.24650 32476168
Wang, S. & Cortes, C. J. Interactions with PDZ proteins diversify voltage-gated calcium channel signaling. J. Neurosci. Res. 99(1), 332–348 (2021).32476168 10.1002/jnr.24650
38. Zhang H Association of CaV1.3 L-type calcium channels with Shank J. Neurosci. 2005 25 5 1037 1049 10.1523/JNEUROSCI.4554-04.2005 15689539
Zhang, H. et al. Association of CaV1.3 L-type calcium channels with Shank. J. Neurosci. 25(5), 1037–1049 (2005).15689539 10.1523/JNEUROSCI.4554-04.2005
39. Pym E Shank is a dose-dependent regulator of Cav1 calcium current and CREB target expression Elife 2017 6 e18931 10.7554/eLife.18931 28477407
Pym, E. et al. Shank is a dose-dependent regulator of Cav1 calcium current and CREB target expression. Elife 6, e18931 (2017).28477407 10.7554/eLife.18931
40. Gao L Shank promotes action potential repolarization by recruiting BK channels to calcium microdomains Elife 2022 11 e75140 10.7554/eLife.75140 35266450
Gao, L. et al. Shank promotes action potential repolarization by recruiting BK channels to calcium microdomains. Elife 11, e75140 (2022).35266450 10.7554/eLife.75140
41. Tonikian R A specificity map for the PDZ domain family PLoS Biol. 2008 6 9 e239 10.1371/journal.pbio.0060239 18828675
Tonikian, R. et al. A specificity map for the PDZ domain family. PLoS Biol. 6(9), e239 (2008).18828675 10.1371/journal.pbio.0060239
42. Graf ER RIM promotes calcium channel accumulation at active zones of the Drosophila neuromuscular junction J. Neurosci. 2012 32 47 16586 16596 10.1523/JNEUROSCI.0965-12.2012 23175814
Graf, E. R. et al. RIM promotes calcium channel accumulation at active zones of the Drosophila neuromuscular junction. J. Neurosci. 32(47), 16586–16596 (2012).23175814 10.1523/JNEUROSCI.0965-12.2012
43. Kaeser PS RIM proteins tether Ca2+ channels to presynaptic active zones via a direct PDZ-domain interaction Cell 2011 144 2 282 295 10.1016/j.cell.2010.12.029 21241895
Kaeser, P. S. et al. RIM proteins tether Ca2+ channels to presynaptic active zones via a direct PDZ-domain interaction. Cell 144(2), 282–295 (2011).21241895 10.1016/j.cell.2010.12.029
44. Kushibiki Y RIMB-1/RIM-binding protein and UNC-10/RIM redundantly regulate presynaptic localization of the voltage-gated calcium channel in Caenorhabditis elegans J. Neurosci. 2019 39 44 8617 8631 10.1523/JNEUROSCI.0506-19.2019 31530643
Kushibiki, Y. et al. RIMB-1/RIM-binding protein and UNC-10/RIM redundantly regulate presynaptic localization of the voltage-gated calcium channel in Caenorhabditis elegans. J. Neurosci. 39(44), 8617–8631 (2019).31530643 10.1523/JNEUROSCI.0506-19.2019
45. Biederer T Südhof TC Mints as adaptors: Direct binding to neurexins and recruitment of munc18 J. Biol. Chem. 2000 275 51 39803 39806 10.1074/jbc.C000656200 11036064
Biederer, T. & Südhof, T. C. Mints as adaptors: Direct binding to neurexins and recruitment of munc18. J. Biol. Chem. 275(51), 39803–39806 (2000).11036064 10.1074/jbc.C000656200
46. Maximov A Südhof TC Bezprozvanny I Association of neuronal calcium channels with modular adaptor proteins J. Biol. Chem. 1999 274 35 24453 24456 10.1074/jbc.274.35.24453 10455105
Maximov, A., Südhof, T. C. & Bezprozvanny, I. Association of neuronal calcium channels with modular adaptor proteins. J. Biol. Chem. 274(35), 24453–24456 (1999).10455105 10.1074/jbc.274.35.24453
47. Bezprozvanny I Maximov A Classification of PDZ domains FEBS Lett. 2001 509 3 457 462 10.1016/S0014-5793(01)03214-8 11749973
Bezprozvanny, I. & Maximov, A. Classification of PDZ domains. FEBS Lett. 509(3), 457–462 (2001).11749973 10.1016/S0014-5793(01)03214-8
48. Long J-F Autoinhibition of X11/Mint scaffold proteins revealed by the closed conformation of the PDZ tandem Nat. Struct. Mol. Biol. 2005 12 8 722 728 10.1038/nsmb958 16007100
Long, J.-F. et al. Autoinhibition of X11/Mint scaffold proteins revealed by the closed conformation of the PDZ tandem. Nat. Struct. Mol. Biol. 12(8), 722–728 (2005).16007100 10.1038/nsmb958
49. Gardezi SR Li Q Stanley EF Inter-channel scaffolding of presynaptic CaV2.2 via the C terminal PDZ ligand domain Biol. Open 2013 2 5 492 498 10.1242/bio.20134267 23789098
Gardezi, S. R., Li, Q. & Stanley, E. F. Inter-channel scaffolding of presynaptic CaV2.2 via the C terminal PDZ ligand domain. Biol. Open 2(5), 492–498 (2013).23789098 10.1242/bio.20134267
50. Borg J-P The phosphotyrosine interaction domains of X11 and FE65 bind to distinct sites on the YENPTY motif of amyloid precursor protein Mol. Cell. Biol. 1996 16 11 6229 6241 10.1128/MCB.16.11.6229 8887653
Borg, J.-P. et al. The phosphotyrosine interaction domains of X11 and FE65 bind to distinct sites on the YENPTY motif of amyloid precursor protein. Mol. Cell. Biol. 16(11), 6229–6241 (1996).8887653 10.1128/MCB.16.11.6229
51. McLoughlin DM Miller CC The intracellular cytoplasmic domain of the Alzheimer's disease amyloid precursor protein interacts with phosphotyrosine-binding domain proteins in the yeast two-hybrid system FEBS Lett. 1996 397 2–3 197 200 10.1016/S0014-5793(96)01128-3 8955346
McLoughlin, D. M. & Miller, C. C. The intracellular cytoplasmic domain of the Alzheimer’s disease amyloid precursor protein interacts with phosphotyrosine-binding domain proteins in the yeast two-hybrid system. FEBS Lett. 397(2–3), 197–200 (1996).8955346 10.1016/S0014-5793(96)01128-3
52. Tanahashi H Tabira T X11L2, a new member of the X11 protein family, interacts with Alzheimer's β-amyloid precursor protein Biochem. Biophys. Res. Commun. 1999 255 3 663 667 10.1006/bbrc.1999.0265 10049767
Tanahashi, H. & Tabira, T. X11L2, a new member of the X11 protein family, interacts with Alzheimer’s β-amyloid precursor protein. Biochem. Biophys. Res. Commun. 255(3), 663–667 (1999).10049767 10.1006/bbrc.1999.0265
53. Okamoto M Amyloid precursor protein associates independently and collaboratively with PTB and PDZ domains of mint on vesicles and at cell membrane Neuroscience 2001 104 3 653 665 10.1016/S0306-4522(01)00124-5 11440799
Okamoto, M. et al. Amyloid precursor protein associates independently and collaboratively with PTB and PDZ domains of mint on vesicles and at cell membrane. Neuroscience 104(3), 653–665 (2001).11440799 10.1016/S0306-4522(01)00124-5
54. Biederer T Regulation of APP-dependent transcription complexes by Mint/X11s: Differential functions of Mint isoforms J. Neurosci. 2002 22 17 7340 7351 10.1523/JNEUROSCI.22-17-07340.2002 12196555
Biederer, T. et al. Regulation of APP-dependent transcription complexes by Mint/X11s: Differential functions of Mint isoforms. J. Neurosci. 22(17), 7340–7351 (2002).12196555 10.1523/JNEUROSCI.22-17-07340.2002
55. Lau K-F X11α and X11β interact with presenilin-1 via their PDZ domains Mol. Cell. Neurosci. 2000 16 5 557 565 10.1006/mcne.2000.0898 11083918
Lau, K.-F. et al. X11α and X11β interact with presenilin-1 via their PDZ domains. Mol. Cell. Neurosci. 16(5), 557–565 (2000).11083918 10.1006/mcne.2000.0898
56. Ho A Liu X Südhof TC Deletion of Mint proteins decreases amyloid production in transgenic mouse models of Alzheimer's disease J. Neurosci. 2008 28 53 14392 14400 10.1523/JNEUROSCI.2481-08.2008 19118172
Ho, A., Liu, X. & Südhof, T. C. Deletion of Mint proteins decreases amyloid production in transgenic mouse models of Alzheimer’s disease. J. Neurosci. 28(53), 14392–14400 (2008).19118172 10.1523/JNEUROSCI.2481-08.2008
57. Ho A Genetic analysis of Mint/X11 proteins: Essential presynaptic functions of a neuronal adaptor protein family J. Neurosci. 2006 26 50 13089 13101 10.1523/JNEUROSCI.2855-06.2006 17167098
Ho, A. et al. Genetic analysis of Mint/X11 proteins: Essential presynaptic functions of a neuronal adaptor protein family. J. Neurosci. 26(50), 13089–13101 (2006).17167098 10.1523/JNEUROSCI.2855-06.2006
58. Snidal CA Molecular characterization of an SV capture site in the mid-region of the presynaptic CaV2.1 calcium channel C-terminal Front. Cell. Neurosci. 2018 12 127 10.3389/fncel.2018.00127 29867360
Snidal, C. A. et al. Molecular characterization of an SV capture site in the mid-region of the presynaptic CaV2.1 calcium channel C-terminal. Front. Cell. Neurosci. 12, 127 (2018).29867360 10.3389/fncel.2018.00127
59. Smith CL Novel cell types, neurosecretory cells, and body plan of the early-diverging metazoan Trichoplax adhaerens Curr. Biol. 2014 24 14 1565 1572 10.1016/j.cub.2014.05.046 24954051
Smith, C. L. et al. Novel cell types, neurosecretory cells, and body plan of the early-diverging metazoan Trichoplax adhaerens. Curr. Biol. 24(14), 1565–1572 (2014).24954051 10.1016/j.cub.2014.05.046
60. Moroz LL Kohn AB Independent origins of neurons and synapses: Insights from ctenophores Philos. Trans. R. Soc. B: Biol. Sci. 2016 371 1685 20150041 10.1098/rstb.2015.0041
Moroz, L. L. & Kohn, A. B. Independent origins of neurons and synapses: Insights from ctenophores. Philos. Trans. R. Soc. B: Biol. Sci. 371(1685), 20150041 (2016).10.1098/rstb.2015.0041
61. Jones P InterProScan 5: Genome-scale protein function classification Bioinformatics 2014 30 9 1236 1240 10.1093/bioinformatics/btu031 24451626
Jones, P. et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 30(9), 1236–1240 (2014).24451626 10.1093/bioinformatics/btu031
62. Wu X Structural Basis for the High-Affinity Interaction between CASK and Mint1 Structure 2020 28 6 664 673.e3 10.1016/j.str.2020.04.001 32348748
Wu, X. et al. Structural Basis for the High-Affinity Interaction between CASK and Mint1. Structure 28(6), 664-673.e3 (2020).32348748 10.1016/j.str.2020.04.001
63. Borg J-P Identification of an evolutionarily conserved heterotrimeric protein complex involved in protein targeting J. Biol. Chem. 1998 273 48 31633 31636 10.1074/jbc.273.48.31633 9822620
Borg, J.-P. et al. Identification of an evolutionarily conserved heterotrimeric protein complex involved in protein targeting. J. Biol. Chem. 273(48), 31633–31636 (1998).9822620 10.1074/jbc.273.48.31633
64. Butz S Okamoto M Südhof TC A tripartite protein complex with the potential to couple synaptic vesicle exocytosis to cell adhesion in brain Cell 1998 94 6 773 782 10.1016/S0092-8674(00)81736-5 9753324
Butz, S., Okamoto, M. & Südhof, T. C. A tripartite protein complex with the potential to couple synaptic vesicle exocytosis to cell adhesion in brain. Cell 94(6), 773–782 (1998).9753324 10.1016/S0092-8674(00)81736-5
65. Mukherjee K Neuron-specific protein interactions of Drosophila CASK-β are revealed by mass spectrometry Front. Mol. Neurosci. 2014 7 58 10.3389/fnmol.2014.00058 25071438
Mukherjee, K. et al. Neuron-specific protein interactions of Drosophila CASK-β are revealed by mass spectrometry. Front. Mol. Neurosci. 7, 58 (2014).25071438 10.3389/fnmol.2014.00058
66. Zhang Z CASK modulates the assembly and function of the Mint1/Munc18-1 complex to regulate insulin secretion Cell Discov. 2020 6 1 1 17 10.1038/s41421-020-00216-3 31934347
Zhang, Z. et al. CASK modulates the assembly and function of the Mint1/Munc18-1 complex to regulate insulin secretion. Cell Discov. 6(1), 1–17 (2020).31934347 10.1038/s41421-020-00216-3
67. Tomita S Interaction of a neuron-specific protein containing PDZ domains with Alzheimer's amyloid precursor protein J. Biol. Chem. 1999 274 4 2243 2254 10.1074/jbc.274.4.2243 9890987
Tomita, S. et al. Interaction of a neuron-specific protein containing PDZ domains with Alzheimer’s amyloid precursor protein. J. Biol. Chem. 274(4), 2243–2254 (1999).9890987 10.1074/jbc.274.4.2243
68. Teber, I. et al. Rab6 interacts with the mint3 adaptor protein. Biol. Chem. 386(7), 671–677 (2005).
69. Zhang Z Sequence-specific recognition of the internalization motif of the Alzheimer's amyloid precursor protein by the X11 PTB domain EMBO J. 1997 16 20 6141 6150 10.1093/emboj/16.20.6141 9321393
Zhang, Z. et al. Sequence-specific recognition of the internalization motif of the Alzheimer’s amyloid precursor protein by the X11 PTB domain. EMBO J. 16(20), 6141–6150 (1997).9321393 10.1093/emboj/16.20.6141
70. Wong FK Synaptic vesicle tethering and the CaV2.2 distal C-terminal Front. Cell. Neurosci. 2014 8 71 10.3389/fncel.2014.00071 24639630
Wong, F. K. et al. Synaptic vesicle tethering and the CaV2.2 distal C-terminal. Front. Cell. Neurosci. 8, 71 (2014).24639630 10.3389/fncel.2014.00071
71. Selkoe DJ The cell biology of β-amyloid precursor protein and presenilin in Alzheimer's disease Trends in cell biology 1998 8 11 447 453 10.1016/S0962-8924(98)01363-4 9854312
Selkoe, D. J. The cell biology of β-amyloid precursor protein and presenilin in Alzheimer’s disease. Trends in cell biology 8(11), 447–453 (1998).9854312 10.1016/S0962-8924(98)01363-4
72. Ganot P Structural molecular components of septate junctions in cnidarians point to the origin of epithelial junctions in eukaryotes Mol. Boil. Evol. 2015 32 1 44 62 10.1093/molbev/msu265
Ganot, P. et al. Structural molecular components of septate junctions in cnidarians point to the origin of epithelial junctions in eukaryotes. Mol. Boil. Evol. 32(1), 44–62 (2015).10.1093/molbev/msu265
73. Guzman, C. et al. Neuronal and non-neuronal functions of the synaptic cell adhesion molecule neurexin in Nematostella vectensis. Nat. Commun. 15(1), 6495 (2024).
74. McLoughlin DM The neuronal adaptor protein X11α interacts with the copper chaperone for SOD1 and regulates SOD1 activity J. Biol. Chem. 2001 276 12 9303 9307 10.1074/jbc.M010023200 11115513
McLoughlin, D. M. et al. The neuronal adaptor protein X11α interacts with the copper chaperone for SOD1 and regulates SOD1 activity. J. Biol. Chem. 276(12), 9303–9307 (2001).11115513 10.1074/jbc.M010023200
75. Tomita S PDZ domain-dependent suppression of NF-κB/p65-induced Aβ42 production by a neuron-specific X11-like protein J. Biol. Chem. 2000 275 17 13056 13060 10.1074/jbc.C000019200 10777610
Tomita, S. et al. PDZ domain-dependent suppression of NF-κB/p65-induced Aβ42 production by a neuron-specific X11-like protein. J. Biol. Chem. 275(17), 13056–13060 (2000).10777610 10.1074/jbc.C000019200
76. Kamm K Trichoplax genomes reveal profound admixture and suggest stable wild populations without bisexual reproduction Sci. Rep. 2018 8 1 1 11 10.1038/s41598-018-29400-y 29311619
Kamm, K. et al. Trichoplax genomes reveal profound admixture and suggest stable wild populations without bisexual reproduction. Sci. Rep. 8(1), 1–11 (2018).29311619 10.1038/s41598-018-29400-y
77. Scheufler C Structure of TPR domain–peptide complexes: Critical elements in the assembly of the Hsp70–Hsp90 multichaperone machine Cell 2000 101 2 199 210 10.1016/S0092-8674(00)80830-2 10786835
Scheufler, C. et al. Structure of TPR domain–peptide complexes: Critical elements in the assembly of the Hsp70–Hsp90 multichaperone machine. Cell 101(2), 199–210 (2000).10786835 10.1016/S0092-8674(00)80830-2
78. Tao X A gating charge transfer center in voltage sensors Science 2010 328 5974 67 73 10.1126/science.1185954 20360102
Tao, X. et al. A gating charge transfer center in voltage sensors. Science 328(5974), 67–73 (2010).20360102 10.1126/science.1185954
79. Sakarya O Evolutionary expansion and specialization of the PDZ domains Mol. Boil. Evol. 2010 27 5 1058 1069 10.1093/molbev/msp311
Sakarya, O. et al. Evolutionary expansion and specialization of the PDZ domains. Mol. Boil. Evol. 27(5), 1058–1069 (2010).10.1093/molbev/msp311
80. Skelton NJ Origins of PDZ domain ligand specificity: Structure determination and mutagenesis of the Erbin PDZ domain J. Biol. Chem. 2003 278 9 7645 7654 10.1074/jbc.M209751200 12446668
Skelton, N. J. et al. Origins of PDZ domain ligand specificity: Structure determination and mutagenesis of the Erbin PDZ domain. J. Biol. Chem. 278(9), 7645–7654 (2003).12446668 10.1074/jbc.M209751200
81. Duquesne AE Solution structure of the second PDZ domain of the neuronal adaptor X11α and its interaction with the C-terminal peptide of the human copper chaperone for superoxide dismutase J. Biomol. NMR 2005 32 3 209 218 10.1007/s10858-005-7333-1 16132821
Duquesne, A. E. et al. Solution structure of the second PDZ domain of the neuronal adaptor X11α and its interaction with the C-terminal peptide of the human copper chaperone for superoxide dismutase. J. Biomol. NMR 32(3), 209–218 (2005).16132821 10.1007/s10858-005-7333-1
82. Jumper J Highly accurate protein structure prediction with AlphaFold Nature 2021 596 7873 583 589 10.1038/s41586-021-03819-2 34265844
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596(7873), 583–589 (2021).34265844 10.1038/s41586-021-03819-2
83. Ashkenazy H ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules Nucleic Acids Res. 2016 44 W1 W344 W350 10.1093/nar/gkw408 27166375
Ashkenazy, H. et al. ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 44(W1), W344–W350 (2016).27166375 10.1093/nar/gkw408
84. Wright PE Dyson HJ Intrinsically disordered proteins in cellular signalling and regulation Nat. Rev. Mol. Cell Biol. 2015 16 1 18 29 10.1038/nrm3920 25531225
Wright, P. E. & Dyson, H. J. Intrinsically disordered proteins in cellular signalling and regulation. Nat. Rev. Mol. Cell Biol. 16(1), 18–29 (2015).25531225 10.1038/nrm3920
85. Warner JF NvERTx: A gene expression database to compare embryogenesis and regeneration in the sea anemone Nematostella vectensis Development 2018 145 10 dev162887 10.1242/dev.162867
Warner, J. F. et al. NvERTx: A gene expression database to compare embryogenesis and regeneration in the sea anemone Nematostella vectensis. Development 145(10), dev162887 (2018).10.1242/dev.162867
86. Sebé-Pedrós A Cnidarian cell type diversity and regulation revealed by whole-organism single-cell RNA-Seq Cell 2018 173 6 1520 1534.e20 10.1016/j.cell.2018.05.019 29856957
Sebé-Pedrós, A. et al. Cnidarian cell type diversity and regulation revealed by whole-organism single-cell RNA-Seq. Cell 173(6), 1520-1534.e20 (2018).29856957 10.1016/j.cell.2018.05.019
87. Weir K A molecular filter for the cnidarian stinging response Elife 2020 9 e57578 10.7554/eLife.57578 32452384
Weir, K. et al. A molecular filter for the cnidarian stinging response. Elife 9, e57578 (2020).32452384 10.7554/eLife.57578
88. Havrilak JA Characterization of NvLWamide-like neurons reveals stereotypy in Nematostella nerve net development Dev. Boil. 2017 431 2 336 346 10.1016/j.ydbio.2017.08.028
Havrilak, J. A. et al. Characterization of NvLWamide-like neurons reveals stereotypy in Nematostella nerve net development. Dev. Boil. 431(2), 336–346 (2017).10.1016/j.ydbio.2017.08.028
89. Admoni Y TATA binding protein (TBP) promoter drives ubiquitous expression of marker transgene in the adult sea anemone Nematostella vectensis Genes 2020 11 9 1081 10.3390/genes11091081 32948049
Admoni, Y. et al. TATA binding protein (TBP) promoter drives ubiquitous expression of marker transgene in the adult sea anemone Nematostella vectensis. Genes 11(9), 1081 (2020).32948049 10.3390/genes11091081
90. Cole AG Updated single cell reference atlas for the starlet anemone Nematostella vectensis Front. Zool. 2024 21 1 8 10.1186/s12983-024-00529-z 38500146
Cole, A. G. et al. Updated single cell reference atlas for the starlet anemone Nematostella vectensis. Front. Zool. 21(1), 8 (2024).38500146 10.1186/s12983-024-00529-z
91. Dolphin AC Lee A Presynaptic calcium channels: Specialized control of synaptic neurotransmitter release Nat. Rev. Neurosci. 2020 21 4 213 229 10.1038/s41583-020-0278-2 32161339
Dolphin, A. C. & Lee, A. Presynaptic calcium channels: Specialized control of synaptic neurotransmitter release. Nat. Rev. Neurosci. 21(4), 213–229 (2020).32161339 10.1038/s41583-020-0278-2
92. Senatore A Zhorov BS Spafford JD Cav3 T-type calcium channels Wiley Interdisciplinary Reviews: Membrane Transport and Signaling 2012 1 4 467 491
Senatore, A., Zhorov, B. S. & Spafford, J. D. Cav3 T-type calcium channels. Wiley Interdisciplinary Reviews: Membrane Transport and Signaling 1(4), 467–491 (2012).
93. Smith CL Evolutionary insights into T-type Ca2+ channel structure, function, and ion selectivity from the Trichoplax adhaerens homologue J. Gen. Physiol. 2017 149 4 483 510 10.1085/jgp.201611683 28330839
Smith, C. L. et al. Evolutionary insights into T-type Ca2+ channel structure, function, and ion selectivity from the Trichoplax adhaerens homologue. J. Gen. Physiol. 149(4), 483–510 (2017).28330839 10.1085/jgp.201611683
94. Ichikawa DM A multireporter bacterial 2-hybrid assay for the high-throughput and dynamic assay of PDZ domain-peptide interactions ACS Synth. Biol. 2019 8 5 918 928 10.1021/acssynbio.8b00499 30969105
Ichikawa, D. M. et al. A multireporter bacterial 2-hybrid assay for the high-throughput and dynamic assay of PDZ domain-peptide interactions. ACS Synth. Biol. 8(5), 918–928 (2019).30969105 10.1021/acssynbio.8b00499
95. Teyra J Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities Protein Sci. 2020 29 2 433 442 10.1002/pro.3759 31654425
Teyra, J. et al. Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities. Protein Sci. 29(2), 433–442 (2020).31654425 10.1002/pro.3759
96. Moroz LL Romanova DY Alternative neural systems: What is a neuron?(Ctenophores, sponges and placozoans) Front. Cell Dev. Biol. 2022 10 1071961 10.3389/fcell.2022.1071961 36619868
Moroz, L. L. & Romanova, D. Y. Alternative neural systems: What is a neuron?(Ctenophores, sponges and placozoans). Front. Cell Dev. Biol. 10, 1071961 (2022).36619868 10.3389/fcell.2022.1071961
97. Laursen L Divergent evolution of a protein–protein interaction revealed through ancestral sequence reconstruction and resurrection Mol. Biol. Evol. 2021 38 1 152 167 10.1093/molbev/msaa198 32750125
Laursen, L. et al. Divergent evolution of a protein–protein interaction revealed through ancestral sequence reconstruction and resurrection. Mol. Biol. Evol. 38(1), 152–167 (2021).32750125 10.1093/molbev/msaa198
98. Wheeler LC Conservation of specificity in two low-specificity proteins Biochemistry 2018 57 5 684 695 10.1021/acs.biochem.7b01086 29240404
Wheeler, L. C. et al. Conservation of specificity in two low-specificity proteins. Biochemistry 57(5), 684–695 (2018).29240404 10.1021/acs.biochem.7b01086
99. Anderson DP Correction: Evolution of an ancient protein function involved in organized multicellularity in animals Elife 2016 5 e14311 10.7554/eLife.14311 26790885
Anderson, D. P. et al. Correction: Evolution of an ancient protein function involved in organized multicellularity in animals. Elife 5, e14311 (2016).26790885 10.7554/eLife.14311
100. Gao M Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis Protein Sci. 2020 29 11 2226 2244 10.1002/pro.3947 32914530
Gao, M. et al. Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis. Protein Sci. 29(11), 2226–2244 (2020).32914530 10.1002/pro.3947
101. Putnam NH Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization Science 2007 317 5834 86 94 10.1126/science.1139158 17615350
Putnam, N. H. et al. Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science 317(5834), 86–94 (2007).17615350 10.1126/science.1139158
102. Koo EH Squazzo SL Evidence that production and release of amyloid beta-protein involves the endocytic pathway J. Biol. Chem. 1994 269 26 17386 17389 10.1016/S0021-9258(17)32449-3 8021238
Koo, E. H. & Squazzo, S. L. Evidence that production and release of amyloid beta-protein involves the endocytic pathway. J. Biol. Chem. 269(26), 17386–17389 (1994).8021238 10.1016/S0021-9258(17)32449-3
103. Ring S The secreted β-amyloid precursor protein ectodomain APPsα is sufficient to rescue the anatomical, behavioral, and electrophysiological abnormalities of APP-deficient mice J. Neurosci. 2007 27 29 7817 7826 10.1523/JNEUROSCI.1026-07.2007 17634375
Ring, S. et al. The secreted β-amyloid precursor protein ectodomain APPsα is sufficient to rescue the anatomical, behavioral, and electrophysiological abnormalities of APP-deficient mice. J. Neurosci. 27(29), 7817–7826 (2007).17634375 10.1523/JNEUROSCI.1026-07.2007
104. Hampel H The amyloid-β pathway in Alzheimer’s disease Mol. Psychiatry 2021 26 10 5481 5503 10.1038/s41380-021-01249-0 34456336
Hampel, H. et al. The amyloid-β pathway in Alzheimer’s disease. Mol. Psychiatry 26(10), 5481–5503 (2021).34456336 10.1038/s41380-021-01249-0
105. Sullivan SE Mint proteins are required for synaptic activity-dependent amyloid precursor protein (APP) trafficking and amyloid β generation J. Biol. Chem. 2014 289 22 15374 15383 10.1074/jbc.M113.541003 24742670
Sullivan, S. E. et al. Mint proteins are required for synaptic activity-dependent amyloid precursor protein (APP) trafficking and amyloid β generation. J. Biol. Chem. 289(22), 15374–15383 (2014).24742670 10.1074/jbc.M113.541003
106. Lin AY A rare autism-associated MINT2/APBA2 mutation disrupts neurexin trafficking and synaptic function Sci. Rep. 2019 9 1 6024 10.1038/s41598-019-42635-7 30988517
Lin, A. Y. et al. A rare autism-associated MINT2/APBA2 mutation disrupts neurexin trafficking and synaptic function. Sci. Rep. 9(1), 6024 (2019).30988517 10.1038/s41598-019-42635-7
107. Servián-Morilla E Proteolytic processing of neurexins by presenilins sustains synaptic vesicle release J. Neurosci. 2018 38 4 901 917 10.1523/JNEUROSCI.1357-17.2017 29229705
Servián-Morilla, E. et al. Proteolytic processing of neurexins by presenilins sustains synaptic vesicle release. J. Neurosci. 38(4), 901–917 (2018).29229705 10.1523/JNEUROSCI.1357-17.2017
108. Senatore A Raiss H Le P Physiology and evolution of voltage-gated calcium channels in early diverging animal phyla: Cnidaria, Placozoa, Porifera and Ctenophora Front. Physiol. 2016 7 481 10.3389/fphys.2016.00481 27867359
Senatore, A., Raiss, H. & Le, P. Physiology and evolution of voltage-gated calcium channels in early diverging animal phyla: Cnidaria, Placozoa, Porifera and Ctenophora. Front. Physiol. 7, 481 (2016).27867359 10.3389/fphys.2016.00481
109. Bullock TH Neuromuscular facilitation in scyphomedusae J. Cell. Comp. Physiol. 1943 22 3 251 272 10.1002/jcp.1030220306
Bullock, T. H. Neuromuscular facilitation in scyphomedusae. J. Cell. Comp. Physiol. 22(3), 251–272 (1943).10.1002/jcp.1030220306
110. Kerfoot P Neuromuscular transmission in the jellyfish Aglantha digitale Journal of experimental biology 1985 116 1 1 25 10.1242/jeb.116.1.1 2865323
Kerfoot, P. et al. Neuromuscular transmission in the jellyfish Aglantha digitale. Journal of experimental biology 116(1), 1–25 (1985).2865323 10.1242/jeb.116.1.1
111. Ryan JF Chiodin M Where is my mind? How sponges and placozoans may have lost neural cell types Philos. Trans. R. Soc. B: Biol. Sci. 2015 370 1684 20150059 10.1098/rstb.2015.0059
Ryan, J. F. & Chiodin, M. Where is my mind? How sponges and placozoans may have lost neural cell types. Philos. Trans. R. Soc. B: Biol. Sci. 370(1684), 20150059 (2015).10.1098/rstb.2015.0059
112. Ho A A role for Mints in transmitter release: Mint 1 knockout mice exhibit impaired GABAergic synaptic transmission Proc. Natl. Acad. Sci. 2003 100 3 1409 1414 10.1073/pnas.252774899 12547917
Ho, A. et al. A role for Mints in transmitter release: Mint 1 knockout mice exhibit impaired GABAergic synaptic transmission. Proc. Natl. Acad. Sci. 100(3), 1409–1414 (2003).12547917 10.1073/pnas.252774899
113. Gross GG X11/Mint genes control polarized localization of axonal membrane proteins in vivo J. Neurosci. 2013 33 19 8575 8586 10.1523/JNEUROSCI.5749-12.2013 23658195
Gross, G. G. et al. X11/Mint genes control polarized localization of axonal membrane proteins in vivo. J. Neurosci. 33(19), 8575–8586 (2013).23658195 10.1523/JNEUROSCI.5749-12.2013
114. Rongo C LIN-10 is a shared component of the polarized protein localization pathways in neurons and epithelia Cell 1998 94 6 751 759 10.1016/S0092-8674(00)81734-1 9753322
Rongo, C. et al. LIN-10 is a shared component of the polarized protein localization pathways in neurons and epithelia. Cell 94(6), 751–759 (1998).9753322 10.1016/S0092-8674(00)81734-1
115. Najle SR Stepwise emergence of the neuronal gene expression program in early animal evolution Cell 2023 186 21 4676 4693.e29 10.1016/j.cell.2023.08.027 37729907
Najle, S. R. et al. Stepwise emergence of the neuronal gene expression program in early animal evolution. Cell 186(21), 4676-4693.e29 (2023).37729907 10.1016/j.cell.2023.08.027
116. Thul PJ Lindskog C The human protein atlas: A spatial map of the human proteome Protein Sci. 2018 27 1 233 244 10.1002/pro.3307 28940711
Thul, P. J. & Lindskog, C. The human protein atlas: A spatial map of the human proteome. Protein Sci. 27(1), 233–244 (2018).28940711 10.1002/pro.3307
117. Armon S Ultrafast epithelial contractions provide insights into contraction speed limits and tissue integrity Proc. Natl. Acad. Sci. 2018 115 44 E10333 E10341 10.1073/pnas.1802934115 30309963
Armon, S. et al. Ultrafast epithelial contractions provide insights into contraction speed limits and tissue integrity. Proc. Natl. Acad. Sci. 115(44), E10333–E10341 (2018).30309963 10.1073/pnas.1802934115
118. Kim J Rewiring of PDZ domain-ligand interaction network contributed to eukaryotic evolution PLoS Genetics 2012 8 2 e1002510 10.1371/journal.pgen.1002510 22346764
Kim, J. et al. Rewiring of PDZ domain-ligand interaction network contributed to eukaryotic evolution. PLoS Genetics 8(2), e1002510 (2012).22346764 10.1371/journal.pgen.1002510
119. Songyang Z Recognition of unique carboxyl-terminal motifs by distinct PDZ domains Science 1997 275 5296 73 77 10.1126/science.275.5296.73 8974395
Songyang, Z. et al. Recognition of unique carboxyl-terminal motifs by distinct PDZ domains. Science 275(5296), 73–77 (1997).8974395 10.1126/science.275.5296.73
120. Wong, Y. Y. et al. Transcriptome profiling of Trichoplax adhaerens highlights its digestive epithelium and a rich set of genes for fast electrogenic and slow neuromodulatory cellular signaling. Res. Sq. preprint (2019).
121. Eitel M Comparative genomics and the nature of placozoan species PLoS Biol. 2018 16 7 e2005359 10.1371/journal.pbio.2005359 30063702
Eitel, M. et al. Comparative genomics and the nature of placozoan species. PLoS Biol. 16(7), e2005359 (2018).30063702 10.1371/journal.pbio.2005359
122. Coordinators NR Database resources of the national center for biotechnology information Nucleic Acids Res. 2016 44 D1 D7 D19 10.1093/nar/gkv1290 26615191
Coordinators, N. R. Database resources of the national center for biotechnology information. Nucleic Acids Res. 44(D1), D7–D19 (2016).26615191 10.1093/nar/gkv1290
123. Edgar RC MUSCLE: Multiple sequence alignment with high accuracy and high throughput Nucleic acids research 2004 32 5 1792 1797 10.1093/nar/gkh340 15034147
Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic acids research 32(5), 1792–1797 (2004).15034147 10.1093/nar/gkh340
124. Capella-Gutiérrez S Silla-Martínez JM Gabaldón T trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses Bioinformatics 2009 25 15 1972 1973 10.1093/bioinformatics/btp348 19505945
Capella-Gutiérrez, S., Silla-Martínez, J. M. & Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15), 1972–1973 (2009).19505945 10.1093/bioinformatics/btp348
125. Minh BQ IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era Mol. Biol. Evol. 2020 37 5 1530 1534 10.1093/molbev/msaa015 32011700
Minh, B. Q. et al. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37(5), 1530–1534 (2020).32011700 10.1093/molbev/msaa015
126. Waterhouse AM Jalview Version 2—a multiple sequence alignment editor and analysis workbench Bioinformatics 2009 25 9 1189 1191 10.1093/bioinformatics/btp033 19151095
Waterhouse, A. M. et al. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25(9), 1189–1191 (2009).19151095 10.1093/bioinformatics/btp033
127. Rice P Longden I Bleasby A EMBOSS: The European molecular biology open software suite Trends Genetics 2000 16 6 276 277 10.1016/S0168-9525(00)02024-2
Rice, P., Longden, I. & Bleasby, A. EMBOSS: The European molecular biology open software suite. Trends Genetics 16(6), 276–277 (2000).10.1016/S0168-9525(00)02024-2
128. Baldwin, S. Compute Canada: Advancing computational research. J. Phys. Conf. Ser. 341(1), 012001 (2012).
129. Goddard TD UCSF ChimeraX: Meeting modern challenges in visualization and analysis Protein Sci. 2018 27 1 14 25 10.1002/pro.3235 28710774
Goddard, T. D. et al. UCSF ChimeraX: Meeting modern challenges in visualization and analysis. Protein Sci. 27(1), 14–25 (2018).28710774 10.1002/pro.3235
130. Kelley LA The Phyre2 web portal for protein modeling, prediction and analysis Nat. Protoc. 2015 10 6 845 858 10.1038/nprot.2015.053 25950237
Kelley, L. A. et al. The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10(6), 845–858 (2015).25950237 10.1038/nprot.2015.053
131. Emsley P Cowtan K Coot: Model-building tools for molecular graphics Acta Crystallogr. Sect. D: Biol. Crystallogr. 2004 60 12 2126 2132 10.1107/S0907444904019158 15572765
Emsley, P. & Cowtan, K. Coot: Model-building tools for molecular graphics. Acta Crystallogr. Sect. D: Biol. Crystallogr. 60(12), 2126–2132 (2004).15572765 10.1107/S0907444904019158
132. Ramachandran S Automated minimization of steric clashes in protein structures Proteins Struct. Funct. Bioinform. 2011 79 1 261 270 10.1002/prot.22879
Ramachandran, S. et al. Automated minimization of steric clashes in protein structures. Proteins Struct. Funct. Bioinform. 79(1), 261–270 (2011).10.1002/prot.22879
133. Micallef L Rodgers P euler APE: Drawing area-proportional 3-Venn diagrams using ellipses PLoS ONE 2014 9 7 e101717 10.1371/journal.pone.0101717 25032825
Micallef, L. & Rodgers, P. euler APE: Drawing area-proportional 3-Venn diagrams using ellipses. PLoS ONE 9(7), e101717 (2014).25032825 10.1371/journal.pone.0101717
134. Schindelin J Fiji: An open-source platform for biological-image analysis Nat. Methods 2012 9 7 676 682 10.1038/nmeth.2019 22743772
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9(7), 676–682 (2012).22743772 10.1038/nmeth.2019
135. Senatore A Boone AN Spafford JD Optimized transfection strategy for expression and electrophysiological recording of recombinant voltage-gated ion channels in HEK-293T cells JoVE 2011 47 e2314
Senatore, A., Boone, A. N. & Spafford, J. D. Optimized transfection strategy for expression and electrophysiological recording of recombinant voltage-gated ion channels in HEK-293T cells. JoVE 47, e2314 (2011).
136. Persikov AV Deep sequencing of large library selections allows computational discovery of diverse sets of zinc fingers that bind common targets Nucleic Acids Res. 2014 42 3 1497 1508 10.1093/nar/gkt1034 24214968
Persikov, A. V. et al. Deep sequencing of large library selections allows computational discovery of diverse sets of zinc fingers that bind common targets. Nucleic Acids Res. 42(3), 1497–1508 (2014).24214968 10.1093/nar/gkt1034
