How nuclear receptors transition between active and inactive forms: An energetic perspective.

Nuclear receptors regulate transcriptional programs in response to the binding of natural and synthetic ligands. These ligands modulate the receptor by inducing dynamic changes in the ligand binding domain that shift the C-terminal helix (H12) between active and inactive conformations. Despite decades of study, many questions persist regarding the nature of the inactive state and how ligands shift receptors between different states. Here, we use molecular dynamics (MD) simulations to investigate the timescale and energetic landscape of the conformational transition between inactive and active forms of progesterone receptor (PR) bound to a partial agonist. We observe that the microsecond timescale is insufficient to observe any transitions; only at millisecond timescales achieved via accelerated MD simulations do we find the inactive PR switches to the active state. Energetic analysis reveals that both active and inactive PR states represent energy minima separated by a barrier that can be traversed. In contrast, little or no transition is observed between active and inactive states when an agonist or antagonist is bound, confirming that ligand identity plays a key role in defining the energy landscape of nuclear receptor conformations.

Structure and assembly of the diiron cofactor in the heme-oxygenase-like domain of the N-nitrosourea-producing enzyme SznF.

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates Nδ and Nω' of Nω-methyl-l-arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme's C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase-like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase-like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O2-reactive Fe2(II/II) complex and the structural basis for cofactor instability-a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon-hydrogen and carbon-carbon bonds to install olefins. Among ∼9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ∼25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

Identification of an Evolutionarily Conserved Allosteric Network in Steroid Receptors.

Allosteric pathways in proteins describe networks comprising amino acid residues which may facilitate the propagation of signals between distant sites. Through inter-residue interactions, dynamic and conformational changes can be transmitted from the site of perturbation to an allosteric site. While sophisticated computational methods have been developed to characterize such allosteric pathways linking specific sites on proteins, few attempts have been made to apply these approaches toward identifying new allosteric sites. Here, we use molecular dynamics simulations and suboptimal path analysis to discover new allosteric networks in steroid receptors with a focus on evolutionarily conserved pathways. Using modern receptors and a reconstructed ancestral receptor, we identify networks connecting several sites to the activation function surface 2 (AF-2), the site of coregulator recruitment. One of these networks is conserved across the entire family, connecting a predicted allosteric site located between helices 9 and 10 of the ligand-binding domain. We investigate the basis of this conserved network as well as the importance of this site, discovering that the site lies in a region of the ligand-binding domain characterized by conserved inter-residue contacts. This study suggests an evolutionarily importance of the helix 9-helix 10 site in steroid receptors and identifies an approach that may be applied to discover previously unknown allosteric sites in proteins.

Allosteric regulation of thioesterase superfamily member 1 by lipid sensor domain binding fatty acids and lysophosphatidylcholine.

Nonshivering thermogenesis occurs in brown adipose tissue to generate heat in response to cold ambient temperatures. Thioesterase superfamily member 1 (Them1) is transcriptionally up-regulated in brown adipose tissue upon exposure to the cold and suppresses thermogenesis in order to conserve energy reserves. It hydrolyzes long-chain fatty acyl-CoAs that are derived from lipid droplets, preventing their use as fuel for thermogenesis. In addition to its enzymatic domains, Them1 contains a C-terminal StAR-related lipid transfer (START) domain with unknown ligand or function. By complementary biophysical approaches, we show that the START domain binds to long-chain fatty acids, products of Them1's enzymatic reaction, as well as lysophosphatidylcholine (LPC), lipids shown to activate thermogenesis in brown adipocytes. Certain fatty acids stabilize the START domain and allosterically enhance Them1 catalysis of acyl-CoA, whereas 18:1 LPC destabilizes and inhibits activity, which we verify in cell culture. Additionally, we demonstrate that the START domain functions to localize Them1 near lipid droplets. These findings define the role of the START domain as a lipid sensor that allosterically regulates Them1 activity and spatially localizes it in proximity to the lipid droplet.

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.

Cutting in-line with iron: ribosomal function and non-oxidative RNA cleavage.

Divalent metal cations are essential to the structure and function of the ribosome. Previous characterizations of the ribosome performed under standard laboratory conditions have implicated Mg2+ as a primary mediator of ribosomal structure and function. Possible contributions of Fe2+ as a ribosomal cofactor have been largely overlooked, despite the ribosome's early evolution in a high Fe2+ environment, and the continued use of Fe2+ by obligate anaerobes inhabiting high Fe2+ niches. Here, we show that (i) Fe2+ cleaves RNA by in-line cleavage, a non-oxidative mechanism that has not previously been shown experimentally for this metal, (ii) the first-order in-line rate constant with respect to divalent cations is >200 times greater with Fe2+ than with Mg2+, (iii) functional ribosomes are associated with Fe2+ after purification from cells grown under low O2 and high Fe2+ and (iv) a small fraction of Fe2+ that is associated with the ribosome is not exchangeable with surrounding divalent cations, presumably because those ions are tightly coordinated by rRNA and deeply buried in the ribosome. In total, these results expand the ancient role of iron in biochemistry and highlight a possible new mechanism of iron toxicity.

Iron mediates catalysis of nucleic acid processing enzymes: support for Fe(II) as a cofactor before the great oxidation event.

Life originated in an anoxic, Fe2+-rich environment. We hypothesize that on early Earth, Fe2+ was a ubiquitous cofactor for nucleic acids, with roles in RNA folding and catalysis as well as in processing of nucleic acids by protein enzymes. In this model, Mg2+ replaced Fe2+ as the primary cofactor for nucleic acids in parallel with known metal substitutions of metalloproteins, driven by the Great Oxidation Event. To test predictions of this model, we assay the ability of nucleic acid processing enzymes, including a DNA polymerase, an RNA polymerase and a DNA ligase, to use Fe2+ in place of Mg2+ as a cofactor during catalysis. Results show that Fe2+ can indeed substitute for Mg2+ in catalytic function of these enzymes. Additionally, we use calculations to unravel differences in energetics, structures and reactivities of relevant Mg2+ and Fe2+ complexes. Computation explains why Fe2+ can be a more potent cofactor than Mg2+ in a variety of folding and catalytic functions. We propose that the rise of O2 on Earth drove a Fe2+ to Mg2+ substitution in proteins and nucleic acids, a hypothesis consistent with a general model in which some modern biochemical systems retain latent abilities to revert to primordial Fe2+-based states when exposed to pre-GOE conditions.

Structural Analysis of the Glucocorticoid Receptor Ligand-Binding Domain in Complex with Triamcinolone Acetonide and a Fragment of the Atypical Coregulator, Small Heterodimer Partner.

The synthetic glucocorticoids (GCs) dexamethasone, mometasone furoate, and triamcinolone acetonide are pharmaceutical mainstays to treat chronic inflammatory diseases. These drugs bind to the glucocorticoid receptor (GR), a ligand-activated transcription factor and member of the nuclear receptor superfamily. The GR is widely recognized as a therapeutic target for its ability to counter proinflammatory signaling. Despite the popularity of GCs in the clinic, long-term use leads to numerous side effects, driving the need for new and improved drugs with less off-target pharmacology. X-ray crystal structures have played an important role in the drug-design process, permitting the characterization of robust structure-function relationships. However, steroid receptor ligand-binding domains (LBDs) are inherently unstable, and their crystallization requires extensive mutagenesis to enhance expression and crystallization. Here, we use an ancestral variant of GR as a tool to generate a high-resolution crystal structure of GR in complex with the potent glucocorticoid triamcinolone acetonide (TA) and a fragment of the small heterodimer partner (SHP). Using structural analysis, molecular dynamics, and biochemistry, we show that TA increases intramolecular contacts within the LBD to drive affinity and enhance stability of the receptor-ligand complex. These data support the emerging theme that ligand-induced receptor conformational dynamics at the mouth of the pocket play a major role in steroid receptor activation. This work also represents the first GR structure in complex with SHP, which has been suggested to play a role in modulating hepatic GR function.

Structure and Dynamics of the Liver Receptor Homolog 1-PGC1α Complex.

Peroxisome proliferator-activated gamma coactivator 1-α (PGC1α) regulates energy metabolism by directly interacting with transcription factors to modulate gene expression. Among the PGC1α binding partners is liver receptor homolog 1 (LRH-1; NR5A2), an orphan nuclear hormone receptor that controls lipid and glucose homeostasis. Although PGC1α is known to bind and activate LRH-1, mechanisms through which PGC1α changes LRH-1 conformation to drive transcription are unknown. Here, we used biochemical and structural methods to interrogate the LRH-1-PGC1α complex. Purified, full-length LRH-1, as well as isolated ligand binding domain, bound to PGC1α with higher affinity than to the coactivator, nuclear receptor coactivator-2 (Tif2), in coregulator peptide recruitment assays. We present the first crystal structure of the LRH-1-PGC1α complex, which depicts several hydrophobic contacts and a strong charge clamp at the interface between these partners. In molecular dynamics simulations, PGC1α induced correlated atomic motion throughout the entire LRH-1 activation function surface, which was dependent on charge-clamp formation. In contrast, Tif2 induced weaker signaling at the activation function surface than PGC1α but promoted allosteric signaling from the helix 6/ß-sheet region of LRH-1 to the activation function surface. These studies are the first to probe mechanisms underlying the LRH-1-PGC1α interaction and may illuminate strategies for selective therapeutic targeting of PGC1α-dependent LRH-1 signaling pathways.

Crystal Structures of the Nuclear Receptor, Liver Receptor Homolog 1, Bound to Synthetic Agonists.

Liver receptor homolog 1 (NR5A2, LRH-1) is an orphan nuclear hormone receptor that regulates diverse biological processes, including metabolism, proliferation, and the resolution of endoplasmic reticulum stress. Although preclinical and cellular studies demonstrate that LRH-1 has great potential as a therapeutic target for metabolic diseases and cancer, development of LRH-1 modulators has been difficult. Recently, systematic modifications to one of the few known chemical scaffolds capable of activating LRH-1 failed to improve efficacy substantially. Moreover, mechanisms through which LRH-1 is activated by synthetic ligands are entirely unknown. Here, we use x-ray crystallography and other structural methods to explore conformational changes and receptor-ligand interactions associated with LRH-1 activation by a set of related agonists. Unlike phospholipid LRH-1 ligands, these agonists bind deep in the pocket and do not interact with residues near the mouth nor do they expand the pocket like phospholipids. Unexpectedly, two closely related agonists with similar efficacies (GSK8470 and RJW100) exhibit completely different binding modes. The dramatic repositioning is influenced by a differential ability to establish stable face-to-face π-π-stacking with the LRH-1 residue His-390, as well as by a novel polar interaction mediated by the RJW100 hydroxyl group. The differing binding modes result in distinct mechanisms of action for the two agonists. Finally, we identify a network of conserved water molecules near the ligand-binding site that are important for activation by both agonists. This work reveals a previously unappreciated complexity associated with LRH-1 agonist development and offers insights into rational design strategies.

Simulations reveal the flexible domain architecture of full-length nuclear receptor complexes.

Nuclear receptors are multidomain transcription factors whose full-length quaternary architecture is poorly described and understood. Most nuclear receptors bind DNA as heterodimers or homodimers, which could encompass a variety of arrangements of the individual domains. Only a handful of experimental structures currently exist describing these architectures. Given that domain interactions and protein-DNA interactions within transcriptional complexes are tightly linked to function, understanding the arrangement of nuclear receptor domains on DNA is of utmost importance. Here, we employ modeling and molecular dynamics (MD) simulations to describe the structure of the full-length farnesoid X receptor (FXR) and retinoid X receptor alpha (RXR) heterodimer bound to DNA. Using over 100 microseconds of atomistic MD simulations, we characterize the dynamic behavior of eight FXR-RXR-DNA complexes, showing that these complexes support a range of quaternary architectures. We reveal the role of DNA binding and the hinge linkers in diversifying domain arrangements, roles that have been hard to appreciate previously due to experimental limitations in studying the flexible hinge. These studies provide a much-needed framework that will enable the field to obtain a complete understanding of nuclear receptor quaternary architectures.

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations.

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Nuclear receptor interdomain communication is mediated by the hinge with ligand specificity.

Nuclear receptors are ligand-induced transcription factors that bind directly to target genes and regulate their expression. Ligand binding initiates conformational changes that propagate to other domains, allosterically regulating their activity. The nature of this interdomain communication in nuclear receptors is poorly understood, largely owing to the difficulty of experimentally characterizing full-length structures. We have applied computational modeling approaches to describe and study the structure of the full length farnesoid X receptor (FXR), approximated by the DNA binding domain (DBD) and ligand binding domain (LBD) connected by the flexible hinge region. Using extended molecular dynamics simulations (> 10 microseconds) and enhanced sampling simulations, we provide evidence that ligands selectively induce domain rearrangement, leading to interdomain contact. We use protein-protein interaction assays to provide experimental evidence of these interactions, identifying a critical role of the hinge in mediating interdomain contact. Our results illuminate previously unknown aspects of interdomain communication in FXR and provide a framework to enable characterization of other full length nuclear receptors.

Accelerating the discovery of alkyl halide-derived natural products using halide depletion.

Even in the genomic era, microbial natural product discovery workflows can be laborious and limited in their ability to target molecules with specific structural features. Here we leverage an understanding of biosynthesis to develop a workflow that targets the discovery of alkyl halide-derived natural products by depleting halide anions, a key biosynthetic substrate for enzymatic halogenation, from microbial growth media. By comparing the metabolomes of bacterial cultures grown in halide-replete and deficient media, we rapidly discovered the nostochlorosides, the products of an orphan halogenase-encoding gene cluster from Nostoc punctiforme ATCC 29133. We further found that these products, a family of unusual chlorinated glycolipids featuring the rare sugar gulose, are polymerized via an unprecedented enzymatic etherification reaction. Together, our results highlight the power of leveraging an understanding of biosynthetic logic to streamline natural product discovery.

Ancient and modern mechanisms compete in progesterone receptor activation.

The progesterone receptor (PR) belongs to the steroid receptor family of ligand-regulated transcription factors, controlling genes important for development, metabolism, and reproduction. Understanding how diverse ligands bind and modulate PR activity will illuminate the design of ligands that control PR-driven signaling pathways. Here, we use molecular dynamics simulations to investigate how PR dynamics are altered by functionally diverse ligands. Using a library of 33 steroidal ligands that range from inactive to EC50 < 0.1 nM, we reveal an unexpected evolutionary basis for the wide gamut of activation. While other oxosteroid receptors employ an evolutionarily conserved mechanism dependent on a hydrogen bond between the receptor and ligand, extant PR has evolved a preference for activation that is not reliant on this polar interaction. We demonstrate that potent ligands utilize the modern PR mechanism while weaker ligands coopt the defunct ancestral mechanism by forming hydrogen bonds with Asn719. Based on their structures and dynamic signatures, ligands partition into four classes (inactive, weak, moderate and high potency) that interact distinctly with the PR binding pocket. Further, we use luciferase reporter assays and PR mutants to probe the roles of pocket residues in mediating distinct PR mechanisms. This combination of MD simulations and in vitro studies provide insight into how the evolutionary history of PR shapes its response to diverse ligands.

The nuclear receptor LRH-1 discriminates between ligands using distinct allosteric signaling circuits.

Nuclear receptors (NRs) are transcription factors that regulate essential biological processes in response to cognate ligands. An important part of NR function involves ligand-induced conformational changes that recruit coregulator proteins to the activation function surface (AFS), ~15 Å away from the ligand binding pocket. Ligands must communicate with the AFS to recruit appropriate coregulators and elicit different transcriptional outcomes, but this communication is poorly understood. These studies illuminate allosteric communication networks underlying activation of liver receptor homolog-1 (LRH-1), a NR that regulates development, metabolism, cancer progression and intestinal inflammation. Using >100 microseconds of all-atom molecular dynamics simulations involving 69 LRH-1 complexes, we identify distinct signaling circuits used by active and inactive ligands for AFS communication. Inactive ligands communicate via strong, coordinated motions along paths through the receptor to the AFS. Activating ligands disrupt the "inactive" circuit by inducing connectivity elsewhere. Ligand-contacting residues in helix 7 help mediate the switch between circuits, suggesting new avenues for developing LRH-1-targeted therapeutics. We also elucidate aspects of coregulator signaling, showing that localized, destabilizing fluctuations are induced by inappropriate ligand-coregulator pairings. These studies have uncovered novel features of LRH-1 allostery, and the quantitative approach used to analyze many simulations provides a framework to study allosteric signaling in other receptors.

The nuclear receptor LRH-1 discriminates between ligands using distinct allosteric signaling circuits.

Nuclear receptors (NRs) are transcription factors that regulate essential biological processes in response to cognate ligands. An important part of NR function involves ligand-induced conformational changes that recruit coregulator proteins to the activation function surface (AFS), ~15 Å away from the ligand-binding pocket. Ligands must communicate with the AFS to recruit appropriate coregulators and elicit different transcriptional outcomes, but this communication is poorly understood. These studies illuminate allosteric communication networks underlying activation of liver receptor homolog-1 (LRH-1), a NR that regulates development, metabolism, cancer progression, and intestinal inflammation. Using >100 µs of all-atom molecular dynamics simulations involving 74 LRH-1 complexes, we identify distinct signaling circuits used by active and inactive ligands for AFS communication. Inactive ligands communicate via strong, coordinated motions along paths through the receptor to the AFS. Activating ligands disrupt the "inactive" circuit and induce connectivity with a second allosteric site. Ligand-contacting residues in helix 7 help mediate the switch between circuits, suggesting new avenues for developing LRH-1-targeted therapeutics. We also elucidate aspects of coregulator signaling, showing that localized, destabilizing fluctuations are induced by inappropriate ligand-coregulator pairings. These studies have uncovered novel features of LRH-1 allostery, and the quantitative approach used to analyze many simulations provides a framework to study allosteric signaling in other receptors.

Bile acids and the gut microbiota: metabolic interactions and impacts on disease.

Despite decades of bile acid research, diverse biological roles for bile acids have been discovered recently due to developments in understanding the human microbiota. As additional bacterial enzymes are characterized, and the tools used for identifying new bile acids become increasingly more sensitive, the repertoire of bile acids metabolized and/or synthesized by bacteria continues to grow. Additionally, bile acids impact microbiome community structure and function. In this Review, we highlight how the bile acid pool is manipulated by the gut microbiota, how it is dependent on the metabolic capacity of the bacterial community and how external factors, such as antibiotics and diet, shape bile acid composition. It is increasingly important to understand how bile acid signalling networks are affected in distinct organs where the bile acid composition differs, and how these networks impact infectious, metabolic and neoplastic diseases. These advances have enabled the development of therapeutics that target imbalances in microbiota-associated bile acid profiles.

Two enhancer binding proteins activate σ54-dependent transcription of a quorum regulatory RNA in a bacterial symbiont.

To colonize a host, bacteria depend on an ensemble of signaling systems to convert information about the various environments encountered within the host into specific cellular activities. How these signaling systems coordinate transitions between cellular states in vivo remains poorly understood. To address this knowledge gap, we investigated how the bacterial symbiont Vibrio fischeri initially colonizes the light organ of the Hawaiian bobtail squid Euprymna scolopes. Previous work has shown that the small RNA Qrr1, which is a regulatory component of the quorum-sensing system in V. fischeri, promotes host colonization. Here, we report that transcriptional activation of Qrr1 is inhibited by the sensor kinase BinK, which suppresses cellular aggregation by V. fischeri prior to light organ entry. We show that Qrr1 expression depends on the alternative sigma factor σ54 and the transcription factors LuxO and SypG, which function similar to an OR logic gate, thereby ensuring Qrr1 is expressed during colonization. Finally, we provide evidence that this regulatory mechanism is widespread throughout the Vibrionaceae family. Together, our work reveals how coordination between the signaling pathways underlying aggregation and quorum-sensing promotes host colonization, which provides insight into how integration among signaling systems facilitates complex processes in bacteria.