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

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

ABEEM/MM浮动电荷力场应用于血红素结构的研究

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

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

应用ABEEMσπ/MM模型对重组人纤溶酶原Kringle 1结构域与配体分子的对接计算

应用原子-键电负性均衡方法的浮动电荷分子力场(ABEEMσπ/MM),对重组人纤溶酶原Kringle1结构域(K1Pg)与配体ε-Aminocaproic acid(EACA),trans-4-(Aminomethyl)cyclohexane-1-carboxylic acid(AMCHA),L-Lysine(Lys),7-Aminoheptanoic acid(7-AHA)和Benzylamine进行了半柔性对接计算.用ABEEMσπ/MM模型优化得到的复合物K1Pg/EACA和K1Pg/AMCHA的结构很接近实验晶体结构.此外通过对5种配体与K1Pg结合能的计算,得出5种配体与K1Pg结合能的大小顺序为AMCHAEACA7-AHALysBenzy-lamine,与实验中测得5种配体与K1pg的平衡结合常数K值大小顺序相一致.     

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

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

碱金属与碱土金属卤化物电荷分布和偶极矩的ABEEMσπ模拟

应用原子键电负性均衡方法模型(ABEEMσπ模型),通过大量量子化学计算,拟合确定了含碱金属和碱土金属的卤化物体系的ABEEMσπ参数.ABEEMσπ模型计算得到的电荷分布与从头算计算的电荷分布有很好的一致性.同时计算了模型分子的偶极矩,与实验值具有很好的一致性.研究表明ABEEMσπ模型可以很好地应用于含碱金属和碱土金属的卤化物的结构和性质的分析.     

鸟嘌呤四链体中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+最容易沿α方向移出.     

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

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

致癌性胺离子反应过程中电荷和化学位移的计算

应用原子-键电负性均衡方法中的σπ模型(ABEEMσπ)计算了由致癌性胺离子所参与的反应过程中的电荷分布,所计算出的电荷分布可以和从头算很好的相关联,并且所需要的时间也大大的缩短;同时应用从头计算程序计算致癌性胺离子反应过程中的NMR化学位移.结果表明,在反应过程中电荷的变化和NMR化学位移的变化有很好的对应关系.     

碱基与氮甲基乙酰胺相互作用的从头算和ABEEMσπ研究

选取氮甲基乙酰胺(N-Methylacetamide,NMA)和DNA碱基为研究对象,应用ABEEMσπ/MM方法研究了NMA与碱基之间形成的二聚体的几何构型、电荷分布及结合能等性质,并进行了相应的从头算(abinitio)MP2水平研究,在B3LYP/6-311 ++G(d,p)水平下对结构进行了优化.将2种方法的研究结果进行比较发现,所获得的构型、电荷及结合能之间具有很好的一致性,NMA与DNA碱基相互作用的强度大小依次为:鸟嘌呤＞胸腺嘧啶＞胞嘧啶＞腺嘌呤,且证明了ABEEMσπ/MM模型的参数具有很好的可转移性.     

Alumina-supported iron catalyst - long-term study of the light product distribution in Fischer-Tropsch synthesis

Several models have been proposed to describe the carbon number product distribution and mechanism in Fischer-Tropsch synthesis (FTS). However, these models have not fully explained the product distribution and mechanism owing to the wide range and complexity of hydrocarbons in FTS. This study was conducted based on the Yao and Anderson-Schulz-Flory (ASF) carbon number product distribution models for light (C₁–C₆) hydrocarbon products of a Fe/Al₂O₃ catalyst. The product distribution based on the molar ratio of olefin to paraffin (O/P) and the neighboring olefins was also studied in order to better understand the mechanism in FTS and C₂ olefin deviation during FTS.Two sets of experiments (A and B) with different reaction conditions were conducted in microtubular fixed-bed reactors on the Fe/Al₂O₃ catalyst for 2249 h and 360 h, respectively. We found that the α values from the Yao and ASF carbon number product distribution models are relatively similar. The α values from the Yao carbon number product distribution plots are relatively constant, irrespective of the reaction conditions.Interestingly, it was also found that the secondary reactions of the C₂ olefin by re-adsorption to produce paraffins and long-chain olefins are dependent on the CO conversion and the reaction temperature during the FTS. Also, the product distribution of the neighboring olefins during the reduction condition gave a similar trend to what was observed for other reaction conditions. This result confirmed what was observed in the Yao and ASF carbon number product distribution of the olefins.     

Macroscopic solvation of molecules in excited states: An MCSCF model including solvent polarization effects-I

An MCSCF model including the effects of solvent polarization is developed. The model is applied within the limitations of INDO approximations to look into the dominant effects of solvent polarization on the electronic structure in the excited states of a model system (e.g.nπ ^* states of H₂CO). Important features of macroscopic solvation-induced reorganization of electron density and some consequence thereof are noted.     

A study of the influence of impurities when discriminating between the terminal and penultimate copolymerization models

We investigate the effect of the presence of reactive impurities during copolymerization on mathematical model discrimination and determine if impurities adversely affect the model discrimination process, depending on the criteria and response(s) used. Nonlinear mechanistic process models are employed in a multiresponse statistical discrimination scenario. The investigation also considers, via several case studies, the extremely interesting question of whether impurity effects could lead to incorrect choices about the copolymerization mechanism itself (terminal versus penultimate model). The answers to the above questions are that reactive impurities do affect the model discrimination process and, what is more, they may mislead researchers as to the choice of the copolymerization mechanism. The latter unexpected result is further illustrated with experimental data from the literature. © 2000 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 38: 2319–2332, 2000     

量子化学模型方法在机械力化学中的应用

机械力化学作为一种无需溶剂的绿色化学技术得到广泛关注。然而,机械力化学反应机制需要从原子和分子尺度上深入理解力诱导的化学反应。在过去的20年中,量子化学模型方法在机械力化学机理研究中得到广泛应用,高精度量化计算可得到外力下变形分子的几何结构、能量、过渡态等诸多性质。本文介绍了目前机械力化学领域的主流量子化学模型的基本原理,同时也关注了这些模型方法在软件上的具体实现,并借助典型的案例阐述了量子化学模型在解释机械力化学机理中的作用与价值。     

非线性和时变参数时间序列模型及其在水质分析中的应用

时间序列分析是根据观测值建立数学模型，研究数据的内在规律，现有文献主要是介绍线性时间序列模型。研究表明化学分析数据很多涉及非线性时间序列，且具有时变参数特性。本文在研究线性、非线性时间序列模型基础上，提出一种具有时变参数特性的非线性时间序列模型。该模型用于加酸调ｐＨ的循环冷却水系统进行预报，可使ｐＨ值极差降低３～６倍，对保证水质稳定具有重要意义。     

蛋白质离子交换平衡模型

本文综述了蛋白质的离子交换平衡模型的研究进展，包括Langmuir模型、质量作用模型以及建立在双电层理论基础上平行板模型和表面过剩模型，分析比较了各模型的特点和存在的问题，对今后的研究方向进行讨论。     

Nonlinear Dependence of Fish Bioconcentration on n-Octanol/Water Partition Coefficient

Abstract

The log-log relationship between the bioconcentration tendency of organic chemicals in fish and the n?octanol/water partition coefficients breaks down for very hydrophobic compounds. The use of parabolic and bilinear models allows this problem to be overcome. The QSAR equation log BCF = 0.910 log P - 1.975 log (6.8 10^?7 P + 1) - 0.786 (n = 154; r = 0.950; s = 0.347; F = 463.51) was found to be a good predictor of bioconcentration in fish.     

Measurement and correlation of the solubility of genistin in eleven organic solvents from T = (283.2 to 323.2) K

The solubilities of genistin in pure solvents including tetrahydrofuran, acetone, ethyl acetate, acetonitrile, isopropanol, n-butyl alcohol, methanol, cyclohexane, n-hexane, chloroform and ethanol were determined by the high performance liquid chromatography (HPLC) analysis method at T = (283.2, 293.2, 303.2, 313.2 and 323.2) K. The fusion enthalpy of genistin was estimated by the group contribution method. The solubility data of genistin were correlated by the simplified thermodynamic model, the modified Apelblat model, λh model and NRTL (Non-Random Two-Liquid) model. The calculated values by all models were in good agreement with the experimental values and however, the NRTL model could give better correlation results than other three models.     

Simplistic modeling approach to heterogeneous dilute-acid hydrolysis of cellulose microcrystallites

The classic kinetic model for cellulose hydrolysis is often referred to as pseudo-homogeneous, a term revealing the insight that the process is actually heterogeneous. During the past 10–15 yr, the shortcomings of this model have been demonstrated in various studies and the interest in the heterogeneous aspects has increased. The present work presents a simplistic model in which the intrinsic, heterogeneous hydrolysis and transport rates are coupled by the assumption of a constant glucosidic surface concentration. The mechanisms affecting these two rates are largely unknown, but the model serves as a guideline for further exploration of the process.