[논문 리뷰] Discovery and analysis of biochemical subnetwork hierarchies

Petter Holme, Mikael Huss|ArXiv.org|2003. 09. 23.

Bioinformatics and Genomic Networks참고 문헌 11인용 수 8

한 줄 요약

이 논문은 이원형 대사 네트워크에서 중심성 중심성의 개념을 활용하여 생화학적 경로 내 계층적 하위망을 탐지하는 그래프 이론적 방법을 제안한다. 높은 중심성 중심성을 가진 반응을 반복적으로 제거함으로써 계층적 계층도를 생성하며, 이는 핵심 기질 계층을 드러낸다. 주요 발견은 강한 피드백 루프(순환 기반 기준)를 가진 하위망만 자율성을 나타내며, *T. pallidum*에서 CoA 및 옻산염 경로가 자율 조절 모듈로 확인된다.

ABSTRACT

The representation of a biochemical network as a graph is the coarsest level of description in cellular biochemistry. By studying the network structure one can draw conclusions on the large scale organisation of the biochemical processes. We describe methods how one can extract hierarchies of subnetworks, how these can be interpreted and further deconstructed to find autonomous subnetworks. The large-scale organisation we find is characterised by a tightly connected core surrounded by increasingly loosely connected substrates.

연구 동기 및 목표

• 대규모 생화학적 네트워크 내 계층적 하위망을 탐지하기 위한 객관적 그래프 이론적 방법을 개발한다.
• 전통적인 기능적 하위망의 한계를 극복하기 위해 전역 네트워크 구조와 동역학을 통합한다.
• 피드백 루프를 통한 내부 조절을 유지하는 자율 하위망을 식별한다.
• 표준 커뮤니티 탐지 기준(예: 진입/진출 차수 균형)이 생화학적 네트워크에서 의미 있는지 테스트한다.
• 생물학적으로 관련성이 있는 하위망을 식별하기 위한 피드백 루프 기반 기준을 제안하고 검증한다.

제안 방법

• 기질(S)과 반응(R)으로 구성된 방향성 이원형 그래프로 생화학적 네트워크를 표현하며, 방향성 간선으로 연결한다.
• 반응의 효과적 중심성 중심성 $ c_B(r) = C_B(r)/k_{\text{in}}(r) $ 를 계산한다. 여기서 $ C_B(r) $ 는 반응 r를 통과하는 최단 경로의 수이다.
• Girvan–Newman 알고리즘을 적용하여 효과적 중심성 중심성이 가장 높은 반응을 반복적으로 제거함으로써 계층적 계층도를 생성한다.
• 약한 커뮤니티를 $ \sum_{s\in S'} K_{\text{in}}(s) > \sum_{s\in S'} K_{\text{out}}(s) $ 를 통해 식별하며, 하위망 내부로의 순방향 피드백을 나타낸다.
• 강한 커뮤니티를 $ K_{\text{in}}(s) > K_{\text{out}}(s) $ 를 통해 식별하며, 모든 $ s \in S' $ 에 대해 고유한 내부 피드백을 나타낸다.
• 피드백 루프 기준을 도입한다: 하위망에 대해 $ L(S') \leq \Lambda |S'| $ 를 만족하며, 여기서 $ L(S') $ 는 길이 >3인 기본 순환 상의 기질 수를 세며, $ \Lambda \in \{0, 1/2, 1\} $ 이다.

실험 결과

연구 질문

• RQ1 효율적 중심성 중심성과 같은 전역 네트워크 성질이 생화학적 하위망의 계층적 조직을 드러내는가?
• RQ2 왜 표준 커뮤니티 탐지 기준(예: 진입/진출 차수 균형)이 생화학적 네트워크에서 의미 있는 하위망을 식별하지 못하는가?
• RQ3 자율적이고 자율 조절되는 하위망과 느슨하게 연결된 하위망을 구분하는 구조적 특징은 무엇인가?
• RQ4 피드백 루프를 정량적으로 어떻게 활용하여 생물학적으로 관련성이 있는 하위망을 정의할 수 있는가?
• RQ5 *T. pallidum* 및 기타 생물에서 하위망이 내부 조절 자율성을 얼마나 잘 유지하는가?

주요 결과

• *T. pallidum* 대사 네트워크의 계층적 계층도는 밀접하게 연결된 핵심 영역을 중심으로 하여 점점 느슨하게 연결된 기질들이 둘러싸고 있음을 드러낸다.
• 표준 커뮤니티 기준(Eq. 3 및 Eq. 4)은 생화학적 네트워크에서 비현실적인 결과를 낳는다. 피드백 루프가 최상위 수준 이외에는 흔하지 않기 때문이다.
• $ \Lambda = 1 $ 조건을 만족하는 하위망(모든 기질이 피드백 순환 상에 있음)은 희귀하지만 생물학적으로 의미가 있으며, CoA 및 옻산염 경로가 자율 모듈로 식별된다.
• 확장된 옻산염 하위망은 $ \Lambda = 1/2 $ 를 만족하지만 $ \Lambda = 1 $ 은 만족하지 않아 부분적인 피드백 자율성을 나타낸다.
• 이 방법은 WIT 데이터베이스 내 43개 생물에서 비트리비얼하고 기능적으로 일관된 하위망을 성공적으로 식별하여 강건성을 확인한다.
• 계층도 상단의 가장 큰 클러스터는 항상 약한 커뮤니티 조건(Eq. 3)을 만족하며 순방향 피드백이 있음을 나타내지만, 강한 조건(Eq. 4)은 만족하지 않는다.