Exploring GPCR–Lipid Interactions by Molecular Dynamics Simulations

Apr 23, 2018 - *Phone: +91-20-2590-2408. E-mail: [email protected]., *Phone: +91-40-2719-2578. E-mail: [email protected]. ACS Editors' Choice - Thi...
0 downloads 4 Views 15MB Size
Subscriber access provided by Warwick University Library

Exploring GPCR-lipid Interactions by Molecular Dynamics Simulations: Excitements, Challenges and the Way Forward Durba Sengupta, Xavier Prasanna, Madhura Mohole, and Amitabha Chattopadhyay J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.8b01657 • Publication Date (Web): 23 Apr 2018 Downloaded from http://pubs.acs.org on April 24, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Invited Feature Article

Exploring GPCR-lipid Interactions by Molecular Dynamics Simulations: Excitements, Challenges and the Way Forward

Durba Sengupta,†,§,* Xavier Prasanna,† Madhura Mohole†,§ and Amitabha Chattopadhyay‡,*



CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411 008, India § Academy of Scientific and Innovative Research, India ‡ CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad 500 007, India

*Address correspondence to Amitabha Chattopadhyay, Tel: +91-40-2719-2578, E-mail: [email protected]; or Durba Sengupta, Tel: +91-20-2590-2408, E-mail: [email protected]

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 35

Abstract G protein-coupled receptors (GPCRs) are seven transmembrane receptors that mediate a large number of cellular responses and are important drug targets. One of the current challenges in GPCR biology is to analyze the molecular signatures of receptor-lipid interactions, and their subsequent effects on GPCR structure, organization and function. Molecular dynamics simulation studies have been successful in predicting molecular determinants of receptor-lipid interactions. In particular, predicted cholesterol interaction sites appear to correspond well with experimentally determined binding sites and estimated time scales of association. In spite of several success stories, the methodologies in molecular dynamics simulations are still emerging. In this Feature Article, we provide a comprehensive overview of coarse-grain and atomistic molecular dynamics simulations of GPCR-lipid interaction in the context of experimental observations. In addition, we discuss the effect of secondary and tertiary structural constraints in coarse-grain simulations in the context of functional dynamics and structural plasticity of GPCRs.

We envision that this

comprehensive overview will help resolve differences in computational studies and provide a way forward.

2

ACS Paragon Plus Environment

Page 3 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

GPCR-lipid Interactions: Where are We Now ?

G protein-coupled receptors (GPCRs) are transmembrane proteins that connect the outside of the cell to inner milieu and act as main hubs in several cell signaling networks.1-3 The major function of GPCRs is information transfer across the cellular plasma membrane. Consequently, they represent major (~50%) drug targets in all clinical areas.4,5

These

receptors are composed of seven transmembrane helices, interconnected by extracellular and intracellular loops.6 The orthosteric ligand binding site of representative GPCRs has been reported to be at the lumen of the transmembrane helices.7 Based on the structural and functional characterization of the receptor, it has been suggested that a dynamic coupling exists between the ligand binding site and the G-protein coupling site that transduces the signal upon ligand binding.8,9 Receptors exhibiting basal level activity such as the β2adrenergic receptor display large dynamics, even in the ligand-free state, that directly relates to their function.9 As a result, GPCR flexibility and dynamics play a critical role in receptor function.10,11 In other words, GPCRs are functional due to their conformational flexibility. In spite of the importance of receptor dynamics in GPCR function, it poses a challenge in structural characterization of GPCRs.

For ease of crystallization, GPCRs are heavily

engineered12-14 in several domains of the receptor, including the extracellular and intracellular loops that modulate ligand binding and function.15 More importantly, crystallization conditions usually consist of non-native, non-lamellar lipid environments.16,17 Consequently, the static structures resolved by crystallography provide important structural information, but may fall short in revealing information on inherent receptor dynamics. A critical aspect of GPCRs is their cross-talk with membrane lipids. The serotonin1A receptor is one of the first GPCRs that was reported to exhibit cholesterol-dependent functional modulation, both in terms of ligand binding and G-protein coupling.18,19

A

variety of experiments using diverse approaches showed that membrane cholesterol modulates the function, stability, organization and oligomerization of the receptor.20-24 The list of GPCRs exhibiting cholesterol-dependent activity and signaling is constantly increasing.

We previously compiled a total of ~15 GPCRs that exhibited cholesterol3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

dependent function.19,25 Subsequently, we have shown cholesterol-dependent signaling for αand β-adrenergic receptors,26 and T2R4 receptor27 (a representative member of the bitter taste receptor family). According to current estimate, more than 50 GPCRs are believed to exhibit cholesterol-dependent function. A comprehensive analysis of such dependence is in progress (Sarkar, P.; Chattopadhyay, A. unpublished observations) and beyond the scope of this article.

Interestingly, specific cholesterol binding sites in GPCR crystal structures

have been reported by x-ray studies,28 while NMR studies suggest that these interactions occur at nanosecond and microsecond timescales.29 Importantly, it has been suggested that GPCRs exhibit cholesterol-dependent conformational dynamics, although molecular details of the conformational dynamics are yet to be resolved.3 Additional membrane components, such as phospholipids and sphingolipids have been reported to modulate ligand binding in GPCRs.30-32 In fact, phospholipids have been suggested to represent a special class of allosteric modulators for GPCRs.33 To gain insight into GPCR-membrane interactions, it would therefore be important to take into account the functional dynamics of GPCRs. In spite of the considerable progress made toward understanding GPCR structure and dynamics by experimental methods in the last few years,34 the information obtained is still limited in terms of spatial and temporal resolution.35 Recent molecular dynamics simulations have made large strides toward complementing experimental data by providing a comprehensive molecular picture of membrane proteins in general and GPCRs in particular.8,35,36

Receptor activation represents an important aspect where atomistic

simulations have helped uncover molecular details in close collaboration with experimental observations. We suggest a recent review by Dror and co-workers to interested readers for an overview of the activation mechanisms.8 However, studies on the interaction of GPCRs with membrane lipids are relatively sparse. In this review, we focus on coarse-grain simulations of GPCR-lipid interactions and how these interactions modulate GPCR structure and organization. In the process, we provide a comprehensive overview of the major success stories in addressing GPCR-membrane interactions using computational approaches.

4

ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Emerging Success Stories: Cholesterol Hot-spots in β2-adrenergic Receptor Compare Well between Experiments and Simulations The β2-adrenergic receptor is an important member of the GPCR family and serves as an excellent prototype for monitoring GPCR structure and function. The receptor is involved in muscle relaxation following activation37 and its dysfunction is associated with cardiac diseases and asthma.38,39

The structural characterization of the β2-adrenergic receptor by

high resolution crystallography provided an important step toward understanding molecular level details of GPCR function.12,13

Subsequently, Kobilka and co-workers were able to

resolve two cholesterol molecules at the groove formed by transmembrane helices I, II, III and IV. This site in the crystal structure of the β2-adrenergic receptor was termed the cholesterol consensus motif (CCM).28 The main residues comprising this site are shown in Figure 1. The hydroxyl moiety of the cholesterol is found to interact with an Arg residue, while the sterol ring interacts with Trp and Ile residues on transmembrane helix IV. Although the CCM site was resolved experimentally, GPCRs may behave differently in cubic and lamellar lipidic mesophases17 and the CCM could be specific to membrane lipid environment. To analyze cholesterol sites for GPCRs in membrane bilayers, we previously performed coarse-grain simulations of the β2-adrenergic receptor embedded in 1-palmitoyl-2oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers at varying cholesterol concentration using the MARTINI coarse-grain force-field.40 The simulations were able to identify the CCM site as a high occupancy site. A representative snapshot from the simulations is shown in Figure 1b. The CCM site was identified by calculating the maximum occupancy of cholesterol, i.e., the maximum time a cholesterol molecule is bound to each site of the receptor, during the course of the simulations. A high cholesterol occupancy was observed at transmembrane helix IV in multiple simulations at three different cholesterol concentrations, thereby indicating that a cholesterol molecule consistently interacted with the receptor at the CCM site. The close comparison between the simulations and the experimentally-resolved site (see Figure 1) is noteworthy and points to a high predictive power of the simulation 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

methodology used. Interestingly, a large dynamics of cholesterol was observed at the CCM site during the course of the simulation. The cholesterol molecule at the CCM site interacted with several adjacent sites in transmembrane helix IV. These sites are represented in Figure 1(c-f). In particular, a site closer to the center of the helix, that displayed prominent ring stacking interactions with the Trp residue of the CCM, was observed to be favorable for cholesterol interaction. The high entropy exhibited by the cholesterol at this site would contribute to the total free energy of interaction at this site. An important aspect of β2-adrenergic receptor-cholesterol interaction is the presence of multiple occupancy sites. The spatial density of cholesterol, averaged over the trajectory, is shown in Figure 2 with the top view of the receptor superimposed. For infinitely long sampling, the total occupancy of cholesterol corresponds to the average spatial density. Several cholesterol sites could be distinguished, with those in the upper leaflet marked e1-3, the inner leaflet i2-5 and a high density site at the center of the bilayer marked m1. The site i3 corresponds to the CCM site discussed above. However, an inter-site hopping between sites i3 and m1 was also observed. The interaction sites were found to be consistent with the study of Cang and co-workers who used atomistic simulations to identify cholesterol interaction sites in the β2-adrenergic receptor.41 Importantly, coarse-grain simulations were able to reproduce most sites identified by atomistic simulations, although the site marked i1 was not observed. Recent work from the Vattulainen group has confirmed the sites i2, i4/5 and e2, although the remaining sites were not identified as high density sites.42 In spite of these differences, the main occupancy sites determined from two independent simulation techniques are similar. These observations point toward force-field independent results and the high accuracy of the coarse-grain methodology. The time scales of cholesterol interaction from simulations were observed to be microseconds for the high occupancy sites and nanoseconds for the low occupancy sites. The microsecond time scale dynamics of high interaction sites, although consistent with nonannular sites, appeared to be faster than previously suggested.43 Interestingly, an NMR study analyzing the interaction of membrane cholesterol with the β2-adrenergic receptor was able to identify two distinct time scales of associations that the authors assigned to the nanosecond 6

ACS Paragon Plus Environment

Page 6 of 35

Page 7 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

and microsecond time scales.29

The time scales reported in the simulations should be

considered as a qualitative description due to the nature of the coarse-grain force-field. Nonetheless, these results emphasize high predictive power of the simulation methodology.

Comparison of Cholesterol Interaction Sites in the Serotonin1A Receptor and β2-adrenergic Receptor The serotonin1A receptor is one of the most comprehensively studied seven transmembrane domain GPCRs in the serotonin receptor family.44,45 As a consequence of its indispensable role in neurological functions, the serotonin1A receptor has emerged as a major drug target in the development of therapeutics against neuropsychiatric disorders such as anxiety, depression, schizophrenia and Parkinson’s disease.46 The cholesterol occupancy sites of the serotonin1A receptor have been calculated at varying cholesterol concentration, using a similar coarse-grain methodology.47 The serotonin1A receptor shares a high sequence similarity with the β2-adrenergic receptor.48 In fact, the serotonin1A receptor was initially discovered as an ‘orphan’ β-adrenergic receptor and was later identified as a serotonin receptor (‘deorphanized’).49 It was not a priori clear whether the GPCR-lipid interactions of the two receptors would differ or not. Coarse-grain simulations were able to identify a large number of cholesterol interaction sites in the serotonin1A receptor.47 Interestingly, the CCM site at transmembrane helix IV could not be distinguished as a high-occupancy cholesterol site. Multiple sites were observed at transmembrane helices I, II, V, VI and VII. Of these, the site at transmembrane helix V corresponds to a putative cholesterol recognition motif, the CRAC motif.50 The site comprises of -L/V-(X)1-5-Y-(X)1-5-R/K- residues, that has also been identified in other cholesterol binding proteins. Subsequently, simulations confirmed the stochasticity of cholesterol interactions.51 To analyze the spatial distribution of cholesterol, we re-analyzed our simulations51 and plotted the spatial distribution function of cholesterol around the receptor for the upper and lower leaflets separately (see Figure 3). As expected, a high cholesterol density site was observed at the lower leaflet at transmembrane helix V near the CRAC site of the serotonin1A 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

receptor. Similarly, transmembrane helix I was found to have a high cholesterol density in its vicinity in the lower leaflet. In the upper leaflet, transmembrane helix VI exhibited high cholesterol density. Surprisingly, we observed high cholesterol density near transmembrane helix IV, but this was not identified as a high occupancy site. It is possible that due to the high dynamics at this region, it has a lower maximum occupancy for cholesterol and longer sampling could help resolve the site. Of these sites, transmembrane helix V was identified as a cholesterol hot-spot by all-atom simulations.52 The close comparison of cholesterol sites from atomistic and coarse-grain simulations clearly demonstrates that these cholesterol interaction sites are force-field independent. In the next step, we analyzed the effect of the cholesterol association sites on the dimerization of the receptor.40,51 Receptor association has been suggested to be modulated by direct effects via cholesterol interaction sites or indirect effects of membrane bilayer properties.53 Simulations indicate that a combination of these effects could modulate the association of the receptors.40,51,53-60

In general, the main sites of receptor association

correspond to those with high membrane perturbations, but these are modulated by specific cholesterol interactions. For example, in case of the β2-adrenergic receptor, high cholesterol occupancy at transmembrane helix IV (from the CCM sites i3 and m1) of the receptor blocks it from the dimer interface at high cholesterol concentrations.40 On the other hand, the cholesterol interaction sites of comparable energetics in the serotonin1A receptor appear to induce greater plasticity in receptor association and as a consequence, several dimer interfaces were identified for the serotonin1A receptor.51 Interestingly, the flexibility of the dimer interface itself was observed to be correlated with the presence of cholesterol at the dimer interface. With increasing computational resources, the next step would be to analyze the molecular details of the role of membrane cholesterol in the oligomerization of these receptors.

Phospholipid Interaction Sites in the Serotonin1A Receptor Match Experimental Phospholipid Sites in Related GPCRs 8

ACS Paragon Plus Environment

Page 8 of 35

Page 9 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

A related question that arises is whether there are any phospholipid interaction sites. A recent crystal structure of the adenosine2A receptor revealed multiple electron density sites corresponding to phospholipid interaction sites.61 Prominent amongst them is the site at the cleft formed by transmembrane helices I and VII (Figure 4). In a series of coarse-grain simulations, we were able to identify the same site as a phospholipid binding site in both β2adrenergic receptor59 and the serotonin1A receptor.51,52 The site was identified as a high occupancy phospholipid site using coarse-grain simulations in β2-adrenergic receptor,59 and both coarse-grain51 and atomistic simulations52 in the serotonin1A receptor. Interestingly, in the serotonin1A receptor, the site appears to be dependent on cholesterol concentration.51 In the absence of cholesterol, or at low cholesterol concentrations, this site is the highest occupancy site, but another high occupancy site was observed at transmembrane helix V with increasing cholesterol concentration.51 Although further work is needed to understand the role of these phospholipid sites, it is encouraging that simulations are able to accurately predict important phosopholipid sites.

Serotonin1A Receptor Binds GM1 Ganglioside via a Sphingolipid Binding Domain

Sphingolipids are essential components of the cell membrane and are recognized as diverse and dynamic regulators of a multitude of cellular processes. Sphingolipids have been shown to regulate the function of GPCRs, such as the serotonin1A receptor.32,62-66

In

particular, ligand binding was shown to be sphingolipid-dependent, but the underlying molecular mechanism was not clear.

Similar to phospholipids and cholesterol, the

modulatory effect of sphingolipids on GPCR function could be a result of direct or indirect interaction, or a combination of both.67 Direct interactions are supported by the fact that several sphingolipid-dependent membrane proteins appear to have a consensus ‘sphingolipid binding domain’ (SBD) that consists of a characteristic combination of aromatic, basic and turn-inducing residues.68,69 The SBD motif has been reported in serotonin receptors and shown to be evolutionarily conserved.68 Interestingly, the motif is present at one of the 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 35

extracellular loops, and not in the transmembrane region of the serotonin1A receptor. However, experiments with isolated SBD peptide did not show significant binding in model membranes,69 thereby highlighting the relevance of the overall ‘context’. In addition, it has been recently reported that representative GPCRs such as the serotonin1A receptor contains a putative sphingolipid-binding motif (SBM), an intrinsic feature of the receptor and is conserved throughout natural evolution.70 We investigated the interaction of GM1 ganglioside with the serotonin1A receptor using coarse-grain molecular dynamics simulations.71 Our results showed that GM1 binds to the predicted SBD in the extracellular loop 1 region of the serotonin1A receptor (see Figure 5). The sugar moiety of GM1 ganglioside was found to interact with the aromatic residue W102, and flanking residues, K101 and T103. An interesting aspect of the interaction site is that though it is independent of cholesterol, the presence of cholesterol allows a closer interaction of GM1 with the receptor. The importance of cholesterol has been previously reported for the interactions of the HIV-1 gp120 glycoprotein with a related glycosphingolipid (Gb3).72 In addition, we observed that the interaction of the serotonin1A receptor with GM1 stabilizes a “flip-out” conformation of W120 such that it points away from the central lumen of the receptor. This conformation is dependent on the cholesterolmodulated GM1 distribution around the receptor. A similar orientation of aromatic residues has been reported in the binding subunit of Shiga-like toxin.73

Current Methodologies: Advantages and Limitations

Molecular dynamics simulations have made enormous strides in the last few years in terms of both accuracy and computational power. In this review, we have highlighted coarsegrain simulation studies on GPCR-lipid interactions that closely match current experimental data and have high predictive power. Despite these successes, several strategies have been used by various groups, and the choice of the methodology lacks consensus.

One of the

main questions that arises is whether a coarse-grain force-field can have the accuracy of 10

ACS Paragon Plus Environment

Page 11 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

atomistic force-fields, despite a lower chemical resolution. Even when a coarse-grain forcefield is chosen, parameter sets used for different applications vary. We discuss below the advantages and limitations of each of these methods in order to compare their suitability for exploring multiple aspects of GPCR-lipid interactions.

Distinguishing Between Atomistic and Coarse-grain Simulations

In conjunction to coarse-grain simulation studies, several atomistic simulations of GPCR-lipid interactions have been reported.74 Atomistic simulations have been able to leapfrog the nanosecond time scale toward the microsecond regime with increasing computational power.

Atomistic simulations have been able to identify important

interactions of GPCRs with lipids,75,76 and a few have been successful in mapping out multiple cholesterol interaction sites.41,42,77,78

Interestingly, a competition between

cholesterol and the ligand at the orthosteric ligand binding site has been reported.79 Interactions with membrane lipids has been shown to result in differences in the overall conformational dynamics of transmembrane domains and loop regions of GPCRs.42,52,80 Several atomistic force-fields, including CHARMM,41,52 OPLS42 and Desmond force-fields78 have been used for these studies. Detailed comparisons between these force-fields have not been reported, although most cholesterol sites that have been identified appear to be forcefield independent. In addition, a rigorous comparison between GPCR-lipid interactions predicted using atomistic simulations and coarse-grain simulations is still missing. We have highlighted a few cases in this review where comparisons can be directly made. In general, atomistic simulations have inherently an increased chemical resolution and as a result a higher predictive power, but the sampling of these interaction sites is limited. The majority of atomistic simulations are able to sample at the microsecond time regime, and it remains difficult for these simulations to comprehensively sample multi microsecond time scale events.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Coarse-grain force-fields, especially the MARTINI force-field, are being increasingly used to analyze GPCR-lipid interactions.81,82

The main advantage of coarse-grain

simulations is the increased sampling efficiency.

Simulation trajectories of multiple

microsecond time scales are easily achievable with current computational resources. The MARTINI force-field is one of the most widely used coarse-grain force-fields mainly due to its transferability and predictive power.83

The force-field was parameterized based on

solvation and partitioning free energies and has been successful in analyzing facets of several membrane proteins.82,83 In particular, the protein-lipid interactions in model systems have been reported to be accurately represented.82 The main limitations of the MARTINI forcefield are simplistic electrostatics, and reduced structural flexibility due to the use of positional restraints on the secondary structural elements.82,83 This is of critical disadvantage in systems where large secondary structural dynamics, such as helix folding/unfolding takes place. The way forward could perhaps be an efficient combination of coarse-grain and atomistic approaches. Overall, the MARTINI force-field appears to be well suited to study GPCR-lipid interactions and has been reported to match previous atomistic simulations, as well as experimental observations.

To Use or Not to Use Artificial Tertiary Structural Constraints in Coarse-grain Simulations? An important methodological difference within the MARTINI force-field is the use of the elastic network model, that is a set of constraints to hold the protein tertiary structure constant. The elastic network model was parameterized by Periole and co-workers, and allows only the harmonic protein dynamics.60,84 Anharmonic modes as well as large-scale deviations from the starting structure (such as helix repositioning) are not permitted by the elastic network. Originally, the MARTINI force-field was parameterized without the elastic network,81,82 and several membrane protein studies have been reported without the use of the elastic network.85-88 A large number of GPCR simulations have been reported without the underlying elastic network40,47,51 or with only an early version of the elastic network.60 12

ACS Paragon Plus Environment

Page 12 of 35

Page 13 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

However, large conformational variations have also been reported in the absence of the elastic network in certain GPCRs.89 Interestingly, several coarse-grain studies described above that are shown to reproduce binding sites determined from experiments or atomistic simulations, were performed without the elastic network. In this regard, it is important to note that the flexibility of GPCRs has been reported to be critical for its function,8,10,11,90,91 and as a first estimate, a more flexible structure is probably a better model for GPCRs. Additionally, the MARTINI coarse-grain force-field has been shown to reproduce transmembrane helix packing92 and association energetics93-95 in single transmembrane helical proteins. Helix packing in the larger GPCRs has not been critically tested, but helixhelix interactions predicted by the MARTINI force-field without considering tertiary structural restraints was shown to reproduce experimental structural data.96 Taken together, the MARTINI force-field appears to have high predictive power and is uniquely positioned to analyze molecular details of GPCR-lipid interactions and its effect on GPCR organization.

The Way Forward

Although simulations have been used to analyze several facets of GPCR-lipid interactions, standardized methodologies are still emerging for different systems. Detailed thermodynamic and kinetic analysis of lipid binding sites that can be achieved by simulations with improved computational resources would provide an unparalleled insight into GPCRlipid interactions. Currently, it is important to compare results obtained from computational approaches with experimental data and test the robustness of the results. As a cautionary note, each receptor should perhaps be tested individually for correct use of the simulation parameters. Another important step could be to match aspects of atomistic and coarse-grain simulations, and a close comparison with experimental data, whenever possible. An emerging concept from molecular dynamics simulations is the large dynamics in GPCR-lipid interactions. For example, several cholesterol interaction sites have been proposed, such as the CRAC50 site, the CCM28 site, and more recently an aromatic signature 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 35

binding mode.97 The plasticity of the binding sites and the dynamics within a site or between neighboring sites increases the repertoire of cholesterol interaction sites. As suggested before,47 stochasticity at the cholesterol sites, as well as competition with other lipid molecules, could provide diversity in cholesterol occupancy sites. This diversity plays a critical role in fine tuning GPCR function in varying membrane microenvironments. However, the link between cholesterol binding and receptor activation has not yet been fully understood. The multiple interaction sites identified point toward a conformational modulation in both the ligand-bound and unbound states of the receptor. Increased receptor flexibility in the absence of cholesterol was previously reported by us by experimental21 and computational52 approaches, and also by Vattulainen and co-workers.42

In addition,

Weinstein/Khelashvili and co-workers earlier reported a link between conformational dynamics upon ligand binding and cholesterol dynamics associated with the serotonin2A receptor.76

Clearly, further studies are necessary to unravel the detailed molecular

mechanism underlying the link between activation of GPCRs and cholesterol interaction sites in GPCRs. In conclusion, the computational frameworks described here match several aspects of current experimental data, and make testable predictions.

As the spatial and temporal

resolution of experimental methods improve, it would be interesting to test the computational predictions to novel experimental data. In particular, mass spectrometric methods are being increasingly used to analyze protein-lipid interactions and could provide an improved resolution of these interactions.98 A closer match between experimental and computational studies would help uncover new aspects of GPCR-lipid interactions and provide a comprehensive understanding of receptor function.

14

ACS Paragon Plus Environment

Page 15 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Acknowledgments

D.S. and A.C. gratefully acknowledge the project (EMR/2016/002294) from the Science and Engineering Research Board (Govt. of India). A.C. gratefully acknowledges support from SERB Distinguished Fellowship (Department of Science and Technology, Govt. of India). A.C. is a Distinguished Visiting Professor at Indian Institute of Technology (Bombay), and Adjunct Professor at Tata Institute of Fundamental Research (Mumbai), RMIT University (Melbourne, Australia) and Indian Institute of Science Education and Research (Kolkata). We thank members of our research groups for critically reading the manuscript.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

References

1.

Pierce, K. L.; Premont, R. T.; Lefkowitz, R. J. Seven-transmembrane Receptors. Nat. Rev. Mol. Cell Biol. 2002, 3, 639-650.

2.

Rosenbaum, D. M.; Rasmussen, S. G. F.; Kobilka, B. K. The Structure and Function of G-protein-coupled Receptors. Nature 2009, 459, 356-363.

3.

Chattopadhyay, A. GPCRs: Lipid-dependent Membrane Receptors that Act as Drug Targets. Adv. Biol. 2014, 2014, 143023.

4.

Jacobson, K. A. New Paradigms in GPCR Drug Discovery. Biochem. Pharmacol. 2015, 98, 541-555.

5.

Cooke, R. M.; Brown, A. J. H.; Marshall, F. H.; Mason, J. S. Structures of G Proteincoupled Receptors Reveal New Opportunities for Drug Discovery. Drug Discov. Today 2015, 20, 1355-1364.

6.

Granier, S.; Kobilka, B. A New Era of GPCR Structural and Chemical Biology. Nat. Chem. Biol. 2012, 8, 670-673.

7.

Venkatakrishnan, A. J.; Deupi, X.; Lebon, G.; Tate, C. G.; Schertler, G. F.; Babu, M. M. Molecular Signatures of G-protein-coupled Receptors. Nature 2013, 494, 185-194.

8.

Latorraca, N. R.; Venkatakrishnan, A. J.; Dror, R. O. GPCR Dynamics: Structures in Motion. Chem. Rev. 2017, 117, 139-155.

9.

Nygaard, R.; Zou, Y.; Dror, R. O.; Mildorf, T. J.; Arlow, D. H.; Manglik, A.; Pan, A. C.; Liu, C. W.; Fung, J. J.; Bokoch, M. P. et al. The Dynamic Process of β2-Adrenergic Receptor Activation. Cell 2013, 152, 532-542.

10.

Deupi, X.; Kobilka, B. K. Energy Landscapes as a Tool to Integrate GPCR Structure, Dynamics, and Function. Physiology 2010, 25, 293-303.

11.

Unal, H.; Karnik, S. S. Domain Coupling in GPCRs: The Engine for Induced Conformational Changes. Trends Pharmacol. Sci. 2012, 33, 79-88.

12.

Cherezov, V.; Rosenbaum, D. M.; Hanson, M. A.; Rasmussen, S. G. F.; Thian, F. S.; Kobilka, T. S.; Choi, H.-J.; Kuhn, P.; Weis, W. I.; Kobilka, B. K. et al. High-resolution

16

ACS Paragon Plus Environment

Page 16 of 35

Page 17 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Crystal Structure of an Engineered Human β2-Adrenergic G Protein-coupled Receptor. Science 2007, 318, 1258-1265. 13.

Rosenbaum, D. M.; Cherezov, V.; Hanson, M. A.; Rasmussen, S. G. F.; Thian, F. S.; Kobilka, T. S.; Choi, H.-J.; Yao, X.-J.; Weis, W. I.; Stevens, R. C. et al. GPCR Engineering Yields High-resolution Structural Insights into β2-Adrenergic Receptor Function. Science 2007, 318, 1266-1273.

14.

Wacker, D.; Wang, S.; McCorvy, J. D.; Betz, R. M.; Venkatakrishnan, A. J.; Levit, A.; Lansu, K.; Schools, Z. L.; Che, T.; Nichols, D. E. et al. Crystal Structure of an LSDbound Human Serotonin Receptor. Cell 2017, 168, 377-389.

15.

Wheatley, M.; Wootten, D.; Conner, M. T.; Simms, J.; Kendrick, R.; Logan, R. T.; Poyner, D. R.; Barwell, J. Lifting the Lid on GPCRs: The Role of Extracellular Loops. Br. J. Pharmacol. 2012, 165, 1688-1703.

16.

Caffrey, M. A Comprehensive Review of the Lipid Cubic Phase or in meso Method for Crystallizing Membrane and Soluble Proteins and Complexes. Acta Crystallogr. F Struct. Biol. Commun. 2015, F71, 3-18.

17.

Khelashvili, G.; Albornoz, P. B. C.; Johner, N.; Mondal, S.; Caffrey, M.; Weinstein, H. Why GPCRs Behave Differently in Cubic and Lamellar Lipidic Mesophases. J. Am. Chem. Soc. 2012, 134, 15858-15868.

18.

Pucadyil, T. J.; Chattopadhyay, A. Role of Cholesterol in the Function and Organization of G-protein Coupled Receptors. Prog. Lipid Res. 2006, 45, 295-333.

19.

Paila, Y. D.; Chattopadhyay, A. Membrane Cholesterol in the Function and Organization of G-protein Coupled Receptors. Subcell. Biochem. 2010, 51, 439-466.

20. Pucadyil, T. J.; Chattopadhyay, A. Cholesterol Modulates Ligand Binding and G-protein Coupling to Serotonin1A Receptors from Bovine Hippocampus. Biochim. Biophys. Acta 2004, 1663, 188-200. 21. Saxena, R.; Chattopadhyay, A. Membrane Cholesterol Stabilizes the Human Serotonin1A Receptor. Biochim. Biophys. Acta 2012, 1818, 2936-2942.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22.

Paila, Y. D.; Kombrabail, M.; Krishnamoorthy, G.; Chattopadhyay, A. Oligomerization of the Serotonin1A Receptor in Live Cells: A Time-resolved Fluorescence Anisotropy Approach. J. Phys. Chem. B 2011, 115, 11439-11447.

23.

Ganguly, S.; Clayton, A. H. A.; Chattopadhyay, A. Organization of Higher-order Oligomers of the Serotonin1A Receptor Explored Utilizing Homo-FRET in Live Cells. Biophys. J. 2011, 100, 361-368.

24.

Chakraborty, H.; Jafurulla, M.; Clayton, A. H. A.; Chattopadhyay, A. Exploring Oligomeric State of the Serotonin1A Receptor utilizing Photobleaching Image Correlation Spectroscopy: Implications for Receptor Function. Faraday Discuss. (in press; DOI: 10.1039/C7FD00192D).

25.

Paila, Y. D.; Chattopadhyay, A. The Function of G-protein Coupled Receptors and Membrane Cholesterol: Specific or General Interaction? Glycoconj. J. 2009, 26, 711720.

26.

Paila, Y. D.; Jindal, E.; Goswami, S. K.; Chattopadhyay, A. Cholesterol Depletion Enhances Adrenergic Signaling in Cardiac Myocytes. Biochim. Biophys. Acta 2011, 1808, 461-465.

27.

Pydi, S. P.; Jafurulla, M.; Wai, L.; Bhullar, R. P.; Chelikani, P.; Chattopadhyay, A. Cholesterol Modulates Bitter Taste Receptor Function. Biochim. Biophys. Acta 2016, 1858, 2081-2087.

28.

Hanson, M. A.; Cherezov, V.; Griffith, M. T.; Roth, C. B.; Jaakola, V.-P.; Chien, E. Y. T.; Velasquez, J.; Kuhn, P.; Stevens, R. C. A Specific Cholesterol Binding Site is Established by the 2.8 Å Structure of the Human β2-Adrenergic Receptor. Structure 2008, 16, 897-905.

29.

Gater, D. L.; Saurel, O.; Iordanov, I.; Liu, W.; Cherezov, V.; Milon, A. Two Classes of Cholesterol Binding Sites for the β2AR Revealed by Thermostability and NMR. Biophys. J. 2014, 107, 2305-2312.

30.

Oates, J.; Watts, A. Uncovering the Intimate Relationship Between Lipids, Cholesterol and GPCR Activation. Curr. Opin. Struct. Biol. 2011, 21, 802-807.

18

ACS Paragon Plus Environment

Page 18 of 35

Page 19 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

31.

Jafurulla, M; Chattopadhyay, A. Membrane Lipids in the Function of Serotonin and Adrenergic Receptors. Curr. Med. Chem. 2013, 20, 47-55.

32.

Jafurulla, M.; Chattopadhyay, A. Sphingolipids in the Function of G Protein-coupled Receptors. Eur. J. Pharmacol. 2015, 763, 241-246.

33.

Dawaliby, R.; Trubbia, C.; Delporte, C.; Masureel, M.; Van Antwerpen, P.; Kobilka, B. K.; Govaerts, C. Allosteric Regulation of G Protein-coupled Receptor Activity by Phospholipids. Nat. Chem. Biol. 2016, 12, 35-39.

34. Ghosh, E.; Nidhi, K.; Shukla, A. K. SnapShot: GPCR-ligand Interactions. Cell 2014, 159, 1712-1712.e1. 35.

Sengupta, D.; Kumar, G. A.; Prasanna, X.; Chattopadhyay, A. Experimental and Computational Approaches to Study Membranes and Lipid-protein Interactions. In Computational Biophysics of Membrane Proteins; Domene, C., Ed.; Royal Society of Chemistry: London, UK, 2017; pp 137-160.

36.

Sengupta, D.; Chattopadhyay, A. Molecular Dynamics Simulations of GPCRcholesterol Interaction: An Emerging Paradigm. Biochim. Biophys. Acta 2015, 1848, 1775-1782.

37.

Ferro, A.; Coash, M.; Yamamoto, T.; Rob, J.; Ji, Y.; Queen, L. Nitric Oxide-dependent β2-Adrenergic Dilatation of Rat Aorta is Mediated through Activation of Both Protein Kinase A and Akt. Br. J. Pharmacol. 2004, 143, 397-403.

38.

Milano, C. A.; Allen, L. F.; Rockman, H. A.; Dolber, P. C.; McMinn, T. R.; Chien, K. R.; Johnson, T. D.; Bond, R. A.; Lefkowitz, R. Enhanced Myocardial Function in Transgenic Mice Overexpressing the β2-Adrenergic Receptor. Science 1994, 264, 582586.

39.

Liggett, S. B. The Pharmacogenetics of β2-Adrenergic Receptors: Relevance to Asthma. J. Allergy Clin. Immunol. 2000, 105, S487-S492.

40.

Prasanna, X.; Chattopadhyay, A.; Sengupta, D. Cholesterol Modulates the Dimer Interface of the β2-Adrenergic Receptor via Cholesterol Occupancy Sites. Biophys J. 2014, 106, 1290-1300.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

41.

Cang, X.; Du, Y.; Mao, Y.; Wang, Y.; Yang, H.; Jiang, H. Mapping the Functional Binding Sites of Cholesterol in β2-Adrenergic Receptor by Long-time Molecular Dynamics Simulations. J. Phys. Chem. B 2013, 117, 1085-1094.

42.

Manna, M.; Niemelä, M.; Tynkkynen, J.; Javanainen, M.; Kulig, W.; Müller, D. J.; Rog, T.; Vattulainen, I. Mechanism of Allosteric Regulation of β2-Adrenergic Receptor by Cholesterol. eLife 2016, 5, e18432.

43.

Paila, Y. D.; Tiwari, S.; Chattopadhyay, A. Are Specific Nonannular Cholesterol Binding Sites Present in G-protein Coupled Receptors? Biochim. Biophys. Acta 2009, 1788, 295-302.

44.

Pucadyil, T. J.; Kalipatnapu, S.; Chattopadhyay, A. The Serotonin1A Receptor: A Representative Member of the Serotonin Receptor Family. Cell. Mol. Neurobiol. 2005, 25, 553-580.

45.

Kalipatnapu, S.; Chattopadhyay, A. Membrane Organization and Function of the Serotonin1A Receptor. Cell. Mol. Neurobiol. 2007, 27, 1097-1116.

46.

Fiorino, F.; Severino, B.; Magli, E.; Ciano, A.; Caliendo, G.; Santagada, V.; Frecentese, F.; Perissutti, E. 5-HT1A Receptor: An Old Target as a New Attractive Tool in Drug Discovery from Central Nervous System to Cancer. J. Med. Chem. 2014, 57, 44074426.

47.

Sengupta, D.; Chattopadhyay, A. Identification of Cholesterol Binding Sites in the Serotonin1A Receptor. J. Phys. Chem. B 2012, 116, 12991-12996.

48.

Paila, Y. D.; Tiwari, S.; Sengupta, D.; Chattopadhyay, A. Molecular Modeling of the Human Serotonin1A Receptor: Role of Membrane Cholesterol in Ligand Binding of the Receptor. Mol. Biosyst. 2011, 7, 224-234.

49.

Fargin, A.; Raymond, J. R.; Lohse, M. J.; Kobilka, B. K.; Caron, M. G.; Lefkowitz, R. J. The Genomic Clone G-21 which Resembles a β-Adrenergic Receptor Sequence Encodes the 5-HT1A Receptor. Nature 1988, 335, 358-360.

50.

Jafurulla, M.; Tiwari, S.; Chattopadhyay, A. Identification of Cholesterol Recognition Amino Acid Consensus (CRAC) Motif in G-protein Coupled Receptors. Biochem. Biophys. Res. Commun. 2011, 404, 569-573. 20

ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

51.

Prasanna, X.; Sengupta, D.; Chattopadhyay, A. Cholesterol-dependent Conformational Plasticity in GPCR Dimers. Sci. Rep. 2016, 6, 31858.

52.

Patra, S. M.; Chakraborty, S.; Shahane, G.; Prasanna, X.; Sengupta, D.; Maiti, P. K.; Chattopadhyay, A. Differential Dynamics of the Serotonin1A Receptor in Membrane Bilayers of Varying Cholesterol Content Revealed by All Atom Molecular Dynamics Simulation. Mol. Membr. Biol. 2015, 32, 127-137.

53.

Mondal, S.; Khelashvili, G.; Johner, N.; Weinstein, H. How the Dynamic Properties and Functional Mechanisms of GPCRs are Modulated by Their Coupling to the Membrane Environment. Adv. Exp. Med. Biol. 2014, 796, 55-74.

54.

Mondal, S.; Johnston, J. M.; Wang, H.; Khelashvili, G.; Filizola, M.; Weinstein, H. Membrane Driven Spatial Organization of GPCRs. Sci. Rep. 2013, 3, 2909.

55.

Provasi, D.; Boz, M. B.; Johnston, J. M.; Filizola, M. Preferred Supramolecular Organization and Dimer Interfaces of Opioid Receptors from Simulated Selfassociation. PLoS Comput. Biol. 2015, 11, e1004148.

56.

Guixà-González, R.; Javanainen, M.; Gómez-Soler, M.; Cordobilla, B.; Domingo, J. C.; Sanz, F.; Pastor, M.; Ciruela, F.; Martinez-Seara, H.; Selent, J. Membrane Omega-3 Fatty Acids Modulate the Oligomerisation Kinetics of Adenosine A2A and Dopamine D2 Receptors. Sci. Rep. 2016, 6, 19839.

57.

Gahbauer, S.; Böckmann, R. A. Membrane-mediated Oligomerization of G Protein Coupled Receptors and its Implications for GPCR Function. Front. Physiol. 2016, 7, 494.

58.

Pluhackova, K.; Gahbauer, S.; Kranz, F.; Wassenaar, T. A.; Böckmann, R. A. Dynamic Cholesterol-conditioned Dimerization of the G Protein Coupled Chemokine Receptor Type 4. PLoS Comput Biol. 2016, 12, e1005169.

59.

Prasanna, X.; Chattopadhyay, A.; Sengupta, D. Role of Lipid-mediated Effects in β2Adrenergic Receptor Dimerization. Adv. Exp. Med. Biol. 2015, 842, 247-261.

60. Periole, X.; Huber, T.; Marrink, S.-J.; Sakmar, T. P. G Protein-coupled Receptors Self Assemble in Dynamics Simulations of Model Bilayers. J. Am. Chem. Soc. 2007, 129, 10126-10132. 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

61.

Liu, W.; Chun, E.; Thompason, A. A.; Chubukov, P.; Xu, F.; Katritch, V.; Han, G. W.; Roth, C. B.; Heitman, L. H.; IJzerman, A. P. et al. Structural Basis for Allosteric Regulation of GPCRs by Sodium Ions. Science 2012, 337, 232-236.

62.

Singh, P.; Paila, Y. D.; Chattopadhyay, A. Role of Glycosphingolipids in the Function of Human Serotonin1A Receptors. J. Neurochem. 2012, 123, 716-724.

63.

Jafurulla, M.; Pucadyil, T. J.; Chattopadhyay, A. Effect of Sphingomyelinase Treatment on Ligand Binding Activity of Human Serotonin1A Receptors. Biochim. Biophys. Acta 2008, 1778, 2022-2025.

64.

Paila, Y. D.; Ganguly, S.; Chattopadhyay, A. Metabolic Depletion of Sphingolipids Impairs Ligand Binding and Signaling of Human Serotonin1A Receptors. Biochemistry 2010, 49, 2389-2397.

65.

Singh, P.; Chattopadhyay, A. Removal of Sphingomyelin Headgroup Inhibits the Ligand Binding Function of Hippocampal Serotonin1A Receptors. Biochem. Biophys. Res. Commun. 2012, 419, 321-325.

66.

Jafurulla, M.; Bandari, S.; Pucadyil, T. J.; Chattopadhyay, A. Sphingolipids Modulate the Function of Human Serotonin1A Receptors: Insights From Sphingolipid-deficient Cells. Biochim. Biophys. Acta 2017, 1859, 598-604.

67.

Fantini, J. How Sphingolipids Bind and Shape Proteins: Molecular Basis of Lipidprotein Interactions in Lipid Shells, Rafts and Related Biomembrane Domains. Cell. Mol. Life Sci. 2003, 60, 1027-1032.

68.

Mahfoud, R.; Garmy, N.; Maresca, M.; Yahi, N.; Puigserver, A.; Fantini, J. Identification of a Common Sphingolipid-binding Domain in Alzheimer, Prion, and HIV-1 Proteins. J. Biol. Chem. 2002, 277, 11292-11296.

69.

Chattopadhyay, A.; Paila, Y. D.; Shrivastava, S.; Tiwari, S.; Singh, P.; Fantini, J. Sphingolipid-binding Domain in the Serotonin1A Receptor. Adv. Exp. Med. Biol. 2012, 749, 279-293.

70.

Shrivastava, S.; Jafurulla, M.; Tiwari, S.; Chattopadhyay, A. Identification of Sphingolipid-binding Motif in G Protein-coupled Receptors. Adv. Exp. Med. Biol. (in press). 22

ACS Paragon Plus Environment

Page 22 of 35

Page 23 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

71.

Prasanna, X.; Jafurulla, M.; Sengupta, D.; Chattopadhyay, A. The Ganglioside GM1 Interacts with the Serotonin1A Receptor via the Sphingolipid Binding Domain. Biochim. Biophys. Acta 2016, 1858, 2818-2826.

72.

Mahfoud, R.; Mylvaganam, M.; Lingwood, C. A.; Fantini, J. A Novel Soluble Analog of the HIV-1 Fusion Cofactor, Globotriaosylceramide (Gb3), Eliminates the Cholesterol Requirement for High Affinity gp120/Gb3 Interaction. J. Lipid Res. 2002, 43, 1670-1679.

73.

Ling, H.; Boodhoo, A.; Hazes, B.; Cummings, M. D.; Armstrong, G. D.; Brunton, J. L.; Read, R. J. Structure of the Shiga-Like Toxin I B-pentamer Complexed with an Analogue of its Receptor Gb3. Biochemistry 1998, 37, 1777-1788.

74.

Grossfield, A. Recent Progress in the Study of G Protein-coupled Receptors with Molecular Dynamics Computer Simulations. Biochim. Biophys. Acta 2011, 1808, 18681878.

75.

Khelashvili, G.; Grossfield, A.; Feller, S. E.; Pitman, M. C.; Weinstein, H. Structural and Dynamic Effects of Cholesterol at Preferred Sites of Interaction with Rhodopsin Identified From Microsecond Length Molecular Dynamics Simulations. Proteins 2009, 76, 403–417.

76.

Shan, J.; Khelashvili, G.; Mondal, S.; Mehler, E. L.; Weinstein, H. Ligand-dependent Conformations and Dynamics of the Serotonin 5-HT2A Receptor Determine its Activation and Membrane-driven Oligomerization Properties. PLoS Comput. Biol. 2012, 8, e1002473.

77.

Rouviere, E.; Arnarez, C.; Yang, L.; Lyman, E. Identification of Two New Cholesterol Interaction Sites on the A2A Adenosine Receptor. Biophys. J. 2017, 113, 2415-2424..

78.

Lee, J. Y.; Lyman, E. Predictions for Cholesterol Interaction Sites on the A2A Adenosine Receptor. J. Am. Chem. Soc. 2012, 134, 16512-16515.

79.

Guixà-González, R.; Albasanz, J. L.; Rodriguez-Espigares, I.; Pastor, M.; Sanz, F.; Martí-Solano, M.; Manna, M.; Martinez-Seara, H.; Hildebrand, P. W.; Martín, M. et al. Membrane Cholesterol Access into a G-protein-coupled Receptor. Nat. Commun. 2017, 8, 14505. 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

80.

Kharche, S.; Joshi, M.; Sengupta, D.; Chattopadhyay, A. Membrane-induced Organization and Dynamics of the N-terminal Domain of Chemokine Receptor CXCR1: Insights from Atomistic Simulations. Chem. Phys. Lipids 2018, 210, 142-148.

81.

Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; de Vries, A. H. The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. B 2007, 111, 7812-7824.

82.

Monticelli, L.; Kandasamy, S. K.; Periole, X.; Larson, R. G.; Tieleman, D. P.; Marrink, S. J. The MARTINI Coarse-grained Force Field: Extension to Proteins. J. Chem. Theory Comput. 2008, 4, 819-834.

83.

Marrink, S. J.; Tieleman, D. P. Perspective on the Martini Model. Chem. Soc. Rev. 2013, 42, 6801-6822.

84.

Periole, X.; Cavalli, M.; Marrink, S. J.; Ceruso, M. A. Combining an Elastic Network with a Coarse-grained Molecular Force Field: Structure, Dynamics, and Intermolecular Recognition. J. Chem. Theory Comput. 2009, 5, 2531-2543.

85.

Schäfer, L. V.; de Jong, D. H.; Holt, A.; Rzepiela, A. J.; de Vries, A. H.; Poolman, B.; Killian, J. A.; Marrink, S. J. Lipid Packing Drives the Segregation of Transmembrane Helices into Disordered Lipid Domains in Model Membranes. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 1343-1348.

86.

Baaden, M.; Marrink, S. J. Coarse-grain Modelling of Protein–protein Interactions. Curr. Opin. Struct. Biol. 2013, 23, 878-886.

87.

Louhivuori, M.; Risselada, H. J.; van der Giessen, E.; Marrink, S. J. Release of Content Through Mechano-sensitive Gates in Pressurized Liposomes. Proc. Nat. Acad. Sci. U. S. A. 2010, 107, 19856-19860.

88.

Sengupta, D.; Rampioni, D.; Marrink, S. J. Simulations of the c-Subunit of ATPSynthase Reveal Helix Rearrangements. Mol. Membr. Biol. 2009, 26, 422-434.

89.

Horn, J. N.; Kao, T.-C.; Grossfield, A. Coarse-grained Molecular Dynamics Provides Insight into the Interactions of Lipids and Cholesterol with Rhodopsin. Adv. Exp. Med. Biol. 2014, 796, 75-94.

24

ACS Paragon Plus Environment

Page 24 of 35

Page 25 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

90.

Huster, D. Solid-state NMR Spectroscopy to Study Protein-lipid Interactions. Biochim. Biophys. Acta 2014, 1841, 1146-1160.

91.

Rasmussen, S. G.; DeVree, B. T.; Zou, Y.; Kruse, A. C.; Chung, K. Y.; Kobilka, T. S.; Thian, F. S.; Chae, P. S.; Pardon, E.; Calinski, D. et al. Crystal Structure of the β2 Adrenergic Receptor-Gs Protein Complex. Nature 2011, 477, 549-555.

92.

Wassenaar, T. A.; Pluhackova, K.; Moussatova, A.; Sengupta, D.; Marrink, S. J.; Tieleman, D. P.; Bockmann, R. A. High-throughput Simulations of Dimer and Trimer Assembly of Membrane Proteins. The DAFT Approach. J. Chem. Theory Comput. 2015, 11, 2278-2291.

93.

Sengupta, D.; Marrink, S. J. Lipid-mediated Interactions Tune the Association of Glycophorin A Helix and its Disruptive Mutants in Membranes. Phys. Chem. Chem. Phys. 2010, 12, 12987-12996.

94.

Prasanna, X.; Praveen, P. J.; Sengupta, D. Sequence Dependent Lipid-mediated Effects Modulate the Dimerization of ErbB2 and its Associative Mutants. Phys. Chem. Chem. Phys. 2013, 15, 19031-19041.

95.

Pawar, A.; Deshpande, S.; Gopal, S. M.; Wassenaar, T. A.; Athale, C. A.; Sengupta, D. Thermodynamic and Kinetic Characterization of Transmembrane Helix Association. Phys. Chem. Chem. Phys. 2015, 17, 1390-1398.

96.

Altwaijry, N. A.; Baron, M.; Wright, D. W.; Coveney, P. V.; Townsend-Nicholson, A. An Ensemble-based Protocol for the Computational Prediction of Helix-helix Interactions in G Protein-coupled Receptors Using Coarse-grained Molecular Dynamics. J. Chem. Theory Comput. 2017, 13, 2254-2270.

97.

Mohole, M.; Prasanna, X.; Sengupta, D.; Chattopadhyay, A. Molecular Signatures of Cholesterol Interaction with Serotonin Receptors. Adv. Exp. Med. Biol. (in press).

98.

Gupta, K.; Donlan, J. A.; Hopper, J. T.; Uzdavinys, P.; Landreh, M.; Struwe, W. B.; Drew, D.; Baldwin, A. J.; Stansfeld, P. J.; Robinson, C.V. The Role of Interfacial Lipids in Stabilizing Membrane Protein Oligomers. Nature 2017, 541, 421-424.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure Legends

Figure 1. The β2-adrenergic receptor and the CCM cholesterol interaction site. (a) A schematic representation of the β2-adrenergic receptor with each helix colored individually. The high occupancy cholesterol site on transmembrane helix IV corresponding to the CCM site (b) resolved crystallographically, and (c) predicted from coarse-grain simulations. (d-g) Dynamics and translation of the cholesterol molecule at transmembrane helix IV sampled during simulations. The transmembrane helix is colored blue, the highlighted side chains are colored white, the bound cholesterol molecule is shown in magenta with the polar head group (-OH) shown in cyan. Reproduced with permission from Ref. 40.

Figure 2. Cholesterol interaction sites of the β2-adrenergic receptor. The high density cholesterol sites predicted from coarse-grain simulations on the (a) upper leaflet, (b) lower leaflet and (c) center of the bilayer marked as e1-3, i2-5 and m1, respectively. The top view of the receptor is superimposed on the plot for clarity.

Figure 3. Cholesterol occupancy sites at the transmembrane helices of the serotonin1A receptor. The high density cholesterol sites predicted from coarse-grain simulations on the (a) upper leaflet, and (b) lower leaflet are shown.

The top view of the receptor is

superimposed on the plot for clarity.

Figure 4. Phospholipid interaction sites in GPCRs. The main interaction site at the groove of transmembrane helices I and VII as observed in (a) crystal structure of A2A adenosine receptor (b) coarse-grain simulations of β2-adrenergic receptor and (c) atomistic simulations of the serotonin1A receptor. The transmembrane helices I and VII are depicted in blue and green respectively. The lipid molecule is shown in gray and the headgroups are colored blue and red.

26

ACS Paragon Plus Environment

Page 26 of 35

Page 27 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5.

Association of GM1 with the serotonin1A receptor. Contact maps of the

extracellular loop 1 (residues 97-109) of the serotonin1A receptor with the GM1 molecule in (a) presence and (b) absence of cholesterol. The GM1 molecule was represented by the headgroup beads (HG1-HG17) and the sphingosine backbone beads (AM1 and AM2). Reproduced with permission from Ref. 71. (c) The coarse-grain representation of GM1 mapped to the underlying structure.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic

28

ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Author Biographies

Dr. Durba Sengupta is a faculty member at the National Chemical Laboratory in Pune, India. She received her B.Sc. degree in Chemistry from St. Stephen's College, University of Delhi and M.Sc. in Biotechnology from Indian Institute of Technology-Bombay, Mumbai. Subsequently, she received her Ph.D. from University of Heidelberg and was a Postdoctoral Fellow at University of Groningen. The central objective of her research is to understand self organization and self assembly in cell membranes using multi-scale simulation approaches. Currently, the main focus of her research group is to analyze functional dynamics of membrane receptors, in particular GPCRs.

Dr. Xavier Prasanna received his B.Sc. and M.Sc. (Microbiology) from St. Joseph's College, Bangalore. He received his Ph.D. from National Chemical Laboratory, Pune. His research interests are in using multiscale simulations to examine the effect of membrane components on higher order organization of G protein-coupled receptors.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Madhura Mohole received her B.Sc. from Department of Biotechnology, M.E.S. Abasaheb Garware College, Pune, and M.Sc in Bioinformatics from Bioinformatics Centre, University of Pune. She is currently pursuing her Ph.D. at National Chemical Laboratory, Pune. Her research includes understanding GPCR structure, organization and function using computational methods.

Prof. Amitabha Chattopadhyay received B.Sc. with Honors in Chemistry from St. Xavier’s College (Calcutta) and M.Sc. in Chemistry from Indian Institute of Technology, Kanpur. He obtained his Ph.D. from the State University of New York at Stony Brook, and was a Postdoctoral Fellow at the University of California, Davis. He subsequently joined the Centre for Cellular and Molecular Biology in Hyderabad, where he is currently a SERB Distinguished Fellow. Prof. Chattopadhyay’s work is focused on the role of membrane lipids on the function of G protein-coupled receptors and its implications in health and disease using experimental and simulation approaches. A translational extension of this work has been on the role of host membrane lipids on the entry of intracellular pathogens into host cells. In addition, his group pioneered the development and application of the wavelengthselective fluorescence approach as a novel tool to monitor organization and dynamics of probes and proteins in membranes and micelles. Prof. Chattopadhyay was awarded the prestigious TWAS (The World Academy of Sciences) Prize, Shanti Swarup Bhatnagar Award and Ranbaxy Research Award. He is an elected Fellow of TWAS, the Royal Society of Biology, Royal Society of Chemistry, and all the Indian Academies of Science. He has served on the editorial board of a large number of reputed journals.

30

ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1 Sengupta et al.

(a)

(b)

(c)

W158 I154

W158 I154 R151

R151 (d)

(e)

(f)

ACS Paragon Plus Environment

(g)

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 35

Figure 2 Sengupta et al. 


(a)

(b)

ACS Paragon Plus Environment

(c)

Page 33 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3 Sengupta et al.

(a)

(b)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 35

Figure 4 Sengupta et al.

(a)

(b)

ACS Paragon Plus Environment

(c)

Page 35 of 35

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5 Sengupta et al.

ACS Paragon Plus Environment