环puckering坐标分析工具:gmx_puckering

类别:    标签: gmx   阅读次数:   版权: (CC) BY-NC-SA

学过化学的人大致都记得, 环己烷是六元环, 它有两种主要的构象, 船式和椅式. 这两种构象之间是可以相互转化的. 我还记得, 上课时老师还拿模型给我们演示过, 怎么通过扭转碳原子将一种构象转换为另一种. 后来我自己好像也试过, 转换是可以转换过去, 但这个过程中各个原子具体怎么动的, 我却缺乏直观的想象. 那时也曾想过, 或许应该能用一些参数来表征这种转换过程吧. 至于到底要用什么样的参数, 如何表征, 因为和自己的工作关系不大, 也就没有上心去查查证思考了.

现在我知道了, 这个问题早在我出生之前就已经解决了, 尽管我上大学时还对其一无所知. 对五元环或六元环的构象进行分析时, 常用的一类坐标是puckering坐标(或可译为折皱/褶皱/摺皱/扭折坐标, 就是不平整, 皱皱巴巴, 咕咕嘟嘟之类的意思). 这类坐标最初由Cremer和Pople提出(D. Cremer, J. A. Pople; J. Am. Chem. Soc. 97(6):1354-1358, 1975; 10.1021/ja00839a011), 所以在文献中常被称为CP坐标. 具体的含义和算法可阅读原文以及一些引用文献(如Lucian Chan, Geoffrey R. Hutchison, Garrett M. Morris; J. Chem. Inf. Model. 61(2):743-755, 2021; 10.1021/acs.jcim.0c01144).

现在, 我们很容易对各种含环的体系进行分子动力学模拟, 其中的糖类分子近年来很受关注. 而糖类分子中常见的是五元环的呋喃糖, 六元环的吡喃糖. 葡萄糖大概是吡喃糖中最为人熟知的吧. 要对模拟轨迹中糖环的构象进行分析, 我们可以计算CP坐标, 然后根据CP坐标, 将环的构象归结到IUPAC定义的38种正则构象之一.

既然想计算CP坐标, 那就要先调研一下, 看有没有现成的工具, 实在没有的话, 再自己想办法. 搜索一圈, 发现

这类坐标的另一个用途是反过来, 根据相应的坐标生成指定构象的分子结构, 或者用于描述环的转变过程, 并计算其自由能形貌图(或曰景观图). 这方面的论文较多, 随意列出几篇吧:

看了一圈之后, 选择很多, 但各有利弊, 可行是可行的, 完美是不可能完美的. 完美的工具只能出自自己之手, 所以只得撸起袖子, 将g_puckering这个年久失修的工具拾捡起来, 擦拭擦拭, 看看能不能修复一下, 物尽其用. 一则可以达到最终目的, 二则可以熟悉如何基于GROMACS已有框架编写分析工具. 断断续续, 磕磕绊绊磕磕, 也算完成了. 因为这个工具是基于GROMACS2019.6的, 所以就唤它作gmx_puckering.

待改进

◆本文地址: , 转载请注明◆
◆评论问题: https://jerkwin.herokuapp.com/category/3/博客, 欢迎留言◆


前一篇: Eigenmath函数说明
后一篇: gmx_puckering使用示例

访问人次(2015年7月 9日起): | 最后更新: 2022-12-17 03:04:57 UTC | 版权所有 © 2008 - 2022 Jerkwin