Molecular basis of ligand-dependent Nurr1-RXRα activation

  1. Xiaoyu Yu
  2. Jinsai Shang
  3. Douglas J Kojetin  Is a corresponding author
  1. Skaggs Graduate School of Chemical and Biological Sciences at Scripps Research, United States
  2. Department of Integrative Structural and Computational Biology, Scripps Research and UF Scripps Biomedical Research, United States
  3. Department of Molecular Medicine, Scripps Research and UF Scripps Biomedical Research, United States

Abstract

Small molecule compounds that activate transcription of Nurr1-retinoid X receptor alpha (RXRα) (NR4A2-NR2B1) nuclear receptor heterodimers are implicated in the treatment of neurodegenerative disorders, but function through poorly understood mechanisms. Here, we show that RXRα ligands activate Nurr1-RXRα through a mechanism that involves ligand-binding domain (LBD) heterodimer protein-protein interaction (PPI) inhibition, a paradigm distinct from classical pharmacological mechanisms of ligand-dependent nuclear receptor modulation. NMR spectroscopy, PPI, and cellular transcription assays show that Nurr1-RXRα transcriptional activation by RXRα ligands is not correlated with classical RXRα agonism but instead correlated with weakening Nurr1-RXRα LBD heterodimer affinity and heterodimer dissociation. Our data inform a model by which pharmacologically distinct RXRα ligands (RXRα homodimer agonists and Nurr1-RXRα heterodimer selective agonists that function as RXRα homodimer antagonists) operate as allosteric PPI inhibitors that release a transcriptionally active Nurr1 monomer from a repressive Nurr1-RXRα heterodimeric complex. These findings provide a molecular blueprint for ligand activation of Nurr1 transcription via small molecule targeting of Nurr1-RXRα.

Editor's evaluation

This is a fundamental study of the activation process of Nurr1, an orphan nuclear receptor that may be a significant target for the treatment of neurodegenerative disorders. Nurr1 functions as a monomer, but may also heterodimerize with RXRα which represses Nurr1 transcriptional activation. The authors provide compelling evidence for Nurr1 activation through ligand-induced dissociation of an inactive Nurr1-RXRα heterodimer. These data will be important for biochemists and cell biologists working on regulatory / activation mechanisms of nuclear hormone receptors.

https://doi.org/10.7554/eLife.85039.sa0

Introduction

Nurr1 (nuclear receptor related 1 protein; NR4A2) is a nuclear receptor (NR) transcription factor that is essential for the development, regulation, and maintenance of several important aspects of mammalian brain development and homeostasis. Nurr1 is critical in the development of dopaminergic neurons that are critical for control of movement that degenerates in Parkinson’s disease (PD) (Decressac et al., 2013; Jiang et al., 2005; Zetterström et al., 1997). Recent studies show that Nurr1 is an important factor in the regulation of neuroinflammation and accumulation of amyloid beta that occurs in the pathogenesis of Alzheimer’s disease (AD) (Jeon et al., 2020; Moon et al., 2019). These and other studies implicate small molecule activation of Nurr1 as a potential therapeutic strategy in aging-associated neurodegenerative and dementia disorders characterized by a loss of neuron function (Moutinho et al., 2019).

Although NRs are considered to be ligand-dependent transcription factors, targeting Nurr1 activity with small molecule compounds has remained challenging. Nurr1 contains the conserved NR domain architecture (Weikum et al., 2018) including an N-terminal activation function-1 (AF-1) domain, a central DNA-binding domain (DBD), and a C-terminal ligand-binding domain (LBD). The classical model of ligand-induced NR transcriptional activation occurs via ligand binding to an orthosteric pocket located in the internal core of the LBD, which stabilizes the activation function-2 (AF-2 surface) in an active conformation resulting in high affinity coactivation interaction to the NR and recruitment of coactivator complexes resulting in chromatin remodeling and increased gene expression. However, crystal structures have revealed properties of Nurr1 LBD that suggest it may function in an atypical ligand-independent manner. The structures show no classical orthosteric ligand-binding pocket volume as well as a reversed ‘charge clamp’ AF-2 surface that prevents the LBD from interacting with typical NR transcriptional coregulator proteins (Wang et al., 2003). However, solution-state structural studies indicate the Nurr1 LBD is dynamic and can likely expand to bind endogenous ligands (de Vera et al., 2019). Although endogenous and synthetic ligands have been reported to bind and/or regulate Nurr1 activity (Bruning et al., 2019; de Vera et al., 2016; Kim et al., 2015; Rajan et al., 2020; Willems and Merk, 2022), some reported Nurr1 ligands function independent of binding to its LBD (Munoz-Tello et al., 2020). These and other observations have directed efforts to discover small molecule modulators of Nurr1 activity through other mechanisms.

An alternative way to modulate Nurr1 transcription activity is through ligands that target its NR heterodimer-binding partner, retinoid X receptor alpha (RXRα; NR2B1). Heterodimerizaton with RXRα represses Nurr1 transcription on monomeric DNA response elements called NBREs (Aarnisalo et al., 2002; Forman et al., 1995) that are present in the promoter regions of genes that regulate dopaminergic signaling including tyrosine hydroxylase (Iwawaki et al., 2000; Kim et al., 2003). The observation that classical pharmacological RXRα agonists enhance Nurr1-RXRα transcription on NBREs (Aarnisalo et al., 2002; Forman et al., 1995; Perlmann and Jansson, 1995; Wallen-Mackenzie et al., 2003) inspired the discovery of RXRα-binding ligands that display biased activation of Nurr1-RXRα heterodimers over other NR-RXRα heterodimers or RXRα homodimers (Giner et al., 2015; Morita et al., 2005; Scheepstra et al., 2017; Spathis et al., 2017; Sundén et al., 2016). However, it remains poorly understood how classical RXRα ligands (pharmacological agonists and antagonists) and Nurr1-RXRα selective agonists impact the structure and function of Nurr1-RXRα.

Here, we tested a set of RXRα ligands in structure-function assays to determine their mechanism of action in regulating Nurr1-RXRα transcriptional activity. Unexpectedly, our studies show that ligand activation of Nurr1-RXRα transcription is not associated with classical pharmacological RXRα agonism, which is defined by a ligand-induced increase in coactivator binding to the RXRα LBD resulting in transcriptional activation. Instead, our data support a model whereby Nurr1-RXRα activating ligands weaken Nurr1-RXRα LBD heterodimer affinity via an allosteric protein-protein interaction (PPI) inhibition mechanism and release a transcriptionally active Nurr1 monomer from the repressive Nurr1-RXRα heterodimer.

Results

RXRα LBD is sufficient for repression of Nurr1 transcription

To confirm the published observation that RXRα represses Nurr1 transcription (Aarnisalo et al., 2002; Forman et al., 1995), we performed a transcriptional reporter assay where SK-N-BE(2) neuronal cells were transfected with a full-length Nurr1 expression plasmid, with or without full-length or domain-truncation RXRα expression plasmids, along with a 3xNBRE-luciferase plasmid containing three copies of the monomeric NBRE DNA-binding response element sequence upstream of luciferase gene (Figure 1a). Cotransfected RXRα expression plasmid repressed Nurr1 transcription, and RXRα truncation constructs show that the RXRα LBD is both necessary and sufficient for repression of Nurr1 transcription (Figure 1b). A slight decrease in Nurr1 transcription is observed with the ΔLBD construct though the effect is not statistically significant. These findings implicate an LBD-driven PPI mechanism with a negligible or minor role for the RXRα DBD, which is consistent with published studies that showed full-length RXRα does not bind to monomeric NBRE sequences via the RXRα DBD, but does interact with Nurr1 that is bound to NBRE DNA sequences via a PPI (Sacchetti et al., 2002). Furthermore, RXR ligands can activate transcription of a Gal4 DBD-Nurr1 LBD fusion protein, but not a Gal4 DBD-Nurr1 LBD mutant (Nurr1dim) that cannot heterodimerize with RXRα (Aarnisalo et al., 2002; Wallen-Mackenzie et al., 2003), which further implicates the RXRα LBD in repression of Nurr1 LBD-mediated transcription.

Contribution of retinoid X receptor alpha (RXRα) domains on repressing Nurr1 transcription.

(a) General scheme of the cellular transcriptional reporter assay. (b) 3xNBRE-luciferase assay performed in SK-N-BE(2)-C cells; see Figure 1—source data 1 for data plotted. Data are normalized to empty vector control (n=9 replicates), shown as a box and whiskers plot with boundaries of the box representing the 25th percentile and the 75th percentile, and representative of two or more independent experiments. Statistical testing was performed and p-values were calculated using the Brown-Forsythe and Welch multiple comparisons test of the FL Nurr1 + RXRα constructs conditions relative to FL Nurr1 control condition.

Figure 1—source data 1

Nurr1+ retinoid X receptor alpha (RXRα) truncated construct luciferase reporter data.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig1-data1-v2.zip

RXRα ligands display graded Nurr1-RXRα transcriptional activation

We assembled a set of 14 commercially available RXRα-binding ligands (Figure 2) described in the literature as pharmacological RXRα agonists and antagonists; a mixed activity PPARγ-RXR selective modulator (LG100754) that antagonizes RXR homodimers but agonizes PPAR-RXR and RAR-RXR heterodimers (Lala et al., 1996) and two RXRα-binding compounds (BRF110 and HX600) described as selective agonists of Nurr1-RXRα heterodimers (Morita et al., 2005; Spathis et al., 2017). To determine how the ligands influence Nurr1-RXRα transcription, we performed a transcriptional reporter assay where SK-N-BE(2) neuronal cells were transfected with full-length Nurr1 and RXRα expression plasmids and the 3xNBRE-luc plasmid then treated with compound or DMSO control (Figure 3a). Compounds reported as RXRα agonists increase Nurr1-RXRα transcription (Figure 3b), and among the most efficacious agonist ligands includes endogenous metabolite 9-cis-retinoic acid (9cRA) and bexarotene, the latter of which was reported to display biased activation of Nurr1-RXR heterodimers over RXR homodimers (McFarland et al., 2013). RXRα antagonists showed relatively no change or slightly decreased transcription of Nurr1-RXRα. The two Nurr1-RXRα selective agonists, BRF110 and HX600, displayed the highest activity of all compounds tested. Taken together, the transcriptional reporter data obtained indicates this RXRα ligand set influences Nurr1-RXRα transcription via a graded activation mechanism.

Retinoid X receptor alpha (RXRα) ligands used in this study.

Grouped by pharmacological phenotype, the set includes ligands that classically activate (agonists) or block (antagonists) activation of RXRα homodimers; a mixed activity modulator (LG100754) that antagonizes RXRα homodimers and activates PPARγ-RXRα and RAR-RXRα heterodimers; and two selective activators of Nurr1-RXRα heterodimers (BRF110 and HX600).

Effect of retinoid X receptor alpha (RXRα)-binding ligands on Nurr1-RXRα transcription.

(a) General scheme of the Nurr1-RXRα/3xNBRE-luciferase cellular transcriptional reporter assay. (b) Nurr1-RXRα/3xBNRE-luciferase transcriptional reporter assay performed in SK-N-BE(2)-C cells treated with RXRα ligand (1 µM) or DMSO (dotted line); see Figure 3—source data 1 for data plotted. Data are normalized to DMSO (n=9 replicates), represent the mean ± s.d., and representative of two or more independent experiments. Statistical testing was performed and p-values were calculated using the Brown-Forsythe and Welch multiple comparisons test relative to DMSO control treated condition.

Figure 3—source data 1

Retinoid X receptor alpha (RXRα) ligand treated Nurr1-RXRα/3xNBRE3-luciferase reporter data.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig3-data1-v2.zip

Nurr1-RXRα activation is not correlated with pharmacological RXRα agonism

Agonist binding to the NR LBD stabilizes an active conformation that facilitates coactivator protein interaction at the AF-2 surface via a ‘charge clamp’ mechanism (Savkur and Burris, 2004) that is important for binding LXXLL-containing motifs present within coactivator proteins resulting in an increase in transcription (Kojetin and Burris, 2013). The Nurr1 LBD does not interact with canonical coregulator proteins because it contains a reversed ‘charge clamp’ in its AF-2 surface (Wang et al., 2003), implicating ligand-dependent coactivator interaction with RXRα in the mechanism of Nurr1-RXRα activation.

To determine if there is a correlation between Nurr1-RXRα transcriptional agonism and increased coactivator recruitment to RXRα that results in transcriptional activation within the RXRα ligand set, we performed a time-resolved fluorescence resonance energy transfer (TR-FRET) biochemical assay (Figure 4a) to assess how the compounds affect interaction between the RXRα LBD and a coregulator peptide derived from a cognate coactivator protein PGC-1α (Delerive et al., 2002; Figure 4b). We also performed a transcriptional reporter assay to determine if there is a correlation between RXRα homodimer-mediated transcription (Figure 4c) and Nurr1-RXRα transcription for the RXRα ligand set. We cotransfected HEK293T cells with a full-length RXRα expression plasmid along with a plasmid containing three copies of the dimeric RXR DNA-binding response element sequence upstream of luciferase gene (3xDR1-luc) then treated the cells with RXRα ligand or DMSO control (Figure 4d). In this analysis, we did not include two RXRα ligands, danthron and rhein, which are colored and cause interference in TR-FRET experiments.

Compound profiling for pharmacological retinoid X receptor alpha (RXRα) agonism and correlation to Nurr1-RXRα agonism.

(a) General scheme of the RXRα ligand-binding domain (LBD) time-resolved fluorescence resonance energy transfer (TR-FRET) coactivator peptide interaction assay. (b) TR-FRET ratio measured in the presence of DMSO (dotted line) or compound (2–4 µM); see Figure 4—source data 1 for data plotted. Data are normalized to DMSO control (n=3 biological replicates), represent the mean ± s.d., representative of two or more independent experiments. Statistical testing was performed and p-values were calculated using the Brown-Forsythe and Welch multiple comparisons test relative to the DMSO control treated condition. (c) General scheme of the RXRα/3xDR1-luciferase cellular transcriptional reporter assay. (d) RXRα/3xDR1-luciferase transcriptional reporter assay performed in HEK293T cells treated with compound (1 µM) or DMSO control (dotted line); see Figure 4—source data 2 for data plotted. Data normalized to DMSO (n=6 replicates), represent the mean ± s.d., and representative of two or more independent experiments. Statistical testing was performed and p-values were calculated using the Brown-Forsythe and Welch multiple comparisons test relative to DMSO control treated condition. (e) Correlation plot of RXRα LBD TR-FRET data vs. Nurr1-RXRα cellular transcription data. (f) Correlation plot of RXRα transcriptional reporter data vs. Nurr1-RXRα cellular transcription data. (g) Correlation plot of RXRα transcriptional reporter data vs. RXRα LBD TR-FRET data. Pearson (rp) and Spearman (rs) correlation coefficients and statistical significance testing are reported above the correlation plots.

Figure 4—source data 1

Retinoid X receptor alpha (RXRα) ligand treated RXRα ligand-binding domain (LBD) time-resolved fluorescence resonance energy transfer (TR-FRET) coactivator interaction data.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig4-data1-v2.zip
Figure 4—source data 2

Retinoid X receptor alpha (RXRα) ligand treated RXRα/3xDR1-luciferase reporter data.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig4-data2-v2.zip

We previously showed that ligands displaying graded pharmacological PPARγ agonism show a strong correlation between coactivator peptide recruitment to the LBD in a TR-FRET biochemical assay and cellular transcription of full-length PPARγ (Shang et al., 2019). However, no significant correlation was observed between coactivator peptide recruitment to the RXRα LBD and Nurr1-RXRα transcription for the RXRα ligand set (Figure 4e). Furthermore, no significant correlation is observed between RXRα homodimer transcription and Nurr1-RXRα transcription (Figure 4f). In contrast, a strong correlation is observed between RXRα LBD TR-FRET data and full-length RXRα transcription (Figure 4g) where ligands that increase PGC-1α coactivator peptide interaction to the RXRα LBD activate transcription of full-length RXRα in cells, similar to what we observed for PPARγ agonists with graded agonist transcriptional activity (Shang et al., 2019).

The two selective Nurr1-RXRα activating compounds, BRF110 and HX600, function as a pharmacological antagonist and weak/partial graded agonist in the RXRα homodimer TR-FRET and transcriptional reporter assay, respectively. These findings suggest that the mechanism by which the RXRα ligand set influences Nurr1-RXRα transcription may occur independent of coregulator recruitment to the RXRα LBD and RXRα-mediated transcription—indicating that ligand-dependent RXRα homodimer modulation and Nurr1-RXRα heterodimer modulation may function through distinct mechanisms.

Nurr1-RXRα activation is correlated with weakening LBD heterodimer affinity

Although ligand binding to NRs typically influences coregulator interaction to the AF-2 surface in the LBD, studies have reported that ligand binding can also weaken or strengthen NR LBD homodimerization and heterodimerization and confer selectivity (Kilu et al., 2021; Powell and Xu, 2008; Rehó et al., 2020; Tamrazi et al., 2002). To determine how the RXRα ligand set influences Nurr1-RXRα LBD heterodimerization affinity, we performed isothermal titration calorimetry (ITC) studies where we titrated Nurr1 LBD into apo/ligand-free or ligand-bound RXRα LBD (Figure 5a) and fitted the data to a homodimer competition model that incorporates apo-RXRα LBD homodimerization affinity (see Materials and methods) and dissociation of RXRα LBD homodimers to obtain Nurr1-RXRα LBD heterodimerization affinity (Table 1 and Figure 5b). One limitation of this ITC data analysis is the assumption that the RXRα LBD homodimer affinity is the same for apo RXRα LBD and the various ligand-bound states. Determining ligand-bound RXRα LBD homodimer KD values using ITC via the dimer dissociation dilution (McPhail and Cooper, 1997) would likely be confounded by ligand dissociation. It is also possible that RXRα ligands could change RXRα LBD homodimer affinity, or potentially change the RXRα LBD dimer equilibrium toward a tetrameric form where the Nurr1 LBD ITC titration could include a component that dissociates a ligand-bound RXRα LBD homotetramer (Chen et al., 1998; Gampe et al., 2000; Zhang et al., 2011a). Despite these limitations, a significant correlation between Nurr1-RXRα LBD heterodimerization affinity and Nurr1-RXRα transcription where RXRα ligands that weaken heterodimerization affinity show higher Nurr1-RXRα transcription (Figure 5c).

Compound profiling for effects on Nurr1-retinoid X receptor alpha (RXRα) ligand-binding domain (LBD) heterodimer affinity and correlation to Nurr1-RXRα agonism.

(a) General scheme of the Nurr1-RXRα LBD isothermal titration calorimetry (ITC) experiment. (b) Nurr1-RXRα LBD heterodimer affinities (log M) in the presence of DMSO (dotted line) or compound determined from the fit of the ITC data (n=2 replicates, except n=1 for IRX4204 and PA425) and represent the mean ± s.d.; see Table 1 and Figure 5—source data 1 for data plotted. Statistical testing was performed and p-values were calculated using ordinary one-way ANOVA relative to DMSO control treated condition. (c) Correlation plot of ITC determined Nurr1-RXRα LBD heterodimer KD log M values vs. Nurr1-RXRα cellular transcription data. (d) Correlation plot of ITC determined Nurr1-RXRα LBD heterodimer KD log M values vs. fitted binding enthalpy (ΔH, which is the ΔHAB component of the homodimer competition model; see Materials and methods section for details). (e) Correlation plot of ITC determined Nurr1-RXRα LBD heterodimer KD log M values vs. calculated binding entropy (TΔS). Pearson (rp) and Spearman (rs) correlation coefficients and statistical significance testing are reported above the correlation plot.

Figure 5—source data 1

Raw isothermal titration calorimetry (ITC) thermograms and fitted data.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig5-data1-v2.zip
Table 1
Fitted and calculated isothermal titration calorimetry (ITC) Nurr1-RXRα binding affinity and thermodynamic parameters.
Ligandlog KD (M)ΔH (kcal/mol)TΔS (kcal/mol)
Averages.d.nReplicatesAverages.d.nReplicatesAverages.d.nReplicates
DMSO–5.3350.1182–5.418–5.2522.9000.24023.0702.73010.1790.401210.4629.895
BRF110–4.3760.0202–4.390–4.3625.3650.36125.1105.62011.3350.334211.09911.571
HX600–4.6620.0282–4.681–4.6427.6954.32024.64010.75014.0544.283211.02617.083
9cRA–5.0260.0602–4.984–5.0685.1100.96225.7904.43013.5393.105212.58911.344
Bexarotene–4.6190.0442–4.588–4.6505.5551.33624.6106.50011.8571.396210.87012.844
LG100268–5.0440.1062–5.119–4.9695.1350.14825.2405.03012.0170.294212.22411.809
CD3254–4.8600.0992–4.930–4.7904.8650.38925.1404.59011.4950.524211.86611.125
SR11237–5.0530.0472–5.087–5.0204.9400.96225.6204.26011.8341.026212.56011.109
UVI3003–4.5650.0212–4.580–4.5506.1252.63827.9904.26012.3532.666214.23810.467
LG100754–5.2220.0262–5.204–5.2404.2700.22624.1104.43011.3940.261211.20911.579
IRX4204–4.6730.0001–4.6738.4600.00018.46014.8350.000114.835
Rhein–5.5540.0772–5.499–5.6092.7800.39623.0602.50010.3570.290210.56310.152
HX531–6.1320.1332–6.226–6.0382.2050.02122.1902.22010.5710.160210.68410.457
Danthron–5.5240.0432–5.494–5.5542.9200.17023.0402.80010.4560.111210.53510.378
PA425–5.4170.0001–5.4174.1200.00014.12011.5100.000111.510

To gain additional insight into the thermodynamic mechanism of ligand-induced Nurr1-RXRα LBD heterodimer dissociation, we analyzed the thermodynamic parameters from the fitted ITC data (Table 1). The change in binding enthalpy (ΔH) upon heterodimerization generally shows an endothermic (positive) profile for all conditions, whereas the change in entropy (TΔS) component shows a favorable (positive) profile for all conditions. Correlation analysis reveals a more significant correlation between the Nurr1-RXRα LBD heterodimerization binding affinity and the ΔH binding component (Figure 5d) compared to the ΔS binding component (Figure 5e). Moreover, compounds that decrease Nurr1-RXRα LBD heterodimer binding affinity show an increasing (less favorable) ΔH component and an increasing (more favorable) ΔS component.

Nurr1-RXRα activation is correlated with dissociation of Nurr1 LBD monomer

To obtain structural insight into the consequence of ligand-induced weakening of Nurr1-RXRα LBD heterodimer and the relationship to transcriptional activation of Nurr1-RXRα, we performed protein NMR structural footprinting analysis. We collected 2D [1H,15N]-TROSY-HSQC NMR data of 15N-labeled Nurr1 LBD in the monomer form or heterodimerized with RXRα LBD. Among the well-resolved Nurr1 LBD peaks in the 2D NMR data that show clear and discernible shifts going from monomeric Nurr1 LBD to the Nurr1-RXRα LBD heterodimer (Figure 6—figure supplement 1) include the NMR peak for Thr411 (Figure 6a). Addition of compounds that activate Nurr1-RXRα transcription and weaken Nurr1-RXRα LBD heterodimer affinity result in the appearance of two Thr411 NMR peaks with chemical shift values corresponding to the heterodimer and Nurr1 monomer populations. In contrast, compounds that display little to no effect on Nurr1-RXRα transcription and do not significantly alter Nurr1-RXRα LBD heterodimer affinity show either a single NMR peak corresponding to the heterodimer population or a lower Nurr1 monomer population on average. These data indicate that Nurr1-RXRα activating ligands perturb the Nurr1-RXRα LBD heterodimer conformational ensemble toward a monomeric Nurr1 LBD population in slow exchange on the NMR time scale.

Figure 6 with 3 supplements see all
Compound profiling for effects on Nurr1-retinoid X receptor alpha (RXRα) ligand-binding domain (LBD) conformational properties in solution and correlation to Nurr1-RXRα agonism.

(a) 2D [1H,15N]-TROSY HSQC data of 15N-labeled Nurr1 LBD heterodimerized with unlabeled RXRα LBD in the presence of RXRα ligands focused on the NMR peak of Thr411. The upper inset shows an overlay of 15N-labeled Nurr1 LBD monomer (200 µM) vs. 15N-labeled Nurr1 LBD-unlabeled RXRα LBD heterodimer (1:2 molar ratio) to demonstrate the shift of the Thr411 peak between monomer (m) and heterodimer (hd) forms; see Figure 6—figure supplement 1 for full spectral overlays. (b) Simulated 1H NMR lineshape analysis of Nurr1 LBD residue Thr411 showing the influence of ligand-induced weakening of Nurr1-RXRα LBD heterodimerization affinity; see Figure 6—source code 1 for calculation input files. (c) NMR estimated Nurr1 LBD monomer populations from the 2D NMR data (n=2 replicates, except n=1 for IRX4204 and PA425) and lines above the replicate values represent the mean ± s.d.; see Figure 6—source data 1 for data plotted. Statistical testing was performed and p-values were calculated using ordinary one-way ANOVA relative to apo/DMSO control treated condition. (d) Correlation plot of Nurr1-RXRα cellular transcription data vs. NMR estimated Nurr1 LBD monomer populations. (e) Correlation plot of ITC determined Nurr1-RXRα LBD heterodimer KD log M values vs. NMR estimated Nurr1 LBD monomer populations. Pearson (rp) and Spearman (rs) correlation coefficients and statistical significance testing are reported above the correlation plots. (f) Analytical size exclusion chromatography (SEC) analysis of Nurr1-RXRα in the presence of RXRα ligands (solid colored lines) relative to Nurr1 LBD monomer (dotted black line) and Nurr1-RXRα LBD heterodimer (solid black line); see Figure 6—figure supplement 2 for all SEC data organized by ligand.

Figure 6—source code 1

Input files for NMR LineShapeKin simulated NMR data analysis in MATLAB (two input files and one readme file).

https://cdn.elifesciences.org/articles/85039/elife-85039-fig6-code1-v2.zip
Figure 6—source data 1

Retinoid X receptor alpha (RXRα) ligand treated Nurr1-RXRα ligand-binding domain (LBD) NMR-observed monomer species.

https://cdn.elifesciences.org/articles/85039/elife-85039-fig6-data1-v2.zip

To gain additional insight into how ligand-dependent changes in Nurr1-RXRα LBD heterodimerization affinity influences the NMR data, we performed NMR lineshape simulations. Simulated 1D spectra of the 1HN dimension of the Thr411 NMR peaks were calculated using a binding model that accounts for an RXRα LBD homodimer component that dissociates into monomers that are capable of binding Nurr1 LBD (see Materials and methods for details). The simulations show that weakening Nurr1-RXRα LBD heterodimerization affinity, for example upon binding an RXRα ligand, shifts the equilibrium from the Nurr1-RXRα LBD heterodimeric (hd) population to Nurr1 LBD monomer (m) population (Figure 6b), consistent with our experimentally collected 2D NMR data.

We calculated the relative population of monomeric Nurr1 LBD dissociated from the Nurr1-RXRα LBD heterodimer by measuring the peak intensity of the monomeric (m) and heterodimer (hd) 2D NMR peaks for Thr411 from two replicate measurements using different batches of protein (Figure 6c). One limitation of this NMR analysis is that the monomeric and heterodimeric NMR peak intensities may over- or underestimate the relative population sizes given that NMR peak lineshapes of the individual monomeric and heterodimer states (molecular sizes effects), and the rate of chemical exchange between these states, can also affect peak lineshapes as well as binding kinetics and affinity. Despite this limitation, a significant correlation is observed between NMR-observed Nurr1 LBD monomeric populations and Nurr1-RXRα transcription (Figure 6d) and ITC-determined Nurr1-RXRα LBD heterodimerization affinity (Figure 6e), indicating a role for dissociation of monomeric Nurr1 LBD from Nurr1-RXRα LBD heterodimer in the mechanism of action of the RXRα ligand set. Notably, one of the replicate NMR measurements shows a more significant correlation to Nurr1-RXRα transcription, which may have at least two origins: errors in concentration between two different batches of protein used in the replicate measurements, and two ligands were not included in one of the replicate measurements—excluding these ligands from the analysis of replicate #2 measurements improves the correlation with rp = 0.77 (p=0.0021) and rs = 0.75 (p=0.0042).

To corroborate our ligand-dependent heterodimer dissociation findings, we performed size exclusion chromatography (SEC) experiments to determine how the compounds influence the oligomeric state of the Nurr1-RXRα LBD heterodimer (Figure 6g and Figure 6—figure supplement 2). The Nurr1-RXRα selective agonist, BRF110, shows the largest effect in dissociating a Nurr1 LBD monomer species in the SEC and NMR studies and is overall the most efficacious compound in activating Nurr1-RXRα transcription. The other Nurr1-RXRα selective compound, HX600, and the classical RXRα agonists dissociate a Nurr1 LBD monomer species although to a lower degree compared to BRF110. Most pharmacological RXRα antagonists we studied do not dissociate a Nurr1 LBD monomer species except for UVI3003, which displays a unique SEC profile that includes a monomeric Nurr1 LBD species and at least two species with a larger hydrodynamic radius (Figure 6—figure supplement 2). One higher order species is consistent with the SEC profile of an RXRα LBD homodimer or Nurr1-RXRα LBD heterodimer, and another with a larger hydrodynamic radius that is more compact than RXRα homotetramer, suggesting it may be a UVI3003-bound RXRα LBD dimer population. Notably, our structure-function profiling data indicate UVI3003 functions similar to the Nurr1-RXRα selective ligands (BRF110 and HX600) in that it activates Nurr1-RXRα transcription while antagonizing RXRα homodimers.

The SEC data provide insight into the potential allosteric mechanism(s) by which compounds in the RXRα ligand set influence Nurr1-RXRα heterodimerzation and transcription via effects on RXRα LBD. In the absence of ligand, RXRα exists as a mixture of homodimers and homotetramers (Kersten et al., 1995a), which can be stabilized as dimers upon agonist binding (Kersten et al., 1995b) or tetramers upon antagonist binding (Zhang et al., 2011b). In our SEC data, the most efficacious compounds not only free a Nurr1 LBD monomeric population, but also increase a higher order species with an elution volume consistent with a RXRα LBD homotetramer population. SDS-PAGE analysis confirms this observation and reveals that the higher order species observed when BRF110 is added to Nurr1-RXRα LBD heterodimers contains RXRα LBD only (Figure 6—figure supplement 3). Furthermore, these compounds change the appearance of species in the SEC data with elution properties similar to Nurr1-RXRα heterodimer and RXRα homodimers, suggesting that they may change the homodimer/heterodimer complex conformation and/or species equilibrium. Thus, the allosteric mechanism may involve ligand-induced stabilization of RXRα oligomeric species that inhibit or are less favorable for Nurr1 heterodimerization.

PCA reveals data features associated with Nurr1-RXRα agonism

Finally, we performed principal component analysis (PCA) using all the experimental data collected as input to determine the data features most associated with transcriptional activation of the 3xNBRE-luciferase reporter by Nurr1-RXRα in an unbiased manner (Figure 7a). Two of the ligands, danthron and rhein, were not included in the analysis because we did not have RXRα LBD TR-FRET and RXRα homodimer transcriptional reporter data (vide infra). Approximately 90% of the variance in the data can be explained by the two PCs (Figure 7b). Several features are notable in the plot of the PCA scores and loadings, the latter of which represent correlation between the experimental data and PCs. The ligands cluster in groups consistent with their known pharmacological phenotypes. The loadings show that Nurr1-RXRα transcription, ITC-determined Nurr1-RXRα heterodimer affinity, and NMR-determined Nurr1 monomeric populations are clustered together and clustered along with the Nurr1-RXRα selective ligands, meaning these features of the data are positively correlated with ligand-dependent Nurr1-RXRα heterodimer dissociation. In contrast, RXRα LBD TR-FRET coactivator interaction and RXRα LBD homodimer transcription is positively correlated with classical RXRα agonists, whereas RXRα antagonists are found opposite to this loading trajectory. Furthermore, the loading projections show that the features of classical RXRα agonism and selective Nurr1-RXRα agonism are uncorrelated since the loadings trajectories associated with these features, and their associated ligand clusters, form a perpendicular (~90°) angle.

Principal component analysis (PCA) of experimental data reveals Nurr1-retinoid X receptor alpha (RXRα) agonism and classical RXRα agonism are uncorrelated.

(a) 2D biplot containing the loadings and PC scores of the first two PCs. PC scores of the ligands are colored according to pharmacological activity (see legend) whereas the loadings (experimental variables) are colored light blue. (b) A proportion of variance plot revealing the amount of variance described by each PC individually and cumulatively.

Discussion

Although Nurr1 is a promising drug target for aging-associated neurodegenerative disorders, uncertainty in the druggability of its LBD has motivated research into alternative approaches to influence Nurr1 transcription. Nurr1-RXRα heterodimer activation via RXRα-binding ligands has emerged as a promising modality to develop neuroprotective therapeutic agents or adjuvants to other therapies. Several RXRα ligands that increase Nurr1-RXRα activation have progressed to clinical trials of PD and AD including bexarotene and IRX4204 (McFarland et al., 2013; Wang et al., 2016). Understanding how these and other Nurr1-RXRα activating ligands function on the molecular level may provide a blueprint to develop more efficacious and potent compounds.

Classical pharmacological modulation of NR transcription by agonists, antagonists, or inverse agonists that activate, block, or repress transcription typically function by enhancing coactivator protein interaction, blocking coregulator interaction, or enhancing corepressor interaction, respectively (Gronemeyer et al., 2004). Our studies demonstrate that Nurr1-RXRα activation occurs through an LBD PPI inhibition heterodimer dissociation mechanism that is distinct from the classical pharmacological properties of NRs ligands. This activation mechanism explains why pharmacological RXRα agonists (e.g. bexarotene and IRX4204) and Nurr1-RXRα selective compounds that function as pharmacological RXRα antagonists (e.g. BRF110 and HX600) both activate Nurr1-RXRα transcription. Although the LBD heterodimer dissociation mechanism is atypical, several studies have reported that pharmacological NR ligands can influence dimer formation in addition to their classical pharmacological properties in regulating coregulator protein interaction (Kilu et al., 2021; Powell and Xu, 2008; Rehó et al., 2020; Tamrazi et al., 2002). Our Nurr1-RXRα LBD heterodimer dissociation findings add to the repertoire of functional ligand targeting mechanisms of NRs that include classical pharmacological ligands, ligand-modulated posttranslational modification (Choi et al., 2011), PROTAC degrader molecules (Flanagan and Neklesa, 2019), and ligand-modulated NR phase separation/biomolecular condensates that include covalent targeting of the disordered N-terminal AF-1 domain (Basu et al., 2022; Xie et al., 2022).

Our data in support of a Nurr1-RXRα LBD heterodimer dissociation activation mechanism along with other published studies inform a model for activation of Nurr1-RXRα transcription by RXRα ligands (Figure 8). Nurr1 displays high transcriptional activity on monomeric NBRE sites as a monomer, and RXRα heterodimerization represses Nurr1 activity (Aarnisalo et al., 2002; Forman et al., 1995). RXRα does not bind to NBRE sequences but does interact with Nurr1 bound to NBRE DNA sequences via a protein-protein tethering interaction (Sacchetti et al., 2002) and recruits corepressor proteins that are released upon binding RXRα ligand (Lammi et al., 2008). Our RXRα truncation studies indicate that RXRα represses Nurr1 transcription through an RXRα LBD-dependent mechanism, pinpointing the Nurr1 and RXRα LBDs as the likely mode of tethering. Pharmacological agonists of RXRα, which stabilize an active RXRα LBD conformation (De Lera et al., 2007), as well as Nurr1-RXRα selective compounds (Morita et al., 2005; Spathis et al., 2017) that function as pharmacological antagonists of RXRα, dissociate tethered Nurr1-RXRα heterodimers leaving a monomeric Nurr1 bound to monomeric NBRE sites that has higher transcriptional activity than Nurr1-RXRα.

Data-informed model for activation of Nurr1-retinoid X receptor alpha (RXRα) transcription by RXRα ligands.

In contrast to the ligand-induced heterodimer dissociation activation mechanism, ligands that stabilize/strengthen Nurr1-RXRα heterodimerization or do not significantly affect Nurr1-RXRα dissociation—such as the RXRα antagonists in our study—may show more complex regulatory mechanisms. In addition to influencing Nurr1-RXRα interaction, antagonist ligands bound to RXRα, which is itself tethered to Nurr1 on NBRE genomic elements, could influence transcription of NBRE-driven gene transcription via classical pharmacological modulation mechanisms. RXRα-bound ligands will influence the conformation of and coregulator protein binding affinity to the RXRα AF-2 surface, which will in turn regulate recruitment of coregulator (coactivator or corepressor) proteins via interaction with the RXRα AF-2 surface.

In summary, our study shows that RXRα ligands can function as allosteric PPI inhibitors that bind to the canonical orthosteric ligand-binding pocket within the RXRα LBD and influence and Nurr1-mediated transcription via Nurr1-RXRα heterodimer dissociation. Our findings provide a molecular understanding into the mechanism of action of Nurr1-RXRα activation that may influence the design and development of new compounds with therapeutic efficacy in PD, AD, and other aging-associated neurodegenerative disorders.

Materials and methods

Ligands, plasmids, and other reagents

Request a detailed protocol

All ligands were obtained from commercial vendors including Cayman Chemicals, Sigma, Axon Medchem, MedChemExpress, or Tocris Bioscience; BRF110 (CAS 2095489-35-1), HX600 (CAS 172705-89-4), 9cRA (CAS 5300-03-8), Bexarotene (CAS 153559-49-0), LG100268 (CAS 153559-76-3), CD3254 (CAS 196961-43-0), SR11237 (CAS 146670-40-8), UVI3003 (CAS 847239-17-2), LG100754 (CAS 180713-37-5), IRX4204 (CAS 220619-73-8), Rhein (CAS 478-43-3), HX531 (CAS 188844-34-0), Danthron (CAS 117-10-2), and PA425 (CAS 457657-34-0). FITC-labeled LXXLL-containing peptide derived from human PGC-1α (137-155; EAEEPSLLKKLLLAPANTQ) was synthesized by LifeTein with an N-terminal FITC label and an amidated C-terminus for stability. Bacterial expression plasmids included human Nurr1 (NR4A2) LBD (residues 353–598) (de Vera et al., 2016) and human RXRα (NR2B1) LBD (residues 223–462) (Kojetin et al., 2015) that were inserted into a pET45b(+) plasmid (Novagen) as a TEV-cleavable N-terminal hexahistidine(6xHis)-tag fusion protein. Luciferase reporter plasmids included a 3xNBRE-luciferase plasmid containing three copies of the NGFI-B response element corresponding to the monomeric binding site for Nurr1 de Vera et al., 2016; Wilson et al., 1991; and a 3xDR1-luciferase containing three copies of the optimal direct repeat 1 (DR1) binding site for RXRα (Hughes et al., 2014; Subauste et al., 1994). Mammalian expression plasmids included full-length human Nurr1 (residues 1–598) in pcDNA3.1 plasmid (de Vera et al., 2016) and full-length human RXRα (residues 1–462) in pCMV-Sport6 plasmid (Zhang et al., 2011c). To clone the RXRα ΔLBD construct (residues 1–226), site directed mutagenesis and PCR was used to insert a stop codon before the start of the LBD using the full-length RXRα expression plasmid the following primers: forward primer, CAGCAGCGCCTAAGAGGACATG; reverse primer, CATGTCCTCTTAGGCGCTGCTG. To clone the ΔNTD RXRα construct (residues 127–462), site directed mutagenesis and PCR was used to add a XhoI cut site and an ATG before at the end of the NTD using the following primers: forward primer, CCACCCCTCGAGAAACATGG; reverse primer, CCATGTTTCTCGAGGGGTGG—then restriction enzymes (XhoI and HindIII) were used to cut out the ΔNTD construct (DBD-hinge-LBD) for insertion using T4 ligase into pcDNA3.1 empty vector that had been linearized with XhoI and HindIII. To clone the RXRα LBD construct (200–462), Xhol was used to cut pcDNA3.1 and subsequently Gibson assembly was used to clone the LBD into the linearized vector using the gBlock sequence provided in Supplementary file 1.

Cell lines for mammalian cell culture

Request a detailed protocol

All cells were obtained from and authenticated by ATCC and determined to be mycoplasma free. HEK293T (#CRL-11268) and SK-N-BE(2)-C (#CRL-2268) cells were cultured according to ATCC guidelines at low passage number (less than 10 passages; typically passages 2–4). HEK293T cells were grown at 37°C and 5% CO2 in DEME (Gibco) supplemented with 10% fetal bovine serum (Gibco) and 100 units/mL of penicillin, 100 μg/mL of streptomycin, and 0.292 mg/mL of glutamine (Gibco) until 90–95% confluence in T-75 flasks prior to subculture or use. SK-N-BE(2)-C were grown at 37°C and 5% CO2 in a media containing 1:1 mixture of EMEM (ATCC) and F12 medium (Gibco) supplemented with 10% fetal bovine serum (Gibco) until 90–95% confluence in T-75 flasks prior to subculture or use.

Protein expression and purification

Request a detailed protocol

Proteins were expressed in Escherichia coli BL21(DE3) cells in autoinduction ZY media (without NMR isotopes) or M9 minimal media (using 15NH4Cl for NMR isotopic labeling). For autoinduction expression, cells were grown for 5 hr at 37°C, then 18 hr at 22°C, then centrifuged for harvesting. For M9 expression, cells were grown at 37°C and induced with 1.0 mM isopropyl β-D-thiogalactoside at OD (600 nm) of 0.6, grown for an additional 18 hr at 22°C, and then centrifuged for harvesting. Cell pellets were lysed using sonication and proteins were purified using Ni-NTA affinity chromatography and gel filtration/SEC. TEV protease was used to cleave the 6xHis-tag for all experiments except protein used for TR-FRET. The purified proteins were verified by SDS-PAGE, then stored in a buffer consisting of 20 mM potassium phosphate (pH 7.4), 50 mM potassium chloride, and 0.5 mM EDTA. All studies that used RXRα LBD protein used pooled SEC fractions for the apo-homodimeric form except for analytical SEC of the apo-tetrameric form.

Transcriptional reporter assays

Request a detailed protocol

SK-N-BE(2)-C cells were seeded in 9.5 cm2 cell culture well (Corning) at 0.5 million for transfection using Lipofectamine 2000 (Thermo Fisher Scientific) and Opti-MEM with an empty vector control or full-length Nurr1 expression plasmid (1 μg), with or without full-length or different RXRα truncation constructs (1 μg), and 3xNBRE-luc (1 μg). HEK293T cells were seeded in 9.5 cm2 cell culture well (Corning) at 0.5 million for transfection using Lipofectamine 2000 (Thermo Fisher Scientific) and Opti-MEM with empty vector control or full-length RXRα expression plasmid (1 μg), and 3xDR1-luc (1 μg). After incubation for 16–20 hr, cells were transferred to white 384-well cell culture plates (Thermo Fisher Scientific) at 0.5 million cells/mL in 20 μL total volume/well. After a 4 hr incubation, cells were treated with 20 μL of vehicle control (DMSO) or 1 μM ligand. After a final 16–20 hr incubation, cells were harvested with 20 μL Britelite Plus (PerkinElmer), and luminescence was measured on a BioTek Synergy Neo multimode plate reader. The luminescence readouts were normalized to cells transfected with the empty vector (truncation construct assay) or cells treated with DMSO (assays with or without ligand treatment). Data were plotted using GraphPad Prism. In Figure 1b, statistical testing was performed and p-values were calculated in GraphPad Prism using the Brown-Forsythe and Welch multiple comparisons test (does not assume equal s.d. values among compared conditions) relative to full-length Nurr1 condition. In Figures 3b and 4d, statistical testing was performed and p-values were calculated in GraphPad Prism using the Brown-Forsythe and Welch multiple comparisons test (does not assume equal s.d. values among compared conditions) relative to DMSO control treated condition. Data are representative of two or more independent experiments (n=6 or 9 biological replicates).

TR-FRET coactivator interaction assay

Request a detailed protocol

Assays were performed in 384-well black plates (Greiner) using 22.5 μL final well volume. Each well contained 4 nM 6xHis-RXRα LBD, 1 nM LanthaScreen Elite Tb-anti-His antibody (Thermo Fisher #PV5895), and 400 nM FITC-labeled PGC1α peptide in a buffer containing 20 mM potassium phosphate (pH 7.4), 50 mM KCl, 5 mM TCEP, and 0.005% Tween 20. Ligand stocks were prepared via serial dilution in DMSO and added to wells (10 µM final concentration) in triplicate. The plates were read using BioTek Synergy Neo multimode plate reader after incubation at 4°C for at least 2 hr. The Tb donor was excited at 340 nm, and its emission was measured at 495 nm; the emission of the acceptor FITC was measured at 520 nm. Data were plotted using GraphPad Prism as TR-FRET ratio (measurement at 520 nm/measurement at 495 nm) at a fixed ligand concentration (2–4 µM depending on the starting compound stock concentration). In Figure 4b, statistical testing was performed and p-values were calculated in GraphPad Prism using the Brown-Forsythe and Welch multiple comparisons test (does not assume equal s.d. values among compared conditions) relative to the DMSO control treated condition. Data are representative of two or more independent experiments (n=3 biological replicates) except for one RXRα ligand (IRX4204) where only one experiment (n=3 biological replicates) was performed.

Isothermal titration calorimetry

Request a detailed protocol

Experiments were performed using a MicroCal iTC200. All experiments were solvent-matched and contained 0.25% DMSO (ligand vehicle) final concentration. RXRα LBD homodimerization affinity was measured using a dimer dissociation dilution method (McPhail and Cooper, 1997) by titrating 1 mM RXRα LBD from the syringe into the ITC buffer in the sample cell at 25°C with a 60 s delay between the 2 μL injections using a mixing speed of 1200 rpm for a total of 20 injections. Data were fitted using the dissociation model in MicroCal Origin 6.0 (MicroCal User Manual, section 7.6). RXRα LBD homodimer affinity (KD) and enthalpy were determined to be 16.3 µM and –13.10 kcal/mol, respectively; this information was used in analysis of the interactions between Nurr1 LBD and ligand-bound RXRα LBD to approximate the dissociation step required for Nurr1 LBD to heterodimerize with RXRα LBD. ITC experiments of Nurr1 LBD to RXRα LBD (apo or ligand-bound) was performed by titrating Nurr1 LBD to RXRα LBD in a 10:1 ratio (500 or 1000 μM, and 50 or 100 µM, respectively) incubated with two equivalent of vehicle (DMSO) or ligand, and the run was performed using the same experimental design. Experiments were performed at 25°C in duplicate, except for two RXRα ligands (IRX4204 and PA425) where only one experiment was performed; and one experimental replicate for HX600 was performed at 15°C, which gave a similar KD value to the data collected at 25°C. NITPIC software (Keller et al., 2012) was used to calculate baselines, integrate curves, prepare experimental data for fitting in SEDPHAT (Brautigam et al., 2016), which was used to obtain binding affinities and thermodynamic parameter measurements using a homodimer competition model (A + B + C <-> AB + C <-> AC + B; competing B and C for A, where A and C is RXRα LBD monomer and B is Nurr1 LBD monomer). Final figures were exported using GUSSI (Brautigam, 2015). In Figure 5b, statistical testing was performed and p-values were calculated in GraphPad Prism using ordinary one-way ANOVA relative to DMSO control treated condition.

NMR spectroscopy

Request a detailed protocol

Two-dimensional [1H,15N]-TROSY-HSQC NMR experiments were performed at 298 K on a Bruker 700 MHz NMR instrument equipped with a QCI cryoprobe. Samples were prepared in a buffer containing 20 mM potassium phosphate (pH 7.4), 50 mM potassium chloride, 0.5 mM EDTA, and 10% D2O. Experiments were collected using 200 μM 15N-labeled Nurr1 LBD with or without two molar equivalents of unlabeled RXRα LBD in the absence or presence of RXRα ligands. Data were processed and analyzed using NMRFx (Norris et al., 2016) using published NMR chemical shift assignments for the Nurr1 LBD (Michiels et al., 2010) that we validated and used previously (de Vera et al., 2019; Munoz-Tello et al., 2020). Relative Nurr1 LBD monomer populations were estimated by the relative peak intensities of the monomeric (Imonomer) and heterodimer (Iheterodimer) species using the following equation:

Imonomer_population=ImonomerImonomer+Iheterodimer

NMR peak intensity errors for the monomeric (Emonomer) and heterodimeric (Eheterodimer) states were calculated in NMRFx and propagated in the Imonomer_population calculation using the following equation:

Emonomer_population=(Imonomer(Iheterodimer+Imonomer)2)+(Iheterodimer(Imonomer+Iheterodimer)2)

Two replicates per condition were performed, except for two RXRα ligands (IRX4204 and PA425) where only one experiment was performed; average and s.d. calculations of the replicate Imonomer_population and Emonomer_population were performed in GraphPad Prism. In Figure 6c, statistical testing was performed and p-values were calculated in GraphPad Prism using ordinary one-way ANOVA relative to apo/DMSO control treated condition. We distinguished the two replicate measurements in the analysis, which were performed using different protein and ligand stock solutions, because one batch of samples (replicate 1) showed a more significant correlation to transcription than the other batch of samples.

NMR lineshape simulations

Request a detailed protocol

NMR lineshape analysis was performed using NMR LineShapeKin version 4 (Kovrigin, 2012) and MATLAB R2022a via NMRbox (Maciejewski et al., 2017). 1D NMR lineshapes for Nurr1 LBD residue T411 (1HN dimension), with or without two molar equivalents of RXRα LBD, were simulated using the U_L2 model (ligand binding to a receptor competing with ligand dimerization) where the ligand (RXRα LBD) exists as a homodimer that is dissociated into monomers upon binding to the receptor (15N-labeled Nurr1 LBD). Simulations were performed using several experimentally defined parameters: ITC-determined apo-RXRα LBD homodimer affinity (16 µM); difference in chemical shift in hertz (150 Hz) between the monomer and heterodimer NMR peaks for T411 measured in 2D [1H,15N]-TROSY HSQC NMR data (w0 for R = –150 Hz; w0 for RL = 0 Hz); a receptor concentration of 200 µM (Rtotal = 2e-3); and ligand added at ±2 molar equivalents (LRratio = 1e-3 and 2.0); and relative relaxation rates for the monomer and heterodimer estimated to give peak full width half height (FWHH) line widths in 2D [1H,15N]-TROSY HSQC NMR data (for R, FWHH = 20, R2=10; for RL, FWHH = 40, R2=20 Hz); otherwise default parameters were used.

Analytical SEC

Request a detailed protocol

To prepare Nurr1-RXRα LBD heterodimer for analytical SEC analysis, purified Nurr1 LBD and RXRα LBD (each 700 μM in 2.5 mL) were incubated together at 4°C in their storage buffer containing 20 mM potassium phosphate (pH 7.4), 50 mM potassium chloride, and 0.5 mM EDTA for 16 hr, injected (5 mL total) into HiLoad 16/600 Superdex 75 pg (Cytiva) connected to an AKTA FPLC system (GE Healthcare Life Science). Gel filtration was performed using the same buffer to purify Nurr1-RXRα LBD heterodimer for analysis, collecting 2 mL fraction to isolate heterodimer population into 15 samples of 0.5 mL (at a protein concentration of 300 μM) for analytical gel filtration. Additionally, purified Nurr1 LBD and RXRα LBD (pooled gel filtration fractions consisting of the homodimeric and homotetrameric species) were also used. Samples containing RXRα LBD (homodimer or heterodimer with Nurr1 LBD) were incubated with two molar equivalents of each ligand for 16 hr at 4°C before injection onto Superdex 75 Increase 10/300 GL (Cytiva) connected to the same AKTA FPLC system. The UV chromatograms were exported in CSV format and plotted using an in-house Jupyter notebook Python script using matplotlib and seaborn packages.

Correlation and other statistical analyses

Request a detailed protocol

Correlation plots were performed using GraphPad Prism to calculate Pearson (rp) and Spearman (rs) correlation coefficients and two-tailed p-values between two experimental measurements per plot. Statistical testing of control to variable conditions was performed using one-way ANOVA testing as detailed above. PCA was performed in GraphPad Prism; all experimentally determined data were used as variables (method = standardize, PCs selected based on parallel analysis at the 95% percentile level with 1000 simulations). p-Value statistical significance shorthand conforms to GraphPad Prism standards: n.s., p>0.05; *, p≤0.05; **, p≤0.01; ***, p≤0.001; and ****, p≤0.0001.

Appendix 1

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Escherichia coli)BL21(DE3)Sigma-AldrichCMC0014Electrocompetent cells
Cell line (Homo sapiens)Human embryonic kidney epithelialATCCCRL-11268
Cell line (Homo sapiens)SK-N-BE(2) neuroblastomaATCCCRL-2271
Chemical compound, drugBRF110Sigma-AldrichCAS 2095489-35-1
Chemical compound, drugHX600Axon MedchemCAS 172705- 89-4
Chemical compound, drug9-cis-Retinoic acidCayman ChemicalsCAS 5300-03-8
Chemical compound, drugBexaroteneCayman ChemicalsCAS 153559- 49-0
Chemical compound, drugLG100268Cayman ChemicalsCAS 153559-76-3
Chemical compound, drugCD3254Cayman ChemicalsCAS 196961-43-0
Chemical compound, drugSR11237Tocris BioscienceCAS 146670-40-8
Chemical compound, drugUVI3003Cayman ChemicalsCAS 847239-17-2
Chemical compound, drugLG100754Cayman ChemicalsCAS 180713-37-5
Chemical compound, drugIRX4204MedChemExpressCAS 220619-73-8
Chemical compound, drugRheinSigma-AldrichCAS 478-43- 3
Chemical compound, drugHX531Cayman ChemicalsCAS 188844-34-0
Chemical compound, drugDanthronSigma-AldrichCAS 117-10-2
Chemical compound, drugPA425Tocris BioscienceCAS 457657-34-0
Peptide, recombinant proteinFITC-PGC1αLifeTeinAmino acid sequence: EAEEPSLLKKLLLAPANTQ, with an N-terminal FITC label and an amidated C-terminus.
Recombinant DNA reagentNurr1-ligand binding domain (LBD) (plasmid)de Vera et al., 2016Bacteria expression plasmid
Recombinant DNA reagentRXRα-ligand binding domain (LBD) (plasmid)Kojetin et al., 2015Bacteria expression plasmid
Recombinant DNA reagentpET45b(+)Novagen71327-3
Transfected construct (Photinus pyralis)3xNBRE-luciferase plasmidde Vera et al., 2016Sanger sequenced
Transfected construct (Photinus pyralis)3xDR1-luciferase plasmidHughes et al., 2014Mammalian expression plasmid, Sanger sequencedThe 3xPPRE-luciferase reporter plasmid in the referenced paper was used in our study.
Transfected construct (human)Full-length human Nurr1de Vera et al., 2016Mammalian expression plasmid, Sanger sequenced
Transfected construct (human)Full-length human RXRα in pCMV-Sport6 vectorZhang et al., 2011cMammalian expression plasmid, Sanger sequencedWe obtained this plasmid from Griffin lab at UF Scripps Institute (see referenced paper).
Recombinant DNA reagentpcDNA3.1 empty vectorThermo Fisher ScientificV790-20
Sequence-based reagentRXRα-ΔLBD-FThis paperPCR primer ordered from SigmaCAGCAGCGCCTAAGAGGACATG
Sequence-based reagentRXRα-ΔLBD-RThis paperPCR primer ordered from SigmaCATGTCCTCTTAGGCGCTGCTG
Sequence-based reagentΔNTD-RXRα-FThis paperPCR primer ordered from SigmaCCACCCCTCGAGAAACATGG
Sequence-based reagentΔNTD-RXRα-RThis paperPCR primer ordered from SigmaCCATGTTTCTCGAGGGGTGG
Gene (human)Nurr1 (NR4A2)UniprotFull length: residues 1–598; LBD: residues 353–598
Gene (human)RXRα (NR2B1)UniprotFull length: residues 1–462; LBD: 223–462
Sequence-based reagentRestriction enzymes, ligase for cloningNEBXhoI, HindIII, T4 ligase
Commercial assay or kitGibson assemblyNBEE2611L
Commercial assay or kitBritelite plus Reporter Gene Assay SystemPerkin Elmer6066769
Sequence-based reagentRXRα-LBDThis papergBlock for Gibson assemblyCCAGCACAGTGGCGGCCGCATGAAGCGGGAAGCCGTGCAGGAGGAGCGGCAGCGTGGCAAGGACCGGAACGAGAATGAGGTGGAGTCGACCAGCAGCGCCAACGAGGACATGCCGGTGGAGAGGATCCTGGAGGCTGAGCTGGCCGTGGAGCCCAAGACCGAGACCTACGTGGAGGCAAACATGGGGCTGAACCCCAGCTCGCCGAACGACCCTGTCACCAACATTTGCCAAGCAGCCGACAAACAGCTTTTCACCCTGGTGGAGTGGGCCAAGCGGATCCCACACTTCTCAGAGCTGCCCCTGGACGACCAGGTCATCCTGCTGCGGGCAGGCTGGAATGAGCTGCTCATCGCCTCCTTCTCCCACCGCTCCATCGCCGTGAAGGACGGGATCCTCCTGGCCACCGGGCTGCACGTCCACCGGAACAGCGCCCACAGCGCAGGGGTGGGCGCCATCTTTGACAGGGTGCTGACGGAGCTTGTGTCCAAGATGCGGGACATGCAGATGGACAAGACGGAGCTGGGCTGCCTGCGCGCCATCGTCCTCTTTAACCCTGACTCCAAGGGGCTCTCGAACCCGGCCGAGGTGGAGGCGCTGAGGGAGAAGGTCTATGCGTCCTTGGAGGCCTACTGCAAGCACAAGTACCCAGAGCAGCCGGGAAGGTTCGCTAAGCTCTTGCTCCGCCTGCCGGCTCTGCGCTCCATCGGGCTCAAATGCCTGGAACATCTCTTCTTCTTCAAGCTCATCGGGGACACACCCATTGACACCTTCCTTATGGAGATGCTGGAGGCGCCGCACCAAATGACTTGATCGAGTCTAGAGGGCCCG
Commercial assay or kitLanthaScreen Elite Tb-anti-His antibodyThermo Fisher#PV5895
Software, algorithmNITPIC softwareKeller et al., 2012Baseline calculation, curve integration
Software, algorithmSEDPHATBrautigam et al., 2016Estimation of binding affinity and thermodynamic parameter measurements
Software, algorithmGUSSIBrautigam, 2015Plot ITC figures
OtherNMR chemical shift assignment of Nurr1 LBDMichiels et al., 2010; de Vera et al., 2019; Munoz-Tello et al., 2020BMRB16541Published NMR peak assignment from Biological Magnetic Resonance Data Bank
Software, algorithmNMRFxNorris et al., 2016NMR data process and analysis
Software, algorithmNMR LineShapeKin version 4Kovrigin, 2012NMR lineshape analysis
Software, algorithmMATLAB R2022a via NMRboxMaciejewski et al., 2017NMR lineshape analysis
Software, algorithmPearson and Spearman correlation analysisGraphPad PrismCorrelation analysis
Software, algorithmPrincipal component analysis (PCA);GraphPad PrismCorrelation analysis
Software, algorithmANOVA multiple comparison testGraphPad PrismStatistical testing

Data availability

Raw ITC thermograms and fitted data are provided as Figure 5-source data 1. Input files for NMR LineShapeKin simulated NMR data analysis in MATLAB are provided as Figure 6-source code 1 (zip file including two input files and one readme file). Raw data used for graphical plots are provided as Figure 1-source data 1 (Nurr1 + RXRα truncated construct luciferase reporter data), Figure 3-source data 1 (RXRα ligand treated Nurr1-RXRα/3xNBRE3-luciferase reporter data), Figure 4-source data 1 (RXRα ligand treated RXRα LBD TR-FRET), Figure 4-source data 2 (RXRα ligand treated RXRα/3xDR1-luciferase reporter data), and Figure 6-source data 1 (RXRα ligand treated Nurr1-RXRα LBD NMR-observed monomer species). All other data generated or analyzed during this study are included in the manuscript and supporting files.

References

Decision letter

  1. Volker Dötsch
    Senior and Reviewing Editor; Goethe University, Germany
  2. Sean W Fanning
    Reviewer; Loyola University Chicago, United States

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Molecular basis of ligand-dependent Nurr1-RXRα activation" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Volker Dötsch as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Sean W Fanning (Reviewer #2).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The RXRɑ antagonist HX531 has little to no effect on Nurr1-RXRɑ reporter activity like the other RXRɑ antagonists Rhein and Danthron (Figure 3b) but appears to selectively induce a >5-fold enhancement in the apparent affinity of LBD heterodimerization (Figure 5b) that seems inconsistent with the proposed model. The authors address a potential caveat with their analysis of the ITC data on lines 68-77 of page 4 but don't comment on any specific ligands or data. As a result, the correlation analyses presented in Figure 5c-d, although largely compelling, may be oversimplified. This is also true in Figure 6, where the correlation between reporter activity and Nurr1 monomerization assessed by NMR (Figure 6c) is much weaker than that with the ITC data (Figure 6e). Can the authors comment on the overall rigor and suitability of this simple correlative approach to justify these apparent discrepancies?

2) Related: In the correlation studies the authors include danthron and rhein as compounds in the Nurr1/RXR assay but not in the RXR PPI and RXR CTF assay. It is not clear if they are included in the correlation graphs in Figure 4 either; however, if they are excluded there should be a reason indicated. How do these compounds behave in the other assays?

3) Can the authors clarify how many independent times ITC data were collected for each ligand? The legend for Figure 5 states n = 1, but the methods state that at least two independent experiments were done. If so, the KD values in Figures 5 and 6 should reflect all of the data unless this is justified explicitly.

4) Related: Asterisks are used to indicate significance but the p-values assigned to them are not discussed in the figure legends nor in the methods. Also, in Figure 5b, error bars are not shown for the Kd values, and significance isn't shown between treatments. Adding in this statistical analysis would help judge whether dimerization differences are meaningful between ligands.

5) Since ITC is being used, information on changes to the heats, Gibbs free energy, enthalpy, and entropy should be available. Some discussion on whether changes to the observed binding affinities are ethically or entropically driven would be of interest to the biothermodynamics community.

6) Is there a way to estimate the error of measuring Nurr1 monomer populations by NMR in Figure 6a, b, c, and e? This could also be useful in making correlations in Figure 6.

7) Response at 1 μM was shown but no dose-response studies were performed to show differences in potency. These could be potentially correlated with the later studies to see if there's a relationship between transcriptional activities and dimerization affinity. Are such dose-response curves available?

8) Binding studies were only done for RXRalpha and the coregulator peptide. It would have been interesting to learn if the coregulator peptide binds the RXR-Nurr1 heterodimer. Comments?

9) HX531 appears to significantly enhance heterodimerization, is there any speculation as to why this is?

10) Should state in the cell culture section whether STR profiling was done to confirm cell line identity during the study, similar to mycoplasma testing.

11) In Figure 1, one of the key constructs (the δ LBD) shows high variability in its activity. Although it is not statistically different than the FL-Nurr1 control the variability that ranges from the FL-Nurr1 control to the FL-Nurr1/Fl-RARa group. Clearly, the LBD is playing a role, but my concern is that it might not be quite so clean. Is it possible that there is a contribution from the RXR DBD in this assay?

12) On page 4 line 43 the authors state,"…the two selective Nurr1-RXRa activating compounds BRF110 and HX600 function as pharmacological RXRa antagonists in these assays – a clear example that pharmacological RXRa modulation of Nurr-1/RXRa modulation is distinct". It isn't clear what assays the authors are referring to where these two compounds behaved as antagonists. This paragraph is prior to the presentation of the heterodimer assay where heterodimer formation is inhibited and in the previously presented assays these compounds were either activators (Nurr1/RXR), inactive (RXR coactivator PPI), or inactive/weak activators (RXR). Please clarify.

Reviewer #1 (Recommendations for the authors):

The RXRɑ antagonist HX531 has little to no effect on Nurr1-RXRɑ reporter activity like the other RXRɑ antagonists Rhein and Danthron (Figure 3b) but appears to selectively induce a >5-fold enhancement in the apparent affinity of LBD heterodimerization (Figure 5b) that seems inconsistent with the proposed model. The authors address a potential caveat with their analysis of the ITC data on lines 68-77 of page 4 but don't comment on any specific ligands or data. As a result, the correlation analyses presented in Figure 5c-d, although largely compelling, may be oversimplified. This is also true in Figure 6, where the correlation between reporter activity and Nurr1 monomerization assessed by NMR (Figure 6c) is much weaker than that with the ITC data (Figure 6e). Can the authors comment on the overall rigor and suitability of this simple correlative approach to justify these apparent discrepancies? I don't think more work necessarily needs to be done, but perhaps a more explicit description of caveats and possible implications for their model would be helpful.

Can the authors clarify how many independent times ITC data were collected for each ligand? The legend for Figure 5 states n = 1, but the methods state that at least two independent experiments were done. If so, the KD values in Figures 5 and 6 should reflect all of the data unless this is justified explicitly.

Is there a way to estimate the error of measuring Nurr1 monomer populations by NMR in Figure 6a, b, c, and e? This could also be useful in making correlations in Figure 6.

Reviewer #2 (Recommendations for the authors):

Congratulations on a compelling and high-impact study. I have a few recommendations that I believe will enhance the significance of this study.

Introduction

There's a typo in the first sentence of paragraph 2.

* are used for significance but the p-values assigned to them are not discussed in the figure legends nor in the methods that I could find.

Results

The second section discussing RXR ligands and Nurr1 transcriptional activity.

– Response at 1 μM was shown but no dose-response studies were performed to show differences in potency. These could be potentially correlated with the later studies to see if there's a relationship between transcriptional activities and dimerization affinity.

Third section looking at RXR coregulator binding.

– Binding studies were only done for RXRalpha and the coregulator peptide. It would have been interesting to learn if the coregulator peptide binds the RXR-Nurr1 heterodimer.

Fourth section using ITC to study heterodimerization.

In figure 5b, error bars are not shown for the Kd values, and significance isn't shown between treatments. Adding in this statistical analysis would help judge whether dimerization differences are meaningful between ligands.

– HX531 appears to significantly enhance heterodimerization, is there any speculation as to why this is?

– Since ITC is being used, you should have information on changes to the heats, Gibbs free energy, enthalpy, and entropy. Some discussion on whether changes to the observed binding affinities are enthalpically or entropically driven would be of interest to the biothermodynamics community.

Materials and methods

Should state in the cell culture section whether STR profiling was done to confirm cell line identity during your study, similar with mycoplasma testing.

Overall: Are there any inferences for this allostery that could be made based on existing RXRalpha co-crystal structures? An analysis of RXR allostery at the end could point to the mechanism of Nurr1 specificity. This would be a very interesting future study.

Reviewer #3 (Recommendations for the authors):

In Figure 1, one of the key constructs (the δ LBD) shows high variability in its activity. Although it is not statistically different than the FL-Nurr1 control the variability that ranges from the FL-Nurr1 control to the FL-Nurr1/Fl-RARa group. Clearly, the LBD is playing a role, but my concern is that it might not be quite so clean. Is it possible that there is a contribution from the RXR DBD in this assay? Otherwise, what do the authors believe causes this variability in this particular group?

In the correlation studies the authors include danthron and rhein as compounds in the Nurr1/RXR assay but not in the RXR PPI and RXR CTF assay. It is not clear if they are included in the correlation graphs in Figure 4 either; however, if they are excluded there should be a reason indicated. How do these compounds behave in the other assays?

On page 4 line 43 the authors state,"…the two selective Nurr1-RXRa activating compounds BRF110 and HX600 function as pharmacological RXRa antagonists in these assays – a clear example that pharmacological RXRa modulation of Nurr-1/RXRa modulation is distinct". It isn't clear to me what assays the authors are referring to where these two compounds behaved as antagonists. This paragraph is prior to the presentation of the heterodimer assay where heterodimer formation is inhibited and in the previously presented assays these compounds were either activators (Nurr1/RXR), inactive (RXR coactivator PPI), or inactive/weak activators (RXR).

https://doi.org/10.7554/eLife.85039.sa1

Author response

Essential revisions:

1) The RXRɑ antagonist HX531 has little to no effect on Nurr1-RXRɑ reporter activity like the other RXRɑ antagonists Rhein and Danthron (Figure 3b) but appears to selectively induce a >5-fold enhancement in the apparent affinity of LBD heterodimerization (Figure 5b) that seems inconsistent with the proposed model. The authors address a potential caveat with their analysis of the ITC data on lines 68-77 of page 4 but don't comment on any specific ligands or data. As a result, the correlation analyses presented in Figure 5c-d, although largely compelling, may be oversimplified. This is also true in Figure 6, where the correlation between reporter activity and Nurr1 monomerization assessed by NMR (Figure 6c) is much weaker than that with the ITC data (Figure 6e). Can the authors comment on the overall rigor and suitability of this simple correlative approach to justify these apparent discrepancies?

We do not believe that the RXRα antagonist data is inconsistent with the ligand-induced activation model and have additional text to the manuscript discussion and a new Results section titled “PCA reveals data features associated with Nurr1-RXRα agonism” where we used principal component analysis (PCA) of all of the experimental data and a new discussion paragraph (second to last) describing how Nurr1-RXRα agonism/heterodimer dissociation and antagonism/heterodimer stabilization likely function via different regulatory mechanisms.

The reviewer also questions the significance of the correlation analysis performed in our study. In our original submission, we used a linear correlation analysis by plotting two sets of experimentally determined data, fitting the data via linear regression, and reporting R2 values as a description of correlation. However, since linear fits are not appropriate in cases where both variables are experimentally measured, in this revised manuscript we instead performed Pearson (rp) and Spearman (rs) correlation analyses and report the associated p-values that report on whether the correlation coefficients (rp or rs) are statistically significant (see grey boxes above plots). Furthermore, the new PCA analysis we added also enhances the rigor of this analysis by analyzing all of the experimental data together in an unbiased manner.

2) Related: In the correlation studies the authors include danthron and rhein as compounds in the Nurr1/RXR assay but not in the RXR PPI and RXR CTF assay. It is not clear if they are included in the correlation graphs in Figure 4 either; however, if they are excluded there should be a reason indicated. How do these compounds behave in the other assays?

We did not include danthron and rhein in the RXRα LBD coactivator peptide interaction TR-FRET assay (Figure 4b) because the compounds are colored and would interfere with the assay. We ran the RXRα/3xDR1-luciferase cellular transcriptional reporter assay (Figure 4d) at the same time as the TR-FRET assay, and therefore we did not include danthron and rhein since no comparison to the TR-FRET data could be made. Danthron and rhein data from the Nurr1-RXRα/3xNBRE-luciferase cellular transcriptional reporter assay are not included in the correlation plots in Figure 4e-g.

Danthron and rhein do not appear to cause interference in the 3xNBRE-luciferase cellular reporter assay we used to measure Nurr1-RXRα transcription, nor in the ITC, NMR, or SEC experiments that measure Nurr1-RXRα LBD interaction; and were therefore included in correlation plots among these experiments.

3) Can the authors clarify how many independent times ITC data were collected for each ligand? The legend for Figure 5 states n = 1, but the methods state that at least two independent experiments were done. If so, the KD values in Figures 5 and 6 should reflect all of the data unless this is justified explicitly.

In this revision the Figures now show data points for duplicate ITC measurements for all but two RXRα ligands (PA425 and IRX4204).

4) Related: Asterisks are used to indicate significance but the p-values assigned to them are not discussed in the figure legends nor in the methods. Also, in Figure 5b, error bars are not shown for the Kd values, and significance isn't shown between treatments. Adding in this statistical analysis would help judge whether dimerization differences are meaningful between ligands.

As requested by the reviewer, we performed statistical testing and p-value analysis for ligand-treated conditions to control conditions (apo/DMSO) for all experimentally measured values in all relevant figures and ensured that all plots include error bars.

5) Since ITC is being used, information on changes to the heats, Gibbs free energy, enthalpy, and entropy should be available. Some discussion on whether changes to the observed binding affinities are ethically or entropically driven would be of interest to the biothermodynamics community.

As requested by the reviewer, we included the measured and fitted ITC thermodynamic parameters to Table 1 (new in the revised manuscript), perform correlation analyses between binding affinity and enthalpy/entropy, and include a discussion of the data and analyses.

6) Is there a way to estimate the error of measuring Nurr1 monomer populations by NMR in Figure 6a, b, c, and e? This could also be useful in making correlations in Figure 6.

We updated Figure 6 to include peak intensity errors calculated by the software we used for NMR analysis; additional details can be found in the methods section.

7) Response at 1 μM was shown but no dose-response studies were performed to show differences in potency. These could be potentially correlated with the later studies to see if there's a relationship between transcriptional activities and dimerization affinity. Are such dose-response curves available?

In this study we are interested in the functional efficacy of the ligand when completely bound to RXRα and therefore performed studies using a single concentration of ligand that is near/above their reported RXRα LBD binding affinity. It is possible there could be a correlation between ligand potency and function, or heterodimerization with Nurr1 could influence ligand affinity/potency vs. RXRα homodimers, which could be apparent in ligand dose-response experiments; we plan to investigate this in future studies.

8) Binding studies were only done for RXRalpha and the coregulator peptide. It would have been interesting to learn if the coregulator peptide binds the RXR-Nurr1 heterodimer. Comments?

The Nurr1 LBD contains a reversed “charge clamp” in its AF-2 surface—a proper charge clamp is important for binding LXXLL-containing motifs present within transcriptional coactivator proteins—and as such Nurr1 does not interact with canonical nuclear receptor coregulator proteins. We updated the first paragraph of the Results section titled “Nurr1-RXRα activation is not correlated with pharmacological RXRα agonism” and included three references to introduce this concept.

If were to perform the TR-FRET coactivator peptide interaction assay to determine how RXRα ligands influence coactivator peptide interaction with Nurr1-RXRα, we anticipate the results would be confounded by several different mechanisms including ligand-dependent Nurr1-RXRα dissociation; ligand-dependent changes in coactivator peptide interaction with RXRα; and ligand-dependent changes in RXRα oligomerization as our size exclusion chromatography (SEC) data indicate ligands differentially affect the hydrodynamic solution properties of RXRα LBD.

9) HX531 appears to significantly enhance heterodimerization, is there any speculation as to why this is?

Our size exclusion chromatography (SEC) data (Figure 6—figure supplement 2) shows that HX531 decreases the hydrodynamic radius of the RXRα LBD (shifts the SEC peak to the right compared to apo-RXRα LBD), which could indicate HX531 changes the RXRα LBD monomer-dimer equilibrium towards a monomeric state (or more generally towards a homodimer conformation that is more compact). Because Nurr1 heterodimerization with RXRα competes for RXRα LBD homodimerization, ligands that favor a monomeric RXRα LBD state would result in higher affinity heterodimerization with Nurr1 LBD. This is an interesting point that we plan to investigate in future studies.

10) Should state in the cell culture section whether STR profiling was done to confirm cell line identity during the study, similar to mycoplasma testing.

We updated the methods cell culture section to state that cells were authenticated by ATCC and cultured at low passage number (less than 10 passages).

11) In Figure 1, one of the key constructs (the δ LBD) shows high variability in its activity. Although it is not statistically different than the FL-Nurr1 control the variability that ranges from the FL-Nurr1 control to the FL-Nurr1/Fl-RARa group. Clearly, the LBD is playing a role, but my concern is that it might not be quite so clean. Is it possible that there is a contribution from the RXR DBD in this assay?

In the first paragraph of the Results section, we updated the description of these data in the Results section and added p-values for the statistical comparison of the RXRα constructs (full-length and truncation) vs. full-length Nurr1. We also added additional information and two more references that implicate a Nurr1 LBD-RXRα LBD interaction in the mechanism by which RXRα represses Nurr1 transcription.

12) On page 4 line 43 the authors state,"…the two selective Nurr1-RXRa activating compounds BRF110 and HX600 function as pharmacological RXRa antagonists in these assays – a clear example that pharmacological RXRa modulation of Nurr-1/RXRa modulation is distinct". It isn't clear what assays the authors are referring to where these two compounds behaved as antagonists. This paragraph is prior to the presentation of the heterodimer assay where heterodimer formation is inhibited and in the previously presented assays these compounds were either activators (Nurr1/RXR), inactive (RXR coactivator PPI), or inactive/weak activators (RXR). Please clarify.

We revised this section of the manuscript removing this confusing statement and added a new section at the end of the Results section titled “Nurr1-RXRα activation is not correlated with pharmacological RXRα agonism” to more concisely make the point that ligand-dependent RXRα homodimer modulation and Nurr1-RXRα heterodimer modulation may function through distinct mechanisms.

Reviewer #1 (Recommendations for the authors):

The RXRɑ antagonist HX531 has little to no effect on Nurr1-RXRɑ reporter activity like the other RXRɑ antagonists Rhein and Danthron (Figure 3b) but appears to selectively induce a >5-fold enhancement in the apparent affinity of LBD heterodimerization (Figure 5b) that seems inconsistent with the proposed model. The authors address a potential caveat with their analysis of the ITC data on lines 68-77 of page 4 but don't comment on any specific ligands or data. As a result, the correlation analyses presented in Figure 5c-d, although largely compelling, may be oversimplified. This is also true in Figure 6, where the correlation between reporter activity and Nurr1 monomerization assessed by NMR (Figure 6c) is much weaker than that with the ITC data (Figure 6e). Can the authors comment on the overall rigor and suitability of this simple correlative approach to justify these apparent discrepancies? I don't think more work necessarily needs to be done, but perhaps a more explicit description of caveats and possible implications for their model would be helpful.

Addressed in our response to Essential revisions #1.

Can the authors clarify how many independent times ITC data were collected for each ligand? The legend for Figure 5 states n = 1, but the methods state that at least two independent experiments were done. If so, the KD values in Figures 5 and 6 should reflect all of the data unless this is justified explicitly.

Addressed in our response to Essential revisions #3.

Is there a way to estimate the error of measuring Nurr1 monomer populations by NMR in Figure 6a, b, c, and e? This could also be useful in making correlations in Figure 6.

Addressed in our response to Essential revisions #6.

Reviewer #2 (Recommendations for the authors):

Congratulations on a compelling and high-impact study. I have a few recommendations that I believe will enhance the significance of this study.

Introduction

There's a typo in the first sentence of paragraph 2.

We fixed it, thanks!

* are used for significance but the p-values assigned to them are not discussed in the figure legends nor in the methods that I could find.

We added the p-values to the figures and described the statistical testing in the figure legends and Materials and methods section.

Results

The second section discussing RXR ligands and Nurr1 transcriptional activity.

– Response at 1 μM was shown but no dose-response studies were performed to show differences in potency. These could be potentially correlated with the later studies to see if there's a relationship between transcriptional activities and dimerization affinity.

Addressed in our response to Essential revisions #7.

Third section looking at RXR coregulator binding.

– Binding studies were only done for RXRalpha and the coregulator peptide. It would have been interesting to learn if the coregulator peptide binds the RXR-Nurr1 heterodimer.

Addressed in our response to Essential revisions #8.

Fourth section using ITC to study heterodimerization.

In figure 5b, error bars are not shown for the Kd values, and significance isn't shown between treatments. Adding in this statistical analysis would help judge whether dimerization differences are meaningful between ligands.

We used GraphPad Prism to calculate Pearson (rp) and Spearman (rs) correlation coefficients and two-tailed p-values to all appropriate figure panels.

– HX531 appears to significantly enhance heterodimerization, is there any speculation as to why this is?

Addressed in our response to Essential revisions #1.

– Since ITC is being used, you should have information on changes to the heats, Gibbs free energy, enthalpy, and entropy. Some discussion on whether changes to the observed binding affinities are enthalpically or entropically driven would be of interest to the biothermodynamics community.

Addressed in our response to Essential revisions #5.

Materials and methods

Should state in the cell culture section whether STR profiling was done to confirm cell line identity during your study, similar with mycoplasma testing.

Addressed in our response to Essential revisions #10.

Overall: Are there any inferences for this allostery that could be made based on existing RXRalpha co-crystal structures? An analysis of RXR allostery at the end could point to the mechanism of Nurr1 specificity. This would be a very interesting future study.

We added a new paragraph to the end of the Results section “Nurr1-RXRα activation is correlated with dissociation of Nurr1 LBD monomer” that describes how our data along with previously reported studies inform on the potential allosteric mechanism by which RXRα ligands influence Nurr1-RXRα heterodimer dissociation. We agree that further studies are warranted to explore this further.

Reviewer #3 (Recommendations for the authors):

In Figure 1, one of the key constructs (the δ LBD) shows high variability in its activity. Although it is not statistically different than the FL-Nurr1 control the variability that ranges from the FL-Nurr1 control to the FL-Nurr1/Fl-RARa group. Clearly, the LBD is playing a role, but my concern is that it might not be quite so clean. Is it possible that there is a contribution from the RXR DBD in this assay? Otherwise, what do the authors believe causes this variability in this particular group?

Addressed in our response to Essential revisions #11.

In the correlation studies the authors include danthron and rhein as compounds in the Nurr1/RXR assay but not in the RXR PPI and RXR CTF assay. It is not clear if they are included in the correlation graphs in Figure 4 either; however, if they are excluded there should be a reason indicated. How do these compounds behave in the other assays?

Addressed in our response to Essential revisions #2.

On page 4 line 43 the authors state,"…the two selective Nurr1-RXRa activating compounds BRF110 and HX600 function as pharmacological RXRa antagonists in these assays – a clear example that pharmacological RXRa modulation of Nurr-1/RXRa modulation is distinct". It isn't clear to me what assays the authors are referring to where these two compounds behaved as antagonists. This paragraph is prior to the presentation of the heterodimer assay where heterodimer formation is inhibited and in the previously presented assays these compounds were either activators (Nurr1/RXR), inactive (RXR coactivator PPI), or inactive/weak activators (RXR).

Addressed in our response to Essential revisions #12.

https://doi.org/10.7554/eLife.85039.sa2

Article and author information

Author details

  1. Xiaoyu Yu

    1. Skaggs Graduate School of Chemical and Biological Sciences at Scripps Research, Jupiter, United States
    2. Department of Integrative Structural and Computational Biology, Scripps Research and UF Scripps Biomedical Research, Jupiter, United States
    Present address
    Department of Biochemistry, Vanderbilt University School of Medicine, Nashville, United States
    Contribution
    Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0549-9560
  2. Jinsai Shang

    Department of Integrative Structural and Computational Biology, Scripps Research and UF Scripps Biomedical Research, Jupiter, United States
    Present address
    School of Basic Medical Sciences, Guangzhou Laboratory, Guangzhou Medical University, Guangzhou, China
    Contribution
    Data curation, Formal analysis, Validation, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  3. Douglas J Kojetin

    1. Department of Integrative Structural and Computational Biology, Scripps Research and UF Scripps Biomedical Research, Jupiter, United States
    2. Department of Molecular Medicine, Scripps Research and UF Scripps Biomedical Research, Jupiter, United States
    Present address
    1. Center for Structural Biology, Vanderbilt University School of Medicine, Nashville, United States
    2. Vanderbilt Institute of Chemical Biology, Vanderbilt University School of Medicine, Nashville, United States
    3. Department of Biochemistry, Vanderbilt University School of Medicine, Nashville, United States
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    douglas.kojetin@vanderbilt.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8058-6168

Funding

National Institute on Aging (R01AG070719)

  • Douglas J Kojetin

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by National Institutes of Health (NIH) grant R01AG070719 from the National Institute on Aging (NIA).

Senior and Reviewing Editor

  1. Volker Dötsch, Goethe University, Germany

Reviewer

  1. Sean W Fanning, Loyola University Chicago, United States

Version history

  1. Preprint posted: November 8, 2022 (view preprint)
  2. Received: November 18, 2022
  3. Accepted: April 26, 2023
  4. Accepted Manuscript published: April 27, 2023 (version 1)
  5. Version of Record published: June 12, 2023 (version 2)

Copyright

© 2023, Yu et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 853
    Page views
  • 167
    Downloads
  • 2
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

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. Xiaoyu Yu
  2. Jinsai Shang
  3. Douglas J Kojetin
(2023)
Molecular basis of ligand-dependent Nurr1-RXRα activation
eLife 12:e85039.
https://doi.org/10.7554/eLife.85039

Further reading

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Daniel Muñoz-Reyes, Levi J McClelland ... Maria Jose Sanchez-Barrena
    Research Article

    The Neuronal Calcium Sensor 1, an EF-hand Ca2+ binding protein, and Ric-8A coregulate synapse number and probability of neurotransmitter release. Recently, the structures of Ric-8A bound to Ga have revealed how Ric-8A phosphorylation promotes Ga recognition and activity as a chaperone and guanine nucleotide exchange factor. However, the molecular mechanism by which NCS-1 regulates Ric-8A activity and its interaction with Ga subunits is not well understood. Given the interest in the NCS-1/Ric-8A complex as a therapeutic target in nervous system disorders, it is necessary to shed light on this molecular mechanism of action at atomic level. We have reconstituted NCS-1/Ric-8A complexes to conduct a multimodal approach and determine the sequence of Ca2+ signals and phosphorylation events that promote the interaction of Ric-8A with Ga. Our data show that the binding of NCS-1 and Ga to Ric-8A are mutually exclusive. Importantly, NCS-1 induces a structural rearrangement in Ric-8A that traps the protein in a conformational state that is inaccessible to Casein Kinase II-mediated phosphorylation, demonstrating one aspect of its negative regulation of Ric-8A-mediated G-protein signaling. Functional experiments indicate a loss of Ric-8A GEF activity towards Ga when complexed with NCS-1, and restoration of nucleotide exchange activity upon increasing Ca2+ concentration. Finally, the high-resolution crystallographic data reported here define the NCS-1/Ric-8A interface and will allow the development of therapeutic synapse function regulators with improved activity and selectivity.

    1. Biochemistry and Chemical Biology
    2. Cell Biology
    Riham Ayoubi, Joel Ryan ... Carl Laflamme
    Research Advance

    Antibodies are critical reagents to detect and characterize proteins. It is commonly understood that many commercial antibodies do not recognize their intended targets, but information on the scope of the problem remains largely anecdotal, and as such, feasibility of the goal of at least one potent and specific antibody targeting each protein in a proteome cannot be assessed. Focusing on antibodies for human proteins, we have scaled a standardized characterization approach using parental and knockout cell lines (Laflamme et al., 2019) to assess the performance of 614 commercial antibodies for 65 neuroscience-related proteins. Side-by-side comparisons of all antibodies against each target, obtained from multiple commercial partners, have demonstrated that: (i) more than 50% of all antibodies failed in one or more applications, (ii) yet, ~50–75% of the protein set was covered by at least one high-performing antibody, depending on application, suggesting that coverage of human proteins by commercial antibodies is significant; and (iii) recombinant antibodies performed better than monoclonal or polyclonal antibodies. The hundreds of underperforming antibodies identified in this study were found to have been used in a large number of published articles, which should raise alarm. Encouragingly, more than half of the underperforming commercial antibodies were reassessed by the manufacturers, and many had alterations to their recommended usage or were removed from the market. This first study helps demonstrate the scale of the antibody specificity problem but also suggests an efficient strategy toward achieving coverage of the human proteome; mine the existing commercial antibody repertoire, and use the data to focus new renewable antibody generation efforts.