DIELECTRIC CONTINUUM MODELING FOR THE DE NOVO DESIGNED CALCIUM CHANNEL CALC6.3: STRUCTURE PREPARATION, MESHGENERATION, AND CALCIUM SELECTIVITY

기술 요약: 데 노보(De Novo) 설계된 칼슘 채널 CalC6 3에 대한 유전체 연속체 모델링

문제 정의
본 연구는 유전체 연속체 모델을 사용하여 이온 채널의 칼슘 선택성을 계산적으로 모델링하는 데 따르는 과제를 다룬다. 고해외 전압 의존성 칼슘(CaV) 채널의 고해상도 구조는 존재하지만, 현재 열린 상태(open state)의 실험적으로 결정된 구조는 부재하여 이온 투과에 대한 직접적인 시뮬레이션을 제한하고 있다. 템플릿 기반 모델링이 존재하기는 하나, 열린 상태에 대한 실험적 검증이 부족하다. 반면, CalC6 3와 같은 데 노보 설계된 채널은 실험적으로 해상된 열린 구멍 구조(PDB ID: 9DZW)와 입증된 칼슘 선택성을 제공한다. 그러나 비균질 크기 수정 포아송-볼츠만(nuSMPBic) 모델을 CalC6 3에 적용하는 데에는 다음과 같은 중대한 기술적 장애물이 존재한다:

1. 구조 준비: 원래의 PDB 파일은 계산적으로 생성된 수소 원자를 포함하고 있으며 175개의 무거운 원자(heavy atoms)가 누락되어 있어, 정제 및 수리가 필요하다.
2. 메시 생성: 채널이 좁고 복잡한 선택성 필터를 가지고 있기 때문에 고품질의 삼각형 분자 표면 메시를 생성하는 것이 어렵다. 표준적인 표면 생성 방법들은 외부 표면의 매끄러움을 달성하기 위해 기공의 기하학적 구조를 왜곡하며, 이는 선택성 필터의 물리적 정확성을 저해한다.
3. 워크플로 통합: 원시 PDB 데이터로부터 정전기 및 이온 농도 방정식의 수치 해를 구하는 과정까지의 워크플로를 자동화하는 통합 소프트웨어 패키지가 필요하다.

방법론
저자들은 CalC6 3에 대한 전체 시뮬레이션 워크플로를 간소화하는 통합 소프트웨어 패키지를 개발하였다. 방법론은 네 가지 단계로 진행된다:

1. 구조 준비: 원본 9DZW.pdb 파일을 처리하기 위한 파이썬 스크립트가 작성되었다. 여기에는 계산적으로 생성된 수소 원자와 CRYST1 레코드를 제거하고, 구조를 변환하여 기공 축이 z축에 정렬되도록 하며, 채널을 원점에 배치하는 과정이 포함되었다. 수리된 구조는 PDB2PQR을 통해 PQR 파일로 변환되었다. 이 단계에서는 누락된 175개의 무거운 원자를 재구성하고, pH 7.0에서의 양성자화 상태를 할당하며, 수소 결합을 최적화하고, AMBER 포스 필드에 기반하여 원자 전하와 반경을 할당하였다.
2. 표면 메시 생성: 삼각형 분자 표면 메시를 생성하기 위해 GICMSG(Gaussian Ion Channel Molecular Surface Generation) 패키지가 사용되었다. 단일 감쇠 파라미터를 사용하는 기존 방식(예: TMSmesh)과 달리, GICMSG는 두 개의 공간적으로 변화하는 파라미터 함수($d_p$ 및 $d_m$)를 사용하여 외부 표면의 매끄러움과 채널 기공의 기하학적 구조를 독립적으로 제어한다. 두 가지 케이스가 테스트되었다: 기본 파라미터를 사용한 Case 1과, 외부 표면은 더 매끄럽게 만들면서 기공의 기하학적 구조는 보존하도록 튜닝된 파라미터를 사용한 Case 2이다.
3. 부피 메시 구축: 생성된 표면 메시는 단백질, 막, 용매 영역을 포함하는 시뮬레이션 박스($\Omega$)를 위한 인터페이스 적합(interface-fitted) 비정형 사면체 메시를 구축하는 데 사용되었다. 박스 크기는 채널의 반데르발스 표면을 포함하도록 설정되었다.
4. 수치 해 도출: 통합 패키지는 메시 생성을 nuSMPBic 유한 요소 솔버와 결합하였다. 수치 실험은 0.1 M NaCl과 0.1 M CaCl2(벌크 농도: 0.1 M Ca²⁺, 0.1 M Na⁺, 0.3 M Cl⁻)를 포함하는 이온 용액에 대해 수행되었다. 수화 반경은 BioNumbers 데이터베이스를 기반으로 할당되었다. 솔버는 비선형 유한 요소 시스템을 해결함으로써 정전기 전위와 이온 농도 분포를 계산하였다.

주요 기여

• 통합 소프트웨어 패키지: 본 논문은 원시 PDB 데이터(9DZW)로부터 정전기 전위 및 이온 농도의 시각화에 이르는 워크플로를 자동화하는 통합 도구를 제시한다. 이는 수동 개입을 줄이고 데 노보 설계된 채널의 연속체 모델링 준비 과정을 표준화한다.
• 기공 보존 메시 생성: GICMSG의 적용을 통해, 유한 요소 계산을 위한 충분한 외부 표면의 매끄러움을 유지하면서도 복잡한 좁은 선택성 필터의 기하학적 구조를 보존하는 고품질의 삼각형 메시를 성공적으로 생성하였다.
• 연속체 모델의 검증: nuSMPBic 모델을 실험적으로 해상된 CalC6 3 구조에 적용했을 때, 기하학적 구조나 전하 분포에 대한 임의적인 단순화 없이도 실험적으로 관찰된 칼슘 선택성을 재현할 수 있음을 입증하였다.

결과

• 칼슘 선택성: 수치 실험 결과, Na⁺가 더 작은 수화 반경을 가짐에도 불구하고, 여섯 개의 Glu 잔기로 형성된 좁은 선택성 필터 내에 Ca²⁺가 우선적으로 축적됨을 확인하였다. Ca²⁺ 농도는 필터의 세포 내 입구 근처에서 약 0.8 mol/L에 도달하였으며, 이는 벌크 농도인 0.1 mol/L를 크게 상회하는 수치이다. Na⁺는 기공으로부터 대부분 배제되었다.
• 메시 매끄러움 및 솔버 성능: Case 1(거친 외부 표면)과 Case 2(매끄러운 외부 표면)를 비교했을 때, 두 케이스 모두 기공의 기하학적 구조를 보존하고 유사한 이온 농도 프로파일을 산출하였으나, Case 2가 솔버 성능을 크게 향상시키는 것으로 나타났다. Case 2는 Case 1(151회 반복)에 비해 72회의 반복만으로 수렴하였으며, 총 계산 시간은 단일 CPU 코어(Mac Studio, M4 Max)에서 약 49초에서 약 44초로 감소하였다. 이는 외부 분자 표면을 매끄럽게 하는 것이 중요한 기공의 물리 법칙을 왜곡하지 않으면서도 유한 요소 솔버의 수렴성을 높인다는 것을 시사한다.
• 이온 분포: 시뮬레이션은 뚜렷한 분포 패턴을 보여주었다: Ca²⁺는 기공과 세포 내 입구 근처에 축적되었고, Na⁺는 배제되었으며, Cl⁻는 이온 전하, 크기 및 채널의 정전기적 환경의 결합된 효과를 반영하여 주로 세포 외 입구 근처와 벌크 용매에 축적되었다.

의의
본 논문은 CalC6 3가 칼슘 선택성을 조사하기 위한 가치 있는 모델 시스템이라고 주장한다. 이는 전압 감지 영역이 없는 데 비해 필수적인 선택성 필터는 유지하고 있어, 천연 CaV 채널보다 단순하면서도 실험적으로 결정된 열린 구멍 구조를 제공하기 때문이다. 본 연구는 유전체 연속체 모델, 특히 nuSMPBic 접근법이 현실적인 구조적 및 유전체 환경 내에서 칼슘 선택성의 평형 특징을 포착할 수 있음을 보여준다. 통합 패키지는 (~142,000개의 미지수를 ~44초 만에 해결하는) 계산 효율적인 도구를 제공하며, 생리적 이온 농도가 낮은 거대 분자 역학(coarse-grained MD) 처리를 위한 기초가 될 수 있다. 저자들은 이 패키지가 현재 CalC6 3에 적용되고 있으나, 다른 이온 채널 단백질을 시뮬레이션하도록 적응 가능하도록 설계되었다고 언급하였다.