Application of Inductive Logic Programming to Structure-Based Drug Design

Enot, D. and King, R. D. (2003) Application of Inductive Logic Programming to Structure-Based Drug Design. 7th European Conference on Principles and Practice of Knowledge Discovery in Databases (PKDD '03). Springer LNAI 2838 p156-167

Structure based drug design with inductive logic programming

David P. Enot and Ross D. King (2003). Structure based drug design with inductive logic programming. The ACS National Meeting Spring 2003, New Orleans

Ligand- and Structure-Based Drug Design Strategies and PPAR delta/alpha Selectivity

Peroxisome-proliferator-activated receptors are a class of nuclear receptors with three subtypes: a, ? and d. Their main function is regulating gene transcription related to lipid and carbohydrate metabolism. Currently, there are no peroxisome-proliferator-activated receptors d drugs being marketed. In this work, we studied a data set of 70 compounds with a and d activity. Three partial least square models were created, and molecular docking studies were performed to understand the main reasons for peroxisome-proliferator-activated receptors d selectivity. The obtained results showed that some molecular descriptors (log P, hydration energy, steric and polar properties) are related to the main interactions that can direct ligands to a particular peroxisome-proliferator-activated receptors subtype.

Structure-based drug design studies on a series of aldolase inhibitors

Human African trypanosomiasis, also known as sleeping sickness, is a major cause of death in Africa, and for which there are no safe and effective treatments available. The enzyme aldolase from Trypanosoma brucei is an attractive, validated target for drug development. A series of alkyl‑glycolamido and alkyl-monoglycolate derivatives was studied employing a combination of drug design approaches. Three-dimensional quantitative structure-activity relationships (3D QSAR) models were generated using the comparative molecular field analysis (CoMFA). Significant results were obtained for the best QSAR model (r2 = 0.95, non-cross-validated correlation coefficient, and q2 = 0.80, cross-validated correlation coefficient), indicating its predictive ability for untested compounds. The model was then used to predict values of the dependent variables (pKi) of an external test set,the predicted values were in good agreement with the experimental results. The integration of 3D QSAR, molecular docking and molecular dynamics simulations provided further insight into the structural basis for selective inhibition of the target enzyme.

Structure- and ligand-based drug design approaches for neglected tropical diseases

Drug discovery has moved toward more rational strategies based on our increasing understanding of the fundamental principles of protein-ligand interactions. Structure( SBDD) and ligand-based drug design (LBDD) approaches bring together the most powerful concepts in modern chemistry and biology, linking medicinal chemistry with structural biology. The definition and assessment of both chemical and biological space have revitalized the importance of exploring the intrinsic complementary nature of experimental and computational methods in drug design. Major challenges in this field include the identification of promising hits and the development of high-quality leads for further development into clinical candidates. It becomes particularly important in the case of neglected tropical diseases (NTDs) that affect disproportionately poor people living in rural and remote regions worldwide, and for which there is an insufficient number of new chemical entities being evaluated owing to the lack of innovation and R&D investment by the pharmaceutical industry. This perspective paper outlines the utility and applications of SBDD and LBDD approaches for the identification and design of new small-molecule agents for NTDs.

Gut hormone GPCRs: structure, function, drug discovery

Crystallization and determination of the high resolution three-dimensional structure of the β2-adrenergic receptor in 2007 was followed by structure elucidation of a number of other receptors, including those for neurotensin and glucagon. These major advances foster the understanding of structure-activity relationship of these receptors and structure-based rational design of new ligands having more predictable activity. At present, structure determination of gut hormone receptors in complex with their ligands (natural, synthetic) and interacting signalling proteins, for example, G-proteins, arrestins, represents a challenge which promises to revolutionize gut hormone endocrinonology.

Crystal structure of human purine nucleoside phosphorylase complexed with acyclovir

In human, purine nucleoside phosphorylase (HsPNP) is responsible for degradation of deoxyguanosine and genetic deficiency of this enzyme leads to profound T-cell mediated immunosuppression. PNP is therefore a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. This work reports the first crystallographic Study of human PNP complexed with acyclovir (HsPNP:Acy). Acyclovir is a potent clinically useful inhibitor of replicant herpes simplex virus that also inhibits human PNP but with a relatively lower inhibitory activity (K-i=90muM). Analysis of the structural differences among the HsPNP:Acy complex, PNP apoenzyme, and HsPNP:Immucillin-H provides explanation for inhibitor binding, refines the purine-binding site, and can be used for future inhibitor design. (C) 2003 Published by Elsevier B.V.

Crystal structure of human PNP complexed with guanine

Purine nucleoside phosphorylase (PNP) catalyzes the phosphorolysis of the N-ribosidic bonds of purine nucleosides and deoxynucleosides. PNP is a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. More recently, the 3-D structure of human PNP has been refined to 2.3 Angstrom resolution, which allowed a redefinition of the residues involved in the substrate-binding sites and provided a more reliable model for structure-based design of inhibitors. This work reports crystallographic study of the complex of Human PNP:guanine (HsPNP:Gua) solved at 2.7 Angstrom resolution using synchrotron radiation. Analysis of the structural differences among the HsPNP:Gua complex, PNP apoenzyme, and HsPNP:immucillin-H provides explanation for inhibitor binding, refines the purine-binding site, and can be used for future inhibitor design. (C) 2003 Elsevier B.V. All rights reserved.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Molecular modelling assisted design and synthesis of cyclothialidine derivatives as DNA gyrase inhibitors

Since cyclothialidine was discovered as the most active DNA gyrase inhibitor in 1994, enormous efforts have been devoted to make it into a commercial medicine by a number of pharmaceutical companies and research groups worldwide. However, no serious breakthrough has been made up to now. An essential problem involved with cyclothialidine is that though it demonstrated the potent inhibition of DNA gyrase, it showed little activity against bacteria. This probably is attributable to its inability to penetrate bacterial cell walls and membranes. We applied the TSAR programme to generate a QSAR equation to the gram-negative organisms. In that equation, LogP is profoundly indicated as the key factor influencing the cyclothialidine activity against bacteria. However, the synthesized new analogues have failed to prove that. In the structure based drug design stage, we designed a group of open chain cyclothialidine derivatives by applying the SPROUT programme and completed the syntheses. Improved activity is found in a few analogues and a 3D pharmacophore of the DNA gyrase B is proposed to lead to synthesis of the new derivatives for development of potent antibiotics.

Relieving the first bottleneck in the drug discovery pipeline:using array technologies to rationalize membrane protein production

The slow down in the drug discovery pipeline is, in part, owing to a lack of structural and functional information available for new drug targets. Membrane proteins, the targets of well over 50% of marketed pharmaceuticals, present a particular challenge. As they are not naturally abundant, they must be produced recombinantly for the structural biology that is a prerequisite to structure-based drug design. Unfortunately, however, obtaining high yields of functional, recombinant membrane proteins remains a major bottleneck in contemporary bioscience. While repeated rounds of trial-and-error optimization have not (and cannot) reveal mechanistic details of the biology of recombinant protein production, examination of the host response has provided new insights. To this end, we published an early transcriptome analysis that identified genes implicated in high-yielding yeast cell factories, which has enabled the engineering of improved production strains. These advances offer hope that the bottleneck of membrane protein production can be relieved rationally.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.