Evolutionary selection across the nuclear hormone receptor superfamily with a focus on the NR1I subfamily (vitamin D,pregnane X,and constitutive androstane receptors)

Background

The nuclear hormone receptor (NR) superfamily complement in humans is composed of 48 genes with diverse roles in metabolic homeostasis, development, and detoxification. In general, NRs are strongly conserved between vertebrate species, and few examples of molecular adaptation (positive selection) within this superfamily have been demonstrated. Previous studies utilizing two-species comparisons reveal strong purifying (negative) selection of most NR genes, with two possible exceptions being the ligand-binding domains (LBDs) of the pregnane X receptor (PXR, NR1I2) and the constitutive androstane receptor (CAR, NR1I3), two proteins involved in the regulation of toxic compound metabolism and elimination. The aim of this study was to apply detailed phylogenetic analysis using maximum likelihood methods to the entire complement of genes in the vertebrate NR superfamily. Analyses were carried out both across all vertebrates and limited to mammals and also separately for the two major domains of NRs, the DNA-binding domain (DBD) and LBD, in addition to the full-length sequences. Additional functional data is also reported for activation of PXR and the vitamin D receptor (VDR; NR1I1) to gain further insight into the evolution of the NR1I subfamily.

Results

The NR genes appear to be subject to strong purifying selection, particularly in the DBDs. Estimates of the ratio of the non-synonymous to synonymous nucleotide substitution rates (the ω ratio) revealed that only the PXR LBD had a sub-population of codons with an estimated ω ratio greater than 1. CAR was also unusual in showing high relative ω ratios in both the DBD and LBD, a finding that may relate to the recent appearance of the CAR gene (presumably by duplication of a pre-mammalian PXR gene) just prior to the evolution of mammals. Functional analyses of the NR1I subfamily show that human and zebrafish PXRs show similar activation by steroid hormones and early bile salts, properties not shared by sea lamprey, mouse, or human VDRs, or by Xenopus laevis PXRs.

Conclusion

NR genes generally show strong sequence conservation and little evidence for positive selection. The main exceptions are PXR and CAR, genes that may have adapted to cross-species differences in toxic compound exposure.

     

Functional evolution of mammalian odorant receptors

The mammalian odorant receptor (OR) repertoire is an attractive model to study evolution, because ORs have been subjected to rapid evolution between species, presumably caused by changes of the olfactory system to adapt to the environment. However, functional assessment of ORs in related species remains largely untested. Here we investigated the functional properties of primate and rodent ORs to determine how well evolutionary distance predicts functional characteristics. Using human and mouse ORs with previously identified ligands, we cloned 18 OR orthologs from chimpanzee and rhesus macaque and 17 mouse-rat orthologous pairs that are broadly representative of the OR repertoire. We functionally characterized the in vitro responses of ORs to a wide panel of odors and found similar ligand selectivity but dramatic differences in response magnitude. 87% of human-primate orthologs and 94% of mouse-rat orthologs showed differences in receptor potency (EC50) and/or efficacy (dynamic range) to an individual ligand. Notably dN/dS ratio, an indication of selective pressure during evolution, does not predict functional similarities between orthologs. Additionally, we found that orthologs responded to a common ligand 82% of the time, while human OR paralogs of the same subfamily responded to the common ligand only 33% of the time. Our results suggest that, while OR orthologs tend to show conserved ligand selectivity, their potency and/or efficacy dynamically change during evolution, even in closely related species. These functional changes in orthologs provide a platform for examining how the evolution of ORs can meet species-specific demands.     

Ligand specificity and evolution of liver X receptors

Liver X receptors (LXRs) are key regulators of lipid and cholesterol metabolism in mammals. Little is known, however, about the function and evolution of LXRs in non-mammalian species. The present study reports the cloning of LXRs from African clawed frog (Xenopus laevis), Western clawed frog (Xenopus tropicalis), and zebrafish (Danio rerio), and their functional characterization and comparison with human and mouse LXRs. Additionally, an ortholog of LXR in the chordate invertebrate Ciona intestinalis was cloned and functionally characterized. Ligand specificities of the frog and zebrafish LXRs were very similar to LXRalpha and LXRbeta from human and mouse. All vertebrate LXRs studied were activated robustly by the synthetic ligands T-0901317 and GW3965 and by a variety of oxysterols. In contrast, Ciona LXR was not activated by T-0901317 or GW3965 but was activated by a limited number of oxysterols, as well as some androstane and pregnane steroids. Pharmacophore analysis, homology modeling, and docking studies of Ciona LXR predict a receptor with a more restricted ligand-binding pocket and less intrinsic disorder in the ligand-binding domain compared to vertebrate LXRs. The results suggest that LXRs have a long evolutionary history, with vertebrate LXRs diverging from invertebrate LXRs in ligand specificity.     

Evolution of pharmacologic specificity in the pregnane X receptor

Background  

The pregnane X receptor (PXR) shows the highest degree of cross-species sequence diversity of any of the vertebrate nuclear hormone receptors. In this study, we determined the pharmacophores for activation of human, mouse, rat, rabbit, chicken, and zebrafish PXRs, using a common set of sixteen ligands. In addition, we compared in detail the selectivity of human and zebrafish PXRs for steroidal compounds and xenobiotics. The ligand activation properties of the Western clawed frog (Xenopus tropicalis) PXR and that of a putative vitamin D receptor (VDR)/PXR cloned in this study from the chordate invertebrate sea squirt (Ciona intestinalis) were also investigated.     

The cholesterol-raising factor from coffee beans, cafestol, as an agonist ligand for the farnesoid and pregnane X receptors

Cafestol, a diterpene present in unfiltered coffee brews such as Scandinavian boiled, Turkish, and cafetière coffee, is the most potent cholesterol-elevating compound known in the human diet. Several genes involved in cholesterol homeostasis have previously been shown to be targets of cafestol, including cholesterol 7alpha-hydroxylase (CYP7A1), the rate-limiting enzyme in bile acid biosynthesis. We have examined the mechanism by which cafestol elevates serum lipid levels. Changes in several lipid parameters were observed in cafestol-treated APOE3Leiden mice, including a significant increase in serum triglyceride levels. Microarray analysis of these mice identified alterations in hepatic expression of genes involved in lipid metabolism and detoxification, many of which are regulated by the nuclear hormone receptors farnesoid X receptor (FXR) and pregnane X receptor (PXR). Further studies demonstrate that cafestol is an agonist ligand for FXR and PXR, and that cafestol down-regulates expression of the bile acid homeostatic genes CYP7A1, sterol 12alpha-hydroxylase, and Na(+)-taurocholate cotransporting polypeptide in the liver of wild-type but not FXR null mice. Cafestol did not affect genes known to be up-regulated by FXR in the liver of wild-type mice, but did increase expression of the positive FXR-target genes intestinal bile acid-binding protein and fibroblast growth factor 15 (FGF15) in the intestine. Because FGF15 has recently been shown to function in an enterohepatic regulatory pathway to repress liver expression of bile acid homeostatic genes, its direct induction in the gut may account for indirect effects of cafestol on liver gene expression. PXR-dependent gene regulation of cytochrome P450 3A11 and other targets by cafestol was also only seen in the intestine. Using a double FXR/PXR knockout mouse model, we found that both receptors contribute to the cafestol-dependent induction of intestinal FGF15 gene expression. In conclusion, cafestol acts as an agonist ligand for both FXR and PXR, and this may contribute to its impact on cholesterol homeostasis.     

Irx4 forms an inhibitory complex with the vitamin D and retinoic X receptors to regulate cardiac chamber-specific slow MyHC3 expression

The slow myosin heavy chain 3 gene (slow MyHC3) is restricted in its expression to the atrial chambers of the heart. Understanding its regulation provides a basis for determination of the mechanisms controlling chamber-specific gene expression in heart development. The observed chamber distribution results from repression of slow MyHC3 gene expression in the ventricles. A binding site, the vitamin D response element (VDRE), for a heterodimer of vitamin D receptor (VDR) and retinoic X receptor alpha (RXR alpha) within the slow MyHC3 promoter mediates chamber-specific expression of the gene. Irx4, an Iroquois family homeobox gene whose expression is restricted to the ventricular chambers at all stages of development, inhibits AMHC1, the chick homolog of quail slow MyHC3, gene expression within developing ventricles. Repression of the slow MyHC3 gene in ventricular cardiomyocytes by Irx4 requires the VDRE. Unlike VDR and RXR alpha, Irx4 does not bind directly to the VDRE. Instead two-hybrid and co-immunoprecipitation assays show that Irx4 interacts with the RXR alpha component of the VDR/RXR alpha heterodimer and that the amino terminus of the Irx4 protein is required for its inhibitory action. These observations indicate that the mechanism of atrial chamber-specific expression requires the formation of an inhibitory protein complex composed of VDR, RXR alpha, and Irx4 that binds at the VDRE inhibiting slow MyHC3 expression in the ventricles.     

Activation of pregnane X receptor disrupts glucocorticoid and mineralocorticoid homeostasis

The pregnane X receptor (PXR) was isolated as a xenobiotic receptor that regulates responses to various xenobiotic agents. In this study, we show that PXR plays an important endobiotic role in adrenal steroid homeostasis. Activation of PXR by genetic (transgene) or pharmacological (ligand, such as rifampicin) markedly increased plasma concentrations of corticosterone and aldosterone, the respective primary glucocorticoid and mineralocorticoid in rodents. The increased levels of corticosterone and aldosterone were associated with activation of adrenal steroidogenic enzymes, including cytochrome P450 (CYP)11a1, CYP11b1, CYP11b2, and 3beta-hydroxysteroid dehydrogenase. The PXR-activating transgenic mice also exhibited hypertrophy of the adrenal cortex, loss of glucocorticoid circadian rhythm, and lack of glucocorticoid responses to psychogenic stress. Interestingly, the transgenic mice had normal pituitary secretion of ACTH and the corticosterone-suppressing effect of dexamethasone was intact, suggesting a functional hypothalamus-pituitary-adrenal axis despite a severe disruption of adrenal steroid homeostasis. The ACTH-independent hypercortisolism in the PXR-activating transgenic mice is reminiscent of the pseudo-Cushing's syndrome in patients. The glucocorticoid effect appears to be PXR specific, as the activation of constitutive androstane receptor in transgenic mice had little effect. We propose that PXR is a potential endocrine disrupting factor that may have broad implications in steroid homeostasis and drug-hormone interactions.     

Coactivator binding promotes the specific interaction between ligand and the pregnane X receptor

The pregnane X receptor (PXR) detects the presence of a wide variety of endogenous and xenobiotic compounds, and is a master regulator of the expression of genes central to drug metabolism and excretion. We present the 2.0A crystal structure of the human PXR ligand-binding domain (LBD) in complex with the cholesterol-lowering compound SR12813 and a 25 amino acid residue fragment of the human steroid receptor coactivator-1 (SRC-1) containing one LXXLL motif. PXR crystallizes as a homodimer in the asymmetric unit in this structure and possesses a novel alpha2 helix adjacent to its ligand-binding cavity. The SRC-1 peptide forms two distinct helices and binds adjacent to the ligand-dependent transactivation AF-2 helix on the surface of PXR. In contrast with previous PXR structures, in which SR12813 bound in multiple orientations, the small SR12813 agonist in this structure binds in a single, unique orientation within the receptor's ligand-binding pocket and contacts the AF-2 helix. Thermal denaturation studies reveal that the SR12813 ligand and SRC-1 coactivator peptide each stabilize the LBD of PXR, and that together they exert an additive effect on the stability of the receptor. These results indicate that the binding of coactivator to the surface of PXR limits the ability of this promiscuous receptor to "breathe" and helps to trap a single, active conformation of SR12813. They further reveal that specificity is required for PXR activation.     

Expression of vitamin D receptors in the superior cervical ganglia of rats

We investigated the role of vitamin D in the sympathetic nervous system including the distribution of vitamin D receptors (VDR), 1α-hydroxylase and 24-hydroxylase (CYP24) in neuronal subpopulations and satellite glia in the superior cervical ganglia (SCGs) of rats using immunohistochemistry. VDR immunoreactivity was observed in the cytoplasm and nucleus of nearly all neurons in the SCG. Intensity of VDR fluorescence was significantly greater in the cytoplasm of neuropeptide Y (NPY) negative somata compared to NPY positive neurons. Immunoreactivity for 1α-hydroxylase also was observed in the cytoplasm of all neurons of the SCG, but the intensity of fluorescence was less in the nuclei. To the contrary, the immunoreactivity for CYP24 was stronger in the nuclei, although it was present at lower intensity also in the cytoplasm of neurons. VDR and 1α-hydroxylase immunofluorescence was observed in many non-neuron cells, except satellite glial cells, which exhibited weak CYP24 immunofluorescence. Expression of VDRs and key metabolizing enzymes indicated the importance of vitamin D in the autonomic nervous system and the ability of sympathetic neurons to activate and deactivate vitamin D for its autocrine and paracrine roles.     

DNA bending is induced by binding of vitamin D receptor-retinoid X receptor heterodimers to vitamin D response elements.

The ability of vitamin D receptor-retinoid X receptor (VDR-RXR) heterodimers to induce a DNA bend upon binding to various vitamin D response elements (VDRE) has been investigated by circular permutation and phasing analysis. Recombinant rat VDR expressed in the baculovirus system and purified recombinant human RXR beta have been used. The VDREs were from 1,25-dihydroxyvitamin D3 (1,25-[OH]2D3) enhanced genes (rat osteocalcin, rOC; mouse osteopontin, mOP, and rat 1,25-dihydroxyvitamin D3-24-hydroxylase, r24-OHase), and a 1,25-(OH)2D3 repressed gene (human parathyroid hormone, hPTH). As shown by circular permutation analysis, VDR-RXR induced a distortion in DNA fragments containing various VDREs. Calculated distortion angles were similar in magnitude (57 degrees, 56 degrees, 61 degrees, and 59 degrees, respectively for rOC, mOP, r24-Ohase, and hPTH). The distortions took place with or without a 1,25-(OH)2D3 ligand. The centers of the apparent bend were found in the vicinity of the midpoint of all VDREs, except for rOC VDRE which was found 4 bp upstream. Phasing analysis was performed with DNA fragments containing mOP VDRE and revealed that VDR-RXR heterodimers induced a directed bend of 26 degrees, not influenced by the presence of hormone. In this study we report that similar to other members of the steroid and thyroid nuclear receptor superfamily, VDR-RXR heterodimers induce DNA bending.     

CYP3A4 and pregnane X receptor humanized mice

Marked species differences exist in P450 expression and activities. In order to produce mouse models that can be used to more accurately predict human drug and carcinogen metabolism, P450- and xenobiotic receptor humanized mice are being prepared using bacterial artificial chromosomes (BAC) and P1 phage artificial chromosomes (PAC) genomic clones. In some cases, transgenic mice carrying the human genes are bred with null-mice to produce fully humanized mice. Mice expressing human CYP1A1, CYP1A2, CYP2E1, CYP2D6, CYP3A4, and CYP3A7 were generated and characterized. Studies with the CYP3A4-humanized (hCYP3A4) mouse line revealed new information on the physiological function of this P450 and its role in drug metabolism in vivo. With this mouse line, CYP3A4, under certain circumstances, was found to alter the serum levels of estrogen resulting in deficient lactation and low pup survival as a result of underdeveloped mammary glands. This hCYP3A4 mouse established the importance of intestinal CYP3A4 in the pharmacokinetics of orally administered drugs. The hCYP3A4 mice were also used to establish the mechanisms of potential gender differences in CYP3A4 expression (adult female > adult male) that could account for human gender differences in drug metabolism and response. The pregnane X receptor (PXR) is also involved in induction of drug metabolism through its target genes including CYP3A4. Since species differences exist in ligand specificity between human and mice, a PXR-humanized mouse (hPXR) was produced that responds to human PXR activators such as rifampicin but does not respond to the rodent activator pregnenalone 16alpha-carbonitrile.     

Are the receptors of 1,25-dihydroxyvitamin D3 vitamin K dependent?

Alimentary deficiency of vitamin K in rats causes a decrease in the level of in vivo occupied nuclear 1,25 (OH)2D3 receptors in small intestinal mucosa and an 2-2.5-fold increase in the ability of cytosolic 1,25 (OH)2D3-receptor complexes to bind to heterologous DNA. The 1,25 (OH)2D3 binding by the receptors is thereby unaffected. Preincubation of kidney and intestinal cytosol of rats with the secondary K-avitaminosis induced by vitamin K antagonist with the microsomal vitamin K-dependent gamma-carboxylation system sharply decreases the binding of the 1.25 (OH)2D3-receptor complexes to DNA. In rats treated with the vitamin K antagonist in combination with a low calcium diet, the subsequent maintenance on a high calcium diet does not cause, in contrast with vitamin K-repleted animals, a sharp decrease of the level of the in vivo occupied 1,25 (OH)2D3 receptors. In vitro Ca2+ cations decrease the binding of the 1,25 (OH)2D3-receptor complexes to DNA only in vitamin K-repleted rats (ED50 = 2.5 x 10(-6) M). The existence of a vitamin K-dependent Ca-sensitive mechanism regulating the binding of the 1,25 (OH)2D3 receptor to DNA has been postulated for the first time.     

Functional characterization and evolution of PTH/PTHrP receptors: insights from the chicken

ABSTRACT: BACKGROUND: The parathyroid hormone (PTH)-family consists of a group of structurally related factors that regulate calcium and bone homeostasis and are also involved in development of organs such as the heart, mammary gland and immune system. They interact with specific members of family 2 B1 G-protein coupled receptors (GPCRs), which have been characterised in teleosts and mammals. Two PTH/PTHrP receptors, PTH1R and PTH2R exist in mammals and in teleost fish a further receptor PTH3R has also been identified. Recently in chicken, PTHfamily members involved in calcium transport were characterized and specific PTHRs are suggested to exist although they have not yet been isolated or functionally characterized. The aim of this study is to further explore the evolution and function of the vertebrate PTH/PTHrP system through the isolation, phylogenetic analysis and functional characterization of the chicken receptors. RESULTS: Two PTHRs were isolated in chicken and sequence comparison and phylogenetic analysis indicate that the chicken receptors correspond to PTH1R and PTH3R, which emerged prior to the teleost/tetrapod divergence since they are present in cartilaginous fish. The vertebrate PTH2R receptor and its ligand TIP39 have been lost from bird genomes. Chicken PTH1R and PTH3R have a divergent and widespread tissue expression and are also evident in very early embryonic stages of development. Receptor stimulation studies using HEK293 cells stably expressing the chicken PTH1R and PTH3R and monitoring cAMP production revealed they are activated by chicken 1-34 N-terminal PTH-family peptides in a dose dependent manner. PTH-L and PTHrP were the most effective peptides in activating PTH1R (EC50 = 7.7 nM and EC50 = 22.7 nM, respectively). In contrast, PTH-L (100 nM) produced a small cAMP accumulation on activation of PTH3R but PTHrP and PTH (EC50 = 2.5 nM and EC50 = 22.1 nM, respectively) readily activated the receptor. PTHrP also stimulated intracellular Ca2+ accumulation on activation of PTH1R but not PTH3R. CONCLUSION: Two PTHR homologues of the vertebrate PTH1R and PTH3R were isolated and functionally characterized in chicken. Their distinct pattern of expression during embryo development and in adult tissues, together with their ligand preference, suggests that they have acquired specific functions, which have contributed to their maintenance in the genome. PTH2R and its activating ligand, TIP39, are absent from bird genomes. Nonetheless identification of putative PTH2R and TIP39 in the genome of an ancient agnathan, lamprey, suggests the PTH/PTHrP ligand and receptor family was already present in an early basal paraphyletic group of vertebrates and during the vertebrate radiation diverged via gene/genome duplication and deletion events. Knowledge of the role PTH/PTHrP system in early vertebrates will help to establish evolution of function.