[논문 리뷰] A simplified derivation of the Linear Noise Approximation

Edward W. Wallace|arXiv (Cornell University)|2010. 04. 24.

Gene Regulatory Network Analysis참고 문헌 11인용 수 9

한 줄 요약

이 논문은 화학 반응 시스템과 같은 이산 확률적 과정에 대한 선형 소음 근사(LNA)의 간소화된 유도를 제시한다. 기존의 마스터 방정식과 포커-플랑크 방정식 대신 확률적 차분 방정식과 화학 랑주인 방정식을 사용함으로써, 소음 편미분에 대한 확률적 미분 방정식을 직접 도출함으로써 기존의 시스템 크기 전개와 동일한 결정론적 및 선형 소음 역학을 회복하면서도, 일반 미분 방정식과 확률적 미분 방정식에 익숙한 청중에게 더 직관적인 접근을 제공한다.

ABSTRACT

Given a discrete stochastic process, for example a chemical reaction system or a birth and death process, we often want to find a continuous stochastic approximation so that the techniques of stochastic differential equations may be brought to bear. One powerful and useful way to do this is the system size expansion of van Kampen to express a trajectory as a small stochastic perturbation to a deterministic trajectory, using a small parameter related to the volume of the system in question. This is usually pursued only up to first order, called the Linear Noise Approximation. The usual derivation of this proceeds via the master equation of the discrete process and derives a Fokker-Planck equation for the stochastic perturbation, both of which are equations for the evolution of probability distributions. Here we present a derivation using stochastic difference equations for the discrete process and leading, via the chemical Langevin equation of Gillespie, directly to a stochastic differential equation for the stochastic perturbation. The new derivation, which does not yield the full system size expansion, draws more explicitly on the intuition of ordinary differential equations so may be more easily digestible for some audiences.

연구 동기 및 목표

• 마스터 방정식과 포커-플랑크 형식에 익숙하지 않은 연구자들에게 선형 소음 근사를 더 접근하기 쉽게 유도하기 위해.
• 기존의 마스터 방정식 기반 시스템 크기 전개를 대체하여, 확률적 차분 방정식과 화학 랑주인 방정식에 기반한 유도를 제공하기 위해.
• 일반 미분 방정식과 확률적 미분 방정식에 익숙한 청중을 대상으로 LNA 유도 과정을 더 직관적으로 만들기 위해.
• 특히 화학 랑주인 근사에 필요한 시간스케일과 시스템 크기 가정을 포함하여, LNA의 유효 조건을 명확히 하기 위해.

제안 방법

• 반응 속도와 시간스케일 τ 동안의 전이 사건을 이용해 이산 확률적 과정을 기술하며, 이때 반응 수가 많고 속도는 약간 변하지 않는 조건을 가정한다.
• τ-리프 근사를 사용하여 종의 농도 증가량을 독립적인 정규분포 확률변수의 합으로 모델링함으로써 화학 랑주인 방정식을 도출한다.
• ansatz K = Ωx + √Ωξ 를 도입하여 시스템을 결정론적 성분 x 와 1/√Ω 스케일의 확률적 편미분 ξ 로 분해한다.
• 결정론적 해 x 주변에서 이동항과 산란항에 대한 테일러 전개를 수행하고, 1/√Ω 에서의 주요 항들만 유지한다.
• 이동항이 결정론적 속도 함수의 자코비안에 의해 지배되고, 산란항이 산란 계수에 비례하는 선형 랑주인 방정식을 ξ 에 대해 도출한다.
• 이토 해석을 사용하여 표준 LNA 제작과의 일관성을 확보하며, 특히 내부 소음 과정에 대해 유의미하다.

실험 결과

연구 질문

• RQ1 마스터 방정식이나 포커-플랑크 방정식을 사용하지 않고도 선형 소음 근사를 도출할 수 있는가?
• RQ2 시스템 크기 전개를 확률적 미분 방정식의 관점에서 어떻게 재구성할 수 있으며, 이는 접근성을 향상시키는가?
• RQ3 τ-리프 근사는 화학 랑주인 방정식을 통해 직접 LNA를 유도하는 데 어떤 역할을 하는가?
• RQ4 간소화된 유도 과정이 기존의 시스템 크기 전개와 동일한 결과를 도출하는 조건은 무엇인가?
• RQ5 왜 일반적으로 LNA는 일阶까지 확장되지 않으며, 더 직관적인 유도 과정은 이러한 절단을 정당화할 수 있는가?

주요 결과

• 확률적 차분 방정식과 화학 랑주인 방정식을 사용한 간소화된 유도 과정은 기존의 시스템 크기 전개와 동일한 결정론적 및 선형 소음 역학을 도출한다.
• 확률적 편미분 ξ 는 결정론적 속도 함수의 자코비안에 의해 지배되는 이동항과 총 전이 속도의 제곱근에 비례하는 소음항을 가진 선형 랑주인 방정식을 만족한다.
• 이 유도 과정은 백색 소음 입력이 상호 독립적인 이유가 기본 전이가 동시에 하나의 집단만 영향을 미치기 때문임을 명확히 드러낸다.
• 결정론적 시스템의 고정점에서, 소음의 진폭은 더 단순한 형태에 따라 정방향 또는 역방향 전이 속도 중 하나로 동일하게 표현될 수 있다.
• 화학 랑주인 방정식이 성립하기 위한 조건—특히 반응 속도의 역수보다 충분히 큰 τ—은 LNA의 유효성에도 필수적이다.
• 이 방법은 일반 미분 방정식과 확률적 미분 방정식에 익숙한 청중을 대상으로 LNA에 도달하는 데 더 직관적이고 접근하기 쉬운 경로를 제공하며, 확률 분포의 진화 방정식에 대한 지식이 없어도 된다.