자율신경계가 조절하는 심혈관 시스템의 수학적 모델링

초록

본 논문은 심혈관계의 자율신경 조절을 통합적으로 기술하는 수학적 모델을 제시한다. 72명의 피험자를 대상으로 500여 회의 임상 자율신경 검사(심전도, 동맥맥압, 호흡, 피부전도 등)를 수행하고, 시간·주파수 영역 분석과 제한 최적화를 통해 모델 파라미터를 추정하였다. 결과는 메이어파가 교감성 프리모터 뉴런의 리듬성 방출에 기인함을 뒷받침하고, 일부 피험자에서 관찰된 미주신경 매개 빈맥이 미주신경에 의한 혈관활성 물질 분비와 연관될 수 있음을 시사한다. 또한 연령에 따라 감소하는 부교감성 톤과, 휴식 시 교감성 톤이 프리모터 뉴런 및 화학수용체 활성에 의해 유지된다는 새로운 지표를 제시한다.

상세 분석

이 연구는 심혈관계 자율조절을 다층 구조의 수학적 프레임워크로 구현한 점이 가장 큰 강점이다. 먼저, 체성·내장·화학·기계 수용체 등 다양한 구심성 감각 뉴런을 입력으로 설정하고, 이들 신호가 시냅스 가중치와 시간 지연을 포함한 미세 조절 메커니즘을 통해 전뇌의 자율 전구조(pre‑motor) 뉴런에 전달되는 과정을 미분 방정식 계열로 모델링하였다. 전구조 뉴런은 교감·부교감 두 갈래로 분리되어 각각 전신성(심장, 혈관) 및 국소성(피부, 호흡) 효과기를 향해 전두엽성 전구신경(pre‑ganglionic)과 말초신경(post‑ganglionic) 신호를 발산한다. 특히, 신경전달물질(노르에피네프린, 아세틸콜린, 혈관활성 물질 등)의 합성·방출·재흡수·분해 과정을 Michaelis‑Menten 형태와 일차 반감기 모델로 결합함으로써 효과기에서의 혈관 저항·심박출량·심박동 간격 변화를 정량화하였다.

실험 설계는 연령대별(11‑82세) 72명을 대상으로 표준 자율신경 검사(베타 차단제, Valsalva, 깊은 호흡, 체위 변화 등)를 500회 이상 수행하고, ECG, 혈압 파형, 호흡, GSR, 피부 온도 등 다중 생리 신호를 동시에 기록하였다. 수집된 데이터는 고해상도 보간과 잡음 제거 후, 시간 영역(심박동 변이, 혈압 변동)과 주파수 영역(저주파·고주파 파워, 메이어파 스펙트럼)으로 분해되었다. 이러한 통계량을 제약조건(생리학적 범위, 신경전달물질 보존 법칙 등)으로 하는 비선형 최적화 문제에 입력함으로써 모델 파라미터(시냅스 가중치, 뉴런 발화 주기, 전송 지연, 물질 반감기 등)를 추정하였다.

시뮬레이션 결과는 두 가지 주요 현상을 설명한다. 첫째, 메이어파(≈0.1 Hz)는 교감성 프리모터 뉴런이 자체적인 펄스 발생기(pacemaker) 역할을 하여 리듬성 방출을 일으키는 것으로 재현되었으며, 이는 기존의 혈압‑피드백 이론을 보완한다. 둘째, 일부 피험자에서 Valsalva와 같은 미주신경 자극 시 관찰되는 급성 빈맥은 미주신경 말단에서 혈관활성 물질(예: NO, 아세틸콜린)의 급증이 심장 전도 속도를 역으로 촉진한다는 가설과 일치한다.

또한, 모델은 휴식 시 교감·부교감 톤을 정량화하는 지표를 제공한다. 연령에 따라 부교감 톤이 감소하고, 교감 톤은 상대적으로 유지되지만 화학수용체(특히 저산소·고이산화탄소 감지)의 활성도가 증가함을 보여준다. 이는 심혈관 질환 위험이 높은 고령 집단에서 교감성 과잉이 어떻게 발생하는지를 메커니즘적으로 설명한다.

한계점으로는 비침습적 측정에 의존해 신경 활동을 간접 추정한다는 점, 파라미터 식별에 사용된 최적화가 지역 최소값에 머물 가능성, 그리고 약물 차단 실험이 부족해 모델 검증이 제한적이라는 점을 들 수 있다. 향후 연구에서는 교감·부교감 차단제 투여 실험, 동물 모델의 직접 신경 전기 기록, 그리고 인공신경망 기반 파라미터 추정 기법을 결합해 모델의 일반화와 정확성을 높일 필요가 있다.