三句代码绘制:分子动力学模拟 自由能景观形貌图
- 2026-02-06 19:31:36
点击上方蓝色，关注稳妥少年
author: "稳妥少年" date: "2026/01/21"
分子动力学模拟与自由能景观
蛋白质是生命活动的“机器”，这些纳米级的小机器并非静止不动，而是时刻在发生着复杂而精密的构象变化。要真正理解它们如何工作，单靠静态结构图片远远不够——我们需要“看到”它们在溶液中的动态过程。
分子动力学模拟通过计算机对蛋白质及其周围成千上万个水分子、离子进行物理建模，我们能够模拟它们在一定温度下随时间变化的运动轨迹。这一技术就像一台“计算显微镜”，让我们在原子尺度上观察蛋白质如何折叠、如何与药物分子结合、如何响应外界信号。模拟结果为我们理解疾病机理、设计新型药物提供了前所未有的动态视角。
然而，一个模拟往往包含数十万原子的运动数据，信息量巨大且高度复杂。我们该如何从中提取关键的运动模式，并理解蛋白质构象转变背后的能量规律呢？
从噪声到信号：主成分分析（PCA）降维提取主导运动
面对庞杂的原子运动轨迹，科学家常常借助主成分分析这一强大的统计工具。简单来说，PCA 能在一大堆看似混乱的运动中，找出几个最能代表整体构象变化的“主成分”，也就是最核心的运动模式。
具体分析流程如下：
计算协方差矩阵
gmx covar -s md_0_1.tpr -f fit.xtc -o eigenvalues.xvg -v eigenvectors.trr -xpma covapic.xpm这一步骤计算所有原子运动的相互关系，并提取特征向量（即主成分方向）与特征值（代表该方向运动的幅度）。
投影至主成分空间
gmx anaeig -s md_0_1.tpr -f fit.xtc -v eigenvectors.trr -first 1 -last 1 -proj pc1.xvg我们将蛋白质-小分子复合物的轨迹投影到最重要的两个主成分（PC1 和 PC2）上，将原本高维的运动简化到最能代表构象变化的二维平面上。
绘制能量地形图：自由能景观揭示稳定状态与转换路径
得到主成分投影后，我们就能进一步构建自由能景观。你可以将其想象成一个“能量地形图”：横纵坐标是 PC1 和 PC2，代表不同的构象状态；颜色高度则代表自由能的高低。
使用以下命令，我们能够将构象分布转化为自由能表面：
gmx sham -tsham 310 -nlevels 100 -f pc_output.xvg -ls pca12_gibbs.xpm -g pca12_gibbs.log -lsh pca12_enthalpy.xpm -lss pca12_entropy.xpm
通过可视化工具（如 xpm_show），我们可以选择不同的插值方式（如 spline16）来平滑并三维展示这张图谱：
dit xpm_show -f pca12_gibbs.xpm -m 3d
dit xpm_show -f pca12_gibbs.xpm -ip spline16
为什么自由能景观如此重要？
自由能景观不是一张普通的图表，它回答了蛋白质动态行为中几个最核心的问题：
稳定构象在哪里？
图中自由能的低谷（蓝色区域）对应蛋白质最稳定或亚稳定的构象状态。构象转变如何发生？
连接不同能量低谷的“鞍部”或路径，揭示了构象转换的可能途径与能量壁垒。小分子如何影响蛋白质？
比较结合小分子前后的自由能景观，可直观看出药物或配体是否稳定了特定构象、是否改变了转换屏障，从而从能量角度解释其作用机制。
通过将复杂的动态轨迹转化为直观的能量地形，自由能景观让我们得以“俯瞰”蛋白质在构象空间中的行走路径，理解其功能实现的能量基础。这不仅在基础研究中意义重大，更为理性药物设计、酶工程改造提供了关键的动态与热力学依据。
在计算模拟的助力下，我们正一步步揭开生命分子机器的动态秘密，让设计药物、改造蛋白质变得更加精准而高效。
相关链接：
直通GROMACS分子动力学模拟【蛋白质-小分子、蛋白质-蛋白质体系】
稳妥少年·药智未来 生信学习资料合订本（分子对接、动力学模拟、转录组分析、单细胞分析、分析流程），转发学习参与课程抽奖
以上就是本期全部内容，祝大家学习愉快，收获满满，请给稳妥少年点一个赞、并转发给新的小伙伴吧！！！
稳妥少年建立了第一个讨论群（稳妥少年_共勉学习交流群），特此欢迎新老朋友进群交流学习，让学习变简单！
交流已超200+，无法直接扫码加入，大家如果需要的可以添加微信，拉您进群的。添加微信时，一定要注明来意哈。
稳妥少年——单细胞分析实战
本课程系统介绍了单细胞RNA测序数据分析的全流程方法学体系。从原始数据读取开始，详细讲解了三种常见格式（标准10X输出、HDF5文件和稀疏矩阵）的处理策略，并统一转换为Seurat对象进行分析。质控阶段通过检测基因数、UMI总数和线粒体基因百分比等指标过滤低质量细胞，确保数据可靠性。标准化环节采用对数归一化和高变基因识别方法，通过主成分分析降维，使用JackStraw检验和肘部图确定最佳维度。
细胞聚类分析采用共享最近邻算法和Louvain聚类，在UMAP降维空间可视化细胞亚群。基于经典标记基因表达模式，对细胞类型进行系统注释。差异表达分析采用多层次策略，包括整体水平比较、细胞类型特异性分析和功能富集评估。细胞比例分析通过统计检验评估实验组间细胞组成的差异。
高级分析模块涵盖细胞间通讯网络构建（CellChat）、发育轨迹推断（Monocle）、代谢通路活性评估（scMetabolism）以及基因集富集分析（GSEA）。特别整合了GSVA通路活性评分、AddModuleScore模块化评分等创新方法，全面解析细胞功能状态。可视化方面采用多类型图表（火山图、热图、轨迹图、网络图等）呈现分析结果，所有方法均经过标准化流程验证，确保分析结果的可靠性和可重复性，为疾病机制研究和生物标志物发现提供系统解决方案。
【淘宝】限时满500减200 https://e.tb.cn/h.73SuB1NblkWjvxx?tk=QJ1vfvkGE3A HU293 「稳妥少年，单细胞分析实战课」
点击链接直接打开 或者 淘宝搜索直接打开
咸鱼
【闲鱼】https://m.tb.cn/h.73SGjJy?tk=tq8efvkucB3 HU071 「我在闲鱼发布了【【稳妥少年-单细胞分析实战全流程视频课】】」
点击链接直接打开
稳妥少年——生信技能课
本课程系统涵盖了生物信息学在基因组数据挖掘与分析中的核心流程与方法。课程从GEO、TCGA等公共数据库的数据检索与下载入门，进而学习差异表达分析（使用limma、DESeq2等工具）及结果可视化。在此基础上，深入进行功能富集分析（GO、KEGG）和蛋白质互作网络构建，并利用Cytoscape筛选关键基因模块。课程还整合了多种机器学习算法进行生物标志物筛选，并教授诊断列线图、神经网络模型等临床预测模型的构建与评估。此外，课程内容扩展至基因集富集分析（GSEA/GSVA）、免疫细胞浸润分析、基因的亚细胞与染色体定位、表型关联分析、表达相关性分析、单基因预后分析以及肿瘤免疫检查点分析等一系列高级专题，旨在提供从基础到前沿的完整生物信息学分析能力。
淘宝
【淘宝】限时满500减200 https://e.tb.cn/h.7X8RTTPSng5IcEN?tk=iJN5fGQL1mn MF278 「稳妥少年一生信技能课」
点击链接直接打开 或者 淘宝搜索直接打开
咸鱼
【闲鱼】https://m.tb.cn/h.72Tg8o6?tk=F34xfGQNQAt HU293 「我在闲鱼发布了【本课程《生物信息学全流程实战技能课》系统讲授疾病研究，尤其是】」
点击链接直接打开
TCGA肿瘤预后+单细胞+虚拟基因敲除
本研究系统整合了多组学数据与生物信息学方法，旨在探究多发性骨髓瘤的分子特征、预后标志及免疫微环境。首先，从TCGA与GEO数据库获取并标准化了批量转录组数据，通过差异表达分析鉴定了疾病相关的关键基因，并利用GO/KEGG富集分析阐释其功能。基于此，采用Cox回归与LASSO算法构建并验证了一个多基因预后风险模型，该模型在独立队列中展现出稳健的预测效能。进一步分析揭示了风险评分与临床病理特征的独立性，并通过列线图实现了个体化生存预测。
为深入解析肿瘤微环境，本研究应用ssGSEA评估了免疫细胞浸润水平，并分析了其与风险模型及免疫检查点表达的相关性。同时，利用单细胞RNA测序技术，在细胞分辨率下刻画了骨髓瘤微环境的细胞组成图谱，识别了主要免疫细胞亚群，并量化了其在不同临床状态下的比例变化。通过CellChat工具，系统推断了疾病组与对照组间细胞通讯网络的差异。此外，采用scTenifoldKnk算法对髓系细胞亚群进行了虚拟基因敲除分析，从网络层面探索了关键基因的调控功能。
综上，本研究通过从宏观群体水平到微观单细胞层面的多层次分析，构建了多发性骨髓瘤的预后预测体系，并深入揭示了其肿瘤免疫微环境的特征与潜在调控机制，为理解疾病生物学及开发治疗策略提供了新的见解和数据支持。
淘宝
【淘宝】https://e.tb.cn/h.74PYbf9cYTjXRw9?tk=o9KTfxHXS3R CZ005 「稳妥少年【肿瘤预后+单细胞+基因敲除】全流程视频课」
点击链接直接打开 或者 淘宝搜索直接打开
咸鱼
【闲鱼】https://m.tb.cn/h.7fqO1d0?tk=KEbifxHcGH5 CA381 「我在闲鱼发布了【肿瘤预后+单细胞+基因敲除全流程视频课】」
点击链接直接打开
稳妥少年【NHANES数据库使用教程】
这门课是面向临床研究的“流水线”实战培训，核心是教你把NHANES数据库的原始文件变成一张张能用于论文的图表和表格。它会手把手带你走通全流程：从用R代码自动下载、合并不同年份的数据模块，到严格清洗变量、筛选研究对象；从正确运用复杂的抽样权重进行基线统计与加权逻辑回归，到使用限制性立方样条（RCS）探索变量间的非线性关系；最后，指导你完成单因素/多因素分析，并绘制出可直接用于发表的森林图。整个过程紧扣实际研究需求，重点解决如何正确处理权重、合并数据、校正混杂因素等具体技术难题，目标是让你能独立完成一项从数据到结果的完整分析。
淘宝
【淘宝】https://e.tb.cn/h.7WtLo3unmwB0dUI?tk=7XlqfEbwdDS MF937 「稳妥少年—NHANES数据库使用 全流程视频课」
点击链接直接打开 或者 淘宝搜索直接打开
咸鱼
【闲鱼】https://m.tb.cn/h.73y4RJk?tk=yhJffEbvOgU MF287 「我在闲鱼发布了【稳妥少年—NHANES数据库使用】」
点击链接直接打开
稳妥少年网络毒理学课程
网络毒理学分析流程详细内容
（1）疾病转录组数据准备（2）毒物毒理学评估（3）毒物靶点基因的获取（4）疾病靶点基因的获取（5）差异表达基因的获取（6）候选基因的获取（7）GO与KEGG富集分析（8）毒物-基因-通路-疾病网络构建（9）蛋白质相互作用网络构建和核心基因筛选（10）表达水平验证筛选关键基因（11）相关性和功能相似性分析（12）转录因子-基因-miRNA调控网络构建（13）关键基因富集分析（14）GeneMANIA分析、亚细胞定位和染色体定位分析（15）分子对接（16）分子动力学模拟
课程购买链接
淘宝
【淘宝】https://e.tb.cn/h.SyH6B9i86ptdeNZ?tk=MiTYfLWdlC2 CZ005 「稳妥少年【网络毒理学全流程视频课】」
点击链接直接打开 或者 淘宝搜索直接打开
咸鱼
【闲鱼】https://m.tb.cn/h.SyhI8J3?tk=1BqmfLWUyIH CZ028 「我在闲鱼发布了【稳妥少年【网络毒理学全流程视频课】】」
点击链接直接打开
分子对接、动力学模拟课程推荐
课程一：《分子对接零基础教学课程上线 | ADT/Vina/Maestro全流程详解》
介绍：
本课程系统讲解从原理到实操的全流程操作，配套三大常用工具教学，助你快速掌握药物设计核心技能。
课程主要内容包括：分子对接原理、全流程梳理、软件下载、蛋白质小分子数据库介绍、ADT受配体结果预处理、对接口袋预测及设置、vina\qvina分子对接计算、结果可视化相互作用力分析、Maestro软件安装介绍、Maestro分子对接计算、高通量分子对接、Maestro蛋白质-蛋白质对接计算
咸鱼平台购买链接
【闲鱼】https://m.tb.cn/h.6lrKXb4?tk=EaHSVelRlCx HU293 「我在闲鱼发布了【【分子对接零基础教学课程上线 | ADT/Vina/Maes】」
点击链接直接打开
淘宝购买链接：
【淘宝】https://e.tb.cn/h.6lrv1KbfzG5szI9?tk=KQ3OVelmPFZ CZ193 「【分子对接零基础教学课程上线 | ADT/Vina/Maestro全流程详解】」
点击链接直接打开 或者 淘宝搜索直接打开
课程二：《分子动力学全流程视频课程》
介绍：
本课程系统讲解从原理到实操的全流程操作，配套三大常用工具教学，助你快速掌握药物设计核心技能。
课程总时长4小时左右，从动力学模拟原理解读、GROMACS软件安装、受配体拆分预处理、拓扑文件构建、模拟体系搭建、动力学模拟以及下游分析、结果图的绘制和解读、文章发表的撰写方法；可以说是全流程一课通！！！
咸鱼平台购买链接
【闲鱼】https://m.tb.cn/h.69TGRAA?tk=bzDYVde3bcC HU293 「我在闲鱼发布了【稳妥少年，从2022年开始主要涉及分子对接和分子动力学模拟领】」
点击链接直接打开
淘宝购买链接：
【淘宝】https://e.tb.cn/h.69hWyklEydyhI5G?tk=hvzpVdVOI6F HU293 「稳妥少年《分子动力学全流程视频课程》时长4小时左右，一课通！」
点击链接直接打开 或者 淘宝搜索直接打开
优惠活动：现在同时购买两套课程，赠送稳妥少年自研的《高通量分子对接工具HTMDv3》，赠送9113个小分子库，以及2101个重要活性分子库，助力大家进行药研相关工作！！！
工具介绍：高通量分子对接软件（HTMDv3）介绍及教学——稳妥少年
祝大家学习愉快！也感谢大家对公众号的支持，继续努力！！！今天分享的内容可能存在不足，希望感兴趣的小伙伴进行批评指正。同时也希望大家可以一键三连（点赞、关注和转发“稳妥少年”）
