Electrostatic versus nonelectrostatic effects in DNA sequence discrimination by divalent ions Mg2+ and Mn2+

Mg2+ and Mn2+ ions are critical to the functioning of phosphoryl transfer enzymes, such as restriction endonucleases. Although these ions play similar roles in the chemical steps, they govern substrate specificity via modulating sequence discrimination by up to a factor of 10(5) with Mg2+ and only up to a factor of 10 with Mn2+. To explain whether such diversity originates in fundamental differences in the electronic structures of the nucleobase-hydrated-metal ion complexes, structures and interaction energies were determined at the density functional (DFT) and second-order M?ller-Plesset (MP2) levels of theory. Although both metal ions favor identical binding sites, Mn2+ complexes exhibit greater distortions from the ideal octahedral geometry and larger variability than the corresponding Mg2+ systems. In inner-shell complexes, with direct contact between the metal and the nucleobase, Mg2+ is preferred over Mn2+ in the gas phase, due primarily to nonelectrostatic effects. The interaction energies of the two metal ions are more similar in the outer-shell complexes, likely due to reduced charge transfer between the hydrated metal ion and the base moieties. Inclusion of solvation effects can amplify the relative nucleobase preferences of Mg2+ and Mn2+, indicating that bulk hydration modulates the balance between electrostatic and nonelectrostatic terms. In most cases, the base substitutions in solution are facilitated more by Mn2+ than by Mg2+. Electrostatic properties of the environment were demonstrated to have a major influence on the nucleobase preferences of the two metal ions. Overall, quantum chemical calculations suggest that the contrasting selectivity of Mg2+ and Mn2+ cofactors toward nucleobases derives from the larger flexibility of the Mn2+ complexes accompanied by the excessive polarization and charge-transfer effects as well as less favorable solvation.     

金属离子对HIV-1整合酶与硫氮硫扎平抑制剂作用模式影响的分子模拟研究

HIV-1整合酶（IN）通过依赖金属离子的两步反应将病毒DNA整合入宿主细胞过程中。结合于HIV-1上的金属离子个数的变化直接影响整合酶与抑制剂之间的结合。本工作用同源模建方法搭建了每条单链核心区具有两个Mg2+ 的（2Mg-IN-Core）和具有一个Mg2+ 的HIV-1 IN二聚体核心区模型（1Mg-IN-Core）。分子对接分别得到它们与硫氮硫扎平类化合物能量较低的复合物结构,把对接结果进行了比较。研究发现：当整合酶中结合的Mg2+个数改变时,它与抑制剂的结合模式也会发生很大的变化;抑制剂能够特异的且稳定的与2Mg-IN-Core模型的活性位点结合;同时与ASP64和GLU152螯合的那个Mg2+离子对于硫氮硫扎平抑制剂与整合酶上的结合有很大的影响。2Mg-IN-Core模型与抑制剂的复合物平均结构进行了2000 ps的 分子动力学模拟,分析发现同时与ASP64及ASP116螯合的Mg2+与IN蛋白形成了四个稳定的螯合键;同时与ASP64及GLU152螯合的Mg2+可与IN结合、也可与抑制剂形成稳定的配位键,这个Mg2+对IN与硫氮硫扎平抑制剂之间的结合有较大影响。     

Influence of Mg2+ and Cd2+ on the interaction between sparfloxacin and calf thymus DNA

Mg(2+) and Cd(2+) have different binding capacity to sparfloxacin, and have different combination modes with calf thymus DNA. Selecting these two different metal ions, the influence of them on the binding constants between SPFX and calf thymus DNA, as well as the related mechanism have been studied by using absorption and fluorescence spectroscopy. The result shows that Cd(2+) has weak binding capacity to SPFX in the SPFX-Cd(2+) binary system, but can decrease the binding between SPFX and DNA obviously in SPFX-DNA-Cd(2+) ternary system. Mg(2+) has strong binding capacity to SPFX. It can increase the binding between SPFX and DNA at concentrations <0.01 mM, and decrease the binding between them at concentrations >0.01 mM. Referring to the different modes of Mg(2+) and Cd(2+) binding to DNA, the mechanism of the influence of metal ions on the binding between SPFX and DNA has been proposed. SPFX can directly bind to DNA by chelating DNA base sites. If a metal ion at certain concentration mainly binds to DNA bases, it can decrease the binding constants between SPFX and DNA through competing with SPFX. While if a metal ion at certain concentration mainly binds to phosphate groups of DNA, it can increase the binding constants by building a bridge between SPFX and DNA. The influence direction of metal ions on the binding between quinolone and DNA relays on their binding ratio of affinity for bases to phosphate groups on DNA. Our result supports Palumbo's conclusion that the binding between SPFX and the phosphate groups is the precondition for the combination between SPFX and DNA, which is stabilized through stacking interactions between the condensed rings of SPFX and DNA bases.     

Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study

A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.     

Active calcium treatment transport via coupling between the enzymatic and the ionophoric sites of Ca2+ + Mg2+-ATPase.

The 20K dalton fragment of Ca2+ + Mg2+-ATPase obtained from th tryptically digested sarcoplasmic reticulum has been further purified using Bio-Gel P-100. This removed low-molecular-weight UV-absorbing and positive Lowry-reacting contaminants. The ionophoric activity of the 20K fragment in both oxidized cholesterol and phosphatidylcholine:cholesterol membranes is unaltered by this further purification. The 20K selectivity sequence in phosphatidylcholine:cholesterol membrane is Ba2+ greater than Ca2+ greater than Sr2+ greater than Mn2+ Mg2+. Digestion of intact sarcoplasmic reticulum vesicles with trypsin, which results in the dissection of the hydrolytic site (30K) from the ionophoric site (20K), is shown to disrupt energy transduction between ATP hydrolysis and calcium transport. This further implicates the 20K dalton fragment as a calcium transport site. These data and previous evidence are discussed in terms of a proposed model for the ATPase molecular structure and the mechanisms of cation transport in sarcoplasmic reticulum.     

Factors governing the substitution of La3+ for Ca2+ and Mg2+ in metalloproteins: a DFT/CDM study

Trivalent lanthanide cations are extensively being used in biochemical experiments to probe various dication-binding sites in proteins; however, the factors governing the binding specificity of lanthanide cations for these binding sites remain unclear. Hence, we have performed systematic studies to evaluate the interactions between La3+ and model Ca2+ - and Mg2+ -binding sites using density functional theory combined with continuum dielectric methods. The calculations reveal the key factors and corresponding physical bases favoring the substitution of trivalent lanthanides for divalent Ca2+ and Mg2+ in holoproteins. Replacing Ca2+ or Mg2+ with La3+ is facilitated by (1) minimizing the solvent exposure and the flexibility of the metal-binding cavity, (2) freeing both carboxylate oxygen atoms of Asp/Glu side chains in the metal-binding site so that they could bind bidentately to La3+, (3) maximizing the number of metal-bound carboxylate groups in buried sites, but minimizing the number of metal-bound carboxylate groups in solvent-exposed sites, and (4) including an Asn/Gln side chain for sites lined with four Asp/Glu side chains. In proteins bound to both Mg2+ and Ca2+, La3+ would prefer to replace Ca2+, as compared to Mg2+. A second Mg2+-binding site with a net positive charge would hamper the Mg2+ --> La3+ exchange, as compared to the respective mononuclear site, although the La3+ substitution of the first native metal is more favorable than the second one. The findings of this work are in accord with available experimental data.     

Fragmentation of doubly-protonated Pro-His-Xaa tripeptides: Formation of b 2 2+ ions

When ionized by electrospray from acidic solutions, the tripeptides Pro-His-Xaa (Xaa=Gly, Ala, Leu) form abundant doubly-protonated ions, [M+2H]²⁺. Collision-induced dissociation (CID) of these doubly-protonated species results, in part, in formation of b ₂ ²⁺ ions, which fragment further by loss of CO to form a ₂ ²⁺ ions; the latter fragment by loss of CO to form the Pro and His iminium [immonium is commonly used in peptide MS work] ions. Although larger doubly-charged b ions are known, this represents the first detailed study of b ₂ ²⁺ ions in CID of small doubly protonated peptides. The most abundant CID products of the studied doubly-protonated peptides arise mainly in charge separation involving two primary fragmentation channels, formation of the b ₂ /y ₁ pair and formation of the a ₁ /y ₂ pair. Combined molecular dynamics and density functional theory calculations are used to gain insight into the structures and fragmentation pathways of doubly-protonated Pro-His-Gly including the energetics of potential protonation sites, backbone cleavages, post-cleavage charge-separation reactions and the isomeric structures of b ₂ ²⁺ ions. Three possible structures are considered for the b ₂ ²⁺ ions: the oxazolone, diketopiperazine, and fused ring isomers. The last is formed by cleavage of the His-Gly amide bond on a pathway that is initiated by nucleophilic attack of one of the His side-chain imidazole nitrogens. Our calculations indicate the b ₂ ²⁺ ion population is dominated by the oxazolone and/or fused ring isomers.     

Ag+与Hg2+对G-四链体DNA结构的破坏及其对构筑DNA逻辑门的应用

本工作利用圆二色光谱研究了Ag⁺与Hg²⁺对4种代表性G-四链体DNA结构的破坏作用。结果表明Ag⁺可能通过与碱基G螯合从而破坏G-四链体结构;Hg²⁺能通过形成T-Hg²⁺-T碱基对,及其他方式破坏G-四链体结构。含巯基(-SH)的半胱氨酸与Ag⁺与Hg²⁺可以发生较强的配位作用,从而使被Ag⁺与Hg²⁺破坏后的G-四链体DNA结构得以回复。基于此,一个新颖的Ag⁺/Hg²⁺-半胱氨酸-DNA逻辑门得以构筑。     

Influence of Mg(2+) on the guanine-cytosine tautomeric equilibrium: simulations of the induced intermolecular proton transfer

Metallic ions are essential for stabilizing the nucleic acid structure, and are also involved in the majority of RNA and DNA biological functions. However, at large concentrations metals may play an opposite role by promoting alterations in the genetic code (mutagenicity). To contribute to the understanding of this effect, theoretical tools are used to investigate the influence of the magnesium dication on the guanine-cytosine (GC) base pair structure and stability. To this end, a fully hydrated Mg(2+) cation is inserted in two models: an isolated GC base pair, and a more realistic DNA model corresponding to a hydrated double-stranded trimer. Calculations performed with a hybrid ONIOM approach reveal that the Mg(2+) cation coordination to the GC base pair alters drastically the natural tautomeric equilibria in DNA by promoting single proton transfer. Nevertheless, the generated rare tautomer will have a limited impact on the total spontaneous mutation due to the low back-reaction barrier allowing a quick return to the canonical form. Additionally, it is demonstrated that the major effects of biological environment arise from the hydration and stacking influence, whereas the impact of phosphate groups is minor.     

复合物Ca+-RNA碱基的结构与稳定性

Gas-phase metal ion affinities and optimized structures of RNA nucleic acid bases for the Ca+ were determined at a density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of RNA nucleic acid and Ca+ on the different binding sites, were considered. For Cytosine, the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base tautomers. As to thymine, the bond energy of the ion with the most stable tautomer of the free nucleic acid base is the weakest among the three tautomer’s complexes, and that of the ion with least stable tautomer of the free nucleic acid base is the strongest . Uracil is similar to thymine. The two kinds of relation, bond energy and total energy for the complex, are in disagreement, as the metal affinities of RNA bases for the Ca+ depend on binding sites, and total energy of complex (Ca+-RNA base) relies on all atoms and their relative positions in the complex.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

土壤颗粒表面电场作用下固-液界面Mg2+-K+与Ca2+-K+交换动力学的比较研究

通过恒流法研究了不同表面电场作用下Mg2+、Ca2+吸附动力学. 结果发现: (1)实验初期阶段是强静电力作用下的零级动力学过程和一定反应时间后的弱静电力作用下的一级动力学过程, 且零级速率过程和一级速率过程之间存在明显的转折点; (2)不同电解质构成中Ca2+的吸附速率明显快于Mg2+的, 平衡吸附量也大于Mg2+的, 且Ca2+在土壤颗粒表面的覆盖度比Mg2+在土壤颗粒表面的覆盖度高; (3)离子的相对有效电荷系数与土壤颗粒表面电场作用的不同是各体系中Ca2+、Mg2+吸附动力学有差别的根本原因; (4)根据离子吸附的理论模型可以分别计算出速率系数、平衡吸附量、离子在土壤颗粒表面的覆盖度以及固定液的体积, 这些参数可以定量评估土壤颗粒表面电场对离子吸附动力学的影响.     

Effects of Carboxylic Acids on Calcite Formation in the Presence of Mg2+ Ions

The effects of seven carboxylic acids on calcite formation in the presence of Mg2+ ions, whose molar concentration ratio Mg2+/Ca2+ = 0.5 exclusively induced aragonite precipitation in the absence of carboxylic acids, were studied using a double diffusion technique. The presence of carboxylic acids, acrylic acid, maleic acid, tartaric acid, malonic acid, malic acid, succinic acid, and citric acid in the gel medium favored the formation of magnesian calcite relative to the amount of the additives. Induction time and the positions of the first precipitation were measured to analyze the behavior of crystallization based on the equivalency rule. The formation of magnesian calcite was also studied with the help of Avrami's equation (solid-state model for transformation). The results of applying this equation suggested that aragonite transformed into calcite through a solid-to-solid process. The formation of magnesian calcite was interpreted as the following process: aragonite nuclei, formed owing to Mg2+ ions at the initial stage of CaCO3 crystallization, transformed into calcite nuclei through a solid-to-solid process while their growth was inhibited by the adsorption of carboxylic acids. The magnesian calcite crystals grew on crystal seeds of calcite formed from aragonite nuclei. Copyright 1999 Academic Press.     

Effects of bridge ions, DNA species, and developing temperature on flat-lying DNA monolayers

Recently, we have successfully constructed flat-lying DNA monolayers on a mica surface (J. Phys. Chem. B 2006, 110, 10792-10798). In this work, the effects of various factors including bridge ions, DNA species, and developing temperature on the configuration of DNA monolayers have been investigated by atomic force microscopy (AFM) in detail. AFM results show that the species of bridge ions and developing temperature play a crucial role during the formation process. For example, the divalent cation Zn2+ resulted in many DNA chains stuck side by side in the monolayers due to the strong interactions between it and DNA's bases or the mica surface. Most DNA chain's conglutinations disappeared when the developing temperature was higher than 40 degrees C. Cd2+ and Ca2+ produced more compact DNA monolayers with some obvious aggregations, especially for the DNA monolayers constructed by using Ca2+ as the bridge ion. Co2+ produced well-ordered, flat-lying DNA monolayers similar to that of Mg2+. Furthermore, it was found that the flat-lying DNA monolayers could still form on a mica surface when plasmid DNA pBR 322 and linear DNA pBR 322/Pst I were used as the DNA source. Whereas, it was hard to form DNA monolayers on a (3-aminopropyl)triethoxysilane-mica surface because the strong interactions between DNA and substrate prevented the lateral movement of DNA molecules. These results suggested that the appropriate interactions between divalent cations and DNA or mica surface were important for the formation of flat-lying DNA monolayers. The obtained information is a necessary supplement to our previous studies on the formation kinetics of such monolayers and may be useful for practical application of the monolayers and further theoretical studies.     

Parameterization of Ca+2-protein interactions for molecular dynamics simulations

Molecular dynamics simulations of Ca+2 ions near protein were performed with three force fields: GROMOS96, OPLS-AA, and CHARMM22. The simulations reveal major, force-field dependent, inconsistencies in the interaction between the Ca+2 ions with the protein. The variations are attributed to the nonbonded parameterizations of the Ca+2-carboxylates interactions. The simulations results were compared to experimental data, using the Ca+2-HCOO- equilibrium as a model. The OPLS-AA force field grossly overestimates the binding affinity of the Ca+2 ions to the carboxylate whereas the GROMOS96 and CHARMM22 force fields underestimate the stability of the complex. Optimization of the Lennard-Jones parameters for the Ca+2-carboxylate interactions were carried out, yielding new parameters which reproduce experimental data.     

The tetracycline: Mg2+ complex: a molecular mechanics force field

Tetracycline (Tc) is an important antibiotic, which binds specifically to the ribosome and several proteins, in the form of a Tc-:Mg2+ complex. To model Tc:protein and Tc:RNA interactions, we have developed a molecular mechanics force field model of Tc, which is consistent with the CHARMM force field for proteins and nucleic acids. We used structures from the Cambridge Crystallographic Data Base to identify the main Tc conformations that are likely to be present in solution and in biomolecular complexes. A conformational search was also done, using the MM3 force field to perform simulated annealing of Tc. Several resulting, low-energy structures were optimized with an ab initio model and used in developing the new Tc force field. Atomic charges and Lennard-Jones parameters were derived from a supermolecule ab initio approach. We considered the ab initio energies and geometries of a probe water molecule interacting with Tc at 36 different positions. We considered both a neutral and a zwitterionic Tc form, with and without bound Mg2+. The final rms deviation between the ab initio and force field energies, averaged over all forms, was just 0.35 kcal/mol. The model also reproduces the ab initio geometry and flexibility of Tc. As further tests, we did simulations of a Tc crystal, of Tc:Mg2+ and Tc:Ca2+ complexes in aqueous solution, and of a solvated complex between Tc:Mg2+ and the Tet repressor protein (TetR). With slight, ad hoc adjustments, the model can reproduce the experimental, relative, Tc binding affinities of Mg2+ and Ca2+. It performs well for the structure and fluctuations of the Tc:Mg2+:TetR complex. The model should therefore be suitable to investigate the interactions of Tc with proteins and RNA. It provides a starting point to parameterize other compounds in the large Tc family.     

二价金属离子对鲑鱼精DNA热稳定性的影响

二价金属离子与 DNA之间的作用在 DNA复制、转录以及新陈代谢过程中起到重要的作用 ,因此越来越受到关注 [1,2 ] .L uck等 [3 ]指出 ,碱土金属与 DNA分子的磷酸基团发生作用使 DNA的热稳定性升高 ,而二价过渡金属离子主要是与碱基作用而使 DNA热稳定性降低 . Eichhorn等 [4 ] 用变温紫外光谱法研究了在 DNA磷酸根离子的浓度较小 (0～ 4.0 )时 ,多种二价金属离子对小牛胸腺 DNA热稳定性的影响 ,发现 Mg2 + 及过渡金属离子 Mn2 + ,Co2 + 和 Ni2 + 都有利于 DNA的稳定 .最近 ,Duguid等 [5]用差示扫描量热法 (DSC)研究了在 c(M2 +…     

Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide

A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.     

Cu2+对腺嘧啶-胸腺嘧啶碱基对阴离子质子转移的影响

采用B3LYP/DZP++方法研究了腺嘌呤-胸腺嘧啶(A-T)碱基对阴离子(AT)-的单质子转移机理以及金属离子Cu2+对(AT)-碱基对质子转移的影响.(AT)-碱基对的单质子转移路径是由胸腺嘧啶N25位上的质子H26沿分子间的氢键N25-H26…N10转移到腺嘌呤的N10位.金属Cu2+可通过络合作用分别吸附在(AT)-碱基对O24、O28、N4、N13上,从而影响(AT)-碱基对中质子转移过程.Cu2+络合作用在胸腺嘧啶(T)的O24、O28上时,发生了从胸腺嘧啶到腺嘌呤方向上的单质子转移反应;而作用在腺嘌呤(A)的N4、N13上时,得到了双质子转移的稳定产物.     

金属离子Na+,Li+及Mg2+与C1O4-和NO3-形成的无水离子缔合物种的从头算量子化学研究

采用量子化学计算方法在B3 LYP/6-311++G**水平下对Na+,Li+和Mg2+与ClO4-和NO3-形成的离子缔合物种的结构以及v1-频率进行了研究,并将结果与SO42-和上述3种阳离子形成的物种进行了对比.在缔合物种结构方面,当阳离子数目≤2时,与SO42-体系相似,ClO4-和NO3-主要与阳离子形成双齿缔合结构,而当阳离子数目＞2时,特别是具有2个正电荷的Mg2+离子数目较多时,由于阳离子间的斥力更大,与阳离 子结合能力较弱的ClO4-和NO3-较难与其形成稳定的离子团簇,而在SO42-体系中,则易形成单齿缔合结构.在v1-频率的变化趋势方面,3种阴离子形成的缔合物种大体相同,说明无水离子团簇的频率变化主要受阳离子性质和缔合结构影响.虽然阴离子性质也有部分影响,但不占主要地位.