Carlos A. Montanari

Inflation of correlation in the pursuit of drug-likeness

Drug-likeness is a frequently invoked, although not always precisely defined, concept in drug discovery. Opinions on drug-likeness are to a large extent shaped by the relationships that are observed between surrogate measures of drug-likeness (e.g. aqueous solubility; permeability; pharmacological promiscuity) and fundamental physicochemical properties (e.g. lipophilicity; molecular size). This article draws on examples from the literature to highlight approaches to data analysis that exaggerate trends in data and the term correlation inflation is introduced in the context of drug discovery. Averaging groups of data points prior to analysis is a common cause of correlation inflation and results from analysis of binned continuous data should always be treated with caution.

Discovery of novel Trypanosoma cruzi glyceraldehyde-3-phosphate dehydrogenase inhibitors

Based on its essential role in the life cycle of Trypanosoma cruzi, the glycolytic enzyme glyceraldehyde 3-phosphate dehydrogenase (GAPDH) has been considered a promising target for the development of novel chemotherapeutic agents for the treatment of Chagas disease. In the course of our research program to discover novel inhibitors of this trypanosomatid enzyme, we have explored a combination of structure and ligand-based virtual screening techniques as a complementary approach to a biochemical screening of natural products using a standard biochemical assay. Seven natural products, including anacardic acids, flavonoid derivatives, and one glucosylxanthone were identified as novel inhibitors of T. cruzi GAPDH. Promiscuous inhibition induced by nonspecific aggregation has been discarded as specific inhibition was not reversed or affected in all cases in the presence of Triton X-100, demonstrating the ability of the assay to find authentic inhibitors of the enzyme. The structural diversity of this series of promising natural products is of special interest in drug design, and should therefore be useful in future medicinal chemistry efforts aimed at the development of new GAPDH inhibitors having increased potency.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: implications to tetramer stability and ligand-binding.

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

Novel insights for dihydroorotate dehydrogenase class 1A inhibitors discovery.

The enzyme dihydroorotate dehydrogenase (DHODH) has been suggested as a promising target for the design of trypanocidal agents. We report here the discovery of novel inhibitors of Trypanosoma cruzi DHODH identified by a combination of virtual screening and ITC methods. Monitoring of the enzymatic reaction in the presence of selected ligands together with structural information obtained from X-ray crystallography analysis have allowed the identification and validation of a novel site of interaction (S2 site). This has provided important structural insights for the rational design of T. cruzi and Leishmania major DHODH inhibitors. The most potent compound (1) in the investigated series inhibits TcDHODH enzyme with Kiapp value of 19.28 μM and possesses a ligand efficiency of 0.54 kcal mol(-1) per non-H atom. The compounds described in this work are promising hits for further development.

Ligand efficiency metrics considered harmful

Ligand efficiency metrics are used in drug discovery to normalize biological activity or affinity with respect to physicochemical properties such as lipophilicity and molecular size. This Perspective provides an overview of ligand efficiency metrics and summarizes thermodynamics of protein–ligand binding. Different classes of ligand efficiency metric are critically examined and the study concludes with suggestions for alternative ways to account for physicochemical properties when prioritizing and optimizing leads.

Non-peptidic Cruzain Inhibitors with Trypanocidal Activity Discovered by Virtual Screening and In Vitro Assay

A multi-step cascade strategy using integrated ligand- and target-based virtual screening methods was developed to select a small number of compounds from the ZINC database to be evaluated for trypanocidal activity. Winnowing the database to 23 selected compounds, 12 non-covalent binding cruzain inhibitors with affinity values (K i) in the low micromolar range (3–60 µM) acting through a competitive inhibition mechanism were identified. This mechanism has been confirmed by determining the binding mode of the cruzain inhibitor Nequimed176 through X-ray crystallographic studies. Cruzain, a validated therapeutic target for new chemotherapy for Chagas disease, also shares high similarity with the mammalian homolog cathepsin L. Because increased activity of cathepsin L is related to invasive properties and has been linked to metastatic cancer cells, cruzain inhibitors from the same library were assayed against it. Affinity values were in a similar range (4–80 µM), yielding poor selectivity towards cruzain but raising the possibility of investigating such inhibitors for their effect on cell proliferation. In order to select the most promising enzyme inhibitors retaining trypanocidal activity for structure-activity relationship (SAR) studies, the most potent cruzain inhibitors were assayed against T. cruzi-infected cells. Two compounds were found to have trypanocidal activity. Using compound Nequimed42 as precursor, an SAR was established in which the 2-acetamidothiophene-3-carboxamide group was identified as essential for enzyme and parasite inhibition activities. The IC50 value for compound Nequimed42 acting against the trypomastigote form of the Tulahuen lacZ strain was found to be 10.6±0.1 µM, tenfold lower than that obtained for benznidazole, which was taken as positive control. In addition, by employing the strategy of molecular simplification, a smaller compound derived from Nequimed42 with a ligand efficiency (LE) of 0.33 kcal mol−1 atom−1 (compound Nequimed176) is highlighted as a novel non-peptidic, non-covalent cruzain inhibitor as a trypanocidal agent candidate for optimization.

ClogPalk: a method for predicting alkane/water partition coefficient

Alkane/water partition coefficients (Palk) are less familiar to the molecular design community than their 1-octanol/water equivalents and access to both data and prediction tools is much more limited. A method for predicting alkane/water partition coefficient from molecular structure is introduced. The basis for the ClogPalk model is the strong (R2xa0=xa00.987) relationship between alkane/water partition coefficient and molecular surface area (MSA) that was observed for saturated hydrocarbons. The model treats a molecule as a perturbation of a saturated hydrocarbon molecule with the same MSA and uses increments defined for functional groups to quantify the extent to which logPalk is perturbed by the introduction each functional group. Interactions between functional groups, such as intramolecular hydrogen bonds are also parameterized within a perturbation framework. The functional groups and interactions between them are specified substructurally in a transparent and reproducible manner using SMARTS notation. The ClogPalk model was parameterized using data measured for structurally prototypical compounds that dominate the literature on alkane/water partition coefficients and then validated using an external test set of 100 alkane/water logP measurements, the majority of which were for drugs.

Multimilligram enantioresolution of sulfoxide proton pump inhibitors by liquid chromatography on polysaccharide-based chiral stationary phase.

The enantiomers of sulfoxide proton pump inhibitors--omeprazole, lansoprazole, rabeprazole and Ro 18-5364--were enantiomerically separated by liquid chromatography at multimilligram scale on a polysaccharide-based chiral stationary phase using normal and polar organic conditions as mobile phase. The values of the recovery and production rate were significant for each enantiomer; better results were achieved using a solid-phase injection system. However, this system was applied just for the enantiomeric separation of omeprazole to demonstrate the applicability of this injection mode at milligram scale. The chiroptical characterization of the compounds was performed using a polarimeter and a circular dichroism detector. The higher enantiomeric purity obtained for the isolated enantiomers suggests that the methods here described should be considered as a simple and rapid way to obtain enantiomeric pure standards for analytical purpose.

Kinetic mechanism and catalysis of Trypanosoma cruzi dihydroorotate dehydrogenase enzyme evaluated by isothermal titration calorimetry.

Trypanosoma cruzi dihydroorotate dehydrogenase (TcDHODH) catalyzes the oxidation of l-dihydroorotate to orotate with concomitant reduction of fumarate to succinate in the de novo pyrimidine biosynthetic pathway. Based on the important need to characterize catalytic mechanism of TcDHODH, we have tailored a protocol to measure TcDHODH kinetic parameters based on isothermal titration calorimetry. Enzymatic assays lead to Michaelis-Menten curves that enable the Michaelis constant (K(M)) and maximum velocity (V(max)) for both of the TcDHODH substrates: dihydroorotate (K(M)=8.6+/-2.6 microM and V(max)=4.1+/-0.7 microMs(-1)) and fumarate (K(M)=120+/-9 microM and V(max)=6.71+/-0.15 microMs(-1)). TcDHODH activity was investigated using dimethyl sulfoxide (10%, v/v) and Triton X-100 (0.5%, v/v), which seem to facilitate the substrate binding process with a small decrease in K(M). Arrhenius plot analysis allowed the determination of thermodynamic parameters of activation for substrates and gave some insights into the enzyme mechanism. Activation entropy was the main contributor to the Gibbs free energy in the formation of the transition state. A factor that might contribute to the unfavorable entropy is the hindered access of substrates to the TcDHODH active site where a loop at its entrance regulates the open-close channel for substrate access.

Pharmacophore Modeling for Anti-Chagas Drug Design Using the Fragment Molecular Orbital Method

Background Chagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives. Methodology/Principal Findings Intermolecular interactions in the complexes of TcDHODH with orotate, oxonate, and 43 orotate derivatives were analyzed by FMO calculation at the MP2/6-31G level. The results indicated that the orotate moiety, which is the base fragment of these compounds, interacts with the Lys43, Asn67, and Asn194 residues of TcDHODH and the cofactor flavin mononucleotide (FMN), whereas functional groups introduced at the orotate 5-position strongly interact with the Lys214 residue. Conclusions/Significance FMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs.