纯粹是1个突发奇想的问题最近看到了几篇公众号看到了CRISPR拿诺奖的也开始做AI蛋白质设计了想想看利用分子动力学MD模拟能否研究DNA碱基变异对基因编辑如CRISPR-Cas/碱基编辑器效率的影响。模拟引擎的话选gromacs或amber生态都没有太大问题一、 力场选择Force Field Selection基因编辑复合体包含蛋白质、RNA和DNA三组分这一点可以参考基因编辑维基百科https://en.wikipedia.org/wiki/CRISPR_gene_editinghttps://en.wikipedia.org/wiki/Genome_editing对力场的“多组分兼容性”要求极高。AMBER 软件生态首选推荐推荐力场组合ff19SB蛋白质 OL3RNA OL15 或 bsc1DNA。优势AMBER 在核酸力场尤其是 DNA 的 bsc1 和 RNA 的 OL3的参数修正上处于领先地位。它能准确模拟大范围核酸非特异性扭曲、碱基堆积以及错配碱基的动态构象非常适合研究碱基变异。GROMACS 软件生态推荐力场组合CHARMM36m蛋白质与核酸其实就是全原子力场。优势GROMACS 的计算速度极快。CHARMM36m 对核酸的主链扭角和核糖构象进行了优化同样能够胜任此类模拟但在处理某些极其特殊的碱基错配时核酸刚性可能略高。二、 核心分析如何验证机制又如何照例是通过上述力场运行出分子动力学轨迹Trajectory后进行后续深入模块分析1. 结构稳定性与刚性分析常用分析工具RMSD均方根偏差、RMSF均方根涨落。具体能做什么计算整个 Cas蛋白-gRNA-DNA 复合体随时间变化的位移以及每个氨基酸/核苷酸的柔性。分别验证什么验证变异是否破坏整体结合如果靶区碱基变异导致系统 RMSD 持续飙升且无法收敛说明该变异彻底破坏了编辑复合体的结构稳定性。验证局部松动若变异位点附近的 RNA/DNA 链 RMSF 显著增大说明变异导致了局部“脱靶”或解离趋势揭示效率下降的动力学原因。2. 核酸构象与空间几何分析常用分析工具Curves 或 3DNA 软件定量分析核酸几何参数。具体能做什么提取模拟轨迹中 DNA 的扭曲角Twist、倾斜角Tilt、滚转角Roll以及沟宽Groove width。分别验证什么验证“非典型变变异”的空间位阻碱基变异如嘌呤变嘧啶会改变 DNA 双螺旋的局部几何形状。通过此分析可验证变异是否导致 DNA 沟槽变窄从而阻止了 Cas 蛋白特异性氨基酸侧链如精氨酸 R、赖氨酸 K的插入和识别。3. 结合自由能计算热力学验证常用分析工具MM-GBSA / MM-PBSA分子力学-广义玻尔兹曼/泊松-玻尔兹曼表面积法。具体能做什么定量计算 Cas蛋白-gRNA 复合物与目标 DNA 之间的总结合自由能Δ G b i n d \Delta G_{bind}ΔGbind并将其分解到每一个氨基酸和碱基上能量贡献分解。分别验证什么验证编辑效率的定量差异比较野生型 DNA 与变异型 DNA 的Δ G \Delta GΔG。Δ G \Delta GΔG越负结合越牢固。如果变异型Δ G \Delta GΔG显著升高亲和力下降直接解释了实验中编辑效率低的原因。识别关键相互作用丢失能量分解能具体找出是哪个残基与变异碱基之间的氢键或静电作用消失了为设计改造 Cas 蛋白以逆转效率下降提供靶点。4. 催化活性中心的动态网络分析常用分析工具氢键网络分析Hydrogen Bond Analysis、全原子接触图Contact Map。具体能做什么监测碱基编辑器中的脱氨酶如 ABE/CBE或 Cas 剪切位点与待编辑单碱基之间的距离涨落及氢键寿命。分别验证什么验证催化活性是否被抑制即使复合物能结合如果碱基变异导致待编辑碱基偏离了脱氨酶的催化活性口袋例如催化距离从正常的3.5 A ˚ 3.5\text{ \AA}3.5A˚增大到了6 A ˚ 6\text{ \AA}6A˚以上就能验证该变异是通过破坏催化构象而非破坏结合来降低编辑效率的。三、 进阶验证增强采样Enhanced Sampling由于常规分子动力学微秒级别很难观测到完整的 DNA解链、gRNA与DNA链置换R-loop形成 过程通常我们会使用增强采样技术深入伞状采样Umbrella Sampling或元动力学Metadynamics验证什么强制将 gRNA 与变异 DNA 一步步拉开或结合计算完整的解离自由能能垒Energy Barrier。这能从动力学速率Kinetic Rate层面彻底验证碱基变异是如何在时间尺度上阻碍编辑复合体正确组装的。