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

Probing the energy landscape of the lipid interactions of the serotonin1A receptor

G protein-coupled receptors (GPCRs) are lipid-regulated transmembrane proteins that play a central role in cell signaling and pharmacology. Although the role of membrane lipids in GPCR function is well established, the underlying GPCR-lipid interactions have not been thermodynamically characterized due to the complexity of these interactions. In this work, we estimate the energetics and dynamics of lipid association from coarse-grain simulations of the serotonin1A receptor embedded in a complex membrane. We show that lipids bind to the receptor with varying energetics of 1-4 kT, and timescales of 1-10 μs. The most favorable energetics and longest residence times are observed for cholesterol, glycosphingolipid GM1, phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids. Multi-exponential fitting of the contact probability suggests distinct dynamic regimes, corresponding to ps, ns and μs timescales, that we correlate with the annular, intermediate and non-annular lipid sites. The timescales of lipid binding correspond to high barrier heights, despite their relatively weaker energetics. Our results highlight that GPCR-lipid interactions are driven by both thermodynamic interactions and the dynamical features of lipid binding.

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.

The density-threshold affinity: Calculating lipid binding affinities from unbiased coarse-grained molecular dynamics simulations

Many membrane proteins are sensitive to their local lipid environment. As structural methods for membrane proteins have improved, there is growing evidence of direct, specific binding of lipids to protein surfaces. Unfortunately the workhorse of understanding protein-small molecule interactions, the binding affinity for a given site, is experimentally inaccessible for these systems. Coarse-grained molecular dynamics simulations can be used to bridge this gap, and are relatively straightforward to learn. Such simulations allow users to observe spontaneous binding of lipids to membrane proteins and quantify localized densities of individual lipids or lipid fragments. In this chapter we outline a protocol for extracting binding affinities from these localized distributions, known as the "density threshold affinity." The density threshold affinity uses an adaptive and flexible definition of site occupancy that alleviates the need to distinguish between "bound'' lipids and bulk lipids that are simply diffusing through the site. Furthermore, the method allows "bead-level" resolution that is suitable for the case where lipids share binding sites, and circumvents ambiguities about a relevant reference state. This approach provides a convenient and straightforward method for comparing affinities of a single lipid species for multiple sites, multiple lipids for a single site, and/or a single lipid species modeled using multiple forcefields.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

Structure and function meet at the nicotinic acetylcholine receptor-lipid interface

The nicotinic acetylcholine receptor (nAChR) is a transmembrane protein that mediates fast intercellular communication in response to the endogenous neurotransmitter acetylcholine. It is the best characterized and archetypal molecule in the superfamily of pentameric ligand-gated ion channels (pLGICs). As a typical transmembrane macromolecule, it interacts extensively with its vicinal lipid microenvironment. Experimental evidence provides a wealth of information on receptor-lipid crosstalk: the nAChR exerts influence on its immediate membrane environment and conversely, the lipid moiety modulates ligand binding, affinity state transitions and gating of ion translocation functions of the receptor protein. Recent cryogenic electron microscopy (cryo-EM) studies have unveiled the occurrence of sites for phospholipids and cholesterol on the lipid-exposed regions of neuronal and electroplax nAChRs, confirming early spectroscopic and affinity labeling studies demonstrating the close contact of lipid molecules with the receptor transmembrane segments. This new data provides structural support to the postulated "lipid sensor" ability displayed by the outer ring of M4 transmembrane domains and their modulatory role on nAChR function, as we postulated a decade ago. Borrowing from the best characterized nAChR, the electroplax (muscle-type) receptor, and exploiting new structural information on the neuronal nAChR, it is now possible to achieve an improved depiction of these sites. In combination with site-directed mutagenesis, single-channel electrophysiology, and molecular dynamics studies, the new structural information delivers a more comprehensive portrayal of these lipid-sensitive loci, providing mechanistic explanations for their ability to modulate nAChR properties and raising the possibility of targetting them in disease.

Polyunsaturated fatty acids inhibit a pentameric ligand-gated ion channel through one of two binding sites

Polyunsaturated fatty acids (PUFAs) inhibit pentameric ligand-gated ion channels (pLGICs) but the mechanism of inhibition is not well understood. The PUFA, docosahexaenoic acid (DHA), inhibits agonist responses of the pLGIC, ELIC, more effectively than palmitic acid, similar to the effects observed in the GABA_A receptor and nicotinic acetylcholine receptor. Using photo-affinity labeling and coarse-grained molecular dynamics simulations, we identified two fatty acid binding sites in the outer transmembrane domain (TMD) of ELIC. Fatty acid binding to the photolabeled sites is selective for DHA over palmitic acid, and specific for an agonist-bound state. Hexadecyl-methanethiosulfonate modification of one of the two fatty acid binding sites in the outer TMD recapitulates the inhibitory effect of PUFAs in ELIC. The results demonstrate that DHA selectively binds to multiple sites in the outer TMD of ELIC, but that state-dependent binding to a single intrasubunit site mediates DHA inhibition of ELIC.

Spontaneous lipid binding to the nicotinic acetylcholine receptor in a native membrane

The nicotinic acetylcholine receptor (nAChR) and other pentameric ligand-gated ion channels are native to neuronal membranes with an unusual lipid composition. While it is well-established that these receptors can be significantly modulated by lipids, the underlying mechanisms have been primarily studied in model membranes with few lipid species. Here, we use coarse-grained molecular dynamics simulation to probe specific binding of lipids in a complex quasi-neuronal membrane. We ran a total of 50 μs of simulations of a single nAChR in a membrane composed of 36 species of lipids. Competition between multiple lipid species produces a complex distribution. We find that overall, cholesterol selects for concave inter-subunit sites and polyunsaturated fatty acids select for convex M4 sites, while monounsaturated and saturated lipids are unenriched in the nAChR boundary. We propose the "density-threshold affinity" as a metric calculated from continuous density distributions, which reduces to a standard affinity in two-state binding. We find that the density-threshold affinity for M4 weakens with chain rigidity, which suggests that flexible chains may help relax packing defects caused by the conical protein shape. For any site, PE headgroups have the strongest affinity of all phospholipid headgroups, but anionic lipids still yield moderately high affinities for the M4 sites as expected. We observe cooperative effects between anionic headgroups and saturated chains at the M4 site in the inner leaflet. We also analyze affinities for individual anionic headgroups. When combined, these insights may reconcile several apparently contradictory experiments on the role of anionic phospholipids in modulating nAChR.

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.

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

Direct binding of phosphatidylglycerol at specific sites modulates desensitization of a ligand-gated ion channel

Pentameric ligand-gated ion channels (pLGICs) are essential determinants of synaptic transmission, and are modulated by specific lipids including anionic phospholipids. The exact modulatory effect of anionic phospholipids in pLGICs and the mechanism of this effect are not well understood. Using native mass spectrometry, coarse-grained molecular dynamics simulations and functional assays, we show that the anionic phospholipid, 1-palmitoyl-2-oleoyl phosphatidylglycerol (POPG), preferentially binds to and stabilizes the pLGIC, Erwinia ligand-gated ion channel (ELIC), and decreases ELIC desensitization. Mutations of five arginines located in the interfacial regions of the transmembrane domain (TMD) reduce POPG binding, and a subset of these mutations increase ELIC desensitization. In contrast, a mutation that decreases ELIC desensitization, increases POPG binding. The results support a mechanism by which POPG stabilizes the open state of ELIC relative to the desensitized state by direct binding at specific sites.

L-3,3',5-triiodothyronine and pregnenolone sulfate inhibit Torpedo nicotinic acetylcholine receptors

The nicotinic acetylcholine receptor (nAChR) is an excitatory pentameric ligand-gated ion channel (pLGIC), homologous to the inhibitory γ-aminobutyric acid (GABA) type A receptor targeted by pharmaceuticals and endogenous sedatives. Activation of the GABA_A receptor by the neurosteroid allopregnanolone can be inhibited competitively by thyroid hormone (L-3,3’,5-triiodothyronine, or T3), but modulation of nAChR by T3 or neurosteroids has not been investigated. Here we show that allopregnanolone inhibits the nAChR from Torpedo californica at micromolar concentrations, as do T3 and the anionic neurosteroid pregnenolone sulfate (PS). We test for the role of protein and ligand charge in mediated receptor inhibition by varying pH in a narrow range around physiological pH. We find that both T3 and PS become less potent with increasing pH, with remarkably similar trends in IC₅₀ when T3 is neutral at pH < 7.3. After deprotonation of T3 (but no additional deprotonation of PS) at pH 7.3, T3 loses potency more slowly with increasing pH than PS. We interpret this result as indicating the negative charge is not required for inhibition but does increase activity. Finally, we show that both T3 and PS affect nAChR channel desensitization, which may implicate a binding site homologous to one that was recently indicated for accelerated desensitization of the GABA_A receptor by PS.

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.

Alzheimer's Disease as a Membrane Disorder: Spatial Cross-Talk Among Beta-Amyloid Peptides, Nicotinic Acetylcholine Receptors and Lipid Rafts

Biological membranes show lateral and transverse asymmetric lipid distribution. Cholesterol (Chol) localizes in both hemilayers, but in the external one it is mostly condensed in lipid-ordered microdomains (raft domains), together with saturated phosphatidyl lipids and sphingolipids (including sphingomyelin and glycosphingolipids). Membrane asymmetries induce special membrane biophysical properties and behave as signals for several physiological and/or pathological processes. Alzheimer’s disease (AD) is associated with a perturbation in different membrane properties. Amyloid-β (Aβ) plaques and neurofibrillary tangles of tau protein together with neuroinflammation and neurodegeneration are the most characteristic cellular changes observed in this disease. The extracellular presence of Aβ peptides forming senile plaques, together with soluble oligomeric species of Aβ, are considered the major cause of the synaptic dysfunction of AD. The association between Aβ peptide and membrane lipids has been extensively studied. It has been postulated that Chol content and Chol distribution condition Aβ production and posterior accumulation in membranes and, hence, cell dysfunction. Several lines of evidence suggest that Aβ partitions in the cell membrane accumulate mostly in raft domains, the site where the cleavage of the precursor AβPP by β- and γ- secretase is also thought to occur. The main consequence of the pathogenesis of AD is the disruption of the cholinergic pathways in the cerebral cortex and in the basal forebrain. In parallel, the nicotinic acetylcholine receptor has been extensively linked to membrane properties. Since its transmembrane domain exhibits extensive contacts with the surrounding lipids, the acetylcholine receptor function is conditioned by its lipid microenvironment. The nicotinic acetylcholine receptor is present in high-density clusters in the cell membrane where it localizes mainly in lipid-ordered domains. Perturbations of sphingomyelin or cholesterol composition alter acetylcholine receptor location. Therefore, Aβ processing, Aβ partitioning, and acetylcholine receptor location and function can be manipulated by changes in membrane lipid biophysics. Understanding these mechanisms should provide insights into new therapeutic strategies for prevention and/or treatment of AD. Here, we discuss the implications of lipid-protein interactions at the cell membrane level in AD.

Essay on Biomembrane Structure

Of all the macromolecular assemblies of life, the least understood is the biomembrane. This is especially true in regard to its atomic structure. Ideas on biomembranes, developed in the last 200 years, culminated in the fluid mosaic model of the membrane. In this essay, I provide a historical outline of how we arrived at our current understanding of biomembranes and the models we use to describe them. A selection of direct experimental findings on the nano-scale structure of biomembranes is taken up to discuss their physical nature, and special emphasis is put on the surprising insights that arise from atomic scale descriptions.