单嘧磺隆除草剂的晶体构象-活性构象转换的密度泛函理论研究

应用密度泛函理论在B3LYP/6-31G(d)水平上对新型除草剂单嘧磺隆绕脲桥部分两个C—N键的内旋转势能面进行计算, 然后对势能面上的驻点进行构型优化和过渡态搜索, 得到单嘧磺隆4种稳定构象和构象转换过程所涉及的8个过渡态结构. 研究结果表明, 单嘧磺隆晶体构象-活性构象转换过程中涉及4种稳定构象和8条转换途径, 脲桥部分NH基团与嘧啶环上N原子所形成的分子内氢键对于构象的稳定性及转换过程起着十分重要的作用. 应用极化连续介质溶剂模型(PCM)在B3LYP/6-31＋＋G(d, p)水平下进行溶剂化效应计算, 结果表明单嘧磺隆从晶体构象转换成活性构象主要是在水相中进行的.     

杯[4]吡咯在不同溶液中的分子动力学模拟研究

运用分子动力学(MD)模拟方法对杯[4]吡咯与不同溶剂之间的相互作用能、杯[4]吡咯在不同溶液中的构象变化以及杯[4]吡咯与溶剂分子之间的氢键相互作用进行了计算研究.模拟发现,杯[4]吡咯与不同溶剂间的相互作用能受溶剂分子偶极矩和杯[4]吡咯-溶剂分子间氢键相互作用影响.杯[4]吡咯在不同溶液中的构象发生翻转的主导因素是杯[4]吡咯与溶剂分子间形成氢键相互作用,溶剂分子的偶极矩不是杯[4]吡咯发生构象转化的主要因素.     

多肽抑制剂抑制淀粉质多肽42构象转换的分子动力学模拟和结合自由能计算

应用分子动力学模拟和结合自由能计算方法研究了多肽抑制剂KLVFF、VVIA和LPFFD抑制淀粉质多肽42 (Aβ42)构象转换的分子机理. 结果表明, 三种多肽抑制剂均能够有效抑制Aβ42的二级结构由α-螺旋向β-折叠的构象转换. 另外, 多肽抑制剂降低了Aβ42分子内的疏水相互作用, 减少了多肽分子内远距离的接触, 有效抑制了Aβ42的疏水塌缩, 从而起到稳定其初始构象的作用. 这些抑制剂与Aβ42之间的疏水和静电相互作用(包括氢键)均有利于它们抑制Aβ42的构象转换. 此外, 抑制剂中的带电氨基酸残基可以增强其和Aβ42之间的静电相互作用(包括氢键), 并降低抑制剂之间的聚集, 从而大大增强对Aβ42构象转换的抑制能力. 但脯氨酸的引入会破坏多肽的线性结构, 从而大大降低其与Aβ42 之间的作用力. 上述分子模拟的结果揭示了多肽抑制剂KLVFF、VVIA和LPFFD抑制Aβ42构象转换的分子机理, 对于进一步合理设计Aβ的高效短肽抑制剂具有非常重要的理论指导意义.     

槲皮素与尿嘧啶和胸腺嘧啶氢键的结合位点

深入理解药物分子和核酸碱基间的相互作用机制对合理设计研发新型高效药物有重要意义. 本文运用密度泛函理论B3LYP方法对核酸碱基尿嘧啶和胸腺嘧啶与药物分子槲皮素间的氢键相互作用位点进行了研究. 使用B3LYP/6-31G(d)方法优化得到了30个稳定的氢键复合物结构, 使用B3LYP/6-311++G(3df,2p)方法计算了这些复合物的结合能. 研究结果表明, 槲皮素可以使用5个不同的结合位点与尿嘧啶或胸腺嘧啶形成氢键复合物, 尿嘧啶或胸腺嘧啶可以使用3个不同的结合位点与槲皮素形成氢键复合物. 当槲皮素的结合位点固定时, 槲皮素与尿嘧啶的位点u1或胸腺嘧啶的位点t1形成的氢键作用最强, 与位点u2或位点t2形成的氢键强度最弱; 当尿嘧啶或胸腺嘧啶的作用位点固定时, 二者与槲皮素的位点qu1 形成的氢键作用最强, 与位点qu5 作用强度次之, 与位点qu3的作用强度最弱. 分子中原子(AIM)和自然键轨道(NBO)分析计算结果表明, 轨道作用在氢键中起重要作用.     

3,4,5-三羟基苯甲酸二聚体氢键结构性质

分子间相互作用在生物和材料等科学中发挥着关键作用,研究分子间相互作用的本质意义重大。氢键是分子间相互作用的一种主要形式,在确定分子构象和晶体结构以及生物分子尤其是核酸和蛋白质的结构功能中起着重要作用[1-3]。苯甲酸衍生物广泛存在于生物大分子内,与生物活性离子通过氢键作用等改变生物活性分子的活性功能,研究苯甲酸衍生物分子间氢键相互作用对于了解生物体内的化学现象具有重要意义。研究表明菱角的抗肿瘤作用明显,实验上已经从菱角中成功提取了活性单体化合物:3,4,5-三羟基苯甲酸二聚体[4],理论研究标题化合物的氢键结构与氢…     

单嘧磺隆除草剂分子印迹聚合物的识别特性研究

为了研究印迹聚合物作为固相提取剂在环境分析中的应用,合成了以除草剂单嘧磺隆为印迹分子的印迹聚合物,通过平衡吸附和液相色谱的方法分析了印迹聚合物的识别特性,并研究了不同色谱条件对单嘧磺隆保留因子的影响研究结果表明,印迹聚合物对印迹分子具有很好的亲和性和特定的选择性,可作为固相提取剂,在单嘧磺隆的残留分析中,对土壤提取样品进行有效的富集与净化     

紫杉醚与αβ微管蛋白的分子对接

采用Insight II/Affinity对紫杉醚与αβ微管蛋白进行分子对接, 共得到10个对接构象. 应用密度泛函B3LYP/6-31G 方法计算对接口袋构象的结合能, 筛选出结合能达-190.53 kJ·mol-1的最优对接构象5. 通过构象分析建立紫杉醚与受体结合的作用模型, 结果表明, 在活性口袋的底部紫杉醚与受体间的作用主要是疏水作用, 而在活性口袋的顶部两者间主要是氢键作用. 氢键作用位置可分为A和B两个作用区, 其中A区有3个氢键, 由C13侧链分别与受体的ASP26和ARG369作用形成; B区也有3个氢键, 是由紫杉醚母环上的极性基团分别与受体的THR276、ARG278和GLN282作用产生的. 紫杉醚与αβ微管蛋白间形成的6个氢键可以有效地将紫杉醚固定在活性口袋中.     

核酸碱基与甘氨酸二肽间氢键作用的最佳位点

优化得到了碱基腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶与甘氨酸二肽分子形成的28个氢键复合物的稳定结构并计算了结合能, 探讨了五种碱基与甘氨酸二肽分子间氢键作用的最佳位点. 本文研究发现: 每种碱基均可以通过不同位点与二肽分子形成氢键复合物, 腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶分别最倾向使用A3、T1、U1、G3及C1位点与甘氨酸二肽分子形成氢键复合物; 碱基分子某位点的质子化反应焓变越负所形成的氢键复合物越稳定, 去质子化反应焓变越小所形成的氢键复合物越稳定; 由氢键复合物的结合能计算得到的稳定性次序与由碱基分子质子化和去质子化反应焓变推得的稳定性次序一致.     

4,5,6-三取代嘧啶苯磺酰脲类化合物的生物活性、分子对接与3D-QSAR关系研究

用柔性分子对接方法(FlexX)将15个4,5,6-三取代嘧啶苯磺酰脲化合物以及3个不含5-位取代嘧啶苯磺酰脲化合物(分别为4,6-双取代嘧啶和4-取代嘧啶)和乙酰羟酸合成酶(AHAS)活性口袋进行了对接, 对接程序预测的抑制剂和酶之间的相互作用能与抑制活性之间有一定的相关性, 相关系数为0.660. 然后采用比较分子相似性指数分析(CoMSIA)对27个新型4,5,6-三取代嘧啶苯磺酰脲类化合物的除草活性进行三维定量构效关系(3D-QSAR)研究. 建立了三维定量构效关系CoMSIA模型, 立体场、静电场和氢键的贡献分别为47.3%, 32.8%, 19.9%. 交叉验证系数q2值为0.520. 根据CoMSIA模型的立体场、静电场、氢键给体场三维等值线图不仅直观地解释了结构与活性的关系, 并且与用FlexX预测的结合模式相一致, 证明了我们预测的结合模式是可靠的, 为进一步设计高活性的标题化合物提供较好的理论指导.     

DFT方法研究基于脲及硫脲衍生物的阴离子受体

采用密度泛函理论(DFT)模拟了两个基于脲和硫脲衍生物的受体分子对卤素阴离子的识别过程.结构优化表明基于脲衍生物的受体分子1最稳定构象为"反反"构象,分子内部形成稳定的C~α-H…O=C分子内氢键;而基于硫脲衍生物的受体分子2,不能形成分子内氢键,最稳定构象为"反顺"构象.受体1、2与卤素阴离子F~-、Cl~-可形成稳定的双氢键复合物,在此过程中,受体2经历了由"反顺"构象到"反反"构象的异构化过程.结构和能量分析表明,1、2受体分子与F~-离子间的氢键强度远大于其与Cl~-离子间的氢键;另一方面,受体2与阴离子间的氢键明显强于受体1,这是由于硫脲基N-H键具有更强的酸性.此外,对受体分子、氢键复合物及去质子化产物的吸收光谱计算结果表明,受体与F~-离子作用可产生明显的吸收光谱红移,而与Cl~-离子的作用对光谱影响较小.     

Conformation and location of membrane-bound salinomycin-sodium complex deduced from NMR in isotropic bicelles

An ionophore antibiotic salinomycin was studied in a membrane environment consisting of isotropic bicelles, a better model for biological membranes than micelles, and its conformation and topological orientation in bicelles was determined. 2D NMR measurements and restrained conformational search revealed that salinomycin-sodium salt in bicelles adopts an open conformation in which the orientation of the E-ring is significantly different from that in crystal and solution structures. This conformational alteration breaks an intramolecular hydrogen bond between 28-OH and 1-O, dislocates the ether oxygen of the E-ring from a coordinated position to the sodium ion observed in the crystal, and consequently weakens the complexation between salinomycin and the sodium ion. Paramagnetic relaxation experiments using doxyl-phospholipids reveal that salinomycin is embedded shallowly in bicelles, with both terminals being closer to the water interface and the olefin portion facing the bicelle interior. Measurements of intermolecular NOEs between salinomycin and phospholipids further supported this orientation. Weaker complexation with sodium ion and positional preference in the membrane polar region may facilitate the catch-and-release of metal ions at the polar/nonpolar interface of bilayers. On the basis of these findings, a model for salinomycin-assisted transport of metal ions across biological membranes is proposed.     

An NMR investigation of the importance of intramolecular hydrogen bonding in determining the conformational equilibrium of ethylene glycol in solution

Although conformational analysis by NMR of ethylene glycol indicates generally strong preferences for the gauche conformation in solvents ranging from water to chloroform, the bulk of the NMR evidence indicates that intramolecular hydrogen bonding between the hydroxyl groups is unlikely to be a significant factor in determining that preference, except possibly in fairly non-polar solvents. The 'gauche effect' is clearly very important, especially in aqueous solution.     

Solvent effects on the polymerization of multi-armed vinyl biphenyl derivatives containing pendant oligo(oxyethylene)cyclotriphosphazenes

Polymerization of vinyl biphenyl derivatives containing a pendant oligo(oxyethylene)cyclotriphosphazene (VBMEP, ? (OCH₂CH₂)nOCH₃, n = 1; VBDEP, n = 2; VBTEP, n = 3) was carried out in various solvents. The conversions of these monomers increased with increasing β values, solvent hydrogen bond acceptor abilities, indicating that the hydrogen bond formation is the most important factor in the polymerization. ¹³CNMR study showed that the reactivity of the monomer is influenced by the hydrogen bond interaction. In ethanol, the kinetic orders of monomer and initiator concentrations for the polymerization of VBDEP were different from those in 1,2-dichloroethane (DCE), which suggest the predominant occurrence of primary radical termination. The intrinsic viscosity of poly(VBDEP) with M?_n = 22 000 in DCE was two times higher than that in ethanol, and plots of intrinsic viscosity versus conversion of VBDEP gave a straight line. The results suggest that the polymer chains in ethanol are in a coiled conformation, whereas in DCE they are in a relatively extended structure, and that the propagation is affected by the conformational change. These behaviors originated from the hydrogen bond formation between polymers and solvents are discussed. The copolymerization of styrene with multiarmed monomers and the properties of polycascade polymers obtained are also described. © 1993 John Wiley & Sons, Inc.     

酶联免疫分析法探测Cry1Ab蛋白在不同介质中的构象变化

采用酶联免疫分析（ELISA）方法,根据构象变化后的蛋白与抗体结合能力下降从而导致ELISA测定值降低的原理,探测了转cry1Ab基因水稻表达的Cry1Ab蛋白在不同溶液介质中的热致构象变化行为,以及不同有机溶剂及溶液pH值对该蛋白构象变化的影响程度。实验表明,Cry1Ab蛋白在不同条件下的构象变化程度可以灵敏地通过ELISA方法检测。在不同的介质中,CrylAb蛋白的热致构象变化程度不同。在Na2SO4介质中,该蛋白具有较高的热稳定性;SDS的存在,可以促进该蛋白的构象变化。常温下,25％（V／V）的有机溶剂乙腈、异丙醇、甲醇、乙醇均能使该蛋白的构象发生转变,其中以乙腈最为显著。醇类溶剂对Cry1Ab蛋白的构象影响程度随疏水性增大而增大;溶液pH值也对该蛋白的构象变化产生影响。pH在8-10之间,该蛋白构象能保持稳定;酸或过碱性的溶液均能使蛋白构象偏离原始状态,从而引起ELISA测定值的降低。另外,腐殖酸能在一定时间内保持Cry1Ab蛋白构象的稳定性。     

Combination of molecular dynamics method and 3D-RISM theory for conformational sampling of large flexible molecules in solution

We have developed an algorithm for sampling the conformational space of large flexible molecules in solution, which combines the molecular dynamics (MD) method and the three-dimensional reference interaction site model (3D-RISM) theory. The solvent-induced force acting on solute atoms was evaluated as the gradient of the solvation free energy with respect to the solute-atom coordinates. To enhance the computation speed, we have applied a multiple timestep algorithm based on the RESPA (Reversible System Propagator Algorithm) to the combined MD/3D-RISM method. By virtue of the algorithm, one can choose a longer timestep for renewing the solvent-induced force compared with that of the conformational update. To illustrate the present MD/3D-RISM simulation, we applied the method to a model of acetylacetone in aqueous solution. The multiple timestep algorithm succeeded in enhancing the computation speed by 3.4 times for this model case. Acetylacetone possesses an intramolecular hydrogen-bonding capability between the hydroxyl group and the carbonyl oxygen atom, and the molecule is significantly stabilized due to this hydrogen bond, especially in gas phase. The intramolecular hydrogen bond was kept intact during almost entire course of the MD simulation in gas phase, while in the aqueous solutions the bond is disrupted in a significant number of conformations. This result qualitatively agrees with the behavior on a free energy barrier lying upon the process for rotating a torsional degree of freedom of the hydroxyl group, where it is significantly reduced in aqueous solution by a cancellation between the electrostatic interaction and the solvation free energy.     

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.     

Conformation, Inversion Barrier, and Solvent-Induced Conformational Shift in Exo- and Endo/Exo-Calix[4]arenes

Calixarenes 4a and 4b having hydroxyl groups in endo and exo positions and the ethanediyl-bridged exo-calixarene 5a were synthesized by a stepwise strategy. Single-crystal X-ray structures were obtained for 4a and for the exo-calixarene 3d, showing the molecules to exist in the 1,2-alternate conformation which is also found for 4a,b in solution. The inversion barriers of 4a and 4b (10.3 and 10.8 kcal mol(-1)) are similar to that determined for the endo-dihydroxycalixarene 12, indicating that the additional intramolecular hydrogen bond between the exo OH groups does not decrease the flexibility of the molecule. In CDCl(3) solution exo-calixarene 5a adopts a 1,2-alternate conformation with the methyl group at the bridge located in an axial position, while in DMSO-d(6) the conformation adopted is the partial cone. Similar solvent-induced conformational shifts were found for the exo-calixarenes 3b and 3d. MM3 calculations predict that the cone form is the lowest energy conformation of 4 and the exo-calixarenes 3 and 5. The calculations suggest that the conformational preferences of the methyl group at the bridge for either the axial or equatorial positions are in large part determined by the repulsive steric interactions with the hydroxyl groups. The inversion barrier of 4b is satisfactorily reproduced by calculations, which indicate that the rotation of the exo rings is less energetically demanding than the rotation of the endo rings.     

Conformational preferences and cis-trans isomerization of azaproline residue

The conformational study of N-acetyl-N'-methylamide of azaproline (Ac-azPro-NHMe, the azPro dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the effects of the replacement of the backbone CHalpha group by the nitrogen atom on the conformational preferences and prolyl cis-trans isomerization in the gas phase and in solution (chloroform and water). The incorporation of the Nalpha atom into the prolyl ring results in the different puckering, backbone population, and barriers to prolyl cis-trans isomerization from those of Ac-Pro-NHMe (the Pro dipeptide). In particular, the azPro dipeptide has a dominant backbone conformation D (beta2) with the cis peptide bond preceding the azPro residue in both the gas phase and solution. This may be ascribed to the favorable electrostatic interaction or intramolecular hydrogen bond between the prolyl nitrogen and the amide hydrogen following the azPro residue and to the absence of the unfavorable interactions between electron lone pairs of the acetyl carbonyl oxygen and the prolyl Nalpha. This calculated higher population of the cis peptide bond is consistent with the results from X-ray and NMR experiments. As the solvent polarity increases, the conformations B and B* with the trans peptide bond become more populated and the cis population decreases more, which is opposite to the results for the Pro dipeptide. The conformation B lies between conformations D and A (alpha) and conformation B* is a mirror image of the conformation B on the phi-psi map. The barriers to prolyl cis-trans isomerization for the azPro dipeptide increase with the increase of solvent polarity, and the cis-trans isomerization proceeds through only the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond for the azPro dipeptide in the gas phase and in solution, as seen for the Pro dipeptide. The pertinent distance d(N...H-NNHMe) and the pyramidality of imide nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure and the lower rotational barriers for the azPro dipeptide than those for the Pro dipeptide in the gas phase and in solution.