DNA小沟结合药物DB818的分子动力学模拟和结合自由能分析

采用分子动力学模拟了DNA小沟结合芳香二脒药物DB818形成的复合物. 通过5 ns的模拟研究表明: DB818药物分子可紧密结合在DNA的AATTC小沟区域, 和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物. 由于噻吩硫原子的弱电负性, 使DB818能够以更大的伸展程度与DNA的小沟结合, 形成更强的结合力. DB818苯并咪唑的氮原子能够与DNA 7位和19位T碱基上的氧原子形成两个稳定的氢键, 同时, DB818末端氨基氮原子分别与DNA 的20位T碱基的氧原子和9位C碱基的氧原子形成两个氢键. 另外, 运用MM_PBSA方法计算了DB293-DNA和DB818-DNA复合物的结合自由能, 计算结合能与实验值能较好的吻合, 通过比较其结合自由能, 从热力学能量角度说明了DB818有较大的熵值与较小的焓值贡献, 从而与DNA小沟结合的结合力比DB293强. 本文在分子水平上提供了DB818直接与双螺旋DNA相互作用的结构及复合物的动态变化情况, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

DNA小沟结合二脒:水分子的识别作用

采用分子动力学模拟了DB921-DNA复合物, 通过7 ns的模拟研究表明: DB921一端的氨基氮原子与一个水分子形成氢键, 同时, 水分子又与DNA的5位A碱基的氮原子形成一个氢键. 水分子在DB921与DNA小沟结合中起了桥连的作用, 使得直线型的芳香二脒化合物DB921通过水桥与DNA小沟结合, 水分子诱导DB921分子与DNA的小沟域构型相适应, 与DNA小沟域的AATTC碱基有较强的结合作用. 在分子水平上提供了DB921与双螺旋DNA相互作用的结构及复合物的动态变化情况, 指出水分子在DNA小沟结合二脒化合物中的识别作用, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.     

Fapy-G对碱基氢键复合物影响的理论研究

应用量子化学方法对2,6-二氨基-4-羟基-5-甲酰胺嘧啶(Fapy-G)与正常碱基作用形成的20种氧化碱基对的多种性质进行了理论分析,碱基G的C8位被氧化后N7和N9位分别增加一个H原子,使其由氢键受体变成氢键供体,N7,N9及O6原子所带的电荷变负,同时O6作为氢键受体的能力增强.与碱基单体相比,碱基对中形成氢键的受体原子所带的电荷平均减小0.05 e;供体H原子所带的电荷平均增大0.02 e.Fapy-G分子中六元环上受体N原子参与形成氢键时,环的呼吸振动模式和N与对位C的振动模式的振动频率蓝移;与氢键相关的振动频率红移.所有氧化碱基对中,NH…N比NH…O氢键作用更强,且在NH…N氢键中,在六元环上的供体N原子形成的氢键比在氨基或开环上的供体N原子形成的氢键强.Fapy-G与碱基A作用结合能区域顺序为1>2>4>3,与碱基T(R)作用区域顺序为3=4>1>2;水溶液使Fapy-G与碱基C作用的结合能减弱程度最大,结合能达到41.84~58.58 k J/mol,且使碱基对结合能力次序发生改变.     

N2O…HOCl复合物结构和性质的理论研究

分别在DFT-B3LYP和MP2/6-311++G**水平上求得HOCl + N2O复合物势能面上的六种(S1, S2, S3, S4, S5和S6)和四种(S1, S2, S4和S5)构型. 频率分析表明,其中的S1和S3为过渡态,其它为稳定构型. 在复合物S3, S5 和S6中,HOCl 单体的σ*(5O-6H)作为质子供体,与N2O单体中作为质子受体的3O原子相互作用,形成氢键结构,而在氢键复合物S2中, 质子受体为N2O单体中的端1N原子;复合物S1中,HOCl分子的σ*(5O-4Cl)作为质子供体与N2O分子中的端1N原子(质子受体)相互作用,形成卤键结构,而复合物S4中的卤键结构的质子受体为N2O分子中的端3O原子. 经B3LYP/6-311++G**水平上的计算,考虑了基组重叠误差(BSSE)校正的单体间相互作用能在-1.56 ~ -8.73 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

单核Zn(Ⅱ)/Ni(Ⅱ)含氮配体配合物的合成、晶体结构及性质

合成了2个新的配合物[Zn(BPP)2(H2O)4](2,6-NDS)·0.5H2O(1)和[Ni(phen)2(H2O)2](A-2,5-DSA)·3H2O(2)(2,6-NDS=2,6-萘二磺酸根,A-2,5-DSA=苯氨-2,5-二磺酸根,BPP=1,3-二(4-吡啶基)丙烷,phen=1,10-邻菲咯啉),用X-射线单晶衍射结构分析方法测定了配合物的晶体结构。配合物1是单核分子,Zn2+离子与2个1,3-二(4-吡啶基)丙烷的2个N原子及4个水分子配位,形成单核配位阳离子。相邻配位阳离子通过配位水分子与氮原子的氢键作用联接成一维双螺旋阳离子链。双螺旋阳离子链与未配位的2,6-萘二磺酸根阴离子通过氢键作用形成二维超分子网。配合物2是单核分子,Ni2+离子与2个1,10-邻菲咯啉分子中的4个N原子及2个水分子配位,形成单核配位阳离子。配位阳离子与游离的水分子及苯氨-2,5-二磺酸根阴离子通过氢键作用构筑成二维超分子网。     

双铂核药物与DNA作用的理论研究

利用分子力学和量子化学方法研究了双铂核药物［{trans-PtCl(NH₃)₂}₂(μ-NH₂(CH₂)_nNH₂)］²⁺与寡聚DNA片段d(ATATG*TACATAT)·d(ATGTG*TACATAT)复合物的几何构型和电子结构. 计算结果表明,Pt配合物与DNA中碱基G的N7原子形成了较强的配位键,并与O6原子之间存在较强的静电作用,使药物与DNA产生稳定作用,药物中的烃链的伸缩性使得DNA在键合药物后其构型并未发生大的变化. 同时,铂配合物中配体NH₃上的H与其邻近的鸟嘌呤的O6,DNA中磷酸根上的O以及与其邻近的碱基T上的O或N等电负性较大的原子间形成的氢键及弱氢键也是影响Pt配合物与DNA键合及其几何结构变化的重要因素. 这些化学键和氢键是药物分子能够对DNA进行识别的重要基础. 因此,可以认为药物结合后所引起DNA的变形较小可能是药物与顺铂产生不同的抗癌机理的主要原因.     

胞嘧啶与一氧化碳复合物的结构与性质

在B3LYP/6-311+G**水平上对胞嘧啶…CO复合物体系进行了理论计算, 发现了6个能量极小的复合物. 其结合方式是CO的C或O原子与胞嘧啶的N—H键形成氢键, 最稳定的复合物的结合能为-8.72 kJ·mol-1. CO的C原子与胞嘧啶的结合具有更强的优势, C原子结合的复合物中CO的键长缩短, 而O结合的复合物中CO键长伸长. 同时, 对复合物的振动分析发现, 在C原子结合的复合物中CO的伸缩频率蓝移, 而O结合的复合物中CO伸缩频率是红移的.     

Schiff碱N－氧化吡啶－2－甲醛缩氨基脲及其铜配合物的合成、晶体结构及谱学性质的研究

合成了Schiff碱N 氧化吡啶 2 甲醛缩氨基脲 (PNOS)及其配合物 [Cu(PNOS) (NO3 ) 2 ],并用单晶X射线衍射法测定了配体和配合物结构 .PNOS晶体中通过传统氢键形成双层二维网状结构 ,再由非传统氢键自组装成三维网状结构 .配合物 [Cu(PNOS) (NO3 ) 2 ]中的铜为六配位 ,畸变八面体结构 ,Schiff碱 (PNOS)通过N 氧化吡啶N—O的O原子 ,亚胺基CN的N原子 ,及羰基CO的O原子与铜配位 ;一个硝基以单齿配体形式与铜配位 ,另一个则以双齿配体形式配位 .配合物分子通过经典氢键相互作用 ,形成单层二维网状结构 ,再通过非经典氢键作用 ,自组装成双层二维网状结构     

O原子与HNCO反应机理的量子化学及电子密度拓扑研究

采用MP2(Full),B3LYP,QCISD/MP2和CCSD(T)/MP2方法在6-311G(d,p)水平上对O原子与HNCO反应的微观机理进行了理论研究．采用MP2(Full)和B3LYP对反应位能面上的各驻点进行几何构型的全优化,振动分析和IRC计算证实了中间体和过渡态的真实性和相互连接关系．四种方法计算得到了四个反应通道的反应活化能．研究表明,O原子进攻HNCO中的H原子为反应的主要反应通道,该通道中形成了两个分子复合物,电子密度拓扑分析表明这两个分子复合物均为氢键复合物．     

Molecular dynamics simulation of thymine glycol-lesioned DNA reveals specific hydration at the lesion

One nanosecond molecular dynamics (MD) simulation of a thymine glycol (TG)-lesioned part of human lymphoblast AG9387 was performed to determine structural changes in DNA molecule caused by the presence of a lesion. These changes can be significant for proper recognition of lesions by a repair enzyme. Thymine glycol is the DNA oxidative lesion formed by addition of OH radicals to C5 and C6 atoms of the thymine base. This lesion is known as causing Cockayne Syndrome-inherited genetic disorder. Distribution of water molecules in a hydration shell around the DNA molecule was analyzed for its contribution to the recognition of the TG lesion by the repair enzyme. The results of MD simulation show there is a specific DNA structural configuration formed at the lesion. After 500 ps the DNA is bent in a kink at the TG site. This change dislocates the glycosyl bond at C5' to a position closer to the DNA surface, and thus its atoms are more exposed to the surrounding water shell. The increased number of water molecules that are close to the TG site indicates that the glycosyl bond may be easily contacted by the repair enzyme. In addition, the higher number of water molecules at the TG site substantiates the importance of water-mediated hydrogen bonds created between the repair enzyme and the DNA upon formation of the complex. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1723-1731, 2001     

Structure and dynamics of poly(T) single-strand DNA: implications toward CPD formation

The formation of cyclobutane pyrimidine dimers between adjacent thymines by UV radiation is thought to be the first event in a cascade leading to skin cancer. Recent studies showed that thymine dimers are fully formed within 1 ps of UV irradiation, suggesting that the conformation at the moment of excitation is the determining factor in whether a given base pair dimerizes. MD simulations on the 50 ns time scale are used to study the populations of reactive conformers that exist at any given time in T18 single-strand DNA. Trajectory analysis shows that only a small percentage of the conformations fulfill distance and dihedral requirements for thymine dimerization, in line with the experimentally observed quantum yield of 3%. Plots of the pairwise interactions in the structures predict hot spots of DNA damage where dimerization in the ssT18 is predicted to be most favored. The importance of hairpin formation by intra-strand base pairing for distinguishing reactive and unreactive base pairs is discussed in detail. The data presented thus explain the structural origin of the results from the ultrafast studies of thymine dimer formation.     

Independent generation of 5-(2'-deoxycytidinyl)methyl radical and the formation of a novel cross-link lesion between 5-methylcytosine and guanine

Reactive oxygen species (ROS) can damage DNA. Although a number of single nucleobase lesions induced by ROS have been structurally characterized, only a few intrastrand cross-link lesions have been identified and characterized, and all of them involve adjacent thymine and guanine or adenine. In mammalian cells, the cytosines at CpG sites are methylated. On the basis of the similar reactivity of 5-methylcytosine and thymine toward hydroxyl radical and the similar orientation of adjacent thymine guanine (TG) and 5-methylcytosine guanine (mCG) in B-DNA, we predict that the cross-link lesion, which was identified in TG and has a covalent bond formed between the 5-methyl carbon atom of T and the C8 carbon atom of G, should also form at mCG site. Here, we report for the first time the independent generation of 5-(2'-deoxycytidinyl)methyl radical, and our results demonstrate that this radical can give rise to the predicted novel intrastrand cross-link lesion in dinucleoside monophosphates d(mCG) and d(GmC). Furthermore, we show that the cross-link lesion can also form in d(mCG) from gamma irradiation under anaerobic conditions.     

Effect of branching architecture on glass transition behavior of hyperbranched copolystyrenes: the experiment and simulation studies

By controlling the feed ratio of CMS/styrene and the polymerization time, a series of hyperbranched copolystyrenes(HBCPS) were synthesized with comparable weight-averaged molecular weights(Mw) but different degree of branching(DB) through atom transfer radical self-condensing vinyl copolymerization(ATR-SCVCP) with Cu Br/2,2?-bipyridyl as the catalyst. The resulting HBCPS samples were used to investigate the effect of branching architecture on their glass transition behavior. With the DB increased, the glass transition temperatures(Tg) of HBCPS samples measured by DMA and DSC both decreased. Their spin-lattice relaxation times(1H T1r) of protons displayed the same downtrend with increasing DB. Besides, a correlation between the Tgs and the DB was well established by all-atom molecular dynamics(MD) simulations. The values of MD-determined Tgs are little higher than the corresponding experimental ones. However, the dependence of Tgs on DB is in good agreement with the experimental results, i.e., Tg decreases both in experiments and simulations with increasing DB.     

The excited states of adenine and thymine nucleoside and nucleotide in aqueous solution: a comparative study by time-dependent DFT calculations

The effect substitutions at nitrogen atom 1 of thymine and nitrogen atom 9 of adenine have on lowest energy excited electronic states has been studied by means of time-dependent PBE0 calculations in aqueous solution. In agreement with the experimental indications, the vertical excitation energy of the bright state of 1,methyl-thymine, thymine nucleoside and thymine nucleotide is red-shifted with respect to that of thymine. Deoxyribose and deoxyribose-phosphate substituents affect mainly the lowest energy dark state of adenine and thymine, slightly increasing their oscillator strength. The excited states of 9, methyl-adenine and 1, methyl-thymine have also been studied by using the recently developed M052X, CAM-B3LYP and LC-ωPBE density functionals. The computed VEE are in good agreement with those obtained by using PBE0, which, however, provides values closer to the experimental band maximum.     

The N6-methyl group of adenine further increases the BI stability of DNA compared to C5-methyl groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.     

Hydrogen-abstracted adenine-thymine radicals with interesting transferable properties

The formation of radicals on DNA bases through various pathways can lead to harmful structural alterations. Such processes are of interest for preventing alteration of healthy DNA and, conversely, to develop more refined methods for inhibiting the replication of unwanted mutagenic DNA. In the present work, we explore theoretically the energetic and structural properties of the nine possible neutral radicals formed via hydrogen abstraction from the adenine-thymine base pair. The lowest energy radical is formed by loss of a hydrogen atom from the methyl group of thymine. The next lowest energy radicals, lying 8 and 9 kcal mol-1 higher than the global minimum, are those in which hydrogens are removed from the two nitrogens that would join the base pair to 2-deoxyribose in double-stranded DNA. The other six radicals lie between 16 and 32 kcal mol-1 higher in energy. Unlike the guanine-cytosine base pair, adenine-thymine (A-T) exhibits only minor structural changes upon hydrogen abstraction, with all A-T derived radicals maintaining planarity. Moreover, the energetic ordering for the radicals of the two isolated bases (adenine and thymine) is preserved upon formation of the base pair, though with a wider spread of energies. Even more significantly, the energetic interleaving of the (A-H)*-T and A-(T-H)* radicals is correctly predicted from the X-H bond dissociation energies of the isolated adenine and thymine. This suggests that the addition of the hydrogen-bonded complement base only marginally affects the bond energies.     

Quantum-chemical study on uracil and thymine nitrosonium complexes

Quantum-chemical calculations at the RI-MP2/L1 level of theory showed that the most energetically favorable complexes of uracil and thymine with nitrosonium cation are those of n-type with NO⁺ coordination at the nitrogen or oxygen atom. A correlation was found between the experimental and calculated affinities of the dioxo tautomer of thymine for nitrosonium ion ( $ A_{NO^ + } $ ). A linear relation was revealed between $ A_{NO^ + } $ values for structurally similar tautomers of uracil and thymine.