Search

Cyclin-Dependent Kinases-2 (CDK2) are members of the serine/threonine protein kinases family. They play an important role in the regulation events of the eukaryotic cell division cycle, especially during the G1 to S phase transition. Experimental evidence indicates that excessive expression of CDK2s should cause abnormal cell cycle regulation. Therefore, since a long time, CDK2s have been considered potential therapeutic targets for cancer therapy. In this work, we collected one-hundred and forty-nine complexes of inhibitors bound in the CDK2-ATP pocket submitting to short MD simulations (10ns) and free energy calculation by means of MM-GBSA. The calculate ∆G values have been compared with experimental data (Ki, Kd, and IC50). Information collected on short MD simulations of protein-ligand complexes has been used to perform molecular modeling approaches that incorporates flexibility into structure-based pharmacophore modeling (Common Hits Approach, CHA,1 and Molecular dYnamics SHAred PharmacophorE, MYSHAPE2 approach) and constraints docking, to enrich the hits list of virtual screening. Short simulations proved to be exhaustive to examine the crucial ligand-protein interactions within the complexes.

Adenosine receptors (AR) belong to the superfamily of G-protein-coupled receptors (GPCRs). They are divided into four subtypes (A1, A2A, A2B, and A3) [1], and can be distinguished on the basis of their distinct molecular structures, distinct tissues distribution, and selectivity for adenosine analogs [2,3]. The hA3R, the most recently identified adenosine receptor, is involved in a variety of intracellular signaling pathways and physiological functions [4]. Expression of A3R was reported to be elevated in cancerous tissues [5], and A3 antagonists have been proposed for therapeutic treatments of cancer. The recent literature availability of crystal structure of hA2A adenosine receptor (PDB code: 3EML) provided us a new template for A3R homology modeling. The validation of the obtained structure model was performed by inspecting the Ramachandran plot (Fig. 1). The modeled protein was optimized using nanosecond scale molecular dynamics simulation. One hundred twenty two active and selective compounds were docked into the obtained model using Induced Fit Docking [6] and used as training set to generate pharmacophore models by means PHASE [7]. Energy-optimized pharmacophore mapping was performed; to each pharmacophore feature site was assigned an energetic value as the sum of the GLIDE XP contributions of the atoms included in the site. This pharmacophore model addresses the prevalent features to be used for the search of new inhibitors. Therefore it was employed as template to screen the ZINC database in the attempt to find new potent and selective human A3R antagonists. [1] B.B. Fredholm, A.P. Ijzerman, K.A. Jacobson, K.N. Klotz, J. Linden, Pharmacol. Rev., 53, 2001, 527. [2] P.G. Baraldi, R. Romagnoli, D. Preti, F. Fruttarolo, M.D. Carrion, M.A. Tabrizi, Curr. Med. Chem., 13, 2006, 3467. [3] M.D.Okusa, Am. J. Physiol. Renal Physiol., 282, 2002, F10. [4] A.Ochaion, S. Bar Yehuda, S. Cohen, F. Barer, R. Patoka, L. Del-Valle, G. Perez-Liz, J. Ophir, R. Galili-Mosberg, T. Reitblat, H. Amital, P. Fishman, Ann. Rheum. Dis., 66, 2007, 446. [5] L.Madi, A. Ochaion, L. Rath-Wolfson, S. Bar-Yehuda, A. Erlanger, G. Ohana, A. Harish, O. Merimski, F. Barer, P. Fishman, Clin. Cancer Res., 10, 2004, 4472. [6] Schrödinger Suite 2009 Induced Fit Docking protocol; Glide version 5.5, Schrödinger, LLC, New York, NY, 2009; Prime version 2.1, Schrödinger, LLC, New York, NY, 2009. [7] Phase, version 3.1, Schrödinger, LLC, New York, NY, 2009.

Zinc(ii) complexes are effective and selective nucleic acid-binders and strongly fluorescent molecules in the low energy range, from the visible to the near infrared. These two properties have often been exploited to quantitatively detect nucleic acids in biological samples, in both in vitro and in vivo models. In particular, the fluorescent emission of several zinc(ii) complexes is drastically enhanced or quenched by the binding to nucleic acids and/or upon visible light exposure, in a different fashion in bulk solution and when bound to DNA. The twofold objective of this perspective is (1) to review recent utilisations of zinc(ii) complexes as selective fluorescent probes for nucleic acids and (2) to highlight their novel potential applications as diagnostic tools based on their photophysical properties. This journal is

We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hit-list. The score for a particular molecule is calculated based on the number of representative pharmacophore models which classified it as active. Hence, the method is called common hits approach (CHA). The steps between the MD simulation and the final hit-list are performed automatically and without user interaction. We test the performance of CHA for virtual screening using screening databases with active and inactive compounds for 40 protein-ligand systems. The results of the CHA are compared to the (i) median screening performance of all representative pharmacophore models of protein-ligand systems, as well as to the virtual screening performance of (ii) a random classifier, (iii) the pharmacophore model derived from the experimental structure in the PDB, and (iv) the representative pharmacophore model appearing most frequently during the MD simulation. For the 34 (out of 40) protein-ligand complexes, for which at least one of the approaches was able to perform better than a random classifier, the highest enrichment was achieved using CHA in 68% of the cases, compared to 12% for the PDB pharmacophore model and 20% for the representative pharmacophore model appearing most frequently. The availabilithy of diverse sets of different pharmacophore models is utilized to analyze some additional questions of interest in 3D pharmacophore-based virtual screening.

In efforts to find agents with improved biological activity against cancer cells, recent years have seen an increased interest in the study of small molecules able to bind the deoxyribonucleic acid (DNA) when it assumes secondary structures known as G-quadruplexes (G4s) preferring them over the B form. Currently, several compounds reported in literature have already shown to be good candidates as G4s DNA stabilizers. Even though some specific features for the G4s affinity are known, such as a π-delocalized system able to stack at the top/end of a G-tetrad and positively charged substituents able to interact with the grooves, it is not clear yet what kind of structural features affect more the G4 arrangement. This is mainly due to the structure heterogeneity of both the G4 stabilizer compounds and the DNA G4s isoforms. In this review, we aim to classify some known G4 binders by analyzing them from a new perspective surprisingly never approached up to date: the symmetry features. Molecular symmetry could be responsible for the specific binding mode to the G4-DNA but could also be crucial in determining different isoform affinity. We propose to classify the G4s stabilizers in five main point group symmetry classes. This classification could be useful to design new ligands able to stabilize a specific G-quadruplex isoform, in order to increase the selectivity of new potential anticancer G-quadruplex targeting drugs, a goal yet highly sought by researchers.

The affinity of three square-planar nickel(II) (1), copper(II) (2) and zinc(II) (3) Schiff-base complexes for wild-type human telomeric (h-Telo) and protooncogene c-myc G-quadruplex (G4) DNA was investigated by UV-visible absorption spectroscopy and circular dichroism. DNA-binding constants (Kb) were determined by spectrophotometric titrations for both G4-DNA and B-DNA. The results obtained point out that the three metal complexes selectively bind G4-DNA with higher affinity, up to two orders of magnitude, with respect to B-DNA. The nickel(II) complex 1 was found to be the most effective G4-DNA stabilizer and the Kb values decrease in the order 1 > 2 ≈ 3. Innovative computational investigations, consisting of molecular dynamics (MD) simulations followed by density functional theory/molecular mechanics (DFT/MM) calculations, provide atomistic support for the interpretation of the binding mechanism to G4-DNA by end stacking and also of the experimental affinity order. Interestingly, 1 is able to induce G4-DNA formation of h-Telo sequences, also in the absence of K+ cations. This last result is nicely confirmed and highlighted by polymerase chain reaction (PCR) stop assays, which show the ability of the title compounds to induce and stabilize G4 structures inhibiting the amplification of PCR products. Finally, compounds 1–3 showed concentration and time-dependent cytotoxicity towards HeLa and MCF-7 human cancer cell lines, inducing significant effects on cell cycle distribution with G2/M arrest in HeLa cells and G0/G1 arrest in MCF-7 cells. Overall, the PCR inhibition and anticancer activity of the three compounds decreases in the same order 1 > 2 ≈ 3, in excellent correlation with the G4-DNA-binding affinity, implying that G4-DNA is the biotarget for their biological activity.

Heat shock protein 60 kDa (Hsp60) is a chaperone classically believed to be involved in assisting the correct folding of other mitochondrial proteins. Hsp60 also plays a role in cytoprotection against cell stressors, displaying for example, antiapoptotic potential. Despite the plethora of studies devoted to the mechanism of Hsp60's function, especially in prokaryotes, fundamental issues still remain unexplored, including the definition of its role in cancer. Key questions still unanswered pertain to the differences in structure-function features that might exist between the well-studied prokaryotic GroEL and the largely unexplored eukaryotic Hsp60 proteins. In this article we discuss these differences in sequence, structure, and roles of Hsp60, focusing on the human ortholog with the view of devising compounds to block its ability to favour tumor-cell growth and survival. Compounds currently known to directly or indirectly affect Hsp60 functions, such as protein folding, HIF-1α accumulation, or Hsp60-induced cell proliferation, are discussed along with strategies that might prove effective for developing Hsp60-targeting drugs for anticancer therapy.

The chemometric protocol VLAK was applied to predict improvement of the biological activity of pyrrolo-pyrimidine derivatives as anticancer agents, by using the NCI ACAM Database as depository of antitumor drugs with a known mechanism of action. Among the selected compounds two of these showed a good increase in the antitumor activity. These new pyrrolo-pyrimidine compounds were demonstrated effective against the full panels of NCI DTP tumour human cell lines. The derivative 8-[3-(piperidino)propyl]-4,10-dimethyl-9-phenyl-6-(methylsulfanyl)-3,4-dihydropyrimido[1,2-c]pyrrolo[3,2-e]pyrimidin-2(8H)-one reveled efficacious against the leukemia subpanel, in particular the RPMI cell line resulted the most sensitive (pGI50 = 6.68). Moreover the derivative 7-(3-Chloropropyl)-9-methyl-5-(methylsulfanyl)-8-phenyl-3H-imidazo[1,2-c]pyrrolo[3,2-e]pyrimidin-2(7H)-one showed a good antitumor activity against the leukemia subpanel with a low cytotoxic activity, above all against the HCT11 human tumour cell line. The VLAK protocol revealed a good method to design new molecules with good antitumor activity, starting from low active compounds. Moreover this protocol focused on the pyrrolo-pyrimidine derivatives as useful starting point for further development to obtain more potent antitumor agents.

p53 is a powerful anti-tumoral molecule frequently inactivated by mutations or deletions in cancer. However, half of all human tumors expresses wild-type p53, and its activation, by antagonizing its negative regulator Mdm2, might offer a new strategy for therapeutic protocol. In this work, we present a molecular dynamics study on Mdm2 structure bound to two different known inhibitors with the aim to investigate the structural transitions between apo-Mdm2 and Mdm2-inhibitor complexes. We tried to gain information about conformational changes binding a benzodiazepine derivative inhibitor with respect the known nutlin and the apo form. The conformational changes alter the size of the cleft and were mainly in the linker regions, suggesting that the overall dynamic nature of Mdm2 is related to dynamic movements in these regions.

Parasitic diseases caused by protozoarian agents are still relevant today more than ever. Recently, we synthesized several polycondensed diazine derivatives by means 1,3-dipolar cycloaddition reactions. A broad selection of these compounds were submitted to in vitro biological screening against Plasmodium falciparum , Leishmania infantum , Trypanosoma brucei , and Trypanosoma cruzi , resulting active at micromolar level. Induced Fit Docking/MM-GBSA studies were performed giving interesting indications about the probable mechanism of action of the most active compounds

The IKK kinases family represents a thrilling area of research because of its importance in regulating the activity of NF-kB transcription factors. The discovery of the central role played by IKK-β in the activation of transcription in response to apoptotic or inflammatory stimuli allowed to considerate its modulation as a promising tool for the treatment of chronic inflammation and cancer. To date, several IKK-β inhibitors have been discovered and tested. In this work, an analysis of the interactions between different classes of inhibitors and their biological target was performed, through the application of Molecular Docking and Pharmacophore/3D-QSAR approaches to a set of 141 inhibitors included in the Binding Database. In order to overcome the difficulty due to the lack of crystallographic data for IKK-β, a homology model of this protein has been built and validated. The results allowed to study in depth the structural bases for the interaction of each family of inhibitors and provided clues for further modifications, with the aim of improving the activity and selectivity of designed drugs targeting this enzyme.

The influence of the structural symmetry of the 2-pi double-reactive-sites component in the 1,3-dipolar cycloaddition reactions, involving nitrilimines as dipoles, was investigated. the experimental data showed that the loss of the symmetry leads to the formation of the monocycloadduct in good yields.

The reactivity of asymmetric benzo-condensed diazines in the 1,3-dipolar cycloaddition reactions with nitrilimines was investigated. The results demonstrated that, at variance with the symmetric quinoxaline, a certain grade of diastereoselectivity emerged. Moreover in the case of the 5-methylquinoxaline and quinazoline a mono-cycloadduct was obtained.

One of the major problems in the fight against cancer is drug-resistance, which, at a molecular level, can be acquired through mutations able to deactivate apoptosis. In particular, proteins in the Bcl-2 family are central regulators of programmed cell death, and members that inhibit apoptosis, such as Bcl-xl and Bcl-2, are overexpressed in many tumours. The development of new inhibitors of these proteins as potential anticancer therapeutics represents a new frontier. In this work, we carried out an in-silico screening of compounds from a free database of more than 2 million structures (ZINC database), which allowed us to identify 17 sulfonamide derivatives as new potential inhibitors; these are currently undergoing biological evaluation.

Various derivatives of the heterocyclic system 1,12,12a,12b-tetrahydrobis-1,2,4-triazolo[4,3-a:3′,4′-c]quinoxaline of pharmaceutical interest have been obtained by double site- and regio-selective 1,3-dipolar cycloaddition of arylnitrilimines to quinoxalines. No evidence for the formation of mono-adducts was obtained, at variance with literature reports. Specific studies to establish the exact stereochemistry of the bis-cycloadducts were undertaken.

In this paper we describe a comparative analysis between multivariate and docking methods in the study of the drug resistance to the reverse transcriptase and the protease inhibitors. In our early papers we developed a simple but efficient method to evaluate the features of compounds that are less likely to trigger resistance or are effective against mutant HIV strains, using the multivariate statistical procedures PCA and DA. In the attempt to create a more solid background for the prediction of susceptibility or resistance, we carried out a comparative analysis between our previous multivariate approach and molecular docking study. The intent of this paper is not only to find further support to the results obtained by the combined use of PCA and DA, but also to evidence the structural features, in terms of molecular descriptors, similarity, and energetic contributions, derived from docking, which can account for the arising of drug-resistance against mutant strains.

New classes of pyrrolopyrimidodiazepinone derivatives (PPDs) and their precursors were studied in silico for their ability to form stable complex with DNA fragment. In the docking studies two binding modes can be envisaged: groove mode and intercalating mode. In the case of the best ligands the docking studies revealed a common binding mode with the chromophore intercalated between GC base pairs whereas the side chain lies close to the minor groove. Synthetic approach to the PPD ring systems is discussed.

2-Amino-N-[2-(3,4-dihydroxy-phenyl)-ethyl]-3-phenyl-propionamide (DA-PHEN) has been previously synthesized to obtain a potential prodrug capable of release dopamine (DA) into CNS. However, DA-PHEN could act per se as a dopaminergic drug. In this study, the permeability transport (Pe), obtained by parallel artificial permeability assay (PAMPA), indicated a low passive transcellular transport (Pe = 0.32 ± 0.01 × 10(-6 )cm/s). Using the Caco-2 cell system, the Papp AP-BL in absorptive direction (3.36 ± 0.02 × 10(-5 )cm/s) was significantly higher than the Papp BL-AP in secretive direction (1.75 ± 0.07 × 10(-5 )cm/s), suggesting a polarized transport. The efflux ratio (Papp AP-BL/Papp BL-AP = 0.52 ± 0.02) indicated a low affinity of DA-PHEN to efflux carriers. The forced swim test highlighted a reduction of immobility time in both pre-test and test sessions (p 0.0001), with an exacerbation in the number of headshakes and divings in the pretest (p 0.0001). Morris water maze strengthened the hypothesis that DA-PHEN induces adaptive responses to environmental challenges which are involved on cognitive functions (DA-PHEN versus CTR: escape latency; p 0.001; distance swum p 0.001, time spent on target quadrant p 0.001), without any change in locomotor activity for the administered dose. The molecular docking revealed the interaction of DA-PHEN with the identified D1 site mapping human brain receptor.

The title indolo-triazolo-pyrimidines were obtained from 3-azidoindoles and can be used as models for the design of DNA-interactive compounds. Hetero-domino reaction of azidoindoles/pyrroles and acetonitriles constitutes the synthetic entry to annelated 1,2,3-triazolo[1,5- a ]pyrimidines.

The aim of this work was devoted to develop a method to predict Delta G values for a series of minor groove binders. Starting from a matrix of docking dataset for 10 minor groove binders (known and not) to 20 DNA fragments, with various sequences, it was possible to analyze the interaction modes and to calculate the Delta G value for new derivatives through MADoSPRO procedure. The method allowed, through the QSPR analysis, to characterize the type of interactions in such complexes, that was demonstrated to be related to quantum chemical and electrostatic descriptors, in agreement with the information available in literature on the structural requirements of specific minor groove ligands. Moreover it was possible also to design new ligands with different lengths and enedine spacers to modulate the binding affinity with various nucleotide sequences. It was observed that in some cases it was possible to increase the affinity toward GC base-pairs notoriously not favoured for most of groove binders.

New benzothieno[3,2-d]-1,2,3-triazines, together with precursors triazenylbenzo[b]thiophenes, were designed, synthesized and screened as anticancer agents. The structural features of these compounds prompted us to investigate their DNA binding capability through UV–vis absorption titrations, circular dichroism, and viscometry, pointing out the occurrence of groove-binding. The derivative 3-(4-methoxy-phenyl)benzothieno[3,2-d]-1,2,3-triazin-4(3H)-one showed the highest antiproliferative effect against HeLa cells and was also tested in cell cycle perturbation experiments. The obtained results assessed for the first time the anticancer activity of benzothieno[3,2-d]-1,2,3-triazine nucleus, and we related it to its DNA-binding properties.

Body weight control is a mechanism finely regulated by several hormonal, metabolic, and nervous pathways. The leptin receptor (Ob-R) is crucial for energy homeostasis and regulation of food uptake. Leptin is a 16 kDa hormone that is mainly secreted by fat cells into the bloodstream, and under normal circumstances, circulating levels are proportionate to the fat body mass. Sensing of elevated leptin levels by the hypothalamic neurocircutry activates a negative feedback loop resulting in reduced food intake and increased energy expenditure. Decreased concentrations lead to opposite effects. Therefore rational design of leptin agonists constitute an appealing challenge in the battle against obesity. In this study, we performed protein-protein docking among the re-built crystal structure of leptin and leptin binding domain (LBD). The obtained complex was used as a starting point to carry out nanosecond-scale molecular dynamics simulations to characterize the key regions in terms of physical-chemical features involved in the protein-protein interaction (dynamic site mapping filtered by means multivariate analysis) and used to carry out a HTVS. The main goal of this study was to suggest guidelines for the rational drug design of new agonists of leptin. Identified hits could be a consistent starting point to carry out in vitro testing.

The in silico COMPARE analysis was performed on 8-[3-(piperidino)propyl]-4,10-dimethyl-9-phenyl-6-(methylsulfanyl)-3,4-dihydropyrimido[1,2-c]pyrrolo[3,2-e]pyrimidin-2(8H)-one, a compound with promising antiproliferative activity, previously synthetized and screened against a panel of 60 human tumor cell lines. The results evidenced that this compound matches the biological properties of Chromomycin A3 and Actinomycin D, known drugs with high DNA binding affinity. Prompted by such results, a thorough spectroscopic investigation of its DNA aqueous solutions was performed, with the aim to verify its DNA-binding properties. DNA groove-binding interaction was assigned by UV-vis spectrophotometric and circular dichroism (CD) titrations. In agreement with the spectroscopic studies, the test compound has achieved significant cell-cycle perturbation with higher accumulation of cells in G0/G1 phase. Further structure modifications of the tested compound, has led to its analogue, which showed DNA base-pairs intercalating ability. Molecular modeling studies on the annelated pyrrolo[3,2-e]pyrimidines tested were in agreement with the biological evaluation.

The competitive binding of nickel(II), copper(II) and zinc(II) complexes toward B- and G4-DNA was addressed through spectroscopic titrations and rationalized by computational investigations, consisting of molecular dynamics simulations followed by density functional theory/molecular mechanics (DFT/MM) calculations [1]. The experimental DNA binding studies clearly highlight the selectivity of the compounds, in particular the nickel(II) complex, toward G4-DNA from both h-Telo and c-myc. Moreover, the compounds show biological activity against HeLa and MCF-7 cancer cell lines. Remarkably, the experimental DNA-binding affinity trend of the three metal complexes, obtained from the DNA-binding constants as ΔG°=−RT ln(Kb), is reproduced by the Gibbs formation free energy calculated by DFT/MM for the DNA-binding complexes, in the implicit water solution.

Nickel(II), copper(II) and zinc(II) complexes of N2O2 tetradentate Schiff base ligands strongly interact with B-DNA, usually by groove-binding and/or by intercalation [1]. It has been also shown that the presence of aromatic substituents on the N,N’ bridge make them suitable G-quadruplex binders [2]. In this context, we have recently investigated the binding toward duplex and G-quadruplex DNA of nickel(II), copper(II) and zinc(II) complexes of N,N’-bis-5-(triethyl ammonium methyl)-salicylidene-2,3-naphthalendiiminato) (see Figure), by spectroscopic and computational methods [3,4]. The compounds show also biological activity against human cancer cell lines. Different substituents are presently considered on the N,N’-bridge, in order to increase their selectivity towards telomeric and oncogene promoter DNA sequences, rather than toward B-DNA, and to target G-quadruplexes on the grooves rather than by π-π stacking.

Chaperonins mediate protein folding in a cavity formed by multisubunit rings. The human CCT has eight non-identical subunits and the His147Arg mutation in one subunit, CCT5, causes neuropathy. Knowledge is scarce on the impact of this and other mutations upon the chaperone's structure and functions. To make progress, experimental models must be developed. We used an archaeal mutant homolog and demonstrated that the His147Arg mutant has impaired oligomeric assembly, ATPase activity, and defective protein homeostasis functions. These results establish for the first time that a human chaperonin gene defect can be reproduced and studied at the molecular level with an archaeal homolog. The major advantage of the system, consisting of rings with eight identical subunits, is that it amplifies the effects of a mutation as compared with the human counterpart, in which just one subunit per ring is defective. Therefore, the slight deficit of a non-lethal mutation can be detected and characterized.

Recently, NiII and ZnII metal complexes of the ligand Salpyrim have been synthesized and characterized. Their affinity for wild-type h-Telo G-quadruplex DNA and for calf thymus DNA was investigated by UV absorption spectroscopy, circular dichroism and viscometry. The data collectively suggest that both complexes bind effectively to G-quadruplexes by direct end-stacking, stabilizing the oligonucleotide secondary structure. The two complexes are also typical B-DNA intercalators. Remarkably, their binding constants, Kb, with the G4s structures are about 10 fold higher than those with B-DNA, highlighting the selectivity. Experiments to evaluate the biological activity of the two complexes against MCF7 and HeLa cancer cell lines are currently ongoing.

1,3-Dipolar cycloaddition reactions can be considered a powerful synthetic tool in the building of heterocyclic rings, with applications in different fields. In this review we focus on the synthesis of biologically active compounds possessing the 1,2,3-triazole core through 1,3-dipolar cycloaddition reactions. The 1,2,3-triazole skeleton can be present as a single disubstituted ring, as a linker between two molecules, or embedded in a polyheterocycle. The cycloaddition reactions are usually catalysed by copper or ruthenium. Domino reactions can be achieved through dipolarophile anion formation, generally followed by cyclisation. The variety of attainable heterocyclic structures gives an illustration of the importance of the 1,2,3-triazole core in medicinal chemistry.

In the last years, the introduction of side chains in different molecular compounds takes increasing importance. As recently reported in literature, heterocyclic scaffolds with poor biological activity , if properly decorated with selected side chains, can improve their anticancer activity against a large spectrum of human tumor cell lines. For example, the annelated Pyrrolo[3,4-e]Pyrimidines and Pyrrolo[3,2-e]Pyrimidines, opportunely functionalized with a large number of side chains, showed a good increase in the antitumor activity with respect to the starting core structure. New compound thus obtained showed antiproliferative activity against all the human tumor cells, generally in the low micromolar concentration range (1). Furthermore, this increase in activity is confirmed in the new linear derivatives Thieno[3,2-d][1,2,3]Triazolo [1,5- α]Pyrimidines (2) and their angular isomer Thieno[2,3-e][1,2,3]Triazolo[1,5- α]Pyrimidines (3). Another confirmation is provided by the Indolo [3,2-e][1,2,3]Triazolo [1,5- α]Pyrimidines, that resulted inactive as antitumor agent, but when a functionalization with side chains takes place, the anticancer activity has interestingly appeared (4). Here, with the aim to explore the importance of the addition of side chains on smaller molecular scaffolds, of easier synthetic access and decoration, we focused our attention on a new series of substituted “chained” pyrrole derivatives in order to evaluate their anticancer potential.

We present a thorough characterization of the interaction of novel nickel(II) (1), copper(II) (2) and zinc(II) (3) Schiff base complexes with native calf thymus DNA (ct-DNA), in buffered aqueous solution at pH 7.5. UV-vis absorption, circular dichroism (CD) and viscometry titrations provided clear evidence of the intercalative mechanism of the three square-planar metal complexes, allowing us to determine the intrinsic DNA-binding constants (K(b)), equal to 1.3 × 10(7), 2.9 × 10(6), and 6.2 × 10(5) M(-1) for 1, 2 and 3, respectively. Preferential affinity, of one order of magnitude, toward AT compared to GC base pair sequences was detected by UV-vis absorption titrations of 1 with [poly(dG-dC)]2 and [poly(dA-dT)]2. Structural details of the intercalation site of the three metal complexes within [dodeca(dA-dT)]2 were obtained by molecular dynamics (MD) simulations followed by density functional theory/molecular mechanics (DFT/MM) calculations. The calculations revealed that three major intermolecular interactions contribute to the strong affinity between DNA and the three metal complexes: (1) the electrostatic attraction between the two positively charged triethylammoniummethyl groups of the metal complexes and the negatively charged phosphate groups of the DNA backbone; (2) the intercalation of the naphthalene moiety within the four nitrogen bases of the intercalation site; (3) the metal coordination by exocyclic donor atoms of the bases, specifically the carbonyl oxygen and amine nitrogen atoms. Remarkably, the Gibbs formation free energy calculated for the intercalation complexes of 1, 2 and 3 with [dodeca(dA-dT)]2 in the implicit water solution is in agreement with the experimental Gibbs free energy values obtained from the DNA-binding constants as ΔG° = -RT ln(K(b)). In particular, the DNA-binding affinity trend, 123, is reproduced. Finally, the first shell coordination distances calculated for the intercalation complex 3/[dodeca(dA-dT)]2 are in excellent agreement with the experimental distances extracted from the extended X-ray absorption fine structure (EXAFS) spectrum of the corresponding 3/ct-DNA solutions. The latter results provided the first evidence of metal ion coordination by native DNA in aqueous solution.

The binding of small molecules with biological targets is associated to interesting chemical and biological properties of the resulting supramolecular systems. We have recently reported on the synthesis and characterization of cationic first row transition metal complexes and the study of their DNA binding properties, in aqueous solutions at neutral pH, essentially investigated by viscosimetry and spectroscopic techniques such as circular dichroism, absorption and fluorescence in the UV-visible wavelength range. Of course, such procedure cannot furnish atomic level details of the molecule-DNA interaction. Computational Chemistry may provide support for the interpretation of experimental data on an atomistic level (Fig.1). For example, we have recently shown that Molecular Dynamics (MD) simulations, followed by quantum mechanics/molecular mechanics (QM/MM) calculation, provided detailed structural informations and binding energies of the complexes between nickel(II), copper(II), zinc(II) metallointercalators with nucleic acids in the canonical B conformation [1]. We are presently applying such complementary experimental and computational approach to the interaction of small molecules with G-quadruplex (G4) DNA. The latter is a non-canonical conformation recently observed in human cells [2], and it has been proposed as a target for a novel class of anticancer drugs. Recently we have performed MD simulation to have an insight into the molecular recognition process of small organic ligands and other biological targets, such as mRNA and proteins [3].

In this work the new protocol BIOlogical Target Assignation (BIOTA) for the prediction of the biological target from molecular structures is proposed. BIOTA is based on the Principal Components Analysis (PCA) application on a matrix of ligands versus molecular descriptors. The application of BIOTA could allow to hypothesize the mechanism of action of a candidate drug prior to its biological evaluation or to repurpose old drugs. The protocol can be fine-tuned by choosing opportune targets (biological or not) and molecular descriptors, and it can be useful in every fields in with it is possible to collect set of compounds with known properties. The robustness of the protocol depends from different factors: the correctness of biological data, the optimization of the molecular structures and their molecular descriptors calculation, the selection of the biological targets. The application of BIOTA to a new class of Hsp90 inhibitors was able to predict quite correctly their affinity for Hsp90.

The OB-receptor or leptin receptor (LR) is crucial for energy homeostasis and regulation of food uptake. Leptin is a 16 kDa hormone that is mainly secreted by fat cells into the bloodstream. Under normal circumstances, circulating leptin levels are proportionate to the fat body mass. Sensing of elevated leptin levels by the hypothalamic neuro-circuitry activates a negative feedback loop resulting in reduced food intake and increased energy expenditure. Decreased leptin concentrations lead to opposite effects. Therefore, rational design of leptin agonists/antagonists could be an appealing challenge in the battle against obesity. The Leptin/LR interactions have been studied in several works by means of different molecular modelling approaches, spreading from homology modelling to manual docking. No small molecules have ever been proposed as agonists of the Ob receptor but researchers’ efforts focused only on leptin-related synthetic peptides as receptor antagonists and on peptidomimetics. In this review we try to track a timeline of obtained in silico information to clarify the mechanism of interaction between leptin and its receptor, together to summarize the more recent efforts to propose new drugs usable in anti-obesity therapy. Final considerations could be useful starting points for the rational drug design of new lead compounds.

An unusual Dimroth rearrangement occurring in the reaction leading to annelated thieno[2,3-e][1,2,3]triazolo[1,5-a]pyrimidine core allowed the isolation of the linear isomer thieno[3,2-d][1,2,3]triazolo[1,5-a]pyrimidine. By decorating the linear isomer with the same chains that improved the biological activity of the angular isomers, new annelated thieno[3,2-d][1,2,3]triazolo[1,5-a]pyrimidines were designed and synthesized. They were selected by the Developmental Therapeutics Program (DTP) of the National Cancer Institute (NCI) for the anticancer screening against a panel of 60 human tumor cell lines. The biological results showed that the new derivatives exhibited strong antiproliferative activity up to nanomolar concentration. In vivo screenings of the most active compound, the N-[2-(1H-imidazol-4-yl)ethyl]-4-(3-phenyl-10-oxo-4,10-dihydrobenzothieno[3,2-d][1,2,3]triazolo[1,5-a]pyrimidin-4-yl)butanamide, showed its low toxicity and high potency.

First observed in cells exposed to high temperatures, Heat shock proteins (Hsps) are nowadays considered the most important cell “chaperone” complexes over-­expressed in response to a number of cell stress stimuli.1 The chaperone activity is the main function of the eukaryotic Heat shock protein 60 kDa (Hsp60), involved in the capture and refold of unfolded or misfolded proteins. Additional roles in signal transduction,2 senescence activation3 and apoptosis4 have been ascribed to cytosolic Hsp60. During the carcinogenIc process, in vivo studies demonstrated increased levels of human Hsp60 in several organs, such as uterine exocervix,5 large bowel,6 and prostate.6 In this context, our study aims to elucidate the structural details of the interaction between Hsp60 and novel designed antagonists able to specifically block this chaperonine. A preliminary virtual screening of 24 million molecules, available in the Zinc database, was carried out on the ATP-­binding site of a bacterial Hsp60 monomer, the coordinates of which were taken from Protein Data Bank (ID: 1WE3), figure 1. Compounds with high affinity were further refined by other in silico protocols previously and successfully applied by us in the study of several biological targets.7The analysis of virtual screening results highlights the N-­{5-­[1H-­imidazol-­4-­yl-­methyl)-­ amino]-­benzofuran-­3-­yl}-­benzamides of type 1 as interesting series for the inhibition of Hsp60 ATP-­binding site. Selected compounds were prepared in excellent yields, following appropriate synthetic pathways. All compounds are currently tested in order to proof they potential anticancer activity as modulator of Hsp60 function in tumor cells.

The reaction mechanism of 1,3-dipolar cycloadditions of both symmetric and unsymmetric benzo-condensed diazines with a nitrilimine dipole, to give two different mono- and bis-cycloadducts, in tetrahydrofuran (THF) solution, was studied by DFT calculatio ns. The results obtained show that each 1,3-dipolar cycloaddition reaction always proceeds by a two steps mechanism, in which the first intermediate shows only one covalent bond between the beta carbon of the nitrilimine and the aromatic nitrogen of the diazine molecule. The structure and energy content of the two transition states of the two cycloaddition steps, in the case of the unsymmetric benzo-condensed diazine, nicely explains why the product of the bis-cycloaddi tions is exclusively observed and why the product of a mono-cycloaddition is not isolated for the symmetric reaction pathway.