Conformations and reaction energetics are important for understanding the interactions be- tween biomolecules and metal ions. In this work, we report a systematic ab initio study on the conformations and metal ion aff...Conformations and reaction energetics are important for understanding the interactions be- tween biomolecules and metal ions. In this work, we report a systematic ab initio study on the conformations and metal ion affinities of glutamine （Gln） binding with alkali and alkaline earth metal ions. An efficient and reliable method of searching low energy conformations of metalated Gln is proposed and applied to the complexes of Gln.Gln.M^＋/＋＋（M＋/＋＋=Li^＋、Na^＋、K^＋、Rb^＋、Cs^＋、Be^＋＋、Mg^＋＋、Ca^＋＋、Sr^＋＋ and Ba^＋＋）.. In addition to all conformers known in literatures, many new important conformations are located, demonstrating the power of the new method and the necessity of the conformational search performed here. The metal coordination modes, relative energies, dipole moments, and equilibrium distributions of all important conformations of Gln.M^＋/＋＋ are calculated by the methods of B3LYP, BHandHLYP, and MP2. IR spectra and metalation enthalpies and free energies are also presented and compared with the available experiments. The results form an extensive database for systematic examination of the metalation properties of Gln.展开更多
An extensive computational study on the conformations of gaseous dipeptide glycine- arginine, GlyArg, has been performed. A large number of trail structures were generated by systematically sampling the potential ener...An extensive computational study on the conformations of gaseous dipeptide glycine- arginine, GlyArg, has been performed. A large number of trail structures were generated by systematically sampling the potential energy surface （PES） of OlyArg. The trial structures were successively optimized with the methods of PM3, HF/3-21G＊, BHandHLYP/6-31G＊ and BHandHLYP/6-311＋＋G＊＊ in order to reliably find the low energy conformations. The eonformational energies were finally determined with the methods of BHandHLYP, cam- B3LYP, B97D, and MP2 using the basis set of 6-311＋＋G（3df,3pd）. The results establish firmly that gaseous GlyArg exists primarily in its canonical form, in sharp contrast with ArgGly that adopts the zwitterionic form. Important data such as the rotational constants, dipole moments, vertical ionization energies, temperature distributions and IR spectra of the low energy conformers are represented for the understanding of the future experiments. Moreover, considering the global minima of all amino acids and many dipeptides, combined with the hydrophobicities of amino acids, a model predicting whether the global minimum configuration of a dipeptide is canonical or zwitterionic is developed.展开更多
Unrestrained molecular dynamics （MD） simulations have been carded out to characterize the stability of DNA conformations and the dynamics of A-DNA^B-DNA conformational transitions in aqueous RbC1 solutions. The PARM...Unrestrained molecular dynamics （MD） simulations have been carded out to characterize the stability of DNA conformations and the dynamics of A-DNA^B-DNA conformational transitions in aqueous RbC1 solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbC1 concentration on the dynamics of the A^B conformational tran- sition in the DNA duplex d（CGCGAATTCGCG）2. Canonical A- and B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was moni- tored for 3.0 ns and the distances between the C5＇ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simula- tions. In addition, increasing the RbC1 concentration in aqueous solution hindered the A^B conformational transition and the transition in aqueous RbC1 solution was faster than that in aqueous NaC1 solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A^B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The ru- bidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help ex- plain the effect of RbC1 concentration on the dynamics of the A^B conformational transition.展开更多
基金ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China （No.11074233 and No.11374272） and the Specialized Research Fund for the Doctoral Program of Higher Education （No.20113402110038 and No.20123402110064）
文摘Conformations and reaction energetics are important for understanding the interactions be- tween biomolecules and metal ions. In this work, we report a systematic ab initio study on the conformations and metal ion affinities of glutamine （Gln） binding with alkali and alkaline earth metal ions. An efficient and reliable method of searching low energy conformations of metalated Gln is proposed and applied to the complexes of Gln.Gln.M^＋/＋＋（M＋/＋＋=Li^＋、Na^＋、K^＋、Rb^＋、Cs^＋、Be^＋＋、Mg^＋＋、Ca^＋＋、Sr^＋＋ and Ba^＋＋）.. In addition to all conformers known in literatures, many new important conformations are located, demonstrating the power of the new method and the necessity of the conformational search performed here. The metal coordination modes, relative energies, dipole moments, and equilibrium distributions of all important conformations of Gln.M^＋/＋＋ are calculated by the methods of B3LYP, BHandHLYP, and MP2. IR spectra and metalation enthalpies and free energies are also presented and compared with the available experiments. The results form an extensive database for systematic examination of the metalation properties of Gln.
基金This work was supported by the National Natu- ral Science Foundation of China （No.11374272） and the Specialized Research Fund for the Doctoral Pro- gram of Higher Education （No.20113402110038 and No.20123402110064）.
文摘An extensive computational study on the conformations of gaseous dipeptide glycine- arginine, GlyArg, has been performed. A large number of trail structures were generated by systematically sampling the potential energy surface （PES） of OlyArg. The trial structures were successively optimized with the methods of PM3, HF/3-21G＊, BHandHLYP/6-31G＊ and BHandHLYP/6-311＋＋G＊＊ in order to reliably find the low energy conformations. The eonformational energies were finally determined with the methods of BHandHLYP, cam- B3LYP, B97D, and MP2 using the basis set of 6-311＋＋G（3df,3pd）. The results establish firmly that gaseous GlyArg exists primarily in its canonical form, in sharp contrast with ArgGly that adopts the zwitterionic form. Important data such as the rotational constants, dipole moments, vertical ionization energies, temperature distributions and IR spectra of the low energy conformers are represented for the understanding of the future experiments. Moreover, considering the global minima of all amino acids and many dipeptides, combined with the hydrophobicities of amino acids, a model predicting whether the global minimum configuration of a dipeptide is canonical or zwitterionic is developed.
基金support from the National Natural Science Foundation of China (21176132)the Specialized Research Fund for the Doctoral Program of Higher Education (2010000211024)
文摘Unrestrained molecular dynamics （MD） simulations have been carded out to characterize the stability of DNA conformations and the dynamics of A-DNA^B-DNA conformational transitions in aqueous RbC1 solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbC1 concentration on the dynamics of the A^B conformational tran- sition in the DNA duplex d（CGCGAATTCGCG）2. Canonical A- and B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was moni- tored for 3.0 ns and the distances between the C5＇ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simula- tions. In addition, increasing the RbC1 concentration in aqueous solution hindered the A^B conformational transition and the transition in aqueous RbC1 solution was faster than that in aqueous NaC1 solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A^B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The ru- bidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help ex- plain the effect of RbC1 concentration on the dynamics of the A^B conformational transition.
