S4 Helix Research Articles

Phosphoinositides regulation of hyperpolarization-activated cyclic-nucleotide gated (HCN) channels.

Several ion channels are regulated by protein-lipid interactions. Hyperpolarization-activated cyclic-nucleotide gated (HCN) channels specifically, have been shown to be regulated by phosphoinositides (PIPs). PIPs enhance HCN activation via a rightward shift in voltage-dependency and has importance applications in neuronal and cardiac function. Using computational and electrophysiological approaches, we aim to examine putative sites for HCN-PIP interactions. Computational docking and coarse-grained simulations indicate that PIP binding to HCN1 and HCN4 channels are not well coordinated, but rather occurs over a broad surface of charged residues primarily in the HCN-domain, S2 and S3 helices that can be loosely organized in 2 or 3 overlapping clusters. Thus, HCN-PIP interactions are more resembling of electrostatic interactions that occur in myristoylated alanine-rich C kinase substrate (MARCKS) proteins, than the specifically coordinated interactions that occur in pleckstrin homology domains (PH domains) or ion channels such as inward rectifier potassium (Kir) channels. Our results also indicate that phosphatidylinositol (PI) interactions with HCN1 are even lower affinity, explaining why unphosphorylated PI have no effect on HCN1 activation unlike phosphorylated PIPs. Experimental validation of putative binding sites are addressed by electrophysiological approaches.

Structural insights for the modulation of KCNQ1 channel by ML277.

KCNQ1 is a voltage-gated potassium ion channel protein which co-assembles with its accessory regulatory KCNE subunits, KCNE1-5, in order to fulfill its diverse physiological roles in various tissues. In the heart, KCNQ1 helps maintain a normal QT interval in the electrocardiogram, and loss of function mutations can lead to long QT syndrome type 1 (LQT1), seizures, and sudden death. Activators of KCNQ1 may be a promising strategy to oppose the effect of LQT mutations. In this work, we determined the cryo-EM structure of Xenopus KCNQ1 bound to Ca2+/calmodulin with or without the KCNQ1 activator, ML277. Our results showed that a binding site for ML277 lies in the transmembrane domain of the channel which is lined by residues of the S4-S5 linker, and the S5 and S6 helices of two adjacent subunits. MD simulations, along with point mutations, support the binding site location and showed that ML277 binds almost identically in open and closed channels. Comparative structural analyses highlighted important residues for the selectivity of ML277 towards KCNQ1, among other KCNQ subfamily members, and provided an explanation for the limited efficacy of ML277 in the presence of KCNE3. This work paves the way for a better understanding of KCNQ1 regulation by small molecules and may lead to the development of novel and potent activators of KCNQ1, providing a framework for therapeutic intervention targeting KCNQ1 loss of function pathologies.

Ionic interactions between gating charges and countercharges in voltage-sensing domain I independently regulate kinetics and voltage-dependence of Cav1.1 gating.

Voltage-gated calcium channels control a variety of processes in excitable cells. However, the exact mechanisms regulating kinetics and voltage-dependence of channel activation are not fully understood. Voltage-gated activation is determined by four distinct voltage-sensing domains (VSD I-IV) coupled to a common pore. Each VSD consists of four transmembrane helices (S1-S4) with S4 containing four to five gating charges. Upon membrane depolarization, consecutive interactions of these gating charges with negatively charged countercharges are believed to facilitate an upward movement of the S4 helices, leading to the opening of the pore and activation of the channel. Previous studies linked naturally occurring mutations of the innermost gating charge R4 (R174W) and its negative countercharge (E100) in CaV1.1 to muscle disease. To investigate the contribution of VSD I and the roles of its gating- and countercharges in channel gating, we combined structure-guided site-directed mutagenesis with patch-clamp analysis in dysgenic myotubes (CaV1.1-null). As E100 in helix S2 together with D126 in helix S3 form the highly conserved charge-transfer center of VSDs, we included this second negative countercharge (D126) in our analysis. While mutation of R174A resulted in a strong right-shift of voltage dependence, charge neutralizing mutations of the countercharges E100 and D126 both lead to a left-shift of voltage-dependence and to a slowing of channel activation. The double mutant of both countercharges resulted in an additive effect on activation kinetics, but left the voltage dependence of activation unaffected compared to the single mutants. Our findings indicate that VSD I substantially contributes to the regulation of both kinetics and voltage-dependence of activation; that the same ionic interactions are critical for both properties; but that the molecular mechanisms governing the two gating properties are functionally separate.

A mutational and modeling study of a binding pocket for IKs activators between KCNQ1 and KCNE1.

The IKs current formed by the co-assembly of KCNQ1 and KCNE1 plays an important role in cardiac repolarization. Loss-of-function mutations in KCNQ1 and KCNE1 are linked to LQT syndrome. Mefenamic acid (Mef) an NSAID, and DIDS a stilbene, are known to enhance IKs currents. The present study explores in detail the role of previously identified K41-KCNE1 in Mef action as well as uncovering the KCNQ1/KCNE1 binding site(s) of Mef and DIDS. In silico experiments of drug docking on a model of IKs channel and molecular dynamic (MD) simulations were used to map critical KCNQ1 and KCNE1 residues, followed by site-directed mutagenesis and functional analysis using whole cell patch clamp. Docking procedures and MD simulations show that Mef and DIDS bind between extracellular KCNE1 residues, specifically to the subunits of the external S6 helix; neighbouring S1; and the pore turret of the opposing subunit. Mef transforms the slowly activating WT IKs current (4:4 EQ) into a current with an almost linear waveform, which causes a left shift in activation V1/2 and inhibits tail current decay at −40 mV. These effects are prevented by the K41C-KCNE1 mutation in a stoichiometric manner. MD simulations and mutagenesis scanning identified three critical residues for Mef binding: W323-KCNQ1; Y46-KCNE1; and K41-KCNE1 in which Mef effects are greatly diminished. The effect of DIDS on WT-EQ are less pronounced than Mef, and K41C-KCNE1 was unable to fully prevent effects of DIDS. Our results show that both activators bind in the same general channel region but depend on different critical residues. These findings provide insights into a common drug-induced binding site, suggesting that the chemical structure of either compound could serve as a starting point for the development of IKs agonists.

PIP2 binding at the voltage sensor domain facilitates KCNQ1 VSD activation and gating.

The cardiac KCNQ1 channel is involved in cell hyperpolarisation, and loss-of-function mutations are associated with Long QT syndrome (LQTS), which can cause sudden cardiac arrest. The KCNQ1 channel is gated by the changes in transmembrane voltage and crucially requires the binding of two phosphatidylinositol-4,5-bisphosphate (PIP2) molecules to open the channel. The recent structure of KCNQ1 has revealed a PIP2 binding site at the voltage-sensor domain (VSD). In addition, previous modelling and electrophysiological studies have proposed another PIP2 binding site at the pore domain. Using enhanced sampling molecular dynamics simulation, we show that PIP2 binding at the VSD facilitates its activation, whereas PIP2 interaction of the pore facilitates VSD-pore coupling by straightening the S6 helix. Our simulation proposes that PIP2 binding at the VSD acts as a counter charge pausing point during the voltage-dependent activation process, lowering the energy barrier during the S4 translocation. Together, our work highlights the importance of PIP2 at the VSD as an activator to the VSD activation, in addition to our previous knowledge as the key coupler between the VSD and the pore.

Activation-pathway transitions in human voltage-gated proton channels revealed by a non-canonical fluorescent amino acid.

Voltage-dependent gating of the voltage-gated proton channels (HV1) remains poorly understood, partly because of the difficulty of obtaining direct measurements of voltage sensor movement in the form of gating currents. To circumvent this problem, we have implemented patch-clamp fluorometry in combination with the incorporation of the fluorescent non-canonical amino acid Anap to monitor channel opening and movement of the S4 segment. Simultaneous recording of currents and fluorescence signals allows for direct correlation of these parameters and investigation of their dependence on voltage and the pH gradient (ΔpH). We present data that indicate that Anap incorporated in the S4 helix is quenched by an aromatic residue located in the S2 helix and that motion of the S4 relative to this quencher is responsible for fluorescence increases upon depolarization. The kinetics of the fluorescence signal reveal the existence of a very slow transition in the deactivation pathway, which seems to be singularly regulated by ΔpH. Our experiments also suggest that the voltage sensor can move after channel opening and that the absolute value of the pH can influence the channel opening step. These results shed light on the complexities of voltage-dependent opening of human HV1 channels.

Functional Characterization of Four Known Cav2.1 Variants Associated with Neurodevelopmental Disorders.

Cav2.1 channels are expressed throughout the brain and are the predominant Ca2+ channels in the Purkinje cells. These cerebellar neurons fire spontaneously, and Cav2.1 channels are involved in the regular pacemaking activity. The loss of precision of the firing pattern of Purkinje cells leads to ataxia, a disorder characterized by poor balance and difficulties in performing coordinated movements. In this study, we aimed at characterizing functional and structural consequences of four variations (p.A405T in I-II loop and p.R1359W, p.R1667W and p.S1799L in IIIS4, IVS4, and IVS6 helices, respectively) identified in patients exhibiting a wide spectrum of disorders including ataxia symptoms. Functional analysis using two major Cav2.1 splice variants (Cav2.1+e47 and Cav2.1-e47) in Xenopus laevis oocytes, revealed a lack of effect upon A405T substitution and a significant loss-of-function caused by R1359W, whereas R1667W and S1799L caused both channel gain-of-function and loss-of-function, in a splice variant-dependent manner. Structural analysis revealed the loss of interactions with S1, S2, and S3 helices upon R1359W and R1667W substitutions, but a lack of obvious structural changes with S1799L. Computational modeling suggests that biophysical changes induced by Cav2.1 pathogenic mutations might affect action potential frequency in Purkinje cells.

An α-π transition in S6 shapes the conformational cycle of the bacterial sodium channel NavAb.

Voltage-gated sodium channels play an important role in electrical signaling in excitable cells. In response to changes in membrane potential, they cycle between nonconducting and conducting conformations. With recent advances in structural biology, structures of sodium channels have been captured in several distinct conformations, which are thought to represent different functional states. However, it has been difficult to capture the intrinsically transient open state. We recently showed that a proposed open state of the bacterial sodium channel NavMs was not conductive and that a conformational change involving a transition to a π-helix in the pore-lining S6 helix converted this structure into a conducting state. However, the relevance of this structural feature in other sodium channels, and its implications for the broader gating cycle, remained unclear. Here, we propose a comparable open state of another class of bacterial channel from Aliarcobacter butzleri (NavAb) with characteristic pore hydration, ion permeation, and drug binding properties. Furthermore, we show that a π-helix transition can lead to pore opening and that such a conformational change blocks fenestrations in the inner helix bundle. We also discover that a region in the C-terminal domain can undergo a disordering transition proposed to be important for pore opening. These results support a role for a π-helix transition in the opening of NavAb, enabling new proposals for the structural annotation and drug modulation mechanisms in this important sodium channel model.

Molecular Modeling of Cardiac Sodium Channel with Mexiletine.

A sodium channel blocker mexiletine (MEX) is used to treat chronic pain, myotonia and some arrhythmias. Mutations in the pore domain (PD) of voltage-gated sodium channels differently affect tonic block (TB) and use-dependent block (UDB) by MEX. Previous studies identified several MEX-sensing residues in the hNav1.5 channel and demonstrated that the channel block by MEX increases with activation of the voltage-sensing domain III (VSDIII), whereas MEX stabilizes the activated state of VSDIII. Structural rationales for these observations are unclear. Here, Monte Carlo (MC) energy minimizations were used to dock MEX and its more potent analog, Thio-Me2, into the hNav1.5 cryo-EM structure with activated VSDs and presumably inactivated PD. Computations yielded two ensembles of ligand binding poses in close contacts with known MEX-sensing residues in helices S6III, S6IV and P1IV. In both ensembles, the ligand NH3 group approached the cation-attractive site between backbone carbonyls at the outer-pore bottom, while the aromatic ring protruded ether into the inner pore (putative UDB pose) or into the III/IV fenestration (putative TB pose). In silico deactivation of VSDIII shifted helices S4-S5III, S5III, S6III and S6IV and tightened the TB site. In a model with activated VSDIII and three resting VSDs, MC-minimized energy profile of MEX pulled from the TB site towards lipids shows a deep local minimum due to interactions with 11 residues in S5III, P1III, S6III and S6IV. The minimum may correspond to an interim binding site for MEX in the hydrophobic path to the TB site along the lipid-exposed sides of repeats III and IV where 15 polar and aromatic residues would attract cationic blockers. The study explains numerous experimental data and suggests the mechanism of allosteric modification of the MEX binding site by VSDIII.

De novo KCNA6 variants with attenuated KV 1.6 channel deactivation in patients with epilepsy.

Mutations in the genes encoding neuronal ion channels are a common cause of Mendelian neurological diseases. We sought to identify novel de novo sequence variants in cases with early infantile epileptic phenotypes and neurodevelopmental anomalies. Following clinical diagnosis, we performed whole exome sequencing of the index cases and their parents. Identified channel variants were expressed in Xenopus oocytes and their functional properties assessed using two-electrode voltage clamp. We identified novel de novo variants in KCNA6 in four unrelated individuals variably affected with neurodevelopmental disorders and seizures with onset in the first year of life. Three of the four identified mutations affect the pore-lining S6 α-helix of KV 1.6. A prominent finding of functional characterization in Xenopus oocytes was that the channel variants showed only minor effects on channel activation but slowed channel closure and shifted the voltage dependence of deactivation in a hyperpolarizing direction. Channels with a mutation affecting the S6 helix display dominant effects on channel deactivation when co-expressed with wild-type KV 1.6 or KV 1.1 subunits. This is the first report of de novo nonsynonymous variants in KCNA6 associated with neurological or any clinical features. Channel variants showed a consistent effect on channel deactivation, slowing the rate of channel closure following normal activation. This specific gain-of-function feature is likely to underlie the neurological phenotype in our patients. Our data highlight KCNA6 as a novel channelopathy gene associated with early infantile epileptic phenotypes and neurodevelopmental anomalies.

Molecular Dynamics Simulations of the Cardiac Ryanodine Receptor Type 2 (RyR2) Gating Mechanism.

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca2+ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.

Abstract 12815: High Throughput Variant Identification in the S3-S5 Helices of KCNH2 (hERG)

Introduction: LQTS is a multi-factorial cardiac repolarization disorder characterized by a prolonged corrected QT interval (QTc) that predisposes individuals to sudden death. Variants in KCNH2 account for ~30% of LQTS cases. Variants in a ‘hotspot’ domain of KCNH2 are of concern, as variant location in KCNH2 is strongly associated with increased risk of cardiac events. Here, we determined the trafficking phenotype of all possible variants in the S3 and S4 helices (Voltage sensor) and S5 helix and p-loop (pore domain), both ‘hotspot’ domains, of KCNH2 . Method: We used a HA tagged KCNH2 construct containing mCherry to enable massively parallel characterization of channel trafficking of missense KCNH2 variants. We studied a total of 3315/3381 possible variants for residues 477-637. Results: Our preliminary data revealed a heterogenous distribution of trafficking defects (range: 0-157% of WT). Some regions (e.g. 562-566) were largely intolerant to substitution (trafficking scores ~0% for all variants) whereas other regions (e.g. 509-521) were very permissive of substitution (trafficking ranged from 78-146% of WT). For 51 clinically identified variants, we compared trafficking results with peak tail current density data and found that 4/6 were correctly classified as functionally normal and 42/44 correctly classified as loss of function. Conclusion: Massively parallel trafficking assay can identify KCNH2 defects and illuminate important channel biophysics. Massively parallel trafficking assays will enable prospective and accurate variant assessment for any missense variant in KCNH2 and subsequently other proteins where misfolding is a major loss-of-function mechanism.

Voltage-sensor movements in the Eag Kv channel under an applied electric field

Voltage-dependent ion channels regulate the opening of their pores by sensing the membrane voltage. This process underlies the propagation of action potentials and other forms of electrical activity in cells. The voltage dependence of these channels is governed by the transmembrane displacement of the positive charged S4 helix within their voltage-sensor domains. We use cryo-electron microscopy to visualize this movement in the mammalian Eag voltage-dependent potassium channel in lipid membrane vesicles with a voltage difference across the membrane. Multiple structural configurations show that the applied electric field displaces S4 toward the cytoplasm by two helical turns, resulting in an extended interfacial helix near the inner membrane leaflet. The position of S4 in this down conformation is sterically incompatible with an open pore, thus explaining how movement of the voltage sensor at hyperpolarizing membrane voltages locks the pore shut in this kind of voltage-dependent K+ (Kv) channel. The structures solved in lipid bilayer vesicles detail the intricate interplay between Kv channels and membranes, from showing how arginines are stabilized deep within the membrane and near phospholipid headgroups, to demonstrating how the channel reshapes the inner leaflet of the membrane itself.

Unexpected expansion of the voltage-gated proton channel family.

Voltage-gated ion channels, whose first identified function was to generate action potentials, are divided into subfamilies with numerous members. The family of voltage-gated proton channels (HV ) is tiny. To date, all species found to express HV have exclusively one gene that codes for this unique ion channel. Here we report the discovery and characterization of three proton channel genes in the classical model system of neural plasticity, Aplysia californica. The three channels (AcHV 1, AcHV 2, and AcHV 3) are distributed throughout the whole animal. Patch-clamp analysis confirmed proton selectivity of these channels but they all differed markedly in gating. AcHV 1 gating resembled HV in mammalian cells where it is responsible for proton extrusion and charge compensation. AcHV 2 activates more negatively and conducts extensive inward proton current, properties likely to acidify the cytosol. AcHV 3, which differs from AcHV 1 and AcHV 2 in lacking the first arginine in the S4 helix, exhibits proton selective leak currents and weak voltage dependence. We report the expansion of the proton channel family, demonstrating for the first time the expression of three functionally distinct proton channels in a single species.

Molecular Determinants for the High-Affinity Blockade of Human Ether-à-go-go-Related Gene K + Channel by Tolterodine.

Tolterodine is a first-line antimuscarinic drug used to treat overactive bladder. Adverse cardiac effects including tachycardia and palpitations have been observed, presumably because of its inhibition of the human ether-à-go-go-related gene (hERG) K + channel. However, the molecular mechanism of hERG channel inhibition by tolterodine is largely unclear. In this study, we performed molecular docking to identify potential binding sites of tolterodine in hERG channel, and two-microelectrode voltage-clamp to record the currents of hERG and its mutants expressed in Xenopus oocytes. The results of computational modeling demonstrated that phenylalanine at position 656 (F656) and tyrosine at position 652 (Y652) on the S6 helix of hERG channel are the most favorable binding residues of tolterodine, which was validated by electrophysiological recordings on Y652A and F656A hERG mutants. The Y652A and F656A mutations decreased inhibitory potency of tolterodine 345-fold and 126-fold, respectively. The Y652A mutation significantly altered the voltage dependence of channel inhibition by tolterodine. For both the wild-type and the mutant channels, tolterodine reduced the currents in a time-dependent manner, and the blockade occurred with the channel activated. Tolterodine did not interfere with hERG channel deactivation, whereas channel inactivation greatly impaired its blocking effect. The inhibition of hERG channel by tolterodine is independent of its action on muscarinic acetylcholine receptors. In conclusion, tolterodine is an open-state blocker of hERG K + channel with nanomolar potency. Y652 and F656, 2 aromatic residues on the inner S6 helix, are responsible for the high-affinity binding of tolterodine to hERG channel.

Scutellarein attenuates atopic dermatitis by selectively inhibiting transient receptor potential vanilloid 3 channels.

Atopic dermatitis (AD) is one of the most common chronic inflammatory cutaneous diseases with unmet clinical needs. As a common ingredient found in several medicinal herbs with efficacy on cutaneous inflammatory diseases, Scutellarein (Scu) has been shown to possess anti-inflammatory and anti-proliferative activities. We aimed to evaluate the therapeutic efficacy of Scu against AD and its underlying molecular mechanism. Efficacy of Scu on AD was evaluated in 2,4-dinitrofluorobenzene (DNFB) and carvacrol-induced dermatitis mouse models. Cytokine mRNA and serum IgE levels were examined using qPCR and ELISA, respectively. Voltage clamp recordings were used to measure currents mediated by transient receptor potential (TRP) channels. In silico docking, site-direct mutagenesis, and covalent modification were used to explore the binding pocket of Scu on TRPV3. Subcutaneous administration of Scu efficaciously suppresses DNFB and carvacrol-induced pruritus, epidermal hyperplasia and skin inflammation in wild type mice but has no additional benefit in Trpv3 knockout mice in the carvacrol model. Scu is a potent and selective TRPV3 channel allosteric negative modulator with an apparent affinity of 1.18 μM. Molecular docking coupled with site-direct mutagenesis and covalent modification of incorporated cysteine residues demonstrate that Scu targets the cavity formed between the pore helix and transmembrane helix S6. Moreover, Scu attenuates endogenous TRPV3 activity in human keratinocytes and inhibits carvacrol-induced proliferative and proinflammatory responses. Collectively, these data demonstrate that Scu ameliorates carvacrol-induced skin inflammation by directly inhibiting TRPV3, and TRPV3 represents a viable therapeutic target for AD treatment.

Vaccine-elicited murine antibody WS6 neutralizes diverse beta-coronaviruses by recognizing a helical stem supersite of vulnerability

Immunization with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike elicits diverse antibodies, but it is unclear if any of the antibodies can neutralize broadly against other beta-coronaviruses. Here, we report antibody WS6 from a mouse immunized with mRNA encoding the SARS-CoV-2 spike. WS6 bound diverse beta-coronavirus spikes and neutralized SARS-CoV-2 variants, SARS-CoV, and related sarbecoviruses. Epitope mapping revealed WS6 to target a region in the S2 subunit, which was conserved among SARS-CoV-2, Middle East respiratory syndrome (MERS)-CoV, and hCoV-OC43. The crystal structure at 2 Å resolution of WS6 revealed recognition to center on a conserved S2 helix, which was occluded in both pre- and post-fusion spike conformations. Structural and neutralization analyses indicated WS6 to neutralize by inhibiting fusion and post-viral attachment. Comparison of WS6 with other recently identified antibodies that broadly neutralize beta-coronaviruses indicated a stem-helical supersite—centered on hydrophobic residues Phe1148, Leu1152, Tyr1155, and Phe1156—to be a promising target for vaccine design.

Structural and electrophysiological basis for the modulation of KCNQ1 channel currents by ML277

The KCNQ1 ion channel plays critical physiological roles in electrical excitability and K+ recycling in organs including the heart, brain, and gut. Loss of function is relatively common and can cause sudden arrhythmic death, sudden infant death, epilepsy and deafness. Here, we report cryogenic electron microscopic (cryo-EM) structures of Xenopus KCNQ1 bound to Ca2+/Calmodulin, with and without the KCNQ1 channel activator, ML277. A single binding site for ML277 was identified, localized to a pocket lined by the S4-S5 linker, S5 and S6 helices of two separate subunits. Several pocket residues are not conserved in other KCNQ isoforms, explaining specificity. MD simulations and point mutations support this binding location for ML277 in open and closed channels and reveal that prevention of inactivation is an important component of the activator effect. Our work provides direction for therapeutic intervention targeting KCNQ1 loss of function pathologies including long QT interval syndrome and seizures.

Common Structural Pattern for Flecainide Binding in Atrial-Selective Kv1.5 and Nav1.5 Channels: A Computational Approach.

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Its treatment includes antiarrhythmic drugs (AADs) to modulate the function of cardiac ion channels. However, AADs have been limited by proarrhythmic effects, non-cardiovascular toxicities as well as often modest antiarrhythmic efficacy. Theoretical models showed that a combined blockade of Nav1.5 (and its current, INa) and Kv1.5 (and its current, IKur) ion channels yield a synergistic anti-arrhythmic effect without alterations in ventricles. We focused on Kv1.5 and Nav1.5 to search for structural similarities in their binding site (BS) for flecainide (a common blocker and widely prescribed AAD) as a first step for prospective rational multi-target directed ligand (MTDL) design strategies. We present a computational workflow for a flecainide BS comparison in a flecainide-Kv1.5 docking model and a solved structure of the flecainide-Nav1.5 complex. The workflow includes docking, molecular dynamics, BS characterization and pattern matching. We identified a common structural pattern in flecainide BS for these channels. The latter belongs to the central cavity and consists of a hydrophobic patch and a polar region, involving residues from the S6 helix and P-loop. Since the rational MTDL design for AF is still incipient, our findings could advance multi-target atrial-selective strategies for AF treatment.

Differential Activation of TRPM8 by the Stereoisomers of Menthol.

The stereoisomers of menthol elicit cooling sensation to various levels. Though the high-resolution three-dimensional structures of the menthol receptor, the transient receptor potential melastatin 8 (TRPM8) ion channels, have been revolved in different states, the menthol-bound state structure is not determined and how the stereoisomers of menthol interact with TRPM8 remains largely elusive. Taking advantage of the identical atom composition but distinct spatial orientation of chemical groups in menthol stereoisomers, we performed thermodynamic mutant cycle analysis (TMCA) with patch-clamp recordings to probe the interaction between these ligands and TRPM8. By comparing (−)-menthol with (+)-neoisomenthol or (+)-neomenthol, we observed that the isopropyl or hydroxyl group in menthol interacts with the S4 or S3 helix in TRPM8, respectively. These interactions were also corroborated in our molecular docking of the stereoisomers, though the predicted structural details in the interactions of these ligands with TRPM8 residues are different. Therefore, we suggest similar molecular mechanisms of TRPM8 activation by the stereoisomers of menthol, while the binding configuration of an individual stereoisomer is varied.