Binding Kinetics and Pathways of Ligands to GPCRs

12-HETE/GPR31 induces endothelial dysfunction in diabetic retinopathy

12-hydroxyeicosatetraenoic acid (12-HETE), a major metabolite of arachidonic acid, is converted by 12/15-lipoxygenase and implicated in diabetic retinopathy (DR). Our previous study demonstrated a positive correlation between 12-HETE and the prevalence of DR. However, reasons for the increased production of 12-HETE are unclear, and the underlying mechanisms through which 12-HETE promotes DR are unknown. This study aimed to elucidate the correlation between 12-HETE and DR onset, investigate potential mechanisms through which 12-HETE promotes DR, and seek explanations for the increased production of 12-HETE in diabetes. We conducted a prospective cohort study, which revealed that higher serum 12-HETE levels could induce DR. Additionally, G protein-coupled receptor 31 (GPR31), a high-affinity receptor for 12-HETE, was expressed in human retinal microvascular endothelial cells (HRMECs). 12-HETE/GPR31-mediated HRMEC inflammation occurred via the p38 MAPK pathway. 12-HETE levels were significantly higher in the retina of mice with high-fat diet (HFD)- and streptozotocin (STZ)-induced diabetes than in those with only STZ-induced diabetes and healthy controls. They were positively correlated with the levels of inflammatory cytokines in the retina, indicating that HFD could induce increased 12-HETE synthesis in patients with diabetes in addition to hyperglycemia. Conclusively, 12-HETE is a potential risk factor for DR. The 12-HETE/GPR31 axis plays a crucial role in HRMEC dysfunction and could be a novel target for DR prevention and control. Nevertheless, further research is warranted to provide comprehensive insights into the complex underlying mechanisms of 12-HETE in DR.

Know your molecule: pharmacological characterization of drug candidates to enhance efficacy and reduce late-stage attrition

Despite advances in chemical, computational and biological sciences, the rate of attrition of drug candidates in clinical development is still high. A key point in the small-molecule discovery process that could provide opportunities to help address this challenge is the pharmacological characterization of hit and lead compounds, culminating in the selection of a drug candidate. Deeper characterization is increasingly important, because the 'quality' of drug efficacy, at least for G protein-coupled receptors (GPCRs), is now understood to be much more than activation of commonly evaluated pathways such as cAMP signalling, with many more 'efficacies' of ligands that could be harnessed therapeutically. Such characterization is being enabled by novel assays to characterize the complex behaviour of GPCRs, such as biased signalling and allosteric modulation, as well as advances in structural biology, such as cryo-electron microscopy. This article discusses key factors in the assessments of the pharmacology of hit and lead compounds in the context of GPCRs as a target class, highlighting opportunities to identify drug candidates with the potential to address limitations of current therapies and to improve the probability of them succeeding in clinical development.

A Pharmacological perspective on the temporal properties of sweeteners

Several non-caloric sweeteners exhibit a delay in sweetness onset and a sweetness linger after sampling. These temporal properties are thought to be the result of non-specific interactions with cell membranes and proteins in the oral cavity. Data and analysis presented in this report also support the potential involvement of receptor affinity and binding kinetics to this phenomenon. In general, affected sweeteners exhibit distinctly higher binding affinity compared to carbohydrate sweeteners, which do not have temporal issues. In addition, binding kinetic simulations illustrate much slower receptor binding association and dissociation kinetics for a set of non-caloric sweeteners presenting temporal issues, in comparison to carbohydrate sweeteners. So, the higher affinity of some non-caloric sweeteners, dictating lower use levels, and affecting binding kinetics, could contribute to their delay and linger in sweetness perception. Simple pharmacology principles could explain, at least in part, some of the temporal issues of sweeteners.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

Activity Models of Key GPCR Families in the Central Nervous System: A Tool for Many Purposes

G protein-coupled receptors (GPCRs) are targets of many drugs, of which ∼25% are indicated for central nervous system (CNS) disorders. Drug promiscuity affects their efficacy and safety profiles. Predicting the polypharmacology profile of compounds against GPCRs can thus provide a basis for producing more precise therapeutics by considering the targets and the anti-targets in that family of closely related proteins. We provide a tool for predicting the polypharmacology of compounds within prominent GPCR families in the CNS: serotonin, dopamine, histamine, muscarinic, opioid, and cannabinoid receptors. Our in-house algorithm, "iterative stochastic elimination" (ISE), produces high-quality ligand-based models for agonism and antagonism at 31 GPCRs. The ISE models correctly predict 68% of CNS drug-GPCR interactions, while the "similarity ensemble approach" predicts only 33%. The activity models correctly predict 56% of reported activities of DrugBank molecules for these CNS receptors. We conclude that the combination of interactions and activity profiles generated by screening through our models form the basis for subsequent designing and discovering novel therapeutics, either single, multitargeting, or repurposed.

Identification, Screening, and Comprehensive Evaluation of Novel DPP-IV Inhibitory Peptides from the Tilapia Skin Gelatin Hydrolysate Produced Using Ginger Protease

PURPOSE

Inhibition of dipeptidyl peptidase-IV (DPP-IV) is an effective therapy for treating type II diabetes (T2D) that has been widely applied in clinical practice. We aimed to evaluate the DPP-IV inhibitory properties of ginger protease hydrolysate (GPH) and propose a comprehensive approach to screen and evaluate DPP-IV inhibitors.

METHODS

We evaluated the in vitro inhibitory properties of fish skin gelatin hydrolysates produced by five proteases, namely, neutral protease, alkaline protease, bromelain, papain, and ginger protease, toward DPP-IV. We screened the most potent DPP-IV inhibitory peptide (DIP) using liquid chromatography-tandem mass spectrometry (LC-MS/MS) coupled with in silico analysis. Next, surface plasmon resonance (SPR) technology was innovatively introduced to explore the interactions between DPP-IV and DIP, as well as the IC₅₀. Furthermore, we performed oral administration of DIP in rats to study its in vivo absorption.

RESULTS

GPH displayed the highest degree of hydrolysis (20.37%) and DPP-IV inhibitory activity (65.18%). A total of 292 peptides from the GPH were identified using LC-MS/MS combined with de novo sequencing. Gly-Pro-Hyp-Gly-Pro-Pro-Gly-Pro-Gly-Pro (GPXGPPGPGP) was identified as the most potent DPP-IV inhibitory peptide after in silico screening (Peptide Ranker and molecular docking). Then, the in vitro study revealed that GPXGPPGPGP had a high inhibitory effect on DPP-IV (IC₅₀: 1012.3 ± 23.3 μM) and exhibited fast kinetics with rapid binding and dissociation with DPP-IV. In vivo analysis indicated that GPXGPPGPGP was not absorbed intact but partially, in the form of dipeptides and tripeptides.

CONCLUSION

Overall, the results suggested that GPH would be a natural functional food for treating T2D and provided new ideas for searching and evaluating potential antidiabetic compounds. The obtained GPXGPPGPGP can be structurally optimized for in-depth evaluation in animal and cellular experiments.

A Single Point Mutation Blocks the Entrance of Ligands to the Cannabinoid CB2 Receptor via the Lipid Bilayer

Molecular dynamic (MD) simulations have become a common tool to study the pathway of ligand entry to the orthosteric binding site of G protein-coupled receptors. Here, we have combined MD simulations and site-directed mutagenesis to study the binding process of the potent JWH-133 agonist to the cannabinoid CB₂ receptor (CB₂R). In CB₂R, the N-terminus and extracellular loop 2 fold over the ligand binding pocket, blocking access to the binding cavity from the extracellular environment. We, thus, hypothesized that the binding pathway is a multistage process consisting of the hydrophobic ligand diffusing in the lipid bilayer to contact a lipid-facing vestibule, from which the ligand enters an allosteric site inside the transmembrane bundle through a tunnel formed between TMs 1 and 7 and finally moving from the allosteric to the orthosteric binding cavity. This pathway was experimentally validated by the Ala282^7.36Phe mutation that blocks the entrance of the ligand, as JWH-133 was not able to decrease the forskolin-induced cAMP levels in cells expressing the mutant receptor. This proposed ligand entry pathway defines transient binding sites that are potential cavities for the design of synthetic modulators.

In-cell NMR: Why and how?

NMR spectroscopy has been applied to cells and tissues analysis since its beginnings, as early as 1950. We have attempted to gather here in a didactic fashion the broad diversity of data and ideas that emerged from NMR investigations on living cells. Covering a large proportion of the periodic table, NMR spectroscopy permits scrutiny of a great variety of atomic nuclei in all living organisms non-invasively. It has thus provided quantitative information on cellular atoms and their chemical environment, dynamics, or interactions. We will show that NMR studies have generated valuable knowledge on a vast array of cellular molecules and events, from water, salts, metabolites, cell walls, proteins, nucleic acids, drugs and drug targets, to pH, redox equilibria and chemical reactions. The characterization of such a multitude of objects at the atomic scale has thus shaped our mental representation of cellular life at multiple levels, together with major techniques like mass-spectrometry or microscopies. NMR studies on cells has accompanied the developments of MRI and metabolomics, and various subfields have flourished, coined with appealing names: fluxomics, foodomics, MRI and MRS (i.e. imaging and localized spectroscopy of living tissues, respectively), whole-cell NMR, on-cell ligand-based NMR, systems NMR, cellular structural biology, in-cell NMR… All these have not grown separately, but rather by reinforcing each other like a braided trunk. Hence, we try here to provide an analytical account of a large ensemble of intricately linked approaches, whose integration has been and will be key to their success. We present extensive overviews, firstly on the various types of information provided by NMR in a cellular environment (the "why", oriented towards a broad readership), and secondly on the employed NMR techniques and setups (the "how", where we discuss the past, current and future methods). Each subsection is constructed as a historical anthology, showing how the intrinsic properties of NMR spectroscopy and its developments structured the accessible knowledge on cellular phenomena. Using this systematic approach, we sought i) to make this review accessible to the broadest audience and ii) to highlight some early techniques that may find renewed interest. Finally, we present a brief discussion on what may be potential and desirable developments in the context of integrative studies in biology.

Pharmacology of Kappa Opioid Receptors: Novel Assays and Ligands

The present study investigated the in vitro pharmacology of the human kappa opioid receptor using multiple assays, including calcium mobilization in cells expressing chimeric G proteins, the dynamic mass redistribution (DMR) label-free assay, and a bioluminescence resonance energy transfer (BRET) assay that allows measurement of receptor interaction with G protein and β-arrestin 2. In all assays, dynorphin A, U-69,593, and [D-Pro10]dyn(1-11)-NH2 behaved as full agonists with the following rank order of potency [D-Pro10]dyn(1-11)-NH2 > dynorphin A ≥ U-69,593. [Dmt1,Tic2]dyn(1-11)-NH2 behaved as a moderate potency pure antagonist in the kappa-β-arrestin 2 interaction assay and as low efficacy partial agonist in the other assays. Norbinaltorphimine acted as a highly potent and pure antagonist in all assays except kappa-G protein interaction, where it displayed efficacy as an inverse agonist. The pharmacological actions of novel kappa ligands, namely the dynorphin A tetrameric derivative PWT2-Dyn A and the palmitoylated derivative Dyn A-palmitic, were also investigated. PWT2-Dyn A and Dyn A-palmitic mimicked dynorphin A effects in all assays showing similar maximal effects but 3–10 fold lower potency. In conclusion, in the present study, multiple in vitro assays for the kappa receptor have been set up and pharmacologically validated. In addition, PWT2-Dyn A and Dyn A-palmitic were characterized as potent full agonists; these compounds are worthy of further investigation in vivo for those conditions in which the activation of the kappa opioid receptor elicits beneficial effects e.g. pain and pruritus.

In-Cell Structural Biology by NMR: The Benefits of the Atomic Scale

In-cell structural biology aims at extracting structural information about proteins or nucleic acids in their native, cellular environment. This emerging field holds great promise and is already providing new facts and outlooks of interest at both fundamental and applied levels. NMR spectroscopy has important contributions on this stage: It brings information on a broad variety of nuclei at the atomic scale, which ensures its great versatility and uniqueness. Here, we detail the methods, the fundamental knowledge, and the applications in biomedical engineering related to in-cell structural biology by NMR. We finally propose a brief overview of the main other techniques in the field (EPR, smFRET, cryo-ET, etc.) to draw some advisable developments for in-cell NMR. In the era of large-scale screenings and deep learning, both accurate and qualitative experimental evidence are as essential as ever to understand the interior life of cells. In-cell structural biology by NMR spectroscopy can generate such a knowledge, and it does so at the atomic scale. This review is meant to deliver comprehensive but accessible information, with advanced technical details and reflections on the methods, the nature of the results, and the future of the field.

Two entry tunnels in mouse TAAR9 suggest the possibility of multi-entry tunnels in olfactory receptors

Orthosteric binding sites of olfactory receptors have been well understood for ligand-receptor interactions. However, a lack of explanation for subtle differences in ligand profile of olfactory receptors even with similar orthosteric binding sites promotes more exploration into the entry tunnels of the receptors. An important question regarding entry tunnels is the number of entry tunnels, which was previously believed to be one. Here, we used TAAR9 that recognizes important biogenic amines such as cadaverine, spermine, and spermidine as a model for entry tunnel study. We identified two entry tunnels in TAAR9 and described the residues that form the tunnels. In addition, we found two vestibular binding pockets, each located in one tunnel. We further confirmed the function of two tunnels through site-directed mutagenesis. Our study challenged the existing views regarding the number of entry tunnels in the subfamily of olfactory receptors and demonstrated the possible mechanism how the entry tunnels function in odorant recognition.

Biased, Bitopic, Opioid-Adrenergic Tethered Compounds May Improve Specificity, Lower Dosage and Enhance Agonist or Antagonist Function with Reduced Risk of Tolerance and Addiction

This paper proposes the design of combination opioid-adrenergic tethered compounds to enhance efficacy and specificity, lower dosage, increase duration of activity, decrease side effects, and reduce risk of developing tolerance and/or addiction. Combinations of adrenergic and opioid drugs are sometimes used to improve analgesia, decrease opioid doses required to achieve analgesia, and to prolong the duration of analgesia. Recent mechanistic research suggests that these enhanced functions result from an allosteric adrenergic binding site on opioid receptors and, conversely, an allosteric opioid binding site on adrenergic receptors. Dual occupancy of the receptors maintains the receptors in their high affinity, most active states; drops the concentration of ligand required for full activity; and prevents downregulation and internalization of the receptors, thus inhibiting tolerance to the drugs. Activation of both opioid and adrenergic receptors also enhances heterodimerization of the receptors, additionally improving each drug's efficacy. Tethering adrenergic drugs to opioids could produce new drug candidates with highly desirable features. Constraints-such as the locations of the opioid binding sites on adrenergic receptors and adrenergic binding sites on opioid receptors, length of tethers that must govern the design of such novel compounds, and types of tethers-are described and examples of possible structures provided.

Kinetics and Mechanism of Fentanyl Dissociation from the μ-Opioid Receptor

Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid overdose is resuscitated by the use of naloxone, which competitively binds and antagonizes the μ-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study, we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H297^6.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the Nδ-protonated H297^6.52, fentanyl can descend as much as 10 Å below the level of the conserved D147^3.32 before escaping the receptor and has a calculated residence time τ of 38 s. In contrast, with the Nϵ- and doubly protonated H297^6.52, the calculated τ are 2.6 and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental τ of ∼4 min for fentanyl and the role of TM6 in mOR activation, the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform evaluation of strategies for opioid overdose reversal. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics.

Identification and Functional Analysis of G Protein-Coupled Receptors in 20-Hydroxyecdysone Signaling From the Helicoverpa armigera Genome

G protein-coupled receptors (GPCRs) are the largest family of membrane receptors in animals and humans, which transmit various signals from the extracellular environment into cells. Studies have reported that several GPCRs transmit the same signal; however, the mechanism is unclear. In the present study, we identified all 122 classical GPCRs from the genome of Helicoverpa armigera, a lepidopteran pest species. Twenty-four GPCRs were identified as upregulated at the metamorphic stage by comparing the transcriptomes of the midgut at the metamorphic and feeding stages. Nine of them were confirmed to be upregulated at the metamorphic stage. RNA interference in larvae revealed the prolactin-releasing peptide receptor (PRRPR), smoothened (SMO), adipokinetic hormone receptor (AKHR), and 5-hydroxytryptamine receptor (HTR) are involved in steroid hormone 20-hydroxyecdysone (20E)-promoted pupation. Frizzled 7 (FZD7) is involved in growth, while tachykinin-like peptides receptor 86C (TKR86C) had no effect on growth and pupation. Via these GPCRs, 20E regulated the expression of different genes, respectively, including Pten (encoding phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase), FoxO (encoding forkhead box O), BrZ7 (encoding broad isoform Z7), Kr-h1 (encoding Krüppel homolog 1), Wnt (encoding Wingless/Integrated) and cMyc, with hormone receptor 3 (HHR3) as their common regulating target. PRRPR was identified as a new 20E cell membrane receptor using a binding assay. These data suggested that 20E, via different GPCRs, regulates different gene expression to integrate growth and development.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from τRAMD Simulations

There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of time scales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the τ-Random Acceleration Molecular Dynamics (τRAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the β2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

Neuronal metabotropic glutamate receptor 8 protects against neurodegeneration in CNS inflammation

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system with continuous neuronal loss. Treatment of clinical progression remains challenging due to lack of insights into inflammation-induced neurodegenerative pathways. Here, we show that an imbalance in the neuronal receptor interactome is driving glutamate excitotoxicity in neurons of MS patients and identify the MS risk-associated metabotropic glutamate receptor 8 (GRM8) as a decisive modulator. Mechanistically, GRM8 activation counteracted neuronal cAMP accumulation, thereby directly desensitizing the inositol 1,4,5-trisphosphate receptor (IP3R). This profoundly limited glutamate-induced calcium release from the endoplasmic reticulum and subsequent cell death. Notably, we found Grm8-deficient neurons to be more prone to glutamate excitotoxicity, whereas pharmacological activation of GRM8 augmented neuroprotection in mouse and human neurons as well as in a preclinical mouse model of MS. Thus, we demonstrate that GRM8 conveys neuronal resilience to CNS inflammation and is a promising neuroprotective target with broad therapeutic implications.

Binding pathway determines norepinephrine selectivity for the human β1AR over β2AR

Beta adrenergic receptors (βARs) mediate physiologic responses to the catecholamines epinephrine and norepinephrine released by the sympathetic nervous system. While the hormone epinephrine binds β₁AR and β₂AR with similar affinity, the smaller neurotransmitter norepinephrine is approximately tenfold selective for the β₁AR. To understand the structural basis for this physiologically important selectivity, we solved the crystal structures of the human β₁AR bound to an antagonist carazolol and different agonists including norepinephrine, epinephrine and BI-167107. Structural comparison revealed that the catecholamine-binding pockets are identical between β₁AR and β₂AR, but the extracellular vestibules have different shapes and electrostatic properties. Metadynamics simulations and mutagenesis studies revealed that these differences influence the path norepinephrine takes to the orthosteric pocket and contribute to the different association rates and thus different affinities.

Dopamine D2 Receptor Agonist Binding Kinetics-Role of a Conserved Serine Residue

The forward (k_on) and reverse (k_off) rate constants of drug–target interactions have important implications for therapeutic efficacy. Hence, time-resolved assays capable of measuring these binding rate constants may be informative to drug discovery efforts. Here, we used an ion channel activation assay to estimate the k_ons and k_offs of four dopamine D₂ receptor (D₂R) agonists; dopamine (DA), p-tyramine, (R)- and (S)-5-OH-dipropylaminotetralin (DPAT). We further probed the role of the conserved serine S193^5.42 by mutagenesis, taking advantage of the preferential interaction of (S)-, but not (R)-5-OH-DPAT with this residue. Results suggested similar k_offs for the two 5-OH-DPAT enantiomers at wild-type (WT) D₂R, both being slower than the k_offs of DA and p-tyramine. Conversely, the k_on of (S)-5-OH-DPAT was estimated to be higher than that of (R)-5-OH-DPAT, in agreement with the higher potency of the (S)-enantiomer. Furthermore, S193^5.42A mutation lowered the k_on of (S)-5-OH-DPAT and reduced the potency difference between the two 5-OH-DPAT enantiomers. Kinetic K_ds derived from the k_off and k_on estimates correlated well with EC₅₀ values for all four compounds across four orders of magnitude, strengthening the notion that our assay captured meaningful information about binding kinetics. The approach presented here may thus prove valuable for characterizing D₂R agonist candidate drugs.

Deciphering the Agonist Binding Mechanism to the Adenosine A1 Receptor

Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A1R) was determined in the active, Gi protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A1R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N6-cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A1R agonists.

Establishment of a NanoBiT-Based Cytosolic Ca2+ Sensor by Optimizing Calmodulin-Binding Motif and Protein Expression Levels

Cytosolic Ca2+ levels ([Ca2+]c) change dynamically in response to inducers, repressors, and physiological conditions, and aberrant [Ca2+]c concentration regulation is associated with cancer, heart failure, and diabetes. Therefore, [Ca2+]c is considered as a good indicator of physiological and pathological cellular responses, and is a crucial biomarker for drug discovery. A genetically encoded calcium indicator (GECI) was recently developed to measure [Ca2+]c in single cells and animal models. GECI have some advantages over chemically synthesized indicators, although they also have some drawbacks such as poor signal-to-noise ratio (SNR), low positive signal, delayed response, artifactual responses due to protein overexpression, and expensive detection equipment. Here, we developed an indicator based on interactions between Ca2+-loaded calmodulin and target proteins, and generated an innovative GECI sensor using split nano-luciferase (Nluc) fragments to detect changes in [Ca2+]c. Stimulation-dependent luciferase activities were optimized by combining large and small subunits of Nluc binary technology (NanoBiT, LgBiT:SmBiT) fusion proteins and regulating the receptor expression levels. We constructed the binary [Ca2+]c sensors using a multicistronic expression system in a single vector linked via the internal ribosome entry site (IRES), and examined the detection efficiencies. Promoter optimization studies indicated that promoter-dependent protein expression levels were crucial to optimize SNR and sensitivity. This novel [Ca2+]c assay has high SNR and sensitivity, is easy to use, suitable for high-throughput assays, and may be useful to detect [Ca2+]c in single cells and animal models.

Discovery and Characterization of BAY 1214784, an Orally Available Spiroindoline Derivative Acting as a Potent and Selective Antagonist of the Human Gonadotropin-Releasing Hormone Receptor as Proven in a First-In-Human Study in Postmenopausal Women

The growth of uterine fibroids is sex hormone-dependent and commonly associated with highly incapacitating symptoms. Most treatment options consist of the control of these hormonal effects, ultimately blocking proliferative estrogen signaling (i.e., oral contraceptives/antagonization of human gonadotropin-releasing hormone receptor [hGnRH-R] activity). Full hGnRH-R blockade, however, results in menopausal symptoms and affects bone mineralization, thus limiting treatment duration or demanding estrogen add-back approaches. To overcome such issues, we aimed to identify novel, small-molecule hGnRH-R antagonists. This led to the discovery of compound BAY 1214784, an orally available, potent, and selective hGnRH-R antagonist. Altering the geminal dimethylindoline core of the initial hit compound to a spiroindoline system significantly improved GnRH-R antagonist potencies across several species, mandatory for a successful compound optimization in vivo. In a first-in-human study in postmenopausal women, once daily treatment with BAY 1214784 effectively lowered plasma luteinizing hormone levels by up to 49%, at the same time being associated with low pharmacokinetic variability and good tolerability.

Application of the Nested Enzyme-Within-Enterocyte (NEWE) Turnover Model for Predicting the Time Course of Pharmacodynamic Effects

The gut wall consists of many biological elements, including enterocytes. Rapid turnover, a prominent feature of the enterocytes, has generally been ignored in the development of enterocyte-targeting drugs, although it has a comparable rate to other kinetic rates. Here, we investigated the impact of enterocyte turnover on the pharmacodynamics of enterocyte-targeting drugs by applying a model accounting for turnover of enterocytes and target proteins. Simulations showed that the pharmacodynamics depend on enterocyte lifespan when drug-target affinity is strong and half-life of target protein is long. Interindividual variability of enterocyte lifespan, which can be amplified by disease conditions, has a substantial impact on the variability of response. However, our comprehensive literature search showed that the enterocyte turnover causes a marginal impact on currently approved enterocyte-targeting drugs due to their relatively weak target affinities. This study proposes a model-informed drug development approach for selecting enterocyte-targeting drugs and their optimal dosage regimens.

In Silico Design, Synthesis and Evaluation of Novel Series of Benzothiazole- Based Pyrazolidinediones as Potent Hypoglycemic Agents

Background:

The discovery of novel ligand binding domain (LBD) of peroxisome proliferator- activated receptor γ (PPARγ) has recently attracted attention to few research groups in order to develop more potent and safer antidiabetic agents.

Objective:

This study is focused on docking-based design and synthesis of novel compounds combining benzothiazole and pyrazolidinedione scaffold as potential antidiabetic agents.

Methods:

Several benzothiazole-pyrazolidinedione hybrids were synthesized and tested for their in vivo anti-hyperglycemic activity. Interactions profile of title compounds against PPARγ was examined through molecular modelling approach.

Results:

All tested compounds exhibited anti-hyperglycemic activity similar or superior to the reference drug Rosiglitazone. Introducing chlorine atom and alkyl group at position-6 and -5 respectively on benzothiazole core resulted in enhancing the anti-hyperglycemic effect. Docking study revealed that such groups demonstrated favorable hydrophobic interactions with novel LBD Ω- pocket of PPARγ protein.

Conclusion:

Among the tested compounds, N-(6-chloro-5-methylbenzo[d]thiazol-2-yl-4-(4((3,5- dioxopyrazolidin-4-ylidene)methyl)phenoxy)butanamide 5b was found to be the most potent compound and provided valuable insights to further develop novel hybrids as anti-hyperglycemic agents.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Detailed In Vitro Pharmacological Characterization of Clinically Tested Negative Allosteric Modulators of the Metabotropic Glutamate Receptor 5

Negative allosteric modulation of the metabotropic glutamate 5 (mGlu₅) receptor has emerged as a potential strategy for the treatment of neurologic disorders. Despite the success in preclinical studies, many mGlu₅ negative allosteric modulators (NAMs) that have reached clinical trials failed due to lack of efficacy. In this study, we provide a detailed in vitro pharmacological characterization of nine clinically and preclinically tested NAMs. We evaluated inhibition of l-glutamate-induced signaling with Ca²⁺ mobilization, inositol monophosphate (IP₁) accumulation, extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation, and real-time receptor internalization assays on rat mGlu₅ expressed in HEK293A cells. Moreover, we determined association rates (k_on) and dissociation rates (k_off), as well as NAM affinities with [³H]methoxy-PEPy binding experiments. k_on and k_off values varied greatly between the nine NAMs (34- and 139-fold, respectively) resulting in long receptor residence times (>400 min) for basimglurant and mavoglurant, medium residence times (10-30 min) for AZD2066, remeglurant, and (RS)-remeglurant, and low residence times (<10 mins) for dipraglurant, F169521, F1699611, and STX107. We found that all NAMs inhibited l-glutamate-induced mGlu₅ receptor internalization, generally with a similar potency to IP₁ accumulation and ERK1/2 phosphorylation, whereas Ca²⁺ mobilization was less potently inhibited. Operational model of allosterism analyses revealed that dipraglurant and (RS)-remeglurant were biased toward (affinity) receptor internalization and away (cooperativity) from the ERK1/2 phosphorylation pathway, respectively. Our study is the first to measure mGlu₅ NAM binding kinetics and negative allosteric modulation of mGlu₅ receptor internalization and adds significant new knowledge about the molecular pharmacology of a diverse range of clinically relevant NAMs. SIGNIFICANCE STATEMENT: The metabotropic glutamate 5 (mGlu₅) receptor is important in many brain functions and implicated in several neurological pathologies. Negative allosteric modulators (NAMs) have shown promising results in preclinical models but have so far failed in human clinical trials. Here we provide the most comprehensive and comparative molecular pharmacological study to date of nine preclinically/clinically tested NAMs at the mGlu₅ receptor, which is also the first study to measure ligand binding kinetics and negative allosteric modulation of mGlu₅ receptor internalization.

Dopamine Receptor Partial Agonists for the Treatment of Bipolar Disorder

Bipolar disorder is a chronic, disabling, and costly illness with frequent relapses and recurrences, high rates of co-morbid conditions, and poor adherence to treatment. Mood stabilizers and antipsychotics are the cornerstones of treatment. Dopamine receptor partial agonists are a novel class of antipsychotic agents with original pharmacodynamic properties. Among them, two have been approved by the US Food and Drug Administration for the treatment of bipolar disorder. Aripiprazole (oral formulation) has been approved as monotherapy for the treatment of manic/mixed episodes in adult and pediatric populations and for maintenance treatment in adults, and as adjunctive treatment to mood stabilizers, for the acute treatment of manic/mixed episodes and for maintenance in adults. An intramuscular formulation of aripiprazole has been approved for the treatment of agitation in mania and a long-acting injectable formulation has been approved as maintenance treatment. In the USA, cariprazine has been approved as monotherapy for the acute treatment of manic/mixed as well as bipolar depressive episodes. Brexpiprazole is not yet approved to treat bipolar disorder. The evidence supporting these indications is reviewed via an analysis of clinical registration trials as well as additional studies, on the basis of a systematic literature search. Further studies dealing with other aspects of bipolar illness are also presented. Aripiprazole and cariprazine are efficacious and generally well tolerated agents that have shown cost effectiveness, and may therefore enrich our therapeutic armamentarium for bipolar illness. Brexpiprazole, which displays an overall promising tolerability profile, deserves further efficacy studies.

Membranes as the third genetic code

Biological membranes and their compositions influence cellular function, age and disease states of organisms. They achieve this by effecting the outcome of bound enzymes/proteins and carbohydrate moieties. While the membrane-bound carbohydrates give rise to antigenicity, membranes impact the eventual outcome of protein structures that would function even outside their enclosure. This is achieved by membrane modulation of translational and post-translational protein folding. Thus, the eventual 3D structures and functions of proteins might not be solely dependent on their primary amino acid sequences and surrounding environments. The 3D protein structures would also depend on enclosing membrane properties such as fluidity, other intrinsic and extrinsic proteins and carbohydrate functionalities. Also, membranes moderate DNA activities with consequences on gene activation-inactivation mechanisms. Consequently, membranes are almost indispensable to the functioning of other cell compositions and serve to modulate these other components. Besides, membrane lipid compositions are also moderated by nutrition and diets and the converse is true. Thus, it could be argued that membranes are the third genetic codes. Suggestively, membranes are at the center of the interplay between nature and nurture in health and disease states.

Crystal structure of the M5 muscarinic acetylcholine receptor

The 5 subtypes of the muscarinic acetylcholine receptors (mAChRs) are expressed throughout the central and peripheral nervous system where they play a vital role in physiology and pathologies. Recently, the M₅ mAChR subtype has emerged as an exciting drug target for the treatment of drug addiction. We have determined the atomic structure of the M₅ mAChR bound to the clinically used inverse agonist tiotropium. The M₅ mAChR structure now allows for a full comparison of all 5 mAChR subtypes and reveals that small differences in the extracellular loop regions can mediate orthosteric and allosteric ligand selectivity. Together, these findings open the door for future structure-based design of selective drugs that target this therapeutically important class of receptors.

The human M₅ muscarinic acetylcholine receptor (mAChR) has recently emerged as an exciting therapeutic target for treating a range of disorders, including drug addiction. However, a lack of structural information for this receptor subtype has limited further drug development and validation. Here we report a high-resolution crystal structure of the human M₅ mAChR bound to the clinically used inverse agonist, tiotropium. This structure allowed for a comparison across all 5 mAChR family members that revealed important differences in both orthosteric and allosteric sites that could inform the rational design of selective ligands. These structural studies, together with chimeric swaps between the extracellular regions of the M₂ and M₅ mAChRs, provided structural insight into kinetic selectivity, where ligands show differential residency times between related family members. Collectively, our study provides important insights into the nature of orthosteric and allosteric ligand interaction across the mAChR family that could be exploited for the design of selective drugs.

Mathematical modeling of the glucagon challenge test

A model for the homeostasis of glucose through the regulating hormones glucagon and insulin is described. It contains a subsystem that models the internalization of the glucagon receptor. Internalization is a mechanism in cell signaling, through which G-protein coupled receptors are taken from the surface of the cell to the endosome. The model is used to interpret data from a glucagon challenge test in which subjects have been under treatment with a novel glucagon receptor anti-sense drug which is aimed at reducing the number of receptors in the liver. It is shown how the receptor internalization results in tolerance of the blood glucose concentration to glucagon-induced hyperglycemia. We quantify the reduction of the number of receptors using the model and the data before and after treatment.

Internalization of Influenza Virus and Cell Surface Proteins Monitored by Site-Specific Conjugation of Protease-Sensitive Probes

Commonly used methods to monitor internalization of cell surface structures involve application of fluorescently or otherwise labeled antibodies against the target of interest. Genetic modification of the protein of interest, for example through creation of fusions with fluorescent or enzymatically active protein domains, is another approach to follow trafficking behavior. The former approach requires indirect methods, such as multiple rounds of cell staining, to distinguish between a target that remains surface-disposed and an internalized and/or recycled species. The latter approach necessitates the creation of fusions whose behavior may not accurately reflect that of their unmodified counterparts. Here, we report a method for the characterization of protein internalization in real time through sortase-mediated, site-specific labeling of single-domain antibodies or viral proteins with a newly developed, cathepsin-sensitive quenched-fluorophore probe. Quenched probes of this type have been used to measure enzyme activity in complex environments and for different cell types, but not as a sensor of protein movement into living cells. This approach allows a quantitative assessment of the movement of proteins into protease-containing endosomes in real time in living cells. We demonstrate considerable variation in the rate of endosomal delivery for different cell surface receptors. We were also able to characterize the kinetics of influenza virus delivery to cathepsin-positive compartments, showing highly coordinated arrival in endosomal compartments. This approach should be useful for identifying proteins expressed on cells of interest for targeted endosomal delivery of payloads, such as antibody-drug conjugates or antigens that require processing.

The steroid hormone 20-hydroxyecdysone binds to dopamine receptor to repress lepidopteran insect feeding and promote pupation

Holometabolous insects stop feeding at the final larval instar stage and then undergo metamorphosis; however, the mechanism is unclear. In the present study, using the serious lepidopteran agricultural pest Helicoverpa armigera as a model, we revealed that 20-hydroxyecdysone (20E) binds to the dopamine receptor (DopEcR), a G protein-coupled receptor, to stop larval feeding and promote pupation. DopEcR was expressed in various tissues and its level increased during metamorphic molting under 20E regulation. The 20E titer was low during larval feeding stages and high during wandering stages. By contrast, the dopamine (DA) titer was high during larval feeding stages and low during the wandering stages. Injection of 20E or blocking dopamine receptors using the inhibitor flupentixol decreased larval food consumption and body weight. Knockdown of DopEcR repressed larval feeding, growth, and pupation. 20E, via DopEcR, promoted apoptosis; and DA, via DopEcR, induced cell proliferation. 20E opposed DA function by repressing DA-induced cell proliferation and AKT phosphorylation. 20E, via DopEcR, induced gene expression and a rapid increase in intracellular calcium ions and cAMP. 20E induced the interaction of DopEcR with G proteins αs and αq. 20E, via DopEcR, induced protein phosphorylation and binding of the EcRB1-USP1 transcription complex to the ecdysone response element. DopEcR could bind 20E inside the cell membrane or after being isolated from the cell membrane. Mutation of DopEcR decreased 20E binding levels and related cellular responses. 20E competed with DA to bind to DopEcR. The results of the present study suggested that 20E, via binding to DopEcR, arrests larval feeding and promotes pupation.

The steroid hormone 20-hydroxyecdysone (20E) represses insect larval feeding and promotes metamorphosis; however, the mechanism is unclear. The dopamine receptor plays important roles in animal motor function and reward-motivated behavior. Using the serious lepidopteran agricultural pest Helicoverpa armigera as a model, we revealed that 20E binds to DopEcR to block the dopamine pathway and initiates the 20E pathway. Dopamine (DA) binds to the dopamine receptor (DopEcR), a G protein-coupled receptor (GPCR), to regulate cell proliferation, larval feeding, and growth. However, 20E competes with DA to bind to DopEcR, which represses larval feeding and triggers the 20E-pathway, leading to metamorphosis. The results suggested that 20E, via binding to DopEcR, stops larval feeding and promotes pupation, which presented an example of the steroid hormone regulating dopamine receptor and behavior. Our study showed that GPCRs can bind 20E and function as 20E cell membrane receptors.

Computational framework for predictive PBPK-PD-Tox simulations of opioids and antidotes

The primary goal of this work was to develop a computational tool to enable personalized prediction of pharmacological disposition and associated responses for opioids and antidotes. Here we present a computational framework for physiologically-based pharmacokinetic (PBPK) modeling of an opioid (morphine) and an antidote (naloxone). At present, the model is solely personalized according to an individual's mass. These PK models are integrated with a minimal pharmacodynamic model of respiratory depression induction (associated with opioid administration) and reversal (associated with antidote administration). The model was developed and validated on human data for IV administration of morphine and naloxone. The model can be further extended to consider different routes of administration, as well as to study different combinations of opioid receptor agonists and antagonists. This work provides the framework for a tool that could be used in model-based management of pain, pharmacological treatment of opioid addiction, appropriate use of antidotes for opioid overdose and evaluation of abuse deterrent formulations.

Short- and Medium-Term Exposure to Ocean Acidification Reduces Olfactory Sensitivity in Gilthead Seabream

The effects of ocean acidification on fish are only partially understood. Studies on olfaction are mostly limited to behavioral alterations of coral reef fish; studies on temperate species and/or with economic importance are scarce. The current study evaluated the effects of short- and medium-term exposure to ocean acidification on the olfactory system of gilthead seabream (Sparus aurata), and attempted to explain observed differences in sensitivity by changes in the protonation state of amino acid odorants. Short-term exposure to elevated PCO₂ decreased olfactory sensitivity to some odorants, such as L-serine, L-leucine, L-arginine, L-glutamate, and conspecific intestinal fluid, but not to others, such as L-glutamine and conspecific bile fluid. Seabream were unable to compensate for high PCO₂ levels in the medium term; after 4 weeks exposure to high PCO₂, the olfactory sensitivity remained lower in elevated PCO₂ water. The decrease in olfactory sensitivity in high PCO₂ water could be partly attributed to changes in the protonation state of the odorants and/or their receptor(s); we illustrate how protonation due to reduced pH causes changes in the charge distribution of odorant molecules, an essential component for ligand-receptor interaction. However, there are other mechanisms involved. At a histological level, the olfactory epithelium contained higher densities of mucus cells in fish kept in high CO₂ water, and a shift in pH of the mucus they produced to more neutral. These differences suggest a physiological response of the olfactory epithelium to lower pH and/or high CO₂ levels, but an inability to fully counteract the effects of acidification on olfactory sensitivity. Therefore, the current study provides evidence for a direct, medium term, global effect of ocean acidification on olfactory sensitivity in fish, and possibly other marine organisms, and suggests a partial explanatory mechanism.

A Metadynamics-Based Protocol for the Determination of GPCR-Ligand Binding Modes

G protein-coupled receptors (GPCRs) are a main drug target and therefore a hot topic in pharmaceutical research. One important prerequisite to understand how a certain ligand affects a GPCR is precise knowledge about its binding mode and the specific underlying interactions. If no crystal structure of the respective complex is available, computational methods can be used to deduce the binding site. One of them are metadynamics simulations which have the advantage of an enhanced sampling compared to conventional molecular dynamics simulations. However, the enhanced sampling of higher-energy states hampers identification of the preferred binding mode. Here, we present a novel protocol based on clustering of multiple walker metadynamics simulations which allows identifying the preferential binding mode from such conformational ensembles. We tested this strategy for three different model systems namely the histamine H₁ receptor in combination with its physiological ligand histamine, as well as the β2 adrenoceptor with its agonist adrenaline and its antagonist alprenolol. For all three systems, the proposed protocol was able to reproduce the correct binding mode known from the literature suggesting that the approach can more generally be applied to the prediction of GPCR ligand binding in future.

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Peeking at G-protein-coupled receptors through the molecular dynamics keyhole

Molecular dynamics is a state of the art computational tool for the investigation of biophysics phenomenon at a molecular scale, as it enables the modeling of dynamic processes, such as conformational motions, molecular solvation and ligand binding. The recent advances in structural biology have led to a bloom in published G-protein-coupled receptor structures, representing a solid and valuable resource for molecular dynamics studies. During the last decade, indeed, a plethora of physiological and pharmacological facets of this membrane protein superfamily have been addressed by means of molecular dynamics simulations, including the activation mechanism, allosterism and, very recently, biased signaling. Here, we try to recapitulate some of the main contributions that molecular dynamics has recently produced in the field.

Drug-receptor kinetics and sigma-1 receptor affinity differentiate clinically evaluated histamine H3 receptor antagonists

The histamine H3 receptor is a G protein-coupled receptor (GPCR) drug target that is highly expressed in the CNS, where it acts as both an auto- and hetero-receptor to regulate neurotransmission. As such, it has been considered as a relevant target in disorders as varied as Alzheimer's disease, schizophrenia, neuropathic pain and attention deficit hyperactivity disorder. A range of competitive antagonists/inverse agonists have progressed into clinical development, with pitolisant approved for the treatment of narcolepsy. Given the breadth of compounds developed and potential therapeutic indications, we assessed the comparative pharmacology of six investigational histamine H3 agents, including pitolisant, using native tissue and recombinant cells. Whilst all of the compounds tested displayed robust histamine H3 receptor inverse agonism and did not differentiate between the main H3 receptor splice variants, they displayed a wide range of affinities and kinetic properties, and included rapidly dissociating (pitolisant, S 38093-2, ABT-239) and slowly dissociating (GSK189254, JNJ-5207852, PF-3654746) agents. S 38093-2 had the lowest histamine H3 receptor affinity (pKB values 5.7-6.2), seemingly at odds with previously reported, potent in vivo activity in models of cognition. We show here that at pro-cognitive and anti-hyperalgesic/anti-allodynic doses, S 38093-2 preferentially occupies the mouse sigma-1 receptor in vivo, only engaging the histamine H3 receptor at doses associated with wakefulness promotion and neurotransmitter (histamine, ACh) release. Furthermore, pitolisant, ABT-239 and PF-3654746 also displayed appreciable sigma-1 receptor affinity, suggesting that this property differentiates clinically evaluated histamine H3 receptor antagonists and may play a role in their efficacy.

Pharmacological Assays for Investigating the NOP Receptor

The nociceptin/orphanin FQ (N/OFQ) peptide receptor (NOP) is a G protein-coupled receptor involved in the regulation of several physiological functions and pathological conditions. Thus, researchers from academia and industry are pursuing NOP to discover and study novel pharmacological entities. In a multidisciplinary effort of pharmacologists, medicinal chemists, and molecular and structural biologists the mechanisms of NOP activation and inhibition have been, at least partially, disentangled. Here, we review the in vitro methodologies employed, which have contributed to our understanding of this target. We hope this chapter guides the reader through the mostly established assay platforms to investigate NOP pharmacology, and gives some hints taking advantage from what has already illuminated the function of other GPCRs. We analyzed the pharmacological results obtained with a large panel of NOP ligands investigated in several assays including receptor binding, stimulation of GTPγS binding, decrease of cAMP levels, calcium flux stimulation via chimeric G proteins, NOP/G protein and NOP/β-arrestin interaction, label-free assays such as dynamic mass redistribution, and bioassays such as the electrically stimulated mouse vas deferens.

Identification of Chemotype Agonists for Human Resolvin D1 Receptor DRV1 with Pro-Resolving Functions

Resolution of acute inflammation is governed, in part, by specialized pro-resolving mediators, including lipoxins, resolvins, protectins, and maresins. Among them, resolvin D1 (RvD1) exhibits potent pro-resolving functions via activating human resolvin D1 receptor (DRV1/GPR32). RvD1 is a complex molecule that requires challenging organic synthesis, diminishing its potential as a therapeutic. Therefore, we implemented a high-throughput screening of small-molecule libraries and identified several chemotypes that activated recombinant DRV1, represented by NCGC00120943 (C1A), NCGC00135472 (C2A), pMPPF, and pMPPI. These chemotypes also elicited rapid impedance changes in cells overexpressing recombinant DRV1. With human macrophages, they each stimulated phagocytosis of serum-treated zymosan at concentrations comparable with that of RvD1, the endogenous DRV1 ligand. In addition, macrophage phagocytosis of live E. coli was significantly increased by these chemotypes in DRV1-transfected macrophages, compared with mock-transfected cells. Taken together, these chemotypes identified by unbiased screens act as RvD1 mimetics, exhibiting pro-resolving functions via interacting with human DRV1.

Binding of histamine to the H1 receptor-a molecular dynamics study

Binding of histamine to the G-protein coupled histamine H₁ receptor plays an important role in the context of allergic reactions; however, no crystal structure of the resulting complex is available yet. To deduce the histamine binding site, we performed unbiased molecular dynamics (MD) simulations on a microsecond time scale, which allowed to monitor one binding event, in which particularly the residues of the extracellular loop 2 were involved in the initial recognition process. The final histamine binding pose in the orthosteric pocket is characterized by interactions with Asp107^3.32, Tyr108^3.33, Thr194^5.43, Asn198^5.46, Trp428^6.48, Tyr431^6.51, Phe432^6.52, and Phe435^6.55, which is in agreement with existing mutational data. The conformational stability of the obtained complex structure was subsequently confirmed in 2 μs equilibrium MD simulations, and a metadynamics simulation proved that the detected binding site represents an energy minimum. A complementary investigation of a D107A mutant, which has experimentally been shown to abolish ligand binding, revealed that this exchange results in a significantly weaker interaction and enhanced ligand dynamics. This finding underlines the importance of the electrostatic interaction between the histamine ammonium group and the side chain of Asp107^3.32 for histamine binding.

Preassembled GPCR signaling complexes mediate distinct cellular responses to ultralow ligand concentrations

G protein-coupled receptors (GPCRs) are the largest class of cell surface signaling proteins, participate in nearly all physiological processes, and are the targets of 30% of marketed drugs. Typically, nanomolar to micromolar concentrations of ligand are used to activate GPCRs in experimental systems. We detected GPCR responses to a wide range of ligand concentrations, from attomolar to millimolar, by measuring GPCR-stimulated production of cyclic adenosine monophosphate (cAMP) with high spatial and temporal resolution. Mathematical modeling showed that femtomolar concentrations of ligand activated, on average, 40% of the cells in a population provided that a cell was activated by one to two binding events. Furthermore, activation of the endogenous β2-adrenergic receptor (β2AR) and muscarinic acetylcholine M3 receptor (M3R) by femtomolar concentrations of ligand in cell lines and human cardiac fibroblasts caused sustained increases in nuclear translocation of extracellular signal-regulated kinase (ERK) and cytosolic protein kinase C (PKC) activity, respectively. These responses were spatially and temporally distinct from those that occurred in response to higher concentrations of ligand and resulted in a distinct cellular proteomic profile. This highly sensitive signaling depended on the GPCRs forming preassembled, higher-order signaling complexes at the plasma membrane. Recognizing that GPCRs respond to ultralow concentrations of neurotransmitters and hormones challenges established paradigms of drug action and provides a previously unappreciated aspect of GPCR activation that is quite distinct from that typically observed with higher ligand concentrations.

Link between a high k on for drug binding and a fast clinical action: to be or not to be?

Review articles on binding kinetics essentially focus on drugs that dissociate slowly from their target since this is required for the successful treatment of many pathophysiological conditions. Recently, the therapeutic benefit of a high k on (i.e. the second order association rate constant) has also been linked to fast association and to a fast clinical action. Other studies, however, called this assertion into question since additional factors, like the dosing paradigm and the binding mechanism, are important as well. The still ongoing reticence about integrating binding kinetics in lead optimization programs motivated us to critically review the link between the drug's kinetic rate constants and their in vitro and in vivo target occupancy profile, with special focus on k on. The presented simulations tally with a positive link between a drug's effective/observed association rate (which is quite easy to determine in vitro) and the swiftness of its clinical action. On the other hand, the simulations show that the k on-concept should not be confounded with the effective association process since increasing this parameter only enhances the drug's in vitro and in vivo association under certain conditions: the binding mechanism should be suitable, rebinding (and thus the factors within the target's micro-environment that favour this mechanism) should not be too prominent and the dosage should not be kept in par with the drug's affinity. Otherwise, increasing k on could be ineffective or even be counter-productive.

Exploring a new ligand binding site of G protein-coupled receptors

Identifying a target ligand binding site is an important step for structure-based rational drug design as shown here for G protein-coupled receptors (GPCRs), which are among the most popular drug targets. We applied long-time scale molecular dynamics simulations, coupled with mutagenesis studies, to two prototypical GPCRs, the M3 and M4 muscarinic acetylcholine receptors. Our results indicate that unlike synthetic antagonists, which bind to the classic orthosteric site, the endogenous agonist acetylcholine is able to diffuse into a much deeper binding pocket. We also discovered that the most recently resolved crystal structure of the LTB4 receptor comprised a bound inverse agonist, which extended its benzamidine moiety to the same binding pocket discovered in this work. Analysis on all resolved GPCR crystal structures indicated that this new pocket could exist in most receptors. Our findings provide new opportunities for GPCR drug discovery.

Binding kinetics of cariprazine and aripiprazole at the dopamine D3 receptor

The dissociation behaviours of aripiprazole and cariprazine at the human D₂ and D₃ receptor are evaluated. A potential correlation between kinetics and in vivo profiles, especially cariprazine’s action on negative symptoms in schizophrenia, is investigated. The binding kinetics of four ligands were indirectly evaluated. After the receptor preparations were pre-incubated with the unlabelled ligands, the dissociation was initiated with an excess of [³H]spiperone. Slow dissociation kinetics characterizes aripiprazole and cariprazine at the D₂ receptor. At the D₃ receptor, aripiprazole exhibits a slow monophasic dissociation, while cariprazine displays a rapid biphasic behaviour. Functional ß-arrestin assays and molecular dynamics simulations at the D₃ receptor confirm a biphasic binding behaviour of cariprazine. This may influence its in vivo action, as the partial agonist could react rapidly to variations in the dopamine levels of schizophrenic patients and the ligand will not quantitatively dissociate from the receptor in one single step. With these findings novel agents may be developed that display rapid, biphasic dissociation from the D₃R to further investigate this effect on in vivo profiles.

Dynamic States of the Ligand-Free Class A G Protein-Coupled Receptor Extracellular Side

G protein-coupled receptors (GPCRs) make up the largest family of drug targets. The second extracellular loop (ECL2) and extracellular end of the third transmembrane helix (TM3) are basic structural elements of the GPCR ligand binding site. Currently, the disulfide bond between the two conserved cysteines in the ECL2 and TM3 is considered to be a basic GPCR structural feature. This disulfide bond has a significant effect on receptor dynamics and ligand binding. Here, molecular dynamics simulations and experimental results show that the two cysteines are distant from one another in the highest-population conformational state of ligand-free class A GPCRs and do not form a disulfide bond, indicating that the dynamics of the GPCR extracellular side are different from our conventional understanding. These surprising dynamics should have important effects on the drug binding process. On the basis of the two distinct ligand-free states, we suggest two kinetic processes for binding of ligands to GPCRs. These results challenge our commonly held beliefs regarding both GPCR structural features and ligand binding.

Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.