在生命科學的廣闊領域中，蛋白質與小分子配體之間的相互作用扮演著至關重要的角色。這些相互作用不僅影響著生物體內的各種生命活動，如信號傳導、代謝調控和藥物作用等，同時也是藥物設計和開發的核心內容。因此，深入理解并模擬這些相互作用過程，對于推動生命科學研究和藥物研發具有重要意義。

本教程旨在為讀者提供一套完整的蛋白質與小分子配體相互作用模擬的流程和方法。通過本教程的學習，您將能夠掌握從蛋白質與小分子配體的結構準備、相互作用模擬到結果分析的全流程，從而能夠自主進行相關的模擬研究。

在本教程中，我們將首先介紹蛋白質與小分子配體相互作用的基本原理和模擬的基本概念，為讀者奠定理論基礎。隨后，我們將詳細闡述模擬的具體步驟，包括結構準備（如蛋白質結構預測、小分子結構優化等）、相互作用模擬（如分子對接、分子動力學模擬等）以及結果分析（如相互作用能計算、軌跡分析等）。在每個步驟中，我們都會結合具體的案例和實例，詳細解釋操作步驟和注意事項，幫助讀者更好地理解和掌握。

具體流程：

一、預處理復合物

1. 蛋白質及配體結構獲取

在本教程中，我們將使用T4溶菌酶L99A/M102Q(PDB ID：3HTB)為例，從PDB蛋白數據庫 (RCSB PDB)下載其晶體結構，去掉晶體水，PO4和 BME。

1. 蛋白及配體力場獲取

只有在力場的.rtp文件中存在構建塊的條目時，拓撲才能自動組裝。而JZ4配體在 GROMACS 提供的任何力場中都不是一個可識別的實體，因此我們將分兩步準備系統拓撲：1)用pdb2gmx準備蛋白質拓撲；2)使用外部工具準備配體拓撲。

2.1 使用pdb2gmx準備蛋白質拓撲

本教程使用的力場為amber14sb.ff，選擇默認的水模型TIP3P，然后為封端選擇“NH3+”和“COO-”，獲得力場文件及完整坐標文件。

2.2 使用外部工具獲得配體拓撲

本教程使用的力場為amber14sb.ff，因此使用GAFF工具生成JZ4配體的top文件。1)使用Avogadro軟件為配體添加氫原子，同時輸出JZ4.com文件，修改Gaussian設置，獲得Gaussian輸入文件JZ4.gjf，進行 Gaussian 優化。2)利用 AmberTools 計算電荷。3)利用 parmchk 檢查成鍵相缺失。4)利用LEaP生成 Amber 格式力場，文件內容見下圖。5)利用acpype將Amber格式轉換為Gormacs格式的GAFF力場文件及坐標文件。

1. 組合蛋白質和配體，生成蛋白質-配體復合物

二、定義盒子，添加溶劑及離子

三、能量最小化

四、限制復合物及體系平衡

1. 限制復合物：通過genrestr創建位置限制文件，定義位置限制。

2. NVT平衡

3. NPT平衡

五、成品模擬

六、分析

1. 執行energy模塊計算蛋白質-配體相互作用

2. 執行rms模塊，計算RMSD

最后，需要相關培訓或者項目合作歡迎通過公眾號“320科技工作室”聯系我們。