무지개송어 에스트로젠 수용체 동형체 결합 친화도 계산 연구

초록

본 연구는 분자동역학 시뮬레이션과 진화유전학 분석을 통해 무지개송어(Oncorhynchus mykiss)의 네 종류의 에스트로젠 수용체(ERα1, ERα2, ERβ1, ERβ2)와 호르몬 17β‑에스트라디올(E2) 사이의 결합 친화도를 비교한다. 결과는 ERα1이 ERα2보다 E2에 더 높은 친화도를 보이며, ERα2의 리간드 결합 도메인이 완화된 선택압을 받는 반면, ERβ1과 ERβ2는 유사한 친화도와 강한 정화 선택을 유지함을 보여준다. 이는 ERα2가 새로운 기능(네오펑셔널라이제이션)으로 전환될 가능성을 시사한다.

상세 분석

본 논문은 무지개송어의 에스트로젠 수용체(ER) 동형체가 진화 과정에서 어떻게 기능적 분화를 겪었는지를 정량적으로 평가하기 위해 두 가지 핵심 접근법을 결합하였다. 첫 번째는 구조 기반의 분자동역학(MD) 시뮬레이션으로, 각 ER 동형체의 리간드 결합 도메인(LBD)을 동종 모델링한 뒤 100 ns 수준의 시뮬레이션을 수행하고, MM‑PBSA(또는 MM‑GBSA) 방법을 이용해 17β‑에스트라디올(E2)과의 결합 자유에너지(ΔG_bind)를 추정하였다. 시뮬레이션 설정은 CHARMM36 힘장, TIP3P 물 모델, 310 K, 1 atm 조건을 사용했으며, 시스템 평형화 후 50 ns 구간을 분석에 활용하였다. 결과는 ERα1‑E2 복합체가 평균 ΔG_bind ≈ ‑12.3 kcal/mol 로 가장 강한 결합을 보인 반면, ERα2‑E2는 ‑9.1 kcal/mol 로 유의하게 약하였다. 반면 ERβ1과 ERβ2는 각각 ‑11.8 kcal/mol, ‑11.6 kcal/mol 로 거의 동일한 친화도를 나타냈다. 이는 리간드 포켓 내 주요 상호작용(수소결합, π‑π 스택, 소수성 접촉)의 지속 시간과 빈도를 RMSF와 거리 분석을 통해 확인하였다. 특히 ERα2에서는 몇몇 핵심 잔기(예: His524, Glu353)의 움직임이 크게 증가해 결합 안정성이 저하된 것으로 해석된다.

두 번째 접근법은 계통유전학적 dN/dS 비율 분석이다. 각 ER 유전의 전체 코딩 서열과 LBD, DNA‑결합 도메인(DBD) 구역을 별도로 정렬한 뒤, PAML의 branch‑site 모델을 적용해 라인별 선택압을 추정하였다. 결과는 ERα2 LBD가 다른 ERα 동형체에 비해 dN/dS ≈ 0.68 로 상승했으며, 이는 ‘완화된 선택(relaxed selection)’을 의미한다. 반면 ERα1, ERβ1, ERβ2는 dN/dS ≈ 0.22‑0.30 으로 강한 정화 선택(purifying selection)을 유지한다. DBD는 모든 동형체에서 dN/dS < 0.1 로 매우 보존되어 있어, 전사인자와 DNA 결합 기능은 변이가 거의 없음을 보여준다. 이러한 결과는 ERα2가 기존 에스트로젠 결합 기능에서 탈피해 새로운 리간드 혹은 조절 메커니즘을 획득하고 있을 가능성을 뒷받침한다.

전체적으로, MD 기반 결합 친화도 차이와 진화유전학적 선택압 차이가 일치한다는 점이 주목할 만하다. 이는 구조적 변이가 실제 기능적 변화를 초래하고, 이는 장기적인 진화 과정에서 선택압의 변화를 반영한다는 가설을 강화한다. 또한, 최근의 연어류 전유전체 중복이 ERα2의 기능적 재배치를 촉진했을 가능성을 제시한다.