In Silico Analysis of Pathogenic Missense Single Nucleotide Polymorphisms in the Bovine IGF-1 Gene

소의 게놈을 살아있는 생명체를 만들고 운영하기 위한 지침서가 담긴 거대하고 오래된 도서관이라고 상상해 보십시오. 이 도서관 내부에서 A, C, G, T는 알파벳이며, 이 글자들이 적절하게 배열될 때 동물을 살리고, 성장시키고, 우유를 생산하게 하는 작은 분자 기계인 단백질의 레시피를 써 내려갑니다. 하지만 때때로 요리책의 오타처럼, 레시피 속의 글자 하나가 다른 글자로 바뀌기도 합니다. 유전학의 세계에서는 이를 단일 염기 다형성(Single Nucleotide Polymorphism), 즉 SNP라고 부릅니다. 이러한 오타 대부분은 국 요리법에서 "소금(salt)"을 "설트(selt)"로 바꾸는 것처럼 해롭지 않습니다. 국 맛은 여전히 괜찮을 테니까요. 그러나 어떤 오타는 "설탕(sugar)"을 "소금(salt)"으로 바꾸어 요리를 망쳐버리는 것처럼 치명적입니다.

소의 도서관에서 가장 중요한 레시피 중 하나는 IGF-1이라는 호르몬에 관한 것입니다. IGF-1을 건설 현장의 현장 소장이라고 생각해보십시오. 현장 소장은 근육을 만들고, 뼈를 키우며, 우유 생산에 필요한 유선(젖샘)을 발달시키라는 명령을 노동자들에게 외칩니다. 만약 레시피의 오타 때문에 현장 소장의 지시 사항이 엉망이 된다면, 건설 현장은 멈춰 서거나 건물이 삐뚤게 지어질 수도 있습니다. 낙농업자들에게 고장 난 현장 소장은 적은 양의 우유와 적은 수의 송아지를 의미하며, 이는 경제에 나쁜 소식입니다. 과학자들은 이러한 유전적 오타가 존재한다는 사실을 오래전부터 알고 있었지만, 수백만 개의 해로운 글자들 사이에서 실제로 기계를 고장 내는 특정 오타를 찾아내는 것은 거대한 시계 속에서 단 하나의 고장 난 톱니바퀴를 찾는 것과 같습니다. 여기서 "인 실리코(in silico)" 분석의 힘이 발휘됩니다. 연구자들은 축사에서 직접 소를 테스트하는 대신, 슈퍼컴퓨터를 사용하여 이러한 오타의 영향을 시뮬레이션합니다. 이는 마치 디지털 탐정이 되어 실제 현실에서 문제가 발생하기 전에 어떤 유전적 오류가 가장 큰 문제를 일으킬지 예측하는 것과 같습니다.

이 디지털 탐정 이야기에서, 이란 로레스탄 대학교의 연구팀은 소의 IGF-1 유전자에서 가장 심각한 오타를 찾아내기 위해 나섰습니다. 그들은 먼저 66가지의 서로 다른 미스센스 SNP(missense SNP)—이는 실제로 하나의 아미노산(단백질의 구성 단위)을 다른 것으로 바꾸어 IGF-1 현장 소장의 형태를 변화시킬 수 있는 특정 오타를 말합니다—의 목록을 수집했습니다. 연구진은 정교한 컴퓨터 도구 세트를 사용하여 소음 속에서 위험 신호를 걸러내는 필터 역할을 수행했습니다.

먼저, 그들은 이 목록을 전문가 심사위원단 역할을 하는 PredictSNP라는 도구에 통과시켰습니다. 66명의 용의자 중 18명이 "해로운(deleterious)" 것으로 분류되었는데, 이는 해당 변이가 해를 끼칠 가능성이 높음을 의미합니다. 연구진은 이 18개를 대상으로 이들이 단백질의 진화 역사에 얼마나 잘 부합하는지 자세히 조사했습니다. ConSurf라는 도구를 사용한 결과, 이 나쁜 지점들이 진화적으로 고도로 보존된 영역, 즉 생존에 필수적이기 때문에 수백만 년 동안 똑같이 유지되어 온 단백질 부위에 위치하고 있음을 발견했습니다. 이는 자동차의 스티어링 휠(핸들) 부분에서 오타가 발생한 것을 발견한 것과 같습니다. 그곳을 건드리면 안 된다는 것을 이미 알고 있는 것이죠.

다음으로, 연구팀은 이러한 돌연변이가 단백질의 물리적 안정성에 어떤 영향을 미칠지 시뮬레이션했습니다. 그들은 디지털 저울(I-Mutant2.0 및 MUpro)을 사용하여 변이가 단백질을 흔들리게 하거나 분해되게 만드는지 확인했습니다. 결과는 거의 모든 18가지 변이가 튼튼한 강철 빔을 약한 플라스틱 빔으로 교체한 것처럼 단백질을 덜 안정적으로 만든다는 것을 보여주었습니다. 또한 MEDUSA라는 도구를 사용하여 단백질의 유연성을 점검했는데, 이 변이들이 구조의 딱딱한 부분에 박혀 있어 단백질이 제 기능을 수행하는 것을 더욱 어렵게 만든다는 것을 발견했습니다.

조사는 MutPred2라는 도구를 통해 더욱 좁혀졌습니다. 이 도구는 돌연변이가 단백질의 분자적 메커니즘을 어떻게 방해할지 예측합니다. 이 단계에서 특히 위험한 7개의 "슈퍼 빌런(super-villains)"이 식별되었습니다. 마지막으로, 연구진은 Project HOPE를 사용하여 단백질의 3D 구조를 살펴보았습니다. 이 도구는 정상적인 단백질과 변이된 단백질의 크기, 전하, 그리고 친수성(분자가 물을 얼마나 좋아하거나 싫어하는지)을 비교했습니다. 그들은 L105R 변이가 가장 파괴적이라는 것을 발견했습니다. 이 변이에서는 중성이고 중간 크기인 아미노산(류신, Leucine)이 크고 양전하를 띤 아미노산(아르기닌, Arginine)으로 바뀌었습니다. 연구진은 이 변화가 단백질 구조에 "혹(bumps)"을 만들 것이라고 시뮬레이션했습니다. 왜냐하면 새로운 조각이 너무 크고, 양전하가 분자의 다른 부분들을 밀어내어 단백을 결합시키는 섬세한 소수성 상호작용을 방해하기 때문입니다.

이 L105R 변이의 실제 영향을 확인하기 위해, 팀은 단백질-단백질 도킹 시뮬레이션을 수행했습니다. 그들은 정상적인 IGF-1 단백질과 변이된 버전이 그들의 파트너인 IGF-1 수용체와 어떻게 악수를 나누는지 모델링했습니다. 시뮬레이션에서 정상 단백질은 25개의 수소 결합을 통해 강력한 악수를 형성했습니다. 그러나 변이된 L105R 단백질은 겨우 21개의 수소 결합만을 형성할 수 있었습니다. 이러한 연결의 감소는 변이된 단백질이 수용체와 효과적으로 결합하는 데 어려움을 겪을 것임을 시사하며, 이는 성장과 우유 생산에 필요한 통신 라인을 끊어버릴 가능성이 있음을 의미합니다.

본 논문은 이러한 발견이 아직 살아있는 소를 대상으로 테스트된 것이 아니라 컴퓨터 시뮬레이션에 기반한 것이지만, L105R 변이가 소의 IGF-1 단백질 기능을 방해할 가능성이 가장 높은 후보라는 점을 결론짓습니다. 저자들은 이 특정 변이가 육종 프로그램에서 분자 표지(molecular marker)로 활용될 수 있으며, 이를 통해 농부들이 이 유전적 "오타"를 가질 가능성이 낮고 생산성이 높은 소를 선택하는 데 도움을 줄 수 있다고 제안합니다. 이 연구는 유제품 생산 문제를 해결했다고 주장하는 것이 아니라, 유전적 기능 장애의 매우 개연성 높은 용의자를 제시함으로써 향 향후 실제 테스트를 위한 길을 열어주고 있습니다.