Anionic lipid and cholesterol interactions with alpha4beta2 nAChR: insights from MD simulations

State-dependent binding of cholesterol and an anionic lipid to the muscle-type Torpedo nicotinic acetylcholine receptor

The ability of the Torpedo nicotinic acetylcholine receptor (nAChR) to undergo agonist-induced conformational transitions requires the presence of cholesterol and/or anionic lipids. Here we use recently solved structures along with multiscale molecular dynamics simulations to examine lipid binding to the nAChR in bilayers that have defined effects on nAChR function. We examine how phosphatidic acid and cholesterol, lipids that support conformational transitions, individually compete for binding with phosphatidylcholine, a lipid that does not. We also examine how the two lipids work synergistically to stabilize an agonist-responsive nAChR. We identify rapidly exchanging lipid binding sites, including both phospholipid sites with a high affinity for phosphatidic acid and promiscuous cholesterol binding sites in the grooves between adjacent transmembrane α-helices. A high affinity cholesterol site is confirmed in the inner leaflet framed by a key tryptophan residue on the MX α-helix. Our data provide insight into the dynamic nature of lipid-nAChR interactions and set the stage for a detailed understanding of the mechanisms by which lipids facilitate nAChR function at the neuromuscular junction.

Structural Elucidation of Ivermectin Binding to α7nAChR and the Induced Channel Desensitization

The α7 nicotinic acetylcholine receptor (α7nAChR) mediates signaling in the central nervous system and cholinergic anti-inflammatory pathways. Ivermectin is a positive allosteric modulator of a full-length α7nAChR and an agonist of the α7nAChR construct containing transmembrane (TMD) and intracellular (ICD) domains, but structural insights of the binding have not previously been determined. Here, combining nuclear magnetic resonance as a primary experimental tool with Rosetta comparative modeling and molecular dynamics simulations, we have revealed details of ivermectin binding to the α7nAChR TMD + ICD and corresponding structural changes in an ivermectin-induced desensitized state. Ivermectin binding was stabilized predominantly by hydrophobic interactions from interfacial residues between adjacent subunits near the extracellular end of the TMD, where the inter-subunit gap was substantially expanded in comparison to the apo structure. The ion-permeation pathway showed a profile distinctly different from the resting-state profile but similar to profiles of desensitized α7nAChR. The ICD also exhibited structural changes, including reorientation of the MX and h3 helices relative to the channel axis. The resulting structures of the α7nAChR TMD + ICD in complex with ivermectin provide opportunities for discovering new modulators of therapeutic potential and exploring the structural basis of cytoplasmic signaling under different α7nAChR functional states.

State-dependent protein-lipid interactions of a pentameric ligand-gated ion channel in a neuronal membrane

Pentameric ligand-gated ion channels (pLGICs) are receptor proteins that are sensitive to their membrane environment, but the mechanism for how lipids modulate function under physiological conditions in a state dependent manner is not known. The glycine receptor is a pLGIC whose structure has been resolved in different functional states. Using a realistic model of a neuronal membrane coupled with coarse-grained molecular dynamics simulations, we demonstrate that some key lipid-protein interactions are dependent on the receptor state, suggesting that lipids may regulate the receptor's conformational dynamics. Comparison with existing structural data confirms known lipid binding sites, but we also predict further protein-lipid interactions including a site at the communication interface between the extracellular and transmembrane domain. Moreover, in the active state, cholesterol can bind to the binding site of the positive allosteric modulator ivermectin. These protein-lipid interaction sites could in future be exploited for the rational design of lipid-like allosteric drugs.

Sequential purification and characterization of Torpedo californica nAChR-DC supplemented with CHS for high-resolution crystallization studies

Over the past 10 years we have been developing a multi-attribute analytical platform that allows for the preparation of milligram amounts of functional, high-pure, and stable Torpedo (muscle-type) nAChR detergent complexes for crystallization purpose. In the present work, we have been able to significantly improve and optimize the purity and yield of nicotinic acetylcholine receptors in detergent complexes (nAChR-DC) without compromising stability and functionality. We implemented new methods in the process, such as analysis and rapid production of samples for future crystallization preparations. Native nAChR was extracted from the electric organ of Torpedo californica using the lipid-like detergent LysoFos Choline 16 (LFC-16), followed by three consecutive steps of chromatography purification. We evaluated the effect of cholesteryl hemisuccinate (CHS) supplementation during the affinity purification steps of nAChR-LFC-16 in terms of receptor secondary structure, stability and functionality. CHS produced significant changes in the degree of β-secondary structure, these changes compromise the diffusion of the nAChR-LFC-16 in lipid cubic phase. The behavior was reversed by Methyl-β-Cyclodextrin treatment. Also, CHS decreased acetylcholine evoked currents of Xenopus leavis oocyte injected with nAChR-LFC-16 in a concentration-dependent manner. Methyl-β-Cyclodextrin treatment do not reverse functionality, however column delipidation produced a functional protein similar to nAChR-LFC-16 without CHS treatment.

A Review on the Receptor-ligand Molecular Interactions in the Nicotinic Receptor Signaling Systems

Nicotine is regarded as the main active addictive ingredient in tobacco products driving continued tobacco abuse behavior (smoking) to the addiction behavior, whereas nicotinic acetylcholine receptors (nAChR) is the crucial effective apparatus or molecular effector of nicotine and acetylcholine and other similar ligands. Many nAChR subunits have been revealed to bind to either neurotransmitters or exogenous ligands, such as nicotine and acetylcholine, being involved in the nicotinic receptor signal transduction. Therefore, the nicotinic receptor signalling molecules and the receptor-ligand molecular interactions between nAChRs and their ligands are universally regarded as crucial mediators of cellular functions and drug targets in medical treatment and clinical diagnosis. Given numerous endeavours have been made in defining the roles of nAChRs in response to nicotine and other addictive drugs, this review focuses on studies and reports in recent years on the receptor-ligand interactions between nAChR receptors and ligands, including lipid-nAChR and protein-nAChR molecular interactions, relevant signal transduction pathways and their molecular mechanisms in the nicotinic receptor signalling systems. All the references were carefully retrieved from the PubMed database by searching key words "nicotine", "acetylcholine", "nicotinic acetylcholine receptor(s)", "nAChR*", "protein and nAChR", "lipid and nAChR", "smok*" and "tobacco". All the relevant referred papers and reports retrieved were fully reviewed for manual inspection. This effort intend to get a quick insight and understanding of the nicotinic receptor signalling and their molecular interactions mechanisms. Understanding the cellular receptor-ligand interactions and molecular mechanisms between nAChRs and ligands will lead to a better translational and therapeutic operations and outcomes for the prevention and treatment of nicotine addiction and other chronic drug addictions in the brain's reward circuitry.

Untangling Direct and Domain-Mediated Interactions Between Nicotinic Acetylcholine Receptors in DHA-Rich Membranes

At the neuromuscular junction (NMJ), the nicotinic acetylcholine receptor (nAChR) self-associates to give rise to rapid muscle movement. While lipid domains have maintained nAChR aggregates in vitro, their specific roles in nAChR clustering are currently unknown. In the present study, we carried out coarse-grained molecular dynamics simulations (CG-MD) of 1-4 nAChR molecules in two membrane environments: one mixture containing domain-forming, homoacidic lipids, and a second mixture consisting of heteroacidic lipids. Spontaneous dimerization of nAChRs was up to ten times more likely in domain-forming membranes; however, the effect was not significant in four-protein systems, suggesting that lipid domains are less critical to nAChR oligomerization when protein concentration is higher. With regard to lipid preferences, nAChRs consistently partitioned into liquid-disordered domains occupied by the omega-3 ([Formula: see text]-3) fatty acid, docosahexaenoic acid (DHA); enrichment of DHA boundary lipids increased with protein concentration, particularly in homoacidic membranes. This result suggests dimer formation blocks access of saturated chains and cholesterol, but not polyunsaturated chains, to boundary lipid sites.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Boundary lipids of the nicotinic acetylcholine receptor: Spontaneous partitioning via coarse-grained molecular dynamics simulation

Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant gain-of-function upon addition of cholesterol to reconstitution mixtures, and cholesterol affects the organization of nAChRs within domain-forming membranes, but whether nAChR partitions to cholesterol-rich liquid-ordered ("raft" or l_o) domains or cholesterol-poor liquid-disordered (l_do) domains is unknown. We use coarse-grained molecular dynamics simulations to observe spontaneous interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA) lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol (DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed spontaneously entering deep "non-annular" cavities in the nAChR TMD, particularly at the subunit interface and the β subunit center, facilitated by the low amino acid density in the cryo-EM structure of nAChR in a native membrane. Cholesterol was highly enriched in the annulus around the TMD, but this effect extended over (at most) 5-10 Å. In domain-forming ternary mixtures containing PUFAs, the presence of a single receptor did not significantly affect the likelihood of domain formation. nAChR partitioned to any cholesterol-poor l_do domain that was present, regardless of whether the l_do or l_o domain lipids had PC or PE headgroups. Enrichment of PUFAs among boundary lipids was positively correlated with their propensity for demixing from cholesterol-rich phases. Long n-3 chains (tested here with Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular embedded sites, partially displacing cholesterol and completely displacing DPPC, and occupying sites even deeper within the bundle. Shorter n-6 chains were far less effective at displacing cholesterol from non-annular sites.

Lipid interaction sites on channels, transporters and receptors: Recent insights from molecular dynamics simulations

Lipid molecules are able to selectively interact with specific sites on integral membrane proteins, and modulate their structure and function. Identification and characterization of these sites are of importance for our understanding of the molecular basis of membrane protein function and stability, and may facilitate the design of lipid-like drug molecules. Molecular dynamics simulations provide a powerful tool for the identification of these sites, complementing advances in membrane protein structural biology and biophysics. We describe recent notable biomolecular simulation studies which have identified lipid interaction sites on a range of different membrane proteins. The sites identified in these simulation studies agree well with those identified by complementary experimental techniques. This demonstrates the power of the molecular dynamics approach in the prediction and characterization of lipid interaction sites on integral membrane proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Competing Lipid-Protein and Protein-Protein Interactions Determine Clustering and Gating Patterns in the Potassium Channel from Streptomyces lividans (KcsA)

There is increasing evidence to support the notion that membrane proteins, instead of being isolated components floating in a fluid lipid environment, can be assembled into supramolecular complexes that take part in a variety of cooperative cellular functions. The interplay between lipid-protein and protein-protein interactions is expected to be a determinant factor in the assembly and dynamics of such membrane complexes. Here we report on a role of anionic phospholipids in determining the extent of clustering of KcsA, a model potassium channel. Assembly/disassembly of channel clusters occurs, at least partly, as a consequence of competing lipid-protein and protein-protein interactions at nonannular lipid binding sites on the channel surface and brings about profound changes in the gating properties of the channel. Our results suggest that these latter effects of anionic lipids are mediated via the Trp(67)-Glu(71)-Asp(80) inactivation triad within the channel structure and its bearing on the selectivity filter.

The influence of cholesterol on membrane protein structure, function, and dynamics studied by molecular dynamics simulations

The plasma membrane, which encapsulates human cells, is composed of a complex mixture of lipids and embedded proteins. Emerging knowledge points towards the lipids as having a regulating role in protein function. Furthermore, insight from protein crystallography has revealed several different types of lipids intimately bound to membrane proteins and peptides, hereby possibly pointing to a site of action for the observed regulation. Cholesterol is among the lipid membrane constituents most often observed to be co-crystallized with membrane proteins, and the cholesterol levels in cell membranes have been found to play an essential role in health and disease. Remarkably little is known about the mechanism of lipid regulation of membrane protein function in health as well as in disease. Herein, we review molecular dynamics simulation studies aimed at investigating the effect of cholesterol on membrane protein and peptide properties. This article is part of a Special Issue entitled: Lipid-protein interactions.

A predicted binding site for cholesterol on the GABAA receptor

Modulation of the GABA type A receptor (GABAAR) function by cholesterol and other steroids is documented at the functional level, yet its structural basis is largely unknown. Current data on structurally related modulators suggest that cholesterol binds to subunit interfaces between transmembrane domains of the GABAAR. We construct homology models of a human GABAAR based on the structure of the glutamate-gated chloride channel GluCl of Caenorhabditis elegans. The models show the possibility of previously unreported disulfide bridges linking the M1 and M3 transmembrane helices in the α and γ subunits. We discuss the biological relevance of such disulfide bridges. Using our models, we investigate cholesterol binding to intersubunit cavities of the GABAAR transmembrane domain. We find that very similar binding modes are predicted independently by three approaches: analogy with ivermectin in the GluCl crystal structure, automated docking by AutoDock, and spontaneous rebinding events in unbiased molecular dynamics simulations. Taken together, the models and atomistic simulations suggest a somewhat flexible binding mode, with several possible orientations. Finally, we explore the possibility that cholesterol promotes pore opening through a wedge mechanism.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

New Insights into the Mechanisms of Action of Cotinine and its Distinctive Effects from Nicotine

Tobacco consumption is far higher among a number of psychiatric and neurological diseases, supporting the notion that some component(s) of tobacco may underlie the oft-reported reduction in associated symptoms during tobacco use. Popular dogma holds that this component is nicotine. However, increasing evidence support theories that cotinine, the main metabolite of nicotine, may underlie at least some of nicotine's actions in the nervous system, apart from its adverse cardiovascular and habit forming effects. Though similarities exist, disparate and even antagonizing actions between cotinine and nicotine have been described both in terms of behavior and physiology, underscoring the need to further characterize this potentially therapeutic compound. Cotinine has been shown to be psychoactive in humans and animals, facilitating memory, cognition, executive function, and emotional responding. Furthermore, recent research shows that cotinine acts as an antidepressant and reduces cognitive-impairment associated with disease and stress-induced dysfunction. Despite these promising findings, continued focus on this potentially safe alternative to tobacco and nicotine use is lacking. Here, we review the effects of cotinine, including comparisons with nicotine, and discuss potential mechanisms of cotinine-specific actions in the central nervous system which are, to date, still being elucidated.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.

Disclosure of cholesterol recognition motifs in transmembrane domains of the human nicotinic acetylcholine receptor

Cholesterol influences ion-channel function, distribution and clustering in the membrane, endocytosis, and exocytic sorting of the nicotinic acetylcholine receptor (AChR). We report the occurrence of a cholesterol recognition motif, here coined “CARC”, in the transmembrane regions of AChR subunits that bear extensive contact with the surrounding lipid, and are thus optimally suited to convey cholesterol-mediated signaling from the latter. Three cholesterol molecules could be docked on the transmembrane segments of each AChR subunit, rendering a total of 15 cholesterol molecules per AChR molecule. The CARC motifs contribute each with an energy of interaction between 35 and 52 kJ.mol⁻¹, adding up to a total of about 200 kJ.mol⁻¹ per receptor molecule, i.e. ∼40% of the lipid solvation free energy/ AChR molecule. The CARC motif is remarkably conserved along the phylogenetic scale, from prokaryotes to human, suggesting that it could be responsible for some of the above structural/functional properties of the AChR.

Membrane interface composition drives the structure and the tilt of the single transmembrane helix protein PMP1: MD studies

PMP1, a regulatory subunit of the yeast plasma membrane H(+)-ATPase, is a single transmembrane helix protein. Its cytoplasmic C-terminus possesses several positively charged residues and interacts with phosphatidylserine lipids as shown through both (1)H- and (2)H-NMR experiments. We used all-atom molecular dynamics simulations to obtain atomic-scale data on the effects of membrane interface lipid composition on PMP1 structure and tilt. PMP1 was embedded in two hydrated bilayers, differing in the composition of the interfacial region. The neutral bilayer is composed of POPC (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylcholine) lipids and the negatively charged bilayer is composed of POPC and anionic POPS (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylserine) lipids. Our results were consistent with NMR data obtained previously, such as a lipid sn-2 chain lying on the W28 aromatic ring and in the groove formed on one side of the PMP1 helix. In pure POPC, the transmembrane helix is two residues longer than the initial structure and the helix tilt remains constant at 6 ± 3°. By contrast, in mixed POPC-POPS, the initial helical structure of PMP1 is stable throughout the simulation time even though the C-terminal residues interact strongly with POPS headgroups, leading to a significant increase of the helix tilt within the membrane to 20 ± 5°.