Multi-objective evolutionary design of adenosine receptor ligands.

A novel multiobjective evolutionary algorithm (MOEA) for de novo design was developed and applied to the discovery of new adenosine receptor antagonists. This method consists of several iterative cycles of structure generation, evaluation, and selection. We applied an evolutionary algorithm (the so-called Molecule Commander) to generate candidate A1 adenosine receptor antagonists, which were evaluated against multiple criteria and objectives consisting of high (predicted) affinity and selectivity for the receptor, together with good ADMET properties. A pharmacophore model for the human A1 adenosine receptor (hA1AR) was created to serve as an objective function for evolution. In addition, three support vector machine models based on molecular fingerprints were developed for the other adenosine receptor subtypes (hA2A, hA2B, and hA3) and applied as negative objective functions, to aim for selectivity. Structures with a higher evolutionary fitness with respect to ADMET and pharmacophore matching scores were selected as input for the next generation and thus developed toward overall fitter ("better") compounds. We finally obtained a collection of 3946 unique compounds from which we derived chemical scaffolds. As a proof-of-principle, six of these templates were selected for actual synthesis and subsequently tested for activity toward all adenosine receptors subtypes. Interestingly, scaffolds 2 and 3 displayed low micromolar affinity for many of the adenosine receptor subtypes. To further investigate our evolutionary design method, we performed systematic modifications on scaffold 3. These modifications were guided by the substitution patterns as observed in the set of generated compounds that contained scaffold 3. We found that an increased affinity with appreciable selectivity for hA1AR over the other adenosine receptor subtypes was achieved through substitution of the scaffold; compound 3a had a Ki value of 280 nM with approximately 10-fold selectivity with respect to hA2AR, while 3g had a 1.6 µM affinity for hA1AR with negligible affinity for the hA2A, hA2B, and hA3 receptor subtypes.

A new generation of adenosine receptor antagonists: from di- to trisubstituted aminopyrimidines.

New adenosine receptor ligands were designed as hybrid structures between previously synthesized substituted dicyanopyridines and aminopyrimidines, yielding two series of cyano-substituted diphenylaminopyrimidines. We were interested in assessing the effect of this substitution pattern on both affinity and intrinsic activity, as the dicyanopyridines comprised both agonists and inverse agonists, whereas the original aminopyrimidines were exclusively inverse agonists. It was found that the new compounds were generally selective for adenosine A(1) receptors, although affinity for the adenosine A(2A) receptor was also noticed for some of the compounds. In a cAMP second messenger assay the compounds behaved as inverse agonists rather than agonists. Among the more A(1) receptor-selective compounds were 5 (LUF6048), 27 (LUF6040) and 53 (LUF6056) with K(i) values of 8.1, 1.2 and 5.7nM, respectively.

Designing active template molecules by combining computational de novo design and human chemist's expertise.

We used a new software tool for de novo design, the "Molecule Evoluator", to generate a number of small molecules. Explicit constraints were a relatively low molecular weight and otherwise limited functionality, for example, low numbers of hydrogen bond donors and acceptors, one or two aromatic rings, and a small number of rotatable bonds. In this way, we obtained a collection of scaffold- or templatelike molecules rather than fully "decorated" ones. We asked medicinal chemists to evaluate the suggested molecules for ease of synthesis and overall appeal, allowing them to make structural changes to the molecules for these reasons. On the basis of their recommendations, we synthesized eight molecules with an unprecedented (not patented) yet simple structure, which were subsequently tested in a screen of 83 drug targets, mostly G protein-coupled receptors. Four compounds showed affinity for biogenic amine targets (receptor, ion channel, and transport protein), reflecting the training of the medicinal chemists involved. Apparently the generation of leadlike solutions helped the medicinal chemists to select good starting points for future lead optimization, away from existing compound libraries.

2,6,8-trisubstituted 1-deazapurines as adenosine receptor antagonists.

In this study we developed a refined pharmacophore model for antagonists of the human adenosine A1 receptor, based on features of known pyrimidine and purine derivatives. The adoption of these updated criteria assisted us in synthesizing a series of 1-deazapurines with consistently high affinity as inverse agonists for the adenosine A1 receptor. These 1-deazapurines (otherwise known as 3H-imidazo[4,5-b]pyridines) were substituted at their 2- and 6-positions, yielding a series with five of the derivatives displaying Ki values in the subnanomolar range. The most potent of these, compound 10 (LUF 5978), displayed an affinity of 0.55 nM at the human adenosine A1 receptor with >300-fold and 45-fold selectivity toward A2A and A3 receptors, respectively. Compound 14 (LUF 5981, Ki = 0.90 nM) appeared to have the best overall selectivity with respect to adenosine A2A (>200-fold) and A3 (700-fold) receptors.

2,6-disubstituted and 2,6,8-trisubstituted purines as adenosine receptor antagonists.

Purines have long been exploited as adenosine receptor antagonists. The substitution pattern about the purine ring has been well investigated, and certain criteria have become almost a prerequisite for good affinity at the adenosine A(1) receptor. The adaptation of the pharmacophore and the initial series of pyrimidines developed in an earlier publication resulted in a series of purines with an entirely new substitution pattern. One compound in particular, 8-cyclopentyl-2,6-diphenylpurine (31, LUF 5962) has been shown to be very promising with an affinity of 0.29 nM at the human adenosine A(1) receptor.

Structure-activity relationships of new 1H-imidazo[4,5-c]quinolin-4-amine derivatives as allosteric enhancers of the A3 adenosine receptor.

1H-Imidazo[4,5-c]quinolin-4-amine derivatives have been synthesized as allosteric modulators of the human A3 adenosine receptor (AR). Structural modifications were made at the 4-amino and 2 positions. The compounds were tested in both binding and functional assays, and many were found to be allosteric enhancers of the action of A3AR agonists by several different criteria. First, a potentiation of the maximum efficacy of the agonist Cl-IB-MECA was observed for numerous derivatives. Also, a number of these compounds decreased the rate of dissociation of the agonist [125I]I-AB-MECA from the A3AR. Most prominently, compound 43 (LUF6000) was found to enhance agonist efficacy in a functional assay by 45% and decrease dissociation rate similarly without influencing agonist potency. The structural requirements for allosteric enhancement at the A3AR were distinct from the requirements to inhibit equilibrium binding. Thus, we have prepared allosteric enhancers of the human A3AR that have an improved allosteric effect in comparison to the inhibition of equilibrium binding at the orthosteric site.

Synthesis and biological evaluation of a new series of 2,3,5-substituted [1,2,4]-thiadiazoles as modulators of adenosine A1 receptors and their molecular mechanism of action.

We synthesized two series (7a-i and 8a-i) of 2,3,5-substituted [1,2,4]-thiadiazole analogues of SCH-202676 (7a, 2,3-diphenyl-5-N-methylimino-2H-[1,2,4]-thiadiazole) with emphasis on the N-imino substituent. Compounds 7a-g,i and 8a-g at a final concentration of 1 microM significantly inhibited [(3)H]CCPA (2-chloro-N(6)- cyclopentyladenosine) agonist binding to human A(1) adenosine receptors. At the same concentration, all compounds appeared to increase [(3)H]DPCPX (1,3-dipropyl-8-cyclopentylxanthine) antagonist binding. Compound 7a and LUF5855 (7g) were selected for further characterization and studied in both equilibrium and kinetic radioligand binding experiments. The results suggest a nonstoichiometric interaction with the receptor. Further bioanalytical procedures (HPLC and MS) provided proof for an unusual receptor interaction in which 7a and 7g upon incubation were transformed into their corresponding thioureas 5a and 5g. We suggest that the thiadiazoles are sulfhydryl modifying agents rather than allosteric modulators, as they appear to reversibly modify the sulfhydryl groups of cysteine residues in cell membrane preparations.

A series of ligands displaying a remarkable agonistic-antagonistic profile at the adenosine A1 receptor.

Adenosine receptor agonists are usually variations of the natural ligand, adenosine. The ribose moiety in the ligand has previously been shown to be of great importance for the agonistic effects of the compound. In this paper, we present a series of nonadenosine ligands selective for the adenosine A(1) receptor with an extraordinary pharmacological profile. 2-Amino-4-benzo[1,3]dioxol-5-yl-6-(2-hydroxyethylsulfanyl)pyridine-3,5-dicarbonitrile (70, LUF 5853) shows full agonistic behavior comparable with the reference compound CPA, while also displaying comparable receptor binding affinity (K(i) = 11 nM). In contrast, compound 58 (2-amino-4-(3-trifluoromethylphenyl)-6-(2-hydroxyethylsulfanyl)pyridine-3,5-dicarbonitrile, LUF 5948) has a binding affinity of 14 nM and acts as an inverse agonist. Also present within this same series are compounds that show neutral antagonism of the adenosine A(1) receptor, for example compound 65 (2-amino-4-(4-difluoromethoxyphenyl)-6-(2-hydroxyethylsulfanyl)pyridine-3,5-dicarbonitrile, LUF 5826).

Evaluation of (4-Arylpiperidin-1-yl)cyclopentanecarboxamides As High-Affinity and Long-Residence-Time Antagonists for the CCR2 Receptor.

Animal models suggest that the chemokine ligand 2/CC-chemokine receptor 2 (CCL2/CCR2) axis plays an important role in the development of inflammatory diseases. However, CCR2 antagonists have failed in clinical trials because of a lack of efficacy. We previously described a new approach for the design of CCR2 antagonists by the use of structure-kinetics relationships (SKRs). Herein we report new findings on the structure-affinity relationships (SARs) and SKRs of the reference compound MK-0483, its diastereomers, and its structural analogues as CCR2 antagonists. The SARs of the 4-arylpiperidine group suggest that lipophilic hydrogen-bond-accepting substituents at the 3-position are favorable. However, the SKRs suggest that a lipophilic group with a certain size is desired [e.g., 3-Br: Ki =2.8 nM, residence time (t(res))=243 min; 3-iPr: Ki =3.6 nM, t(res) =266 min]. Alternatively, additional substituents and further optimization of the molecule, while keeping a carboxylic acid at the 3-position, can also prolong t(res); this was most prominently observed in MK-0483 (Ki =1.2 nM, t(res) =724 min) and a close analogue (Ki =7.8 nM) with a short residence time.

When structure-affinity relationships meet structure-kinetics relationships: 3-((Inden-1-yl)amino)-1-isopropyl-cyclopentane-1-carboxamides as CCR2 antagonists.

Chemokine ligand 2 (CCL2) mediates chemotaxis of monocytes to inflammatory sites via interaction with its G protein-coupled receptor CCR2. Preclinical animal models suggest that the CCL2-CCR2 axis has a critical role in the development and maintenance of inflammatory disease states (e.g., multiple sclerosis, atherosclerosis, insulin resistance, restenosis, and neuropathic pain), which can be treated through inhibition of the CCR2 receptor. However, in clinical trials high-affinity inhibitors of CCR2 have often demonstrated a lack of efficacy. We have previously described a new approach for the design of high-affinity CCR2 antagonists, by taking their residence time (RT) on the receptor into account. Here, we report our findings on both structure-affinity relationship (SAR) and structure-kinetic relationship (SKR) studies for a series of 3-((inden-1-yl)amino)-1-isopropyl-cyclopentane-1-carboxamides as CCR2 antagonists. SAR studies showed that this class of compounds tolerates a vast diversity of substituents on the indenyl ring with only small changes in affinity. However, the SKR is affected greatly by minor modifications of the structure. The combination of SAR and SKR in the hit-to-lead process resulted in the discovery of a new high-affinity and long-residence-time CCR2 antagonist (compound 15a, Ki = 2.4 nM; RT = 714 min).

New, non-adenosine, high-potency agonists for the human adenosine A2B receptor with an improved selectivity profile compared to the reference agonist N-ethylcarboxamidoadenosine.

The adenosine A(2B) receptor is the least well characterized of the four known adenosine receptor subtypes because of the absence of potent, selective agonists. Here, we present five non-adenosine agonists. Among them, 2-amino-4-(4-hydroxyphenyl)-6-(1H-imidazol-2-ylmethylsulfanyl)pyridine-3,5-dicarbonitrile, 17, LUF5834, is a high-efficacy partial agonist with EC(50) = 12 nM and 45-fold selectivity over the adenosine A(3) receptor but lacking selectivity versus the A(1) and A(2A) subtypes. Compound 18, LUF5835, the 3-hydroxyphenyl analogue, is a full agonist with EC(50) = 10 nM.

2,4,6-trisubstituted pyrimidines as a new class of selective adenosine A1 receptor antagonists.

Adenosine receptor antagonists usually possess a bi- or tricyclic heteroaromatic structure at their core with varying substitution patterns to achieve selectivity and/or greater affinity. Taking into account molecular modeling results from a series of potent adenosine A1 receptor antagonists, a pharmacophore was derived from which we show that a monocyclic core can be equally effective. To achieve a compound that may act at the CNS we propose imposing a restriction related to its polar surface area (PSA). In consequence, we have synthesized two novel series of pyrimidines, possessing good potency at the adenosine A1 receptor and desirable PSA values. In particular, compound 30 (LUF 5735) displays excellent A1 affinity (Ki = 4 nM) and selectivity (< or =50% displacement of 1 muM concentrations of the radioligand at the other three adenosine receptors) and has a PSA value of 53 A2.

Binding kinetics of ZM241385 derivatives at the human adenosine A2A receptor.

Classical drug design and development rely mostly on affinity- or potency-driven structure-activity relationships (SAR). Thus far, a given compound's binding kinetics have been largely ignored, the importance of which is now being increasingly recognized. In the present study, we performed an extensive structure-kinetics relationship (SKR) study in addition to a traditional SAR analysis at the adenosine A2A receptor (A2A R). The ensemble of 24 A2A R compounds, all triazolotriazine derivatives resembling the prototypic antagonist ZM241385 (4-(2-((7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-yl)amino)ethyl)phenol), displayed only minor differences in affinity, although they varied substantially in their dissociation rates from the receptor. We believe that such a combination of SKR and SAR analyses, as we have done with the A2A R, will have general importance for the superfamily of G protein-coupled receptors, as it can serve as a new strategy to tailor the interaction between ligand and receptor.

Agonists for the adenosine A1 receptor with tunable residence time. A Case for nonribose 4-amino-6-aryl-5-cyano-2-thiopyrimidines.

We report the synthesis and evaluation of previously unreported 4-amino-6-aryl-5-cyano-2-thiopyrimidines as selective human adenosine A1 receptor (hA1AR) agonists with tunable binding kinetics, this without affecting their nanomolar affinity for the target receptor. They show a very diverse range of kinetic profiles (from 1 min (compound 52) to 1 h (compound 43)), and their structure-affinity relationships (SAR) and structure-kinetics relationships (SKR) were established. When put in perspective with the increasing importance of binding kinetics in drug discovery, these results bring new evidence of the consequences of affinity-only driven selection of drug candidates, that is, the potential elimination of slightly less active compounds that may display preferable binding kinetics.

Strategies to reduce HERG K+ channel blockade. Exploring heteroaromaticity and rigidity in novel pyridine analogues of dofetilide.

Drug-induced blockade of the human ether-a-go-go-related gene K(+) channel (hERG) represents one of the major antitarget concerns in pharmaceutical industry. SAR studies of this ion channel have shed light on the structural requirements for hERG interaction but most importantly may reveal drug design principles to reduce hERG affinity. In the present study, a novel library of neutral and positively charged heteroaromatic derivatives of the class III antiarrhythmic agent dofetilide was synthesized and assessed for hERG affinity in radioligand binding and manual patch clamp assays. Structural modifications of the pyridine moiety, side chain, and peripheral aromatic moieties were evaluated, thereby revealing approaches for reducing hERG binding affinity. In particular, we found that the extra rigidity imposed close to the positively charged pyridine moiety can be very efficient in decreasing hERG affinity.

Dual-point competition association assay: a fast and high-throughput kinetic screening method for assessing ligand-receptor binding kinetics.

The concept of ligand-receptor binding kinetics is emerging as an important parameter in the early phase of drug discovery. Since the currently used kinetic assays are laborious and low throughput, we developed a method that enables fast and large format screening. It is a so-called dual-point competition association assay, which measures radioligand binding at two different time points in the absence or presence of unlabeled competitors. Specifically, this assay yields the kinetic rate index (KRI), which is a measure for the binding kinetics of the unlabeled ligands screened. As a prototypical drug target, the adenosine A(1) receptor (A(1)R) was chosen for assay validation and optimization. A screen with 35 high-affinity A(1)R antagonists yielded seven compounds with a KRI value above 1.0, which indicated a relatively slow dissociation from the target. All other compounds had a KRI value below or equal to 1.0, predicting a relatively fast dissociation rate. Several compounds were selected for follow-up kinetic quantifications in classical kinetic assays and were shown to have kinetic rates that corresponded to their KRI values. The dual-point assay and KRI value may have general applicability at other G-protein-coupled receptors, as well as at drug targets from other protein families.

Removal of human ether-à-go-go related gene (hERG) K+ channel affinity through rigidity: a case of clofilium analogues.

Cardiotoxicity is a side effect that plagues modern drug design and is very often due to the off-target blockade of the human ether-à-go-go related gene (hERG) potassium channel. To better understand the structural determinants of this blockade, we designed and synthesized a series of 40 derivatives of clofilium, a class III antiarrhythmic agent. These were evaluated in radioligand binding and patch-clamp assays to establish structure-affinity relationships (SAR) for this potassium channel. Efforts were especially focused on studying the influence of the structural rigidity and the nature of the linkers composing the clofilium scaffold. It was shown that introducing triple bonds and oxygen atoms in the n-butyl linker of the molecule greatly reduced affinity without significantly modifying the pKa of the essential basic nitrogen. These findings could prove useful in the first stages of drug discovery as a systematic way of reducing the risk of hERG K(+) channel blockade-induced cardiotoxicity.

Structure-kinetic relationships--an overlooked parameter in hit-to-lead optimization: a case of cyclopentylamines as chemokine receptor 2 antagonists.

Preclinical models of inflammatory diseases (e.g., neuropathic pain, rheumatoid arthritis, and multiple sclerosis) have pointed to a critical role of the chemokine receptor 2 (CCR2) and chemokine ligand 2 (CCL2). However, one of the biggest problems of high-affinity inhibitors of CCR2 is their lack of efficacy in clinical trials. We report a new approach for the design of high-affinity and long-residence-time CCR2 antagonists. We developed a new competition association assay for CCR2, which allows us to investigate the relation of the structure of the ligand and its receptor residence time [i.e., structure-kinetic relationship (SKR)] next to a traditional structure-affinity relationship (SAR). By applying combined knowledge of SAR and SKR, we were able to re-evaluate the hit-to-lead process of cyclopentylamines as CCR2 antagonists. Affinity-based optimization yielded compound 1 with good binding (Ki = 6.8 nM) but very short residence time (2.4 min). However, when the optimization was also based on residence time, the hit-to-lead process yielded compound 22a, a new high-affinity CCR2 antagonist (3.6 nM), with a residence time of 135 min.

Understanding of molecular substructures that contribute to hERG K+ channel blockade: synthesis and biological evaluation of E-4031 analogues.

Cardiotoxicity is a common side effect of a large variety of drugs that is often caused by off-target human ether-à-go-go-related gene (hERG) potassium channel blockade. In this study, we designed and synthesized a series of derivatives of the class III antiarrhythmic agent E-4031. These compounds where evaluated in a radioligand binding assay and automated patch clamp assay to establish structure-activity relationships (SAR) for their inhibition of the hERG K(+) channel. Structural modifications of E-4031 were made by altering the peripheral aromatic moieties with a series of distinct substituents. Additionally, we synthesized several derivatives with a quaternary nitrogen and modified the center of the molecule by introduction of an additional nitrogen and deletion of the carbonyl oxygen. Some modifications caused a great increase in affinity for the hERG K(+) channel, while other seemingly minor changes led to a strongly diminished affinity. Structures with quaternary amines carrying an additional aromatic moiety were found to be highly active in radioligand binding assay. A decrease in affinity was achieved by introducing an amide functionality in the central scaffold without directly interfering with the pK(a) of the essential basic amine. The knowledge gained from this study could be used in early stages of drug discovery and drug development to avoid or circumvent hERG K(+) channel blockade, thereby reducing the risk of cardiotoxicity, related arrhythmias and sudden death.

Determination of different putative allosteric binding pockets at the lutropin receptor by using diverse drug-like low molecular weight ligands.

The lutropin/choriogonadotrophin receptor (LHCGR) is a family A G protein-coupled receptor (GPCR) which binds the endogenous hormone-ligands at the large extracellular domain. In contrast, several drug-like low-molecular-weight ligands (LMWs) have been reported to interact allosterically within the seven transmembrane domain (7TMD) of the LHCGR. Here, we were interested to study the putative allosteric LHCGR binding region with focus on the determination of two pockets for LMW ligands. A library of compounds was screened for their ability to modify the binding of an allosteric radiolabeled LMW agonist [³H]Org 43553. Further experimental and computational studies revealed that the putative binding pocket for a newly identified allosteric enhancer (LUF5419) and a previously described allosteric inhibitor (LUF5771) are overlapping and that this site is different from the Org 43553 binding site. The present study showed that these compounds are useful tools to reveal details on different allosteric binding sites located within the 7TMD of the LHCGR.