应用ABEEMσπ确定大分子体系的电荷分布

以密度泛函理论（DFT）和电负性均衡原理为基础,明确处理了双键的结构,发展建立原子一键电负性均衡方法中的σπ模型（ABEEMσπ）,本模型将双键划分为一个σ键区域和四个π键区域（每个双键原子各有2个π键区域）,其中,σ电荷中心位于两成键原子之间共价半径之比处;π电荷中心垂直于双键所在平面,置于双键原子上下两侧。本文给出分子中各部分有效电负性的精密公式,以及参数确定方法,并应用该方法简捷快速地计算了C18N3O3H27C15NO2H21,以促甲状腺激素放激素TRH等有机和生物大分子体系的电荷分布,得到的结束可以和相应的从头算项媲美。     

应用ab initio方法和ABEEMσπ模型研究含铝金属酶

应用原子-键电负性均衡方法中的σπ模型(ABEEMσπ模型),通过大量量子化学计算,拟合确定了含铝金属酶体系的ABEEMσπ参数.将这些参数应用到含铝金属酶大分子体系的电荷分布及Fukui函数的计算,结果显示,ABEEMσπ模型计算得到的电荷分布及Fukui函数与从头算和实验结论均有很好的一致性.还进一步计算分析了1L3R酶与丝氨酸结合前后的分子各区域的电荷分布,结果表明,Al 3+是1L3R酶的活性中心,根据结合后分子的Fukui函数可以得出丝氨酸会使1L3R酶的活性降低.另外,通过比较两者结合前后Al 3+的广义Fukui函数,证明了广义Fukui函数可用于该体系分子间反应活性的比较,同时也说明利用ABEEMσπ模型来预测含铝金属酶的抑制剂是可行的.     

应用ABEEM_(σπ)模型研究非金属硼化物的电荷分布

通过大量量子化学计算,拟合确定了非金属硼化物体系模型分子的ABEEMσπ参数.将这些参数应用到ABEEMσπ模型中计算非金属硼化物模型分子的电荷分布,计算结果显示,ABEEMσπ模型计算得到的电荷分布与从头算计算的电荷分布都有很好的一致性.还计算了测试分子的电荷分布从而验证了ABEEMσπ参数的可转移性,电荷分布变化的规律是在BF3,BCl3,BBr3,BI3分子中B原子电荷逐渐减小,F,Cl,Br,I原子电荷逐渐增大.ABEEMσπ模型可以很好地应用于非金属硼化物的结构和性质的分析.     

应用ABEEM σπ模型研究含氮杂环化合物的电荷分布

应用在密度泛函理论和电负性均衡原理基础上发展的原子-键电负性均衡方法的σπ模型(即ABEEM σπ模型),并结合自编程序,对5种典型的含氮杂环化合物的电荷分布进行了计算. 计算结果表明,所得的电荷分布可以和从头计算结果有非常好的一致性,进一步验证了ABEEM σπ模型的合理性和可靠性     

并行程序实现ABEEM σπ模型电荷分布计算

对以密度泛函理论和电负性均衡原理为基础发展的原子-键电负性均衡方法中的σπ模型(ABEEM σπ模型)进行分析得出,利用该模型计算分子体系的电荷分布时,最耗费时间的部分是求解线性方程组.根据解线性方程组的串行程序,我们提出在并行环境下不带平方根的Cholesky分解方法.结果表明,利用改编后的并行程序能快速而准确地计算分子的电荷分布,算法随着矩阵规模的增大,并行效率也随之增高,即分子体系越大结果越理想.因此本算法适用于大规模问题的计算.     

应用ABEEMσπ模型研究氨基酸和多肽分子的电荷分布

在密度泛函理论和电负性均衡原理的理论框架下 ,发展了直接用于计算原子及键电荷分布的原子 键电负性均衡方法的σπ模型 ,本文应用此模型计算了氨基酸和一些较大多肽分子的电荷分布 .     

应用ABEEMσπ模型确定大分子体系的电荷分布

以密度泛函理论(DFT)和电负性均衡原理为基础,明确处理了双键的结构,发展建立原子一键电负性均衡方法中的σπ模型(ABEFMσπ).本模型将双键划分为一个σ键区域和四个π键区域(每个双键原子各有2个π键区域),其中,σ电荷中心位于两成键原子之间共价半径之比处;π电荷中心垂直于双键所在平面,置于双键原子上下两侧.本文给出分子中各部分有效电负性的精密公式,以及参数确定方法.并应用该方法简捷快速地计算了C     

组氨酸二肽-水团簇的Ab Initio和ABEEMσπ力场的研究

使用密度泛函B3LYP/6-311++G(d,p)方法对组氨酸二肽与水团簇的结构进行优化,在MP2/aug-cc-pVDZ水平下计算了这些体系的结合能,同时考虑了基组重叠误差(BSSE)和零点能(ZPE)校正.应用ABEEMσπ浮动电荷分子力场优化了组氨酸二肽与水分子所形成的团簇结构,计算了氢键键长和氢键键角,同时计算了组氨酸二肽与1~6个水分子所形成的团簇His(H2O)n(n=1~6)的结合能,探讨了氢键的协同效应.将ABEEMσπ浮动电荷分子力场、OPLSAA和AMBER力场所得的结果与从头算方法的结果进行了比较,ABEEMσπ的结果好于OPLS-AA和AMBER力场的,可与从头算方法所得到的结果相媲美.     

应用ABEEMσπ模型快速计算分子静电势

分子静电势是量子化学的理论指标,可以用来判断分子的某些性质,尤其是分子的反应活性.由于静电势由分子本身的性质决定,不同分子在其周围空间各点产生的静电势不同,因此可以采用静电势来描述分子本身的性质.建立了一种新的快速计算分子静电势的方法,并以HF、H2O、NH3和CH4小分子为例计算了他们的静电势.其中的电荷采用原子-键电负性均衡方法σπ模型（ABEEMσπ）计算的结果,结构由MP2/6-311＋＋G（d,p）方法优化所得.基于杨等人建立的内禀特征轮廓理论,绘制了相应的分子静电势分布图.结果与从头算符合得很好,为大分子静电势的计算打下基础.     

ABEEMσπ/MM模型中静电相互作用能的并行化

ABEEMσπ/MM模型程序中,计算静电相互作用能非常耗费机时.针对原串行程序中多个循环相互嵌套的求解部分,进行循环带状划分并行化处理.经测试表明,利用新编制的并行程序进行动力学模拟,并行加速比以线性趋势提高、求解静电相互作用能速度大幅度加快、尤其是针对原子数较多的分子体系效果比较理想.利用36个CPU,对于位点数为10 000左右的蛋白质体系,进行1ns的动力学模拟,至少可以节省1年左右的时间,明显地提高了研究蛋白质体系性质的效率.     

应用ABEEMσπ/MM方法计算小分子肽的总能量

把分子力场(MM)和原子-键电负性均衡方法(ABEEMσπ)协调融合在一起,采用ABEEMσπ/MM方法的分子能量计算方法,选取直链烷烃(n=1~10)、丙酮、丁酮、3-戊酮、N-甲基乙酰胺(NMA)、丙氨酸二肽和若干种氨基酸作为模型分子,计算模型分子能量,从而得到不同类型原子价态能量参数.应用这些原子价态能量参数,计算小分子肽的能量.把所得结果与从头算(MP2/6-311++G(d,p))方法计算的结果相比,得出小分子肽能量的绝对偏差小于9.000 0 kcal/mol,相对偏差小于19.000 0×10-6,且均方根偏差3.450 0 kcal/mol,相对均方根偏差6.919 0×10-6.以上结果表明,ABEEMσπ/MM方法计算的分子能量结果与从头算(MP2/6-311++G(d,p))方法计算的结果具有较好的一致性,此外,用ABEEMσπ/MM方法计算分子能量要远远快于从头算(MP2/6-311++G(d,p))方法.     

ABEEMσπ／MM 模型下8种类型氢键的氢键拟合函数

在ABEEMσπ/MM浮动电荷模型下,为了防止电荷过度极化,提高力场的准确性,引入了氢键拟合函数klp,H(Rlp,H)在氢键相互作用区域代替总校正系数k来描述氢键相互作用.选取了8个模型小分子二聚体,通过拟合在氢键两个区域之间不同距离时从头计算结果的结合能,来拟合ABEEMσπ/MM方法下的氢键拟合函数.应用该氢键拟合函数在ABEEMσπ/MM方法下计算小分子二聚体的结合能,与从头算方法得到的结合能有很好的一致性.同时对具有同样类型氢键的其他二聚体的结合能进行了计算,与从头算结果相近,说明这些函数具有良好的可转移性.     

用DFT和ABEEMσπ对不同形状单层石墨烯结构研究

选取完整的单层矩形(41个碳环)石墨烯构型和完整单层正六边形(37个碳环)石墨烯构型,采用B3LYP/6-31G(d,p)方法进行几何优化,计算完整石墨烯构型的HOMO和LUMO分布、能隙及电荷分布.研究结果显示,完整的矩形石墨烯的HOMO和LUMO只分布在锯齿形边缘,能隙约0.15 eV;正六边形石墨烯的HOMO和LUMO广泛地分布在整个片层上,能隙约2.0 eV.石墨烯的几何构型包含扶手椅形边缘和锯齿形边缘,然后定义了2种描述C原子的ABEEM参数,即石墨烯边缘碳原子、内部碳原子,以及与这2种碳原子相关的键.使用ABEEM参数计算出的电荷分布与HF/STO-3G方法所计算的电荷分布相比较,得到的线性相关系数都达到了0.99以上.认为所获得的碳原子和π键ABEEM参数是合理的,可以应用于大石墨烯分子体系.     

分子中电荷分布和分子电负性的预测

使用原子－键电负性均衡模型（ＡＢＥＥＭ），计算了一些任意选取的有机大分子的电荷分布和分子电负性。计算结果很好地重复了相应的从头计算结果，但却更加快捷、省时，展示了原子－键电负性均衡模型的广阔应用前景。     

应用ABEEM方法计算含金属离子Ga3+蛋白的电荷分布

首次开发ABEEM方法应用于含金属离子Ga3+蛋白体系的研究.ABEEM方法将分子电荷分解到了原子区域、σ键区域、π键区域和孤对电子区域.对金属离子Ga3+与蛋白之间的相互作用采用成键模型,Ga3+与配体原子之间有键电荷分布.通过蛋白晶体数据库的搜索,总结出Ga3+离子和蛋白相互作用的模型分子,确定了相关的新的电荷参数...     

稠杂环化合物ABEEM参数的拟合及简单应用

应用以密度泛函理论和电负性均衡原理为基础的原子-键电负性均衡方法中的σ-π模型（ABEEMσ-π）,通过大量的量子化学计算,利用最小二乘法,并结合我们自编的程序,拟合确定了稠杂环化合物的ABEEMσ-π各参数.进而,将这些参数应用到研究叶酸及维生素B2等稠杂环药物分子体系的电荷分布,结果显示由ABEEMσ-π模型获得的电荷分布与量子化学方法的计算结果有很好的一致性,线性相关系数均达到0.98以上.这说明我们拟合的参数是正确的,也进一步验证了ABEEMσ-π模型的合理性和可靠性.     

ABEEMσπ/MM模型对蛋白质G_A88/G_B88水溶液动力学性质的研究

应用ABEEMσπ/MM（σπ水平的原子与键电负性均衡方法融合进分子力学）浮动电荷模型以及显性ABEEM-7P水模型,对GA88和GB88两个蛋白质分子进行了分子动力学模拟.分析了2个蛋白质的动力学性质,包括蛋白质的回旋半径、疏水表面积和亲水表面积、各类原子位置的均方根偏差以及氢键分布.通过对比水溶液和真空下2个蛋白质的回旋半径,表明该模型很好地体现了蛋白质的＂电致紧缩＂现象;对疏水表面积和亲水表面积的计算表明,GB88中残基与溶剂的相互作用更强一些;非氢原子位置的均方根偏差及氢键分布情况与实验结构相比较表明,ABEEMσπ/MM浮动电荷模型模拟的GA88和GB88的结构与实验结构有很好的一致性,进而说明该模型的合理性和参数的可转移性.     

应用ABEEMσπ/GB-SA方法估算4个带电二肽的溶剂化自由能

采用B3LYP/6-311++G(d,p)理论水平优化获得体系的稳定几何结构.在相同的理论水平下,使用SMD溶剂模型,计算这些分子在水溶剂中的溶剂化自由能,并以其作为标准值.将各分子的稳定结构作为初始结构,在298 K、NVT系综下,应用ABEEMσπ可极化分子力场,对氨基酸侧链分子的水溶液进行分子动力学模拟,模拟的时间为1ns.平衡后的结构作为训练集,以调节和确定ABEEMσπ的相关参数,最终得到与从头计算相一致的溶剂化自由能.使用相同的参数,计算带电氨基酸二肽的溶剂化自由能,ABEEMσπ的计算结果与从头计算结果具有良好的一致性,表明该参数具有可转移性,这为研究多肽和蛋白质分子的溶剂化自由能奠定了基础.     

应用原子键电负性均衡方法σπ模型研究顺式丙烯酸的分子内氢键和构象

在密度泛函理论框架下发展的原子一键电负性均衡方法σπ模型(ABEEMσπ),在探讨有机分子体系的结构和性质等方面取得了令人满意的成果．本利用该模型计算了顺式丙烯酸不同构象的氢键值．根据氢键能与构象稳定性的关系,预测并解释了构象的稳定性,并与从头计算结果相比较,取得了很好的一致性．     

应用ABEEM/MM方法研究尿素与缬氨酸二肽间的相互作用

对尿素与二肽的相互作用的研究可为探讨尿素致蛋白质变性的机理提供基础.应用ABEEM/MM方法,对尿素和缬氨酸二肽间的相互作用进行了理论研究.结果 表明,尿素对缬氨酸二肽的几何结构影响很大,使其骨架二面角发生了明显的扭转.此外,用该方法计算的电荷分布和键长、键角、二面角等结构参数,都与从头算的结果有很好的一致性.