Comprehensive study of nuclear receptor DNA binding provides a revised framework for understanding receptor specificity

Drug-induced cis-regulatory elements in human hepatocytes affect molecular phenotypes associated with adverse reactions

Genomic variation drives phenotypic diversity, including individual differences in drug response. While coding polymorphisms linked to drug efficacy and adverse reactions are well characterized, the contribution of noncoding regulatory elements remains underexplored. Using CAGE (Cap Analysis of Gene Expression), profiling transcription initiations of mRNAs and enhancer RNAs, we identify candidate cis-regulatory elements (CREs) and assessed their activities simultaneously in HepG2 cells expressing the drug-responsive transcription factor pregnane X receptor (PXR). Comparison with GWAS data reveals strong enrichment of the drug-induced CREs near variants associated with bilirubin and vitamin D levels. Among those bound by PXR in primary hepatocytes, we identify enhancers of UGT1A1, TSKU, and CYP24A1 and functional alleles that alter regulatory activities. We also find that TSKU influences expression of vitamin D-metabolizing enzymes. This study expands the landscape of PXR-mediated regulatory elements and uncovers noncoding variants impacting drug response, providing insights into the genomic basis of adverse drug reactions.

Genomic regions occupied by both RARα and VDR are involved in the convergence and cooperation of retinoid and vitamin D signaling pathways

Retinoic acid receptors (RARs) and the vitamin D receptor (VDR) regulate distinct but overlapping gene sets in multiple cell types. The abundance and characteristics of regulatory regions, occupied by both RARs and VDR are largely unexplored. We used global approaches (ChIP-seq, RNA-seq, and ATAC-seq) and bioinformatics tools to map and characterize common binding regions of RARα and VDR in differentiated human THP-1 cells. We found that the cistromes of ligand-activated RARα and VDR largely overlapped, and their agonists (AM580 and calcitriol) co-regulated several genes, often cooperatively. Common binding regions were frequently (but not exclusively) annotated with co-regulated genes and exhibited increased MED1 occupancy upon ligand stimulation, suggesting their involvement in gene regulation. Chromatin accessibility was typically higher in the common regions than in regions occupied exclusively by RARα or VDR. DNA response elements for RARα (DR1/2/5) and VDR (DR3) were enriched in the common regions, albeit the co-occurrence of the two types of canonical motifs was low (8.4%), suggesting that "degenerate" DR1/2/5 and DR3 motifs or other sequences could mediate the binding. In summary, common binding regions of RARα and VDR are at the crossroads of the retinoid and vitamin D pathways, playing important roles in their convergence and cooperation.

Coupled Heterogeneity to Dimeric Site-Specific Binding by the POU-Family Transcription Factor OCT2

POU-family transcription factors regulate metazoan gene expression via a bipartite DNA-binding domain consisting of two covalently linked helix-turn-helix subdomains, POUS and POUH. POU factors bind as dimers to DNA half-sites to form complexes with a variable quaternary structure. To enhance the knowledge of the physical chemistry of dimeric POU/DNA recognition, we carried out a crystallographic and titration analysis of the cooperative homodimer formed by the POU factor OCT2 and an optimized palindromic DNA site known as MORE. The data evidence strong heterogeneity in the binding and formation of secondary complexes in site-specific DNA recognition by OCT2 at thermodynamic equilibrium. These secondary complexes are strictly contingent to the site-specific complex, detectable at subsaturating OCT2 concentrations, and cooperate with nonspecific binding to guide the affinity of the site-specific complex. Modulation with salt and poly[d(I-C)] unmasks the compensation between nonspecific DNA depleting unbound OCT2 on the one hand while driving specific binding by intermolecular transfer of OCT2 via secondary complexes on the other. Molecular dynamics simulations extend a mechanism, previously proposed for POU monomers, in which the two subdomains dynamically cross-link DNA strands to form supramolecular dimeric POU/DNA complexes at equilibrium.

The farnesoid X receptor activates transcription independently of RXR at non-canonical response elements

The farnesoid X receptor (FXR) is a nuclear receptor (NR) known to obligately heterodimerize with the retinoid X receptor (RXR). FXR is expressed as four isoforms (α1-α4) that drive transcription from IR-1 (inverted repeat-1) response elements (REs). Recently, we found that FXR isoforms α2/α4 also activate transcription from non-canonical ER-2 (everted repeat-2) REs, mediating most metabolic effects of general FXR activation. Here, we explored molecular determinants of regulation by FXRα2 from ER-2 REs through quantitative interaction proteomics, site-directed mutagenesis and transcriptomics. We discovered FXRα2 binds to and activates ER-2 elements in vitro and in reporter assays independently of RXR. Genome-wide binding analysis in mouse liver revealed higher ER-2 motif enrichment in FXR sites lacking RXR. Abrogation of FXRα2:RXR heterodimerization abolished IR-1, but preserved ER-2 transactivation. Transcriptome-wide, RXR overexpression inhibited 25% of FXRα2 targets in HepG2. These genes were specifically activated by the heterodimerization-deficient mutant FXRα2L434R, enriched for ER-2 motifs at their promoters, and involved in lipid metabolism and ammonia detoxification. In conclusion, RXR acts as a molecular switch, inhibiting FXRα2 activation from ER-2 while enhancing it from canonical IR-1 REs. Our results showcase FXR as the first NR with isoform-specific RXR-independent REs, highlighting a new layer of regulation and complexity for RXR-heterodimerizing NRs.

Hormone response elements for the thyroid receptor-α include specific distal 5'-flanking DNA

Optimal gene transcription is achieved through precise interactions between transcription factors and their DNA binding sites. We provide evidence that conserved distally located 5'-flanking sequences interact directly with the intrinsically disordered amino-terminal region of the thyroid receptor-α (TRα) to control transcriptional activity. Simulated modeling and dynamics with multiple ChIP-seq-derived sequences consistently reveal specific lysine/arginine-DNA minor groove interactions. The impact of these interactions is to distort DNA structural conformations, and these are also revealed with atomic force microscopy. The importance of the 5'-flanking DNA is further emphasized with reporter gene assays and comparisons with canonical response elements. Overall, the study reveals the inadequacy of current definitions of the DNA hormone response element (HRE) and suggests that future descriptions of the HRE include the conserved distal DNA sequences. The broad impact of this study is further underscored by the common occurrence of Lys/Arg-rich motifs within the intrinsically disordered regions of nuclear receptors.

Nr5a2 is dispensable for zygotic genome activation but essential for morula development

Early embryogenesis is driven by transcription factors (TFs) that first activate the zygotic genome and then specify the lineages constituting the blastocyst. Although the TFs specifying the blastocyst's lineages are well characterized, those playing earlier roles remain poorly defined. Using mouse models of the TF Nr5a2, we show that Nr5a2-/- embryos arrest at the early morula stage and exhibit altered lineage specification, frequent mitotic failure, and substantial chromosome segregation defects. Although NR5A2 plays a minor but measurable role during zygotic genome activation, it predominantly acts as a master regulator at the eight-cell stage, controlling expression of lineage-specifying TFs and genes involved in mitosis, telomere maintenance, and DNA repair. We conclude that NR5A2 coordinates proliferation, genome stability, and lineage specification to ensure correct morula development.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here, we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced histone 3 lysine 27 acetylation (H3K27ac). Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data support clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

LXR-dependent enhancer activation regulates the temporal organization of the liver's response to refeeding leading to lipogenic gene overshoot

Transitions between the fed and fasted state are common in mammals. The liver orchestrates adaptive responses to feeding/fasting by transcriptionally regulating metabolic pathways of energy usage and storage. Transcriptional and enhancer dynamics following cessation of fasting (refeeding) have not been explored. We examined the transcriptional and chromatin events occurring upon refeeding in mice, including kinetic behavior and molecular drivers. We found that the refeeding response is temporally organized with the early response focused on ramping up protein translation while the later stages of refeeding drive a bifurcated lipid synthesis program. While both the cholesterol biosynthesis and lipogenesis pathways were inhibited during fasting, most cholesterol biosynthesis genes returned to their basal levels upon refeeding while most lipogenesis genes markedly overshoot above pre-fasting levels. Gene knockout, enhancer dynamics, and ChIP-seq analyses revealed that lipogenic gene overshoot is dictated by LXRα. These findings from unbiased analyses unravel the mechanism behind the long-known phenomenon of refeeding fat overshoot.

Long Non-Coding RNAs, Nuclear Receptors and Their Cross-Talks in Cancer-Implications and Perspectives

Long non-coding RNAs (lncRNAs) play key roles in various epigenetic and post-transcriptional events in the cell, thereby significantly influencing cellular processes including gene expression, development and diseases such as cancer. Nuclear receptors (NRs) are a family of ligand-regulated transcription factors that typically regulate transcription of genes involved in a broad spectrum of cellular processes, immune responses and in many diseases including cancer. Owing to their many overlapping roles as modulators of gene expression, the paths traversed by lncRNA and NR-mediated signaling often cross each other; these lncRNA-NR cross-talks are being increasingly recognized as important players in many cellular processes and diseases such as cancer. Here, we review the individual roles of lncRNAs and NRs, especially growth factor modulated receptors such as androgen receptors (ARs), in various types of cancers and how the cross-talks between lncRNAs and NRs are involved in cancer progression and metastasis. We discuss the challenges involved in characterizing lncRNA-NR associations and how to overcome them. Furthering our understanding of the mechanisms of lncRNA-NR associations is crucial to realizing their potential as prognostic features, diagnostic biomarkers and therapeutic targets in cancer biology.

Nuclear Receptors and the Hidden Language of the Metabolome

Nuclear hormone receptors (NHRs) are a family of ligand-regulated transcription factors that control key aspects of development and physiology. The regulation of NHRs by ligands derived from metabolism or diet makes them excellent pharmacological targets, and the mechanistic understanding of how NHRs interact with their ligands to regulate downstream gene networks, along with the identification of ligands for orphan NHRs, could enable innovative approaches for cellular engineering, disease modeling and regenerative medicine. We review recent discoveries in the identification of physiologic ligands for NHRs. We propose new models of ligand-receptor co-evolution, the emergence of hormonal function and models of regulation of NHR specificity and activity via one-ligand and two-ligand models as well as feedback loops. Lastly, we discuss limitations on the processes for the identification of physiologic NHR ligands and emerging new methodologies that could be used to identify the natural ligands for the remaining 17 orphan NHRs in the human genome.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

Mitotic bookmarking redundancy by nuclear receptors in pluripotent cells

Mitotic bookmarking transcription factors (TFs) are thought to mediate rapid and accurate reactivation after mitotic gene silencing. However, the loss of individual bookmarking TFs often leads to the deregulation of only a small proportion of their mitotic targets, raising doubts on the biological significance and importance of their bookmarking function. Here we used targeted proteomics of the mitotic bookmarking TF ESRRB, an orphan nuclear receptor, to discover a large redundancy in mitotic binding among members of the protein super-family of nuclear receptors. Focusing on the nuclear receptor NR5A2, which together with ESRRB is essential in maintaining pluripotency in mouse embryonic stem cells, we demonstrate conjoint bookmarking activity of both factors on promoters and enhancers of a large fraction of active genes, particularly those most efficiently reactivated in G1. Upon fast and simultaneous degradation of both factors during mitotic exit, hundreds of mitotic targets of ESRRB/NR5A2, including key players of the pluripotency network, display attenuated transcriptional reactivation. We propose that redundancy in mitotic bookmarking TFs, especially nuclear receptors, confers robustness to the reestablishment of gene regulatory networks after mitosis.

PAPerFly: Partial Assembly-based Peak Finder for ab initio binding site reconstruction

BACKGROUND

The specific recognition of a DNA locus by a given transcription factor is a widely studied issue. It is generally agreed that the recognition can be influenced not only by the binding motif but by the larger context of the binding site. In this work, we present a novel heuristic algorithm that can reconstruct the unique binding sites captured in a sequencing experiment without using the reference genome.

RESULTS

We present PAPerFly, the Partial Assembly-based Peak Finder, a tool for the binding site and binding context reconstruction from the sequencing data without any prior knowledge. This tool operates without the need to know the reference genome of the respective organism. We employ algorithmic approaches that are used during genome assembly. The proposed algorithm constructs a de Bruijn graph from the sequencing data. Based on this graph, sequences and their enrichment are reconstructed using a novel heuristic algorithm. The reconstructed sequences are aligned and the peaks in the sequence enrichment are identified. Our approach was tested by processing several ChIP-seq experiments available in the ENCODE database and comparing the results of Paperfly and standard methods.

CONCLUSIONS

We show that PAPerFly, an algorithm tailored for experiment analysis without the reference genome, yields better results than an aggregation of ChIP-seq agnostic tools. Our tool is freely available at https://github.com/Caeph/paperfly/ or on Zenodo ( https://doi.org/10.5281/zenodo.7116424 ).

Retinoid X Receptor agonists as selective modulators of the immune system for the treatment of cancer

Upon heterodimerizing with other nuclear receptors, retinoid X receptors (RXR) act as ligand-dependent transcription factors, regulating transcription of critical signaling pathways that impact numerous hallmarks of cancer. By controlling both inflammation and immune responses, ligands that activate RXR can modulate the tumor microenvironment. Several small molecule agonists of these essential receptors have been synthesized. Historically, RXR agonists were tested for inhibition of growth in cancer cells, but more recent drug discovery programs screen new molecules for inhibition of inflammation or activation of immune cells. Bexarotene is the first successful example of an effective therapeutic that molecularly targets RXR; this drug was approved to treat cutaneous T cell lymphoma and is still used as a standard of care treatment for this disease. No additional RXR agonists have yet achieved FDA approval, but several promising novel compounds are being developed. In this review, we provide an overview of the multiple mechanisms by which RXR signaling regulates inflammation and tumor immunity. We also discuss the potential of RXR-dependent immune cell modulation for the treatment or prevention of cancer and concomitant challenges and opportunities.

NR5A2 connects zygotic genome activation to the first lineage segregation in totipotent embryos

Zygotic genome activation (ZGA) marks the beginning of the embryonic program for a totipotent embryo, which gives rise to the inner cell mass (ICM) where pluripotent epiblast arises, and extraembryonic trophectoderm. However, how ZGA is connected to the first lineage segregation in mammalian embryos remains elusive. Here, we investigated the role of nuclear receptor (NR) transcription factors (TFs), whose motifs are highly enriched and accessible from the 2-cell (2C) to 8-cell (8C) stages in mouse embryos. We found that NR5A2, an NR TF strongly induced upon ZGA, was required for this connection. Upon Nr5a2 knockdown or knockout, embryos developed beyond 2C normally with the zygotic genome largely activated. However, 4-8C-specific gene activation was substantially impaired and Nr5a2-deficient embryos subsequently arrested at the morula stage. Genome-wide chromatin binding analysis showed that NR5A2-bound cis-regulatory elements in both 2C and 8C embryos are strongly enriched for B1 elements where its binding motif is embedded. NR5A2 was not required for the global opening of its binding sites in 2C embryos but was essential to the opening of its 8C-specific binding sites. These 8C-specific, but not 2C-specific, binding sites are enriched near genes involved in blastocyst and stem cell regulation, and are often bound by master pluripotency TFs in blastocysts and embryonic stem cells (ESCs). Importantly, NR5A2 regulated key pluripotency genes Nanog and Pou5f1/Oct4, and primitive endoderm regulatory genes including Gata6 among many early ICM genes, as well as key trophectoderm regulatory genes including Tead4 and Gata3 at the 8C stage. By contrast, master pluripotency TFs NANOG, SOX2, and OCT4 targeted both early and late ICM genes in mouse ESCs. Taken together, these data identify NR5A2 as a key regulator in totipotent embryos that bridges ZGA to the first lineage segregation during mouse early development.

High-throughput data and modeling reveal insights into the mechanisms of cooperative DNA-binding by transcription factor proteins

Cooperative DNA-binding by transcription factor (TF) proteins is critical for eukaryotic gene regulation. In the human genome, many regulatory regions contain TF-binding sites in close proximity to each other, which can facilitate cooperative interactions. However, binding site proximity does not necessarily imply cooperative binding, as TFs can also bind independently to each of their neighboring target sites. Currently, the rules that drive cooperative TF binding are not well understood. In addition, it is oftentimes difficult to infer direct TF-TF cooperativity from existing DNA-binding data. Here, we show that in vitro binding assays using DNA libraries of a few thousand genomic sequences with putative cooperative TF-binding events can be used to develop accurate models of cooperativity and to gain insights into cooperative binding mechanisms. Using factors ETS1 and RUNX1 as our case study, we show that the distance and orientation between ETS1 sites are critical determinants of cooperative ETS1-ETS1 binding, while cooperative ETS1-RUNX1 interactions show more flexibility in distance and orientation and can be accurately predicted based on the affinity and sequence/shape features of the binding sites. The approach described here, combining custom experimental design with machine-learning modeling, can be easily applied to study the cooperative DNA-binding patterns of any TFs.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

Overview of Ursolic Acid Potential for the Treatment of Metabolic Disorders, Autoimmune Diseases, and Cancers via Nuclear Receptor Pathways

Nuclear receptors (NRs) form a family of druggable transcription factors that are regulated by ligand binding to orchestrate multifaceted physiological functions, including reproduction, immunity, metabolism, and growth. NRs represent attractive and valid targets for the management and treatment of a vast array of ailments. Pentacyclic triterpenes (PTs) are ubiquitously distributed natural products in medicinal and aromatic plants, of which ursolic acid (UA) is an extensively studied member, due to its diverse bio-pertinent activities against different cancers, inflammation, aging, obesity, diabetes, dyslipidemia, and liver injury. In fact, PTs share a common lipophilic structure that resembles NRs' endogenous ligands. Herein, we present a review of the literature on UA's effect on NRs, showcasing the resulting health benefits and potential therapeutic outcomes. De facto, UA exhibited numerous pharmacodynamic effects on PPAR, LXR, FXR, and PXR, resulting in remarkable anti-inflammatory, anti-hyperlipidemic, and hepatoprotective properties, by lowering lipid accumulation in hepatocytes and mitigating non-alcoholic steatohepatitis (NASH) and its subsequent liver fibrosis. Furthermore, UA reversed valproate and rifampicin-induced hepatic lipid accumulation. Additionally, UA showed great promise for the treatment of autoimmune inflammatory diseases such as multiple sclerosis and autoimmune arthritis by antagonizing RORγ. UA exhibited antiproliferative effects against skin, prostate, and breast cancers, partially via PPARα and RORγ pathways. Herein, for the first time, we explore and provide insights into UA bioactivity with respect to NR modulation.

Harnessing Oleanolic Acid and Its Derivatives as Modulators of Metabolic Nuclear Receptors

Nuclear receptors (NRs) constitute a superfamily of ligand-activated transcription factors with a paramount role in ubiquitous physiological functions such as metabolism, growth, and reproduction. Owing to their physiological role and druggability, NRs are deemed attractive and valid targets for medicinal chemists. Pentacyclic triterpenes (PTs) represent one of the most important phytochemical classes present in higher plants, where oleanolic acid (OA) is the most studied PTs representative owing to its multitude of biological activities against cancer, inflammation, diabetes, and liver injury. PTs possess a lipophilic skeleton that imitates the NRs endogenous ligands. Herein, we report a literature overview on the modulation of metabolic NRs by OA and its semi-synthetic derivatives, highlighting their health benefits and potential therapeutic applications. Indeed, OA exhibited varying pharmacological effects on FXR, PPAR, LXR, RXR, PXR, and ROR in a tissue-specific manner. Owing to these NRs modulation, OA showed prominent hepatoprotective properties comparable to ursodeoxycholic acid (UDCA) in a bile duct ligation mice model and antiatherosclerosis effect as simvastatin in a model of New Zealand white (NZW) rabbits. It also demonstrated a great promise in alleviating non-alcoholic steatohepatitis (NASH) and liver fibrosis, attenuated alpha-naphthol isothiocyanate (ANIT)-induced cholestatic liver injury, and controlled blood glucose levels, making it a key player in the therapy of metabolic diseases. We also compiled OA semi-synthetic derivatives and explored their synthetic pathways and pharmacological effects on NRs, showcasing their structure-activity relationship (SAR). To the best of our knowledge, this is the first review article to highlight OA activity in terms of NRs modulation.

Natural compounds targeting nuclear receptors for effective cancer therapy

Human nuclear receptors (NRs) are a family of forty-eight transcription factors that modulate gene expression both spatially and temporally. Numerous biochemical, physiological, and pathological processes including cell survival, proliferation, differentiation, metabolism, immune modulation, development, reproduction, and aging are extensively orchestrated by different NRs. The involvement of dysregulated NRs and NR-mediated signaling pathways in driving cancer cell hallmarks has been thoroughly investigated. Targeting NRs has been one of the major focuses of drug development strategies for cancer interventions. Interestingly, rapid progress in molecular biology and drug screening reveals that the naturally occurring compounds are promising modern oncology drugs which are free of potentially inevitable repercussions that are associated with synthetic compounds. Therefore, the purpose of this review is to draw our attention to the potential therapeutic effects of various classes of natural compounds that target NRs such as phytochemicals, dietary components, venom constituents, royal jelly-derived compounds, and microbial derivatives in the establishment of novel and safe medications for cancer treatment. This review also emphasizes molecular mechanisms and signaling pathways that are leveraged to promote the anti-cancer effects of these natural compounds. We have also critically reviewed and assessed the advantages and limitations of current preclinical and clinical studies on this subject for cancer prophylaxis. This might subsequently pave the way for new paradigms in the discovery of drugs that target specific cancer types.

The Endogenous Retinoic Acid Receptor Pathway Is Exploited by Mycobacterium tuberculosis during Infection, Both In Vitro and In Vivo

Retinoic acid (RA) is a fundamental vitamin A metabolite involved in regulating immune responses through the nuclear RA receptor (RAR) and retinoid X receptor. While performing experiments using THP-1 cells as a model for Mycobacterium tuberculosis infection, we observed that serum-supplemented cultures displayed high levels of baseline RAR activation in the presence of live, but not heat-killed, bacteria, suggesting that M. tuberculosis robustly induces the endogenous RAR pathway. Using in vitro and in vivo models, we have further explored the role of endogenous RAR activity in M. tuberculosis infection through pharmacological inhibition of RARs. We found that M. tuberculosis induces classical RA response element genes such as CD38 and DHRS3 in both THP-1 cells and human primary CD14+ monocytes via a RAR-dependent pathway. M. tuberculosis-stimulated RAR activation was observed with conditioned media and required nonproteinaceous factor(s) present in FBS. Importantly, RAR blockade by (4-[(E)-2-[5,5-dimethyl-8-(2-phenylethynyl)-6H-naphthalen-2-yl]ethenyl]benzoic acid), a specific pan-RAR inverse agonist, in a low-dose murine model of tuberculosis significantly reduced SIGLEC-F+CD64+CD11c+high alveolar macrophages in the lungs, which correlated with 2× reduction in tissue mycobacterial burden. These results suggest that the endogenous RAR activation axis contributes to M. tuberculosis infection both in vitro and in vivo and reveal an opportunity for further investigation of new antituberculosis therapies.

Vitamin A and retinoid signaling in the kidneys

Vitamin A (VA, retinol) and its metabolites (commonly called retinoids) are required for the proper development of the kidney during embryogenesis, but retinoids also play key roles in the function and repair of the kidney in adults. Kidneys filter 180-200 liters of blood per day and each kidney contains approximately 1 million nephrons, which are often referred to as the 'functional units' of the kidney. Each nephron consists of a glomerulus and a series of tubules (proximal tubule, loop of Henle, distal tubule, and collecting duct) surrounded by a network of capillaries. VA is stored in the liver and converted to active metabolites, most notably retinoic acid (RA), which acts as an agonist for the retinoic acid receptors ((RARs α, β, and γ) to regulate gene transcription. In this review we discuss some of the actions of retinoids in the kidney after injury. For example, in an ischemia-reperfusion model in mice, injury-associated loss of proximal tubule (PT) differentiation markers occurs, followed by re-expression of these differentiation markers during PT repair. Notably, healthy proximal tubules express ALDH1a2, the enzyme that metabolizes retinaldehyde to RA, but transiently lose ALDH1a2 expression after injury, while nearby myofibroblasts transiently acquire RA-producing capabilities after injury. These results indicate that RA is important for renal tubular injury repair and that compensatory mechanisms exist for the generation of endogenous RA by other cell types upon proximal tubule injury. ALDH1a2 levels also increase in podocytes, epithelial cells of the glomeruli, after injury, and RA promotes podocyte differentiation. We also review the ability of exogenous, pharmacological doses of RA and receptor selective retinoids to treat numerous kidney diseases, including kidney cancer and diabetic kidney disease, and the emerging genetic evidence for the importance of retinoids and their receptors in maintaining or restoring kidney function after injury. In general, RA has a protective effect on the kidney after various types of injuries (eg. ischemia, cytotoxic actions of chemicals, hyperglycemia related to diabetes). As more research into the actions of each of the three RARs in the kidney is carried out, a greater understanding of the actions of vitamin A is likely to lead to new insights into the pathology of kidney disorders and the development of new therapies for kidney diseases.

Hidden modes of DNA binding by human nuclear receptors

Human nuclear receptors (NRs) are a superfamily of ligand-responsive transcription factors that have central roles in cellular function. Their malfunction is linked to numerous diseases, and the ability to modulate their activity with synthetic ligands has yielded 16% of all FDA-approved drugs. NRs regulate distinct gene networks, however they often function from genomic sites that lack known binding motifs. Here, to annotate genomic binding sites of known and unexamined NRs more accurately, we use high-throughput SELEX to comprehensively map DNA binding site preferences of all full-length human NRs, in complex with their ligands. Furthermore, to identify non-obvious binding sites buried in DNA-protein interactomes, we develop MinSeq Find, a search algorithm based on the MinTerm concept from electrical engineering and digital systems design. The resulting MinTerm sequence set (MinSeqs) reveal a constellation of binding sites that more effectively annotate NR-binding profiles in cells. MinSeqs also unmask binding sites created or disrupted by 52,106 single-nucleotide polymorphisms associated with human diseases. By implicating druggable NRs as hidden drivers of multiple human diseases, our results not only reveal new biological roles of NRs, but they also provide a resource for drug-repurposing and precision medicine.

Prediction of cooperative homeodomain DNA binding sites from high-throughput-SELEX data

Homeodomain proteins constitute one of the largest families of metazoan transcription factors. Genetic studies have demonstrated that homeodomain proteins regulate many developmental processes. Yet, biochemical data reveal that most bind highly similar DNA sequences. Defining how homeodomain proteins achieve DNA binding specificity has therefore been a long-standing goal. Here, we developed a novel computational approach to predict cooperative dimeric binding of homeodomain proteins using High-Throughput (HT) SELEX data. Importantly, we found that 15 of 88 homeodomain factors form cooperative homodimer complexes on DNA sites with precise spacing requirements. Approximately one third of the paired-like homeodomain proteins cooperatively bind palindromic sequences spaced 3 bp apart, whereas other homeodomain proteins cooperatively bind sites with distinct orientation and spacing requirements. Combining structural models of a paired-like factor with our cooperativity predictions identified key amino acid differences that help differentiate between cooperative and non-cooperative factors. Finally, we confirmed predicted cooperative dimer sites in vivo using available genomic data for a subset of factors. These findings demonstrate how HT-SELEX data can be computationally mined to predict cooperativity. In addition, the binding site spacing requirements of select homeodomain proteins provide a mechanism by which seemingly similar AT-rich DNA sequences can preferentially recruit specific homeodomain factors.

Zinc cluster transcription factors frequently activate target genes using a non-canonical half-site binding mode

Gene expression changes are orchestrated by transcription factors (TFs), which bind to DNA to regulate gene expression. It remains surprisingly difficult to predict basic features of the transcriptional process, including in vivo TF occupancy. Existing thermodynamic models of TF function are often not concordant with experimental measurements, suggesting undiscovered biology. Here, we analyzed one of the most well-studied TFs, the yeast zinc cluster Gal4, constructed a Shea-Ackers thermodynamic model to describe its binding, and compared the results of this model to experimentally measured Gal4p binding in vivo. We found that at many promoters, the model predicted no Gal4p binding, yet substantial binding was observed. These outlier promoters lacked canonical binding motifs, and subsequent investigation revealed Gal4p binds unexpectedly to DNA sequences with high densities of its half site (CGG). We confirmed this novel mode of binding through multiple experimental and computational paradigms; we also found most other zinc cluster TFs we tested frequently utilize this binding mode, at 27% of their targets on average. Together, these results demonstrate a novel mode of binding where zinc clusters, the largest class of TFs in yeast, bind DNA sequences with high densities of half sites.

Structural basis of the farnesoid X receptor/retinoid X receptor heterodimer on inverted repeat DNA

Farnesoid X receptor (FXR) is a ligand-activated transcription factor known as bile acid receptor (BAR). FXR plays critical roles in various biological processes, including metabolism, immune inflammation, liver regeneration and liver carcinogenesis. FXR forms a heterodimer with the retinoid X receptor (RXR) and binds to diverse FXR response elements (FXREs) to exert its various biological functions. However, the mechanism by which the FXR/RXR heterodimer binds the DNA elements remains unclear. In this study, we aimed to use structural, biochemical and bioinformatics analyses to study the mechanism of FXR binding to the typical FXRE, such as the IR1 site, and the heterodimer interactions in the FXR-DBD/RXR-DBD complex. Further biochemical assays showed that RAR, THR and NR4A2 do not form heterodimers with RXR when bound to the IR1 sites, which indicates that IR1 may be a unique binding site for the FXR/RXR heterodimer. Our studies may provide a further understanding of the dimerization specificity of nuclear receptors.

Effects of all-trans and 9-cis retinoic acid on differentiating human neural stem cells in vitro

Cyanobacterial blooms are known sources of environmentally-occurring retinoid compounds, including all-trans and 9-cis retinoic acids (RAs). The developmental hazard for aquatic organisms has been described, while the implications for human health hazard assessment are not yet sufficiently characterized. Here, we employ a human neural stem cell model that can differentiate in vitro into a mixed culture of neurons and glia. Cells were exposed to non-cytotoxic 8-1000 nM all-trans or 9-cis RA for 9-18 days (DIV13 and DIV22, respectively). Impact on biomarkers was analyzed on gene expression (RT-qPCR) and protein level (western blot and proteomics) at both time points; network patterning (immunofluorescence) on DIV22. RA exposure significantly concentration-dependently increased gene expression of retinoic acid receptors and the metabolizing enzyme CYP26A1, confirming the chemical-specific response of the model. Expression of thyroid hormone signaling-related genes remained mostly unchanged. Markers of neural progenitors/stem cells (PAX6, SOX1, SOX2, NESTIN) were decreased with increasing RA concentrations, though a basal population remained. Neural markers (DCX, TUJ1, MAP2, NeuN, SYP) remained unchanged or were decreased at high concentrations (200-1000 nM). Conversely, (astro-)glial marker S100β was increased concentration-dependently on DIV22. Together, the biomarker analysis indicates an RA-dependent promotion of glial cell fates over neural differentiation, despite the increased abundance of neural protein biomarkers during differentiation. Interestingly, RA exposure induced substantial changes to the cell culture morphology: while low concentrations resulted in a network-like differentiation pattern, high concentrations (200-1000 nM RA) almost completely prevented such network patterning. After functional confirmation for implications in network function, such morphological features could present a proxy for network formation assessment, an apical key event in (neuro-)developmental Adverse Outcome Pathways. The described application of a human in vitro model for (developmental) neurotoxicity to emerging environmentally-relevant retinoids contributes to the evidence-base for the use of differentiating human in vitro models for human health hazard and risk assessment.

Parallelized multidimensional analytic framework applied to mammary epithelial cells uncovers regulatory principles in EMT

A proper understanding of disease etiology will require longitudinal systems-scale reconstruction of the multitiered architecture of eukaryotic signaling. Here we combine state-of-the-art data acquisition platforms and bioinformatics tools to devise PAMAF, a workflow that simultaneously examines twelve omics modalities, i.e., protein abundance from whole-cells, nucleus, exosomes, secretome and membrane; N-glycosylation, phosphorylation; metabolites; mRNA, miRNA; and, in parallel, single-cell transcriptomes. We apply PAMAF in an established in vitro model of TGFβ-induced epithelial to mesenchymal transition (EMT) to quantify >61,000 molecules from 12 omics and 10 timepoints over 12 days. Bioinformatics analysis of this EMT-ExMap resource allowed us to identify; -topological coupling between omics, -four distinct cell states during EMT, -omics-specific kinetic paths, -stage-specific multi-omics characteristics, -distinct regulatory classes of genes, -ligand-receptor mediated intercellular crosstalk by integrating scRNAseq and subcellular proteomics, and -combinatorial drug targets (e.g., Hedgehog signaling and CAMK-II) to inhibit EMT, which we validate using a 3D mammary duct-on-a-chip platform. Overall, this study provides a resource on TGFβ signaling and EMT.

A structural perspective of liver X receptors

Liver X receptors α and β are members of the nuclear receptor family, which comprise a flexible N-terminal domain, a DNA binding domain, a hinge linker, and a ligand binding domain. Liver X receptors are important regulators of cholesterol and lipid homeostasis by controlling the transcription of numerous genes. Key to their transcriptional role is synergetic interaction among the domains. DNA binding domain binds on DNA; ligand binding domain is a crucial switch to control the transcription activity through conformational change caused by ligand binding. The Liver X receptors form heterodimers with retinoid X receptor and then the liganded heterodimer may recruit other necessary transcription components to form an active transcription complex.

Interactions governing transcriptional activity of nuclear receptors

The key players in transcriptional regulation are transcription factors (TFs), proteins that bind specific DNA sequences. Several mechanisms exist to turn TFs ‘on’ and ‘off’, including ligand binding which induces conformational changes within TFs, subsequently influencing multiple inter- and intramolecular interactions to drive transcriptional responses. Nuclear receptors are a specific family of ligand-regulated TFs whose activity relies on interactions with DNA, coregulator proteins and other receptors. These multidomain proteins also undergo interdomain interactions on multiple levels, further modulating transcriptional outputs. Cooperation between these distinct interactions is critical for appropriate transcription and remains an intense area of investigation. In this review, we report and summarize recent findings that continue to advance our mechanistic understanding of how interactions between nuclear receptors and diverse partners influence transcription.

Hormonal Crosstalk Between Thyroid and Breast Cancer

Differentiated thyroid cancer and breast cancer account for a significant portion of endocrine-related malignancies and predominately affect women. As hormonally responsive tissues, the breast and thyroid share endocrine signaling. Breast cells are responsive to thyroid hormone signaling and are affected by altered thyroid hormone levels. Thyroid cells are responsive to sex hormones, particularly estrogen, and undergo protumorigenic processes upon estrogen stimulation. Thyroid and sex hormones also display significant transcriptional crosstalk that influences oncogenesis and treatment sensitivity. Obesity-related adipocyte alterations-adipocyte estrogen production, inflammation, feeding hormone dysregulation, and metabolic syndromes-promote hormonal alterations in breast and thyroid tissues. Environmental toxicants disrupt endocrine systems, including breast and thyroid homeostasis, and influence pathologic processes in both organs through hormone mimetic action. In this brief review, we discuss the hormonal connections between the breast and thyroid and perspectives on hormonal therapies for breast and thyroid cancer. Future research efforts should acknowledge and further explore the hormonal crosstalk of these tissues in an effort to further understand the prevalence of thyroid and breast cancer in women and to identify potential therapeutic options.

Plasmodium falciparum Drug Resistance Genes pfmdr1 and pfcrt In Vivo Co-Expression During Artemether-Lumefantrine Therapy

Background: Artemisinin-based combination therapies (ACTs) are the global mainstay treatment of uncomplicated Plasmodium falciparum infections. PfMDR1 and PfCRT are two transmembrane transporters, associated with sensitivity to several antimalarials, found in the parasite food vacuole. Herein, we explore if their relatedness extends to overlapping patterns of gene transcriptional activity before and during ACT administration.

Methods: In a clinical trial performed in Tanzania, we explored the pfmdr1 and pfcrt transcription levels from 48 patients with uncomplicated P. falciparum malaria infections who underwent treatment with artemether-lumefantrine (AL). Samples analyzed were collected before treatment initiation and during the first 24 h of treatment. The frequency of PfMDR1 N86Y and PfCRT K76T was determined through PCR-RFLP or direct amplicon sequencing. Gene expression was analyzed by real-time quantitative PCR.

Results: A wide range of pre-treatment expression levels was observed for both genes, approximately 10-fold for pfcrt and 50-fold for pfmdr1. In addition, a significant positive correlation demonstrates pfmdr1 and pfcrt co-expression. After AL treatment initiation, pfmdr1 and pfcrt maintained the positive co-expression correlation, with mild downregulation throughout the 24 h post-treatment. Additionally, a trend was observed for PfMDR1 N86 alleles and higher expression before treatment initiation.

Conclusion:pfmdr1 and pfcrt showed significant co-expression patterns in vivo, which were generally maintained during ACT treatment. This observation points to relevant related roles in the normal parasite physiology, which seem essential to be maintained when the parasite is exposed to drug stress. In addition, keeping the simultaneous expression of both transporters might be advantageous for responding to the drug action.

Teratogenicity of retinoids detected in surface waters in zebrafish embryos and its predictability by in vitro assays

Retinoids are newly detected compounds in aquatic ecosystems associated with cyanobacterial water blooms. Their potential health risks are only scarcely described despite numerous detections of all-trans retinoic acid (ATRA) and its derivatives in the environment. Besides the known teratogen ATRA there is only little or no information about their potency and namely their effects in vivo. We characterize ATRA and 8 other retinoids reported to occur in the environment for their bioactivity and teratogenicity using four in vitro reporter gene assays and zebrafish (Danio rerio) embryotoxicity assay. Our results document the ability of these compounds to interfere with retinoid signalling and cause teratogenicity at environmentally relevant levels with EC₅₀ values at nM (hundreds of ng/L) levels and teratogenic indexes ranging from 2.8 (9cis retinoic acid) to 15.8 (retinal). The relative potency of individual compounds for teratogenicity ranged from 0.059 (retinal) to 0.96 (5,6-epoxy ATRA) when compared to ATRA. An environmentally relevant mixture of retinoids was tested showing good predictability of teratogenicity from the in vitro activities and additive toxicity of the mixture. The high teratogenicity of the newly described compounds associated with cyanobacteria presents a concern for developmental stages due to high conservation of the retinoid signalling across vertebrates.

PPARγ phase separates with RXRα at PPREs to regulate target gene expression

Peroxisome proliferator-activated receptor (PPAR)-γ is a key transcription activator controlling adipogenesis and lipid metabolism. PPARγ binds PPAR response elements (PPREs) as the obligate heterodimer with retinoid X receptor (RXR) α, but exactly how PPARγ orchestrates the transcriptional response is unknown. This study demonstrates that PPARγ forms phase-separated droplets in vitro and solid-like nuclear condensates in cell, which is intriguingly mediated by its DNA binding domain characterized by the zinc finger motif. Furthermore, PPARγ forms nuclear condensates at PPREs sites through phase separation to compartmentalize its heterodimer partner RXRα to initiate PPARγ-specific transcriptional activation. Finally, using an optogenetic approach, the enforced formation of PPARγ/RXRα condensates leads to preferential enrichment at PPREs sites and significantly promotes the expression of PPARγ target genes. These results define a novel mechanism by which PPARγ engages the phase separation principles for efficient and specific transcriptional activation.

Gene repression through epigenetic modulation by PPARA enhances hepatocellular proliferation

Peroxisome proliferator-activated receptor α (PPARA) is a key mediator of lipid metabolism and inflammation. Activation of PPARA in rodents causes hepatocyte proliferation, but the underlying mechanism is poorly understood. This study focused on genes repressed by PPARA and analyzed the mechanism by which PPARA promotes hepatocyte proliferation in mice. Activation of PPARA by agonist treatment was autoregulated, and induced expression of the epigenetic regulator UHRF1 via activation of the newly described PPARA target gene E2f8, which, in turn, regulates Uhrf1. UHRF1 strongly repressed the expression of CDH1 via methylation of the Cdh1 promoter marked with H3K9me3. Repression of CDH1 by PPARA activation was reversed by PPARA deficiency or knockdown of E2F8 or UHRF1. Furthermore, a forced expression of CDH1 inhibited expression of the Wnt signaling target genes such as Myc after PPARA activation, and suppressed hepatocyte hyperproliferation. These results demonstrate that the PPARA-E2F8-UHRF1-CDH1 axis causes epigenetic regulation of hepatocyte proliferation.

•

PPARA activation induces the UHRF1 expression via novel PPARA target gene E2f8

•

Induction of UHRF1 by PPARA activation represses Cdh1 gene marked with H3K9me3

•

CDH1 suppresses hepatocyte proliferation after PPARA activation

•

Autoinduction of PPARA by agonist enhances cell proliferation via E2F8-UHRF1-CDH1

Mechanisms of Feedback Regulation of Vitamin A Metabolism

Vitamin A is an essential nutrient required throughout life. Through its various metabolites, vitamin A sustains fetal development, immunity, vision, and the maintenance, regulation, and repair of adult tissues. Abnormal tissue levels of the vitamin A metabolite, retinoic acid, can result in detrimental effects which can include congenital defects, immune deficiencies, proliferative defects, and toxicity. For this reason, intricate feedback mechanisms have evolved to allow tissues to generate appropriate levels of active retinoid metabolites despite variations in the level and format, or in the absorption and conversion efficiency of dietary vitamin A precursors. Here, we review basic mechanisms that govern vitamin A signaling and metabolism, and we focus on retinoic acid-controlled feedback mechanisms that contribute to vitamin A homeostasis. Several approaches to investigate mechanistic details of the vitamin A homeostatic regulation using genomic, gene editing, and chromatin capture technologies are also discussed.

CASCADE: high-throughput characterization of regulatory complex binding altered by non-coding variants

Non-coding DNA variants (NCVs) impact gene expression by altering binding sites for regulatory complexes. New high-throughput methods are needed to characterize the impact of NCVs on regulatory complexes. We developed CASCADE (Customizable Approach to Survey Complex Assembly at DNA Elements), an array-based high-throughput method to profile cofactor (COF) recruitment. CASCADE identifies DNA-bound transcription factor-cofactor (TF-COF) complexes in nuclear extracts and quantifies the impact of NCVs on their binding. We demonstrate CASCADE sensitivity in characterizing condition-specific recruitment of COFs p300 and RBBP5 (MLL subunit) to the CXCL10 promoter in lipopolysaccharide (LPS)-stimulated human macrophages and quantify the impact of all possible NCVs. To demonstrate applicability to NCV screens, we profile TF-COF binding to ~1,700 single-nucleotide polymorphism quantitative trait loci (SNP-QTLs) in human macrophages and identify perturbed ETS domain-containing complexes. CASCADE will facilitate high-throughput testing of molecular mechanisms of NCVs for diverse biological applications.

Bray et al. develop CASCADE, a method to profile transcription factor (TF)-cofactor (COF) complexes binding to DNA. They demonstrate the approach by profiling complex binding across the CXCL10 cytokine promoter and to ~1,700 single-nucleotide polymorphisms (SNPs). They anticipate that CASCADE can be applied to diverse biological systems to examine regulatory complex binding to DNA.

Tbx5 drives Aldh1a2 expression to regulate a RA-Hedgehog-Wnt gene regulatory network coordinating cardiopulmonary development

The gene regulatory networks that coordinate the development of the cardiac and pulmonary systems are essential for terrestrial life but poorly understood. The T-box transcription factor Tbx5 is critical for both pulmonary specification and heart development, but how these activities are mechanistically integrated remains unclear. Here using Xenopus and mouse embryos, we establish molecular links between Tbx5 and retinoic acid (RA) signaling in the mesoderm and between RA signaling and sonic hedgehog expression in the endoderm to unveil a conserved RA-Hedgehog-Wnt signaling cascade coordinating cardiopulmonary (CP) development. We demonstrate that Tbx5 directly maintains expression of aldh1a2, the RA-synthesizing enzyme, in the foregut lateral plate mesoderm via an evolutionarily conserved intronic enhancer. Tbx5 promotes posterior second heart field identity in a positive feedback loop with RA, antagonizing a Fgf8-Cyp regulatory module to restrict FGF activity to the anterior. We find that Tbx5/Aldh1a2-dependent RA signaling directly activates shh transcription in the adjacent foregut endoderm through a conserved MACS1 enhancer. Hedgehog signaling coordinates with Tbx5 in the mesoderm to activate expression of wnt2/2b, which induces pulmonary fate in the foregut endoderm. These results provide mechanistic insight into the interrelationship between heart and lung development informing CP evolution and birth defects.

PPARγ attenuates hepatic inflammation and oxidative stress of non‑alcoholic steatohepatitis via modulating the miR‑21‑5p/SFRP5 pathway

Inflammation and oxidative stress are key steps in the progression of non‑alcoholic steatohepatitis (NASH). Intervention in these two processes will therefore benefit NASH treatment. Peroxisome proliferator‑activated receptor γ (PPARγ), as a multiple functional transcription factor, has been reported to be involved in the prevention of NASH progression. However, the mechanism by which PPARγ prevents NASH remains to be elucidated. The present study demonstrated that the level of PPARγ was inversely correlated with that of microRNA (miRNA/miRs)‑21‑5p in both mice and humans with NASH. Activation of PPARγ inhibited lipid droplet accumulation, hepatic inflammation and oxidative stress by downregulating miR‑21‑5p in an in vitro model. Luciferase reporter and chromatin immunoprecipitation assays demonstrated that PPARγ suppressed transcriptional activity of miR‑21‑5p and bound to miR‑21‑5p promoter region. Furthermore, PPARγ downregulated miR‑21‑5p while miR‑21‑5p upregulated secreted frizzled‑related protein 5 (SFRP5) by targeting the 3'‑UTR of its mRNA. In vivo experiments revealed that PPARγ repressed inflammation and oxidative stress and miR‑21‑5p expression while increased SFRP5 level in a NASH mouse model. In summary, PPARγ attenuates inflammation and oxidative stress in NASH by modulating the miR‑21‑5p/SFRP5 pathway, thus holding promise of a new target for NASH treatment.

Ligands and DNA in the allosteric control of retinoid receptors function

Retinoids are a family of compounds that include both vitamin A (all-trans retinol) and its naturally occurring metabolites such as retinoic acids (e.g. all-trans retinoic acid) as well as synthetic analogs. They are critically involved in the regulation of a wide variety of essential biological processes, such as embryogenesis and organogenesis, apoptosis, reproduction, vision, and the growth and differentiation of normal and neoplastic cells in vertebrates. The ability of these small molecules to control the expression of several hundred genes through binding to nuclear ligand-dependent transcription factors accounts for most of their functions. Three retinoic acid receptor (RARα,β,γ) and three retinoid X receptor (RXRα,β,γ) subtypes form a variety of RXR-RAR heterodimers that have been shown to mediate the pleiotropic effects of retinoids through the recruitment of high-molecular weight co-regulatory complexes to response-element DNA sequences found in the promoter region of their target genes. Hence, heterodimeric retinoid receptors are multidomain entities that respond to various incoming signals, such as ligand and DNA binding, by allosteric structural alterations which are the basis of further signal propagation. Here, we provide an overview of the current state of knowledge with regard to the structural mechanisms by which retinoids and DNA response elements act as allosteric effectors that may combine to finely tune RXR-RAR heterodimers activity.

Competition for DNA binding between paralogous transcription factors determines their genomic occupancy and regulatory functions

Most eukaryotic transcription factors (TFs) are part of large protein families, with members of the same family (i.e., paralogous TFs) recognizing similar DNA-binding motifs but performing different regulatory functions. Many TF paralogs are coexpressed in the cell and thus can compete for target sites across the genome. However, this competition is rarely taken into account when studying the in vivo binding patterns of eukaryotic TFs. Here, we show that direct competition for DNA binding between TF paralogs is a major determinant of their genomic binding patterns. Using yeast proteins Cbf1 and Pho4 as our model system, we designed a high-throughput quantitative assay to capture the genomic binding profiles of competing TFs in a cell-free system. Our data show that Cbf1 and Pho4 greatly influence each other's occupancy by competing for their common putative genomic binding sites. The competition is different at different genomic sites, as dictated by the TFs' expression levels and their divergence in DNA-binding specificity and affinity. Analyses of ChIP-seq data show that the biophysical rules that dictate the competitive TF binding patterns in vitro are also followed in vivo, in the complex cellular environment. Furthermore, the Cbf1-Pho4 competition for genomic sites, as characterized in vitro using our new assay, plays a critical role in the specific activation of their target genes in the cell. Overall, our study highlights the importance of direct TF-TF competition for genomic binding and gene regulation by TF paralogs, and proposes an approach for studying this competition in a quantitative and high-throughput manner.

Retinoic acid signaling is critical during the totipotency window in early mammalian development

Totipotent cells hold enormous potential for regenerative medicine. Thus, the development of cellular models recapitulating totipotent-like features is of paramount importance. Cells resembling the totipotent cells of early embryos arise spontaneously in mouse embryonic stem (ES) cell cultures. Such '2-cell-like-cells' (2CLCs) recapitulate 2-cell-stage features and display expanded cell potential. Here, we used 2CLCs to perform a small-molecule screen to identify new pathways regulating the 2-cell-stage program. We identified retinoids as robust inducers of 2CLCs and the retinoic acid (RA)-signaling pathway as a key component of the regulatory circuitry of totipotent cells in embryos. Using single-cell RNA-seq, we reveal the transcriptional dynamics of 2CLC reprogramming and show that ES cells undergo distinct cellular trajectories in response to RA. Importantly, endogenous RA activity in early embryos is essential for zygotic genome activation and developmental progression. Overall, our data shed light on the gene regulatory networks controlling cellular plasticity and the totipotency program.

Cellular Responses to the Efferocytosis of Apoptotic Cells

The rapid and efficient phagocytic clearance of apoptotic cells, termed efferocytosis, is a critical mechanism in the maintenance of tissue homeostasis. Removal of apoptotic cells through efferocytosis prevents secondary necrosis and the resultant inflammation caused by the release of intracellular contents. The importance of efferocytosis in homeostasis is underscored by the large number of inflammatory and autoimmune disorders, including atherosclerosis and systemic lupus erythematosus, that are characterized by defective apoptotic cell clearance. Although mechanistically similar to the phagocytic clearance of pathogens, efferocytosis differs from phagocytosis in that it is immunologically silent and induces a tissue repair response. Efferocytes face unique challenges resulting from the internalization of apoptotic cells, including degradation of the apoptotic cell, dealing with the extra metabolic load imposed by the processing of apoptotic cell contents, and the coordination of an anti-inflammatory, pro-tissue repair response. This review will discuss recent advances in our understanding of the cellular response to apoptotic cell uptake, including trafficking of apoptotic cell cargo and antigen presentation, signaling and transcriptional events initiated by efferocytosis, the coordination of an anti-inflammatory response and tissue repair, unique cellular metabolic responses and the role of efferocytosis in host defense. A better understanding of how efferocytic cells respond to apoptotic cell uptake will be critical in unraveling the complex connections between apoptotic cell removal and inflammation resolution and maintenance of tissue homeostasis.

Danthron ameliorates obesity and MAFLD through activating the interplay between PPARα/RXRα heterodimer and adiponectin receptor 2

Obesity and associated metabolic associated fatty liver diseases (MAFLD) are strongly associated with dysfunction of glucose and lipid metabolism. AMPKα and PPARα are key regulators in the lipid and glucose homeostasis, indicating that novel agents to activate them are promising therapeutic approaches for metabolic syndrome. Noticeably, as a natural anthraquinone derivative extracted from rhubarb, danthron can activate AMPKα in vitro. However, the protective effect of danthron on obesity and associated MAFLD in vivo, as well as the underlying mechanism remains unknown. In this study, obesity and associated MAFLD was induced in C57BL/6J mice by high fat diet (HFD), which were subjected to evaluations on the parameters of systematic metabolism. Simultaneously, the molecular mechanism of danthron on lipid metabolism was investigated in 3T3-L1-derived adipocytes and HepG2 cells in vitro. In vivo, danthron significantly attenuated the obesity and MAFLD by enhancing hepatic fatty acid oxidation, decreasing lipid synthesis, and promoting mitochondrial homeostasis. Mechanistically, danthron significantly promoted combination of RXRα and PPARα, enhanced the binding of RXRα/PPARα heterodimer to the promoter of adiponectin receptor 2 (AdipoR2), by which activating the AMPKα and PPARα pathway. Moreover, PPARα and AdipoR2 can interplay in a loop style. Collectively, this study demonstrates that danthron can substantially ameliorate obesity and associated hepatic steatosis via AdipoR2-mediated dual PPARα/AMPKα activation, which suggests that danthron might be a novel therapeutic approach for inhibition of obesity and hepatic steatosis.

Interaction of BDE-47 with nuclear receptors (NRs) based on the cytotoxicity: In vitro investigation and molecular interaction

Polybrominated diphenyl ethers (PBDEs) are endocrine-disrupting chemicals that possess neuroendocrine and reproductive toxicity to humans and disturb thyroid hormone homeostasis, neurobehavior, and development. The most predominant congener of PBDEs in humans and other organisms is 2,2',4,4'-tetrabromodiphenyl ether (BDE-47); however, the molecular mechanisms underlying its cytotoxicity remain largely unknown. Here, we evaluated the toxic effect and underlying mechanism of nuclear receptors (NRs) induced by BDE-47 in SK-N-SH human neuroblastoma cells. The CCK-8 cell viability assay showed that the proliferation of human SK-N-SH cells exposed to BDE-47 was significantly inhibited in time- and dose-dependent manners, and flow cytometry showed that cell cycle was arrested at the S phase after BDE-47 exposure. Moreover, compared with the control group, the expression of retinoic acid receptor alpha (RXRα), pregnane X receptor (PXR), thyroid hormone receptors (TRs), and peroxisome proliferator-activated receptors (PPARs) at the mRNA and protein levels was significantly increased, as determined by quantitative PCR and western blot analysis, demonstrating that BDE-47 activated the NRs in vitro. Moreover, BDE-47 could bind to all four NRs in the affinity order of PPARγ > PXR > TRβ > RXRα under molecular dynamics. Because RXR is the promiscuous dimerization partner for a large number of NRs, ZDock was used to calculate its interaction with other three NRs. Taking the number of hydrogen bonds and ZDock scores into account, the rank of docking ability between RXRα and the NRs was PXR > TRβ > PPARγ. Further analysis of the interaction between BDE-47 and dimerized-NRs, the affinity order was RXRα > TRβ > PXR > PPARγ via Glide. The results of this study demonstrated that BDE-47 interfered the cross-talk among NRs, especially the promiscuous RXRα, which might be critical for the harmonized re-adjustment of cytotoxicity and biological regulation. Our findings provide a better understanding of the mechanisms underlying toxic effects and intermolecular interaction induced by BDE-47.

High-Throughput Analysis of the Cell and DNA Site-Specific Binding of Native NF-κB Dimers Using Nuclear Extract Protein-Binding Microarrays (NextPBMs)

Nuclear factor-kappa B (NF-κB) transcription factors coordinate gene expression in response to a broad array of cellular signals. In vertebrates, there are five NF-κB proteins (c-Rel, RelA/p65, RelB, p50, and p52) that can form various dimeric combinations exhibiting both common and dimer-specific DNA-binding specificity. In this chapter, we describe the use of the nuclear extract protein-binding microarray (nextPBM), a high-throughput method to characterize the DNA binding of transcription factors present in cell nuclear extracts. NextPBMs allow for sensitive analysis of the DNA binding of NF-κB dimers and their interactions with cell-specific cofactors.

G ATA2 mediates the negative regulation of the prepro-thyrotropin-releasing hormone gene by liganded T3 receptor β2 in the rat hypothalamic paraventricular nucleus

Thyroid hormone (T3) inhibits thyrotropin-releasing hormone (TRH) synthesis in the hypothalamic paraventricular nucleus (PVN). Although the T3 receptor (TR) β2 is known to mediate the negative regulation of the prepro-TRH gene, its molecular mechanism remains unknown. Our previous studies on the T3-dependent negative regulation of the thyrotropin β subunit (TSHβ) gene suggest that there is a tethering mechanism, whereby liganded TRβ2 interferes with the function of the transcription factor, GATA2, a critical activator of the TSHβ gene. Interestingly, the transcription factors Sim1 and Arnt2, the determinants of PVN differentiation in the hypothalamus, are reported to induce expression of TRβ2 and GATA2 in cultured neuronal cells. Here, we confirmed the expression of the GATA2 protein in the TRH neuron of the rat PVN using immunohistochemistry with an anti-GATA2 antibody. According to an experimental study from transgenic mice, a region of the rat prepro-TRH promoter from nt. -547 to nt. +84 was able to mediate its expression in the PVN. We constructed a chloramphenicol acetyltransferase (CAT) reporter gene containing this promoter sequence (rTRH(547)-CAT) and showed that GATA2 activated the promoter in monkey kidney-derived CV1 cells. Deletion and mutation analyses identified a functional GATA-responsive element (GATA-RE) between nt. -357 and nt. -352. When TRβ2 was co-expressed, T3 reduced GATA2-dependent promoter activity to approximately 30%. Unexpectedly, T3-dependent negative regulation was maintained after mutation of the reported negative T3-responsive element, site 4. T3 also inhibited the GATA2-dependent transcription enhanced by cAMP agonist, 8-bromo-cAMP. A rat thyroid medullary carcinoma cell line, CA77, is known to express the preproTRH mRNA. Using a chromatin immunoprecipitation assay with this cell line where GATA2 expression plasmid was transfected, we observed the recognition of the GATA-RE by GATA2. We also confirmed GATA2 binding using gel shift assay with the probe for the GATA-RE. In CA77 cells, the activity of rTRH(547)-CAT was potentiated by overexpression of GATA2, and it was inhibited in a T3-dependent manner. These results suggest that GATA2 transactivates the rat prepro-TRH gene and that liganded TRβ2 interferes with this activation via a tethering mechanism as in the case of the TSHβ gene.

Structural basis for DNA recognition and allosteric control of the retinoic acid receptors RAR-RXR

Retinoic acid receptors (RARs) as a functional heterodimer with retinoid X receptors (RXRs), bind a diverse series of RA-response elements (RAREs) in regulated genes. Among them, the non-canonical DR0 elements are bound by RXR–RAR with comparable affinities to DR5 elements but DR0 elements do not act transcriptionally as independent RAREs. In this work, we present structural insights for the recognition of DR5 and DR0 elements by RXR–RAR heterodimer using x-ray crystallography, small angle x-ray scattering, and hydrogen/deuterium exchange coupled to mass spectrometry. We solved the crystal structures of RXR–RAR DNA-binding domain in complex with the Rarb2 DR5 and RXR–RXR DNA-binding domain in complex with Hoxb13 DR0. While cooperative binding was observed on DR5, the two molecules bound non-cooperatively on DR0 on opposite sides of the DNA. In addition, our data unveil the structural organization and dynamics of the multi-domain RXR–RAR DNA complexes providing evidence for DNA-dependent allosteric communication between domains. Differential binding modes between DR0 and DR5 were observed leading to differences in conformation and structural dynamics of the multi-domain RXR–RAR DNA complexes. These results reveal that the topological organization of the RAR binding element confer regulatory information by modulating the overall topology and structural dynamics of the RXR–RAR heterodimers.