Virtual screening strategies in drug discovery: a critical review

Virtual screening (VS) is a powerful technique for identifying hit molecules as starting points for medicinal chemistry. The number of methods and softwares which use the ligand and target-based VS approaches is increasing at a rapid pace. What, however, are the real advantages and disadvantages of the VS technology and how applicable is it to drug discovery projects? This review provides a comprehensive appraisal of several VS approaches currently available. In the first part of this work, an overview of the recent progress and advances in both ligand-based VS (LBVS) and structure-based VS (SBVS) strategies highlighting current problems and limitations will be provided. Special emphasis will be given to in silico chemogenomics approaches which utilize annotated ligand-target as well as protein-ligand interaction databases and which could predict or reveal promiscuous binding and polypharmacology, the knowledge of which would help medicinal chemists to design more potent clinical candidates with fewer side effects. In the second part, recent case studies (all published in the last two years) will be discussed where the VS technology has been applied successfully. A critical analysis of these case studies provides a good platform in order to estimate the applicability of various VS strategies in the new lead identification and optimization.

New 2-(aryloxy)-3-phenylpropanoic acids as peroxisome proliferator-activated receptor α/γ dual agonists able to upregulate mitochondrial carnitine shuttle system gene expression

The preparation of a series of 2-(aryloxy)-3-phenylpropanoic acids, resulting from the introduction of different substituents into the biphenyl system of the previously reported peroxisome proliferator-activated receptor α/γ (PPARα/γ) dual agonist 1, allowed the identification of new ligands with higher potency on PPARα and fine-tuned moderate PPARγ activity. For the most promising stereoisomer (S)-16, X-ray and calorimetric studies in PPARγ revealed, at high ligand concentration, the presence of two molecules simultaneously bound to the receptor. On the basis of these results and docking experiments in both receptor subtypes, a molecular explanation was provided for its different behavior as a full and partial agonist of PPARα and PPARγ, respectively. The effects of (S)-16 on mitochondrial acylcarnitine carrier and carnitine-palmitoyl-transferase 1 gene expression, two key components of the carnitine shuttle system, were also investigated, allowing the hypothesis of a more beneficial pharmacological profile of this compound compared to the less potent PPARα agonist fibrates currently used in therapy.

CDC25 phosphatase inhibitors: an update

The cell division cycle 25 (CDC25) family of proteins is a group of highly conserved dual-specificity phosphatases. They are key regulators of normal cell division and the cell response to DNA damage, and play a fundamental role in transitions between cell cycle phases during normal cell division, via the activation of CdK/cyclin complexes. Their abnormal expression, detected in a number of tumors, often correlated with a poor clinical prognosis, implies that their dysregulation is involved in malignant transformation. Thus, inhibition of these proteins represents an attractive therapeutic target in oncology, as evidenced from many patents and papers published on the subject in recent years. Hence, this review aims to provide an overview of recent developments in the field of CDC25 phosphatase inhibitor design since 2008.

Frontal affinity chromatography with MS detection of the ligand binding domain of PPARγ receptor: ligand affinity screening and stereoselective ligand-macromolecule interaction

In this study we report the development of new chromatographic tools for binding studies based on the gamma isoform ligand binding domain (LBD) of peroxisome proliferator-activated receptor (PPARγ) belonging to the nuclear receptor superfamily of ligand-activated transcription factors. PPARγ subtype plays important roles in the functions of adipocytes, muscles, and macrophages with a direct impact on type 2 diabetes, dyslipidemia, atherosclerosis, and cardiovascular disease. In order to set up a suitable immobilization chemistry, the LBD of PPARγ receptor was first covalently immobilized onto the surface of aminopropyl silica particles to create a PPARγ-Silica column for zonal elution experiments and then onto the surface of open tubular (OT) capillaries to create PPARγ-OT capillaries following different immobilization conditions. The capillaries were used in frontal affinity chromatography coupled to mass spectrometry (FAC-MS) experiments to determine the relative binding affinities of a series of chiral fibrates. The relative affinity orders obtained for these derivatives were consistent with the EC(50) values reported in literature. The optimized PPARγ-OT capillary was validated by determining the K(d) values of two selected compounds. Known the role of stereoselectivity in the binding of chiral fibrates, for the first time a detailed study was carried out by analysing two enantioselective couples on the LBD-PPARγ capillary by FAC and a characteristic two-stairs frontal profile was derived as the result of the two saturation events. All the obtained data indicate that the immobilized form of PPARγ-LBD retained the ability to specifically bind ligands.

[1,2,4]Triazino[4,3-a]benzimidazole acetic acid derivatives: a new class of selective aldose reductase inhibitors

Acetic acid derivatives of [1,2,4]triazino[4,3-a]benzimidazole (TBI) were synthesized and tested in vitro and in vivo as a novel class of aldose reductase (ALR2) inhibitors. Compound 3, (10-benzyl[1,2,4]triazino[4,3-a]benzimidazol-3,4(10H)-dion-2-yl)acetic acid, displayed the highest inhibitory activity (IC(50) = 0.36 microM) and was found to be effective in preventing cataract development in severely galactosemic rats when administered as an eyedrop solution. All the compounds investigated were selective for ALR2, since none of them inhibited appreciably aldehyde reductase, sorbitol dehydrogenase, or glutathione reductase. The activity of 3 was lowered by inserting various substituents on the pendant phenyl ring, by shifting the acetic acid moiety from the 2 to the 3 position of the TBI nucleus, or by cleaving the TBI system to yield benzimidazolylidenehydrazines as open-chain analogues. A three-dimensional model of human ALR2 was built, taking into account the conformational changes induced by the binding of inhibitors such as zopolrestat, to simulate the docking of 3 into the enzyme active site. The theoretical binding mode of 3 was fully consistent with the structure-activity relationships in the TBI series and will guide the design of novel ALR2 inhibitors.

Structure-based design, synthesis, and biological evaluation of conformationally restricted novel 2-alkylthio-6-[1-(2,6-difluorophenyl)alkyl]-3,4-dihydro-5-alkylpyrimidin-4(3H)-ones as non-nucleoside inhibitors of HIV-1 reverse transcriptase

5-Alkyl-2-(alkylthio)-6-(2,6-difluorobenzyl)-3,4-dihydropyrimidin-4(3H)-ones (S-DABOs, 2) have been recently described as a new class of human immunodeficiency virus type 1 (HIV-1) non-nucleoside reverse transcriptase (RT) inhibitors (NNRTIs) active at nanomolar concentrations (Mai, A. et al. J. Med. Chem. 1999, 42, 619-627). In pursuing our lead optimization efforts, we designed novel conformationally restricted S-DABOs, 3, featuring a methyl at the benzylic carbon (Y = Me) and at the pyrimidine 5-position (R = Me). Conformational analyses and docking simulations suggested that the presence of both methyls would significantly reduce conformational flexibility without compromising, in the R enantiomers, the capability of fitting into the RT non-nucleoside binding pocket. To develop structure-activity relationships, we prepared several congeners of type 3 belonging to the thymine (R = Me) and uracil (R = H) series, featuring various 2-alkylthio side chains (X = Me, i-Pr, n-Bu, i-Bu, s-Bu, c-pentyl, and c-hexyl) and aryl moieties different from the 2,6-difluorophenyl (Ar = phenyl, 2,6-dichlorophenyl, 1-naphthyl). Moreover, alpha-ethyl derivatives (Y = Et) were included in the synthetic project in addition to alpha-methyl derivatives (Y = Me). All of the new compounds were evaluated for their cytotoxicity and anti-HIV-1 activity in MT-4 cells, and some of them were assayed against highly purified recombinant wild-type HIV-1 RT using homopolymeric template primers. The results were expressed as CC(50) (cytotoxicity), EC(50) (anti-HIV-1 activity), SI (selectivity, given by the CC(50)/EC(50) ratio), and IC(50) (RT inhibitory activity) values. In the 2,6-difluorobenzylthymine (R = Me) series, methylation of the benzylic carbon improved anti-HIV-1 and RT inhibitory activities together with selectivity. Compound 3w (Ar = 2,6-F(2)-Ph, R = Y = Me, X = c-pentyl) turned out the most potent and selective among the S-DABOs reported to date (CC(50) > 200 microM, EC(50) = 6 nM, IC(50) = 5 nM, and SI > 33 333). Assays performed on the pure enantiomer (+)-3w, much more active than (-)-3w, yielded the following results: CC(50) > 200 microM, EC(50) = 2 nM, IC(50) = 8 nM, and SI > 100 000, under conditions wherein MKC-442 was less active and selective (CC(50) > 200 microM, EC(50) = 30 nM, IC(50) = 40 nM, SI > 6666). The 2,6-difluorophenylethylthymines (R = Me) were generally endowed with higher potency compared with the uracil counterparts (R = H). In the 2,6-difluorophenyl series the best and the least performant 2-alkylthio side chains were the 2-c-pentylthio and the 2-methylthio, respectively. When the methyl at the benzylic carbon was replaced by an ethyl, activity was retained or decreased slightly, thus suggesting that the dimensions of the cavity within the RT hosting this substituent would not be compatible with groups larger than ethyl. Aryl moieties different from the 2,6-difluorophenyl (phenyl, 1-naphthyl, 2,6-dichlorophenyl) were generally detrimental to activity, consistent with a favorable electronic effect exerted by the 2,6-fluorines on a putative charge-transfer interaction between the aromatic moieties of the inhibitor and Tyr188.

Novel N-(arylalkyl)indol-3-ylglyoxylylamides targeted as ligands of the benzodiazepine receptor: synthesis, biological evaluation, and molecular modeling analysis of the structure-activity relationships

A series of N-(arylalkyl)indol-3-ylglyoxylylamides (4-8) was synthesized as ligands of the benzodiazepine receptor (BzR) and tested for their ability to displace [(3)H]flumazenil from bovine brain membranes. The new compounds, bearing a branched (4) or a geometrically constrained benzyl/phenylethyl amide side chain (5-8), represent the continuation of our research on N-benzylindol-3-ylglyoxylylamides 1 (Da Settimo et al., 1996), N'-phenylindol-3-ylglyoxylohydrazides 2 (Da Settimo et al., 1998), and N-(indol-3-ylglyoxylyl)alanine derivatives 3 (Primofiore et al., 1989). A few indoles belonging to the previously investigated benzylamides 1 and phenylhydrazides 2 were synthesized and tested to enrich the SARs in these two series. The affinities and the GABA ratios of selected compounds for clonal mammalian alpha(1)beta(2)gamma(2), alpha(3)beta(2)gamma(2), and alpha(5)beta(3)gamma(2) BzR subtypes were also determined. It was hypothesized that the reduced flexibility of indoles 4-8 would both facilitate the mapping of the BzR binding cleft and increase the chances of conferring selectivity for the considered receptor subtypes. In the series of indoles 4, the introduction of a methyl group on the benzylic carbon with the R configuration improved affinity of the 5-substituted (5-Cl and 5-NO(2)) derivatives, whereas it was detrimental for their 5-unsubtituted (5-H) counterparts. All S enantiomers were less potent than the R ones. Replacement of the methyl with hydrophilic substituents on the benzylic carbon lowered affinity. The isoindolinylamide side chain was tolerated if the 5-position was unsubstituted (K(i) of 5a = 123 nM), otherwise affinity was abolished (5b, c). All the 2-indanylamides 6 and (S)-1-indanylamides 8 were devoid of any appreciable affinity. The 5-Cl and 5-NO(2) (R)-1-indanylamides 7b (K(i) 80 nM) and 7c (K(i) 28 nM) were the most potent among the indoles 5-8 geometrically constrained about the side chain. The 5-H (R)-1-indanylamide 7a displayed a lower affinity (K(i) 675 nM). The SARs developed from the new compounds, together with those collected from our previous studies, confirmed the hypothesis of different binding modes for 5-substituted and 5-unsubstituted indoles, suggesting that the shape of the lipophilic pocket L(1) (notation in accordance with Cook's BzR topological model) is asymmetric and highlighted the stereoelectronic and conformational properties of the amide side chain required for high potency. Several of the new indoles showed selectivity for the alpha(1)beta(2)gamma(2) subtype compared with the alpha(3)beta(2)gamma(2) and alpha(5)beta(3)gamma(2) subtypes (e.g.: 4t and 7c bind to these three BzR isoforms with K(i) values of 14 nM, 283 nM, 239 nM, and 9 nM, 1960 nM, 95 nM, respectively). The GABA ratios close to unity exhibited by all the tested compounds on each BzR subtype were predictive of an efficacy profile typical of antagonists.

3-Aryl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-ones: a new class of selective A1 adenosine receptor antagonists

Radioligand binding assays using bovine cortical membrane preparations and biochemical in vitro studies revealed that various 3-aryl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one (ATBI) derivatives, previously reported by us as ligands of the central benzodiazepine receptor (BzR) (Primofiore, G.; et al. J. Med. Chem. 2000, 43, 96-102), behaved as antagonists at the A1 adenosine receptor (A1AR). Alkylation of the nitrogen at position 10 of the triazinobenzimidazole nucleus conferred selectivity for the A1AR vs the BzR. The most potent ligand of the ATBI series (10-methyl-3-phenyl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one 12) displayed a Ki value of 63 nM at the A1AR without binding appreciably to the adenosine A2A and A3 nor to the benzodiazepine receptor. Pharmacophore-based modeling studies in which 12 was compared against a set of well-established A1AR antagonists suggested that three hydrogen bonding sites (HB1 acceptor, HB2 and HB3 donors) and three lipophilic pockets (L1, L2, and L3) might be available to antagonists within the A1AR binding cleft. According to the proposed pharmacophore scheme, the lead compound 12 engages interactions with the HB2 site (via the N2 nitrogen) as well as with the L2 and L3 sites (through the pendant and the fused benzene rings). The results of these studies prompted the replacement of the methyl with more lipophilic groups at the 10-position (to fill the putative L1 lipophilic pocket) as a strategy to improve A1AR affinity. Among the new compounds synthesized and tested, the 3,10-diphenyl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one (23) was characterized by a Ki value of 18 nM which represents a 3.5-fold gain of A1AR affinity compared with the lead 12. A rhodopsin-based model of the bovine adenosine A1AR was built to highlight the binding mode of 23 and two well-known A1AR antagonists (III and VII) and to guide future lead optimization projects. In our docking simulations, 23 receives a hydrogen bond (via the N1 nitrogen) from the side chain of Asn247 (corresponding to the HB1 and HB2 sites) and fills the L1, L2, and L3 lipophilic pockets with the 10-phenyl, 3-phenyl, and fused benzene rings, respectively.

Computer-assisted design, synthesis and biological evaluation of novel pyrrolyl heteroaryl sulfones targeted at HIV-1 reverse transcriptase as non-nucleoside inhibitors

Three pyrrolyl heteroaryl sulfones (ethyl 1-[(1H-benzimidazol-2(3H)one-5-yl)sulfonyl]-1H-pyrrole-2-carboxyla te, ethyl 1-[(1H-benzimidazol-5(6)-yl)sulfonyl]-1H-pyrrole-2-carboxylate and ethyl 1-[(1H-benzotriazol-5(6)-yl)sulfonyl]-1H-pyrrole-2-carboxylate) were designed as novel HIV-1 reverse transcriptase non-nucleoside inhibitors using structure-based computational methods. Although these compounds were inactive in the cell-based assay, they inhibited the target enzyme with micromolar potency (IC50s = 2 microM, 3 microM and 9 microM, respectively).

Synthesis, molecular modeling, and opioid receptor affinity of 9, 10-diazatricyclo[4.2.1.1(2,5)]decanes and 2,7-diazatricyclo[4.4.0. 0(3,8)]decanes structurally related to 3,8-diazabicyclo[3.2. 1]octanes

Various lines of evidence, including molecular modeling studies, imply that the endoethylenic bridge of 3,8-diazabicyclo[3.2. 1]octanes (DBO, 1) plays an essential role in modulating affinity toward mu opioid receptors. This hypothesis, together with the remarkable analgesic properties observed for N(3) propionyl, N(8) arylpropenyl derivatives (2) and of the reverted isomers (3), has prompted us to insert an additional endoethylenic bridge on the piperazine moiety in order to identify derivatives with increased potency toward this receptor class. In the present report, we describe the synthesis of the novel compounds 9,10-diazatricyclo[4.2. 1.1(2,5)]decane (4) and 2,7-diazatricyclo[4.4.0.0(3,8)]decane (5), as well as the representative derivatives functionalized at the two nitrogen atoms by propionyl and arylpropenyl groups (6a-e, 7a-d). Opioid receptor binding assays revealed that, among the compounds tested, the N-propionyl-N-cinnamyl derivatives 6a and 7a exhibited the highest mu-receptor affinity, and remarkably, compound 7a displayed in vivo (mice) an analgesic potency 6-fold that of morphine.

Modeling of kappa-opioid receptor/agonists interactions using pharmacophore-based and docking simulations

The interaction of the kappa-opioid receptor with arylacetamide and benzomorphan derivatives acting as agonists was modeled through pharmacophore-based and docking calculations. Potentially bioactive conformations of representative ligands (U-50,488 and its benzo-fused analogues 4 and 6 for arylacetamides and MPCB for benzomorphans) were identified by systematic conformational analysis and docked into a 3D model of the kappa-receptor. The obtained complexes, refined by energy-minimization and molecular dynamics, were evaluated for their consistency with structure-activity relationships and site-directed mutagenesis data. The following interactions are hypothesized to govern the ligand-receptor recognition process: (i) a salt bridge between the Asp138 carboxylate and the protonated nitrogen of the bound agonist; (ii) a hydrogen bond donated by the Tyr312 hydroxyl to the carbonyl oxygen of arylacetamides and MPCB; (iii) hydrophobic interactions established by the dichlorophenyl moiety of arylacetamides and the pendant phenyl ring of MPCB with the surrounding side chains of Tyr312, Leu224, Leu295, and Ala298; (iv) a pi-stacking contact between the Tyr312 side chain and the phenyl ring of arylacetamides; (v) a hydrogen bond linking the His291 imidazole ring to the phenolic hydroxy group featured by typical benzomorphans and the arylacetamides 4 and 6.

Structure-activity relationship studies on potential non-nucleoside DABO-like inhibitors of HIV-1 reverse transcriptase

Using 2,6-dichloro-4-aminopyrimidine, a number of uracil and cytosine derivatives with both arylthio and alkoxy moieties were prepared. These novel pyrimidines share chemical similarities with DABOs and HEPTs, two classes of non-nucleoside human immunodeficiency virus type 1 (HIV-1) reverse transcriptase inhibitors (NNRTIs), which have been widely studied of late. All new derivatives were tested in MT-4 cells to explore their potential in vivo anti-HIV activity. Like other NNRTIs, they selectively inhibit HIV-1 but not HIV-2. The majority of test derivatives were found to have low potency and were sometimes more cytotoxic than zidovudine and emivirine (formerly MKC-442), used here as reference drugs. Uracil and cytosine derivatives bearing a sec-butoxy chain and a methyl-substituted benzenesulphonyl moiety were the most potent. Enzyme assays proved that these derivatives target RT. Structure-activity relationship studies established a correlation between the anti-HIV-1 activity and the meta substitution on the phenyl ring; furthermore, oxidation of sulphide to sulphone significantly increased the potency of certain derivatives.

Solution conformation of a potent cyclic analogue of tuftsin: low- temperature nuclear magnetic resonance study in a cryoprotective mixture

Tuftsin, a linear tetrapeptide (Thr-Lys-Pro-Arg), corresponding to the sequence 289-292 of the heavy chain of leukokinin, has been the object of intensive SAR studies during the past 30 years, owing to its numerous biological activities and to the possibility of generating a novel anticancer drug. A cyclic tuftsin analogue, c-[T-K-P-R-G], has biological activity 50 times higher than that of the parent linear peptide. Here we present a conformational study of c-[T-K-P-R-G] based on NMR data in a cryoprotective DMSO/water mixture. The preferred conformation is a type VIa turn centered on the K-P residues. The orientation of the side chains of the two basic residues (K and R) may represent the essential feature of the bioactive conformation of tuftsin. A possible role of tuftsin as a DNA binding motif is suggested by the similarity of the bioactive conformation of c-[T-K-P-R-G] and of the beta-turn conformation proposed by Suzuki for the [T,S]-P-K-R motif.

Urachal adenocarcinoma: Case report

The Authors describe a case of urachal adenocarcinoma. The patient was a 36-year-old man who had already undergone segmentai (cupula) cystectomy at another hospital. The result of the histological test was: urachal adenocarcinoma. Afterwards a C.T. scan was taken which was negative. Then the patient had a cystoscopy with biopsy (negative) in our hospital, and underwent a new segmentai cystectomy releasing the urachus. The result of the second histological test was urachal adenocarcinoma. The Authors underline the rareness of this pathology, and the importance of total or segmentai cystectomy for treatment of the above.