HEPT类逆转录酶抑制剂的三维定量构效关系

利用比较分子力场分析(CoMFA)方法对32个HEPT类HIV-1逆转录酶抑制剂(RTIs)的三维定量构效关系(3D-QSAR)进行了分析,建立了HIV-1逆转录酶抑制剂的3种3D-QSAR模型,发现影响其生物活性的主要因素为立体场因素,这与HIV-1RT的非底物结合部位(NNBS)的疏水性环境相吻合.进一步分析表明,适当长度的1-位侧链对保持化合物的抗病毒活性致关重要;增大5-位取代基的体积可增强生物活性;在1-位苄氧甲基的对位引入大体积基团有利于提高活性.同时考察立体场、静电场与生物活性的关系,表明,CoMFA模型为最佳预测模型,其交叉验证系数R_CV²=0.870,传统相关系数R²=0.986,标准偏差SE=0.146,F=294.546.用此模型预测了检验组3个HEPT类化合物的-lgEC₅₀,R_pred²=0.850,表明模型具有很好的预测能力,可为HEPT类HIV-1逆转录酶抑制剂的结构优化提供理论指导.     

芳香噻嗪类衍生物作为选择性醛糖还原酶抑制剂的分子对接和基于受体的三维定量构效关系研究

芳香噻嗪类衍生物被证明是一类选择性较好的高活性醛糖还原酶抑制剂(ARIs).本文对44个芳香噻嗪类化合物进行了分子对接(docking)和三维定量构效关系(3D-QSAR)研究,并探索了此类化合物与醛糖还原酶(ALr²)的作用机理.醛糖还原酶与醛还原酶(ALR1)活性位点的叠加结果显示, ALr²中残基Leu 300和Cys298的存在是化合物1m具有高选择性的原因.分别建立了比较分子场分析方法(CoMFA, q² = 0.649, r² =0.934; q²:交叉验证相关系数, r²:非交叉验证相关系数)和比较分子相似性指数分析方法(CoMSIA, q² = 0.746, r² = 0.971)模型,并对影响此类化合物生物活性的结构进行了鉴定.结果显示,两个模型均具有较高预测能力,并通过测试集中的7个化合物进行了验证,其中CoMFA模型和CoMSIA模型的预测相关系数(r_Pred²)分别为0.748和0.828. 3D-QSAR模型中的三维等值线图表明,在化合物1m的苄基环上C3和C4位置以及苯并噻嗪母核上C5和C7位置进行改进可能对生物活性的提高有利,此预测与我们前期报道的苯并噻嗪母核C7位改进结果一致.本文所建3D-QSAR模型能够在理性设计具有更高生物活性的新型ARIs中发挥重要作用.     

基于分子对接和QSAR方法预测B-Raf II型抑制剂活性

B-Raf激酶在促分裂素原活化蛋白激酶(MAPK)信号转导通路中起着重要作用,已被确定为癌症治疗非常有吸引力的靶标.新型高效B-Raf抑制剂的开发成为癌症治疗的一个热门研究领域.本文以结构多样的B-Raf II型抑制剂为研究对象,联合应用分子对接和定量构效关系(QSAR)模型研究其定量构效关系去探讨抑制活性的起源.两个主题作为研究重点:生物活性构象和描述符.首先对分子对接方法(Glide、Gold、LigandFit、Cdocker和Libdock)进行准确性评价,后将研究的对象分子对接到B-Raf活性位点并获得生物活性构象.基于准确的对接结果,计算得到16个打分评价函数和21个能量描述符,以此构建定量构效关系模型. QSAR结果表明模型具有高度精确的拟合和强的预测能力(模型M1: r² = 0.852, r_(CV)² = 0.790, r_pre² = 0.864;模型M2: r² = 0.738, r_(CV)² = 0.812, r_pre² = 0.8605).同时探讨了对抑制活性有重要影响的描述符,结果表明打分评价函数(G_Score, -ECD, Dock_Score, PMF)与能量描述符(S(hb_ext), DE(int), Emodel)对抑制活性影响非常大.通过虚拟筛选和QSAR模型理论预测,一些新的具有潜在抑制活性的化合物作为B-Raf II型抑制剂被获得.上述信息对于进一步设计新颖高效的B-Raf II型抑制剂提供了有用的指导.     

取代嘧啶衍生物抗前列腺癌活性的3D-QSAR研究

运用比较分子力场分析(CoMFA)方法,建立18种取代嘧啶衍生物抗前列腺癌活性(pM)的三维定量构效关系。训练集中15个化合物用于建立预测模型,测试集13个化合物(含10号模板分子和新设计的9个分子)作为模型验证。建立的CoMFA模型的交叉验证系数(R_ev²)、非交叉验证系数(R²)分别为0.344、0.935,说明所建模型具有较强的鲁棒性和良好的预测能力。该模型中立体场、静电场贡献率依次为71.6%、28.6%。影响取代嘧啶衍生物抗前列腺癌活性的主要因素是取代基的疏水作用和空间位阻,其次是取代基的库仑力、氢键及配位作用。基于此研究结果,设计了9个新化合物,其抗前列腺癌活性有待医学实验验证。     

氟喹诺酮C-3均三唑衍生物抗增殖活性的QSAR研究

基于拓扑化学理论,用原子类型的电性拓扑状态指数(E_A)描述了48个氟喹诺酮C-3均三唑衍生物分子的化学微环境。基于E_A和最佳变量子集回归,建立上述化合物对人白血病细胞(HL60)的抗增殖活性(pH)的定量构效关系(QSAR)模型。其最优四元(E₈、E₁₃、E₁₉、E₄₂)模型的判定系数(R²)和逐一剔除法交叉验证系数(R_cv²)分别为0.853和0.813。经R²、R_cv²、F_IT、A_IC、F、V_IF等统计指标检验,所建模型具有良好的稳定性、相关性和预测能力。结果显示影响氟喹诺酮C-3均三唑衍生物对人白血病细胞(HL60)的抗增殖活性的主要因素是与靶标形成氢键,以及配位和疏水作用。     

脒类KARI酶抑制剂的分子对接和3D-QSAR研究

对设计合成的20个单脒类化合物的水稻KARI酶体外抑制活性和体内除草活性进行了分子对接和三维定量构效关系研究. 前者采用AutoDock3.0方法, 研究发现化合物活性变化趋势与分子对接计算结果基本一致, 通过分析化合物9与KARI酶活性氨基酸残基结合模式发现, 残基Glu319, Asp315, Glu496, Gly253, Met254, Cys517等对氢键和静电相互作用以及疏水作用都有重要贡献; 后者研究采用比较分子力场分析(CoMFA)方法, 结果表明立体场和静电场对活性的贡献分别为67.8%和32.2%, 交叉验证系数r_cv²＝0.774, 非交叉验证r²＝0.999, F＝1593.134, 标准偏差s＝0.036, 所建立的3D-QSAR模型对化合物除草活性具有较好的预测能力. 两种方法研究结果为进一步设计合成更高活性的KARI酶抑制剂提供了指导.     

电拓扑状态预测多氯二苯并噻吩及噻吩砜化合物的气相色谱保留指数

以原子类型电拓扑状态指数(ETSI)有效表征了135 个多氯二苯并噻吩(PCDT)和135 个多氯二苯并噻吩砜(PCDTO₂)的分子结构, 应用基于预测的变量选择与模型化(VSMP)方法建立PCDT 和PCDTO₂ 化合物在DB-5 气相色谱柱上的气相色谱保留指数(RI)与分子结构(ETSI)的定量相关模型, 模型的相关系数r² 分别为0.9939 和0.9729, LOO 交叉验证相关系数 q² 分别为0.9921 和0.9692. 为验证模型稳定性和预测能力, 应用17 个PCDT 和PCDTO₂ 训练集样本构建的QSRR 模型的r² 分别为0.9959 和0.9783, LOO 交叉验证相关系数 q² 分别为0.9921 和0.9740, 说明模型具有良好的稳定性. 以此模型预测外部8 个检验集及110 个预测集的RI 值, 8 个检验集样本的结果表明训练集模型具有良好预测能力.     

3C-like蛋白酶抑制剂的构效关系、分子对接和分子动力学

3C-like蛋白酶是中东呼吸综合征冠状病毒（MERS-CoV）等其它冠状病毒的繁殖过程中极为重要的蛋白酶。它已成为人类在抗冠状病毒领域中的研究热点。本文基于计算生物学方法对与MERS-CoV同属的蝙蝠冠状病毒HKU4（HKU4-CoV）的43个肽类3C-like蛋白酶抑制剂分子,建立三维定量构效关系（3D-QSAR）模型。在基于配体叠合的基础上,发现比较分子相似性指数分析法（CoMSIA）中的四个场组合（位阻场、静电场、氢键供体场与氢键受体场）为最优的模型（Q²=0.522,R_ncv²=0.996,R_pre²=0.904;Q²：交叉验证相关系数,R_ncv²：非交叉验证相关系数,R_pre²：验证集分子的预测值相关系数）,并借助该模型通过分子对接（docking）与分子动力学（MD）方法阐明了配受体结合作用。实验结果表明：（1）基于最优的CoMSIA模型基础上的三维等势图形象地说明了分子基团的位阻作用、静电作用、氢键供体与氢键受体作用对分子生物活性的影响;（2）分子对接研究结果显示了疏水性以及结晶水、氨基酸His166和Glu169在配体和受体结合过程中产生重要作用;（3）分子动力学模拟进一步验证了分子对接模型的可靠性,并发现了两个新的关键氨基酸Ser24与Gln192,它们与配体产生了两个较强的氢键。此外,根据这些结果,一些新的具有潜在抑制活性的肽类化合物作为3C-like蛋白酶抑制剂被获得。以上结果能够帮助深入了解3C-like蛋白酶与肽类抑制剂的作用机理,并且能够为今后的抗MERS-CoV药物设计提供有价值的参考。     

磷改性失活TS-1高效催化C5=裂解制备C2=/C3=反应的研究

发展C₅⁼催化裂解高效制乙烯和丙烯(C₂⁼/C₃⁼)技术具有重要的研究意义及工业应用价值,其核心在于高效催化剂的开发.本文发展了以失活钛硅分子筛TS-1(De-TS-1)为催化剂的C₅⁼催化裂解过程,研究发现,P改性后的De-TS-1显示出优异的催化反应性能,C₂⁼+C₃⁼选择性与单程寿命分别达77.1%和213 h,其综合催化性能优于经典的ZSM-5催化剂.进一步研究表明,De-TS-1的Br?nsted性质符合C₅⁼裂解高效裂解的基本要求,P改性可进一步降低De-TS-1中Br?nsted酸的酸密度和酸强度,从而促进了目标产物C₂⁼/C₃⁼生成的主反应,抑制了C₅     

电性拓扑态指数预测PCDFs光解半衰期

以电性拓扑态指数(E-state indices)为参数, 应用偏最小二乘法(PLS)对48种多氯代二苯并呋喃(PCDFs)在云杉(Picea abies)针叶和飞灰表面的光解半衰期(t_1/2)进行模拟分析, 分别获得5变量和4变量定量结构-性质相关(QSPR)模型. 模型交叉验证相关系数(Q_cum²)分别为0.758和0.701, 标准偏差(SD)分别为0.047和0.049, 模型具有较高的预测能力和可靠性. 随着氯原子取代数目的增多, PCDFs与活性氢供体接触的“位阻效应”增大, PCDFs分子中苯环的静电荷密度随之降低, 以致高氯代PCDFs较难发生光解作用.     

拉曼光谱法研究Na+, Mg2+//SO42-, Cl-, H2O四元体系中SO42-离子缔合结构特征及其变化规律

复杂水盐体系存在稳态和介稳固液相平衡以及复杂的成盐规律。为了解固液平衡状态下液相的结构特征,本文采用拉曼光谱技术并结合高斯-洛伦兹去卷积分峰拟合程序对Na⁺, Mg²⁺//SO₄^2-, Cl^-, H₂O四元体系及其二元和三元子体系中ν₁-SO₄^2-的离子缔合结构特征进行了分析。研究结果表明:SO₄^2-在Na₂SO₄-H₂O体系存在自由态SO₄^2-和SO₄^2-离子簇两种结构,在MgSO₄-H₂O, MgSO₄-MgCl₂-H₂O及Na⁺, Mg²⁺//SO₄^2-, Cl^-, H₂O等含镁体系中,还有Mg²⁺-H₂O-SO₄^2-和Mg²⁺-OSO₃^2-两种缔合结构。在二元和三元体系中ν₁-SO₄^2-的离子缔合结构以自由态SO₄^2-为主,随着SO₄^2-离子总浓度的变化,上述四种结构所占比例会发生规律性变化。Na⁺, Mg²⁺//SO₄^2-, Cl^-, H₂O四元体系在NaCl减少及等温蒸发过程中,自由态SO₄^2-结构比例逐步降低, Mg²⁺和SO₄^2-相结合形成Mg²⁺-H₂O-SO₄^2-或Mg²⁺-OSO₃^2-结构的机会增多,在复盐区还会形成SO₄^2-离子簇结构。由此判断溶液结构的适应性变化是导致介稳现象的重要原因。进一步的相关分析表明:SO₄^2-的浓度和耶涅克指数J与ν₁-SO₄^2-峰的峰强度和峰面积存在正相关关系, Mg²⁺浓度是影响ν₁-SO₄^2-峰中四种缔合结构的比例发生变化的主要因素。     

Hydrogen peroxide photoproduction by the semicarbazide—tris(2,2′-bipyridine)ruthenium(II)—oxygen system

A light-driven system consisting of tris(2,2′-bipyridine)ruthenium(II) (Ru(bpy)₃²⁺) as the photosensitizer, semicarbazide as the electron donor and molecular oxygen as the electron acceptor has been employed for hydrogen peroxide production. The efficiency of this photosystem markedly depends on pH: while the peroxide yield is almost negligible at acid, neutral or slightly alkaline pH, it reaches significant values at high hydroxide concentrations, the initial rate of H₂O₂ formation drastically increasing from pH 12 to pH 14. In 1 M NaOH solutions containing Ru(bpy)₃²⁺ and semicarbazide at optimum concentrations, the number of catalytic cycles (or turnover number) undergone by the ruthenium complex over the complete course of the photochemical reaction is as high as 1.1 × 10⁴.

Spectrofluorometric and laser flash photolysis techniques were used to study the primary photochemical reactions involving the excited state of the ruthenium complex as well as the photochemically generated species Ru(bpy)₃³⁺ and Ru(bpy)₃⁺. It is proposed that at pH 14 a sequence of reactions leading to O₂ photoreduction by electrons from semicarbazide takes place, with the concomitant formation of H₂O₂; the excited state of Ru(bpy)₃²⁺ appears to react via oxidative quenching by oxygen rather than via reductive quenching by semicarbazide. At neutral pH, in contrast, there is no H₂O₂ formation owing to the fact that semicarbazide is unable to reduce (Ru(bpy)₃³⁺ to Ru(bpy)₃²⁺, although the photoexcited ruthenium complex is quenched equally by oxygen.     

Ion chromatography of organic-rich natural waters from peatlands: I. Cl, NO2, Br, NO3, HPO4, SO4 and oxalate

Organic-rich natural waters from peat bogs in continental (Switzerland) and maritime (Shetland Islands, Scotland) areas were analysed for Cl⁻, NO₂⁻, Br⁻, NO₃⁻, HPO₄²⁻, SO₄²⁻ and oxalate using ion chromatography. These anions can be determined simultaneously in the surface and pore water samples from the continental bogs using a 250-μl injection loop. Using this loop, the detection limits were ca. 5 ng/g for the monovalent anions and SO₄²⁻ and 10 ng/g for HPO₄²⁻ and oxalate. An organics-removal cartridge (Dionex OnGuard P) was used to remove humic materials. These cartridges did not significantly affect the measured concentrations of anions in blind standards. Analyses of deionized water treated with these cartridges are not significantly different from those for untreated deionized water. For the maritime bogs, the relatively high concentrations of Cl⁻ (more than 100μ/g in many samples) and SO₄²⁻ (up to 50 μg/g) require two separate determinations for complete analyses. A 10-μl injection loop was used to determine Cl⁻, Br⁻ and SO₄²⁻. A 250-μl injection loop was used to measure NO₂⁻, NO₃⁻, HPO ₄²⁻ and oxalate. In each instance a Dionex OnGuard P cartridge was used to remove humic materials. In addition, a chloride-removal cartridge (Dionex OnGuard AG) was used to remove Cl⁻ when the larger injection loop was used. This cartridge has no significant effect on the measurement of HPO_4-²⁻ at concentrations of 20 ng/g. In each of the bog water chromatograms there were usually a number of unknown peaks. These are probably due mainly to organic anions.     

The vibration-rotation spectrum and anharmonic potential of H3

Some qualitative effects of anharmonicity on the spectra of H₃⁺ and D₃⁺ between low vibrational levels are described. Using large-basis vibration-rotation calculations with a Morse-based discrete variable representation for the vibrations and a symmetric-top basis for the rotations, new spectra of H₃⁺ and D₃⁺ have been assigned. This procedure was assisted by adjusting eight coefficients for H₃⁺ and six coefficients for D₃⁺ in the Meyer-Botschwina-Burton ab initio potential function, and eventually 621 new and old lines of H₃⁺ to levels up to 3ν₂ and 529 new and old lines of D₃⁺ to levels up to 2ν₂ have been fitted with standard deviations of 0.118 and 0.059 cm⁻¹, respectively. An attempt is made to compare five different potential energy functions for the H₃⁺ system, two ab initio and three adjusted to fit spectra of H₃⁺ or D₃⁺, by expanding them by the same procedure in the same variables. For extension of the present work to higher vibrational levels, more accurate boundary behaviour at linear configurations will be required, and some aspects of the use of hyperspherical coordinates are discussed.     

Testing the uncertainty principle using (n,γ) photon beams and pulsed neutron beams

Nuclear resonance scattering of a photon beam based on the (n,γ) reaction was used for measuring the mean square zero-point linear momenta p_a² and p_c² of the atoms of a metallic Zn single crystal. A similar measurement was carried out using a pulsed neutron source. These were combined with corresponding values of x_a² and x_c² obtained using the Mössbauer effect. In both the cases the measured values of p_a² and p_c² were found to conform to the requirements of the Uncertainty principle for a real solid to within 3%.     

Flexible regulation of C3=/C2= ratio in methanol-to-hydrocarbons by delicate control of acidity of ZSM-5 catalyst

A very wide range of the C₃⁼/C₂⁼ ratio from 0.72 to 7.56 with high C₂⁼+C₃⁼ selectivity of around 66% in the methanol-to-hydrocarbons process can be realized over ZSM-5 catalyst in a fixed-bed reactor. We firstly conduct a single factor experiment of acidity, demonstrating that the acidity control of MTH catalyst is crucial to adjusting light olefins selectivity. Weak Brønsted acid sites favor to high C₃⁼ selectivity (59.0%) due to the suppression of the conversion reactions from the alkene-based to arene-based cycle, while Lewis acid sites conduce to high C₂⁼ selectivity (39.6%) due to the promotion of the conversion reactions for the aromatics formation and steric constraints of Lewis acid sites making the aromatics crack more efficiently to C₂⁼.     

VO1-4+阳离子与n-CmH2m+2 (m = 3, 5, 7)烷烃的反应性研究：氧含量和碳链长度的影响

Vanadium oxides are one of the most important heterogeneous catalysts that are widely used to oxidize hydrocarbon molecules into value-added chemicals. In order to reveal the mechanisms and the nature of active sites, numerous experimental and theoretical studies have been reported on the reactivities of gas-phase vanadium oxide clusters toward small molecules. However, there has been very limited research on the chemical reactivity changes associated with the oxygen contents of vanadium oxides and the carbon chain lengths of alkanes. In this work, the reactions of vanadium oxide ions VO₁₋₄⁺ with alkanes (n-C_mH_2m+2, m = 3, 5, 7) were systematically investigated by time-of-flight mass spectrometry and the reactions of VO₁₋₃⁺ with pentane were further studied by density functional theory calculations. Experimental results show that in the reactions of VO⁺, VO₃⁺, and VO₄⁺ with n-C₅H₁₂, in addition to the major adsorption processes, the activation of the C―H and C―C bonds of n-C₅H₁₂ was observed. The activation of both the bonds was observed experimentally during the reaction of VO₂⁺ with n-C₅H₁₂ with large branching ratios. Among the vanadium oxide cations studied, VO₂⁺ shows the strongest oxidizability and the generation of lighter alkanes and alkenes dominates the reactions; VO⁺ is more reactive than VO₃⁺. VO₄⁺ pocesses only one η²-O₂ unit. Due to the weak bond between VO₂⁺ and η²-O₂, the η²-O₂ unit is released in VO₄⁺/n-C₅H₁₂ system leading to the formation of VO₂⁺; thus VO₄⁺ cations reflect some reactivity of VO₂⁺. Although the oxidation states in the vanadium oxide clusters increase from +Ⅲ in VO⁺ to +Ⅴ in VO₂⁺ and +Ⅳ in VO₃⁺, the reactivity does not gradually increase. Moreover, the reactivity of the mononuclear vanadium oxide cations also does not exhibit a gradually increasing trend with the increase in oxygen content. Based on the observed reactivity trend, the adsorption channels gradually become weak as the carbon chain lengths increase; meanwhile, the dehydrogenation and C―C bond activation channels gradually become obvious and some oxygen transfer products appear. Therefore, much lighter fragments of alkanes/alkenes will be obtained if linear alkanes with more carbon atoms were reacted with VO₁₋₄⁺. The theoretical results are generally consistent with those obtained from the experiments. The various reaction channels and versatile reactivity of the mononuclear vanadium oxide cations investigated in this study not only offer new insights into gas-phase reactions but also shed light on the processes occurring on the surfaces of the corresponding condensed-phase catalysts.     

SO42-/ZrO2和SO42-/Al2O3-ZrO2催化剂上的正戊烷反应

研究了添加Al对SO₄^2-/ZrO₂超强酸样品的晶化、比表面、硫含量、超强酸性和正戊烷反应性能的影响,考察了活化温度、反应温度、Al含量和载Pt对催化剂活性和选择性的影响.SO₄^2-/Al₂O₃ZrO₂催化剂的酸强度与SO₄^2-/ZrO₂基本相当,但超强酸位比后者多,未载Pt时正戊烷反应活性和稳定性明显高于后者.负载Pt后,正戊烷异构化选择性和稳定性大大提高,但Pt/SO₄^2-/Al₂O₃ZrO₂催化剂的反应活性与Pt/SO₄^2-/ZrO₂相近,Al的促进作用不明显.     

大黄素的电化学氧化还原机理研究

本文利用循环伏安法(CV)、红外光谱循环伏吸法(CVA)和导数循环伏吸法(DCVA)研究大黄素(Q)在乙腈溶剂中的电子转移机理.Q的还原过程中阴离子自由基Q^?-会结合中性分子Q生成二聚物Q₂^?-.Q₂^?-在更负的电位下进一步还原为Q₂^2-.当扫描范围为-0.2 ~ -2.0 V时,经过一个循环伏安过程,在扫描结束物质并没有回到反应物Q,而是Q₂^2-. Q₂^2-会继续发生电化学反应,经历两步一电子过程,分别生成Q₂^3-. 和Q₂^4-.,对应CV图中峰C₃和C₄.当扫描范围扩大至1.0~-2.0 V时,在更正的电位下,观察到两个新的氧化峰A₁和A₂,该范围内的三圈扫描结果表明,在扫描结束物质重新氧化回到Q.当扫描范围缩小至0.3 ~ -1.4 V,A2峰随着扫描圈数的增加而增大,与A2峰对应的氧化产物Q₂^?-在溶液中不断积累.A1峰对应于Q₂^?-氧化回到Q.