Analysing an allelic series of rare missense variants of CACNA1I in a Swedish schizophrenia cohort

Using automated patch clamp electrophysiology platforms in ion channel drug discovery: an industry perspective

INTRODUCTION

Automated patch clamp (APC) is now well established as a mature technology for ion channel drug discovery in academia, biotech and pharma companies, and in contract research organizations (CRO), for a variety of applications including channelopathy research, compound screening, target validation and cardiac safety testing.

AREAS COVERED

Ion channels are an important class of drugged and approved drug targets. The authors present a review of the current state of ion channel drug discovery along with new and exciting developments in ion channel research involving APC. This includes topics such as native and iPSC-derived cells in ion channel drug discovery, channelopathy research, organellar and biologics in ion channel drug discovery.

EXPERT OPINION

It is our belief that APC will continue to play a critical role in ion channel drug discovery, not only in 'classical' hit screening, target validation and cardiac safety testing, but extending these applications to include high throughput organellar recordings and optogenetics. In this way, with advancements in APC capabilities and applications, together with high resolution cryo-EM structures, ion channel drug discovery will be re-invigorated, leading to a growing list of ion channel ligands in clinical development.

Acute ketamine induces neuronal hyperexcitability and deficits in prepulse inhibition by upregulating IL-6

Acute ketamine administration results in psychotic symptoms similar to those observed in schizophrenia and is regarded as a pharmacological model of schizophrenia. Accumulating evidence suggests that patients with schizophrenia show increased IL-6 levels in the blood and cerebrospinal fluid and that IL-6 levels are associated with the severity of psychotic symptoms. In the present study, we found that a single ketamine exposure led to increased expression of IL-6 and IL-6Rα, decreased dendritic spine density, increased expression and currents of T-type calcium channels, and increased neuron excitability in the hippocampal CA1 area 12 h after exposure. Acute ketamine administration also led to impaired prepulse inhibition (PPI) 12 h after administration. Additionally, we found that the expression of signaling molecules IKKα/β, NF-κB, JAK2, and STAT3 was upregulated 12 h after a single ketamine injection. The decreases in dendritic spine density, the increases in calcium currents and neuron excitability, and the impairments in PPI were ameliorated by blocking IL-6 or IL-6Rα. Our findings show that blocking IL-6 or its receptor may protect hippocampal neurons from hyperexcitability, thereby ameliorating ketamine-induced psychotic effects. Our study provides additional evidence that targeting IL-6 and its receptor is a potential strategy for treating psychotic symptoms in acute ketamine-induced psychosis and schizophrenia.

The T-type calcium channelosome

T-type calcium channels perform crucial physiological roles across a wide spectrum of tissues, spanning both neuronal and non-neuronal system. For instance, they serve as pivotal regulators of neuronal excitability, contribute to cardiac pacemaking, and mediate the secretion of hormones. These functions significantly hinge upon the intricate interplay of T-type channels with interacting proteins that modulate their expression and function at the plasma membrane. In this review, we offer a panoramic exploration of the current knowledge surrounding these T-type channel interactors, and spotlight certain aspects of their potential for drug-based therapeutic intervention.

The Thalamocortical Mechanism Underlying the Generation and Regulation of the Auditory Steady-State Responses in Awake Mice

The auditory steady-state response (ASSR) is a cortical oscillation induced by trains of 40 Hz acoustic stimuli. While the ASSR has been widely used in clinic measurement, the underlying neural mechanism remains poorly understood. In this study, we investigated the contribution of different stages of auditory thalamocortical pathway-medial geniculate body (MGB), thalamic reticular nucleus (TRN), and auditory cortex (AC)-to the generation and regulation of 40 Hz ASSR in C57BL/6 mice of both sexes. We found that the neural response synchronizing to 40 Hz sound stimuli was most prominent in the GABAergic neurons in the granular layer of AC and the ventral division of MGB (MGBv), which were regulated by optogenetic manipulation of TRN neurons. Behavioral experiments confirmed that disrupting TRN activity has a detrimental effect on the ability of mice to discriminate 40 Hz sounds. These findings revealed a thalamocortical mechanism helpful to interpret the results of clinical ASSR examinations.Significance Statement Our study contributes to clarifying the thalamocortical mechanisms underlying the generation and regulation of the auditory steady-state response (ASSR), which is commonly used in both clinical and neuroscience research to assess the integrity of auditory function. Combining a series of electrophysiological and optogenetic experiments, we demonstrate that the generation of cortical ASSR is dependent on the lemniscal thalamocortical projections originating from the ventral division of medial geniculate body to the GABAergic interneurons in the granule layer of the auditory cortex. Furthermore, the thalamocortical process for ASSR is strictly regulated by the activity of thalamic reticular nucleus (TRN) neurons. Behavioral experiments confirmed that dysfunction of TRN would cause a disruption of mice's behavioral performance in the auditory discrimination task.

Molecular and Cellular Context Influences SCN8A Variant Function

Pathogenic variants in SCN8A , which encodes the voltage-gated sodium (Na V ) channel Na V 1.6, are associated with neurodevelopmental disorders including epileptic encephalopathy. Previous approaches to determine SCN8A variant function may be confounded by the use of a neonatal-expressed alternatively spliced isoform of Na V 1.6 (Na V 1.6N), and engineered mutations to render the channel tetrodotoxin (TTX) resistant. In this study, we investigated the impact of SCN8A alternative splicing on variant function by comparing the functional attributes of 15 variants expressed in two developmentally regulated splice isoforms (Na V 1.6N, Na V 1.6A). We employed automated patch clamp recording to enhance throughput, and developed a novel neuronal cell line (ND7/LoNav) with low levels of endogenous Na V current to obviate the need for TTX-resistance mutations. Expression of Na V 1.6N or Na V 1.6A in ND7/LoNav cells generated Na V currents that differed significantly in voltage-dependence of activation and inactivation. TTX-resistant versions of both isoforms exhibited significant functional differences compared to the corresponding wild-type (WT) channels. We demonstrated that many of the 15 disease-associated variants studied exhibited isoform-dependent functional effects, and that many of the studied SCN8A variants exhibited functional properties that were not easily classified as either gain- or loss-of-function. Our work illustrates the value of considering molecular and cellular context when investigating SCN8A variants.

Prophecy or empiricism? Clinical value of predicting versus determining genetic variant functions

The recent explosion of epilepsy genetic testing has created challenges for interpretation of gene variants. Assessments of the functional consequences of genetic variants either by predictive or experimental strategies can contribute to estimating pathogenicity, but there is no consensus on which approach is best. The Special Interest Group on Epilepsy Genetics hosted a session during the Annual American Epilepsy Society Meeting in December 2022 to discuss this topic. The session featured a debate of the relative advantages and limitations of predicting (prophecy) versus experimentally determining (empiricism) variant function using ion channel gene variants as examples. This commentary summarizes these discussions.

Classification of missense variants in the N-methyl-d-aspartate receptor GRIN gene family as gain- or loss-of-function

Advances in sequencing technology have generated a large amount of genetic data from patients with neurological conditions. These data have provided diagnosis of many rare diseases, including a number of pathogenic de novo missense variants in GRIN genes encoding N-methyl-d-aspartate receptors (NMDARs). To understand the ramifications for neurons and brain circuits affected by rare patient variants, functional analysis of the variant receptor is necessary in model systems. For NMDARs, this functional analysis needs to assess multiple properties in order to understand how variants could impact receptor function in neurons. One can then use these data to determine whether the overall actions will increase or decrease NMDAR-mediated charge transfer. Here, we describe an analytical and comprehensive framework by which to categorize GRIN variants as either gain-of-function (GoF) or loss-of-function (LoF) and apply this approach to GRIN2B variants identified in patients and the general population. This framework draws on results from six different assays that assess the impact of the variant on NMDAR sensitivity to agonists and endogenous modulators, trafficking to the plasma membrane, response time course and channel open probability. We propose to integrate data from multiple in vitro assays to arrive at a variant classification, and suggest threshold levels that guide confidence. The data supporting GoF and LoF determination are essential to assessing pathogenicity and patient stratification for clinical trials as personalized pharmacological and genetic agents that can enhance or reduce receptor function are advanced. This approach to functional variant classification can generalize to other disorders associated with missense variants.

The Action Potential Clamp Technique as a Tool for Risk Stratification of Sinus Bradycardia Due to Loss-of-Function Mutations in HCN4: An In Silico Exploration Based on In Vitro and In Vivo Data

These days, in vitro functional analysis of gene variants is becoming increasingly important for risk stratification of cardiac ion channelopathies. So far, such risk stratification has been applied to SCN5A, KCNQ1, and KCNH2 gene variants associated with Brugada syndrome and long QT syndrome types 1 and 2, respectively, but risk stratification of HCN4 gene variants related to sick sinus syndrome has not yet been performed. HCN4 is the gene responsible for the hyperpolarization-activated 'funny' current If, which is an important modulator of the spontaneous diastolic depolarization underlying the sinus node pacemaker activity. In the present study, we carried out a risk classification assay on those loss-of-function mutations in HCN4 for which in vivo as well as in vitro data have been published. We used the in vitro data to compute the charge carried by If (Qf) during the diastolic depolarization phase of a prerecorded human sinus node action potential waveform and assessed the extent to which this Qf predicts (1) the beating rate of the comprehensive Fabbri-Severi model of a human sinus node cell with mutation-induced changes in If and (2) the heart rate observed in patients carrying the associated mutation in HCN4. The beating rate of the model cell showed a very strong correlation with Qf from the simulated action potential clamp experiments (R2 = 0.95 under vagal tone). The clinically observed minimum or resting heart rates showed a strong correlation with Qf (R2 = 0.73 and R2 = 0.71, respectively). While a translational perspective remains to be seen, we conclude that action potential clamp on transfected cells, without the need for further voltage clamp experiments and data analysis to determine individual biophysical parameters of If, is a promising tool for risk stratification of sinus bradycardia due to loss-of-function mutations in HCN4. In combination with an If blocker, this tool may also prove useful when applied to human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) obtained from mutation carriers and non-carriers.

Scanning mutagenesis of the voltage-gated sodium channel NaV1.2 using base editing

It is challenging to apply traditional mutational scanning to voltage-gated sodium channels (NaVs) and functionally annotate the large number of coding variants in these genes. Using a cytosine base editor and a pooled viability assay, we screen a library of 368 guide RNAs (gRNAs) tiling NaV1.2 to identify more than 100 gRNAs that change NaV1.2 function. We sequence base edits made by a subset of these gRNAs to confirm specific variants that drive changes in channel function. Electrophysiological characterization of these channel variants validates the screen results and provides functional mechanisms of channel perturbation. Most of the changes caused by these gRNAs are classifiable as loss of function along with two missense mutations that lead to gain of function in NaV1.2 channels. This two-tiered strategy to functionally characterize ion channel protein variants at scale identifies a large set of loss-of-function mutations in NaV1.2.

Contribution of Axon Initial Segment Structure and Channels to Brain Pathology

Brain channelopathies are a group of neurological disorders that result from genetic mutations affecting ion channels in the brain. Ion channels are specialized proteins that play a crucial role in the electrical activity of nerve cells by controlling the flow of ions such as sodium, potassium, and calcium. When these channels are not functioning properly, they can cause a wide range of neurological symptoms such as seizures, movement disorders, and cognitive impairment. In this context, the axon initial segment (AIS) is the site of action potential initiation in most neurons. This region is characterized by a high density of voltage-gated sodium channels (VGSCs), which are responsible for the rapid depolarization that occurs when the neuron is stimulated. The AIS is also enriched in other ion channels, such as potassium channels, that play a role in shaping the action potential waveform and determining the firing frequency of the neuron. In addition to ion channels, the AIS contains a complex cytoskeletal structure that helps to anchor the channels in place and regulate their function. Therefore, alterations in this complex structure of ion channels, scaffold proteins, and specialized cytoskeleton may also cause brain channelopathies not necessarily associated with ion channel mutations. This review will focus on how the AISs structure, plasticity, and composition alterations may generate changes in action potentials and neuronal dysfunction leading to brain diseases. AIS function alterations may be the consequence of voltage-gated ion channel mutations, but also may be due to ligand-activated channels and receptors and AIS structural and membrane proteins that support the function of voltage-gated ion channels.

Shared genetic loci and causal relations between schizophrenia and obsessive-compulsive disorder

Based on the clinical overlap between schizophrenia (SCZ) and obsessive-compulsive disorder (OCD), both disorders may share neurobiological substrates. In this study, we first analyzed recent large genome-wide associations studies (GWAS) on SCZ (n = 53,386, Psychiatric Genomics Consortium Wave 3) and OCD (n = 2688, the International Obsessive-Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and the OCD Collaborative Genetics Association Study (OCGAS)) using a conjunctional false discovery rate (FDR) approach to evaluate overlap in common genetic variants of European descent. Using a variety of biological resources, we functionally characterized the identified genomic loci. Then we used two-sample Mendelian randomization (MR) to estimate the bidirectional causal association between SCZ and OCD. Results showed that there is a positive genetic correlation between SCZ and OCD (r_g = 0.36, P = 0.02). We identified that one genetic locus (lead SNP rs5757717 in an intergenic region at CACNA1I) was jointly associated with SCZ and OCD (conjFDR = 2.12 × 10^-2). Mendelian randomization results showed that variants associated with increased risk for SCZ also increased the risk of OCD. This study broadens our understanding of the genetic architectures underpinning SCZ and OCD and suggests that the same molecular genetic processes may be responsible for shared pathophysiological and clinical characteristics between the two disorders.

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. As a consequence, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion-channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated. We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion-channel families. We collected and curated 3049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12 546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures. We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis centre and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain versus loss-of-channel function. In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional readouts and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.

Reversal of hyperactive higher-order thalamus attenuates defensiveness in a mouse model of PTSD

Posttraumatic stress disorder (PTSD) is a highly prevalent and debilitating psychiatric disease often accompanied by severe defensive behaviors, preventing individuals from integrating into society. However, the neural mechanisms of defensiveness in PTSD remain largely unknown. Here, we identified that the higher-order thalamus, the posteromedial complex of the thalamus (PoM), was overactivated in a mouse model of PTSD, and suppressing PoM activity alleviated excessive defensive behaviors. Moreover, we found that diminished thalamic inhibition derived from the thalamic reticular nucleus was the major cause of thalamic hyperactivity in PTSD mice. Overloaded thalamic innervation to the downstream cortical area, frontal association cortex, drove abnormal defensiveness. Overall, our study revealed that the malfunction of the higher-order thalamus mediates defensive behaviors and highlighted the thalamocortical circuit as a potential target for treating PTSD-related overreactivity symptoms.

Epigenome-wide DNA methylation in leukocytes and toenail metals: The normative aging study

BACKGROUND

Environmental metal exposures have been associated with multiple deleterious health endpoints. DNA methylation (DNAm) may provide insight into the mechanisms underlying these relationships. Toenail metals are non-invasive biomarkers, reflecting a medium-term time exposure window.

OBJECTIVES

This study examined variation in leukocyte DNAm and toenail arsenic (As), cadmium (Cd), lead (Pb), manganese (Mn), and mercury (Hg) among elderly men in the Normative Aging Study, a longitudinal cohort.

METHODS

We repeatedly collected samples of blood and toenail clippings. We measured DNAm in leukocytes with the Illumina HumanMethylation450 K BeadChip. We first performed median regression to evaluate the effects of each individual toenail metal on DNAm at three levels: individual cytosine-phosphate-guanine (CpG) sites, regions, and pathways. Then, we applied a Bayesian kernel machine regression (BKMR) to assess the joint and individual effects of metal mixtures on DNAm. Significant CpGs were identified using a multiple testing correction based on the independent degrees of freedom approach for correlated outcomes. The approach considers the effective degrees of freedom in the DNAm data using the principal components that explain >95% variation of the data.

RESULTS

We included 564 subjects (754 visits) between 1999 and 2013. The numbers of significantly differentially methylated CpG sites, regions, and pathways varied by metals. For example, we found six significant pathways for As, three for Cd, and one for Mn. The As-associated pathways were associated with cancer (e.g., skin cancer) and cardiovascular disease, whereas the Cd-associated pathways were related to lung cancer. Metal mixtures were also associated with 47 significant CpG sites, as well as pathways, mainly related to cancer and cardiovascular disease.

CONCLUSIONS

This study provides an approach to understanding the potential epigenetic mechanisms underlying observed relations between toenail metals and adverse health endpoints.

CaV3.3 Channelopathies

CaV3.3 is the third member of the low-voltage-activated calcium channel family and the last to be recognized as disease gene. Previously, CACNA1I, the gene encoding CaV3.3, had been described as schizophrenia risk gene. More recently, de novo missense mutations in CACNA1I were identified in patients with variable degrees of neurodevelopmental disease with and without epilepsy. Their functional characterization indicated gain-of-function effects resulting in increased calcium load and hyperexcitability of neurons expressing CaV3.3. The amino acids mutated in the CaV3.3 disease variants are located in the vicinity of the channel's activation gate and thus are classified as gate-modifying channelopathy mutations. A persistent calcium leak during rest and prolonged calcium spikes due to increased voltage sensitivity of activation and slowed kinetics of channel inactivation, respectively, may be causal for the neurodevelopmental defects. The prominent expression of CaV3.3 in thalamic reticular nucleus neurons and its essential role in generating the rhythmic thalamocortical network activity are consistent with a role of the mutated channels in the etiology of epileptic seizures and thus suggest T-type channel blockers as a viable treatment option.

CACNA1C (CaV1.2) and other L-type calcium channels in the pathophysiology and treatment of psychiatric disorders: Advances from functional genomics and pharmacoepidemiology

A role for voltage-gated calcium channels (VGCCs) in psychiatric disorders has long been postulated as part of a broader involvement of intracellular calcium signalling. However, the data were inconclusive and hard to interpret. We review three areas of research that have markedly advanced the field. First, there is now robust genomic evidence that common variants in VGCC subunit genes, notably CACNA1C which encodes the L-type calcium channel (LTCC) Ca_V1.2 subunit, are trans-diagnostically associated with psychiatric disorders including schizophrenia and bipolar disorder. Rare variants in these genes also contribute to the risk. Second, pharmacoepidemiological evidence supports the possibility that calcium channel blockers, which target LTCCs, might have beneficial effects on the onset or course of these disorders. This is especially true for calcium channel blockers that are brain penetrant. Third, long-range sequencing is revealing the repertoire of full-length LTCC transcript isoforms. Many novel and abundant CACNA1C isoforms have been identified in human and mouse brain, including some which are enriched compared to heart or aorta, and predicted to encode channels with differing functional and pharmacological properties. These isoforms may contribute to the molecular mechanisms of genetic association to psychiatric disorders. They may also enable development of therapeutic agents that can preferentially target brain LTCC isoforms and be of potential value for psychiatric indications.

Investigation of CACNA1I Cav3.3 Dysfunction in Hemiplegic Migraine

Familial hemiplegic migraine (FHM) is a severe neurogenetic disorder for which three causal genes, CACNA1A, SCN1A, and ATP1A2, have been implicated. However, more than 80% of referred diagnostic cases of hemiplegic migraine (HM) are negative for exonic mutations in these known FHM genes, suggesting the involvement of other genes. Using whole-exome sequencing data from 187 mutation-negative HM cases, we identified rare variants in the CACNA1I gene encoding the T-type calcium channel Cav3.3. Burden testing of CACNA1I variants showed a statistically significant increase in allelic burden in the HM case group compared to gnomAD (OR = 2.30, P = 0.00005) and the UK Biobank (OR = 2.32, P = 0.0004) databases. Dysfunction in T-type calcium channels, including Cav3.3, has been implicated in a range of neurological conditions, suggesting a potential role in HM. Using patch-clamp electrophysiology, we compared the biophysical properties of five Cav3.3 variants (p.R111G, p.M128L, p.D302G, p.R307H, and p.Q1158H) to wild-type (WT) channels expressed in HEK293T cells. We observed numerous functional alterations across the channels with Cav3.3-Q1158H showing the greatest differences compared to WT channels, including reduced current density, right-shifted voltage dependence of activation and inactivation, and slower current kinetics. Interestingly, we also found significant differences in the conductance properties exhibited by the Cav3.3-R307H and -Q1158H variants compared to WT channels under conditions of acidosis and alkalosis. In light of these data, we suggest that rare variants in CACNA1I may contribute to HM etiology.