Comparative pharmacophore modeling of human adenosine receptor A1 and A3 antagonists

Adenosine receptors are promising therapeutic targets in drug discovery. In this study, three-dimensional pharmacophore models of human adenosine receptor A1 and A3 antagonists were developed based on 26 and 23 diverse compounds, respectively. The best A1 pharmacophore model (A 1 _Hopy1) consists of four features: one hydrogen bond donor, one hydrophobic point and two ring aromatics, while the best A 3 pharmacophore model (A3 _Hopy1) also has four features: one hydrogen bond acceptor, one hydrophobic point and two ring aromatics. The correlation coefficients were 0.840 for A 1 test set with 146 diverse compounds and 0.827 for A3 test set with 238 diverse compounds. In the simulated virtual screening experiments, high enrichment factors of 6.51 and 6.90 were obtained for A 1 _Hopy1 and A3 _Hopy1 models, respectively. Moreover, two models also showed high subtype-selectivity in the simulated virtual screening experiments. These results could be helpful for the discovery of novel potent and selective A 1 and A3 antagonists.     

A combined structure-based pharmacophore modeling and 3D-QSAR study on a series of N-heterocyclic scaffolds to screen novel antagonists as human DHFR inhibitors

Structural Chemistry - Dihydrofolate reductase (DHFR) is an essential enzyme that participates in folate metabolism and purine and thymidylate synthesis in cell proliferation. It converts...     

3D-pharmacophores of flavonoid binding at the benzodiazepine GABA(A) receptor site using 4D-QSAR analysis

4D-QSAR analysis was applied to a training set of 38 flavonoids where affinity constants, Ki, to the GABA(A) benzodiazepine receptor site, BzR, were determined. It was found that the -logKi values of the compounds are highly dependent on the size and electrostatics character of the substituents at the R(3') and R(6) positions of the flavonoid scaffold. Polar negative groups correctly embedded in the R(3') and/or R(6) substituents are predicted to increase -logKi values. A planar conformation of the flavonoid scaffold was found not to be a requirement for the flavonoids to be active. A test set of four compounds was used to evaluate the predictivity of the 4D-QSAR models.     

A cation-pi binding interaction with a tyrosine in the binding site of the GABAC receptor

GABA(C) (rho) receptors are members of the Cys-loop superfamily of neurotransmitter receptors, which includes nicotinic acetylcholine (nACh), 5-HT(3), and glycine receptors. As in other members of this family, the agonist binding site of GABA(C) receptors is rich in aromatic amino acids, but while other receptors bind agonist through a cation-pi interaction to a tryptophan, the GABA(C) binding site has tyrosine at the aligning positions. Incorporating a series of tyrosine derivatives at position 198 using unnatural amino acid mutagenesis reveals a clear correlation between the cation-pi binding ability of the side chain and EC(50) for receptor activation, thus demonstrating a cation-pi interaction between a tyrosine side chain and a neurotransmitter. Comparisons among four homologous receptors show variations in cation-pi binding energies that reflect the nature of the cationic center of the agonist.     

Substituted pyrazolo[3,4-b]pyridines as human A1 adenosine antagonists: developments in understanding the receptor stereoselectivity

A(1) adenosine receptor antagonists have been proposed to possess an interesting range of potential therapeutic applications. We have already reported the synthesis and the biological characterization of a family of pyrazolo[3,4-b]pyridine derivatives as A(1) adenosine ligands endowed with an antagonistic profile. In the present work, we report the LC separation of enantiomers of our most active A(1) antagonists together with the determination of their absolute configuration by means of X-ray crystal structure analysis. Biological assays confirmed a different activity for the two enantiomers, with the R one showing the higher human A(1)AR affinity. We also developed a homology model of this receptor subtype in order to suggest a binding disposition of the ligands into the hA(1)AR. All of the obtained data suggest that the compound's chirality plays a key role in A(1) affinity.     

Delineation of agonist binding to the human histamine H4 receptor using mutational analysis, homology modeling, and ab initio calculations

A three-dimensional homology model of the human histamine H 4 receptor was developed to investigate the binding mode of a series of structurally diverse H 4-agonists, i.e. histamine, clozapine, and the recently described selective, nonimidazole agonist VUF 8430. Mutagenesis studies and docking of these ligands in a rhodopsin-based homology model revealed two essential points of interactions in the binding pocket, i.e. Asp3.32 and Glu5.46 (Ballesteros-Weinstein numbering system). It is postulated that Asp3.32 interacts in its anionic state, whereas Glu5.46 interacts in its neutral form. The hypothesis was tested with the point mutations D3.32N and E5.46Q. For the D3.32N no binding affinity toward any of the ligands could be detected. This is in sharp contrast to the E5.46Q mutant, which discriminates between various ligands. The affinity of histamine-like ligands was decreased approximately a 1000-fold, whereas the affinity of all other ligands remained virtually unchanged. The proposed model for agonist binding as well as ab initio calculations for histamine and VUF 8430 explain the observed differences in binding to the H 4R mutants. These studies provide a molecular understanding for the action of a variety of H 4 receptor-ligands. The resulting H 4 receptor model will be the basis for the development of new H 4 receptor-ligands.     

Ligand-, structure- and pharmacophore-based molecular fingerprints: a case study on adenosine A1, A2A, A2B, and A3 receptor antagonists

FLAP fingerprints are applied in the ligand-, structure- and pharmacophore-based mode in a case study on antagonists of all four adenosine receptor (AR) subtypes. Structurally diverse antagonist collections with respect to the different ARs were constructed by including binding data to human species only. FLAP models well discriminate ??active?? (=highly potent) from ??inactive?? (=weakly potent) AR antagonists, as indicated by enrichment curves, numbers of false positives, and AUC values. For all FLAP modes, model predictivity slightly decreases as follows: A_2BR?>?A_2AR?>?A₃R?>?A₁R antagonists. General performance of FLAP modes in this study is: ligand-?>?structure-?>?pharmacophore- based mode. We also compared the FLAP performance with other common ligand- and structure-based fingerprints. Concerning the ligand-based mode, FLAP model performance is superior to ECFP4 and ROCS for all AR subtypes. Although focusing on the early first part of the A_2A, A_2B and A₃ enrichment curves, ECFP4 and ROCS still retain a satisfactory retrieval of actives. FLAP is also superior when comparing the structure-based mode with PLANTS and GOLD. In this study we applied for the first time the novel FLAPPharm tool for pharmacophore generation. Pharmacophore hypotheses, generated with this tool, convincingly match with formerly published data. Finally, we could demonstrate the capability of FLAP models to uncover selectivity aspects although single AR subtype models were not trained for this purpose.     

A neoceptor approach to unraveling microscopic interactions between the human A2A adenosine receptor and its agonists

Strategically mutated neoceptors, e.g., with anionic residues in TMs 3 and 7 intended for pairing with positively charged amine-modified nucleosides, were derived from the antiinflammatory A(2A) adenosine receptor (AR). Adenosine derivatives functionalized at the 5', 2, and N(6) positions were synthesized. The T88D mutation selectively enhanced the binding of the chain-length-optimized 5'-(2-aminoethyl)uronamide but not 5'-(2-hydroxyethyl)uronamide, suggesting a critical role of the positively charged amine. Combination of this modification with the N(6)-(2-methylbenzyl) group enhanced affinity at the Q89D- and N181D- but not the T88D-A(2A)AR. Amino groups placed near the 2- or N(6)-position only slightly affected the binding to mutant receptors. The 5'-hydrazide MRS3412 was 670- and 161-fold enhanced, in binding and functionally, respectively, at the Q89D-A(2A)AR compared to the wild-type. Thus, we identified and modeled pairs of A(2A)AR-derived neoceptor-neoligand, which are pharmacologically orthogonal with respect to the native species.     

On the stability of 2-aminoselenophene-3-carboxylates: potential dual-acting selenium-containing allosteric enhancers of A1 adenosine receptor binding

Ethyl-2-amino-4,5,6,7-tetrahydro-1-benzoselenophene-3-carboxylate (4), has been prepared as a potential dual-acting selenium-containing allosteric enhancer of adenosine A(1)A receptor binding utilising a modified Gewald reaction. While preliminary testing indicated that 4 is a superior enhancer of A(1)AR binding than its thiophene counterpart, its instability under mildly acidic conditions is cause for concern. X-Ray crystallography, together with DFT calculations, provide evidence that the decomposition of 4 involves the ring-opening of selenophenium ion (12b) followed by the loss of elemental selenium through a radical chain process.     

Direct NMR detection of the binding of functional ligands to the human sweet receptor, a heterodimeric family 3 GPCR

We present a robust method for monitoring the binding of ligands to the heterodimeric (T1R2+T1R3) human sweet receptor (a family 3 GPCR receptor). The approach utilizes saturation transfer difference (STD) NMR spectroscopy with receptor proteins expressed on the surface of human epithelial kidney cells. The preparation investigated by NMR can contain either live cells or membranes isolated from these cells containing the receptor. We have used this approach to confirm the noncompetitive binding of alitame and cyclamate to the receptor and to determine that greatly reduced receptor binding affinity compared to wild-type brazzein explains the lack of sweetness of brazzein mutant A16C17. This approach opens new avenues for research on the mechanism of action of the sweet receptor and for the design of new noncalorigenic sweeteners.     

Revealing the ionization ability of binding site I of human serum albumin using 2-(2'-hydroxyphenyl)benzoxazole as a pH sensitive probe

The ability of site I of human serum albumin (HSA) to bind medium sized molecules is important for the distribution, metabolism, and efficacy of many drugs. Herein, we show that this binding site has the ionization ability that may alter the drug structure during the process of its delivery. We reveal this ability by employing 2-(2'-hydroxyphenyl)benzoxazole (HBO) as a pH sensitive probe. Binding of HBO in site I is studied here at physiological pH 7.2 using steady-state and lifetime spectroscopic measurements, molecular docking and molecular dynamics (MD) simulation methods. The complex photophysics of HBO and the unique fluorescence signature of its anionic form indicate that, upon binding with HSA, the molecule exists in equilibrium between the anionic and the syn-keto forms. The position of HBO inside the binding site was determined experimentally by measuring the fluorescence quenching of W214, the sole tryptophan residue in HSA. The ionization degree of HBO inside the binding site was estimated to be close to the ionization degree of HBO in an aqueous solution of pH 10. This was concluded by comparing the fluorescence behavior of bound HBO to that of HBO in different solvents and in aqueous solutions of different pH values. Molecular docking and MD simulations show that HBO binds in site I close to W214, confirming the experimental results, and pinpoint the dominant role of hydrophobic interactions in the binding site. The formation of the anionic form is proposed to be due to through-space interaction between the OH group of HBO and both R222 and I290 with a binding mode similar to that of warfarin in site I. Comparison of the results with those of HBO mixed with key amino acids in solution indicates the importance of through-space interaction in the formation of the anion, similar to enzymatic reactions.     

A 3D QSAR approach to the search for geometrical similarity in a series of nonpeptide angiotensin II receptor antagonists

Summary A 3D QSAR methodology based on the combined use of conformational analysis and chemometrics was applied to perform a comparative analysis of the 3D conformational features of 13 nonpeptide angiotensin II receptor antagonists showing different levels of binding affinity. Conformational analysis by using a molecular mechanics MM2 method was carried out for each of these structures to obtain conformational minima. These minima were described by ten interatomic distances which define the relative spatial disposition of five significant atoms belonging to relevant functional groups present in all the 13 molecules. The structure-activity relationship between the interatomic distances and the biological activity was then assessed by using chemometric methods (cluster analysis, principal component analysis, classification methods). With our indirect approach based on the search for geometrical similarity it was possible, even though structural information on the receptor active site was lacking, to identify the 3D geometrical requirements for the binding affinity of nonpeptide angiotensin II receptor inhibitors.     

Molecular modeling and mutagenesis reveals a tetradentate binding site for Zn2+ in GABA(A) alphabeta receptors and provides a structural basis for the modulating effect of the gamma subunit

Gamma-aminobutyric acid type A receptors (GABA(A)-R) containing alpha1beta2gamma2 subunits are weakly inhibited by Zn2+, whereas receptors containing only the alpha1beta2 subunits are strongly inhibited. We built homology models of the ion pores of alpha1beta2 and alpha1beta2gamma2 GABA(A)-R using coordinates of the nicotinic acetylcholine receptor as a template. Threading the GABA(A)-R beta2 sequence onto this template placed the 17' histidine and the 20' glutamate residues at adjacent locations in the mouth of the pore, such that a nearly ideal tetradentate site for Zn2+ was formed from two histidine and two glutamate residues between adjacent beta subunits in the alpha1beta2 GABA(A)-R. Following optimization with CHARMM, the distance between the alpha-carbons of the adjacent histidine residues was approximately 9.2 A, close to the ideal distance for a Zn2+ binding site. Loss of inhibition by Zn2+ in alpha1beta2gamma2 GABA(A)-R can be explained by the geometry of these residues in the arrangement alpha1beta2gamma2alpha1beta2, in which the nearest C-alpha-C-alpha distance between the histidine residues is 15.5 A, too far apart for an energetically optimal Zn2+ binding site. We then mutated the gamma subunit at the 17' and/or 20' positions. Zn2+ inhibition was not restored in alpha1beta2gamma2 (I282H) receptors. A novel finding is that the modeling shows the native 20' lysine in gamma2 can compete with Zn2+ for binding to the inserted 17' histidine. Sensitivity to Zn2+ was restored in the double mutant receptor, alpha1beta2gamma2 (I282H; K285E), in which the competition with lysine was removed and a more favorable Zn2+ binding site was formed.     

Polyaniline immobilized on polycaprolactam nanofibers as a sorbent in electrochemically controlled solid-phase microextraction coupled with HPLC for the determination of angiotensin II receptor antagonists in human blood plasma

In this research, electrospun polycaprolactam nanofibers were collected on a fine stainless steel mesh sheet without a binder, and a layer of conductive polyaniline was chemically deposited on the nanofibers. The polyaniline immobilized on the polycaprolactam nanofibers provided high electrical conductivity, acceptable mechanical stability, and a large surface area. This assembly was then used as a working electrode in electrochemically controlled solid-phase microextraction (EC-SPME), a fast and environmentally friendly method. The polymer layers were characterized by SEM and FTIR techniques. Significant factors affecting the EC-SPME efficiency were investigated, including the desorption conditions, the sorbent used, the pH of the sample solution, the extraction voltage, the extraction time, and the ionic strength. Under the optimum conditions, the limits of detection and quantification for the target analytes were 0.9–1.8 μg L⁻¹ and 3.0–6.1 μg L⁻¹, respectively. The linear dynamic range was 5–2000 μg L⁻¹, with R² > 0.993. The method was coupled with HPLC analysis and applied to the determination of angiotensin ΙΙ receptor antagonists (ARA-ΙΙs) in human plasma, and relative recoveries of 91.1–104.3% with RSDs of ≤8.3% were obtained.

     

High-throughput approaches to the quantitative analysis of ketoconazole, a potent inhibitor of cytochrome P450 3A4, in human plasma

Ketoconazole, an imidazole-piperazine compound, is an orally active antimycotic agent. In addition, ketoconazole is a specific inhibitor of cytochrome P450 3A4. As about 60% of oxidized drugs are biotransformed by this isoform, the potential effect of a concomitant administration of ketoconazole on drug disposition may be of interest during drug development. The present paper describes three different approaches (methods A, B, and C) to attain high-throughput sample preparation and analysis in the quantification of ketoconazole in human plasma. Method A consisted of acetonitrile precipitation in a 96-well plate, transfer of the supernatant via a Tomtec Quadra 96 Model 320, and subsequent injection onto a 50 x 4.6 mm (i.d.) Develosil Combi-RP-5 column (packed with C30 bonded silica particles). Method B consisted of an identical sample preparation to method A with the exception that a Michrom Magic Bullet(trade mark) column, 2.0 --> 0.50 mm (i.d., tapered bore) x25 mm length, was used. Lastly, in method C, a turbulent-flow chromatography (TurboFlow LC/APCI-MS/MS) module was used for the direct analysis of ketoconazole in human plasma. A Sciex API 3000 was used in methods A and B, while a Micromass Quattro LC was employed in method C. Based on the values obtained for the calibrator (standard) and quality control samples, all three protocols yielded satisfactory accuracy, precision, and reduced manual sample preparation time.     

LHMDS mediated tandem acylation-cyclization of 2-aminobenzenecarbonitriles with 2-benzymidazol-2-yl acetates: a short and efficient route to the synthesis of 4-amino-3-benzimidazol-2-ylhydroquinolin-2-ones

We herein describe the discovery of a mild, one-pot tandem acylation-cyclization for the synthesis of 4-amino-3-benzimidazol-2-ylhydroquinolin-2-ones from 2-aminobenzenecarbonitriles and ethyl 2-benzimidazol-2-yl acetates.     

Binding of the Zn2+ ion to ferric uptake regulation protein from E. coli and the competition with Fe2+ binding: a molecular modeling study of the effect on DNA binding and conformational changes of Fur

The three dimensional structure of Ferric uptake regulation protein dimer from E. coli, determined by molecular modeling, was docked on a DNA fragment (iron box) and Zn²⁺ ions were added in two steps. The first step involved the binding of one Zn²⁺ ion to what is known as the zinc site which consists of the residues Cys 92, Cys 95, Asp 137, Asp141, Arg139, Glu 140, His 145 and His 143 with an average metal-Nitrogen distance of 2.5 Å and metal-oxygen distance of 3.1–3.2 Å. The second Zn²⁺ ion is bound to the iron activating site formed from the residues Ile 50, His 71, Asn 72, Gly 97, Asp 105 and Ala 109. The binding of the second Zn²⁺ ion strengthened the binding of the first ion as indicated by the shortening of the zinc-residue distances. Fe²⁺, when added to the complex consisting of 2Zn²⁺/Fur dimer/DNA, replaced the Zn²⁺ ion in the zinc site and when a second Fe²⁺ was added, it replaced the second zinc ion in the iron activating site. The binding of both zinc and iron ions induced a similar change in Fur conformations, but shifted residues closer to DNA in a different manner. This is discussed along with a possible role for the Zn²⁺ ion in the Fur dimer binding of DNA in its repressor activity. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.