Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations

Aberrant lipid accumulation and retinal pigmental epithelium dysfunction in PRCD-deficient mice

Progressive Rod-Cone Degeneration (PRCD) is an integral membrane protein found in photoreceptor outer segment (OS) disc membranes and its function remains unknown. Mutations in Prcd are implicated in Retinitis pigmentosa (RP) in humans and multiple dog breeds. PRCD-deficient models exhibit decreased levels of cholesterol in the plasma. However, potential changes in the retinal cholesterol remain unexplored. In addition, impaired phagocytosis observed in these animal models points to potential deficits in the retinal pigment epithelium (RPE). Here, using a Prcd -/- murine model we investigated the alterations in the retinal cholesterol levels and impairments in the structural and functional integrity of the RPE. Lipidomic and immunohistochemical analyses show a 5-fold increase in the levels of cholesteryl esters (C.Es) and accumulation of neutral lipids in the PRCD-deficient retina, respectively, indicating alterations in total retinal cholesterol. Longitudinal fundus and spectral domain optical coherence tomography (SD-OCT) examinations showed focal lesions and RPE hyperreflectivity. Strikingly, the RPE of Prcd -/- mice exhibited age-related pathological features such as neutral lipid deposits, lipofuscin accumulation, Bruch's membrane (BrM) thickening and drusenoid focal deposits, mirroring an Age-related Macular Degeneration (AMD)-like phenotype. We propose that the extensive lipofuscin accumulation likely impairs lysosomal function, leading to the defective phagocytosis observed in Prcd -/- mice. Our findings support the dysregulation of retinal cholesterol homeostasis in the absence of PRCD. Further, we demonstrate that progressive photoreceptor degeneration in Prcd -/- mice is accompanied by progressive structural and functional deficits in the RPE, which likely exacerbates vision loss over time.

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.

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.

Insights of ligand binding in modeled h5-HT1A receptor: homology modeling, docking, MM-GBSA, screening and molecular dynamics

The pharmacologically characterized receptor subtype of the serotonin family, the 5HT1A receptor is implicated in the pathophysiology and treatment of depression and anxiety-related disorders. Being the most extensively targeted receptor for developing novel antidepressants and anxiolytics, a near-ideal theoretical model can aid in high-throughput screening of promising drug candidates. However, the design of potential drug candidates suffers owing to a lack of complete structural information. In this work, homology models of 5-HT1A receptor are generated using two distinct alignments (CW and PSTA) and model building methods (KB and EB). The developed models are validated for virtual screening using a ligand dataset of agonists and antagonists. The best-suited model was efficient in discriminating agonist/antagonist binding. Correlation plots between pKi and docking (R²_agonist≥ 0.6, R²_antagonist≥ 0.7) and MM-GBSA dG bind values (R²_agonist≥ 0.5, R²_antagonist≥ 0.7) revealed optimum corroboration between in vitro and in silico outcomes, which further suggested the usefulness of the developed model for the design of high-affinity probes for the neurological disorders.Communicated by Ramaswamy H. Sarma.

Lysine 101 in the CRAC Motif in Transmembrane Helix 2 Confers Cholesterol-Induced Thermal Stability to the Serotonin1A Receptor

G protein-coupled receptors (GPCRs) constitute the largest class of membrane proteins that transduce signals across the plasma membrane and orchestrate a multitude of physiological processes within cells. The serotonin_1A receptor is a crucial neurotransmitter receptor in the GPCR family involved in a multitude of neurological, behavioral and cognitive functions. We have previously shown, using a combination of experimental and simulation approaches, that membrane cholesterol acts as a key regulator of organization, dynamics, signaling and endocytosis of the serotonin_1A receptor. In addition, we showed that membrane cholesterol stabilizes the serotonin_1A receptor against thermal deactivation. In the present work, we explored the molecular basis of cholesterol-induced thermal stability of the serotonin_1A receptor. For this, we explored the possible role of the K101 residue in a cholesterol recognition/interaction amino acid consensus (CRAC) motif in transmembrane helix 2 in conferring the thermal stability of the serotonin_1A receptor. Our results show that a mutation in the K101 residue leads to loss in thermal stability of the serotonin_1A receptor imparted by cholesterol, independent of membrane cholesterol content. We envision that our results could have potential implications in structural biological advancements of GPCRs and design of thermally stabilized receptors for drug development.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β2-adrenergic receptor (β2AR) and the cholecystokinin (CCK) receptor subfamily. The β2AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β2AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK1R shows functional sensitivity to membrane cholesterol while CCK2R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β2AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-Tm) shows that mutation of a conserved CRAC sequence on TM7 of the β2AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.

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.

Modeling and Structural Characterization of the Sweet Taste Receptor Heterodimer

Sweet taste receptor, a heterodimer belonging to the class C G-protein coupled receptor (GPCR) family and composed of the T1R2 and T1R3 subunits, is responsible for the perception of natural sugars, sweet proteins, various d-amino acids, as well as artificial sweeteners. Despite the critical importance of the sweet receptor not only in mediating gustation but also in its role in the food industry, the architecture of the T1R2-T1R3 complex and the mechanism by which extracellular stimuli induce conformational changes that are propagated to the intracellular milieu, i.e., the signal transduction pathway, remain largely unknown. Here, we constructed and characterized a full-length structural model of the T1R2-T1R3 receptor, including both the transmembrane (TM) and extracellular (EC) domains of the heterodimer, using comparative modeling and extensive all-atom molecular dynamics simulations. Several heterodimer interfaces were first examined for the TM domain, and conformational changes occurring at the intracellular side and associated with the receptor's activation were characterized. From the analysis on the simulated data, putative allosteric binding sites for ligands, ions, and cholesterol were proposed. Also, insights into the protein interface of the TM domain upon activation are provided. The EC domain of the heterodimer, including both the Venus flytrap and cysteine-rich domains, was also investigated. Several important intersubunit interactions located at regions responsible for the receptor's proper function were observed, which resemble those recently identified in other class C GPCR members. Integration of the results from the TM and EC domains facilitates the generation of a full-length T1R2-T1R3 receptor. These findings along with the full-length structural model of the T1R2-T1R3 receptor provide a structural framework that may assist in understanding the mechanistic details associated with the receptor activation process for the sweet T1R2-T1R3 receptor as well as other members of the same family.

Reconstruction of apo A2A receptor activation pathways reveal ligand-competent intermediates and state-dependent cholesterol hotspots

G-protein coupled receptors (GPCRs) play a pivotal role in transmitting signals at the cellular level. Structural insights can be exploited to support GPCR structure-based drug discovery endeavours. Despite advances in GPCR crystallography, active state structures are scarce. Molecular dynamics (MD) simulations have been used to explore the conformational landscape of GPCRs. Efforts have been made to retrieve active state conformations starting from inactive structures, however to date this has not been possible without using an energy bias. Here, we reconstruct the activation pathways of the apo adenosine receptor (A2A), starting from an inactive conformation, by applying adaptive sampling MD combined with a goal-oriented scoring function. The reconstructed pathways reconcile well with experiments and help deepen our understanding of A2A regulatory mechanisms. Exploration of the apo conformational landscape of A2A reveals the existence of ligand-competent states, active intermediates and state-dependent cholesterol hotspots of relevance for drug discovery. To the best of our knowledge this is the first time an activation process has been elucidated for a GPCR starting from an inactive structure only, using a non-biased MD approach, opening avenues for the study of ligand binding to elusive yet pharmacologically relevant GPCR states.

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.

Interfacial Binding Sites for Cholesterol on G Protein-Coupled Receptors

A docking procedure is described that allows the transmembrane surface of a G protein-coupled receptor (GPCR) to be swept rapidly for potential binding sites for cholesterol at the bilayer interfaces on the two sides of the membrane. The procedure matches 89% of the cholesterols resolved in published GPCR crystal structures, when cholesterols likely to be crystal packing artifacts are excluded. Docking poses are shown to form distinct clusters on the protein surface, the clusters corresponding to "greasy hollows" between protein ridges. Docking poses depend on the angle of tilt of the GPCR in the surrounding lipid bilayer. It is suggested that thermal motion could alter the optimal binding pose for a cholesterol molecule, with the range of binding poses within a cluster providing a guide to the range of thermal motions likely for a cholesterol within a binding site.

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.

Peeking at G-protein-coupled receptors through the molecular dynamics keyhole

Molecular dynamics is a state of the art computational tool for the investigation of biophysics phenomenon at a molecular scale, as it enables the modeling of dynamic processes, such as conformational motions, molecular solvation and ligand binding. The recent advances in structural biology have led to a bloom in published G-protein-coupled receptor structures, representing a solid and valuable resource for molecular dynamics studies. During the last decade, indeed, a plethora of physiological and pharmacological facets of this membrane protein superfamily have been addressed by means of molecular dynamics simulations, including the activation mechanism, allosterism and, very recently, biased signaling. Here, we try to recapitulate some of the main contributions that molecular dynamics has recently produced in the field.

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.

State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes

Membranes are known to have modulatory effects on G protein-coupled receptors (GPCRs) via specific lipid interactions. However, the mechanisms of such modulations in physiological conditions and how they influence GPCR functions remain unclear. Here we report coarse-grained molecular dynamics simulations on the Adenosine A2a receptor in different conformational states embedded in an in vivo-mimetic membrane model. Nine lipid interaction sites were revealed. The strength of lipid interactions with these sites showed a degree of dependence on the conformational states of the receptor, suggesting that these lipids may regulate the conformational dynamics of the receptor. In particular, we revealed a dual role of PIP₂ on A2aR activation that involves both stabilization of the characteristic outward tilt of TM6 and enhancement of A2aR-mini-Gs association. Our results demonstrated that the bound lipids allosterically regulate the functional properties of GPCRs. These protein-lipid interactions provide a springboard for design of allosteric modulators of GPCRs.

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.

Structural Changes Due to Antagonist Binding in Ligand Binding Pocket of Androgen Receptor Elucidated Through Molecular Dynamics Simulations

When a small molecule binds to the androgen receptor (AR), a conformational change can occur which impacts subsequent binding of co-regulator proteins and DNA. In order to accurately study this mechanism, the scientific community needs a crystal structure of the Wild type AR (WT-AR) ligand binding domain, bound with antagonist. To address this open need, we leveraged molecular docking and molecular dynamics (MD) simulations to construct a structure of the WT-AR ligand binding domain bound with antagonist bicalutamide. The structure of mutant AR (Mut-AR) bound with this same antagonist informed this study. After molecular docking analysis pinpointed the suitable binding orientation of a ligand in AR, the model was further optimized through 1 μs of MD simulations. Using this approach, three molecular systems were studied: (1) WT-AR bound with agonist R1881, (2) WT-AR bound with antagonist bicalutamide, and (3) Mut-AR bound with bicalutamide. Our structures were very similar to the experimentally determined structures of both WT-AR with R1881 and Mut-AR with bicalutamide, demonstrating the trustworthiness of this approach. In our model, when WT-AR is bound with bicalutamide, Val716/Lys720/Gln733, or Met734/Gln738/Glu897 move and thus disturb the positive and negative charge clumps of the AF2 site. This disruption of the AF2 site is key for understanding the impact of antagonist binding on subsequent co-regulator binding. In conclusion, the antagonist induced structural changes in WT-AR detailed in this study will enable further AR research and will facilitate AR targeting drug discovery.

Identification of Two New Cholesterol Interaction Sites on the A2A Adenosine Receptor

By mole, cholesterol is the most abundant component of animal cell plasma membranes. Many membrane proteins have been shown to be functionally dependent on cholesterol, several of which have also been shown to bind cholesterol at well-defined locations on their membrane-facing surface. In this work, a combination of coarse-grained "Martini" and all-atom simulations are used to identify two, to our knowledge, new cholesterol-binding sites on the A_2A adenosine receptor, a G-protein-coupled receptor that is a target for the treatment of Parkinson's disease. One of the sites is also observed to bind cholesterol in several recent, high-resolution crystal structures of the protein, and in the simulations, interacts with cholesterol only when bound to the inverse agonist ZM241385. Cataloguing cholesterol-binding sites is a vital step in the effort to understand cholesterol-dependent function of membrane proteins. Given that cholesterol content in plasma membranes varies with cell type and on administration of widely prescribed pharmaceuticals, such as statins, understanding cholesterol-dependent function is an important step toward exploiting membrane compositional variation for therapeutic purposes.

Membrane cholesterol effect on the 5-HT2A receptor: Insights into the lipid-induced modulation of an antipsychotic drug target

The serotonin 5-hydroxytryptamine 2A (5-HT_2A ) receptor is a G-protein-coupled receptor (GPCR) relevant for the treatment of CNS disorders. In this regard, neuronal membrane composition in the brain plays a crucial role in the modulation of the receptor functioning. Since cholesterol is an essential component of neuronal membranes, we have studied its effect on the 5-HT_2A receptor dynamics through all-atom MD simulations. We find that the presence of cholesterol in the membrane increases receptor conformational variability in most receptor segments. Importantly, detailed structural analysis indicates that conformational variability goes along with the destabilization of hydrogen bonding networks not only within the receptor but also between receptor and lipids. In addition to increased conformational variability, we also find receptor segments with reduced variability. Our analysis suggests that this increased stabilization is the result of stabilizing effects of tightly bound cholesterol molecules to the receptor surface. Our finding contributes to a better understanding of membrane-induced alterations of receptor dynamics and points to cholesterol-induced stabilizing and destabilizing effects on the conformational variability of GPCRs.

Computational studies of G protein-coupled receptor complexes: Structure and dynamics

G protein-coupled receptors (GPCRs) are ubiquitously expressed transmembrane proteins associated with a wide range of diseases such as Alzheimer's, Parkinson, schizophrenia, and also implicated in in several abnormal heart conditions. As such, this family of receptors is regarded as excellent drug targets. However, due to the high number of intracellular signaling partners, these receptors have a complex interaction networks and it becomes challenging to modulate their function. Experimentally determined structures give detailed information on the salient structural properties of these signaling complexes but they are far away from providing mechanistic insights into the underlying process. This chapter presents some of the computational tools, namely molecular dynamics, molecular docking, and molecular modeling and related analyses methods that have been used to complement experimental findings.

Impact of Lipid Composition and Receptor Conformation on the Spatio-temporal Organization of μ-Opioid Receptors in a Multi-component Plasma Membrane Model

The lipid composition of cell membranes has increasingly been recognized as playing an important role in the function of various membrane proteins, including G Protein-Coupled Receptors (GPCRs). For instance, experimental and computational evidence has pointed to lipids influencing receptor oligomerization directly, by physically interacting with the receptor, and/or indirectly, by altering the bulk properties of the membrane. While the exact role of oligomerization in the function of class A GPCRs such as the μ-opioid receptor (MOR) is still unclear, insight as to how these receptors oligomerize and the relevance of the lipid environment to this phenomenon is crucial to our understanding of receptor function. To examine the effect of lipids and different MOR conformations on receptor oligomerization we carried out extensive coarse-grained molecular dynamics simulations of crystal structures of inactive and/or activated MOR embedded in an idealized mammalian plasma membrane composed of 63 lipid types asymmetrically distributed across the two leaflets. The results of these simulations point, for the first time, to specific direct and indirect effects of the lipids, as well as the receptor conformation, on the spatio-temporal organization of MOR in the plasma membrane. While sphingomyelin-rich, high-order lipid regions near certain transmembrane (TM) helices of MOR induce an effective long-range attractive force on individual protomers, both long-range lipid order and interface formation are found to be conformation dependent, with a larger number of different interfaces formed by inactive MOR compared to active MOR.

The μ-opioid receptor (MOR) is an important pharmaceutical target in the treatment of pain. In order to develop novel pain therapies, devoid of the serious side-effects of present opioid analgesics, we need to understand the fundamentals of how MOR works on the molecular level. While some studies suggest that oligomers of MOR could play a role in signaling, how MOR forms dimers, which interfaces form, and the exact role of oligomers in MOR function remain unclear. While research has shown that the membrane environment can affect membrane protein function, most previous computational work to study oligomerization has been performed in a very simple membrane. Here, we use molecular dynamics simulations of MOR in a heterogeneous plasma membrane model (comprising 63 lipid types) to investigate how the presence of the protein modulates its lipid environment, affecting species distribution and sculpting characteristic order and thickness profiles around the receptors. Such modulations, in turn, induce long-range interactions between the proteins and favor the formation of specific dimeric conformations.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Cholesterol-dependent Conformational Plasticity in GPCR Dimers

The organization and function of the serotonin_1A receptor, an important member of the GPCR family, have been shown to be cholesterol-dependent, although the molecular mechanism is not clear. We performed a comprehensive structural and dynamic analysis of dimerization of the serotonin_1A receptor by coarse-grain molecular dynamics simulations totaling 3.6 ms to explore the molecular details of its cholesterol-dependent association. A major finding is that the plasticity and flexibility of the receptor dimers increase with increased cholesterol concentration. In particular, a dimer interface formed by transmembrane helices I-I was found to be sensitive to cholesterol. The modulation of dimer interface appears to arise from a combination of direct cholesterol occupancy and indirect membrane effects. Interestingly, the presence of cholesterol at the dimer interface is correlated with increased dimer plasticity and flexibility. These results represent an important step in characterizing the molecular interactions in GPCR organization with potential relevance to therapeutic interventions.

Differential dynamics of the serotonin1A receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin1A receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin1A receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3 μs, to analyze the effect of cholesterol on the structure and dynamics of the serotonin1A receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin1A receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.

Cholesterol Corrects Altered Conformation of MHC-II Protein in Leishmania donovani Infected Macrophages: Implication in Therapy

Background

Previously we reported that Kala-azar patients show progressive decrease in serum cholesterol as a function of splenic parasite burden. Splenic macrophages (MΦ) of Leishmania donovani (LD) infected mice show decrease in membrane cholesterol, while LD infected macrophages (I-MΦ) show defective T cell stimulating ability that could be corrected by liposomal delivery of cholesterol. T helper cells recognize peptide antigen in the context of class II MHC molecule. It is known that the conformation of a large number of membrane proteins is dependent on membrane cholesterol. In this investigation we tried to understand the influence of decreased membrane cholesterol in I-MΦ on the conformation of MHC-II protein and peptide-MHC-II stability, and its bearing on the antigen specific T-cell activation.

Methodology/Principal Findings

MΦ of CBA/j mice were infected with Leishmania donovani (I-MΦ). Two different anti-A^κ mAbs were used to monitor the status of MHC-II protein under parasitized condition. One of them (11.5–2) was conformation specific, whereas the other one (10.2.16) was not. Under parasitized condition, the binding of 11.5–2 decreased significantly with respect to the normal counterpart, whereas that of 10.2.16 remained unaltered. The binding of 11.5–2 was restored to normal upon liposomal delivery of cholesterol in I-MΦ. By molecular dynamics (MD) simulation studies we found that there was considerable conformational fluctuation in the transmembrane domain of the MHC-II protein in the presence of membrane cholesterol than in its absence, which possibly influenced the distal peptide binding groove. This was evident from the faster dissociation of the cognate peptide from peptide-MHC complex under parasitized condition, which could be corrected by liposomal delivery of cholesterol in I-MΦ.

Conclusion

The decrease in membrane cholesterol in I-MΦ may lead to altered conformation of MHC II, and this may contribute to a faster dissociation of the peptide. Furthermore, liposomal delivery of cholesterol in I-MΦ restored its normal antigen presenting function. This observation brings strength to our previous observation on host directed therapeutic application of liposomal cholesterol in experimental visceral leishmaniasis.

The disease visceral leishmaniasis is caused by the protozoan parasite Leishmania donovani (LD). One of the hallmarks of the disease is immune suppression. The parasites replicate within the macrophages and dendritic cells, and such cells are known as antigen presenting cells (APCs). APCs present peptide to T-helper cells in association with the transplantation antigen-II (MHC-II). The infected macrophages show decrease in membrane cholesterol leading to increase in membrane fluidity. The membrane cholesterol is important for maintaining conformation of membrane proteins. Here we show that conformation of MHC-II protein is altered in parasitized macrophages, which results faster dissociation of peptide from peptide-MHC-II complex as compared to normal counterpart. The conformational change in MHC-II protein is also supported by molecular dynamic simulation studies, as there is considerable structural fluctuation of MHC-II peptide binding domain in presence and absence of cholesterol. This observation indicated that cholesterol is important for maintaining conformation of MHC-II protein and stability of the peptide-MHC complex. Thus, Leishmania parasites by modulating membrane cholesterol influence above processes leading to defective T-cell stimulation in leishmaniasis. The above defects displayed by infected macrophages could be corrected by liposomal delivery of cholesterol indicating a possible therapeutic role of liposomal cholesterol in infection.

Cholesterol in the rod outer segment: A complex role in a "simple" system

The rod outer segment (ROS) of retinal photoreceptor cells consists of disk membranes surrounded by the plasma membrane. It is a relatively uncomplicated system in which to investigate cholesterol distribution and its functional consequences in biologically relevant membranes. The light sensitive protein, rhodopsin is the major protein in both membranes, but the lipid compositions are significantly different in the disk and plasma membranes. Cholesterol is high in the ROS plasma membrane. Disk membranes are synthesized at the base of the ROS and are also high in cholesterol. However, cholesterol is rapidly depleted as the disks are apically displaced. During this apical displacement the disk phospholipid fatty acyl chains become progressively more unsaturated, which creates an environment unfavorable to cholesterol. Membrane cholesterol has functional consequences. The high cholesterol found in the plasma membrane and in newly synthesized disks inhibits the activation of rhodopsin. As disks are apically displaced and cholesterol is depleted rhodopsin becomes more responsive to light. This effect of cholesterol on rhodopsin activation has been shown in both native and reconstituted membranes. The modulation of activity can be at least partially explained by the effect of cholesterol on bulk lipid properties. Cholesterol decreases the partial free volume of the hydrocarbon region of the bilayer and thereby inhibits rhodopsin conformational changes required for activation. However, cholesterol binds to rhodopsin and may directly affect the protein also. Furthermore, cholesterol stabilizes rhodopsin to thermal denaturation. The membrane must provide an environment that allows rhodopsin conformational changes required for activation while also stabilizing the protein to thermal denaturation. Cholesterol thus plays a complex role in modulating the activity and stability of rhodopsin, which have implications for other G-protein coupled receptors.

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.

The membrane as the gatekeeper of infection: Cholesterol in host-pathogen interaction

The cellular plasma membrane serves as a portal for the entry of intracellular pathogens. An essential step for an intracellular pathogen to gain entry into a host cell therefore is to be able to cross the cell membrane. In this review, we highlight the role of host membrane cholesterol in regulating the entry of intracellular pathogens using insights obtained from work on the interaction of Leishmania and Mycobacterium with host cells. The entry of these pathogens is known to be dependent on host membrane cholesterol. Importantly, pathogen entry is inhibited either upon depletion (or complexation), or enrichment of membrane cholesterol. In other words, an optimum level of host membrane cholesterol is necessary for efficient infection by pathogens. In this overall context, we propose a general mechanism, based on cholesterol-induced conformational changes, involving cholesterol binding sites in host cell surface receptors that are implicated in this process. A therapeutic strategy targeting modulation of membrane cholesterol would have the advantage of avoiding the commonly encountered problem of drug resistance in tackling infection by intracellular pathogens. Insights into the role of host membrane cholesterol in pathogen entry would be instrumental in the development of novel therapeutic strategies to effectively tackle intracellular pathogenesis.

Applying high-performance computing in drug discovery and molecular simulation

In recent decades, high-performance computing (HPC) technologies and supercomputers in China have significantly advanced, resulting in remarkable achievements. Computational drug discovery and design, which is based on HPC and combines pharmaceutical chemistry and computational biology, has become a critical approach in drug research and development and is financially supported by the Chinese government. This approach has yielded a series of new algorithms in drug design, as well as new software and databases. This review mainly focuses on the application of HPC to the fields of drug discovery and molecular simulation at the Chinese Academy of Sciences, including virtual drug screening, molecular dynamics simulation, and protein folding. In addition, the potential future application of HPC in precision medicine is briefly discussed.

Characterizing rhodopsin signaling by EPR spectroscopy: from structure to dynamics

Electron paramagnetic resonance (EPR) spectroscopy, together with spin labeling techniques, has played a major role in the characterization of rhodopsin, the photoreceptor protein and G protein-coupled receptor (GPCR) in rod cells. Two decades ago, these biophysical tools were the first to identify transmembrane helical movements in rhodopsin upon photo-activation, a critical step in the study of GPCR signaling. EPR methods were employed to identify functional loop dynamics within rhodopsin, to measure light-induced millisecond timescale changes in rhodopsin conformation, to characterize the effects of partial agonists on the apoprotein opsin, and to study lipid interactions with rhodopsin. With the emergence of advanced pulsed EPR techniques, the stage was set to determine the amplitude of structural changes in rhodopsin and the dynamics in the rhodopsin signaling complexes. Work in this area has yielded invaluable information about mechanistic properties of GPCRs. Using EPR techniques, receptors are studied in native-like membrane environments and the effects of lipids on conformational equilibria can be explored. This perspective addresses the impact of EPR methods on rhodopsin and GPCR structural biology, highlighting historical discoveries made with spin labeling techniques, and outlining exciting new directions in the field.

Binding of the Antagonist Caffeine to the Human Adenosine Receptor hA2AR in Nearly Physiological Conditions

Lipid composition may significantly affect membrane proteins function, yet its impact on the protein structural determinants is not well understood. Here we present a comparative molecular dynamics (MD) study of the human adenosine receptor type 2A (hA(2A)R) in complex with caffeine--a system of high neuro-pharmacological relevance--within different membrane types. These are POPC, mixed POPC/POPE and cholesterol-rich membranes. 0.8-μs MD simulations unambiguously show that the helical folding of the amphipathic helix 8 depends on membrane contents. Most importantly, the distinct cholesterol binding into the cleft between helix 1 and 2 stabilizes a specific caffeine-binding pose against others visited during the simulation. Hence, cholesterol presence (~33%-50% in synaptic membrane in central nervous system), often neglected in X-ray determination of membrane proteins, affects the population of the ligand binding poses. We conclude that including a correct description of neuronal membranes may be very important for computer-aided design of ligands targeting hA(2A)R and possibly other GPCRs.

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.

Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins

The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.

Coarse-grained molecular dynamics provides insight into the interactions of lipids and cholesterol with rhodopsin

Protein function is a complicated interplay between structure and dynamics, which can be heavily influenced by environmental factors and conditions. This is particularly true in the case of membrane proteins, such as the visual receptor rhodopsin. It has been well documented that lipid headgroups, polyunsaturated tails, and the concentration of cholesterol in membranes all play a role in the function of rhodopsin. Recently, we used all-atom simulations to demonstrate that different lipid species have preferential interactions and possible binding sites on rhodopsin's surface, consistent with experiment. However, the limited timescales of the simulations meant that the statistical uncertainty of these results was substantial. Accordingly, we present here 32 independent 1.6 μs coarse-grained simulations exploring lipids and cholesterols surrounding rhodopsin and opsin, in lipid bilayers mimicking those found naturally. Our results agree with those found experimentally and in previous simulations, but with far better statistical certainty. The results demonstrate the value of combining all-atom and coarse-grained models with experiment to provide a well-rounded view of lipid-protein interactions.

Structure and dynamics of G-protein coupled receptors

G-protein coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell-to-cell communication. They also form the largest superfamily of drug targets. Hence detailed studies of the three dimensional structure and dynamics are critical to understanding the functional role of GPCRs in signal transduction pathways, and for drug design. In this chapter we compare the features of the crystal structures of various biogenic amine receptors, such as β1 and β2 adrenergic receptors, dopamine D3 receptor, M2 and M3 muscarinic acetylcholine receptors. This analysis revealed that conserved residues are located facing the inside of the transmembrane domain in these GPCRs improving the efficiency of packing of these structures. The NMR structure of the chemokine receptor CXCR1 without any ligand bound, shows significant dynamics of the transmembrane domain, especially the helical kink angle on the transmembrane helix6. The activation mechanism of the β2-adrenergic receptor has been studied using multiscale computational methods. The results of these studies showed that the receptor without any ligand bound, samples conformations that resemble some of the structural characteristics of the active state of the receptor. Ligand binding stabilizes some of the conformations already sampled by the apo receptor. This was later observed in the NMR study of the dynamics of human β2-adrenergic receptor. The dynamic nature of GPCRs leads to a challenge in obtaining purified receptors for biophysical studies. Deriving thermostable mutants of GPCRs has been a successful strategy to reduce the conformational heterogeneity and stabilize the receptors. This has lead to several crystal structures of GPCRs. However, the cause of how these mutations lead to thermostability is not clear. Computational studies are beginning to shed some insight into the possible structural basis for the thermostability. Molecular Dynamics simulations studying the conformational ensemble of thermostable mutants have shown that the stability could arise from both enthalpic and entropic factors. There are regions of high stress in the wild type GPCR that gets relieved upon mutation conferring thermostability.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Dynamic PIP2 interactions with voltage sensor elements contribute to KCNQ2 channel gating

The S4 segment and the S4-S5 linker of voltage-gated potassium (Kv) channels are crucial for voltage sensing. Previous studies on the Shaker and Kv1.2 channels have shown that phosphatidylinositol-4,5-bisphosphate (PIP2) exerts opposing effects on Kv channels, up-regulating the current amplitude, while decreasing the voltage sensitivity. Interactions between PIP2 and the S4 segment or the S4-S5 linker in the closed state have been highlighted to explain the effects of PIP2 on voltage sensitivity. Here, we show that PIP2 preferentially interacts with the S4-S5 linker in the open-state KCNQ2 (Kv7.2) channel, whereas it contacts the S2-S3 loop in the closed state. These interactions are different from the PIP2-Shaker and PIP2-Kv1.2 interactions. Consistently, PIP2 exerts different effects on KCNQ2 relative to the Shaker and Kv1.2 channels; PIP2 up-regulates both the current amplitude and voltage sensitivity of the KCNQ2 channel. Disruption of the interaction of PIP2 with the S4-S5 linker by a single mutation decreases the voltage sensitivity and current amplitude, whereas disruption of the interaction with the S2-S3 loop does not alter voltage sensitivity. These results provide insight into the mechanism of PIP2 action on KCNQ channels. In the closed state, PIP2 is anchored at the S2-S3 loop; upon channel activation, PIP2 interacts with the S4-S5 linker and is involved in channel gating.

Non-covalent binding of membrane lipids to membrane proteins

Polar lipids and membrane proteins are major components of biological membranes, both cell membranes and membranes of enveloped viruses. How these two classes of membrane components interact with each other to influence the function of biological membranes is a fundamental question that has attracted intense interest since the origins of the field of membrane studies. One of the most powerful ideas that driven the field is the likelihood that lipids bind to membrane proteins at specific sites, modulating protein structure and function. However only relatively recently has high resolution structure determination of membrane proteins progressed to the point of providing atomic level structure of lipid binding sites on membrane proteins. Analysis of X-ray diffraction, electron crystallography and NMR data over 100 specific lipid binding sites on membrane proteins. These data demonstrate tight lipid binding of both phospholipids and cholesterol to membrane proteins. Membrane lipids bind to membrane proteins by their headgroups, or by their acyl chains, or binding is mediated by the entire lipid molecule. When headgroups bind, binding is stabilized by polar interactions between lipid headgroups and the protein. When acyl chains bind, van der Waals effects dominate as the acyl chains adopt conformations that complement particular sites on the rough protein surface. No generally applicable motifs for binding have yet emerged. Previously published biochemical and biophysical data link this binding with function. This Article is Part of a Special Issue Entitled: Membrane Structure and Function: Relevance in the Cell's Physiology, Pathology and Therapy.

Cholesterol lowering drug may influence cellular immune response by altering MHC II function

Major histocompatibility complex class II (MHC II) expressed on the surface of antigen-presenting cells (APCs) displays peptides to CD4⁺ T cells. Depletion of membrane cholesterol from APCs by methyl β-cyclodextrin treatment compromises peptide-MHC II complex formation coupled with impaired binding of conformational antibody, which binds close to the peptide binding groove of MHC II. Interestingly, the total cell surface of MHC II remains unaltered. These defects can be corrected by restoring membrane cholesterol. In silico docking studies with a three-dimensional model showed the presence of a cholesterol binding site in the transmembrane domain of MHC II (TM-MHC-II). From the binding studies it was clear that cholesterol, indeed, interacts with the TM-MHC-II and alters its conformation. Mutation of cholesterol binding residues (F240, L243, and F246) in the TM-MHC-II decreased the affinity for cholesterol. Furthermore, transfection of CHO cells with full-length mutant MHC II, but not wild-type MHC II, failed to activate antigen-specific T cells coupled with decreased binding of conformation-specific antibodies. Thus, cholesterol-induced conformational change of TM-MHC-II may allosterically modulate the peptide binding groove of MHC II leading to T cell activation.