{UTF-8}
********************************************************
	.thalidomid
********************************************************~
		Simulation: Peter Bützer
	|

K=
	kirs/kisr
	~	Dmnl [0,?]
	~	Gleichgewichtskonstante carakterisiert nur die Lage, nicht aber die \
		Geschwindigkeit der Gleichgewichtseinstellung
	|

Br=
	krr*R1*(1-Bv)
	~	mg/Hour [0,?]
	~	Biologisch nicht verfügbar, geht durch den Verdauungstrakt ohne \
		aufgenommen zu werden
	|

Bs=
	ksr*S1*(1-Bv)
	~	mg/Hour [0,?]
	~	Biologisch nicht verfügbar, geht durch den Verdauungstrakt ohne \
		aufgenommen zu werden
	|

Bv=
	0.4
	~	Dmnl [0,0.6]
	~	Biologische Verfügbarkeit, experimentell ermittelt
	|

HWZrs=
	2.5
	~	Hour [2.5,4.1]
	~	krs=ln(0.25)/5=0.2773, HWZrs=ln(2)/0.2773=2.5; andere Angaben HWZ = 4.1 h, \
		Diese Halbwertszeit wird durch das keto-enol-Gleichgewicht und die \
		enzymatische Isomerisierung bestimmt
	|

Rr=
	krr*R1*Bv
	~	mg/Hour [0,?]
	~	Die Resorption ist abhängig von der vorhandenen Menge (1.Ordnung), ergibt \
		zusammen mit der Elimination die Bateman-Funktion
	|

S1= INTEG (
	-Bs-Sr,
		Dosis*"RS-S"*0.4)
	~	mg [0,?]
	~	Menge S-Isomer im Verdauungstrakt; Biologische Verfügbarkeit: 40% --> 0.4
	|

Sr=
	ksr*S1*Bv
	~	mg/Hour [0,?]
	~	Die Resorption ist abhängig von der vorhandenen Menge (1.Ordnung), ergibt \
		zusammen mit der Elimination die Bateman-Funktion
	|

R1= INTEG (
	-Br-Rr,
		Dosis*(1-"RS-S")*0.4)
	~	mg [0,?]
	~	Menge R-Isomer  im Verdauungstrakt; Biologische Verfügbarkeit: 40% --> 0.4
	|

Dosis=
	100
	~	mg [50,400,50]
	~		|

R2= INTEG (
	+Rr+RSi-Re+SRi,
		0)
	~	mg [0,?]
	~	Menge R-Isomer  im Blut
	|

S2= INTEG (
	+Sr-RSi-Se-SRi,
		0)
	~	mg [0,?]
	~	Menge S-Isomer  im Blut
	|

kisr=
	LN(2)/HWZsr
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Inversion S-> R
	|

HWZsr=
	3.75
	~	Hour [3.75,5.8]
	~	ksr=ln(0.25)/7.5=0.1848, HWZsr=ln(2)/0.1848=3.75; andere Angaben: HWZ= 5.8 \
		h
	|

SRi=
	kisr*S2
	~	mg/Hour [0,?]
	~	Die Isomerisierungsgeschwindigkeit ist abhängig von der vorhandenen Menge \
		(1.Ordung)
	|

krr=
	LN(2)/HWZr
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Resorption
	|

kse=
	LN(2)/HWZe
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Elimination des S-Isomers
	|

HWZr=
	0.75
	~	Hour [0.5,2]
	~	Der maximale Plasmaspiegel wird in 2.9 - 5.7 Stunden erreicht (Mittelwert: \
		4h); daraus ergibt sich für die Resorption eine  HWZr=0.75 h
	|

ksr=
	LN(2)/HWZr
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Resorption
	|

HWZe=
	12
	~	Hour [8,15]
	~	Experimentelle Angabe, ist ein Mittelwert
	|

kre=
	LN(2)/HWZe
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Elimination des R-Isomers
	|

kirs=
	LN(2)/HWZrs
	~	1/Hour [0,?]
	~	Reaktionsgeschwindigkeits-Konstante der Inversion R-> S
	|

"RS-S"=
	1
	~	Dmnl [0.5,1]
	~	Anteil S-Isomer: Reines S-Isomer=1, reines Racemat 0.5, reines R-Isomer=0
	|

RSi=
	-kirs*R2
	~	mg/Hour [0,?]
	~	Die Isomerisierungsgeschwindigkeit ist abhängig von der vorhandenen Menge \
		(1.Ordung)
	|

Re=
	kre*R2
	~	mg/Hour [0,?]
	~	Die Eliminationsgeschwindigkeit ist abhängig von der vorhandenen Menge \
		(1.Ordnung), ergibt zusammen mit der Resorption die Bateman-Funktion
	|

Se=
	kse*S2
	~	mg/Hour [0,?]
	~	Die Eliminationsgeschwindigkeit ist abhängig von der vorhandenen Menge \
		(1.Ordnung), ergibt zusammen mit der Resorption die Bateman-Funktion
	|

********************************************************
	.Control
********************************************************~
		Simulation Control Parameters
	|

FINAL TIME  = 72
	~	Hour
	~	The final time for the simulation.
	|

INITIAL TIME  = 0
	~	Hour
	~	The initial time for the simulation.
	|

SAVEPER  = 
        TIME STEP
	~	Hour [0,?]
	~	The frequency with which output is stored.
	|

TIME STEP  = 0.1
	~	Hour [0,?]
	~	The time step for the simulation.
	|

\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*View 1
$192-192-192,0,Times New Roman|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|96,96,100,0
10,1,R1,313,117,32,16,3,3,0,4,0,0,0,0,0-0-0,255-192-128,|12||0-0-0
10,2,S1,314,217,32,16,3,3,0,4,0,0,0,0,0-0-0,192-255-255,|12||0-0-0
10,3,R2,470,117,32,16,3,3,0,4,0,0,0,0,0-0-0,255-192-128,|12||0-0-0
10,4,S2,473,220,32,16,3,3,0,4,0,0,0,0,0-0-0,192-255-255,|12||0-0-0
1,5,7,3,36,0,0,22,2,0,0,-1--1--1,|12||0-0-0,1|(449,146)|
1,6,7,4,68,0,0,22,2,0,0,-1--1--1,|12||0-0-0,1|(449,187)|
11,7,444,449,165,6,5,33,3,0,0,4,0,0,0
10,8,RSi,471,165,16,9,40,131,0,0,-1,0,0,0
1,9,11,3,4,0,0,22,0,0,0,-1--1--1,,1|(417,117)|
1,10,11,1,100,0,0,22,0,0,0,-1--1--1,,1|(365,117)|
11,11,396,391,117,5,6,34,3,0,0,1,0,0,0
10,12,Rr,391,133,15,10,40,3,0,0,-1,0,0,0
1,13,15,4,4,0,0,22,0,0,0,-1--1--1,,1|(419,217)|
1,14,15,2,100,0,0,22,0,0,0,-1--1--1,,1|(366,217)|
11,15,316,392,217,5,6,34,3,0,0,1,0,0,0
10,16,Sr,392,233,14,10,40,3,0,0,-1,0,0,0
12,17,48,630,117,10,8,0,3,0,0,-1,0,0,0
1,18,20,17,4,0,0,22,0,0,0,-1--1--1,,1|(589,117)|
1,19,20,3,100,0,0,22,0,0,0,-1--1--1,,1|(525,117)|
11,20,48,553,117,5,6,34,3,0,0,1,0,0,0
10,21,Re,553,133,10,10,40,3,0,0,-1,0,0,0
12,22,48,631,218,10,8,0,3,0,0,-1,0,0,0
1,23,25,22,4,0,0,22,0,0,0,-1--1--1,,1|(589,218)|
1,24,25,4,100,0,0,22,0,0,0,-1--1--1,,1|(526,218)|
11,25,48,552,218,5,6,34,3,0,0,1,0,0,0
10,26,Se,552,234,10,10,40,3,0,0,-1,0,0,0
1,27,3,20,1,0,0,0,0,64,0,-1--1--1,,1|(521,82)|
1,28,4,25,1,0,0,0,0,64,0,-1--1--1,,1|(521,194)|
1,29,3,7,1,0,0,0,0,64,0,-1--1--1,,1|(418,130)|
1,30,1,11,1,0,0,0,0,64,0,-1--1--1,,1|(364,82)|
1,31,2,15,1,0,0,0,0,64,0,-1--1--1,,1|(362,193)|
10,32,kirs,393,164,16,8,8,131,0,0,0,0,0,0
10,33,krr,340,55,14,10,8,3,0,0,0,0,0,0
10,34,ksr,332,262,26,10,8,131,0,0,0,0,0,0
10,35,kse,491,269,13,10,8,3,0,0,0,0,0,0
10,36,kre,514,53,12,10,8,3,0,0,0,0,0,0
1,37,36,20,1,0,0,0,0,64,0,-1--1--1,,1|(553,78)|
1,38,35,26,1,0,0,0,0,64,0,-1--1--1,,1|(528,268)|
1,39,34,16,1,0,0,0,0,64,0,-1--1--1,,1|(378,265)|
1,40,33,11,1,0,0,0,0,64,0,-1--1--1,,1|(389,80)|
1,41,32,7,1,0,0,0,0,64,0,-1--1--1,,1|(414,149)|
10,42,HWZrs,321,166,35,14,8,131,0,0,0,0,0,0
1,43,42,32,1,0,0,0,0,64,0,-1--1--1,,1|(356,148)|
10,44,"RS-S",206,183,18,10,8,3,0,0,0,0,0,0
10,45,HWZe,453,70,26,10,8,3,0,0,0,0,0,0
1,46,45,36,1,0,0,0,0,64,0,-1--1--1,,1|(480,54)|
1,47,45,35,1,0,0,0,0,64,0,-1--1--1,,1|(447,230)|
10,48,HWZr,122,166,22,10,8,3,0,0,0,0,0,0
1,49,48,33,1,0,0,0,0,64,0,-1--1--1,,1|(188,73)|
1,50,48,34,1,0,0,0,0,64,0,-1--1--1,,1|(219,266)|
10,51,HWZsr,612,178,25,10,8,3,0,0,0,0,0,0
10,52,kisr,561,183,14,10,8,3,0,0,0,0,0,0
1,53,51,52,1,0,0,0,0,64,0,-1--1--1,,1|(591,157)|
1,54,56,3,4,0,0,22,0,0,0,-1--1--1,,1|(494,148)|
1,55,56,4,100,0,0,22,0,0,0,-1--1--1,,1|(494,188)|
11,56,476,494,168,6,5,33,3,0,0,4,0,0,0
10,57,SRi,514,168,14,10,40,3,0,0,-1,0,0,0
1,58,52,56,1,0,0,0,0,64,0,-1--1--1,,1|(534,156)|
1,59,4,56,1,0,0,0,0,64,0,-1--1--1,,1|(530,187)|
10,60,Dosis,214,86,20,10,8,3,0,0,0,0,0,0
12,61,0,736,506,266,163,3,188,0,0,1,0,0,0
Vergleich
12,62,48,206,118,10,8,0,3,0,0,-1,0,0,0
1,63,65,62,4,0,0,22,0,0,0,-1--1--1,,1|(229,118)|
1,64,65,1,100,0,0,22,0,0,0,-1--1--1,,1|(267,118)|
11,65,48,248,118,6,8,34,3,0,0,1,0,0,0
10,66,Br,248,137,10,11,40,3,0,0,-1,0,0,0
12,67,48,206,217,10,8,0,3,0,0,-1,0,0,0
1,68,70,67,4,0,0,22,0,0,0,-1--1--1,,1|(229,217)|
1,69,70,2,100,0,0,22,0,0,0,-1--1--1,,1|(268,217)|
11,70,48,249,217,6,8,34,3,0,0,1,0,0,0
10,71,Bs,249,236,10,11,40,3,0,0,-1,0,0,0
1,72,1,65,1,0,0,0,0,64,0,-1--1--1,,1|(281,149)|
1,73,2,70,1,0,0,0,0,64,0,-1--1--1,,1|(275,188)|
10,74,Bv,173,142,11,11,8,3,0,0,0,0,0,0
1,75,74,65,1,0,0,0,0,64,0,-1--1--1,,1|(212,143)|
1,76,74,11,1,0,0,0,0,64,0,-1--1--1,,1|(272,152)|
1,77,74,15,1,0,0,0,0,64,0,-1--1--1,,1|(317,183)|
1,78,74,70,1,0,0,0,0,64,0,-1--1--1,,1|(237,178)|
1,79,33,65,1,0,0,0,0,64,0,-1--1--1,,1|(275,78)|
1,80,34,70,1,0,0,0,0,64,0,-1--1--1,,1|(272,255)|
1,81,60,1,0,0,0,0,0,64,1,-1--1--1,,1|(250,97)|
1,82,44,1,0,0,0,0,0,64,1,-1--1--1,,1|(248,156)|
1,83,60,2,0,0,0,0,0,64,1,-1--1--1,,1|(256,142)|
1,84,44,2,0,0,0,0,0,64,1,-1--1--1,,1|(246,194)|
10,85,K,593,140,8,11,8,3,0,0,0,0,0,0
1,86,32,85,1,0,0,0,0,64,0,-1--1--1,,1|(472,139)|
1,87,52,85,1,0,0,0,0,64,0,-1--1--1,,1|(559,148)|
///---\\\
:GRAPH GRAPH
:TITLE Vergleich
:X-AXIS Time
:SCALE
:VAR S2|S-Form
:Y-MIN 0
:Y-MAX 0.5
:DATASET S-Form
:SCALE
:VAR R2|Racemat
:Y-MIN 0
:Y-MAX 0.5
:DATASET Racemat

:GRAPH Vergleich
:TITLE Vergleich
:X-AXIS Time
:X-LABEL Stunden
:SCALE
:VAR R2|R-Enantiomer
:UNITS mg
:LINE-WIDTH 2
:VAR S2|S-Enantiomer
:UNITS mg
:LINE-WIDTH 2
:L<%^E!@
1:Current.vdf
9:Current
23:0
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
4:Time
5:Br
24:0
25:71.9995
26:71.9995