기술 요약: 아릴말로네이트 탈카복실화 효소의 조상들이 보여주는 높은 활성, 안정성 및 입체 선택성

문제 정의

Bordetella bronchiseptica 유래 박테리아 아릴말로네이트 탈카복실화 효소(BbAMDase)는 $\alpha$-카이랄 카르복실산(비스테로이드성 항염증제(NSAIDs) 및 카이랄 빌딩 블록 포함)을 생성하는 데 귀중한 바이오 촉매이다. 그러나 이 효소의 광범어적인 산업적 적용은 세 가지 주요 제한 사항으로 인해 저해된다:

1. 낮은 안정성: 이 효소는 열적 및 운전 안정성이 낮아, 30 °C에서 반감기($t_{1/2}$)가 1.2 시간에 불과하며 40 °C 이상의 온도에서 상당한 활성 손실을 보인다.
2. 제한된 기질 범위: 효소는 활성 부위의 입체적 제약으로 인해 $\alpha$-에틸 치환 말로네이트를 처리하는 데 어려움을 겪으며, 더 작은 치환기에 대해서는 활성이 감소한다.
3. 최적화되지 않은 입체 선택성: 일반적으로 선택성을 보이지만, 입체적으로 유사한 치환기(예: 2-메틸-2-비닐말로닉 산)를 가진 기질에 대해서는 거울상 이성질체 과잉률(ee)이 감소하며, (R)에서 (S)로 입체 선택성을 반전시키기 위해 광범위한 공학적 설계가 필요하다.

부위 특이적 돌연변이 도입을 통한 전통적인 단백질 공학은 안정성 또는 선택성 향상을 위한 돌연변이를 도입할 때 촉매 효율을 저해하는 경우가 많다. 저자들은 조상 서열 재구성(Ancestral Sequence Reconstruction, ASR)이 이러한 특성들을 회복하거나 강화할 수 있는 적절한 서열 맥락을 제공할 수 있다는 가설을 세웠다.

방법론

본 연구는 계통 발생학, 계산 생물학 및 생화학적 특성 분석을 결wise 결합한 다단계 접근 방식을 채택하였다:

1. 계통 발생 분석 및 ASR:

   • 222개의 AMDase 서열과 4개의 외집단(outgroup) 서열(cis-trans 말레산 이성질화 효소, MI) 데이터셋을 정렬하였다.
   • WAG 모델을 사용하여 최대 가능도(Maximum Likelihood, ML) 계통수를 구축하였다.
   • 10개의 조상 노드(N5, N10, N31, N87, N107, N109, N131, N151, N164, N200)를 선택하여 서로 다른 클레이드와 계층적 수준을 포괄하도록 하였다.
   • 갭(gap)을 별개의 상태로 취급하여 phangorn을 이용한 주변부 재구성(marginal reconstruction) 방식으로 조상 서열을 재구성하였다.

2. 생화학적 특성 분석:

   • 조상 유전자를 E. coli BL21 (DE3)에서 발현시켰다.
   • 안정성: 차등 주사 형광 분석법(Differential Scanning Fluorimetry, DSF)을 통해 열 안정성을 평가하여 녹는 온도($T_m$)를 결정하였다. 운전 안정성(반감기)은 30 °C 및 45 °C에서 측정하였다.
   • 활성: 말론산 패널(1a–6a, 프로카이랄 및 $\alpha$-에틸 기질 포함)에 대한 비특이적 활성을 측정하였다. 유망한 변이체들에 대해 운동학적 파라미터($K_M$, $k_{cat}$)를 계산하였다.
   • 입체 선택성: (R)-선택적 조상들과 (S)-선택적 변이체들(촉매 시스테인 G74C/C188G 교체 및 소수성 포켓 돌연변이 도입을 통해 생성됨)에 대해 거울상 이성질체 과잉률(ee)을 카이랄 GC로 분석하였다.

3. 계산 분석:

   • 서열 분석: 서열의 열 안정성, 용해도 및 표면 전하 특성을 예측하기 위해 ProtParam, SCMTPP, Protein-Sol 도구를 사용하였다.
   • 분자 역학(MD): 야생형 BbAMDase와 가장 안정적인 조상인 N131에 대해 300 K 및 360 K에서 10 $\times$ 0.5 $\mu$s MD 시뮬레이션을 수행하였다. 용매 접근 가능 표면적(SASA), 수화 껍질, 수소 결합 및 염교(salt bridge) 분석에 집중하였다.
   • 메커니즘 연구: 탈카복실화의 입체 화학적 경로를 추적하기 위해 동위원소 표지된 4a (pro-R-13C)를 사용하였다.

4. 준비 규모 적용:

   • N131을 사용하여 2-메틸-2-페닐 말로네이트(2a)의 그램(gram) 규모 탈카복실화를 수행하여 시공간 수율(space-time yield)과 총 전환수(total turnover number, TTN)를 평가하였다.

주요 결과

1. 향상된 안정성 및 활성

• 열 안정성: 조상 N131은 62 °C의 $T_m$을 나타냈다(야생형보다 7 °C 높음). 더욱 놀랍게도, 45 °C에서 N131은 BbAMDase(1.3 시간)에 비해 294배 증가한 반감기(381 시간)를 보였다. 조상 N5 또한 72배의 개선을 보였다.
• 촉매 효율: N131은 2-메틸-2-페닐 말로네이트(2a)에 대해 야생형 BbAMDase(25.1 U·mg⁻¹ vs. 14.3 U·mg⁻¹)보다 1.8배 높은 비활성을 나타냈다. 운동학적 분석 결과, N131은 이 기질에 대해 현저히 낮은 $K_M$(0.3 mM vs. 25.5 mM)과 더 높은 $k_{cat}$(33.9 s⁻¹ vs. 29 s⁻¹)을 가졌다.
• 기질 범위: N131은 입체적으로 까다로운 $\alpha$-에틸 기질(3a)에 대해 완전한 전환을 달면서, BbAMDase는 부분적 전환(15%)만을 달성했다.

2. 개선된 입체 선택성

• 자연적 선택성: 조상 N31은 2-메틸-2-비닐말로닉 산으로부터 (R)-4b를 99.7% ee로 생성하였으며, 이는 야생형의 96.8% ee보다 유의미하게 향상된 수치이다.
• 역전된 선택성: N5의 엔지니어링된 (S)-선택적 변이체(특히 N5 ICPLLG)는 (S)-2b에 대해 >99.8% ee를 달성하여, 엔지니어링된 BbAMDase ICPLLG(89.0% ee)를 능가하였다.
• 메커니즘: 동위원소 표지는 N31의 높은 선택성이 pro-S 탈카복실화 경로의 억제로부터 기인함을 확인하였으며, 이는 전이 상태 간의 에너지 차이를 야생형(2.5 kcal·mol⁻¹)에 비해 3.9 kcal·mol⁻¹로 증가시켰다.

3. 안정성의 구조적 근거 (N131)

계산 분석은 N131의 탁월한 안정성을 야생형에 비해 6% 증가한 표면 글루타메이트(Glu) 잔기 덕분이라고 설명하였다.

• 수화 껍질: MD 시뮬레이션은 이러한 표면 Glu 잔기들이 단백질 표면을 보호하는 보호적이고 질서 정연한 수화 껍질을 형성함을 보여주었다.
• 염교(Salt Bridges): 고온(360 K)에서 N131은 더 많은 분자 내 수소 결합을 유지하고 열 유도 염교를 형성한 반면, 야생형은 이러한 상호작용을 상실하였다.
• 용해도: 증가된 표면 극성은 소수성 SASA 비율을 감소시켜 열적 응집을 방지하였다.

4. 준비 규모 성능

그램 규모 반응에서 N131은 기질 2a를 31분 만에 거의 완결 수준으로 전환하였다(야생형의 느린 전환 대비).

• 생산성: N131은 엔자임 1g당 시간당 96.2g의 생성물이라는 특정 생산성을 달성하였으며, 이는 BbAMDase보다 3.4배 높은 수치이다.
• 총 전환수 (TTN): N131은 $3.0 \times 10^8$의 TTN을 달축하였으며, 이는 BbAMDase 대비 633배 향상된 결과이다.

의의 및 주장

본 논문은 조상 서열 재구성이 전통적인 지향적 진화(directed evolution)로는 달성하기 어려운 자연적으로 진화된 바이오 촉매의 개선을 발견하는 강력한 전략임을 주장한다.

• 안정성과 활성의 디커플링: 본 연구는 조상 효소가 활성을 희생하지 않고도 더 높은 열 안정성과 더 높은 촉매 활성을 동시에 보유할 수 있음을 입증하여, 안정성이 종종 활성을 저해한다는 통념에 도전한다.
• 메커니즘적 통찰: 수화 차폐를 통한 열 안정성의 핵심 동력으로서 표면 Glu 잔기를 식별한 것은 단백질 공학을 위한 새로운 설계 원칙을 제공한다.
• 실용적 유용성: N131 조상은 더 높은 TTN과 빠른 반응 시간을 통해 생산 비용과 환경 발자국을 크게 줄일 수 있는, 광학적으로 순수한 $\alpha$-카이랄 카르복실산의 산업적 합성을 위한 우수한 바이오 촉매로 제시된다.
• 공학적 가이드: 계통 발생 분석은 이전에 합리적 설계(예: 소수성 포켓 내 돌연변이)를 통해 확인된 핵심 돌연변이들이 진화 계통도상에서 자연적으로 발생함을 보여주었으며, 이는 ASR이 안정성 및 선택성 과제에 대한 "자연적인" 해결책을 강조함으로써 향후 공학 캠페인을 안내할 수 있음을 시사한다.

저자들은 모든 재구성된 조상이 개선을 보이는 것은 아니지만, N131과 N31의 성공적인 식별은 ASR이 효소 패밀리의 촉매 다양성과 견고성을 확장하는 강력한 방법임을 검증한다고 결론지었다.