Mutation-induced conformational remodeling explains gain- and loss-of-function phenotypes in Euglena gracilis phytoene synthase

기술 요약: Euglena gracilis 파이토엔 합성효소의 돌연변이 유도 구조적 재구성

문제 제기
파이토엔 합성효소(CrtB)는 카로티노이드 생합성의 속도 제한 단계인 첫 번째 단계를 촉매하며, 이는 대사 공학의 핵심적인 표적이 된다. Euglena gracilis에서의 실험적 스크리닝을 통해 특정 돌연변이(Y197A, E121H, L349F)가 상반된 기능 획득(카로티노이드 생산 증가) 및 기능 상실(생산 감소) 표현형을 유발한다는 사실이 확인되었으나, 이러한 차이를 유도하는 구조적 메커니즘은 여전히 불분명하다. 정적인 구조 모델과 실험적 관찰만으로는 단일 잔기 치환이 어떻게 효소 활성을 변화시키는지 설명하기에 부족하다. 본 연구는 특정 돌연변이가 변동된 카로티노이드 생산량과 연결되는 동적 구조적 근거에 대한 이해의 공백을 다룬다.

연구 방법
저자들은 야생형 효소와 세 가지 변이체(Y197A, E121H, L349F)의 구조적 역동성을 조사하기 위해 E. gracilis CrtB의 아포(apo, 기질 결합 전) 형태에 대한 분자 동역학(MD) 시뮬레이션을 수행하였다. 계산적 접근 방식에는 다음이 포함된다:

• 에너지 분석: 돌연변이에 따른 에너지적 섭동 및 안정성($\Delta\Delta G$)을 평가하기 위한 MM-GBSA 및 MM-PBSA 계산.
• 정전기적 특성 규명: 등전점(pI)의 변화와 pH 7.0에서의 순전하를 결정하기 위한 PROPKA 분석.
• 전역 구조 지표: 전역적 조밀함과 형태적 샘플링을 평가하기 위한 용매 접근 가능 표면적(SASA) 및 회전 반경($R_g$) 분석.
• 국소적 역동성: 잔기별 가동성 변화를 매핑하기 위한 제곱평균제곱근 변동(RMSF) 및 $\Delta$RMSF 분석.
• 활성 부위 구조: Asp-풍부 모티프(Mg²⁺ 배위에 필수적) 및 소수성 플랩(hydrophobic flap) 내의 용매 접근성, 결합 포켓 부피 및 수소 결합 네트워크에 대한 상세 조사.

주요 결과
시뮬레이션 결과, 돌연변이들은 전역적인 구조 변화를 크게 일으키지 않으면서도 국소적 상호작용 네트워크와 활성 부위 구조에 영향을 미치는 독특한 형태적 에너지 지형의 재구성을 유도함을 밝혔다:

1. 에너지 및 정전기적 변화:

   • E121H는 가장 큰 형태적 역동성 재분배를 보였으며, 양(+)의 $\Delta\Delta G$ 값은 야생형 단백질 대비 에너지적 불안정성을 나타냈다. 또한 가장 높은 예측 pI와 중성 pH에서 약간의 양전하를 나타냈다.
   • Y197A는 중립에 가깝거나 유리한 에너지 변화를 보였다.
   • L349F는 완만한 에너지 섭동을 유도했으나, PROPKA 폴딩 자유 에너지 예측에서 뚜렷한 불안정성을 보였다.

2. 형태적 샘플링 및 유연성:

   • 기능 획득 변이체 (Y197A, E121H): 두 변이 모두 상대적으로 접근 가능한 소수성 플랩을 유지했다. E121H는 가장 넓은 $R_g$ 분포를 보였으며, Y197A와 E121H는 야생형보다 넓은 SASA 분포를 보여, 전체적인 조밀함은 유지하면서도 향상된 형태적 샘플링을 시사했다.
   • 기능 상실 변이체 (L349F): 전역적 SASA 및 $R_g$ 분포에서 야생형과 매우 유사했으나, 소수성 플랩 포켓의 급격하고 국소적인 수축(약 35% 부피 감소)을 보였다.

3. 국소적 상호작용 네트워크 및 활성 부위 재구성:

   • 소수성 플랩: L349F 돌연변이는 플랩 포켓 부피와 용매 접근성의 유의미한 감소를 초래했으며, 이는 더 닫힌 형태와 일치한다. 반면, 기능 획득 변이체들은 더 높은 플랩 노출도를 유지했다.
   • Asp-풍부 모티프: 모든 돌연변이는 Asp-풍부 1 영역에서 수소 결합 점유율 증가를 보였다. 특히, L349F는 낮은 촉매 성능에도 불구하고 가장 높은 점유율을 보였는데, 이는 촉매 효율이 최대의 국소적 안정화보다는 안정성과 유연성 사이의 균형에 의존함을 시사한다.
   • L349F의 경직화: L349F 돌연변이의 구조 분석 결과, Tyr127 잔기가 더 경직된 형태를 취하며 Arg154 및 Glu158과의 $\pi$-양이온 및 $\pi$-음이온 상호작용을 포함하는 지속적인 상호작용 네트워크를 형성함을 발견했다. 이는 야생형 및 기능 획득 변이체에서 Tyr127이 형태적으로 유연했던 것과 대조적이다.

의의 및 주장
본 논문은 돌연변이에 의한 촉매 중심 주변의 형태적 재구성이 실험적으로 관찰된 카로티노이드 생산량 차이와 연결된다는 구조적 프레임워크를 제공한다고 주장한다. 저자들은 상반된 표현형이 단순히 돌연변이 지점에서의 국소적 구조 섭동 때문이 아니라, 기능적으로 중요한 모티프(특히 소수성 플랩 및 인접한 Asp-풍부 1 모티프) 주변의 형태적 샘플링 재분배에서 기인한다고 단언한다.

연구는 L349F의 기능 상실 표현형이 소수성 플랩의 "두드러진 수축"과 "경직된 상호작용 네트워크"로 특징지어지는 반면, 기능 획득 변이체들은 "상대적으로 접근 가능한 소수성 플랩"과 더 큰 유연성을 유지한다는 점을 강조한다. 이러한 결과는 동적 구조 조직이 CrtB 기능의 결정적인 요소임을 시사한다. 저자들은 이러한 통찰이 파이토엔 합성효소의 합리적 공학을 위한 기초를 제공한다고 결론지었으나, 이러한 고유한 형태적 변화가 기질 인식 및 촉매 작용에 어떻게 영향을 미치는지 완전히 규명하기 위해서는 기질(GGPP)과 Mg²⁺ 이온을 포함한 향후 시뮬레이션이 필요함을 명시적으로 언급하였다.