Center for Genomics and Bioinformatics (A.S., W.W.), Karolinska Institutet, 113 37 Stockholm, Sweden; Centre for Molecular Medicine and Therapeutics (W.W.W.), British Columbia Childrens and Womens Hospitals, and Department of Medical Genetics (W.W.W.), University of British Columbia, Vancouver, Canada V5Z 4H4
Address all correspondence and requests for reprints to: Wyeth W. Wasserman, Centre for Molecular Medicine and Therapeutics, 3018-950 West 28th Avenue, Vancouver, British Columbia, Canada V5Z 4H4. E-mail: wyeth{at}cmmt.ubc.ca.
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The nuclear receptor (NR) family constitutes a medically important group of TFs, involved in critical biological processes in all multicellular metazoans (4). The NR proteins are unique in directly linking regulation of gene transcription to the concentrations of key physiological signaling molecules, such as estrogen, testosterone, and glucocorticoids (4, 5). Thus, NRs constitute established drug targets for combating disease (6), for instance in the treatment of diabetes [peroxisome proliferator-activated receptor (PPAR) ] (7), cardiovascular disease (estrogen receptor, PPAR
, and ligand X receptor) (8, 9, 10), and osteoporosis and autoimmune diseases (vitamin D receptor) (11). NR proteins are composed, in most cases, of several domains, including a DNA binding domain, a ligand-binding domain, and one or more activation domains (12). With certain exceptions, NRs operate as dimerseither as homodimers or in heterodimeric combinations (4, 5). The members of the NR family have been categorized based on function and evolution. The classification scheme suggested by Laudet et al. (13) divides NRs into seven groups (NR0-NR6). Whereas most classes exhibit similar properties in DNA binding specificity, the NR0 and NR3 groups deserve special treatment. The NR0 group describes unusual NRs lacking either the ligand binding domain or DNA binding domain (13, 14). Because both domains are required for DNA binding by dimerization, these factors either cannot bind DNA or bind as monomers (15). The NR3 group (also known as nuclear steroid receptors) differs significantly in terms of both DNA binding specificity and functional ligands (16). The binding distinction of the NR0 and NR3 subclasses is not further addressed in this work. The NR12 and NR46 subclasses are collectively referred to as nuclear hormone receptors (NHRs). NHR dimers bind to regulatory sequences composed of two half-sites that are separated by variable spacing and can occur in different orientations (12, 17). Half-sites have the consensus sequence AGGTCA, but most factors tolerate considerable variation from the consensus within functional binding sites. Combinations include the direct repeat (DR) AGGTCAnxAGGTCA (or the reverse complement TGACCTnxTGACCT), the inverted repeat (IR) AGGTCAnxTGACCT, and the everted repeat (ER) TGACCTnxAGGTCA (5), where the number of spacer nucleotides nx can vary from zero to eight nucleotides. For clarity, when we refer to a specific configuration, the repeat abbreviation (DR, ER, or IR), is followed by the number of spacers between half-sites (for example, a direct repeat with a 5-bp spacer is referred to as DR5). A few NHRs (e.g. retinoic acid receptor-related orphan receptor, nerve growth factor-induced clone B, and steroidogenic factor-1) are known to function in vivo as monomers, where the single half-site has an adenosine/thymidine-rich 5' extension (4). This situation is not addressed in this work because single half-sites can be adequately modeled with well-established profile models termed position-weight matrices.
![]() |
Computational Prediction of NHR Target Sites |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
HMMs are widely used in bioinformatics in problems ranging from the prediction of protein domains to the analysis of gene structure (23, 24, 25, 26). Briefly, an HMM model consists of a set of states, where each state can emit symbols (nucleotides in this case) based on a probability distribution. The emission probability for a certain nucleotide is specific for each state. States are connected in a chain-like structure, where the probability of moving from one state to another is termed a transition probability. Any specific path through the states for a given sequence will have a defined probability (effectively the product of all emission probabilities for respective nucleotides and all transitions probabilities for each move between states). It is possible to computationally calculate the optimal route through the chain for a given sequence (the Viterbi algorithm), or the total probability of the model producing a certain sequence (the Forward algorithm). Excellent textbooks and reviews are available discussing HMMs in detail (23, 24).
Based on an extensive collection of verified binding sites, we have developed an HMM model capable of simultaneous recognition and classification of NHR binding sites. Analysis of the models performance indicates high sensitivity for classification of real sites, while maintaining selectivity comparable to weight matrix binding profiles for TFs outside the NR class. Application of the model to the genome of the pufferfish (Fugu rubripes) (27), identifies a high cross-regulatory potential between NR encoding genes. The bioinformatics method is implemented in a user-friendly web interface, available at http://mordor.cgb.ki.se/NHR-scan.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Model Architecture
The model architecture was chosen based on the biochemical characteristics of NHR binding and the volume of training data. Briefly, the model is a parallel classification-type HMM combined with circular local alignment architecture (23) (Fig. 1). From the background state B, it is possible to move to each type of match state: DR, ER, and IR, of which each consists of two half-site models separated by a spacer-state classifier. Direct repeats can be viewed in two directions; thus, a reverse complement version is required for an accurate description. In the case of ER and IR sites, one must consider the binding site characteristics closely. It would not be illogical to argue that a reverse complement model is duplicative, based on the palindromic pattern. However, imperfect palindromes are frequently observed, resulting in a de facto orientation. Therefore, for each repeat class, we allow for two match states; one in each strand direction. As the default setting, the combined probability j of entering all match states is internally divided equally among the six alternatives (three structures on each strand). This parameter can be adjusted.
|
Model Parameter Estimation
Advanced machine learning algorithms, such as the Baum-Welch algorithm (28), are often employed for estimation of model parameters (i.e. emission and transition probabilities). However, if each nucleotide in the input sequences is labeled with respective state, emission parameter estimation is produced by simply counting each nucleotide found in each state in the input data (a maximum likelihood procedure). Transition probabilities can be estimated analogously (23). Because all nucleotides in the input sequences can be labeled based on the literature-derived data, training of emission and transmission probabilities is straightforwardcomparable to constructing common binding site profiles by counting nucleotides in each column of an alignment. A pseudocount is often added to emission and transition counts to avoid overfitting the model to small data sets. We used a Laplaces rule pseudocount procedure (i.e. in each state, add one to the count of observed nucleotides) (23) because the background nucleotide distribution in promoters is approximately uniform as judged from the Eukaryotic Promoter Database (29).
Because DNA strands can be described in two directions, and it is not given that all experimentally confirmed sites are reported in the same strand, we used a Gibbs Sampling-based alignment algorithm [ANN-Spec (30)] to orient each reference site within each class (DR, IR, and ER). Because ANN-Spec does not produce gapped alignments, concatenated half-sites were aligned (i.e. the regions of variable spacing were removed). Gap state emissions (but not transitions) were estimated from a pool consisting of all gap-labeled nucleotides to avoid overfitting sparsely populated states.
Scoring Algorithms
To identify and classify candidate binding sites, the Viterbi algorithm was used to find the most probable state path through the model with a given sequence. Given the model, the Viterbi algorithm is applied to identify the most probable chain of states (labels) that are consistent with the observed sequence (23). A disadvantage with the Viterbi algorithm is that predicted binding sites cannot overlap. The alternative is to use the Forward algorithm on a fixed window that slides along the input sequence. The latter approach is computationally expensive but produces a probability score for the given subsequence and the model, corresponding to prediction scores using profile models. We have implemented both methods in the web application, allowing the user to select between the faster Viterbi approach for long sequences and the slower Forward approach for detailed analysis.
Model Evaluation
Most models for binding site prediction are limited by inadequate supplies of confirmed binding sites for training and verification. We have evaluated the model sensitivity with the standard cross-validation (jack-knife test) procedure and using two additional independent data sets (Fig. 1, A and B), as described in Materials and Methods.
In a cross-validation test, one binding site is removed from the training set and subsequently assessed with the model created with the remaining sites. The procedure is repeated iteratively to cover every site in the training set. In any sensitivity test for an HMM, the outcome is strongly dependent on the probability j of entering match states initially. In the cross-validation test (Fig. 1A), a plateau is reached at approximately 85% correctly labeled nucleotides for DR and ER sites and approximately 75% for IR sites, whereafter less stringent values of j(>0.04) do not improve sensitivity substantially.
To verify that the model is not overfitted to the original training data, we tested the method on two additional large data sets: response elements identified from in vitro selection studies (31, 32) and literature (33), respectively (Fig. 2B). The large majority of sites in both sets can be identified with low settings of j. Although almost all the in vitro sites are identified at stringent thresholds, the literature-derived sites are more varied (detectable with less-stringent settings). Nevertheless, the scores for the binding sites in both test sets are clearly distinct from the distribution of scores observed for the sequences in the randomly generated negative set (Fig. 1B
). As in the previous test, a sensitivity plateau is reached after j > 0.04.
|
The Limitations of Phylogenetic Footprinting for Detection of NHR Response Elements
TFBS prediction is limited by poor selectivity (i.e. most predicted sites are not functional) (19, 34), due to the limited size of binding sites and the considerable sequence variation in binding sites tolerated by most TFs (reviewed in Refs. 18 and 19).
Cross-species comparison (phylogenetic footprinting) is a viable method for addressing the selectivity problem (2). Such approaches involve the identification and alignment of promoter regions for orthologous genes in two separate species (for instance, human and mouse), followed by analysis with TF binding profiles within preferentially conserved regions. The success of alignment-based phylogenetic footprinting is limited by two important alignment problems:
1) Sites have to be situated in regions that can be readily aligned. Regulatory regions located far from promoters (distal regions) tend to lack long continuous segments of conservation that are required for detection by global alignment algorithms, even if the binding sites in themselves are conserved. This is an important problem for the analysis of NHR TFBS because many of the sites are distally positioned (35, 36, 37).
2) Although regulatory properties of essential systems can remain intact over long periods of evolution, alignment-based phylogenetic footprinting lends itself best to moderate evolutionary distances such as that between human and mouse.
Alignment-based analysis of gene promoters separated by longer evolutionary distances is often problematic because alignment methods often are unable to generate reliable alignments over the entire regulatory sequence due to the low level of sequence similarity (2).
Analysis of NHR Binding Sites in the F. rubripes Genome
The pufferfish F. rubripes (27) genome is a valuable reference, both for its compact size (450 Mb) and its phylogenetic distance from humans (450 My). The F. rubripes and human genomes display remarkable similarity in the content and order of genes (38). In diverse studies, NR genes, in particular the ligand-binding and DNA-binding region, are remarkably well conserved over evolution in vertebrates (6, 39, 40). Successful regulatory analyses have been enabled by sequence comparisons between F. rubripes and mammals (41, 42, 43), including cases that identified NR binding sites (40).
Because the genes encoding NRs are well represented in F. rubripes (44), and in light of the importance of NRs in the regulation of critical genes, the F. rubripes genome is a powerful platform for computational prediction of regulatory elements.
To motivate the use of the compact F. rubripes genome as a platform for prediction of NHR binding sites, we explored whether genes regulated by NHRs in man are likely to be similarly regulated in the pufferfish. We investigated whether a binding site configuration (e.g. DR1) active in human and rodent genes is more abundant than expected by chance in the promoters of homologous F. rubripes genes (as described in Materials and Methods). Homologs to NHR-regulated human and rodent genes were identified in the F. rubripes based on EnsEMBL orthology mappings. DR1 sites are significantly overrepresented in these F. rubripes promoters (P = 0.0015, cumulative binomial distribution (as explained in Materials and Methods), at cutoff 0.02). DR3, DR5, ER1, and ER6 configured sites were also significantly overrepresented (P < 0.05). The other site configurations, which are not populous in the reference collection, were not found to be significantly overrepresented.
Case Example: Evidence for Cross/Autoregulation of NHR Genes
Numerous NR genes in mammals are autoregulated or cross-regulated by other NRs (45, 46, 47, 48, 49). To demonstrate an application of the NHR binding site prediction model in a genome-scale application, we investigated the properties of cross/autoregulation. As discussed above, for genome analysis, the compact F. rubripes genome is a preferred target for analysis.
Occurrence of NHR binding sites in upstream sequences of NR genes was compared with a background promoter set. Despite the inclusion of steroid NR genes (whose binding specificity is not addressed by the model) and the possible presence of atypical NRs (for instance, pseudogenes), we observed a clear overrepresentation (P < 0.05, cumulative binomial distribution) of certain site configurations: most notably the direct repeats DR1, DR3, DR4, DR5, and DR8. In a further analysis step, we investigated whether the six subclasses (NR16) of NRs had different cross/autoregulatory potential (Fig. 3). The same statistical analysis was performed for the six different classes of NRs found in the set. In general, cross/autoregulatory potential remains high in most groups; however, certain differences in site distributions are evident. The number of NR1, NR2, and NR3 class genes is sufficient to support analysis of site frequency (Fig. 3A
). Although the regulatory potential of NR1 class is composed of a wide range of site configurations (DR4, DR5, DR8, ER1, and ER6 are all significantly overrepresented), the NR2 class clearly has a higher tendency for direct repeats (DR1 and DR3 are significantly overrepresented) (Fig. 3B
). Interestingly, the NR3 (steroid receptors) class has a defined cross-regulatory potential, with an abundance of direct repeats that are unlikely to be bound by NR3 members (DR1 and DR3 are highly overrepresented).
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Performance of the HMM Model Framework
Sensitivity
NHR-scan accurately describes the properties of NR binding sites. Even a small probability j of entering the match states confers high predictive sensitivity (Fig. 2A). At higher settings of j (>0.1) (i.e. less stringent), more than 90% of reference binding sites are identified, although some are mislabeled (data not shown). The latter is possibly due to 1) the few base pair changes required to transfer a weak site of one configuration type to another, and 2) potentially erroneous classifications of sites in the literature. Misclassifications were only observed for IR and ER sites, which, although consistent with the greater variation tolerated in these site formats, is most likely due the limited number of validated sites. The IR and ER site collections are each 12 sites deep, compared with 83 DR sites. The classification and detection of IR sites are less reliable than those produced by the DR and ER classifiers (Fig. 2A
). This is entirely consistent with the sequence variability observed in one of the half-sites in the IR modelthe IR sites are less informative than the others (Fig. 1
). Because the data depth is low, we cannot estimate whether this is a sampling bias or whether it reflects a general property of IR sites. Clearly, both IR and ER models would benefit from incorporation of more functionally validated sites. Despite the limited data, the sensitivity tests using independent data sets confirm that the model is not overfitted.
Selectivity
The frequency of predicted binding sites is directly dependent on the user-specified probability of entering the match states. For predictive purposes, a structure for a NR binding site (i.e. a given spacing and orientation) is highly specific and thus generates few false predictions. Because all known half-site architectures are allowed by NHR-scan, predictions will be more frequent than those produced by a rigid model. For example, at threshold j = 0.01 we predict, on average, one site approximately every 2000 bp, while classifying approximately 80% of the labels correctly (Fig. 2A). Although the selectivity is higher than a corresponding analysis with a typical matrix model for TF binding sites (which produce predictions at an average rate of
1/500 bp), the model is in itself not sufficiently selective for meaningful analysis of entire genomes.
F. rubripes NRs Have a High Cross-Regulatory Potential
We have identified site configurations that occur preferentially in F. rubripes NR promoters, indicating a considerable potential for cross/autoregulation of these genes. Furthermore, classes of receptors (13) exhibit differential preferences in terms of site configurationsmost notably difference between the NR1, NR2, and NR3 classes. Because the NR3 class (steroid receptors) is functionally distinct from NR1 and NR2 (both in terms of ligand affinities and binding site preferences), it is not surprising that the regulatory sequences display distinct characteristics. The NR3 class proteins preferentially bind a distinct site configuration, consisting of an IR3 site with a different half-site consensus (16). It is possible that some of the detected DR3 sites are nuclear steroid receptor binding sites; however, initial tests of the NHR-scan model on 43 validated nuclear steroid receptor sites at the same settings only detected six sites (data not shown). Supporting the hypothesis, several findings indicate that the different classes of receptors often are involved in cross-signaling (51, 52). For instance, in human, thyroid hormone receptor ß is known to be autoregulated by a proximal DR4 response element (47). Interestingly, thyroid hormone ß receptor (here labeled THB2_HUMAN) has one of the highest potentials for cross/autoregulation within the test set. In addition, recent studies indicate that the thyroid receptor ß gene is cross-regulated by retinoid X receptor (47). We interpret the abundance of NHR binding sites to indicate potential for considerable cross/autoregulatory networks within the NR genes.
Comparison to Similar Methods
The only comparable published resource for NHR TFBS prediction is NUBI-scan (53). This service is based on an algorithm that scans sequences with a single half-site model [(using profile-based prediction methodology (18)] and subsequently evaluates pairs of half-sites with a chosen configuration (for example DR4) by the product of half-site scores. Thus, the user must directly specify possible site configurations for analysis to proceed. Large-scale sensitivity and selectivity performance data has not been published for NUBI-scan. Several characteristics distinguish NHR-scan from NUBI-scan, including the depth of the training data (the NUBI-scan model is based on 11 reference sites), the model framework, and the assessment of predictions.
In addition to obvious algorithmic differences, the statistical methods for evaluating hits are different. NHR-scan uses established HMM scoring algorithms, whereas NUBI-scan reports predictions exceeding 67 Z scores. These Z scores are derived by comparison to a distribution of all scores generated by the input sequence; NUBI-scan results are dependent on the length and composition of the sequence analyzed
Roelet et al. (54) have developed a prediction method for the CCAAT-binding transcription factor/nuclear factor I transcription factor that binds DNA with two half-sites with a variable spacing (but not orientation), based on a generalized profile concept (54). Although not stated as a HMM, the algorithm has certain similarities to NHR-scan but is focused on another biological problem.
Future Directions
Models for Each Type of NR Class (or Individual Factors).
NHR-scan models the binding sites for all NHRs and is thus classified as a meta-model. To model individual factors, or subclasses of NRs, would be a logical extension of the present work. However, the body of experimentally validated in vivo sites is insufficient for most factors. Promising developments are underway to increase the collections of known sites, both through laboratory methods for chromatin immunoprecipitation of TFs bound to active sites (55) and through systematic evaluations of sites in vitro (56).
Input Sequence Weighting.
Because the binding sites in the reference collection cannot be regarded as fully representative (for example, the collection contains nine ligand X receptor sites but only three testicular orphan receptor-4 sites), implanting a system to weight the contribution of each site could have utility. Such a system could be based on the similarity between the binding domains of the NHRs or between target genes and modified by the number of sites available. It is hard to evaluate weighting procedures given the data limitations. Because a weighting system would produce a more generalized model, sensitivity is likely to decline, and we would have no means of measuring potential benefits.
Regulatory Modules.
The focus on isolated binding sites of a single TF is biologically unrealistic. This is reflected in the limited selectivity of binding site predictions. In the recent years, cis-regulatory module detection (i.e. detection of functional clusters of binding sites) has proven a potent method for delineating regulatory regions on a genome scale (57, 58). NHR-scan can conceptually be integrated in the mathematical frameworks employed for such module discrimination (which usually are based on profile representations). For example, the liver regulatory model by Krivan and Wasserman (22), which includes a profile model for the HNF4 NHR, would likely be enhanced with an HMM model describing the variable spacing of HNF4 binding sites.
Background Model.
Background models for describing sequence characteristics can dramatically impact the performance of predictive algorithms. NHR-scan presently uses a simple nucleotide distribution. It is known that the intricacy of genomic DNA requires additional model complexity for rigorous representations, for instance as an n-state Markov chain (30, 59). This concept could be extended to include different background Markov chains for various genomic features: repeat regions, coding regions, and regions with regulatory potential (60).
Conclusion
We have presented a HMM framework for predicting NHR binding sites in genomic sequences. The model is sensitive, while maintaining respectable selectivity. A genome-scale application of the model highlighted the high cross-regulatory potential of F. rubripes NRs. The NHR-scan model provided through a user-friendly Internet interface will positively impact efforts to identify functional NHR binding sites.
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Algorithms
For scoring purposes, we used the general HMM (GHMM) C++ library (http://ghmm.org/), where the model is described as an XML object using the XMLIO library (http://xmlio.sourceforge.net/). The combination of libraries enables efficient central processing unit usage and human-readable models and outputs. It should be noted that the model can be used by any HMM software package that allows for custom state architectures.
Cross-Comparison (Jack-Knife) Assessment of Sensitivity
Both sensitivity and selectivity are dependent on both the quality of the individual match states and the probability of entering a match state. For these reasons, we performed the test with varying probabilities of entering match states. The value j describes the combined probability of entering any match state; thus, the probability of entering any of the match chain is j/6 because each match state is composed of a forward- and a reverse-strand chain. We measured performance as the fraction of correctly labeled nucleotides (background, DR, IR, ER, or spacer).
Assessment of Sensitivity Using Independent Data Sets
We collected two sets of independent binding site data: 1) 134 in vitro selected binding sites from the JASPAR database (31) and other sources (32); and 2) 69 known HNF-4 binding sites from literature sources, collected by Ellrot et al. (33),
These data sets consists of binding sites of length approximately 1520 bp. Unlike the training set, we have no information about the exact positions of binding sites, and thus the half-site configuration for a given site is not known. Thus, we evaluate sensitivity by measuring how many of the sequences are predicted to contain a site, given a cutoff j. Each site was extended with three randomly generated nucleotides because the model for technical reasons requires that the first and last nucleotide in any analyzed sequence must be labeled background.
Selectivity Assessment
For selectivity assessment, we analyzed the vertebrate representative subset of the Eukaryotic Promoter Database (29) for NHR sites, using different thresholds as above. Because most promoters are not expected to harbor a functional NHR site, we measure selectivity as the frequency with which we predict predicted sites (average number of base pairs between two sites).
Coregulation F. rubripes/Human and Rodent
Homologs to human and rodent NHR-targeted genes were located in F. rubripes genome (assembly 17.2.1) using Swissprot (61) identifiers and EnsEMBL (50) gene annotations. For each gene, we retrieved F. rubripes promoter sequence corresponding to 1700 to +300 relative to the predicted transcription start site. For an F. rubripes promoter to be classified as positive, we required that one or more sites of the same type and half-site configuration as found in the reference human/rodent gene were predicted, using a set cutoff. For clarity, this measure does not correspond to the quantity of sites in a single promoter because a promoter can only be positive (i.e. containing one or more sites) or negative (containing no sites). With a given cutoff j, for each site configuration, a certain number of positives are reported out of the number of genes suspected to be regulated by that particular site type. For example, nine genes in the reference set (human/rodent) are regulated by DR1 sites. Setting the cutoff j to 0.02, we classify seven of the corresponding nine F. rubripes promoters as positive.
We estimated the chance of classifying a single promoter as positive by analyzing nonoverlapping promoters of the same size from the vertebrate representative subset of the Eukaryotic Promoter Database (29) (in total 1,072,239 bp). With such an analysis, we have an assessment of the chance p of randomly drawing a positive promoter using a certain type of site and configuration. The chance of drawing r positive sites out of n trials (corresponding to the number of reference genes regulated by the same site configuration) follows a binomial distribution:
![]() |
![]() |
NR Cross/Autoregulatory Potential
Using EnsEMBL (50) annotation, we retrieved the 57 genes containing the characteristic C4H4 NR DNA binding domain [INTERPRO (62) domain IPR001628]. This is comparable to the 68 genes identified using protein similarity searches in another study (44). As shown previously, we analyzed the region 1700 to +300 relating to the TSS, using the same statistical method. We divided the NRs into subclasses (as defined by Ref. 13), based on annotation or sequence similarity [(Blastp (63)]. Four of the 57 sequences could not be classified reliably. No promoter pair displayed any noteworthy similarity, based on ClustalW (64) analysis. Additional data and resources, including color versions of images, training data, and model framework in GHMM XML format, are available at the NHR-scan web page http://mordor.cgb.ki.se/NHR-scan.
![]() |
ACKNOWLEDGMENTS |
---|
![]() |
FOOTNOTES |
---|
First Published Online November 24, 2004
Abbreviations: DR, Direct repeat; ER, everted repeat; GHMM, general HMM; HMM, Hidden Markov Model; IR, inverted repeat; NHR, nuclear hormone receptor; NR, nuclear receptor; PPAR, peroxisome proliferator-activated receptor; TF, transcription factor; TFBS, transcription factor binding site.
Received for publication March 9, 2004. Accepted for publication November 17, 2004.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |