Homo-oligomerization of the human adenosine A2A receptor is driven by the intrinsically disordered C-terminus

  1. Khanh Dinh Quoc Nguyen
  2. Michael Vigers
  3. Eric Sefah
  4. Susanna Seppälä
  5. Jennifer Paige Hoover
  6. Nicole Star Schonenbach
  7. Blake Mertz
  8. Michelle Ann O'Malley  Is a corresponding author
  9. Songi Han  Is a corresponding author
  1. Department of Chemistry and Biochemistry, University of California, Santa Barbara, United States
  2. Department of Chemical Engineering, University of California, Santa Barbara, United States
  3. C. Eugene Bennett Department of Chemistry, West Virginia University, United States
8 figures, 1 table and 2 additional files

Figures

Figure 1 with 2 supplements
Method for collecting size-exclusion chromatography (SEC) data and assessing A2AR oligomerization.

The SEC data is recorded every second as absorbance at 280 nm. The baseline is corrected to ensure uniform fitting and integration across the peaks. The areas under the curve, resulting from a multiple-Gaussian curve fit, express the population of each oligomeric species. The reported standard errors of integration are within a 95% confidence interval and are calculated from the variance of the fit, not experimental errors. The levels of high-molecular-weight oligomer and dimer are expressed relative to the monomeric population in arbitrary units. A representative calculation defining the oligomer levels is given in the box.

Figure 1—figure supplement 1
The purity and identity of A2AR are confirmed with total protein stain and western blot.

(A) Representative total protein stain (upper panel) and western blot (lower panel) of A2AR-WT during purification. Positive ([+] ctrl) and negative ([–] ctrl) controls consist of 5 OD cell lysate of Saccharomyces cerevisiae BJ5464 cells expressing and not expressing A2AR WT, respectively. ‘IMAC FT’ indicates the flow-through from IMAC step. ‘XAC inactive’ and ‘XAC active’ indicate the fractions that do not and do bind to XAC during the ligand-affinity chromatography step. (B) Representative western blot of A2AR-WT during size-exclusion chromatography (SEC) separation. The fractions are matched to the distinct oligomeric peaks in the SEC chromatogram. Each lane on the blot is from 0.5 mL fractions eluted from a Superdex 200 10/300 GL (GE Healthcare) column. MagicMark protein ladder (LC5602) is used as the molecular weight standard.

Figure 1—figure supplement 1—source data 1

Raw representative total protein stain of A2AR-WT during purification.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data1-v2.tif.zip
Figure 1—figure supplement 1—source data 2

Labeled representative total protein stain of A2AR-WT during purification.

Positive ([+] ctrl) and negative ([–] ctrl) controls consist of 5 OD cell lysate of Saccharomyces cerevisiae BJ5464 cells expressing and not expressing A2AR WT, respectively. ‘IMAC FT’ indicates the flow-through from IMAC step. ‘XAC inactive’ and ‘XAC active’ indicate the fractions that do not and do bind to XAC during the ligand-affinity chromatography step. MagicMark protein ladder (LC5602) is used as the molecular weight standard.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data2-v2.tif.zip
Figure 1—figure supplement 1—source data 3

Raw representative western blot of A2AR-WT during purification.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data3-v2.tif.zip
Figure 1—figure supplement 1—source data 4

Labeled representative western blot of A2AR-WT during purification.

Positive ([+] ctrl) and negative ([–] ctrl) controls consist of 5 OD cell lysate of Saccharomyces cerevisiae BJ5464 cells expressing and not expressing A2AR WT, respectively. ‘IMAC FT’ indicates the flow-through from IMAC step. ‘XAC inactive’ and ‘XAC active’ indicate the fractions that do not and do bind to XAC during the ligand-affinity chromatography step. MagicMark protein ladder (LC5602) is used as the molecular weight standard.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data4-v2.tif.zip
Figure 1—figure supplement 1—source data 5

Raw representative western blot of A2AR-WT during size-exclusion chromatography separation.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data5-v2.tif.zip
Figure 1—figure supplement 1—source data 6

Labeled representative western blot of A2AR-WT during size-exclusion chromatography separation.

Each lane on the blot is from 0.5 mL fractions eluted from a Superdex 200 10/300 GL (GE Healthcare) column. MagicMark protein ladder (LC5602) is used as the molecular weight standard.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp1-data6-v2.tif.zip
Figure 1—figure supplement 2
Size-exclusion chromatographic traces and data distribution of all A2AR variants used in the main text of this study.

(A) Curve fitting using OriginLab of all A2AR variants used in the main text of this study, listed by the order they appear. By default, each oligomeric peak is fitted with one curve using Gaussian distribution and displayed by different color shades, with the high-molecular-weight (HMW) oligomer eluted first (dark orange), followed by the dimer (lighter orange), followed by the monomer (lightest orange). However, the HMW oligomer peak in some cases cannot be fitted with one curve and thus is fitted with two curves instead. This discrepancy can be explained by variation in HMW oligomerization order among the variants. The identity of each peak is confirmed with western blotting. The value and error from the curve fitting of each peak are given in Supplementary file 1. (B) Data distribution of all variants used in this study in comparison to five experimental replicates of A2AR-WT. The C-terminally truncated mutants are represented by different shades of green in increasing darkness corresponding to the increased length of the C-terminus, with the lightest shade representing the mutant with the shortest C-terminus (A316ΔC) and the darkest shade for the mutant with the longest C-terminus (P395ΔC). The levels of dimer and HMW oligomer are expressed relative to the monomeric population in arbitrary unit, with reported errors calculated from the variance of the fit, not experimental variation. There are significant variations in the dimer and HMW oligomer levels among the WT replicates, stemming from experimental errors. These variations are mitigated when the two parameters are added as the data distribution becomes more uniform. Also, the oligomerization levels of the WT replicates are consistently higher than the mutated and truncated variants.

Figure 1—figure supplement 2—source data 1

Raw size-exclusion chromatography data of five experimental replicates of A2AR-WT.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig1-figsupp2-data1-v2.xlsx
Residue C394 helps stabilize A2AR oligomerization via disulfide bonds.

(A) The effect of C394X substitutions on A2AR oligomerization. The levels of dimer (dark colors) and high-molecular-weight oligomer (light colors) are expressed relative to the monomeric population in arbitrary units, with reported errors calculated from the variance of the fit, not experimental variation. (B) Line densitometry of western blot bands on size-exclusion chromatography (SEC)-separated dimeric populations of A2AR-WT and Q372ΔC with and without 5 mM TCEP. The level of dimer is expressed relative to the monomeric population in arbitrary units similarly to the SEC analysis. MagicMark protein ladder (LC5602) is used as the molecular weight standard.

Figure 2—source data 1

Raw western blot of size-exclusion chromatography-separated dimeric populations of A2AR-WT with and without 5 mM TCEP.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig2-data1-v2.tif.zip
Figure 2—source data 2

Raw western blot of size-exclusion chromatography-separated dimeric populations of A2AR-WT with and without 5 mM TCEP.

MagicMark protein ladder (LC5602) is used as the molecular weight standard.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig2-data2-v2.tif.zip
Figure 2—source data 3

Raw western blot of size-exclusion chromatography-separated dimeric populations of A2AR-Q372ΔC with and without 5 mM TCEP.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig2-data3-v2.tif.zip
Figure 2—source data 4

Raw western blot of size-exclusion chromatography-separated dimeric populations of A2AR-Q372ΔC with and without 5 mM TCEP.

MagicMark protein ladder (LC5602) is used as the molecular weight standard.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig2-data4-v2.tif.zip
Figure 2—source data 5

Raw size-exclusion chromatography data of A2AR-WT and C394X variants.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig2-data5-v2.xlsx
Truncating the C-terminus systematically affects A2AR oligomerization.

(A) Depiction of where the truncation points are located on the C-terminus, with region 354–359 highlighted (in black) showing critical residues. (B) The levels of dimer and high-molecular-weight (HMW) oligomer are expressed relative to the monomeric population as an arbitrary unit and plotted against the residue number of the truncation sites, with reported errors calculated from the variance of the fit, not experimental variation. Region 354–359 is emphasized (in black and gray) due to a drastic change in the dimer and HMW oligomer levels. (C) The dependence of A2AR oligomerization on three consecutive charged residues 355ERR357. The substitution of residues 355ERR357 to 355AAA357 is referred to as the ERR:AAA mutations. The levels of dimer and HMW oligomer are expressed relative to the monomeric population as an arbitrary unit, with reported errors calculated from the variance of the fit, not experimental variation.

Figure 3—source data 1

Raw size-exclusion chromatography data of A2AR-WT and C-terminally truncated ΔC variants.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig3-data1-v2.xlsx
Non-bonded interactions of the extended C-terminus of A2AR play a critical role in stabilization of the dimeric interface.

(A) Dimer configurations from cluster analysis in GROMACS of the 394-residue variant identify two major clusters involving either (1) the C-terminus of one protomer and the C-terminus, ICL2, and ICL3 of the second protomer or (2) the C-terminus of one protomer and ICL2, ICL3, and ECL2 of the second protomer. Spheres: residues forming intermolecular electrostatic contacts. (B) Average number of residues that form electrostatic contacts as a function of sequence length of A2AR. (C) Average number of residues that form hydrogen bonds as a function of sequence length of A2AR. The criteria for designating inter-A2AR contacts as electrostatic interactions or hydrogen bonds are described in detail in Materials and methods.

Figure 4—source data 1

Detailed data regarding the multiple interfaces of A2AR and the network of non-bonded interactions that stabilize these interfaces.

(A) Dimer configurations from cluster analysis in GROMACS of the 394-residue variant. (B) Average number of residues that form electrostatic contacts as a function of sequence length of A2AR. (C) Average number of residues that form hydrogen bonds as a function of sequence length of A2AR.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig4-data1-v2.xlsx
Figure 5 with 1 supplement
The effects of ionic strength on the oligomerization of various A2AR variants reveal the involvement of depletion interactions.

The levels of (A) dimer and (B) high-molecular-weight oligomer + dimer are expressed relative to the monomeric population as an arbitrary unit and plotted against ionic strength, with reported errors calculated from the variance of the fit, not experimental variation. NaCl concentration is varied to achieve ionic strengths of 0.15, 0.45, and 0.95 M.

Figure 5—source data 1

Raw size-exclusion chromatography data of various A2AR variants under different ionic strengths of 0.15, 0.45, and 0.95 M.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig5-data1-v2.xlsx
Figure 5—figure supplement 1
The dimer/oligomerization of A2AR is a thermodynamic process where the dimer and high-molecular-weight oligomer once formed are kinetically trapped.

(A) Size-exclusion chromatography (SEC) chromatograms of the consecutive rounds of SEC performed on A2AR-WT and Q372ΔC. The first rounds of SEC are to separate the dimer/oligomer population and the monomer population, while the second rounds of SEC are performed on these SEC-separated populations to assess their stability and reversibility. The total oligomer level is expressed relative to the monomeric population in arbitrary units. (B) Energy diagram depicting A2AR oligomerization progress. The monomer needs to overcome an activation barrier (EA), driven by depletion interactions, to form the dimer/oligomer. Once formed, the dimer/oligomer populations are kinetically trapped by disulfide linkages.

Figure 5—figure supplement 1—source data 1

Raw size-exclusion chromatography (SEC) data of the consecutive rounds of SEC performed on A2AR-WT and Q372ΔC.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig5-figsupp1-data1-v2.xlsx
Figure 6 with 1 supplement
The A2AR C-terminus is prone to aggregation.

(A) Absorbance at 450 nm of the A2AR C-terminus in solution, with NaCl and GdnHCl concentrations varied to achieve ionic strengths 0–4 M. Inset: the solution at ionic strength 4 M achieved with NaCl. The Hofmeister series is provided to show the ability of cations to salt-out (blue) or salt-in (red) proteins. (B) SYPRO orange fluorescence of solutions containing the A2AR C-terminus as the temperature was varied from 20°C to 70°C (gray). The change in fluorescence, measured in relative fluorescence unit (RFU), was calculated by taking the first derivative of the fluorescence curve (black).

Figure 6—source data 1

Detailed data showing the propensity of A2AR C-terminus to aggregate.

(A) Absorbance at 450 nm of the A2AR C-terminus in solution, with NaCl and GdnHCl concentrations varied to achieve ionic strengths 0–4 M. (B) SYPRO orange fluorescence of solutions containing the A2AR C-terminus as the temperature was varied from 20°C to 70°C (gray). The change in fluorescence, measured in relative fluorescence unit (RFU), was calculated by taking the first derivative of the fluorescence values.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig6-data1-v2.xlsx
Figure 6—figure supplement 1
The C-terminus of A2AR can form non-polar contacts.

(A) Hydropathy plot against A2AR residue number showing the hydrophobicity of A2AR C-terminus, scored with ProtScale using method described by Kyte and Doolittle, window size of 3. Positive scores represent hydrophobicity and negative scores hydrophilicity. (B) The non-polar residues in A2AR C-terminus. (C) Average number of residues that form non-polar contacts as a function of sequence length of A2AR. The criteria for designating inter-A2AR contacts as non-polar interactions are described in detail in Materials and methods.

Figure 6—figure supplement 1—source data 1

Detailed data showing the ability of A2AR C-terminus to form non-polar contacts.

(A) Hydropathy plot against A2AR residue number showing the hydrophobicity of A2AR. C-terminus, scored with ProtScale using method described by Kyte and Doolittle, window size of 3. Positive scores represent hydrophobicity and negative scores hydrophilicity. (C) Average number of residues that form non-polar contacts as a function of sequence length of A2AR.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig6-figsupp1-data1-v2.xlsx
Figure 7 with 1 supplement
Visualizing A2AR dimeric interface and observing conformational changes of the TM7 using MD simulations.

(A) Representative snapshot of A2AR-C394ΔC dimers shows salt bridge formation between a sample trajectory. The insets are closeups of the salt bridges, which can be both intra- and intermolecular. The last inset shows a network of salt bridges with the charged cluster 355ERR357 involved. (B) Helical tilt angles for TM7 helix in A2AR as a function of protein length. Systematic truncations of the C-terminus lead to rearrangement of the heptahelical bundle. The participation of the C-terminus in A2AR dimerization increases the tilting of the TM7 domain, which is in closest proximity to the C-terminus.

Figure 7—source data 1

MD simulations data used to visualize A2AR dimeric interface and observe the conformational changes of the TM7.

(A) List of all C-terminal residue pairs of A2AR-C394ΔC dimers engaging in electrostatic interactions. (B) Helical tilt angles for TM7 helix in A2AR as a function of protein length.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig7-data1-v2.xlsx
Figure 7—figure supplement 1
Helical tilt angles for TM1–6 helices in A2AR as a function of protein length.

Systematic truncations of the C-terminus lead to rearrangement of the heptahelical bundle, propagated to the entire receptor and is especially pronounced in helices proximal to the C-terminus, that is, TM1, TM2, TM7. For almost all TM helices, a noticeable shift in tilt angle occurs upon modeling the full-length (394 residues) variant. This behavior is fundamentally different from the conventional model of G protein-coupled receptor (GPCR) activation, in which TM 1, 2, 4, and 7 remain rigid, with TM5 and TM6 undergoing an outward tilt/rotation to enable binding to the cognate G protein. Relaxation of the heptahelical bundle (i.e., an increase in helical tilt) as a function of protein length and dimerization could potentially be critical to our understanding of the activation mechanism of A2AR as past studies have overwhelmingly focused on activation of the monomer.

Figure 7—figure supplement 1—source data 1

Helical tilt angles for TM1–6 helices in A2AR as a function of protein length.

https://cdn.elifesciences.org/articles/66662/elife-66662-fig7-figsupp1-data1-v2.xlsx
Author response image 1
Decorrelation time calculated from representative trajectories of each A2AR truncated system.

Tables

Key resources table
Reagent type
(species)
or resource
DesignationSource or
reference
IdentifiersAdditional
information
Recombinant DNA reagentpITy (plasmid)Parekh et al., 1996
Strain, strain background (Saccharomyces cerevisiae)BJ5464Robinson Lab – Carnegie Mellon University
Strain, strain background (Escherichia coli)BL21 (DE3)Sigma, St. Louis, MO, USA#CMC0014
Chemical compound, drugDDMAnatrace, Maumee, OH, USA#D310
Chemical compound, drugCHAPSAnatrace, Maumee, OH, USA#C216
Chemical compound, drugCHSAnatrace, Maumee, OH, USA#CH210
Chemical compound, drugXanthine amine congenerSigma, St. Louis, MO, USA#X103
Chemical compound, drugTheophyllineSigma, St. Louis, MO, USA#T1633
Commercial assay, kitAffigel 10 resinBioRad, Hercules, CA, USA#1536099
Commercial assay, kitTricorn Superdex 200 10/300 GL columnGE Healthcare, Pittsburgh, PA, USA#17-5175-01
AntibodyAnti-A2AR, clone 7F6-G5-A2 (Mouse monoclonal)Millipore, Burlington, MA, USA#05-717(1:500) dilution
AntibodyAnti-Mouse IgG H&L DyLight 550 (Goat monoclonal)Abcam, Cambridge, MA, USA#ab96880(1:600) dilution
Software, algorithmMODELLER 9.23Eswar et al., 2006
Software, algorithmmartinize.py scriptde Jong et al., 2013
Software, algorithmELNeDyn elastic networkPeriole et al., 2009
Software, algorithmMARTINI coarse-grained force field v2.2Monticelli et al., 2008
Software, algorithmGROMACS 2016Abraham et al., 2015
Software, algorithmbackward.py scriptWassenaar et al., 2014
Software, algorithmLINCSHess et al., 1997
Software, algorithmCHARMM36 and TIP3P force fieldsBest et al., 2012; Jorgensen et al., 1983
Software, algorithmLOOSRomo and Grossfield, 2009
Software, algorithmVMDHumphrey et al., 1996

Additional files

Supplementary file 1

Results from curve fitting using OriginLab and calculations of the high-molecular-weight (HMW) oligomer and dimer levels for all A2AR variants used in the main text of this study.

The variants are grouped by the order they appear and numbered corresponding to Figure 1—figure supplement 2. The levels of dimer and HMW oligomer are expressed relative to the monomeric population in arbitrary units as monomer-equivalent concentration ratios. The errors are calculated from the variance of the fit, not experimental variations, and are within 95% confidence interval. Only the WT replicates are represented with standard deviation as experimental variations (last row; n = 5; mean ± SD).

https://cdn.elifesciences.org/articles/66662/elife-66662-supp1-v2.docx
Transparent reporting form
https://cdn.elifesciences.org/articles/66662/elife-66662-transrepform-v2.docx

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. Khanh Dinh Quoc Nguyen
  2. Michael Vigers
  3. Eric Sefah
  4. Susanna Seppälä
  5. Jennifer Paige Hoover
  6. Nicole Star Schonenbach
  7. Blake Mertz
  8. Michelle Ann O'Malley
  9. Songi Han
(2021)
Homo-oligomerization of the human adenosine A2A receptor is driven by the intrinsically disordered C-terminus
eLife 10:e66662.
https://doi.org/10.7554/eLife.66662