X-ray structures of progesterone receptor ligand binding domain in its agonist state reveal differing mechanisms for mixed profiles of 11β-substituted steroids

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.

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.

Determination of Molecule Category of Ligands Targeting the Ligand-Binding Pocket of Nuclear Receptors with Structural Elucidation and Machine Learning

The mechanism of transcriptional activation/repression of the nuclear receptors (NRs) involves two main conformations of the NR protein, namely, the active (agonistic) and inactive (antagonistic) conformations. Binding of agonists or antagonists to the ligand-binding pocket (LBP) of NRs can regulate the downstream signaling pathways with different physiological effects. However, it is still hard to determine the molecular type of a LBP-bound ligand because both the agonists and antagonists bind to the same position of the protein. Therefore, it is necessary to develop precise and efficient methods to facilitate the discrimination of agonists and antagonists targeting the LBP of NRs. Here, combining structural and energetic analyses with machine-learning (ML) algorithms, we constructed a series of structure-based ML models to determine the molecular category of the LBP-bound ligands. We show that the proposed models work robustly and with high accuracy (ACC > 0.9) for determining the category of molecules derived from docking-based and crystallized poses. Furthermore, the models are also capable of determining the molecular category of ligands with dual opposite functions on different NRs (i.e., working as an agonist in one NR target, whereas functioning as an antagonist in another) with reasonable accuracy. The proposed method is expected to facilitate the determination of the molecular properties of ligands targeting the LBP of NRs with structural interpretation.

Identification of ELK1 interacting peptide segments in the androgen receptor

Prostate cancer (PCa) growth requires tethering of the androgen receptor (AR) to chromatin by the ETS domain transcription factor ELK1 to coactivate critical cell proliferation genes. Disruption of the ELK1-AR complex is a validated potential means of therapeutic intervention in PCa. AR associates with ELK1 by co-opting its two ERK docking sites, through the amino-terminal domain (A/B domain) of AR. Using a mammalian two-hybrid assay, we have now functionally mapped amino acids within the peptide segments 358-457 and 514-557 in the A/B domain as required for association with ELK1. The mapping data was validated by GST (glutathione S-transferase)-pulldown and BRET (bioluminescence resonance energy transfer) assays. Comparison of the relative contributions of the interacting motifs/segments in ELK1 and AR to coactivation of ELK1 by AR suggested a parallel mode of binding of AR and ELK1 polypeptides. Growth of PCa cells was partially inhibited by deletion of the upstream segment in AR and nearly fully inhibited by deletion of the downstream segment. Our studies have identified two peptide segments in AR that mediate functional association of AR with its two docking sites in ELK1. Identification of the ELK1 recognition sites in AR should enable further structural studies of the ELK1-AR interaction and rational design of small molecule drugs to disrupt this interaction.

Structural overview and perspectives of the nuclear receptors, a major family as the direct targets for small-molecule drugs

The nuclear receptors (NRs) are an evolutionarily related family of transcription factors, which share certain common structural characteristics and regulate the expressions of various genes by recognizing different response elements. NRs play important roles in cell differentiation, proliferation, survival and apoptosis, rendering them indispensable in many physiological activities including growth and metabolism. As a result, dysfunctions of NRs are closely related to a variety of diseases, such as diabetes, obesity, infertility, inflammation, the Alzheimer's disease, cardiovascular diseases, prostate and breast cancers. Meanwhile, small-molecule drugs directly targeting NRs have been widely used in the treatment of above diseases. Here we summarize recent progress in the structural biology studies of NR family proteins. Compared with the dozens of structures of isolated DNA-binding domains (DBDs) and the striking more than a thousand of structures of isolated ligand-binding domains (LBDs) accumulated in the Protein Data Bank (PDB) over thirty years, by now there are only a small number of multi-domain NR complex structures, which reveal the integration of different NR domains capable of the allosteric signal transduction, or the detailed interactions between NR and various coregulator proteins. On the other hand, the structural information about several orphan NRs is still totally unavailable, hindering the further understanding of their functions. The fast development of new technologies in structural biology will certainly help us gain more comprehensive information of NR structures, inspiring the discovery of novel NR-targeting drugs with a new binding site beyond the classic LBD pockets and/or a new mechanism of action.

Identification of steroidal saponins from Tribulus terrestris and their in silico docking studies

Tribulus terrestris is known to possess many pharmacological properties, most notably its anticancer activities, owing to its rich steroidal saponin contents. Even though many reports are available elucidating the anticancer potential of the herb, we, for the very first time have attempted to isolate and identified the active compound present in seed crude saponin extract and confers its in silico docking ability with various cellular targets proteins. High performance thin layer chromatography confirms the presence of active saponins in leaf and seed saponin extracts which were further fractionated by silica gel column chromatography. Fractions collected were assessed for cytotoxicity on human breast cancer cells. High resolution liquid chromatography and mass spectroscopy was employ to identify the active components present in fraction with highest cytotoxicity. Intriguingly, Nautigenin type of steriodal saponin was identified to present in the active fraction of seed extract (SF11) and the identified compound was further analyzed for its in silico docking interaction using PyRx AutodockVina. Docking studies revealed the high binding affinity of Nuatigenin at significant sites with apoptotic proteins Bcl-2, Bax, caspase-3, caspase-8, p53 and apoptosis inducing factor along with cell surface receptors estrogen receptor, projesterone receptor, epidermal growth factor receptor, and human epidermal growth factor receptor-2. Thus, the conclusions were drawn that saponin fraction of Tribulus terrestis possesses active compounds having anticancer property and specifically, Nuatigenin saponin can be considered as an important therapeutic drug for the breast cancer treatment.

Functional Effects In Silico Prediction for Androgen Receptor Ligand-Binding Domain Novel I836S Mutation

Over 1000 mutations are described in the androgen receptor (AR) gene. Of those, about 600 were found in androgen insensitivity syndrome (AIS) patients, among which 400 mutations affect the ligand-binding domain (LBD) of the AR protein. Recently, we reported a novel missense mutation c.2507T>G I836S (ClinVarID: 974911) in a patient with complete AIS (CAIS) phenotype. In the present study, we applied a set of computational approaches for the structural analysis of the ligand-binding domains in a wild-type and mutant AR to evaluate the functional impact of the novel I836S mutation. We revealed that the novel I836S substitution leads to a shorter existence time of the ligand's gating tunnel and internal cavity, occurring only in the presence of S836 phosphorylation. Additionally, the analysis of phosphorylation of the 836 mutant residues explained the negative impact on AR homodimerization, since monomer surface changes indirectly impacted the binding site. Our analyses provide evidence that I836S causes disruptions of AR protein functionality and development of CAIS clinical features in patients.

Unraveling the Allosteric Cross-Talk between the Coactivator Peptide and the Ligand-Binding Site in the Glucocorticoid Receptor

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known allosteric cross-talk between the ligand and coregulator binding sites within the GR ligand-binding domain that is crucial for the control of the functional response. However, the molecular mechanisms underlying such an allosteric control remain elusive. Here, molecular dynamics (MD) simulations, bioinformatic analysis, and biophysical measurements are integrated to capture the structural and dynamic features of the allosteric cross-talk within the GR. We identified a network of evolutionarily conserved residues that enables the allosteric signal transduction, in agreement with experimental data. MD simulations clarify how such a network is dynamically interconnected and offer a mechanistic explanation of how different peptides affect the intensity of the allosteric signal. This study provides useful insights to elucidate the GR allosteric regulation, ultimately providing a foundation for designing novel drugs.

EC313-a tissue selective SPRM reduces the growth and proliferation of uterine fibroids in a human uterine fibroid tissue xenograft model

Uterine fibroids (UFs) are associated with irregular or excessive uterine bleeding, pelvic pain or pressure, or infertility. Ovarian steroid hormones support the growth and maintenance of UFs. Ulipristal acetate (UPA) a selective progesterone receptor (PR) modulator (SPRM) reduce the size of UFs, inhibit ovulation and lead to amenorrhea. Recent liver toxicity concerns with UPA, diminished enthusiasm for its use and reinstate the critical need for a safe, efficacious SPRM to treat UFs. In the current study, we evaluated the efficacy of new SPRM, EC313, for the treatment for UFs using a NOD-SCID mouse model. EC313 treatment resulted in a dose-dependent reduction in the fibroid xenograft weight (p < 0.01). Estradiol (E2) induced proliferation was blocked significantly in EC313-treated xenograft fibroids (p < 0.0001). Uterine weight was reduced by EC313 treatment compared to UPA treatment. ER and PR were reduced in EC313-treated groups compared to controls (p < 0.001) and UPA treatments (p < 0.01). UF specific desmin and collagen were markedly reduced with EC313 treatment. The partial PR agonism and no signs of unopposed estrogenicity makes EC313 a candidate for the long-term treatment for UFs. Docking studies have provided a structure based explanation for the SPRM activity of EC313.

Ligand induced dissociation of the AR homodimer precedes AR monomer translocation to the nucleus

The androgen receptor (AR) regulates male sexual development. We have now investigated AR homodimerization, hormone-dependent monomerization and nuclear translocation in PC-3 and COS-1 cells, by utilizing mutations associated with the androgen insensitivity syndrome: Pro767Ala, Phe765Leu, Met743Val and Trp742Arg. AR wild type (WT) was expressed as a homodimer in the cytoplasm, while none of these mutants formed homodimers. Unlike AR WT which responded to 1 nM dihydrotestosterone (DHT) to dissociate and translocate into the nucleus, AR Pro767Ala and Phe765Leu mutants remain as the monomer in the cytoplasm. In the crystal structure of the AR LBD homodimer, Pro767 and Phe765 reside closely on a loop that constitutes the dimer interface; their sidechains interact with the Pro767 of the other monomer and with the DHT molecule in the ligand-binding pocket. These observations place Phe765 at a position to facilitate DHT binding to Pro767 and lead to dissociation of the AR homodimer in the cytoplasm. This Pro-Phe Met relay may constitute a structural switch that mediates androgen signaling and is conserved in other steroid hormone receptors.

Structural Reassignment of rel-(3'Z,3R,6R,7R,3a'R,6'R)-3,8-Dihydrodiligustilide and the Activity of Diligustilide and 3,8-Dihydro- and 3,8,7',7a'-Tetrahydrodiligustilides as Progestins

Several phthalides were semisynthesized, including a 3,8-dihydrodiligustilide with progesterone-like activity, previously isolated from Ligusticum chuanxiong, the structure of which was earlier assigned to a semisynthetic product with nonidentical spectroscopic constants. The structure of this natural phthalide was reassigned with a proposal of its absolute configuration. Phthalides acted as progestins in cell viability assays, immunofluorescence microscopy, and docking analysis. Therefore, the structures for natural and semisynthetic phthalides with potential use in hormone-related therapies were reassigned.

Ligand Binding Induces Agonistic-Like Conformational Adaptations in Helix 12 of Progesterone Receptor Ligand Binding Domain

Progesterone receptor (PR) is a member of the nuclear receptor (NR) superfamily and plays a vital role in the female reproductive system. The malfunction of it would lead to several types of cancers. The understanding of conformational changes in its ligand binding domain (LBD) is valuable for both biological function studies and therapeutically intervenes. A key unsolved question is how the binding of a ligand (agonist, antagonist, or a selective modulator) induces conformational changes of PR LBD, especially its helix 12. We applied molecular dynamics (MD) simulations to explore the conformational adaptations of PR LBD with or without a ligand or the co-repressor peptides binding. From the simulations, both the agonist progesterone (P4) and the selective PR modulator (SPRM) asoprisnil induces agonistic-like helix 12 conformations (the "closed" states) in PR LBD and the complex of LBD-SPRM is less stable, comparing to the agonist-liganded PR LBD. The results, therefore, explain the partial agonism of the SPRM, which could induce weak agonistic effects in PR. We also found that co-repressor peptides could be stably associated with the LBD and stabilize the LBD in a "semi-open" state for helix 12. These findings would enhance our understanding of PR structural and functional relationships and would also be useful for future structure and knowledge-based drug discovery.

Progestin therapy to prevent preterm birth: History and effectiveness of current strategies and development of novel approaches

In the 1930s the "progestin" hormone produced by the corpus luteum was isolated and found to be a Δ4-keto-steroid. It was aptly named progesterone (P4) and in the following 30 years the capacity of P4 and derivatives to prevent preterm birth (PTB) was examined. Outcomes of multiple small studies suggested that progestin prophylaxis beginning at mid-gestation decreases the risk for PTB. Subsequent larger trials found that prophylaxis with weekly intramuscular injections of 17α-hydroxyprogesterone caproate (17HPC) beginning at mid-gestation decreased PTB risk in women with a history of PTB. Other trials found that daily vaginal P4 prophylaxis, also beginning at mid-gestation decreased PTB risk in women with a short cervix. Currently, prophylaxis with 17HPC (in women with a history of PTB) or vaginal P4 (in women with a short cervix) are used to prevent PTB. Recent advances in understanding the molecular biology of P4 signaling in uterine cells is revealing novel progestin-based targets for PTB prevention. One possibility is to use selective P4 receptor (PR) modulators (SPRMs) to boost PR anti-inflammatory activity that blocks labor, while simultaneously preventing PR phosphorylation that causes loss of P4/PR anti-inflammatory activity. This may be achieved by SPRMs that induce a specific PR conformation that prevents site-specific serine phosphorylation that inhibits anti-inflammatory activity. Further advances in understanding how P4 promotes uterine quiescence and how its labor blocking actions are withdrawn to trigger parturition will reveal novel therapeutic targets to more effectively prevent PTB.

Exploring Flexibility of Progesterone Receptor Ligand Binding Domain Using Molecular Dynamics

Progesterone receptor (PR), a member of nuclear receptor (NR) superfamily, plays a vital role for female reproductive tissue development, differentiation and maintenance. PR ligand, such as progesterone, induces conformation changes in PR ligand binding domain (LBD), thus mediates subsequent gene regulation cascades. PR LBD may adopt different conformations upon an agonist or an antagonist binding. These different conformations would trigger distinct transcription events. Therefore, the dynamics of PR LBD would be of general interest to biologists for a deep understanding of its structure-function relationship. However, no apo-form (non-ligand bound) of PR LBD model has been proposed either by experiments or computational methods so far. In this study, we explored the structural dynamics of PR LBD using molecular dynamics simulations and advanced sampling tools in both ligand-bound and the apo-forms. Resolved by the simulation study, helix 11, helix 12 and loop 895–908 (the loop between these two helices) are quite flexible in antagonistic conformation. Several residues, such as Arg899 and Glu723, could form salt-bridging interaction between helix 11 and helix 3, and are important for the PR LBD dynamics. And we also propose that helix 12 in apo-form PR LBD, not like other NR LBDs, such as human estrogen receptor α (ERα) LBD, may not adopt a totally extended conformation. With the aid of umbrella sampling and metadynamics simulations, several stable conformations of apo-form PR LBD have been sampled, which may work as critical structural models for further large scale virtual screening study to discover novel PR ligands for therapeutic application.

Crystalline arrays of side chain modified bile acids derivatives. Two novel self-assemblies based on π-π and belly-to-belly interactions

Crystalline derivatives of side chain modified bile acids were efficiently prepared from the naturally occurring steroids by palladium-catalyzed cross coupling reaction as a key step. The solvent-free crystalline bile acids derivatives 2b-e are readily accessed by slow evaporation from selected solvents. A variety of steroidal scaffolds were found and elucidated by SXRD studies. The crystal packing of the title compounds are dominated by hydrogen-bonding interactions established between differently positioned acetyl protecting groups, which in the case of 2b and 2e take advantage of the facial amphiphilicity producing two novel steroidal supramolecular self-assemblies combining π-π and strong facial interactions. Thus, these crystalline arrays of side chain modified bile acids represent promising scaffolds for research and implementation in biomolecular materials or inclusion phenomena.

Progesterone-targeted magnetic resonance imaging probes

Determination of progesterone receptor (PR) status in hormone-dependent diseases is essential in ascertaining disease prognosis and monitoring treatment response. The development of a noninvasive means of monitoring these processes would have significant impact on early detection, cost, repeated measurements, and personalized treatment options. Magnetic resonance imaging (MRI) is widely recognized as a technique that can produce longitudinal studies, and PR-targeted MR probes may address a clinical problem by providing contrast enhancement that reports on PR status without biopsy. Commercially available MR contrast agents are typically delivered via intravenous injection, whereas steroids are administered subcutaneously. Whether the route of delivery is important for tissue accumulation of steroid-modified MRI contrast agents to PR-rich tissues is not known. To address this question, modification of the chemistry linking progesterone with the gadolinium chelate led to MR probes with increased water solubility and lower cellular toxicity and enabled administration through the blood. This attribute came at a cost through lower affinity for PR and decreased ability to cross the cell membrane, and ultimately it did not improve delivery of the PR-targeted MR probe to PR-rich tissues or tumors in vivo. Overall, these studies are important, as they demonstrate that targeted contrast agents require optimization of delivery and receptor binding of the steroid and the gadolinium chelate for optimal translation in vivo.

Recent advances in structure of progestins and their binding to progesterone receptors

The role of progesterone in women's cancers as well as the knowledge of the progesterone receptor (PR) structure has prompted the design of different therapies. The aim of this review is to describe the basic structure of PR agonists and antagonists as well as the recent treatments for illness associated with the progesterone receptor. The rational design for potent and effective drugs for the treatment of female cancer must consider the structural changes of the androgen and progestogen skeleton which are an indicator of their activity as progestins or antiprogestins. The presence of a hydroxyl group at C-17 in the progesterone skeleton brings about a loss of progestational activity whereas acetylation induces a progestational effect. The incorporation of an ethynyl functional group to the testosterone framework results in a loss of androgenic activity with a concomitant enhancement of the progestational effect. On the other hand, an ester function at C-3 of dehydroepiandrosterone skeleton induces partial antagonism to the PR.

Progesterone receptor targeting with radiolabelled steroids: an approach in predicting breast cancer response to therapy

Steroid receptors have demonstrated to be potentially useful biological targets for the diagnosis and therapy follow-up of hormonally responsive cancers. The over-expression of these proteins in human cancer cells as well as their binding characteristics provides a favourable mechanism for the localization of malignant tumours. The need for newer and more selective probes to non-invasively assess steroid receptor expression in hormone-responsive tumours has encouraged the synthesis and the biological evaluation of several steroidal derivatives labelled with positron and gamma emitters. The physiological effects of the steroid hormone progesterone are mediated by the progesterone receptor (PR). Since PR expression is stimulated by the oestrogen receptor (ER), PR status has been considered as a biomarker of ER activity and its value for predicting and monitoring therapeutic efficacy of hormonal therapy has been studied. Imaging of PR-expressing breast cancer patients under hormonal therapy may be advantageous, since the response to therapy can be more accurately predicted after quantification of both ER and PR status. Thus, ligands for PR targeting, although much less explored than ER ligands, have gained some importance lately as potential PET and SPECT tumour imaging agents. In this review, we present a brief survey of explored approaches for progesterone targeting using radiolabelled progestins as potential clinical probes to predict responsiveness to breast cancer therapy. This article is part of a Special Issue entitled "Synthesis and biological testing of steroid derivatives as inhibitors".

A new strategy for selective targeting of progesterone receptor with passive antagonists

Currently available progesterone (P4) receptor (PR) antagonists, such as mifepristone (RU486), lack specificity and display partial agonist properties, leading to potential drawbacks in their clinical use. Recent x-ray crystallographic studies have identified key contacts involved in the binding of agonists and antagonists with PR opening the way for a new rational strategy for inactivating PR. We report here the synthesis and characterization of a novel class of PR antagonists (APRn) designed from such studies. The lead molecule, the homosteroid APR19, displays in vivo endometrial anti-P4 activity. APR19 inhibits P4-induced PR recruitment and transactivation from synthetic and endogenous gene promoters. Importantly, it exhibits high PR selectivity with respect to other steroid hormone receptors and is devoid of any partial agonist activity on PR target gene transcription. Two-hybrid and immunostaining experiments reveal that APR19-bound PR is unable to interact with either steroid receptor coactivators 1 and 2 (SRC1 and SCR2) or nuclear receptor corepressor (NcoR) and silencing mediator of retinoid acid and thyroid hormone receptor (SMRT), in contrast to RU486-PR complexes. APR19 also inhibits agonist-induced phosphorylation of serine 294 regulating PR transcriptional activity and turnover kinetics. In silico docking studies based on the crystal structure of the PR ligand-binding domain show that, in contrast to P4, APR19 does not establish stabilizing hydrogen bonds with the ligand-binding cavity, resulting in an unstable ligand-receptor complex. Altogether, these properties highly distinguish APR19 from RU486 and likely its derivatives, suggesting that it belongs to a new class of pure antiprogestins that inactivate PR by a passive mechanism. These specific PR antagonists open new perspectives for long-term hormonal therapy.