Computationally guided modifications of CviUPO to improve catalytic activity

기술 요약: 촉매 활성 개선을 위한 CviUPO의 계산 유도 변형 연구

문제 제기
비특이적 퍼옥시게네이스(UPO)는 수소 퍼옥사이드($H_2O_2$)만을 공기질(co-substrate)로 사용하여 포화 탄화수소를 선택적으로 산화할 수 있는 유망한 바이오 촉매이다. 그러나 이들의 산업적 응용은 퍼옥사이드에 의한 효소 불활성화에 의해 저해되는데, 이는 활성 포켓 내의 아미노산을 산화시키는 하이드록실 라디칼에 의해 유발되는 것으로 보인다. 따라서 이러한 문제를 완화하고 성능을 향상시키기 위해 표적 효소 공학이 필수적이지만, 핵심 활성 중심 잔기가 촉매 주기에 미치는 구체적인 영향은 여전히 잘 알려져 있지 않다. 본 연구는 Collariella virescens UPO(CviUPO)의 주요 활성 중심 잔기 변형이 에너지 장벽과 촉매 효율에 어떻게 영향을 미치는지 규명하고자 하는 필요성을 다룬다.

방법론
저자들은 CviUPO의 두 가지 특정 변형을 조사하기 위해 분자 동역학(MD) 시뮬레이션과 양자 역학/분자 역학(QM/MM) 계산을 결ิด합한 시너지 효과가 있는 계산 워크플로우를 채택하였다:

1. Glu162 $\to$ Asp162: 프로톤 추상화(proton abstraction)에 관여하는 글루탐산을 더 짧은 아스파르트산으로 교체하여, 포켓 부피의 증가가 활성을 개선하는지 테스트함.
2. Cys19 $\to$ His19: 헤임(heme) 고정용 시스테인을 히스티딘으로 교체하여, 헤임 표백(heme bleaching)에 대한 저항성 또는 퍼옥시다제와 유사한 거동으로의 전환 가능성을 조사함.

본 연구에서는 다음의 계산 기술들을 활용하였다:

• 고전적 MD 시뮬레이션: 야생형(native), Asp 변형형, His 변형형 CviUPO에 대해 120 ns 시뮬레이션을 수행하여 구조적 안정성, 활성 포켓 부피 및 수화 역학을 평가하였다. 기질(ethylbenzene/ETBE)의 안정성은 도킹 후 100 ns MD 시뮬레이션을 통해 테스트되었다.
• QM/MM 시뮬레이션: ORCA 소프트웨어를 사용하였으며, B3LYP 범함수(def2-TZVP 기저 집합)와 CHARMM36 힘의 장(force field)을 사용하여 화합물 I(반응성 철-옥소 복합체) 형성까지의 촉매 주기를 모델링하였다. 여기에는 저스핀(LS), 중간 스핀(IS), 고스핀(HS) 상태에 대한 기하 구조 최적화 및 진동수 분석이 포함되었다.
• Nudged Elastic Band (NEB) 시뮬레이션: CI-NEB 계산을 수행하여 최소 에너지 경로를 매핑하고 주요 중간체 간(초기 상태 $\to$ Compound 0 $\to$ Compound I)의 에너지 장벽을 식별하였다.
• 완화된 에너지 표면 스캔(Relaxed Energy Surface Scans): 활성 포켓 내에 기질이 존재하는 상태에서 야생형 및 His 변형 효소를 대상으로 수행되어, 촉매 주기의 두 번째 단계(Compound I $\to$ Compound II $\to$ Product)를 조사하였다.

주요 결과

• 동적 거동 및 수화: 세 가지 효소 변이체 모두 120 ns 동안 안정성을 유지하였다. Asp 및 His 변형은 야생형 효소에 비해 활성 포켓 부피가 약간 작아졌으나, 포켓 내 평균 물 분자 밀도는 모든 변이체에서 유사하게 나타났다. 특히, His 변형 효소는 포켓 부피의 시간적 변화가 더 크게 나타났으며, 이는 잠재적인 구조적 불안정성을 시사한다.
• Compound I 형성 (야생형 vs. 변형형):
  • Asp 변형: NEB 시뮬레이션 결과, 초기 상태에서 Compound 0로 가는 에너지 장벽이 현저히 증가(약 3 eV)한 것으로 나타났다. 이는 더 짧은 아스파르트산 잔기로 인해 $H_2O_2$로부터 프로톤이 이동해야 하는 거리가 멀어져 탈양성자화(deprotonation)를 방해하기 때문인 것으로 분석된다. 결과적으로, Asp 변형은 촉매 활성에 해로운 것으로 예측된다.
  • His 변형: NEB 결과, 초기 상태와 Compound 0 사이에 에너지 장벽이 없는 것으로 나타나 페릭 하이드로페록소(ferric hydroperoxo) 복합체의 형성이 용이함을 시사했다. 그러나 Compound 0와 Compound I 사이의 장벽은 야생형보다 높았다(1.25–1.63 eV). 또한, 최종 상태(Compound I)가 Compound 0보다 에너지적으로 낮지 않아, 촉매 주기가 국소 최솟값(local minimum)에 갇힐 가능성을 제기했다.
• 스핀 상태: 변형된 효소들의 경우, 계산된 LS, IS, HS 상태의 에너지가 거의 동일하게 나타났으며, 이는 국부적인 전자 환경의 작은 변화가 선호되는 전자 배치에 큰 영향을 미친다는 것을 의미한다. 반면, 야생형 효소는 스핀 상태 간의 구분이 더 뚜렷했으며, 일반적으로 IS 상태가 가장 낮은 에너지 값을 보였다.
• 기질 상호작용: MD 시뮬레이션은 ETBE가 야생형 및 His 변형 효소의 활성 포켓 내에서 안정적이지만, Asp 변형 변이체에서는 불안정함을 확인하였다.
• 기질과의 반응 경로: His 변형 효소에 대한 완화된 표면 스캔 결과, Compound II 형성을 위한 장벽이 야생형(1.37 eV)에 비해 더 낮음(0.67 eV)을 보여주었다. 그러나 이 스캔은 하이드록실 그룹의 직접적인 리바운드(rebound) 대신 물 분자가 생성물 형성을 촉진할 수 있다는 메커니즘적 편차를 시사했다. 이는 His 변형이 효소의 메커니즘을 퍼옥시게네이스에서 퍼옥시다제로 변화시킬 가능성을 제기하지만, 이것이 확정적으로 확인되지는 않았다.

의의 및 주장
본 논문은 계산 시뮬레이션이 CviUPO 촉매 주기의 근본적인 과정을 성공적으로 규명하고 활성 포켓의 수화 및 스핀 상태의 결정적 역할을 강조했음에도 불구하고, 이러한 변형이 전체적인 촉매 활성을 증가시키는지 결정하는 데에는 불충분하다고 결론지었다.

• Asp 변형은 프로톤 전달과 관련된 높은 에너지 장벽으로 인해 계산적으로 불리한 것으로 예측된다.
• His 변형은 복잡한 트레이드오프(trade-off)를 보여준다: 페릭 하이드로페록소 복합체의 초기 형성은 용이하게 하지만, 주기 후반부에 더 높은 장벽을 도입하며 반응 메커니즘을 완전히 바꿀 가능성이 있다.
  저자들은 시뮬레이션만으로는 기질 특이성을 확정하거나 His 변형이 효소를 퍼옥시다제로 전환하는지 여부를 확정할 수 없음을 강조하며, 효율적인 효소 공학을 위해서는 이론적 예측과 실험적 검증을 연결하는 시너지 접근 방식이 필요하다고 주장한다. His 변형 효소의 발현을 위한 예비 실험적 시도는 낮은 효소량과 낮은 헤임 로딩(heme loading)으로 인해 어려움을 겪었으며, 이는 실험적 검증의 필요성을 더욱 뒷받침한다.

본 연구는 고수준의 양자 화학 계산이 분자 동역학과 통합될 때 전자적 상호작용과 에너지 지형에 대한 깊은 통찰을 제공할 수 있음을 보여주는 동시에, 순수 이론적 예측의 한계를 인정하며 향다의 합리적 설계(rational design)를 가이드하는 데 기여한다.