Saen-Oon Suwipa, Hannongbua Supa, Wolschann Peter
Chemistry Department, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand.
J Chem Inf Comput Sci. 2003 Sep-Oct;43(5):1412-22. doi: 10.1021/ci0340299.
The conformational analysis of the HIV-1 reverse transcriptase inhibitor, (+)-(s)-4,5,6,7-tetrahydro-9-chloro-5-methyl-6-(3-methyl-2-butenyl)imidazol[4,5,1-jk][1,4]benzodiazepin-2(1H)-thione or 9-Cl TIBO, has been investigated using high level of calculations, ab initio, and DFT theory. The potential energy surface as the function of two important rotatable dihedral angles of the 9-Cl TIBO side chain was generated by the Hartree-Fock method at the 3-21G basis set. Eight pronounced local minima were found to exist within an energy difference of less than 10 kJ/mol. The energy barriers between the different local minima are lower than 15 kJ/mol. A second derivative (frequency) analysis showed that all conformers are stable at this level of theory. These structures were used as starting points for full geometry optimizations at the HF/6-31G and B3LYP/6-31G levels of theory to obtain the absolute geometries and structural information. The comparisons of calculated conformers with the bound conformer in the X-ray structure were sequentially considered. Additionally, to obtain some structural information and to correlate between calculated structures and the structure in solution, NMR chemical shift calculations were also performed on all eight local minimum structures at B3LYP/6-311++G level, using the GIAO approach. The calculated (1)H NMR and (13)C NMR chemical shifts for the lowest energetic conformer give the greatest correspondence with the experimental results.
已使用高水平计算、从头算和密度泛函理论(DFT)对HIV-1逆转录酶抑制剂(+)-(s)-4,5,6,7-四氢-9-氯-5-甲基-6-(3-甲基-2-丁烯基)咪唑并[4,5,1-jk][1,4]苯并二氮杂䓬-2(1H)-硫酮或9-Cl TIBO进行了构象分析。通过Hartree-Fock方法在3-21G基组下生成了9-Cl TIBO侧链两个重要可旋转二面角的势能面。发现在能量差小于10 kJ/mol的范围内存在八个明显的局部极小值。不同局部极小值之间的能垒低于15 kJ/mol。二阶导数(频率)分析表明,在该理论水平下所有构象异构体都是稳定的。这些结构被用作在HF/6-31G和B3LYP/6-31G理论水平下进行全几何优化的起点,以获得绝对几何结构和结构信息。依次考虑了计算得到的构象异构体与X射线结构中结合构象的比较。此外,为了获得一些结构信息并关联计算结构与溶液中的结构,还使用GIAO方法在B3LYP/6-311++G水平下对所有八个局部极小值结构进行了NMR化学位移计算。能量最低的构象异构体的计算(1)H NMR和(13)C NMR化学位移与实验结果的对应度最高。