手機版
化工儀器網(wǎng)手機版
化工儀器網(wǎng)小程序
官方微信
公眾號：chem17
掃碼關(guān)注視頻號
科研助攻 | 分子動力學模擬：藥物設(shè)計的超能力！ | MedChemExpress (MCE)
Hey, 準備好解鎖藥物設(shè)計的最·前·沿科技--分子動力學模擬！其能夠?qū)崟r觀察分子如何與疾病靶點“對話”，預(yù)測它們的每一次親密接觸！這究竟是如何實現(xiàn)的？一文解析~
01
分子模擬的簡要介紹
在介紹分子動力學模擬方法之前，我們首先來認識一下分子模擬的概念，分子模擬是一種計算化學方法，它通過計算機程序模擬分子或分子體系的行為，以研究其結(jié)構(gòu)、性質(zhì)和相互作用。這種方法能夠在原子或分子水平上提供對化學過程和物理現(xiàn)象的深入理解，是實驗和理論化學的重要補充。
分子模擬主要包括分子力學 (Molecular Mechanics, MM)、蒙特卡洛 (Monte Carlo, MC) 模擬、分子動力學 (Molecular Dynamics, MD) 模擬、量子化學計算、第·一性原理分子動力學 (Ab Initio Molecular Dynamics, AIMD)、粗粒化模型 (Coarse-Grained Models) 以及自由能計算等方法 (圖 1)，其中應(yīng)用較為廣泛的是分子動力學方法和蒙特卡洛方法。與蒙特卡洛模擬相比，分子動力學模擬在多個方面展現(xiàn)出其獨·特的優(yōu)勢，尤其在處理分子體系的動態(tài)行為和時間分辨問題方面。
圖 1. 結(jié)合分子動力學模擬的虛擬篩選工作流程[1]。
02
分子動力學模擬的作用原理
分子動力學模擬的作用原理是利用經(jīng)典力學的牛頓運動定律，結(jié)合統(tǒng)計力學的概念，通過數(shù)值方法模擬原子或分子在給定力場作用下的運動和相互作用。
在模擬開始時，設(shè)定分子體系的初始構(gòu)型和速度，這些初始條件通?；趯嶒炗^測或理論預(yù)測。隨后，模擬中的每個原子或分子根據(jù)其所受的力場作用，通過數(shù)值積分方法 (如 Verlet 或 Beeman 算法，圖 2) 來更新其位置和速度，這些力場包括成鍵相互作用和非鍵相互作用，可以通過參數(shù)化的勢函數(shù)來描述。在每個時間步長中，系統(tǒng)的狀態(tài) (包括原子的位置和速度) 都會被更新，從而模擬出原子在時間上的動態(tài)行為。時間步長的選擇對于保證模擬的準確性和穩(wěn)定性至關(guān)重要，它必須足夠小以捕捉到體系中的關(guān)鍵物理現(xiàn)象。在模擬過程中，可以對系統(tǒng)的溫度和壓強進行控制，以模擬不同的熱力學條件。通過收集和分析模擬過程中的數(shù)據(jù)，可以得到分子體系的結(jié)構(gòu)、動力學和熱力學性質(zhì)，如均方位移、相關(guān)函數(shù)和自由能等。
圖 2. 分子動力學基本算法[2]。
03
分子動力學在藥物研發(fā)中的應(yīng)用
分子動力學 (MD) 模擬技術(shù)不僅能夠揭示蛋白質(zhì)與配體之間復(fù)雜的能量相互作用，還為我們提供了深入理解生物大分子動力學結(jié)構(gòu)的豐富信息。MD 模擬已經(jīng)成為現(xiàn)代藥物發(fā)現(xiàn)過程中不·可或缺的一部分，其應(yīng)用廣泛且成效顯著。以下是 MD 模擬在藥物發(fā)現(xiàn)中的一些具體應(yīng)用實例。
分子動力學模擬揭示抗病毒化合物與 PrPC “熱點”區(qū)域的結(jié)合、穩(wěn)定機制
PrPC 的結(jié)構(gòu)轉(zhuǎn)變到 PrPSc (病理型朊蛋白) 是導致傳染性海綿狀腦病 (即朊疾病) 的主要原因，因此開發(fā)能抑制這一轉(zhuǎn)變的有效抗朊病毒·藥物是一個有希望的治療方向。
Shuangyan Zhou 等人使用分子動力學模擬方法，研究了幾種已知能特異性結(jié)合 PrPC“熱點”區(qū)域并有效穩(wěn)定其天然結(jié)構(gòu)的抗朊病毒化合物的穩(wěn)定化機制。結(jié)果發(fā)現(xiàn)抗朊病毒化合物通過三種機制穩(wěn)定 PrPC：(1) 穩(wěn)定 α2 的柔性 C 末端和疏水核心，如 BMD42-29 和 GN8；(2) 僅穩(wěn)定疏水核心，如 J1 和 GJP49；(3) 通過高結(jié)合親和力穩(wěn)定 PrPC 的整體結(jié)構(gòu)，如 NPR-056。這些相互作用關(guān)鍵殘基的確定對于未來 PrPSc 的藥物設(shè)計和發(fā)現(xiàn)非常重要。
圖 3. PrPC 與抗原化合物 BMD42-29、GN8、J1、NPR-056 和 GJP49 之間的詳細相互作用[3]。
分子動力學模擬助力 GPCR 藥物設(shè)計
設(shè)計 G 蛋白偶聯(lián)受體 (GPCR) 變構(gòu)調(diào)節(jié)劑是現(xiàn)代藥物研究中的一個活躍領(lǐng)域，但由于缺乏這些藥物的結(jié)合模式和分子機制的知識，這一直是一個挑戰(zhàn)。
Ron O. Dror等人利用分子動力學模擬揭示了 M2 毒蕈堿乙酰膽堿受體 (M2 受體) 的幾種變構(gòu)調(diào)節(jié)劑的結(jié)合位點、結(jié)合構(gòu)象和特異性藥物-受體相互作用。圖 4 還顯示了經(jīng)典配體結(jié)合的正、負變構(gòu)調(diào)節(jié)機制，包括兩個結(jié)合位點的耦合構(gòu)象變化和這些位點上配體之間的靜電相互作用。這項研究不僅增進了對 GPCR 變構(gòu)調(diào)節(jié)機制的理解，也為開發(fā)新型藥物提供了重要信息。
圖 4. 正、負變構(gòu)調(diào)節(jié)機制[4]。
分子動力學模擬闡明 INSTIs 交叉耐藥分子機制
HIV-1 整合酶是治療 HIV-1 感染的重要藥物靶點。然而，對 INSTIs (HIV 整合酶鏈轉(zhuǎn)移抑制劑) 的交叉耐藥性問題日益嚴重，特別是 E138K/Q148K 雙突變體對幾種重要的 INSTIs 表現(xiàn)出高水平的耐藥性。
Weiwei Xue 等人基于同源建模構(gòu)建的 HIV-1 整合酶四聚體結(jié)構(gòu)，使用分子動力學 (MD) 模擬和 RIN 分析來研究交叉耐藥突變對 INSTIs 的影響。MD 模擬顯示，INSTIs 與 HIV-1 整合酶活性位點的結(jié)合模式受到 E138K/Q148K 突變的影響，導致活性位點的結(jié)構(gòu)重排 (圖 5A、B)。通過結(jié)合自由能計算，發(fā)現(xiàn) Pro145 殘基在 INSTIs 與 HIV-1 整合酶結(jié)合中起到重要作用，特別是其與 INSTIs 的疏水相互作用 (圖 5C、D)。該研究不僅為理解 HIV-1 整合酶抑制劑的交叉耐藥機制提供了重要信息，而且為未來合理設(shè)計新的 INSTIs 以克服交叉耐藥性提供了理論基礎(chǔ)。
圖 5. 分子動力學模擬在闡明 HIV-1 鏈轉(zhuǎn)移抑制劑交叉耐藥分子機制中的應(yīng)用[5]。
04
小結(jié)
今天，小 M 帶大家深入了解了分子動力學模擬技術(shù)在藥物研發(fā)領(lǐng)域的應(yīng)用前景及核心原理。作為一種創(chuàng)新的技術(shù)，分子動力學模擬技術(shù)將在未來的藥物研發(fā)中發(fā)揮更加關(guān)鍵的作用~
| 產(chǎn)品推薦 |
| Virtual Screening 虛擬篩選 (Virtual Screening,VS) 是基于小分子數(shù)據(jù)庫開展的活性化合物篩選。利用小分子化合物與藥物靶標間的分子對接運算，虛擬篩選可快速從幾十至上百萬分子中，遴選出具有成藥性的活性化合物，大大降低實驗篩選化合物數(shù)量，縮短研究周期，降低藥物研發(fā)的成本。 |
| MCE 50K Diversity Library (HY-L901) MCE 50K Diversity Library 由 50,000 種類藥化合物組成。依據(jù)谷本相似性 (Tanimoto Coefficient) 及聚類算法 (Bemis-Murcko) 對上百萬化合物進行篩選以確保結(jié)構(gòu)多樣性。本多樣性庫具備新穎性、類藥性，化合物結(jié)構(gòu)類型多樣、化學空間豐富，庫中化合物可重復(fù)供應(yīng)，是新藥研發(fā)的有力工具，可以廣泛地應(yīng)用于高通量篩選 (HTS) 和高內(nèi)涵篩選 (HCS)。 |
| 分子動力學模擬 分子動力學 (Molecular Dynamics, MD) 模擬是一種基于牛頓力學，綜合了物理、數(shù)學和化學等多種學科的計算機模擬方法，用于研究分子體系的運動和相互作用、預(yù)測分子體系的行為和結(jié)構(gòu)性質(zhì)。 |
參考詳情：
[1] Liu X, Shi D, Zhou S, Liu H, Liu H, Yao X. Molecular dynamics simulations and novel drug discovery. Expert Opin Drug Discov. 2018 Jan;13(1):23-37.
[2] Go?i JR, Orozco M, Gelpí JL. Molecular dynamics simulations: advances and applications. Adv Appl Bioinform Chem. 2015 Nov 19;8:37-47.
[3] Zhou S, Liu X, An X, Yao X, Liu H. Molecular Dynamics Simulation Study on the Binding and Stabilization Mechanism of Antiprion Compounds to the "Hot Spot" Region of PrPC. ACS Chem Neurosci. 2017 Nov 15;8(11):2446-2456.
[4] Dror RO, Green HF, Valant C, Borhani DW, Valcourt JR, Pan AC, Arlow DH, Canals M, Lane JR, Rahmani R, Baell JB, Sexton PM, Christopoulos A, Shaw DE. Structural basis for modulation of a G-protein-coupled receptor by allosteric drugs. Nature. 2013 Nov 14;503(7475):295-9.
[5] Xue W, Jin X, Ning L, Wang M, Liu H, Yao X. Exploring the molecular mechanism of cross-resistance to HIV-1 integrase strand transfer inhibitors by molecular dynamics simulation and residue interaction network analysis. J Chem Inf Model. 2013 Jan 28;53(1):210-22.
相關(guān)產(chǎn)品
免責聲明
- 凡本網(wǎng)注明“來源：化工儀器網(wǎng)”的所有作品，均為浙江興旺寶明通網(wǎng)絡(luò)有限公司-化工儀器網(wǎng)合法擁有版權(quán)或有權(quán)使用的作品，未經(jīng)本網(wǎng)授權(quán)不得轉(zhuǎn)載、摘編或利用其它方式使用上述作品。已經(jīng)本網(wǎng)授權(quán)使用作品的，應(yīng)在授權(quán)范圍內(nèi)使用，并注明“來源：化工儀器網(wǎng)”。違反上述聲明者，本網(wǎng)將追究其相關(guān)法律責任。
- 本網(wǎng)轉(zhuǎn)載并注明自其他來源（非化工儀器網(wǎng)）的作品，目的在于傳遞更多信息，并不代表本網(wǎng)贊同其觀點和對其真實性負責，不承擔此類作品侵權(quán)行為的直接責任及連帶責任。其他媒體、網(wǎng)站或個人從本網(wǎng)轉(zhuǎn)載時，必須保留本網(wǎng)注明的作品第一來源，并自負版權(quán)等法律責任。
- 如涉及作品內(nèi)容、版權(quán)等問題，請在作品發(fā)表之日起一周內(nèi)與本網(wǎng)聯(lián)系，否則視為放棄相關(guān)權(quán)利。
采購中心