Multiple Binding Sites Contribute to the Mechanism of Mixed Agonistic and Positive Allosteric Modulators of the Cannabinoid CB1 Receptor

The cannabinoid CB1 receptor (CB1R) is an abundant metabotropic G‐protein‐coupled receptor that has been difficult to address therapeutically because of CNS side effects exerted by orthosteric drug candidates. Recent efforts have focused on developing allosteric modulators that target CB1R. Compounds from the recently discovered class of mixed agonistic and positive allosteric modulators (Ago‐PAMs) based on 2‐phenylindoles have shown promising functional and binding properties as CB1R ligands. Here, we identify binding modes of both the CP 55,940 agonist and GAT228, a 2‐phenylindole allosteric modulator, by using our metadynamics simulation protocol, and quantify their affinity and cooperativity by atomistic simulations. We demonstrate the involvement of multiple adjunct binding sites in the Ago‐PAM characteristics of the 2‐phenylindole modulators and explain their ability to compete with orthosteric agonists at higher concentrations. We validate these results experimentally by showing the contribution of multiple sites on the allosteric binding of ZCZ011, another homologous member of the class, together with the orthosteric agonist.     

Universality of the Sodium Ion Binding Mechanism in Class A G‐Protein‐Coupled Receptors

The allosteric modulation of G‐protein‐coupled receptors (GPCRs) by sodium ions has received significant attention as crystal structures of several receptors show Na⁺ ions bound to the inactive conformations at the conserved Asp^2.50. To date, structures from 24 families of GPCRs have been determined, though mechanistic insights into Na⁺ binding to the allosteric site are limited. We performed hundreds‐of‐microsecond long simulations of 18 GPCRs and elucidated their Na⁺ binding mechanism. In class A GPCRs, the Na⁺ ion binds to the conserved residue 2.50 whereas in class B receptors, it binds at 3.43b, 6.53b, and 7.49b. Using Markov state models, we obtained the free energy profiles and kinetics of Na⁺ binding to the allosteric site, which reveal a conserved mechanism of Na⁺ binding for GPCRs and show the residues that act as major barriers for ion diffusion. Furthermore, we also show that the Na⁺ ion can bind to GPCRs from the intracellular side when the allosteric site is inaccessible from the extracellular side.     

Allostery Inhibition of BACE1 by Psychotic and Meroterpenoid Drugs in Alzheimer’s Disease Therapy

In over a century since its discovery, Alzheimer’s disease (AD) has continued to be a global health concern due to its incurable nature and overwhelming increase among older people. In this paper, we give an overview of the efforts of researchers towards identifying potent BACE1 exosite-binding antibodies and allosteric inhibitors. Herein, we apply computer-aided drug design (CADD) methods to unravel the interactions of some proposed psychotic and meroterpenoid BACE1 allosteric site inhibitors. This study is aimed at validating the allosteric potentials of these selected compounds targeted at BACE1 inhibition. Molecular docking, molecular dynamic (MD) simulations, and post-MD analyses are carried out on these selected compounds, which have been experimentally proven to exhibit allosteric inhibition on BACE1. The SwissDock software enabled us to identify more than five druggable pockets on the BACE1 structural surface using docking. Besides the active site region, a melatonin derivative (compound 1) previously proposed as a BACE1 allostery inhibitor showed appreciable stability at eight different subsites on BACE1. Refinement with molecular dynamic (MD) simulations shows that the identified non-catalytic sites are potential allostery sites for compound 1. The allostery and binding mechanism of the selected potent inhibitors show that the smaller the molecule, the easier the attachment to several enzyme regions. This finding hereby establishes that most of these selected compounds failed to exhibit strong allosteric binding with BACE1 except for compound 1. We hereby suggest that further studies and additional identification/validation of other BACE1 allosteric compounds be done. Furthermore, this additional allosteric site investigation will help in reducing the associated challenges with designing BACE1 inhibitors while exploring the opportunities in the design of allosteric BACE1 inhibitors.     

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1−SUMO1 complex.     

Mutant-Selective Allosteric EGFR Degraders are Effective Against a Broad Range of Drug-Resistant Mutations

Targeting epidermal growth factor receptor (EGFR) through an allosteric mechanism provides a potential therapeutic strategy to overcome drug-resistant EGFR mutations that emerge within the ATP binding site. Here, we develop an allosteric EGFR degrader, DDC-01-163, which can selectively inhibit the proliferation of L858R/T790M (L/T) mutant Ba/F3 cells while leaving wildtype EGFR Ba/F3 cells unaffected. DDC-01-163 is also effective against osimertinib-resistant cells with L/T/C797S and L/T/L718Q EGFR mutations. When combined with an ATP-site EGFR inhibitor, osimertinib, the anti-proliferative activity of DDC-01-163 against L858R/T790M EGFR-Ba/F3 cells is enhanced. Collectively, DDC-01-163 is a promising allosteric EGFR degrader with selective activity against various clinically relevant EGFR mutants as a single agent and when combined with an ATP-site inhibitor. Our data suggests that targeted protein degradation is a promising drug development approach for mutant EGFR.     

Mutant‐Selective Allosteric EGFR Degraders are Effective Against a Broad Range of Drug‐Resistant Mutations

Targeting epidermal growth factor receptor (EGFR) through an allosteric mechanism provides a potential therapeutic strategy to overcome drug‐resistant EGFR mutations that emerge within the ATP binding site. Here, we develop an allosteric EGFR degrader, DDC‐01‐163, which can selectively inhibit the proliferation of L858R/T790M (L/T) mutant Ba/F3 cells while leaving wildtype EGFR Ba/F3 cells unaffected. DDC‐01‐163 is also effective against osimertinib‐resistant cells with L/T/C797S and L/T/L718Q EGFR mutations. When combined with an ATP‐site EGFR inhibitor, osimertinib, the anti‐proliferative activity of DDC‐01‐163 against L858R/T790M EGFR‐Ba/F3 cells is enhanced. Collectively, DDC‐01‐163 is a promising allosteric EGFR degrader with selective activity against various clinically relevant EGFR mutants as a single agent and when combined with an ATP‐site inhibitor. Our data suggests that targeted protein degradation is a promising drug development approach for mutant EGFR.     

牛血清白蛋白-Cu2＋结合过程中盐酸小檗碱的变构效应

用核磁共振(1H NMR)、圆二色谱(CD)、荧光光谱(FS)以及紫外光谱(UV)技术考察了中药有效成分盐酸小檗碱(BC)对牛血清白蛋白(BSA)-Cu2＋结合过程的变构效应, 得到分别表征BSA内源荧光猝灭、BSA-Cu2＋复合物稳定性以及Cu2＋在BSA分子上的结合位点发生变构的定量效应参数βQ (βA和βn)和效率参数γQ (γA和γn). 结果表明, BC对Cu2＋猝灭BSA内源荧光呈负变构效应(0＜βQ＜1), 而对BSA-Cu2＋复合物稳定性以及Cu2＋在BSA分子上的结合位点呈正变构效应(βA＞1, βn＞1); 变构效应随BC浓度增加而增强, BC对BSA-Cu2＋复合物稳定性的变构效率明显高于其对荧光猝灭和结合位点的变构; BSA分子构象转变是变构效应的主要原因.     

Computational Study of SARS-CoV-2 RNA Dependent RNA Polymerase Allosteric Site Inhibition

The COVID-19 pandemic has caused millions of fatalities since 2019. Despite the availability of vaccines for this disease, new strains are causing rapid ailment and are a continuous threat to vaccine efficacy. Here, molecular docking and simulations identify strong inhibitors of the allosteric site of the SARS-CoV-2 virus RNA dependent RNA polymerase (RdRp). More than one hundred different flavonoids were docked with the SARS-CoV-2 RdRp allosteric site through computational screening. The three top hits were Naringoside, Myricetin and Aureusidin 4,6-diglucoside. Simulation analyses confirmed that they are in constant contact during the simulation time course and have strong association with the enzyme’s allosteric site. Absorption, distribution, metabolism, excretion and toxicity (ADMET) data provided medicinal information of these top three hits. They had good human intestinal absorption (HIA) concentrations and were non-toxic. Due to high mutation rates in the active sites of the viral enzyme, these new allosteric site inhibitors offer opportunities to drug SARS-CoV-2 RdRp. These results provide new information for the design of novel allosteric inhibitors against SARS-CoV-2 RdRp.     

Allosteric Modulation of the CB1 Cannabinoid Receptor by Cannabidiol—A Molecular Modeling Study of the N-Terminal Domain and the Allosteric-Orthosteric Coupling

The CB₁ cannabinoid receptor (CB₁R) contains one of the longest N termini among class A G protein-coupled receptors. Mutagenesis studies suggest that the allosteric binding site of cannabidiol (CBD) involves residues from the N terminal domain. In order to study the allosteric binding of CBD to CB₁R we modeled the whole N-terminus of this receptor using the replica exchange molecular dynamics with solute tempering (REST2) approach. Then, the obtained structures of CB₁R with the N terminus were used for ligand docking. A natural cannabinoid receptor agonist, Δ⁹-THC, was docked to the orthosteric site and a negative allosteric modulator, CBD, to the allosteric site positioned between extracellular ends of helices TM1 and TM2. The molecular dynamics simulations were then performed for CB₁R with ligands: (i) CBD together with THC, and (ii) THC-only. Analyses of the differences in the residue-residue interaction patterns between those two cases allowed us to elucidate the allosteric network responsible for the modulation of the CB₁R by CBD. In addition, we identified the changes in the orthosteric binding mode of Δ⁹-THC, as well as the changes in its binding energy, caused by the CBD allosteric binding. We have also found that the presence of a complete N-terminal domain is essential for a stable binding of CBD in the allosteric site of CB₁R as well as for the allosteric-orthosteric coupling mechanism.     

A Computational Approach Applied to the Study of Potential Allosteric Inhibitors Protease NS2B/NS3 from Dengue Virus

Dengue virus (DENV) is a danger to more than 400 million people in the world, and there is no specific treatment. Thus, there is an urgent need to develop an effective method to combat this pathology. NS2B/NS3 protease is an important biological target due it being necessary for viral replication and the fact that it promotes the spread of the infection. Thus, this study aimed to design DENV NS2B/NS3pro allosteric inhibitors from a matrix compound. The search was conducted using the Swiss Similarity tool. The compounds were subjected to molecular docking calculations, molecular dynamics simulations (MD) and free energy calculations. The molecular docking results showed that two compounds, ZINC000001680989 and ZINC000001679427, were promising and performed important hydrogen interactions with the Asn152, Leu149 and Ala164 residues, showing the same interactions obtained in the literature. In the MD, the results indicated that five residues, Lys74, Leu76, Asn152, Leu149 and Ala166, contribute to the stability of the ligand at the allosteric site for all of the simulated systems. Hydrophobic, electrostatic and van der Waals interactions had significant effects on binding affinity. Physicochemical properties, lipophilicity, water solubility, pharmacokinetics, druglikeness and medicinal chemistry were evaluated for four compounds that were more promising, showed negative indices for the potential penetration of the Blood Brain Barrier and expressed high human intestinal absorption, indicating a low risk of central nervous system depression or drowsiness as the the side effects. The compound ZINC000006694490 exhibited an alert with a plausible level of toxicity for the purine base chemical moiety, indicating hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms. All of the compounds selected in this study showed a synthetic accessibility (SA) score lower than 4, suggesting the ease of new syntheses. The results corroborate with other studies in the literature, and the computational approach used here can contribute to the discovery of new and potent anti-dengue agents.     

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

     

The Antagonist pGlu-βGlu-Pro-NH2 Binds to an Allosteric Site of the Thyrotropin-Releasing Hormone Receptor

After we identified pGlu-βGlu-Pro-NH₂ as the first functional antagonist of the cholinergic central actions of the thyrotropin-releasing hormone (TRH, pGlu-His-Pro-NH₂), we became interested in finding the receptor-associated mechanism responsible for this antagonism. By utilizing a human TRH receptor (hTRH-R) homology model, we first refined the active binding site within the transmembrane bundle of this receptor to enhance TRH’s binding affinity. However, this binding site did not accommodate the TRH antagonist. This directed us to consider a potential allosteric binding site in the extracellular domain (ECD). Searches for ECD binding pockets prompted the remodeling of the extracellular loops and the N-terminus. We found that different trajectories of ECDs produced novel binding cavities that were then systematically probed with TRH, as well as its antagonist. This led us to establish not only a surface-recognition binding site for TRH, but also an allosteric site that exhibited a selective and high-affinity binding for pGlu-βGlu-Pro-NH₂. The allosteric binding of this TRH antagonist is more robust than TRH’s binding to its own active site. The findings reported here may shed light on the mechanisms and the multimodal roles by which the ECD of a TRH receptor is involved in agonist and/or antagonist actions.     

Small Molecule Tools to Study Cellular Target Engagement for the Intracellular Allosteric Binding Site of GPCRs

A conserved intracellular allosteric binding site (IABS) has recently been identified at several G protein-coupled receptors (GPCRs). Ligands targeting the IABS, so-called intracellular allosteric antagonists, are highly promising compounds for pharmaceutical intervention and currently evaluated in several clinical trials. Beside co-crystal structures that laid the foundation for the structure-based development of intracellular allosteric GPCR antagonists, small molecule tools that enable an unambiguous identification and characterization of intracellular allosteric GPCR ligands are of utmost importance for drug discovery campaigns in this field. Herein, we discuss recent approaches that leverage cellular target engagement studies for the IABS and thus play a critical role in the evaluation of IABS-targeted ligands as potential therapeutic agents.     

New Promise and Opportunities for Allosteric Kinase Inhibitors

Drugs that function through allosteric inhibition of kinase signaling represent a promising approach for the targeted discovery of therapeutics. The majority of developed allosteric kinase inhibitors are characterized as type III and IV inhibitors that show good kinome selectivity but generally lack the subtype selectivity of same kinase family. Recently allosteric inhibitors have been developed that bind outside the catalytic kinase domain with high selectivity for specific kinase subtypes. Allosteric inhibitors that bind to the pseudokinase domain of pseudokinase or the extracellular domain of receptor tyrosine kinases are reviewed. We also review recent developments in the field of allosteric kinase inhibitors including examples of proteolysis targeting chimeras, and highlight the unique binding modes for each type of inhibitors and address future opportunities in this area.     

Computational Investigation on the Allosteric Modulation of Androgen Receptor

Androgens have similar structures with different biological activities. To identify molecular determinants responsible for the activity difference, we have docked six steroidal androgens to the binding site or the surface of androgen receptor by using molecular docking with computational investigation. The energy was calculated respectively based on the QM (quantum mechanics) and MM (molecular mechanics) methods. The result shows that the allosteric modulation of androgen receptor plays an important role in the binding process between androgens and receptor. The open state receptor is less stable than the close state one, but the latter is more favorable for binding with androgens. It is worthy of note that when the androgen receptors binding or without binding with androgen are in close state, they are difficult to return to their open state. This phenomenon is an exception of the well known two-state model theory in which the two states are reversible. Whether the internal of close state androgen receptor has a combination of androgen or not, the androgen receptor surface can be combined with another androgen, and their surface binding energies could be very close. The result is consistent with the experimental observations, but this phenomenon of continuous combination from open state is also an exception of the two-state model theory.     

Parameterization of Org27569: An allosteric modulator of the cannabinoid CB1 G protein‐coupled receptor

The cannabinoid CB₁ receptor is a class A G protein‐coupled receptor (GPCR) that is the most widely expressed GPCR in the brain. Many GPCRs contain allosteric binding sites for endogenous and/or synthetic ligands, which are topographically distinct from the agonist‐binding site that is known as the orthosteric site. While both endogenous and synthetic ligands that act at the CB₁ orthosteric site have been known for some time, compounds that act at a CB₁ allosteric site have only recently been discovered. The most studied of these is 5‐chloro‐3‐ethyl‐1H‐indole‐2‐carboxylic acid [2‐(4‐piperidin‐1‐ylphenyl)ethyl]amide (Org27569). Because allosteric ligands are thought to act through conformational changes in the receptor that are transmitted from the allosteric to the orthosteric site, computational studies of the structural and dynamic interactions of Org27569 with the CB₁ receptor are crucial to achieve a molecular level understanding of the basis of action of this important new class of compounds. To date, such computational studies have not been possible due to the lack of a complete set of molecular mechanics force field parameters for Org27569. Here, we present the development of missing CHARMM force field parameters for Org27569 using previously published methods and the validation and application of these new parameters using normal mode analysis and molecular dynamics simulations combined with experimental infrared measurements. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Structure-activity analysis of the purine binding site of human liver glycogen phosphorylase

Human liver glycogen phosphorylase (HLGP) catalyzes the breakdown of glycogen to maintain serum glucose levels and is a therapeutic target for diabetes. HLGP is regulated by multiple interacting allosteric sites, each of which is a potential drug binding site. We used surface plasmon resonance (SPR) to screen for compounds that bind to the purine allosteric inhibitor site. We determined the affinities of a series of compounds and solved the crystal structures of three representative ligands with K(D) values from 17-550 microM. The crystal structures reveal that the affinities are partly determined by ligand-specific water-mediated hydrogen bonds and side chain movements. These effects could not be predicted; both crystallographic and SPR studies were required to understand the important features of binding and together provide a basis for the design of new allosteric inhibitors targeting this site.     

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

Allostery, which is one of the most direct and efficient methods to fine-tune protein functions, has gained increasing recognition in drug discovery. However, there are several challenges associated with the identification of allosteric sites, which is the fundamental cornerstone of drug design. Previous studies on allosteric site predictions have focused on communication signals propagating from the allosteric sites to the orthosteric sites. However, recent biochemical studies have revealed that allosteric coupling is bidirectional and that orthosteric perturbations can modulate allosteric sites through reversed allosteric communication. Here, we proposed a new framework for the prediction of allosteric sites based on reversed allosteric communication using a combination of computational and experimental strategies (molecular dynamics simulations, Markov state models, and site-directed mutagenesis). The desirable performance of our approach was demonstrated by predicting the known allosteric site of the small molecule MDL-801 in nicotinamide dinucleotide (NAD⁺)-dependent protein lysine deacetylase sirtuin 6 (Sirt6). A potential novel cryptic allosteric site located around the L116, R119, and S120 residues within the dynamic ensemble of Sirt6 was identified. The allosteric effect of the predicted site was further quantified and validated using both computational and experimental approaches. This study proposed a state-of-the-art computational pipeline for detecting allosteric sites based on reversed allosteric communication. This method enabled the identification of a previously uncharacterized potential cryptic allosteric site on Sirt6, which provides a starting point for allosteric drug design that can aid the identification of candidate pockets in other therapeutic targets.

Using reversed allosteric communication, we performed MD simulations, MSMs, and mutagenesis experiments, to discover allosteric sites. It reproduced the known allosteric site for MDL-801 on Sirt6 and uncovered a novel cryptic allosteric Pocket X.     

Allosteric Regulation of Enzymatic Reactions in a Transparent Inorganic Sol-Gel Material

Glutamate dehydrogenase is encapsulated in a transparent porous silicate matrix by using sol-gel techniques. The inorganic polymer is formed around the enzyme (MW > 300,000 D). The enzyme is active in the material, catalyzes the reaction of L-glutamate to 2-oxoglutarate and follows Michaelis-Menten kinetics. The allosteric regulators ADP and GTP inhibit or activate the reaction; at pH 6, GTP acts as a strong activator and ADP acts as an inhibitor. This system involves a complex series of interactions; the co-enzyme NAD⁺ is required for catalysis, large-scale conformational changes accompany the binding of the substrate and coenzyme to the enzyme, the activators/inhibitors must bind to the enzyme to regulate the reactions, and the substrates and products must diffuse through the matrix to and from the binding site. The influence of the unique matrix on the complex enzymatic system is discussed.     

Equilibrium dialysis of metal-serum albumin (II) Allosteric effect in Ni(II)-serum albumin systems

Detailed studies were carried out on equilibrium dialysis of the binding of Ni²⁺ ion to human serum albumin (HSA) and bovine serum albumin (BSA). The successive stability constants were obtained by the leastsquares fitting. The eight binding sites found for both Ni(II)-HSA and Ni(II)-BSA systems can be divided into two different sets; and for both systems, there exist two identical prior binding sites where the bound Ni²⁺ ions can be considered as allosteric effectors, which induce the allosteric effect in accordance with the model proposed by Monodet al. As indicated by allosteric parameters, the ability of R-state to bind Ni²⁺ ionions is ca. 100 times as much as that of T-state, and the conformation of HSA is markedly tenser than that of BSA. Project supported by the National Natural Science Foundation of China.