理论线性溶解能参数在苯砜基羧酸酯类化合物急性毒性QSAR研究中的应用

采用密度泛函理论方法,在B3LYP/6-31G(d)理论水平下,计算了45种苯砜基羧酸酯化合物的量子化学参数.经多元线性回归分析,得到描述此类化合物对发光菌急性毒性的模型:-lgEC50=3.02 6.24EHOMO-0.091μ-0.006P 1.22q(1)-6.67q(10),其中R=0.92,r2adj=0.82,F=42.0,q2=0.79.通过对模型进行分析,得到如下结论:苯环和酯基取代基的电负性越大,分子体积越小,毒性越大.该研究为探讨此类化合物急性毒性的机理奠定了理论基础.     

含有咪唑生色团系列树型有机分子二阶非线性光学性质的密度泛函理论研究

应用密度泛函理论(DFT)B3LYP/6-31G*方法计算研究了系列树型含有咪唑生色团的有机分子的结构和非线性光学性质.计算结果表明:该系列分子具有A- -D- -A(A:受体,D:给体)结构,分子基态的偶极矩、极化率、二阶NLO系数( )随共轭链的增长及吸电子基的增强而增大;同时,前线轨道能级差值越小此类分子的二阶极化率总有效值( )越大.计算的吸收光谱显示此系列树型分子在低能区域247.79nm-419.87nm都有一个最强吸收,并且均是最高占据轨道与最低空轨道之间的跃迁.     

苯衍生物二阶非线性光学性质的量子化学研究

本文采用引入外场微扰的CNDO/S—CI方法,计算了一系列苯衍生物的分子二阶非线性光学系数,探讨了取代基的电子性质、取代位置及取代基数目对分子二阶非线性光学系数的影响.结果表明;取代基的供电能力越强,分子的二阶非线性光学系数越大.D,A对位二取代苯和1—D,2.4—A,A三取代苯具有较高的非线性光学效应,苯衍生物中取代基数目为3时分子具有最佳的非线性光学效应.文中还对上述结论给出了初步的理论解释.     

水杨醛缩胺类Schiff碱光致变色性质

用时间相关光谱法研究了水杨醛缩胺类Schiff碱光致变色的性质,发现分子共轭体系大小和分子骨架刚性等因素对光致变色性质存在规律性的影响.共轭体系越大,发生光致变色的速率越小;分子骨架的刚性越大,化合物越难发生光致变色;反之,则光致变色较容易进行.     

醇(酚)羟基的硅烷化保护与去保护研究进展

醇(酚)羟基的硅烷化保护是一类重要的有机合成手段,目的在于使羟基稳定化,消除或减轻其引起的副反应。保护基中硅原子连接的基团空间位阻越小,该保护基的反应活性越大,生成的相应硅醚的稳定性则越差,在弱酸或弱碱的条件下即可脱除;硅原子连接的基团越大,该保护基的反应活性则越小,硅醚化反应越难发生,需要借助催化剂才能进行。本文尽可能全面地论述了有机硅烷保护基的类型,如三甲基硅烷、三乙基硅烷、叔丁基二甲基硅烷、三异丙基硅烷、苯基取代硅烷和桥型硅烷等,并讨论了其在不同环境下的活性及稳定性。     

13位取代苦参碱衍生物抗癌活性的定量构效关系研究

采用密度泛函理论B3LYP方法计算了13个13位取代苦参碱衍生物的电子结构,研究了化合物结构与抑制人肝癌细胞HepG2抗癌活性的定量构效关系(QSAR).结果表明:(1)13位取代的苦参碱类衍生物的最低空轨道能ELUMO越低,最低空轨道与最高占据轨道的能隙ΔE越小,化合物抗癌活性越高;(2)分子的能量Etotal、面积S以及体积V越大,其极化度P越大,活性越大;(3)分子的油水分配系数logP越大,活性越大,即分子的疏水性增大活性增强.综合得到了显著性较好的QSAR方程:-lgIC50=97.008-11.759ΔE+818.602QC2-2.132×10-4Etotal,可用于预测该类衍生物抑制人肝癌细胞HepG2的活性并进行分子设计.     

氟代十二顶点碳硼烷结构和非线性光学性质的理论研究

采用密度泛函理论(DFT)方法,计算并分析了四氟取代十二顶点碳硼烷及其衍生物的结构和非线性光学(NLO)性质.结果表明,改变四氟碳硼烷取代基的共轭性或给吸电子能力,会使分子中碳硼笼原子间距离发生改变.碳硼烷取代基的给吸电子能力越强,其偶极矩越大.分子极化率随取代基共轭性和体积的增加而增大.引入强吸电子基或增加四氟碳硼烷取代基的共轭性,可使其二阶NLO响应明显增强.通过分析分子的电子光谱和对应的分子轨道组成可知,第一超极化率最大的分子4a’发生碳硼笼到并苯取代基的电荷转移.     

药物水溶解度的QSAR研究和VolSurf参数

药物的水溶解度与其吸收密切相关。本文利用一种新的计算方法,VolSurf,预测药物的水溶解度并测定有利于药物水溶解度的主要分子特征。被测化合物包括26个结构不同的药物,通过偏最小二乘分析法,对药物水溶解度实验值与分子特征进行相关,得到较好的模型(r2=0.90,q2=0.77)。将化合物分为训练集和预测集进行相关分析,结果表明以18个化合物所建立的训练集模型对其余8个化合物有较好的预测能力,预测的标准偏差(SDEP)为0.59。参数分析表明分子与水相互作用的3个局部能量最小值越小,且它们之间的距离越大,对其水溶解度越有利;亲水性占主导因素的分子有高的水溶解度;分子的疏水性越强,在水中的溶解性越弱;大分子的溶解度较小分子溶解度低。     

8-羟基喹啉银(铂)金属配合物电子光谱与非线性光学性质

运用密度泛函理论B3LYP方法对8-羟基喹啉(银、铂)(AgQ、PtQ₂)金属配合物及其衍生物的非线性光学性质进行理论计算研究. 结果表明: 引入取代基使铂配合物的最强吸收波长产生较大红移. 最低能量跃迁吸收来自最高占据分子轨道(HOMO)到最低空分子轨道(LUMO)的d→π^*和π→π^*跃迁, 属于金属配体电荷转移(MLCT)与配体配体电荷转移(LLCT). 金属银和铂掺杂8-羟基喹啉使其三阶非线性光学系数γ值明显增大, 并且在配合物上引入―Ph, ―PhOCH₃, ―PhF₂, ―PhF₅基团将进一步增大γ值. 引入基团的供电子性越强, γ值增大的幅度越大, 引入基团的吸电子性越强, γ值增大的幅度越小.     

C(4)取代紫杉醇类似物的定量构效关系研究

用量子化学B3LYP／6-31G方法计算了23个C（4）取代紫杉醇类似物的结构,用遗传算法（GFA）对能量、电性、拓扑及热力学等类型的278个结构描述符进行筛选,并回归建立其抑制人体结肠癌细胞HCT-116活性的定量构效关系（QSAR）．QSAR方程含分子体积Vm、分子分支度指数CHI-O、分子中带正电荷原子的溶剂可积面积与其所带电荷之积的加和值Jurs-PPSA-3以及分子表面积S4个结构描述符．方程的拟合相关系数的平方R^0及交叉验证系数Q^2分别为0．956和0．913,所得QSAR具有可信的预报能力．由优化后的几何构型知,C（4）取代基、C（13）侧链和2-OBz三基团共同形成疏水腔,C（4）取代基的改变影响C（13）侧链的电子结构．C（13）连接的18号O原子的负电荷越大、3’位连接的NHBz基团的极性越小活性越高;C（4）取代基若为吸电子基对活性不利;适当增大分子体积、表面积和疏水性,保持一定的分支度对活性有利．     

CoMFA of the acute toxicity of phenylsulfonyl carboxylates to Vibrio fischeri

From the Comparative Molecular Field Analysis (CoMFA) method, the paper describes a three-dimentional quantitative structure-activity relationship (3D-QSAR) model for the acute toxicity logEC50 (15min-EC50 in micromol l(-1)) of 56 phenylsulfonyl carboxylates on Vibrio fischeri. The achievement of a high leave-one-out (LOO) cross-validated correlation coefficient q2 of 0.790 with four optimum components indicates the significance of the correlation of the steric and electrostatic fields with the biological activities. The key features in the CoMFA contour maps are critical to trace the important properties and gain insight into the toxic mechanism of the tested phenylsulfonyl carboxylates.     

磺酰脲类化合物除草活性的QSAR研究

采用密度泛函理论方法, 在B3LYP/6-31G(d)水平下, 计算了23种磺酰类化合物的分子极化率及分子骨架中各原子的Milliken电荷. 提出了一种新的QSAR建模方法, 并据此对其中18种化合物进行多元线性回归分析, 建立了除草活性的预测模型(R=0.96, R2=0.92, r2adj=0.88, F=26.26, q2=0.71, p<0.01, SE=0.36), 对剩余五种化合物进行预测, 结果吻合. 该模型从化合物的亲水性、分子几何特征的角度对如何提高磺酰脲类化合物的除草活性进行了分析, 并对提高化合物除草活性的方法做出预测: 提高苯环和嘧啶环取代基的亲水性, 增加N13周围的电子云密度, 为苯环接入较小的取代基团, 在嘧啶环上接入较大取代基团都可提高化合物的除草活性. 预测结果与3D-QSAR方法的预测结果一致.     

Fluxional processes in diamagnetic and paramagnetic allyl dicarbonyl and 2-methylallyl dicarbonyl molybdenum histidinato complexes as revealed by spectroscopic data and density functional calculations

This work describes a detailed study on the structure and dynamics of pseudooctahedral low-valent complexes of the type [Mo(His-N(epsilon)-R)(eta-2-R'-allyl)(CO)(2)] (His=N(delta),N,O-L-histidinate; R=H, R'=H (1); R=C(2)H(4)CO(2)Me, R'=H (2); R=H, R'=Me (3); R=C(2)H(4)CO(2)Me, R'=Me (4)). These diamagnetic 18-electron complexes were comprehensively characterized spectroscopically and by X-ray crystallography. In the solid state, the (substituted) allyl ligand is in an endo position in all compounds, but it is trans to the His-N(delta) atom in 1 and 2, whereas it is trans to the carboxylate O atom for the 2-Me-allyl compounds 3 and 4. In solution, both isomers are present in a solvent-dependent equilibrium. The third isomer (allyl trans to His-NH(2)) is not spectroscopically observed in solution. This is in agreement with the results from density functional (DFT) computations (BPW 91 functional) for 1 and 3, which predict a considerably higher energy (+6.3 and +5.9 kJ mol(-1), respectively) for this isomer. A likely path for isomerization is calculated, which is consistent with the activation energy determined by variable temperature NMR measurements. At least for 3, the preferred path involves several intermediates and a rotation of the 2-Me-allyl ligand. For the paramagnetic 17-electron congeners, DFT predicts the exo isomer of 3(+) with the 2-Me-allyl ligand trans to the carboxylate O atom to be by far the most stable isomer. For 1(+), an endo-exo equilibrium between the isomers with the allyl ligand trans to the carboxylate O atom is suggested. These suggestions are confirmed by EPR spectroscopy on the electrochemically generated species, which show signals for one- (4) and two- (2) metal-containing compounds. The appearance of the EPR spectra may be rationalized by inspection of the SOMOs from DFT calculations of the species in question. The notion of a metal-centered oxidation is also substantiated by IR spectroelectrochemistry and by UV/Vis spectra of the 17-electron complexes. Upon depleting the metal of electron density, the stretching vibrations of the carbonyl ligands shift more than 100 cm(-1) to higher wavenumbers, and the carbonyl vibration of the metal-coordinated carboxylate shifts by about 50 cm(-1). A color change from yellow to green upon oxidation is observed visually and quantified by the appearance of a new band at 622 nm (2(+)) and 546 nm (4(+)), respectively.     

弹性竿模型下超螺旋DNA分子的构象研究

采用弹性竿模型 (Elasticrodmodel) ,用MonteCarlo方法对DNA分子的构象进行研究 .通过计算发现 ,DNA分子的能量是由弯曲势能EB 和扭转势能ET 两部分组成 ,通常EB 比ET 大一至两个数量级 .同时给出了均方回转半径与链长之间的关系为〈R2g〉 =1 1 69 5 -3 5×n +0 0 2 5×n2 ,它体现了DNA分子结构的特点 .验证了公式Lk=Wr+Tw ,得出Lk与Wr比较接近的结论 ,考虑DNA分子的构型 ,意味着DNA分子容易被弯曲而不易被扭转 ,但随着连接系数的增加 ,DNA被扭转的几率也在增加 .这为分析DNA分子的结构特征提供了一种新方法     

Molecular Modeling and Design of Arylthioindole Derivatives as Tubulin Inhibitors

Three-dimensional quantitative structure activity relationship （3D-QSAR） and docking studies of a series of arylthioindole derivatives as tubulin inhibitors against human breast cancer cell line MCF-7 have been carried out. An optimal 3D-QSAR model from the comparative molecular field analysis （CoMFA） for training set with significant statistical quality （R2=0.898） and predictive ability （q2=0.654） was established. The same model was further applied to predict pIC50 values of the compounds in test set, and the resulting predictive correlation coefficient R2（pred） reaches 0.816, further showing that this CoMFA model has high predictive ability. Moreover, the appropriate binding orientations and conformations of these compounds interacting with tubulin are located by docking study, and it is very interesting to find the consistency between the CoMFA field distribution and the 3D topology structure of active site of tubulin. Based on CoMFA along with docking results, some important factors improving the activities of these compounds were discussed in detail and were summarized as follows： the substituents R3-R5 （on the phenyl ring） with higher electronegativity, the substituent R6 with higher eleetropositivity and bigger bulk, the substituent R7 with smaller bulk, and so on. In addition, five new compounds with higher activities have been designed. Such results can offer useful theoretical references for experimental works.     

A study of trialkyltin β-aryl-β-triphenylgermyl propionates

Twenty new compounds of the form Ph₃GeCHArCH₂COOSnR₃ (R = n-Bu, cyclohexyl; Ar = substituted phenyl) have been synthesized. Their structures were characterized by IR and ¹¹⁹Sn and ¹H NMR spectroscopy. The compounds are five-coordinated carboxylate bridged polymers when R = n– Bu; when R = cyclohexyl (Cy) they are four-coordinate. ¹¹⁹Sn NMR measurements of chemical shift for the two series of compounds have shown that there is a good linear relationship for the chemical shift of ¹¹⁹Sn NMR between the tributyltin and tricyclohexyltin propionates, viz. δ¹¹⁹Sn(Bu₃Sn) = 1.0474 δ ¹¹⁹Sn(Cy₃Sn) + 95.8076, n = 5, r = 0.993. The structure of one compound was determined by X-ray diffraction. It exists as a monomeric four-coordinated species in a distorted tetrahedronal geometry.