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.

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.