带电组氨酸侧链与DNA碱基间非键作用强度的理论研究

采用MP2方法和6-31+G(d,p)基组优化得到了带有一个正电荷的组氨酸侧链与4个DNA碱基间形成的18个氢键复合物的气相稳定结构, 从文献中获取了组氨酸侧链与DNA碱基间形成的12个堆积和T型复合物的气相稳定结构, 使用包含基组重叠误差(BSSE)校正的MP2方法和aug-cc-pVTZ基组及密度泛函理论M06-2X-D3方法和aug-cc-pVDZ基组计算了这些复合物的结合能. 研究结果表明, 包含BSSE校正的M06-2X-D3方法和aug-cc-pVDZ基组能够给出较准确的结合能; 气相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用明显强于堆积作用和T型作用, 组氨酸侧链最易通过离子氢键与胞嘧啶C和鸟嘌呤G作用形成氢键复合物, 组氨酸与胞嘧啶C和鸟嘌呤G间的T型作用强于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用; 水相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用仍明显强于堆积作用和T型作用, 组氨酸侧链更易与胞嘧啶C和鸟嘌呤G相互作用形成氢键复合物, 但是最强的组氨酸侧链与胞嘧啶C间的T型作用明显弱于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用, 说明水相条件下组氨酸侧链与DNA碱基间主要通过离子氢键作用形成氢键复合物.     

精氨酸侧链和核酸碱基间离子氢键作用强度分析

采用MP2/6-31+G(d,p)方法优化得到了22个由精氨酸侧链与碱基尿嘧啶、 胸腺嘧啶、 胞嘧啶、 鸟嘌呤及腺嘌呤形成的氢键复合物的气相稳定结构, 使用包含BSSE校正的MP2/aug-cc-pVTZ方法计算得到了复合物的气相结合能, 通过MP2/6-31+G(d,p)方法和PCM模型优化得到了复合物的水相稳定结构, 采用MP2/aug-cc-pVTZ方法和PCM模型计算得到了复合物的水相结合能. 研究发现, 精氨酸侧链与碱基间的离子氢键作用强度与单体间电荷转移量、 氢键临界点电子密度及二阶作用稳定化能密切相关. 与中性氢键相比, 离子氢键作用具有更显著的共价作用成分. 研究还发现, 精氨酸侧链和碱基间形成的氢键复合物的稳定性次序可以通过氢键受体碱基分子上氧原子和氮原子的质子化反应焓变进行预测, 质子化反应焓变越负, 形成的氢键复合物越稳定.     

槲皮素与尿嘧啶和胸腺嘧啶氢键的结合位点

深入理解药物分子和核酸碱基间的相互作用机制对合理设计研发新型高效药物有重要意义. 本文运用密度泛函理论B3LYP方法对核酸碱基尿嘧啶和胸腺嘧啶与药物分子槲皮素间的氢键相互作用位点进行了研究. 使用B3LYP/6-31G(d)方法优化得到了30个稳定的氢键复合物结构, 使用B3LYP/6-311++G(3df,2p)方法计算了这些复合物的结合能. 研究结果表明, 槲皮素可以使用5个不同的结合位点与尿嘧啶或胸腺嘧啶形成氢键复合物, 尿嘧啶或胸腺嘧啶可以使用3个不同的结合位点与槲皮素形成氢键复合物. 当槲皮素的结合位点固定时, 槲皮素与尿嘧啶的位点u1或胸腺嘧啶的位点t1形成的氢键作用最强, 与位点u2或位点t2形成的氢键强度最弱; 当尿嘧啶或胸腺嘧啶的作用位点固定时, 二者与槲皮素的位点qu1 形成的氢键作用最强, 与位点qu5 作用强度次之, 与位点qu3的作用强度最弱. 分子中原子(AIM)和自然键轨道(NBO)分析计算结果表明, 轨道作用在氢键中起重要作用.     

核酸碱基与甘氨酸二肽间氢键作用的最佳位点

优化得到了碱基腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶与甘氨酸二肽分子形成的28个氢键复合物的稳定结构并计算了结合能, 探讨了五种碱基与甘氨酸二肽分子间氢键作用的最佳位点. 本文研究发现: 每种碱基均可以通过不同位点与二肽分子形成氢键复合物, 腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶分别最倾向使用A3、T1、U1、G3及C1位点与甘氨酸二肽分子形成氢键复合物; 碱基分子某位点的质子化反应焓变越负所形成的氢键复合物越稳定, 去质子化反应焓变越小所形成的氢键复合物越稳定; 由氢键复合物的结合能计算得到的稳定性次序与由碱基分子质子化和去质子化反应焓变推得的稳定性次序一致.     

在遗传字母序列中新添两个字母

在此以前遗传字母序列中仅有四个字母(或碱基),即腺嘌呤(A),鸟嘌呤(G)、胸腺嘧啶(u或DNA的T)和胞嘧啶(C)。这些碱基通过氨基和羰基间的氢键使双链螺旋型RNA中的两条互补链结合在一起。其中,A与u形成两个氢键,G与C形成三个氢键,从而组成两个碱基对。生物世界就是靠这两个“碱基对”来运载遗传信息。最近,瑞士生物化学家Steven A.Benner证实其它碱基对亦可能运载遗传信息。例如,将鸟嘌呤或胞嘧啶中的羰基和氨基互相异位,可得到异嘌呤或异胞嘧啶。含这些非标准碱基的核苷酸能通过三个氢键作用形成所谓的Watson-Crick碱基对。     

碱基腺嘌呤与多肽酰胺间氢键作用的最佳位点

正确理解核酸碱基和蛋白质多肽间的作用机制有助于人们利用这些生物分子有效地进行分子设计,进而制备具有特殊纳米结构和功能的生物分子材料.本文优化得到了碱基腺嘌呤与N-甲基乙酰胺、甘氨酸二肽、丙氨酸二肽形成的20个氢键复合物的结构并计算了结合能,探讨了腺嘌呤与多肽酰胺间氢键作用的最佳位点.研究发现:腺嘌呤可以使用两个不同位点(A1位点和A2位点)与N-甲基乙酰胺形成N―H…N型或者N―H…O=C型氢键复合物,腺嘌呤使用A1位点与N-甲基乙酰胺形成的N―H…N型氢键复合物更稳定;二肽分子可以使用主链上两个不同位点(丙氨酸的Ala7位点和Ala5位点或者甘氨酸的Gly7位点和Gly5位点)与腺嘌呤形成含有N―H…N和N―H…O=C两条氢键的复合物,二肽分子使用Ala7或Gly7位点与腺嘌呤形成的氢键复合物更稳定;腺嘌呤与多肽间的氢键作用强于其与N-甲基乙酰胺的作用.基于分子中的原子理论与自然键轨道计算结果分析了氢键作用的本质.     

水溶液中氨基酸侧链对G∶C碱基对间氢键影响的理论研究

利用量子化学方法研究了气相和水溶液下,氨基酸侧链与鸟嘌呤和胞嘧啶间的氢键作用.应用B3LYP/6-31+G(d,p)方法优化复合物几何结构,使用MP2/aug-cc-p VDZ方法进行复合物能量、自然键轨道(NBO)电荷和二阶稳定化能的计算.结果表明,水溶液可使氨基酸侧链与碱基或碱基对之间氢键键能显著减小;带电复合物气相和水溶液氢键键能之差范围为50.63~146.48 k J/mol,中性为0.17~24.94 k J/mol;电荷的转移量与氢键键能成正比,电荷转移量越多,复合物越稳定;二阶稳定化能与氢键键长成反比,与电荷转移量成正比,且气相与水溶液氢键二阶稳定化能之比约为两相的电荷转移量之比.水溶液对该类体系中氢键作用具有明显影响.     

取代基对腺嘌呤与胸腺嘧啶氢键复合物体系结合能的影响

优化得到了17个取代胸腺嘧啶与腺嘌呤形成的氢键复合物的结构, 并计算了这些复合物的结合能, 探讨了胸腺嘧啶上不同取代基对结合能的影响. 结果表明, CF₃取代的胸腺嘧啶与腺嘌呤间的结合能大于胸腺嘧啶与腺嘌呤间的结合能, 这可能是屈氟尿苷具有阻止病毒及肿瘤扩散功能的原因所在. SO₃H, CN和NO₂取代的胸腺嘧啶与腺嘌呤间具有更大的结合能, 表明这3个基团取代的胸腺嘧啶也可能具有潜在的抗肿瘤作用. 分子中原子理论与自然键轨道分析表明, 在所有体系中, 氢键N—H…N最强, N—H…O=C次之, C—H…O=C最弱, 轨道作用在氢键作用中占有重要地位.     

氨基酸侧链与氧化鸟嘌呤碱基对体系的氢键作用

应用量子化学方法,分别在气相和水溶液中对氨基酸侧链与氧化鸟嘌呤碱基对(8-oxo-G∶C)形成的三体复合物的氢键键能、几何结构、电荷分布及二阶稳定化能进行了研究.结果表明,水溶液的存在削弱了复合物中的氢键强度,电荷分布变化明显,水溶液中形成氢键位点的电荷变化量约为气相中的10倍,而几何结构变化不明显、对于酶与DNA之间的相互作用的研究需在水溶液中进行.水溶液对带电三体复合物中8-oxo-G∶C与氨基酸侧链间的氢键有较大影响,键能平均减小了69.23 k J/mol,不带电复合物仅减小了3.60k J/mol.水溶液中三体复合物中8-oxo-G∶C间的氢键受侧链的影响不大,且与侧链带电与否无关,带电复合物和不带电复合物的氢键强度分别减小了24.57和30.05 k J/mol,且二阶稳定化能越大,其对应的氢键键长越短.     

O6-鸟嘌呤和O4-胸腺嘧啶甲基化与DNA碱基异常氢键作用的理论研究

用密度泛函B3LYP方法在6-311＋G**基组水平上对顺(cis-)、反式(anti-)O⁶-甲基鸟嘌呤(O⁶-MeG)和O⁴-甲基胸腺嘧啶(O⁴-MeT)与DNA碱基(腺嘌呤A、鸟嘌呤G)的非Watson-Crick氢键二聚体进行了优化. 在MP2/cc-pVXZ (X＝D,T)//B3LYP/6-311＋G**水平上, 采用完全基组外推方法校正了氢键二聚体的相互作用能, 并用完全均衡校正法(CP)校正了基组重叠误差(BSSE). 此外, 在B3LYP/6-311＋G**水平上计算了各氢键碱基对的全电子波函数, 并用分子中的原子理论(AIM)和电子密度拓扑方法分析了碱基间的弱相互作用. 计算结果显示, 甲基化使碱基对间的氢键作用模式发生了明显的扭转和不同程度的位移, 碱基间的电子密度分布和氢键作用能明显减小, 甲基化对O⁶-MeG和O⁴-MeT与DNA碱基间的氢键作用是去稳定化的, 这种影响主要来自于大体积的甲基的空间效应和给电子效应, 且对顺式的影响明显大于反式. 计算结果与文献给出的实验结论基本一致.     

A method for molecular correlation energy calculations. Application to the determination of dissociation energies of diatomic and polyatomic molecules

A simple and economical method for molecular correlation energy calculations is developed. In this method, the internal part of the correlation energy is calculated by means of a CI in a minimal basis set and the non-internal part (semi-internal and all-external) is evaluated using an original atoms-in-molecule method. It is successfully applied to the determination of dissociation energies of some diatomic (H₂, NH, C₂, CN, N₂, CO, NO, O₂, F₂) and polyatomic (H₂O, N₂O, CO₂, N₃H, CH₂N₂, CH₂CO, C₂N₂) molecules. The results are compared to those obtained using very elaborate variational methods.Aspirant du Fonds National Belge de la Recherche Scientifique.     

Partitioning of atomization energy

The present work provides a technique for partitioning the atomization energy of a molecule into diatomic contributions. The method is largely based on the redistribution of the kinetic energy term in Mayer's energy partitioning and uses free‐atom energies as a reference. The comparison of Mayer's original method, the alternative Ichikawa–Yoshida approach, and the new atomization energy partitioning (AEP) shows that the new approach has advantages in describing trends in diatomic energies in molecules with triple bonds, as well as for hydrogen bonds. The proposed AEP is a viable alternative to Mayer's energy partitioning method. © 2007 Wiley Periodicals, Inc. Int. J. Quantum Chem, 2008     

Interaction energy contours of molecules as probed by a test atom

Potential energy contour curves for some simple molecules interacting with a test atom (a Li atom in the ground state) are calculated by theab initio method. Relative hardness of the repulsive molecular surface is found to be anisotropic. This indicates that the outer boundary surface of a molecule changes its form depending on the extent of applied external forces. The attractive nature of the molecular surface is demonstrated to be local. Surroundings of electronegative atoms having nonbonded electron pairs are found to be attractieve to the test atom having both electron-donating and electron-accepting abilities.     

A representation of the polarization term in the interaction energy between a molecule and a point-like charge

The interaction energy of a molecule M with a point-like charge q can be partitioned into simpler contributions, two of which can be expressed in terms of the charge distribution _M of the sole M. The first term, qV(r), represents the interaction of q with the undistorted charge _M ⁰ of M while the second q ² P(r) gives the additional contributions due to the polarization of _M ⁰ under the influence of the charge q placed at the point r. In this paper we investigate the possibility of getting an inexpensive and sufficiently accurate analytical representation of P(r) over the whole space outside the van der Waals volume of M.     

Synthesis of Imines via Reactions of Benzyl Alcohol with Amines Using Half‐Sandwich (η6‐p‐cymene) Ruthenium(II) Complexes Stabilised by 2‐aminofluorene Derivatives

A new class of half‐sandwich (η⁶‐p‐cymene) ruthenium(II) complexes supported by 2‐aminofluorene derivatives [Ru(η⁶‐p‐cymene)(Cl)(L)] ( L  = 2‐(((9H‐fluoren‐2‐yl)imino)methyl)phenol ( L ¹ ), 2‐(((9H‐fluoren‐2‐yl)imino)methyl)‐3‐methoxyphenol ( L ² ), 1‐(((9H‐fluoren‐2‐yl)imino)methyl)naphthalene‐2‐ol ( L ³ ) and N‐((1H‐pyrrol‐2‐yl)methylene)‐9H‐fluorene‐2‐amine ( L ⁴ )) were synthesized. All compounds were fully characterized by analytical and spectroscopic techniques (IR, UV–Vis, NMR) and also by mass spectrometry. The solid state molecular structures of the complexes [Ru(η⁶‐p‐cymene)(Cl)(L²)], [Ru(η⁶‐p‐cymene)(Cl)(L³)] and [Ru(η⁶‐p‐cymene)(Cl)(L⁴)] revealed that the 2‐aminofluorene and p‐cymene moieties coordinate to ruthenium(II) in a three‐legged piano‐stool geometry. The synthesized complexes were used as catalysts for the dehydrogenative coupling of benzyl alcohol with a range of amines (aliphatic, aromatic and heterocyclic). The reactions were carried out under thermal heating, ultrasound and microwave assistance, using solvent or solvent free conditions, and the catalytic performance was optimized regarding the solvent, the type of base, the catalyst loading and the temperature. Moderately high to very high isolated yields were obtained using [Ru(η⁶‐p‐cymene)(Cl)(L⁴)] at 1 mol%. In general, microwave irradiation produced better yields than the other two techniques irrespective of the nature of the substituents.     

能源转换演示装置的设计

针对高中化学教材中提到的新能源问题,探索设计了可将太阳能转化为电能,电能转化为化学能,化学能转化为电能,电能转化为热能的能源转换演示装置。通过该实验的演示,可以让学生直观地感受到几种能量之间的转化过程,可供教学参考和借鉴。     

Influence of multivalent ions on power production from mixing salt and fresh water with a reverse electrodialysis system

Reverse electrodialysis is a membrane-based technique for production of sustainable electricity from controlled mixing of a diluted electrolyte solution (e.g., river water) and a concentrated electrolyte solution (e.g., sea water). Reverse electrodialysis has been investigated with pure sodium chloride solutions. In practice, however, in most cases also other ions are present in both feed solutions. In the present paper, the effect of multivalent ions on the performance of a reverse electrodialysis stack was investigated. Results show that, besides a higher stack resistance in presence of multivalent ions, especially the presence of multivalent ions in the dilute solution has a lowering effect on the stack voltage. This can be explained by an observed transport of these ions from the diluted electrolyte solution to the concentrated electrolyte solution. In order to prevent or hamper this transport against the activity gradient, monovalent-selective membranes can be used. This shows indeed better results with respect to the stack voltage. Therefore, it would be beneficial to use monovalent-selective membranes in reverse electrodialysis, especially in the case of a relatively high content of multivalent ions in the dilute (i.e., in the first stages of the installation where the sodium chloride content in the dilute is still relatively low).     

Efficient calculation of two-dimensional adiabatic and free energy maps: Application to the isomerization of the C13C14 and C15N16 bonds in the retinal of bacteriorhodopsin

Accurate calculation of potential energy and free-energy profiles along reaction coordinates of biological processes such as enzymatic reactions or conformational changes is fundamental to the obtention of theoretical insight into protein function. We describe here the practical implementation of the Automatic Map Refinement Procedure (AMRP) and two-dimensional Weighted Histogram Analysis Method (WHAM) for efficient computation of adiabatic potential energy and free-energy maps, respectively. Methods for efficiently sampling configuration space with high-energy barriers and for removing hysteresis in the case of periodic reaction coordinates are presented. The application of these techniques to the isomerization of the C13C14 and C15N16 bonds in the retinal of bacteriorhodopsin is described. In dark-adapted bacteriorhodopsin (bR), the retinal moiety exists in two conformers, all-trans and (13,15)cis, with the latter making ≃67% of the population. This experimental free energy difference is reproduced here to within k_BT. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1644–1658, 1999     

Electron impact ionization of amines using crossed electron-molecular beams

Ionization energy functions of triethyl amine, tri-n-propyl amine, N,N-dimethyl-benzenamine and N,N-diethyl-benzenamine, and appearance energy functions of major products of dissociative ionization have been measured using a high sensitivity crossed electron beam-molecular beam apparatus incorporating a quadrupole mass spectrometer and an on-line laboratory computer. The data obtained is useful in understanding the role of dopants in gas laser plasmas.     

太阳能光解水制氢的研究进展

本文概述了利用光催化技术催化分解水制氢的反应机理和研究进展。结合作者的最近研究,重点描述了TiO2及过渡金属氧化物,层状金属氧化物以及某些能利用可见光的光催化材料的结构和光催化特性,阐述了核课题的意义和今后的研究方向。