EP1456646A2 - Computersystem und verfahren zur berechnung von adme-eigenschaften - Google Patents
Computersystem und verfahren zur berechnung von adme-eigenschaftenInfo
- Publication number
- EP1456646A2 EP1456646A2 EP02783091A EP02783091A EP1456646A2 EP 1456646 A2 EP1456646 A2 EP 1456646A2 EP 02783091 A EP02783091 A EP 02783091A EP 02783091 A EP02783091 A EP 02783091A EP 1456646 A2 EP1456646 A2 EP 1456646A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- substance
- properties
- adme
- database
- molecular properties
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/40—Searching chemical structures or physicochemical data
Landscapes
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Complex Calculations (AREA)
Abstract
Die Erfindung betrifft ein Computersystem und ein Verfahren zur Berechnung einer ADME-Eigenschaft einer Substanz mit folgenden Schritten: a. Eingabe von molekularen Eigenschaften der Substanz in ein biophysikalisches Modell, wobei das biophysikalische Modell einen Zusammenhang zwischen den molekularen Eigenschaften und der ADME-Eigenschaft beschreibt, b. Ausgabe der ADME-Eigenschaft.
Description
Computersystem und Verfahren zur Berechnung von ADME-Eigenschaften
Die Erfindung betrifft ein Computersystem und ein Verfahren zur Berechnung von ADME-Eigenschaften einer Substanz, insbesondere für eine Substanz mit pharma- kologischer Wirkung oder einer Substanz für Pflanzenschutzanwendungen.
Die Wirksamkeit von Wirkstoffen wird bestimmt durch deren Wechselwirkung mit dem molekularen, biologischen Target, sowie durch die Konzentration am Zielort. Beide Größen werden in der Regel von unterschiedlichen molekularen Parametern bestimmt und können daher in gewissen Grenzen unabhängig voneinander optimiert werden. Während die intrinsische biochemische Wirkung in in-vitro Tests in einem sehr f ühen Forschungsstadium für großen Anzahlen von Substanzen bestimmt werden kann, lässt sich die Konzentration im Wirkort nur über Versuche am ganzen Organismus (Tier, Pflanze oder Pilz) untersuchen. Letzteres bedeutet, dass die Informationen durch die Aufwändigkeit der Versuche erst in späten Forschungsstadien durchgeführt werden können, und daher für die ersten Optimierungszyklen nicht zur Verfügung stehen.
Dies hat in den letzten Jahren dazu geführt, dass in der pharmazeutischen und
Pflanzenschutz-Industrie nach alternativen Wegen gesucht wurde, frühzeitig Informationen über das ADME- Verhalten (Absorption, Distribution, Metabolismus, Exkretion) von Wirkstoffen zu gewinnen. Da große Teile des ADME- Verhaltens durch einfach messbare physikochemische Eigenschaften, bzw. aus der chemischen Struktur berechenbare Größen beeinflusst werden, hat sich inzwischen die
Vorgehensweise etabliert, solche Größen mit hohem Durchsatz experimentell zu bestimmen, bzw. zu berechnen [H. van de Waterbeemd, D. A. Smith, K. Beaumont, D. K. Walker, J. Med. Chem. 44, 1-21 (2001)].
Typische Eigenschaften die dabei üblicherweise berücksichtigt werden, sind z. B.
Lipophilie, Wasserlöslichkeit, Permeabilitäten über künstliche Membranen oder
Zellschichten, Molekülgewicht und Anzahlen bestimmter Stakturmerkmale wie Wasserstoffdonoren und -akzeptoren. Zur Beurteilung der Substanzen dient dann in der Regel die Einhaltung bestimmter Grenzen, welche üblicherweise aus Erfahrungswerten oder aus der statistischen Verteilung der Eigenschaften für Handelsprodukte gewonnen werden [C. A. Lipinski, F. Lombardo, B. W. Dominy, P. J. Feeney, Adv.
Drug Delivery Rev. 23, 3-25 (1997) und C. M. Tice, Pest Management Sei. 57, 3-16 (2001)].
Nachteil dieses Verfahrens ist, dass starre Grenzen für einzelne nur indirekt relevante Eigenschaften betrachtet werden. Die eigentlich wichtigen ADME-Eigenschaften hängen aber in der Regel von mehreren dieser Größen gleichzeitig ab, was dazu führt, auch die tolerierbaren Grenzen der Einzelgrößen von deren Wert abhängig werden, und daher absolute Grenzwerte nur sehr grob festzulegen sind.
Aus US-A-5.901,069 ist ein Computerprogrammprodukt zur zumindest teilweise automatischen Berechnung von molekularen Größen basierend auf einer Substanzbibliothek bekannt. Dieses Verfahren ermöglicht jedoch nicht die Berechnung von ADME-Eigenschaften.
Aus US-A-5 680 590 ist ein Simulationssystem zur Berechnung physiologischer
Daten hinsichtlich pharmakokinetischer und pharmakodynamischer Parameter bekannt. Nachteilig bei diesem Simulationssystem ist, dass die Berechnung nicht direkt basierend auf molekularen Eigenschaften einer zu bewertenden Substanz beruht. Dieses Modell ist nur für Ausbildungszwecke geeignet und erlaubt nicht die Berechnung der ADME-Eigenschaften einer neuen Substanz.
Der Erfindung liegt die Aufgabe zu Grunde ein verbessertes Verfahren zur Berechnung einer ADME-Eigenschaft einer Substanz sowie ein entsprechendes Computerprogramm und Computersystem zu schaffen.
Die der Erfindung zu Grunde liegende Aufgabe wird mit den Merkmalen der unabhängigen Patentansprüche jeweils gelöst.
Bevorzugte Ausführungsformen der Erfindung sind in den abhängigen Patent- ansprächen angegeben.
Die Erfindung kann vorteilhafterweise sowohl für bereits synthetisierte Substanzen als auch für „virtuelle Strukturen" von Substanzbibliotheken verwendet werden. Für die Berechnung der ADME-Eigenschaften werden erfindungsgemäß biophysi- kaiische Zusammenhänge zwischen molekularen und physikochemischen Eigenschaften der Substanzen und den betrachteten ADME-Eigenschaften verwendet. Die mathematische Beschreibung erfolgt dann entweder durch eine analytische Formel oder bei komplexeren Zusammenhängen durch eine numerische Simulation.
Ein besonderer Vorteil der Erfindung liegt darin ,dass die Berechnung von ADME-
Eigenschaften unmittelbar basierend auf den molekularen Eigenschaften der zu untersuchenden Substanz erfolgt. Hierzu wird ein biophysikalisches Modell verwendet, in welches die molekularen Eigenschaften als Eingangsgrößen eingegeben werden. Das biophysikalische Modell stellt einen Zusammenhang zwischen diesen molekularen Eigenschaften und der oder den ADME-Eigenschaften her. Auf diese Art und Weise erhält man eine unmittelbare Aussage zu den ADME- Eigenschaften und nicht lediglich zu korrelierenden Surrogat-Größen, wie das im Stand der Technik der Fall ist.
Entsprechend einer bevorzugten Ausführungsform der Erfindung werden als molekulare Eigenschaften die Lipophilie, die Bindung an Proteine und die Molekülgröße verwendet. Beispielsweise erlaubt dies die Modellierung der Absorption von Wirkstoffen in einem menschlichen, tierischen oder pflanzlichen Organismus. Die Lipophilie wird dabei durch Verteilungskoeffizienten zwischen einer lipoiden Phase (z.B. Oktanol, Speiseöl, Hexan, Phospholipid-Membranen) und
Wasser beschrieben. Als Maß für die Molekülgröße können z.B. das Molgewicht oder das Molvolumen dienen.
Da diese in der Regel von einem Permeationsprozess bestimmt wird, und die Permeabilität bekanntermaßen von der Lipophilie und der Größe des permeierenden
Moleküls abhängt, können aus aus diesen Eigenschaften Schlüsse auf die Absorptionsrate gezogen werden. Da die Abhängigkeiten aber gegenläufig sind, können zum Beispiel mit steigender Lipophilie höhere Molgewichte toleriert werden, um die gleiche Absorptionsraten zu erreichen. Dieser Fakt kann aber im Stand der Technik auch durch feste Grenzwerte für optimale Lipophilie und Molekülgröße nicht ausreichend berücksichtigt werden.
Die Erfindung erlaubt dagegen unmittelbar von den molekularen Eigenschaften auf die ADME-Eigenschaften zu schließen, so dass eine Berechnung mit verbesserter Genauigkeit möglich ist.
Nach einer bevorzugten Ausführungsform der Erfindung handelt es sich bei dem biophysikalischen Modell um ein physiologisch basiertes pharmakokinetisches Modell. Für die Untersuchung der ADME-Eigenschaften einer Substanz bei- spielsweise im menschlichen Körper umfasst das Modell zumindest die für die
Untersuchung wesentlichen Organe, wie z. B. die Lunge, Leber und Nieren sowie den Blutkreislauf. Die verschiedenen Teilmodelle der Organe werden durch Massenerhaltungsgleichungen miteinander verknüpft.
Nach einer weiteren bevorzugten Ausführungsform der Erfindung werden die
Massenerhaltungsgleichungen in der Form eines Differenzialgleichungssystems ausgedrückt, wobei die Eingangsgrößen des Differenzialgleichungssystems unmittelbar aus den berechneten molekularen Eigenschaften gewonnen werden.
Nach einer weiteren bevorzugten Ausführungsform der Erfindung erfolgt die
Eingabe der chemischen Sfrukiur einer zu untersuchenden Substanz in eine Daten-
bank. Die chemische Struktur kann für diesen Zweck beispielsweise auch in der Form eines Deskriptors oder eines sogenannten Fingerprints dargestellt werden.
Die Eingabe der chemischen Struktur kann dabei dezentral von einem Client- Computer erfolgen. Der Client-Computer befindet sich beispielsweise unmittelbar an dem Arbeitsplatz eines Chemikers zur Eingabe von neuen chemischen Strukturen, für die die ADME-Eigenschaften vorherbestimmt werden sollen.
Die Abfrage dieser Datenbank erfolgt dann zyklisch, beispielsweise durch einen Server-Computer. Sobald die Eingabe einer neuen chemischen Struktur durch den
Server-Computer festgestellt worden ist, werden die molekularen Eigenschaften dieser chemischen Struktur durch einen entsprechenden Programmaufruf automatisch berechnet.
Alternative oder zusätzlich kann der Server-Computer auf eine weitere Datenbank zugreifen, in der experimentell ermittelte molekulare Eigenschaften der Substanz abgespeichert sind. Nach der Ermittlung der molekularen Eigenschaften werden diese in das biophysikalische Modell durch eine weiteren. Programmaufruf eingegeben, so dass die Berechnung der gewünschten ADME-Eigenschaften automatisch erfolgt.
Vorzugsweise wird dieser Vorgang für verschiedene Substanzen, beispielsweise desselben chemischen Grundkörpers, wiederholt durchgeführt. Die Ergebnisse der ADME-Berechnungen werden dann in einer strukturierten Form ausgegeben, beispielsweise nach dem Wert einer bestimmten ADME-Eigenschaft sortiert oder nach einem gewichteten Index von ADME-Eigenschaften sortiert.
Nach einer weiteren bevorzugten Ausführungsform der Erfindung wird ein statistisches Verfahren zur Berechnung der molekularen Eigenschaften aus der chemischen Struktur verwendet, wie zum Beispiel ein QSAR oder HQSAR Ver- fahren, oder ein auf einem neuronalen Netz basierenden Verfahren. Solche Verfahren zur Bestimmung von molekularen Eigenschaften - beispielsweise aus einem
Deskriptor der chemischen Struktur - sind an sich aus dem Stand der Technik bekannt.
Für Pharma-Anwendungen ermöglicht die Erfindung insbesondere die Berechnung folgender ADME-Eigenschaften:
- absorbierter Anteil einer Dosis nach oraler Applikation (Fraction absorbed),
- Konzentrations-Zeit-Kurven in der Portalvene nach oraler Applikation,
- freie Fraktion im Plasma,
- Organ/Plasma- Verteilungskoeffizient,
- Blut/Plasma- Verteilungskoeffizient
- Konzentrations-Zeit-Kurven im Blutplasma und in Organen nach oraler oder intravenöser Applikation.
- Aus den Konzentrations-Zeit-Kurven abgeleitete übliche pharmakokinetische Parameter wie z.B. Maximale Konzentration, Zeitpunkt der max. Konzentration, Verteilungsvolumen, Halbwertzeiten.
Für Pflanzenschutzanwendungen ermöglicht die Erfindung insbesondere die Berechnung der folgenden ADME-Eigenschaften:
- Kennzahl für Geschwindigkeit der Aufnahme in das Blatt nach Spritzapplikation,
- Kennzahl für Geschwindigkeit der Aufnahme in ein Insekt über die Kutikula,
- Kennzahl für Geschwindigkeit der Verteilung in der Pflanze nach Blattapplikation (Phloemmobilität),
- Kennzahl für Geschwindigkeit der Verteilung in der Pflanze nach Wurzel- applikation(Xylemmobilität),
- Konzentrationsverteilung zwischen den Organen eines Insekts.
Die berechneten Daten werden entweder direkt in einer Datenbank abgelegt, oder in einer Tabelle ausgegeben. Die so vorliegenden Daten bilden die Grundlage für ein
Ranking der betrachteten Substanzen oder Strukturen und die Auswahl der Kandidaten für weitere Optimierungen anhand dieses Rankings.
Für ein Rahking sind zunächst diejenigen Eigenschaften auszuwählen, die für die angestrebte Anwendung des zu optimierenden Wirkstoffs entscheidend sind. Die
Auswertung erfolgt manuell mit einer geeigneten Datenauswerte- und Visuali- sierungssof ware. Dabei werden durch Ranking die Substanzen oder Strukturen gesucht, die in der Eigenschaftsverteilung im optimalen Bereich liegen ( z.B. die 10 Substanzen mit höchster Absorption nach oraler Applikation). Werden mehr als eine Eigenschaft beim Ranking berücksichtigt, so kann ein Index berechnet werden, der dies beinhaltet (im einfachsten Fall die Summe der Werte aller Größen), dabei können auch Wichtungen der Eigenschaften nach Relevanz vorgenommen werden.
Alternativ zur manuellen Auswertung besteht auch die Möglichkeit Ausgabemasken für die Daten zu verwenden, die eine automatische Bewertung vornehmen, indem sie für jede Struktur beim Heranziehen der Daten überprüfen, ob bestimmte Werte oder auch Kombinationen von Werten (Indizes) in einem für die betrachtete Anwendung optimalen Bereich liegen (Anzeige, z. B. über Ampelfarben). Die projektspezifischen . Regeln werden dann in der Auswerte Maske hinterlegt (z.B. in einem Tabellen- kalkulationsprogramm) .
Im Weiteren werden die bevorzugten Ausführungsformen der Erfindung mit Bezugnahme auf die Zeichnungen näher erläutert. Es zeigen:
Figur 1 ein Flussdiagramm zur Berechnung von ADME-Eigenschaften aus den molekularen Eigenschaften einer Substanz,
Figur 2 eine Ausführungsform eines biophysikalischen Modells zur
Herstellung eines Zusammenhangs zwischen den molekularen Eigenschaften und den ADME-Eigenschaften,
Figur 3 eine Tabelle zur Berechnung von Organ/Blut- Verteilungskoeffizienten für das Modell der Figur 2,
Figur 4 eine Ausführungsform eines erfindungsgemäßen Computersystems.
Die Figur 1 zeigt ein Flussdiagramm für die Berechnung von ADME-Eigenschaften einer Substanz. In dem Schritt 1 wird die molekulare Eigenschaft einer Substanz ermittelt. Dies kann experimentelle erfolgen, wenn die Substanz bereits synthetisiert worden ist. Ferner kann die Bestimmung der molekularen Eigenschaften auch durch eine Berechnung erfolgen.
Hierzu sind an sich aus dem Stand der Technik Verfahren bekannt, welche es erlauben, aus der chemischen Struktur molekulare Eigenschaften zu berechnen. Beispiel für solche Verfahren sind QSAR, HQSAR und neuronale Netze. Als Eingangsgröße für solche Berechnungs verfahren dienen Deskriptoren oder
Fingerprints der chemischen Struktur der zu untersuchenden Substanz.
Ferner ist es möglich, in dem Schritt 1 sowohl auf experimentell ermittelte molekulare Eigenschaften als auch auf durch Berechnung ermittelte molekulare Eigenschaften der Substanz zuzugreifen. Auf diese Art und Weise können sich
experimentelle Verfahren mit den Berechnungs verfahren zur Bestimmung der molekularen Eigenschaften ergänzen.
Die in dem Schritt 1 ermittelten molekularen Eigenschaften der Substanz werden in dem Schritt 2 in ein biophysikalisches Modell eingegeben. Das biophysikalische
Modell stellt einen Zusammenhang her zwischen den molekularen Größen und den interessierenden ADME-Eigenschaften. Hierbei kann es sich beispielsweise um ein physiologisch basiertes pharmako-kinetisches Modell handeln. Eine Ausfüh- rungsform eines solchen biophysikalischen Modells wird unten mit Bezugnahme auf die Figuren 2 und 3 näher erläutert.
In dem Schritt 3 werden die ADME-Eigenschaften aus dem biophysikalischen Modell ausgegeben. Von besonderem Vorteil ist hierbei, dass die ADME-Eigenschaften unmittelbar aus den molekularen Eigenschaften bestimmt werden, und zwar ohne Zwischenschaltung von interpretationsbedürftigen Surrogat-Größen. Dies erlaubt einen vollautomatischen Ablauf für die Berechnung der ADME-Eigenschaften.
Die Figur 2 zeigt ein biophysikalisches Modell 4 eines Warmblüters, beispielsweise eines Menschen. Das biophysikalische Modell 4 beinhaltet eine Anzahl von Teil-
Modellen für diejenigen Organe, die für die Verteilung der Substanz im Körper am relevantesten sind. In dem Beispiel der Figur 2 sind dies das Teil-Modell 5 für Lunge, Teil-Modell 6 für die Leber, Teil-Modell 7 für die Nieren und Teil-Modelle 8 jeweils für verschiedene andere Organe X.
Die Teil-Modelle 5, 1, 8 werden durch venöses Blut 9 und arterielles Blut 10 miteinander „verschaltet". Das venöse Blut 9 tritt in das Teil-Modell 5 für die Lunge ein, wo es in arterielles Blut 10 überführt wird. Das arterielle Blut 10 gelangt dann in die weiteren Teil-Modelle 6, 1 und 8, von wo es jeweils wieder als venöses Blut 9 austritt. Die verschiedenen Teil-Modelle der Organe werden also durch das venöse
Blut 9 und das arterielle Blut 10 „parallel geschaltet".
Ferner beinhaltet das biophysikalische Modell 4 ein Ausscheidungs-Modell 11 für die Teil-Modelle 6 und 7, das heißt, für die Leber und die Nieren.
In dem betrachteten Beispielsfall der Figur 2 ist das biophysikalische Modell für die Modellierung der zeitlichen Konzentrationsverteilung einer Dosis der zu untersuchenden Substanz vorgesehen. Die Dosis wird dabei über das venöse Blut 9 zugeführt, das heißt, beispielsweise über die Portalvene nach oraler Applikation oder durch Injektion in eine Vene.
Für die Durchflussrate Qum e von venösem Blut 9 durch die Lunge 5 kann auf eine entsprechende Tabelle von Erfahrungswerten zugegriffen werden. Ebenso kann für die Durchflussraten von arteriellem Blut 10 durch die weiteren Organteil-Modelle 6, 7, 8 auf entsprechende Erfahrungswerte für die Durchflussraten O -ber, Queren und Qx zugegriffen werden.
Bei den Größen Cx handelt es sich um die Konzentration der Substanz in dem betreffenden Organ X zu einem bestimmten Zeitpunkt. Der Parameter Kx bezeichnet den Verteilungskoeffizient der Substanz zwischen Blut und dem Organ X im Gleichgewichtszustand. Die Parameter CL e er und CLNieren bezeichnen die intrinsische Ausscheidung der Leber bzw. der Nieren.
Basierend auf dem biophysikalischen Modell 4 lässt sich für jedes Organ X eine Massengleichgewichtsbeziehung durch eine Differenzialgleichung der folgenden Form aufstellen:
dC C
V - = 0 - C - O •~ -
Vχ dt üx ar ü* Kx
Vx = Volumen des Organs X Cx = Konzentration der Substanz in dem Organ X
Qx = Durchflussrate von Blut durch das Organ X fu = Anteil der Substanz, die nicht in Plasma gebunden ist
Car = Konzentration der Substanz, die das Organ über das arterielle Blut erreicht Kx - Verteilungskoeffizient der Substanz zwischen Blut und Organ X im Gleich- gewichtszustand
Der Parameter fu berechnet sich aus dem Kehrwert des Verteilungskoeffizienten der Substanz im Gleichgewicht zwischen Blutplasma und Wasser.
Die entsprechenden Differenzialgleichungen für die Leber und die Nieren beinhalten einen zusätzlichen Ter , der die Ausscheidung der Substanz beschreibt. Eine solche Differenzialgleichung ist unten für die Nieren angegeben; für die Leber gilt Entsprechendes:
Vj = Volumen der Nieren
C i = Konzentration der Substanz in den Nieren
Qki = Durchflussrate des Bluts durch die Nieren fu = Anteil der Substanz, die nicht in Plasma gebunden ist
Car = Konzentration der Substanz, die die Nieren über das arterielle Blut erreicht
Kw = Verteilungskoeffizient der Substanz zwischen Blut und Nieren im
Gleichgewichtszustand
CLw = intrinsische Ausscheidung der Nieren
Daraus lässt sich eine Gleichung für das venöse Blut aufstellen, und zwar durch
Aufsummierung aller „Output" Konzentrationen der verschiedenen Organe und der
CLk. intravenös zugefügten Dosis der Substanz. Der Term ist abhängig vom
Lipophilie-Wert der Substanz, kann also auch aus einer molekularen Eigenschaft bestimmt werden.
Die Differentialgleichung für die Lunge stellt eine Verbindung zwischen dem venösen und dem arteriellen Blut her. Die entsprechenden Gleichungen sind unten angegeben:
~ = ∑Qx -Q, - vs +dose . l(t)
V dClu - Q . r -O .^i dt Klυ
Vve = Volumen des venösen Bluts
Cve = Konzentration der Substanz in der Lunge zugeführtem venösen Blut
Qx = Durchflussraten von Blut durch das Organ X fu = Anteil der Substanz, die nicht in Plasma gebunden ist Cx = Konzentration der Substanz in dem Organ X
Kx = Verteilungskoeffizient der Substanz zwischen Blut und Organ X im Gleichgewichtszustand
Qιu = Blutdurchfluss durch die Lunge = Σ aller Blutdurchflüsse
Viu = Lungenvolumen Ciu = Konzentration der Substanz in den Lungen
Kiu = Verteilungskoeffizient der Substanz im Gleichgewicht zwischen Blut und Lungen
Dosis = Dosis der Substanz, die dem venösem Blut über die Zeit t gemäß einer Eingabefunkti on I(t) zugefügt wird.
Zur Lösung des hieraus resultierenden Differenzialgleichungssystems des biophysikalischen Modells 4 ist die Kenntnis der verschiedenen Verteilungskoeffizienten im Gleichgewichtszustand erforderlich. Diese lassen sich
aus molekularen Eigenschaften der Substanz, welche zuvor experimentell ermittelt oder berechnet worden sind, bestimmen.
Im Wege der Berechnung lassen sich für eine Substanz die Verteilungskoeffizienten im Gleichgewicht zwischen Fett und Wasser (KFett) und zwischen Protein und Wasser (Kprotein) bestimmen. Die genannten Verteilungskoeffizienten werden entweder rechnerisch oder experimentell ermittelt; sowohl für die rechnerische als auch die experimentelle Ermittlung gibt es an sich aus dem Stand der Technik bekannte Verfahren.
Ferner werden für die Berechnung die Organzusammensetzungen aus den Bestandteilen Wasser, Fett und Protein herangezogen. Diese können aus der Tabelle gemäß Figur 3 ermittelt werden. Die Tabelle der Figur 3 gibt für verschiedene Organe die Gesamtmasse des betreffenden Organs sowie die absoluten Anteile in Gramm von Wasser, Fett und Protein in dem Organ an.
Daraus lassen sich Fwasser, Fpett und Fprotein bestimmen, wobei
F asser = Gesamtmasse von Wasser in dem Organ / Gesamtmasse des Organs
FFett = Gesamtmasse von Fett in dem Organ / Gesamtmasse des Organs
Fprotein = Gesamtmasse von Protein in dem Organ / Gesamtmasse des Organs
Daraus lassen sich die Verteilungskoeffizienten der Substanz zwischen einem Organ und Wasser im Gleichgewichtszustand (Koi-gan /Wasser) berechnen:
Kθrgan/Wasser= F asser+ KFett' FFett+ Kprotein' Fprotem (1)
Daraus lassen sich wiederum die Organ/Blut oder Organ/Plasma
Verteilungskoeffizienten berechnen:
Kθrgan Blut — Korgan/Wasser/ Kfilut/Wasser (2)
-^Organ/Plasma — -"^Organ/Wasser ' piasma Water ( )
Die Koeffizienten Kßiut/wasser und Kpιasma/Wasser werden ebenso nach Formel (1) berechnet.
Mit Hilfe des biophysikalischen Modells 4 (vgl. Figur 2) ist es also möglich unmittelbar aus den molekularen Eigenschaften der Substanz KFett und Kprotein die Konzentration der Substanz in einem bestimmten Organ X (Cx) sowie die
Konzentration der Substanz in dem arteriellem Blut (Car) durch Lösung der oben genannten Gleichungssysteme zu bestimmen wie auch den zeitlichen Verlauf dieser Konzentrationen.
Die Figur 4 zeigt ein Blockdiagramm eines Computersystems für die Berechnung von ADME-Eigenschaften einer Substanz.
Eine Datei 12 beinhaltet die chemische Struktur der Substanz, beispielsweise in Form eines sogenannten Deskriptors oder eines Fingerprints. Die Datei 12 wird von einem Chemiker per Hand eingegeben oder es handelt sich um einen Teil einer
Substanzbibliothek chemischer Strukturen deren ADME-Eigenschaften zu bestimmen sind.
Die Datei 12 wird in eine Datenbank 13 eingegeben. Die Datenbank 13 dient zur Speicherung von chemische Stnikturen beschreibende Dateien 12. Die Datei 13 wird zyklisch von einem Programm 14 abgefragt, und zwar dahingehend, ob eine neue Datei 12 in dem Zeitraum zwischen der vorhergehenden Abfrage eingegeben worden ist.
Beispielsweise kann die Eingabe der Datei 12 von einem Client-Computer aus vorgenommen werden. Die Datenbank 13 befindet sich beispielsweise auf einem
S erver-Computer sowie auch das Programm 14, welches die Datenbank 13 zyklisch abfragt. Auf diese Art und Weise lässt sich ein verteiltes System realisieren.
Wenn das Programm 14 feststellt, dass eine neue Datei 12 zwischenzeitlich in die Datenbank 13 eingegeben worden ist, so startet das Programm 14 automatisch ein
Programm 15 zur Berechnung der molekularen Eigenschaften der durch die Datei 12 beschriebenen Substanz. Die von dem Programm 15 berechneten molekularen Eigenschaften werden in einer Datei zwischengespeichert oder in einer Datenbank 16 gespeichert.
Nach Beendigung der Berechnung der molekularen Eigenschaften und deren Speicherung in der Datenbank 16 startet das Programm 15 automatisch ein Programm 17 für die Berechnung von einer oder mehreren ADME-Eigenschaften der Substanz. Dazu greift das Programm 17 auf die Datenbank 16 zu, um die zuvor von dem Programm 15 berechneten molekularen Eigenschaften der Substanz abzurufen.
Alternativ oder zusätzlich greift das Programm 17 auf eine Datenbank 18 zu, welche weitere experimentell ermittelte molekulare Eigenschaften der Substanz beinhaltet. Dies setzt voraus, dass die Substanz zuvor synthetisiert worden ist, um experimentell ermittelte molekulare Eigenschaften der Substanz in einer Datei 19 in die Datenbank
18 eingeben zu können.
Die von dem Programm 17 basierend auf den in der Datenbank 16 und / oder der Datenbank 18 gespeicherten molekularen Eigenschaften berechneten ADME- Eigenschaften werden in einer Datenbank 20 gespeichert. Auf die Datenbank 20 greift ein Programm 21 zu, um eine strukturierte Ausgabe zu erzeugen. Dies kann in der Form einer tabellarischen Ausgabe in Form eines Spread Sheets erfolgen. Ferner kann die Ausgabe sortiert nach bestimmten ADME-Eigenschaften erfolgen.
Bezugszeichenl iste
biophysikalisches Modell 4 Teil-Modell 5
Teil-Modell 6
Teil-Modell 7
Teil-Modell 8 venöses Blut 9 arterielles Blut 10
Ausscheidungsmodell 11
Datei 12
Datenbank 13
Programm 14 Programm 15
Datenbank 16
Programm 17
Datenbank 18
Datei 19 Datenbank 20
Programm 21
Claims
1. Verfahren zur Berechnung von ADME-Eigenschaften einer Substanz mit folgenden Schritten:
Eingabe von molekularen Eigenschaften der Substanz in ein biophysikalisches Modell, wobei das biophysikalische Modell einen Zusammenhang zwischen den molekularen Eigenschaften und der ADME-Eigenschaft beschreibt,
Ausgabe der ADME-Eigenschaft.
2. Verfahren nach Anspruch 1 mit folgenden weiteren Schritten:
- Eingabe der chemischen Struktur der Substanz in eine Datenbank,
zyklisches Abfragen der Datenbank auf eine Eingabe einer chemischen Struktur,
- automatische Berechnung von molekularen Eigenschaften der
Substanz, basierend auf der chemischen Struktur, nachdem die Eingabe in die Datenbank durch das zyklische Abfragen festgestellt worden ist,
- automatische Eingabe der molekularen Eigenschaften in das biophysikalische Modell.
3. Verfahren nach Anspruch 1 oder 2 mit folgenden weiteren Schritten:
- Zugriff auf eine zweite Datenbank mit weiteren experimentell ermittelten molekularen Eigenschaften der Substanz, automatische Eingabe der berechneten und der weiteren molekularen Eigenschaften in das biophysikalische Modell.
4. Verfahren nach Anspruch 1, 2 oder 3, wobei es sich bei den biophysikalischen Modell um ein physiologisch basiertes pharma- kokinetisches Modell handelt, wobei das Modell mehrere Organe umfasst und für jedes der Organe ein Teil-Modell beinhaltet, wobei die Teil-Modelle durch ein Blutflussmodell miteinander verbunden sind und der Zusam- - menhang durch Massengleichgewichtsbeziehungen beschrieben wird.
5. Verfahren nach einem der vorhergehenden Ansprüche 1 bis 4, wobei die molekularen Eigenschaften der Substanz aus der chemischen Struktur, vorzugsweise in Form eines Deskriptors oder eines Fingerprints, durch ein statistisches Verfahren, vorzugsweise QSAR oder HQSAR, oder mittels eines neuronalen Netzes berechnet werden.
6. Verfahren nach einem der vorhergehenden Ansprüche 1 bis 5, wobei es sich um eine Substanz mit einer pharmakologischen Wirkung auf den mensch- liehen und/oder tierischen Körper handelt und eine oder mehrere der folgenden ADME-Eigenschaften berechnet werden:
absorbierter Anteil einer Dosis der Substanz nach oraler Applikation,
- Konzentrations-Zeit-Kurvenverlauf in der Portalvene nach oraler
Applikation der Substanz,
freie Fraktion im Plasma,
- Organ - Plasma -Verteilungskoeffizienten, Konzentrations-Zeit-Kurven im Blutplasma und in den Organen nach oraler oder intravenöser Applikation der Substanz.
7. Verfahren nach einem der vorhergehenden Ansprüche 1 bis 6, wobei es sich bei der Substanz um eine Substanz für eine Pflanzenschutzanwendung handelt, und eine oder mehrere der folgenden ADME-Eigenschaften bereclmet werden:
a. Kennzahl für die Geschwindigkeit der Aufnahme in ein Blatt nach . einer Spritzapplikation,
b. Kennzahl für die Geschwindigkeit der Aufhahme in ein Insekt über die Kutikula,
c. Kennzahl für die Geschwindigkeit der Verteilung in der Pflanze nach
Blattapplikation (Phloemmobilität),
d. Kennzahl für die Geschwindigkeit der Verteilung in der Pflanze nach Wurzelapplikation (Xylemmobilität),
e. Konzentrationsverteilung zwischen den Organen eines Insekts.
8. Verfahren nach einem der vorhergehenden Ansprüche 1 bis 7 mit folgenden weiteren Schritten:
a. wiederholte Berechnung der ADME-Eigenschaft verschiedener Substanzen, insbesondere von verschiedenen Substanzen desselben chemischen Grundkörpers,
b. Sortierung der Substanzen nach dem Wert der ADME-Eigenschaft.
9. Verfahren nach Anspruch 8, wobei jeweils mehrere ADME-Eigenschaften für jede der Substanzen berechnet werden und die Sortierung der Substanzen durch die Bildung eines Index aus zwei oder mehreren der ADME-Eigenschaften erfolgt.
10. Verfahren nach Anspruch 9, wobei der Index durch eine gewichtete Summe der zumindest zwei ADME-Eigenschaften gebildet wird.
11. Verfahren nach Anspruch 8, 9 oder 10, wobei die Ausgabe der ADME- Eigenschaften der Substanzen in tabellarischer Form, vorzugsweise in einer relationalen Tabelle, erfolgt.
12. Computerprogramm zur Durchführung eines Verfahren nach einem der vorhergehenden Ansprüche 1 bis 11.
13. Computersystem zur Durchführung eines Verfahrens nach einem der vorhergehenden Ansprüche 1 bis 11.
14. Computersystem nach Anspruch 13 mit
a. einer ersten Datenbank (13) zur Eingabe der chemischen Struktur einer Substanz,
b. Mitteln (14) zum zyklischen Abfragen der Datenbank auf eine Eingabe einer chemischen Struktur,
c. Mitteln (15) zur automatischen Berechnung von molekularen Eigenschaften der Substanz basierend auf der chemischen Struktur, nachdem durch das zyklische Abfragen der Datenbank die Eingabe der chemischen Struktur erkannt worden ist, d. Mitteln (17) zur automatischen Berechnung einer oder mehrerer der ADME-Eigenschaften.
15. Computersystem nach Anspruch 13 oder 14 mit einer zweiten Datenbank (18) zur Speicherung von weiteren experimentell ermittelten molekularen
Eigenschaften der Substanz, wobei die automatische Berechnung der ADME- Eigenschaft basierend auf den berechneten und/oder den weiteren experimentell ermittelten molekularen Eigenschaften erfolgt.
16. Computersystem nach einem der vorhergehenden Ansprüche 13, 14 oder 15 mit Mitteln (21) zur Ausgabe der berechneten ADME-Eigenschaften verschiedener Substanzen in einer strukturierten Form, vorzugsweise in Tabellenform, insbesondere in Form einer relttionalen Tabelle.
17. Computersystem nach einem der vorhergehenden Ansprüche 13 bis 16 mit
Mitteln zur Berechnung eines Index aus den ADME-Eigenschaften zur Sortierung der Substanzen.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE10160270A DE10160270A1 (de) | 2001-12-07 | 2001-12-07 | Computersystem und Verfahren zur Berechnung von ADME-Eigenschaften |
| DE10160270 | 2001-12-07 | ||
| PCT/EP2002/014150 WO2003048720A2 (de) | 2001-12-07 | 2002-11-25 | Computersystem und verfahren zur berechnung von adme-eigenschaften |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| EP1456646A2 true EP1456646A2 (de) | 2004-09-15 |
Family
ID=7708460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| EP02783091A Withdrawn EP1456646A2 (de) | 2001-12-07 | 2002-11-25 | Computersystem und verfahren zur berechnung von adme-eigenschaften |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US7765092B2 (de) |
| EP (1) | EP1456646A2 (de) |
| JP (1) | JP2005512176A (de) |
| CN (1) | CN1659435A (de) |
| AU (1) | AU2002346846A1 (de) |
| CA (1) | CA2469388A1 (de) |
| DE (1) | DE10160270A1 (de) |
| RU (1) | RU2004120787A (de) |
| WO (1) | WO2003048720A2 (de) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10256315A1 (de) * | 2002-12-03 | 2004-06-17 | Bayer Ag | Computersystem und Verfahren zur Berechnung eines pharmakokinetischen Verhaltens einer chemischen Substanz in Insekten |
| DE10350525A1 (de) * | 2003-10-29 | 2005-06-09 | Bayer Technology Services Gmbh | Verfahren zur Visualisierung der ADME-Eigenschaften chemischer Substanzen |
| DE102004010516A1 (de) | 2004-03-04 | 2005-09-22 | Bayer Technology Services Gmbh | Verbessertes Verfahren zur zeitlichen Dosierung von Arzneistoffen |
| WO2006065950A2 (en) * | 2004-12-16 | 2006-06-22 | Pharmix Corporation | Modeling biological effects of molecules using molecular property models |
| DE102005028080A1 (de) | 2005-06-17 | 2006-12-21 | Bayer Technology Services Gmbh | Verfahren zur zeitlich gesteuerten intravenösen Verabreichung des Narkosemittels Propofol |
| US8380730B2 (en) * | 2008-10-09 | 2013-02-19 | International Business Machines Corporation | Program invocation from a query interface to parallel computing system |
| DE102010060311B3 (de) * | 2010-11-02 | 2012-02-02 | Freie Universität Berlin | Verfahren zur Unterstützung klinischer Studien |
| EP3236956A1 (de) | 2014-12-23 | 2017-11-01 | Bosteels, Arnaud | Kombination von remifentanil und propofol |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5751605A (en) * | 1996-08-15 | 1998-05-12 | Tripos, Inc. | Molecular hologram QSAR |
| US20020035459A1 (en) * | 1998-09-14 | 2002-03-21 | George M. Grass | Pharmacokinetic-based drug design tool and method |
| ATE295541T1 (de) * | 1999-06-18 | 2005-05-15 | Biacore Ab | Verfahren und vorrichtung zur untersuchung von wirksstoffskandidaten und zur bestimmung ihrer pharmakokinetischen parametern |
| EP1167969A2 (de) * | 2000-06-14 | 2002-01-02 | Pfizer Inc. | Verfahren und System zur Voraussage von pharmakokinetischen Eigenschaften |
| CA2416787A1 (en) * | 2000-07-28 | 2002-02-07 | Lion Bioscience Ag | System and method for predicting adme/tox characteristics of a compound |
| WO2002010741A2 (en) * | 2000-07-28 | 2002-02-07 | Lion Bioscience Ag | Regional intestinal permeability model |
| US20020169561A1 (en) * | 2001-01-26 | 2002-11-14 | Benight Albert S. | Modular computational models for predicting the pharmaceutical properties of chemical compunds |
-
2001
- 2001-12-07 DE DE10160270A patent/DE10160270A1/de not_active Ceased
-
2002
- 2002-11-25 CN CN028278100A patent/CN1659435A/zh active Pending
- 2002-11-25 EP EP02783091A patent/EP1456646A2/de not_active Withdrawn
- 2002-11-25 RU RU2004120787/15A patent/RU2004120787A/ru not_active Application Discontinuation
- 2002-11-25 AU AU2002346846A patent/AU2002346846A1/en not_active Abandoned
- 2002-11-25 CA CA002469388A patent/CA2469388A1/en not_active Abandoned
- 2002-11-25 WO PCT/EP2002/014150 patent/WO2003048720A2/de not_active Application Discontinuation
- 2002-11-25 US US10/497,759 patent/US7765092B2/en not_active Expired - Fee Related
- 2002-11-25 JP JP2003549866A patent/JP2005512176A/ja active Pending
Non-Patent Citations (1)
| Title |
|---|
| See references of WO03048720A2 * |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005512176A (ja) | 2005-04-28 |
| US20050119832A1 (en) | 2005-06-02 |
| US7765092B2 (en) | 2010-07-27 |
| CN1659435A (zh) | 2005-08-24 |
| DE10160270A1 (de) | 2003-06-26 |
| WO2003048720A2 (de) | 2003-06-12 |
| RU2004120787A (ru) | 2005-04-10 |
| WO2003048720A3 (de) | 2004-04-22 |
| CA2469388A1 (en) | 2003-06-12 |
| AU2002346846A1 (en) | 2003-06-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1671251B1 (de) | Verfahren zur simulation der wechselwirkung von chemischen substanzen mit lebenden organismen | |
| DE60037116T2 (de) | System zur behandlung eines patienten mit einem medikament zur optimierung der therapie und zur verhinderung eines negativen nebeneffektes | |
| WO2003048720A2 (de) | Computersystem und verfahren zur berechnung von adme-eigenschaften | |
| DE602004010154T2 (de) | Verfahren zur bestimmung der auswirkung eines aus mehreren komponenten bestehenden naturproduktgemischs auf das biologische profil einer krankheit in einer gruppe lebender systeme sowie entwicklung und qualitätskontrolle von medizin auf naturproduktbasis | |
| WO2021165386A1 (de) | Vorhersage von formulierungseigenschaften | |
| DE102018207557A1 (de) | Bestimmung eines Dehnungsgrades von Haar mit einem NIR-Sensor | |
| EP1426763B1 (de) | Computersystem und Verfahren zur Berechnung eines pharmakokinetischen Verhaltens einer chemischen Substanz in Insekten | |
| DE10350525A1 (de) | Verfahren zur Visualisierung der ADME-Eigenschaften chemischer Substanzen | |
| DE102018207558A1 (de) | Ermitteln von äußeren Schäden von Haaren | |
| LU102506A1 (de) | Ein Verfahren, ein System, ein Speichermedium und Programm zum Detektieren der Aromastilmerkmale des rauchgehärteten Tabaks | |
| Flache et al. | Multi-Agenten-Modelle | |
| DE2919517C2 (de) | ||
| EP1451750B1 (de) | Verfahren zur identifikation von pharmakophoren | |
| Holländer et al. | Current reference values of organ weights and indices for forensic practice, part 1: Brain and heart | |
| DE19781642B4 (de) | Verfahren zum Bestimmen eines die maximale Sauerstoffaufnahme beschreibenden Wertes eines zu beurteilenden Lebenswesens | |
| Wilde et al. | Pilot study on quality assurance in cardiac surgery | |
| Utschig | Analysis of Single Tree Productivity and Growing Space on Long Term Trial Plots–Methods, Program Development and Initial Results– | |
| DE10233022A1 (de) | Verfahren und Anlage zur Lösung von Aufgaben der adaptiven Chemie | |
| EP0598741B1 (de) | Verfahren zum generieren eines referenzstrings | |
| WO1994007197A1 (de) | Verfahren zur bearbeitung eines anwenderprogramms auf einem parallelrechnersystem | |
| WO2021255289A1 (de) | Verfahren zur bestimmung individueller pflegemittelformulierungen | |
| Schmidt et al. | Chronic pain patient’s acceptance and satisfaction with their analgesic medication: Development and validation of the QAMPAS questionnaire | |
| DE102005029437B4 (de) | Verfahren und Vorrichtung zum computergestützten Auffinden von ähnlichen Molekülen | |
| Meister | Evaluation verschiedener Imputationsverfahren zur Aufbereitung großer Datenbestände am Beispiel der SrV-Studie von 2013 | |
| Gyimesi et al. | HASH (0x556da14d4e68) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
| AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SK TR |
|
| AX | Request for extension of the european patent |
Extension state: AL LT LV MK RO SI |
|
| 17P | Request for examination filed |
Effective date: 20041022 |
|
| 17Q | First examination report despatched |
Effective date: 20041223 |
|
| STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN |
|
| 18D | Application deemed to be withdrawn |
Effective date: 20060202 |