손박사의 도킹과 분자동역학 첫걸음

45,000개에서 하나를 골랐다.
그 하나를 얼마나 믿어야 하나

가상탐색은 깔때기입니다. 싼 방법으로 많이 버리고, 비싼 방법은 남은 소수에만 씁니다. 그런데 앞 단계에서 잘못 버린 것은 절대 돌아오지 않고, 앞 단계가 가장 부정확합니다.
이 20강은 ZINC → LigPrep → HTVS → ADMET → Vina → SP → XP → GROMACS → MM-PBSA 를 한 논문의 실제 값으로 따라가며, 도킹 점수 −20.66 kcal/mol 이 왜 자유에너지가 아닌지, 왜 탄닌이 1등을 했는지, 484 kcal/mol 이 왜 있을 수 없는 값인지를 직접 계산합니다.

🐍 Python🔗 도킹🧊 GROMACS ⚖️ MM-PBSA💊 ADMET🚩 위양성🇰🇷 한국어

1강 시작하기 →

이렇게 배웁니다

🔗깔때기를 한 단계씩

45,000 → 7,779 → 300 → 50 → 10. 각 단계에서 무엇을 버렸고 왜 버렸는지를 하나씩 따집니다.

📄진짜 논문의 진짜 숫자

Table 2 의 gscore −20.66, ADMET 의 logP 3.43, MD 의 200 ns / 20 ns — 논문에 실린 값을 그대로 씁니다.

🚩속는 법을 먼저

점수를 크기로 나눠 보면 1등이 꼴찌가 됩니다. 탄닌·PAINS·부호 오류·범위 밖의 값 — 직접 계산해 봐야 안 속습니다.

🛠️만지며 익히기

강마다 실험실 위젯이 하나씩 — 컷오프를 밀고, 유연성을 올리고, 궤적을 늘이고, MM-PBSA 항을 켜고 꺼 봅니다.

전체 커리큘럼

45,000개에서 하나로 — 깔때기의 모양

Nath 2017 「Ligand Preparation」~「Protein-Ligand Docking」

ZINC — 살 수 있는 분자들의 목록

「the ligands were downloaded from the ZINC database」

리간드 준비 — 분자 하나가 여럿이 되는 이유

「the ligands were prepared and optimized using LigPrep」

도킹이 푸는 두 문제 — 자세와 점수

어디에 어떻게 놓을 것인가, 그리고 그게 얼마나 좋은가

kcal/mol 인데 자유에너지가 아니다

Table 2 「Glide gscore」 −20.66 kcal/mol 을 어떻게 읽나

무엇을 움직이게 둘 것인가 — 강체·유연 리간드·유연 단백질

「In XP docking, flexibility of the interacting regions were considered」

HTVS → SP → XP — 깔때기를 세 번 좁힌다

「a cut off docking score value of −7 kcal/mol」

Vina 의 10.9 와 Glide 의 −20.66 — 견줄 수 없는 두 눈금

「Only ligands with an affinity value of 10.9 and above were considered」

Fig. 7 읽는 법 — 2D 상호작용 그림에 없는 것

Fig. 7 「Hydrogen bond interactions within the Protein and ligand 15 complex」

약물다움 — 리핀스키 다섯 규칙과 그 예외

Fig. 6 「Parameters used to screen ligands based on drug likeliness」

ADMET — 붙는 것과 약이 되는 것은 다르다

FAF-Drugs3 (Mobyle server) · 7,779개가 남은 자리

탄닌이 1등을 했다 — 위양성의 교과서적 사례

「Beta-1,2,3,4,6-Penta-O-Galloyl-D-Glucopyranose … a naturally occurring tannin」

분자동역학이 도킹에 없는 것을 준다 — 시간

GROMACS · 단백질 200 ns · 복합체 20 ns

상자 짓기 — 물·이온·힘장

「SPC/E water model and AMBER99SB-ILDN protein force-field … 3 Cl atoms were added」

에너지 최소화와 평형화 — 왜 먼저 눌러 주나

「energetically minimized for 0.1 ns and for a maximum force ≥ 10.0 kJ/mol」

RMSD 그림 읽는 법 — 이 논문은 평평해지지 않았다

Fig. 8a · 「the deviation … continues to increase for the full simulation length」

RMSF 와 경첩 잔기 — 어디가 흔들리나

Fig. 8b·c · 최대 요동 209·210·212–215·256·344·345 · 경첩 195·221·274·310·322·332·361

수소결합 수 · SASA · 2차 구조 — 나머지 그림들

「higher numbers of hydrogen bonds between the protein residues than with non-protein residues」

【실습】 MM-PBSA — 200~484 kcal/mol 이 왜 이상한가

「the top hundred ligands based on binding energy with energy in the range of 200–484 kcal/mol」

정리 — 인실리코 결과를 실험으로 넘기는 문턱

무엇이 있으면 이 논문이 강해졌을까