DFT-based conformational analysis of negatively-charged cytosine-containing artificial Threose Nucleic Acid (TNA) nucleotides

기술 요약: 음전하를 띤 사이토신 함유 인공 TNA 뉴클레오타이드의 DFT 기반 형태 분석

문제 제기
트레오스 핵산(Threose Nucleic Acid, TNA)은 구조적 단순성(4탄소 L-트레오스 당)과 RNA 및 DNA와의 안정적인 하이브리드 형성 능력 덕분에 유력한 생명 기원 유전 고분자 후보로 꼽히지만, 그 형태적 지형(conformational landscape)에 대한 포괄적인 계산 연구는 부족한 실정이다. 기존의 실험 연구들은 평균화된 구조 데이터를 제공하지만, TNA 단량체(monomer)가 점유하는 형태적 평형의 복잡성을 포착하지 못한다. 구체적으로, 음전하(인산기 도입)와 사이토신 함유 TNA 단량체의 특정 치환 패턴이 푸라노스 고리(furanose ring) 및 엑소사이클릭(exocyclic) 그룹의 형태적 다양성에 어떻게 영향을 미치는지 이해할 필요가 있다. 본 연구는 이러한 음전하를 띤 TNA 단량체에 대한 이론적 이해의 공백을 메우기 위해, 이전에 연구된 중성 유도체들과 비교 분석한다.

방법론
저자들은 $\alpha$-L-트레오푸라노스 당, C1' 위치의 사이토신 염기, 그리고 C3' 위치의 메톡시인산기를 특징으로 하는 음전하를 띤 TNA 단량체에 대해 포괄적인 계산 연구를 수행하였다. 견고성을 확보하기 위해 다음과 같이 다양한 수준의 이론을 적용한 DFT(밀도 범함수 이론)를 사용하였다:

• 기상(Gas Phase): HF/6-3G* 수준에서 구조를 사전 최적화한 후, B3LYP/6-31+G**, B3LYP/6-311++G**, 및 M062X/6-31+G에서 전체 최적화를 수행하였다. DFT 결과를 벤치마킹하기 위해 MP2/aug-cc-pVDZ//B3LYP/6-311++G 수준에서 단일 지점 에너지(single-point energy) 계산을 수행하였다.
• 수용액상(Aqueous Phase): 생물학적 환경을 모델링하기 위해, 기상에서의 최솟값들을 M062X/6-31+G** 수준의 SMD 용매화 모델을 사용하여 재최적화하였다.
• 분석: 분자 내 상호작용, 특히 수소 결합을 특성화하기 위해 NBO(Natural Bond Orbital) 및 QTAIM(Quantum Theory of Atoms in Molecules) 분석을 활용하였다. 푸라노스 고리의 형태(위상각 $P$ 및 퍼커링 진폭 $\nu_m$)를 정의하기 위해 Altona와 Sundaralingam(AS) 의 가상 회전(pseudorotation) 시스템을 사용하여 형태적 선호도를 분석하였다.

주요 기여 및 결과

1. 방법론적 평가:
   본 연구는 B3LYP와 M062X 범함수의 성능을 MP2 벤치마크와 비교하였다. 저자들은 M062X 범함수가 B3LYP에 비해 이러한 음이온 형태의 형태 에너지를 질적 및 양적으로 더 잘 설명한다는 것을 발견하였다. M062X 결과는 낮은 에너지 영역에서 MP2 데이터와 더 적은 산포도와 더 높은 상관관계를 보였으며, 반면 B3LYP는 대부분의 로타머(rotamer)에 대해 에너지를 과대평가하는 경향이 있었다.

2. 기상 형태 지형:

• B3LYP 결과: Conformer 3를 전역 최솟값(global minimum)으로 식별하였으며, 이는 평형 혼합물의 대부분(>95%)을 차지한다. 이 형태는 high-anti 염기 배향, ap 2'-OH 그룹, 그리고 -sc 인산 잔기를 특징으로 한다. 푸라노스 고리는 2T1 형태($P = 138.1^\circ$)를 취한다. 이는 인산기의 산소와 2'-OH 수소 사이의 상대적으로 강한 분자 내 수소 결합에 의해 안정화된다.
• M062X 결과: 더 복잡한 다성분 평형을 드러냈다. Conformer 1이 전역 최솟값으로 식별되었으며, 전체 혼합물의 50% 이상을 차지하며, Conformer 3는 약 27% 존재한다. Conformer 1은 anti 염기 배향, +sc 2'-OH, 그리고 ap 인산 잔기를 특징으로 한다. 푸라노스 고리는 가상 회전 바퀴의 북동쪽 영역(A-형 DNA의 특징)에 위치한 3E 형태($P = 14.1^\circ$)를 취한다.
• 안정화: NBO 및 QTAIM 분석은 Conformer 1이 2~3개의 상대적으로 약한 부적절한(improper) 수소 결합에 의해 안정화됨을 나타냈다. 이와 대조적으로, Conformer 3는 두 개의 상호작용, 그 중 하나는 부분적으로 공유 결합 성격을 띠는 강한 수소 결합에 의해 안정화된다.

3. 수용액 환경 효과:
   수용액 환경(SMD 모델)으로의 전이는 이전에 연구된 비전하 TNA 시스템에서는 관찰되지 않았던 급격한 형태적 안정성 변화를 일으켰다.

• Conformer 12가 수용액에서의 전역 최솟값으로 등장하였으며, 평형 분포의 **65.14%**를 차지하였다.
• 주목할 점은, Conformer 12는 기상에서의 Conformer 1을 안정화했던 분자 내 수소 결합에 의해 안정화되지 않는다는 것이다. 이러한 변화는 인산 잔기의 회전(기상 Conformer 1의 $\gamma$ 비틀림 각 ~165°에서 Conformer 12의 ~-68°로의 변화)에 의해 구동되며, 이는 분자 내 상호작용을 깨뜨리는 대신 용매화를 선호하게 만든다.
• 이전에 지배적이었던 기상 형태들(1 및 3)은 수용액에서 그 분포가 크게 감소하였다(각각 <5% 및 ~2.6%).

의의
본 논문은 음전하의 존재가 중성 대응체와 비교하여 TNA 단량체의 형태적 선호도를 유의미하게 변화시킨다고 주장한다. 본 연구는 다음을 강조한다:

• 음이온성 핵산 유사체의 형태적 평형을 정확하게 예측하기 위해서는 계산 범함수의 선택(M062X vs. B3LYP)이 매우 중요하다.
• 용매화 효과는 매우 강력하여, 기상에서 발견된 형태들의 안정성 순서를 완전히 역전시킬 수 있다.
• 수용액 환경에서 가장 안정적인 형태(Conformer 12)는 기상 최솟값에서 보이는 분자 내 수소 결합 네트워크가 결여되어 있으며, 이는 용매 상호작용이 생물학적 맥락에서 TNA의 구조를 결정하는 데 지배적인 역할을 한다는 것을 시사한다.

이러한 발견은 생명의 기원 및 생명공학(예: 압타머 및 안티센스 약물)에 대한 TNA의 응용 가능성을 규명하는 데 필수적인 구조적 행동에 대한 이론적 토대를 제공한다. 본 연구는 인공 핵산의 구조-활성 관계를 이해하기 위해 상세한 형태적 선호도에 대한 지식이 선행 조건임을 강조한다.