Search

Biased agonists, which selectively stimulate certain signaling pathways controlled by a G protein-coupled receptor (GPCR), hold great promise as drugs that maximize efficacy while minimizing dangerous side effects. Biased agonists of the μ-opioid receptor (μOR) are of particular interest as a means to achieve analgesia through G protein signaling without dose-limiting side effects such as respiratory depression and constipation. Rational structure-based design of biased agonists remains highly challenging, however, because the ligand-mediated interactions that are key to activation of each signaling pathway remain unclear. We identify several compounds for which the R- and S-enantiomers have distinct bias profiles at the μOR. These compounds serve as excellent comparative tools to study bias because the identical physicochemical properties of enantiomer pairs ensure that differences in bias profiles are due to differences in interactions with the μOR binding pocket. Atomic-level simulations of compounds at μOR indicate that R- and S-enantiomers adopt different poses that form distinct interactions with the binding pocket. A handful of specific interactions with highly conserved binding pocket residues appear to be responsible for substantial differences in arrestin recruitment between enantiomers. Our results offer guidance for rational design of biased agonists at μOR and possibly at related GPCRs.

The development of G protein-biased agonists for the μ-opioid receptor (MOR) offers a clear drug discovery rationale for improved analgesia and reduced side-effects of opiate pharmacotherapy. However, our understanding of the molecular mechanisms governing ligand bias is limited, which hinders our ability to rationally design biased compounds. We have investigated the role of MOR binding site residues W320 and Y328 in controlling bias, by receptor mutagenesis. The pharmacology of a panel of ligands in a cAMP and a β-arrestin2 assay were compared between the wildtype and mutated receptors, with bias factors calculated by operational analysis using ΔΔlog(τ/KA) values. [3H]diprenorphine competition binding was used to estimate affinity changes. Introducing the mutations W320A and Y328F caused changes in pathway bias, with different patterns of change between ligands. For example, DAMGO increased relative β-arrestin2 activity at the W320A mutant, whilst its β-arrestin2 response was completely lost at Y328F. In contrast, endomorphin-1 gained activity with Y328F but lost activity at W320A, in both pathways. For endomorphin-2 there was a directional shift from cAMP bias at the wildtype towards more β-arrestin2 bias at W320A. We also observe clear uncoupling between mutation-driven changes in function and binding affinity. These findings suggest that the mutations influenced the balance of pathway activation in a ligand-specific manner, thus identifying residues in the MOR binding pocket that govern ligand bias. This increases our understanding of how ligand/receptor binding interactions can be translated into agonist-specific pathway activation.

Diseases such as chronic pain with complex etiologies are unlikely to respond to single, target-specific therapeutics but rather require intervention at multiple points within a perturbed disease system. Such approaches are being enabled by the rise of computational methods to identify key points of intervention and by new screening techniques that focus on a relevant condition or phenotype, rather than a specific target. Here we apply an in silico network pharmacology approach to identify small-molecule compounds with the potential to selectively disrupt the structure of a chronic-pain specific disease network, which we validate using a novel phenotypic screen that recapitulates key aspects of neuronal and pain biology by measuring changes in neuronal excitability in native sensory neurons. The combination of network pharmacology with a phenotypic screen is a powerful approach; we show that hit rates increase from 26% to 42%. This represents a rational approach to the discovery of compounds with a poly-pharmacology based therapeutic value, which will be vital for the discovery of treatments for complex disease.

Summary Small conductance potassium (SK) ion channels define neuronal firing rates by conducting the after-hyperpolarization current. They are key targets in developing therapies where neuronal firing rates are dysfunctional, such as in epilepsy, Parkinson's, and amyotrophic lateral sclerosis (ALS). Here, we characterize a binding pocket situated at the intracellular interface of SK2 and calmodulin, which we show to be shared by multiple small-molecule chemotypes. Crystallization of this complex revealed that riluzole (approved for ALS) and an analog of the anti-ataxic agent (4-chloro-phenyl)-[2-(3,5-dimethyl-pyrazol-1-yl)-pyrimidin-4-yl]-amine (CyPPA) bind to and allosterically modulate via this site. Solution-state nuclear magnetic resonance demonstrates that riluzole, NS309, and CyPPA analogs bind at this bipartite pocket. We demonstrate, by patch-clamp electrophysiology, that both classes of ligand interact with overlapping but distinct residues within this pocket. These data define a clinically important site, laying the foundations for further studies of the mechanism of action of riluzole and related molecules.

The discovery and selection of a highly potent and selective Na

Selective modulators of the γ-amino butyric acid (GABA

Potassium channels are of paramount physiological and pathological importance and therefore constitute significant drug targets. One of the keys to rationalize the way drugs modulate ion channels is to understand the ability of such small molecules to access their respective binding sites, from which they can exert an activating or inhibitory effect. Many computational studies have probed the energetics of ion permeation, and the mechanisms of voltage gating, but little is known about the role of fenestrations as possible mediators of drug entry in potassium channels. To explore the existence, structure, and conformational dynamics of transmembrane fenestrations accessible by drugs in potassium channels, molecular dynamics simulation trajectories were analyzed from three potassium channels: the open state voltage-gated channel Kv1.2, the G protein-gated inward rectifying channel GIRK2 (Kir3.2), and the human two-pore domain TWIK-1 (K2P1.1). The main results of this work were the identification of the sequence identity of four main lateral fenestrations of similar length and with bottleneck radius in the range of 0.9-2.4 Å for this set of potassium channels. It was found that the fenestrations in Kv1.2 and Kir3.2 remain closed to the passage of molecules larger than water. In contrast, in the TWIK-1 channel, both open and closed fenestrations are sampled throughout the simulation, with bottleneck radius shown to correlate with the random entry of lipid membrane molecules into the aperture of the fenestrations. Druggability scoring function analysis of the fenestration regions suggests that Kv and Kir channels studied are not druggable in practice due to steric constraining of the fenestration bottleneck. A high (>50%) fenestration sequence identity was found in each potassium channel subfamily studied, Kv1, Kir3, and K2P1. Finally, the reported fenestration sequence of TWIK-1 compared favorably with another channel, K2P channel TREK-2, reported to possess open fenestrations, suggesting that K2P channels could be druggable via fenestrations, for which we reported atomistic detail of the fenestration region, including the flexible residues M260 and L264 that interact with POPC membrane in a concerted fashion with the aperture and closure of the fenestrations.

Supplemental Digital Content is available in the text., Background: Elabela/toddler (ELA) is a critical cardiac developmental peptide that acts through the G-protein–coupled apelin receptor, despite lack of sequence similarity to the established ligand apelin. Our aim was to investigate the receptor pharmacology, expression pattern, and in vivo function of ELA peptides in the adult cardiovascular system, to seek evidence for alteration in pulmonary arterial hypertension (PAH) in which apelin signaling is downregulated, and to demonstrate attenuation of PAH severity with exogenous administration of ELA in a rat model. Methods: In silico docking analysis, competition binding experiments, and downstream assays were used to characterize ELA receptor binding in human heart and signaling in cells expressing the apelin receptor. ELA expression in human cardiovascular tissues and plasma was determined using real-time quantitative polymerase chain reaction, dual-labeling immunofluorescent staining, and immunoassays. Acute cardiac effects of ELA-32 and [Pyr1]apelin-13 were assessed by MRI and cardiac catheterization in anesthetized rats. Cardiopulmonary human and rat tissues from PAH patients and monocrotaline- and Sugen/hypoxia-exposed rats were used to show changes in ELA expression in PAH. The effect of ELA treatment on cardiopulmonary remodeling in PAH was investigated in the monocrotaline rat model. Results: ELA competed for binding of apelin in human heart with overlap for the 2 peptides indicated by in silico modeling. ELA activated G-protein– and β-arrestin–dependent pathways. We detected ELA expression in human vascular endothelium and plasma. Comparable to apelin, ELA increased cardiac contractility, ejection fraction, and cardiac output and elicited vasodilatation in rat in vivo. ELA expression was reduced in cardiopulmonary tissues from PAH patients and PAH rat models, respectively. ELA treatment significantly attenuated elevation of right ventricular systolic pressure and right ventricular hypertrophy and pulmonary vascular remodeling in monocrotaline-exposed rats. Conclusions: These results show that ELA is an endogenous agonist of the human apelin receptor, exhibits a cardiovascular profile comparable to apelin, and is downregulated in human disease and rodent PAH models, and exogenous peptide can reduce the severity of cardiopulmonary remodeling and function in PAH in rats. This study provides additional proof of principle that an apelin receptor agonist may be of therapeutic use in PAH in humans.

Supplemental Digital Content is available in the text., [Pyr1]apelin-13 is an endogenous vasodilator and inotrope but is downregulated in pulmonary hypertension and heart failure, making the apelin receptor an attractive therapeutic target. Agonists acting at the same G-protein–coupled receptor can be engineered to stabilize different conformational states and function as biased ligands, selectively stimulating either G-protein or β-arrestin pathways. We used molecular dynamics simulations of apelin/receptor interactions to design cyclic analogues and identified MM07 as a biased agonist. In β-arrestin and internalization assays (G-protein–independent), MM07 was 2 orders of magnitude less potent than [Pyr1]apelin-13. In a G-protein–dependent saphenous vein contraction assay, both peptides had comparable potency (pD2:[Pyr1]apelin-13 9.93±0.24; MM07 9.54±0.42) and maximum responses with a resulting bias for MM07 of ≈350- to 1300-fold for the G-protein pathway. In rats, systemic infusions of MM07 (10-100nmol) caused a dose-dependent increase in cardiac output that was significantly greater than the response to [Pyr1]apelin-13. Similarly, in human volunteers, MM07 produced a significant dose-dependent increase in forearm blood flow with a maximum dilatation double that is seen with [Pyr1]apelin-13. Additionally, repeated doses of MM07 produced reproducible increases in forearm blood flow. These responses are consistent with a more efficacious action of the biased agonist. In human hand vein, both peptides reversed an established norepinephrine constrictor response and significantly increased venous flow. Our results suggest that MM07 acting as a biased agonist at the apelin receptor can preferentially stimulate the G-protein pathway, which could translate to improved efficacy in the clinic by selectively stimulating vasodilatation and inotropic actions but avoiding activating detrimental β-arrestin–dependent pathways.

Metabolism of xenobiotic and endogenous compounds is frequently complex, not completely elucidated, and therefore often ambiguous. The prediction of sites of metabolism (SoM) can be particularly helpful as a first step toward the identification of metabolites, a process especially relevant to drug discovery. This paper describes a reactivity approach for predicting SoM whereby reactivity is derived directly from the ground state ligand molecular orbital analysis, calculated using Density Functional Theory, using a novel implementation of the average local ionization energy. Thus each potential SoM is sampled in the context of the whole ligand, in contrast to other popular approaches where activation energies are calculated for a predefined database of molecular fragments and assigned to matching moieties in a query ligand. In addition, one of the first descriptions of molecular dynamics of cytochrome P450 (CYP) isoforms 3A4, 2D6, and 2C9 in their Compound I state is reported, and, from the representative protein structures obtained, an analysis and evaluation of various docking approaches using GOLD is performed. In particular, a covalent docking approach is described coupled with the modeling of important electrostatic interactions between CYP and ligand using spherical constraints. Combining the docking and reactivity results, obtained using standard functionality from common docking and quantum chemical applications, enables a SoM to be identified in the top 2 predictions for 75%, 80%, and 78% of the data sets for 3A4, 2D6, and 2C9, respectively, results that are accessible and competitive with other recently published prediction tools.

Fibroblast growth factors (FGFs) are recognized targets for the development of therapies against angiogenesis-driven diseases, including cancer. The formation of a ternary complex with the transmembrane tyrosine kinase receptors (FGFRs), and heparan sulphate proteoglycans (HSPGs) is required for FGF2 pro-angiogenic activity. Here by using a combination of techniques including Nuclear Magnetic Resonance, Molecular Dynamics, Surface Plasmon Resonance and cell-based binding assays we clarify the molecular mechanism of inhibition of an angiostatic small molecule, sm27, mimicking the endogenous inhibitor of angiogenesis, thrombospondin-1. NMR and MD data demonstrate that sm27 engages the heparin-binding site of FGF2 and induces long-range dynamics perturbations along FGF2/FGFR1 interface regions. The functional consequence of the inhibitor binding is an impaired FGF2 interaction with both its receptors, as demonstrated by SPR and cell-based binding assays. We propose that sm27 antiangiogenic activity is based on a twofold-direct and allosteric-mechanism, inhibiting FGF2 binding to both its receptors.