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.

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.

Channel plasmon-polariton modes in V grooves filled with dielectric

We investigated the effect of dielectric filling in a V groove on the propagation parameters of channel plasmon-polariton (CPP) modes. In particular, existence conditions and critical groove angles, mode localization, field structure, dispersion, and propagation distances of CPP modes are analyzed as functions of dielectric permittivity inside the groove. It is demonstrated that increasing dielectric permittivity in the groove results in a rapid increase of mode localization near the tip of the groove and increase of both the critical angles that determine a range of groove angles for which CPP modes can exist. Detailed analysis of the field structure has demonstrated that the maximum of the field in a CPP mode is typically reached at a small distance from the tip of the groove. The effect of rounded tip is also investigated.

Stochastic models and simulation of ion channel dynamics

The behaviour of ion channels within cardiac and neuronal cells is intrinsically stochastic in nature. When the number of channels is small this stochastic noise is large and can have an impact on the dynamics of the system which is potentially an issue when modelling small neurons and drug block in cardiac cells. While exact methods correctly capture the stochastic dynamics of a system they are computationally expensive, restricting their inclusion into tissue level models and so approximations to exact methods are often used instead. The other issue in modelling ion channel dynamics is that the transition rates are voltage dependent, adding a level of complexity as the channel dynamics are coupled to the membrane potential. By assuming that such transition rates are constant over each time step, it is possible to derive a stochastic differential equation (SDE), in the same manner as for biochemical reaction networks, that describes the stochastic dynamics of ion channels. While such a model is more computationally efficient than exact methods we show that there are analytical problems with the resulting SDE as well as issues in using current numerical schemes to solve such an equation. We therefore make two contributions: develop a different model to describe the stochastic ion channel dynamics that analytically behaves in the correct manner and also discuss numerical methods that preserve the analytical properties of the model.

On the oscillatory behavior of transient Rayleigh Benard convection of air for 2D channel flow at moderate Rayleigh numbers

Unsteady numerical simulation of Rayleigh Benard convection heat transfer from a 2D channel is performed. The oscillatory behavior is attributed to recirculation of ascending and descending flows towards the core of the channel producing organized rolled motions. Variation of the parameters such as Reynolds number, channel outlet flow area and inclination of the channel are considered. Increasing Reynolds number (for a fixed Rayleigh number), delays the generation of vortices. The reduction in the outflow area leads to the later and the less vortex generation. As the time progresses, more vortices are generated, but the reinforced mean velocity does not let the eddies to enter the core of the channel. Therefore, they attach to the wall and reduce the heat transfer area. The inclination of the channel (both positive and negative) induces the generated vortices to get closer to each other and make an enlarged vortex.

Cerebral cortex activation mapping upon electrical muscle stimulation by 32-channel time-domain functional near-infrared spectroscopy

The application of different EMS current thresholds on muscle activates not only the muscle but also peripheral sensory axons that send proprioceptive and pain signals to the cerebral cortex. A 32-channel time-domain fNIRS instrument was employed to map regional cortical activities under varied EMS current intensities applied on the right wrist extensor muscle. Eight healthy volunteers underwent four EMS at different current thresholds based on their individual maximal tolerated intensity (MTI), i.e., 10 % < 50 % < 100 % < over 100 % MTI. Time courses of the absolute oxygenated and deoxygenated hemoglobin concentrations primarily over the bilateral sensorimotor cortical (SMC) regions were extrapolated, and cortical activation maps were determined by general linear model using the NIRS-SPM software. The stimulation-induced wrist extension paradigm significantly increased activation of the contralateral SMC region according to the EMS intensities, while the ipsilateral SMC region showed no significant changes. This could be due in part to a nociceptive response to the higher EMS current intensities and result also from increased sensorimotor integration in these cortical regions.

Comparison between electrically-evoked and voluntary wrist movements on sensorimotor and prefrontal cortical activation : A multi-channel time domain fNIRS study

Neuromuscular electrical stimulation (NMES) has been consistently demonstrated to improve skeletal muscle function in neurological populations with movement disorders, such as poststroke and incomplete spinal cord injury (Vanderthommen and Duchateau, 2007). Recent research has documented that rapid, supraspinal central nervous system reorganisation/neuroplastic mechanisms are also implicated during NMES (Chipchase et al., 2011). Functional neuroimaging studies have shown NMES to activate a network of sub-cortical and cortical brain regions, including the sensorimotor (SMC) and prefrontal (PFC) cortex (Blickenstorfer et al., 2009; Han et al., 2003; Muthalib et al., 2012). A relationship between increase in SMC activation with increasing NMES current intensity up to motor threshold has been previously reported using functional MRI (Smith et al., 2003). However, since clinical neurorehabilitation programmes commonly utilise NMES current intensities above the motor threshold and up to the maximum tolerated current intensity (MTI), limited research has determined the cortical correlates of increasing NMES current intensity at or above MTI (Muthalib et al., 2012). In our previous study (Muthalib et al., 2012), we assessed contralateral PFC activation using 1-channel functional near infrared spectroscopy (fNIRS) during NMES of the elbow flexors by increasing current intensity from motor threshold to greater than MTI and showed a linear relationship between NMES current intensity and the level of PFC activation. However, the relationship between NMES current intensity and activation of the motor cortical network, including SMC and PFC, has not been clarified. Moreover, it is of scientific and clinical relevance to know how NMES affects the central nervous system, especially in comparison to voluntary (VOL) muscle activation. Therefore, the aim of this study was to utilise multi-channel time domain fNIRS to compare SMC and PFC activation between VOL and NMESevoked wrist extension movements.

Two-channel waveguide modulator based on surface wave propagating in a semiconductor-metal thin-film structure

The results of theoretical investigations of two-channel waveguide modulator based on Surface Wave (SW) propagation are presented. The structure studied consists of two n-type semiconductor waveguide channels separated from each other by a dielectric gap and coated by a metal. The SW propagates at the semiconductor-metal interface across an external magnetic field which is parallel to the interface. An external dc voltage is applied to the metal surface of one channel to provide a small phase shift between two propagating modes. In a coupled mode approximation, two possible regimes of operation of the structure, namely as a directional coupler and as an electro-optical modulator, are considered. Our results suggest new applications in millimeter and submillimeter wave solid-state electronics and integrated optics.

UWB channel measurement and data transfer analysis for multiuser Infostation applications

In this paper, the results of the time dispersion parameters obtained from a set of channel measurements conducted in various environments that are typical of multiuser Infostation application scenarios are presented. The measurement procedure takes into account the practical scenarios typical of the positions and movements of the users in the particular Infostation network. To provide one with the knowledge of how much data can be downloaded by users over a given time and mobile speed, data transfer analysis for multiband orthogonal frequency division multiplexing (MB-OFDM) is presented. As expected, the rough estimate of simultaneous data transfer in a multiuser Infostation scenario indicates dependency of the percentage of download on the data size, number and speed of the users, and the elapse time.

Towards a single crystal Raman spectrum of kaolinite at 77 K

The Raman spectra at 77 K of the hydroxyl stretching of kaolinite were obtained along the three axes perpendicular to the crystal faces. Raman bands were observed at 3616, 3658 and 3677 cm−1 together with a distinct band observed at 3691 cm−1 and a broad profile between 3695 and 3715 cm−1. The band at 3616 cm−1 is assigned to the inner hydroxyl. The bands at 3658 and 3677 cm−1 are attributed to the out-of-phase vibrations of the inner surface hydroxyls. The Raman spectra of the in-phase vibrations of the inner-surface hydroxyl-stretching region are described in terms of transverse and longitudinal optic splitting. The band at 3691 cm−1 is assigned to the transverse optic and the broad profile to the longitudinal optic mode. This splitting remained even at liquid nitrogen temperature. The transverse optic vibration may be curve resolved into two or three bands, which are attributed to different types of hydroxyl groups in the kaolinite.

Micropolar flow in a porous channel with high mass transfer

Two dimensional flow of a micropolar fluid in a porous channel is investigated. The flow is driven by suction or injection at the channel walls, and the micropolar model due to Eringen is used to describe the working fluid. An extension of Berman's similarity transform is used to reduce the governing equations to a set of non-linear coupled ordinary differential equations. The latter are solved for large mass transfer via a perturbation analysis where the inverse of the cross-flow Reynolds number is used as the perturbing parameter. Complementary numerical solutions for strong injection are also obtained using a quasilinearisation scheme, and good agreement is observed between the solutions obtained from the perturbation analysis and the computations.