Định luật tác dụng khối lượng (Law of Mass Action) và Mô hình ức chế
Trong mô hình hóa toán sinh (mathematical biology), đặc biệt là động học dược lý (pharmacodynamics) và động học enzyme, Định luật Tác dụng khối lượng (Law of Mass Action) là nền tảng số một để chuyển hóa các hiện tượng tương tác sinh học thành các phương trình vi phân (mathematical models).
Bài viết này sẽ đi từ trực giác cơ bản của định luật này đến việc thiết lập phương trình cho tỷ lệ ức chế trong trạng thái dừng (steady state) – một dạng cấu trúc thường xuyên xuất hiện trong các luận án tại OUCRU.
1. Trực giác từ phản ứng cơ bản $A + B \xrightarrow{k_+} C$
Hãy tưởng tượng một không gian nơi các phân tử loại $A$ và $B$ di chuyển ngẫu nhiên (chuyển động Brown). Để tạo ra một phân tử $C$, một hạt $A$ và một hạt $B$ phải va chạm trực tiếp với nhau.
Xác suất để một va chạm xảy ra tỷ lệ thuận với mật độ không gian (nồng độ) của $A$ và mật độ của $B$. Do đó, tốc độ hình thành của $C$ (rate of change) tỷ lệ thuận với tích số của $A$ và $B$:
$$ \frac{dC}{dt} = k_+ \cdot A \cdot B $$Trong đó:
- Các biến $A, B, C$ biểu diễn nồng độ (concentration/density) tại một thời điểm $t$.
- Hằng số $k_+$ (rate constant) đại diện cho xác suất một vụ va chạm thực sự dẫn đến phản ứng thành công.
Đây chính là cốt lõi của Định luật Tác dụng khối lượng: Tốc độ của một phản ứng hóa học tỷ lệ thuận với tích nồng độ của các chất tham gia.
2. Mô hình hóa tỷ lệ ức chế: $N_a + n \cdot CA \rightleftharpoons N_i$
Giả sử ta có một hệ thống thụ thể (receptor) hoặc kênh ion. Gọi $N$ là tổng mật độ các thụ thể này. Tại bất kỳ thời điểm nào, thụ thể có thể ở một trong hai trạng thái động lực học (dynamical states):
- $N_a$: Mật độ thụ thể đang hoạt động (Active).
- $N_i$: Mật độ thụ thể bị ức chế/bất hoạt (Inactive).
Áp dụng Định luật Tác dụng khối lượng, ta có thể viết phương trình vi phân mô tả tốc độ thay đổi của trạng thái bị ức chế $N_i$:
- Flux đi tới (Inactivation flux): $k_+ \cdot N_a \cdot (CA)^n$ (Lưu ý: Vì cần $n$ hạt $CA$ cùng va chạm nên nồng độ $CA$ được nhân lên $n$ lần, tạo thành lũy thừa bậc $n$).
- Flux quay lui (Activation flux): $k_- \cdot N_i$
3. Tỷ lệ ức chế ở Trạng thái dừng (Steady State)
Trong nhiều bài toán dược lý học, các phản ứng gắn kết ở mức phân tử xảy ra cực kỳ nhanh so với sự tiến triển của bệnh lý. Do đó, ta thường áp dụng giả thiết trạng thái tựa dừng (quasi-steady-state assumption), tức là cho $\frac{dN_i}{dt} \approx 0$.
$$ k_+ N_a (CA)^n = k_- N_i $$$$ k_+ (N - N_i) (CA)^n = k_- N_i $$$$ k_+ N (CA)^n - k_+ N_i (CA)^n = k_- N_i $$$$ k_+ N (CA)^n = N_i (k_- + k_+ (CA)^n) $$$$ N_i = N \frac{k_+ (CA)^n}{k_- + k_+ (CA)^n} $$Chia cả tử và mẫu cho $k_+$, và định nghĩa Hằng số phân ly (Dissociation constant) $K_d = \frac{k_-}{k_+}$, ta có được tỷ lệ ức chế (fraction of inhibition):
$$ \frac{N_i}{N} = \frac{(CA)^n}{K_d + (CA)^n} $$4. Kết luận
Công thức cuối cùng mà ta vừa dẫn ra chính là Phương trình Hill (Hill Equation) kinh điển.
- Khi nồng độ thuốc/chất ức chế $(CA)$ rất nhỏ, tỷ lệ ức chế tiến về 0.
- Khi $(CA)$ rất lớn (tiến tới vô cực), tỷ lệ ức chế bão hòa tiến về 1 ($100%$ thụ thể bị khóa).
- Tham số $n$ (Hệ số Hill - Hill coefficient) quyết định độ dốc của đường cong. Khi $n > 1$, hệ thống có tính hiệp đồng (cooperativity), nghĩa là việc một phân tử $CA$ gắn vào sẽ làm tăng khả năng gắn của các phân tử tiếp theo, tạo ra một sự chuyển trạng thái cực kỳ đột ngột và “sắc nét”.
Việc hiểu rõ cách xây dựng phương trình Hill từ những viên gạch cơ bản của Định luật Tác dụng khối lượng là yêu cầu bắt buộc khi xây dựng các mô hình toán sinh (mathematical models) thực tế.