Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A₁AR, A_2AAR, A_2BAR and A₃AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A₁AR and A₃AR preferentially couple to the G_i/o proteins, while the A_2AAR and A_2BAR prefer coupling to the G_s proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.     

Binding of Androgen- and Estrogen-Like Flavonoids to Their Cognate (Non)Nuclear Receptors: A Comparison by Computational Prediction

Flavonoids are plant bioactives that are recognized as hormone-like polyphenols because of their similarity to the endogenous sex steroids 17β-estradiol and testosterone, and to their estrogen- and androgen-like activity. Most efforts to verify flavonoid binding to nuclear receptors (NRs) and explain their action have been focused on ERα, while less attention has been paid to other nuclear and non-nuclear membrane androgen and estrogen receptors. Here, we investigate six flavonoids (apigenin, genistein, luteolin, naringenin, quercetin, and resveratrol) that are widely present in fruits and vegetables, and often used as replacement therapy in menopause. We performed comparative computational docking simulations to predict their capability of binding nuclear receptors ERα, ERβ, ERRβ, ERRγ, androgen receptor (AR), and its variant AR^T877A and membrane receptors for androgens, i.e., ZIP9, GPRC6A, OXER1, TRPM8, and estrogens, i.e., G Protein-Coupled Estrogen Receptor (GPER). In agreement with data reported in literature, our results suggest that these flavonoids show a relevant degree of complementarity with both estrogen and androgen NR binding sites, likely triggering genomic-mediated effects. It is noteworthy that reliable protein–ligand complexes and estimated interaction energies were also obtained for some suggested estrogen and androgen membrane receptors, indicating that flavonoids could also exert non-genomic actions. Further investigations are needed to clarify flavonoid multiple genomic and non-genomic effects. Caution in their administration could be necessary, until the safe assumption of these natural molecules that are largely present in food is assured.     

Molecular Modeling of Histamine Receptors—Recent Advances in Drug Discovery

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer’s disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014–2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.     

The Structural Determinants for α1-Adrenergic/Serotonin Receptors Activity among Phenylpiperazine-Hydantoin Derivatives

Several studies confirmed the reciprocal interactions between adrenergic and serotoninergic systems and the influence of these phenomena on the pathogenesis of anxiety. Hence, searching for chemical agents with a multifunctional pharmacodynamic profile may bring highly effective therapy for CNS disorders. This study presents a deep structural insight into the hydantoin-arylpiperazine group and their serotonin/α-adrenergic activity. The newly synthesized compounds were tested in the radioligand binding assay and the intrinsic activity was evaluated for the selected derivatives. The computer-aided SAR analysis enabled us to answer questions about the influence of particular structural fragments on selective vs. multifunctional activity. As a result of the performed investigations, there were two leading structures: (a) compound 12 with multifunctional adrenergic-serotonin activity, which is a promising candidate to be an effective anxiolytic agent; (b) compound 14 with high α_1A/α_1D affinity and selectivity towards α_1B, which is recommended due to the elimination of probable cardiotoxic effect. The structural conclusions of this work provide significant support for future lead optimization in order to achieve the desired pharmacodynamic profile in searching for new CNS-modulating agents.     

Stalis: A Computational Method for Template-Based Ab Initio Ligand Design

Proteins interact with small molecules through specific molecular recognition, which is central to essential biological functions in living systems. Therefore, understanding such interactions is crucial for basic sciences and drug discovery. Here, we present S tructure t emplate-based a b initio li gand design s olution (Stalis), a knowledge-based approach that uses structure templates from the Protein Data Bank libraries of whole ligands and their fragments and generates a set of molecules (virtual ligands) whose structures represent the pocket shape and chemical features of a given target binding site. Our benchmark performance evaluation shows that ligand structure-based virtual screening using virtual ligands from Stalis outperforms a receptor structure-based virtual screening using AutoDock Vina, demonstrating reliable overall screening performance applicable to computational high-throughput screening. However, virtual ligands from Stalis are worse in recognizing active compounds at the small fraction of a rank-ordered list of screened library compounds than crystal ligands, due to the low resolution of the virtual ligand structures. In conclusion, Stalis can facilitate drug discovery research by designing virtual ligands that can be used for fast ligand structure-based virtual screening. Moreover, Stalis provides actual three-dimensional ligand structures that likely bind to a target protein, enabling to gain structural insight into potential ligands. Stalis can be an efficient computational platform for high-throughput ligand design for fundamental biological study and drug discovery research at the proteomic level. © 2019 Wiley Periodicals, Inc.     

High‐Throughput Kinetic Analysis for Target‐Directed Covalent Ligand Discovery

Cysteine‐reactive small molecules are used as chemical probes of biological systems and as medicines. Identifying high‐quality covalent ligands requires comprehensive kinetic analysis to distinguish selective binders from pan‐reactive compounds. Quantitative irreversible tethering (qIT), a general method for screening cysteine‐reactive small molecules based upon the maximization of kinetic selectivity, is described. This method was applied prospectively to discover covalent fragments that target the clinically important cell cycle regulator Cdk2. Crystal structures of the inhibitor complexes validate the approach and guide further optimization. The power of this technique is highlighted by the identification of a Cdk2‐selective allosteric (type IV) kinase inhibitor whose novel mode‐of‐action could be exploited therapeutically.     

分子钳人工受体研究进展

分子识别是生物体系的基本特征, 并在生命活动中起中心作用. 利用合成的人工受体与适当底物间的分子识别以建立化学模型或化学仿生体系对生命过程中的分子识别现象进行模拟研究是生物有机化学和超分子化学前沿富于挑战的课题之一. 按照不同的隔离基, 综述了分子钳人工受体的研究进展.     

Computational Prediction and Experimental Validation of the Unique Molecular Mode of Action of Scoulerine

Scoulerine is a natural compound that is known to bind to tubulin and has anti-mitotic properties demonstrated in various cancer cells. Its molecular mode of action has not been precisely known. In this work, we perform computational prediction and experimental validation of the mode of action of scoulerine. Based on the existing data in the Protein Data Bank (PDB) and using homology modeling, we create human tubulin structures corresponding to both free tubulin dimers and tubulin in a microtubule. We then perform docking of the optimized structure of scoulerine and find the highest affinity binding sites located in both the free tubulin and in a microtubule. We conclude that binding in the vicinity of the colchicine binding site and near the laulimalide binding site are the most likely locations for scoulerine interacting with tubulin. Thermophoresis assays using scoulerine and tubulin in both free and polymerized form confirm these computational predictions. We conclude that scoulerine exhibits a unique property of a dual mode of action with both microtubule stabilization and tubulin polymerization inhibition, both of which have similar affinity values.     

A Critical Evaluation of Terpenoid Signaling at Cannabinoid CB1 Receptors in a Neuronal Model

In addition to phytocannabinoids, cannabis contains terpenoids that are claimed to have a myriad of effects on the body. We tested a panel of five common cannabis terpenoids, myrcene, linalool, limonene, α-pinene and nerolidol, in two neuronal models, autaptic hippocampal neurons and dorsal root ganglion (DRG) neurons. Autaptic neurons express a form of cannabinoid CB1 receptor-dependent retrograde plasticity while DRGs express a variety of transient receptor potential (TRP) channels. Most terpenoids had little or no effect on neuronal cannabinoid signaling. The exception was nerolidol, which inhibited endocannabinoid signaling. Notably, this is not via inhibition of CB1 receptors but by inhibiting some aspect of 2-arachidonoylglycerol (2-AG) production/delivery; the mechanism does not involve reducing the activity of the 2-AG-synthesizing diacylglycerol lipases (DAGLs). Nerolidol was also the only terpenoid that activated a sustained calcium response in a small (7%) subpopulation of DRGs. In summary, we found that only one of five terpenoids tested had notable effects on cannabinoid signaling in two neuronal models. Our results suggest that a few terpenoids may indeed interact with some components of the cannabinoid signaling system and may therefore offer interesting insights upon further study.     

Binding affinity of substituted ureido‐benzenesulfonamide ligands to the carbonic anhydrase receptor: A theoretical study of enzyme inhibition

The binding properties of a series of benzenesulfonamide inhibitors (4‐substituted‐ureido‐benzenesulfonamides, UBSAs) of human carbonic anhydrase II (hCA II) enzyme with active site residues have been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. To account for the important docking interactions between the UBSAs ligand and hCA II enzyme, a molecular docking program AutoDock Vina is used. The molecular docking results obtained by AutoDock Vina revealed that the docked conformer has root mean square deviation value less than 1.50 Å compared to X‐ray crystal structures. The inhibitory activity of UBSA ligands against hCA II is found to be in good agreement with the experimental results. The thermodynamic parameters for inhibitor binding show that hydrogen bonding, hydrophilic, and hydrophobic interactions play a major role in explaining the diverse inhibitory range of these derivatives. Additionally, natural bond orbital analysis is performed to characterize the ligand–metal charge transfer stability. The insights gained from this study have great potential to design new hCA‐II inhibitor, 4‐[3‐(1‐p‐Tolyl‐4‐trifluoromethyl‐1H‐pyrazol‐3‐yl)‐ureido]‐benzenesulfonamide, which belongs to the family of UBSA inhibitors and shows similar type of inhibitor potency with hCA II. This work also reveals that a QM/MM model and molecular docking method are computationally feasible and accurate for studying substrate–protein inhibition. © 2013 Wiley Periodicals, Inc.     

甘脲类分子胶囊受体合成方法的研究进展

甘脲类分子胶囊是一类独特的由非共价键弱相互作用力组装而成的人工受体, 在分子识别、自组装、分子微反应器、智能材料等方面得到了广泛而深入的研究. 系统综述了网球型(Tennis ball)、垒球型(Softball)、足球型(Football)及圆盘型(Jelly Doughnut)等分子胶囊的合成方法.     

新型希夫碱型钳形人工受体的微波合成及其对中性分子的识别性能

以对苯二甲酸为起始原料,设计并合成了8个新型的芳环钳形人工受体(3a～3h),其结构经1H NMR,IR,ESI-MS和元素分析表征。利用UV-Vis滴定法初步考察了3a～3h对中性分子(邻苯二胺,间苯二胺,对苯二胺和二苯甲酮)的识别性能,结果表明3a～3h对其具有良好的识别性能。     

An Inside Job: Molecular Determinants for Postsynaptic Localization of Nicotinic Acetylcholine Receptors

Nicotinic acetylcholine receptors (nAChRs) mediate fast synaptic transmission at neuromuscular and autonomic ganglionic synapses in the peripheral nervous system. The postsynaptic localization of muscle ((α1)₂β1γδ) and neuronal ((α3β4)₂β4) nicotinic receptors at these synapses is mediated by interactions between the nAChR intracellular domains and cytoplasmic scaffolding proteins. Recent high resolution structures and functional studies provide new insights into the molecular determinants that mediate these interactions. Surprisingly, they reveal that the muscle nAChR binds 1–3 rapsyn scaffolding molecules, which dimerize and thereby form an interconnected lattice between receptors. Moreover, rapsyn binds two distinct sites on the nAChR subunit cytoplasmic loops; the MA-helix on one or more subunits and a motif specific to the β subunit. Binding at the latter site is regulated by agrin-induced phosphorylation of βY390, and increases the stoichiometry of rapsyn/AChR complexes. Similarly, the neuronal nAChR may be localized at ganglionic synapses by phosphorylation-dependent interactions with 14-3-3 adaptor proteins which bind specific motifs in each of the α3 subunit cytoplasmic loops. Thus, postsynaptic localization of nAChRs is mediated by regulated interactions with multiple scaffolding molecules, and the stoichiometry of these complexes likely helps regulate the number, density, and stability of receptors at the synapse.     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.     

Macrocyclic Receptors for Identification and Selective Binding of Substrates of Different Nature

Molecular recognition of host/guest molecules represents the basis of many biological processes and phenomena. Enzymatic catalysis and inhibition, immunological response, reproduction of genetic information, biological regulatory functions, the effects of drugs, and ion transfer—all these processes include the stage of structure recognition during complexation. The goal of this review is to solicit and publish the latest advances in the design and sensing and binding abilities of porphyrin-based heterotopic receptors with well-defined geometries, the recognition ability of which is realized due to ionic, H-bridge, charge transfer, hydrophobic, and hydrophilic interactions. The dissection of the considered low-energy processes at the molecular scale expands our capabilities in the development of effective systems for controlled recognition, selective delivery, and prolonged release of substrates of different natures (including drugs) to their sites of functioning.     

The Interplay of Cholesterol and Ligand Binding in hTSPO from Classical Molecular Dynamics Simulations

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol–binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N–terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C–terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150–X–Y152–X(3)–R156 and R135–X(2)–Y138–X(2)–L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC–resembling motif that we observed in TM I (residues L17–X(2)–F20–X(3)–R24) and to CRAC in TM V. We expect that the CRAC–like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.     

Identification of Multiple Druggable Secondary Sites by Fragment Screening against DC-SIGN

DC-SIGN is a cell-surface receptor for several pathogenic threats, such as HIV, Ebola virus, or Mycobacterium tuberculosis. Multiple attempts to develop inhibitors of the underlying carbohydrate–protein interactions have been undertaken in the past fifteen years. Still, drug-like DC-SIGN ligands are sparse, which is most likely due to its hydrophilic, solvent-exposed carbohydrate-binding site. Herein, we report on a parallel fragment screening against DC-SIGN applying SPR and a reporter displacement assay, which complements previous screenings using ¹⁹F NMR spectroscopy and chemical fragment microarrays. Hit validation by SPR and ¹H–¹⁵N HSQC NMR spectroscopy revealed that although no fragment bound in the primary carbohydrate site, five secondary sites are available to harbor drug-like molecules. Building on key interactions of the reported fragment hits, these pockets will be targeted in future approaches to accelerate the development of DC-SIGN inhibitors.