Transcript κινητικη αναπτυξης μικροβιων και παραγωγης μεταβολικων προϊοντων

ΚΙΝΗΤΙΚΗ ΑΝΑΠΤΥΞΗΣ
ΜΙΚΡΟΒΙΩΝ ΚΑΙ ΠΑΡΑΓΩΓΗΣ
ΜΕΤΑΒΟΛΙΚΩΝ ΠΡΟΪΟΝΤΩΝ
Μικροβιακή ανάπτυξη
Μέτρηση βιομάζας με
φασματοφωτομετρία
πηγή φωτός
φωτοηλεκτρικό κύτταρο
οπτική
πυκνότητα
βιομάζα (mg)
Τριβλεία Petri
Ανάπτυξη καλλιέργειας
ln (αρ/ml)
IV
III
II
I
Χρόνος
V
VI
ΓΕΝΙΚΕΣ ΕΝΝΟΙΕΣ ΚΙΝΗΤΙΚΗΣ
aA+ bB + cC+... pP + qQ+ rR+...
• Η αντίδραση μπορεί να αντιπροσωπεύει μία
συνολική αντίδραση στην οποία περίπτωση,
όπως είδαμε, οι στοιχειομετρικοί συντελεστές
δεν είναι σταθεροί.
• Όταν δεν μεταβάλλονται αισθητά, μπορούμε να
τους θεωρήσουμε σταθερούς και να χειριστούμε
την αντίδραση σαν να ήταν απλή.
Ρυθμός αντίδρασης
1 d d 
r=
=
V dt
dt
*
• V ο όγκος του χώρου που λαμβάνει χώρα η
αντίδραση
• ξ η έκταση της αντίδρασης
• ξ*=ξ/V η ειδική έκταση
Η τελευταία ισότητα ισχύει μόνο για συστήματα
σταθερού όγκου.
Για κάθε ουσία Ni
1 dN i
dci
= ri =  i r =
V dt
dt
Ο
•
•
•
ρυθμός εν γένει εξαρτάται από:
την θερμοκρασία,
την πίεση
την σύσταση
Εκτός από τις ουσίες που συμπεριλαμβάνονται στην αντίδραση,
υπάρχουν και άλλες ουσίες οι οποίες αν και παρούσες ως
ιχνοστοιχεία μπορούν να επηρεάζουν σημαντικά τον ρυθμό
αντίδρασης.
Μακριά από την ισορροπία:
r = k (ci)
• φ(ci) συνάρτηση των συγκεντρώσεων ci
• και k μία "σταθερά" που ονομάζεται κινητική σταθερά.
• Η k είναι σταθερά μόνο όσο αφορά τις συγκεντρώσεις
των αντιδρώντων και προϊόντων αλλά εξαρτάται από την
θερμοκρασία με την σχέση Arrhenius:
k = Ae
-E
RT
•R η σταθερά των αερίων
•Τ η απόλυτη θερμοκρασία
•Α μία σταθερά που ονομάζεται προεκθετικός συντελεστής
Αν ο ρυθμός μπορεί να γραφεί στην
μορφή:
A
B
r = k c A cB ...
• τα βi είναι οι τάξεις της αντίδρασης ως προς τις
ουσίες i
• Τέτοιου είδους εκφράσεις παίρνουμε για τις απλές
αντιδράσεις  κινητική δρώντων μαζών
• Τότε η τάξη ως προς μία ουσία συμπίπτει με τον
αριθμό των μορίων που παίρνουν μέρος στην
αντίδραση (μοριακότητα).
• Δηλαδή για κινητική δρώντων μαζών: βΑ=a, βB=b…
• Όταν δεν παρατηρείται κινητική δρώντων μαζών, τότε
μπορούμε να συμπεράνουμε ότι πρόκειται για συνολική
αντίδραση
• Το σύνολο των απλών αντιδράσεων που οδηγούν σε κάποια
συνολική αντίδραση ονομάζεται μηχανισμός της συνολικής
αντίδρασης.
• Η ύπαρξη πολλαπλών μηχανισμών συμβατών με την συνολική
αντίδραση δυσχεραίνει το πρόβλημα της διαλεύκανσης του
σωστού μηχανισμού.
• Δεν υπάρχει γενική μέθοδος μέτρησης του ρυθμού
αντίδρασης. Αν η στοιχειομετρία είναι γνωστή (και σταθερή)
αρκεί η παρακολούθηση μίας μόνο ουσίας ,μια και για κάθε i:
1 dci
r=
 i dt
Μέθοδοι επεξεργασίας κινητικών
δεδομένων
Διακρίνουμε τρεις γενικές κατηγορίες:
• διαφορικές μεθόδους
• ολοκληρωτικές μεθόδους
• μεθόδους που βασίζονται στην απλοποίηση
της παράστασης του ρυθμού
(χρησιμοποιώντας π.χ. περίσσεια κάποιων
αντιδρώντων)
Δύο Διαφορικές μέθοδοι
• Μέθοδος Α: παραγώγιση των δεδομένων από
μία και μοναδική πειραματική καμπύλη.
• Μέθοδος B: Με βάση πολλαπλές μετρήσεις του
αρχικού ρυθμού αντίδρασης αλλάζοντας κάθε
φορά την αρχική συγκέντρωση κάποιων
αντιδρώντων.
Μέθοδος Α
• με δεδομένα τα ci συναρτήσει του t
• υποθέτουμε τη μορφή της φ(ci)
• παραγωγίζοντας γραφικά το διάγραμμα της
συγκέντρωσης συναρτήσει του χρόνου γίνεται
προσδιορισμός του ρυθμού συναρτήσει του
χρόνου.
• υπολογίζοντας την παράσταση φ(ci) συναρτήσει
του χρόνου και κατασκευάζοντας το διάγραμμα
του ρυθμού r ως προς φ(ci) βλέπουμε αν τα
σημεία κείνται σε ευθεία.
η υπόθεση για την μορφή
εξάρτησης ήταν σωστή και η κλίση
της ευθείας μας δίνει την κινητική
σταθερά.
υποθέτουμε διαφορετική μορφή
της φ(ci) και ξαναδοκιμάζουμε
μέχρι να φέρουμε τα δεδομένα
σε ευθεία γραμμή.
Μέθοδος πολλαπλών αρχικών
συγκεντρώσεων
ro = k (cio)
Ολοκληρωτικές μέθοδοι
• Ολοκληρώνουμε μία
έκφραση για τον ρυθμό

που έχουμε αρχικά
1 d
d
= k ( ci )  kt = 

υποθέσει για να
V( ci )
V dt
0
προβλέψουμε τις
συγκεντρώσεις
συναρτήσει του χρόνου
• Εάν το διάγραμμα του
ψ(ci) συναρτήσει του t
δώσει ευθεία γραμμή, η
έκφραση που υποτέθηκε
είναι σωστή, και η κλίση
μάς δίνει την κινητική
σταθερά.
• Εάν όχι τότε υποθέτουμε
νέα έκφραση και να
επαναλαμβάνουμε τη
διαδικασία.
 (ci )
Οι ρυθμοί
(α) ο ρυθμός ανάπτυξης βιομάζας rX=dx/dt,
(β) o ρυθμός κατανάλωσης υποστρώματος rS=-d[S]/dt,
(γ) o ρυθμός σχηματισμού προϊόντος rP=d[P]/dt, και
(δ) o ρυθμός παραγωγής θερμότητας rH=dH/dt.
Ειδικοί ρυθμοί (ανά g κυτταρικής μάζας)
1 dx
=
x dt
1 d[S]
qS = x dt
1 d[P]
qP =
x dt
1 dH
qH =
x dt
Ταξινόμηση μοντέλων μικροβιακής ανάπτυξης
ΜΗ ΔΟΜΗΜΕΝΑ
ΔΟΜΗΜΕΝΑ
ΜΗ
Ο κυτταρικός πληθυσμός είναι
Διακρίνουμε πολλά
ΚΑΤΑΝΕΜΗΜΕΝΑ
μία μεταβλητή ενός και μόνο
συστατικά που απαρτίζουν
συστατικού
την κυτταρική μάζα
Όλα τα κύτταρα έχουν την ίδια
Όλα τα κύτταρα έχουν την
ηλικία.
ίδια ηλικία.
Θεωρούμε την κυτταρική μάζα
Διακρίνουμε πολλά
ως μία μεταβλητή
συστατικά που απαρτίζουν
ΚΑΤΑΝΕΜΗΜΕΝΑ
Διακρίνουμε κατανομή ηλικίας
κυττάρων
την κυτταρική μάζα
Διακρίνουμε κατανομή
ηλικίας κυττάρων
Εξισορροπημένη ανάπτυξη
«Μέσο»
ΜΗ ΔΟΜΗΜΕΝΑ
κύτταρο
ΜΗ ΚΑΤΑΝΕΜΗΜΕΝΑ


Ο κυτταρικός
ΔΟΜΗΜΕΝΑ

πληθυσμός μία
συστατικά που
μεταβλητή ενός και
απαρτίζουν την
μόνο συστατικού
κυτταρική μάζα
όλα τα κύτταρα έχουν

την ίδια ηλικία.
ΚΑΤΑΝΕΜΗΜΕΝΑ


Διακρίνουμε πολλά
Θεωρούμε την
όλα τα κύτταρα έχουν
την ίδια ηλικία.

Διακρίνουμε πολλά
κυτταρική μάζα ως μία
συστατικά που
μεταβλητή
απαρτίζουν την
κυτταρική μάζα
Διακρίνουμε κατανομή
ηλικίας κυττάρων

Διακρίνουμε κατανομή
ηλικίας κυττάρων
Καμπύλη ανάπτυξης
ln (αρ/ml)
IV
III
II
I
Χρόνος
V
VI
I. φάση
καθυστέρησης
II. φάση επιτάχυνσης
III.φάση εκθετικής
ανάπτυξης
IV. φάση
επιβράδυνσης
V. στάσιμη φάση
VI. φάση θανάτου ή
αποδόμησης
ΜΗ ΔΟΜΗΜΕΝΑ - ΜΗ
ΚΑΤΑΝΕΜΗΜΕΝΑ ΜΟΝΤΕΛΑ
ΜΟΝΤΕΛΑ ΑΝΑΠΤΥΞΗΣ ΜΙΑΣ
ΜΕΤΑΒΛΗΤΗΣ
dx
= f(x)
rX =
dt
Malthus
dx
= x
f(x) =  x 
dt
• ο ειδικός ρυθμός ανάπτυξης μ είναι
σταθερός.
• Ολοκληρώνοντας την εξίσωση με αρχική
συνθήκη x(0)=xo:
t
x(t)= xo e
Ανεξέλεγκτη ανάπτυξη!
• Ενώ μπορεί να περιγράψει την ανάπτυξη
στην εκθετική φάση, δεν μπορεί να
προβλέψει άλλες φάσεις.
• το μοντέλο αυτό είναι ανίκανο να
προβλέψει την επήρεια στον ρυθμό
ανάπτυξης οιασδήποτε παραμέτρου όπως
συγκεντρώσεις θρεπτικών συστατικών
κ.λ.π.
• Από την εξίσωση μπορούμε να
προβλέψουμε τον χρόνο διπλασιασμού
της βιομάζας:
ln 2
tD =

Παράδειγμα: αποδόμηση
σουλφαμεθοξαζόλης
350
250
300
200
250
Biomass (mg/L)
Sulfamethoxazole (mg/L)
300
150
100
50
200
150
100
50
0
0
2
4
6
Time (days)
8
10
12
0
0
2
4
6
Time (days)
8
10
12
Malthus model
x
50
70
140
150
180
200
225
250
ln(x)
3,9
4,2
4,9
5,0
5,2
5,3
5,4
5,5
dX
 X
dt
μ=0.1894 d-1
7
6
5
ln(x)
time (day)
1
2
5
6
7
8
9
10
4
3
2
1
0
0
5
10
Time (d)
15
Λογιστικό ή σιγμοειδές
Verhulst 1844, Pearl και Reed το 1920
βx
Malthus
dx
= kx(1-  x)
dt
λογιστική
kt
xo e
x=
1-  xo (1- ekt )
• Προβλέπει εκθετική αλλά και στάσιμη
Λογιστικό
dx
= kx(1-  x)
dt
k=0,328 d-1
β=0,032 l/mg
Πρόβλεψη λογιστικού μοντέλου
Ανάπτυξη μέχρι 315 mg/l σε 22 περίπου d
Όμως στάσιμη φάση στις 10 d
Volterra
t
dx
= kx(1-  x)+ k o|  x(r)dr|
dt
0
• περιλαμβάνει έναν επί πλέον όρο "μνήμης"
αναγνωρίζοντας ότι οι παρελθούσες συνθήκες
επηρεάζουν τον τρέχοντα ρυθμό ανάπτυξης είτε
θετικά (ko>0) είτε αρνητικά (ko<0)
Τα μοντέλα μιας μεταβλητής:
• είναι όλα εμπειρικά στην φύση τους
• δεν προβλέπουν φάση καθυστέρησης
• δεν περιγράφουν τις επιδράσεις των
διαφόρων παραγόντων στον ρυθμό
ανάπτυξης
• δεν κάνουν καμία χρήση των γνώσεων της
βιοχημείας και μικροβιολογίας.
ΜΟΝΤΕΛΑ ΑΝΑΠΤΥΞΗΣ ΜΕ
ΠΕΡΙΟΡΙΣΜΟ ΥΠΟΣΤΡΩΜΑΤΟΣ
r X =  ([S],x)x
• Διάφορα μοντέλα ανάλογα με τη
συνάρτηση
Το μοντέλο Monod
μ max S
μ
KS  S
όπου:
S: περιοριστικό υπόστρωμα (mg/L)
Ks: σταθερά κορεσμού (ίση με την συγκέντρωση του
υποστρώματος στην οποία ο ρυθμός ανάπτυξης είναι
ίσος με το ήμισυ του μέγιστου ειδικού ρυθμού
ανάπτυξης)
μmax: μέγιστος ειδικός ρυθμός ανάπτυξης
300
350
250
300
Biomass (mg/L)
Sulfamethoxazole (mg/L)
Παράδειγμα σουλφομεθοξαζόλης
200
150
100
250
200
150
100
50
50
0
0
0
2
4
6
8
10
12
0
Time (days)
2
4
6
8
10
12
Time (days)
Kinetic model
Kinetic constants
Sulf   X
dX
  max
dt
K S  Sulf 
μmax

Sulf   X
d Sulf 
  max
dt
Y x / s K S  Sulf 
Yx/s
= 0.21 d-1
Ks = 21.15 mg/L
= 0.95 m g biom ass
m g sulfamethoxazole
Εναλλακτικές εκφράσεις
(α) Teissier:
[S]
 =  max (1 - e- K )
S
(β) Moser:
 max [S ]
n
=
n
K S + [S ]
(γ) Contois:
=
 max [S]
Bx + [S]
Εναλλακτικές εκφράσεις
(δ) Powell:
=
 max [S]
( K S + K D )+ [S]
(ε) Blackman (ή μοντέλο δύο φάσεων):

[S]
=
,[S] < 2K
 max 2K

 max
= 1,[S] > 2K
(στ) Dabes:
 - KS
[S] =   +
 max - 
Konak
d[

 max
d[S]
]
= k  max
p-1
(1 -

 max
p
)
• Το μοντέλο αυτό είναι υπό μορφή
κανονικής διαφορικής εξίσωσης
Γενικευμένη εξίσωση Kargi και Shuler
d  rel
m
p
= K(  rel ) (1-  rel )
d[S]
μrel=μ/μmax
ΜΟΝΤΕΛΟ
K
m
p
MONOD
1/KS
0
2
TEISSIER
1/KS
0
1
MOSER
n/(KS)(1/n)
1-1/n
1+1/n
CONTOIS
1/(Bx)
0
2
KONAK
k μmax
p
p
Παρεμπόδιση υποστρώματος
μmax=0,5 d-1
Κs=30 mg/l
KΙ=100 mg/l
• Andrews
=
 max
[S]
K
S
1+
+
[S] K I
Πολλαπλά υποστρώματα
• Διαδοχική χρήση (διαυξική ανάπτυξη εναλλακτικών
υποστρωμάτων)
• Ταυτόχρονη χρήση (συμπληρωματικών
υποστρωμάτων)
• Οι περιπτώσεις της ταυτόχρονης χρήσης
διακρίνονται περαιτέρω επί τη βάσει του αν ο
περιορισμός είναι ταυτόχρονος ή όχι.
Ταυτόχρονος περιορισμός
[ S1 ]
[ S2 ]
 =  max
K S1 +[ S 1 ] K S2 +[ S 2 ]
Μη ταυτόχρονος περιορισμός
 max1[ S 1 ] [ S 1 ] [ S 2 ]
=
,
<
K S1 + [ S 1 ] K S1
K S2
 max2[ S 2 ] [ S 1 ] [ S 2 ]
=
,
>
K S2 + [ S 2 ] K S1
K S2
[S2]/KS2
[S2]/KS2
μ/μmax
μ/μmax
0,7
0,5
[S1]/KS1
(α) ταυτόχρονος περιορισμός
0,7
0,5
[S1]/KS1
(β) μη ταυτόχρονος περιορισμός
Διαδοχική χρήση
 =  ([ S 1 ])+  ([ S 2 ]) f  .
(α) Imanaka:
f  =
[ S2 ]
[ S1 ] + [ S 2 ]
(β) Moser:
f  =
1
[ S1 ]
1+
K
(γ) Bergter and Knorre:
df 
1
KS
=
(
- f  )
dt
t L2 K S + [ S 1 ]
Επαύξηση
• Στην περίπτωση επαύξησης του ρυθμού
από την παρουσία κάποιας ένωσης SE
ακολουθώντας τους Tsao και Hanson
έχουμε την ακόλουθη τροποποίηση του
μέγιστου ειδικού ρυθμού ανάπτυξης:
 max =  maxo +
 maxe [ S E ]
K SE + [ S E ]
Στην περίπτωση της παρεμπόδισης από κάποια ουσία Ι, συνήθως
χρησιμοποιείται μία από τις ακόλουθες εκφράσεις:
 max =  maxo (1- K I [I])
 max =  maxo e-K [I]
I
 max =
 maxo
[I]
1+
KI
Συχνά η παρεμπόδιση προέρχεται από κάποιο μεταβολικό προϊόν P.
Συνήθως και για την παρεμπόδιση προϊόντος χρησιμοποιούνται οι
πιο πάνω εκφράσεις.
Μοντέλα καθυστέρησης
• O’Neil and
Lyberatos
• Wang and
Stephanopoulos
dz
=  ([S] - z)
dt
 max [S]
d
=(
- )
dt
K S + [S]
ΣΥΝΤΗΡΗΣΗ ΚΑΙ ΕΝΔΟΓΕΝΗΣ
ΜΕΤΑΒΟΛΙΣΜΟΣ
• κατανάλωση υποστρώματος όχι μόνο για
ανάπτυξη αλλά και για συντήρηση με ρυθμό
ανάλογο της συγκέντρωσης της ζώσας βιομάζας:
r m = mS xv
• όπου xv η συγκέντρωση της ζώσας βιομάζας.
• Η τιμή της σταθεράς ms εξαρτάται από τις
συνθήκες.
• όσο πιο αντίξοες οι συνθήκες, τόσο μεγαλύτερη
και η σταθερά
Ενδογενής μεταβολισμός ή
ενδογενής αναπνοή
• Η ενέργεια που απαιτείται για την διατήρηση ζώντων
οργανισμών και γίνεται εις βάρος της βιομάζας
• δεν επηρεάζεται το ισοζύγιο του υποστρώματος αλλά
από τον ρυθμό ανάπτυξης βιομάζας αφαιρούμε ένα
ρυθμό ενδογενούς μεταβολισμού:
r s = k e xv
Ταξινόμηση μεταβολικών
προϊόντων με βάση την κινητική
(α) Προϊόντα συσχετιζόμενα με την ανάπτυξη, με ρυθμό
παραγωγής ανάλογο του ρυθμού ανάπτυξης:
rP =  rX
(β) Προϊόντα μη συσχετιζόμενα με την ανάπτυξη, με ρυθμό
παραγωγής ανάλογο της συγκέντρωσης της βιομάζας:
r P =  xv
(γ) Προϊόντα μικτής κινητικής:
r P =  r X +  xv
ΜΟΝΤΕΛΑ ΓΙΑ ΤΗΝ ΧΡΗΣΗ
ΥΠΟΣΤΡΩΜΑΤΟΣ
n
r
r
X
Pi
=
+
+ mS x v
rS

Y X/S i=1 Y Pi/S
• όπου P1,...,Pn είναι τα εξωκυτταρικά προϊόντα, και YPi/S
είναι οι συντελεστές απόδοσης που συσχετίζουν το
υπόστρωμα με το προϊόν i.
• Υποθέτοντας αμελητέα χρήση για συντήρηση και ότι όλα
τα προϊόντα είναι συσχετιζόμενα με την ανάπτυξη:
rX
rS =
Y
Πότε ισχύει το κάτωθι;
rX =
 max [S]
K S + [S]
x
1  max [S]
x
rS =
Y K S + [S]
Όταν: 1. όλα τα προϊόντα είναι σχετιζόμενα με την ανάπτυξη
2. η συντήρηση και ο ενδογενής μεταβολισμός
μπορούν να αγνοηθούν
Φάση καθυστέρησης
Dean και Hinshelwood
*
c aV
a
a


tL=
a"
No+
a
• ο χρόνος καθυστέρησης tL είναι ο χρόνος μετά την μεταφορά
της μαγιάς στο νέο θρεπτικό μέσο στο οποίο κάποιο κρίσιμο
ένζυμο φτάνει σε κάποια κρίσιμη συγκέντρωση c*
–
–
–
–
V: ο λόγος των όγκων του προηγούμενου μέσου προς τον όγκο του νέου
Νο : ο αριθμός των παλαιών κυττάρων ανά μονάδα όγκου στο νέο μέσο
a: η συγκέντρωση του ενζύμου στο προηγούμενο μέσο
a’: ο ρυθμός παραγωγής του ενζύμου από τα "παλαιά" κύτταρα ανά
κύτταρο
– a" : ο ρυθμός παραγωγής του ενζύμου από νέα κύτταρα
Προκειμένου να ελαχιστοποιήσουμε τον χρόνο καθυστέρησης
χρειάζεται μεγάλη μαγιά (No) και μεγάλος λόγος όγκων (V).
Στην πράξη προσπαθούμε να χρησιμοποιούμε όγκο μαγιάς περίπου
το ένα δέκατο του όγκου του νέου μέσου.
Στάσιμη φάση
• Για την στάσιμη φάση κατ' αρχήν διακρίνουμε
δύο πιθανούς λόγους διακοπής της ανάπτυξης.
– ο πρώτος είναι η εξάντληση κάποιου θρεπτικού
συστατικού
– ο δεύτερος η συσσώρευση κάποιου τοξικού
προϊόντος.
• Στην πρώτη περίπτωση η τελική συγκέντρωση
της βιομάζας θα είναι απλώς:
x s = xo + Y X/S S o
Στάσιμη φάση
• Στη δεύτερη περίπτωση πρέπει να βρούμε την
εξάρτηση του ειδικού ρυθμού ανάπτυξης από το
τοξικό προϊόν και να εκφράσουμε τον ρυθμό
παραγωγής του τοξικού προϊόντος συναρτήσει του
υποστρώματος και της βιομάζας.
• Η ανάπτυξη θα διακοπεί όταν ο ειδικός ρυθμός
ανάπτυξης μηδενιστεί, δηλαδή όταν η συγκέντρωση
φτάσει κάποιο συγκεκριμένο επίπεδο.
• Επιλύνοντας τα ισοζύγια από κοινού βρίσκουμε τις
τελικές συγκεντρώσεις του υποστρώματος και της
βιομάζας.
Στάσιμη φάση
• Για να μπορέσουμε να διακρίνουμε
πειραματικά τον λόγο της διακοπής της
ανάπτυξης, δεν έχουμε παρά να αραιώσουμε
ένα δείγμα με νερό.
• Αν παρατηρηθεί περαιτέρω ανάπτυξη ο λόγος
είναι η συσσώρευση τοξίνης, ενώ αν δεν
παρατηρηθεί πρόκειται σίγουρα για
εξάντληση υποστρώματος.
Φάση θανάτου
-k d t
x = xs e
• Η σταθερά kd είναι συνάρτηση της
θερμοκρασίας και προσδιορίζεται εύκολα από
την κλίση διαγράμματος του λογαρίθμου της
βιομάζας συναρτήσει του χρόνου.
ΔΟΜΗΜΕΝΑ ΜΗ
ΚΑΤΑΝΕΜΗΜΕΝΑ ΜΟΝΤΕΛΑ
• Όταν
– ενδιαφέρει η εξέλιξη κάποιου συγκεκριμένου
κυτταρικού συστατικού, ή
– απαιτείται η αναγνώριση των κύριων συστατικών του
κυττάρου για την επαρκή περιγραφή της συνολικής
κινητικής ανάπτυξης (όταν δηλαδή η παραδοχή της
εξισορροπημένης ανάπτυξης δεν ευσταθεί)
• πρέπει να επιστρατεύσουμε τις γνώσεις μας για
την βιοχημεία και φυσιολογία του κυττάρου για
την ανάπτυξη κατάλληλων δομημένων μοντέλων.
Πως συνδέεται η συγκέντρωση κάποιου
συγκεκριμένου ενδοκυτταρικού παράγοντα με τον
ρυθμό ανάπτυξης του κυττάρου που αναπτύσσεται
σε κάποιο μέσο όγκου VR?
• Ας υποθέσουμε ότι το σύστημα είναι κλειστό.
• Έστω cj η συγκέντρωση του συστατικού j μέσα στο
κύτταρο, σε moles j ανά μονάδα κυτταρικού όγκου.
• Αν η συγκέντρωση της βιομάζας είναι x (σε g/l) τότε
η συνολική βιομάζα είναι VRx
• Αν η πυκνότητα ρc των κυττάρων θεωρείται
σταθερή, τότε ο ολικός όγκος της βιολογικής
φάσεως (του συνόλου των κυττάρων) είναι (VRx)/ρc.
• Αρα ο συνολικός αριθμός των moles του j μέσα στο
μέσο είναι (VRxcj)/ρc.
Αν ο καθαρός ρυθμός σχηματισμού του j ανά μονάδα
χρόνου ανά μονάδα όγκου μέσα στο κύτταρο είναι rfj τότε
το ισοζύγιο του j μπορεί να γραφεί ως ακολούθως:
d V R x cj
VRx
(
)=
r fj


dt
c
c
Θεωρώντας τον όγκο και την πυκνότητα σταθερές:
dx
dc j
= x r fj 
cj +x
dt
dt
1 dx
dc j
= r fj cj 
dt
x dt
dc j
= r fj -  c j
dt
dc j
= r fj -  c j
dt
• Ο πρώτος όρος είναι ο καθαρός ρυθμός μεταβολής
λόγω αντίδρασης
• Ο δεύτερος όρος αντιπροσωπεύει την μείωση που
παρατηρείται εξ αιτίας της αραίωσης που προκύπτει
από την αύξηση του κυτταρικού όγκου λόγω
ανάπτυξης.
• Η σχέση αυτή είναι ιδιαίτερα σημαντική μια και
χρησιμοποιείται για την κατασκευή οιουδήποτε
δομημένου μοντέλου.
ΔΙΑΜΕΡΙΣΜΑΤΙΚΑ ΔΟΜΗΜΕΝΑ
ΜΟΝΤΕΛΑ
• Τα διαμερισματικά μοντέλα χωρίζουν την
βιομάζα σε επί μέρους τμήματα.
• Έτσι π.χ. μπορούμε να διακρίνουμε ένα
συνθετικό μέρος βιομάζας (βασικά RNA) και ένα
δομικό μέρος (DNA και πρωτείνες).
• Εναλλακτικά μπορούμε να διακρίνουμε την
βιομάζα σε ένα αφομοιωτικό τμήμα και ένα
συνθετικό τμήμα, κ.ο.κ.
• Ο διαχωρισμός σε τμήματα έχει κάποια βιοχημική
και δομική βάση, αλλά δεν αντιστοιχεί επακριβώς
σε κάποιο πραγματικό μηχανισμό.
Williams
• Η βιομάζα αποτελείται από δύο τμήματα 1 και 2.
Το πρώτο είναι συνθετικό και το δεύτερο δομικό.
• Το συνθετικό μέρος παράγεται με την
ταυτόχρονη κατανάλωση υποστρώματος και
διέπεται από συντελεστή απόδοσης Υ που
δίνεται σε g μέρους 1 / g υποστρώματος
• Το δεύτερο δομικό μέρος βιομάζας παράγεται με
ρυθμό ανάλογο του γινομένου ρ1ρ2 όπου ρi η
συγκέντρωση του μέρους i μέσα στο κύτταρο
• Αναγκαία και ικανή συνθήκη για διχοτόμηση
αποτελεί ο διπλασιασμός του μέρους 2
Για ένα καλώς αναδευόμενο αντιδραστήρα διαλείποντος
έργου (κλειστό ως προς μάζα), οι εξισώσεις που
περιγράφουν
(συγκέντρωση
την
μεταβολή
των
υποστρώματος),
βιομάζας) και ρi, i=1,2, είναι:
dx
= k 1[S] x
dt
d[S]
1
= - k 1[S] x
dt
Y
d 1
= k 1[S](  1 +  2 )- k 2  1  2 -   1
dt
d 2
= k2 1  2 -   2
dt
μεταβλητών
x
[S]
(συγκέντρωση
• το μοντέλο υποθέτει ότι ο ειδικός ρυθμός ανάπτυξης
είναι ανάλογος της συγκέντρωσης υποστρώματος
(μ=k1[S]) και ότι η πυκνότητα της βιομάζας ρκ=ρ1+ρ2
είναι σταθερή.
• Ουσιαστικά οι δύο πρώτες εξισώσεις είναι ισοδύναμες με
το μοντέλο δύο φάσεων μια και μπορούν να επιλυθούν
ανεξάρτητα από τις δύο τελευταίες.
• Επίσης πολλαπλασιάζοντας την πρώτη με Υ και
προσθέτοντας με την δεύτερη βλέπουμε ότι το άθροισμα
Υ[S]+x είναι σταθερό και ίσο με Υ[Sο]+xο.
• Αντικαθιστώντας το [S] συναρτήσει του x καταλήγουμε
ουσιαστικά στην λογιστική εξίσωση των Verhulst και
Pearl. Ορίζοντας f2=ρ2/ρκ (το κλάσμα του τμήματος 2)
παίρνουμε:
df 2
= - k 1 [S] f 2 + k 2   f 2 (1- f 2 )
dt
• Από την παραδοχή για διχοτόμηση
φαίνεται ότι ο αριθμός των κυττάρων είναι
ανάλογος του f2x
• επομένως το μέγεθος των κυττάρων είναι
ανάλογο του x/(f2x)=1/f2
Καλλιέργεια σε αντιδραστήρα διαλείποντος
έργου με μαγιά που προέρχεται από στάσιμη
φάση (ρ1=0)
Σχετικό μέγεθος 1/f2-1
Βιομάζα x/(Y.So)
Υπόστρωμα S/So
Αριθμός f2.x/(Yso)
K1=0,005
K2=0,1
Y=0,5
So=270
Xo=5
Προβλέπει ορισμένα ποιοτικά χαρακτηριστικά που
συνήθως παρατηρούνται σε πραγματικές
καλλιέργειες:
(α) Ύπαρξη φάσης καθυστέρησης κατά την οποία
αυξάνεται το μέγεθος των κυττάρων χωρίς αύξηση του
αριθμού τους.
(β) Ύπαρξη φάσης εκθετικής ανάπτυξης κατά την οποία και
το μέσο κυτταρικό μέγεθος παίρνει την μεγαλύτερη του
τιμή.
(γ) Μεταβολή της σύστασης των κυττάρων κατά τον κύκλο
ανάπτυξης.
(δ) Ύπαρξη στάσιμης φάσης με μικρά σε μέγεθος κύτταρα.
(ε) Ανομοιόμορφη μεταβολή του αριθμού και της μάζας
των κυττάρων.
ΜΗΧΑΝΙΣΤΙΚΑ ΔΟΜΗΜΕΝΑ ΜΟΝΤΕΛΑ
• Τα μηχανιστικά μοντέλα κάνουν άμεση χρήση των
γνώσεων που έχουμε για τον μεταβολισμό ή/και την
γενετική του κυττάρου.
• τα μοντέλα αυτά περιλαμβάνουν μόνο στοιχεία που:
– κρίνονται ως απαραίτητα για την περιγραφή της κινητικής
(μια και παίζουν καθοριστικό ρόλο στην παρατηρούμενη
κινητική και στοιχειομετρία) ή
– συνδέονται άμεσα με την κινητική επί μέρους
συστατικών που μας ενδιαφέρουν (π.χ. προϊόντων).
• καθοριστικό ρόλο παίζει η δυνατότητα μέτρησης
των μεταβλητών που περιλαμβάνονται στο μοντέλο.
Παραγωγή κάποιας πρωτείνης P που
λαμβάνει υπόψη του τις διεργασίες της
μεταγραφής και μετάφρασης:
Αν ο γόνος που έχει την πληροφορία για την παρασκευή της
πρωτείνης είναι G, το ισοζύγιο μάζας για το μήνυμα RNA είναι:
d[mRNA]
= k p  [G] - k d [mRNA] -  [mRNA]
dt
όπου kp η κινητική σταθερά μεταγραφής, kd η κινητική
σταθερά της απενεργοποίησης/ καταστροφής του mRNA.
Ο συντελεστής αποτελεσματικότητας η συμπεριλαμβάνεται για
να περιγράψει πιθανές τροποποιήσεις του ρυθμού μεταγραφής.
Το ισοζύγιο μάζας για την πρωτείνη μπορεί να
γραφεί ως:
d[P]
= k q  [mRNA] - k e [P] -  [P]
dt
Εδώ ξ είναι ο συντελεστής αποδοτικότητας της μετάφρασης.
Αν υποθέσουμε εξισορροπημένη ανάπτυξη,
οι χρονικές παράγωγοι στο αριστερό μέρος των
εξισώσεων μηδενίζονται, και οι αλγεβρικές
εξισώσεις που προκύπτουν δίνουν:
k p k q [G]
[P] =
( k d +  )( k e +  )
• Η σχέση αυτή δίνει την συγκέντρωση του
πρωτεϊνικού προϊόντος συναρτήσει των
ρυθμών ανάπτυξης, μεταγραφής και
μετάφρασης
ΚΙΝΗΤΙΚΗ ΑΠΟΣΤΕΙΡΩΣΗΣ
• Αποστείρωση ενός μέσου ονομάζουμε την εξόντωση ή
απομάκρυνση όλων των μικροοργανισμών με ταυτόχρονη
απενεργοποίηση όλων των ιών.
• Η αποστείρωση είναι μία ιδιαίτερα βασική διεργασία που
συνήθως προηγείται της κύριας βιοτεχνολογικής διεργασίας
που περιλαμβάνει ανάπτυξη μικροοργανισμών (αποστείρωση
θρεπτικού μέσου).
Στις φαρμακοβιομηχανίες και στις βιομηχανίες
τροφίμων χρησιμοποιούνται μία ή περισσότερες
από τις ακόλουθες βασικές μεθόδους
αποστείρωσης:
(α) υγρή θέρμανση (ατμοί)
(β) ξερή θέρμανση (κλίβανοι)
(γ) μικροβιοκτόνοι παράγοντες όπως χλώριο, οξείδιο του
αιθυλενίου και φορμαλδεύδη
(δ) ακτινοβολία
(ε) διήθηση
Μία κατάλληλη επεξεργασία φέρνει την επιθυμητή
απομάκρυνση μικροοργανισμών χωρίς να προκαλέσει βλάβες
στα προϊόντα.
Η πιο συνηθισμένη διεργασία είναι η χρήση υπέρθερμων
ατμών.
Η κινητική της θερμικής αποστείρωσης είναι πρώτης τάξης ως προς την
συγκέντρωση των ζώντων οργανισμών, δηλαδή περιγράφεται από την
εξίσωση:
dn
= - kd n
dt
η οποία δίνει:
-k d t
n(t)= no e
όπου no η αρχική συγκέντρωση μικροοργανισμών.
Η σταθερά kd αυξάνει με την θερμοκρασία σύμφωνα με την σχέση τύπου
Arrhenius:
k d = k do e
Ed
RT
Κινητική αποστείρωσης
Επιζόν κλάσμα πληθυσμού
1
10-1
50οC
10-2
55oC
10-3
60oC
10-4
2
4
6
8
10 χρόνος (min)
Όταν ο αριθμός μειωθεί σημαντικά μόνο μία
στοχαστική περιγραφή του
φαινομένου της αποστείρωσης είναι στην πραγματικότητα δικαιολογημένη.
Αν υποθέσουμε ότι:
 η τύχη κάθε οργανισμού είναι ανεξάρτητη από την τύχη των άλλων
 αρχικά έχουμε ολικό αριθμό Νο
 η πιθανότητα ένας οργανισμός να παραμένει ζωντανός μετά από χρόνο t είναι e-kt
τότε η πιθανότητα να έχουμε Ν οργανισμούς εν ζωή μετά από χρόνο t είναι:
P N (t)=
No!
N
( - N)
( e-kt ) (1- e-kt ) N o
( N o - N)!N!
Επομένως η πιθανότητα να έχουμε εξοντώσει όλους τους οργανισμούς σε χρόνο t
είναι:
No
Po (t)= 1-(1- e-kt )
Αυτή η πιθανότητα είναι χρήσιμη για περιπτώσεις που απαιτείται σχεδόν τέλεια
αποστείρωση.
Μοντέλο δύο φάσεων
• Πολλοί οργανισμοί δεν εξουδετερώνονται
άμεσα, αλλά μετατρέπονται πρώτα σε
σπόρια, τα οποία εξουδετερώνονται στη
συνέχεια με άλλη κινητική.
• Σε αυτή την περίπτωση έχουμε:
ΝR  NS  ND
με σταθερές KR και KS αντίστοιχα
dN R
=  K R NR
dt
dN S
= K R NR  K S NS
dt
Η επίλυση με ΝR(0)=No και ΝS(0)=0 δίνει:
N K R   K st K s  K Rt 
 e 

=
e
No K S 
KR

όπου Ν=ΝR +NS
ΠΡΟΒΛΗΜΑ
Κατά τη θερμική απενεργοποίηση σπορίων βακτηρίου στους 121ºC προέκυψαν τα
δεδομένα του Πίνακα που παρατίθεται.
Ποιο πρότυπο ακολουθεί η θερμική καταστροφή του μικροοργανισμού;
Εκτιμείστε την προσαρμογή των δεδομένων στο προτεινόμενο μοντέλο και υπολογίστε
τις σχετικές σταθερές.
Χρόνος (sec)
Αριθμός ζωντανών σπορίων
(Ν)
Χρόνος (sec)
Αριθμός ζωντανών σπορίων (Ν)
0
1
2
3
4
5
6
7
100
99
98
95
91
88
83
79
8
9
10
15
20
30
40
50
75
70
66
46
31
13
5
2
Επίλυση
Χρόνος (sec)
0
1
2
3
4
5
6
7
8
9
10
15
20
30
40
50
Ν/Νο
1,000
0,990
0,980
0,950
0,910
0,880
0,830
0,790
0,750
0,700
0,660
0,460
0,310
0,130
0,050
0,020
ln(Ν/No)
0,000
-0,010
-0,020
-0,051
-0,094
-0,128
-0,186
-0,236
-0,288
-0,357
-0,416
-0,777
-1,171
-2,040
-2,996
-3,912
Χρόνος (sec)
ln(N/No)
0
10
20
30
40
50
60
0
-0,5
-1
-1,5
-2
-2,5
-3
-3,5
-4
-4,5
• Η προσαρμογή ευθείας στα σημεία του γραφήματος δίνει
συντελεστή συσχέτισης R2=0,9875
• ενώ από την κλίση της ευθείας υπολογίζεται ο ειδικός
ρυθμός θερμικής καταστροφής k=0,080 sec-1.
Στη συνέχεια εξετάζεται η προσαρμογή του μοντέλου θερμικής καταστροφής δύο σταδίων
Επιλύοντας την εξίσωση
N K R   K st K s  K Rt 
 e 

=
e
No K S 
KR

χρησιμοποιώντας αλγόριθμο μη γραμμικής βελτιστοποίησης προκύπτει συντελεστής συσχέτισης
R2=0.9999, ενώ οι τιμές των σταθερών υπολογίζονται ίσες με:
KS = 0.098 σπόρια s-1
KR = 0.152 σπόρια s-1
1,0
0,8
N/No
0,6
0,4
0,2
0,0
0
10
20
30
40
50
× ñüí ï ò (sec)
Από τα δύο πρότυπα που εξετάστηκαν, το μοντέλο δύο σταδίων εμφανίζει τον
υψηλότερο συντελεστή συσχέτισης, επομένως μπορεί να περιγράψει πιο
ικανοποιητικά τη θερμική απενεργοποίηση των σπορίων του βακτηρίου.
ΑΠΟΛΥΜΑΝΣΗ
• Απολύμανση ονομάζουμε την επιλεκτική μείωση
οργανισμών που μπορεί να μη συνοδεύεται από πλήρη
εξόντωση σαπροφυτικών οργανισμών και ενδόσπορων.
• Η απολύμανση απαντάται στην πράξη σε ιατρικές
εφαρμογές καθώς και στην επεξεργασία νερού και
λυμάτων.
• Για τον καθαρισμό του νερού το πιο συνηθισμένο
απολυμαντικό μέσο είναι το χλώριο.
• Η κινητική της απολύμανσης είναι ανάλογη της κινητικής
θερμικής αποστείρωσης (κινητική πρώτης τάξης ως προς
την συγκέντρωση οργανισμών).
Κινητική απολύμανσης
• Η εξάρτηση της κινητικής από την συγκέντρωση
του απολυμαντικού παράγοντα δίνεται συνήθως
εμπειρικά από την σχέση:
Cnt=σταθερά
όπου:
t ο χρόνος που απαιτείται για να επιτευχθεί μία
δεδομένη ποσοστιαία μείωση του πληθυσμού
C η συγκέντρωση του απολυμαντικού παράγοντα.
π.χ. αν t είναι ο χρόνος για 99% μείωση του πληθυσμού
E.coli σε διάλυμα υποχλωριώδους οξέος, η σταθερά είναι
0,24 ενώ n=0,86.