Refine
Document Type
- Doctoral Thesis (6)
Has Fulltext
- yes (6)
Is part of the Bibliography
- no (6)
Keywords
- 5-LO (1)
- 5-Lipoxygenase (1)
- Apoptose (1)
- Apoptosis (1)
- B-Lymphocytes (1)
- B-Lymphozyten (1)
- Caspase-6 (1)
- Entzündung (1)
- In vitro (1)
- In vitro assay (1)
Institute
- Pharmazie (6) (remove)
Als zellulärer Sensor für Gallensäuren und als Regulator zahlreicher metabolischer und inflamma-torischer Gene stellt der nukleäre Farnesoid X Rezeptor (FXR) ein vielversprechendes neues Wirkstofftarget dar. Die Aktivierung von FXR mit natürlichen oder synthetischen Liganden führte in vitro und in vivo zu zahlreichen wünschenswerten Effekten wie gesteigerter Insulinfreisetzung, verringerter Insulinresistenz oder verbessertem Lipidprofil. Daneben stellt die Aktivierung von FXR ein Prinzip zur Behandlung von Lebererkrankungen wie nicht-alkoholischer Fettleber und primärer billiärer Zirrhose dar, das mit dem FXR-Agonisten Obeticholsäure bereits in klinischen Studien überprüft wird. Existierende synthetische FXR-Liganden sind Fettsäure- bzw. Gallensäuremimetika und imitieren die physiologischen FXR-Agonisten. Die meisten synthetischen FXR-Liganden sind jedoch aufgrund von Toxizität, geringer Selektivität oder schlechter Bioverfügbarkeit nicht zur wie-teren klinischen Entwicklung geeignet. Sie stellen außerdem vornehmlich vollagonistische FXR-Ligan-den dar, doch die klinischen Erfahrungen mit Liganden anderer nukleärer Rezeptoren wie den Peroxi-somen Proliferator-aktivierten Rezeptoren (PPAR) oder den Estrogenrezeptoren (ER) haben gezeigt, dass eine zu starke Aktivierung eines Ligand-aktivierten Transkriptionsfaktors Risiken erheblicher Nebenwirkungen bergen kann. Eine Möglichkeit, dieser Gefahr vorzubeugen, bietet die Entwicklung partialagonistischer FXR-Liganden, die den Rezeptor nur mit moderater Amplitude aktivieren.
In dieser Arbeit wurde ausgehend von der in einem virtuellen Screening identifizierten Leitstruktur 1, durch medizinisch chemische Optimierung und Studien zu den Struktur-Wirkungs-Beziehungen (SAR) ein potenter und selektiver FXR-Partialagonist entwickelt. Die drei Molekülteile der Leitstruktur 1 (azide Kopfgruppe, zentraler Anthranilamidkörper und Acylsubstituent) wurden einzeln hinsichtlich ihrer Potenz an FXR untersucht und optimiert. In der Untersuchung der SAR des Acylsubstituenten zeigten sich ein 2-Naphthoyl- und ein 4-tert-Butylbenzoylsubstituent der in 1 enthaltenen 4-Methylbenzoylgruppe überlegen. Unter Beibehaltung des 2-Naphthoylsubstituenten wurde hierauf durch selektive Methylierung bzw. Reduktion die Notwendigkeit beider Amidbindungen der Subs-tanzklasse für Aktivität an FXR nachgewiesen. Durch Erweiterung der aziden Kopfgruppe um einen zusätzlichen aromatischen Ring gelang eine weitere Potenzsteigerung, die sich durch eine Methyl-gruppe an 6-Position dieses neu eingeführten Ringes noch erhöhen ließe. Andere Substituenten am aromatischen Ring der Kopfgruppe führten dagegen an keiner Position zu einer Aktivitätsverbes-serung. Der Austausch der freien Carbonsäure durch metabolisch stabilere Bioisostere wie ein Methylketon oder ein Nitril stellte sich als ohne Aktivitätsverlust möglich heraus, wobei das Tetrazol als klassisches Carbonsäurebioisoster eine Ausnahme mit geringerer Potenz bildete. Die entscheiden-de Steigerung der Aktivität der Acylanthranilamide an FXR resultierte aus der Einführung eines zusätzlichen Substituenten in 4-Position des zentralen aromatischen Ringes, wobei eine Methoxy-gruppe zur größten Potenz führte. Das resultierende Anthranilamid 2 stellt einen hochpotenten FXR-Partialagonisten mit einem EC50-Wert von 8±3 nM in einem flFXR-Reportergenassay bei 18±1% Maximalaktivierung dar und ist der Leitstruktur somit um mehr als einen Faktor 1000 überlegen.
Die optimierte Verbindung 2 wurde aufgrund ihrer großen Potenz ausführlich in vitro pharma-kologisch charakterisiert. Dabei stellte sich die Substanz als metabolisch sehr stabil, moderat löslich in Wasser und gemessen an ihrer hohen Aktivität an FXR als wenig toxisch heraus. Darüber hinaus er-wies sich 2 als selektiv für FXR über den membranständigen G-Protein-gekoppelten Gallen-säurerezeptor TGR5 (Faktor >1000) sowie über die nukleären Rezeptoren PPARα (>1000), PPARγ (~375) und PPARδ (>1000). Bei der Quantifizierung seiner Effekte auf die FXR-Targetgene SHP, CYP7A1, BSEP, OSTα und IBABP durch qRT-PCR übte 2 im Bereich 0,1 µM bis 10 µM einen konzen-trationsunabhängigen partialagonistischen Effekt von etwa 40% des Effektes des physiologischen FXR-Agonisten Chenodeoxycholsäure (CDCA) aus. Mit der Verbindung 2 wurde somit ein hochpotenter, selektiver und metabolisch stabiler FXR-Partialagonist entwickelt und charakterisiert, der sich für künftige in vitro und in vivo Studien zu partieller FXR-Aktivierung empfehlen kann.
In silico Methoden spielen in der Wirkstoffentwicklung eine immer größere Bedeutung. Sie können eine Größe Hilfe in der Analyse des Targets oder beim Screening von neuen Liganden sein. Ihre Stärken liegen vorallem in der Zeit- und Kostenreduzierung während einer
Wirkstoffentwicklung.
Ziel der Arbeit war die Entwicklung neuer COX-2 Liganden mit Hilfe von in silico Methoden. Weil von der mCOX-2 keine Kristallstruktur in der PDB publiziert war, begann die Arbeit mit der Modellierung der mCOX-2. Dafür wurde aus der Sequenz der hCOX-2 aus UniProt mit der ID P35354 mit Hilfe der Kristallstruktur 3LN1 ein Homologie Modell entwickelt und im Anschluss über eine Validierungsmethode, den Ramachandran Plot, analysiert. Der Ramachandran Plot zeigte, dass 93.7% der Aminosäuren in favorisierten Regionen, 6.1% in
erlaubten Regionen, 0.2% in geduldeten Regionen und 0.0% in unerlaubten Regionen lagen. Mit diesem Modell wurde eine MD-Simulation durchgeführt, um die Energie des Modells zu
minimieren.
Die neuen Verbindungen wurden über drei verschiedene Ansätze designt. Im ersten Ansatz wurde die Software DOGS verwendet. Dabei handelt es sich um ein de novo Design Programm, welches nicht nur neue Verbindungen entwickelt, sondern auch deren Syntheseweg vorschlägt. Die vorgeschlagenen Verbindungen wurden über eine Docking-Studie analysiert, wobei die Verbindungen aus Abbildung 15 identifiziert werden konnten. Verbindung 22 wurde ohne weitere Variationen synthetisiert. Die Verbindungen 71 und 86 wurden aus der modellierten Verbindung 87, welche von DOGS vorgeschlagen wurde, weiterentwickelt. Dabei wurde Verbindung 71 als ein Fragment von Verbindung 85 entwickelt. Verbindung 86 wurde direkt aus Verbindung 87 entwickelt, wobei einige Variationen durchgenommen wurde. Hierbei sollte vorallem die Form von Verbindung 87 beibehalten werden.
Literatur verwendet, um ausgehend von den Verbindungen APHS und ASS neue COX-2 Inhibitoren zu entwickeln (siehe Abb. 16). Dabei wurden mehrere Verbindungen designt, wovon Verbindung 3 als ein leichter Inhibitor identifiziert werden konnte. 3 enthält keine für COX-Inhibitoren typische polare Gruppe, besitzt dafür aber eine Acetylgruppe, die gemäß in silico Untersuchungen in der Lage sein könnte, Ser530 in COX-2 zu acetylieren.
In der letzten Studie wurden mit Hilfe eines Fragment-basierten Designs neue Verbindungen entwickelt, wobei das Benzensulfonamid von Celecoxib aus der Kristallstruktur 1PTH extrahiert und mit kleinen Fragmente verknüpft wurde, welche zuvor über eine Docking-Studie analysiert wurden. Hieraus entwickelte sich Verbindung 35, die in einer kleinen SAR-Studie zu 70 optimiert werden konnte. Dabei konnte das Sulfonamid, welches typisch für Coxibe ist, gegen eine Carbonsäure ausgetauscht werden (69). Erst durch eine Vergrößerung der Verbindung um einen Benzyl-Rest am sekundären Amin von 69 führte zur potenten Verbindung 70.
Zusammenfassend konnten in dieser Arbeit fünf neue COX-2 Inhibitoren als Leitstrukturen entwickelt werden. Dabei kamen fortschrittliche in silico Methoden wie die De-Novo Design Software DOGS aber auch rationale Designmethoden zum Einsatz. Beide Methoden boten Vor- und Nachteile und haben jeweils zu guten Ergebnissen geführt. Bei der Entwicklung der vielversprechendsten Leitstrukturen 70 und 71 wurden die Vorteile beider Ansätze kombiniert.
Die vorliegende, in kumulativer Schreibweise verfasste Arbeit erläutert die Entwicklung, Charakterisierung und Optimierung zweier unterschiedlicher Leitstrukturen, die als Agonisten von Peroxisomen Proliferator-aktivierten Rezeptoren (PPAR) und gleichsam als duale Inhibitoren der mikrosomalen Prostaglandin E2 Synthase-1 (mPGES-1) und der 5-Lipoxygenase (5-LO) wirken. Chemisch betrachtet sind dies zum ersten die Gruppe der alpha-n-Hexyl-Pirinixinsäurederivate und zum zweiten die Gruppe der 2-(Phenylthio)-hexansäurederivate. Die Publikation zur Synthese und in vitro-pharmakologischen Charakterisierung der alpha-n-Hexyl-Pirinixinsäurederivate an PPAR (Zettl et al., QSAR & Combinatorial Science, 28:576–586, 2009) enthält einerseits die strukturelle Optimierung durch Variation der Aryl-Substitution des zentralen Pyrimidinringes der Leitstruktur und andererseits die durch Docking-Verfahren gestützte Untersuchung des Einflusses der Stereochemie auf die PPAR-Aktivierung. Letztlich konnte durch die Einführung von Biphenyl-Substituenten eine Verbesserung insbesondere der PPARalpha-Aktivität gegenüber der als strukturellen Referenz dienenden alpha-n-Hexyl-Pirinixinsäure (Rau et al., Archiv der Pharmazie, 341:191–195, 2008) erreicht werden. Mit Hilfe von präparativer enantioselektiver HPLC wurde eine ausgewählte Verbindung in ihre beiden Enantiomere getrennt. Deren in vitro-pharmakologische Charakterisierung ergab, dass das (R)-Enantiomer insbesondere bei PPARalpha als Eutomer fungiert. Dieses Ergebnis konnte mit Hilfe von Docking-Studien weiter untermauert werden. Hierbei wurde deutlich, dass die Besetzung der linken proximalen Bindetasche der PPARalpha-Liganden-Bindungs-Domäne durch den alpha-n-Hexyl-Rest lediglich im Fall einer (R)-Konfiguration optimal erfolgen kann. Die Synthese und die in vitro-pharmakologische Charakterisierung der Substanzklasse der 2-(Phenylthio)-hexansäurederivate an PPAR sind in Zettl et al., Bioorganic & Medicinal Chemistry Letters, 19: 4421-4426, 2009 zusammengefasst. Bei der Analyse der Struktur-Wirkungs-Beziehungen erwies sich die Leitstruktur als hochaktiv und sehr robust. Je nach Substitutionsmuster des lipophilen Molekülteils wurden potente selektive PPARalpha-Agonisten wie auch PPARalpha-präferenzielle duale PPARalpha/gamma-Agonisten dargestellt. Durch die Synthese von Kohlenstoff-Analoga und alpha-unsubstituierten Verbindungen wurde des Weiteren der Einfluss des Schwefelatoms und des n-Butylrestes in alpha-Position zur Carbonsäure auf die PPAR-Aktivität untersucht. Hierbei konnte gezeigt werden, dass beide Strukturelemente einen großen Beitrag zur hohen PPARalpha-Aktivität der Leitstruktur leisten. Wie auch bei den alpha-n-Hexyl-Pirinixinsäurederivaten wurde eine ausgewählte Verbindung in ihre Enantiomere getrennt und der Einfluss des Stereozentrums in alpha-Position zur Carbonsäure untersucht. Das Ergebnis bestätigte die Resultate der vorangegangenen Studie: Das (R)-Enantiomer wirkte als Eutomer, wobei der stereochemische Einfluss bei PPARalpha besonders deutlich war. Ausgewählte Synthesen und die in vitro-pharmakologische Charakterisierung von Pirinixinsäurederivaten an mPGES-1, 5-LO sowie der Cyclooxygenase (COX) sind in Koeberle und Zettl et al., Journal of Medicinal Chemistry, 51:8068–8076, 2009 publiziert. Die Arbeit beinhaltet eine umfassende Reihe an Pirinixinsäurederivaten mit Strukturvariationen in alpha-Position zur Carbonsäure und im Aryl-Substitutionsmuster des Pyrimidinringes. Hinsichtlich der alpha-Substitution zeigte sich, dass für Alkylreste eine Kettenlänge von mindestens 6 Kohlenstoffatomen für einen dualen Wirkmechanismus erforderlich ist. Als Leitstruktur für duale mPGES-1/5-LO-Inhibitoren ergab sich somit alpha-n-Hexyl-substituierte Pirinixinsäure, deren Aryl-Substitutionsmuster am zentralen Pyrimidin weiter optimiert wurde. Als vorteilhaft erwies sich die Substitution mit Biphenylresten, wodurch die Darstellung von niedrig mikromolar aktiven dualen mPGES-1/5-LO-Inhibitoren gelang. Bei der Analyse der Strukur-Wirkungs-Beziehungen von unterschiedlichen Biphenylresten zeigte sich eine hohe strukturelle Toleranz hinsichtlich der dualen inhibitorischen Aktivität an der mPGES-1 und der 5-LO. Somit stellen die alpha-n-Hexyl-Pirinixinsäurederivate die ersten publizierten dualen mPGES-1/5-LO-Inhibitoren dar.
The goal of this thesis was the development, evaluation and application of novel virtual screening approaches for the rational compilation of high quality pharmacological screening libraries. The criteria for a high quality were a high probability of the selected molecules to be active compared to randomly selected molecules and diversity in the retrieved chemotypes of the selected molecules to be prepared for the attrition of single lead structures. For the latter criterion the virtual screening approach had to perform “scaffold hopping”. The first molecular descriptor that was explicitly reported for that purpose was the topological pharmacophore CATS descriptor, representing a correlation vector (CV) of all pharmacophore points in a molecule. The representation is alignment-free and thus renders fast screening of large databases feasible. In a first series of experiments the CATS descriptor was conceptually extended to the three-dimensional pharmacophore-pair CATS3D descriptor and the molecular surface based SURFCATS descriptor. The scaling of the CATS3D descriptor, the combination of CATS3D with different similarity metrics and the dependence of the CATS3D descriptor on the threedimensional conformations of the molecules in the virtual screening database were evaluated in retrospective screening experiments. The “scaffold hopping” capabilities of CATS3D and SURFCATS were compared to CATS and the substructure fingerprint MACCS keys. Prospective virtual screening with CATS3D similarity searching was applied for the TAR RNA and the metabotropic glutamate receptor 5 (mGlur5). A combination of supervised and unsupervised neural networks trained on CATS3D descriptors was applied prospectively to compile a focused but still diverse library of mGluR5 modulators. In a second series of experiments the SQUID fuzzy pharmacophore model method was developed, that was aimed to provide a more general query for virtual screening than the CATS family descriptors. A prospective application of the fuzzy pharmacophore models was performed for TAR RNA ligands. In a last experiment a structure-/ligand-based pharmacophore model was developed for taspase1 based on a homology model of the enzyme. This model was applied prospectively for the screening for the first inhibitors of taspase1. The effect of different similarity metrics (Euc: Euclidean distance, Manh: Manhattan distance and Tani: Tanimoto similarity) and different scaling methods (unscaled, scaling1: scaling by the number of atoms, and scaling2: scaling by the added incidences of potential pharmacophore points of atom pairs) on CATS3D similarity searching was evaluated in retrospective virtual screening experiments. 12 target classes of the COBRA database of annotated ligands from recent scientific literature were used for that purpose. Scaling2, a new development for the CATS3D descriptor, was shown to perform best on average in combination with all three similarity metrics (enrichment factor ef (1%): Manh = 11.8 ± 4.3, Euc = 11.9 ± 4.6, Tani = 12.8 ± 5.1). The Tanimoto coefficient was found to perform best with the new scaling method. Using the other scaling methods the Manhattan distance performed best (ef (1%): unscaled: Manh = 9.6 ± 4.0, Euc = 8.1 ± 3.5, Tani = 8.3 ± 3.8; scaling1: Manh = 10.3 ± 4.1, Euc = 8.8 ± 3.6, Tani = 9.1 ± 3.8). Since CATS3D is independent of an alignment, the dependence of a “receptor relevant” conformation might also be weaker compared to other methods like docking. Using such methods might be a possibility to overcome problems like protein flexibility or the computational expensive calculation of many conformers. To test this hypothesis, co-crystal structures of 11 target classes served as queries for virtual screening of the COBRA database. Different numbers of conformations were calculated for the COBRA database. Using only a single conformation already resulted in a significant enrichment of isofunctional molecules on average (ef (1%) = 6.0 ± 6.5). This observation was also made for ligand classes with many rotatable bonds (e.g. HIV-protease: 19.3 ± 6.2 rotatable bonds in COBRA, ef (1%) = 12.2 ± 11.8). On average only an improvement from using the maximum number of conformations (on average 37 conformations / molecule) to using single conformations of 1.1 fold was found. It was found that using more conformations actives and inactives equally became more similar to the reference compounds according to the CATS3D representations. Applying the same parameters as before to calculate conformations for the crystal structure ligands resulted in an average Cartesian RMSD of the single conformations to the crystal structure conformations of 1.7 ± 0.7 Å. For the maximum number of conformations, the RMSD decreased to 1.0 ± 0.5 Å (1.8 fold improvement on average). To assess the virtual screening performance and the scaffold hopping potential of CATS3D and SURFACATS, these descriptors were compared to CATS and the MACCS keys, a fingerprint based on exact chemical substructures. Retrospective screening of ten classes of the COBRA database was performed. According to the average enrichment factors the MACCS keys performed best (ef (1%): MACCS = 17.4 ± 6.4, CATS = 14.6 ± 5.4, CATS3D = 13.9 ± 4.9, SURFCATS = 12.2 ± 5.5). The classes, where MACCS performed best, consisted of a lower average fraction of different scaffolds relative to the number of molecules (0.44 ± 0.13), than the classes, where CATS performed best (0.65 ± 0.13). CATS3D was the best performing method for only a single target class with an intermediate fraction of scaffolds (0.55). SURFCATS was not found to perform best for a single class. These results indicate that CATS and the CATS3D descriptors might be better suited to find novel scaffolds than the MACCS keys. All methods were also shown to complement each other by retrieving scaffolds that were not found by the other methods. A prospective evaluation of CATS3D similarity searching was done for metabotropic glutamate receptor 5 (mGluR5) allosteric modulators. Seven known antagonists of mGluR5 with sub-micromolar IC50 were used as reference ligands for virtual screening of the 20,000 most drug-like compounds – as predicted by an artificial neural network approach – of the Asinex vendor database (194,563 compounds). Eight of 29 virtual screening hits were found with a Ki below 50 µM in a binding assay. Most of the ligands were only moderately specific for mGluR5 (maximum of > 4.2 fold selectivity) relative to mGluR1, the most similar receptor to mGluR5. One ligand exhibited even a better Ki for mGluR1 than for mGluR5 (mGluR5: Ki > 100 µM, mGluR1: Ki = 14 µM). All hits had different scaffolds than the reference molecules. It was demonstrated that the compiled library contained molecules that were different from the reference structures – as estimated by MACCS substructure fingerprints – but were still considered isofunctional by both CATS and CATS3D pharmacophore approaches. Artificial neural networks (ANN) provide an alternative to similarity searching in virtual screening, with the advantage that they incorporate knowledge from a learning procedure. A combination of artificial neural networks for the compilation of a focused but still structurally diverse screening library was employed prospectively for mGluR5. Ensembles of neural networks were trained on CATS3D representations of the training data for the prediction of “mGluR5-likeness” and for “mGluR5/mGluR1 selectivity”, the most similar receptor to mGluR5, yielding Matthews cc between 0.88 and 0.92 as well as 0.88 and 0.91 respectively. The best 8,403 hits (the focused library: the intersection of the best hits from both prediction tasks) from virtually ranking the Enamine vendor database (ca. 1,000,000 molecules), were further analyzed by two self-organizing maps (SOMs), trained on CATS3D descriptors and on MACCS substructure fingerprints. A diverse and representative subset of the hits was obtained by selecting the most similar molecules to each SOM neuron. Binding studies of the selected compounds (16 molecules from each map) gave that three of the molecules from the CATS3D SOM and two of the molecules from the MACCS SOM showed mGluR5 binding. The best hit with a Ki of 21 µM was found in the CATS3D SOM. The selectivity of the compounds for mGluR5 over mGluR1 was low. Since the binding pockets in the two receptors are similar the general CATS3D representation might not have been appropriate for the prediction of selectivity. In both SOMs new active molecules were found in neurons that did not contain molecules from the training set, i. e. the approach was able to enter new areas of chemical space with respect to mGluR5. The combination of supervised and unsupervised neural networks and CATS3D seemed to be suited for the retrieval of dissimilar molecules with the same class of biological activity, rather than for the optimization of molecules with respect to activity or selectivity. A new virtual screening approach was developed with the SQUID (Sophisticated Quantification of Interaction Distributions) fuzzy pharmacophore method. In SQUID pairs of Gaussian probability densities are used for the construction of a CV descriptor. The Gaussians represent clusters of atoms comprising the same pharmacophoric feature within an alignment of several active reference molecules. The fuzzy representation of the molecules should enhance the performance in scaffold hopping. Pharmacophore models with different degrees of fuzziness (resolution) can be defined which might be an appropriate means to compensate for ligand and receptor flexibility. For virtual screening the 3D distribution of Gaussian densities is transformed into a two-point correlation vector representation which describes the probability density for the presence of atom-pairs, comprising defined pharmacophoric features. The fuzzy pharmacophore CV was used to rank CATS3D representations of molecules. The approach was validated by retrospective screening for cyclooxygenase 2 (COX-2) and thrombin ligands. A variety of models with different degrees of fuzziness were calculated and tested for both classes of molecules. Best performance was obtained with pharmacophore models reflecting an intermediate degree of fuzziness. Appropriately weighted fuzzy pharmacophore models performed better in retrospective screening than CATS3D similarity searching using single query molecules, for both COX-2 and thrombin (ef (1%): COX-2: SQUID = 39.2., best CATS3D result = 26.6; Thrombin: SQUID = 18.0, best CATS3D result = 16.7). The new pharmacophore method was shown to complement MOE pharmacophore models. SQUID fuzzy pharmacophore and CATS3D virtual screening were applied prospectively to retrieve novel scaffolds of RNA binding molecules, inhibiting the Tat-TAR interaction. A pharmacophore model was built up from one ligand (acetylpromazine, IC50 = 500 µM) and a fragment of another known ligand (CGP40336A), which was assumed to bind with a comparable binding mode as acetylpromazine. The fragment was flexible aligned to the TAR bound NMR conformation of acetylpromazine. Using an optimized SQUID pharmacophore model the 20,000 most druglike molecules from the SPECS database (229,658 compounds) were screened for Tat-TAR ligands. Both reference inhibitors were also applied for CATS3D similarity searching. A set of 19 molecules from the SQUID and CATS3D results was selected for experimental testing. In a fluorescence resonance energy transfer (FRET) assay the best SQUID hit showed an IC50 value of 46 µM, which represents an approximately tenfold improvement over the reference acetylpromazine. The best hit from CATS3D similarity searching showed an IC50 comparable to acetylpromazine (IC50 = 500 µM). Both hits contained different molecular scaffolds than the reference molecules. Structure-based pharmacophores provide an alternative to ligand-based approaches, with the advantage that no ligands have to be known in advance and no topological bias is introduced. The latter is e.g. favorable for hopping from peptide-like substrates to drug-like molecules. A homology model of the threonine aspartase taspase1 was calculated based on the crystal structures of a homologous isoaspartyl peptidase. Docking studies of the substrate with GOLD identified a binding mode where the cleaved bond was situated directly above the reactive N-terminal threonine. The predicted enzyme-substrate complex was used to derive a pharmacophore model for virtual screening for novel taspase1 inhibitors. 85 molecules were identified from virtual screening with the pharmacophore model as potential taspase1- inhibitors, however biochemical data was not available before the end of this thesis. In summary this thesis demonstrated the successful development, improvement and application of pharmacophore-based virtual screening methods for the compilation of molecule-libraries for early phase drug development. The highest potential of such methods seemed to be in scaffold hopping, the non-trivial task of finding different molecules with the same biological activity.
Das Enzym 5-Lipoxygenase (5-LO) spielt eine essentielle Rolle in der Biosynthese der Leukotriene, bioaktiver Metabolite der Arachidonsäure (AA), die an einer Vielzahl entzündlicher und allergischer Erkrankungen beteiligt sind. Die 5-LO wird bevorzugt in Zellen myeloiden Ursprungs wie Granulozyten, Monozyten oder B-Lymphozyten exprimiert. In die Regulation der zellulären 5-LO-Aktivität in der Epstein-Barr Virus-transformierten B-lymphozytären Zelllinie BL41-E95-A sind Caspasen, Aspartat-spezifische Cysteinproteasen, involviert. Das Passagieren von BL41-E95-A führt zu einer Erhöhung der Proliferationsrate der B-Lymphozyten sowie zu einem deutlichen Verlust der 5-LO-Aktivität, der mit dem Auftreten eines 62 kDa-Spaltproduktes der 5-LO und einer signifikanten Aktivitätserhöhung der Caspase-8 und -6 korreliert. Isolierte humane 5-LO wird durch rekombinante Caspase-6 zwischen Asp170 und Ser171 zu einem 58 kDa-Fragment in vitro gespalten, wobei das Tetrapeptid VEID170 innerhalb der 5-LO als Erkennungsmotiv für den Angriff der Caspase-6 dient. In einigen weiteren untersuchten Zelllinien wie Mono Mac 6 (MM6), RBL-1, PMNL oder HeLa, die nicht den B-Lymphozyten angehören, konnte die 5-LO-Spaltung weder durch das Passagieren von Zellen noch durch die Behandlung mit diversen proapoptotischen Agentien ausgelöst werden. Laut Ergebnissen aus in vitro-Untersuchungen scheinen 5-LO-positive HeLa- bzw. MM6-Zellen einen Faktor zu exprimieren, der die 5-LO direkt oder indirekt vor dem Angriff der Caspase-6 und anschließender Prozessierung schützt. Die in den BL41-E95-A-Zellen beobachtete Aktivierung der Caspasen mit anschließender Prozessierung der 5-LO lässt sich durch zwei Pflanzeninhaltsstoffe supprimieren, das Hyperforin (HP) aus Johanniskraut-Extrakten und das Myrtucommulon (MC) aus Myrte-Blättern. Beide Verbindungen scheinen in B-Lymphozyten zu einer Hemmung der Caspasen-Aktivierung zu führen. Nichtsdestotrotz führt die Behandlung der B-Lymphozyten mit HP bzw. MC zu einem apoptotischen Tod der Zellen. Offensichtlich wird dabei ein (unbekannter) einzigartiger Mechanismus der Apoptose-Induktion ausgelöst. In der vorliegenden Arbeit konnte zum ersten Mal eine potente Apoptose-induzierende Wirkung des natürlich vorkommenden Myrtucommulons auf Krebszelllinien gezeigt werden. In allen getesteten Krebszelllinien führte Myrtucommulon zum Zelltod, wobei die HL-60-Zellen mit einem IC50-Wert von 3,26 ± 0,51 µM MC am sensitivsten gegenüber MC-Einfluss waren. Zusätzlich konnte in HL-60- und MM6-Zellen nach MC-Behandlung neben einer erhöhten Caspasen-Aktivität und PARP-Spaltung ein signifikanter DNA-Abbau detektiert werden. Von besonderer Bedeutung ist die Tatsache, dass die zytotoxische MC-Wirkung eine bemerkenswerte Selektivität für entartete Zelllinien zu besitzen scheint und gegenüber nicht-transfizierten Zellen minimal ist.
The goal of this thesis was to gain further insight into the binding behavior of ligands in the heptahelical domain (HD) of group I metabotropic glutamate receptors (mGluRs). This was realized by the establishment of strategies for the detection and optimization of molecules acting as non-competitive antagonists of group I mGluRs (mGluR1/5). These strategies should guarantee high diversity in the retrieved chemotypes of the detected compounds not resembling original reference molecules (“scaffold-hopping”). The detection of new scaffolds, in turn, was divided into two approaches: First the development of pharmacological assays to screen compounds at a certain target for bioactivity (here: affinity towards the allosteric recognition site of mGluR1 and mGluR5), and second the evaluation of computer assisted methods for the identification of virtual hits to be screened afterwards on the pharmacological assays established before. Promising molecules should be optimized with respect to activity/affinity and selectivity, their binding mode investigated and, finally, compared to existing lead compounds. Initially, membrane based binding assays for the HD of mGlu1 and mGlu5 receptors with enhanced throughput (shifting from 24-well plates to 96-well plates) were set up. For the mGluR1 assay the potent antagonist EMQMCM exhibited high affinity towards the binding site (Ki ~3nM), which is in accordance with published data from Mabire et al. (functional IC50 3nM). For mGluR5 the reference antagonist MPEP binds with high affinity to the receptor (binding IC50 13.8nM), which confirmed earlier findings from Anderson et al. (binding IC50 15nM). In another series of experiments the properties of rat cerebellar (mGluR1) and corticalmembranes (mGluR5) as well as of radiotracers were investigated by means of binding saturation studies and kinetic experiments. Furthermore, the influence of the solvent DMSO, necessary for compound screening of lipophilic substances, on positive and negative controls was evaluated. As the precise architecture of the HD of mGluR1 is still not known our efforts in identifying new ligands for this receptor focused on the ligand-based approach. All computer assisted methods that were applied to virtually screen large compound collections and to retrieve potential hits (“activity-enriched subsets”) acting at the heptahelical domain of mGluR1 relied on the existence of a valid dataset of reference molecules. This was realized by an initial compilation of a mGluR reference data collection comprising in total 357 entries predominantly negative but also some positive allosteric modulators for mGluR1 and mGluR5. In the next step a pharmacophore model for non-competitive mGluR1 antagonists was constructed. It was based upon six selective, potent and structurally diverse ligands. Prospective virtual screening was performed using the CATS atom-pair descriptor. The Asinex Gold-Collection was screened for each seed compound and some of the most similar compounds (according to the CATS descriptor) were ordered and tested forbinding affinity and functional activity at mGluR1. A high hit rate of approximately 26% (IC50 < 15 micro M) was yielded confirming the applicability of this method. One compound exerted functional activity below one micro molar (IC50-value of C-07:362nM ± 0.03). Moreover, non-linear principal component analysis was employed. Again the Asinex vendor database served as test database and was filtered by the pharmacophore model for mGluR1 established before. Test molecules that were adjacently located with mGluR1 antagonist references were selected. 15 compounds were tested on mGluR1 in binding and functional assays and three of them exhibited functional activity (IC50) below 15 micro M. The most potent molecule P-06 revealed an IC50-value of 1.11 micro M (± 0.41). The COBRA database comprising 5,376 structurally diverse bioactive molecules affecting various targets was encoded with the CATS descriptor and used for training two selforganizing maps (SOM). The encoded mGluR reference data collection was projected onto this map according to the SOM algorithm. This projection allowed to clearly distinguish between antagonists of mGluR1 and mGluR5 subtype. 28 compounds were ordered and tested on activity and affinity for mGluR1. They exhibited functional activity down to the sub-micro molar range (IC50-value of S-08: 744nM ± 0.29) yielding a final hit rate of 46% (<15 micro M). Then, the Asinex collection was screened using the SOM approach. For a predicted target panel including the muscarinic mACh (M1) receptor, the histamine H1-receptor and the dopamine D2/D3 receptors, the tested mGluR ligands exhibited the calculated binding pattern. This virtual screening concept might provide a basis for early recognition of potential sideeffects in lead discovery. We superimposed a set of 39 quinoline derivatives as non-competitive mGluR1 antagonists that were recently published by Mabire and co-workers. A CoMFA model (QSAR) was established and the influence of several side chains on functional activity was investigated. The coumarine derivative C-07 was obtained as a result of similarity searching. Starting from this compound a series of chemical derivatives was synthesized. This led to the discovery of potent (B-28, IC50: 58nM ± 0.008; Ki: 293nM ± 0.022) and selective (rmGluR5 IC50: 28.6 micro M) mGluR1 antagonists. From a homology model of mGluR1 we derived a potential binding mode for coumarines within the allosteric transmembrane region. Potential interacting patterns with amino acids were proposed considering the difference of the binding pockets between rat and human receptors. The proposed binding modes for quinolines (here:EMQMCM) and coumarines (here:B-04) were compared and discussed considering in particular the influence on activity of several side chains of quinolines obtained from the QSAR studies. The present studies demonstrated the applicability of ligand-based virtual screening for non-competitive antagonists of a G-protein coupled receptor, resulting in novel, potent and selective agents.