Sandro Cosconati

Virtual Screening with AutoDock: Theory and Practice.

Importance of the field: Virtual screening is a computer-based technique for identifying promising compounds to bind to a target molecule of known structure. Given the rapidly increasing number of protein and nucleic acid structures, virtual screening continues to grow as an effective method for the discovery of new inhibitors and drug molecules. Areas covered in this review: We describe virtual screening methods that are available in the AutoDock suite of programs and several of our successes in using AutoDock virtual screening in pharmaceutical lead discovery. What the reader will gain: A general overview of the challenges of virtual screening is presented, along with the tools available in the AutoDock suite of programs for addressing these challenges. Take home message: Virtual screening is an effective tool for the discovery of compounds for use as leads in drug discovery, and the free, open source program AutoDock is an effective tool for virtual screening.

Tandem Application of Virtual Screening and NMR Experiments in the Discovery of Brand New DNA Quadruplex Groove Binders

In the past decade, DNA G-quadruplexes have come into the limelight thanks to their biological implications and to their potential druggability in anticancer therapy. In particular, it has been found that small molecules that stabilize G-quadruplex structures are effective inhibitors of telomerase which plays a critical role in tumorigenesis. So far, the quadruplex groove recognition, which is expected to give a higher degree of selectivity over the other DNA structures, has been demonstrated for very few compounds. Thus with the aim of detecting new and structurally diverse groove binders, a structure-based virtual screening campaign has been performed using the X-ray structure of the [d(TGGGGT)](4) quadruplex. Remarkable results were achieved, and six brand new different molecular entities have been found to interact with the groove through NMR experiments. The reported results will certainly stimulate further studies aimed at the design and optimization of new quadruplex-specific groove binders to be applied as anticancer agents and for other diseases.

Sampling protein motion and solvent effect during ligand binding

An exhaustive description of the molecular recognition mechanism between a ligand and its biological target is of great value because it provides the opportunity for an exogenous control of the related process. Very often this aim can be pursued using high resolution structures of the complex in combination with inexpensive computational protocols such as docking algorithms. Unfortunately, in many other cases a number of factors, like protein flexibility or solvent effects, increase the degree of complexity of ligand/protein interaction and these standard techniques are no longer sufficient to describe the binding event. We have experienced and tested these limits in the present study in which we have developed and revealed the mechanism of binding of a new series of potent inhibitors of Adenosine Deaminase. We have first performed a large number of docking calculations, which unfortunately failed to yield reliable results due to the dynamical character of the enzyme and the complex role of the solvent. Thus, we have stepped up the computational strategy using a protocol based on metadynamics. Our approach has allowed dealing with protein motion and solvation during ligand binding and finally identifying the lowest energy binding modes of the most potent compound of the series, 4-decyl-pyrazolo[1,5-a]pyrimidin-7-one.

Design, Synthesis, and Functionalization of Dimeric Peptides Targeting Chemokine Receptor CXCR4

The chemokine receptor CXCR4 is a critical regulator of inflammation and immune surveillance, and it is specifically implicated in cancer metastasis and HIV-1 infection. On the basis of the observation that several of the known antagonists remarkably share a C(2) symmetry element, we constructed symmetric dimers with excellent antagonistic activity using a derivative of a cyclic pentapeptide as monomer. To optimize the binding affinity, we investigated the influence of the distance between the monomers and the pharmacophoric sites in the synthesized constructs. The affinity studies in combination with docking computations support a two-site binding model. In a final step, 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) was introduced as chelator for (radio-)metals, thus allowing to exploit these compounds as a new group of CXCR4-binding peptidic probes for molecular imaging and endoradiotherapeutic purposes. Both the DOTA conjugates and some of their corresponding metal complexes retain good CXCR4 affinity, and one (68)Ga labeled compound was studied as PET tracer.

Conformational control of integrin subtype selectivity in isoDGR peptide motifs: A biological switch

Thumbnail image of graphical abstract The rearrangement of asparagine to isoaspartate (isoD) is responsible for the deactivation of many functional proteins. However, the isoDGR motif, which is optimally presented as a conformationally controlled cyclic pentapeptide, binds selectively to a5s1 integrin (see the docking model) with an affinity comparable to that of the peptidic antitumor agent Cilengitide

Modeling of Cdc25B dual specifity protein phosphatase inhibitors: docking of ligands and enzymatic inhibition mechanism.

The Cdc25 dual specificity phosphatases have central roles in coordinating cellular signalling processes and cell proliferation. It has been reported that an improper amplification or activation of these enzymes is a distinctive feature of a number of human cancers, including breast cancers. Thus, the inhibition of Cdc25 phosphatases might provide a novel approach for the discovery of new and selective antitumor agents. By using the crystal structure of the catalytic domain of Cdc25B, structural models for the interaction of various Cdc25B inhibitors (1–13) with the enzyme were generated by computational docking. The parallel use of two efficient and predictive docking programs, AutoDock and GOLD, allowed mutual validation of the predicted binding poses. To evaluate their quality, the models were validated with known structure–activity relationships and site‐directed mutagenesis data. The results provide an improved basis for structure‐based ligand design and suggest a possible explanation for the inhibition mechanism of the examined Cdc25B ligands. We suggest that the recurring motif of a tight interaction between the inhibitor and the two arginine residues, 482 and 544, is of prime importance for reversible enzyme inhibition. In contrast, the irreversible inhibition mechanism of 1–4 seems to be associated with the close vicinity of the quinone ring and the Cys473 catalytic thiolate. We believe that this extensive study might provide useful hints to guide the development of new potent Cdc25B inhibitors as novel anticancer drugs.

Selective non-nucleoside inhibitors of human DNA methyltransferases active in cancer including in cancer stem cells.

DNA methyltransferases (DNMTs) are important enzymes involved in epigenetic control of gene expression and represent valuable targets in cancer chemotherapy. A number of nucleoside DNMT inhibitors (DNMTi) have been studied in cancer, including in cancer stem cells, and two of them (azacytidine and decitabine) have been approved for treatment of myelodysplastic syndromes. However, only a few non-nucleoside DNMTi have been identified so far, and even fewer have been validated in cancer. Through a process of hit-to-lead optimization, we report here the discovery of compound 5 as a potent non-nucleoside DNMTi that is also selective toward other AdoMet-dependent protein methyltransferases. Compound 5 was potent at single-digit micromolar concentrations against a panel of cancer cells and was less toxic in peripheral blood mononuclear cells than two other compounds tested. In mouse medulloblastoma stem cells, 5 inhibited cell growth, whereas related compound 2 showed high cell differentiation. To the best of our knowledge, 2 and 5 are the first non-nucleoside DNMTi tested in a cancer stem cell line.

Protein Flexibility in Virtual Screening: The BACE-1 Case Study

Simulating protein flexibility is a major issue in the docking-based drug-design process for which a single methodological solution does not exist. In our search of new anti-Alzheimer ligands, we were faced with the challenge of including receptor plasticity in a virtual screening campaign aimed at finding new β-secretase inhibitors. To this aim, we incorporated protein flexibility in our simulations by using an ensemble of static X-ray enzyme structures to screen the National Cancer Institute database. A unified description of the protein motion was also generated by computing and combining a set of grid maps using an energy weighting scheme. Such a description was used in an energy-weighted virtual screening experiment on the same molecular database. Assessment of the enrichment factors from these two virtual screening approaches demonstrated comparable predictive powers, with the energy-weighted method being faster than the ensemble method. The in vitro evaluation demonstrated that out of the 32 tested ligands, 17 featured the predicted enzyme inhibiting property. Such an impressive success rate (53.1%) demonstrates the enhanced power of the two methodologies and suggests that energy-weighted virtual screening is a more than valid alternative to ensemble virtual screening given its reduced computational demands and comparable performance.

Novel N2-substituted pyrazolo[3,4-d]pyrimidine adenosine A3 receptor antagonists: inhibition of A3-mediated human glioblastoma cell proliferation.

Adenosine induces glioma cell proliferation by means of an antiapoptotic effect, which is blocked by cotreatment with selective A(3) AR antagonists. In this study, a novel series of N(2)-substituted pyrazolo[3,4-d]pyrimidines 2a-u was developed as highly potent and selective A(3) AR antagonists. The most performing compounds were derivatives 2a (R(1) = CH(3) and R(2) = COC(6)H(5); K(i) 334, 728, and 0.60 nM at the human A(1), A(2A), and A(3) ARs, respectively) and 2b (R(1) = CH(3) and R(2) = COC(6)H(4)-4-OCH(3); K(i) 1037, 3179, and 0.18 nM at the human A(1), A(2A), and A(3) ARs, respectively), which counteracted the effect of the A(3) AR agonists Cl-IB-MECA and IB-MECA on human glioma U87MG cell proliferation. This effect was concentration-dependent, with IC(50) values comparable to A(3) AR binding affinity values of 2a and 2b, thereby suggesting that their effects were receptor-mediated. Furthermore, the antiproliferative activity of the new compounds was demonstrated to be mediated by the block of A(3) AR agonist activation of intracellular kinases ERK 1/2.

A Conformationally Frozen Peptoid Boosts CXCR4 Affinity and Anti‐HIV Activity

The chemokine receptor subtype CXCR4 belongs to the G-protein coupled receptors (GPCRs) and is, together with itsnatural ligand CXCL12 (or SDF-1), a central part of thesignaling system in the human body. Its functions range fromstem-cell trafficking during embryogenesis, through cardio-vascular, hematopoietic, and brain development, to signalingin the nervous and immune system.