Somatosensory neurons have historically been classified by a variety of approaches, including structural, anatomical, and genetic markers; electrophysiological properties; pharmacological sensitivities; and more recently, transcriptional profile differentiation. These methodologies, used separately, have yielded inconsistent classification schemes. Here, we describe phenotypic differences in response to pharmacological agents as measured by changes in cytosolic calcium concentration for the rapid classification of neurons in vitro; further analysis with genetic markers, whole-cell recordings, and single-cell transcriptomics validated these findings in a functional context. Using this general approach, which we refer to as tripartite constellation analysis (TCA), we focused on large-diameter dorsal-root ganglion (L-DRG) neurons with myelinated axons. Divergent responses to the K-channel antagonist, κM-conopeptide RIIIJ (RIIIJ), reliably identified six discrete functional cell classes. In two neuronal subclasses (L1 and L2), block with RIIIJ led to an increase in [Ca]i. Simultaneous electrophysiology and calcium imaging showed that the RIIIJ-elicited increase in [Ca]i corresponded to different patterns of action potentials (APs), a train of APs in L1 neurons, and sporadic firing in L2 neurons. Genetically labeled mice established that L1 neurons are proprioceptors. The single-cell transcriptomes of L1 and L2 neurons showed that L2 neurons are Aδ–low-threshold mechanoreceptors. RIIIJ effects were replicated by application of the Kv1.1 selective antagonist, Dendrotoxin-K, in several L-DRG subclasses (L1, L2, L3, and L5), suggesting the presence of functional Kv1.1/Kv1.2 heteromeric channels. Using this approach on other neuronal subclasses should ultimately accelerate the comprehensive classification and characterization of individual somatosensory neuronal subclasses within a mixed population.

The molecular sensor of innocuous (painless) cold sensation is well-established to be Transient Receptor Potential cation channel, subfamily M, member 8 (TRPM8). However, the role of Transient Receptor Potential cation channel, subfamily A, member 1 (TRPA1) in noxious (painful) cold sensation has been controversial. We find that TRPA1 channels contribute to the noxious cold sensitivity of mouse somatosensory neurons, independent of TRPM8 channels, and that TRPA1-expressing neurons are largely non-overlapping with TRPM8-expressing neurons in mouse dorsal-root ganglia (DRG). However, relatively few TRPA1-expressing neurons (e.g., responsive to allyl isothiocyanate or AITC, a selective TRPA1 agonist) respond overtly to cold temperature in vitro, unlike TRPM8-expressing neurons, which almost all respond to cold. Using somatosensory neurons from TRPM8 −/− mice and subtype-selective blockers of TRPM8 and TRPA1 channels, we demonstrate that responses to cold temperatures from TRPA1-expressing neurons are mediated by TRPA1 channels. We also identify two factors that affect the cold-sensitivity of TRPA1-expressing neurons: 1) Cold-sensitive AITC-sensitive neurons express relatively more TRPA1 transcripts than cold-insensitive AITC-sensitive neurons and 2) voltage-gated potassium (KV) channels attenuate the cold-sensitivity of some TRPA1-expressing neurons. The combination of these two factors, combined with the relatively weak agonist-like activity of cold temperature on TRPA1 channels, partially explains why few TRPA1-expressing neurons respond to cold. Blocking KV channels also reveals another subclass of noxious cold-sensitive DRG neurons that do not express TRPM8 or TRPA1 channels. Altogether, the results of this study provide novel insights into the cold-sensitivity of different subclasses of somatosensory neurons.

The vast complexity of native heteromeric K+ channels is largely unexplored. Defining the composition and subunit arrangement of individual subunits in native heteromeric K+ channels and establishing their physiological roles is experimentally challenging. Here we systematically explored this “zone of ignorance” in molecular neuroscience. Venom components, such as peptide toxins, appear to have evolved to modulate physiologically relevant targets by discriminating among closely related native ion channel complexes. We provide proof-of-principle for this assertion by demonstrating that κM-conotoxin RIIIJ (κM-RIIIJ) from Conus radiatus precisely targets “asymmetric” Kv channels composed of three Kv1.2 subunits and one Kv1.1 or Kv1.6 subunit with 100-fold higher apparent affinity compared with homomeric Kv1.2 channels. Our study shows that dorsal root ganglion (DRG) neurons contain at least two major functional Kv1.2 channel complexes: a heteromer, for which κM-RIIIJ has high affinity, and a putative Kv1.2 homomer, toward which κM-RIIIJ is less potent. This conclusion was reached by (i) covalent linkage of members of the mammalian Shaker-related Kv1 family to Kv1.2 and systematic assessment of the potency of κM-RIIIJ block of heteromeric K+ channel-mediated currents in heterologous expression systems; (ii) molecular dynamics simulations of asymmetric Kv1 channels providing insights into the molecular basis of κM-RIIIJ selectivity and potency toward its targets; and (iii) evaluation of calcium responses of a defined population of DRG neurons to κM-RIIIJ. Our study demonstrates that bioactive molecules present in venoms provide essential pharmacological tools that systematically target specific heteromeric Kv channel complexes that operate in native tissues.

A pressing need in neurobiology is the comprehensive identification and characterization of neuronal subclasses within the mammalian nervous system. To this end, we used constellation pharmacology as a method to interrogate the neuronal and glial subclasses of the mouse cerebellum individually and simultaneously. We then evaluated the data obtained from constellation-pharmacology experiments by cluster analysis to classify cells into neuronal and glial subclasses, based on their functional expression of glutamate, acetylcholine, and GABA receptors, among other ion channels. Conantokin peptides were used to identify N-methyl-d-aspartate (NMDA) receptor subtypes, which revealed that neurons of the young mouse cerebellum expressed NR2A and NR2B NMDA receptor subunits. Additional pharmacological tools disclosed differential expression of α-amino-3-hydroxy-5-methyl-4-isoxazloepropionic, nicotinic acetylcholine, and muscarinic acetylcholine receptors in different neuronal and glial subclasses. Certain cell subclasses correlated with known attributes of granule cells, and we combined constellation pharmacology with genetically labeled neurons to identify and characterize Purkinje cells. This study illustrates the utility of applying constellation pharmacology to classify neuronal and glial subclasses in specific anatomical regions of the brain.

In a program to identify pain treatments with low addiction potential, we isolated five steroids, conosteroids A–E (1–5), from the hypobranchial gland of the mollusk Conus geographus. Compounds 1–5 were active in a mouse dorsal root ganglion (DRG) assay that suggested that they might be analgesic. A synthetic analogue 6 was used for a detailed pharmacological study. Compound 6 significantly increased the pain threshold in mice in the hot-plate test at 2 and 50 mg/kg. Compound 6 at 500 nM antagonizes type-A γ-aminobutyric acid receptors (GABAARs). In a patch-clamp experiment, out of the six subunit combinations tested, 6 exhibited subtype selectivity, most strongly antagonizing α1β1γ2 and α4β3γ2 receptors (IC50 1.5 and 1.0 μM, respectively). Although the structures of 1–6 differ from those of known neuroactive steroids, they are cell-type-selective modulators of GABAARs, expanding the known chemical space of neuroactive steroids.