yy.vip易游-Gaussian教程 计算分子的比旋光度

　　YYVIP易游·(中国有限公司)官方网站-

　　GIST 结合位点水分析在先导化合物优化中的应用 —— 以 15-PGDH 抑制剂设计为例

　　静电互补性分析：一种高效替代FEP的配体优化策略——以MALT1别构抑制剂为例

　　Flare™ V11 正式发布：新增生成式 AI 工具、约束的蛋白质-蛋白质对接功能、分离计算模式以及众多其他增强功能

　　摘要：本教程介绍了手性分子的旋光度、比旋光度、旋光色散等相关概念。并以甲基环氧乙烷分子为例，详细介绍了如何利用Gaussian程序计算分子的比旋光度，并通过与实验测量值进行对比以确定分子的绝对构型。

　　旋光度（Optical Rotation）是指当一束偏振光经过某类物质时，其相对于自身运动方向旋转的角度。而物质具有的这种使偏振光旋转的性质也称为旋光性，具有旋光性的物质包括晶体，自旋极化分子和手性分子。

　　对于同一种分子，当入射的偏振光波长发生变化时，观察到的旋光度也会相应发生变化，这种变化也称为旋光色散。

　　通常情况下，由于测量时，旋光管，溶液浓度，温度等条件有所不同，因此对于同一种分子，在相同波长的入射偏振光下测得的旋光度也不相同。为了便于比较，通常把在不同条件下测得的旋光度换算成单位长度旋光管和单位浓度下的旋光度值，这个值也称为比旋光度，其单位为degrees/(dm·g/ml)。它是物质的固有属性，不随测量条件的变化而变化。

　　其中，[a]为比旋光度，a为旋光度，C为浓度(单位为克每100毫升)，L为光通过旋光管的光程（单位dm）。在比较比旋光度的时候需要指明测量时的温度和入射光波长。

　　Gaussian程序提供了2个关键字（连同选项）：Polar=OptRot（指示程序计算旋光度），CPHF=RdFreq（指示程序读取特定的波长）完成特定入射波长下旋光度的计算，计算结果即为比旋光度。

　　需要注意的是：由于比旋光度受溶剂的影响较大，因此计算一定要考虑溶剂的影响，并且在较大的基组下进行计算。

　　上述我们提到，手性分子具有旋光性，而对于一对对映体，其旋光的方向正好相反，因此，通过与实验测得的比旋光度相比较，能够确定手性分子的绝对构型。

　　在本教程中，我们以甲基环氧乙烷作为例子，通过计算两种手性分子构型的比旋光度，并与实验值对比，以确定分子的绝对构型。实验数据我们选用气相条件下，在355nm和633nm下测得的比旋光度。

　　使用第一步优化得到的结构作为初始结构，计算S和R型甲基环氧乙烷的比旋光度。计算S和R型甲基环氧乙烷在气相条件下，355nm和633nm入射光下比旋光度的输入文件分别如下。

　　如果上一步优化步骤已经保存了check文件，比如opt.chk,则可以简化为下面的输入文件：

　　利用GaussView查看分子的比旋光度，我们以S构型甲基环氧乙烷比旋光度计算结果为例，步骤如下：

　　在Results下拉菜单中点击ORD选项，将显示不同入射波长下的比旋光度。用鼠标点击图中的点，在下方将显示指定波长下比旋光度。

　　除了用GaussView查看外，还可以直接从输出文件里查看，找[Alpha]行就可以，比如下面这个样子:

　　利用同样的方法我们可以查看每个构型在不同入射频率下的的比旋光度，下表列出了所有计算和实验得到的比旋光度，可以看出在这个例子中，S构型计算得到的比旋光度与实验测得的更为接近，因此我们认为实验测量的分子为S构型，这也与事实一致。

　　需要注意的是：对于柔性分子，需要进行构象搜索，对计算结果进行玻尔兹曼平均，具体方法参见NMR计算教程。