Negative allosteric modulators of cannabinoid receptor 2: protein modeling,binding site identification and molecular dynamics simulations in the presence of an orthosteric agonist

Abstract

Selective activation of the cannabinoid receptor subtype 2 (CB2) shows promise for treating pain, inflammation, multiple sclerosis, cancer, ischemic/reperfusion injury and osteoporosis. Target selectivity and off-target side effects are two major limiting factors for orthosteric ligands, and therefore, the search for allosteric modulators (AMs) is a widely used drug discovery approach. To date, only a limited number of negative CB2 AMs have been identified, possessing only micromolar activity at best, and the CB2 receptor’s allosteric site(s) are not well characterized. Herein, we used computational approaches including receptor modeling, site mapping, docking, molecular dynamics (MD) simulations and binding free energy calculations to predict, characterize and validate allosteric sites within the complex of the CB2 receptor with bound orthosteric agonist CP55,940. After docking of known negative CB2 allosteric modulators (NAMs), dihydro-gambogic acid (DHGA) and trans-β-caryophyllene (TBC) (note that TBC also shows agonist activity), at the predicted allosteric sites, the best total complex with CB2, CP55,940 and NAM was embedded into a hydrated lipid bilayer and subjected to a 200 ns MD simulation. The presence of an AM affected the CB2–CP55,940 complex, altering the relative positioning of the toggle switch residues and promoting a strong π–π interaction between Phe117^3.36 and Trp258^6.48. Binding of either TBC or DHGA to a putative allosteric pocket directly adjacent to the orthosteric ligand reduced the binding free energy of CP55,940, which is consistent with the expected effect of a negative AM. The identified allosteric sites present immense scope for the discovery of novel classes of CB2 AMs.     

Ago-Allosteric Modulation and Other Types of Allostery in Dimeric 7TM Receptors

Conventionally, an allosteric modulator is neutral in respect of efficacy and binds to a receptor site distant from the orthosteric site of the endogenous agonist. However, recently compounds being ago-allosteric modulators have been described i.e., compounds acting both as agonists on their own and as enhancers for the endogenous agonists in both increasing agonist potency and providing additive efficacy—superagonism. The additive efficacy can also be observed with agonists, which are neutral or even negative modulators of the potency of the endogenous ligand. Based on the prevailing dimeric concept for 7TM receptors, it is proposed that the ago-allosteric modulators bind in the orthosteric binding site, but–importantly–in the “other” or allosteric protomer of the dimer. Hereby, they can act both as additive co-agonists, and through intermolecular cooperative effects between the protomers, they may influence the potency of the endogenous agonist. It is of interest that at least some endogenous agonists can only occupy one protomer of a dimeric 7TM receptor complex at a time and thereby they leave the orthosteric binding site in the allosteric protomer free, potentially for binding of exogenous, allosteric modulators. If the allosteric modulator is an agonist, it is an ago-allosteric modulator; if it is neutral, it is a classical enhancer. Molecular mapping in hetero-dimeric class-C receptors, where the endogenous agonist clearly binds only in one protomer, supports the notion that allosteric modulators can act through binding in the “other” protomer. It is suggested that for the in vivo, clinical setting a positive ago-allosteric modulator should be the preferred agonist drug.     

Integrative computational protocol for the discovery of inhibitors of the <Emphasis Type="Italic">Helicobacter pylori</Emphasis> nickel response regulator (NikR)

In order to identify novel inhibitors of the Helicobacter pylori nickel response regulator (HpNikR) an integrative protocol was performed for half a million compounds retrieved from the ZINC database. We firstly implement a structure-based virtual screening to build a library of potential inhibitors against the HpNikR using a docking analysis (AutoDock Vina). The library was then used to perform a hierarchical clustering of docking poses, based on protein-contact footprints calculation from the multiple conformations given by the AutoDock Vina software, and the drug-protein interaction analyses to identify and remove potential promiscuous compounds likely interacting with human proteins, hence causing drug side effects. 250 drug-like compounds were finally proposed as non-promicuous potential inhibitors for HpNikR. These compounds target the DNA-binding sites of HpNikR so that HpNikR-compound binding could be able to mimic key interactions in the DNA-protein recognition process. HpNikR inhibitors with promising potential against H. pylori could also act against other human bacterial pathogens due to the conservation of targeting motif of NikR involved in DNA-protein interaction.     

Allosteric small molecules unveil a role of an extracellular E2/transmembrane helix 7 junction for G protein-coupled receptor activation

G protein-coupled receptors represent the largest superfamily of cell membrane-spanning receptors. We used allosteric small molecules as a novel approach to better understand conformational changes underlying the inactive-to-active switch in native receptors. Allosteric molecules bind outside the orthosteric area for the endogenous receptor activator. The human muscarinic M(2) acetylcholine receptor is prototypal for the study of allosteric interactions. We measured receptor-mediated G protein activation, applied a series of structurally diverse muscarinic allosteric agents, and analyzed their cooperative effects with orthosteric receptor agonists. A strong negative cooperativity of receptor binding was observed with acetylcholine and other full agonists, whereas a pronounced negative cooperativity of receptor activation was observed with the partial agonist pilocarpine. Applying a newly synthesized allosteric tool, point mutated receptors, radioligand binding, and a three-dimensional receptor model, we found that the deviating allosteric/orthosteric interactions are mediated through the core region of the allosteric site. A key epitope is M(2)Trp(422) in position 7.35 that is located at the extracellular top of transmembrane helix 7 and that contacts, in the inactive receptor, the extracellular loop E2. Trp 7.35 is critically involved in the divergent allosteric/orthosteric cooperativities with acetylcholine and pilocarpine, respectively. In the absence of allosteric agents, Trp 7.35 is essential for receptor binding of the full agonist and for receptor activation by the partial agonist. This study provides first evidence for a role of an allosteric E2/transmembrane helix 7 contact region for muscarinic receptor activation by orthosteric agonists.     

Investigation of the influence of external factors on the conformational dynamics of rhodopsin-like receptors by means of molecular dynamics simulation

Abstract

The study reports about the influence of binding of orthosteric ligands on the conformational dynamics of β-2-adrenoreceptor. Using molecular dynamics (MD) simulation, we found that there was a little fraction of active states of the receptor in its apo (ligand-free) ensemble. Analysis of MD trajectories indicated that such spontaneous activation of the receptor is accompanied by the motion in intracellular part of its alpha-helices. Thus, receptor’s constitutive activity directly results from its conformational dynamics. On the other hand, the binding of a full agonist resulted in a significant shift of the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized the receptor in its inactive state. It is likely that the binding of inverse agonists might be a universal way of constitutive activity inhibition in vivo. Our results indicate that ligand binding redistribute pre-existing conformational degrees of freedom (in accordance to the Monod–Wyman–Changeux Model) of the receptor rather than cause induced fit in it. Therefore, the ensemble of biologically relevant receptor conformations is encoded in its spatial structure, and individual conformations from that ensemble might be used by the cell in conformity with the physiological behavior.     

Influence of orthosteric ligand binding on the conformational dynamics of the β-2-adrenergic receptor via essential dynamics sampling simulation

The influence of the binding of orthosteric ligands on the conformational dynamics of the β-2-adrenoreceptor was identified using the molecular dynamics method. It was found that there was a small fraction of active states of the receptor in its apo (ligand free) ensemble. An analysis of the MD trajectories indicated that this spontaneous activation of the receptor was accompanied by the motion of its VI helix. Thus, the receptor’s constitutive activity is a direct result of its conformational dynamics. On the other hand, the binding of the full agonist resulted in a significant shift in the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized receptor in its inactive state. It is likely that the binding of the inverse agonists might be a universal method of the constitutive activity inhibition. Our results indicate that ligand binding redistributes preexisting conformational degrees of freedom (in accordance to the Monod-Wyman-Changeux Model), rather than causes an induced fit. Therefore, the ensemble of the biologically relevant receptors conformations has been encoded in its spatial structure and individual conformations from that ensemble and might be used by the cell according to the physiological behavior.     

Virtual Screening for HIV Protease Inhibitors: A Comparison of AutoDock 4 and Vina

Background

The AutoDock family of software has been widely used in protein-ligand docking research. This study compares AutoDock 4 and AutoDock Vina in the context of virtual screening by using these programs to select compounds active against HIV protease.

Methodology/Principal Findings

Both programs were used to rank the members of two chemical libraries, each containing experimentally verified binders to HIV protease. In the case of the NCI Diversity Set II, both AutoDock 4 and Vina were able to select active compounds significantly better than random (AUC = 0.69 and 0.68, respectively; p<0.001). The binding energy predictions were highly correlated in this case, with r = 0.63 and ι = 0.82. For a set of larger, more flexible compounds from the Directory of Universal Decoys, the binding energy predictions were not correlated, and only Vina was able to rank compounds significantly better than random.

Conclusions/Significance

In ranking smaller molecules with few rotatable bonds, AutoDock 4 and Vina were equally capable, though both exhibited a size-related bias in scoring. However, as Vina executes more quickly and is able to more accurately rank larger molecules, researchers should look to it first when undertaking a virtual screen.     

面向虚拟筛选的GPU加速的分子对接方法

虚拟筛选是在计算机上对化合物分子进行模拟预筛选,找出容易和药物靶标结合的小分子(配体),从而降低实际实验测试次数,提高药物先导化合物的发现效率。常用的分子对接软件可以用于基于结构的虚拟筛选,寻找配体与靶标的最佳的作用模式和结合构象,并通过打分函数来筛选出潜在的配体。现有的对接软件如AutoDock Vina等在分子对接过程中需要耗费大量时间和计算资源,特别是面对大规模分子对接时,过长的筛选时间不能满足应用需求,因此,本文在最高效的QVina2对接软件基础上,提出一种基于GPU的QVina 2并行化方法QVina2-GPU,利用GPU硬件高度并行体系加速分子对接。具体包括增加初始化分子构象数量,以扩展蒙特卡罗的迭代局部搜索中线程的并行规模,增加蒙特卡罗的迭代搜索的广度以减少每次蒙特卡罗迭代搜索深度,并利用Wolfe-Powell准则改进局部搜索算法,提高了对接精度,进一步减少蒙特卡罗迭代搜索深度,最后,在NVIDIA Geforce RTX 3090平台上在公开的配体数据库上验证了QVina2-GPU的性能,实验表明在保证分子对接精度的基础上,我们提出的QVina2-GPU对Qvina2的平均加速比达到5.18倍,最大加速比达到12.28倍。     

Ago-allosteric modulation and other types of allostery in dimeric 7TM receptors

Conventionally, an allosteric modulator is neutral in respect of efficacy and binds to a receptor site distant from the orthosteric site of the endogenous agonist. However, recently compounds being ago-allosteric modulators have been described i.e., compounds acting both as agonists on their own and as enhancers for the endogenous agonists in both increasing agonist potency and providing additive efficacy-superagonism. The additive efficacy can also be observed with agonists, which are neutral or even negative modulators of the potency of the endogenous ligand. Based on the prevailing dimeric concept for 7TM receptors, it is proposed that the ago-allosteric modulators bind in the orthosteric binding site, but-importantly-in the "other" or allosteric protomer of the dimer. Hereby, they can act both as additive co-agonists, and through intermolecular cooperative effects between the protomers, they may influence the potency of the endogenous agonist. It is of interest that at least some endogenous agonists can only occupy one protomer of a dimeric 7TM receptor complex at a time and thereby they leave the orthosteric binding site in the allosteric protomer free, potentially for binding of exogenous, allosteric modulators. If the allosteric modulator is an agonist, it is an ago-allosteric modulator; if it is neutral, it is a classical enhancer. Molecular mapping in hetero-dimeric class-C receptors, where the endogenous agonist clearly binds only in one protomer, supports the notion that allosteric modulators can act through binding in the "other" protomer. It is suggested that for the in vivo, clinical setting a positive ago-allosteric modulator should be the preferred agonist drug.     

Contrasting Properties of α7-Selective Orthosteric and Allosteric Agonists Examined on Native Nicotinic Acetylcholine Receptors

Subtype-selective ligands are important tools for the pharmacological characterisation of neurotransmitter receptors. This is particularly the case for nicotinic acetylcholine receptors (nAChRs), given the heterogeneity of their subunit composition. In addition to agonists and antagonists that interact with the extracellular orthosteric nAChR binding site, a series of nAChR allosteric modulators have been identified that interact with a distinct transmembrane site. Here we report studies conducted with three pharmacologically distinct nicotinic ligands, an orthosteric agonist (compound B), a positive allosteric modulator (TQS) and an allosteric agonist (4BP-TQS). The primary focus of the work described in this study is to examine the suitability of these compounds for the characterisation of native neuronal receptors (both rat and human). However, initial experiments were conducted on recombinant nAChRs demonstrating the selectivity of these three compounds for α7 nAChRs. In patch-clamp recordings on rat primary hippocampal neurons we found that all these compounds displayed pharmacological properties that mimicked closely those observed on recombinant α7 nAChRs. However, it was not possible to detect functional responses with compound B, an orthosteric agonist, using a fluorescent intracellular calcium assay on either rat hippocampal neurons or with human induced pluripotent stem cell-derived neurons (iCell neurons). This is, presumably, due to the rapid desensitisation of α7 nAChR that is induced by orthosteric agonists. In contrast, clear agonist-evoked responses were observed in fluorescence-based assays with the non-desensitising allosteric agonist 4BP-TQS and also when compound B was co-applied with the non-desensitising positive allosteric modulator TQS. In summary, we have demonstrated the suitability of subtype-selective orthosteric and allosteric ligands for the pharmacological identification and characterisation of native nAChRs and the usefulness of ligands that minimise receptor desensitisation for the characterisation of α7 nAChRs in fluorescence-based assays.     

Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.     

Benzodiazepine ligands can act as allosteric modulators of the Type 1 cholecystokinin receptor

The cholecystokinin (CCK₁) receptor is a G protein-coupled receptor important for nutrient homeostasis. The molecular basis of CCK-receptor binding has been debated, with one prominent model suggesting occupation of the same region of the intramembranous helical bundle as benzodiazepines. Here, we used a specific assay of allosteric ligand interaction to probe the mode of binding of devazepide, a prototypic benzodiazepine ligand. Devazepide elicited marked slowing of dissociation of pre-bound CCK, only possible through binding to a topographically distinct allosteric site. This effect was disrupted by chemical modification of a cysteine in the benzodiazepine-binding pocket. Application of an allosteric model to the equilibrium interaction between a series of benzodiazepine ligands and CCK yielded quantitative estimates of each modulator’s affinity for the allosteric site, as well as the degree of negative cooperativity for the interaction between occupied orthosteric and allosteric sites. The allosteric nature of benzodiazepine binding to the CCK₁ receptor provides new opportunities for small molecule drug development.     

Regulation of signal transduction at M2 muscarinic receptor

Muscarinic acetylcholine receptors mediate transmission of an extracellular signal represented by released acetylcholine to neuronal or effector cells. There are five subtypes of closely homologous muscarinic receptors which are coupled by means of heterotrimeric G-proteins to a variety of signaling pathways resulting in a multitude of target cell effects. Endogenous agonist acetylcholine does not discriminate among individual subtypes and due to the close homology of the orthosteric binding site the same holds true for most of exogenous agonists. In addition to the classical binding site muscarinic receptors have one or more allosteric binding sites at extracellular domains. Binding of allosteric modulators induces conformational changes in the receptor that result in subtype-specific changes in orthosteric binding site affinity for both muscarinic agonists and antagonists. This overview summarizes our recent experimental effort in investigating certain aspects of M2 muscarinic receptor functioning concerning i) the molecular determinants that contribute to the binding of allosteric modulators, ii) G-protein coupling specificity and subsequent cellular responses and iii) possible functional assays that exploit the unique properties of allosteric modulators for characterization of muscarinic receptor subtypes in intact tissue. A detailed knowledge of allosteric properties of muscarinic receptors is required to permit drug design that will modulate signal transmission strength of specific muscarinic receptor subtypes. Furthermore, allosteric modulation of signal transmission strength is determined by cooperativity rather than concentration of allosteric modulator and thus reduces the danger of overdose.     

The Activity of GAT107, an Allosteric Activator and Positive Modulator of α7 Nicotinic Acetylcholine Receptors (nAChR), Is Regulated by Aromatic Amino Acids That Span the Subunit Interface

GAT107, the (+)-enantiomer of racemic 4-(4-bromophenyl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-sulfonamide, is a strong positive allosteric modulator (PAM) of α7 nicotinic acetylcholine receptor (nAChR) activation by orthosteric agonists with intrinsic allosteric agonist activities. The direct activation produced by GAT107 in electrophysiological studies is observed only as long as GAT107 is freely diffusible in solution, although the potentiating activity primed by GAT107 can persist for over 30 min after drug washout. Direct activation is sensitive to α7 nAChR antagonist methyllycaconitine, although the primed potentiation is not. The data are consistent with GAT107 activity arising from two different sites. We show that the coupling between PAMs and the binding of orthosteric ligands requires tryptophan 55 (Trp-55), which is located at the subunit interface on the complementary surface of the orthosteric binding site. Mutations of Trp-55 increase the direct activation produced by GAT107 and reduce or prevent the synergy between allosteric and orthosteric binding sites, so that these mutants can also be directly activated by other PAMs such as PNU-120596 and TQS, which do not activate wild-type α7 in the absence of orthosteric agonists. We identify Tyr-93 as an essential element for orthosteric activation, because Y93C mutants are insensitive to orthosteric agonists but respond to GAT107. Our data show that both orthosteric and allosteric activation of α7 nAChR require cooperative activity at the interface between the subunits in the extracellular domain. These cooperative effects rely on key aromatic residues, and although mutations of Trp-55 reduce the restraints placed on the requirement for orthosteric agonists, Tyr-93 can conduct both orthosteric activation and desensitization among the subunits.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

mGluR5: Exploration of Orthosteric and Allosteric Ligand Binding Pockets and Their Applications to Drug Discovery

Since its discovery in 1992, mGluR5 has attracted significant attention and been linked to several neurological and psychiatric diseases. Ligand development was initially focused on the orthosteric binding pocket, but lack of subtype selective ligands changed the focus to the transmembrane allosteric binding pocket. This strategy has resulted in several drug candidates in clinical testing. In the present article we explore the orthosteric and allosteric binding pockets in terms of structure and ligand recognition across the mGluR subtypes and groups, and discuss the clinical potential of ligands targeting these pockets. We have performed binding mode analyses of non- and group-selective orthosteric ligands based on molecular docking in mGluR crystal structures and models. For the analysis of the allosteric binding pocket we have combined data from all mGluR5-mutagenesis studies, collectively reporting five negative allosteric modulators and 47 unique mutations, and compared it to the closest related homolog, mGluR1.     

Homology modeling,docking, and molecular dynamics simulation of the receptor GALR2 and its interactions with galanin and a positive allosteric modulator

Galanin receptor type 2 (GALR2) is a class A G-protein-coupled receptor (GPCR), and it has been reported that orthosteric ligands and positive allosteric modulators (PAMs) of GALR2 could potentially be used to treat epilepsy. So far, the X-ray structure of this receptor has not been resolved, and knowledge of the 3D structure of GALR2 may prove informative in attempts to design novel ligands and to explore the mechanism for the allosteric modulation of this receptor. In this study, homology modeling was used to obtain several GALR2 models using known templates. ProSA-web Z-scores and Ramachandran plots as well as pre-screening against a test dataset of known compounds were all utilized to select the best model of GALR2. Molecular dockings of galanin (a peptide) and a nonpeptide ligand were carried out to choose the (GALR2 model)–galanin complex that showed the closest agreement with the corresponding experimental data. Finally, a 50-ns MD simulation was performed to study the interactions between the GALR2 model and the synthetic and endogenous ligands. The results from docking and MD simulation showed that, besides the reported residues, Tyr160^4.60, Ile105^3.32, Ala274^7.35, and Tyr163^ECL2 also appear to play important roles in the binding of galanin. The potential allosteric binding pockets in the GALR2 model were then investigated via MD simulation. The results indicated that the mechanism for the allosteric modulation caused by PAMs is the binding of the PAM at pocket III, which is formed by galanin, ECL2, TM2, TM3, and ECL1; this results in the disruption of the Na⁺-binding site and/or the Na⁺ ion pathway, leading to GALR2 agonism.