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

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

二极管特性简化模型在染料敏化太阳能电池电流-电压曲线拟合及在光电特性参数评估中的应用

利用瞬态光电导技术研究了在一个大的偏压范围内染料敏化太阳能电池(DSSC)中的电荷收集和电荷复合过程的动力学. 结果表明, 在不同的电压下电荷收集速率远快于电荷复合速率, 用以解释电荷的收集效率几乎不变. 基于这个结果, 简化了DSSC二极管特性模型, 从而实现了对电流-电压(I-V)曲线的直接拟合. 利用这一模型拟合的结果提取出一系列与工作条件下DSSC光电转化过程相关的参数, 可以用以描述包括电荷生成、电荷收集和复合以及DSSC器件的整流特性等关键性质. 将这一拟合方法应用于不同叔丁基吡啶(TBP)浓度电解液的电池, 获得的结果表明, 不同TBP浓度可以导致电荷复合速率有较大的差别, 但对电荷的收集速率影响甚小, 这与I-V曲线拟合的结果非常吻合. 同时研究证明TBP浓度对电池的理想因子(m)的影响较为明显,即高浓度TBP的DSSC对应大的m值以及较慢的电荷复合速率.     

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

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

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

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

有机污染物在碳纳米管吸附的定量结构-性质关系

对一系列共59个芳香性分子进行了HF/6-31G*水平上的结构优化,并在优化结构上进行了分子静电势及其导出参数的计算,应用多元线性回归方法建立了碳纳米管吸附有机污染物的平衡常数与分子结构间的定量关系.结果表明,分子表面静电势参数(Vmin、σ2+和ΣV+ind)结合分子表面积(S)和最低空轨道能级(εLUMO)可以很好地用于构建碳纳米管吸附的定量结构-性质关系(QSPR)模型.模型中引入的参数均具有明确的物理意义,其合理性可以从污染物与碳纳米管或水分子间相互作用的角度进行解释.模型的稳定性和预测能力经"留一法"和Monte Carlo交叉验证法进行了确证.本文亦采用支持向量机(SVM)、最小二乘支持向量机(LSSVM)和高斯过程(GP)等三种方法建立了上述参数与碳纳米管吸附性质的非线性模型.SVM和LSSVM模型表现出强的拟合能力,但预测能力明显不如其他模型.GP模型无论是拟合能力还是预测能力都是最佳,但并没有明显地优于线性模型,说明对本文研究体系而言,其分子结构与性质间的关系主要以线性形式存在.     

用电荷云模型计算氢原子和氢分子体系的能量

一般原子、分子体系的能量是通过解薛定谔方程式得到的,这种方法需要高深的数学推导,对初学者来说是不好理解的.下面介绍一种计算原子、分子体系能量的简单方法——原子分子结构的电荷云模型. 一、氢原子体系的能量氢原子电荷云模型(见图1)是由一个带正电的原子核和围绕着它的半径为R的电荷云组成的.电荷云的电荷分布是均匀的,即在整个球内任一点处的电荷密度ρ(r)=-q(4/3πR~3)~(-1).     

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

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

支持向量机用于多氯代萘毒性的定量构效研究

用偏最小二乘法(PLS)和留一交叉验证从90多个量子化学参数中筛选出极化率、分子量、部分原子上的净电荷、静电势等作为描述符,应用支持向量机(SVM)对20个多氯代萘同系物的三组毒性数据分别建立了定量构效关系模型.所得模型的交叉验证相关系数的平方分别为0.805、0.890、0.936.并将偏最小二乘法建模所得结果与之进行比较,结果表明,SVM预报能力优于PLS.     

紧束缚模型方法及其在DNA分子中的应用

近年来,紧束缚模型方法被广泛应用于计算生物大分子体系.本文从第一性原理出发,根据紧束缚近似的思想,推导出生物大分子体系中的单电子运动方程.在此基础上给出了紧束缚模型方法中所涉及参数(在位能和迁移积分)的计算公式,在理论上完善了紧束缚模型方法.我们将所提出的参数化方法应用于理想B型DNA分子,给出了各种序列组合下的在位能和迁移积分.此外,我们还计算了周期性DNA分子poly(A)-poly(T)和poly(G)-poly(C)中空穴在位能和迁移积分随格点间距离的变化,为改进现有的SSH极化子模型提供了新的思路,有助于DNA中电荷输运的极化子机理的研究.     

结合三维静电势参数研究化合物溶解度与溶剂分子结构间的定量关系

对一系列溶剂分子进行了结构优化和静电势及其导出参数的计算,运用多元线 性回归方法对10种化合物萘、菲、蒽、联苯、苊、六氯苯、苯偶酰、噻吨-9-酮、 二苯砜和敌草隆的溶解性能与溶剂分子的结构参数进行了关联。结果表明：分子表 面最正和最负的静电势V_(s,max)和V_(s,min)、电荷分离度Π以及分子的静电相互 作用趋势量τ这四个三维静电势参数,加上分子的前线轨道能级ε_(HOMO)和ε_ (LUMO)能很好地用于表达这些化合物在不同纯溶剂中的溶解度与溶剂分子结构间的 定量关系。     

Atomic multipoles: Electrostatic potential fit,local reference axis systems,and conformational dependence

Currently, all standard force fields for biomolecular simulations use point charges to model intermolecular electrostatic interactions. This is a fast and simple approach but has deficiencies when the electrostatic potential (ESP) is compared to that from ab initio methods. Here, we show how atomic multipoles can be rigorously implemented into common biomolecular force fields. For this, a comprehensive set of local reference axis systems is introduced, which represents a universal solution for treating atom‐centered multipoles for all small organic molecules and proteins. Furthermore, we introduce a new method for fitting atomic multipole moments to the quantum mechanically derived ESP. This methods yields a 50–90% error reduction compared to both point charges fit to the ESP and multipoles directly calculated from the ab initio electron density. It is shown that it is necessary to directly fit the multipole moments of conformational ensembles to the ESP. Ignoring the conformational dependence or averaging over parameters from different conformations dramatically deteriorates the results obtained with atomic multipole moments, rendering multipoles worse than partial charges. © 2012 Wiley Periodicals, Inc.     

Various atomic charge calculation schemes of CoMFA on HIF‐1 inhibitors of moracin analogs

Selection of appropriate partial charges in a molecule is crucial to derive good quantitative structure–activity relationship models. In this work, several partial atomic charges were assigned and tested in a comparative molecular field analysis (CoMFA) models. Many CoMFA models were generated for a series of hypoxia inducible factor 1 (HIF‐1) inhibitors using various partial atomic charges including charge equalization, Mülliken population analysis (MPA), natural population analysis, and electrostatic potential (ESP)‐derived charges. These atomic charges were investigated at various theoretical levels such as empirical, semiempirical, Hartree–Fock (HF), and density functional theory (DFT). Among them, Merz‐Singh‐Kollman (MK) ESP‐derived charges at the level of HF/6‐31G* gave the highest predictive q² with experimental pIC₅₀ values. With this charge scheme, a detailed analysis of CoMFA model was performed to understand the electrostatic interactions between ligand and receptor. More elaborate charge calculation schemes such as HF and DFT correlated more strongly with activity than empirical or semiempirical schemes. The choice of optimization methods was important. As geometries were fully optimized at the given levels of theory, the aligned structures were different. They differed considerably, especially for the flexible parts. This was likely the source of the substantial variation of q² values, even when the same steric factor was considered without electrostatic parameters. ESP‐derived charges were most appropriate to describe CoMFA electrostatic interactions among MPA, NBA, and ESP charges. Overall q² values vary considerably (0.8–0.5) depending on the charge schemes applied. The results demonstrate the need to consider more appropriate atomic charges rather than default CoMFA charges. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A numerically stable restrained electrostatic potential charge fitting method

Inspired by the idea of charge decomposition in calculation of the dipole preserving and polarization consistent charges (Zhang et al., J. Comput. Chem. 2011, 32, 2127), we have proposed a numerically stable restrained electrostatic potential (ESP)‐based charge fitting method for protein. The atomic charge is composed of two parts. The dominant part is fixed to a predefined value (e.g., AMBER charge), and the residual part is to be determined by restrained fitting to residual ESP on grid points around the molecule. Nonuniform weighting factors as a function of the dominant charge are assigned to the atoms. Because the residual part is several folds to several orders smaller than the dominant part, the impact of ill‐conditioning is alleviated. This charge fitting method can be used in quantum mechanical/molecular mechanical (QM/MM) simulations and similar studies, where QM calculated electronic properties are frequently mapped to partial atomic charges. © 2012 Wiley Periodicals, Inc.     

Conformational preferences for hydroxyl groups in substituted tetrahydropyrans

Quantum mechanical (ab initio and semiempirical) and force field calculations are reported for representative torsion potentials in several tetrahydropyran derivatives. The overall agreement between the various methods is quite good except that the AMBER torsion profiles are sensitive to the choice of atomic point charges. Using electrostatic potential (ESP) derived atomic point charges determined with the STO-3G basis set we find that AMBER is able to match the best quantum mechanical results quite well. However, when the point charges are derived using the 6-31G* basis set we find that scaling the intramolecular electrostatic nonbond interactions is necessary. AM1 does not work very well for these compounds when compared to the ab initio methods and, therefore, should only be used in cases when ab initio calculations would be prohibitive. Based upon our results we feel that any force field that makes use of 6-31G* ESP derived atomic point charges will need to scale intramolecular interactions. Implications of scaling intramolecular interactions to the development of force fields based on 6-31G* ESP derived atomic point charges are discussed. © 1992 by John Wiley & Sons, Inc.     

Building molecular charge distributions from fragments: Application to HIV-1 protease inhibitors

Interaction energies are a function of the molecular charge distribution. In previous work, we found that the set of atomic partial charges giving the best agreement with experimental vacuum dipole moments were from density functional theory calculations using an extended basis set. Extension of such computations to larger molecules requires an atomic partial charge calculation beyond present computational resources. A solution to this problem is the calculation of atomic partial charges for segments of the molecule and reassociation of such fragments to yield partial charges for the entire molecule. Various partitions and reassociation methods for five molecules relevant to HIV-1 protease inhibitors are examined. A useful method of reassociation is introduced in which atomic partial charges for a large molecule are computed by fitting to the combined electrostatic potential calculated from the fragment partial charges. As expected, the best sites for partitions are shown to be carbon—carbon rather than carbon—nitrogen bonds. © 1997 by John Wiley & Sons, Inc.     

Communication: Quantum polarized fluctuating charge model: a practical method to include ligand polarizability in biomolecular simulations

We present a simple and practical method to include ligand electronic polarization in molecular dynamics (MD) simulation of biomolecular systems. The method involves periodically spawning quantum mechanical (QM) electrostatic potential (ESP) calculations on an extra set of computer processors using molecular coordinate snapshots from a running parallel MD simulation. The QM ESPs are evaluated for the small-molecule ligand in the presence of the electric field induced by the protein, solvent, and ion charges within the MD snapshot. Partial charges on ligand atom centers are fit through the multi-conformer restrained electrostatic potential (RESP) fit method on several successive ESPs. The RESP method was selected since it produces charges consistent with the AMBER/GAFF force-field used in the simulations. The updated charges are introduced back into the running simulation when the next snapshot is saved. The result is a simulation whose ligand partial charges continuously respond in real-time to the short-term mean electrostatic field of the evolving environment without incurring additional wall-clock time. We show that (1) by incorporating the cost of polarization back into the potential energy of the MD simulation, the algorithm conserves energy when run in the microcanonical ensemble and (2) the mean solvation free energies for 15 neutral amino acid side chains calculated with the quantum polarized fluctuating charge method and thermodynamic integration agree better with experiment relative to the Amber fixed charge force-field.     

On combining Thole's induced point dipole model with fixed charge distributions in molecular mechanics force fields

The Thole induced point dipole model is combined with three different point charge fitting methods, Merz–Kollman (MK), charges from electrostatic potentials using a grid (CHELPG), and restrained electrostatic potential (RESP), and two multipole algorithms, distributed multipole analysis (DMA) and Gaussian multipole model (GMM), which can be used to describe the electrostatic potential (ESP) around molecules in molecular mechanics force fields. This is done to study how the different methods perform when intramolecular polarizability contributions are self‐consistently removed from the fitting done in the force field parametrization. It is demonstrated that the polarizable versions of the partial charge models provide a good compromise between accuracy and computational efficiency in describing the ESP of small organic molecules undergoing conformational changes. For the point charge models, the inclusion of polarizability reduced the the average root mean square error of ESP over the test set by 4–10%. © 2015 Wiley Periodicals, Inc.     

Clarification of the role of protein in carbonmonoxy myoglobin by investigating electronic states

This article reports the proton tautomerization effects of distal histidine residues in carbonmonoxy myoglobin according to the density functional calculations of the whole protein. The electron eigenstates and electrostatic potential (ESP) distributed around heme and its pocket vary significantly depending on the protonation positions of the distal histidine residues. To investigate the range over which the electronic structures are affected by the proton tautomerization, the quantum mechanics/molecular mechanics (QM/MM) method is applied to probe the QM size to reproduce the atomic partial charges and ESP around the active center. Consequently, we show that these properties converged for the 300 pm QM/MM system in this study. During the analysis, we also find that amino residues such as Phe43, Val68, and Phe138 interact strongly with heme through orbital mixing, indicating that the protein is a medium not only interacting with the reaction center, but also buffering on electrons. © 2013 Wiley Periodicals, Inc.