[논문 리뷰] Application of Fractional Calculus to Reaction-Subdiffusion Processes and Morphogen Gradient Formation

S. B. Yuste, E. Abad외 1명|arXiv (Cornell University)|2010. 06. 14.

Field-Flow Fractionation Techniques인용 수 4

한 줄 요약

이 논문은 분자의 기울기 형성에서 반응-하위분산 과정을 모델링하기 위해 분수적 미적분을 적용하며, 기억 효과와 공간에 따라 변하는 분해를 고려한 CTRW 기반의 분수적 반응-하위분산 방정식을 유도한다. 정적인 분자의 기울기 프로파일은 재반응성의 공간적 의존성과 비정상적 분산 지수 γ에 매우 민감하며, 비균일한 재반응성이 존재할 경우 γ와 분해 파rameter에 따라 프로파일이 감쇠하거나 일정하게 유지되거나 무한히 증가하는 등 다양한 형태를 띈다.

ABSTRACT

It is a well known fact that subdiffusion equations in terms of fractional derivatives can be obtained from Continuous Time Random Walk (CTRW) models with long-tailed waiting time distributions. Over the last years various authors have shown that extensions of such CTRW models incorporating reactive processes to the mesoscopic transport equations may lead to non-intuitive reaction-subdiffusion equations. In particular, one such equation has been recently derived for a subdiffusive random walker subject to a linear (first-order) death process. We take this equation as a starting point to study the developmental biology key problem of morphogen gradient formation, both for the uniform case where the morphogen degradation rate coefficient (reactivity) is constant and for the non-uniform case (position-dependent reactivity). In the uniform case we obtain exponentially decreasing stationary concentration profiles and we study their robustness with respect to perturbations in the incoming morphogen flux. In the non-uniform case we find a rich phenomenology at the level of the stationary profiles. We conclude that the analytic form of the long-time morphogen concentration profiles is very sensitive to the spatial dependence of the reactivity and the specific value of the anomalous diffusion coefficient.

연구 동기 및 목표

• 하위분산 운반과 반응(분해)을 고려한 분수적 미적분을 사용하여 분자의 기울기 형성 모델링.
• 긴 꼬리가 있는 대기 시간 분포를 갖는 연속 시간 랜덤 워크(CTRW)를 통해 기억 효과를 포함한 고전적 반응-확산 모델의 확장.
• 비정상적 분산과 공간적으로 의존하는 분해가 정적인 분자의 농도 프로파일의 형태와 강건성에 미치는 영향 탐구.
• 물리적으로 일관된 분수적 반응-하위분산 방정식을 유도하고, 음수 농도와 같은 비물리적 결과를 피하는 분석.
• 이러한 프로파일이 발달 생물학에 미치는 함의, 특히 외부 자극에 대한 강건성과 기울기 형태에 대한 고려.

제안 방법

• 기타 사망 과정을 포함한 CTRW 모델에서 리만-리우빌 분수 미분을 사용하여 분수적 반응-하위분산 방정식 유도.
• CTRW 형식을 통해 비마르코프 대기 시간 분포를 포함하여 γ < 1 인 하위분산 유도.
• 라플라스 변환 기법을 사용하여 균일 및 비균일 재반응성 경우의 분수 편미분 방정식 해법.
• 상수, 조각별 상수, 지수, 거듭제곱 법칙 감쇠 형태의 다양한 공간적 재반응성 k(x) 형태에 대한 정적 해 분석.
• 해석적 해가 존재하지 않는 경우 임시적 프로파일을 연구하기 위해 라플라스 변환의 수치적 역변환 수행.
• 유도된 방정식을 통해 시간 및 공간에 따라 변화하는 재반응성 k(x,t)에 대한 일반화를 위해 분수 시간 미분과 지수 기억 커널을 활용.

실험 결과

연구 질문

• RQ1 분수 미분으로 모델링된 하위분산은 고전적 확산과 비교해 분자의 기울기 형성의 형태에 어떻게 영향을 미치는가?
• RQ2 공간적으로 변하는 분해(재반응성)는 하위분산 운반에서 정적인 농도 프로파일에 어떤 영향을 미치는가?
• RQ3 비정상적 분산 지수 γ는 다양한 재반응성 함수 하에서 분자의 기울기 프로파일의 장기적 행동에 어떻게 영향을 미치는가?
• RQ4 이 모델은 하위분산 영역에서 비단조화적 또는 발산하는 프로파일을 예측할 수 있으며, 이러한 행동은 어떤 조건에서 발생하는가?
• RQ5 분자의 기울기 프로파일이 유입 흐름 변화에 대한 강건성은 하위분산 매개변수 γ와 분해 속도 k에 어떻게 의존하는가?

주요 결과

• 균일한 재반응성 조건에서는 정적인 분자의 기울기 프로파일이 지수적으로 감쇠되며, 분해 속도 k가 클수록 유입 흐름 변화에 대한 강건성이 증가한다.
• k₁ = 0 인 영역이 존재하는 조각별 상수 재반응성의 경우 불연속적인 농도 프로파일이 나타나며, 정적인 상태에 도달하지 못하며, 이는 일반적인 확산에서는 관찰되지 않는 새로운 효과이다.
• 지수 감쇠 재반응성의 경우 γ = 1 이면 무한대에서 프로파일이 일정해지지만, γ < 1 이면 무한히 증가한다.
• 거듭제곱 법칙 감쇠 재반응성의 경우, 큰 x에서의 정적 프로파일 행동은 감쇠 지수와 γ 간의 상호작용에 따라 결정되며, 감쇠, 일정, 또는 발산하는 프로파일이 될 수 있다.
• 정적 프로파일의 해석적 형태는 재반응성 k(x)의 함수 형태와 비정상적 분산 계수 γ의 값에 매우 민감하여 강한 비보편적 행동을 나타낸다.
• 유도된 분수적 반응-하위분산 방정식은 고전 모델을 일반화하며 기억 효과를 정확히 반영하여 음수 농도와 같은 비물리적 결과를 피한다.