Functional analysis of missense variants in the TRESK (KCNK18) K+ channel

A loss of function mutation in the TRESK K2P potassium channel (KCNK18), has recently been linked with typical familial migraine with aura. We now report the functional characterisation of additional TRESK channel missense variants identified in unrelated patients. Several variants either had no apparent functional effect, or they caused a reduction in channel activity. However, the C110R variant was found to cause a complete loss of TRESK function, yet is present in both sporadic migraine and control cohorts, and no variation in KCNK18 copy number was found. Thus despite the previously identified association between loss of TRESK channel activity and migraine in a large multigenerational pedigree, this finding indicates that a single non-functional TRESK variant is not alone sufficient to cause typical migraine and highlights the genetic complexity of this disorder. Migraine is a common, disabling neurological disorder with a genetic, environmental and in some cases hormonal component. It is characterized by attacks of severe, usually unilateral and throbbing headache, can be accompanied by nausea, vomiting and photophobia and is clinically divided into two main subtypes, migraine with aura (MA) when a migraine is accompanied by transient and reversible focal neurological symptoms and migraine without aura (MO)1. The multifactorial and clinical heterogeneity of the disorder have considerably hindered the identification of common migraine susceptibility genes and most of our current understanding comes from the studies of familial hemiplegic migraine (FHM), a rare monogenic autosomal dominant form of MA2. So far, the three susceptibility genes that have been convincingly identified in FHM families all encode ion channels or transporters: CACNA1A encoding the α1 subunit of the Cav2.1 calcium channel3, SCN1A encoding the Nav1.1 sodium channel4 and ATP1A2 encoding the α2 subunit of the Na+/K+ pump5. It is believed that mutations in these genes may lead to increased efflux of glutamate and potassium in the synapse and thereby cause migraine by rendering the brain more susceptible to cortical spreading depression (CSD)6 which is thought to play a role in initiating a migraine attack7,8. However, these genes have not to date been implicated in common forms of migraine9. Nevertheless, current opinion suggests that typical migraine, like FHM, is also disorder of neuronal excitability, ion homeostasis and neurotransmitter release10,11,12. Mutations in the SLC4A4 gene encoding the sodium-bicarbonate cotransporter NBCe1, have recently been implicated in several different forms of migraine13, and a variety of genes involved in glutamate homeostasis (PGCP, MTDH14 and LRP115) and a cation channel (TRPM8)15 have also recently been implicated in migraine via genome-wide association studies. Ion channels are therefore highly likely to play an important role in the pathogenesis of typical migraine. TRESK (KCNK18), is a member of the two-pore domain (K2P) family of potassium channels involved in the control of cellular electrical excitability16. Regulation of TRESK activity by the calcium-dependent phosphatase calcineurin17, as well as its expression in dorsal root ganglia (DRG)18 and trigeminal ganglia (TG)19,20 has led to a proposed role for this channel in a variety of pain pathways. In a recent study, a frameshift mutation (F139Wfsx24) in TRESK was identified in a large multigenerational pedigree where it co-segregated perfectly with typical MA and a significant genome-wide linkage LOD score of 3.0. Furthermore, functional analysis revealed that this mutation caused a complete loss of TRESK function and that the truncated subunit was also capable of down regulating wild-type channel function. This therefore highlighted KCNK18 as potentially important candidate gene and suggested that TRESK dysfunction might play a possible role in the pathogenesis of familial migraine with visual aura20. Additional screening for KCNK18 mutations in unrelated sporadic migraine and control cohorts also identified a number of other missense variants; R10G, A34V, C110R, S231P and A233V20. The A233V variant was found only in the control cohort, whilst A34V was identified in a single Australian migraine proband for which family samples were not available, but it was not detected in controls. By contrast, the R10G, C110R, and S231P variants were found in both migraineurs and controls in both cohorts. In this study, we have investigated the functional effect of these variants to further probe the potential association of TRESK dysfunction with typical migraine.

Structures of μOconotoxins from Conusmarmoreus. Inhibitors of tetrodotoxin (TTX)-sensitive and TTX-resistant sodium channels in mammalian sensory neurons

The μO-conotoxins are an intriguing class of conotoxins targeting various voltage-dependent sodium channels and molluscan calcium channels. In the current study, we have shown MrVIA and MrVIB to be the first known peptidic inhibitors of the transient tetrodotoxin-resistant (TTX-R) Na+ current in rat dorsal root ganglion neurons, in addition to inhibiting tetrodotoxin-sensitive Na+ currents. Human TTX-R sodium channels are a therapeutic target for indications such as pain, highlighting the importance of the μO-conotoxins as potential leads for drug development. Furthermore, we have used NMR spectroscopy to provide the first structural information on this class of conotoxins. MrVIA and MrVIB are hydrophobic peptides that aggregate in aqueous solution but were solubilized in 50% acetonitrile/water. The three-dimensional structure of MrVIB consists of a small β-sheet and a cystine knot arrangement of the three-disulfide bonds. It contains four backbone “loops” between successive cysteine residues that are exposed to the solvent to varying degrees. The largest of these, loop 2, is the most disordered part of the molecule, most likely due to flexibility in solution. This disorder is the most striking difference between the structures of MrVIB and the known δ- and ω-conotoxins, which along with the μO-conotoxins are members of the O superfamily. Loop 2 of ω-conotoxins has previously been shown to contain residues critical for binding to voltage-gated calcium channels, and it is interesting to speculate that the flexibility observed in MrVIB may accommodate binding to both sodium and molluscan calcium channels.

Involvement of Ca2+-dependent hyperpolarization in sleep duration in mammals

The detailed molecular mechanisms underlying the regulation of sleep duration in mammals are still elusive. To address this challenge, we constructed a simple computational model, which recapitulates the electrophysiological characteristics of the slow-wave sleep and awake states. Comprehensive bifurcation analysis predicted that a Ca2+-dependent hyperpolarization pathway may play a role in slow-wave sleep and hence in the regulation of sleep duration. To experimentally validate the prediction, we generate and analyze 21 KO mice. Here we found that impaired Ca2+-dependent K+ channels (Kcnn2 and Kcnn3), voltage-gated Ca2+ channels (Cacna1g and Cacna1h), or Ca2+/calmodulin-dependent kinases (Camk2a and Camk2b) decrease sleep duration, while impaired plasma membrane Ca2+ ATPase (Atp2b3) increases sleep duration. Pharmacological intervention and whole-brain imaging validated that impaired NMDA receptors reduce sleep duration and directly increase the excitability of cells. Based on these results, we propose a hypothesis that a Ca2+-dependent hyperpolarization pathway underlies the regulation of sleep duration in mammals.

Gonadotropin-induced ovarian cancer cell migration and proliferation require extracellular signal-regulated kinase 1/2 activation regulated by calcium and protein kinase Cδ

The gonadotropin hypothesis proposes that elevated serum gonadotropin levels may increase the risk of epithelial ovarian cancer (EOC). We have studied the effect of treating EOC cell lines (OV207 and OVCAR-3) with FSH or LH. Both gonadotropins activated the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase 1/2 (ERK1/2) pathway and increased cell migration that was inhibited by the MAPK 1 inhibitor PD98059. Both extra- and intracellular calcium ion signalling were implicated in gonadotropin-induced ERK1/2 activation as treatment with either the calcium chelator EGTA or an inhibitor of intracellular calcium release, dantrolene, inhibited gonadotropin-induced ERK1/2 activation. Verapamil was also inhibitory, indicating that gonadotropins activate calcium influx via L-type voltage-dependent calcium channels. The cAMP/protein kinase A (PKA) pathway was not involved in the mediation of gonadotropin action in these cells as gonadotropins did not increase intracellular cAMP formation and inhibition of PKA did not affect gonadotropin-induced phosphorylation of ERK1/2. Activation of ERK1/2 was inhibited by the protein kinase C (PKC) inhibitor GF 109203X as well as by the PKCδ inhibitor rottlerin, and downregulation of PKCδ was inhibited by small interfering RNA (siRNA), highlighting the importance of PKCδ in the gonadotropin signalling cascade. Furthermore, in addition to inhibition by PD98059, gonadotropin-induced ovarian cancer cell migration was also inhibited by verapamil, GF 109203X and rottlerin. Similarly, gonadotropin-induced proliferation was inhibited by PD98059, verapamil, GF 109203X and PKCδ siRNA. Taken together, these results demonstrate that gonadotropins induce both ovarian cancer cell migration and proliferation by activation of ERK1/2 signalling in a calcium- and PKCδ-dependent manner.

Look before you leap : a confidence-based method for selecting species criticality while avoiding negative populations in τ-leaping

The stochastic simulation algorithm was introduced by Gillespie and in a different form by Kurtz. There have been many attempts at accelerating the algorithm without deviating from the behavior of the simulated system. The crux of the explicit τ-leaping procedure is the use of Poisson random variables to approximate the number of occurrences of each type of reaction event during a carefully selected time period, τ. This method is acceptable providing the leap condition, that no propensity function changes “significantly” during any time-step, is met. Using this method there is a possibility that species numbers can, artificially, become negative. Several recent papers have demonstrated methods that avoid this situation. One such method classifies, as critical, those reactions in danger of sending species populations negative. At most, one of these critical reactions is allowed to occur in the next time-step. We argue that the criticality of a reactant species and its dependent reaction channels should be related to the probability of the species number becoming negative. This way only reactions that, if fired, produce a high probability of driving a reactant population negative are labeled critical. The number of firings of more reaction channels can be approximated using Poisson random variables thus speeding up the simulation while maintaining the accuracy. In implementing this revised method of criticality selection we make use of the probability distribution from which the random variable describing the change in species number is drawn. We give several numerical examples to demonstrate the effectiveness of our new method.

Epistatic effects of potassium channel variation on cardiac repolarization and atrial fibrillation risk

Objectives The aim of this study was to evaluate the role of cardiac K+ channel gene variants in families with atrial fibrillation (AF). Background The K+ channels play a major role in atrial repolarization but single mutations in cardiac K+ channel genes are infrequently present in AF families. The collective effect of background K+ channel variants of varying prevalence and effect size on the atrial substrate for AF is largely unexplored. Methods Genes encoding the major cardiac K+ channels were resequenced in 80 AF probands. Nonsynonymous coding sequence variants identified in AF probands were evaluated in 240 control subjects. Novel variants were characterized using patch-clamp techniques and in silico modeling was performed using the Courtemanche atrial cell model. Results Nineteen nonsynonymous variants in 9 genes were found, including 11 rare variants. Rare variants were more frequent in AF probands (18.8% vs. 4.2%, p < 0.001), and the mean number of variants was greater (0.21 vs. 0.04, p < 0.001). The majority of K+ channel variants individually had modest functional effects. Modeling simulations to evaluate combinations of K+ channel variants of varying population frequency indicated that simultaneous small perturbations of multiple current densities had nonlinear interactions and could result in substantial (>30 ms) shortening or lengthening of action potential duration as well as increased dispersion of repolarization. Conclusions Families with AF show an excess of rare functional K+ channel gene variants of varying phenotypic effect size that may contribute to an atrial arrhythmogenic substrate. Atrial cell modeling is a useful tool to assess epistatic interactions between multiple variants.

Association analysis of a highly polymorphic CAG Repeat in the human potassium channel gene KCNN3 and migraine susceptibility

Background Migraine is a polygenic multifactorial disease, possessing environmental and genetic causative factors with multiple involved genes. Mutations in various ion channel genes are responsible for a number of neurological disorders. KCNN3 is a neuronal small conductance calcium-activated potassium channel gene that contains two polyglutamine tracts, encoded by polymorphic CAG repeats in the gene. This gene plays a critical role in determining the firing pattern of neurons and acts to regulate intracellular calcium channels. Methods The present association study tested whether length variations in the second (more 3') polymorphic CAG repeat in exon 1 of the KCNN3 gene, are involved in susceptibility to migraine with and without aura (MA and MO). In total 423 DNA samples from unrelated individuals, of which 202 consisted of migraine patients and 221 non-migraine controls, were genotyped and analysed using a fluorescence labelled primer set on an ABI310 Genetic Analyzer. Allele frequencies were calculated from observed genotype counts for the KCNN3 polymorphism. Analysis was performed using standard contingency table analysis, incorporating the chi-squared test of independence and CLUMP analysis. Results Overall, there was no convincing evidence that KCNN3 CAG lengths differ between Caucasian migraineurs and controls, with no significant difference in the allelic length distribution of CAG repeats between the population groups (P = 0.090). Also the MA and MO subtypes did not differ significantly between control allelic distributions (P > 0.05). The prevalence of the long CAG repeat (>19 repeats) did not reach statistical significance in migraineurs (P = 0.15), nor was there a significant difference between the MA and MO subgroups observed compared to controls (P = 0.46 and P = 0.09, respectively), or between MA vs MO (P = 0.40). Conclusion This association study provides no evidence that length variations of the second polyglutamine array in the N-terminus of the KCNN3 channel exert an effect in the pathogenesis of migraine.

Fluorescence detection of plant extracts that affect neuronal voltage-gated Ca2+ channels

Structurally novel compounds able to block voltage-gated Ca2+ channels (VGCCs) are currently being sought for the development of new drugs directed at neurological disorders. Fluorescence techniques have recently been developed to facilitate the analysis of VGCC blockers in a multi-well format. By utilising the small cell lung carcinoma cell line, NCI-H146, we were able to detect changes in intracellular Ca2+ concentration ([Ca2+]i) using a fluorescence microplate reader. NCI-H146 cells have characteristics resembling those of neuronal cells and express multiple VGCC subtypes, including those of the L-, N- and P-type. We found that K+-depolarisation of fluo-3 loaded NCI-H146 cells causes a rapid and transient increase in fluorescence, which was readily detected in a 96-well plate. Extracts of Australian plants, including those used traditionally as headache or pain treatments, were tested in this study to identify those affecting Ca2+ influx following membrane depolarisation of NCI-H146 cells. We found that E. bignoniiflora, A. symphyocarpa and E. vespertilio caused dose-dependent inhibition of K+-depolarised Ca2+ influx, with IC50 values calculated to be 234, 548 and 209 μg/ml, respectively. This data suggests an effect of these extracts on the function of VGCCs in these cells. Furthermore, we found similar effects using a fluorescence laser imaging plate reader (FLIPR) that allows simultaneous measurement of real-time fluorescence in a multi-well plate. Our results indicate that the dichloromethane extract of E. bignoniiflora and the methanolic extract of E. vespertilio show considerable promise as antagonists of neuronal VGCCs. Further analysis is required to characterise the function of the bioactive constituents in these extracts and determine their selectivity on VGCC subtypes.

Distamycin A affects the stability of NF-kappaB p50-DNA complexes in a sequence-dependent manner

The effect of two different DNA minor groove binding molecules, Hoechst 33258 and distamycin A, on the binding kinetics of NF-κB p50 to three different specific DNA sequences was studied at various salt concentrations. Distamycin A was shown to significantly increase the dissociation rate constant of p50 from the sequences PRDII (5′-GGGAAATTCC-3′) and Ig-κ B (5′-GGGACTTTCC-3′) but had a negligible effect on the dissociation from the palindromic target-κB binding site (5′-GGGAATTCCC-3′). By comparison, the effect of Hoechst 33258 on binding of p50 to each sequence was found to be minimal. The dissociation rates for the protein–DNA complexes increased at higher potassium chloride concentrations for the PRDII and Ig-κB binding motifs and this effect was magnified by distamycin A. In contrast, p50 bound to the palindromic target-κB site with a much higher intrinsic affinity and exhibited a significantly reduced salt dependence of binding over the ionic strength range studied, retaining a KD of less than 10 pM at 150 mM KCl. Our results demonstrate that the DNA binding kinetics of p50 and their salt dependence is strongly sequence-dependent and, in addition, that the binding of p50 to DNA can be influenced by the addition of minor groove-binding drugs in a sequence-dependent manner.

Induction of epithelial-mesenchymal transition (EMT) in breast cancer cells is calcium signal dependent

Signals from the tumor microenvironment trigger cancer cells to adopt an invasive phenotype through epithelial-mesenchymal transition (EMT). Relatively little is known regarding key signal transduction pathways that serve as cytosolic bridges between cell surface receptors and nuclear transcription factors to induce EMT. A better understanding of these early EMT events may identify potential targets for the control of metastasis. One rapid intracellular signaling pathway that has not yet been explored during EMT induction is calcium. Here we show that stimuli used to induce EMT produce a transient increase in cytosolic calcium levels in human breast cancer cells. Attenuation of the calcium signal by intracellular calcium chelation significantly reduced epidermal growth factor (EGF)- and hypoxia-induced EMT. Intracellular calcium chelation also inhibited EGF-induced activation of signal transducer and activator of transcription 3 (STAT3), while preserving other signal transduction pathways such as Akt and extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation. To identify calcium-permeable channels that may regulate EMT induction in breast cancer cells, we performed a targeted siRNA-based screen. We found that transient receptor potential-melastatin-like 7 (TRPM7) channel expression regulated EGF-induced STAT3 phosphorylation and expression of the EMT marker vimentin. Although intracellular calcium chelation almost completely blocked the induction of many EMT markers, including vimentin, Twist and N-cadherin, the effect of TRPM7 silencing was specific for vimentin protein expression and STAT3 phosphorylation. These results indicate that TRPM7 is a partial regulator of EMT in breast cancer cells, and that other calcium-permeable ion channels are also involved in calcium-dependent EMT induction. In summary, this work establishes an important role for the intracellular calcium signal in the induction of EMT in human breast cancer cells. Manipulation of calcium-signaling pathways controlling EMT induction in cancer cells may therefore be an important therapeutic strategy for preventing metastases.

Time-Dependent Density Functional Molecular Orbital and Excited State Calculations on Bis(porphyrinyl)butadiynes in the Monocationic, Neutral, Monoanionic, and Dianionic Oxidation States

We report a theoretical study of the multiple oxidation states (1+, 0, 1−, and 2−) of a meso,meso-linked diporphyrin, namely bis[10,15,20-triphenylporphyrinatozinc(II)-5-yl]butadiyne (4), using Time-Dependent Density Functional Theory (TDDFT). The origin of electronic transitions of singlet excited states is discussed in comparison to experimental spectra for the corresponding oxidation states of the close analogue bis{10,15,20-tris[3‘,5‘-di-tert-butylphenyl]porphyrinatozinc(II)-5-yl}butadiyne (3). The latter were measured in previous work under in situ spectroelectrochemical conditions. Excitation energies and orbital compositions of the excited states were obtained for these large delocalized aromatic radicals, which are unique examples of organic mixed-valence systems. The radical cations and anions of butadiyne-bridged diporphyrins such as 3 display characteristic electronic absorption bands in the near-IR region, which have been successfully predicted with use of these computational methods. The radicals are clearly of the “fully delocalized” or Class III type. The key spectral features of the neutral and dianionic states were also reproduced, although due to the large size of these molecules, quantitative agreement of energies with observations is not as good in the blue end of the visible region. The TDDFT calculations are largely in accord with a previous empirical model for the spectra, which was based simplistically on one-electron transitions among the eight key frontier orbitals of the C4 (1,4-butadiyne) linked diporphyrins.