应用分子图形学、分子力学、量子化学及静电势研究农药分子结构与性能的关系（Ⅴ）──磺酰脲类分子中的化学键及电子结构

用半经验分子轨道方法ＭＮＤＯ及ＣＮＤＯ／２对农药分子结构进行计算，确认磺酰脲超高效除草剂分子中存在３个原子平面，分别形成离域π键，硫原子的空ｄ轨道参与了共轭体系的形成，静电势计算结果表明在分子两侧存在着２个明显的负电势区域，该结构与分子活性密切相关。     

应用分子图形学,分子力学,量子化学及静电势研究农药分…

用半经验分子轨道方法ＭＮＤＯ及ＣＮＤＯ／２对农药分子结构进行计算，确认磺酰脲类超高效除草剂分子中存在３个原子平面，分别形成离域π键，硫原子的空ｄ轨道参与了共轭体系的形成，静电势计算结果表明在分子两侧存在着２个明显的负电势区域，该结构与分子活性密切相关。     

分子间的非共价相互作用：σ-穴键和π-穴键

第四到第七主族元素以及惰性气体元素可以作为电子受体位点形成分子间吸引性的非共价键合作用。这些电子受体位点,即缺电子密度的位点,大多情况下具有正的分子表面静电势,可以将其分为2类。沿着σ共价键键轴延伸方向原子的外围中心具有正的分子表面静电势的区域称之为σ-穴。而垂直于分子σ-骨架平面具有正的分子表面静电势区域称之为π-穴。σ-穴和π-穴与富电性位点之间吸引性的相互作用分别称之为σ-穴键和π-穴键。σ-穴键倾向于180°,而π-穴键倾向于90°。按照元素周期表族的名称,σ-穴键又分为碳素键、氮素键、氧素键、卤键、惰素键等。可见流行的卤键只是σ-穴键的一个子集。π-穴键分类显得复杂些,可简单分为单原子和多原子π-穴键。     

从头计算法研究沙蚕毒和杀虫环分子的几何构型与电子结构

用Hartree-Fock/6-31G^*从头算确定了沙蚕毒和杀虫环分子的几何构型,在全局优化中发现杀虫环分子的椅式和船式两种稳定构象,在二级Moller-plesset微扰理论MP2/6-31^*水平下,椅式较船式稳定27.06kJ/mol.用MP2/6-31G^*波函数计算电子相关校正的分子静电势,以此为基础讨论生物活性与静电势的关系。发现对此二分子,Mulliken布居分析获得的原子净电荷存在问题,本文用Breneman提出的从静电势导出原子净电荷的CHELPG方法计算了原子净电荷。     

预测亲核反应位点方法的比较

预测分子不同位点发生亲核取代反应的活性具有重要的理论和实际意义.目前已提出了许多基于反应物自身电子结构特点的预测方法.本文将碳基化合物、芳香族化合物和吡啶衍生物这3类分子作为测试体系,对14种预测方法的可靠性进行了详细的比较分析.结果表明,体现局部电子软度的方法能很好地预测反应位点,如简缩双描述符方法,但表现静电效应的预测方法,如原子电荷分析和静电势分析,整体表现很差.对于本文中所用的分子体系来讲,简缩双描述符和Hirshfeld电荷分析方法对分子反应位点预测最准确.     

环多肽晶体的浮动电荷极化力场模拟

利用原子键电负性均衡结合分子力场方法(ABEEM/MM)对五种环多肽晶体进行了研究. 与传统力场相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理地体现了分子中的电荷分布. 相对其他极化力场模型, 具有计算量较小的特点. 该模型下计算得到的环多肽分子单元相对实验测得的结构的原子位置、氢键长度和二面角的均方根偏差分别为0.009 nm、0.013 nm和5.16°, 能够很好地重复实验结果. 总体上, 其结果优于或相当于其他力场模型, 适用于对实际蛋白质体系的模拟和研究.     

富勒烯结构B_(80)分子静电势的理论研究

在混合密度泛函B3LYP理论下,用6-31G*基函数对富勒烯结构B80分子的3个异构体(1个具有Ih对称性,2个具有Th对称性)构型进行优化和分子静电势计算.结果表明:3个异构体球内全部为正电势,球外五元环中心所对应的区域都为负电势,B80Ih,Th(A)和Th(B)球外静电势的最大负值分别对应于20个六元环中心的B原子,五元环中心和12个六元环中心的B原子周围,它们组成了化学反应中最可能的活性点.     

分子亲脂-亲水性的量子化学描述(Ⅰ)——分子的亲脂和亲水表面

给出基于分子结构的“启发式”亲脂-亲水势HMLP(Heuristic molecular lipophilicity-hydrophilicity potential)的理论分析和有说服力的算例.用量子化学计算其分子表面的静电势V(r)的分布,通过与周围原子表面静电势的比较,构造表达分子静电势的极性和大小的函数L(r).亲脂势L(r)保留了静电势V(r)描述分子静电作用的能力,并把应用范围扩展到疏水性的描述.HMLP不使用原子的经验参数,但在L(r)的构造中使用了经验的函数形式.经参数化和指标化后,HMLP有望用于蛋白质结构与功能的研究和药物分子配体与生物大分子受体结合自由能的估算.     

基于密度泛函理论研究碱性水溶液中乙内酰脲及其衍生物的结构及性质

乙内酰脲及其衍生物对多种金属离子具有良好的络合效果,是一种有望取代氰化物的绿色电镀络合剂.本文以密度泛函理论(DFT)方法对乙内酰脲及其衍生物的反应活性进行了研究.通过分子极性指数和极性表面积所占总面积百分比研究多种乙内酰脲衍生物在水中的溶解性,理论计算结果与文献报道的实验结果基本一致.通过原子电荷、表面静电势、原子对...     

利用浮动电荷力场对N-甲基乙酰胺水溶液的分子动力学模拟

利用原子键电负性均衡结合分子力场方法(ABEEM/MM)对N-甲基乙酰胺(NMA)分子的水溶液体系进行了分子动力学模拟. 与经典的力场模型相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理体现了分子中的电荷分布. 相对其它极化力场模型, 该模型具有计算量较小的特点. 在该模型下对NMA纯溶液和其水溶液体系进行了分子动力学模拟, 得到的径向分布函数、汽化热和偶极矩等物理量与实验值和其它极化力场方法符合很好, 合理描述了溶质与溶剂之间的静电极化和分子内的电荷转移.     

Exploring the conformation,charge density distribution and the electrostatic properties of galanthamine molecule in the active site of AChE using DFT and AIM theory

Quantum chemical calculations and charge density analysis were carried out to understand the geometry, charge density distribution and the electrostatic properties of isolated galanthamine molecule (form I) and for the same lifted out from the active site (form II) of AChE. The optimized geometry of isolated galanthamine was obtained from a hybrid density functional theory (B3LYP/6‐311G**) calculation. A docking analysis on galanthamine with AChE was performed, and the lowest docked energy structure was selected from the active site of AChE for the further study. A single point energy quantum chemical calculation (B3LYP/6‐311G**) was carried out for the lowest energy structure, which was lifted from the galanthamine–AChE complex from molecular docking analysis. The structural comparison between (I) and (II) helps to understand the conformational modification of the galanthamine molecule in the active site. When the molecule present in the active site, the molecular geometry is seen to be significantly altered, specifically, large changes were observed in the outer core of the molecule while the inner core geometry is intact. The bond topological and electrostatic properties of (I) and (II) were calculated. The dipole moment of the galanthamine molecule also increases from 2.09 to 2.67 D in the process. A large negative electrostatic potential region is found at the vicinity of oxygen and nitrogen atoms of the molecule, which predominantly involve strong hydrophobic and electrostatic interactions with the amino acid residues TRP84, PHE330, GLY118, TYR70, and SER122 present in the active site of AChE. © 2013 Wiley Periodicals, Inc.     

Topological analysis of the electronic charge density in the ethene protonation reaction catalyzed by acidic zeolite

In the present work, the distribution of the electronic charge density in the ethene protonation reaction by a zeolite acid site is studied within the framework of the density functional theory and the atoms in molecules (AIM) theory. The key electronic effects such as topological distribution of the charge density involved in the reaction are presented and discussed. The results are obtained at B3LYP/6-31G(**) level theory. Attention is focused on topological parameters such as electron density, its Laplacian, kinetic energy density, potential energy density, and electronic energy density at the bond critical points (BCP) in all bonds involved in the interaction zone, in the reactants, pi-complex, transition state, and alkoxy product. In addition, the topological atomic properties are determined on the selected atoms in the course of the reaction (average electron population, N(Omega), atomic net charge, q(Omega), atomic energy, E(Omega), atomic volume, v(Omega), and first moment of the atomic charge distribution, M(Omega)) and their changes are analyzed exhaustively. The topological study clearly shows that the ethene interaction with the acid site of the zeolite cluster, T5-OH, in the ethene adsorbed, is dominated by a strong O-H...pi interaction with some degree of covalence. AIM analysis based on DFT calculation for the transition state (TS) shows that the hydrogen atom from the acid site in the zeolitic fragment is connected to the carbon atom by a covalent bond with some contribution of electrostatic interaction and to the oxygen atom by closed shell interaction with some contribution of covalent character. The C-O bond formed in the alkoxy product can be defined as a weaker shared interaction. Our results show that in the transition state, the dominant interactions are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the concentration and accumulation of the charge density along the bond path between the nuclei linked increases.     

Electrostatic features of molecular recognition by cyclic urea mimics of chymotrypsin

The molecular electrostatic-potential pattern was used to investigate the electrostatic features of molecular recognition by two cyclic urea mimics of the active site of α-chymotrypsin. The structures of the mimics were obtained by molecular-mechanics evaluation of the conformational potential-energy surface of the molecules. Calculations were done by using two different atomic point-charge sets in order to assess the effect of charge on the electrostatic potential pattern. The molecules studies were: (1) a “full” mimic of chymotrypsin containing the hydroxyl, imidazole, and carboxylate anion functionalities typical of the active site of the enzyme, and (2) a “partial” mimic with only the hydroxyl and imidazole functional groups. Comparison of the molecular electrostatic-potential patterns of the two mimics in both charge sets showed that the largest differences were due to the structural addition of the carboxylate anion, rather than any particular differences in the choice of atomic point charge. For the full mimic, the pattern was essentially dominated by the negative charge on the carboxylate. Small structural changes which optimized the orientation of the catalytic components had little effect on the electrostatic potential pattern of the molecule. In the absence of the anionic functionality, greater differences were noted in the electrostatic potential pattern of the partial mimic in the two charge sets. The choice of atomic point charge was seen to influence the hydrogen-bonding pattern of the hydroxyl and imidazole moieties, resulting in differences in the spatial orientation of the electrostatic potential minima. In general, both charge sets produced molecular electrostatic-potential patterns which indicated that long-range electrostatic interactions would direct the cationic end of the substrate into the electron-rich binding site. However, specific local features of the electrostatic potential pattern were found to depend on point-charge set through the influence of charge on the hydrogen-bonding pattern.     

18-冠醚-6及其碱金属离子配合物的原子电荷及构象能的分子力学研究

本文推导了电负性电荷与静电势电荷的关系, 并把电负性电荷用于分子力学计算。我们的推导和计算结果表明, 电负性电荷也能很好的反映分子中的静电相互作用, 并用于分子力学计算, 从而我们为用电负性作为分子力学力场参数提供了理论根据。     

Structural character and energetics of tyrosyl radical formation by electron/proton transfers of a covalently linked histidine-tyrosine: a model for cytochrome C oxidase

The structural, energetic, and electronic and IR spectroscopic properties for a model of the cross-linked histidine-tyrosine (His-Tyr) residues as found in cytochrome c oxidase (CcO) are investigated by ab initio methods. The formation of a His-Tyr radical is studied by two paths: proton release followed by electron release and vice versa. The energetics for the proton/electron releases of the Tyr depend modestly on the cross-linked His substituent and, more sensitively, on the charge of the cation attached to the imino N site of the His residue. Protonation of the imino N site significantly increases the electron ionization potential and decreases the proton dissociation energy, making them competitive processes. A positive charge placed at the imino N site, whose value is scanned from zero to one, shows a continuous increase in ionization potential and a decrease in proton dissociation energy, with the +1 limit agreeing well with the protonated imino N site result, indicating a dominant electrostatic effect. The charge populations and the spin density distributions of the His-Tyr model, the radical cation formed by electron ionization, the anion formed by proton dissociation, and the final His-Tyr radical depend sensitively on the substituents, implying a modulation role on the charge transfer between the phenol and imidazole rings, especially for the charged species. His-Tyr and protonated His-Tyr exhibit differences among their respective structural isomers with consequences on their IR absorptions. Small barriers between their pseudo-cis and pseudo-trans rotamers demonstrate the relative flexibility between the two rings, and these may facilitate proton release and charge transfer. The cation effect demonstrates that the cationized cross-linked His-Tyr should be the best candidate to mimic the covalently ring-linked histidine-tyrosine structure in CcO.     

The effects of atomic multipole moments obtained by the potential-derived method on hydrogen bonding

A potential-derived atomic multipole method called the cumulative potential-derived atomic multipole method is developed, with which electrostatic atomic multipole moments are derived by fitting the molecular electric potential in a cumulative way. It is applied to the hydrides of N , O , F , S , Cl , and methanol and the hydrogen-bonded dimers formed between them. The relationship between atomic multipole moments and molecular charge distributions is found. The structures calculated with Buckingham's electrostatic model are in good agreement with experiments. The phenomena of nonlinear structures of most H -bonded complexes—the deviations of symmetry axes of electron donors from H bonds—and correct distinguishing between two alternative structures are attributed to atomic dipole and quadrupole moments. Compared with other methods, this method has a quantitative and qualitative advantage and simple algorithm. The main conclusion is that the atomic multipole moments play a substantial role, although a potential-derived charge model was deemed sufficient previously. © 1993 John Wiley & Sons, Inc.     

A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges

We report an approach for the determination of atomic monopoles of macromolecular systems using connectivity and geometry parameters alone. The method is appropriate also for the calculation of charge distributions based on the quantum mechanically determined wave function and does not suffer from the mathematical instability of other electrostatic potential fit methods.