A Gating Mechanism of the Serotonin 5-HT3 Receptor

Mechanism of hydrophobic gating in the acetylcholine receptor channel pore

Neuromuscular acetylcholine receptors (AChRs) are hetero-pentameric, ligand-gated ion channels. The binding of the neurotransmitter acetylcholine (ACh) to two target sites promotes a global conformational change of the receptor that opens the channel and allows ion conduction through the channel pore. Here, by measuring free-energy changes from single-channel current recordings and using molecular dynamics simulations, we elucidate how a constricted hydrophobic region acts as a "gate" to regulate the channel opening in the pore of AChRs. Mutations of gate residues, including those implicated in congenital myasthenia syndrome, lower the permeation barrier of the channel substantially and increase the unliganded gating equilibrium constant (constitutive channel openings). Correlations between hydrophobicity and the observed free-energy changes, supported by calculations of water densities in the wild-type versus mutant channel pores, provide evidence for hydrophobic wetting-dewetting transition at the gate. The analysis of a coupled interaction network provides insight into the molecular mechanism of closed- versus open-state conformational changes at the gate. Studies of the transition state by "phi"(φ)-value analysis indicate that agonist binding serves to stabilize both the transition and the open state. Intersubunit interaction energy measurements and molecular dynamics simulations suggest that channel opening involves tilting of the pore-lining M2 helices, asymmetric outward rotation of amino acid side chains, and wetting transition of the gate region that lowers the barrier to ion permeation and stabilizes the channel open conformation. Our work provides new insight into the hydrophobic gate opening and shows why the gate mutations result in constitutive AChR channel activity.

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.

Phylogenetic analyses of 5-hydroxytryptamine 3 (5-HT3) receptors in Metazoa

The 5-hydroxytrptamine 3 (5-HT3) receptor is a member of the 'Cys-loop' family and the only pentameric ligand gated ion channel among the serotonin receptors. 5-HT3 receptors play an important role in controlling growth, development, and behaviour in animals. Several 5-HT3 receptor antagonists are used to treat diseases (e.g., irritable bowel syndrome, nausea and emesis). Humans express five different subunits (A-E) enabling a variety of heteromeric receptors to form but all contain 5HT3A subunits. However, the information available about the 5-HT3 receptor subunit occurrence among the metazoan lineages is minimal. In the present article we searched for 5-HT3 receptor subunit homologs from different phyla in Metazoa. We identified more than 1000 5-HT3 receptor subunits in Metazoa in different phyla and undertook simultaneous phylogenetic analysis of 526 5HT3A, 358 5HT3B, 239 5HT3C, 70 5HT3D, and 173 5HT3E sequences. 5-HT3 receptor subunits were present in species belonging to 11 phyla: Annelida, Arthropoda, Chordata, Cnidaria, Echinodermata, Mollusca, Nematoda, Orthonectida, Platyhelminthes, Rotifera and Tardigrada. All subunits were most often identified in Chordata phylum which was strongly represented in searches. Using multiple sequence alignment, we investigated variations in the ligand binding region of the 5HT3A subunit protein sequences in the metazoan lineage. Several critical amino acid residues important for ligand binding (common structural features) are commonly present in species from Nematoda and Platyhelminth gut parasites through to Chordata. Collectively, this better understanding of the 5-HT3 receptor evolutionary patterns raises possibilities of future pharmacological challenges facing Metazoa including effects on parasitic and other species in ecosystems that contain 5-HT3 receptor ligands.

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.

When is a hydrophobic gate not a hydrophobic gate?

The flux of ions through a channel is most commonly regulated by changes that result in steric occlusion of its pore. However, ion permeation can also be prevented by formation of a desolvation barrier created by hydrophobic residues that line the pore. As a result of relatively minor structural changes, confined hydrophobic regions in channels may undergo transitions between wet and dry states to gate the pore closed without physical constriction of the permeation pathway. This concept is referred to as hydrophobic gating, and many examples of this process have been demonstrated. However, the term is also now being used in a much broader context that often deviates from its original meaning. In this Viewpoint, we explore the formal definition of a hydrophobic gate, discuss examples of this process compared with other gating mechanisms that simply exploit hydrophobic residues and/or lipids in steric closure of the pore, and describe the best practice for identification of a hydrophobic gate.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

Using Metadynamics To Explore the Free Energy of Dewetting in Biologically Relevant Nanopores

Water confined within hydrophobic spaces can undergo cooperative dewetting transitions due to slight changes in water density and pressure that push water toward the vapor phase. Many transmembrane protein ion channels contain nanoscale hydrophobic pores that could undergo dewetting transitions, sometimes blocking the flow of ions without physical blockages. Standard molecular dynamics simulations have been extensively applied to study the behavior of water in nanoscale pores, but the large free energy barriers of dewetting often prevent direct sampling of both wet and dry states and quantitative studies of the hydration thermodynamics of biologically relevant pores. Here, we describe a metadynamics protocol that uses the number of waters within the pore as the collective variable to drive many reversible transitions between relevant hydration states and calculate well-converged free energy profiles of pore hydration. By creating model nanopore systems and changing their radius and morphology and including various cosolvents, we quantify how these pore properties and cosolvents affect the dewetting transition. The results reveal that the dewetting free energy of nanoscale pores is determined by two key thermodynamic parameters, namely, the effective surface tension coefficients of water-air and water-pore interfaces. Importantly, while the effect of salt can be fully captured in the water activity dependence, amphipathic cosolvents such as alcohols modify both dry and wet states of the pore and dramatically shift the wet-dry equilibrium. The metadynamics approach could be applied to studies of dewetting transitions within nanoscale pores of proteins and provide new insights into why different pore properties evolved in biological systems.

PharmacoGenetic targeting of a C. elegans essential neuron provides an in vivo screening for novel modulators of nematode ion channel function

Chemical or drug treatments are successfully used to treat parasitic nematode infections that impact human, animal and plant health. Many of these exert their effects through modifying neural function underpinning behaviours essential for parasite viability. Selectivity against the parasite may be achieved through distinct pharmacological properties of the parasite nervous system, as exemplified by the success of the ivermectin which target a glutamate-gated chloride channel found only in invertebrates. Despite the success of the ivermectins, emerging resistance and concerns around eco-toxicity are driving the search for new nematocidal chemicals or drugs. Here, we describe the potential of a 5-HT-gated chloride channel MOD-1, which is involved in vital parasite behaviours with constrained distribution in the invertebrate phyla. This ion channel has potential pharmacophores that could be targeted by new nematocidal chemicals and drugs. We have developed a microtiter based bioassay for MOD-1 pharmacology based on its ectopic expression in the Caenorhabditis elegans essential neuron M4. We have termed this technology 'PhaGeM4' for 'Pharmacogenetic targeting of M4 neuron'. Exposure of transgenic worms harbouring ectopically expressed MOD-1 to 5-HT results in developmental arrest. By additional expression of a fluorescence marker in body wall muscle to monitor growth we demonstrate that this assay is suitable for the identification of receptor agonists and antagonists. Indeed, the developmental progression is a robustly quantifiable bioassay that resolves MOD-1 activation by quipazine, 5-carboxyamidotryptamine and fluoxetine and highlight methiothepin as a potent antagonist. This assay has the intrinsic ability to highlight compounds with optimal bioavailability and furthermore to filter out off-target effects. It can be extended to the investigation of other classes of membrane receptors and modulators of neuronal excitation. This approach based on heterologous modulation of the essential M4 neuron function offers a route to discover new effective and selective anthelmintics potentially less confounded by disruptive environmental impact.

Hydroxy Pentacyclic Triterpene Acid, Kaempferol, Inhibits the Human 5-Hydroxytryptamine Type 3A Receptor Activity

Monoamine serotonin is a major neurotransmitter that acts on a wide range of central nervous system and peripheral nervous system functions and is known to have a role in various processes. Recently, it has been found that 5-HT is involved in cognitive and memory functions through interaction with cholinergic pathways. The natural flavonoid kaempferol (KAE) extracted from Cudrania tricuspidata is a secondary metabolite of the plant. Recently studies have confirmed that KAE possesses a neuroprotective effect because of its strong antioxidant activity. It has been confirmed that KAE is involved in the serotonergic pathway through an in vivo test. However, these results need to be confirmed at the molecular level, because the exact mechanism that is involved in such effects of KAE has not yet been elucidated. Therefore, the objective of this study is to confirm the interaction of KAE with 5-HT_3A through electrophysiological studies at the molecular level using KAE extracted from Cudrania tricuspidata. This study confirmed the interaction between 5-HT_3A and KAE at the molecular level. KAE inhibited 5-HT_3A receptors in a concentration-dependent and voltage-independent manner. Site-directed mutagenesis and molecular-docking studies confirmed that the binding sites D177 and F199 are the major binding sites of human 5-HT_3A receptors of KAE.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.

Structure, Function and Physiology of 5-Hydroxytryptamine Receptors Subtype 3

5-hydroxytryptamine receptor subtype 3 (5-HT₃R) is a pentameric ligand-gated ion channel (pLGIC) involved in neuronal signaling. It is best known for its prominent role in gut-CNS signaling though there is growing interest in its other functions, particularly in modulating non-serotonergic synaptic activity. Recent advances in structural biology have provided mechanistic understanding of 5-HT₃R function and present new opportunities for the field. This chapter gives a broad overview of 5-HT₃R from a physiological and structural perspective and then discusses the specific details of ion permeation, ligand binding and allosteric coupling between these two events. Biochemical evidence is summarized and placed within a physiological context. This perspective underscores the progress that has been made as well as outstanding challenges and opportunities for future 5-HT₃R research.

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.

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.

International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-hydroxytryptamine; Pharmacology and Function

5-HT receptors expressed throughout the human body are targets for established therapeutics and various drugs in development. Their diversity of structure and function reflects the important role 5-HT receptors play in physiologic and pathophysiological processes. The present review offers a framework for the official receptor nomenclature and a detailed understanding of each of the 14 5-HT receptor subtypes, their roles in the systems of the body, and, where appropriate, the (potential) utility of therapeutics targeting these receptors. SIGNIFICANCE STATEMENT: This review provides a comprehensive account of the classification and function of 5-hydroxytryptamine receptors, including how they are targeted for therapeutic benefit.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor Channel

Ion channel proteins form water-filled nanoscale pores within lipid bilayers, and their properties are dependent on the complex behavior of water in a nanoconfined environment. Using a simplified model of the pore of the 5-HT3 receptor (5HT3R) which restrains the backbone structure to that of the parent channel protein from which it is derived, we compare additive with polarizable models in describing the behavior of water in nanopores. Molecular dynamics simulations were performed with four conformations of the channel: two closed state structures, an intermediate state, and an open state, each embedded in a phosphatidylcholine bilayer. Water density profiles revealed that for all water models, the closed and intermediate states exhibited strong dewetting within the central hydrophobic gate region of the pore. However, the open state conformation exhibited varying degrees of hydration, ranging from partial wetting for the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. Water dipole moments calculated using polarizable force fields also revealed that water molecules remaining within dewetted sections of the pore resemble gas phase water. Free energy profiles for Na+ and for Cl- ions within the open state pore revealed more rugged energy landscapes using polarizable force fields, and the hydration number profiles of these ions were also sensitive to induced polarization resulting in a substantive reduction of the number of waters within the first hydration shell of Cl- while it permeates the pore. These results demonstrate that induced polarization can influence the complex behavior of water and ions within nanoscale pores and provides important new insights into their chemical properties.

Activation and Signaling Mechanism Revealed by Cannabinoid Receptor-Gi Complex Structures

Human endocannabinoid systems modulate multiple physiological processes mainly through the activation of cannabinoid receptors CB1 and CB2. Their high sequence similarity, low agonist selectivity, and lack of activation and G protein-coupling knowledge have hindered the development of therapeutic applications. Importantly, missing structural information has significantly held back the development of promising CB2-selective agonist drugs for treating inflammatory and neuropathic pain without the psychoactivity of CB1. Here, we report the cryoelectron microscopy structures of synthetic cannabinoid-bound CB2 and CB1 in complex with G_i, as well as agonist-bound CB2 crystal structure. Of important scientific and therapeutic benefit, our results reveal a diverse activation and signaling mechanism, the structural basis of CB2-selective agonists design, and the unexpected interaction of cholesterol with CB1, suggestive of its endogenous allosteric modulating role.

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

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

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.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

CHAP: A Versatile Tool for the Structural and Functional Annotation of Ion Channel Pores

The control of ion channel permeation requires the modulation of energetic barriers or “gates” within their pores. However, such barriers are often simply identified from the physical dimensions of the pore. Such approaches have worked well in the past, but there is now evidence that the unusual behavior of water within narrow hydrophobic pores can produce an energetic barrier to permeation without requiring steric occlusion of the pathway. Many different ion channels have now been shown to exploit “hydrophobic gating” to regulate ion flow, and it is clear that new tools are required for more accurate functional annotation of the increasing number of ion channel structures becoming available. We have previously shown how molecular dynamics simulations of water can be used as a proxy to predict hydrophobic gates, and we now present a new and highly versatile computational tool, the Channel Annotation Package (CHAP) that implements this methodology.

•

Some ion channels exhibit hydrophobic gating via dewetting of the central pore.

•

This cannot be predicted from their pore radius alone.

•

The hydrophobicity of a pore contributes to this dewetting effect.

•

We have developed a new tool (CHAP) that combines these measurements.

•

CHAP now enables the rapid functional annotation of ion channel structures.

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.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Trans- Cis Proline Switches in a Pentameric Ligand-Gated Ion Channel: How They Are Affected by and How They Affect the Biomolecular Environment

Pentameric ligand-gated ion channels (pLGICs) are important neuroreceptors, embedded in neuronal membranes, that mediate fast synaptic transmission. The molecular details of their working mechanisms have still to be fully unravelled due to their complexity and limited structural information available. Here we focus on a potential molecular switch in a prototypical pLGIC, the serotonin-activated 5-HT₃ receptor, consisting of the trans- cis isomerization of a proline at the interface between the extracellular and transmembrane domain. Mutagenesis electrophysiology experiments previously showed that if such isomerization could not take place, the channel would not open, but the hypothetical role of this mechanism as key to channel gating is still debated. We investigate this switch within the receptor with molecular dynamics and enhanced sampling simulations. We analyze how the isomerization free energy landscape is affected by the receptor environment in comparison to simplified models. Moreover, we reveal how the isomerization, in turn, affects the structural and electrostatic properties of the receptor at the extracellular-transmembrane domain interface, e.g., by tuning the ion selectivity filter.

Cryo-EM reveals two distinct serotonin-bound conformations of full-length 5-HT3A receptor

Serotonin receptor (5-HT_3AR)¹, a cationic pentameric ligand-gated ion channel (pLGIC), is the clinical target for management of nausea and vomiting associated with radiation and chemotherapies². Upon binding, serotonin induces a global conformational change encompassing the ligand-binding extracellular domain (ECD), the transmembrane domain (TMD), and the intracellular domain (ICD), the molecular details of which are unclear. Here, we present two serotonin-bound structures of the full-length 5-HT_3AR in distinct conformations at 3.32 Å and 3.89 Å resolutions that reveal the mechanism underlying channel activation. When compared to Apo-5-HT_3AR, serotonin-bound states underwent a large twisting motion in the ECD and TMD leading to the opening of a 165 Å long permeation pathway. Notably, this motion results in creation of lateral portals for ion permeation at the interface of the TMD and ICD. Combined with molecular dynamics simulations, these structures provide novel insights into conformational coupling across domains and functional modulation.

Approaches for Differentiation and Interconverting GPCR Agonists and Antagonists

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

The C and E subunits of the serotonin 5-HT3 receptor subtly modulate electrical properties of the receptor

Serotonin type 3 (5-hydroxytrptamine-3, 5-HT₃) receptors are ligand-gated cation channels present in both central and peripheral nervous systems. In humans there are five different subunits (A, B, C, D and E) of 5-HT₃ receptors which can form homomeric or heteromeric receptors that may account for discrepancies in patient responses to treatments. The present study commences characterisation of the profiles of human 5-HT₃ receptors containing C and/or E subunits. Recombinant 5-HT₃ receptors were expressed transiently in HEK293T cells and expression was checked via immunocytochemistry staining against each epitope-tagged subunits. Functional characterisation of different combinations of 5-HT₃ receptor complexes was studied via patch clamp whole cell recordings. In this study, increased current was seen in cells containing A and C subunits but only subtle changes were seen in the electrical properties of cells expressing A, AE, or ACE subunits in response to the ligand, 5-HT. Both di- and tri-heteromeric 5-HT₃ receptors were significantly inhibited by the antagonists, ondansetron and palonosetron. Notably, palonosetron exerted stronger and more rapid inhibition on the 5-HT₃ receptor ACE tri-heteromer compared to homomeric and di-heteromeric counterparts. This study demonstrated that the C and E subunits when assembled as simple or complex heteromeric 5-HT₃ receptors may alter efficacies of 5-HT and clinically used antagonists such as ondansetron and palonosetron, and this in turn may have implications for patient responses to therapies.

Ligand-induced action of the W2866.48 rotamer toggle switch in the β2-adrenergic receptor

Studies focused on GPCRs, particularly on the β₂-adrenergic receptor (β₂-AR), have demonstrated the relationship between ligand structure, receptor conformational changes and the corresponding pharmacological outcomes. Herein, we studied the molecular details of the rotameric flip of the W286^6.48 sidechain, i.e. a presumed action switch that has not been reported in native β₂-AR thus far. It is believed that although both the 'active' and 'inactive' conformers of β₂-AR exhibit similar conformations of this switch, it may still play a substantial role in the ligand-induced activation of the receptor. By using both experimental methods (time-resolved fluorescence spectroscopy) and molecular modeling techniques (enhanced-sampling molecular dynamics), we characterized the conformational rearrangements of W286^6.48 in relation to the type of ligand present in the binding cavity and to the conformation of the receptor ('active' vs. 'inactive' β₂-AR). We found that the conformational behaviour of W286^6.48 is correlated with the pharmacological character of the ligand present in the binding cavity but not with the instantaneous conformation of the receptor. Namely, agonists promote the W286^6.48 conformations that facilitate the increase of the solvation within the inner receptor channel. In contrast, antagonists and inverse agonists act toward the decrease of the solvation in the inner channel. This creates an opportunity for using computational methodologies in determining the pharmacological properties of various ligands. The combination of the time-resolved fluorescence spectroscopy technique with the enhanced-sampling molecular dynamics simulations is shown to be a powerful tool for studying the ligand-induced conformational rearrangements in GPCRs.

Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter L-glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The μs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the β1-β2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with L-glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.

Functional Annotation of Ion Channel Structures by Molecular Simulation

Ion channels play key roles in cell membranes, and recent advances are yielding an increasing number of structures. However, their functional relevance is often unclear and better tools are required for their functional annotation. In sub-nanometer pores such as ion channels, hydrophobic gating has been shown to promote dewetting to produce a functionally closed (i.e., non-conductive) state. Using the serotonin receptor (5-HT₃R) structure as an example, we demonstrate the use of molecular dynamics to aid the functional annotation of channel structures via simulation of the behavior of water within the pore. Three increasingly complex simulation analyses are described: water equilibrium densities; single-ion free-energy profiles; and computational electrophysiology. All three approaches correctly predict the 5-HT₃R crystal structure to represent a functionally closed (i.e., non-conductive) state. We also illustrate the application of water equilibrium density simulations to annotate different conformational states of a glycine receptor.

Superagonist, Full Agonist, Partial Agonist, and Antagonist Actions of Arylguanidines at 5-Hydroxytryptamine-3 (5-HT3) Subunit A Receptors

Introduction of minor variations to the substitution pattern of arylguanidine 5-hydroxytryptamine-3 (5-HT₃) receptor ligands resulted in a broad spectrum of functionally-active ligands from antagonist to superagonist. For example, (i) introduction of an additional Cl-substituent(s) to our lead full agonist N-(3-chlorophenyl)guanidine (mCPG, 2; efficacy % = 106) yielded superagonists 7-9 (efficacy % = 186, 139, and 129, respectively), (ii) a positional isomer of 2, p-Cl analog 11, displayed partial agonist actions (efficacy % = 12), and (iii) replacing the halogen atom at the meta or para position with an electron donating OCH₃ group or a stronger electron withdrawing (i.e., CF₃) group resulted in antagonists 13-16. We posit based on combined mutagenesis, crystallographic, and computational analyses that for the 5-HT₃ receptor, the arylguanidines that are better able to simultaneously engage the primary and complementary subunits, thus keeping them in close proximity, have greater agonist character while those that are deficient in this ability are antagonists.

Two-dimensional crystallization of the mouse serotonin 5-HT3A receptor

The mouse serotonin 5-HT_3A receptor is a homo-pentameric ligand-gated ion channel (pLGIC) mediating fast excitatory neurotransmission in the central nervous system. The molecular mechanism of ion permeation of 5-HT_3A receptors triggered by the neurotransmitter serotonin is not yet fully understood. The recent X-ray structure of the mouse serotonin 5-HT_3A receptor in complex with a stabilizing nanobody revealed for the first time the entire structure of a mammalian pLGIC in detergent. Structural information of the receptor in a lipid bilayer however is still limited primarily due to the lack of 2D crystals of the receptor in a lipid bilayer. Here we present our results on the formation and improvement of diffracting 2D crystals of the mouse 5-HT_3A by limited proteolysis and addition of conformational nanobodies.

The Molecular Mechanism of P2Y1 Receptor Activation

Human purinergic G protein-coupled receptor P2Y1 (P2Y1 R) is activated by adenosine 5'-diphosphate (ADP) to induce platelet activation and thereby serves as an important antithrombotic drug target. Crystal structures of P2Y1 R revealed that one ligand (MRS2500) binds to the extracellular vestibule of this GPCR, whereas another (BPTU) occupies the surface between transmembrane (TM) helices TM2 and TM3. We introduced a total of 20 μs all-atom long-timescale molecular dynamic (MD) simulations to inquire why two molecules in completely different locations both serve as antagonists while ADP activates the receptor. Our results indicate that BPTU acts as an antagonist by stabilizing extracellular helix bundles leading to an increase of the lipid order, whereas MRS2500 blocks signaling by occupying the ligand binding site. Both antagonists stabilize an ionic lock within the receptor. However, binding of ADP breaks this ionic lock, forming a continuous water channel that leads to P2Y1 R activation.