[수학 × 생명과학 세특] : 5. 미분방정식을 활용한 약물의 체내 흡수와 반감기 모델링

📑 [미적분×의약학] 미분방정식을 이용한 약물 반감기 모델링과 복용 주기 설계

1. 교과 연결 핵심 개념

• 수학: 미적분(미분방정식, 변수분리법, 정적분), 수학Ⅰ(지수함수와 로그함수, 상용로그와 자연로그)
• 과학: 화학Ⅱ(반응 속도, 1차 반응), 생명과학Ⅱ(세포막을 통한 물질 수송, 신장의 배설 작용)

2. 탐구 동기 및 목적

일상적으로 복용하는 약물의 효과가 시간이 지남에 따라 감소하는 이유를 생명과학적 현상에 머무르지 않고, "현재 체내 약물 농도와 감소 속도가 비례한다"는 미분방정식적 가정을 도입하여 정량화하고자 함. 약물의 종류에 따라 하루 복용 횟수가 다른 원인을 수학적 모델링을 통해 규명하고, 최적의 복용 주기를 설계하는 프로세스를 탐구함.

3. 수학적 유도 프로세스 (고교 과정 중심)

① 미분방정식 세우기 (1차 반응 모델)

체내 약물 제거 속도는 현재 존재하는 약물량 $A$에 비례하므로, 시간 $t$에 따른 약물량의 변화율 $\frac{dA}{dt}$는 다음과 같은 1차 미분방정식으로 표현됩니다.

$$\frac{dA}{dt} = -kA \quad (k > 0, \, k \text{는 제거 속도 상수})$$

② 변수분리와 적분

양변을 $A$와 $t$에 대해 각각 분리합니다.

$$\frac{1}{A} \, dA = -k \, dt$$

양변을 부정적분합니다.

$$\int \frac{1}{A} \, dA = \int -k \, dt \implies \ln|A| = -kt + C_1$$

약물의 양 $A>0$이므로 절댓값을 벗기고 지수형태로 변환합니다.

$$A(t) = e^{-kt + C_1} = e^{C_1} \cdot e^{-kt} = C \cdot e^{-kt}$$

초기 조건 $t=0$일 때의 약물량을 $A_0$라 하면 $A(0) = C = A_0$가 되므로, 최종적인 지수감소 함수식이 도출됩니다.

$$A(t) = A_0 e^{-kt}$$

③ 로그 연산을 통한 반감기($t_{1/2}$) 공식 유도

반감기는 체내 약물량이 처음의 절반($\frac{A_0}{2}$)이 되는 시간입니다. $A(t_{1/2}) = \frac{A_0}{2}$를 대입합니다. $$\frac{A_0}{2} = A_0 e^{-k \cdot t_{1/2}} \implies \frac{1}{2} = e^{-k \cdot t_{1/2}}$$ 양변에 자연로그($\ln$)를 취하면 다음과 같습니다.

$$\ln\left(\frac{1}{2}\right) = -k \cdot t_{1/2} \implies -\ln 2 = -k \cdot t_{1/2}$$

$$t_{1/2} = \frac{\ln 2}{k} \approx \frac{0.693}{k}$$

• 인사이트: 약물의 반감기는 초기 약물량($A_0$)과 무관하며, 오직 체내 제거 속도 상수($k$)에 의해서만 결정됨을 수학적으로 증명할 수 있습니다.

4. 학생 자기주도 활동: 엑셀(Excel)을 활용한 데이터 시각화

고등학생 수준에서 복잡한 실험 장비 없이도 구글 스프레드시트나 엑셀을 활용해 수치적 가독성을 높이는 시각화 활동을 수행할 수 있습니다.

• 조건 설정: 초기 약물량 $A_0 = 100\text{mg}$, 제거 상수 $k = 0.2$
• 수식 적용: $A(t) = 100 \cdot \exp(-0.2t)$

시간 (t, 시간)	계산된 약물량 (A(t), \text{mg})	순간 변화율 (기울기, −0.2A)
0	100.0	-20.0
2	67.0	-13.4
4	44.9	-9.0
6	30.1	-6.0
8	20.2	-4.0
10	13.5	-2.7

📊 엑셀 탐구 분석 포인트 (보고서 작성용)

1. 곡선의 기하학적 특징: 그래프를 그렸을 때 직선이 아닌 아래로 볼록한 지수감소 곡선이 나타납니다.
2. 변화율의 변화: 시간이 지날수록(즉, 체내 약물 농도가 낮아질수록) 순간 변화율(기울기)의 절댓값이 감소합니다. 이는 "체내 잔류량이 적을수록 배설·대사되는 절대적인 양도 줄어든다"는 생리적 현상을 수학적으로 뒷받침합니다.

5. 심화 확장: 다회 복용 주기와 최소 유효 농도(MEC) 모델링

실제 의학 현장에서는 약물을 1회만 투여하지 않고 일정 주기마다 반복 투여합니다. 약효가 나타나는 최소 유효 농도(MEC)와 독성이 나타나는 최소 독성 농도(MTC) 사이의 '유효 혈중 농도'를 유지하는 시뮬레이션을 수행합니다.

• 가정: 반감기가 약 3.47시간($k=0.2$)인 약물 $100\text{mg}$을 8시간 간격으로 반복 복용할 때.
• 수학적 모델링 (수열의 극한 연결):
  • $t=0$ (1회 복용 직후): $A_1 = 100$
  • $t=8$ (2회 복용 직전 남은 양): $100 \cdot e^{-0.2 \times 8} \approx 20.2$
  • $t=8$ (2회 복용 직후 총량): $A_2 = 100 + 100e^{-1.6} \approx 120.2$
  • $t=16$ (3회 복용 직후 총량): $A_3 = 100 + A_2e^{-1.6} = 100 + 100e^{-1.6} + 100e^{-3.2} \approx 124.3$

이를 일반화하면 $n$회 복용 직후의 약물량은 첫째항이 $100$, 공비가 $e^{-1.6} \approx 0.202$인 등비수열의 합이 됩니다. 무한히 복용한다고 가정했을 때 체내 최대 약물량의 수렴값(무한등비급수)은 다음과 같습니다.

$$\lim_{n \to \infty} A_n = \frac{100}{1 - e^{-1.6}} \approx \frac{100}{1 - 0.202} \approx 125.3\text{mg}$$

• 인사이트: 반복 복용하더라도 약물이 체내에 무한히 쌓여 독성을 일으키지 않고, 수학적 극한값에 의해 일정 범위(약 $25.3\text{mg} \sim 125.3\text{mg}$) 내에서 정조 상태(Steady State)를 유지하게 됨을 파악할 수 있습니다.

6. 실생활 연계 및 진로 탐구 방향

• 약학·의예과: 수면제(반감기가 길어 다음 날 부작용 유발 가능성)나 카페인(유전자 변이에 따른 제거 속도 상수 $k$의 개인차)의 특성을 조사하고, 이를 기반으로 환자 맞춤형 약물 투여 설계(Pharmacokinetics)의 중요성을 강조.
• 화공·신소재공학과: 약물이 체내에서 서서히 방출되도록 제형을 제어하는 'DDS(약물 전달 시스템)' 기술과 0차 반응(일정한 속도로 방출) 모델의 수학적 차이점 비교 분석.

7. 학생부종합전형(세특) 기재 문장 추천 예시

[미적분] 약물의 체내 농도 변화율이 현재 농도에 비례하여 감소한다는 생리적 가정을 바탕으로, 1차 미분방정식을 세우고 변수분리법을 통해 약물 농도의 지수감소 함수 모델($A(t)=A_0e^{-kt}$)을 성공적으로 유도함. 양변에 자연로그를 취하는 연산 과정을 통해 약물의 반감기 공식($t_{1/2}=\frac{\ln 2}{k}$)을 도출하고, 반감기가 초기 투여량과 무관한 고유한 특성임을 수학적으로 증명함. 나아가 공학적 도구(Excel)를 활용해 시간별 잔류량 그래프를 시각화하여 순간 변화율의 감소 추이를 분석하였으며, 등비급수의 성질을 응용해 '일정 주기 반복 복용 시 체내 약물 농도의 수렴 한계값'을 정량적으로 계산해 냄. 이 과정을 통해 복용 간격 설정의 의학적 원리를 규명하는 등, 미적분 개념을 실생활의 보건·의약학적 문제 해결에 적용하는 탁월한 수학적 모델링 능력과 학제 간 융합 탐구 역량을 보여줌.