[논문 리뷰] A dynamic p53-mdm2 model with delay kernel

Florin George Horhat, Mihaela Neamţu외 1명|ArXiv.org|2006. 08. 22.

Cancer-related Molecular Pathways참고 문헌 11인용 수 3

한 줄 요약

이 논문은 DNA 손상에 대한 반응에서의 진동 다이내믹스를 연구하기 위해 분포 지연과 약한 커널 함수를 사용한 p53-Mdm2 피드백 루프의 지연 미분 방정식 모델을 개발한다. 허프 브리지테이션 분석을 통해 시스템이 지연 매개변수에 따라 불안정하거나 안정적인 주기적 진동을 보이는 초하위 브리지테이션(subcritical Hopf bifurcation)을 나타냄을 입증하며, 만성세포에서 실험적으로 관찰된 p53 진동에 대한 정교한 다이내믹스 해석을 제공한다.

ABSTRACT

Specific activator and repressor transcription factors which bind to specific regulator DNA sequences, play an important role in gene activity control. Interactions between genes coding such transcripion factors should explain the different stable or sometimes oscillatory gene activities characteristic for different tissues. In this paper, the dynamic P53-Mdm2 interaction model with distributed delays and weak kernel, is investigated. Choosing the delay or the kernel's coefficient as a bifurcation parameter, we study the direction and stability of the bifurcating periodic solutions. Some numerical examples are finally given for justifying the theoretical results.

연구 동기 및 목표

• 더 생물학적으로 현실적인 다이내믹스를 위해 분포 지연과 약한 커널 함수를 사용한 p53-Mdm2 피드백 루프 모델링.
• 허프 브리지테이션 분석을 통해 시스템이 진동 행동를 보이는 조건을 조사.
• 지연 또는 커널 매개변수를 변화시켜 발생하는 주기적 해의 방향, 안정성 및 주기를 결정.
• 생물학적으로 타당한 매개변수 값을 사용한 수치 시뮬레이션을 통해 이론적 결과를 검증.
• 이전의 이산 지연 모델을 개선하기 위해 분포 지연과 더 부드러운 커널 함수를 통합.

제안 방법

• 분포 지연을 포함한 p53 및 Mdm2 mRNA와 단백질 농도에 대한 지연 미분 방정식(DDE) 시스템을 수립.
• 지연 매개변수 $ \tau $ 또는 커널 계수 $ q $ 를 브리지테이션 매개변수로 간주하여 평형 상태의 국소 안정성 분석 및 허프 브리지테이션 조건 규명.
• 약한 커널 함수 $ k_i(s) = \frac{q^{p+1}_i}{p!} s^p e^{-q_i s} $ 를 사용하며, $ p=0 $ 인 경우 약한 커널로 분포 지연을 허용.
• 중앙다양체 이론 및 정규형 방법(Hassard의 방법)을 적용하여 발생하는 주기적 해의 방향, 안정성 및 주기를 계산.
• 시스템 매개변수와 지연에 기반한 핵심 브리지테이션 양수인 $ \mu_2 $, $ \beta_2 $, 및 $ T_2 $ 의 명시적 표현 유도.
• 생물학적으로 타당한 매개변수를 사용한 Maple 9.5 기반 수치 시뮬레이션을 수행하여 파형과 위상도를 시각화.

실험 결과

연구 질문

• RQ1 분포 지연을 가진 p53-Mdm2 시스템이 언제 허프 브리지테이션을 겪는가?
• RQ2 지연 커널의 선택(약한 vs. 이산)이 주기적 해의 안정성과 성격에 어떻게 영향을 미치는가?
• RQ3 브리지테이션 발생 주기적 해의 방향(초하위/초하위)과 안정성(안정/불안정)을 결정하는 요소는 무엇인가?
• RQ4 지연 매개변수에 따라 진동 주기는 어떻게 변화하는가?
• RQ5 이 모델은 이전의 이산 지연 모델보다 실험적으로 관찰된 p53 진동을 더 정확하게 재현할 수 있는가?

주요 결과

• 이산 지연 모델에서 $ \tau_2 = 3 $ 인 경우, 시스템은 $ \tau_0 = 9.54 $ 에서 초하위 허프 브리지테이션을 겪으며, $ \mu_2 = -0.420 $ 이므로 $ \tau > \tau_0 $ 일 때 브리지테이션 해가 존재함.
• 주기적 해는 궤도적으로 불안정함($ \beta_2 = 0.280 $)이며, 지연이 임계값을 초과할수록 주기가 증가함($ T_2 = 0.0005 $).
• 약한 커널 모델에서 $ q_2 = 0.5 $ 인 경우, 허프 브리지테이션은 초하위($ \mu_2 = -0.599 $)이며 주기적 해는 궤도적으로 안정함($ \beta_2 = -0.748 $), 지연이 $ \tau_1 > 32.37 $ 인 경우 주기가 증가함($ T_2 = 0.180 $).
• 모델은 DNA 손상에 대한 반응에서 실험적으로 관찰된 p53 다이내믹스와 일치하는 진동 행동를 성공적으로 재현함.
• 분포 지연과 약한 커널을 사용함으로써 더 부드러운 다이내믹스와 이전의 이산 지연 모델보다 더 생물학적으로 타당한 진동을 도출함.
• 수치 시뮬레이션은 이론적 예측을 확인하였으며, p53 및 Mdm2 단백질 농도에서 특징적인 파형과 위상도 패턴을 보이는 지속적인 진동을 보여줌.