Identification of overlapping but differential binding sites for the high-affinity CXCR3 antagonists NBI-74330 and VUF11211.

CXC chemokine receptor CXCR3 and/or its main three ligands CXCL9, CXCL10, and CXCL11 are highly upregulated in a variety of diseases. As such, considerable efforts have been made to develop small-molecule receptor CXCR3 antagonists, yielding distinct chemical classes of antagonists blocking binding and/or function of CXCR3 chemokines. Although it is suggested that these compounds bind in an allosteric fashion, thus far no evidence has been provided regarding the molecular details of their interaction with CXCR3. Using site-directed mutagenesis complemented with in silico homology modeling, we report the binding modes of two high-affinity CXCR3 antagonists of distinct chemotypes: VUF11211 [(S)-5-chloro-6-(4-(1-(4-chlorobenzyl)piperidin-4-yl)-3-ethylpiperazin-1-yl)-N-ethylnicotinamide] (piperazinyl-piperidine) with a rigid elongated structure containing two basic groups and NBI-74330 [(R)-N-(1-(3-(4-ethoxyphenyl)-4-oxo-3,4-dihydropyrido[2,3-d]pyrimidin-2-yl)ethyl)-2-(4-fluoro-3-(trifluoromethyl)phenyl)-N-(pyridin-3-ylmethyl)acetamide] (8-azaquinazolinone) without any basic group. Here we show that NBI-74330 is anchored in the transmembrane minor pocket lined by helices 2 (W2.60, D2.63), 3 (F3.32), and 7 (S7.39, Y7.43), whereas VUF11211 extends from the minor pocket into the major pocket of the transmembrane domains, located between residues in helices 1 (Y1.39), 2 (W2.60), 3 (F3.32), 4 (D4.60), 6 (Y6.51), and 7 (S7.39, Y7.43). Mutation of these residues did not affect CXCL11 binding significantly, confirming the allosteric nature of the interaction of these small molecules with CXCR3. Moreover, the model derived from our in silico-guided studies fits well with the already published structure-activity relationship data on these ligands. Altogether, in this study, we show overlapping, yet different binding sites for two high-affinity CXCR3 antagonists, which offer new opportunities for the structure-based design of allosteric modulators for CXCR3.

From heptahelical bundle to hits from the Haystack: structure-based virtual screening for GPCR ligands.

This review will focus on the construction, refinement, and validation of G-protein-coupled receptor (GPCR) structural models for the purpose of structure-based virtual screening (SBVS) and ligand design. The review will present a comparative analysis of GPCR crystal structures and their implication on GPCR (homology) modeling. The challenges associated with steps along the modeling workflow will be discussed: the use of experimental anchors to steer the modeling procedure, amino acid sequence alignment and template selection, receptor structure refinement, loop modeling, ligand-binding mode prediction, and virtual screening for novel ligands. An overview of several successful structure-based ligand discovery and design studies shows that receptor models, despite structural inaccuracies, can be efficiently used to find novel ligands for GPCRs. Moreover, the recently solved GPCR crystal structures have further increased the opportunities in structure-based ligand discovery for this pharmaceutically important protein family.

Combining quantum mechanical ligand conformation analysis and protein modeling to elucidate GPCR-ligand binding modes.

SAR beyond protein-ligand interactions: By combining structure-affinity relationships, protein-ligand modeling studies, and quantum mechanical calculations, we show that ligand conformational energies and basicity play critical roles in ligand binding to the histamine H4 receptor, a GPCR that plays a key role in inflammation.

Chemical subtleties in small-molecule modulation of peptide receptor function: the case of CXCR3 biaryl-type ligands.

The G protein-coupled chemokine receptor CXCR3 plays a role in numerous inflammatory events. The endogenous ligands for the chemokine receptors are peptides, but in this study we disclose small-molecule ligands that are able to activate CXCR3. A class of biaryl-type compounds that is assembled by convenient synthetic routes is described as a new class of CXCR3 agonists. Intriguingly, structure-activity relationship and structure-function relationship studies reveal that subtle chemical modifications on the outer aryl ring (e.g., either the size or position of a halogen atom) result in a full spectrum of agonist efficacies on CXCR3. Quantum mechanics calculations and nuclear Overhauser effect spectroscopy NMR studies suggest that the biaryl dihedral angle and the electronic nature of ortho-substituents play an important role in determining agonist efficacies. Compounds 38 (VUF11222) and 39 (VUF11418) are the first reported nonpeptidomimetic agonists on CXCR3, rendering them highly useful chemical tools for detailed assessment of CXCR3 activation as well as for studying downstream CXCR3 signaling.

Regiochemistry of the condensation of 2-aroyl-cyclohexanones and 2-cyanoacetamide: 13C-labeling studies and semiempirical MO calculations.

Hydroxy-aryl-5,6,7,8-tetrahydroisoquinoline-4-carbonitriles represent interesting chemical scaffolds, but synthetic access to these compounds is limited. The reaction of 2-aroyl-cyclohexanones with 2-cyanoacetamide and base in ethanol has been reported to lead to the formation of the tetrahydroisoquinoline isomer. We show that depending on the electronic nature of the para-substituent on the aryl ring, formation of the regioisomeric tetrahydroquinoline isomer can significantly compete. The electron-donating or -withdrawing properties of the para-substituent of the aryl ring determines the ratio of product isomers. A series of 2-aroyl-cyclohexanones, with para-substituents ranging from electron-donating to electron-withdrawing, were reacted with [2-(13)C]-cyanoacetamide. The product ratio and absolute regiochemistry were directly determined by quantitative (13)C, HMBC, and NOESY NMR spectroscopy on the reaction mixtures. A clear relationship between the regioisomeric product ratio and the Hammett sigma values of the substituents is demonstrated. This is explained by the separate in situ yields, which reveal that the pathway leading to the tetrahydroquinoline regioisomer is significantly more sensitive toward the electronic nature of the para-substituent than the pathway leading to the tetrahydroisoquinoline. Semiempirical AM1 molecular orbital calculations on the starting electrophile 2-aroyl-cyclohexanone support a correlation between the energy of the LUMOs and the regioisomeric product ratio. Our results facilitate synthetic access to a range of these interesting synthetic intermediates.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design.

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the ß-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

Molecular determinants of ligand binding modes in the histamine H(4) receptor: linking ligand-based three-dimensional quantitative structure-activity relationship (3D-QSAR) models to in silico guided receptor mutagenesis studies.

The histamine H(4) receptor (H(4)R) is a G protein-coupled receptor (GPCR) that plays an important role in inflammation. Similar to the homologous histamine H(3) receptor (H(3)R), two acidic residues in the H(4)R binding pocket, D(3.32) and E(5.46), act as essential hydrogen bond acceptors of positively ionizable hydrogen bond donors in H(4)R ligands. Given the symmetric distribution of these complementary pharmacophore features in H(4)R and its ligands, different alternative ligand binding mode hypotheses have been proposed. The current study focuses on the elucidation of the molecular determinants of H(4)R-ligand binding modes by combining (3D) quantitative structure-activity relationship (QSAR), protein homology modeling, molecular dynamics simulations, and site-directed mutagenesis studies. We have designed and synthesized a series of clobenpropit (N-(4-chlorobenzyl)-S-[3-(4(5)-imidazolyl)propyl]isothiourea) derivatives to investigate H(4)R-ligand interactions and ligand binding orientations. Interestingly, our studies indicate that clobenpropit (2) itself can bind to H(4)R in two distinct binding modes, while the addition of a cyclohexyl group to the clobenpropit isothiourea moiety allows VUF5228 (5) to adopt only one specific binding mode in the H(4)R binding pocket. Our ligand-steered, experimentally supported protein modeling method gives new insights into ligand recognition by H(4)R and can be used as a general approach to elucidate the structure of protein-ligand complexes.

Identification of a novel allosteric binding site in the CXCR2 chemokine receptor.

We have shown previously that different chemical classes of small-molecule antagonists of the human chemokine CXCR2 receptor interact with distinct binding sites of the receptor. Although an intracellular binding site for diarylurea CXCR2 antagonists, such as N-(2-bromophenyl)-N'-(7-cyano-1H-benzotriazol-4-yl)urea (SB265610), and thiazolopyrimidine compounds was recently mapped by mutagenesis studies, we now report on an imidazolylpyrimidine antagonist binding pocket in the transmembrane domain of CXCR2. Using different CXCR2 orthologs, chimeric proteins, site-directed mutagenesis, and in silico modeling, we have elucidated the binding mode of this antagonist. Our in silico-guided mutagenesis studies indicate that the ligand binding cavity for imidazolylpyrimidine compounds in CXCR2 is located between transmembrane (TM) helices 3 (Phe130(3.36)), 5 (Ser217(5.44), Phe220(5.47)), and 6 (Asn268(6.52), Leu271(6.55)) and suggest that these antagonists enter CXCR2 via the TM5-TM6 interface. It is noteworthy that the same interface is postulated as the ligand entry channel in the opsin receptor and is occupied by lipid molecules in the recently solved crystal structure of the CXCR4 chemokine receptor, suggesting a general ligand entrance mechanism for nonpolar ligands to G protein-coupled receptors. The identification of a novel allosteric binding cavity in the TM domain of CXCR2, in addition to the previously identified intracellular binding site, shows the diversity in ligand recognition mechanisms by this receptor and offers new opportunities for the structure-based design of small allosteric modulators of CXCR2 in the future.

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs.

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ∼75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.

Application of a ligand-based theoretical approach to derive conversion paths and ligand conformations in CYP11B2-mediated aldosterone formation.

The biosynthesis of the mineralocorticoid hormone aldosterone involves a multistep hydroxylation of 11-deoxycorticosterone at the 11- and 18-positions, resulting in the formation of corticosterone and 18-hydroxycorticosterone, the final precursor of aldosterone. Two members of the cytochrome P450 11B family, CYP11B1 and CYP11B2, are known to catalyze these 11- and 18-hydroxylations, however, only CYP11B2 can oxidize 18-hydroxycorticosterone to aldosterone. It is unknown what sequence of hydroxylations leads to the formation of 18-hydroxycorticosterone. In this study we have investigated which of the possible conversion paths towards formation of 18-hydroxycorticosterone and aldosterone are most likely from the ligand perspective. Therefore, we combined quantum mechanical investigations on the steroid conformations of 11-deoxycorticosterone and its ensuing reaction intermediates with Fukui indices calculations to predict the reactivity of their carbon atoms for an attack by the iron-oxygen species. Both F(-) and F(0) were calculated to account for different mechanisms of substrate conversion. We show which particular initial conformations of 11-deoxycorticosterone and which conversion paths are likely to result in the successful synthesis of aldosterone, and thereby may be representative for the mechanism of aldosterone biosynthesis by CYP11B2. Moreover, we found that the most likely path for aldosterone synthesis coincides with the substrate conformation proposed in an earlier publication. To summarize, we show that on a theoretical and strictly ligand-directed basis only a limited number of reaction paths in the conversion of 11-deoxycorticosterone to aldosterone is possible. Despite its theoretical nature, this knowledge may help to understand the catalytic function of CYP11B1 and CYP11B2.

Synthesis, biological evaluation, and molecular modeling of 1-benzyl-1H-imidazoles as selective inhibitors of aldosterone synthase (CYP11B2).

Reducing aldosterone action is beneficial in various major diseases such as heart failure. Currently, this is achieved with mineralocorticoid receptor antagonists, however, aldosterone synthase (CYP11B2) inhibitors may offer a promising alternative. In this study, we used three-dimensional modeling of CYP11B2 to model the binding modes of the natural substrate 18-hydroxycorticosterone and the recently published CYP11B2 inhibitor R-fadrozole as a rational guide to design 44 structurally simple and achiral 1-benzyl-1H-imidazoles. Their syntheses, in vitro inhibitor potencies, and in silico docking are described. Some promising CYP11B2 inhibitors were identified, with our novel lead MOERAS115 (4-((5-phenyl-1H-imidazol-1-yl)methyl)benzonitrile) displaying an IC(50) for CYP11B2 of 1.7 nM, and a CYP11B2 (versus CYP11B1) selectivity of 16.5, comparable to R-fadrozole (IC(50) for CYP11B2 6.0 nM, selectivity 19.8). Molecular docking of the inhibitors in the models enabled us to generate posthoc hypotheses on their binding modes, providing a valuable basis for future studies and further design of CYP11B2 inhibitors.

Fadrozole reverses cardiac fibrosis in spontaneously hypertensive heart failure rats: discordant enantioselectivity versus reduction of plasma aldosterone.

Reversal of cardiac fibrosis is a major determinant of the salutary effects of mineralocorticoid receptor antagonists in heart failure. Recently, R-fadrozole was coined as an aldosterone biosynthesis inhibitor, offering an appealing alternative to mineralocorticoid receptor antagonists to block aldosterone action. The present study aimed to evaluate the effects of R- and S-fadrozole on plasma aldosterone and urinary aldosterone excretion rate and to compare their effectiveness vs. the mineralocorticoid receptor antagonist potassium canrenoate to reverse established cardiac fibrosis. Male lean spontaneously hypertensive heart failure (SHHF) rats (40 wk) were treated for 8 wk by sc infusions of low (0.24 mg/kg.d) or high (1.2 mg/kg.d) doses of R- or S-fadrozole or by potassium canrenoate via drinking water (7.5 mg/kg.d). At the high dose, plasma aldosterone levels were decreased similarly by R- and S-fadrozole, whereas urinary aldosterone excretion rate was reduced only by S-fadrozole. In contrast, whereas at the high dose, R-fadrozole effectively reversed preexistent left ventricular interstitial fibrosis by 50% (vs. 42% for canrenoate), S-fadrozole was devoid of an antifibrotic effect. The low doses of the fadrozole enantiomers did not change cardiac fibrosis or plasma aldosterone but similarly reduced urinary aldosterone excretion rate. In conclusion, R-fadrozole may possess considerable therapeutic merit because of its potent antifibrotic actions in the heart. However, the observed discordance between the aldosterone-lowering and antifibrotic effects of the fadrozole enantiomers raises some doubt about the mechanism by which R-fadrozole diminishes cardiac collagen and about the generality of the concept of lowering aldosterone levels to treat the diseased heart.

Construction of 3D models of the CYP11B family as a tool to predict ligand binding characteristics.

Aldosterone is synthesised by aldosterone synthase (CYP11B2). CYP11B2 has a highly homologous isoform, steroid 11beta-hydroxylase (CYP11B1), which is responsible for the biosynthesis of aldosterone precursors and glucocorticoids. To investigate aldosterone biosynthesis and facilitate the search for selective CYP11B2 inhibitors, we constructed three-dimensional models for CYP11B1 and CYP11B2 for both human and rat. The models were constructed based on the crystal structure of Pseudomonas Putida CYP101 and Oryctolagus Cuniculus CYP2C5. Small steric active site differences between the isoforms were found to be the most important determinants for the regioselective steroid synthesis. A possible explanation for these steric differences for the selective synthesis of aldosterone by CYP11B2 is presented. The activities of the known CYP11B inhibitors metyrapone, R-etomidate, R-fadrazole and S-fadrazole were determined using assays of V79MZ cells that express human CYP11B1 and CYP11B2, respectively. By investigating the inhibitors in the human CYP11B models using molecular docking and molecular dynamics simulations we were able to predict a similar trend in potency for the inhibitors as found in the in vitro assays. Importantly, based on the docking and dynamics simulations it is possible to understand the enantioselectivity of the human enzymes for the inhibitor fadrazole, the R-enantiomer being selective for CYP11B2 and the S-enantiomer being selective for CYP11B1.

Switching from S- to R-selectivity in the Candida antarctica lipase B-catalyzed ring-opening of omega-methylated lactones: tuning polymerizations by ring size.

Novozym 435-catalyzed ring-opening of a range of omega-methylated lactones demonstrates fascinating differences in rate of reaction and enantioselectivity. A switch from S- to R-selectivity was observed upon going from small (ring sizes or=8). This was attributed to the transition from a cisoid to a transoid conformational preference of the ester bond on going from small to large lactones. The S-selectivity of the ring-opening of the small, cisoid lactones was low to moderate, while the R-selectivity of the ring-opening of the large transoid lactones was surprisingly high. The S-selectivity of the ring-opening of the small, cisoid lactones combined with the established R-selectivity of the transesterification of (aliphatic) secondary alcohols prevented polymerization from taking place. Ring-opening of the large, transoid lactones was R-selective with high enantioselectivity. As a result, these lactones could be polymerized, without exception, by straightforward kinetic resolution polymerization, yielding the enantiopure R-polyester with excellent enantiomeric excess (>99%).