Identification of Cholesterol Binding Sites in the Serotonin1A Receptor

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

Metabotropic glutamate receptors (mGluRs) are class C G protein-coupled receptors that function as obligate dimers in regulating neurotransmission and synaptic plasticity in the central nervous system. The mGluR1 subtype has been shown to be modulated by the membrane lipid environment, particularly cholesterol, though the molecular mechanisms remain elusive. In this study, we employed all-atom molecular dynamics simulations to investigate the effects of cholesterol on the conformational dynamics of the mGluR1 seven-transmembrane (7TM) domain in an inactive state model. Simulations were performed with three different cholesterol concentrations (0%, 10%, and 25%) in a palmitoyl-oleoyl phosphatidylcholine (POPC) lipid bilayer system. Our results demonstrate that cholesterol induces conformational changes in the mGluR1 dimer more significantly than in the individual protomers. Notably, cholesterol modulates the dynamics and conformations of the TM1 and TM2 helices at the dimer interface. Interestingly, an intermediate cholesterol concentration of 10% elicits more pronounced conformational changes compared to both cholesterol-depleted (0%) and cholesterol-enriched (25%) systems. Specific electrostatic interaction unique to the 10% cholesterol system further corroborate these conformational differences. Given the high sequence conservation of the 7TM domains across mGluR subtypes, the cholesterol-dependent effects observed in mGluR1 are likely applicable to other members of this receptor family. Our findings provide atomistic insights into how cholesterol modulates the conformational landscape of mGluRs, which could impact their function and signaling mechanisms.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

How Neuromembrane Lipids Modulate Membrane Proteins: Insights from G-Protein-Coupled Receptors (GPCRs) and Receptor Tyrosine Kinases (RTKs)

Lipids play a diverse and critical role in cellular processes in all tissues. The unique lipid composition of nerve membranes is particularly interesting because it contains, among other things, polyunsaturated lipids, such as docosahexaenoic acid, which the body only gets through the diet. The crucial role of lipids in neurological processes, especially in receptor-mediated cell signaling, is emphasized by the fact that in many neuropathological diseases there are significant deviations in the lipid composition of nerve membranes compared to healthy individuals. The lipid composition of neuromembranes can significantly affect the function of receptors by regulating the physical properties of the membrane or by affecting specific interactions between receptors and lipids. In addition, it is worth noting that the ligand-binding pocket of many receptors is located inside the cell membrane, due to which lipids can even modulate the binding of ligands to their receptors. These mechanisms highlight the importance of lipids in the regulation of membrane receptor activation and function. In this article, we focus on two major protein families: G-protein-coupled receptors (GPCRs) and receptor tyrosine kinases (RTKs) and discuss how lipids affect their function in neuronal membranes, elucidating the basic mechanisms underlying neuronal function and dysfunction.

Cholesterol in Class C GPCRs: Role, Relevance, and Localization

G-protein coupled receptors (GPCRs), one of the largest superfamilies of cell-surface receptors, are heptahelical integral membrane proteins that play critical roles in virtually every organ system. G-protein-coupled receptors operate in membranes rich in cholesterol, with an imbalance in cholesterol level within the vicinity of GPCR transmembrane domains affecting the structure and/or function of many GPCRs, a phenomenon that has been linked to several diseases. These effects of cholesterol could result in indirect changes by altering the mechanical properties of the lipid environment or direct changes by binding to specific sites on the protein. There are a number of studies and reviews on how cholesterol modulates class A GPCRs; however, this area of study is yet to be explored for class C GPCRs, which are characterized by a large extracellular region and often form constitutive dimers. This review highlights specific sites of interaction, functions, and structural dynamics involved in the cholesterol recognition of the class C GPCRs. We summarize recent data from some typical family members to explain the effects of membrane cholesterol on the structural features and functions of class C GPCRs and speculate on their corresponding therapeutic potential.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

PI(4,5)P2 and Cholesterol: Synthesis, Regulation, and Functions

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) is the most abundant membrane phosphoinositide and cholesterol is an essential component of the plasma membrane (PM). Both lipids play key roles in a variety of cellular functions including as signaling molecules and major regulators of protein function. This chapter provides an overview of these two important lipids. Starting from a brief description of their structure, synthesis, and regulation, the chapter continues to describe the primary functions and signaling processes in which PI(4,5)P₂ and cholesterol are involved. While PI(4,5)P₂ and cholesterol can act independently, they often act in concert or affect each other's impact. The chapters in this volume on "Cholesterol and PI(4,5)P₂ in Vital Biological Functions: From Coexistence to Crosstalk" focus on the emerging relationship between cholesterol and PI(4,5)P₂ in a variety of biological systems and processes. In this chapter, the next section provides examples from the ion channel field demonstrating that PI(4,5)P₂ and cholesterol can act via common mechanisms. The chapter ends with a discussion of future directions.

Synergistic and Competitive Lipid Interactions in the Serotonin1A Receptor Microenvironment

The interaction of lipids with G-protein-coupled receptors (GPCRs) has been shown to modulate and dictate several aspects of GPCR organization and function. Diverse lipid interaction sites have been identified from structural biology, bioinformatics, and molecular dynamics studies. For example, multiple cholesterol interaction sites have been identified in the serotonin_1A receptor, along with distinct and overlapping sphingolipid interaction sites. How these lipids interact with each other and what is the resultant effect on the receptor is still not clear. In this work, we have analyzed lipid-lipid crosstalk at the receptor of the serotonin_1A receptor embedded in a membrane bilayer that mimics the neuronal membrane composition by long coarse-grain simulations. Using a set of similarity coefficients, we classified lipids that bind at the receptor together as synergistic cobinding, and those that bind individually as competitive. Our results show that certain lipids interact with the serotonin_1A receptor in synergy with each other. Not surprisingly, the ganglioside GM1 and cholesterol show a synergistic cobinding, along with the relatively uncommon GM1-phosphatidylethanolamine (PE) and cholesterol-PE synergy. In contrast, certain lipid pairs such as cholesterol and sphingomyelin appear to be in competition at several sites, despite their coexistence in lipid nanodomains. In addition, we observed intralipid competition between two lipid tails, with the receptor exhibiting increased interactions with the unsaturated lipid tails. We believe our work represents an important step in understanding the diversity of GPCR-lipid interactions and exploring synergistic cobinding and competition in natural membranes.

Labelfree mapping and profiling of altered lipid homeostasis in the rat hippocampus after traumatic stress: Role of oxidative homeostasis

Oxidative and lipid homeostasis are altered by stress and trauma and post-traumatic stress disorder (PTSD) is associated with alterations to lipid species in plasma. Stress-induced alterations to lipid oxidative and homeostasis may exacerbate PTSD pathology, but few preclinical investigations of stress-induced lipidomic changes in the brain exist. Currently available techniques for the quantification of lipid species in biological samples require tissue extraction and are limited in their ability to retrieve spatial information. Raman imaging can overcome this limitation through the quantification of lipid species in situ in minimally processed tissue slices. Here, we utilized a predator exposure and psychosocial stress (PE/PSS) model of traumatic stress to standardize Raman imaging of lipid species in the hippocampus using LC-MS based lipidomics and these data were confirmed with qRT-PCR measures of mRNA expression of relevant enzymes and transporters. Electron Paramagnetic Resonance Spectroscopy (EPR) was used to measure free radical production and an MDA assay to measure oxidized polyunsaturated fatty acids. We observed that PE/PSS is associated with increased cholesterol, altered lipid concentrations, increased free radical production and reduced oxidized polyunsaturated fats (PUFAs) in the hippocampus (HPC), indicating shifts in lipid and oxidative homeostasis in the HPC after traumatic stress.

Structural dynamics of Smoothened (SMO) in the ciliary membrane and its interaction with membrane lipids

The Smoothened receptor (SMO, a 7 pass transmembrane domain, Class F GPCR family protein) plays a crucial role in the Hedgehog (HH) signaling pathway, which is involved in embryonic development and is implicated in various types of cancer throughout the animal kingdom. In the absence of HH signaling, SMO is inhibited by Patched 1 (PTC1; a 12 pass transmembrane domain protein), which is localized in the primary cilia. HH binding leads to the dislocation of PTC1 from the cilia, thus making way for SMO to localize in the primary cilia, as an essential prerequisite for its activation. We have carried out MARTINI coarse-grained molecular dynamics simulations of SMO in POPC and in ciliary membrane models, respectively, to study the interactions of SMO with cholesterol and other lipid molecules in the ciliary membrane, and to gain molecular-level insights into the role of the primary cilia in shaping the functional dynamics of SMO. We are able to identify the interaction of membrane cholesterols with definite sites and domains within SMO and relate them with known cholesterol-binding sequence and structure motifs. We show that cholesterol interactions with the transmembrane domain TMD, unlike those with the cysteine-rich domain (CRD) and the intracellular domain (ICD), are through residues belonging to known cholesterol-binding motifs. Notably, a few persistent interactions of cholesterol with lower TM cholesterol-binding domains are governed by the presence of multiple cholesterol-binding motifs. These analyses have also helped to identify and define a strict cholesterol consensus motif (CCM), which may well steer cholesterol into the hitherto identified binding sites within the TMD of SMO. We have also reported the interaction of phosphatidylinositol 4-phosphate with the intracellular region of transmembrane (TM) helices (TM1, TM3, TM4, and TM5), intracellular loop1, helix8, and Arg/Lys clusters of the ICD. Structural analysis of SMO domains shows significant changes in the CRD and ICD, during the course of the simulation. Further detailed analysis of the dynamics of the TMD reveals the movements of TM5, TM6, and TM7, linked with the helix8, which are possibly involved in shaping the conformational disposition of the ICD. The movement of these TM helices could possibly be a consequence of interactions involving the extracellular domain and extracellular loops. In addition, our analysis also shows that phosphatidylinositol-4-phosphate (PI4P), along with some ICD cholesterols, are implicated in anchoring SMO in the membrane.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

Impact of lipid rafts on transient receptor potential channel activities

Members of the transient receptor potential (TRP) superfamily are cation channels that are expressed in nearly every mammalian cell type and respond as cellular sensors to various environmental stimuli. Light, pressure, osmolarity, temperature, and other stimuli can induce TRP calcium conductivity and correspondingly trigger many signaling processes in cells. Disruption of TRP channel activity, as a rule, harms cellular function. Despite numerous studies, the mechanisms of TRP channel regulation are not yet sufficiently clear, in part, because TRP channels are regulated by a broad set of ligands having diverse physical and chemical features. It is now known that some TRP members are located in membrane microdomains termed lipid rafts. Moreover, interaction between specific raft-associated lipids with channels may be a key regulation mechanism. This review examines recent findings related to the roles of lipid rafts in regulation of TRP channel activity. The mechanistic events of channel interactions with the main lipid raft constituent, cholesterol, are being clarified. Better understanding of mechanisms behind such interactions would help establish the key elements of TRP channel regulation and hence allow control of cellular responses to environmental stimuli.

Lipid–Protein Interactions in Plasma Membrane Organization and Function

Lipid–protein interactions in cells are involved in various biological processes, including metabolism, trafficking, signaling, host–pathogen interactions, and transmembrane transport. At the plasma membrane, lipid–protein interactions play major roles in membrane organization and function. Several membrane proteins have motifs for specific lipid binding, which modulate protein conformation and consequent function. In addition to such specific lipid–protein interactions, protein function can be regulated by the dynamic, collective behavior of lipids in membranes. Emerging analytical, biochemical, and computational technologies allow us to study the influence of specific lipid–protein interactions, as well as the collective behavior of membranes on protein function. In this article, we review the recent literature on lipid–protein interactions with a specific focus on the current state-of-the-art technologies that enable novel insights into these interactions.

Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Caveolin induced membrane curvature and lipid clustering: two sides of the same coin?

Caveolin-1 (cav-1) is a multi-domain membrane protein that is a key player in cell signaling, endocytosis and mechanoprotection. It is the principle component of cholesterol-rich caveolar domains and has been reported to induce membrane curvature. The molecular mechanisms underlying the interactions of cav-1 with complex membranes, leading to modulation of membrane topology and the formation of cholesterol-rich domains, remain elusive. In this study, we aim to understand the effect of lipid composition by analyzing the interactions of cav-1 with complex membrane bilayers comprised of about sixty lipid types. We have performed a series of coarse-grain molecular dynamics simulations using the Martini force-field with a cav-1 protein construct (residue 82-136) that includes the membrane binding domains and a palmitoyl tail. We observe that cav-1 induces curvature in this complex membrane, though it is restricted to a nanometer length scale. Concurrently, we observe a clustering of cholesterol, sphingolipids and other lipid molecules leading to the formation of nanodomains. Direct microsecond timescale interactions are observed for specific lipids such as cholesterol, phosphatidylserine and phosphatidylethanolamine lipid types. The results indicate that there is an interplay between membrane topology and lipid species. Our work is a step toward understanding how lipid composition and organization regulate the formation of caveolae, in the context of endocytosis and cell signaling.

Cholesterol in GPCR Structures: Prevalence and Relevance

Bound cholesterol molecules are emerging as important hallmarks of GPCR structures. In this commentary, we analyze their statistical prevalence and biological relevance.

Cholesterol footprint in high-resolution structures of serotonin receptors: Where are we now and what does it mean?

An emerging feature of several high-resolution GPCR structures is the presence of closely bound cholesterol molecules. In this Perspective, we share the excitement of the recent advancements in GPCR structural biology. We further highlight our laboratory's journey in comprehensively elucidating functional sensitivity of GPCRs (using the serotonin_1A receptor as a representative neurotransmitter GPCR) to membrane cholesterol and validation using a variety of assays and molecular dynamics simulations. Although high-resolution structures of many GPCRs have been reported in the last few years, the structure of the serotoin_1A receptor proved to be elusive for a long time. Very recently the cryo-EM structure of the serotoin_1A receptor displaying 10 bound cholesterol molecules has been reported. We conclude by providing a critical analysis of caveats involved in GPCR structure determination.

In Silico Identification of Cholesterol Binding Motifs in the Chemokine Receptor CCR3

CC motif chemokine receptor 3 (CCR3) is a Class A G protein-coupled receptor (GPCR) mainly responsible for the cellular trafficking of eosinophils. As such, it plays key roles in inflammatory conditions, such as asthma and arthritis, and the metastasis of many deadly forms of cancer. However, little is known about how CCR3 functionally interacts with its bilayer environment. Here, we investigate cholesterol binding sites in silico through Coarse-Grained Molecular Dynamics (MD) and Pylipid analysis using an extensively validated homology model based on the crystal structure of CCR5. These simulations identified several cholesterol binding sites containing Cholesterol Recognition/Interaction Amino Acid Consensus motif (CRAC) and its inversion CARC motifs in CCR3. One such site, a CARC site in TM1, in conjunction with aliphatic residues in TM7, emerged as a candidate for future investigation based on the cholesterol residency time within the binding pocket. This site forms the core of a cholesterol binding site previously observed in computational studies of CCR2 and CCR5. Most importantly, these cholesterol binding sites are conserved in other chemokine receptors and may provide clues to cholesterol regulation mechanisms in this subfamily of Class A GPCRs.

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

Mechanistic Insights into Pore Formation by an α-Pore Forming Toxin: Protein and Lipid Bilayer Interactions of Cytolysin A

Pore forming toxins (PFTs) are the largest class of bacterial toxins playing a central role in bacterial pathogenesis. They are proteins specifically designed to form nanochannels in the membranes of target cells, ultimately resulting in cell death and establishing infection. PFTs are broadly classified as α- and β-PFTs, depending on secondary structures that form the transmembrane channel. A unique feature about this class of proteins is the drastic conformational changes and complex oligomerization pathways that occur upon exposure to the plasma membrane. A molecular understanding of pore formation has implications in designing novel intervention strategies to combat rising antimicrobial resistance, targeted-cancer therapy, as well as designing nanopores for specialized technologies. Central to unraveling the pore formation pathway is the availability of high resolution crystal structures. In this regard, β-toxins are better understood, when compared with α-toxins whose pore forming mechanisms are complicated by an incomplete knowledge of the driving forces for amphiphatic membrane-inserted helices to organize into functional pores. With the publication of the first crystal structure for an α-toxin, cytolysin A (ClyA), in 2009 we embarked on an extensive multiscale study to unravel its pore forming mechanism. This Account represents the collective mechanistic knowledge gained in our laboratories using a variety of experimental and theoretical techniques which include large scale molecular dynamics (MD) simulations, kinetic modeling studies, single-molecule fluorescence imaging, and super-resolution spectroscopy. We reported MD simulations of the ClyA protomer, oligomeric intermediates, and full pore complex in a lipid bilayer and mapped the conformational transitions that accompany membrane binding. Using single-molecule fluorescence imaging, the conformational transition was experimentally verified by analysis of various diffusion states of membrane bound ClyA. Importantly, we have uncovered a hitherto unknown putative cholesterol binding motif in the membrane-inserted helix of ClyA. Distinct binding pockets for cholesterol formed by adjacent membrane-inserted helices are revealed in MD simulations. Cholesterol appears to play a dual role by stabilizing both the membrane-inserted protomer as well as oligomeric intermediates. Molecular dynamics simulations and kinetic modeling studies suggest that the membrane-inserted arcs oligomerize reversibly to form the predominant transmembrane oligomeric intermediates during pore formation. We posit that this mechanistic understanding of the complex action of α-PFTs has implications in unraveling pore assembly across the wider family of bacterial toxins. With emerging antimicrobial resistance, alternate therapies may rely on disrupting pore functionality or oligomerization of these pathogenic determinants utilized by bacteria, and our study includes assessing the potential for dendrimers as pore blockers.

Interplay of cholesterol, membrane bilayers and the AT1R: A cholesterol consensus motif on AT1R is revealed

Hypertension, mediated by the Angiotensin II receptor type 1 (AT1R), is still the major cause of premature death despite the discovery of novel therapeutics, highlighting the importance of an in depth understanding of the drug-AT1R recognition mechanisms coupled with the impact of the membrane environment on the interaction of drugs with AT1R. Herein, we examine the interplay of cholesterol-lipid-candesartan and the AT1R using Molecular Dynamics simulations of a model membrane consisting of 60:40 mol%. DPPC:cholesterol, candesartan and the AT1R, mimicking the physiological cholesterol concentration in sarcolemma membranes. The simulations of the model membrane of 60:40 mol%. DPPC:cholesterol were further validated using DOSY NMR experiments. Interestingly, our results suggest a significant role of cholesterol in the AT1R function imposed through a Cholesterol Consensus Motif (CCM) in the receptor, which could be crucial in the drug binding process. Candesartan diffusion towards AT1R through incorporation into lipid bilayers, appears to be retarded by the presence of cholesterol. However, its direct approach towards AT1R may be facilitated through the mobility induced on the N-terminus by the cholesterol binding on the CCM these novel insights could pave the way towards the development of more potent pharmaceutical agents to combat hypertension more effectively.

Structure, dynamics and lipid interactions of serotonin receptors: excitements and challenges

Serotonin (5-hydroxytryptamine, 5-HT) is an intrinsically fluorescent neurotransmitter found in organisms spanning a wide evolutionary range. Serotonin exerts its diverse actions by binding to distinct cell membrane receptors which are classified into many groups. Serotonin receptors are involved in regulating a diverse array of physiological signaling pathways and belong to the family of either G protein-coupled receptors (GPCRs) or ligand-gated ion channels. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions such as sleep, mood, pain, anxiety, depression, aggression, and learning. Serotonin receptors act as drug targets for a number of diseases, particularly neuropsychiatric disorders. The signaling mechanism and efficiency of serotonin receptors depend on their amazing ability to rapidly access multiple conformational states. This conformational plasticity, necessary for the wide variety of functions displayed by serotonin receptors, is regulated by binding to various ligands. In this review, we provide a succinct overview of recent developments in generating and analyzing high-resolution structures of serotonin receptors obtained using crystallography and cryo-electron microscopy. Capturing structures of distinct conformational states is crucial for understanding the mechanism of action of these receptors, which could provide important insight for rational drug design targeting serotonin receptors. We further provide emerging information and insight from studies on interactions of membrane lipids (such as cholesterol) with serotonin receptors. We envision that a judicious combination of analysis of high-resolution structures and receptor-lipid interaction would allow a comprehensive understanding of GPCR structure, function and dynamics, thereby leading to efficient drug discovery.

Molecular evolution of a collage of cholesterol interaction motifs in transmembrane helix V of the serotonin1A receptor

The human serotonin_1A receptor is a representative member of the superfamily of G protein-coupled receptors (GPCRs) and an important drug target for neurological disorders. Using a combination of biochemical, biophysical and molecular dynamics simulation approaches, we and others have shown that membrane cholesterol modulates the organization, dynamics and function of vertebrate serotonin_1A receptors. Previous studies have shown that the cytoplasmic portion of transmembrane helix V (TM V) and the extramembraneous intracellular loop 3 are critical for G-protein coupling, phosphorylation and desensitization of the receptor. We have recently resolved a collage of putative cholesterol interaction motifs from the amino acid sequence overlapping this region. In this paper, we explore the sequence plasticity of this fragment that may have adapted to altered membrane lipidome, after vertebrates evolved from primordial invertebrates. Since invertebrates have lower levels of membrane cholesterol relative to vertebrates, we compared TM V sequence fragments from invertebrate serotonin₁ receptors with vertebrate orthologs to infer the sequence plasticity in TM V. We report that the average number of cholesterol interaction motifs in TM V for diverse phyla represents an increasing trend that could mirror vertebrate evolution from primordial invertebrates. By statistical modeling, we propose that the collage of cholesterol interaction motifs in TM V of the human serotonin_1A receptor may have evolved from rudimentary collages, reminiscent of primordial invertebrate orthologs. Taken together, we propose that a repertoire of cholesterol-philic nonsynonymous substitutions may have enhanced collage complexity in TM V during vertebrate evolution.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.

Mechanisms of Action of Cassiae Semen for Weight Management: A Computational Molecular Docking Study of Serotonin Receptor 5-HT2C

Overweight and obesity is a growing global health concern. Current management of obesity includes lifestyle intervention, bariatric surgery and medication. The serotonin receptor, 5-HT2C, is known to mediate satiety, appetite and consumption behaviour. Lorcaserin, an appetite control drug, has demonstrated efficacy in appetite control by targeting 5-HT2C but causes undesirable side effects. This study aimed to explore the potential usage of Cassiae semen (CS), a well-known traditional Chinese medicine used to treat obesity. A computational molecular docking study was performed to determine the binding mechanism of CS compounds to the 5-HT2C receptors in both active, agonist-bound and inactive, antagonist-bound conformations. By comparing binding poses and predicted relative binding affinities towards the active or inactive forms of the receptor, we hypothesise that two of the CS compounds studied may be potent agonists which may mimic the appetite suppression effects of lorcaserin: obtusifoliol and cassiaside B2. Furthermore, two ligands, beta-sitosterol and juglanin, were predicted to bind favourably to 5-HT2C outside of the known agonist binding pocket in the active receptor, suggesting that such ligands may serve as positive allosteric modulators of 5-HT2C receptor function. Overall, this study proposed several CS compounds which may be responsible for exerting anti-obesity effects via appetite suppression by 5-HT2C receptor activation.

Cholesterol interaction motifs in G protein-coupled receptors: Slippery hot spots?

G protein-coupled receptors (GPCRs) are cell membrane associated signaling hubs that orchestrate a multitude of cellular functions upon binding to a diverse variety of extracellular ligands. Since GPCRs are integral membrane proteins with seven-transmembrane domain architecture, their function, organization and dynamics are intimately regulated by membrane lipids, such as cholesterol. Cholesterol is an extensively studied lipids in terms of its effects on GPCR structure and function. One of the possible mechanisms underlying modulation of GPCR function by cholesterol is via specific interaction of GPCRs with membrane cholesterol. These interactions of GPCRs with membrane cholesterol are often attributed to structural features of GPCRs that could facilitate their preferential association with cholesterol. In this backdrop, cholesterol interaction motifs represent putative interaction sites on GPCRs that could facilitate cholesterol-sensitive function of these receptors. In this review, we provide an overview of cholesterol interaction motifs found in GPCRs, which have been identified through a combination of crystallography, bioinformatics analysis, and functional studies. In addition, we will highlight, using specific examples, why mere presence of a cholesterol interaction motif at a given site may not directly implicate its role in interaction with membrane cholesterol. We therefore believe that experimental approaches, followed by functional analysis of cholesterol sensitivity of GPCRs, would provide a better understanding of the role played by these motifs in cholesterol-sensitive function. We envision that a comprehensive knowledge of cholesterol interaction sites in GPCRs would allow us to develop a better understanding of GPCR structure-function paradigm, and could be useful in future therapeutics. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Analytical and Computational Methods > Computational Methods Laboratory Methods and Technologies > Macromolecular Interactions, Methods.

Role of cholesterol-mediated effects in GPCR heterodimers

G protein-coupled receptors (GPCRs) are transmembrane receptors that mediate a large number of cellular responses. The organization of GPCRs into dimers and higher-order oligomers is known to allow a larger repertoire of downstream signaling events. In this context, a crosstalk between the adenosine and dopamine receptors has been reported, indicating the presence of heterodimers that are functionally relevant. In this paper, we explored the effect of membrane cholesterol on the adenosine_2A (A_2A) and dopamine D₃ (D₃) receptors using coarse-grain molecular dynamics simulations. We analyzed cholesterol interaction sites on the A_2A receptor and were able to reproduce the sites indicated by crystallography and previous atomistic simulations. We predict novel cholesterol interaction sites on the D₃ receptor that could be important in the reported cholesterol sensitivity in receptor function. Further, we analyzed the formation of heterodimers between the two receptors. Our results suggest that membrane cholesterol modulates the relative population of several co-existing heterodimer conformations. Both direct receptor-cholesterol interaction and indirect membrane effects contribute toward the modulation of heterodimer conformations. These results constitute one of the first examples of modulation of GPCR hetero-dimerization by membrane cholesterol, and could prove to be useful in designing better therapeutic strategies.

Molecular Dynamics Simulations of Adenosine Receptors: Advances, Applications and Trends

:

Adenosine receptors (ARs) are transmembrane proteins that belong to the G protein-coupled receptors (GPCRs) superfamily and mediate the biological functions of adenosine. To date, four AR subtypes are known, namely A1, A2A, A2B and A3 that exhibit different signaling pathways, tissue localization, and mechanisms of activation. Moreover, the widespread ARs and their implication in numerous physiological and pathophysiological conditions had made them pivotal therapeutic targets for developing clinically effective agents.

:

The crystallographic success in identifying the 3D crystal structures of A2A and A1 ARs has dramatically enriched our understanding of their structural and functional properties such as ligand binding and signal transduction. This, in turn, has provided a structural basis for a larger contribution of computational methods, particularly molecular dynamics (MD) simulations, toward further investigation of their molecular properties and designing bioactive ligands with therapeutic potential. MD simulation has been proved to be an invaluable tool in investigating ARs and providing answers to some critical questions. For example, MD has been applied in studying ARs in terms of ligand-receptor interactions, molecular recognition, allosteric modulations, dimerization, and mechanisms of activation, collectively aiding in the design of subtype selective ligands.

:

In this review, we focused on the advances and different applications of MD simulations utilized to study the structural and functional aspects of ARs that can foster the structure-based design of drug candidates. In addition, relevant literature was briefly discussed which establishes a starting point for future advances in the field of drug discovery to this pivotal group of drug targets.

Resolving the conformational dynamics of ErbB growth factor receptor dimers

The combinatorial dimerization of the ErbB growth factor receptors (ErbB¹- ErbB⁴) are critical for their function. Here, we have characterized the conformational dynamics of ErbB transmembrane homo-dimers and hetero-dimers by using a coarse-grain simulation framework. All dimers, except ErbB^4-4 and ErbB^1-4, exhibit at least two conformations. The reported NMR structures correspond to one of these conformations, representing the N-terminal active state in ErbB^1-1 (RH2), ErbB^2-2 (RH1) and ErbB^4-4 (RH) homo-dimers and the LH dimer in ErbB^3-3 homo-dimer, validating the computational approach. Further, we predict a right-handed ErbB^3-3 dimer conformer that warrants experimental testing. The five hetero-dimers that have not yet been experimentally resolved display prominent right-handed dimers associating by the SmXXXSm motif. Our results provide insights into the constitutive signaling of ErbB⁴ after cleavage of the extracellular region. The presence of the inactive-like dimer conformers leading to symmetric kinase domains gives clues on the autoinhibition of the receptor dimers. The dimer states characterized here represent an important step towards understanding the combinatorial cross associations in the ErbB family.

Changes in the function of angiotensin II type 1 receptor due to cholesterol depletion from cell membrane

Blockers of G-protein coupled receptors (GPCRs), angiotensin II (Ang II) type 1 (AT₁) receptor and β₁-adrenergic (Ad) receptor, have been shown to improve the prognosis of cardiovascular disease. Cholesterol molecules in the cell membrane are needed to stabilize GPCRs as well as the cell membrane itself. We determined whether the functions of AT₁ and β₁-Ad receptors were changed by cholesterol depletion from cardiovascular cell membranes. Ang II-induced inositol phosphate production through AT₁ receptor was suppressed by cholesterol depletion from cell membranes using rosuvastatin or methyl-β-cyclodextrin (MβCD), whereas isoproterenol-induced cyclic AMP production through β₁-Ad receptor did not change after cholesterol depletion. In addition, the binding affinities of Ang II and AT₁ receptor blocker after cholesterol depletion were significantly lower than those before depletion. Although AT₁ receptor expression levels did not change after cholesterol depletion, the expression levels of AT₁ receptor that could bind to Ang II significantly decreased after depletion. The changes in the structure of AT₁ receptor due to depletion were confirmed by substituted-cysteine accessibility mapping. In conclusion, Ang II-induced activation of AT₁ receptor is reduced without affecting the function of β₁-Ad receptor after cholesterol depletion from cardiovascular cell membranes.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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.

Interaction of POPC, DPPC, and POPE with the μ opioid receptor: A coarse-grained molecular dynamics study

The μ opioid receptor (μOR), which is part of the G protein-coupled receptors family, is a membrane protein that is modulated by its lipid environment. In the present work, we model μOR in three different membrane systems: POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), and DPPC (1, 2-dipalmitoyl-sn-glycero-3-phosphocholine) through 45 μs molecular dynamics (MD) simulations at the coarse-grained level. Our theoretical studies provide new insights to the lipid-induced modulation of the receptor. Particularly, to characterize how μOR interacts with each lipid, we analyze the tilt of the protein, the number of contacts occurring between the lipids and each amino acid of the receptor, and the μOR-lipid interface described as a network graph. We also analyze the variations in the number and the nature of the protein contacts that are induced by the lipid structure. We show that POPC interacts preferentially with helix 1 (H1) and helices H5-H6, POPE, with H5-H6 and H6-H7, and DPPC, with H4 and H6. We demonstrate how each of the three lipids shape the structure of the μOR.

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.

Influence of Cholesterol and Its Stereoisomers on Members of the Serotonin Receptor Family

Despite the ubiquity of cholesterol within the cell membrane, the mechanism by which it influences embedded proteins remains elusive. Numerous G-protein coupled receptors exhibit dramatic responses to membrane cholesterol with regard to the ligand-binding affinity and functional properties, including the 5-HT receptor family. Here, we use over 25 μs of unbiased atomistic molecular dynamics simulations to identify cholesterol interaction sites in the 5-HT_1B and 5-HT_2B receptors and evaluate their impact on receptor structure. Susceptibility to membrane cholesterol is shown to be subtype dependent and determined by the quality of interactions between the extracellular loops. Charged residues are essential for maintaining the arrangement of the extracellular surface in 5-HT_2B; in the absence of such interactions, the extracellular surface of the 5-HT_1B is malleable, populating a number of distinct conformations. Elevated cholesterol density near transmembrane helix 4 is considered to be conducive to the conformation of extracellular loop 2. Occupation of this site is also shown to be stereospecific, illustrated by differential behavior of nat-cholesterol isomers, ent- and epi-cholesterol. In simulations containing the endogenous agonist, serotonin, cholesterol binding at transmembrane helix 4 biases bound serotonin molecules toward an unexpected binding mode in the extended binding pocket. The results highlight the capability of membrane cholesterol to influence the mobility of the extracellular surface in the 5-HT₁ receptor family and manipulate the architecture of the extracellular ligand-binding pocket.

A collage of cholesterol interaction motifs in the serotonin1A receptor: An evolutionary implication for differential cholesterol interaction

The serotonin_1A receptor is a representative member of the G protein-coupled receptor (GPCR) superfamily and acts as an important drug target. In our previous work, we comprehensively demonstrated that membrane cholesterol is necessary in the organization, dynamics and function of the serotonin_1A receptor. In this context, analysis of high-resolution GPCR crystal structures in general and in silico studies of the serotonin_1A receptor in particular, have suggested the presence of cholesterol interaction sites (hotspots) in various regions of the receptor. In this work, we have identified an evolutionarily conserved collage of four categories of cholesterol interaction motifs associated with transmembrane helix V and the adjacent intracellular loop 3 fragment of the vertebrate serotonin_1A receptor. This collage of motifs represents a total of twenty diverse context-dependent cholesterol interaction configurations. We envision that the gamut of cholesterol interaction sites, characterized by sequence plasticity in cholesterol interaction, could be relevant in receptor-cholesterol interaction in membranes of varying cholesterol content and organization, as found in diverse cell types. We conclude that an evolutionarily conserved mechanism of GPCR-cholesterol interaction allows the serotonin_1A receptor to adapt to diverse membrane cholesterol levels during natural evolution.

Taste Receptors: New Players in Sperm Biology

Taste receptors were first described as sensory receptors located on the tongue, where they are expressed in small clusters of specialized epithelial cells. However, more studies were published in recent years pointing to an expression of these proteins not only in the oral cavity but throughout the body and thus to a physiological role beyond the tongue. The recent observation that taste receptors and components of the coupled taste transduction cascade are also expressed during the different phases of spermatogenesis as well as in mature spermatozoa from mouse to humans and the overlap between the ligand spectrum of taste receptors with compounds in the male and female reproductive organs makes it reasonable to assume that sperm "taste" these different cues in their natural microenvironments. This assumption is assisted by the recent observations of a reproductive phenotype of different mouse lines carrying a targeted deletion of a taste receptor gene as well as the finding of a significant correlation between human male infertility and some polymorphisms in taste receptors genes. In this review, we depict recent findings on the role of taste receptors in male fertility, especially focusing on their possible involvement in mechanisms underlying spermatogenesis and post testicular sperm maturation. We also highlight the impact of genetic deletions of taste receptors, as well as their polymorphisms on male reproduction.

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.

Role of membrane environment and membrane-spanning protein regions in assembly and function of the Class II Major Histocompatibility complex

Class II Major Histocompatibility complex (MHC-II) is a polymorphic heterodimer that binds antigen-derived peptides and presents them on the surface of antigen presenting cells. This mechanism of antigen presentation leads to recognition by CD4 T-cells and T-cell activation, making it a critical element of adaptive immune response. For this reason, the structural determinants of MHC-II function have been of great interest for the past 30 years, resulting in a robust structural understanding of the extracellular regions of the complex. However, the membrane-localized regions have also been strongly implicated in protein-protein and protein-lipid interactions that facilitate Class II assembly, transport and function, and it is these regions that are the focus of this review. Here we describe studies that reveal the strong and selective interactions between the transmembrane domains of the MHC α, and invariant chains which, when altered, have broad reaching impacts on antigen presentation and Class II function. We also summarize work that clearly demonstrates the link between membrane lipid composition (particularly the presence of cholesterol) and MHC-II conformation, subsequent peptide binding, and downstream T-cell activation. We have integrated these studies into a comprehensive view of Class II transmembrane domain biology.

A Critical Analysis of Molecular Mechanisms Underlying Membrane Cholesterol Sensitivity of GPCRs

G protein-coupled receptors (GPCRs) are the largest and a diverse family of proteins involved in signal transduction across biological membranes. GPCRs mediate a wide range of physiological processes and have emerged as major targets for the development of novel drug candidates in all clinical areas. Since GPCRs are integral membrane proteins, regulation of their organization, dynamics, and function by membrane lipids, in particular membrane cholesterol, has emerged as an exciting area of research. Cholesterol sensitivity of GPCRs could be due to direct interaction of cholesterol with the receptor (specific effect). Alternately, GPCR function could be influenced by the effect of cholesterol on membrane physical properties (general effect). In this review, we critically analyze the specific and general mechanisms of the modulation of GPCR function by membrane cholesterol, taking examples from representative GPCRs. While evidence for both the proposed mechanisms exists, there appears to be no clear-cut distinction between these two mechanisms, and a combination of these mechanisms cannot be ruled out in many cases. We conclude that classifying the mechanism underlying cholesterol sensitivity of GPCR function merely into these two mutually exclusive classes could be somewhat arbitrary. A more holistic approach could be suitable for analyzing GPCR-cholesterol interaction.

Interplay between Membrane Curvature and Cholesterol: Role of Palmitoylated Caveolin-1

Caveolin-1 (cav-1) is an important player in cell signaling and endocytosis that has been shown to colocalize with cholesterol-rich membrane domains. Experimental studies with varying cav-1 constructs have suggested that it can induce both cholesterol clustering and membrane curvature. Here, we probe the molecular origin of membrane curvature and cholesterol clustering by cav-1 by using coarse-grain molecular dynamics simulations. We have performed a series of simulations of a functionally important cav-1 construct, comprising the membrane-interacting domains and a C-terminal palmitoyl tail. Our results suggest that cav-1 is able to induce cholesterol clustering in the membrane leaflet to which it is bound as well as the opposing leaflet. A positive membrane curvature is observed upon cav-1 binding in cholesterol-containing bilayers. Interestingly, we observe an interplay between cholesterol clustering and membrane curvature such that cav-1 is able to induce higher membrane curvature in cholesterol-rich membranes. The role of the cav-1 palmitoyl tail is less clear and appears to increase the membrane contacts. Further, we address the importance of the secondary structure of cav-1 domains and show that it could play an important role in membrane curvature and cholesterol clustering. Our work is an important step toward a molecular picture of caveolae and vesicular endocytosis.

Lipid composition and salt concentration as regulatory factors of the anion selectivity of VDAC studied by coarse-grained molecular dynamics simulations

The voltage-dependent anion channel (VDAC) is a mitochondrial outer membrane protein whose fundamental function is to facilitate and regulate the flow of metabolites between the cytosol and the mitochondrial intermembrane space. Using coarse-grained molecular dynamics simulations, we investigated the dependence of VDAC selectivity towards small inorganic anions on two factors: the ionic strength and the lipid composition. In agreement with experimental data we found that VDAC becomes less anion selective with increasing salt concentration due to the screening of a few basic residues that point into the pore lumen. The molecular dynamics simulations provide insight into the regulation mechanism of VDAC selectivity by the composition in the lipid membrane and suggest that the ion distribution is differently modulated by POPE compared to the POPC bilayer. This occurs through the more persistent interactions of acidic residues located at both rims of the β-barrel with head groups of POPE which in turn impact the electrostatic potential and thereby the selectivity of the pore. This mechanism occurs not only in POPE single component membranes but also in a mixed POPE/POPC bilayer by an enrichment of POPE over POPC lipids on the surface of VDAC. Thus we show here that computationally-inexpensive coarse-grained simulations are able to capture, in a semi-quantitative way, essential features of VDAC anion selectivity and could pave the way toward a molecular level understanding of metabolite transport in natural membranes.

A Streamlined, General Approach for Computing Ligand Binding Free Energies and Its Application to GPCR-Bound Cholesterol

The theory of receptor-ligand binding equilibria has long been well-established in biochemistry, and was primarily constructed to describe dilute aqueous solutions. Accordingly, few computational approaches have been developed for making quantitative predictions of binding probabilities in environments other than dilute isotropic solution. Existing techniques, ranging from simple automated docking procedures to sophisticated thermodynamics-based methods, have been developed with soluble proteins in mind. Biologically and pharmacologically relevant protein-ligand interactions often occur in complex environments, including lamellar phases like membranes and crowded, nondilute solutions. Here, we revisit the theoretical bases of ligand binding equilibria, avoiding overly specific assumptions that are nearly always made when describing receptor-ligand binding. Building on this formalism, we extend the asymptotically exact Alchemical Free Energy Perturbation technique to quantifying occupancies of sites on proteins in a complex bulk, including phase-separated, anisotropic, or nondilute solutions, using a thermodynamically consistent and easily generalized approach that resolves several ambiguities of current frameworks. To incorporate the complex bulk without overcomplicating the overall thermodynamic cycle, we simplify the common approach for ligand restraints by using a single distance-from-bound-configuration (DBC) ligand restraint during AFEP decoupling from protein. DBC restraints should be generalizable to binding modes of most small molecules, even those with strong orientational dependence. We apply this approach to compute the likelihood that membrane cholesterol binds to known crystallographic sites on three GPCRs (β₂-adrenergic, 5HT-2B, and μ-opioid) at a range of concentrations. Nonideality of cholesterol in a binary cholesterol:phosphatidylcholine (POPC) bilayer is characterized and consistently incorporated into the interpretation. We find that the three sites exhibit very different affinities for cholesterol: The site on the adrenergic receptor is predicted to be high affinity, with 50% occupancy for 1:10⁹ CHOL:POPC mixtures. The sites on the 5HT-2B and μ-opioid receptor are predicted to be lower affinity, with 50% occupancy for 1:10³ CHOL:POPC and 1:10² CHOL:POPC, respectively. These results could not have been predicted from the crystal structures alone.

Molecular Dynamics Simulations of Kir2.2 Interactions with an Ensemble of Cholesterol Molecules

Cholesterol is a major regulator of multiple types of ion channels, but the specific mechanisms and the dynamics of its interactions with the channels are not well understood. Kir2 channels were shown to be sensitive to cholesterol through direct interactions with "cholesterol-sensitive" regions on the channel protein. In this work, we used Martini coarse-grained simulations to analyze the long (μs) timescale dynamics of cholesterol with Kir2.2 channels embedded into a model membrane containing POPC phospholipid with 30 mol% cholesterol. This approach allows us to simulate the dynamic, unbiased migration of cholesterol molecules from the lipid membrane environment to the protein surface of Kir2.2 and explore the favorability of cholesterol interactions at both surface sites and recessed pockets of the channel. We found that the cholesterol environment surrounding Kir channels forms a complex milieu of different short- and long-term interactions, with multiple cholesterol molecules concurrently interacting with the channel. Furthermore, utilizing principles from network theory, we identified four discrete cholesterol-binding sites within the previously identified cholesterol-sensitive region that exist depending on the conformational state of the channel-open or closed. We also discovered that a twofold decrease in the cholesterol level of the membrane, which we found earlier to increase Kir2 activity, results in a site-specific decrease of cholesterol occupancy at these sites in both the open and closed states: cholesterol molecules at the deepest of these discrete sites shows no change in occupancy at different cholesterol levels, whereas the remaining sites showed a marked decrease in occupancy.

Cholesterol promotes Cytolysin A activity by stabilizing the intermediates during pore formation

Pore-forming toxins (PFTs) form nanoscale pores across target membranes causing cell death. Cytolysin A (ClyA) from Escherichia coli is a prototypical α-helical toxin that contributes to cytolytic phenotype of several pathogenic strains. It is produced as a monomer and, upon membrane exposure, undergoes conformational changes and finally oligomerizes to form a dodecameric pore, thereby causing ion imbalance and finally cell death. However, our current understanding of this assembly process is limited to studies in detergents, which do not capture the physicochemical properties of biological membranes. Here, using single-molecule imaging and molecular dynamics simulations, we study the ClyA assembly pathway on phospholipid bilayers. We report that cholesterol stimulates pore formation, not by enhancing initial ClyA binding to the membrane but by selectively stabilizing a protomer-like conformation. This was mediated by specific interactions by cholesterol-interacting residues in the N-terminal helix. Additionally, cholesterol stabilized the oligomeric structure using bridging interactions in the protomer-protomer interfaces, thereby resulting in enhanced ClyA oligomerization. This dual stabilization of distinct intermediates by cholesterol suggests a possible molecular mechanism by which ClyA achieves selective membrane rupture of eukaryotic cell membranes. Topological similarity to eukaryotic membrane proteins suggests evolution of a bacterial α-toxin to adopt eukaryotic motifs for its activation. Broad mechanistic correspondence between pore-forming toxins hints at a wider prevalence of similar protein membrane insertion mechanisms.

Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins

Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.