应用ABEEM/MM蛋白质力场模型研究半胱氨酸二肽构象的性质

ABEEM/MM蛋白质力场模型是应用于蛋白质体系的原子-键电负性均衡方法(ABEEM)与力场(MM)相结合的浮动电荷模型.该模型能够准确地描述分子在环境变化时的静电极化,并能快速计算气态和溶液多肽的结构和能量.首次应用ABEEM/MM蛋白质力场模型研究半胱氨酸二肽构象的性质,如构象能、氢键等.此外,应用从头计算HF/6-31G**方法对其性质进行计算.ABEEM/MM蛋白质力场模型可以快速准确地得到半胱氨酸二肽分子不同稳定构象的性质,其结果可以和从头计算相媲美.以上研究有助于加深对半胱氨酸二肽构象性质的了解,从而也为进一步验证ABEEM/MM蛋白质力场模型的正确性以及参数的合理性提供可靠的依据.     

微过氧化物酶水溶液的ABEEM/MM动力学模拟

应用原子-键电负性均衡浮动电荷分子力场(ABEEM/MM), 对微过氧化物酶水溶液进行了分子动力学模拟. 研究了水溶液对微过氧化物酶的结构, 血红素的皱裂构象以及轴配体咪唑基的取向的影响. 结果表明, 在水溶液中微过氧化物酶的骨架氨基酸是稳定的, 而血红素的皱裂构象在水分子的作用下趋于平面. 与血红素轴配体咪唑基键连的组氨酸决定着咪唑基的空间取向, 而咪唑基与血红素侧链的丙酸基的静电作用对其取向仅起次要作用.     

ABEEMσπ和OPLS-AA力场对蛋白质GA88/GB88的分子动力学模拟

将原子与键电负性均衡方法融入分子力学方法,即利用ABEEMσπ浮动电荷力场与ABEEM-7P水模型相结合的方法及OPLS-AA固定电荷力场方法,对GA88和GB88蛋白进行了水溶液(温度295 K)和真空中的分子动力学模拟.比较两种方法得到的两个蛋白质的结构与实验结构的均方根偏差,分析了两种方法得到的两个蛋白质的回旋半径、氢键分布、径向分布及电荷分布情况.结果表明,ABEEMσπ和OPLS-AA力场均能正确模拟蛋白质结构,得到的各项偏差值接近,但从各偏差的波动大小可见,ABEEMσπ力场的模拟更稳定;回旋半径模拟很好地体现了蛋白质的"电致紧缩"现象;氢键分布、径向分布及电荷分布表明,与OPLS-AA固定电荷力场相比,ABEEMσπ浮动电荷力场能更好地体现蛋白质和周围水分子的极化效应.     

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

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

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

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

血红素近轴侧基氢键的ABEEM/MM分子动力学模拟

应用ABEEM/MM浮动电荷力场对鲸鱼肌红蛋白及突变体进行了分子动力学模拟. 结果表明, 血红素近轴侧基不存在稳定的双氢键, 该氢键对轴配体咪唑的取向不起决定性作用, 而咪唑的取向与键联的组氨酸有密切联系. 同时表明, 血红素轴配体的柔性与其邻近的氨基酸和咪唑体积有关.     

计算受体和配体结合自由能的新方法

利用ABEEMσπ浮动电荷力场与连续介质模型相结合的方法,计算了受体和配体的结合自由能.将结合自由能分解为真空中的力场作用项、溶剂化能量以及熵效应.由于ABEEMσπ/MM方法充分考虑了外界环境发生变化引起的体系中各个位点之间的电荷极化,因而极大地提高了结合自由能的计算精度.利用该方法计算的2个复合物的结合自由能与实验值的偏差均小于0.5kJ/mol.     

应用ABEEM/MM力场探究Na~+与氮-甲基乙酰胺及水的相互作用

应用量子化学方法 MP2/6-311++G(d,p)//B3LYP/6-31+G(d,p)和原子-键电负性均衡浮动电荷分子力场(ABEEM/MM),对[Na(H_2O)n]+(n=3～8),[Na(NMA)n]+(n=3～8)和[Na(NMA)n(H_2O)m]+(n+m=4,6)(NMA=氮-甲基乙酰胺)体系的结构、结合能和电荷分布进行研究.在计算结果的基础上,构建上述体系的ABEEM/MM可极化势能函数,优选并确定相关参数.结果表明,ABEEM/MM的计算结果与量子化学的计算结果相符:Na+与配体间距离的平均绝对偏差(AAD)小于0. 007 nm,相对均方根偏差(RRMSD)小于3. 5%,夹角的AAD小于2. 4°,RRMSD小于2. 0%,结合能的AAD小于8. 9 k J/mol,RRMSD小于12. 4%;ABEEM/MM电荷分布与量子力学(QM)电荷分布的线性相关系数在0. 97以上.     

ABEEMσп/MM力场模型下的多肽构象

在ABEEM/MM蛋白质浮动电荷力场模型的基础上,加入孤对电子和π电荷位点,从而能够体现多肽和蛋白质分子中一些重要的各向异性极化性质,允许非化学键方向的电子转移和极化.利用从头计算数据拟合模型相关参数.计算得到的小分子团簇结合能、偶极矩、氢键键长等性质与从头计算结果符合很好.该经典极化模型力场能够重复量子场下丙氨酸二肽、丙氨酸四肽、甘氨酸三肽的各稳定构象,其稳定性顺序与精密从头计算结果相一致,其结构和能量性质较以往模型有一定提高,并优于其他力场模型.     

原子-键电负性均衡方法融合进分子力场(ABEEM/MM)应用于蛋白质分子(Crambin)模拟

建立应用于多肽和蛋白质模拟的ABEEM/MM浮动电荷力场.利用该模型和参数,对实际蛋白质分子Crambin(植物种子中的一种小的蛋白质)进行模拟,得到了满意的结果,为其更广泛的应用开辟了道路.     

Study of peptide conformation in terms of the ABEEM/MM method

The ABEEM/MM model (atom-bond electronegativity equalization method fused into molecular mechanics) is applied to study of the polypeptide conformations. The Lennard-Jones and torsional parameters were optimized to be consistent with the ABEEM/MM fluctuating charge electrostatic potential. The hydrogen bond was specially treated with an electrostatic fitting function. Molecular dipole moments, dimerization energies, and hydrogen bond lengths of complexes are reasonably achieved by our model, compared to ab initio results. The ABEEM/MM fluctuating charge model reproduces both the peptide conformational energies and structures with satisfactory accuracy with low computer cost. The transferability is tested by applying the parameters of our model to the tetrapeptide of alanine and another four dipeptides. The overall RMS deviations in conformational energies and key dihedral angles for four di- or tetrapeptide, is 0.39 kcal/mol and 7.7 degrees . The current results agree well with those by the accurate ab initio method, and are comparable to those from the best existing force fields. The results make us believe that our fluctuating charge model can obtain more promising results in protein and macromolecular modeling with good accuracy but less computer cost.     

Sr2+和Ba2+水分子体系结构和结合能的从头算和ABEEM/MM研究

应用从头算方法和ABEEM/MM浮动电荷分子力场, 研究了水合碱土离子团簇Sr²⁺/Ba²⁺(H₂O)_n (n=1-6), 构建了离子-水相互作用的ABEEM/MM势能函数, 获得了水合离子团簇的稳定结构, 计算了结合能. 计算结果表明, ABEEM/MM方法的结果和从头算方法的结果有很好的一致性. 进一步应用ABEEM/MM对Sr²⁺和Ba²⁺水溶液进行了分子动力学模拟. 对Sr²⁺水溶液, 得到的Sr²⁺-水中氧原子的径向分布函数的第一和第二最高峰分别位于0.257和0.464 nm处, 第一和第二水合层的配位水分子数分别为9.2和11.4; 对Ba²⁺水溶液, 得到的Ba²⁺与水中氧原子的径向分布函数的第一和第二最高峰分别位于0.269和0.467 nm处, 第一和第二水合层的配位水分子数分别为9.9和12.4. 这与实验值或其它理论模拟结果有较好的一致性. 对比外层的水分子, 金属离子的极化作用使得溶液中第一水合层中水分子的O―H键长增长, HOH键角减小.     

从头算和ABEEMσπ/MM对(CH_3OH)_n(n=3～12)和[Na(CH_3OH)_n]~+(n=3～6)体系的研究

采用从头算方法(ab initio)和原子-键电负性均衡浮动电荷分子力场方法(ABEEMσπ/MM),对甲醇团簇(CH_3OH)n(n=3~12)和[Na(CH_3OH)_n]~+(n=3~6)体系的结构、电荷分布和结合能进行研究.依据从头算结果构建上述体系的ABEEMσπ/MM浮动电荷势能函数,并确定相关参数.结果表明,ABEEMσπ/MM所获得的结构和结合能等均优于OPLS/AA力场,并与从头算结果相符,其中键长的平均绝对偏差(AAD)小于0.004 nm,键长、键角和结合能的相对均方根偏差(RRMSD)分别小于3.8%,1.7%和6.8%;电荷分布与从头算结果的线性相关系数均大于0.99.     

Theoretical study on the hydration of hydrogen peroxide in terms of ab initio method and atom-bond electronegativity equalization method fused into molecular mechanics

In this paper, the interaction between hydrogen peroxide (HP) and water were systemically studied by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) and ab initio method. The results show that the optimized geometries, interaction energies and dipole moments of hydrated HP clusters HP(H₂O) _n (n = 1–6) calculated by ABEEM/MM model are fairly consistent with the MP2/aug-cc-pVTZ//MP2/aug-cc-pVDZ results. The ABEEM/MM results indicate that n = 4 is the transition state structure from 2D planar structure to 3D network structure. The variations of the average hydrogen bond length with the increasing number of water molecules given by ABEEM/MM model agree well with those of ab initio studies. Moreover, the radial distribution functions (RDFs) of water molecule around HP in HP aqueous solution have been analyzed in detail. It can be confirmed that HP is a good proton donor and poor proton acceptor in aqueous solution by analysis of the RDFs.     

An investigation of crambin and BPTI based on ABEEM/MM model

Molecular dynamics simulation studies on crambin,BPTI(298 K,in vacuo) have been performed by ABEEM/MM method. Some structural properties were discussed.The results show fair consistency with those from X-ray experiment,Moreover, ABEEM/MM model can properly describe the interactions of hydrogen bond of protein systems.     

Modelling of carbohydrate–aromatic interactions: <Emphasis Type="Italic">ab initio </Emphasis>energetics and force field performance

Summary Aromatic amino acid residues are often present in carbohydrate-binding sites of proteins. These binding sites are characterized by a placement of a carbohydrate moiety in a stacking orientation to an aromatic ring. This arrangement is an example of CH/π interactions. Ab initio interaction energies for 20 carbohydrate–aromatic complexes taken from 6 selected ultra-high resolution X-ray structures of glycosidases and carbohydrate-binding proteins were calculated. All interaction energies of a pyranose moiety with a side chain of an aromatic residue were calculated as attractive with interaction energy ranging from −2.8 to −12.3 kcal/mol as calculated at the MP2/6-311+G(d) level. Strong attractive interactions were observed for a wide range of orientations of carbohydrate and aromatic ring as present in selected X-ray structures. The most attractive interaction was associated with apparent combination of CH/π interactions and classical H-bonds. The failure of Hartree–Fock method (interaction energies from +1.0 to −6.9 kcal/mol) can be explained by a dispersion nature of a majority of the studied complexes. We also present a comparison of interaction energies calculated at the MP2 level with those calculated using molecular mechanics force fields (OPLS, GROMOS, CSFF/CHARMM, CHEAT/CHARMM, Glycam/AMBER, MM2 and MM3). For a majority of force fields there was a strong correlation with MP2 values. RMSD between MP2 and force field values were 1.0 for CSFF/CHARMM, 1.2 for Glycam/AMBER, 1.2 for GROMOS, 1.3 for MM3, 1.4 for MM2, 1.5 for OPLS and to 2.3 for CHEAT/CHARMM (in kcal/mol). These results show that molecular mechanics approximates interaction energies very well and support an application of molecular mechanics methods in the area of glycochemistry and glycobiology.     

ABEEMσπ/MM力场模型下的多肽构象

在ABEEM/MM蛋白质浮动电荷力场模型的基础上,加入孤对电子和 电荷位点,从而能够体现多肽和蛋白质分子中一些重要的各向异性极化性质,允许非化学键方向的电子转移和极化。利用从头计算数据拟合模型相关参数。计算得到的小分子团簇结合能、偶极矩、氢键键长等性质与从头计算结果符合很好。该经典极化模型力场能够重复量子场下丙氨酸二肽、丙氨酸四肽、甘氨酸三肽的各稳定构象,其稳定性顺序与精密从头计算结果相一致,其结构和能量性质较以往模型有一定提高,并优于其他力场模型。     

Ab initio and ABEEM/MM fluctuating charge model studies of dimethyl phosphate anion in a microhydrated environment

Dimethyl phosphate (DMP) anion has been used extensively as a model compound to simulate the properties of phosphate group. A 35-point DMP anion potential model is constructed based on the atom-bond electronegativity equalization fluctuating charge molecular force field (ABEEM/MM), and it is employed to study the properties of gas-phase DMP anion and DMP-(H₂O) _n (n = 1–3) clusters. The ABEEM/MM model reproduces well the properties obtained by available experiments and QM calculations, including charge distributions, geometries, and conformational energies of gas-phase DMP-water complexes. Furthermore, molecular dynamics simulation on the DMP anion in aqueous solution based on the ABEEM/MM shows that a remarkable first hydration shell around the nonesterified oxygen atom of DMP anion is formed with a coordination number of 5.2. It is also found that two hydrogen atoms of one water molecule form two hydrogen bonds with two nonesterified oxygen atoms of DMP anion simultaneously. This work could be used as a starting point for us to establish the ABEEM/MM nucleic acid force field.