伞形采样是一种强大的增强采样技术非常适合研究蛋白质-蛋白质相互作用。下面我将为您提供一个详细的操作指南说明如何使用GROMACS进行伞形采样来研究Her2蛋白的聚集过程第一阶段准备系统确定反应坐标对于Her2蛋白聚集最简单的反应坐标是两个蛋白质质心之间的距离也可以选择关键结构域或残基之间的距离确定初始和最终状态初始状态两个Her2蛋白分离最终状态两个Her2蛋白形成稳定复合物第二阶段生成初始构象拉伸模拟创建拉伸模拟的MDP文件pull.mdp; 基本MD参数 integrator md dt 0.002 ; 时间步长(ps) nsteps 500000 ; 总步数(1 ns) nstxout 5000 nstvout 5000 nstfout 5000 nstxtcout 1000 ; 温度和压力耦合 tcoupl v-rescale tc-grps Protein Non-Protein tau-t 0.1 0.1 ref-t 300 300 pcoupl Parrinello-Rahman pcoupltype isotropic tau-p 2.0 ref-p 1.0 compressibility 4.5e-5 ; 非键相互作用 cutoff-scheme Verlet coulombtype PME rcoulomb 1.0 vdwtype Cut-off rvdw 1.0 ; 约束 constraints h-bonds constraint-algorithm LINCS ; 拉伸设置 pull yes pull-ngroups 2 pull-group1-name ProteinA ; 第一个Her2蛋白(根据您的索引文件命名) pull-group2-name ProteinB ; 第二个Her2蛋白 pull-ncoords 1 pull-coord1-type umbrella pull-coord1-geometry distance pull-coord1-groups 1 2 pull-coord1-dim Y Y Y ; 所有方向都拉伸 pull-coord1-start yes ; 从当前距离开始 pull-coord1-rate 0.01 ; nm/ps缓慢拉伸 pull-coord1-k 1000 ; kJ/mol/nm^2, 力常数创建索引文件定义拉伸组# 创建索引文件 gmx make_ndx -f equilibrated_system.gro -o pull.ndx # 在交互提示符下创建两个Her2蛋白的组 # 假设第一个Her2蛋白是残基1-500第二个是残基501-1000 r 1-500 t Protein name 8 ProteinA r 501-1000 t Protein name 9 ProteinB q执行拉伸模拟# 准备TPR文件 gmx grompp -f pull.mdp -c equilibrated_system.gro -p topol.top -n pull.ndx -o pull.tpr -maxwarn 1 # 运行拉伸模拟 gmx mdrun -deffnm pull -pf pullf.xvg -px pullx.xvg从轨迹中提取构象# 分析轨迹提取多个窗口的初始构象 gmx trjconv -s pull.tpr -f pull.xtc -o conf.gro -sep这会生成一系列名为conf0.gro、conf1.gro等的文件代表拉伸过程中的不同构象。第三阶段设置伞形采样窗口准备伞形采样MDP文件umbrella.mdp; 基本MD参数类似于pull.mdp但有一些修改 integrator md dt 0.002 nsteps 250000 ; 每个窗口0.5 ns nstxout 5000 nstvout 5000 nstfout 5000 nstxtcout 1000 ; 温度和压力耦合 (与pull.mdp相同) tcoupl v-rescale tc-grps Protein Non-Protein tau-t 0.1 0.1 ref-t 300 300 pcoupl Parrinello-Rahman pcoupltype isotropic tau-p 2.0 ref-p 1.0 compressibility 4.5e-5 ; 非键相互作用 (与pull.mdp相同) cutoff-scheme Verlet coulombtype PME rcoulomb 1.0 vdwtype Cut-off rvdw 1.0 ; 约束 (与pull.mdp相同) constraints h-bonds constraint-algorithm LINCS ; 伞形采样设置 pull yes pull-ngroups 2 pull-group1-name ProteinA pull-group2-name ProteinB pull-ncoords 1 pull-coord1-type umbrella ; 伞形约束 pull-coord1-geometry distance pull-coord1-groups 1 2 pull-coord1-dim Y Y Y pull-coord1-start no ; 使用参考值 pull-coord1-init WINDOW_DIST ; 将在脚本中替换 pull-coord1-rate 0 ; 不进行拉伸 pull-coord1-k 1000 ; kJ/mol/nm^2创建伞形采样窗口编写一个简单的脚本来设置和运行所有窗口#!/bin/bash # 定义窗口数量和初始距离 WINDOW_START0.5 # 起始距离(nm) WINDOW_END4.0 # 结束距离(nm) WINDOW_STEP0.2 # 窗口间距(nm) # 创建目录存放所有窗口 mkdir -p umbrella_windows # 为每个窗口创建MDP文件和运行模拟 WINDOW_DIST$WINDOW_START WINDOW_NUM0 while (( $(echo $WINDOW_DIST $WINDOW_END | bc -l) )); do # 创建窗口目录 WINDOW_DIRumbrella_windows/window_${WINDOW_NUM} mkdir -p $WINDOW_DIR # 创建MDP文件替换窗口距离 sed s/WINDOW_DIST/$WINDOW_DIST/g umbrella.mdp $WINDOW_DIR/umbrella_${WINDOW_NUM}.mdp # 使用最接近的构象作为起点 # 寻找最接近当前窗口距离的构象 CLOSEST_CONF$(grep WINDOW_DIST pullx.xvg | awk -v target$WINDOW_DIST \ {if(sqrt(($2-target)^2)min || min) {minsqrt(($2-target)^2); file$1}} END {print file}) # 准备TPR文件 gmx grompp -f $WINDOW_DIR/umbrella_${WINDOW_NUM}.mdp -c conf${CLOSEST_CONF}.gro \ -p topol.top -n pull.ndx -o $WINDOW_DIR/umbrella_${WINDOW_NUM}.tpr -maxwarn 1 # 运行这个窗口的模拟 gmx mdrun -deffnm $WINDOW_DIR/umbrella_${WINDOW_NUM} -pf $WINDOW_DIR/pullf_${WINDOW_NUM}.xvg \ -px $WINDOW_DIR/pullx_${WINDOW_NUM}.xvg # 更新窗口距离和计数器 WINDOW_DIST$(echo $WINDOW_DIST $WINDOW_STEP | bc -l) WINDOW_NUM$((WINDOW_NUM 1)) done注意如果您的系统很大可能需要在集群上并行运行这些窗口。第四阶段PMF计算收集所有窗口的数据# 创建tpr文件列表 ls umbrella_windows/window_*/umbrella_*.tpr tpr-files.dat # 创建pullx文件列表 ls umbrella_windows/window_*/pullx_*.xvg pullx-files.dat使用WHAM方法分析数据# 运行WHAM分析 gmx wham -it tpr-files.dat -if pullx-files.dat -o pmf.xvg -hist histo.xvg计算误差估计# 使用bootstrap方法估计误差 gmx wham -it tpr-files.dat -if pullx-files.dat -o pmf_errors.xvg -bs 50第五阶段分析结果PMF曲线分析PMF最小值位置代表最稳定的蛋白质-蛋白质距离从分离到结合的能垒高度反映结合动力学结合能可以从PMF深度估计结构分析分析不同窗口的构象特别是PMF最小值处的构象识别关键的相互作用残基# 提取PMF最小值处的构象 gmx trjconv -s umbrella_windows/window_X/umbrella_X.tpr -f umbrella_windows/window_X/umbrella_X.xtc -o bound_structure.pdb -b 400检验聚集界面分析氢键和疏水接触gmx hbond -f umbrella_windows/window_X/umbrella_X.xtc -s umbrella_windows/window_X/umbrella_X.tpr -num hbnum.xvgHer2蛋白特异性考虑对于Her2蛋白您应特别注意以下方面选择合适的反应坐标Her2蛋白具有多个结构域可能需要考虑特定结构域之间的距离例如可以选择胞外域(ECD)或酪氨酸激酶域(TKD)之间的距离设置合理的窗口范围Her2是大型蛋白窗口距离范围应该足够大例如0.5-8.0 nm窗口间距可以根据距离调整近距离时更密集如0.1 nm远距离时更稀疏如0.3 nm延长采样时间对于大型蛋白复合物每个窗口可能需要1-5 ns的采样时间检查每个窗口的构象是否充分平衡和采样实用提示系统准备使用质心距离作为初始反应坐标简单且有效确保系统在拉伸前已充分平衡技术细节力常数选择约1000-3000 kJ/mol/nm²根据需要调整窗口数量通常15-30个窗口足以覆盖感兴趣的反应坐标范围避免窗口之间有间隙确保良好的直方图重叠计算资源伞形采样非常适合并行计算每个窗口可以独立运行对于Her2蛋白这样的大系统考虑使用GPU加速通过这个详细的操作指南您应该能够成功设置和执行Her2蛋白聚集的伞形采样模拟。如果在实施过程中遇到具体问题请随时提问以获取更针对性的帮助。[ref:1,8,2,4]