Category: Publications

Published 15 July 2022 in Nature Communications (doi 10.1038/s41467-022-29594-w):

Cryo-EM structure of the human Kv3.1 channel reveals gating control by the cytoplasmic T1 domain

Gamma Chi, Qiansheng Liang, Akshay Sridhar, John B Cowgill, Kasim Sader, Mazdak Radjainia, Pu Qian, Pablo Castro-Hartmann, Shayla Venkaya, Nanki Kaur Singh, Gavin McKinley, Alejandra Fernandez-Cid, Shubhashish M M Mukhopadhyay, Nicola A Burgess-Brown, Lucie Delemotte, Manuel Covarrubias, Katharina L Dürr

Kv3 channels have distinctive gating kinetics tailored for rapid repolarization in fast-spiking neurons. Malfunction of this process due to genetic variants in the KCNC1 gene causes severe epileptic disorders, yet the structural determinants for the unusual gating properties remain elusive. Here, we present cryo-electron microscopy structures of the human Kv3.1a channel, revealing a unique arrangement of the cytoplasmic tetramerization domain T1 which facilitates interactions with C-terminal axonal targeting motif and key components of the gating machinery. Additional interactions between S1/S2 linker and turret domain strengthen the interface between voltage sensor and pore domain. Supported by molecular dynamics simulations, electrophysiological and mutational analyses, we identify several residues in the S4/S5 linker which influence the gating kinetics and an electrostatic interaction between acidic residues in α6 of T1 and R449 in the pore-flanking S6T helices. These findings provide insights into gating control and disease mechanisms and may guide strategies for the design of pharmaceutical drugs targeting Kv3 channels.

Read the full publication here.

Published 1 June 2022 in eLife (doi 10.7554/eLife.77672):

Subtype-specific responses of hKv7.4 and hKv7.5 channels to polyunsaturated fatty acids reveal an unconventional modulatory site and mechanism

Damon JA Frampton, Koushik Choudhury, Johan Nikesjö, Lucie Delemotte, Sara I Liin

The KV7.4 and KV7.5 subtypes of voltage-gated potassium channels play a role in important physiological processes such as sound amplification in the cochlea and adjusting vascular smooth muscle tone. Therefore, the mechanisms that regulate KV7.4 and KV7.5 channel function are of interest. Here, we study the effect of polyunsaturated fatty acids (PUFAs) on human KV7.4 and KV7.5 channels expressed in Xenopus oocytes. We report that PUFAs facilitate activation of hKV7.5 by shifting the V50 of the conductance versus voltage (G(V)) curve toward more negative voltages. This response depends on the head group charge, as an uncharged PUFA analogue has no effect and a positively charged PUFA analogue induces positive V50 shifts. In contrast, PUFAs inhibit activation of hKV7.4 by shifting V50 toward more positive voltages. No effect on V50 of hKV7.4 is observed by an uncharged or a positively charged PUFA analogue. Thus, the hKV7.5 channel’s response to PUFAs is analogous to the one previously observed in hKV7.1–7.3 channels, whereas the hKV7.4 channel response is opposite, revealing subtype-specific responses to PUFAs. We identify a unique inner PUFA interaction site in the voltage-sensing domain of hKV7.4 underlying the PUFA response, revealing an unconventional mechanism of modulation of hKV7.4 by PUFAs.

Read the full publication here.

Published 6 April 2022 in Journal of Fluid Mechanics (doi 10.1017/jfm.2022.219):

Nanoscale sheared droplet: volume-of-fluid, phase-field and no-slip molecular dynamics

Uǧis Lācis, Michele Pellegrino, Johan Sundin, Gustav Amberg, Stéphane Zaleski, Berk Hess, Shervin Bagheri

The motion of the three-phase contact line between two immiscible fluids and a solid surface arises in a variety of wetting phenomena and technological applications. One challenge in continuum theory is the effective representation of molecular motion close to the contact line. Here, we characterize the molecular processes of the moving contact line to assess the accuracy of two different continuum two-phase models. Specifically, molecular dynamics simulations of a two-dimensional droplet between two moving plates are used to create reference data for different capillary numbers and contact angles. We use a simple-point-charge/extended water model. This model provides a very small slip and a more realistic representation of the molecular physics than Lennard-Jones models. The Cahn–Hilliard phase-field model and the volume-of-fluid model are calibrated against the drop displacement from molecular dynamics reference data. It is shown that the calibrated continuum models can accurately capture droplet displacement and droplet break-up for different capillary numbers and contact angles. However, we also observe differences between continuum and atomistic simulations in describing the transient and unsteady droplet behaviour, in particular, close to dynamical wetting transitions. The molecular dynamics of the sheared droplet provide insight into the line friction experienced by the advancing and receding contact lines. The presented results will serve as a stepping stone towards developing accurate continuum models for nanoscale hydrodynamics.

Read the full publication here.

Published 23 February 2022 in eLife (doi 10.7554/eLife.74773):

Identification of electroporation sites in the complex lipid organization of the plasma membrane

Lea Rems, Xinru Tang, Fangwei Zhao, Sergio Pérez-Conesa, Ilaria Testa, Lucie Delemotte

The plasma membrane of a biological cell is a complex assembly of lipids and membrane proteins, which tightly regulate transmembrane transport. When a cell is exposed to strong electric field, the membrane integrity becomes transiently disrupted by formation of transmembrane pores. This phenomenon termed electroporation is already utilized in many rapidly developing applications in medicine including gene therapy, cancer treatment, and treatment of cardiac arrhythmias. However, the molecular mechanisms of electroporation are not yet sufficiently well understood; in particular, it is unclear where exactly pores form in the complex organization of the plasma membrane. In this study, we combine coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis to identify how formation of pores depends on the local lipid organization. We show that pores do not form homogeneously across the membrane, but colocalize with domains that have specific features, the most important being high density of polyunsaturated lipids. We further show that knowing the lipid organization is sufficient to reliably predict poration sites with machine learning. Additionally, by analysing poration kinetics with Bayesian survival analysis we show that poration does not depend solely on local lipid arrangement, but also on membrane mechanical properties and the polarity of the electric field. Finally, we discuss how the combination of atomistic and coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis can guide the design of future experiments and help us to develop an accurate description of plasma membrane electroporation on the whole-cell level. Achieving this will allow us to shift the optimization of electroporation applications from blind trial-and-error approaches to mechanistic-driven design.

Read the full publication here.

Published 4 January 2022 in Biophysical Journal (doi 10.1016/j.bpj.2021.12.010):

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Koushik Choudhury, Marina A Kasimova, Sarah McComas, Rebecca J Howard, Lucie Delemotte

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Read the full publication here.

Published December 2021 in Journal of Biological Chemistry (doi 10.1016/j.jbc.2021.101355):

A missense mutation converts the Na⁺,K⁺-ATPase into an ion channel and causes therapy-resistant epilepsy

Sofia Ygberg, Evgeny E Akkuratov*, Rebecca J Howard*, Fulya Taylan, Daniel C Jans*, Dhani R Mahato, Adriana Katz, Paula F Kinoshita, Benjamin Portal, Inger Nennesmo, Maria Lindskog, Steven JD Karlish, Magnus Andersson*, Anna Lindstrand, Hjalmar Brismar*, Anita Aperia*

*Current and former members of the SciLifeLab biophysics community

The ion pump Na⁺,K⁺-ATPase is a critical determinant of neuronal excitability; however, its role in the etiology of diseases of the central nervous system (CNS) is largely unknown. We describe here the molecular phenotype of a Trp931Arg mutation of the Na⁺,K⁺-ATPase catalytic α1 subunit in an infant diagnosed with therapy-resistant lethal epilepsy. In addition to the pathological CNS phenotype, we also detected renal wasting of Mg²⁺. We found that membrane expression of the mutant α1 protein was low, and ion pumping activity was lost. Arginine insertion into membrane proteins can generate water-filled pores in the plasma membrane, and our molecular dynamic (MD) simulations of the principle states of Na⁺,K⁺-ATPase transport demonstrated massive water inflow into mutant α1 and destabilization of the ion-binding sites. MD simulations also indicated that a water pathway was created between the mutant arginine residue and the cytoplasm, and analysis of oocytes expressing mutant α1 detected a nonspecific cation current. Finally, neurons expressing mutant α1 were observed to be depolarized compared with neurons expressing wild-type protein, compatible with a lowered threshold for epileptic seizures. The results imply that Na⁺,K⁺-ATPase should be considered a neuronal locus minoris resistentia in diseases associated with epilepsy and with loss of plasma membrane integrity.

Read the full publication here.

Published 15 November 2021 in Journal of Chemical Information & Modeling (doi 10.1021/acs.jcim.1c00826):

Allosteric effect of nanobody binding on ligand-specific active states of the β2 adrenergic receptor

Yue Chen, Oliver Fleetwood, Sergio Pérez-Conesa, Lucie Delemotte

Nanobody binding stabilizes G-protein-coupled receptors (GPCR) in a fully active state and modulates their affinity for bound ligands. However, the atomic-level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes structural features of the β2AR with an ∼14 Å outward movement of transmembrane helix 6 and a close proximity of transmembrane (TM) helices 5 and 7, and favors the fully active-like conformation of the receptor, independent of ligand binding, in contrast to the conditions under which no intracellular binding partner is bound, in which case the receptor is only stabilized in an intermediate-active state. This activation is supported by the residues located at hotspots located on TMs 5, 6, and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intracellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding triggers tighter and stronger local communication networks between the Nb80 and the ligand-binding sites, primarily involving residues around ICL2 and the intracellular end of TM3, TM5, TM6, as well as ECL2, ECL3, and the extracellular ends of TM6 and TM7. In particular, we identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in conformations modulated by a full agonist, BI167107, and a G-protein-biased partial agonist, salmeterol, involving mainly TM1 and TM2, and TM5, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which should be taken into account in structure-based drug discovery.

Read the full publication here.

Published 15 October 2021 in eLife (doi 10.7554/eLife.68369):

Markov state models of proton- and pore-dependent activation in a pentameric ligand-gated ion channel

Cathrine Bergh, Stephanie A Heusser, Rebecca Howard, Erik Lindahl

Ligand-gated ion channels conduct currents in response to chemical stimuli, mediating electrochemical signaling in neurons and other excitable cells. For many channels the details of gating remain unclear, partly due to limited structural data and simulation timescales. Here, we used enhanced sampling to simulate the pH-gated channel GLIC, and construct Markov state models (MSMs) of gating. Consistent with new functional recordings we report in oocytes, our analysis revealed differential effects of protonation and mutation on free-energy wells. Clustering of closed- versus open-like states enabled estimation of open probabilities and transition rates, while higher-order clustering affirmed conformational trends in gating. Furthermore, our models uncovered state- and protonation-dependent symmetrization. This demonstrates the applicability of MSMs to map energetic and conformational transitions between ion-channel functional states, and how they reproduce shifts upon activation or mutation, with implications for modeling neuronal function and developing state-selective drugs.

Read the full publication here.

Published 14 September 2021 in Proceedings of the National Academy of Sciences of the USA (doi 10.1073/pnas.2108006118):

Probing solution structure of the pentameric ligand-gated ion channel GLIC by small-angle neutron scattering

Marie Lycksell, Urška Rovšnik, Cathrine Bergh, Nicolai T. Johansen, Anne Martel, Lionel Porcar, Lise Arleth, Rebecca J Howard, Erik Lindahl

Pentameric ligand-gated ion channels undergo subtle conformational cycling to control electrochemical signal transduction in many kingdoms of life. Several crystal structures have now been reported in this family, but the functional relevance of such models remains unclear. Here, we used small-angle neutron scattering (SANS) to probe ambient solution-phase properties of the pH-gated bacterial ion channel GLIC under resting and activating conditions. Data collection was optimized by inline paused-flow size-exclusion chromatography, and exchanging into deuterated detergent to hide the micelle contribution. Resting-state GLIC was the best-fit crystal structure to SANS curves, with no evidence for divergent mechanisms. Moreover, enhanced-sampling molecular-dynamics simulations enabled differential modeling in resting versus activating conditions, with the latter corresponding to an intermediate ensemble of both the extracellular and transmembrane domains. This work demonstrates state-dependent changes in a pentameric ion channel by SANS, an increasingly accessible method for macromolecular characterization with the coming generation of neutron sources.

Read the full publication here.

Published online August 31, 2021 in Journal of Investigative Dermatology (doi https://doi.org/10.1016/j.jid.2021.06.037):

The Skin’s Barrier: A Cryo-EM Based Overview of its Architecture and Stepwise Formation

Lars Norlén, Magnus Lundborg, Christian Wennberg, Ali Narangifard, Bertil Daneholt

A major role of the skin is to serve as a barrier toward the environment. The skin’s permeability barrier consists of a lipid structure positioned in the stratum corneum. Recent progress in high-resolution cryo-electron microscopy (cryo-EM) has allowed for elucidation of the architecture of the skin’s barrier and its stepwise formation process representing the final stage of epidermal differentiation. In this review, we present an overview of the skin’s barrier structure and its formation process, as evidenced by cryo-EM.

Read the full publication here.

Published 17 August 2021 in Proceedings of the National Academy of Sciences of the USA (doi 10.1073/pnas.2025320118):

Functional cross-talk between phosphorylation and disease-causing mutations in the cardiac sodium channel Na_v1.5

Iacopo Galleano, Hendrik Harms, Koushik Choudhury, Keith Khoo, Lucie Delemotte, Stephan A Pless

The voltage-gated sodium channel Na_v1.5 initiates the cardiac action potential. Alterations of its activation and inactivation properties due to mutations can cause severe, life-threatening arrhythmias. Yet despite intensive research efforts, many functional aspects of this cardiac channel remain poorly understood. For instance, Na_v1.5 undergoes extensive posttranslational modification in vivo, but the functional significance of these modifications is largely unexplored, especially under pathological conditions. This is because most conventional approaches are unable to insert metabolically stable posttranslational modification mimics, thus preventing a precise elucidation of the contribution by these modifications to channel function. Here, we overcome this limitation by using protein semisynthesis of Na_v1.5 in live cells and carry out complementary molecular dynamics simulations. We introduce metabolically stable phosphorylation mimics on both wild-type (WT) and two pathogenic long-QT mutant channel backgrounds and decipher functional and pharmacological effects with unique precision. We elucidate the mechanism by which phosphorylation of Y1495 impairs steady-state inactivation in WT Na_v1.5. Surprisingly, we find that while the Q1476R patient mutation does not affect inactivation on its own, it enhances the impairment of steady-state inactivation caused by phosphorylation of Y1495 through enhanced unbinding of the inactivation particle. We also show that both phosphorylation and patient mutations can impact Na_v1.5 sensitivity toward the clinically used antiarrhythmic drugs quinidine and ranolazine, but not flecainide. The data highlight that functional effects of Na_v1.5 phosphorylation can be dramatically amplified by patient mutations. Our work is thus likely to have implications for the interpretation of mutational phenotypes and the design of future drug regimens.

Read the full publication here.

Available online 1 July 2021 in Life Science Alliance (doi 10.26508/lsa.202101011):

Dynamic closed states of a ligand-gated ion channel captured by cryo-EM and simulations

Urška Rovšnik, Yuxuan Zhuang, Björn O Forsberg, Marta Carroni, Linnea Yvonnesdotter, Rebecca J Howard, Erik Lindahl 

Ligand-gated ion channels are critical mediators of electrochemical signal transduction across evolution. Biophysical and pharmacological characterization of these receptor proteins relies on high-quality structures in multiple, subtly distinct functional states. However, structural data in this family remain limited, particularly for resting and intermediate states on the activation pathway. Here, we report cryo-electron microscopy (cryo-EM) structures of the proton-activated Gloeobacter violaceus ligand-gated ion channel (GLIC) under three pH conditions. Decreased pH was associated with improved resolution and side chain rearrangements at the subunit/domain interface, particularly involving functionally important residues in the β1–β2 and M2–M3 loops. Molecular dynamics simulations substantiated flexibility in the closed-channel extracellular domains relative to the transmembrane ones and supported electrostatic remodeling around E35 and E243 in proton-induced gating. Exploration of secondary cryo-EM classes further indicated a low-pH population with an expanded pore. These results allow us to define distinct protonation and activation steps in pH-stimulated conformational cycling in GLIC, including interfacial rearrangements largely conserved in the pentameric channel family.

Read the full publication here.

Available online 3 July 2021 in the Journal of Molecular Biology (doi 10.1016/j.jmb.2021.167128):

Elephants in the dark: insights and incongruities in pentameric ligand-gated ion channel models

Rebecca J Howard

The superfamily of pentameric ligand-gated ion channels (pLGICs) comprises key players in electrochemical signal transduction across evolution, including historic model systems for receptor allostery and targets for drug development. Accordingly, structural studies of these channels have steadily increased, and now approach 250 depositions in the protein data bank. This review contextualizes currently available structures in the pLGIC family, focusing on morphology, ligand binding, and gating in three model subfamilies: the prokaryotic channel GLIC, the cation-selective nicotinic acetylcholine receptor, and the anion-selective glycine receptor. Common themes include the challenging process of capturing and annotating channels in distinct functional states; partially conserved gating mechanisms, including remodeling at the extracellular/transmembrane-domain interface; and diversity beyond the protein level, arising from posttranslational modifications, ligands, lipids, and signaling partners. Interpreting pLGIC structures can be compared to describing an elephant in the dark, relying on touch alone to comprehend the many parts of a monumental beast: each structure represents a snapshot in time under specific experimental conditions, which must be understood and integrated with further structure, function, and simulations data to build a comprehensive model, and understand how one channel may fundamentally differ from another.

Read the full publication here.

Available online 19 June 2021 in the Journal of Biological Chemistry (doi 10.1016/j.jbc.2021.100899):

Regulation of a pentameric ligand-gated ion channel by a semi-conserved cationic-lipid binding site

Akshay Sridhar*, Sarah CR Lummis*, Diletta Pasini, Aujan Mehregan, Marijke Brams, Kumiko Kambara, Daniel Bertrand, Erik Lindahl, Rebecca J Howard, Chris Ulens

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA, and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of non-canonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABA_A receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABA_A receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Read the full publication here.

*equal contributions

Available online 6 June 2021 in the Journal of Molecular Biology (doi 10.1016/j.jmb.2021.167095):

Structure and sequence-based computational approaches to allosteric signal transduction: application to electromechanical coupling in voltage-gated ion channels

Ahmad Elbahnsi, Lucie Delemotte

Allosteric signaling underlies the function of many biomolecules, including membrane proteins such as ion channels. Experimental methods have enabled specific quantitative insights into the coupling between the voltage sensing domain and the pore gate of voltage-gated ion channels, located tens of Angstrom apart from one another, as well as pinpointed specific residues and domains that participate in electromechanical signal transmission. Nevertheless, an overall atomic-level resolution picture is difficult to obtain from these methods alone. Today, thanks to the cryo-EM resolution revolution, we have access to high resolution structures of many different voltage-gated ion channels in various conformational states, putting a quantitative description of the processes at the basis of these changes within our close reach. Here, we review computational methods that build on structures to detect and characterize allosteric signaling and pathways. We then examine what has been learned so far about electromechanical coupling between VSD and pore domain using such methods. While no general theory of electromechanical coupling in voltage-gated ion channels integrating results from all these methods is available yet, we outline the types of insights that could be achieved in the near future using the methods that have not yet been put to use in this field of application.

Read the full publication here.

Available online 25 May 2021 in The Journal of Chemical Physics (doi 10.1063/5.0044352):

The accelerated weight histogram method for alchemical free energy calculations

Magnus Lundborg, Jack Lidmar, Berk Hess

The accelerated weight histogram method is an enhanced sampling technique used to explore free energy landscapes by applying an adaptive bias. The method is general and easy to extend. Herein, we show how it can be used to efficiently sample alchemical transformations, commonly used for, e.g., solvation and binding free energy calculations. We present calculations and convergence of the hydration free energy of testosterone, representing drug-like molecules. We also include methane and ethanol to validate the results. The protocol is easy to use, does not require a careful choice of parameters, and scales well to accessible resources, and the results converge at least as quickly as when using conventional methods. One benefit of the method is that it can easily be combined with other reaction coordinates, such as intermolecular distances.

Read the full publication here.

Available online 30 April 2021 in Trends in Biochemical Sciences (doi 10.1016/j.tibs.2021.04.005):

Molecular dynamics simulations of ion channels

Vincenzo Carnevale*, Lucie Delemotte*, Rebecca J Howard*

Propelled by enormous increases in computational power, molecular dynamics (MD) simulations were first reported in 1957 by B.J. Alder and T.E. Wainwright and since then have moved from this proof of concept to routinely investigating the dynamics of complex systems composed of up to tens of millions of atoms. MD simulations are based on the idea that the equations of motion of a multi-particle system can be solved numerically within an acceptable level of accuracy; the resulting trajectory is key for calculating occupancy probabilities of distinct conformational states (sampling). In advanced protocols, enhanced sampling ensures wide exploration of the configurational space. Further, a strength of MD simulations is that by starting from descriptions of the motions of all atoms, several data learning techniques can be used to conceptualize trajectories. In ion channel biophysics, these tools are used to study ion permeation, conformational cycling, drug binding, and lipid–channel interactions.

Read the full publication here.

*equal contributions

Published 22 April 2021 in the Journal of Chemical Physics (doi 10.1063/5.0040649):

Informing NMR experiments with molecular dynamics simulations to characterize the dominant activated state of the KcsA ion channel

Sergio Pérez-Conesa, Eric G. Keeler, Dongyu Zhang, Lucie Delemotte, Ann E. McDermott

As the first potassium channel with an x-ray structure determined, and given its homology to eukaryotic channels, the pH-gated prokaryotic channel KcsA has been extensively studied. Nevertheless, questions related, in particular, to the allosteric coupling between its gates remain open. The many currently available x-ray crystallography structures appear to correspond to various stages of activation and inactivation, offering insights into the molecular basis of these mechanisms. Since these studies have required mutations, complexation with antibodies, and substitution of detergents in place of lipids, examining the channel under more native conditions is desirable. Solid-state nuclear magnetic resonance (SSNMR) can be used to study the wild-type protein under activating conditions (low pH), at room temperature, and in bacteriomimetic liposomes. In this work, we sought to structurally assign the activated state present in SSNMR experiments. We used a combination of molecular dynamics (MD) simulations, chemical shift prediction algorithms, and Bayesian inference techniques to determine which of the most plausible x-ray structures resolved to date best represents the activated state captured in SSNMR. We first identified specific nuclei with simulated NMR chemical shifts that differed significantly when comparing partially open vs fully open ensembles from MD simulations. The simulated NMR chemical shifts for those specific nuclei were then compared to experimental ones, revealing that the simulation of the partially open state was in good agreement with the SSNMR data. Nuclei that discriminate effectively between partially and fully open states belong to residues spread over the sequence and provide a molecular level description of the conformational change.

Read the full publication here.

Published 9 April 2021 in the Journal of General Physiology (doi 10.1085/jgp.202012701):

Cannabidiol inhibits the skeletal muscle Nav1.4 by blocking its pore and by altering membrane elasticity

Mohammad-Reza Ghovanloo, Koushik Choudhury, Tagore S. Bandaru, Mohamed A. Fouda, Kaveh Rayani, Radda Rusinova, Tejas Phaterpekar, Karen Nelkenbrecher, Abeline R. Watkins, Damon Poburko, Jenifer Thewalt, Olaf S. Andersen, Lucie Delemotte, Samuel J. Goodchild, Peter C. Ruben 

Cannabidiol (CBD) is the primary nonpsychotropic phytocannabinoid found in Cannabis sativa, which has been proposed to be therapeutic against many conditions, including muscle spasms. Among its putative targets are voltage-gated sodium channels (Navs), which have been implicated in many conditions. We investigated the effects of CBD on Nav1.4, the skeletal muscle Nav subtype. We explored direct effects, involving physical block of the Nav pore, as well as indirect effects, involving modulation of membrane elasticity that contributes to Nav inhibition. MD simulations revealed CBD’s localization inside the membrane and effects on bilayer properties. Nuclear magnetic resonance (NMR) confirmed these results, showing CBD localizing below membrane headgroups. To determine the functional implications of these findings, we used a gramicidin-based fluorescence assay to show that CBD alters membrane elasticity or thickness, which could alter Nav function through bilayer-mediated regulation. Site-directed mutagenesis in the vicinity of the Nav1.4 pore revealed that removing the local anesthetic binding site with F1586A reduces the block of I_Na by CBD. Altering the fenestrations in the bilayer-spanning domain with Nav1.4-WWWW blocked CBD access from the membrane into the Nav1.4 pore (as judged by MD). The stabilization of inactivation, however, persisted in WWWW, which we ascribe to CBD-induced changes in membrane elasticity. To investigate the potential therapeutic value of CBD against Nav1.4 channelopathies, we used a pathogenic Nav1.4 variant, P1158S, which causes myotonia and periodic paralysis. CBD reduces excitability in both wild-type and the P1158S variant. Our in vitro and in silico results suggest that CBD may have therapeutic value against Nav1.4 hyperexcitability.

Read the full publication here.

Published 16 March 2021 in ACS Pharmacology & Translational Sciences (doi 10.1021/acsptsci.0c00215):

A blueprint for high affinity SARS-CoV-2 Mpro inhibitors from activity-based compound library screening guided by analysis of protein dynamics

Jonas Gossen, Simone Albani, Anton Hanke, Benjamin P Joseph, Cathrine Bergh, Maria Kuzikov, Elisa Costanzi, Candida Manelfi, Paola Storici, Philip Gribbon, Andrea R Beccari, Carmine Talarico, Francesca Spyrakis, Erik Lindahl, Andrea Zaliani, Paolo Carloni, Rebecca C Wade, Francesco Musiani, Daria B Kokh, and Giulia Rossetti

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein’s druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.

Read the full publication here.

Published 8 March 2021 in the Journal of General Physiology (doi 10.1085/jgp.202012676):

Resin-acid derivatives bind to multiple sites on the voltage-sensor domain of the Shaker potassium channel

Malin Silverå Ejneby, Arina Gromova, Nina E Ottosson, Stina Borg, Argel Estrada-Mondragón, Samira Yazdi, Panagiotis Apostolakis, Fredrik Elinder, Lucie Delemotte

Voltage-gated potassium (Kv) channels can be opened by negatively charged resin acids and their derivatives. These resin acids have been proposed to attract the positively charged voltage-sensor helix (S4) toward the extracellular side of the membrane by binding to a pocket located between the lipid-facing extracellular ends of the transmembrane segments S3 and S4. By contrast to this proposed mechanism, neutralization of the top gating charge of the Shaker KV channel increased resin-acid–induced opening, suggesting other mechanisms and sites of action. Here, we explore the binding of two resin-acid derivatives, Wu50 and Wu161, to the activated/open state of the Shaker KV channel by a combination of in silico docking, molecular dynamics simulations, and electrophysiology of mutated channels. We identified three potential resin-acid–binding sites around S4: (1) the S3/S4 site previously suggested, in which positively charged residues introduced at the top of S4 are critical to keep the compound bound, (2) a site in the cleft between S4 and the pore domain (S4/pore site), in which a tryptophan at the top of S6 and the top gating charge of S4 keeps the compound bound, and (3) a site located on the extracellular side of the voltage-sensor domain, in a cleft formed by S1–S4 (the top-VSD site). The multiple binding sites around S4 and the anticipated helical-screw motion of the helix during activation make the effect of resin-acid derivatives on channel function intricate. The propensity of a specific resin acid to activate and open a voltage-gated channel likely depends on its exact binding dynamics and the types of interactions it can form with the protein in a state-specific manner.

Read the full publication here.

Published 28 January 2021 in eLife (doi 10.7554/eLife.60715):

Identification of ligand-specific G-protein coupled receptor states and prediction of downstream efficacy via data-driven modeling

Oliver Fleetwood, Jens Carlsson, Lucie Delemotte

Ligand binding stabilizes different G protein-coupled receptor states via a complex allosteric process that is not completely understood. Here, we have derived free energy landscapes describing activation of the β2 adrenergic receptor bound to ligands with different efficacy profiles using enhanced sampling molecular dynamics (MD) simulations. These reveal shifts towards active-like states at the G protein binding site for receptors bound to partial and full agonists and that the ligands modulate the conformational ensemble of the receptor by tuning protein microswitches. We indeed find an excellent correlation between the conformation of the microswitches close to the ligand binding site and in the transmembrane region and experimentally reported cAMP signaling responses. Dimensionality reduction further reveals the similarity between the unique conformational states induced by different ligands and examining the output of classifiers highlights two distant hotspots governing agonism on transmembrane helices 5 and 7.

Read the full publication here.

Published 19 January 2021 in Scientific Reports (doi 10.1038/s41598-020-80352-8):

Direct detection of SARS-CoV-2 using non-commercial RT-LAMP reagents on heat-inactivated samples

Alisa Alekseenko, Donal Barrett, Yerma Pareja-Sanchez, Rebecca J Howard, Emilia Strandback, Henry Ampah-Korsah, Urška Rovšnik, Silvia Zuniga-Veliz, Alexander Klenov, Jayshna Malloo, Shenglong Ye, Xiyang Liu, Björn Reinius, Simon J. Elsässer, Tomas Nyman, Gustaf Sandh, Xiushan Yin, Vicent Pelechano

RT-LAMP detection of SARS-CoV-2 has been shown to be a valuable approach to scale up COVID-19 diagnostics and thus contribute to limiting the spread of the disease. Here we present the optimization of highly cost-effective in-house produced enzymes, and we benchmark their performance against commercial alternatives. We explore the compatibility between multiple DNA polymerases with high strand-displacement activity and thermostable reverse transcriptases required for RT-LAMP. We optimize reaction conditions and demonstrate their applicability using both synthetic RNA and clinical patient samples. Finally, we validate the optimized RT-LAMP assay for the detection of SARS-CoV-2 in unextracted heat-inactivated nasopharyngeal samples from 184 patients. We anticipate that optimized and affordable reagents for RT-LAMP will facilitate the expansion of SARS-CoV-2 testing globally, especially in sites and settings where the need for large scale testing cannot be met by commercial alternatives.

Read the full publication here.

Published 11 December 2020 in Science Advances (doi 10.1126/sciadv.abd6798):

Calmodulin acts as a state-dependent switch to control a cardiac potassium channel opening

Po Wei Kang,* Annie M Westerlund,* Jingyi Shi, Kelli McFarland White, Alex K Dou, Amy H Cui, Jonathan R Silva, Lucie Delemotte, Jianmin Cui
*contributed equally to this work

Calmodulin (CaM) and phosphatidylinositol 4,5-bisphosphate (PIP₂) are potent regulators of the voltage-gated potassium channel KCNQ1 (K_V7.1), which conducts the cardiac I_Ks current. Although cryo–electron microscopy structures revealed intricate interactions between the KCNQ1 voltage-sensing domain (VSD), CaM, and PIP₂, the functional consequences of these interactions remain unknown. Here, we show that CaM-VSD interactions act as a state-dependent switch to control KCNQ1 pore opening. Combined electrophysiology and molecular dynamics network analysis suggest that VSD transition into the fully activated state allows PIP₂ to compete with CaM for binding to VSD. This leads to conformational changes that alter VSD-pore coupling to stabilize open states. We identify a motif in the KCNQ1 cytosolic domain, which works downstream of CaM-VSD interactions to facilitate the conformational change. Our findings suggest a gating mechanism that integrates PIP₂and CaM in KCNQ1 voltage-dependent activation, yielding insights into how KCNQ1 gains the phenotypes critical for its physiological function.

Read the full publication here.

Published 22 October 2020 in Journal of Investigative Dermatology (doi 10.1016/j.jid.2020.07.040):

Molecular reorganization during formation of the human skin barrier studied in situ

Ali Narangifard, Christian L Wennberg, Lianne den Hollander, Ichiro Iwai, HongMei Han, Magnus Lundborg, Sergej Masich, Erik Lindahl, Bertil Daneholt, Lars Norlén

In vertebrates, skin upholds homeostasis by preventing body water loss. The skin’s permeability barrier is located intercellularly in stratum corneum and consists of stacked lipid lamellae composed of ceramides, cholesterol and free fatty acids. We have combined cryo-EM with molecular dynamics modelling and EM-simulation in our analysis of the lamellae’s formation, a maturation process beginning in stratum granulosum and ending in stratum corneum. Previously, we have revealed the lipid lamellae’s initial- and end-stage molecular organizations. Here, we reveal two cryo-EM patterns representing intermediate stages in the lamellae’s maturation process: a single-band pattern with 2.0-2.5 nm periodicity and a two-band pattern with 5.5-6.0 nm periodicity, that may be derived from lamellar lipid structures with 4.0-5.0 nm and 5.5-6.0 nm periodicity, respectively. Based on the analysis of the data now available on the four maturation stages identified, we can present a tentative molecular model for the complete skin barrier formation process.

Read the full publication here.