Study of ligand-based virtual screening tools in computer-aided drug design

Virtual screening is a central technique in drug discovery today. Millions of molecules can be tested in silico with the aim to only select the most promising and test them experimentally. The topic of this thesis is ligand-based virtual screening tools which take existing active molecules as starting point for finding new drug candidates. One goal of this thesis was to build a model that gives the probability that two molecules are biologically similar as function of one or more chemical similarity scores. Another important goal was to evaluate how well different ligand-based virtual screening tools are able to distinguish active molecules from inactives. One more criterion set for the virtual screening tools was their applicability in scaffold-hopping, i.e. finding new active chemotypes. In the first part of the work, a link was defined between the abstract chemical similarity score given by a screening tool and the probability that the two molecules are biologically similar. These results help to decide objectively which virtual screening hits to test experimentally. The work also resulted in a new type of data fusion method when using two or more tools. In the second part, five ligand-based virtual screening tools were evaluated and their performance was found to be generally poor. Three reasons for this were proposed: false negatives in the benchmark sets, active molecules that do not share the binding mode, and activity cliffs. In the third part of the study, a novel visualization and quantification method is presented for evaluation of the scaffold-hopping ability of virtual screening tools.

A novel Bayesian approach to drug design with simple and complex enzymes: Use with lens aldehyde dehydrogenase and the glyoxalase pathway

Purpose: Acquiring details of kinetic parameters of enzymes is crucial to biochemical understanding, drug development, and clinical diagnosis in ocular diseases. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. Methods: We have developed Bayesian utility functions to minimise kinetic parameter variance involving differentiation of model expressions and matrix inversion. These have been applied to the simple kinetics of the enzymes in the glyoxalase pathway (of importance in posttranslational modification of proteins in cataract), and the complex kinetics of lens aldehyde dehydrogenase (also of relevance to cataract). Results: Our successful application of Bayesian statistics has allowed us to identify a set of rules for designing optimum kinetic experiments iteratively. Most importantly, the distribution of points in the range is critical; it is not simply a matter of even or multiple increases. At least 60 % must be below the KM (or plural if more than one dissociation constant) and 40% above. This choice halves the variance found using a simple even spread across the range.With both the glyoxalase system and lens aldehyde dehydrogenase we have significantly improved the variance of kinetic parameter estimation while reducing the number and costs of experiments. Conclusions: We have developed an optimal and iterative method for selecting features of design such as substrate range, number of measurements and choice of intermediate points. Our novel approach minimises parameter error and costs, and maximises experimental efficiency. It is applicable to many areas of ocular drug design, including receptor-ligand binding and immunoglobulin binding, and should be an important tool in ocular drug discovery.

Advances in drug design based on the amino acid approach: Taurine analogues for the treatment of CNS diseases

Amino acids are well known to be an important class of compounds for the maintenance of body homeostasis and their deficit, even for the polar neuroactive aminoacids, can be controlled by supplementation. However, for the amino acid taurine (2-aminoethanesulfonic acid) this is not true. Due its special physicochemical properties, taurine is unable to cross the blood-brain barrier. In addition of injured taurine transport systems under pathological conditions, CNS supplementation of taurine is almost null. Taurine is a potent antioxidant and anti-inflammatory semi-essential amino acid extensively involved in neurological activities, acting as neurotrophic factor, binding to GABA A/glycine receptors and blocking the excitotoxicity glutamate-induced pathway leading to be a neuroprotective effect and neuromodulation. Taurine deficits have been implicated in several CNS diseases, such as Alzheimer's, Parkinson's, epilepsy and in the damage of retinal neurons. This review describes the CNS physiological functions of taurine and the development of new derivatives based on its structure useful in CNS disease treatment.&; 2012 by the authors; licensee MDPI, Basel, Switzerland.

Anti-inflammatory drug design by molecular hybridization approach

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.

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.

Synthesis of Biologically Active Small Molecules: Different Approaches to Drug Design

In the past years, genome biology had disclosed an ever-growing kind of biological targets that emerged as ideal points for therapeutic intervention. Nevertheless, the number of new chemical entities (NCEs) translated into effective therapies employed in the clinic, still not observed. Innovative strategies in drug discovery combined with different approaches to drug design should be searched for bridge this gap. In this context organic synthetic chemistry had to provide for effective strategies to achieve biologically active small molecules to consider not only as potentially drug candidates, but also as chemical tools to dissect biological systems. In this scenario, during my PhD, inspired by the Biology-oriented Synthesis approach, a small library of hybrid molecules endowed with privileged scaffolds, able to block cell cycle and to induce apoptosis and cell differentiation, merged with natural-like cores were synthesized. A synthetic platform which joined a Domino Knoevenagel-Diels Alder reaction with a Suzuki coupling was performed in order to reach the hybrid compounds. These molecules can represent either antitumor lead candidates, or valuable chemical tools to study molecular pathways in cancer cells. The biological profile expressed by some of these derivatives showed a well defined antiproliferative activity on leukemia Bcr-Abl expressing K562 cell lines. A parallel project regarded the rational design and synthesis of minimally structurally hERG blockers with the purpose of enhancing the SAR studies of a previously synthesized collection. A Target-Oriented Synthesis approach was applied. Combining conventional and microwave heating, the desired final compounds were achieved in good yields and reaction rates. The preliminary biological results of the compounds, showed a potent blocking activity. The obtained small set of hERG blockers, was able to gain more insight the minimal structural requirements for hERG liability, which is mandatory to investigate in order to reduce the risk of potential side effects of new drug candidates.

Conformationally homogeneous heterocyclic pseudotetrapeptides as three-dimensional scaffolds for rational drug design: receptor-selective somatostatin analogues

A would-be amide: A 1,4-disubstituted 1,2,3-triazole was used as a surrogate for a trans amide bond to create a library of 16 diastereomeric pseudotetrapeptides as beta-turn mimetics. High-resolution structural analysis indicated that these scaffolds adopt distinct, rigid, conformationally homogeneous beta-turn-like structures (see example), some of which bind somatostatin receptor subtypes selectively, and some of which show broad-spectrum activity.

Mapping the active site in vasoactive intestinal peptide to a core of four amino acids: Neuroprotective drug design

The understanding of the molecular mechanisms leading to peptide action entails the identification of a core active site. The major 28-aa neuropeptide, vasoactive intestinal peptide (VIP), provides neuroprotection. A lipophilic derivative with a stearyl moiety at the N-terminal and norleucine residue replacing the Met-17 was 100-fold more potent than VIP in promoting neuronal survival, acting at femtomolar–picomolar concentration. To identify the active site in VIP, over 50 related fragments containing an N-terminal stearic acid attachment and an amidated C terminus were designed, synthesized, and tested for neuroprotective properties. Stearyl-Lys-Lys-Tyr-Leu-NH2 (derived from the C terminus of VIP and the related peptide, pituitary adenylate cyclase activating peptide) captured the neurotrophic effects offered by the entire 28-aa parent lipophilic derivative and protected against β-amyloid toxicity in vitro. Furthermore, the 4-aa lipophilic peptide recognized VIP-binding sites and enhanced choline acetyltransferase activity as well as cognitive functions in Alzheimer’s disease-related in vivo models. Biodistribution studies following intranasal administration of radiolabeled peptide demonstrated intact peptide in the brain 30 min after administration. Thus, lipophilic peptide fragments offer bioavailability and stability, providing lead compounds for drug design against neurodegenerative diseases.

Structure of the HIV-1 integrase catalytic domain complexed with an inhibitor: A platform for antiviral drug design

HIV integrase, the enzyme that inserts the viral DNA into the host chromosome, has no mammalian counterpart, making it an attractive target for antiviral drug design. As one of the three enzymes produced by HIV, it can be expected that inhibitors of this enzyme will complement the therapeutic use of HIV protease and reverse transcriptase inhibitors. We have determined the structure of a complex of the HIV-1 integrase core domain with a novel inhibitor, 5ClTEP, 1-(5-chloroindol-3-yl)-3-hydroxy-3-(2H-tetrazol-5-yl)-propenone, to 2.1-Å resolution. The inhibitor binds centrally in the active site of the integrase and makes a number of close contacts with the protein. Only minor changes in the protein accompany inhibitor binding. This inhibitor complex will provide a platform for structure-based design of an additional class of inhibitors for antiviral therapy.

Structural plasticity of the cyclic-cystine-knot framework: implications for biological activity and drug design

The cyclotide family of plant proteins is of interest because of their unique topology, which combines a head-to-tail cyclic backbone with an embedded cystine knot, and because their-remarkable chemical and biological properties make them ideal candidates as grafting templates for biologically active peptide epitopes. The present Study describes the first steps towards exploiting the cyclotide framework by synthesizing and structurally characterizing two grafted analogues of the cyclotide kalata B1. The modified peptides have polar or charged residues substituted for residues that form part of a surface-exposed hydrophobic patch that plays a significant role in the folding and biological activity of kalata B1. Both analogues retain the native cyclotide fold, but lack the undesired haemolytic activity of their parent molecule, kalata B1. This finding confirms the tolerance of the cyclotide framework to residue Substitutions and opens up possibilities for the Substitution of biologically active peptide epitopes into the framework.

The cyclotides and related macrocyclic peptides as scaffolds in drug design

The applicability of linear peptides as drugs is potentially limited by their susceptibility to proteolytic cleavage and poor bioavailability. Cyclotides are macrocyclic cystine-knotted mini-proteins that have a broad range of bioactivities and are exceptionally stable, being resistant to chemical, thermal and enzymatic degradation. The general limitations of peptides as drugs can potentially be overcome by using the cyclotide framework as a scaffold onto which new activities may be engineered. The potential use of cyclotides and related peptide scaffolds for drug design is evaluated herein, with reference to increasing knowledge of the structures and sequence diversity of natural cyclotides and the emergence of new approaches in protein engineering.