ABEEM浮动电荷分子力场在含离子体系中的发展和应用

离子水合及生物分子体系内离子选择性的微观作用机制是人们长期探索的重要课题,其难点在于如何合理精确地描述上述体系内的离子-水、离子-生物分子等各种相互作用.本文主要总结近年来原子-键电负性均衡浮动电荷分子力场（ABEEM/MM）在含离子体系中的发展和应用,包括离子水溶液、金属蛋白、离子-核酸碱基体系的研究.我们优选相关参数,构建上述体系的势能函数,并对气相水合离子团簇、离子水溶液、金属蛋白、离子-核酸碱基体系进行研究,模拟其结构、活性、热力学和动力学等性质.研究和比较结果表明,我们的ABEEM浮动电荷力场总体上优于其它力场方法,其精度可达到或接近高水平从头计算MP2方法.这为进一步探讨生物分子体系内的离子选择性、金属酶及其它含离子体系的结构和性质奠定了基础.     

Valence state parameters of all transition metal atoms in metalloproteins—development of ABEEMσπ fluctuating charge force field

To promote accuracy of the atom‐bond electronegativity equalization method (ABEEMσπ) fluctuating charge polarizable force fields, and extend it to include all transition metal atoms, a new parameter, the reference charge is set up in the expression of the total energy potential function. We select over 700 model molecules most of which model metalloprotein molecules that come from Protein Data Bank. We set reference charges for different apparent valence states of transition metals and calibrate the parameters of reference charges, valence state electronegativities, and valence state hardnesses for ABEEMσπ through linear regression and least square method. These parameters can be used to calculate charge distributions of metalloproteins containing transition metal atoms (Sc‐Zn, Y‐Cd, and Lu‐Hg). Compared the results of ABEEMσπ charge distributions with those obtained by ab initio method, the quite good linear correlations of the two kinds of charge distributions are shown. The reason why the STO‐3G basis set in Mulliken population analysis for the parameter calibration is specially explained in detail. Furthermore, ABEEMσπ method can also quickly and quite accurately calculate dipole moments of molecules. Molecular dynamics optimizations of five metalloproteins as the examples show that their structures obtained by ABEEMσπ fluctuating charge polarizable force field are very close to the structures optimized by the ab initio MP2/6–311G method. This means that the ABEEMσπ/MM can now be applied to molecular dynamics simulations of systems that contain metalloproteins with good accuracy. © 2014 Wiley Periodicals, Inc.     

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.     

Theoretical Study on Co^3＋ in Aqueous Solution in Terms of ABEEM/MM Model

A detailed theoretical investigation on Co^3＋ hydration in aqueous solution has been carded out by means of molecular dynamics （MD） simulations based on the atom-bond electronegativity equalization method fused into molecular mechanics （ABEEM/MM）. The effective Co^3＋ ion-water potential has been constructed by fitting to ab initio structures and binding energies for ionic clusters. And then the ion-water interaction potential was applied in combination with the ABEEM-7P water model to molecular dynamics simulations of single Co^3＋（aq.） solution, managing to reproduce many experimental structural and dynamical properties of the solution. Here, not only the common properties （radial distribution function, angular distribution function and solvation energy） obtained for Co^3＋ in ABEEM-7P water solution were in good agreement with those from the experimental methods and other molecular dynamics simulations but also very interesting properties of charge distributions, geometries of water molecules, hydrogen bond, diffusion coefficients, vibrational spectra are investigated by ABEEM/MM model.     

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

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

Theoretical studies on the hydration of formic acid by ab initio and ABEEMσπ fluctuating charge model

The interaction between formic acid (FA) and water was systemically investigated by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEMσπ/MM) and ab initio methods. The geometries of 20 formic acid–water complexes (FA–water) were obtained using B3LYP/aug-cc-pVTZ level optimizations, and the energies were determined at the MP2/aug-cc-pVTZ level with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) corrections. The ABEEMσπ potential model gives reasonable properties of these clusters when compared with the present ab initio data. For interaction energies, the root mean square deviation is 0.74 kcal/mol, and the linear coefficient reaches 0.993. Next, FA in aqueous solution was also studied. The hydrogen-bonding pattern due to the interactions with water has been analyzed in detail. Furthermore, the ABEEMσπ charges changed when H₂O interacted with the FA molecule, especially at the sites where the hydrogen bonds form. These results show that the ABEEMσπ fluctuating charge model is fine giving the overall characteristic hydration properties of FA–water systems in good agreement with the high-level ab initio calculations.     

Molecular-dynamics simulations of alkaline-earth metal cations in water by atom-bond electronegativity equalization method fused into molecular mechanics

Intermolecular potential for alkaline-earth metal (Be(2+), Mg(2+), and Ca(2+)) cations in water has been derived using the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), and it is consistent with what was previously applied to the hydration study of the monovalent cations. Parameters for the effective interaction between a cation and a water molecule were determined, reproducing the ab initio results. The static, dynamic, and thermodynamic properties of Be(2+)(aq), Mg(2+)(aq), and Ca(2+)(aq) were studied using these potential parameters. Be(2+) requires a more complicated form of the potential function than Mg(2+) and Ca(2+) in order to obtain better fits. Strong influences of the twofold charged cations on the structures of the hydration shells and some other properties of aqueous ionic solutions are discussed and compared with the results of a previous study of monovalent cations in water. At the same time, comparative study of the hydration properties of each cation is also discussed. This work demonstrates that ABEEM/MM provides a useful tool in the exploration of the hydration of double-charged cations in water.     

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

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

New force field parameters for metalloproteins I: Divalent copper ion centers including three histidine residues and an oxygen‐ligated amino acid residue

Transition metal ion complexation with proteins is ubiquitous across such diverse fields as neurodegenerative and cardiovascular diseases and cancer. In this study, the structures of divalent copper ion centers including three histidine and one oxygen‐ligated amino acid residues and the relative binding affinities of the oxygen‐ligated amino acid residues with these metal ion centers, which are debated in the literature, are presented. Furthermore, new force field parameters, which are currently lacking for the full‐length metal‐ligand moieties, are developed for metalloproteins that have these centers. These new force field parameters enable investigations of metalloproteins possessing these binding sites using molecular simulations. In addition, the impact of using the atom equivalence and inequivalence atomic partial charge calculation procedures on the simulated structures of these metallopeptides, including hydration properties, is described. © 2014 Wiley Periodicals, Inc.     

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

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

从头算和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.     

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浮动电荷力场应用于血红素结构的研究

运用拟合的参数, 应用ABEEM/MM浮动电荷力场对血红素分子结构进行了模拟. 结果表明, 该力场与CHARMM力场相比, 能更好地模拟晶体结构. 计算的ruffing构象能与B3LYP/6-31G*计算结果的线性相关系数在0.98以上, 同时表明血红素分子中twist-angle对ruffing构象具有明显影响. ABEEM/MM力场计算的细胞色素c552中血红素分子的电荷分布与CHARMM固定电荷力场的比较, 更准确地反映了血红素分子的电荷分布以及极化现象.     

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.     

Study of water clusters in the n = 2–34 size regime, based on the ABEEM/MM model

Various properties of water clusters in the n = 2–34 size regime with the change of cluster size have been systemically explored based on the newly developed flexible-body and charge-fluctuating ABEEM/MM water potential model. The ABEEM/MM water model is to take ABEEM charges of all atoms, bonds, and lone-pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The computed correlating properties characterizing water clusters (H₂O) _n (n = 2–34) include optimal structures, structural parameters, ABEEM charge distributions, binding energies, hydrogen bonds, dipole moments, and so on. The study of optimal structures shows that the ABEEM/MM model can correctly predict the following important structural features, such as the transition from two-dimensional (from dimer to pentamer) to three-dimensional (for clusters larger than the hexamer) structures at hexamer region, the transition from cubes to cages at dodecamer (H₂O)₁₂, the transition from all-surface (all water molecules on the surface of the cluster) to one water-centered (one water molecule at the center of the cluster, fully solvated) structures at (H₂O)₁₇, the transition from one to two internal molecules in the cage at (H₂O)₃₃, and so on. The first three structural transitions are in good agreement with those obtained from previous work, while the fourth transition is different from that identified by Hartke. Subsequently, a systematic investigation of structural parameters, ABEEM charges, energetic properties, and dipole moments of water clusters with increasing cluster size can provide important reference for describing the objective trait of hydrogen bonds in water cluster system, and also provide a strong impetus toward understanding how the water clusters approach the bulk limit.     

Application of the ABEEM/MM model in studying the properties of the water clusters (H2O) n (n=7−10)

ABEEM/MM model has been applied to compute the various properties characterizing water clusters (H₂O) _n (n = 7?10), such as optimized geometries, the hydrogen bonds number, cluster interaction energies, stabilities, ABEEM charge distributions, dipole moments, structural parameters, and so on, and to describe the transition reflected by the hexamer region from two-dimensional (from dimer to pentamer) to three-dimensional structures (for clusters larger than the hexamer).     

鸟嘌呤四链体中K+移动路径的预测

结合鸟嘌呤四链体(G4)的结构特点,取G4中的一个K+与一片g-tetrad(g-tetrad-K+)为研究对象,和一个K+与两片g-tetrad(1,2-g-tetrads-K+和2,3-g-tetrads-K+)为研究对象.为g-tetrad-K+,1,2-g-tetrads-K+和2,3-g-tetrads-K+三体系中的K+设计了几种移动路径,来探讨G4中K+与g-tetrad的作用.首先分别应用从头计算MP2/6-31G(d,p)和ABEEMσπ/MM(σπ水平的原子键电负性均衡原理融入到分子力学)方法,对g-tetrad-K+中K+移动过程中体系的结构和能量变化进行了详细研究.结果表明ABEEMσπ/MM方法能够准确地描述g-tetrad-K+体系的结构、结合能等性质.为了更接近实际体系,进而用ABEEMσπ/MM方法预测了更大的体系:一个K+和两片g-tetrad.由于体系较大,MP2方法无法计算其结合能.ABEEMσπ/MM方法对g-tetrad-K+,1,2-g-tetrads-K+和2,3-g-tetrads-K+体系中K+移动而产生的结合能变化表明:如果吸收足够的能量,G4中K+最容易沿α方向移出.     

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.