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.

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.

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.

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.

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.

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.

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.

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.

Ligand-induced shifts in conformational ensembles that describe transcriptional activation

Nuclear receptors function as ligand-regulated transcription factors whose ability to regulate diverse physiological processes is closely linked with conformational changes induced upon ligand binding. Understanding how conformational populations of nuclear receptors are shifted by various ligands could illuminate strategies for the design of synthetic modulators to regulate specific transcriptional programs. Here, we investigate ligand-induced conformational changes using a reconstructed, ancestral nuclear receptor. By making substitutions at a key position, we engineer receptor variants with altered ligand specificities. We combine cellular and biophysical experiments to characterize transcriptional activity, as well as elucidate mechanisms underlying altered transcription in receptor variants. We then use atomistic molecular dynamics (MD) simulations with enhanced sampling to generate ensembles of wildtype and engineered receptors in combination with multiple ligands, followed by conformational analysis and correlation of MD-based predictions with functional ligand profiles. We determine that conformational ensembles accurately describe ligand responses based on observed population shifts. These studies provide a platform which will allow structural characterization of physiologically-relevant conformational ensembles, as well as provide the ability to design and predict transcriptional responses in novel ligands.

Differential Modulation of Nuclear Receptor LRH-1 through Targeting Buried and Surface Regions of the Binding Pocket

Liver receptor homologue-1 (LRH-1) is a phospholipid-sensing nuclear receptor that has shown promise as a target for alleviating intestinal inflammation and metabolic dysregulation in the liver. LRH-1 contains a large ligand-binding pocket, but generating synthetic modulators has been challenging. We have had recent success generating potent and efficacious agonists through two distinct strategies. We targeted residues deep within the pocket to enhance compound binding and residues at the mouth of the pocket to mimic interactions made by phospholipids. Here, we unite these two designs into one molecule to synthesize the most potent LRH-1 agonist to date. Through a combination of global transcriptomic, biochemical, and structural studies, we show that selective modulation can be driven through contacting deep versus surface polar regions in the pocket. While deep pocket contacts convey high affinity, contacts with the pocket mouth dominate allostery and provide a phospholipid-like transcriptional response in cultured cells.

A phospholipid mimetic targeting LRH-1 ameliorates colitis

Phospholipids are ligands for nuclear hormone receptors (NRs) that regulate transcriptional programs relevant to normal physiology and disease. Here, we demonstrate that mimicking phospholipid-NR interactions is a robust strategy to improve agonists of liver receptor homolog-1 (LRH-1), a therapeutic target for colitis. Conventional LRH-1 modulators only partially occupy the binding pocket, leaving vacant a region important for phospholipid binding and allostery. Therefore, we constructed a set of molecules with elements of natural phospholipids appended to a synthetic LRH-1 agonist. We show that the phospholipid-mimicking groups interact with the targeted residues in crystal structures and improve binding affinity, LRH-1 transcriptional activity, and conformational changes at a key allosteric site. The best phospholipid mimetic markedly improves colonic histopathology and disease-related weight loss in a murine T cell transfer model of colitis. This evidence of in vivo efficacy for an LRH-1 modulator in colitis represents a leap forward in agonist development.

Structural basis for an unprecedented enzymatic alkylation in cylindrocyclophane biosynthesis

The cyanobacterial enzyme CylK assembles the cylindrocyclophane natural products by performing two unusual alkylation reactions, forming new carbon–carbon bonds between aromatic rings and secondary alkyl halide substrates. This transformation is unprecedented in biology, and the structure and mechanism of CylK are unknown. Here, we report X-ray crystal structures of CylK, revealing a distinctive fusion of a Ca²⁺-binding domain and a β-propeller fold. We use a mutagenic screening approach to locate CylK’s active site at its domain interface, identifying two residues, Arg105 and Tyr473, that are required for catalysis. Anomalous diffraction datasets collected with bound bromide ions, a product analog, suggest that these residues interact with the alkyl halide electrophile. Additional mutagenesis and molecular dynamics simulations implicate Asp440 in activating the nucleophilic aromatic ring. Bioinformatic analysis of CylK homologs from other cyanobacteria establishes that they conserve these key catalytic amino acids, but they are likely associated with divergent reactivity and altered secondary metabolism. By gaining a molecular understanding of this unusual biosynthetic transformation, this work fills a gap in our understanding of how alkyl halides are activated and used by enzymes as biosynthetic intermediates, informing enzyme engineering, catalyst design, and natural product discovery.

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 Mg²⁺ as a primary mediator of ribosomal structure and function. Possible contributions of Fe²⁺ as a ribosomal cofactor have been largely overlooked, despite the ribosome's early evolution in a high Fe²⁺ environment, and the continued use of Fe²⁺ by obligate anaerobes inhabiting high Fe²⁺ niches. Here, we show that (i) Fe²⁺ 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 Fe²⁺ than with Mg²⁺, (iii) functional ribosomes are associated with Fe²⁺ after purification from cells grown under low O₂ and high Fe²⁺ and (iv) a small fraction of Fe²⁺ 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.

The Rho-GEF PIX-1 directs assembly or stability of lateral attachment structures between muscle cells

PIX proteins are guanine nucleotide exchange factors (GEFs) that activate Rac and Cdc42, and are known to have numerous functions in various cell types. Here, we show that a PIX protein has an important function in muscle. From a genetic screen in C. elegans, we found that pix-1 is required for the assembly of integrin adhesion complexes (IACs) at borders between muscle cells, and is required for locomotion of the animal. A pix-1 null mutant has a reduced level of activated Rac in muscle. PIX-1 localizes to IACs at muscle cell boundaries, M-lines and dense bodies. Mutations in genes encoding proteins at known steps of the PIX signaling pathway show defects at muscle cell boundaries. A missense mutation in a highly conserved residue in the RacGEF domain results in normal levels of PIX-1 protein, but a reduced level of activated Rac in muscle, and abnormal IACs at muscle cell boundaries.

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.

Development of the First Low Nanomolar Liver Receptor Homolog-1 Agonist through Structure-guided Design

As a key regulator of metabolism and inflammation, the orphan nuclear hormone receptor, liver receptor homolog-1 (LRH-1), has potential as a therapeutic target for diabetes, nonalcoholic fatty liver disease, and inflammatory bowel diseases (IBD). Discovery of LRH-1 modulators has been difficult, in part due to the tendency for synthetic compounds to bind unpredictably within the lipophilic binding pocket. Using a structure-guided approach, we exploited a newly discovered polar interaction to lock agonists in a consistent orientation. This enabled the discovery of the first low nanomolar LRH-1 agonist, one hundred times more potent than the best previous modulator. We elucidate a novel mechanism of action that relies upon specific polar interactions deep in the LRH-1 binding pocket. In an organoid model of IBD, the new agonist increases expression of LRH-1-controlled steroidogenic genes and promotes anti-inflammatory gene expression changes. These studies constitute major progress in developing LRH-1 modulators with potential clinical utility.

Comparison of the prevalence of adhesions at the time of diagnostic laparoscopy for infertility between patient who had open myomectomy and those who had no previous pelvic-abdominal surgery or pelvic inflammatory disease

Objective

The aim of this study is to determine the prevalence and patterns of adhesions in infertile women with prior open myomectomy compared with women without prior pelvic-abdominal surgery.

Methods

A nested case-control study of infertile women who had diagnostic laparoscopy after open myomectomy between January 2008 and June 2015 in Life Institute for Endoscopy Limited, Nnewi Nigeria, was conducted. At diagnostic laparoscopy, the presence, site, and quality (density) of adhesions was noted and recorded. Women with endometriosis and documented prior pelvic infections were excluded.

Results

Of the 348 women who were eligible, 121 women had prior open myomectomy (study group) while 227 had none (control group). Of the 121 women, adhesions were detected in 97 women; thus, the prevalence of adhesion after open myomectomy was 80.2% versus 20.3% (46/227) in controls (odds ratio [OR] =15.90; 95% confidence interval [CI] =8.86-28.76; P < 0.001). Compared with controls, adhesions were statistically higher in these sites: uterus (68.6% vs. 14.5%), fallopian tubes (77.7% vs. 16.7%), ovaries (62.8% vs. 20.3%), cul-de-sac (66.1% vs. 16.3%), bladder (45.5% vs. 7.9%), and bowel (53.7% vs. 6.2%) (P < 0.001, for all). Cohesive form of adhesions was also statistically higher, 54 (44.6%) vs. 32 (14.1%) (P < 0.05).

Conclusion

In infertile women, post-open myomectomy adhesions have a high prevalence of 80.2% with high predilection in fallopian tubes and uterus, compared with women without prior abdominopelvic surgery at 20.3%. Cohesive form of adhesion predominates.

The first crystal structure of a DNA-free nuclear receptor DNA binding domain sheds light on DNA-driven allostery in the glucocorticoid receptor

The glucocorticoid receptor (GR) is a steroid hormone receptor of the nuclear receptor family that regulates gene expression in response to glucocorticoid hormone signaling. Interaction with specific GR DNA binding sequences causes conformational changes in the GR DNA binding domain (DBD) that result in recruitment of specific sets of co-regulators that determine transcriptional outcomes. We have solved the crystal structure of GR DBD in its DNA-free state, the first such crystal structure from any nuclear receptor. In contrast to previous NMR structures, this crystal structure reveals that free GR DBD adopts a conformation very similar to DNA-bound states. The lever arm region is the most variable element in the free GR DBD. Molecular dynamics of the free GR DBD as well as GR DBD bound to activating and repressive DNA elements confirm lever arm flexibility in all functional states. Cluster analysis of lever arm conformations during simulations shows that DNA binding and dimerization cause a reduction in the number of conformations sampled by the lever arm. These results reveal that DNA binding and dimerization drive conformational selection in the GR DBD lever arm region and show how DNA allosterically controls GR structure and dynamics.

Structural and Dynamics Comparison of Thermostability in Ancient, Modern, and Consensus Elongation Factor Tus

Rationally engineering thermostability in proteins would create enzymes and receptors that function under harsh industrial applications. Several sequence-based approaches can generate thermostable variants of mesophilic proteins. To gain insight into the mechanisms by which proteins become more stable, we use structural and dynamic analyses to compare two popular approaches, ancestral sequence reconstruction (ASR) and the consensus method, used to generate thermostable variants of Elongation Factor Thermo-unstable (EF-Tu). We present crystal structures of ancestral and consensus EF-Tus, accompanied by molecular dynamics simulations aimed at probing the strategies employed to enhance thermostability. All proteins adopt crystal structures similar to extant EF-Tus, revealing no difference in average structure between the methods. Molecular dynamics reveals that ASR-generated sequences retain dynamic properties similar to extant, thermostable EF-Tu from Thermus aquaticus, while consensus EF-Tu dynamics differ from evolution-based sequences. This work highlights the advantage of ASR for engineering thermostability while preserving natural motions in multidomain proteins.

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.

Ribosomal small subunit domains radiate from a central core

The domain architecture of a large RNA can help explain and/or predict folding, function, biogenesis and evolution. We offer a formal and general definition of an RNA domain and use that definition to experimentally characterize the rRNA of the ribosomal small subunit. Here the rRNA comprising a domain is compact, with a self-contained system of molecular interactions. A given rRNA helix or stem-loop must be allocated uniquely to a single domain. Local changes such as mutations can give domain-wide effects. Helices within a domain have interdependent orientations, stabilities and interactions. With these criteria we identify a core domain (domain A) of small subunit rRNA. Domain A acts as a hub, linking the four peripheral domains and imposing orientational and positional restraints on the other domains. Experimental characterization of isolated domain A, and mutations and truncations of it, by methods including selective 2'OH acylation analyzed by primer extension and circular dichroism spectroscopy are consistent with our architectural model. The results support the utility of the concept of an RNA domain. Domain A, which exhibits structural similarity to tRNA, appears to be an essential core of the small ribosomal subunit.

RNA-magnesium-protein interactions in large ribosomal subunit

Some of the magnesium ions in the ribosome are coordinated by multiple rRNA phosphate groups. These magnesium ions link distal sequences of rRNA, primarily by incorporating phosphate groups into the first coordination shell. Less frequently, magnesium interacts with ribosomal proteins. Ribosomal protein L2 appears to be unique by forming specific magnesium-mediated interactions with rRNA. Using optimized models derived from X-ray structures, we subjected rRNA/magnesium/water/rProtein L2 assemblies to quantum mechanical analysis using the density functional theory and natural energy decomposition analysis. The combined results provide estimates of energies of formation of these assemblies, and allow us to decompose the energies of interaction. The results indicated that RNA immobilizes magnesium by multidentate chelation with phosphate, and that the magnesium ions in turn localize and polarize water molecules, increasing energies and specificities of interaction of these water molecules with L2 protein. Thus, magnesium plays subtle, yet important, roles in ribosomal assembly beyond neutralization of electrostatic repulsion and direct coordination of RNA functional groups.