Allosteric Binding Sites On Nuclear Receptors: Focus On Drug Efficacy and Selectivity

NRPreTo: A Machine Learning-Based Nuclear Receptor and Subfamily Prediction Tool

The nuclear receptor (NR) superfamily includes phylogenetically related ligand-activated proteins, which play a key role in various cellular activities. NR proteins are subdivided into seven subfamilies based on their function, mechanism, and nature of the interacting ligand. Developing robust tools to identify NR could give insights into their functional relationships and involvement in disease pathways. Existing NR prediction tools only use a few types of sequence-based features and are tested on relatively similar independent datasets; thus, they may suffer from overfitting when extended to new genera of sequences. To address this problem, we developed Nuclear Receptor Prediction Tool (NRPreTo), a two-level NR prediction tool with a unique training approach where in addition to the sequence-based features used by existing NR prediction tools, six additional feature groups depicting various physiochemical, structural, and evolutionary features of proteins were utilized. The first level of NRPreTo allows for the successful prediction of a query protein as NR or non-NR and further subclassifies the protein into one of the seven NR subfamilies in the second level. We developed Random Forest classifiers to test on benchmark datasets, as well as the entire human protein datasets from RefSeq and Human Protein Reference Database (HPRD). We observed that using additional feature groups improved the performance. We also observed that NRPreTo achieved high performance on the external datasets and predicted 59 novel NRs in the human proteome. The source code of NRPreTo is publicly available at https://github.com/bozdaglab/NRPreTo.

Modeling, molecular dynamics and docking studies of a full-length Echinococcus granulosus 2DBD nuclear receptor

Nuclear receptors are ligand-activated transcription factors capable of regulating the expression of complex gene networks. The family includes seven subfamilies of protein with a wide phylogenetic distribution. A novel subfamily with two DNA-binding domains (2DBDs) has been first reported in Schistosoma mansoni (Platyhelminth, Trematoda). Employing an ab initio protocol and homology modeling methods, the full-length 3D structure of the Eg2DBDα.1 nuclear receptor from Echinococcus granulosus (Platyhelminth, Cestoda) was generated. The model analysis reveals the presence of the conserved three-layered alpha-helical sandwich structure in the ligand binding domain, and a particularly long and flexible hinge region. Molecular dynamics simulations were performed previous to dock a conformational library of fatty acids and retinoic acids. Our results indicate that oleic and linoleic acids are suitable ligands to this receptor. The ligand-protein complex is stabilized mainly by hydrogen bonds and hydrophobic interactions. The fact that 2DBD nuclear receptors have not been identified in vertebrates confers particular interest to these nuclear receptors, not only concerning their structure and function but as targets of new anthelmintic drugs.Communicated by Ramaswamy H. Sarma.

Allosteric Antagonism of the Pregnane X Receptor (PXR): Current-State-of-the-Art and Prediction of Novel Allosteric Sites

The pregnane X receptor (PXR, NR1I2) is a xenobiotic-activated transcription factor with high levels of expression in the liver. It not only plays a key role in drug metabolism and elimination, but also promotes tumor growth, drug resistance, and metabolic diseases. It has been proposed as a therapeutic target for type II diabetes, metabolic syndrome, and inflammatory bowel disease, and PXR antagonists have recently been considered as a therapy for colon cancer. There are currently no PXR antagonists that can be used in a clinical setting. Nevertheless, due to the large and complex ligand-binding pocket (LBP) of the PXR, it is challenging to discover PXR antagonists at the orthosteric site. Alternative ligand binding sites of the PXR have also been proposed and are currently being studied. Recently, the AF-2 allosteric binding site of the PXR has been identified, with several compounds modulating the site discovered. Herein, we aimed to summarize our current knowledge of allosteric modulation of the PXR as well as our attempt to unlock novel allosteric sites. We describe the novel binding function 3 (BF-3) site of PXR, which is also common for other nuclear receptors. In addition, we also mention a novel allosteric site III based on in silico prediction. The identified allosteric sites of the PXR provide new insights into the development of safe and efficient allosteric modulators of the PXR receptor. We therefore propose that novel PXR allosteric sites might be promising targets for treating chronic metabolic diseases and some cancers.

Allosteric binding on nuclear receptors: Insights on screening of non-competitive endocrine-disrupting chemicals

Endocrine-disrupting chemicals (EDCs) can compete with endogenous hormones and bind to the orthosteric site of nuclear receptors (NRs), affecting normal endocrine system function and causing severe symptoms. Recently, a series of pharmaceuticals and personal care products (PPCPs) have been discovered to bind to the allosteric sites of NRs and induce similar effects. However, it remains unclear how diverse EDCs work in this new way. Therefore, we have systematically summarized the allosteric sites and underlying mechanisms based on existing studies, mainly regarding drugs belonging to the PPCP class. Advanced methods, classified as structural biology, biochemistry and computational simulation, together with their advantages and hurdles for allosteric site recognition and mechanism insight have also been described. Furthermore, we have highlighted two available strategies for virtual screening of numerous EDCs, relying on the structural features of allosteric sites and lead compounds, respectively. We aim to provide reliable theoretical and technical support for a broader view of various allosteric interactions between EDCs and NRs, and to drive high-throughput and accurate screening of potential EDCs with non-competitive effects.

A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment

Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 μs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.

Mechanistic Insights into the Allosteric Inhibition of Androgen Receptors by Binding Function 3 Antagonists from an Integrated Molecular Modeling Study

An androgen receptor (AR) is an intensively studied treatment target for castration-resistant prostate cancer that is irresponsive to conventional antiandrogen therapeutics. Binding function 3 (BF3) inhibitors with alternative modes of action have emerged as a promising approach to overcoming antiandrogen resistance. However, how these BF3 inhibitors modulate AR function remains elusive, hindering the development of BF3-targeting agents. Here, we performed an integrated computational study to interrogate the binding mechanism of several known BF3 inhibitors with ARs. Our results show that the inhibitory effect of the BF3 antagonists arises from their allosteric modulation of the activation function (AF2) site, which alters the dynamic coupling between the BF3 and AF2 sites as well as the AF2-coactivator (SRC2-3) interaction. Moreover, the per-residue binding energy analyses reveal the "anchor" role of the linker connecting the phenyl ring and benzimidazole/indole in these BF3 inhibitors. Furthermore, the allosteric driver-interacting residues are found to include both "positive", e.g., Phe673 and Asn833, and "negative" ones, e.g., Phe826, and the differential interactions with these residues provide an explanation why stronger binding does not necessarily result in higher inhibitory activities. Finally, our allosteric communication pathway analyses delineate how the allosteric signals triggered by BF3 binding are propagated to the AF2 pocket through multiple short- and/or long-ranged transmission pathways. Collectively, our combined computational study provides a comprehensive structural mechanism underlying how the selected set of BF3 inhibitors modulate AR function, which will help guide future development of BF3 antagonists.

Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (Mpro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 Mpro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV Mpro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 Mpro inhibitors, as we report it, was not previously conducted.

Computational Assessment of Combination Therapy of Androgen Receptor-Targeting Compounds

The ligand-binding domain of the androgen receptor (AR) is a target for drugs against prostate cancer and offers three distinct binding sites for small molecules. Drugs acting on the orthosteric hormone binding site suffer from resistance mechanisms that can, in the worst case, reverse their therapeutic effect. While many allosteric ligands targeting either the activation function-2 (AF-2) or the binding function-3 (BF-3) have been reported, their potential for simultaneous administration with currently prescribed antiandrogens was disregarded. Here, we report results of 60 μs molecular dynamics simulations to investigate combinations of orthosteric and allosteric AR antagonists. Our results suggest BF-3 inhibitors to be more suitable in combination with classical antiandrogens as opposed to AF-2 inhibitors based on binding free energies and binding modes. As a mechanistic explanation for these observations, we deduced a structural adaptation of helix-12 involved in the formation of the AF-2 site by classical AR antagonists. Additionally, the changes were accompanied by an expansion of the orthosteric binding site. Considering our predictions, the selective combination of AR-targeting compounds may improve the treatment of prostate cancer.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.