1. 서론

1.1 탐구 동기

계성고등학교의 3학기제는 정규 학기 외의 시기에 한 과목을 더 깊게 수강할 수 있도록 마련되어 있어, 평소 교과서에서 짧게 지나갔던 개념을 더 좁고 분명한 문제의식으로 발전시키기에 적절한 환경을 제공한다. 이 제도를 통해 대학 수준의 분자생물학을 수강하면서 가장 크게 달라진 점은 유전 정보의 흐름을 바라보는 시선이었다. 고등학교 생명과학에서는 DNA에서 RNA를 거쳐 단백질이 만들어지는 과정을 비교적 선명한 도식으로 배우지만, 실제 진핵세포에서는 전사 직후 생성된 pre-mRNA가 그대로 번역되는 것이 아니라 여러 단계의 가공을 거쳐 최종적인 발현 산물이 정해진다. 그중에서도 스플라이싱은 인트론을 제거하고 엑손을 연결하는 기본 과정에 그치지 않고, 어떤 엑손을 남기고 어떤 엑손을 제외할 것인지까지 조절하는 핵심 단계라는 점에서 인상적이었다.

특히 같은 유전자로부터 서로 다른 단백질 산물이 만들어지는 대안적 스플라이싱은, 유전자 수보다 훨씬 큰 단백질 다양성이 어떻게 가능해지는지를 설명하는 중요한 원리이다. 처음에는 스플라이싱을 단순한 RNA 정리 과정으로 이해했지만, 관련 내용을 더 읽어 갈수록 실제 결정은 훨씬 복잡하다는 사실이 드러났다. 동일한 pre-mRNA라도 엑손 내부와 인트론 내부의 조절 서열, splice site의 강도, 스플라이소솜 구성요소의 접근성, 그리고 RNA 결합 단백질의 경쟁 관계에 따라 전혀 다른 결과가 나타날 수 있었다. 그 과정에서 가장 자주 등장한 조절 인자가 SR 단백질과 hnRNP였다. SR 단백질은 대체로 스플라이소솜 조립을 돕고 엑손 포함을 촉진하는 방향으로 설명되며, hnRNP는 splice site 인식을 방해하거나 silencer 서열에 결합해 엑손 배제를 유도하는 방향으로 설명된다.

하지만 문헌을 더 자세히 살펴보면 이 둘을 단순히 촉진자와 억제자로 나누는 설명만으로는 충분하지 않았다. 실제 스플라이싱 결과는 결합 위치와 주변 서열 문맥, 약한 splice site의 존재, 세포 유형, 전사 속도에 따라 달라질 수 있고, 같은 단백질도 상황에 따라 다른 방향으로 작동하기 때문이다. 따라서 스플라이싱 조절을 제대로 이해하려면 단백질 이름을 나열하는 수준을 넘어서, 어떤 조건에서 어떤 인자가 스플라이소솜의 인식 과정을 강화하거나 약화시키는지까지 추적해야 한다고 판단하였다. 이러한 문제의식은 생명과학Ⅱ에서 배운 유전자 발현과 RNA 가공 개념을 대학 수준의 분자 조절 메커니즘으로 확장하는 데 적절하다고 보였다.

이에 본 탐구는 SR 단백질과 hnRNP의 길항 작용이 약한 splice site를 가진 cassette exon의 포함과 배제를 어떻게 결정하는지에 초점을 두고, 문헌 분석과 가상 데이터 모델링을 결합한 방식으로 대안적 스플라이싱 조절의 구조를 정량적이고 기전적으로 해석하고자 하였다. 이 주제는 전사 후 조절이 유전자 발현 다양성의 핵심 층위임을 보여줄 뿐 아니라, 비정상적 스플라이싱이 암과 신경계 질환 같은 질병의 분자적 기반이 될 수 있다는 점에서도 탐구 가치가 크다고 보았다.

1.2 탐구 목적

본 탐구의 목적은 진핵세포의 대안적 스플라이싱 과정에서 SR 단백질과 hnRNP가 pre-mRNA의 조절 서열에 결합하여 엑손 포함과 배제를 어떻게 조절하는지 분석하고, 이러한 길항 작용이 스플라이소솜 조립과 splice site 인식에 미치는 영향을 분자생물학적으로 해석하는 데 있다. 이를 위해 먼저 스플라이싱의 기본 원리와 exon definition의 개념을 정리하고, 이후 SR 단백질과 hnRNP의 기능을 비교하며, 마지막으로 가상 데이터 모델을 통해 두 인자의 상대적 활성 비율이 엑손 포함 비율에 미치는 정량적 경향을 분석하였다.

보다 구체적으로는 첫째, SR 단백질과 hnRNP가 각각 어떤 RNA 서열 요소에 결합하며 어떤 방식으로 spliceosome recruitment를 변화시키는지 밝히고자 하였다. 둘째, 강한 splice site와 약한 splice site의 차이를 바탕으로 왜 특정 엑손에서 조절 효과가 더 크게 나타나는지 해석하고자 하였다. 셋째, SR/hnRNP 활성 비율이라는 가상의 독립변수를 설정하여 엑손 포함 비율 PSI의 변화를 수치적으로 모델링함으로써, 스플라이싱 결과가 단순 선형 반응이 아니라 임계 구간을 가지는 비선형 조절 과정일 가능성을 검토하고자 하였다.

1.3 연구 질문

본 탐구의 중심 질문은 다음과 같이 정리할 수 있다. 약한 splice site를 가진 cassette exon에서 SR 단백질과 hnRNP의 상대적 활성이 바뀔 때 엑손 포함 비율은 어떤 방식으로 달라지는가. 또한 그 변화는 단순한 비례 증가인지, 아니면 특정 균형점 부근에서 급격히 달라지는 비선형 조절 양상인지가 핵심 질문이다. 이를 통해 스플라이싱이 자동적 절단 과정이 아니라 RNA 위에서 벌어지는 경쟁적 해석 과정이라는 점을 확인하고자 하였다.

2. 이론적 배경

2.1 전사 후 조절과 대안적 스플라이싱의 의미

진핵세포에서 유전자는 전사된 직후 바로 단백질로 번역되지 않는다. 먼저 pre-mRNA 상태에서 5' cap 형성, polyadenylation, 스플라이싱과 같은 가공 과정을 거친 뒤에야 성숙 mRNA가 만들어진다. 그중 스플라이싱은 인트론을 제거하고 엑손을 연결하는 과정이며, 대안적 스플라이싱은 동일한 pre-mRNA로부터 서로 다른 조합의 엑손을 포함한 여러 mRNA isoform을 생성하게 한다. 이 과정은 하나의 유전자가 하나의 단백질만을 만든다는 단순한 도식을 수정하게 만들며, 조직 특이성이나 발생 단계에 따라 서로 다른 단백질 기능이 나타나는 이유를 설명하는 데 중요한 역할을 한다.

따라서 대안적 스플라이싱은 단순한 RNA 후처리 단계가 아니다. 오히려 DNA에 기록된 정보를 어떤 형태의 단백질 산물로 실현할 것인지를 RNA 수준에서 다시 결정하는 핵심 조절 층위에 가깝다. 같은 유전자를 가진 세포라도 스플라이싱 방식이 달라지면 최종 단백질 구조와 기능이 달라질 수 있으며, 이 점은 유전자 발현 조절을 전사 단계에만 한정해서는 충분히 이해할 수 없음을 보여준다.

2.2 스플라이소솜과 기본 스플라이싱 과정

스플라이싱은 스플라이소솜이라는 대형 리보핵단백질 복합체에 의해 수행된다. 스플라이소솜은 U1, U2, U4, U5, U6 snRNP를 핵심 구성요소로 가지며, pre-mRNA의 특정 부위를 순차적으로 인식하여 반응성 복합체를 형성한다. 일반적으로 5' splice site, branch point, polypyrimidine tract, 3' splice site가 정확히 인식되어야 인트론이 제거될 수 있다. 먼저 U1 snRNP가 5' splice site에 결합하고, 3' 쪽에서는 U2AF와 SF1이 polypyrimidine tract 및 branch point 인근을 인식한다. 이후 U2 snRNP가 branch point에 결합하고, 뒤이어 U4/U6.U5 삼중 복합체가 합류하면서 활성화 단계를 거쳐 촉매 반응이 진행된다.

이 반응은 두 번의 전이 반응으로 요약된다. 첫 번째 반응에서는 branch point의 아데노신이 5' splice site를 공격하여 lariat 구조를 형성하고, 두 번째 반응에서는 상류 엑손의 3' 말단이 3' splice site를 공격하여 두 엑손이 연결된다. 결과적으로 인트론은 고리 모양으로 제거되고 성숙 mRNA가 생성된다. 이처럼 화학적 절단과 연결은 일정한 원리를 따르지만, 실제 생물학적 변동성의 대부분은 어느 splice site가 얼마나 안정적으로 인식되느냐에서 비롯된다. 그래서 스플라이싱은 촉매 반응 자체보다 인식과 선택의 문제로 접근해야 한다.

2.3 강한 splice site와 약한 splice site

모든 엑손이 동일한 확률로 선택되는 것은 아니다. splice site의 염기서열이 합의 서열과 얼마나 유사한지에 따라 스플라이소솜의 기본 인식 효율은 크게 달라진다. 일반적으로 5' splice site나 3' splice site가 보존성이 높아 U1 snRNP나 U2AF가 안정적으로 결합할 수 있는 경우를 강한 splice site라고 부른다. 이런 경우에는 외부 조절 인자의 개입이 상대적으로 적더라도 엑손이 성숙 mRNA에 포함될 가능성이 높다.

반면 약한 splice site는 합의 서열과의 일치도가 낮거나 주변 맥락이 불리하여 기본적인 spliceosome recognition만으로는 안정적인 스플라이싱이 어렵다. 대안적 스플라이싱이 빈번하게 일어나는 cassette exon은 약한 splice site를 가지는 경우가 많고, 이때 엑손 포함 여부는 조절 인자의 존재에 매우 민감해진다. 본 탐구가 약한 splice site를 가진 reporter exon을 가정한 것도 이러한 이유 때문이다. 강한 splice site에서는 SR/hnRNP 비율 변화의 효과가 상대적으로 제한될 수 있지만, 약한 splice site에서는 작은 조절 변화가 PSI의 큰 변동으로 이어질 가능성이 높다.

2.4 ESE, ESS, ISE, ISS와 RNA 서열 기반 조절

pre-mRNA는 단순히 단백질 서열 정보를 담는 분자가 아니다. 동시에 스플라이싱 조절을 위한 서열 신호도 포함하고 있다. 엑손 내부에는 ESE와 ESS가, 인트론 내부에는 ISE와 ISS가 존재할 수 있으며, 이러한 서열 요소는 스플라이싱 촉진 또는 억제 인자의 결합 지점으로 기능한다. ESE는 대체로 SR 단백질의 결합을 유도하여 인접 splice site의 인식을 강화하고, ESS는 hnRNP 계열 단백질의 결합을 유도하여 해당 엑손의 포함을 억제하는 경향이 있다. ISE와 ISS 역시 유사한 방식으로 인트론 내부에서 조절 효과를 낸다.

그러나 이 분류는 작동 경향을 설명할 뿐 절대 법칙은 아니다. 동일한 단백질이라도 결합 위치가 달라지면 결과가 달라지고, 주변에 어떤 다른 서열 요소가 존재하는지에 따라서도 효과가 변한다. 따라서 대안적 스플라이싱은 단일 신호의 결과라기보다 여러 조절 신호가 겹쳐지는 문맥 의존적 산물이라고 이해하는 편이 정확하다.

2.5 SR 단백질의 구조와 기능

SR 단백질은 세린과 아르지닌이 풍부한 도메인을 가지는 RNA 결합 단백질군으로, 대안적 스플라이싱 조절에서 가장 대표적인 촉진 인자로 알려져 있다. 이들은 하나 또는 둘 이상의 RNA recognition motif를 통해 pre-mRNA의 특정 서열에 결합하고, RS domain을 통해 다른 스플라이싱 인자와 단백질-단백질 상호작용을 수행한다. 이러한 구조적 특징 덕분에 SR 단백질은 단순 결합 인자를 넘어서 RNA와 스플라이소솜 구성요소 사이를 연결하는 분자적 가교 역할을 수행할 수 있다.

기능적으로 SR 단백질은 주로 ESE에 결합하여 약한 5' splice site의 U1 snRNP 인식을 강화하거나 3' splice site 쪽의 U2AF 결합을 안정화한다. 이 과정은 결국 엑손 양쪽 경계를 하나의 단위로 인식하게 만드는 exon definition을 촉진한다. 특히 약한 splice site를 가진 cassette exon에서는 기본 결합의 불안정성을 보완하는 역할이 중요하므로, SR 단백질의 농도 변화나 인산화 상태 변화가 엑손 포함률을 크게 바꿀 수 있다.

2.6 hnRNP의 구조와 기능

hnRNP는 heterogeneous nuclear ribonucleoprotein의 약자로, 핵 내 pre-mRNA와 결합하여 다양한 전사 후 조절에 관여하는 단백질군이다. 이들은 구조적으로 RNA 결합 도메인을 가지며, 구성원에 따라 RRM, KH domain, glycine-rich domain 등 서로 다른 부위를 가진다. 스플라이싱 조절에서 hnRNP는 주로 ESS 또는 ISS에 결합하여 스플라이소솜의 접근을 억제하거나 RNA의 2차 구조를 변화시켜 특정 splice site를 물리적으로 덜 노출시키는 방향으로 작용한다.

대표적으로 hnRNP A1은 특정 엑손 주변에 결합하여 exon skipping을 촉진하는 사례가 많이 알려져 있다. 다만 hnRNP 또한 항상 억제자로만 작용하는 것은 아니며, 결합 위치와 주변 인자의 조합에 따라 다른 억제 인자의 접근을 차단하거나 RNA 구조를 유리하게 바꾸어 오히려 엑손 포함에 기여하는 경우도 보고되어 있다. 따라서 hnRNP의 본질은 단순한 반대항이 아니라, 스플라이싱 결과를 문맥 의존적으로 재배열하는 조절자에 가깝다.

2.7 SR 단백질과 hnRNP의 길항 작용

대안적 스플라이싱의 핵심은 특정 엑손이 스플라이소솜에 의해 충분히 인식되는지 여부이며, 이 인식 과정은 RNA 위의 조절 단백질들에 의해 지속적으로 조정된다. SR 단백질이 enhancer 서열에 결합하여 U1 snRNP와 U2AF의 결합을 돕는다면, hnRNP는 silencer 서열에 결합하여 동일 부위 또는 인접 부위의 접근성을 떨어뜨릴 수 있다. 이 두 작용은 독립적이라기보다 같은 pre-mRNA 위에서 경쟁적으로 작동하며 최종 결과를 결정한다.

이 길항 작용은 단순한 양적 대결이 아니라 결합 친화도, 결합 위치, 단백질 간 상호작용, RNA 구조의 노출 정도, 전사 속도 같은 여러 요소가 얽힌 동적 과정이다. 그럼에도 상대적 활성 비율이라는 개념은 전체 경향을 정량적으로 표현하는 데 유용하다. 일반적으로 SR 단백질의 상대적 우세가 커질수록 엑손 포함이 증가하고, hnRNP의 상대적 우세가 커질수록 엑손 배제가 증가할 가능성이 높다고 볼 수 있다.

2.8 exon definition과 질병 연계성

고등 진핵생물의 유전자는 대체로 인트론이 길고 엑손이 짧은 구조를 가지므로, 스플라이소솜은 인트론 전체보다 엑손 경계를 중심으로 인식하는 exon definition 전략을 많이 사용한다. 이때 한 엑손의 양쪽 경계를 안정적으로 묶어 주는 조절 인자의 존재가 중요하며, SR 단백질은 이러한 묶임을 촉진하는 방향으로 작용한다. 반대로 hnRNP가 splice site 근처의 접근을 차단하거나 RNA 구조를 불리하게 바꾸면 exon definition이 실패하여 해당 엑손은 건너뛰어질 수 있다.

이러한 원리가 중요한 이유는 스플라이싱 오류가 곧 기능적 단백질 구성의 오류로 이어질 수 있기 때문이다. 실제로 암에서는 특정 SR 단백질의 과발현이나 hnRNP 발현 이상이 세포 증식이나 침습성과 관련된 isoform 생성을 유도하는 경우가 보고되어 있다. 또한 신경계 질환과 일부 유전질환에서도 exon skipping 또는 비정상적 inclusion이 병인 기전으로 작용한다. 따라서 스플라이싱 조절 인자의 기능을 이해하는 일은 단순한 RNA 가공 개념 정리를 넘어 질병의 분자적 기반을 해석하는 작업과 직접 연결된다.

3. 연구 방법

3.1 연구 설계

본 탐구는 문헌 분석과 가상 데이터 모델링을 결합한 이론 기반 탐구 연구로 설계하였다. 연구 질문은 SR 단백질과 hnRNP의 상대적 활성 변화가 약한 splice site를 가진 cassette exon의 포함 비율에 어떤 영향을 주는가로 설정하였다. 이에 따라 독립변수는 SR/hnRNP 활성 비율 지수로, 종속변수는 exon inclusion 비율을 나타내는 PSI 값으로 설정하였다. 이 설계는 스플라이싱 조절 메커니즘을 직접 실험으로 구현하지 않더라도, 기존에 알려진 기전적 경향을 정량적으로 재구성함으로써 탐구 문제를 더 분명하게 분석할 수 있다는 장점을 가진다.

연구 가설은 두 단계로 구성하였다. 첫째, SR 단백질의 상대적 활성이 증가하고 hnRNP의 상대적 활성이 감소할수록 약한 splice site를 가진 엑손의 PSI는 증가할 것이라고 가정하였다. 둘째, 이러한 증가는 단순 선형이 아니라 일정 임계 비율 근처에서 스플라이소솜 조립 효율이 급격히 증가하는 비선형 패턴을 보일 것이라고 가정하였다. 이 가설은 스플라이싱 조절이 특정 단백질의 존재 유무보다 상대적 균형과 결합 경쟁의 결과라는 이론적 배경에 근거한다.

3.2 자료 선정과 분석 대상

자료는 스플라이싱 조절과 관련된 리뷰 논문과 원저 논문을 중심으로 정리하였다. 분석 대상은 약한 5' splice site를 가진 cassette exon reporter 시스템으로 설정하였다. 이는 실제 연구에서 자주 사용되는 minigene reporter 개념을 단순화한 것이다. 모든 조건에서 pre-mRNA 총 발현량, 세포 종류, 배양 조건, 측정 시점은 동일하게 유지된 것으로 가정하였고, SR 단백질과 hnRNP의 상대적 활성만 변화시키는 방식으로 모델을 구성하였다. 이러한 설정은 독립변수와 종속변수의 관계를 더 명료하게 보기 위한 통제 조건이다.

3.3 가상 실험 조건과 변수

가상 실험 조건은 SR/hnRNP 비율을 0.50, 0.80, 1.00, 1.20, 1.50, 2.00의 여섯 수준으로 나누고, 각 수준마다 네 번의 반복 측정을 수행한 것으로 설정하였다. PSI 값은 비율이 증가할수록 상승하되, 중간 구간에서 민감하게 증가하고 높은 구간에서는 포화에 가까워지는 시그모이드형 경향을 반영하도록 구성하였다. 반복값 사이에는 현실성을 높이기 위해 소규모 오차를 부여하였고, 일부 조건에서는 약한 이상값을 포함시켜 생물학적 측정의 변동성을 반영하였다.

정량적 경향을 설명하기 위한 개념적 모델은 $PSI=20+\dfrac{65}{1+e^{-4(x-1.2)}}+\varepsilon$로 둘 수 있다. 여기서 $x$는 SR/hnRNP 활성 비율이고, $\varepsilon$는 측정 오차를 뜻한다. 이 식은 본 탐구의 실측값이 아니라 스플라이싱 조절이 임계 구간을 가질 수 있다는 점을 설명하기 위한 이론적 틀로 사용하였다.

3.4 분석 방법

각 조건의 평균 PSI와 표준편차를 계산하고, 전체 경향이 양의 상관을 보이는지 확인하였다. 또한 단순 직선적 증가가 아닌 비선형적 변화 여부를 검토하기 위해 낮은 구간, 중간 구간, 높은 구간의 증가 폭을 비교하였다. 이와 함께 중간 비율 영역에서 분산이 커지는지 살펴봄으로써, 스플라이싱 결정이 전이 구간에서 상대적으로 민감해질 가능성도 함께 해석하였다.

4. 결과

4.1 SR/hnRNP 비율에 따른 PSI 원시 데이터

가상 실험에서 SR/hnRNP 활성 비율을 여섯 수준으로 나누어 측정한 결과는 다음과 같이 정리되었다. 각 조건은 동일한 약한 splice site를 가진 cassette exon reporter를 대상으로 네 번 반복 측정한 것으로 설정하였다.

원시 데이터 수준에서 확인되는 가장 기본적인 패턴은 SR/hnRNP 비율이 증가함에 따라 PSI가 전반적으로 상승한다는 점이다. 가장 낮은 비율인 0.50 조건에서는 반복 측정값이 21.8%에서 25.1% 범위에 머물렀고, 가장 높은 비율인 2.00 조건에서는 74.6%에서 83.0% 범위로 상승하였다. 즉, hnRNP 우세 조건에서는 엑손 배제가 상대적으로 우세하고, SR 우세 조건에서는 엑손 포함이 현저히 증가하는 방향성이 일관되게 나타났다.

또한 반복 간 편차가 완전히 동일하지 않다는 점도 확인되었다. 예를 들어 1.00 조건의 48.2%는 같은 조건의 다른 반복에 비해 다소 높은 값이며, 2.00 조건의 74.6%는 같은 집단 평균보다 다소 낮다. 그러나 이러한 값들은 전체 경향을 뒤집는 수준은 아니며, 실제 RT-PCR 기반 정량이나 밴드 해석 과정에서 나타날 수 있는 실험적 변동 범위를 모사한 것으로 해석할 수 있다.

4.2 기초 통계량과 변화 폭

각 조건별 평균 PSI와 표준편차를 정리하면 다음과 같다.

평균 PSI는 23.6%, 32.1%, 42.1%, 54.1%, 67.3%, 79.3% 순으로 증가하였다. 최저 조건과 최고 조건을 비교하면 총 증가 폭은 $79.3-23.6=55.7$%p이다. 또한 최고 조건의 PSI는 최저 조건의 약 $\dfrac{79.3}{23.6}\approx3.36$배에 해당한다. 이 수치는 SR 단백질의 상대적 우세가 약한 splice site를 가진 엑손의 포함 가능성을 매우 크게 변화시킬 수 있음을 보여준다.

표준편차는 1.4에서 4.6 범위로 나타났으며, 1.00 조건에서 가장 크게 관찰되었다. 이는 중간 비율 영역에서 스플라이싱 결정이 상대적으로 불안정하거나 민감하게 반응할 수 있음을 시사한다. 반대로 0.50과 0.80 조건에서는 비교적 낮은 분산이 관찰되어, hnRNP 우세 구간에서는 엑손 배제가 안정적으로 유지될 가능성을 보여준다.

조건 간 변화 폭을 구체적으로 비교하면 0.50에서 0.80 구간은 8.5%p, 0.80에서 1.00 구간은 10.0%p, 1.00에서 1.20 구간은 12.0%p, 1.20에서 1.50 구간은 13.2%p, 1.50에서 2.00 구간은 12.0%p의 증가를 보였다. 특히 1.00에서 1.50 사이의 누적 증가 폭은 $67.3-42.1=25.2$%p로 크게 나타나, 이 구간이 스플라이싱 결과의 민감 구간일 가능성을 뒷받침한다.

4.3 시각화와 경향 분석

평균 PSI를 그래프로 나타내면 SR/hnRNP 비율이 증가함에 따라 엑손 포함 비율이 꾸준히 상승하는 흐름을 더 분명하게 확인할 수 있다.

그래프는 낮은 비율 구간에서는 완만한 상승이 나타나고, 1.00에서 1.50 부근에서는 증가 폭이 커지며, 높은 비율 구간에서는 상승세가 유지되지만 점차 포화에 가까워지는 양상을 보여준다. 이는 스플라이싱 조절이 단순 선형 반응보다는 시그모이드형 또는 완만한 비선형 증가 모델과 더 잘 맞을 수 있음을 시사한다.

4.4 비율 변화에 따른 패턴 정리

데이터를 구간별로 나누어 보면 세 가지 패턴이 확인된다. 첫째, 낮은 비율 구간인 0.50과 0.80에서는 PSI가 각각 23.6%와 32.1%로 낮게 유지되며 엑손 배제가 우세하다. 이는 hnRNP의 억제 작용이 상대적으로 강하여 스플라이소솜이 약한 splice site를 안정적으로 인식하지 못하는 상황과 부합한다. 둘째, 중간 비율 구간인 1.00과 1.20, 1.50에서는 PSI가 42.1%, 54.1%, 67.3%로 비교적 큰 폭으로 상승한다. 이는 SR 단백질과 hnRNP의 균형점 주변에서 작은 비율 변화가 exon definition의 성공률을 크게 바꿀 수 있음을 보여준다. 셋째, 높은 비율 구간인 2.00에서는 PSI가 79.3%까지 상승하지만 증가 속도는 무한정 가팔라지지 않고 일정 수준에서 포화에 접근하는 모습을 보인다.

이러한 패턴은 SR 단백질이 일정 수준 이하일 때는 inclusion이 낮은 상태로 유지되다가, 특정 임계 영역을 지나면 spliceosome recruitment가 급격히 강화되며 inclusion이 빠르게 증가하고, 이후에는 대부분의 전사산이 엑손 포함 isoform으로 전환되면서 증가 폭이 제한되는 구조로 해석할 수 있다. 따라서 본 결과는 스플라이싱 조절이 문턱값을 가지는 조절 시스템이라는 가설과 대체로 일치한다.

5. 논의 및 결론

5.1 수치 결과 해석과 이론 연결

가상 데이터 분석 결과 SR/hnRNP 활성 비율이 증가함에 따라 PSI가 일관되게 상승하였으며, 특히 1.00에서 1.50 사이의 중간 구간에서 증가 폭이 크게 나타났다. 이러한 결과는 SR 단백질이 enhancer 서열에 결합하여 U1 snRNP 및 U2AF의 인식을 강화하고, hnRNP가 silencer 서열에 결합하여 이를 방해한다는 기존 이론과 잘 부합한다. 즉, 대안적 스플라이싱의 결과는 특정 단백질 하나의 존재 여부가 아니라 같은 pre-mRNA 위에서 작동하는 상반된 조절 인자들의 상대적 우세에 의해 결정된다고 해석할 수 있다.

특히 약한 splice site를 가진 엑손을 연구 대상으로 설정했을 때 이러한 경향이 분명해졌다는 점은 중요하다. 강한 splice site는 기본 인식이 안정적이어서 조절 인자의 효과가 제한될 수 있지만, 약한 splice site는 기본적으로 인식 효율이 낮기 때문에 SR 단백질의 작은 증가만으로도 exon definition이 성공할 가능성이 커진다. 반대로 hnRNP가 상대적으로 우세하면 이미 취약한 인식 과정이 더 불안정해져 엑손 배제가 쉽게 일어난다. 따라서 약한 splice site는 대안적 스플라이싱 조절에서 변화 감도가 높은 핵심 분기점으로 이해할 수 있다.

또한 중간 비율 영역에서 표준편차가 상대적으로 크게 나타난 점은 스플라이싱 결정이 이 구간에서 불안정한 전이 상태를 가진다는 해석과 연결될 수 있다. 이는 생물학적 시스템에서 자주 관찰되는 임계 전이 현상과 유사하며, 일정 수준 이상의 촉진 신호가 누적되기 전까지는 결과가 쉽게 흔들릴 수 있음을 시사한다. 스플라이싱을 이처럼 확률적이면서도 조절 가능한 의사결정 과정으로 이해하면 RNA 가공을 훨씬 더 동적인 관점에서 해석할 수 있다.

5.2 비교 분석과 대안 해석

본 결과를 비교 분석하면 SR 우세 조건과 hnRNP 우세 조건의 차이는 단순 평균 차이를 넘어 스플라이싱 시스템의 작동 방식 차이를 보여준다. hnRNP 우세 구간에서는 PSI가 낮고 분산도 비교적 작아 스플라이싱 배제가 안정적으로 유지되는 구조가 나타났다. 반면 SR 우세 구간에서는 PSI가 높아지면서도 일부 변동이 남아 있어, 촉진 상태에서도 항상 완전한 일방향 결정이 이루어지는 것은 아님을 보여준다. 이는 실제 세포에서 스플라이싱이 언제나 100% 포함 또는 100% 배제로 수렴하지 않고, 일정한 isoform 혼합 상태를 유지하는 경우가 많다는 점과도 맞닿아 있다.

다만 본 탐구의 기본 모델은 SR 단백질을 촉진자, hnRNP를 억제자로 설정한 단순화된 구조를 기반으로 한다. 실제 문헌에서는 같은 SR 단백질이나 hnRNP라도 결합 위치와 주변 문맥에 따라 반대 방향의 효과를 나타낼 수 있음이 알려져 있다. 따라서 본 결과를 해석할 때는 SR/hnRNP 비율 자체가 절대적 원인이라고 단정하기보다, 조절 방향의 대표적 경향을 수치화한 모형으로 이해하는 것이 타당하다. 예를 들어 특정 인트론 위치에 결합한 hnRNP가 다른 억제 인자의 접근을 차단하여 오히려 inclusion에 기여하는 경우도 있을 수 있다. 이와 같은 대안 해석은 스플라이싱 조절이 본질적으로 문맥 의존적이라는 점을 다시 확인하게 한다.

5.3 방법론적 한계와 후속 탐구

본 탐구는 문헌 기반 이론 분석과 가상 데이터 모델링을 결합한 형태이므로, 실제 세포에서 측정한 실험값을 직접 사용하지 않았다는 한계를 가진다. 따라서 수치 결과는 실험적 사실의 재현이라기보다, 알려진 기전을 바탕으로 합리적 경향을 모형화한 것이다. 이 점은 탐구의 정량적 구조를 강화하는 장점이 있는 동시에, 실제 세포 환경의 복잡성을 충분히 반영하지 못한다는 제한도 남긴다.

또한 스플라이싱 결과는 SR 단백질과 hnRNP 외에도 RNA 2차 구조, 전사 속도, 크로마틴 상태, 다른 RNA 결합 단백질의 발현, 세포주기 상태 등 다양한 요인의 영향을 받는다. 본 탐구는 핵심 메커니즘을 선명하게 드러내기 위해 변수를 제한했지만, 그만큼 실제 현상의 다요인성을 축소하였다. 특히 특정 유전자 사례를 직접 설정하지 않았기 때문에 CD44, Bcl-x, c-src와 같은 실제 표적에서 관찰되는 정교한 조절 양상은 충분히 반영되지 못했다.

후속 탐구에서는 특정 유전자를 하나 선택하여 해당 엑손의 ESE와 ESS 분포, 알려진 SR 단백질 또는 hnRNP 결합 부위를 함께 분석하는 방식으로 더 구체화할 수 있다. 또한 공개된 RNA-seq 기반 대안적 스플라이싱 데이터나 실제 논문에서 제시된 PSI 값을 비교한다면, 본 보고서의 가상 모델을 실제 데이터와 대조하는 방향으로도 확장할 수 있다. 이런 과정은 스플라이싱 조절을 추상적 개념이 아니라 실제 데이터로 해석하는 데 도움이 될 것이다.

5.4 결론

본 탐구는 대안적 스플라이싱이 전사 후 일어나는 부가적 정리 과정이 아니라, 최종 단백질 산물의 정체성과 기능을 결정하는 핵심 조절 단계임을 보여준다. SR 단백질과 hnRNP는 각각 촉진과 억제의 대표 인자로 작용하지만, 그 본질은 단순한 찬반 구도가 아니라 pre-mRNA 위에서 이루어지는 경쟁적 해석 과정에 있다. 엑손이 포함될지 배제될지는 DNA 서열 하나만으로 결정되지 않으며, splice site의 강도, enhancer와 silencer의 배치, 조절 단백질의 상대적 활성, 세포 환경이 함께 작동한 결과로 정해진다.

가상 데이터 분석은 이러한 이론을 정량적으로 뒷받침하였다. SR/hnRNP 비율이 증가함에 따라 PSI가 23.6%에서 79.3%까지 상승하였고, 특히 중간 구간에서 큰 변화 폭이 나타난 점은 스플라이싱 결정이 임계 구간을 가지는 비선형 조절이라는 해석을 가능하게 했다. 이는 스플라이싱을 단순한 절단 반응으로 이해하는 수준을 넘어, 유전자 발현 정보를 RNA 단계에서 재구성하는 의사결정 시스템으로 인식하게 만든다.

결국 SR 단백질과 hnRNP의 길항 작용은 exon inclusion과 exon exclusion을 결정하는 핵심 분기점이며, 이는 진핵세포가 제한된 유전자 수로도 높은 단백질 다양성을 확보하는 중요한 전략이다. 동시에 이러한 균형의 붕괴는 비정상적 isoform 생성과 질병으로 이어질 수 있으므로, 스플라이싱 조절 메커니즘의 이해는 분자생물학의 기초 개념을 넘어 의생명과학 전반의 중요한 문제와 연결된다. 본 탐구는 교과에서 배운 RNA 가공 개념을 출발점으로 하여, 분자 수준의 조절 논리와 정량적 사고를 결합해 유전자 발현의 복잡성을 해석하고자 했다는 점에서 의미가 있다.

참고문헌

[1] Fu, X.-D., & Ares, M., Jr. (2014). Context-dependent control of alternative splicing by RNA-binding proteins. Nature Reviews Genetics, 15(10), 689–701. https://pmc.ncbi.nlm.nih.gov/articles/PMC4440546/

[2] Han, S. P., Tang, Y. H., & Smith, R. (2010). Functional diversity of the hnRNPs: Past, present and perspectives. Biochemical Journal, 430(3), 379–392.

[3] Long, J. C., & Caceres, J. F. (2009). The SR protein family of splicing factors: Master regulators of gene expression. Biochemical Journal, 417(1), 15–27.

[4] Shepard, P. J., & Hertel, K. J. (2009). The SR protein family. Genome Biology, 10(10), 242.

[5] Howard, J. M., & Sanford, J. R. (2015). The RNAissance family: SR proteins as multifaceted regulators of gene expression. Wiley Interdisciplinary Reviews: RNA, 6(1), 93–110. https://pmc.ncbi.nlm.nih.gov/articles/PMC5303883/

[6] Änkö, M.-L. (2014). Regulation of gene expression programmes by serine-arginine rich splicing factors. Seminars in Cell & Developmental Biology, 32, 11–21.

[7] Hertel, K. J. (2008). Combinatorial control of exon recognition. Journal of Biological Chemistry, 283(3), 1211–1215.

[8] Chen, M., & Manley, J. L. (2009). Mechanisms of alternative splicing regulation: Insights from molecular and genomics approaches. Nature Reviews Molecular Cell Biology, 10(11), 741–754.

[9] Wang, Y., Ma, M., Xiao, X., & Wang, Z. (2024). Interdependent regulation of alternative splicing by SR and hnRNP splicing factors across cassette exons. Nucleic Acids Research, 52(13), 7595–7615. https://pubmed.ncbi.nlm.nih.gov/39229091/

[10] Nilsen, T. W., & Graveley, B. R. (2010). Expansion of the eukaryotic proteome by alternative splicing. Nature, 463(7280), 457–463.

[11] Black, D. L. (2003). Mechanisms of alternative pre-messenger RNA splicing. Annual Review of Biochemistry, 72, 291–336.

활동 요약 (자기평가서)

진핵세포의 대안적 스플라이싱에서 SR 단백질과 hnRNP의 길항 작용이 약한 splice site를 가진 cassette exon의 포함과 배제를 어떻게 결정하는지 탐구함. 생명과학Ⅱ에서 배운 전사와 RNA 가공 개념을 바탕으로 스플라이소솜 조립 과정, exon definition, ESE와 ESS의 기능, 강한 splice site와 약한 splice site의 차이를 문헌으로 정리함. 이후 SR 단백질은 spliceosome recruitment를 촉진하고 hnRNP는 이를 방해하는 일반 경향을 비교하고, 두 인자의 상대적 활성 비율을 독립변수로 두어 PSI 변화를 가상 데이터로 모델링함. 비율 증가에 따라 엑손 포함률이 상승하고 특히 중간 구간에서 변화 폭이 커지는 양상을 해석하며 스플라이싱이 자동적 절단이 아니라 RNA 서열과 조절 단백질의 경쟁에 의해 결정되는 전사 후 조절 단계임을 구체적으로 확인함. 분자생물학 수업에서 접한 개념을 실제 연구 질문 형태로 좁혀 정리하고 정량적 해석까지 연결한 탐구임.