Yurenko Yevgen P, Zhurakivsky Roman O, Ghomi Mahmoud, Samijlenko Svitlana P, Hovorun Dmytro M
Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, vul. Zabolotnoho 150, 03143, Kyiv, Ukraine.
J Phys Chem B. 2007 Jun 7;111(22):6263-71. doi: 10.1021/jp066742h. Epub 2007 May 16.
A comprehensive conformational analysis of isolated 2'-beta-deoxy-6-azacytidine (d6AC), an analogue of therapeutically active 6-azacytidine (6AC), has been performed by means of ab initio calculations at the MP2/6-311++G(2df,pd)//DFT B3LYP/6-31G(d,p) level of theory. Among the 81 conformers located within a 7.83 kcal/mol Gibbs energy range at T = 298.15 K, 38 contain syn-oriented bases with respect to 2'-deoxyribose; the other conformers include anti-oriented bases. Energetic analysis of these conformers shows that conformational equilibrium of isolated d6AC at T = 298.15 K is shifted to syn conformation with a syn/anti ratio estimated as 61.4%:38.6%. As far as the sugar conformation is concerned, 40 conformers contain north (N) (with 0.3 degrees < or = P < or = 40.1 degrees), and the rest possess south (S) (with 157.1 degrees < or = P < or = 207.0 degrees) puckers, where P is the pseudorotational angle of the furanose ring. The S/N occupancy ratio is estimated as 80.2%:19.8% (T = 298.15 K). The two most stable conformers are energetically quasidegenerate and correspond to both C2'-endo/syn conformers differing only by orientation of the O3'H hydroxyl group. They are both stabilized by means of similar intramolecular H-bonds, i.e., O5'H...O2, C2'H2...O2, and C2'H2...O5'. As examined by AIM criteria, from 1 to 3 H-bonds per conformer were identified among 13 possible interactions: O5'H...O2, O5'H...N6, O3'H...O5', O5'H...O3', C1'H...O2, C2'H2...O2, C2'H2...O5', C3'H...O2, C3'H...N6, C5'H1...O2, C5'H2...O2, C5'H1...N6, and C5'H2...N6. The biological effect of d6AC is conceived as an inhibition of replicative DNA polymerase caused by an unusual orientation of the sugar residue against the base in the only A form DNA-like conformer.
通过在MP2/6 - 311++G(2df,pd)//DFT B3LYP/6 - 31G(d,p)理论水平下的从头算计算,对治疗活性6 - 氮杂胞苷(6AC)的类似物——分离的2'-β-脱氧-6 - 氮杂胞苷(d6AC)进行了全面的构象分析。在T = 298.15 K时,吉布斯自由能范围为7.83 kcal/mol内找到的81个构象异构体中,38个含有相对于2'-脱氧核糖呈顺式取向的碱基;其他构象异构体则包含反式取向的碱基。对这些构象异构体的能量分析表明,在T = 298.15 K时,分离的d6AC的构象平衡向顺式构象偏移,顺式/反式比例估计为61.4%:38.6%。就糖的构象而言,40个构象异构体具有北(N)型(0.3°≤P≤40.1°),其余的具有南(S)型(157.1°≤P≤207.0°)褶皱,其中P是呋喃糖环的假旋转角。S/N占有率估计为80.2%:19.8%(T = 298.15 K)。两个最稳定的构象异构体在能量上几乎简并,对应于仅O3'H羟基取向不同的两个C2'-内型/顺式构象异构体。它们都通过类似的分子内氢键,即O5'H...O2、C2'H2...O2和C2'H2...O5'而稳定。根据AIM标准检查,在13种可能的相互作用:O5'H...O2、O5'H...N6、O