摘要:很多场合需要用到两面角分布,比如轴手性分子的立体化学结构确证。本文详述了如何利用CONFLEX构象搜索的结果,获取两面角的分布,方便地对构象进行分类。

作者:肖高铿
时间:2017-11-30

一. 问题

有时候,我们需要绘制Energy Vs Dihedral angle的图以发现两面角分布与能量之间的关系。比如,在轴手性化合物的结构确证中,我们需要(1)根据某个旋转轴的两面角将构象分成不同的两类;(2)分别计算ECD或VCD。

CONFLEX进行构象搜索的优势之一是穷尽的构象搜索,多种方式评估构象能:构象能与吉布斯自由能,并可以设定溶剂环境。本文的目的是:如何利用CONFLEX的计算结果,分析化合物的两面角分布,以便将构象按指定的旋转轴两面角分类。因此,我们需要知道知道每个构象的每个两面角的角度。

二. 解决方式

1. 通过CONFLEX图形界面确认关心的两面角

CONFLEX中两面角可以通过通过4个原子定义,比如I-L-J-K。用CONFLEX的Label命令给出两面角的原子组成:CONFLEX View>Label,如图1所示。

CONFLEX教程 | 两面角的分布-墨灵格的博客

图1. 两面角通过4个原子I-L-J-K定义,比如上图的1-2-3-9。

2. 两面角信息的获取

两面角的信息写入到CONFLEX输出文件ls2里,典型的ls2文件如下:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
!====================================================================================!
 !                                                                                    !
 !  LIST OF CONFORMERS (LIST FORM 2)                                                  !
 !                                                                                    !
 !------------------------------------------------------------------------------------!
 !                                                                                    !
 !  DATE: 2009/07/06    TIME: 23:34:59.04                                             !
 !  clus-BGLU: b-glucose.mol                                                          !
 !  EMPIRICAL FORMULA:  C6H12O6     MW =   180.050                                    !
 !  FORCE FIELD:  MMFF94S(2006-11-24HGTEMP)                                           !
 !  TOTAL NUMBER OF CONFORMERS FOUND:   220                                           !
 !  STERIC ENERGY:  MIN=     61.6760     MAX=     74.6791     AVERAGE=     62.1996    !
 !                                                                                    !
 !------------------------------------------------------------------------------------!
 !                                                                                    !
 !  DEFINITIONS:                                                                      !
 !      NO. = INDEX NUMBER IN THE ORDER OF STERIC ENERGY                              !
 !      ID  = IDENTIFICATION NUMBER OF EACH CONFORMER                                 !
 !      STERIC E = STERIC ENERGY OF EACH CONFORMER (KCAL/MOL)                         !
 !      DELTA E  = RELATIVE STERIC ENERGY FROM THE GLOBAL ENERGY MINIMUM (KCAL/MOL)   !
 !      DISTRIB. = DISTRIBUTION BASED ON STERIC ENERGY (%)                            !
 !                                                                                    !
 !------------------------------------------------------------------------------------!
 !                                                                                    !
 !  DEFINITIONS OF TORSION ANGLES:                                                    !
 !      D 1 =    4 -  1 -  2 - 10                                                     !
 !      D 2 =    1 -  2 - 10 - 11                                                     !
 !      D 3 =    2 - 10 - 11 -  3                                                     !
 !      D 4 =   10 - 11 -  3 -  4                                                     !
 !      D 5 =   11 -  3 -  4 -  1                                                     !
 !      D 6 =    3 -  4 -  1 -  2                                                     !
 !      D 7 =    1 -  2 -  9 - 21                                                     !
 !      D 8 =    4 -  3 -  5 - 16                                                     !
 !      D 9 =    1 -  4 -  8 - 12                                                     !
 !      D10 =   17 -  6 - 11 -  3                                                     !
 !      D11 =   18 -  7 - 10 -  2                                                     !
 !      D12 =    4 -  8 - 12 - 24                                                     !
 !                                                                                    !
 !====================================================================================!
 
 
 ==================================================================================================================================
    NO.   ID     STERIC E   DELTA E   DISTRIB.        D 1     D 2     D 3     D 4     D 5     D 6     D 7     D 8     D 9     D10
 ----------------------------------------------------------------------------------------------------------------------------------
     1  00022     61.6760    0.0000   65.1796       -61.93   56.63  -52.15   51.91  -55.54   61.42  -70.38  172.08   59.84 -171.77
     2  00003     62.6838    1.0078   11.8954       -61.65   57.11  -52.55   51.47  -54.27   60.22  -69.49  175.75  -57.90 -171.45
     3  00001     62.9538    1.2778    7.5419       -61.87   56.00  -52.00   52.52  -56.56   62.07  -67.55  170.26  169.19 -171.29
     4  00028     63.1511    1.4750    5.4062       -61.44   55.54  -51.23   52.00  -56.12   61.56   61.52  171.80   60.57 -173.85
     5  00021     63.2235    1.5475    4.7842       -62.70   56.25  -51.72   52.20  -56.41   62.72  -64.84  171.01   67.77 -171.49
     6  00004     64.0344    2.3583    1.2174       -61.99   56.38  -52.06   51.96  -55.39   61.41  -63.44  172.54   66.33 -171.11
     7  00012     64.1467    2.4707    1.0072       -61.29   56.07  -51.60   51.53  -54.91   60.48   62.53  175.57  -57.41 -173.56
     8  00013     64.3784    2.7024    0.6811       -61.12   57.18  -52.98   51.73  -53.98   59.50  -64.04  175.09  -65.78 -170.20
     9  00002     64.5796    2.9035    0.4851       -61.58   55.20  -51.19   52.49  -56.92   62.12   58.76  170.06  169.78 -173.28
    10  00059     64.7430    3.0670    0.3681       -63.43   58.24  -52.92   50.50  -53.40   61.18  125.31  -44.21  179.82   50.08
    11  00014     64.9437    3.2677    0.2623       -61.14   55.80  -52.39   53.18  -56.71   61.52  -63.37  170.23  -69.00 -170.32
    12  00127     65.1685    3.4925    0.1795       -62.90   56.66  -52.82   53.00  -57.09   63.36  127.02  -76.29  -64.37   51.38
    13  00040     65.2332    3.5572    0.1609       -60.57   56.38  -52.36   51.81  -54.22   59.25   55.53  174.85  -65.42 -172.11
    14  00035     65.2758    3.5998    0.1498       -62.26   55.56  -51.11   52.24  -56.65   62.51   53.57  170.76   67.92 -173.30
    15  00064     65.3378    3.6618    0.1349       -59.88   53.82  -50.63   51.30  -55.00   60.47  -81.62  -44.01  179.20   49.66
    16  00131     65.6085    3.9325    0.0854       -59.19   51.96  -50.32   53.87  -59.01   62.78  -82.45  -77.61  -65.02   51.11
    17  00020     65.6576    3.9816    0.0786       -61.78   55.97  -51.49   51.79  -55.37   61.15   53.65  172.50   64.03 -172.96
    18  00068     65.9370    4.2610    0.0491       -63.15   57.70  -52.76   51.00  -54.19   61.53  118.68  -53.30  175.75   50.20
    19  00057     66.0429    4.3668    0.0410       -63.39   57.46  -52.86   51.79  -55.39   62.68  129.38  -73.35   61.37   53.41
    20  00033     66.0595    4.3834    0.0399       -60.79   55.26  -51.85   53.14  -56.75   61.22   54.07  170.05  -69.67 -172.12
    21  00074     66.2649    4.5889    0.0282       -62.02   55.58  -49.73   49.17  -53.62   60.78  -71.60  -44.87  178.93   54.56
    22  00058     66.4894    4.8134    0.0193       -63.10   57.12  -52.04   51.01  -54.77   62.07  115.35   53.54   62.60   44.72
    23  00082     66.5018    4.8258    0.0189       -59.76   52.86  -50.44   52.76  -57.41   62.18  -83.85  -76.44   60.81   53.39
    24  00088     66.5450    4.8689    0.0176       -63.29   58.11  -50.71   47.18  -50.76   59.70  -69.94  -50.05  178.14 -177.39

在DEFINITIONS OF TORSION ANGLES区块,定义了两面角I-L-J-K及其名称(D1,D2,D3...),每个两面角一行。

根据Delta E,Steric E与两面角的值,很容易就做出Energy vs 两面角分布图。在轴手性化合物研究中,就可以根据的两面角对构象分组,分别计算ECD或VCD。分组很简单了,会Linux的用户可以用shell命令直接搞定;如果不会也没有关系:复制、黏贴到Excel里,然后对两面角排序,再根据ID号选择构象进行计算。

三. 相关主题或教程

  1. CONFLEX教程 | CONFLEX构象搜索时如何约束两面角
  2. http://blog.molcalx.com.cn/2017/06/30/tutorial-conflex-torsion-fix.html

  3. CONFLEX教程 | 构象搜索
  4. http://blog.molcalx.com.cn/2016/04/12/conflex-conformation-search.html

四. 联系我们