Allosteric potentiation of a ligand-gated ion channel is mediated by access to a deep membrane-facing cavity

Discovery of lipid binding sites in a ligand-gated ion channel by integrating simulations and cryo-EM

Ligand-gated ion channels transduce electrochemical signals in neurons and other excitable cells. Aside from canonical ligands, phospholipids are thought to bind specifically to the transmembrane domain of several ion channels. However, structural details of such lipid contacts remain elusive, partly due to limited resolution of these regions in experimental structures. Here, we discovered multiple lipid interactions in the channel GLIC by integrating cryo-electron microscopy and large-scale molecular simulations. We identified 25 bound lipids in the GLIC closed state, a conformation where none, to our knowledge, were previously known. Three lipids were associated with each subunit in the inner leaflet, including a buried interaction disrupted in mutant simulations. In the outer leaflet, two intrasubunit sites were evident in both closed and open states, while a putative intersubunit site was preferred in open-state simulations. This work offers molecular details of GLIC-lipid contacts particularly in the ill-characterized closed state, testable hypotheses for state-dependent binding, and a multidisciplinary strategy for modeling protein-lipid interactions.

Mechanisms of Altered Immune Response in Skin Melanoma

Melanoma, a deadly form of skin cancer, poses significant challenges to the host immune system, allowing tumor cells to evade immune surveillance and persist. This complex interplay between melanoma and the immune system involves a multitude of mechanisms that impair immune recognition and promote tumor progression. This review summarizes the intricate strategies employed by melanoma cells to evade the immune response, including defective immune recognition, immune checkpoint activation, and the role of regulatory T-cells, myeloid-derived suppressor cells, and exosomes in suppressing anti-tumor immunity. Additionally, we discuss potential therapeutic targets aimed at reversing immune evasion in melanoma, highlighting the importance of understanding these mechanisms for developing more effective immunotherapies. Improved insights into the interactions between melanoma and the immune system will aid in the development of novel treatment strategies to enhance anti-tumor immune responses and improve patient outcomes.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

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

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.

Mutational analysis to explore long-range allosteric couplings involved in a pentameric channel receptor pre-activation and activation

Pentameric ligand-gated ion channels (pLGICs) mediate chemical signaling through a succession of allosteric transitions that are yet not completely understood as intermediate states remain poorly characterized by structural approaches. In a previous study on the prototypic bacterial proton-gated channel GLIC, we generated several fluorescent sensors of the protein conformation that report a fast transition to a pre-active state, which precedes the slower process of activation with pore opening. Here, we explored the phenotype of a series of allosteric mutations, using simultaneous steady-state fluorescence and electrophysiological measurements over a broad pH range. Our data, fitted to a three-state Monod-Wyman-Changeux model, show that mutations at the subunit interface in the extracellular domain (ECD) principally alter pre-activation, while mutations in the lower ECD and in the transmembrane domain principally alter activation. We also show that propofol alters both transitions. Data are discussed in the framework of transition pathways generated by normal mode analysis (iModFit). It further supports that pre-activation involves major quaternary compaction of the ECD, and suggests that activation involves principally a reorganization of a ‘central gating region’ involving a contraction of the ECD β-sandwich and the tilt of the channel lining M2 helix.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Regulation and drug modulation of a voltage-gated sodium channel: Pivotal role of the S4-S5 linker in activation and slow inactivation

Voltage-gated sodium channels initiate electric signals in cell communications. The S4–S5 linker between the voltage-sensing and pore modules transmits depolarization signals to trigger channel activation. The mechanisms of this action, however, remain elusive. By combining biophysical and computational approaches, we identify a critical residue, T140, in the S4–S5 linker of the bacterial sodium channel NaChBac, which plays a pivotal role in channel activation and drug modulation of slow inactivation. Specifically, we discovered conformation-dependent drug binding at this site and unveiled a toggling mode of action by T140, which switches interaction partners with different S6 residues to regulate channel activation and slow inactivation. These observations suggest the possibility of conformation-specific drugs targeting the gating machinery of voltage-gated ion channels.

Voltage-gated sodium (Na_V) channels control excitable cell functions. While structural investigations have revealed conformation details of different functional states, the mechanisms of both activation and slow inactivation remain unclear. Here, we identify residue T140 in the S4–S5 linker of the bacterial voltage-gated sodium channel NaChBac as critical for channel activation and drug effects on inactivation. Mutations at T140 either attenuate activation or render the channel nonfunctional. Propofol, a clinical anesthetic known to inhibit NaChBac by promoting slow inactivation, binds to a pocket between the S4–S5 linker and S6 helix in a conformation-dependent manner. Using ¹⁹F-NMR to quantify site-specific binding by saturation transfer differences (STDs), we found strong STDs in inactivated, but not activated, NaChBac. Molecular dynamics simulations show a highly dynamic pocket in the activated conformation, limiting STD buildup. In contrast, drug binding to this pocket promotes and stabilizes the inactivated states. Our results provide direct experimental evidence showing distinctly different associations between the S4–S5 linker and S6 helix in activated and inactivated states. Specifically, an exchange occurs between interaction partners T140 and N234 of the same subunit in activation, and T140 and N225 of the domain-swapped subunit in slow inactivation. The drug action on slow inactivation of prokaryotic Na_V channels seems to have a mechanism similar to the recently proposed “door-wedge” action of the isoleucine-phenylalanine-methionine (IFM) motif on the fast inactivation of eukaryotic Na_V channels. Elucidating this gating mechanism points to a possible direction for conformation-dependent drug development.

A Single Mutation in the Outer Lipid-Facing Helix of a Pentameric Ligand-Gated Ion Channel Affects Channel Function Through a Radially-Propagating Mechanism

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic transmission and are crucial drug targets. Their gating mechanism is triggered by ligand binding in the extracellular domain that culminates in the opening of a hydrophobic gate in the transmembrane domain. This domain is made of four α-helices (M1 to M4). Recently the outer lipid-facing helix (M4) has been shown to be key to receptor function, however its role in channel opening is still poorly understood. It could act through its neighboring helices (M1/M3), or via the M4 tip interacting with the pivotal Cys-loop in the extracellular domain. Mutation of a single M4 tyrosine (Y441) to alanine renders one pLGIC-the 5-HT3A receptor-unable to function despite robust ligand binding. Using Y441A as a proxy for M4 function, we here predict likely paths of Y441 action using molecular dynamics, and test these predictions with functional assays of mutant receptors in HEK cells and Xenopus oocytes using fluorescent membrane potential sensitive dye and two-electrode voltage clamp respectively. We show that Y441 does not act via the M4 tip or Cys-loop, but instead connects radially through M1 to a residue near the ion channel hydrophobic gate on the pore-lining helix M2. This demonstrates the active role of the M4 helix in channel opening.

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

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.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Mechanistic insights into volatile anesthetic modulation of K2P channels

K2P potassium channels are known to be modulated by volatile anesthetic (VA) drugs and play important roles in clinically relevant effects that accompany general anesthesia. Here, we utilize a photoaffinity analog of the VA isoflurane to identify a VA-binding site in the TREK1 K2P channel. The functional importance of the identified site was validated by mutagenesis and biochemical modification. Molecular dynamics simulations of TREK1 in the presence of VA found multiple neighboring residues on TREK1 TM2, TM3, and TM4 that contribute to anesthetic binding. The identified VA-binding region contains residues that play roles in the mechanisms by which heat, mechanical stretch, and pharmacological modulators alter TREK1 channel activity and overlaps with positions found to modulate TASK K2P channel VA sensitivity. Our findings define molecular contacts that mediate VA binding to TREK1 channels and suggest a mechanistic basis to explain how K2P channels are modulated by VAs.

The Biology of General Anesthesia from Paramecium to Primate

General anesthesia serves a critically important function in the clinical care of human patients. However, the anesthetized state has foundational implications for biology because anesthetic drugs are effective in organisms ranging from paramecia, to plants, to primates. Although unconsciousness is typically considered the cardinal feature of general anesthesia, this endpoint is only strictly applicable to a select subset of organisms that are susceptible to being anesthetized. We review the behavioral endpoints of general anesthetics across species and propose the isolation of an organism from its environment - both in terms of the afferent arm of sensation and the efferent arm of action - as a generalizable definition. We also consider the various targets and putative mechanisms of general anesthetics across biology and identify key substrates that are conserved, including cytoskeletal elements, ion channels, mitochondria, and functionally coupled electrical or neural activity. We conclude with a unifying framework related to network function and suggest that general anesthetics - from single cells to complex brains - create inefficiency and enhance modularity, leading to the dissociation of functions both within an organism and between the organism and its surroundings. Collectively, we demonstrate that general anesthesia is not restricted to the domain of modern medicine but has broad biological relevance with wide-ranging implications for a diverse array of species.

Menthol Binding to the Human α4β2 Nicotinic Acetylcholine Receptor Facilitated by Its Strong Partitioning in the Membrane

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.