EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Mapeamento das interações proteína-ligante através de técnicas de RMN de ¹H utilizando detecção do ligante

NMR is a valuable screening tool for the binding of ligands to proteins providing structural information on both protein and ligands and is thus largely applied to drug-discovery. Among the recent NMR techniques to probe weak binding protein-ligand complexes we have critically evaluated the advantages and disadvantages of STD (Saturation Transfer Difference), WaterLOGSY (Water Ligand Observation with Gradient Spectroscopy), NOE pumping and DOSY-NOESY (Diffusion-Ordered NOESY) using a mixture of BSA (bovine serum albumin) plus salicylic acid, caffeine, citric acid, adipic acid and D-glucose.

Síntese, avaliação biológica e modelagem molecular de arilfuranos como inibidores da enzima tripanotiona redutase

Trypanosoma cruzi is a protozoan parasite that causes a severe disease (Chagas'disease) in Central and South America. The currently available chemotherapeutic agents against this disease are still inadequate. The enzyme trypanothione reductase (TR) is considered a validated molecular target for the development of new drugs against this parasite. In this regard, a series of arylfurans based on 2,5-bis-(4-acetamidophenyl)furan was synthesized and tested for their in vitro inhibitory activity against TR. Molecular modeling studies of putative enzyme-inhibitor complexes revealed a possible mechanism of interaction. From synthesized compounds, a benzylaminofuran derivative was found to be more active than the lead compound.

INTERAÇÕES FÁRMACO-RECEPTOR: APLICAÇÕES DE TÉCNICAS COMPUTACIONAIS EM AULA PRÁTICA SOBRE A EVOLUÇÃO DOS INIBIDORES DA ENZIMA CONVERSORA DE ANGIOTENSINA

Teaching classes and events regarding the molecular aspects of drug-receptor interactions is not an easy task. The ligand stereochemistry and the spatial arrangement of the macromolecular targets highly increase the complexity of the process. In this context, the use of alternative and more playful approaches could allow students to gain a more thorough understanding of this important topic in medicinal chemistry. Herein, we describe a practical teaching approach that uses computational strategies as a tool for drug-receptor interaction studies performed for angiotencsin converting enzyme inhibitors (ACEi). Firstly, the students learn how to find the crystallographic structure (enzyme-ligand complex). Then, they proceed to the treatment of crude crystallographic data. Thereafter, they learn how to analyze the positioning of the drug on the active site of the enzyme, looking for regions related to the molecular recognition. At the end of the study, students can summarize the molecular requirements for the interaction and the structure-activity relationships of the studied drugs.

Response to zinc deficiency of two rice lines with contrasting tolerance is determined by root growth maintenance and organic acid exudation rates, and not by zinc-transporter activity

Zinc (Zn)-deficient soils constrain rice (Oryza sativa) production and cause Zn malnutrition. The identification of Zn-deficiency-tolerant rice lines indicates that breeding might overcome these constraints. Here, we seek to identify processes underlying Zn-deficiency tolerance in rice at the physiological and transcriptional levels. A Zn-deficiency-tolerant line RIL46 acquires Zn more efficiently and produces more biomass than its nontolerant maternal line (IR74) at low Zn(ext) under field conditions. We tested if this was the result of increased expression of Zn(2+) transporters; increased root exudation of deoxymugineic acid (DMA) or low-molecular-weight organic acids (LMWOAs); and/or increased root production. Experiments were performed in field and controlled environment conditions. There was little genotypic variation in transcript abundance of Zn-responsive root Zn(2+)-transporters between the RIL46 and IR74. However, root exudation of DMA and LMWOA was greater in RIL46, coinciding with increased root expression of putative ligand-efflux genes. Adventitious root production was maintained in RIL46 at low Zn(ext), correlating with altered expression of root-specific auxin-responsive genes. Zinc-deficiency tolerance in RIL46 is most likely the result of maintenance of root growth, increased efflux of Zn ligands, and increased uptake of Zn-ligand complexes at low Zn(ext); these traits are potential breeding targets.

Perspectives for Novel Mixed Diruthenium-Organic Drugs as Metallopharmaceuticals in Cancer Therapy

Ruthenium compounds have been actively studied as metallodrugs for cancer therapy. Representatives of ruthenium-based antitumor drugs are the classes of ruthenium(III)-chlorido-(N-ligand)complexes, including the drugs namely NAMI-A and KP1019 in clinical trials, and ruthenium(II)-arene organometallics, with some compounds currently undergoing advanced preclinical testing. An alternative approach for tumor-inhibiting metallodrugs is the coordination of metal ions to organic pharmaceuticals. The combination of antitumor-active ruthenium ion with biologically-active pro-ligands in single compounds can result in the enhancement of activity, for example through synergistic effects. In the present article, some developments in the ruthenium-based antitumor drugs field are briefly highlighted and recent studies on mixed diruthenium-organic drugs as metallopharmaceuticals in cancer therapy are described. Novel organic pharmaceuticals-containing diruthenium(II, III)complexes have shown promising antitumor activity for C6 rat glioma - a model for glioblastoma multiforme (GBA).

Spectrophotometric study of binary and ternary systems involving metal ions and benzylidenepyruvates: Equilibria in aqueous solutions

The protonation of 4-dimethylaminobenzylidenepyruvate (DMBP) and 2-chloro-4-dimethylaminobenzylidenepyruvate (2-CI-DMBP) and their complex formation with Mn(II), Co(II), Ni(II), Cu(II), Zn(II), Pb(II), Cd(II) and Al(III) have been studied by potentiometric and spectrophotometric methods at 25 °C and ionic strength 0.500 M, held with sodium perchlorate. The stability order found for 1 :1 complexes of both ligands is Al(III) > Cu(II) > Pb(II) > Ni(II) > Zn(II) > Co(II) > Cd(II) > Mn(II). The stability changes move in the same direction as the pKa of the ligands. The results are compared with literature values reported for metal ion pyruvate systems. Thermodynamic stabilities of ternary complexes formed in Cu(II)-B-L- systems, where B = 2,2′-bipyridyl (bipy), ethylenediamine or glycinate and L = DMBP or 2-CI-DMBP, were also determined. The Cu(bipy)L+ species are more stable than would be expected on purely statistical grounds. The importance of the :t system associated with bipy on the enhanced stability of its mixed ligand complexes is stressed. Analytical applications of the investigated ligands are outlined.

Structure and catalytic mechanism of glucosamine 6-phosphate deaminase from Escherichia coli at 2.1 Å resolution

Background: Glucosamine 6-phosphate deaminase from Escherichia coli is an allosteric hexameric enzyme which catalyzes the reversible conversion of D-glucosamine 6-phosphate into D-fructose 6-phosphate and ammonium ion and is activated by N-acetyl-D-glucosamine 6-phosphate. Mechanistically, it belongs to the group of aldose-ketose isomerases, but its reaction also accomplishes a simultaneous amination/deamination. The determination of the structure of this protein provides fundamental knowledge for understanding its mode of action and the nature of allosteric conformational changes that regulate its function. Results: The crystal structure of glucosamine 6-phosphate deaminase with bound phosphate ions is presented at 2.1 Å resolution together with the refined structures of the enzyme in complexes with its allosteric activator and with a competitive inhibitor. The protein fold can be described as a modified NAD-binding domain. Conclusions: From the similarities between the three presented structures, it is concluded that these represent the enzymatically active R state conformer. A mechanism for the deaminase reaction is proposed. It comprises steps to open the pyranose ring of the substrate and a sequence of general base-catalyzed reactions to bring about isomerization and deamination, with Asp72 playing a key role as a proton exchanger.

Enzymatic transesterification of soybean oil with ethanol using lipases immobilized on highly crystalline PVA microspheres

Polyvinyl alcohol (PVA) microspheres with different degree of crystallinity were used as solid supports for Rhizomucor miehei lipase immobilization, and the enzyme-PVA complexes were used as biocatalysts for the transesterification of soybean oil to fatty acid ethyl esters (FAEE). The amounts of immobilized enzyme on the polymeric supports were similar for both the amorphous microspheres (PVA4) and the high crystalline microspheres (PVA25). However, the enzymatic activity of the immobilized enzymes was depended on the crystallinity degree of the PVA microspheres: enzymes immobilized on the PVA4 microspheres have shown low enzymatic activity (6.13 U mg-1), in comparison with enzymes immobilized on the high crystalline PVA25 microspheres (149.15 U mg-1). A synergistic effect was observed for the enzyme-PVA25 complex during the transesterification reaction of soybean oil to FAEE: transesterification reactions with free enzyme with the equivalent amount of enzyme that were immobilized onto the PVA25 microspheres (5.4 U) have yielded only 20% of FAEE, reactions with the pure highly crystalline microsphere PVA25 have not yielded FAEE, however reactions with the enzyme-PVA25 complexes have yielded 66.3% of FAEE. This synergistic effect of an immobilized enzyme on a polymeric support has not been observed before for transesterification reaction of triacylglycerides into FAEE. Based on ATR-FTIR, 23Na- and 13C-NMR-MAS spectroscopic data and the interaction of the polymeric network intermolecular hydrogen bonds with the lipases residual amino acids a possible explanation for this synergistic effect is provided. © 2013 Elsevier Ltd. All rights reserved.

Complexos de inclusão com amilose: obtenção, caracterização e avaliação como sistemas de liberação controlada de fármacos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Untersuchung der DNA-bindenden Eigenschaften neuer Nukleobasen-gekoppelter Pyrrolcarboxamide, Combilexine und Hetaren[a]carbazol-Derivate: Validierung und Etablierung von in-vitro- und in-silico-Methoden

Da maligne Neoplasien durch Mutationen in Proto-Onko- und/oder Tumorsuppressorgenen ausgelöst werden, stellt die DNA eines der wichtigsten Targets für die Entwicklung neuer Zytostatika dar. Auch bei den im Arbeitskreis Pindur designten und synthetisier-ten Verbindungen der Nukleobasen-gekoppelten Pyrrolcarboxamid-, der Hetaren[a]carbazol- und der Combilexin-Reihe handelt es sich um DNA-Liganden mit potentiell antitumoraktiven Eigenschaf-ten. Die einen dualen Bindemodus aufweisenden Combilexine bestehen aus einem Interkalator (u. a. Naphthalimid, Acridon), der über einen Linker variabler Kettenlänge mit einer rinnenbin-denden, von Netropsin abgeleiteten Bispyrrol-, oder einer bioisosteren Imidazol-, Thiazol- oder Thiophen-pyrrolcarboxamid-struktur verknüpft ist. Das N-terminale Ende der Combilexine wird von einer N,N-Dimethylaminopropyl- oder -ethyl-Seitenkette gebildet. Die DNA-Affinitäten der Liganden wurden mittels Tm-Wert-Messung-en bestimmt. Diese Denaturierungsexperimente wurden sowohl mit poly(dAdT)2- als auch mit Thymus-DNA (~42% GC-Anteil) durchge-führt, um Aussagen zur Stärke und zur Sequenzselektivität der DNA-Bindung machen zu können. Des Weiteren wurden die Bindekon-stanten einiger ausgewählter Vertreter mit Hilfe des Ethidium-bromid-Verdrängungsassays ermittelt; einige Testverbindungen wurden zudem auf potentiell vorhandene, TOPO I-inhibierende Eigenschaften untersucht. Diese biochemischen und biophysika-lischen Tests wurden durch Molecular Modelling-Studien ergänzt, die die Berechnung von molekularen Eigenschaften, die Durch-führung von Konformerenanalysen und die Simulation von DNA-Ligand-Komplexen (Docking) umfassten. Durch Korrelation der in vitro-Befunde mit den in silico-Daten gelang es, vor allem für die Substanzklasse der Combilexine einige richtungweisende Struktur-Wirkungsbeziehungen aufzustellen. So konnte gezeigt werden, dass die Einführung eines Imidazol-Rings in die rinnen-bindende Hetaren-pyrrolcarboxamid-Struktur der Combilexine aufgrund der H-Brücken-Akzeptor-Funktion des sp2-hybridisierten N-Atoms eine Verschiebung der Sequenzselektivität der DNA-Bindung von AT- zu GC-reichen Arealen der DNA bedingt. Zudem erwies sich ein C3-Linker für die Verknüpfung des Naphthalimids mit dem rinnenbindenden Strukturelement als am besten geeignet, während bei den Acridon-Derivaten die Verbindungen mit einem N-terminalen Buttersäure-Linker die höchste DNA-Affinität aufwiesen. Dies ist sehr wahrscheinlich auf die im Vergleich zum Naphthalimid-Molekül geringere y-Achsen-Ausdehnung (bzgl. eines x/y-Koordinatensystems) des Acridons zurückzuführen. Die ermittelten Struktur-Wirkungsbeziehungen können dazu herangezogen werden, das rationale Design neuer DNA-Liganden mit potentiell stärkerer DNA-Bindung zu optimieren.

Estrogen receptor (ER) modulators each induce distinct conformational changes in ER α and ER β

Estrogen receptor (ER) modulators produce distinct tissue-specific biological effects, but within the confines of the established models of ER action it is difficult to understand why. Previous studies have suggested that there might be a relationship between ER structure and activity. Different ER modulators may induce conformational changes in the receptor that result in a specific biological activity. To investigate the possibility of modulator-specific conformational changes, we have applied affinity selection of peptides to identify binding surfaces that are exposed on the apo-ERs α and β and on each receptor complexed with estradiol or 4-OH tamoxifen. These peptides are sensitive probes of receptor conformation. We show here that ER ligands, known to produce distinct biological effects, induce distinct conformational changes in the receptors, providing a strong correlation between ER conformation and biological activity. Furthermore, the ability of some of the peptides to discriminate between different ER α and ER β ligand complexes suggests that the biological effects of ER agonists and antagonists acting through these receptors are likely to be different.

Identification of 5′-deoxyribose phosphate lyase activity in human DNA polymerase γ and its role in mitochondrial base excision repair in vitro

Mitochondria have been proposed to possess base excision repair processes to correct oxidative damage to the mitochondrial genome. As the only DNA polymerase (pol) present in mitochondria, pol γ is necessarily implicated in such processes. Therefore, we tested the ability of the catalytic subunit of human pol γ to participate in uracil-provoked base excision repair reconstituted in vitro with purified components. Subsequent to actions of uracil-DNA glycosylase and apurinic/apyrimidinic endonuclease, human pol γ was able to fill a single nucleotide gap in the presence of a 5′ terminal deoxyribose phosphate (dRP) flap. We report here that the catalytic subunit of human pol γ catalyzes release of the dRP residue from incised apurinic/apyrimidinic sites to produce a substrate for DNA ligase. The heat sensitivity of this activity suggests the dRP lyase function requires a three-dimensional protein structure. The dRP lyase activity does not require divalent metal ions, and the ability to trap covalent enzyme-DNA complexes with NaBH4 strongly implicates a Schiff base intermediate in a β-elimination reaction mechanism.

Conversion of agonist site to metal-ion chelator site in the β2-adrenergic receptor

Previously metal-ion sites have been used as structural and functional probes in seven transmembrane receptors (7TM), but as yet all the engineered sites have been inactivating. Based on presumed agonist interaction points in transmembrane III (TM-III) and -VII of the β2-adrenergic receptor, in this paper we construct an activating metal-ion site between the amine-binding Asp-113 in TM-III—or a His residue introduced at this position—and a Cys residue substituted for Asn-312 in TM-VII. No increase in constitutive activity was observed in the mutant receptors. Signal transduction was activated in the mutant receptors not by normal catecholamine ligands but instead either by free zinc ions or by zinc or copper ions in complex with small hydrophobic metal-ion chelators. Chelation of the metal ions by small hydrophobic chelators such as phenanthroline or bipyridine protected the cells from the toxic effect of, for example Cu2+, and in several cases increased the affinity of the ions for the agonistic site. Wash-out experiments and structure–activity analysis indicated, that the high-affinity chelators and the metal ions bind and activate the mutant receptor as metal ion guided ligand complexes. Because of the well-understood binding geometry of the small metal ions, an important distance constraint has here been imposed between TM-III and -VII in the active, signaling conformation of 7TM receptors. It is suggested that atoxic metal-ion chelator complexes could possibly in the future be used as generic, pharmacologic tools to switch 7TM receptors with engineered metal-ion sites on or off at will.

Molecular determinants of tissue selectivity in estrogen receptor modulators

Interaction of the estrogen receptor/ligand complex with a DNA estrogen response element is known to regulate gene transcription. In turn, specific conformations of the receptor-ligand complex have been postulated to influence unique subsets of estrogen-responsive genes resulting in differential modulation and, ultimately, tissue-selective outcomes. The estrogen receptor ligands raloxifene and tamoxifen have demonstrated such tissue-specific estrogen agonist/antagonist effects. Both agents antagonize the effects of estrogen on mammary tissue while mimicking the actions of estrogen on bone. However, tamoxifen induces significant stimulation of uterine tissue whereas raloxifene does not. We postulate that structural differences between raloxifene and tamoxifen may influence the conformations of their respective receptor/ligand complexes, thereby affecting which estrogen-responsive genes are modulated in various tissues. These structural differences are 4-fold: (A) the presence of phenolic hydroxyls, (B) different substituents on the basic amine, (C) incorporation of the stilbene moiety into a cyclic benzothiophene framework, and (D) the imposition of a carbonyl “hinge” between the basic amine-containing side chain and the olefin. A series of raloxifene analogs that separately exemplify each of these differences have been prepared and evaluated in a series of in vitro and in vivo assays. This strategy has resulted in the development of a pharmacophore model that attributes the differences in effects on the uterus between raloxifene and tamoxifen to a low-energy conformational preference imparting an orthogonal orientation of the basic side chain with respect to the stilbene plane. This three-dimensional array is dictated by a single carbon atom in the hinge region of raloxifene. These data indicate that differences in tissue selective actions among benzothiophene and triarylethylene estrogen receptor modulators can be ascribed to discrete ligand conformations.