Spatial O-GalNAc glycoproteomics atlas: region-specific regulation across the mouse brain

기술 요약: 생쥐 뇌의 공간적 O-GalNAc 당단백질체 지도

문제 정의
뮤신형 O-GalNAc 당쇄화(glycosylation)는 신경 세포 간의 통신, 부착 및 신호 전달을 조절하는 중요한 번역 후 변형이다. 이 과정을 개시하는 20개의 GALNT(GalNAc 전이효소)는 부위 특이적 제어를 제공하지만, 포유류 뇌 전체에 걸친 체계적이고 영역별로 해상도가 높은 정량적 O-GalNAc 당쇄화 지도는 그동안 부재했다. Allen Brain Atlas와 같은 기존 자원들은 전사체학에 집중되어 있어 단백질 당쇄화를 포착하지 못한다. 또한, 분석적 난제들로 인해 희소하고 해부학적으로 정의된 뇌 영역에 당단백질체학을 적용하는 데 어려움이 있었다. 특히, O-GalNAc 당쇄의 불균질성과 Ser/Thr이 풍부한 영역에서의 클러스터링은 표준 고에너지 충돌 해리(HCD) 파편화 하에서 위치 국소화(localization)를 복잡하게 만드는데, HCD에서는 당쇄 모이어티(moiety)가 우선적으로 절단되기 때문이다. 이는 데이터 의존적 획득(DDA) 스펙트럼 라이브러리에서 "위치적 중복성(positional redundancy)"을 초래하며, 단일 전구체 신호가 당쇄 부위 국소화만 다른 상호 배타적인 여러 글리코펩타이드 이성질체에 매칭되게 한다. 이러한 중복성은 검색 공간을 팽창시키고, 통계적 점수 산출을 복잡하게 하며, 다수의 샘플에 대해 일관된 정량을 위해 선호되는 데이터 독립적 획득(DIA) 워크플로우의 확장성을 제한한다.

방법론: Glyco-DIA 2.0
이러한 한계를 해결하기 위해 저자들은 정량적 O-GalNAc 당단백질체학을 위해 설계된 비중복적, DIA 중심 워크플로우인 Glyco-DIA 2.0을 개발하였다.

• 비중복적 스펙트럼 라이브러리 설계: 핵심 혁신은 당쇄 조성과 특정 글리코사이트 국소화를 분리하는 라이브러리 구축 전략이다. 당쇄를 개별 세린(serine) 또는 트레오닌(threonine) 잔기에 할당하는 대신(이는 HCD 하에서 이성질체 중복성을 생성함), 이 워크플로우는 누적 당쇄 질량을 펩타이드 백본의 가상 N-말단 변형으로 취급한다. 이는 "당쇄 차감(glycan-subtraction)" 검색 전략(subHCD)을 통해 달성되는데, 여기서 데이터베이스 검색 전 전구체 질량에서 누적 당쇄 질량을 차감한다. 이를 통해 O-글리코펩타이드 MS/MS 스펙트럼을 실질적으로 비당쇄화된 펩타이드 스펙트럼처럼 해석할 수 있게 하여, 스펙트럼 중복성을 크게 줄이고 식별률을 높인다.
• 샘플 준비: 생후 21일(P21) C57BL/6j 생쥐의 뇌를 10개의 거대 해부학적 영역(후각구, 대뇌 피질, 해마, 시상, 시상하부, 중뇌, 교뇌, 연수, 소뇌, 맥락총)으로 미세 절단하였다. 샘플은 말단 불균질성을 줄이고 렉틴 결합을 강화하기 위해 탈시알산화(desialylation) 과정을 거쳤으며, 이후 트립신 소화 및 Jacalin 기반 약한 친화성 크로마토그래피를 통한 T- 및 Tn-타입 에피토프 농축을 수행하였다.
• 질량 분석:
  • 라이브러리 생성: 긴 기울기의 nanoLC-MS/MS를 DDA 모드로 사용하여 전체 뇌 및 풀링된 구획 샘플로부터 심층 참조 라이브러리를 생성하였다. 전체 뇌 샘플은 고-pH 역상 크로마토그래피를 통해 추가로 분획되었다.
  • 정량화: 9개 뇌 구획(입력 차이로 인해 전체 뇌와 맥락총은 차별 분석에서 제외)에 대해 DIA 질량 분석법을 사용하여 공간적 매핑을 수행하였다.
  • 검증: 구조적 주석을 위한 명확한 글리코사이트 지도를 생성하기 위해 별도의 ETD 기반 워크플로우를 사용하였다.
• 데이터 분석: 데이터 분석에는 Proteome Discoverer DDA 출력을 Spectronaut 호환 비중복 라이브러리로 변환하는 자동화된 R 기반 파이프라인이 사용된다. 정량 데이터는 견고한 통계 분석을 보장하기 위해 분산 안정화 정규화(vsn), 배치 교정(ComBat), 그리고 랜덤 포레스트 기반의 결측치 임퓨테이션(imputation) 과정을 거쳤다.

주요 결과

• 심층 참조 라이브 라이브러리: 본 연구는 13,007개의 O-GalNAc 글리코펩타이드 전구체와 9,769개의 고유 글리코펩타이드 서열을 포함하는 포괄적인 참조 라이브러리를 생성하였으며, 이는 817개의 당단백질로부터 유래되었다. 이는 최근의 데이터셋과 비교했을 때 알려진 생쥐 뇌 O-GalNAc 당단백질체를 520개의 당단백질만큼 확장한 것이다.
• 공간적 지도: Glyco-DIA 2.0을 9개 뇌 영역에 적용하여 534개의 서로 다른 당단백질에서 유래된 4,008개의 고신뢰도 고유 글리코펩타이드 데이터셋을 확보하였다.
• 당쇄 조성: 데이터셋은 Mono-glycosylated Hex(1)HexNAc(1) (T 에피토프, 47.2%) 및 HexNAc (Tn 항원, 8.7%)에 의해 지배되었다. 고차 조성은 단일 복합 당쇄보다는 여러 부위에 걸친 짧은 core-1 구조의 누적 점유율을 나타내는 것으로 결정되었다.
• 영역 특이적 조절: 주성분 분석(PCA)은 명확한 영역 의존적 조직화를 보여주었으며, 시상하부와 후각구의 뚜렷한 클러스터링, 그리고 교뇌, 연수, 소뇌의 관련 클러스터링을 나타냈다.
• 조절의 디커플링(Decoupling): 정량적 단백질체학과의 통합을 통해 당쇄화 패턴이 단백질 풍부도에 의해서만 결정되지 않음을 밝혔다. 연구는 세 가지 조절 상태를 식별하였다:
  • 일치 조절 (Concordant Regulation, CR): 글리코펩타이드와 단백질 풍부도가 함께 변화함.
  • 글리코펩타이드 특이적 조절 (Glycopeptide-Specific Regulation, GSR): 단백질 수준과 독립적으로 당쇄화 점유율 변화가 발생함.
  • 불일치 조절 (Discordant Regulation, DR): 글리코펩타이드와 단백질 풍부도가 반대 방향으로 변화함.
  • VGF, IGFBP5, CX3CL1, CD44, GFRA2, APLP1을 포함한 많은 단백질이 GSR 또는 DR을 나타냈으며, 이는 당쇄화 점유율 수준에서의 복잡한 조절을 시사한다.
• 단백질 내 리모델링: 데이터는 개별 단백질 내의 공간적 리모델링을 드러냈다. 예를 들어, CX3CL1은 중뇌/후각구의 막 근접 영역과 달리 교뇌/연수에서는 N-말단 스템(stalk)에서 우선적인 당쇄화를 보였으며, 이는 영역 특이적 글리코프로포름(glycoproteoform)의 생성을 시사한다.

의의 및 주장
저자들은 본 연구가 생쥐 뇌에 대한 최초의 공간적으로 해상된 O-GalNAc 당쇄화 점유 지도를 구축했다고 주장한다. DIA 워크플로우의 위치적 중복성이라는 분석적 병목 현상을 극복함으로써, Glyco-DIA 2.0은 제한된 샘플 입력량에서도 정량적 조직 특이적 당단백질체학을 위한 확장 가능한 프레임워크를 제공한다.

본 연구는 O-GalNAc 당쇄화가 단백질 발현 수준과는 구별되며 종종 독립적인, 포유류 뇌에서의 추가적인 공간적 분자 조절 계층을 구성함을 입증한다. 이 조절은 단백질 접힘, 수송, 프로테아제 민감도 및 수용체 클러스터링에 영향을 미친다. 저자들은 이러한 영역 특이적 글리코프로포름이 시냅스 조직, 축삭 유도 및 마이엘린 생물학에 차별적인 영향을 미칠 수 있다고 가정한다.

커뮤니티 접근 및 추가 연구를 촉진하기 위해, 저자들은 인터랙티브 웹 리소스인 **GlycoDB (glycodb.org)**를 공개하였다. 이 플랫폼을 통해 사용자는 공간적 당쇄화 패턴을 탐색하고, 영역별 풍부도를 비교하며, 의사 히트맵(pseudoheatmap)을 통해 단백질 내 리모델링을 시각화하고, 구조 데이터(AlphaFold 모델)를 당쇄화 부위와 통합하여 확인할 수 있다. 저자들은 이 리소스가 중추 신경계에서 O-GalNAc 당쇄화의 역할을 조사하기 위한 토대를 제공한다고 강조하면서도, 현재의 지도가 단일 발달 단계(P21) 및 야생형 생쥐로 제한되어 있음을 언급하였다.