2025년 9월 26일
이 기사에서는 여러 리간드의 분자 특징을 통합하여 합의 약리단 모델을 구성하기 위한 프로토콜을 제시합니다. 이 방법은 알려진 리간드 결합 형태를 가진 모든 생물학적 표적을 표적으로 하는 약물 발견 노력에 적용할 수 있으며, 가상 스크리닝 및 합리적인 약물 설계를 위한 주요 상호 작용 특징을 식별할 수 있습니다.
우리의 연구는 대규모 리간드 라이브러리를 사용하여 약리단을 모델링하는 것을 탐구합니다. 이 연구에서 우리는 약물 발견을 안내하는 구조적 패턴을 식별하기 위해 SARS-CoV-2 Mpro에 중점을 둡니다. 약리단 모델링의 최근 발전은 더 크고 다양한 데이터 세트를 사용하여 신약 발견을 위한 핵심 분자 피더를 식별하는 정확도를 향상시킵니다.
분자 역학 시뮬레이션과 머신 러닝은 현재 약리단 모델링을 발전시키고 신약 발견 결과를 개선하는 데 사용됩니다. 주요 과제는 크고 다양한 리간드 데이터 세트를 관리하는 동시에 향후 추출 정확도를 개선하고 계산 비용을 줄이는 것입니다. 우리는 고가의 리간드 라이브러리가 약리단의 미래 정확도를 향상시켜 보다 신뢰할 수 있는 예측과 잠재적인 약물 후보의 더 나은 식별을 가능하게 한다는 것을 보여주었습니다.
시작하려면 PyMOL 소프트웨어를 실행하세요. 소프트웨어에서 모든 단백질 리간드 복합체를 정렬합니다. 정렬된 각 복합체에 대해 리간드 상합체를 추출하고 SDF 형식의 개별 파일로 저장합니다.
각 리간드 파일을 개별적으로 업로드하여 Load Features 옵션을 사용하여 파밍합니다. 로드 후 세션 저장 옵션을 사용하여 약리단 파일을 JSON 형식으로 다운로드합니다. 다운로드한 모든 약리단 JSON 파일을 나중에 Google Colab 환경에서 사용할 수 있도록 단일 폴더에 함께 저장합니다.
웹 브라우저에서 Google Colab을 실행하고 새 노트북을 만듭니다. Conda 및 PyMOL을 설치한 다음 재생 아이콘을 클릭하거나 Shift 및 Enter 키를 눌러 셀을 실행하여 성공적인 실행을 확인합니다. 올바르게 완료되면 셀 위에 녹색 막대가 나타납니다.
이제 CANFAR Python 패키지를 설치하고 필요한 모듈을 가져옵니다. 재생 아이콘을 클릭하여 셀을 실행하여 성공적인 설치를 확인합니다. 실행에 성공하면 확인 메시지가 나타납니다.
다음으로, 약리단 JSON 파일을 저장하기 위한 새 폴더를 만듭니다. 폴더를 열고 그 안을 마우스 오른쪽 버튼으로 클릭한 다음 업로드를 선택하여 JSON 파일을 업로드합니다. 약리단 특징을 구문 분석하고 통합하려면 업로드된 모든 파일에서 특징을 추출하십시오.
재생 아이콘을 클릭하여 셀을 실행하고 각 개별 리간드의 특징을 포함하는 통합 데이터 프레임을 생성합니다. 모든 약리단 설명자를 표시한 다음 셀을 실행하여 클러스터링 패턴을 표시합니다. 이제 약리단 모델을 PyMOL 형식으로 저장합니다.
셀을 실행하여 pse 파일을 생성한 다음, 해당 코드 셀을 실행하여 모델을 JSON 파일로 저장합니다. 제공된 코드를 사용하여 기능 클러스터형 출력 및 덴도그램을 생성합니다. 재생을 클릭하고 내보내기 명령을 실행하여 합의 약리단 데이터를 CSV 파일로 내보냅니다.
가상 스크리닝을 위해 생성된 합의 약리단 모델을 사용합니다. 서로 다른 비공유 억제제와 공결정화된 100개의 주요 프로테아제 또는 Mpro 복합체를 정렬하여 결합 포켓 내의 리간드 다양성을 시각화했습니다. 각 리간드는 정렬된 복합체에서 개별적으로 추출되어 가상 스크리닝을 위한 별도의 파일을 생성했습니다.
추출된 리간드는 pharmit 서버에 업로드되어 약리단 모델링을 위한 구조화된 JSON 파일을 생성하는 데 사용되었습니다. JSON 파일은 1, 450개의 약리단 특징을 산출했으며, 23개의 방향족, 30개의 수소 결합 수용체, 16개의 수소 결합 공여체, 36개의 소수성 및 5개의 음이온 클러스터를 포함하여 110개의 클러스터로 그룹화되었습니다. 합의 모델을 구축하기 위해 가장 밀도가 높은 클러스터만 유지되고 구성원이 부족하여 음이온 클러스터는 제외되었습니다.
최종 합의 약리단 모델은 수소 결합 수용체에 대한 3개의 방향족, 2개의 수소 결합 공여체 및 2개의 소수성 특징으로 구성된 11개의 특징으로 구성되었습니다. 가장 대표적이지 않은 방향족 특징인 aro 1, 2개의 후보 화합물이 확인되었습니다. 화합물 101 267 741의 2차원 화학 구조를 시각화하여 약리단 모델과의 호환성을 확인하였다.
화합물 101 267 741은 Mpro 바인딩 포켓에 깊숙이 들어가 S1 포켓을 향해 퍼지는 참조 리간드 38a보다 S1 및 S2 서브 포켓을 더 효과적으로 차지합니다. 화합물 101 267 741은 11개의 분자간 상호 작용, 7개의 수소 결합 및 4개의 소수성 접촉을 형성했습니다.
본 논문은 여러 리간드로부터 분자적 특징을 통합하여 컨센서스 파마코포어 모델을 구축하는 프로토콜을 제시하며, 특히 SARS-CoV-2 Mpro에 초점을 맞춥니다. 이 방법은 가상 스크리닝 및 합리적 약물 설계를 위한 핵심 상호작용 특징을 식별함으로써 약물 발굴 노력을 강화합니다.
광범위한 리간드 라이브러리를 이용한 컨센서스 약리단 모델링은 초기 단계 신약 개발, 특히 SARS-CoV-2 Mpro와 같은 표적에 대해 더 높은 예측 신뢰도를 제공합니다. 다양한 리간드 결합 컨포메이션을 통합함으로써 모델의 편향을 줄이고 강력한 가상 스크리닝을 지원하여, 리드 화합물 발굴 및 포트폴리오 우선순위 결정에 직접적인 영향을 미칩니다. 이러한 접근 방식은 풍부한 구조적 데이터가 있는 표적에 대해 합리적인 후보 물질 선택 과정을 효율화하여 기업의 R&D 효율성을 높입니다.
본 합의 약리단 모델링 프로토콜은 초기 발견 단계부터 리드 화합물 도출 단계까지 적용 가능하며, 구조 생물학과 정보학을 활용하여 가상 스크리닝 및 후보 물질 우선순위 선정에 필요한 정보를 제공합니다.