Tryptophan-scanning mutagenesis in the S1 domain of mammalian HCN channel reveals residues critical for voltage-gated activation

Cannabidiol potentiates hyperpolarization-activated cyclic nucleotide-gated (HCN4) channels

Cannabidiol (CBD), the main non-psychotropic phytocannabinoid produced by the Cannabis sativa plant, blocks a variety of cardiac ion channels. We aimed to identify whether CBD regulated the cardiac pacemaker channel or the hyperpolarization-activated cyclic nucleotide-gated channel (HCN4). HCN4 channels are important for the generation of the action potential in the sinoatrial node of the heart and increased heart rate in response to β-adrenergic stimulation. HCN4 channels were expressed in HEK 293T cells, and the effect of CBD application was examined using a whole-cell patch clamp. We found that CBD depolarized the V1/2 of activation in holo-HCN4 channels, with an EC50 of 1.6 µM, without changing the current density. CBD also sped activation kinetics by approximately threefold. CBD potentiation of HCN4 channels occurred via binding to the closed state of the channel. We found that CBD's mechanism of action was distinct from cAMP, as CBD also potentiated apo-HCN4 channels. The addition of an exogenous PIP2 analog did not alter the ability of CBD to potentiate HCN4 channels, suggesting that CBD also acts using a unique mechanism from the known HCN4 potentiator PIP2. Lastly, to gain insight into CBD's mechanism of action, computational modeling and targeted mutagenesis were used to predict that CBD binds to a lipid-binding pocket at the C-terminus of the voltage sensor. CBD represents the first FDA-approved drug to potentiate HCN4 channels, and our findings suggest a novel starting point for drug development targeting HCN4 channels.

The effect of alterations of schizophrenia-associated genes on gamma band oscillations

Abnormalities in the synchronized oscillatory activity of neurons in general and, specifically in the gamma band, might play a crucial role in the pathophysiology of schizophrenia. While these changes in oscillatory activity have traditionally been linked to alterations at the synaptic level, we demonstrate here, using computational modeling, that common genetic variants of ion channels can contribute strongly to this effect. Our model of primary auditory cortex highlights multiple schizophrenia-associated genetic variants that reduce gamma power in an auditory steady-state response task. Furthermore, we show that combinations of several of these schizophrenia-associated variants can produce similar effects as the more traditionally considered synaptic changes. Overall, our study provides a mechanistic link between schizophrenia-associated common genetic variants, as identified by genome-wide association studies, and one of the most robust neurophysiological endophenotypes of schizophrenia.

Computational Modeling of Genetic Contributions to Excitability and Neural Coding in Layer V Pyramidal Cells: Applications to Schizophrenia Pathology

Pyramidal cells in layer V of the neocortex are one of the most widely studied neuron types in the mammalian brain. Due to their role as integrators of feedforward and cortical feedback inputs, they are well-positioned to contribute to the symptoms and pathology in mental disorders-such as schizophrenia-that are characterized by a mismatch between the internal perception and external inputs. In this modeling study, we analyze the input/output properties of layer V pyramidal cells and their sensitivity to modeled genetic variants in schizophrenia-associated genes. We show that the excitability of layer V pyramidal cells and the way they integrate inputs in space and time are altered by many types of variants in ion-channel and Ca²⁺ transporter-encoding genes that have been identified as risk genes by recent genome-wide association studies. We also show that the variability in the output patterns of spiking and Ca²⁺ transients in layer V pyramidal cells is altered by these model variants. Importantly, we show that many of the predicted effects are robust to noise and qualitatively similar across different computational models of layer V pyramidal cells. Our modeling framework reveals several aspects of single-neuron excitability that can be linked to known schizophrenia-related phenotypes and existing hypotheses on disease mechanisms. In particular, our models predict that single-cell steady-state firing rate is positively correlated with the coding capacity of the neuron and negatively correlated with the amplitude of a prepulse-mediated adaptation and sensitivity to coincidence of stimuli in the apical dendrite and the perisomatic region of a layer V pyramidal cell. These results help to uncover the voltage-gated ion-channel and Ca²⁺ transporter-associated genetic underpinnings of schizophrenia phenotypes and biomarkers.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Inhibition of hyperpolarization-activated cyclic nucleotide-gated channels by β-blocker carvedilol

BACKGROUND AND PURPOSE

Carvedilol is a clinically effective β-blocker broadly used for treating congestive heart failure (CHF), and several clinical trials have demonstrated that it shows a favourable effect compared with other β-blockers in patients with CHF. The mechanism underlying this beneficial effect of carvedilol compared to other β-blockers is not clearly understood. In addition to β-blockers, inhibitors of hyperpolarization-activated cyclic nucleotide (HCN)-gated channels, which play a critical role in spontaneous rhythmic activity in the heart, have also been proposed to be suitable drugs for reducing heart rate and, therefore, beneficial for treating CHF. In the present study, we investigated the effect of carvedilol on HCN channels.

EXPERIMENTAL APPROACH

Whole-cell patch-clamp recordings were used to assess the effect of carvedilol on currents from wild-type and mutant HCN1, HCN2 and HCN4 channels expressed in CHO cells.

KEY RESULTS

Carvedilol was the only β-blocker tested that showed inhibitory effects on the major sinoatrial HCN channel isoform HCN4. Carvedilol inhibited HCN4 in a concentration-dependent manner with an EC₅₀ of 4.4 μM. In addition, carvedilol also inhibited HCN1 and HCN2 channels. Carvedilol blocked HCN channels by decelerating the rate of channel activation and increasing that of deactivation, and shifted the voltage-dependence of activation leftwards. Our data also showed that carvedilol, unlike other inhibitors of this channel (ivabradine and ZD7288), is not an 'open-channel' inhibitor of HCN4.

CONCLUSIONS AND IMPLICATIONS

Carvedilol is a negative gating modulator of HCN channels. It represents a novel structure for future drug design of HCN channel inhibitors.

Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling

Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na⁺ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca²⁺ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.

A structural and functional comparison of gap junction channels composed of connexins and innexins

Methods such as electron microscopy and electrophysiology led to the understanding that gap junctions were dense arrays of channels connecting the intracellular environments within almost all animal tissues. The characteristics of gap junctions were remarkably similar in preparations from phylogenetically diverse animals such as cnidarians and chordates. Although few studies directly compared them, minor differences were noted between gap junctions of vertebrates and invertebrates. For instance, a slightly wider gap was noted between cells of invertebrates and the spacing between invertebrate channels was generally greater. Connexins were identified as the structural component of vertebrate junctions in the 1980s and innexins as the structural component of pre-chordate junctions in the 1990s. Despite a lack of similarity in gene sequence, connexins and innexins are remarkably similar. Innexins and connexins have the same membrane topology and form intercellular channels that play a variety of tissue- and temporally specific roles. Both protein types oligomerize to form large aqueous channels that allow the passage of ions and small metabolites and are regulated by factors such as pH, calcium, and voltage. Much more is currently known about the structure, function, and structure-function relationships of connexins. However, the innexin field is expanding. Greater knowledge of innexin channels will permit more detailed comparisons with their connexin-based counterparts, and provide insight into the ubiquitous yet specific roles of gap junctions. © 2016 Wiley Periodicals, Inc. Develop Neurobiol 77: 522-547, 2017.

Direct Interaction between the Voltage Sensors Produces Cooperative Sustained Deactivation in Voltage-gated H+ Channel Dimers

The voltage-gated H(+) channel (Hv) is a voltage sensor domain-like protein consisting of four transmembrane segments (S1-S4). The native Hv structure is a homodimer, with the two channel subunits functioning cooperatively. Here we show that the two voltage sensor S4 helices within the dimer directly cooperate via a π-stacking interaction between Trp residues at the middle of each segment. Scanning mutagenesis showed that Trp situated around the original position provides the slow gating kinetics characteristic of the dimer's cooperativity. Analyses of the Trp mutation on the dimeric and monomeric channel backgrounds and analyses with tandem channel constructs suggested that the two Trp residues within the dimer are functionally coupled during Hv deactivation but are less so during activation. Molecular dynamics simulation also showed direct π-stacking of the two Trp residues. These results provide new insight into the cooperative function of voltage-gated channels, where adjacent voltage sensor helices make direct physical contact and work as a single unit according to the gating process.

Functional Effects of Schizophrenia-Linked Genetic Variants on Intrinsic Single-Neuron Excitability: A Modeling Study

BACKGROUND

Recent genome-wide association studies have identified a large number of genetic risk factors for schizophrenia (SCZ) featuring ion channels and calcium transporters. For some of these risk factors, independent prior investigations have examined the effects of genetic alterations on the cellular electrical excitability and calcium homeostasis. In the present proof-of-concept study, we harnessed these experimental results for modeling of computational properties on layer V cortical pyramidal cells and identified possible common alterations in behavior across SCZ-related genes.

METHODS

We applied a biophysically detailed multicompartmental model to study the excitability of a layer V pyramidal cell. We reviewed the literature on functional genomics for variants of genes associated with SCZ and used changes in neuron model parameters to represent the effects of these variants.

RESULTS

We present and apply a framework for examining the effects of subtle single nucleotide polymorphisms in ion channel and calcium transporter-encoding genes on neuron excitability. Our analysis indicates that most of the considered SCZ-related genetic variants affect the spiking behavior and intracellular calcium dynamics resulting from summation of inputs across the dendritic tree.

CONCLUSIONS

Our results suggest that alteration in the ability of a single neuron to integrate the inputs and scale its excitability may constitute a fundamental mechanistic contributor to mental disease, alongside the previously proposed deficits in synaptic communication and network behavior.

Tryptophan Scanning Reveals Dense Packing of Connexin Transmembrane Domains in Gap Junction Channels Composed of Connexin32

Tryptophan was substituted for residues in all four transmembrane domains of connexin32. Function was assayed using dual cell two-electrode voltage clamp after expression in Xenopus oocytes. Tryptophan substitution was poorly tolerated in all domains, with the greatest impact in TM1 and TM4. For instance, in TM1, 15 substitutions were made, six abolished coupling and five others significantly reduced function. Only TM2 and TM3 included a distinct helical face that lacked sensitivity to tryptophan substitution. Results were visualized on a comparative model of Cx32 hemichannel. In this model, a region midway through the membrane appears highly sensitive to tryptophan substitution and includes residues Arg-32, Ile-33, Met-34, and Val-35. In the modeled channel, pore-facing regions of TM1 and TM2 were highly sensitive to tryptophan substitution, whereas the lipid-facing regions of TM3 and TM4 were variably tolerant. Residues facing a putative intracellular water pocket (the IC pocket) were also highly sensitive to tryptophan substitution. Although future studies will be required to separate trafficking-defective mutants from those that alter channel function, a subset of interactions important for voltage gating was identified. Interactions important for voltage gating occurred mainly in the mid-region of the channel and focused on TM1. To determine whether results could be extrapolated to other connexins, TM1 of Cx43 was scanned revealing similar but not identical sensitivity to TM1 of Cx32.

Emerging approaches to probing ion channel structure and function

Ion channels, as membrane proteins, are the sensors of the cell. They act as the first line of communication with the world beyond the plasma membrane and transduce changes in the external and internal environments into unique electrical signals to shape the responses of excitable cells. Because of their importance in cellular communication, ion channels have been intensively studied at the structural and functional levels. Here, we summarize the diverse approaches, including molecular and cellular, chemical, optical, biophysical, and computational, used to probe the structural and functional rearrangements that occur during channel activation (or sensitization), inactivation (or desensitization), and various forms of modulation. The emerging insights into the structure and function of ion channels by multidisciplinary approaches allow the development of new pharmacotherapies as well as new tools useful in controlling cellular activity.

Tryptophan scanning mutagenesis of the first transmembrane domain of the innexin Shaking-B(Lethal)

The channel proteins of gap junctions are encoded by two distinct gene families, connexins, which are exclusive to chordates, and innexins/pannexins, which are found throughout the animal kingdom. Although the relationship between the primary structure and function of the vertebrate connexins has been relatively well studied, there are, to our knowledge, no structure-function analyses of invertebrate innexins. In the first such study, we have used tryptophan scanning to probe the first transmembrane domain (M1) of the Drosophila innexin Shaking-B(Lethal), which is a component of rectifying electrical synapses in the Giant Fiber escape neural circuit. Tryptophan was substituted sequentially for 16 amino acids within M1 of Shaking-B(Lethal). Tryptophan insertion at every fourth residue (H27, T31, L35, and S39) disrupted gap junction function. The distribution of these sites is consistent with helical secondary structure and identifies the face of M1 involved in helix-helix interactions. Tryptophan substitution at several sites in M1 altered channel properties in a variety of ways. Changes in sensitivity to transjunctional voltage (Vj) were common and one mutation (S39W) induced sensitivity to transmembrane voltage (Vm). In addition, several mutations induced hemichannel activity. These changes are similar to those observed after substitutions within the transmembrane domains of connexins.

Peripheral N- and C-terminal domains determine deactivation kinetics of HCN channels

Among four subtypes of mammalian HCN channels, HCN1 has the fastest activation and deactivation kinetics while HCN4 shows the slowest. We previously showed that the activation kinetics are determined mainly by S1, S1-S2, and the S6-cyclic nucleotide binding domain. However, the effects of those regions on the deactivation kinetics were relatively small. Therefore, we investigated the structural basis for deactivation kinetics. Substitution of the core region (from S3 to S6) between HCN1 and HCN4 did not affect deactivation kinetics. This suggests that the peripheral regions (outside of S3 to S6) determine subtype-specific deactivation kinetics. Furthermore, we examined whether peripheral regions determined the deactivation kinetics across species by introducing the core region of DMIH (Drosophila homologue) into both HCN1 and HCN4. The DMIH core with HCN1 activated and deactivated more than threefold faster than that with HCN4. Taken together, the peripheral domains are diversified to create distinct kinetics.