생리학적 조건에서 효소 촉매에 대한 용매 흐름 효과 모델링

초록

본 논문은 멀티입자 충돌 동역학(MPC)과 탄성 네트워크 모델을 결합해, 인산기 전달 효소인 포스포글리세르산 키나아제(PGK)의 촉매 사이클을 물리적 용매 흐름 하에서 시뮬레이션한다. 확률적 첫 통과 시간 밀도와 직접 시뮬레이션 결과를 통합해 전체 반응 주기의 시간 분포를 도출하고, 용매-효소 수소 결합 및 수소 흐름이 효소의 힌지 움직임과 기질 결합·폐쇄 과정에 미치는 영향을 정량적으로 분석한다.

상세 분석

이 연구는 효소 촉매를 ‘확률적 확산‑반응’ 시스템으로 모델링하고, 두 단계의 접근법을 결합한다. 첫 번째는 구형 쉘 안에서 확산 입자가 흡수 경계에 도달하는 첫 통과 시간(first‑passage time, FPT) 밀도를 해석적으로 구해, 기질이 효소 표면에 도달하는 확률적 특성을 정량화한다. 여기서 구의 내부 반경은 효소의 활성 부위 반경, 외부 반경은 용매 내 평균 기질 농도를 반영한다. 두 번째는 탄성 네트워크 모델(ENM)로 단백질 구조를 저해상도화하고, 멀티입자 충돌 동역학(MPC)으로 용매 입자의 수소 흐름과 점성 효과를 재현한다. ENM은 각 아미노산을 질점으로, 인접 질점 사이를 스프링으로 연결해 효소의 대규모 힌지 움직임을 자연스럽게 발생시킨다. MPC는 입자를 격자 셀에 배치하고, 일정 시간 간격마다 충돌 단계에서 무작위 회전 변환을 적용해 유체의 Navier‑Stokes 방정식을 통계적으로 만족한다.

시뮬레이션은 두 가지 결합 방식을 비교한다. (1) 명시적 상호작용: 효소와 기질 사이에 Lennard‑Jones 및 전기적 포텐셜을 부여해 직접적인 용매‑단백질 상호작용을 구현한다. (2) 암시적 결합: 효소·기질을 MPC 충돌 단계에 포함시켜, 용매 흐름이 효소에 전달되는 수소적 효과만을 반영한다. 또한, 수소 흐름이 전혀 없는 ‘무수소’ 시뮬레이션을 대조군으로 설정해 흐름 효과를 정량화한다.

결과는 다음과 같다. 수소 흐름이 존재할 때, 효소의 힌지 회전이 평균 30 % 가속화되고, 기질 결합 후 폐쇄(conformational closure) 단계의 평균 소요 시간이 20 % 감소한다. 특히, 효소 표면 근처에 존재하는 기질의 첫 통과 시간 분포가 비대칭적인 롱테일을 보이며, 흐름이 강할수록 롱테일이 억제돼 전체 사이클 시간이 크게 단축된다. 무수소 경우에는 힌지 움직임이 크게 저해되어, 전체 사이클 평균 시간이 약 2배 이상 늘어난다. 또한, 명시적 상호작용 모델은 수소 흐름에 의한 촉진 효과를 더 크게 나타내어, 실제 세포 내 복합 용매 환경을 보다 정확히 재현한다는 결론을 도출한다.

이러한 발견은 효소 동역학 연구에서 용매 흐름을 무시하면, 특히 대규모 구조 변화를 동반하는 촉매 메커니즘을 과소평가하게 됨을 시사한다. 또한, MPC‑ENM 결합 프레임워크는 실험적으로 측정하기 어려운 미시적 시간 스케일을 접근 가능하게 하여, 효소 설계와 약물 타깃팅에 새로운 계산적 도구를 제공한다.