Sapse A M, Feng W, Fugler-Domenico L, Kabir S, Joseph T, Lown J W
Department of Physiology and Biophysics, Mount Sinai School of Medicine, New York, N.Y.
J Biomol Struct Dyn. 1993 Feb;10(4):709-26. doi: 10.1080/07391102.1993.10508002.
Ab initio calculations (Hartree-Fock) using the 6-31G basis set have been performed on two chiral oligopeptide antitumor antibiotics amidinomycin 5 and noformycin 6. The latter are DNA minor groove binding agents related to the A.T recognizing netropsin 4 and distamycin 3 but, unlike the latter, bear stereocenters (two for 5 and one for 6) that may be expected to affect binding to the B-DNA receptor. Geometry optimized conformations, energies and distribution of electrostatic charges within the molecules were derived. The rotational barrier for bond C3-C6 in 6 was calculated to be ca. 6 kcal.mole-1 and the dipole moment for 6 was 7.69D and for 5 was 5.58D. The ab initio derived parameters of the geometry optimized conformations of the different possible stereoisomeric forms of 5 and 6 were used to interpret their different interactions with the minor groove of DNA at both A.T and G.C sequences and the results were compared with molecular mechanics calculations. The order of binding of the four stereoisomers of 5 at the preferred (A.T)n sequences by both ab initio and molecular mechanics calculations is 1S,3R > RR > RS > SS. The predicted energy differences for complexation with DNA of the other stereoisomers from that of 1S,3R are: RR (4.2%); RS (6.7%) and SS (21.5%). In the case of noformycin the 4R structure binds more effectively than the enantiomer. Considerations of phasing in the computed distances between hydrogen bond donating sites in the DNA-bound antibiotics provide further insight into the binding processes. In the complexes of noformycin 6 the N-N1-N4 and N1-N5 distances (9.05 and 9.15 A respectively for 4R-6 and 9.23 and 9.26 A respectively for 4S-6) are close to the optimum value of 9.1 A for effective binding. In the case of amidinomycin 5 the best agreement with the optimum value occurs with the strongest binding diastereomer 1S,3R (N1-N3 = 8.91, N1-N4 = 9.41 A). The unexpected result, consistent in both ab initio and molecular mechanics treatments, is that, in contrast to the cases of kikumycin 1 and anthelvencin 2, the natural 3S configuration of 5 and 4S of 6 do not confer maximal binding efficiency. This suggests that biogenetic factors in the generation of the oligopeptide antibiotics lead to maximum DNA binding in the cases of kikumycin and anthelvencin but not in the cases of amidinomycin and noformycin.
使用6-31G基组进行了从头算(哈特里-福克)计算,研究对象为两种手性寡肽抗肿瘤抗生素脒霉素5和诺福霉素6。后者是与识别A.T的纺锤菌素4和偏端霉素3相关的DNA小沟结合剂,但与后者不同的是,它们带有立体中心(5有两个,6有一个),预计这些立体中心会影响与B-DNA受体的结合。得出了分子内几何结构优化的构象、能量和静电荷分布。计算得出6中C3-C6键的旋转势垒约为6 kcal·mol⁻¹,6的偶极矩为7.69D,5的偶极矩为5.58D。利用从头算得出的5和6不同可能立体异构体形式的几何结构优化构象参数,来解释它们在A.T和G.C序列处与DNA小沟的不同相互作用,并将结果与分子力学计算结果进行比较。从头算和分子力学计算得出,5的四种立体异构体在优选的(A.T)n序列处的结合顺序为1S,3R > RR > RS > SS。与1S,3R相比,其他立体异构体与DNA络合的预测能量差异为:RR(4.2%);RS(6.7%)和SS(21.5%)。就诺福霉素而言,4R结构比其对映体结合更有效。对DNA结合抗生素中氢键供体位点之间计算距离的相位考虑,为结合过程提供了进一步的见解。在诺福霉素6的络合物中,N-N1-N4和N1-N5距离(4R-6分别为9.05和9.15 Å,4S-6分别为9.23和9.26 Å)接近有效结合的最佳值9.1 Å。就脒霉素5而言,与最佳值的最佳一致性出现在结合力最强的非对映异构体1S,3R中(N1-N3 = 8.91,N1-N4 = 9.41 Å)。在从头算和分子力学处理中都一致的意外结果是,与菊霉素1和抗蠕虫菌素2的情况相反,5的天然3S构型和6的4S构型并未赋予最大结合效率。这表明,寡肽抗生素生成过程中的生物发生因素在菊霉素和抗蠕虫菌素的情况下导致最大的DNA结合,但在脒霉素和诺福霉素的情况下并非如此。