抗癌性喹唑啉类二氢叶酸还原酶抑制剂三维定量构效关系

在对一系列抗癌性7,8-二烃基-1,3-二氨基吡咯-[3,2-f]喹唑啉类二氢叶酸还原酶抑制剂的二维定量构效关系（2D—QSAa）研究基础上,应用比较分子场分析法对该类配合物进行了三维定量构效关系（3D—QSAR）研究．建立了具有良好的统计学性能及预报能力的3D．QSAR模型,非交叉验证相关系数为0．993,交叉验证相关系数为0．619,估算的标准误差0．208,统计方差比193．4．该模型表明立体场因素的影响比静电场因素大很多,此结果与我们已经报道的2D—QSAR模型结果相一致．然而,3D—QSAR模型提供了可视化的立体场、静电场因素对活性的影响．3D—QSAR研究对实验上提出的二氢叶酸还原酶与药物分子的疏水键合作用机理得到了进一步的理论解释．     

抗癌性吲哚喹唑啉衍生物的定量构效关系

用量子化学密度泛函理论(DFT)、分子力学(MM＋)及回归分析方法,对一系列抗癌性吲哚喹唑啉衍生物进行了定量构效关系(QSAR)的研究.通过回归分析,筛选了影响抗癌活性的主要因素,建立了定量构效关系方程.结果表明,化合物的最低未占据分子轨道(LUMO)与最高占据分子轨道(HOMO)之间的能量差(ΔεL-H)、分子的疏水性(lgP)以及环D上的总电荷(ΣQD)和环D上R1取代基的第一个原子的净电荷(QFR1)是影响化合物抗癌活性的主要因素.所得模型对化合物抗癌活性有较好的预测效果. 同时, 与ΔεL-H密切相关的LUMO轨道能量及共轭平面面积对药物的DNA-结合及其活性起着十分重要的作用,可通过选取具有较强的拉电子性质同时又能与本系列化合物的骨架形成更大共轭体系的取代基R1,设计抗癌活性较高的化合物.     

2-(4-取代-苯基)-3-异噻唑啉酮抗菌活性的定量构效关系研究

用密度泛函理论和逐步回归分析方法对22种新合成的2-(4-取代-苯基)-3-异噻唑啉酮化合物进行了结构与抗菌活性研究. 通过对平衡几何构型、前线轨道组成和Mulliken电荷分布的分析得出: S(1), N(2)原子是该类化合物的主要活性部位, 且5位氯取代的化合物抗菌活性较好. 通过回归分析, 筛选了影响抗菌活性的主要因素, 建立了定量构效关系方程. 结果表明, 硫原子的亲核前线电子密度、3D-Balaban指数、S(1)-N(2) Wiberg键级和Schultz分子价拓扑指数是影响异噻唑啉酮类化合物抗大肠杆菌活性的主要因素, 所得模型对化合物抗菌活性有较好的预报能力.     

DFT‐based QSAR study and molecular design of AHMA derivatives as potent anticancer agents

A quantitative structure–activity relationship (QSAR) of 3‐(9‐acridinylamino)‐5‐hydroxymethylaniline (AHMA) derivatives and their alkylcarbamates as potent anticancer agents has been studied using density functional theory (DFT), molecular mechanics (MM+), and statistical methods. In the best established QSAR equation, the energy (E_NL) of the next lowest unoccupied molecular orbital (NLUMO) and the net charges (Q_FR) of the first atom of the substituent R, as well as the steric parameter (MR₂) of subsituent R₂ are the main independent factors contributing to the anticancer activity of the compounds. A new scheme determining outliers by “leave‐one‐out” (LOO) cross‐validation coefficient (q) was suggested and successfully used. The fitting correlation coefficient (R²) and the “LOO” cross‐validation coefficient (q²) values for the training set of 25 compounds are 0.881 and 0.829, respectively. The predicted activities of 5 compounds in the test set using this QSAR model are in good agreement with their experimental values, indicating that this model has excellent predictive ability. Based on the established QSAR equation, 10 new compounds with rather high anticancer activity much greater than that of 34 compounds have been designed and await experimental verification. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

1-取代-2-氨基苯并咪唑衍生物的理论计算及QSAR研究

在B3LYP/6-31＋g(d)计算水平下, 通过运用密度泛函理论(DFT)量子化学方法, 对7种1-取代-2-氨基苯并咪唑衍生物的电子结构特征进行了研究, 获得了它们的电离能、亲合能、化学硬度、化学软度、化学势、电负性、亲电性、分配系数、折射率、极化率、分子体积和分子表面积等参数. 结果表明, 在2-氨基苯并咪唑衍生物的1位引入取代基团, 对化合物的电荷布居和结构性质都有较大的影响|关联了定量结构-活性关系(QSAR), 7种化合物的半数致死量LD₅₀与极化率等有较好的相关性. 这些结果可为1-取代-2-氨基苯并咪唑衍生物的氧化代谢实验研究和毒性机理揭示等提供重要信息.     

3D-QSAR Study on a Series of Indolo[1,2-b]quinazoline Derivatives with Anticancer Activity and their Molecular Design

Indolo[1,2-b]quinazoline derivatives have recently been reported as a type of compound with potential anticancer activity. On the basis of our the published two-dimensional quantitative structure-activity relationship (2D-QSAR) of these compounds, a further study on the three-dimensional quantitative structure-activity relationship (3D-QSAR) was carried out using the method of comparative molecular field analysis (CoMFA). A reasonable, receivable, and an effective 3D-QSAR model has been established, in which the correlation coefficient (r²) and cross-validation coefficient (q²) values are 0.986 and 0.695, respectively, the statistical squared deviation ratio (F) is 114.6, and the standard deviation (SD) is 0.084. The results suggest that the electrostatic effect of substituent R₁ and steric effect of substituent R₂ play a very important role in the improvement of the anticancer activity of these compounds. In this article, some significant conclusions, which are in good agreement with the conclusions obtained using 2D-QSAR, were drawn as follows. (1) The electrostatic effect in the substituent R₁ part plays a major role, and it is very important to make the first atom of R₁ carrying more positive charges in order to improve the anticancer activity of the compounds. (2) The steric effect plays a major role in the substituent R₂ part, and the volume of R₂ should be moderate. Based on the above conclusions, three new molecules of Indolo[1,2-b]quinazoline derivatives with higher anticancer activity have been theoretically designed and are waiting for support from experiment. The QSAR results can offer a theoretical reference for the pharmaceutical synthesis.     

Experimental and computational investigations on the high binding-selectivity of pyrimidine derivatives by a pillar[5]arene

The binding selectivity of an adenine-monofunctionalized pillar[5]arene (H) with a series of pyrimidine derivatives were investigated through ¹H NMR experiments and density functional theory (DFT) study. High binding-selectivity was demonstrated. Typically, H displayed very strong binding strength with 6-(2,4-dioxo-3, 4-dihydropyrimidin-1 (2H)-yl)hexanenitrile (G1) [K_a >10⁵ M^?1], up to about 3000-fold as compared with 1-hexylpyrimidine-2,4(1H, 3H)-dione (G5) [K_a = 31 M^?1]. The strong binding ability of H with G1 was due to the cooperative multiple hydrogen bond, dipole-dipole, C-H···π and π···π interactions. The high binding-selectivity was also verified by calculation results. The calculated interaction energy (ΔE_i) of G1?H was ?12.92 Kcal·mol^?1 while that of G5?H was ?2.85 Kcal·mol^?1.     

DFT Study for a Series of Less-Symmetrical Crown Ethers and Their Complexes with Alkali Metal Cations

A series of ring‐contracted (14‐crown‐5, 17‐crown‐6) and ring‐enlarged (16‐crown‐5, 17‐crown‐5, 19‐crown‐6, 20‐crown‐6) crown ethers and their complexes with alkali‐metal cations Na⁺ and K⁺ had been explored using density functional theory (DFT) at B3LYP/6‐31G* level in order to reveal the effects of the methylene‐chain length in a crown ether. The nucleophilicity of all crown ethers had been investigated by the Fukui functions. The quantum chemistry parameters, such as the energy gap (ΔE), the highest occupied molecular orbital energy (E_HOMO) and the lowest unoccupied molecular orbital energy (E_LUMO) for less‐symmetrical crown ethers and symmetrical frameworks (15‐crown‐5, 18‐crown‐6) had been calculated. In addition, the thermodynamic energies of complexation reactions had also been studied. The results of the DFT calculations show that the methylene‐chain length plays an important role in determining the structure characters of the crown ethers and also strongly influences the properties of the ethers. Some of the calculated results are in a good agreement with the experimental values.     

Theoretical study on conformational features and cation-binding properties of a diquinone calix[4]arene

Theoretical studies of a diquinone calix[4]arene and its interactions with the cations Li⁺, Na⁺, K⁺ and Ag⁺ have been performed. Conformational features and cation-binding properties were evaluated with the restricted hybrid Becke three-parameter exchange functional method using the 6-31G(d) basis set and its relativistic effective core potentials. To model the effect of medium, the polarisable continuum model was also used. Four typical conformations of the parent diquinone calix[4]arene were studied. The calculated results show that the most stable conformers are 1,3-alternate and partial cone in the gas phase and in CH₂Cl₂ solution, respectively. The optimised geometric structures were used to perform natural bond orbital analysis. The two main types of driving force metal–ligand and cation–π interactions are investigated. The calculated binding energy for cations (Li⁺, Na⁺, K⁺ and Ag⁺) is discussed. The calculated results indicate that cone complexes are the most stable.     

Studies on Quinazolines. Part III. Synthesis of Some Fused-Ring Cationic and Mesoionic Derivatives of the 1,3-Thiazolo[3,2-c]quinazoline Ring System

4-Allylthio-2-arylquinazolines 4a–c undergo cyclization by action of bromine to furnish 5-aryl-3-bromomethyl-2,3-dihydrothiazolo[3,2-c]quinazolin-4-ium bromide 5a–c. Compounds 5a–c undergo ring opening by action of water under acid catalysis to afford the corresponding dibromide derivatives 6a–c. Bromination of 3-allyl-2-aryl-4(3H)quinazolinethiones 7a–c leads to 5-aryl-2-bromomethyl-2,3-dihydrothiazolo[3,2-c]quinazolin-4-ium bromides 8a–c. However, 5-aryl-1,3-thiazolo[3,2-c]quinazolinium-3-olate 10a–c were prepared by the cyclodehydration of the corresponding thioglycolic acids 9a–c with Ac₂O.     

Dynamic Covalent Properties of a Novel Indolo[3,2-b]carbazole Diradical

This work describes the synthesis and properties of a dicyanomethylene-substituted indolo[3,2-b]carbazole diradical ICz-CN. This quinoidal system dimerises almost completely to (ICz-CN)₂, which contains two long C(sp³)−C(sp³) σ-bonds between the dicyanomethylene units. The minor open-shell ICz-CN component in the solid-state mixture was identified by EPR spectroscopy. Cyclic voltammetry and UV–visible spectroelectrochemical data, as well as comparison with reference monomer ICz-Br reveal that the nature of the one-electron oxidation of (ICz-CN)₂ at ambient temperature and ICz-CN at elevated temperature is very similar in all these compounds due to the prevailing localization of their HOMO on the ICz backbone. The peculiar cathodic behaviour reflects the co-existence of (ICz-CN)₂ and ICz-CN. The involvement of the dicyanomethylene groups stabilizes the close-lying LUMO and LUMO+1 of (ICz-CN)₂ and especially ICz-CN compared to ICz-Br, resulting in a distinctive cathodic response at low overpotentials. Differently from neutral ICz-CN, its radical anion and dianion are remarkably stable under ambient conditions. The UV/Vis(–NIR) electronic transitions in parent (ICz-CN)₂ and ICz-CN and their different redox forms have been assigned convincingly with the aid of TD-DFT calculations. The σ-bond in neutral (ICz-CN)₂ is cleaved in solution and in the solid-state upon soft external stimuli (temperature, pressure), showing a strong chromism from light yellow to blue–green. Notably, in the solid state, the monomeric diradical species is predominantly formed under high hydrostatic pressure (>1 GPa).     

功能离子液体[NH2p-bim]BF4吸收CO2的密度泛函研究

用密度泛函理论对功能化离子液体[NH2p-bim]BF4吸收CO2的作用机制进行了理论研究。在RB3LYP/6-311 G**的计算水平对离子液体[NH2p-bim]BF4的结构以及CO2与该离子液体反应可能生成的产物进行了全优化,获得了优化结构的振动频率和热力学数据。计算结果表明,离子液体[NH2p-bim]BF4吸收CO2主要是通过离子液体的阳离子[NH2p-bim] 自偶解离的[NHp-bim]与CO2分子结合生成[O2C-NHp-bim],其结合能为238-260 kJ/mol。     

杯[4]二吡咯与卤素阴离子相互作用的密度泛函理论研究

采用密度泛函B3LYP方法, 在6\|31G*/LANL2DZ水平上对杯[4]二吡咯主体分子及其与卤素阴离子形成的复合物进行研究.  结果表明, 杯[4]二吡咯可与卤素阴离子相互作用形成具有较高对称性的复合物, 其8个吡咯NH基团上的H原子均可以和卤素阴离子形成氢键; 杯[4]二吡咯与卤素阴离子结合能力大小的顺序为F- >Cl- >Br- > I-.  振动光谱、 电荷分布以及前线轨道计算结果表明, 杯[4]二吡咯与卤素阴离子相互作用的本质为氢键, 经BSSE校正的结合能与电荷转移程度、 N-H键键长和N-H伸缩振动频率变化呈线性关系.      

Molecular design of a “molecular syringe” mimic for metal cations using a 1,3‐alternate calix[4]arene cavity

The chemically switchable actions well imitate the function of a “molecular syringe,” has been studied in theory using the 1,3‐alternate calix [4]arene bearing a nitrogen‐containing crown cap at one side and a bis(ethoxyethoxy) group at another side by the π‐basic calixtube as a pipette and the crown ring as a rubber cap. The model is characterized by geometry optimization using density functional theory (DFT) at B3LYP/6‐31G level. The obtained optimized structures are used to perform natural bond orbital (NBO) and frequency analysis. The electron‐donating heteroatoms: O and N offer lone pair electrons to the contacting RY* (1‐center Rydberg) or LP* (1‐center valence antibond lone pair) orbitals of K⁺, Ag⁺. The results indicate that when the nitrogen atom in the crown ring is protonated, K⁺ and Ag⁺ will be pushed out to the bis(ethoxyethoxy) side through a π‐basic calixtube. When the nitrogen·H⁺ in the crown ring is deprotonated, K⁺ and Ag⁺ are sucked back to the crown‐capped side again. In the course of the coordination, both the intermolecular electrostatic interactions and the cation‐π interactions between the metal ion and π‐orbitals of the two pairs facing inverted benzene rings play a significant role. It is believed that this prototype of a “molecular syringe” is a novel molecular architecture for the action of metal cations. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Theoretical Study on a 42 T Model of Beta Zeolite by Density Functional Theory Simulation

Density functional theory was applied to study the structure of Beta zeolite. A model cluster containing 41Si atoms, 1 Al atom, 70 O atoms and 29 H atoms was constructed. The model structures were optimized using the Becke＇s three-parameter hybrid method with the Lee-Yang-Parr correlation functional （B3LYP） and the 6-31G basis set applying the Gaussian03 program package. The NMR parameters were calculated to validate the rationality of the model. It was found that in the optimization models, all O-H bond lengths were in range of 0.984-0.985A^°, among which the model with O-H bond length of 0.98478A^° was more stable than the others. The ^1H and ^27Al chemical shifts of the most stable model were 4.03434 and 55.74 ppm, which were pretty consistent with Larry＇ s experimental data of 4.1 and 54 ppm. The relationship between other structure parameters and total relative electric energy has also been found. All the results exhibit that the 42 T （the total number of Si and Al atoms is 42） model has common properties of the standard of zeolite Beta.     

多溴吩噻嗪系列化合物的稳定性和热力学性质的密度泛函理论研究

用Gaussian 03程序, 在B3LYP/6-31G*水平上全优化计算了吩噻嗪和135个多溴吩噻嗪系列化合物(PBPTHs)在298.15 K和101.3 kPa状态时的热力学参数. 设计等键反应, 计算了PBPTHs系列化合物的标准生成热( )和标准生成自由能( ). 同时研究了这些参数与溴原子的取代位置及取代数目(NPBS)之间的关系. 结果表明: 多溴吩噻嗪的热能校正值(Eth)、恒容热容( )、标准熵( )、标准焓( )以及标准自由能( )与NPBS之间有很强的相关性(r2≥0.998). 在相关方程中, 溴原子的取代个数对多溴代吩噻嗪热力学数值的大小有很大影响. 根据 的相对大小, 从理论上求得异构体的相对稳定性顺序.     

密度泛函理论处理激发态与多重态结构研究进展

用密度泛函理论处理分子激发态和电子多重态结构是密度泛函理论发展的重要方面, 研究工作一直很活跃。本文对近年来在这一领域取得的重要进展作简要的综述。     

氮混杂金属卟啉吸收光谱性质的理论研究

采用密度泛函理论在B3LYP/6-31+G(d)水平上研究了4种金属Mg, Ni, Cu, Zn配位的自由卟啉(FBP)及氮混杂卟啉(NECP)的几何结构及分子轨道能级. 采用含时密度泛函理论(TD-DFT)方法计算了金属与2种卟啉配位后在气体条件下的电子吸收光谱, 包括激发能、 吸收波长、 跃迁组成和振子强度.计算结果表明, 与金属配位的FBP(M-FBP)具有D_4h对称性, 分子轨道能级HOMO/HOMO-1和LUMO/LUMO+1因能级相近发生简并, HOMO-LUMO轨道能级差大约3.0 eV, 在Soret带出现较强吸收峰.由于C/N原子位置的改变, 非对称性结构的M-NECP前线轨道组成发生改变, 轨道能级差(HOMO-LUMO)减小至2.6 eV左右, 且能级发生分裂, Soret带出现多个电子吸收谱峰, Q带也出现吸收峰. 本文研究了水、 氯仿和苯3种不同极性溶剂对M-FBP和M-NECP的分子轨道及电子吸收光谱的影响, 结果表明, 随溶剂极性减弱金属配合物的电子吸收光谱发生红移, 并且吸收峰强度增强.     

抗癌性吲哚喹唑啉衍生物3D-QSAR研究及其分子设计

吲哚喹唑啉衍生物是近年来发现的一类具有良好抗癌活性的化合物. 作者在最近报道的二维定量构效关系(2D-QSAR)的基础上, 采用比较分子力场方法(CoMFA)进一步对该系列化合物进行三维定量构效关系(3D-QSAR)研究, 建立了3D-QSAR的CoMFA模型, 其非交叉验证相关系数r2=0.986, 标准偏差SD=0.084, 统计方差比F=114.6, 交叉验证相关系数q2=0.695, 表明该模型合理、可信, 并具有良好的预测能力. 研究结果表明: (1) 取代基R1的部位上静电效应起主要作用, 并且确保取代基R1的第一个原子具有较大的净正电荷, 对提高化合物的抗癌活性十分重要. 这与2D-QSAR研究结果相一致. (2) 取代基R2的部位上立体效应起主要作用, R2的体积大小要适中. 应用这些规律进行了分子设计, 在理论上获得了一些具有较高抗癌活性的新的吲哚喹唑啉衍生物, 并期待实验证实. 该QSAR的研究结果可为实验工作者合成新药提供理论参考.