Papyrus: a large-scale curated dataset aimed at bioactivity predictions.

With the ongoing rapid growth of publicly available ligand-protein bioactivity data, there is a trove of valuable data that can be used to train a plethora of machine-learning algorithms. However, not all data is equal in terms of size and quality and a significant portion of researchers' time is needed to adapt the data to their needs. On top of that, finding the right data for a research question can often be a challenge on its own. To meet these challenges, we have constructed the Papyrus dataset. Papyrus is comprised of around 60 million data points. This dataset contains multiple large publicly available datasets such as ChEMBL and ExCAPE-DB combined with several smaller datasets containing high-quality data. The aggregated data has been standardised and normalised in a manner that is suitable for machine learning. We show how data can be filtered in a variety of ways and also perform some examples of quantitative structure-activity relationship analyses and proteochemometric modelling. Our ambition is that this pruned data collection constitutes a benchmark set that can be used for constructing predictive models, while also providing an accessible data source for research.

Impact of cancer-associated mutations in CC chemokine receptor 2 on receptor function and antagonism.

CC chemokine receptor 2 (CCR2), a G protein-coupled receptor, plays a role in many cancer-related processes such as metastasis formation and immunosuppression. Since âˆ¼ 20 % of human cancers contain mutations in G protein-coupled receptors, ten cancer-associated CCR2 mutants obtained from the Genome Data Commons were investigated for their effect on receptor functionality and antagonist binding. Mutations were selected based on either their vicinity to CCR2's orthosteric or allosteric binding sites or their presence in conserved amino acid motifs. One of the mutant receptors, namely S101P2.63 with a mutation near the orthosteric binding site, did not express on the cell surface. All other studied mutants showed a decrease in or a lack of G protein activation in response to the main endogenous CCR2 ligand CCL2, but no change in potency was observed. Furthermore, INCB3344 and LUF7482 were chosen as representative orthosteric and allosteric antagonists, respectively. No change in potency was observed in a functional assay, but mutations located at F1163.28 impacted orthosteric antagonist binding significantly, while allosteric antagonist binding was abolished for L134Q3.46 and D137N3.49 mutants. As CC chemokine receptor 2 is an attractive drug target in cancer, the negative effect of these mutations on receptor functionality and drugability should be considered in the drug discovery process.

Pan-cancer functional analysis of somatic mutations in G protein-coupled receptors.

G Protein-coupled receptors (GPCRs) are the most frequently exploited drug target family, moreover they are often found mutated in cancer. Here we used a dataset of mutations found in patient samples derived from the Genomic Data Commons and compared it to the natural human variance as exemplified by data from the 1000 genomes project. We explored cancer-related mutation patterns in all GPCR classes combined and individually. While the location of the mutations across the protein domains did not differ significantly in the two datasets, a mutation enrichment in cancer patients was observed among class-specific conserved motifs in GPCRs such as the Class A "DRY" motif. A Two-Entropy Analysis confirmed the correlation between residue conservation and cancer-related mutation frequency. We subsequently created a ranking of high scoring GPCRs, using a multi-objective approach (Pareto Front Ranking). Our approach was confirmed by re-discovery of established cancer targets such as the LPA and mGlu receptor families, but also discovered novel GPCRs which had not been linked to cancer before such as the P2Y Receptor 10 (P2RY10). Overall, this study presents a list of GPCRs that are amenable to experimental follow up to elucidate their role in cancer.

Correlation between human ether-a-go-go-related gene channel inhibition and action potential prolongation.

BACKGROUND AND PURPOSE: Human ether-a-go-go-related gene (hERG; Kv 11.1) channel inhibition is a widely accepted predictor of cardiac arrhythmia. hERG channel inhibition alone is often insufficient to predict pro-arrhythmic drug effects. This study used a library of dofetilide derivatives to investigate the relationship between standard measures of hERG current block in an expression system and changes in action potential duration (APD) in human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). The interference from accompanying block of Cav 1.2 and Nav 1.5 channels was investigated along with an in silico AP model. EXPERIMENTAL APPROACH: Drug-induced changes in APD were assessed in hiPSC-CMs using voltage-sensitive dyes. The IC50 values for dofetilide and 13 derivatives on hERG current were estimated in an HEK293 expression system. The relative potency of each drug on APD was estimated by calculating the dose (D150 ) required to prolong the APD at 90% (APD90 ) repolarization by 50%. KEY RESULTS: The D150 in hiPSC-CMs was linearly correlated with IC50 of hERG current. In silico simulations supported this finding. Three derivatives inhibited hERG without prolonging APD, and these compounds also inhibited Cav 1.2 and/or Nav 1.5 in a channel state-dependent manner. Adding Cav 1.2 and Nav 1.2 block to the in silico model recapitulated the direction but not the extent of the APD change. CONCLUSIONS AND IMPLICATIONS: Potency of hERG current inhibition correlates linearly with an index of APD in hiPSC-CMs. The compounds that do not correlate have additional effects including concomitant block of Cav 1.2 and/or Nav 1.5 channels. In silico simulations of hiPSC-CMs APs confirm the principle of the multiple ion channel effects.

Small molecule absorption by PDMS in the context of drug response bioassays.

The polymer polydimethylsiloxane (PDMS) is widely used to build microfluidic devices compatible with cell culture. Whilst convenient in manufacture, PDMS has the disadvantage that it can absorb small molecules such as drugs. In microfluidic devices like "Organs-on-Chip", designed to examine cell behavior and test the effects of drugs, this might impact drug bioavailability. Here we developed an assay to compare the absorption of a test set of four cardiac drugs by PDMS based on measuring the residual non-absorbed compound by High Pressure Liquid Chromatography (HPLC). We showed that absorption was variable and time dependent and not determined exclusively by hydrophobicity as claimed previously. We demonstrated that two commercially available lipophilic coatings and the presence of cells affected absorption. The use of lipophilic coatings may be useful in preventing small molecule absorption by PDMS.

Kinetic binding and activation profiles of endogenous tachykinins targeting the NK1 receptor.

Ligand-receptor binding kinetics (i.e. association and dissociation rates) are emerging as important parameters for drug efficacy in vivo. Awareness of the kinetic behavior of endogenous ligands is pivotal, as drugs often have to compete with those. The binding kinetics of neurokinin 1 (NK1) receptor antagonists have been widely investigated while binding kinetics of endogenous tachykinins have hardly been reported, if at all. Therefore, the aim of this research was to investigate the binding kinetics of endogenous tachykinins and derivatives thereof and their role in the activation of the NK1 receptor. We determined the binding kinetics of seven tachykinins targeting the NK1 receptor. Dissociation rate constants (koff) ranged from 0.026±0.0029min-1 (Sar9,Met(O2)11-SP) to 0.21±0.015min-1 (septide). Association rate constants (kon) were more diverse: substance P (SP) associated the fastest with a kon value of 0.24±0.046nM-1min-1 while neurokinin A (NKA) had the slowest association rate constant of 0.001±0.0002nM-1min-1. Kinetic binding parameters were highly correlated with potency and maximal response values determined in label-free impedance-based experiments on U-251 MG cells. Our research demonstrates large variations in binding kinetics of tachykinins which correlate to receptor activation. These findings provide new insights into the ligand-receptor interactions of tachykinins and underline the importance of measuring binding kinetics of both drug candidates and competing endogenous ligands.

Getting personal: Endogenous adenosine receptor signaling in lymphoblastoid cell lines.

Genetic differences between individuals that affect drug action form a challenge in drug therapy. Many drugs target G protein-coupled receptors (GPCRs), and a number of receptor variants have been noted to impact drug efficacy. This, however, has never been addressed in a systematic way, and, hence, we studied real-life genetic variation of receptor function in personalized cell lines. As a showcase we studied adenosine A2A receptor (A2AR) signaling in lymphoblastoid cell lines (LCLs) derived from a family of four from the Netherlands Twin Register (NTR), using a non-invasive label-free cellular assay. The potency of a partial agonist differed significantly for one individual. Genotype comparison revealed differences in two intron SNPs including rs2236624, which has been associated with caffeine-induced sleep disorders. While further validation is needed to confirm genotype-specific effects, this set-up clearly demonstrated that LCLs are a suitable model system to study genetic influences on A2AR response in particular and GPCR responses in general.

Persistent GnRH receptor activation in pituitary αT3-1 cells analyzed with a label-free technology.

The gonadotropin-releasing hormone (GnRH) receptor is a drug target for certain hormone-dependent diseases such as prostate cancer. In this study, we examined the activation profiles of the endogenous ligand, GnRH and a well-known marketed analog, buserelin using a label-free assay in pituitary αT3-1 cells with endogenous GnRH receptor expression. This whole cell impedance-based technology allows for the real-time measurement of morphological cellular changes. Both agonists dose-dependently decreased the impedance as a result of GnRH receptor activation with potencies of 9.3 ± 0.1 (pEC50 value, buserelin) and 7.8 ± 0.06 (pEC50 value, GnRH). Subsequently, GnRH receptor activation was completely abolished with a selective Gαq inhibitor, thereby confirming the Gαq-coupling of the GnRH receptor in pituitary αT3-1 cells. Additionally, we observed continued responses after agonist stimulation of αT3-1 cells indicating long-lasting cellular effects. Wash-out experiments demonstrated that the long-lasting effects induced by GnRH were most likely caused by rebinding since over 70% of the original response was abolished after wash-out. In contrast, a long receptor residence time was responsible for the prolonged effects caused by buserelin, with over 70% of the original response remaining after wash-out. In summary, we validated that impedance-based label-free technology is suited for studying receptor-mediated activation in cell lines endogenously expressing the target of interest. Moreover, this real-time monitoring allows the examination of binding kinetics and its influence on receptor activation at a cellular level.

Mass spectrometry-based ligand binding assays on adenosine A1 and A2A receptors.

Conventional methods to measure ligand-receptor binding parameters typically require radiolabeled ligands as probes. Despite the robustness of radioligand binding assays, they carry inherent disadvantages in terms of safety precautions, expensive synthesis, special lab requirements, and waste disposal. Mass spectrometry (MS) is a method that can selectively detect ligands without the need of a label. The sensitivity of MS equipment increases progressively, and currently, it is possible to detect low ligand quantities that are usually found in ligand binding assays. We developed a label-free MS ligand binding (MS binding) assay on the adenosine A(1) and A(2A) receptors (A(1)AR and A(2A)AR), which are well-characterized members of the class A G protein-coupled receptor (GPCR) family. Radioligand binding assays for both receptors are well established, and ample data is available to compare and evaluate the performance of an MS binding assay. 1,3-Dipropyl-8-cyclopentyl-xanthine (DPCPX) and 4-(2-((7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a]-[1,3,5]triazin-5-yl)amino)ethyl)phenol (ZM-241,385) are high-affinity ligands selective for the A(1)AR and A(2A)AR, respectively. To proof the feasibility of MS binding on the A(1)AR and A(2A)AR, we first developed an MS detection method for unlabeled DPCPX and ZM-241,385. To serve as internal standards, both compounds were also deuterium-labeled. Subsequently, we investigated whether the two unlabeled compounds could substitute for their radiolabeled counterparts as marker ligands in binding experiments, including saturation, displacement, dissociation, and competition association assays. Furthermore, we investigated the accuracy of these assays if the use of internal standards was excluded. The results demonstrate the feasibility of the MS binding assay, even in the absence of a deuterium-labeled internal standard, and provide great promise for the further development of label-free assays based on MS for other GPCRs.

Optogenetic activation of intracellular adenosine A2A receptor signaling in the hippocampus is sufficient to trigger CREB phosphorylation and impair memory.

Human and animal studies have converged to suggest that caffeine consumption prevents memory deficits in aging and Alzheimer's disease through the antagonism of adenosine A2A receptors (A2ARs). To test if A2AR activation in the hippocampus is actually sufficient to impair memory function and to begin elucidating the intracellular pathways operated by A2AR, we have developed a chimeric rhodopsin-A2AR protein (optoA2AR), which retains the extracellular and transmembrane domains of rhodopsin (conferring light responsiveness and eliminating adenosine-binding pockets) fused to the intracellular loop of A2AR to confer specific A2AR signaling. The specificity of the optoA2AR signaling was confirmed by light-induced selective enhancement of cAMP and phospho-mitogen-activated protein kinase (p-MAPK) (but not cGMP) levels in human embryonic kidney 293 (HEK293) cells, which was abolished by a point mutation at the C terminal of A2AR. Supporting its physiological relevance, optoA2AR activation and the A2AR agonist CGS21680 produced similar activation of cAMP and p-MAPK signaling in HEK293 cells, of p-MAPK in the nucleus accumbens and of c-Fos/phosphorylated-CREB (p-CREB) in the hippocampus, and similarly enhanced long-term potentiation in the hippocampus. Remarkably, optoA2AR activation triggered a preferential p-CREB signaling in the hippocampus and impaired spatial memory performance, while optoA2AR activation in the nucleus accumbens triggered MAPK signaling and modulated locomotor activity. This shows that the recruitment of intracellular A2AR signaling in the hippocampus is sufficient to trigger memory dysfunction. Furthermore, the demonstration that the biased A2AR signaling and functions depend on intracellular A2AR loops prompts the possibility of targeting the intracellular A2AR-interacting partners to selectively control different neuropsychiatric behaviors.

Kv 11.1 (hERG)-induced cardiotoxicity: a molecular insight from a binding kinetics study of prototypical Kv 11.1 (hERG) inhibitors.

BACKGROUND AND PURPOSE: Drug-induced arrhythmia due to blockade of the Kv 11.1 channel (also known as the hERG K(+) channel) is a frequent side effect. Previous studies have primarily focused on equilibrium parameters, i.e. affinity or potency, of drug candidates at the channel. The aim of this study was to determine the kinetics of the interaction with the channel for a number of known Kv 11.1 blockers and to explore a possible correlation with the affinity or physicochemical properties of these compounds. EXPERIMENTAL APPROACH: The affinity and kinetic parameters of 15 prototypical Kv 11.1 inhibitors were evaluated in a number of [(3) H]-dofetilide binding assays. The lipophilicity (logKW - C8 ) and membrane partitioning (logKW - IAM ) of these compounds were determined by means of HPLC analysis. KEY RESULTS: A novel [(3) H]-dofetilide competition association assay was set up and validated, which allowed us to determine the binding kinetics of the Kv 11.1 blockers used in this study. Interestingly, the compounds' affinities (Ki values) were correlated to their association rates rather than dissociation rates. Overall lipophilicity or membrane partitioning of the compounds were not correlated to their affinity or rate constants for the channel. CONCLUSIONS AND IMPLICATIONS: A compound's affinity for the Kv 11.1 channel is determined by its rate of association with the channel, while overall lipophilicity and membrane affinity are not. In more general terms, our findings provide novel insights into the mechanism of action for a compound's activity at the Kv 11.1 channel. This may help to elucidate how Kv 11.1-induced cardiotoxicity is governed and how it can be circumvented in the future.

Design and synthesis of novel small molecule CCR2 antagonists: evaluation of 4-aminopiperidine derivatives.

A novel N-(2-oxo-2-(piperidin-4-ylamino)ethyl)-3-(trifluoromethyl)benzamide series of human CCR2 chemokine receptor antagonists was identified. With a pharmacophore model based on known CCR2 antagonists a new core scaffold was designed, analogues of it synthesized and structure­affinity relationship studies derived yielding a new high affinity CCR2 antagonist N-(2-((1-(4-(3-methoxyphenyl)cyclohexyl)piperidin-4-yl)amino)-2-oxoethyl)-3-(trifluoromethyl)benzamide.

Structure-activity relationships of pentamidine-affected ion channel trafficking and dofetilide mediated rescue.

BACKGROUND AND PURPOSE: Drug interference with normal hERG protein trafficking substantially reduces the channel density in the plasma membrane and thereby poses an arrhythmic threat. The chemical substructures important for hERG trafficking inhibition were investigated using pentamidine as a model drug. Furthermore, the relationship between acute ion channel block and correction of trafficking by dofetilide was studied. EXPERIMENTAL APPROACH: hERG and K(IR)2.1 trafficking in HEK293 cells was evaluated by Western blot and immunofluorescence microscopy after treatment with pentamidine and six pentamidine analogues, and correction with dofetilide and four dofetilide analogues that displayed different abilities to inhibit IKr . Molecular dynamics simulations were used to address mode, number and type of interactions between hERG and dofetilide analogues. KEY RESULTS: Structural modifications of pentamidine differentially affected plasma membrane levels of hERG and K(IR)2.1. Modification of the phenyl ring or substituents directly attached to it had the largest effect, affirming the importance of these chemical residues in ion channel binding. PA-4 had the mildest effects on both ion channels. Dofetilide corrected pentamidine-induced hERG, but not K(IR)2.1 trafficking defects. Dofetilide analogues that displayed high channel affinity, mediated by pi-pi stacks and hydrophobic interactions, also restored hERG protein levels, whereas analogues with low affinity were ineffective. CONCLUSIONS AND IMPLICATIONS: Drug-induced trafficking defects can be minimized if certain chemical features are avoided or 'synthesized out'; this could influence the design and development of future drugs. Further analysis of such features in hERG trafficking correctors may facilitate the design of a non-blocking corrector for trafficking defective hERG proteins in both congenital and acquired LQTS.

The role of the second and third extracellular loops of the adenosine A1 receptor in activation and allosteric modulation.

The adenosine A1 receptor is a member of the large membrane protein family that signals through G proteins, the G protein-coupled receptors (GPCRs). GPCRs consist of seven transmembrane domains connected by three intracellular and three extracellular loops. Their N-terminus is extracellular, the C-terminal tail is in the cytoplasm. The transmembrane domains in receptor subfamilies that bind the same endogenous ligand, such as dopamine or adenosine, tend to be highly similar. In contrast, the loop regions can vary greatly, both in sequence and in length, and the role these loops have in the activation mechanism of the receptors remains unclear. Here, we investigated the activating role of the second and third extracellular loop of the human adenosine A1 receptor. By means of an (Ala)3 mutagenic scan in which consecutive sets of three amino acids were mutated into alanine residues in EL2 and a classical alanine scan in EL3, we revealed a strong regulatory role for the second extracellular loop (EL2) of the human adenosine A1 receptor. Besides many residues in the second and the third extracellular loops important for adenosine A1 receptor activation, we also identified two residues in EL2, a tryptophan and a glutamate, that affect the influence of the allosteric modulator PD81,723. These results, combined with a comparison of the different receptor loop regions, provide insight in the activation mechanism of this typical class A GPCR and further emphasize the unique pharmacological profile the loops can provide to individual receptors, even within subfamilies of GPCRs.

Structure-activity relationships of trans-substituted-propenoic acid derivatives on the nicotinic acid receptor HCA2 (GPR109A).

Nicotinic acid (niacin) has been used for decades as an antidyslipidemic drug in man. Its main target is the hydroxy-carboxylic acid receptor HCA2 (GPR109A), a G protein-coupled receptor. Other acids and esters such as methyl fumarate also interact with the receptor, which constituted the basis for the current study. We synthesized a novel series of substituted propenoic acids, such as fumaric acid esters, fumaric acid amides and cinnamic acid derivatives, and determined their affinities for the HCA2 receptor. We observed a rather restricted binding pocket on the receptor with trans-cinnamic acid being the largest planar ligand in our series with appreciable affinity for the receptor. Molecular modeling and analysis of the structure-activity relationships in the series suggest a planar trans-propenoic acid pharmacophore with a maximum length of 8 Å and out-of-plane orientation of the larger substituents.

Importance of the extracellular loops in G protein-coupled receptors for ligand recognition and receptor activation.

G protein-coupled receptors (GPCRs) are the major drug target of medicines on the market today. Therefore, much research is and has been devoted to the elucidation of the function and three-dimensional structure of this large family of membrane proteins, which includes multiple conserved transmembrane domains connected by intra- and extracellular loops. In the last few years, the less conserved extracellular loops have garnered increasing interest, particularly after the publication of several GPCR crystal structures that clearly show the extracellular loops to be involved in ligand binding. This review will summarize the recent progress made in the clarification of the ligand binding and activation mechanism of class-A GPCRs and the role of extracellular loops in this process.

Interleukin-6 upregulates neuronal adenosine A1 receptors: implications for neuromodulation and neuroprotection.

The immunological response in the brain is crucial to overcome neuropathological events. Some inflammatory mediators, such as the immunoregulatory cytokine interleukin-6 (IL-6) affect neuromodulation and may also play protective roles against various noxious conditions. However, the fundamental mechanisms underlying the long-term effects of IL-6 in the brain remain unclear. We now report that IL-6 increases the expression and function of the neuronal adenosine A1 receptor, with relevant consequences to synaptic transmission and neuroprotection. IL-6-induced amplification of A1 receptor function enhances the responses to readily released adenosine during hypoxia, enables neuronal rescue from glutamate-induced death, and protects animals from chemically induced convulsing seizures. Taken together, these results suggest that IL-6 minimizes the consequences of excitotoxic episodes on brain function through the enhancement of endogenous adenosinergic signaling.

Population pharmacokinetic modeling of blood-brain barrier transport of synthetic adenosine A1 receptor agonists.

A population pharmacokinetic model is proposed for estimation of the brain distribution clearance of synthetic A1 receptor agonists in vivo. Rats with permanent venous and arterial cannulas in combination with a microdialysis probe in the striatum received intravenous infusions of 8-methylamino-N6-cyclopentyladenosine (MCPA) and 2'-deoxyribose-N6-cyclopentyladenosine (2'-dCPA) (10 mg kg(-1)). The clearance for transport from blood to the brain was estimated by simultaneous analysis of the blood and extracellular fluid concentrations using a compartmental pharmacokinetic model. The proposed pharmacokinetic model consists of three compartments describing the time course of the concentration in blood in combination with three compartments for the brain extracellular fluid concentrations. The blood clearance was 7.4 +/- 0.5 for MCPA and 7.2 +/- 1.4 ml min(-1) for 2'-dCPA. The in vivo microdialysis recoveries determined by the dynamic-no-net-flux method were independent of time with values of 0.21 +/- 0.02 and 0.22 +/- 0.01 for MCPA and 2'-dCPA, respectively. The values of the intercompartmental clearance for the distribution from blood to brain were 1.9 +/- 0.4 versus 1.6 +/- 0.3 mul min(-1) for MCPA and 2'-dCPA, respectively. It is concluded that on basis of the novel six-compartment model precise estimates of the rate of brain distribution are obtained that are independent of eventual differences in systemic exposure. The low brain distribution rates of MCPA and 2'-dCPA were consistent with in vitro tests. Furthermore, a slow elimination from the brain compartment was observed, indicating that the duration of central nervous system effects may be much longer than expected on the basis of the terminal half-life in blood.

Pyrazole derivatives as partial agonists for the nicotinic acid receptor.

Nicotinic acid as a hypolipidemic agent appears unique due to its potential to increase HDL cholesterol levels to a greater extent than other drugs. However, it has some side effects, among which severe skin flushing is the most frequent and often limits patients' compliance. In a search for novel agonists for the recently identified and cloned G protein-coupled nicotinic acid receptor, we synthesized a series of substituted pyrazole-3-carboxylic acids that proved to have substantial affinity for this receptor. The affinities were measured by inhibition of [(3)H]nicotinic acid binding to rat spleen membranes. Potencies and intrinsic activities relative to nicotinic acid were determined by their effects on [(35)S]GTPgammaS binding to rat adipocyte and spleen membranes. Interestingly, most compounds were partial agonists. In particular, 2-diazabicyclo[3,3,0(4,8)]octa-3,8-diene-3-carboxylic acid (4c) and 5-propylpyrazole-3-carboxylic acid (4f) proved active with K(i) values of approximately 0.15 microM and EC(50) values of approximately 6 microM, while their intrinsic activity was only approximately 50% when compared to nicotinic acid. Even slightly more active was 5-butylpyrazole-3-carboxylic acid (4g) with a K(i) value of 0.072 microM, an EC(50) value of 4.12 microM, and a relative intrinsic activity of 75%. Of the aralkyl derivatives, 4q (5-(3-chlorobenzyl)pyrazole-3-carboxylic acid) was the most active with a relatively low intrinsic activity of 39%. Partial agonism of the pyrazole derivatives was confirmed by inhibition of G protein activation in response to nicotinic acid by these compounds. The pyrazoles both inhibited the maximum effect elicited by 100 microM nicotinic acid and concentration dependently shifted nicotinic acid concentration-response curves to the right, pointing to a competitive mechanism of action.