The Promiscuity of Allosteric Regulation of Nuclear Receptors by Retinoid X Receptor

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.

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.

Statistical Analysis of Protein-Ligand Interaction Patterns in Nuclear Receptor RORγ

The receptor RORγ belongs to the nuclear receptor superfamily that senses small signaling molecules and regulates at the gene transcription level. Since RORγ has a high basal activity and plays an important role in immune responses, inhibitors targeting this receptor have been a focus for many studies. The receptor-ligand interaction is complex, and often subtle differences in ligand structure can determine its role as an inverse agonist or an agonist. We examined more than 130 existing RORγ crystal structures that have the same receptor complexed with different ligands. We reported the features of receptor-ligand interaction patterns and the differences between agonist and inverse agonist binding. Specific changes in the contact interaction map are identified to distinguish active and inactive conformations. Further statistical analysis of the contact interaction patterns using principal component analysis reveals a dominant mode which separates allosteric binding vs. canonical binding and a second mode which may indicate active vs. inactive structures. We also studied the nature of constitutive activity by performing a 100-ns computer simulation of apo RORγ. Using constitutively active nuclear receptor CAR as a comparison, we identified a group of conserved contacts that have similar contact strength between the two receptors. These conserved contact interactions, especially a couple key contacts in H11–H12 interaction, can be considered essential to the constitutive activity of RORγ. These protein-ligand and internal protein contact interactions can be useful in the development of new drugs that direct receptor activity.

Discovery of a "Gatekeeper" Antagonist that Blocks Entry Pathway to Retinoid X Receptors (RXRs) without Allosteric Ligand Inhibition in Permissive RXR Heterodimers

Retinoid X receptor (RXR) heterodimers such as PPAR/RXR, LXR/RXR, and FXR/RXR can be activated by RXR agonists alone and are therefore designated as permissive. Similarly, existing RXR antagonists show allosteric antagonism toward partner receptor agonists in these permissive RXR heterodimers. Here, we show 1-(3-(2-ethoxyethoxy)-5,5,8,8-tetramethyl-5,6,7,8-tetrahydronaphthalen-2-yl)-2-(trifluoromethyl)-1H-benzo[d]imidazole-5-carboxylic acid (14, CBTF-EE) as the first RXR antagonist that does not show allosteric inhibition in permissive RXR heterodimers. This compound was designed based on the hypothesis that RXR antagonists that do not induce conformational changes of RXR would not exhibit such allosteric inhibition. CD spectra and X-ray co-crystallography of the complex of 14 and the RXR ligand binding domain (LBD) confirmed that 14 does not change the conformation of hRXR-LBD. The X-ray structure analysis revealed that 14 binds at the entrance of the ligand binding pocket (LBP), blocking access to the LBP and thus serving as a "gatekeeper".

Protein promiscuity in drug discovery, drug-repurposing and antibiotic resistance

Proteins are supposed to bind to their substrates/ligands in a specific manner via their pre-formed binding sites, according to classical biochemistry. In recent years, several types of deviations from this norm have been observed and called promiscuous behavior. Enzymatic promiscuities allow several biochemical functions to be carried out by the same enzyme. The promiscuous activity can also be the origin of "new proteins" via gene duplication. In more recent years, proteins from prokaryotes, eukaryotes and viruses have been found to have intrinsic disorder and lack a preformed binding site. Intrinsic disorder is exploited in regulatory proteins such as those that are involved in transcription and signal transduction. Such proteins function by folding locally while binding to their ligands or interacting with other proteins. These phenomena have also been classified as examples of protein promiscuity and encompass diverse kinds of ligands that can bind to a protein. Given the significant extent of structural homology in many protein families, it is not surprising that ligands also have been found to display promiscuity. Promiscuous behavior of proteins offers both challenges and opportunities to the drug discovery programs such as drug repurposing. Pathogens when exposed to antibiotics exploit protein promiscuity in several ways to develop resistance to the drug. There is increasing evidence now to support that the disorder in proteins is a major tool used by pathogens for virulence and evade drug action by exploiting protein promiscuity. This review provides a holistic view of this multi-faceted phenomenon called protein promiscuity.

Transcriptome analyses of liver in newly-hatched chicks during the metabolic perturbation of fasting and re-feeding reveals THRSPA as the key lipogenic transcription factor

Background

The fasting-refeeding perturbation has been used extensively to reveal specific genes and metabolic pathways that control energy metabolism in the chicken. Most global transcriptional scans of the fasting-refeeding response in liver have focused on juvenile chickens that were 1, 2 or 4 weeks old. The present study was aimed at the immediate post-hatch period, in which newly-hatched chicks were subjected to fasting for 4, 24 or 48 h, then refed for 4, 24 or 48 h, and compared with a fully-fed control group at each age (D1-D4).

Results

Visual analysis of hepatic gene expression profiles using hierarchical and K-means clustering showed two distinct patterns, genes with higher expression during fasting and depressed expression upon refeeding and those with an opposing pattern of expression, which exhibit very low expression during fasting and more abundant expression with refeeding. Differentially-expressed genes (DEGs), identified from five prominent pair-wise contrasts of fed, fasted and refed conditions, were subjected to Ingenuity Pathway Analysis. This enabled mapping of analysis-ready (AR)-DEGs to canonical and metabolic pathways controlled by distinct gene interaction networks. The largest number of hepatic DEGs was identified by two contrasts: D2FED48h/D2FAST48h (968 genes) and D2FAST48h/D3REFED24h (1198 genes). The major genes acutely depressed by fasting and elevated upon refeeding included ANGTPL, ATPCL, DIO2, FASN, ME1, SCD, PPARG, SREBP2 and THRSPA—a primary lipogenic transcription factor. In contrast, major lipolytic genes were up-regulated by fasting or down-regulated after refeeding, including ALDOB, IL-15, LDHB, LPIN2, NFE2L2, NR3C1, NR0B1, PANK1, PPARA, SERTAD2 and UPP2.

Conclusions

Transcriptional profiling of liver during fasting/re-feeding of newly-hatched chicks revealed several highly-expressed upstream regulators, which enable the metabolic switch from fasted (lipolytic/gluconeogenic) to fed or refed (lipogenic/thermogenic) states. This rapid homeorhetic shift of whole-body metabolism from a catabolic-fasting state to an anabolic-fed state appears precisely orchestrated by a small number of ligand-activated transcription factors that provide either a fasting-lipolytic state (PPARA, NR3C1, NFE2L2, SERTAD2, FOX01, NR0B1, RXR) or a fully-fed and refed lipogenic/thermogenic state (THRSPA, SREBF2, PPARG, PPARD, JUN, ATF3, CTNNB1). THRSPA has emerged as the key transcriptional regulator that drives lipogenesis and thermogenesis in hatchling chicks, as shown here in fed and re-fed states.

Regulatory Mechanics of Constitutive Androstane Receptors: Basal and Ligand-Directed Actions

Constitutive androstane receptor (CAR) is a nuclear hormone receptor that primarily functions in sensing and metabolizing xenobiotics. The basal activity of this receptor is relatively high, and CAR is deemed active in the absence of ligand. The (over)activation can promote drug toxicity and tumor growth. Thus, therapeutic treatments seek inverse agonists to inhibit or modulate CAR activities. To advance our understanding of the regulatory mechanisms of CAR, we used computational and experimental approaches to elucidate three aspects of CAR activation and inactivation: (1) ligand-dependent actions, (2) ligand-orthologue specificity, and (3) constitutive activity. For ligand-dependent actions, we examined the ligand-bound simulations and identified two sets of ligand-induced contacts promoting CAR activation via coactivator binding (H11-H12 contact) or inactivation via corepressor binding (H4-H11 contact). For orthologue specificity, we addressed a puzzling fact that murine CAR (mCAR) and human CAR (hCAR) respond differently to the same ligand (CITCO), despite their high sequence homology. We found that the helix H7 of hCAR is responsible for a stronger binding of the ligand CITCO compared to mCAR, hence a stronger CITCO-induced activation. For basal activity, we reported computer-generated unliganded CAR structures and critical mutagenesis (mCAR's V209A and N333D) results of a cell-based transcription assay. Our results reveal that the basal conformation of CAR shares prominent features with the agonist-bound form, and helix HX has an important contribution to the constitutive activity. These findings altogether can be useful for the understanding of constitutively active receptors and the design of drug molecules targeting them.

Characterizing the 3D structure and dynamics of chromosomes and proteins in a common contact matrix framework

Conformational ensembles of biopolymers, whether proteins or chromosomes, can be described using contact matrices. Principal component analysis (PCA) on the contact data has been used to interrogate both protein and chromosome structures and/or dynamics. However, as these fields have developed separately, variants of PCA have emerged. Previously, a variant we hereby term Implicit-PCA (I-PCA) has been applied to chromosome contact matrices and revealed the spatial segregation of active and inactive chromatin. Separately, Explicit-PCA (E-PCA) has previously been applied to proteins and characterized their correlated structure fluctuations. Here, we swapped analysis methods (I-PCA and E-PCA), applying each to a different biopolymer type (chromosome or protein) than the one for which they were initially developed. We find that applying E-PCA to chromosome distance matrices derived from microscopy data can reveal the dominant motion (concerted fluctuation) of these chromosomes. Further, by applying E-PCA to Hi-C data across the human blood cell lineage, we isolated the aspects of chromosome structure that most strongly differentiate cell types. Conversely, when we applied I-PCA to simulation snapshots of proteins, the major component reported the consensus features of the structure, making this a promising approach for future analysis of semi-structured proteins.

Effector-Binding-Directed Dimerization and Dynamic Communication between Allosteric Sites of Ribonucleotide Reductase

Proteins forming dimers or larger complexes can be strongly influenced by their effector-binding status. We investigated how the effector-binding event is coupled with interface formation via computer simulations, and we quantified the correlation of two types of contact interactions: between the effector and its binding pocket and between protein monomers. This was achieved by connecting the protein dynamics at the monomeric level with the oligomer interface information. We applied this method to ribonucleotide reductase (RNR), an essential enzyme for de novo DNA synthesis. RNR contains two important allosteric sites, the s-site (specificity site) and the a-site (activity site), which bind different effectors. We studied these different binding states with atomistic simulation and used their coarse-grained contact information to analyze the protein dynamics. The results reveal that the effector-protein dynamics at the s-site and dimer interface formation are positively coupled. We further quantify the resonance level between these two events, which can be applied to other similar systems. At the a-site, different effector-binding states (ATP vs dATP) drastically alter the protein dynamics and affect the activity of the enzyme. On the basis of these results, we propose a new mechanism of how the a-site regulates enzyme activation.

The crustacean ecdysone cassette: A gatekeeper for molt and metamorphosis

Arthropods have long been utilized as models to explore molecular function, and the findings derived from them can be applied throughout metazoa, including as a basis for medical research. This has led to the adoption of many representative insect models beyond Drosophila, as each lends its own unique perspective to questions in endocrinology and genetics. However, non-insect arthropods are yet to be realised for the potential insight they may provide in such studies. The Crustacea are among the most ancient arthropods from which insects descended, comprising a huge variety of life histories and ecological roles. Of the events in a typical crustacean development, metamorphosis is perhaps the most ubiquitous, challenging and highly studied. Despite this, our knowledge of the endocrinology which underpins metamorphosis is rudimentary at best; although several key molecules have been identified and studied in depth, the link between them is quite nebulous and leans heavily on well-explored insect models, which diverged from the Pancrustacea over 450 million years ago. As omics technologies become increasingly accessible, they bring the prospect of explorative molecular research which will allow us to uncover components and pathways unique to crustaceans. This review reconciles known components of crustacean metamorphosis and reflects on our findings in insects to outline a future search space, with focus given to the ecdysone cascade. To expand our knowledge of this ubiquitous endocrine system not only aids in our understanding of crustacean metamorphosis, but also provides a deeper insight into the adaptive capacity of arthropods throughout evolution.

Characterizing protein conformations by correlation analysis of coarse-grained contact matrices

We have developed a method to capture the essential conformational dynamics of folded biopolymers using statistical analysis of coarse-grained segment-segment contacts. Previously, the residue-residue contact analysis of simulation trajectories was successfully applied to the detection of conformational switching motions in biomolecular complexes. However, the application to large protein systems (larger than 1000 amino acid residues) is challenging using the description of residue contacts. Also, the residue-based method cannot be used to compare proteins with different sequences. To expand the scope of the method, we have tested several coarse-graining schemes that group a collection of consecutive residues into a segment. The definition of these segments may be derived from structural and sequence information, while the interaction strength of the coarse-grained segment-segment contacts is a function of the residue-residue contacts. We then perform covariance calculations on these coarse-grained contact matrices. We monitored how well the principal components of the contact matrices is preserved using various rendering functions. The new method was demonstrated to assist the reduction of the degrees of freedom for describing the conformation space, and it potentially allows for the analysis of a system that is approximately tenfold larger compared with the corresponding residue contact-based method. This method can also render a family of similar proteins into the same conformational space, and thus can be used to compare the structures of proteins with different sequences.

Allosteric pathways in nuclear receptors - Potential targets for drug design

The nuclear receptor family of transcription factor proteins mediates endocrine function and plays critical roles in the development, physiology and pharmacology. Malfunctioning nuclear receptors are associated with several disease states. The functional activity of nuclear receptors is regulated by small molecular hormonal and synthetic molecules. Multiple sources of evidence have identified and distinguished between the different allosteric pathways initiated by ligands, DNA and cofactors such as co-activators and co-repressors. Also, these biophysical studies are attempting to determine how these pathways that regulate co-activator and DNA recognition can control gene transcription. Thus, there is a growing interest in determining the genome-scale impact of allostery in nuclear receptors. Today, it is accepted that a detailed understanding of the allosteric regulatory pathways within the nuclear receptor molecular complex will enable the development of efficient drug therapies in the long term.

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue-residue associations

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue-residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.