Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Serotonin Signaling through Lipid Membranes

Serotonin (5-HT) is a vital modulatory neurotransmitter responsible for regulating most behaviors in the brain. An inefficient 5-HT synaptic function is often linked to various mental disorders. Primarily, membrane proteins controlling the expression and activity of 5-HT synthesis, storage, release, receptor activation, and inactivation are critical to 5-HT signaling in synaptic and extra-synaptic sites. Moreover, these signals represent information transmission across membranes. Although the lipid membrane environment is often viewed as fairly stable, emerging research suggests significant functional lipid-protein interactions with many synaptic 5-HT proteins. These protein-lipid interactions extend to almost all the primary lipid classes that form the plasma membrane. Collectively, these lipid classes and lipid-protein interactions affect 5-HT synaptic efficacy at the synapse. The highly dynamic lipid composition of synaptic membranes suggests that these lipids and their interactions with proteins may contribute to the plasticity of the 5-HT synapse. Therefore, this broader protein-lipid model of the 5-HT synapse necessitates a reconsideration of 5-HT's role in various associated mental disorders.

Jankowski W, Várnai C, Irwin BWJ, Payne MC, Chau PL. Efficiently Differentiating Agonists and Competitive Antagonists for Weak Allosteric Protein-Ligand Interactions with Linear Response Theory

What makes an agonist and a competitive antagonist? In this work, we aim to answer this question by performing parallel tempering Monte Carlo simulations on the serotonin type 3A (5-HT3A) receptor. We use linear response theory to predict conformational changes in the 5-HT3A receptor active site after weak perturbations are applied to its allosteric binding sites. A covariance tensor is built from conformational sampling of its apo state, and a harmonic approximation allows us to substitute the calculation of ligand-induced forces with the binding site's displacement vector. Remarkably, our study demonstrates the feasibility of effectively discerning between agonists and competitive antagonists for multiple ligands, requiring computationally expensive calculations only once per protein.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Recent Insight into Lipid Binding and Lipid Modulation of Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels (pLGICs) play a leading role in synaptic communication, are implicated in a variety of neurological processes, and are important targets for the treatment of neurological and neuromuscular disorders. Endogenous lipids and lipophilic compounds are potent modulators of pLGIC function and may help shape synaptic communication. Increasing structural and biophysical data reveal sites for lipid binding to pLGICs. Here, we update our evolving understanding of pLGIC–lipid interactions highlighting newly identified modes of lipid binding along with the mechanistic understanding derived from the new structural data.

Mutations of the nACh Receptor M4 Helix Reveal Different Phenotypes in Different Expression Systems: Could Lipids be Responsible?

The role of the outermost helix (M4) in the pentameric ligand-gated ion channel (pLGIC) family is currently not fully understood. It is known that M4 is important for receptor assembly, possibly via interactions with neighboring M1 and M3 helices. M4 can also transmit information on the lipid content of the membrane to the gating mechanism, and it may form a link to the extracellular domain via the Cys-loop. Our previous study examining the α4β2 nACh receptor M4 helix using HEK cells indicated M4 here is more sensitive to change than those of other pLGIC. Many of these other studies, however, were performed in Xenopus oocytes. Here we examine the nine previously identified nonfunctional α4β2 nACh receptor M4 mutant receptors using this system. The data reveal that seven of these mutant receptors do function when expressed in oocytes, with only 2, the conserved Asp at the intracellular end of M4 and a Phe in the center, having a similar phenotype (nonfunctional) in both HEK cells and oocytes. The oocyte data are more consistent with studies in other pLGIC and demonstrate the importance of the expression system used. Of the many differences between these two expression systems, we suggest that the different lipid content of the plasma membrane is a possible candidate for explaining these discrepancies.

Conformational decoupling in acid-sensing ion channels uncovers mechanism and stoichiometry of PcTx1-mediated inhibition

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels involved in fast synaptic transmission. Pharmacological inhibition of ASIC1a reduces neurotoxicity and stroke infarct volumes, with the cysteine knot toxin psalmotoxin-1 (PcTx1) being one of the most potent and selective inhibitors. PcTx1 binds at the subunit interface in the extracellular domain (ECD), but the mechanism and conformational consequences of the interaction, as well as the number of toxin molecules required for inhibition, remain unknown. Here, we use voltage-clamp fluorometry and subunit concatenation to decipher the mechanism and stoichiometry of PcTx1 inhibition of ASIC1a. Besides the known inhibitory binding mode, we propose PcTx1 to have at least two additional binding modes that are decoupled from the pore. One of these modes induces a long-lived ECD conformation that reduces the activity of an endogenous neuropeptide. This long-lived conformational state is proton-dependent and can be destabilized by a mutation that decreases PcTx1 sensitivity. Lastly, the use of concatemeric channel constructs reveals that disruption of a single PcTx1 binding site is sufficient to destabilize the toxin-induced conformation, while functional inhibition is not impaired until two or more binding sites are mutated. Together, our work provides insight into the mechanism of PcTx1 inhibition of ASICs and uncovers a prolonged conformational change with possible pharmacological implications.

Tapping into 5-HT3 Receptors to Modify Metabolic and Immune Responses

5-hydroxytryptamine type 3 (5-HT₃) receptors are ligand gated ion channels, which clearly distinguish their mode of action from the other G-protein coupled 5-HT or serotonin receptors. 5-HT₃ receptors are well established targets for emesis and gastrointestinal mobility and are used as adjunct targets in treating schizophrenia. However, the distribution of these receptors is wider than the nervous system and there is potential that these additional sites can be targeted to modulate inflammatory and/or metabolic conditions. Recent progress in structural biology and pharmacology of 5-HT₃ receptors have provided profound insights into mechanisms of their action. These advances, combined with insights into clinical relevance of mutations in genes encoding 5-HT₃ subunits and increasing understanding of their implications in patient's predisposition to diseases and response to the treatment, open new avenues for personalized precision medicine. In this review, we recap on the current status of 5-HT₃ receptor-based therapies using a biochemical and physiological perspective. We assess the potential for targeting 5-HT₃ receptors in conditions involving metabolic or inflammatory disorders based on recent findings, underscoring the challenges and limitations of this approach.

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.

Considerations of Recent All-Atom Lipid Force Field Development

Molecular simulations of biological molecules require an accurate description of molecular interactions through a force field (FF). The focus of this Perspective is on all-atom lipid FFs. Recent additions to the CHARMM36 lipid FF continue to expand a researcher's ability to probe membrane structure and function with a wide variety of biologically important lipids. Currently, there is an effort to reduce the assumptions in all-atom lipid FFs. The inclusion of long-range dispersion interaction through particle-mesh Ewald is allowing for more accurate descriptions of lipid bilayer and monolayer properties without additional computational cost. Soon, simulations with lipid FFs will no longer depend on short-range cutoffs and will accurately represent long-range dispersion. This requires efficient FF parametrization with an automated approach due to FF complexity. In addition, polarizable FFs for lipids will be important for the next generation of simulations that accurately represent how molecule interactions respond to a varied environment.

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.

Asymmetric opening of the homopentameric 5-HT3A serotonin receptor in lipid bilayers

Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop receptor family are key players in fast signal transduction throughout the nervous system. They have been shown to be modulated by the lipid environment, however the underlying mechanism is not well understood. We report three structures of the Cys-loop 5-HT_3A serotonin receptor (5HT₃R) reconstituted into saposin-based lipid bilayer discs: a symmetric and an asymmetric apo state, and an asymmetric agonist-bound state. In comparison to previously published 5HT₃R conformations in detergent, the lipid bilayer stabilises the receptor in a more tightly packed, ‘coupled’ state, involving a cluster of highly conserved residues. In consequence, the agonist-bound receptor conformation adopts a wide-open pore capable of conducting sodium ions in unbiased molecular dynamics (MD) simulations. Taken together, we provide a structural basis for the modulation of 5HT₃R by the membrane environment, and a model for asymmetric activation of the receptor.

Pentameric ligand-gated ion channels (pLGICs) are key players in neurotransmission and have been shown to be modulated by the lipid environment, however the underlying mechanism is not well understood. Here, the authors report structures of the pLGIC 5-HT_3A serotonin receptor reconstituted into lipid bilayer discs and reveal lipid–protein interactions as well as asymmetric activation of the homopentameric receptor.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Tracking Membrane Protein Dynamics in Real Time

Abstract

Membrane proteins govern critical cellular processes and are central to human health and associated disease. Understanding of membrane protein function is obscured by the vast ranges of structural dynamics—both in the spatial and time regime—displayed in the protein and surrounding membrane. The membrane lipids have emerged as allosteric modulators of membrane protein function, which further adds to the complexity. In this review, we discuss several examples of membrane dependency. A particular focus is on how molecular dynamics (MD) simulation have aided to map membrane protein dynamics and how enhanced sampling methods can enable observing the otherwise inaccessible biological time scale. Also, time-resolved X-ray scattering in solution is highlighted as a powerful tool to track membrane protein dynamics, in particular when combined with MD simulation to identify transient intermediate states. Finally, we discuss future directions of how to further develop this promising approach to determine structural dynamics of both the protein and the surrounding lipids.

Graphic Abstract

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.

Spontaneous Membrane Nanodomain Formation in the Absence or Presence of the Neurotransmitter Serotonin

Detailed knowledge on the formation of biomembrane domains, their structure, composition, and physical characteristics is scarce. Despite its frequently discussed importance in signaling, e.g., in obtaining localized non-homogeneous receptor compositions in the plasma membrane, the nanometer size as well as the dynamic and transient nature of domains impede their experimental characterization. In turn, atomistic molecular dynamics (MD) simulations combine both, high spatial and high temporal resolution. Here, using microsecond atomistic MD simulations, we characterize the spontaneous and unbiased formation of nano-domains in a plasma membrane model containing phosphatidylcholine (POPC), palmitoyl-sphingomyelin (PSM), and cholesterol (Chol) in the presence or absence of the neurotransmitter serotonin at different temperatures. In the ternary mixture, highly ordered and highly disordered domains of similar composition coexist at 303 K. The distinction of domains by lipid acyl chain order gets lost at lower temperatures of 298 and 294 K, suggesting a phase transition at ambient temperature. By comparison of domain ordering and composition, we demonstrate how the domain-specific binding of the neurotransmitter serotonin results in a modified domain lipid composition and a substantial downward shift of the phase transition temperature. Our simulations thus suggest a novel mode of action of neurotransmitters possibly of importance in neuronal signal transmission.

High-resolution structures of multiple 5-HT3AR-setron complexes reveal a novel mechanism of competitive inhibition

Serotonin receptors (5-HT3AR) play a crucial role in regulating gut movement, and are the principal target of setrons, a class of high-affinity competitive antagonists, used in the management of nausea and vomiting associated with radiation and chemotherapies. Structural insights into setron-binding poses and their inhibitory mechanisms are just beginning to emerge. Here, we present high-resolution cryo-EM structures of full-length 5-HT3AR in complex with palonosetron, ondansetron, and alosetron. Molecular dynamic simulations of these structures embedded in a fully-hydrated lipid environment assessed the stability of ligand-binding poses and drug-target interactions over time. Together with simulation results of apo- and serotonin-bound 5-HT3AR, the study reveals a distinct interaction fingerprint between the various setrons and binding-pocket residues that may underlie their diverse affinities. In addition, varying degrees of conformational change in the setron-5-HT3AR structures, throughout the channel and particularly along the channel activation pathway, suggests a novel mechanism of competitive inhibition.

Microsecond dynamics in proteins by two-dimensional ESR: Predictions

Two-dimensional electron-electron double resonance (2D-ELDOR) provides extensive insight into molecular motions. Recent developments permitting experiments at higher frequencies (95 GHz) provide molecular orientational resolution, enabling a clearer description of the nature of the motions. In this work, simulations are provided for the example of domain motions within proteins that are themselves slowly tumbling in solution. These show the nature of the exchange cross-peaks that are predicted to develop in real time from such domain motions. However, we find that the existing theoretical methods for computing 2D-ELDOR experiments over a wide motional range begin to fail seriously when applied to very slow motions characteristic of proteins in solution. One reason is the failure to obtain accurate eigenvectors and eigenvalues of the complex symmetric stochastic Liouville matrices describing the experiment when computed by the efficient Lanczos algorithm in the range of very slow motion. Another, perhaps more serious, issue is that these matrices are "non-normal," such that for the very slow motional range even rigorous diagonalization algorithms do not yield the correct eigenvalues and eigenvectors. We have employed algorithms that overcome both these issues and lead to valid 2D-ELDOR predictions even for motions approaching the rigid limit. They are utilized to describe the development of cross-peaks in 2D-ELDOR at 95 GHz for a particular case of domain motion.

Assessing the extent of the structural and dynamic modulation of membrane lipids due to pore forming toxins: insights from molecular dynamics simulations

Infections caused by many virulent bacterial strains are triggered by the release of pore forming toxins (PFTs), which form oligomeric transmembrane pore complexes on the target plasma membrane. The spatial extent of the perturbation to the surrounding lipids during pore formation is relatively unexplored. Using all-atom molecular dynamics simulations, we investigate the changes in the structure and dynamics of lipids in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer in the presence of contrasting PFTs. Cytolysin A (ClyA), an α toxin with its inserted wedge shaped bundle of inserted α helices, induces significant asymmetry across the membrane leaflets in comparison with α hemolysin (AHL), a β toxin. Despite the differences in hydrophobic mismatch and uniquely different topologies of the two oligomers, perturbations to lipid order as reflected in the tilt angle and order parameters and membrane thinning are short ranged, lying within ∼2.5 nm from the periphery of either pore complex, and commensurate with distances typically associated with van der Waals forces. In contrast, the spatial extent of perturbations to the lipid dynamics extends outward to at least 4 nm for both proteins, and the continuous survival probabilities reveal the presence of a tightly bound shell of lipids in this region. Displacement probability distributions show long tails and the distinctly non-Gaussian features reflect the induced dynamic heterogeneity. A detailed profiling of the protein-lipid contacts with tyrosine, tryptophan, lysine and arginine residues shows increased non-polar contacts in the cytoplasmic leaflet for both PFTs, with a higher number of atomic contacts in the case of AHL in the extracellular leaflet due to the mushroom-like topology of the pore complex. The short ranged nature of the perturbations observed in this simple one component membrane suggests inherent plasticity of membrane lipids enabling the recovery of the structure and membrane fluidity even in the presence of these large oligomeric transmembrane protein assemblies. This observation has implications in membrane repair processes such as budding or vesicle fusion events used to mitigate PFT virulence, where the underlying lipid dynamics and fluidity in the vicinity of the pore complex are expected to play an important role.

Functional movements of the GABA type A receptor

We have performed a Monte Carlo simulation of the GABA type A receptor. We have analysed the configurations and developed a correlation tensor method to predict receptor gating.