水合锰(Ⅱ)结构的量子化学和ABEEM/MM研究

应用新一代可极化分子力场——原子-键电负性均衡浮动电荷分子力场ABEEM/MM,结合精密量子化学方法, 构建了精确的Mn²⁺-H₂O 相互作用的势能函数, 确定了相关参数.将该势能函数用于计算[Mn(H₂O)_n]²⁺(n=1~12)的结构和结合能, 得到了与量子化学一致的结果. 进一步对Mn²⁺水溶液进行ABEEM/MM 动力学模拟, 得到的Mn²⁺–O 径向分布函数的第一和第二最高峰分别处于0.218 和0.435 nm 处, 积分得到第一和第二水合层的配位水分子数分别为7.03 和17.74; 对于O–Mn²⁺–O 角度分布函数, 其第一和第二最高峰分别位于80°和140°附近, 这些结果与实验和其他理论方法的结果有很好的一致性. Mn²⁺的极化作用使得第一水合层中水分子的键长明显增长, 键角明显减小; 而Mn²⁺对第二水合层及外层水分子的结构影响较小. 分析体系的电荷分布表明, 与ABEEM-7P 纯水相比, Mn²⁺水溶液中参与形成氢键的氢原子和孤对电子的电荷变化较大, 且Mn²⁺和其邻近的水分子间存在明显的电荷转移.     

原子-键电负性均衡融合进分子力场对CoCl3水溶液的分子动力学模拟研究

298．15K时运用原子-键电负性均衡融合进分子力场（ABEEM／MM）对CoCl3水溶液进行了分子动力学模拟．结果表明，ABEEM／MM力场模型可以获得体系的动态电荷，充分考虑了体系中存在的极化现象，从而及时地调节改变的静电场，真实地模拟了溶液的微观结构；并且CoCl3溶液中钴离子更倾向于和氯离子配位，而使钴离子与水分子的配位能力降低．     

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

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

HIV-1病毒中NCp7锌指区域和配体的电荷分布研究

采用B3LYP/6-311++G(d,p)方法,优化了4个锌指蛋白分子和31个模型分子的几何构型,使用HF/STO-3G计算了优化后体系的Mulliken电荷分布.使用线性回归和最小二乘法,拟合确定了ABEEMσπ参数(包括参考电荷,价态电负性和价态硬度).使用ABEEMσπ参数,计算获得了4个锌指蛋白和11个配体的电荷分布.ABEEMσπ方法所计算的电荷分布与从头计算方法的线性相关方程斜率k接近1.000,截距b接近于0.000,线性相关系数R在0.943 4以上.由此可见,ABEEMσπ方法与从头计算方法相比有很好的一致性.验证了笔者所拟合的ABEEMσπ参数是合理的和可转移的,可应用于类似体系的电荷分布的计算.     

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

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

碱基与氨基酸二肽间的氢键相互作用

对于碱基-蛋白质之间的相互作用,氢键的形成及断裂起着重要的作用.选取7种氨基酸二肽和3种碱基,以它们之间形成单根氢键的体系为研究对象,在ab initioB3LYP/6-311++G(d,p)水平下进行结构优化,MP2/6-311++G(d,p)水平下进行能量的计算.同时应用ABEEMσπ/MM方法优化结构和计算能量.通过改变氢键相互作用区域的静电相互作用参数kH-bond和形成氢键的氢原子与其受体原子的孤对电子之间距离Rlp,H来拟合函数kH-bond(Rlp,H).用所得函数对选取的复合物进行结合能计算,其结果与ab initio方法所得的结果符合的很好.     

组氨酸二肽-水团簇的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浮动电荷模型（原子一键电负性均衡方法融合进分子力场）对气态丝氨酸和苏氨酸残基二肽分子构象进行了初步的研究,与经典的力场模型相比,该方法中的静电势包含了分子内和分子间的静电极化作用,以及分子内电荷转移影响,同时加入了化学键等非原子中心电荷位点,合理的体现了分子中的电荷分布。相对其他极化力场模型,该模型具有计算量较小的特点。结果表明：我们模型计算的两种二肽分子的构象能和关键二面角结构与从头计算结果符合得很好,优于其他力场模型。