Mechanism of allosteric regulation of β2-adrenergic receptor by cholesterol

  1. Moutusi Manna
  2. Miia Niemelä
  3. Joona Tynkkynen
  4. Matti Javanainen
  5. Waldemar Kulig
  6. Daniel J Müller
  7. Tomasz Rog
  8. Ilpo Vattulainen  Is a corresponding author
  1. Tampere University of Technology, Finland
  2. University of Helsinki, Finland
  3. ETH Zürich, Switzerland
  4. University of Southern Denmark, Denmark
10 figures, 2 videos and 2 tables

Figures

Figure 1 with 1 supplement
Conformational dynamics of β2AR.

(A) The distances between the Cα atoms of D1133.32–S2075.46 (distance defined as LL) and R1313.50–E2686.30 (LG) pairs used to measure the fluctuations at the ligand and G-protein binding sites, respectively. (B–C) The conformational distributions of β2AR in membranes with 0 and 10 mol% cholesterol (Chol) as a function of LL and LG. The gray dotted lines represent the corresponding LL and LG values in the inactive crystal structure of β2AR (Hanson et al., 2008). The cartoon diagram shows the fluctuations of LL and LG at the ligand and G-protein binding sites of the receptor, respectively. (D–E) The time evolution of LL (light red) and LG (light blue) in systems with 0 and 10 mol% cholesterol. Corresponding 50-point running averages are shown in dark colors. (F) Standard deviation for the distribution of the distance between the intracellular (IC) (or extracellular (EC)) end of H5 and its average position, and its dependence on whether the given end of H5 is in contact with cholesterol or not; similarly for H6.

https://doi.org/10.7554/eLife.18432.004
Figure 1—figure supplement 1
Conformational distributions of β2AR in lipid bilayers with various cholesterol (Chol) concentrations.

In panels (A–F) the distributions are plotted as a function of LL (distance between the Cα atoms of D1133.32 and S2075.46) at the ligand binding site and LG (distance between the Cα atoms of R1313.50 and E2686.30) at the G protein-binding site. (A–E) Starting from situations where no cholesterol molecules were initially bound to β2AR, distributions are plotted over all independent trajectories of a given system, where the equilibration time (the first 100 ns) was discarded from the analysis. (F) β2AR conformational distribution in control simulations, where cholesterol molecules were initially bound at the eight interaction sites of β2AR predicted by our simulations, but no further cholesterol was in the membrane (total (average) cholesterol concentration 1.9 mol%). Here, there is reason to keep in mind the rapid migration of cholesterols away from the receptor surface when the cholesterol concentration is low (see main text and Figure 5), implying that panel (F) corresponds to cholesterol-rich conditions in the vicinity of the receptor at very short times but to cholesterol-poor conditions at long times.

https://doi.org/10.7554/eLife.18432.005
Figure 2 with 10 supplements
Cholesterol interaction sites on β2AR.

(A–B) 2D number densities of cholesterol (Chol) around β2AR. The data are averaged over all independent trajectories for a given cholesterol concentration (Table 1) and normalized with respect to the maximum density for that particular cholesterol concentration. The intracellular (IC) and extracellular (EC) bilayer leaflets are depicted separately. The major cholesterol interaction sites (IC1, IC2 and EC1) are marked in the density plots. The IC and EC sides of the transmembrane regions (H1–H7) of β2AR are shown in gray scale (the darker the color, the higher is the number density) and numbered accordingly. (C–D) Cartoon representation of three main cholesterol interaction sites in β2AR. IC1 (H1–H4) and IC2 (H5–H6) are located on the intracellular side, and EC1 comprised of two closely placed cholesterols between H5-H6 and H6-ECL3-H7 is located on the extracellular side of β2AR.

https://doi.org/10.7554/eLife.18432.006
Figure 2—figure supplement 1
Residues of β2AR involved in cholesterol binding, and cholesterol interaction sites on β2AR.

Panels (A–B) (top): Cholesterol occupancy time per residue of β2AR described in terms of the normalized time fraction, where a value of one stands for a contact throughout the simulation trajectory and zero means no contact. Results are given for (A) 10 and (B) 40 mol% of cholesterol. The residues of β2AR are defined to be in contact with cholesterol when any non-hydrogen atom of the residue is within ≤0.5 nm of any heavy atom of cholesterol. The data show that there are several hot spots (blue) as cholesterol binding sites. These plots were averaged from all independent simulations for a given cholesterol concentration, where the equilibration time (the first 100 ns of the simulation) was disregarded from the analysis. Panels (C) (bottom): Interaction sites as obtained from our simulations, are shown from two perspectives around the protein. EC and IC stand for extracellular and intracellular, respectively. Interaction sites at the intracellular (IC) side: IC1 (dark green) between helices (H) 1–4, IC2 (red) between H5 and H6, IC3 (magenta) between H3 and H5, and IC4 (orange) between H1 and H8. Interaction sites at the extracellular (EC) side: EC1 comprised of two closely placed cholesterol molecules between H5 and H6 (cyan) and in space surrounded by H6-ECL3-H7 (green), EC2 (purple) between H3 and H4, and EC3 (blue) between H1-H2-ECL1 (where ECL stands for the extracellular loop).

https://doi.org/10.7554/eLife.18432.007
Figure 2—figure supplement 2
Sequence alignment of β2AR orthologues around the cholesterol-binding site IC1.

The residues that play a major role (contact fraction ≥ 0.4, where one stands for maximum contact and zero for no contact) in cholesterol binding are highlighted. Here for IC1, the residues in the cholesterol consensus motif are highlighted in red. Following sequence alignment, shown are the contact fraction per residue (tan bars) and its occurrence in the set of sequences (cyan bars) [one stands for 100% and zero for no occurrence]. The occurrence represents the extent a particular residue is conserved.

https://doi.org/10.7554/eLife.18432.008
Figure 2—figure supplement 3
Sequence alignment of β2AR orthologues around the cholesterol-binding site IC2.

The residues that play a major role (contact fraction ≥ 0.4, where one stands for maximum contact and zero for no contact) in cholesterol binding are highlighted. Following sequence alignment, shown are the contact fraction per residue (tan bars) and its occurrence in the set of sequences (cyan bars) [one stands for 100% and zero for no occurrence]. The occurrence represents the extent a particular residue is conserved.

https://doi.org/10.7554/eLife.18432.009
Figure 2—figure supplement 4
Sequence alignment of β2AR orthologues around the cholesterol-binding site EC1.

The residues that play a major role (contact fraction ≥ 0.4, where one stands for maximum contact and zero for no contact) in cholesterol binding are highlighted. Following sequence alignment, shown are the contact fraction per residue (tan bars) and its occurrence in the set of sequences (cyan bars) [one stands for 100% and zero for no occurrence]. The occurrence represents the extent a particular residue is conserved.

https://doi.org/10.7554/eLife.18432.010
Figure 2—figure supplement 5
Cholesterol density around the receptor at low cholesterol concentrations.

Two-dimensional (2D) averaged and normalized number densities of cholesterol around β2AR shown at low cholesterol concentrations (2 and 5 mol%). The intracellular and extracellular leaflets are depicted separately. The intracellular and extracellular sides of the transmembrane regions of β2AR are shown in gray scale (the darker the color, the higher is the number density), and they are numbered accordingly to show the locations of the individual helices (H1–H7).

https://doi.org/10.7554/eLife.18432.011
Figure 2—figure supplement 6
Structure of cholesterol analogues and properties of sterol-containing bilayers.

(A) The different cholesterol analogues used in the current study. (B–D) Average lipid chain order parameter SCD of DOPC bilayers with different concentrations of cholesterol or cholesterol-analogues. (E–G) Average bilayer thickness in DOPC bilayers with different concentrations of cholesterol or cholesterol-analogues. Error bars for order parameter and thickness are less than 0.02 and 0.005 Å, respectively.

https://doi.org/10.7554/eLife.18432.012
Figure 2—figure supplement 7
Interactions of cholesterol and cholesterol-like molecules with β2AR.

The average interaction energies for van der Waals (vdW) and electrostatic interactions are determined separately. Error bars are in the range of 0.1–1 kJ/mol. The lower panel represents the oxysterol-containing systems, where a fraction of cholesterol is replaced by 4β-OH-Chol and 27-OH-Chol, respectively, resulting in bilayers with 4 mol% oxysterol and 21 mol% cholesterol.

https://doi.org/10.7554/eLife.18432.013
Figure 2—figure supplement 8
Densities of sterols around β2AR.

Normalized 2D average number densities around β2AR: (A–B) CHSA (the deprotonated form of cholesteryl hemisuccinate (CHS)); (C–F) CHS. Densities of sterols in mixed sterol-containing bilayers with other molecules: (G–H) 4β-hydroxy-Chol (4β-OH-Chol); (K–L) 27-hydroxy-Chol (27-OH-Chol). The densities of 4β-OH-Chol and 27-OH-Chol are shown separately: (I–J) 4β-OH-Chol; (M–N) 27-OH-Chol. For descriptions of models, see Table 1. For each system, the intracellular and extracellular bilayer leaflets are depicted separately. The intracellular and extracellular sides of β2AR transmembrane regions (H1–H7) are shown in gray scale (the darker the color, the higher the number density) and numbered accordingly.

https://doi.org/10.7554/eLife.18432.014
Figure 2—figure supplement 9
Conformational distributions of β2AR in lipid bilayers with different cholesterol analogues.

(A–B) Oxysterol-containing systems having 4 mol% of oxysterol (27-OH-Chol or 4β-OH-Chol) and 21% cholesterol. (C–D) DOPC bilayer with 10 mol% and 40 mol% of CHS. Conformational distributions are calculated over all independent trajectories of a given system, where the equilibration time (100 ns) is disregarded from the analysis. The CαD1133.32–CαS2075.46 (defined as LL) and CαR1313.50–CαE2686.30 (LG) distances represent the fluctuations in the ligand and the G protein-binding sites, respectively.

https://doi.org/10.7554/eLife.18432.015
Figure 2—figure supplement 10
IC1 interaction site.

Specific cholesterol binding site in β2AR with the cholesterol consensus motif displayed with side chain positions of the conserved amino acid residues, as found in (A) the crystal structure (ref. 17) and (B) during our simulation. In the simulation snapshot, residues are colored according to their strength of interaction with cholesterol (red represents the weakest and blue represents the strongest interaction). (C) As to the time-dependent distance between H4 and its average position, as the H4 helix fluctuates around its average location, shown here are results for the standard deviation of the distance fluctuations. Data are given for cases, where IC1 is occupied (blue) or unoccupied (orange) by cholesterol.

https://doi.org/10.7554/eLife.18432.016
Figure 3 with 1 supplement
Effect of cholesterol on the active conformation of β2AR.

Cytosolic view of β2AR (A) in the beginning of a simulation (active state) as well as in representative simulation snapshots in (B) a DOPC bilayer and (C) in the presence of 40 mol% cholesterol. The dotted line represents the distance between the Cα atoms of R1313.50–E2686.30 (defined as LG), used to measure the fluctuation at the G protein-binding site. (D) Simulation snapshot (in the presence of 40 mol% cholesterol) showing cholesterol binding at the interaction sites of β2AR. (E) The time evolution of LG in systems with 0 (light red) and 40 mol% cholesterol (light blue). Corresponding 50-point running averages are shown in dark colors (red, blue). (F) 2D number densities of cholesterol around β2AR (cytosolic view). The individual transmembrane helixes of β2AR are numbered and shown in gray scale (as in Figure 2A,B).

https://doi.org/10.7554/eLife.18432.017
Figure 3—figure supplement 1
Conformational distribution of β2AR starting from the active state.

The conformational distributions of β2AR in (left) a DOPC bilayer and (right) a DOPC bilayer with 40 mol% cholesterol (Chol) as a function of LL and LG. The gray dotted lines represent the corresponding LL and LG values in the initial active crystal structure of β2AR. The distribution is averaged over the different replicas of each system.

https://doi.org/10.7554/eLife.18432.018
Figure 4 with 1 supplement
Impact of membrane-mediated effects on the β2AR conformation.

The conformational distribution of β2AR in bilayers composed of (A) long-chain PC-20:0/22:1 c13 lipids and (C) DOPC with 20 mol% pyrene (Pyrene20). (B) 3D-distribution of bilayer thickness in the thicker PC-20:0/22:1 c13 membrane. The receptor is shown as a purple cartoon. (D) 2D number density of pyrene around β2AR.

https://doi.org/10.7554/eLife.18432.019
Figure 4—figure supplement 1
Properties of thick and/or ordered cholesterol-free bilayers.

(A–B) Long-chain PC bilayer properties compared to those of cholesterol-rich and DOPC systems. (A) The average bilayer thickness in several different bilayer systems (see Table 1). (B) The average lipid chain order parameter of a bilayer composed of long-chain PC-20:0/22:1 c13 lipids. Error bars for thickness and order parameter are less than 0.005 Å and 0.02, respectively. Panels (C–F) show results of pyrene-containing bilayer systems. (C) The structure of the polycyclic aromatic hydrocarbon compound pyrene. Bilayer properties of the pyrene-containing membrane: (D) lipid chain order parameter and (E) the average bilayer thickness, in comparison with other cholesterol-containing bilayers. Error bars for thickness are less than 0.005 Å. (F) The probability of finding pyrene at the β2AR surface. The corresponding fit based on exponential decay is shown in red.

https://doi.org/10.7554/eLife.18432.020
Figure 5 with 1 supplement
Binding time of cholesterol.

(A–C) Time-correlation function of cholesterol (Chol) at the three major interaction sites (IC1, IC2, EC1) on the β2AR surface. Initially cholesterol is bound to the site (distance ≤ 0.5 nm) and the correlation function describes the probability that cholesterol remains bound to the given site for increasing time. Data are shown for DOPC-cholesterol membranes with 2, 5, 10, and 40 mol% of cholesterol. (D–E) Schematic representation showing the transition from fast to slow exchange as cholesterol concentration increases from 2 to 40 mol%. Color code: β2AR (blue), DOPC (thin grey lines), cholesterol molecules bound to the interaction sites (purple), and other cholesterol molecules not bound to the receptor (yellow sticks).

https://doi.org/10.7554/eLife.18432.021
Figure 5—figure supplement 1
Interaction of cholesterol with β2AR.

Time development for the distances of cholesterol molecules from the β2AR surface, where these cholesterol molecules were initially bound at the eight binding sites identified in this study (cholesterol-bound, see Table 1; Figure 2—figure supplement 1). Here, EC1-A and EC1-B stand for the two cholesterol molecules in the EC1 binding site. The rest of the membrane was initially cholesterol-free. Shown here are the data based on the three independent repeats.

https://doi.org/10.7554/eLife.18432.022
Author response image 1
Specific cholesterol binding site in β2AR with CCM displayed with side chain positions of conserved amino acid residues, as found in (A) the crystal structure (1) and (B) during our simulation.

In the simulation snapshot, residues are colored according to their strength of interaction with cholesterol (red represents the weakest and blue represents the strongest interaction).

https://doi.org/10.7554/eLife.18432.026
Author response image 2
For the time-dependent distance betweenH4 and its average position, as the H4 helix fluctuates around its average location, shown here are results for the standard deviation of the distance fluctuations.

Data are given for cases, where IC1 is occupied (blue) or unoccupied (orange) by cholesterol.

https://doi.org/10.7554/eLife.18432.027
Author response image 3
Distributions of LL and LG distances from individual trajectories (shown in different colors) for various cholesterol concentrations.
https://doi.org/10.7554/eLife.18432.028
Author response image 4
Area in the 2D histogram visited by the receptor conformations.

The bin edge length was set to 0.1 Å in both dimensions.

https://doi.org/10.7554/eLife.18432.029
Author response image 5
Cytosolic view of β2AR (A) in the beginning of simulation (active state) as well as in representative simulation snapshots in (B) a DOPC bilayer and (C) in the presence of 40 mol% cholesterol.

The dotted line represents the distance between the Cα atoms of R1313.50–E2686.30 (defined as LG) used to measure the fluctuation at the G protein binding site. (D) Simulation snapshot (in the presence of 40 mol% cholesterol) showing cholesterol binding at the interaction sites of β2AR. (E) The time evolution of LG in systems with 0 (light red) and 40 mol% cholesterol (light blue). Corresponding 50-point running averages are shown in dark colors. (F) 2D number densities of cholesterol (Chol) around β2AR (cytosolic view). The transmembrane regions of β2AR are shown in gray scale (the darker the color, the higher is the number density), and they are numbered accordingly to show the locations of the individual helices (H1-H7).

https://doi.org/10.7554/eLife.18432.030

Videos

Video 1
Spontaneous binding/unbinding of cholesterol at the three main cholesterol interaction sites of β2AR during a 2.5-μs simulation with 10 mol% of cholesterol.

Cholesterols interacting at the cholesterol-binding sites are highlighted (yellow at IC1; green at IC2; and blue and red at EC1). Other cholesterols are shown in gray. For clarity, other lipids in a membrane are not shown.

https://doi.org/10.7554/eLife.18432.023
Video 2
Spontaneous binding/unbinding of cholesterol at the three main cholesterol interaction sites of β2AR during a 2.5-μs simulation with 40 mol% of cholesterol.

Cholesterols interacting at the cholesterol-binding interaction sites are highlighted (yellow and green at IC1; red, blue and orange at IC2; and pink, purple and cyan at EC1). Other cholesterols are shown in gray. For clarity, other lipids in a membrane are not shown.

https://doi.org/10.7554/eLife.18432.024

Tables

Table 1

Descriptions of systems simulated: β2AR in bilayers with varying lipid compositions. ‘Chol’ stands for cholesterol.

https://doi.org/10.7554/eLife.18432.003

Systems*

Initial lipid arrangement around β2AR

Lipids

Sterol mol %

No. of repeats

Time (μs)

DOPC

Random

DOPC

0

3

3×2.5

DOPC-active

Random

DOPC

0

3

3×2.5

C

H

O

L

Chol2

Random

DOPC + Chol

2

3

3×2.5

R

A

N

D

O

M

Chol5

Random

DOPC + Chol

5

3

3×2.5

Chol10

Random

DOPC + Chol

10

3

3×2.5

Chol25

Random

DOPC + Chol

25

2

2×2

Chol40

Random

DOPC + Chol

40

3

3×2.5

Chol40-active

Random

DOPC + Chol

40

3

3×2.5

C

H

S

CHS10

Random

DOPC + CHS

10

2

2×2

CHS40

Random

DOPC + CHS

40

2

2×2

CHSA10

[A for anionic]

Random

DOPC + CHSA

10

1

2

CHSA40

Random

DOPC + CHSA

40

1

2

O

X

Y

S

T

E

R

O

L

27-OH-Chol

Random

[16 mol % Chol was randomly replaced by 27-OH-Chol]

DOPC +

Chol +

27-OH-Chol

25

(4 mol% 27-OH-Chol + 21 mol% Chol)

3

2 + 1 + 1

4β-Chol

Random

[16 mol% Chol was randomly replaced by 4β-OH-Chol]

DOPC +

Chol +

4β-OH-Chol

25

(4 mol% 4β-OH-Chol + 21 mol% Chol)

3

1 + 1 + 1

Chol-Bound§

8 cholesterols bound at sites predicted by simulations

DOPC + Chol

1.9

3

3×2.5

B
O
U
N
D

Chol-IC1

2 Chol bound at IC1

DOPC + Chol

<1

2

2×2

CHS-IC1

2 CHS bound at IC1

DOPC + CHS

<1

1

2

CHSA-IC1

2 CHSA bound at IC1

DOPC + CHSA

<1

1

2

PC-20:0–22:1 c13

[Double bond at carbon 13]

Random

PC-20:0–22:1 c13

0

3

3×1.5

Pyrene20

Random

DOPC +

20 mol% pyrene

0

3

3×1.5

  1. *In the DOPC-active and Chol40-active systems, we used the active-state conformation of the receptor as the starting structure; for all the other systems, we used the inactive conformation.

  2. For systems with no sterols initially bound to β2AR, i.e., the systems which started with a random distribution of lipids, a number of different repeat simulations for each lipid composition were performed with different initial lipid arrangements around the receptor. For systems with sterols initially bound to β2AR (seed and BOUND), different replicas were generated with different starting velocities.

  3. Listed are the simulation times of production simulations; the equilibration time of the systems (100 ns) is not included.

  4. §In the Chol-Bound system, eight cholesterol molecules were initially (at time zero of the simulation) bound at eight binding sites predicted by the present simulations, while the rest of the system had no cholesterol at all.

Table 2

Interactions* of sterols at the three high-affinity cholesterol-binding sites.

https://doi.org/10.7554/eLife.18432.025

Cholesterol/Cholesterol analogue

High-affinity cholesterol interaction sites

IC1

IC2

EC1

vdW interaction energy (kJ/mol)

No. of contacts

vdW interaction energy (kJ/mol)

No. of contacts

vdW interaction energy (kJ/mol)

No. of contacts

Cholesterol

−138.04 ± 0.20

141.02 ± 0.22

−95.06 ± 0.12

90.65 ± 0.16

−129.51 ± 0.29

104.38 ± 0.28

CHS

−29.63 ± 0.14

28.78 ± 0.16

−98.75 ± 0.11

96.30 ± 0.16

-

-

27-OH-Chol

−32.17 ± 0.30

34.95 ± 0.33

−22.69 ± 0.23

28.41 ± 0.28

−132.85 ± 0.27

120.20 ± 0.30

4β-OH-Chol

-

-

-

-

−41.80 ± 0.48

33.41 ± 0.42

  1. * Shown are the total van der Waals (vdW) interaction energy and the number of contacts between cholesterol and β2AR, when cholesterol is in the IC1, IC2, or EC1 binding site (and similarly for the cholesterol analogues).

  2. Calculations are based on systems having ≥10 mol% cholesterol. Shown here are the average values over different trajectories.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Moutusi Manna
  2. Miia Niemelä
  3. Joona Tynkkynen
  4. Matti Javanainen
  5. Waldemar Kulig
  6. Daniel J Müller
  7. Tomasz Rog
  8. Ilpo Vattulainen
(2016)
Mechanism of allosteric regulation of β2-adrenergic receptor by cholesterol
eLife 5:e18432.
https://doi.org/10.7554/eLife.18432