Mechanism Of Ligand Binding Research Articles

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins.

Modulating the function of Intrinsically Disordered Proteins (IDPs) with small molecules is of considerable importance given the crucial roles of IDPs in the pathophysiology of numerous diseases. Reported binding affinities for ligands to diverse IDPs vary broadly, and little is known about the detailed molecular mechanisms that underpin ligand efficacy. Molecular simulations of IDP ligand binding mechanisms can help us understand the mode of action of small molecule inhibitors of IDP function, but it is still unclear how binding energies can be modeled rigorously for such a flexible class of proteins. Here, we compare alchemical absolute binding free energy calculations (ABFE) and Markov-State Modeling (MSM) protocols to model the binding of the small molecule 10058-F4 to a disordered peptide extracted from a segment of the oncoprotein c-Myc. The ABFE results produce binding energy estimates that are sensitive to the choice of reference structure. In contrast, the MSM results produce more reproducible binding energy estimates consistent with weak mM binding affinities and transient intermolecular contacts reported in the literature.

Insight into structural properties of viral G protein-coupled receptors and their role in the viral infection: IUPHAR Review 41.

G protein-coupled receptors (GPCRs) are pivotal in cellular signalling and drug targeting. Herpesviruses encode GPCRs (vGPCRs) to manipulate cellular signalling, thereby regulating various aspects of the virus life cycle, such as viral spreading and immune evasion.vGPCRs mimic host chemokine receptors, often with broader signalling and high constitutive activity. This review focuses on the recent advancements in structural knowledge about vGPCRs, with an emphasis on molecular mechanisms of action and ligand binding. The structures of US27 and US28 from human cytomegalovirus (HCMV) are compared to their closest human homologue, CX3CR1. Contrasting US27 and US28, the homotrimeric UL78 structure (HCMV) reveals more distance to chemokine receptors. Open reading frame 74 (ORF74; Kaposi's sarcoma-associated herpesvirus) is compared to CXCRs, whereas BILF1 (Epstein-Barr virus) is discussed as a putative lipid receptor. Furthermore, the rolesof vGPCRs in latency and lytic replication, reactivation, dissemination and immune evasion arereviewed, together with their potential as drug targets for virus infections and virus-related diseases.

Intense Circular Dichroism and Spin Selectivity in AgBiS2 Nanocrystals by Chiral Ligand Exchange.

Chiral semiconducting nanomaterials offer many potential applications in photodetection, light emission, quantum information, and so on. However, it is difficult to achieve a strong circular dichroism (CD) signal in semiconducting nanocrystals (NCs) due to the complexity of chiral ligand surface engineering and multiple, uncertain mechanisms of chiroptical behavior. Here, a chiral ligand exchange strategy with cysteine on the ternary metal chalcogenide AgBiS2 NCs is developed, and a strong, long-lasting CD signal in the near-UV region is achieved. By carefully optimizing the ligand concentration, the CD peaks are observed at 260 and 320nm, respectively, giving insight into the different ligand binding mechanisms influencing the CD signal of AgBiS2 NCs. Using density-functional theory, a large degree of crystal distortion by the bidentate mode of ligand chelation, and efficient ligand-NC electron transfer, synergistically resulting in the strongest CD signal (g-factor over 10-2) observed in chiral ligand-exchanged semiconductor NCs to date, is demonstrated. To demonstrate the effective chiral properties of these AgBiS2 NCs, a spin-filter device with over 86% efficiency is fabricated. This work represents a considerable leap in the field of chiral semiconductor NCs and points toward their future applications.

Synthesis of metallic nanoparticles using biometabolites: mechanisms and applications.

Bio-metabolites have played a crucial role in the recent green synthesis of nanoparticles, resulting in more versatile, safer, and effective nanoparticles. Various primary and secondary metabolites, such as proteins, carbohydrates, lipids, nucleic acids, enzymes, vitamins, organic acids, alkaloids, flavonoids, and terpenes, have demonstrated strong metal reduction and stabilization properties that can be utilized to synthesize nanomaterials and influence their characters. While physical and chemical methods were previously used to synthesize these nanomaterials, their drawbacks, including high energy consumption, elevated cost, lower yield, and the use of toxic chemicals, have led to a shift towards eco-friendly, rapid, and efficient alternatives. Biomolecules act as reducing agents through deprotonation, nucleophilic reactions, transesterification reactions, ligand binding, and chelation mechanisms, which help sequester metal ions into stable metal nanoparticles (NPs). Engineered NPs have potential applications in various fields due to their optical, electronic, and magnetic properties, offering improved performance compared to bulkier counterparts. NPs can be used in medicine, food and agriculture, chemical catalysts, energy harvesting, electronics, etc. This review provides an overview of the role of primary and secondary metabolites in creating effective nanostructures and their potential applications.

Self-organizing maps of unbiased ligand-target binding pathways and kinetics.

The interpretation of ligand-target interactions at atomistic resolution is central to most efforts in computational drug discovery and optimization. However, the highly dynamic nature of protein targets, as well as possible induced fit effects, makes difficult to sample many interactions effectively with docking studies or even with large-scale molecular dynamics (MD) simulations. We propose a novel application of Self-Organizing Maps (SOMs) to address the sampling and dynamic mapping tasks, particularly in cases involving ligand flexibility and induced fit. The SOM approach offers a data-driven strategy to create a map of the interaction process and pathways based on unbiased MD. Furthermore, we show how the preliminary SOM mapping is complementary to kinetic analysis, with the employment of both network-based approaches and Markov state models. We demonstrate the method by comprehensively mapping a large dataset of 640 μs of unbiased trajectories sampling the recognition process between the phosphorylated YEEI peptide and its high-specificity target lck-SH2. The integration of SOM into unbiased simulation protocols significantly advances our understanding of the ligand binding mechanism. This approach serves as a potent tool for mapping intricate ligand-target interactions with unprecedented detail, thereby enhancing the characterization of kinetic properties crucial to drug design.

Ligand-Coupled Conformational Changes in a Cyclic Nucleotide-Gated Ion Channel Revealed by Time-Resolved Transition Metal Ion FRET.

Allosteric regulation is pivotal for the function of most proteins, especially ion channels like the cyclic nucleotide-binding domain (CNBD) channels. This study examines the allosteric mechanism of ligand binding in the C-terminal region of the prokaryotic CNBD ion channel SthK using steady-state and time-resolved tmFRET. We uncovered significant structural and energetic changes induced by ligand binding with the full-agonist cAMP and the weak partial agonist cGMP. Our approach also highlights the effectiveness of using fluorescence lifetimes to reveal conformational heterogeneity and free energy changes in proteins. These findings deepen our understanding of CNBD channel activation overall and lay the groundwork for a more comprehensive characterization of the effects of mutations and pharmacological agents in these channels.

Functional Characterization Supports Multiple Evolutionary Origins of Pheromone Receptors in Bark Beetles.

Chemical communication using pheromones is thought to have contributed to the diversification and speciation of insects. The species-specific pheromones are detected by specialized pheromone receptors (PRs). Whereas the evolution and function of PRs have been extensively studied in Lepidoptera, only a few PRs have been identified in beetles, which limits our understanding of their evolutionary histories and physiological functions. To shed light on these questions, we aimed to functionally characterize potential PRs in the spruce bark beetle Ips typographus ("Ityp") and explore their evolutionary origins and molecular interactions with ligands. Males of this species release an aggregation pheromone comprising 2-methyl-3-buten-2-ol and (4S)-cis-verbenol, which attracts both sexes to attacked trees. Using two systems for functional characterization, we show that the highly expressed odorant receptor (OR) ItypOR41 responds specifically to (4S)-cis-verbenol, with structurally similar compounds eliciting minor responses. We next targeted the closely related ItypOR40 and ItypOR45. Whereas ItypOR40 was unresponsive, ItypOR45 showed an overlapping response profile with ItypOR41, but a broader tuning. Our phylogenetic analysis shows that these ORs are present in a different OR clade as compared to all other known beetle PRs, suggesting multiple evolutionary origins of PRs in bark beetles. Next, using computational analyses and experimental validation, we reveal two amino acid residues (Gln179 and Trp310) that are important for ligand binding and pheromone specificity of ItypOR41 for (4S)-cis-verbenol, possibly via hydrogen bonding to Gln179. Collectively, our results shed new light on the origins, specificity, and ligand binding mechanisms of PRs in beetles.

Modulation of dynamic DNA G-quadruplex structures in the hTERT promoter region by ligands.

Small molecules can inhibit cellular processes such as replication and transcription by binding to the promoter regions that are prone to form G-quadruplexes. However, since G-quadruplexes exist throughout the human genome, the G-quadruplex binders suffer from specificity issues. To tackle this problem, a G-quadruplex binder (Pyridostatin, or PDS) is conjugated with a ligand (Polyamide, or PA) that can specifically recognize DNA sequences flanking the G-quadruplex forming region. The binding mechanism of this hybrid ligand to the hTERT promoter region (hTERT 5-12) is then elucidated using optical tweezers. During mechanical unfolding processes, different intermediate structures of hTERT 5-12 in presence of PDS, PA, or PA-PDS conjugateare observed. These intermediate structures are consistent with two folding patterns of G-quadruplexes in the hTERT 5-12 fragment. While the duplex DNA binder PA facilitates the folding of a hairpin-G-quadruplex structure, the PDS assists the formation of two tandem G-quadruplexes. Both replication stop assay in vitro and dual luciferase assay in vivo established the effectiveness of the PA-PDS conjugate for hTERT 5-12 targeting. We expect such a ligand dependent folding dynamics will provide guidelines to the development of drugs that not only target hTERT expressions, but also other oncogenes via interactions with specific G-quadruplex structures formed in their promotor regions.

Mathematical Modeling of Tumor Immune Interactions: The Role of Anti-FGFR and Anti-PD-1 in the Combination Therapy.

Bladder cancer poses a significant global health burden with high incidence and recurrence rates. This study addresses the therapeutic challenges in advanced bladder cancer, focusing on the competitive mechanisms of ligand or drug binding to receptors. We developed a refined mathematical model that integrates the dynamics of tumor cells and immune responses, particularly targeting fibroblast growth factor receptor 3 (FGFR3) and immune checkpoint inhibitors (ICIs). This study contributes to understanding combination therapies by elucidating the competitive binding dynamics and quantifying the synergistic effects. The findings highlight the importance of personalized immunotherapeutic strategies, considering factors such as drug dosage, dosing schedules, and patient-specific parameters. Our model further reveals that ligand-independent activated-state receptors are the most essential drivers of tumor proliferation. Moreover, we found that PD-L1 expression rate was more important than PD-1 in driving the dynamic evolution of tumor and immune cells. The proposed mathematical model provides a comprehensive framework for unraveling the complexities of combination therapies in advanced bladder cancer. As research progresses, this multidisciplinary approach contributes valuable insights toward optimizing therapeutic strategies and advancing cancer treatment paradigms.

Leveraging coevolutionary insights and AI-based structural modeling to unravel receptor–peptide ligand-binding mechanisms

Secreted signaling peptides are central regulators of growth, development, and stress responses, but specific steps in the evolution of these peptides and their receptors are not well understood. Also, the molecular mechanisms of peptide-receptor binding are only known for a few examples, primarily owing to the limited availability of protein structural determination capabilities to few laboratories worldwide. Plants have evolved a multitude of secreted signaling peptides and corresponding transmembrane receptors. Stress-responsive SERINE RICH ENDOGENOUS PEPTIDES (SCOOPs) were recently identified. Bioactive SCOOPs are proteolytically processed by subtilases and are perceived by the leucine-rich repeat receptor kinase MALE DISCOVERER 1-INTERACTING RECEPTOR-LIKE KINASE 2 (MIK2) in the model plant Arabidopsis thaliana. How SCOOPs and MIK2 have (co)evolved, and how SCOOPs bind to MIK2 are unknown. Using in silico analysis of 350 plant genomes and subsequent functional testing, we revealed the conservation of MIK2 as SCOOP receptor within the plant order Brassicales. We then leveraged AI-based structural modeling and comparative genomics to identify two conserved putative SCOOP-MIK2 binding pockets across Brassicales MIK2 homologues predicted to interact with the "SxS" motif of otherwise sequence-divergent SCOOPs. Mutagenesis of both predicted binding pockets compromised SCOOP binding to MIK2, SCOOP-induced complex formation between MIK2 and its coreceptor BRASSINOSTEROID INSENSITIVE 1-ASSOCIATED KINASE 1, and SCOOP-induced reactive oxygen species production, thus, confirming our in silico predictions. Collectively, in addition to revealing the elusive SCOOP-MIK2 binding mechanism, our analytic pipeline combining phylogenomics, AI-based structural predictions, and experimental biochemical and physiological validation provides a blueprint for the elucidation of peptide ligand-receptor perception mechanisms.

Structural basis of tethered agonism and G protein coupling of protease-activated receptors

Protease-activated receptors (PARs) are a unique group within the G protein-coupled receptor superfamily, orchestrating cellular responses to extracellular proteases via enzymatic cleavage, which triggers intracellular signaling pathways. Protease-activated receptor 1 (PAR1) is a key member of this family and is recognized as a critical pharmacological target for managing thrombotic disorders. In this study, we present cryo-electron microscopy structures of PAR1 in its activated state, induced by its natural tethered agonist (TA), in complex with two distinct downstream proteins, the Gq and Gi heterotrimers, respectively. The TA peptide is positioned within a surface pocket, prompting PAR1 activation through notable conformational shifts. Contrary to the typical receptor activation that involves the outward movement of transmembrane helix 6 (TM6), PAR1 activation is characterized by the simultaneous downward shift of TM6 and TM7, coupled with the rotation of a group of aromatic residues. This results in the displacement of an intracellular anion, creating space for downstream G protein binding. Our findings delineate the TA recognition pattern and highlight a distinct role of the second extracellular loop in forming β-sheets with TA within the PAR family, a feature not observed in other TA-activated receptors. Moreover, the nuanced differences in the interactions between intracellular loops 2/3 and the Gα subunit of different G proteins are crucial for determining the specificity of G protein coupling. These insights contribute to our understanding of the ligand binding and activation mechanisms of PARs, illuminating the basis for PAR1’s versatility in G protein coupling.

Computational Exploration and Molecular Docking Analysis of Mycobacterium Tuberculosis RpfB Protein

The RpfB protein, a vital component of bacterial cell wall metabolism, is particularly significant in the context of Mycobacterium tuberculosis, the causative agent of tuberculosis. Within this pathogen, RpfB plays a pivotal role in cell growth, division, and virulence, contributing to the survival and persistence of M. tuberculosis within the host. As a member of the Resuscitation-Promoting Factor (Rpf) family, RpfB is crucial for reactivating dormant cells, a process essential for M. tuberculosis to establish latent infections and evade host immune responses. Through its enzymatic activity, RpfB facilitates peptidoglycan remodeling and cell wall synthesis in M. tuberculosis, influencing bacterial physiology and pathogenicity. Elucidating the molecular mechanisms by which RpfB operates in M. tuberculosis not only enhances our understanding of bacterial persistence but also unveils potential targets for novel antimicrobial interventions against tuberculosis, a global health threat of immense significance. The vulnerability of Mycobacterium tuberculosis to various drugs and its persistence has stood as a hurdle in the race against eradication of the pathogenicity of the bacteria. Identification of novel antituberculosis compounds is highly demanding as the available drugs are resistant. The ability of the bacteria to surpass the body's defenses and adapt itself to survive for disease reactivation is contributed by secreted proteins called resuscitating promoting factors (Rpfs). These factors aid in virulence and resuscitation from dormancy of the bacteria. Sequence analysis of RpfB was performed and compounds were first screened for toxicity and high-throughput virtual screening eliminating the toxic compounds. To understand the mechanism of ligand binding and interaction, molecular docking was performed for the compounds passing through the filter resulting with better docking studies predicting the possible binding mode of the inhibitors to the protein.

Decoding the κ Opioid Receptor (KOR): Advancements in Structural Understanding and Implications for Opioid Analgesic Development.

The opioid crisis in the United States is a significant public health issue, with a nearly threefold increase in opioid-related fatalities between 1999 and 2014. In response to this crisis, society has made numerous efforts to mitigate its impact. Recent advancements in understanding the structural intricacies of the κ opioid receptor (KOR) have improved our knowledge of how opioids interact with their receptors, triggering downstream signaling pathways that lead to pain relief. This review concentrates on the KOR, offering crucial structural insights into the binding mechanisms of both agonists and antagonists to the receptor. Through comparative analysis of the atomic details of the binding site, distinct interactions specific to agonists and antagonists have been identified. These insights not only enhance our understanding of ligand binding mechanisms but also shed light on potential pathways for developing new opioid analgesics with an improved risk-benefit profile.

Mitoribosome structure with cofactors and modifications reveals mechanism of ligand binding and interactions with L1 stalk

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNAVal. The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome.

Fatty acid binding to the human transport proteins FABP3, FABP4, and FABP5 from a Ligand’s perspective

Fatty acid binding proteins (FABPs) are a family of amphiphilic transport proteins with high diversity in terms of their amino acid sequences and binding preferences. Beyond their main biological role as cytosolic fatty acid transporters, many aspects regarding their binding mechanism and functional specializations in human cells remain unclear. In this work, the binding properties and thermodynamics of FABP3, FABP4 and FABP5 were analyzed under various physical conditions. For this purpose, the FABPs were loaded with fatty acids bearing fluorescence or spin probes as model ligands, comparing their binding affinities via microscale thermophoresis (MST) and continuous-wave electron paramagnetic resonance (CW EPR) spectroscopy. The CW EPR spectra of non-covalently bound 5- and 16-DOXYL stearic acid (5/16-DSA) deliver in-depth information about the dynamics and chemical environments of ligands inside the binding pockets of the FABPs. EPR spectral simulations allow the construction of binding curves, revealing two different binding states (‘intermediately’ and ‘strongly’ bound). The proportion of bound 5/16-DSA depends strongly on the FABP concentration and the temperature, but with remarkable differences between the three isoforms. Additionally, the more dynamic state (‘intermediately bound’) seems to dominate at body temperature with thermodynamic preference. The ligand binding studies were supplemented by aggregation studies via dynamic light scattering (DLS) and bioinformatic analyses. Beyond the remarkably fine-tuned binding properties exhibited by each FABP, which were discernible with our EPR-centered approach, the results of this work attest the power of simple spectroscopic experiments to provide new insights into the ligand binding mechanisms of proteins in general on a molecular level.

Computation-guided transcription factor biosensor specificity engineering for adipic acid detection

Transcription factor (TF)-based biosensors that connect small-molecule sensing with readouts such as fluorescence have proven to be useful synthetic biology tools for applications in biotechnology. However, the development of specific TF-based biosensors is hindered by the limited repertoire of TFs specific for molecules of interest since current construction methods rely on a limited set of characterized TFs. In this study, we present an approach for engineering the specificity of TFs through a computation-based workflow using molecular docking that enables targeted alteration of TF ligand specificity. Using this method, we engineer the LysR family BenM TF to alter its specificity from its cognate ligand cis,cis-muconic acid to adipic acid through a single amino acid substitution identified by our computational workflow. When implemented in a cell-free system, the engineered biosensor shows higher ligand sensitivity, expanding the potential applications of this circuit. We further investigate ligand binding through molecular dynamics to analyze the substitution, elucidating the impact of modulating a single amino acid position on the mechanism of BenM ligand binding. This study represents the first application of biomolecular modeling methods for altering BenM specificity and for gaining insights into how mutations influence the structural dynamics of BenM. Such methods can potentially be applied to other TFs to alter specificity and analyze the dynamics responsible for these changes, highlighting the applicability of computational tools for informing experiments. In addition, our developed adipic acid biosensor can be applied for the identification and engineering of enzymes to produce adipic acid.

αMI-domain of integrin Mac-1 binds the cytokine pleiotrophin using multiple mechanisms

The integrin Mac-1 (αMβ2, CD11b/CD18, CR3) is an adhesion receptor expressed on macrophages and neutrophils. Mac-1 is also a promiscuous integrin that binds a diverse set of ligands through its αMI-domain. However, the binding mechanism of most ligands remains unclear. We have characterized the interaction of αMI-domain with the cytokine pleiotrophin (PTN), a protein known to bind αMI-domain and induce Mac-1-mediated cell adhesion and migration. Our data show that PTN's N-terminal domain binds a unique site near the N- and C-termini of the αMI-domain using a metal-independent mechanism. However, a stronger interaction is achieved when an acidic amino acid in a zwitterionic motif in PTN's C-terminal domain chelates the divalent cation in the metal ion-dependent adhesion site of active αMI-domain. These results indicate that αMI-domain can bind ligands using multiple mechanisms and that the active αMI-domain has a preference for motifs containing both positively and negatively charged amino acids.

Mechanism of antagonist ligand binding to REV-ERBα.

REV-ERBα, a therapeutically promising nuclear hormone receptor, plays a crucial role in regulating various physiological processes such as the circadian clock, inflammation, and metabolism. However, the availability of chemical probes to investigate the pharmacology of this receptor is limited, with SR8278 being the only identified synthetic antagonist. Moreover, no X-ray crystal structures are currently available that demonstrate the binding of REV-ERBα to antagonist ligands. This lack of structural information impedes the development of targeted therapeutics. To address this issue, we employed Gaussian accelerated molecular dynamics (GaMD) simulations to investigate the binding pathway of SR8278 to REV-ERBα. For comparison, we also used GaMD to observe the ligand binding process of STL1267, for which an X-ray structure is available. GaMD simulations successfully captured the binding of both ligands to the receptor's orthosteric site and predicted the ligand binding pathway and important amino acid residues involved in the antagonist SR8278 binding. This study highlights the effectiveness of GaMD in investigating protein-ligand interactions, particularly in the context of drug recognition for nuclear hormone receptors.

Ligand binding properties of three odorant-binding proteins in striped flea beetle Phyllotreta striolata towards two phthalate esters.

Odorant-binding proteins (OBPs) initiate insect olfactory perception and mediate specific binding and selection of odorants via uncertain binding mechanisms. We characterized the binding characteristics of four OBPs from the striped flea beetle Phyllotreta striolata (SFB), a major cruciferous crop pest. Tissue expression analysis revealed that the two ABPII OBPs (PstrOBP12 and PstrOBP19) were highly expressed mainly in the antenna, whereas the two minus-C OBPs (PstrOBP13 and PstrOBP16) showed a broad expression pattern. Competitive binding assays of cruciferous plant volatiles showed that PstrOBP12, PstrOBP16 and PstrOBP19 had very strong binding capacities for only two phthalate esters (Ki < 20 μM), and PstrOBP13 specifically bound to four aromatic volatiles (Ki < 11 μM). Fluorescence quenching assays displayed that two phthalate esters bound to three PstrOBPs via different quenching mechanisms. PstrOBP12/PstrOBP16-diisobutyl phthalate and PstrOBP19-bis(6-methylheptyl) phthalate followed static quenching, while PstrOBP12/PstrOBP16-bis(6-methylheptyl) phthalate and PstrOBP19-diisobutyl phthalate followed dynamic quenching. Homology modelling and molecular docking displayed that PstrOBP12-diisobutyl phthalate was driven by H-bonding and van der Waals interactions, while PstrOBP16-diisobutyl phthalate and PstrOBP19-bis(6-methylheptyl) phthalate followed hydrophobic interactions. Finally, behavioural activity analysis demonstrated that phthalate esters exhibited different behavioural activities of SFB at different doses, with low doses attracting and high doses repelling. Overall, we thus revealed the different binding properties of the three PstrOBPs to two phthalate esters, which was beneficial in shedding light on the ligand-binding mechanisms of OBPs.

Structural Perspective of NR4A Nuclear Receptor Family and Their Potential Endogenous Ligands.

There are 48 nuclear receptors in the human genome, and many members of this superfamily have been implicated in human diseases. The NR4A nuclear receptor family consisting of three members, NR4A1, NR4A2, and NR4A3 (formerly annotated as Nur77, Nurr1, and NOR1, respectively), are still orphan receptors but exert pathological effects on immune-related and neurological diseases. We previously reported that prostaglandin A1 (PGA1) and prostaglandin A2 (PGA2) are potent activators of NR4A3, which bind directly to the ligand-binding domain (LBD) of the receptor. Recently, the co-crystallographic structures of NR4A2-LBD bound to PGA1 and PGA2 were reported, followed by reports of the neuroprotective effects of these possible endogenous ligands in mouse models of Parkinson's disease. Based on these structures, we modeled the binding structures of the other two members (NR4A1 and NR4A3) with these potential endogenous ligands using a template-based modeling method, and reviewed the similarity and diversity of ligand-binding mechanisms in the nuclear receptor family.