JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

化學動力學與藥物動力學

一杯正在反應的分子,和一劑進入血液的藥物,竟服從同一組微分方程——而一顆藥丸的效果為何會沿著一條你能預測的曲線上升、達峰、再消退。

質量作用定律:把化學變成一條 ODE

想像一個燒瓶裡的分子在溶液中相互碰撞。當合適的夥伴相撞時,反應就發生,所以反應速率應與它們相遇的頻率成正比——而這個頻率又與每種反應物的擁擠程度成正比。這就是質量作用定律:反應速率正比於各反應物濃度的乘積。把濃度寫成方括號,例如 [A],於是這個化學動力學模型幾乎自動寫成一組 ODE。

取最簡單的情形:A 衰變成 B。每個 A 每秒都有固定的轉變機率,所以消失速率正比於現存的 A 量:[A]' = -k [A]。到現在你已見過這條方程上百次了——它就是純粹的指數衰減,與支配放射性原子核和冷卻咖啡的那條可分離方程完全相同。其解為 [A](t) = [A]0 e^(-k t),而化學家的半衰期 ln(2)/k,正是穿上實驗袍的物理學家半衰期。

反應的級數,不過是其速率定律中濃度的總次方,而它徹底改變了方程的性質。一階步驟給出線性的 A' = -k A;像 A + A → 產物這樣的二階步驟則給出非線性的 A' = -k A^2,其解 1/[A] = 1/[A]0 + k t 衰減得慢得多,只像 1/t 那樣。速率定律的形態,直接從化學讀出,正是溫和的指數消退與頑固的代數尾巴之間的分野。

連串、平衡,以及一個拒絕安定下來的反應

大多數真實的化學是一條連鎖:A → B → C,其中 B 由 A 而生,再被消耗成 C。這給出一對耦合方程,[A]' = -k1 [A] 與 [B]' = k1 [A] - k2 [B]——正是你學過用特徵值求解的那種一階系統。中間物種 B 最有意思:當 A 充裕時它上升,A 耗盡後它又下降,畫出一個只有單一峰值的小丘。記住這幅圖,因為過一會兒 B 就會變成你血液中的藥物,而那個小丘就是你實際感受到的劑量。

當反應可逆,A ⇌ B 時,正向與逆向速率彼此較勁:[A]' = -kf [A] + kr [B]。系統停止變化,並非在一切都反應完時,而是在兩個速率達到平衡時——這是一個平衡,此時 kf [A] = kr [B],給出平衡常數 K = kf/kr = [B]/[A]。這其實是一個偽裝的吸引點:從任何初始混合出發,濃度都會滑向那個平衡點並停留,正是你在相線上研究過的穩定平衡行為。

房室:化學如何變成醫學

現在把燒瓶換成一具身體。藥物動力學問的是身體對藥物做了什麼——它如何被吸收、分佈與清除——而它用一個房室模型來回答。訣竅是假裝身體是少數幾個充分混合的水箱。最簡單的是單一房室:血液,內有藥量 A(t),以正比於現存量的速率把藥物流失給腎臟與肝臟。這又是 A' = -k A,同樣的一階衰減,因為清除機制每小時清掉固定的比例。

然而一顆藥丸不會一下子全部出現在血液中;它在腸道溶解,再滲透過去。所以我們用兩個房室——腸道 G 與血液 C——耦合方式恰如上面的 A → B → C 連鎖,帶有吸收速率 ka 與清除速率 ke。這就是經典的藥物動力學模型,而它的血藥濃度曲線,正是我們在中間物種 B 看到的那個單峰小丘。

G' = -ka G,            G(0) = D   (whole dose in gut)
C' =  ka G - ke C,     C(0) = 0   (none in blood yet)

solution (ka != ke):
C(t) = (D ka)/(ka - ke) * ( e^(-ke t) - e^(-ka t) )

rises from 0, peaks at  t_max = ln(ka/ke)/(ka - ke),  then fades like e^(-ke t)
兩房室口服劑量模型及其閉式血藥濃度——吸收把曲線推高,清除再把它帶回低處。

把那個解當故事來讀。早期 e^(-ka t) 項又大又負,所以隨著藥物被吸收,C 攀升;稍後那個快速項消亡,只剩 e^(-ke t),於是 C 沿著由半衰期 ln(2)/ke 設定的緩慢清除時間尺度衰減。中間的峰值就是你實際感受到的劑量;把它推得太高就有中毒之虞,讓它垂得太低藥物就失效。整套給藥的藝術——多少、多久一次——就是在讀這條曲線。

穩態、重複給藥與靜脈滴注

真實的病人每八小時服一顆藥,而非只服一次。每劑都在上一劑殘留之上疊加一個新的小丘,而由於方程是線性的,疊加原理讓你能直接把各個反應相加。藥物一劑一劑地累積,但不會無止境:清除量隨著上升的濃度同步增長,直到經過幾個半衰期後,你每個間隔加入的量恰等於你清除的量。血藥濃度於是在某個穩態附近趨於平台,在峰值與谷值之間溫和振盪。

一條穩定的靜脈滴注讓同樣的想法更乾淨。現在輸入是一個常數速率 R,所以 C' = R - ke C——一條帶常數驅動的一階線性方程。其解就是課本上的暫態加穩態的拆分:一個會消退的暫態 e^(-ke t),與一個留存的穩態平台 C_ss = R/ke。藥物朝 R/ke 攀升並趨平,而這趨近的速度完全由時間常數 1/ke 決定——約四個時間常數就能落到目標值的百分之幾以內。

模型隱藏了什麼,以及我們如何檢驗它

對這些虛構要誠實。身體不是一個攪拌槽,所以充分混合的房室是刻意的簡筆畫;有些藥物需要第三個房室來容納像脂肪那樣緩慢平衡的組織,而即便如此那仍是近似。質量作用的速率常數 k 也不是自然界的常數——它們隱藏了溫度、pH 與酵素,而在高劑量下清除機制會飽和,使得一階清除崩解為非線性的米氏動力學。認清一個模型的極限不是悲觀;它是建模中保障其安全的那一部分。

那我們憑什麼相信這些呢?我們去量測。在數個時刻抽血,再把 ka、ke 與劑量比例擬合到資料上——模型擬合與驗證。一個漂亮的診斷直接從數學中掉出來:吸收完成後,ln C(t) 對 t 是一條斜率為 -ke 的直線,所以把尾段畫在對數刻度上,就能用肉眼讀出清除速率。當擬合好的曲線能追上它未曾見過的新量測時,這模型才掙得了一點點信任——也僅僅是一點點。

退一步,欣賞這份統一。一瓶分子、腸道裡的一顆藥丸、一份放射性樣本、一個正在充電的 RC 電路——它們全是同一小撮一階線性 ODE 換上了不同的戲服。這正是整道階梯安靜的收穫:一旦你認得 A' = -k A 及其耦合的表親,你就能像讀一個自由落體那樣流利地讀懂化學與醫學,因為在底層,它們從來就不是不同的方程。