金属丝桃蒽酮结构变化的理论研究

用MMX和AM1方法对金属桃蒽酮(HYP)结构变化过程的生成热进行了计算.发现:1.HYP1,2位OH可越过20kJ/mol左右的势,绕C--O键旋转而形成分子内氢键,并估算出键能约为10kJ/mol;2.HYP其他四个OH也可进行类似的构象变化,势垒在26kJ/mol左右,相应分子内氢键键能约为20kJmol;3.HYP分子内氢传递产生的异构体在能量上不稳定,进而在基础上探讨了1,2位OH分子内氢键的形成对HYP光敏活性的影响。     

痂囊腔菌素A分子构型和分子内氢键的理论研究

用半经验量子化学AM1方法对天然苝醌化合物痂囊腔菌素A(EA)的分子构型和分子内氢键进行了研究;从EA可能的64种构型中选择16种进行了计算.结果表明,EA的X射线晶体结构对应的构型是II型左旋(A)a,a型(II-L-A-a,a);小的生成热差值可以使得异构体间的转换容易进行,有利于发生分子内质子传递.尽管采用AM1方法计算得到的EA各种构象的平面性有所差别,但都很接近晶体的平面性.此外,EA分子内氢键键能的平均值为22.9kJ/mol;II型的氢键键能比I型的大,(9,10)位的氢键键能比(3,4)位的大;EA的平面性是由苝醌环上的侧链取代所决定,而与分子内是否存在氢键无关.     

磷酰化丝氨酸形成六配位磷中间体的理论研究

用MNDO方法对磷酰化丝氨酸仿生化反应机理中六配位磷中间体的形成过程进行了研究.磷酰化丝氨酸(1)形成分子内磷酸-羧酸分子内混酐的五配位磷中间体(2)后,其酸性质子解离,分子经过具有氢桥键结构的过渡态,使氨基酸侧链羟基上的氢通过氢键作用向磷上的O₁进行转移,然后再经过构型由三角双锥向八面体的转变,形成六配位磷中间体(3).氢桥键的存在使反应过渡态能量降低,其相对能量为148.5kJ/mol.理论计算较成功的解释了六配位磷中间体的形成机理以及磷酰化丝氨酸仿生化反应中羧基和侧链羟基共同参与的实验结果.     

SBA-15分子筛羟基结构性质及其催化活性的理论研究

采用密度泛函方法对SBA-15分子筛簇模型化合物进行了计算模拟,重点研究了簇模型化合物羟基的几何构型和电子结构性质。从键级、前线分子轨道、静电势和质子化能等方面探讨了表面羟基的酸碱及氧化还原性质。分子筛模型化合物的静电势图显示SBA-15分子筛表面孤立羟基的H原子为分子筛表面的L酸性位,而氢键羟基的O原子则为分子筛表面的B碱性位。分子筛的前线分子轨道研究发现,表面羟基是分子筛表面的氧化还原活性位。计算得到的羟基SiO-H键级及形成的氢键键级分别在0.677 5～0.710 5和0.055 7～0.092 6范围,计算羟基OH的质子化能在1 471～1 589 kJ.mol-1范围。考察分子筛表面的质子化能显示未参与形成氢键的羟基H质子具有较强的B酸性。     

苝醌类光敏剂构象的分子力学研究

用分子力学方法,计算了苝醌类光敏剂竹红菌甲素(HA),乙素(HB)等不同构象的生成热及酚羟基质子解离前后生成热的变化.发现:(1)1(1)HA,HB4种构象导构体的生成热相近,因而室温下即可相互转化;(2)2(2)HA,HB酚羟基质子解离能力有差别,HA强于HB.(3)3(3)HA,HB存在两个分子内氢键,键能均在8kJ/mol左右,而且Ⅰ型构象键能低于Ⅱ型构象,HA低于HB;(4)4(4)当酚羟基键从形成氢键位置扭转180°左右时存在一能量低点,提示这种构象也是可能存在的;(5)5(5)HBMC的酚羟基氢与相邻的氮原子之问也形成分子内氢键,这是HBMC仍有光敏活性的结构基础.     

铑(Ⅰ)催化烯烃碳酰化反应机理:手性聚稠环过程的对映选择性

采用密度泛函理论方法研究了[Rh(R,R-DIOP)]+[DIOP=(1R,2R)-1,2-O-异丙叉-1,2-二醚氧基-1,2-双(二苯基磷基)乙烷]催化下苯并环丁酮手性聚稠环过程在气相、四氢乙呋喃(THF)及水中的反应机理.计算结果表明,在气相中反应容易进行,经TS2形成六元环的过程为决速步骤,但产物无明显的对映选择性.在THF中,S-和R-通道的C—C键活化能垒仅由79.5和69.3 kJ/mol提高到80.2和88.8 kJ/mol,未改变反应的实质;Rh与2个C原子的配位明显弱于气相,相对于催化剂和反应物自由能之和,S-和R-通道的反应总能垒分别提高到63.8和68.1 kJ/mol.对于S-通道,溶剂THF使经TS2的能垒升至112.0 kJ/mol,仍为整个过程的决速步骤;然而对R-通道,溶剂使经TS1形成五元环过程的能垒升至91.5 kJ/mol,但使经TS2的能垒由78.9 kJ/mol降至77.7 kJ/mol,IM1→TS1成为决速步骤.在以水为溶剂时,经TS1形成五元环的过程成为2个通道的决速步骤,其在S-和R-通道中的能垒分别为102.5和94.9 kJ/mol.因此,溶剂改变了反应的决速步骤及能垒.3种方法均预测R-通道为主反应路径,但THF能明显增加产物的对映选择性.采用自然键轨道(NBO)电荷分析了反应过程中电荷的变化.     

苄嘧磺隆印迹聚合物的波谱分析及吸附性能研究

以苄嘧磺隆为模板分子, α-甲基丙烯酸为功能单体, 三甲氧基丙烷三甲基丙烯酸酯为交联剂, 二氯甲烷为致孔剂, 在20 ℃温度下, 采用紫外引发沉淀聚合方法制备了苄嘧磺隆分子印迹聚合物. 紫外光谱和核磁共振氢谱实验提示了聚合前模板分子与功能单体之间的相互作用主要是分子间氢键, 分子间氢键相互作用能和双氢键的键距由Hyperchem 7.0和Gaussian 03W软件分别计算为: －28.6163 kJ/mol和0.179～0.181 nm. 制备的印迹聚合物在高效液相色谱和吸附动力学实验中都表现出对模板分子有较强的吸附作用.     

氯化1-乙基-3-甲基烷基咪唑分子结构及其氢键作用能的密度泛函研究

采用量子化学密度范函方法计算研究了氯化1-乙基-3-甲基咪唑盐（[EMIM]C1）正负离子间相互作用的平衡构型和Cl^-离子从不同方位与咪唑阳离子的结合能．计算发现Cl^-可在咪唑环平面上结合形成离子键，其结合能为．368．97kJ／mol，还可与环上的三个H原子形成氢键，其氢键键长和结合能分别为0．2000nm／-378．03kJ／mol、0．2339nm／-344，69kJ／mol和0．2301nm／-346，01kJ／mol．结合能包括氢键作用能和正负离子的静电作用能，后者为主要贡献．从理论上展示了该离子液体的氢键结构．     

儿茶素-胞嘧啶分子间相互作用的密度泛函研究

采用密度泛函理论的B3LYP方法,在6-31+G*基组水平上研究了儿茶素-胞嘧啶分子间相互作用机制,得到稳定的儿茶素-胞嘧啶复合物11个.计算结果表明氢键对于复合物的稳定性起着重要的作用,并且当复合物形成2个或更多的氢键时,氢键的类型及强度共同决定着复合物的稳定性.我们还应用了分子中的原子(AIM)理论和自然键轨道(NBO)理论对这11种复合物中氢键的性质和特征进行了分析.通过研究发现,所有的氢键复合物进行基组重叠误差(BSSE)校正后的相互作用能为-17.35～-43.27kJ/mol,相互作用能主要由氢键所贡献.振动分析显示,氢键的形成使相对应键的对称伸缩振动频率减小,说明这些复合物中形成的氢键都是正常的红移型氢键,与实验结果相一致.     

苏丹红Ⅱ与肌红蛋白相互作用的分子模拟与实验

用荧光光谱法研究了苏丹红Ⅱ与肌红蛋白(Mb)之间的相互作用, 实验结果表明,二者结合位点数近似为1, 结合常数K=3.84×107L/mol, 有很强的相互作用. 用分子柔性对接技术模拟确定了它们之间的作用位点、作用力类型及相互作用能. 理论计算的结果表明,苏丹红Ⅱ和Mb相互作用的势能为-9419.9 kJ/mol, 静电能为-7468.8 kJ/mol, 范德华能为-1951.0 kJ/mol. 苏丹红Ⅱ与Mb中His64残基形成氢键, 苏丹红Ⅱ也能与疏水氨基酸残基, 如能产生内源荧光的Phe33、Phe43、Phe106 和 Phe138等发生作用, 这与苏丹红Ⅱ能使Mb荧光猝灭的实验结果是一致的.     

Theoretical elucidation of structure-activity relationship of flavonoid antioxidants

AM1 method was employed to calculate flavonoid antioxidants, and the results obtained are as follows. Firstly, flavonoid hydroxyls at ortho position were more active than the hydroxyls at meta position in scavenging oxygen-free radicals, which resulted from the facts that (ⅰ) the former were stabilized by forming intramolecular hydrogen bond and (ⅱ) ortho benzoquinone formed in the former structures through resonance, which resulted in large percentage of distribution of spin density on ortho oxygen and low internal energy. Secondly, electron-attracting effect of ring C of chromone-flavonoids showed some passive effects on hydroxyls of ring A, making the OH less active. As ring C had little effect on ring B and hydroxyls of ring B in most flavonoids were at ortho position, the rule summarized from experiments showing that hydroxyls of ring B were more active in scavenging oxygen-free radicals was elucidated.     

Investigation of the influence of hydroxy groups on the radical scavenging ability of polyphenols

Recently, O-H bond dissociation enthalpies (BDEs) have been successfully used to express the free radical scavenging ability of polyphenolic antioxidants. In this work, the BDEs of phenol, catechol, resorcinol, hydroquinone, pyrogallol, phloroglucinol, 1,2,4-benzenetriol, and 5-hydroxypyrogallol have been calculated at B3LYP/6-311G++(3df, 3pd) and used to elucidate the effect of OH groups. Increasing the number of OH groups in the adjacent (vicinal) position decreases the BDE of phenols. Increasing the number of O-H groups in the alternative position C(1,3) as in resorcinol and C(1,3,5) as in phloroglucinol does not show any notable change in the BDEs when compared to that of OH in C(1) as in phenol. 5-Hydroxypyrogallol has the smallest BDE (250.3 kJ mol(-1)) followed by pyrogallol (289.4 kJ mol(-1)), then 1,2,4-benzenetriol (294.8 kJ mol(-1)), and then catechol (312.8 kJ mol(-1)). Overall, our results indicated that the presence of ortho and para hydroxy groups reduces the BDEs. An intramolecular hydrogen bond (IHB) develops due to the ortho arrangement of OH's and plays a dominant role in decreasing the BDEs. This key study on phenols showed that the reactive order of OH position in the benzene ring is the following: 5-hydroxypyrogallol > pyrogallol > 1,2,4-benzenetriol > catechol > hydroquinone > phenol approximately resorcinol approximately phloroglucinol.     

Intramolecular hydrogen bonding and anion binding of N-benzamido-N'-benzoylthioureas

N-(p-Dimethylamino)benzoyl-N'-phenylthiourea as an N-acylthiourea is known to be unable to bind anions due to a strong intramolecular hydrogen bond (IHB). We show here that by inserting an amido group in the N'-phenyl side the newly designed N-benzamido-N'-benzoylthioureas, despite this IHB too, bind strongly to anions with binding constants on the order of 10(6)-10(7) mol(-1) L. Results suggest that potential anion receptors or organocatalysts could be developed on the basis of this framework with a wide structural diversity.     

Solvent and substituent effects on the conformational equilibria and intramolecular hydrogen bonding of 4-substituted-2-hydroxybenzaldehydes

By applying the B3LYP/6-31G(d) method with the SCIPCM model on seven 4X substituted 2-hydroxybenzaldehydes, some structural characteristics related with their conformational equilibria and intramolecular hydrogen bonds have been clarified. The compounds are almost completely under the planar conformation characterized by a strong intramolecular hydrogen bond, which decreases in those solvents that possess a higher hydrogen bond donating capability and polarity. The substituents exert a marked influence on the conformational equilibrium constants and the strength of the IHB. Moreover, the excellent Hammett-type equations obtained support the proposed conformational reactions to quantify the IHB in the o-hydroxybenzaldehydes studied.     

The presentation of an approach for estimating the intramolecular hydrogen bond strength in conformational study of β-Aminoacrolein

All the plausible conformations of β-aminoacrolein (AMAC) have been investigated by the Bekes-Lee-Yang-Parr (B3LYP) nonlocal density functional with extended 6-311++G** basis set for studying the stability order of conformers and the various possibilities of intramolecular hydrogen bonding formation. In general the ketoamine (KA) conformers of AMAC, by mean average, are more stable than the corresponding enolimine (EI) and ketoimine (KI) analogues and this stability is mainly due to the π-electron resonance in these conformers that established by NH2 functional group. The contribution of resonance to the stability of chelated KA conformers is about 75.6 kJ/mol, which is greater than that of the hydrogen bond energy (E_HB=35.0 kJ/mol). The relative decreasing order of the various hydrogen bond energies was found to be: O–HN_imine(strong)>N_amine-HO_keto (normal)>N_imine-HO_hydroxyl (weak) > N_imine-HO_keto (weak). Hydrogen bond energies for all systems were obtained from the method that we called related rotamers method (RRM). The topological properties of the electron density contributions for various type of intramolecular hydrogen bond have been analyzed in term of the Bader theory of atoms in molecules (AIM). The results of these calculations support the previous calculations, which obtained by the related rotamer methods.     

Nitro-substituted 2-phenyl- and 2-(2-tosylaminophenyl)-4H-3,1-benzoxazin-4-one

Nitro-substituted 2-phenyl- and 2-(2-tosylaminophenyl)-4H-3,1-benzoxazin-4-ones were synthesized. The UV, IR, and luminescence spectra were studied. The position of the nitro group affects the strength of the intramolecular hydrogen bond (IHB). The luminescence properties of nitro-substituted 2-(2-tosylaminophenyl)-4H-3,1-benzoxazin-4-ones are associated with the strength of the IHB. The luminescence maximum is shifted to the short-wave region as the IHB becomes stronger.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 5, pp. 616–621, May, 1972.     

Gas-phase structure of protonated histidine and histidine methyl ester: combined experimental mass spectrometry and theoretical ab initio study

Gas-phase H/D exchange experiments with CD3OD and D2O and quantum chemical ab initio G3(MP2) calculations were carried out on protonated histidine and protonated histidine methyl ester in order to elucidate their bonding and structure. The H/D exchange experiments show that both ions have three equivalent fast hydrogens and one appreciably slower exchangeable hydrogen assigned to the protonated amino group participating in a strong intramolecular hydrogen bond (IHB) with the nearest N(sp2) nitrogen of the imidazole fragment and to the distal ring NH-group, respectively. It is taken for granted that the proton exchange in the IHB is much faster than the H/D exchange. Unlike in other protonated amino acids (glycine, proline, phenylalanine, tyrosine, and tryptophan) studied earlier, the exchange rate of the carboxyl group in protonated histidine is slower than that of the amino group. The most stable conformers and the enthalpies of neutral and protonated histidine and its methyl ester are calculated at the G3(MP2) level of theory. It is shown that strong intramolecular hydrogen bonding between the amino group and the imidazole ring nitrogen sites is responsible for the stability and specific properties of the protonated histidine. It is found that the proton fluctuates between the amino and imidazole groups in the protonated form across an almost vanishing barrier. Proton affinity (PA) of histidine calculated by the G3(MP2) method is 233.2 and 232.4 kcal mol(-1) for protonation at the imidazole ring and at the amino group nitrogens, respectively, which is about 3-5 kcal mol(-1) lower than the reported experimental value.     

Topological study of intramolecular hydrogen bonding in beta-hydroxyethylperoxy radical and beta-hydroxyethoxy radical along its dissociation pathway

The presence of intramolecular hydrogen bond (IHB) in beta-hydroxyethylperoxy and beta-hydroxyethoxy radicals was investigated using the QTAIM topological study on UB3LYP/6-311++G(2d,2p) charge densities. Only one of the two conformers of beta-hydroxyethylperoxy radical which were previously considered to present IHB displays a bond critical point (BCP) associated to an IHB. Furthermore, the atomic energies and electron populations indicate no evidence of IHB in the second conformer. Nevertheless, very small differences in molecular energies were obtained using several one-step and multi-step methods (G3, G3B3) between both conformers. No BCP is found between the hydroxyl hydrogen and the oxygen in the most stable conformer or in the transition state for the dissociation path of beta-hydroxyethoxy radical. However, a BCP is formed in the last steps of this path, thereby yielding H-bonded products.     

二氟尼柳/水滑石插层组装结构、氢键及水合特性的分子动力学模拟

采用分子动力学方法模拟二氟尼柳插层水滑石(DIF/LDHs)的超分子结构, 研究复合材料主客体间形成的氢键以及水合膨胀特性.结果表明, 当水分子总数与DIF分子总数之比Nw≤3时, 层间距dc保持基本恒定, 约1.80 nm; 当Nw≥4时, 层间距逐渐增大, 且符合dc=1.2611Nw+13.63线性方程. 随着水分子个数增加, 水合能驻UH逐渐增大. 当Nw≤16时, 由于⊿UH<-41.84 kJ·mol-1, LDHs-DIF可以持续吸收水, 从而使材料层间距不断膨胀. 但当Nw≥24时, ⊿UH>-41.84 kJ·mol-1, 此时LDHs-DIF层间不能再进一步水合, 因此LDHs-DIF在水环境中膨胀具有一定的限度. 水滑石层间存在复杂的氢键网络. DIF/LDHs水合过程中, 水分子首先同步与层板和阴离子构成氢键; 当阴离子趋于饱和后, 水分子继续与层板形成氢键, 并逐步发生L-W型氢键取代L-A型氢键, 驱使阴离子向层间中央移动, 与层板发生隔离; 最后水分子在水滑石羟基表面形成有序结构化水层.