Integrative In Silico Analysis of Synonymous Variants in ABCG2 and PPARA: Implications for Valproate Resistance in Epilepsy

기술적 요약: ABCG2 및 PPARA 내 동의 치환 변이(Synonymous Variants)에 대한 통합적 인 실리코(In Silico) 분석

문제 제기
뇌전증 환자의 약 3분의 1은 항경련제, 특히 발프로산(VPA)에 대해 내성을 보입니다. ABCG2 수송체와 그 전사 조절 인자인 PPARA는 VPA 내성과 관련이 있는데, 구체적으로는 VPA가 PPARA를 상향 조절하고, 이것이 다시 ABCG2 발현을 증가시켜 뇌 내 약물 생체 이용률을 제한하는 경로를 통해 작용합니다. 그러나 이러한 과정에서 동의 단일 염기 다형성(sSNP)이 기여하는 바는 여전히 명확히 밝혀지지 않았습니다. 역사적으로 '침묵'하며 기능적으로 중립적인 것으로 간주되어 온 sSNP는, 최근 pre-mRNA 스플라이싱 변화, mRNA 이차 구조 변화, 번역 역학 조절 등의 메커니즘을 통해 유전자 발현을 교란할 수 있는 능력이 있는 것으로 인식되고 있습니다. 본 연구는 ABCG2 및 PPARA의 sSNP가 어떻게 기계적으로 개인 간의 VPA 반응 차이에 기여할 수 있는지에 대한 이해의 공백을 다룹니다.

방법론
저자들은 dbSNP 데이터베이스(GRCh38 참조)에서 추출한 ABCG2의 262개 및 PPARA의 109개 sSNP에 대해 통합적인 in silico 분석을 수행했습니다. 워크플로우는 다음 네 가지 조절 계층에 걸쳐 기능적 영향을 평가했습니다:

1. 스플라이싱 영향 (Splicing Impact): 변이가 스플계 부위(acceptor/donor)의 변화(수용체/공여체 획득 또는 손실)를 예측하도록 SpliceAI(심층 신경망)를 사용하여 분석하였으며, 높은 신뢰도 임계값($\Delta \ge 0.5$)과 민감도 임계값($\Delta \ge 0.2$)을 적용했습니다. 유망한 후보들은 Human Splicing Finder(HSF)를 통해 엑손 스플라이싱 인핸서(ESE), 실런서(ESS) 및 암호적 스플라이스 부위(cryptic splice sites)에 미치는 영향을 평가하여 추가 검증되었습니다.
2. RNA 이차 구조 (RNA Secondary Structure): RNAsnp(긴 전사체를 위한 Mode 2 및 국소 창을 위한 Mode 1)를 사용하여 국소적 구조 변화를 예측했습니다. 유의성은 p-값($\le 0.1$)과 최소 자유 에너지 변화($\Delta\Delta G$)로 결정되었습니다.
3. 번역 역학 (Translational Kinetics):
   • 코돈 사용 편향 (Codon Usage Bias): 상대적 동의 코돈 사용 비율(RSCU)을 계산하여 코돈 선호도의 변화를 식별했습니다.
   • 번역 효율성 (Translational Efficiency): GtRNAdb의 tRNA 유전자 복제수(tGCN)와 효모 유래 선택적 제약 파라미터($s_{ij}$)를 적용한 tRNA 적응 지수(tAI) 프레임워크를 사용하여 디코딩 효율성($\Delta W_i$)의 변화를 추정했습니다.
   • 디코딩 분류 (Decoding Classification): 변이를 Watson-Crick, wobble, 그리고 혼합 디코딩 메커니즘 간의 전이로 분류했습니다.
4. 공번역 접힘 (Co-Translational Folding): 단백질 접힘에 영향을 미칠 수 있는 것으로 알려진 민감한 영역, 즉 초기 N-말단 영역(코돈 120)과 도메인 경계로부터 2050 코돈 하류 지역을 매핑했습니다.
5. 통계 분석: 코돈 사용 편향 변화와 번역 효율성 사이의 관계를 평가하기 위해 피어슨 상관 분석을 수행했습니다.

주요 결과
분석 결과, 두 유전자 모두에서 기능적 결과가 예측되는 sSNP 하위 집단이 확인되었습니다:

• 스플라이싱 교란: SpliceAI는 ABCG2에서 16개, PPARA에서 3개의 유의미한 sSNP를 식별했습니다. 주목할 만한 변이는 ABCG2 c.1134A>G(공여체 획득/손실 이중 발생) 및 c.786A>T(암호적 공여체 생성 및 심각한 ESE/ESS 교란)입니다. PPARA의 경우, c.333C>T 및 c.126C>T가 암호적 공여체 부위를 생성할 것으로 예측되었습니다. HSF 분석은 이러한 변이들이 종종 조절 모티프를 교란하여 엑손 건너뛰기(exon skipping)나 프레임시프트(frameshift)를 유발할 수 있음을 확인했습니다.
• mRNA 이차 구조: ABCG2에서는 20개(7.6%)의 변이가 유의미한 구조적 교란($p \le 0.1$)을 보였습니다. 가장 불안정화시킨 변이는 c.480G>A($\Delta\Delta G = +6.6$ kcal/mol)였으며, c.1044A>G는 유의미한 안정화($\Delta\Delta G = -4.0$ kcal/mol)를 보였습니다. PPARA에서는 4개의 변이가 유의성 임계값을 충족했으며, c.1323G>A는 240 뉴클레오타이드에 걸친 장거리 헬릭스 상호작용을 방해하는 전역적 구조 재배열을 일으켰습니다.
• 번역 효율성 및 코돈 편향:
  • RSCU: ABCG2에서 29개 변이는 선호도가 증가(비율 >2)했고, 56개는 선호도가 감소(비율 <0.5)했습니다. PPARA에서도 유사한 경향이 관찰되었습니다.
  • 효율성 ($\Delta W_i$): ABCG2의 경우, 45.8%의 변이가 번역 효율성을 감소시킨 반면, 40.8%는 증가시켰습니다. PPARA에서는 51.3%가 효율성을 감소시켰습니다.
  • 상관관계: 두 유전자 모두에서 코돈 사용 편향 변화($\Delta$RSCU)와 번역 효율성 변화($\Delta W_i$) 사이에 중간 정도의 양의 상관관계가 발견되었습니다 (ABCG2 $r=0.478$, PPARA $r=0.567$).
• 공번역 접힘: 민감한 접힘 영역에서 변이가 확인되었습니다. ABCG2의 경우 초기 N-말단 영역에서 11개, 도메인 경계 하류에서 25개가 발견되었으며, PPARA는 각각 11개와 8개가 발견되었습니다.

의의 및 주장
본 논문은 ABCG2 및 PPARA의 동의 변이가 기능적으로 중립적인 것이 아니라, 여러 분자적 메커니즘을 통해 발프로산 내성에 영향을 미칠 수 있다고 주장합니다:

1. 스플라이싱 충실도 (Splicing Fidelity): 스플라이싱 조절 요소의 교란은 비정상적인 전사체 이소형(isoform)을 유발하여, 잠재적으로 비기능적 수송체나 수용체를 생성할 수 있습니다.
2. mRNA 안정성: 구조적 불안정화는 전사체의 분해를 가속화하여 단백질 양을 줄일 수 있는 반면, 안정화는 발현 수준을 변화시킬 수 있습니다.
3. 번역 역학: 코돈 사용 및 디코딩 효율성의 변화는 리보솜의 신장 속도를 변화시켜, 공번역 단백질 접힘을 방해하는 리보솜 정체(ribosome pausing)를 유발할 수 있습니다.

저자들은 PPARA의 sSNP가 VPA에 반응하여 ABCG2를 상향 조절하는 능력을 저해하거나, ABCG2의 sSNP가 직접적으로 수송체 기능이나 안정성을 변화시킬 수 있다는 기계적 프레임워크를 제안합니다. 이러한 섭동(perturbations)은 혈액-뇌 장벽에서의 발프로산 약동학에 영향을 미치는 방향으로 수렴됩니다.

결론
본 연구는 ABCG2 및 PPARA의 생물학적으로 의미 있는 일련의 동의 변이들이 발프로산 내성에 기여하기에 충분한 예측된 섭동을 나타낸다고 결론짓습니다. 저자들은 약물 유전체학 스크리닝의 정밀도를 높이기 위해 sSNP 프로파일링을 체계적으로 통합할 것을 권고합니다. 이들은 이러한 발견이 아미노산 보존성만을 근거로 중립성을 가정하는 대신, 구조적, 스플라이싱 및 번역 차원을 모두 고려하는 변이별 평가 프레임워크를 지지한다는 점을 강조합니다. 또한, 이러한 결과가 유전자형-표현형 관계를 확증하기 위한 향후 실험적 검증과 임상 연관 연구를 필요로 한다는 점을 명시하고 있습니다.