Density functional theory study on hydrogen‐bonded complexes of adenine with polyformamide molecules

The structures of hydrogen‐bonded complexes A–F_n (n = 2–7) of adenine with polyformamide molecules have been fully optimized at B3LYP/6‐31G(d) basis set level. All the formamide molecules prefer to be N? H proton donor rather than C? H proton donor and are favorably bound to the five‐numbered moiety of adenine. A displacement of formamide molecules to one side of adenine mean plane has happened with an increasing number of formamide molecules. An obvious effect of hydrogen‐bonding cooperativity can be seen during the complex process. The most interesting geometrical change of adenine upon the complex is the shortening of the bond C4? N6 resulting from the strengthening of the conjugation between the π system of the adenine ring and the lone pair of the nitrogen atom. An existence of weak N? H···π bonding interaction between the π system of adenine and N? H bond of F7 is found and further conformed by an natural bond orbital analysis specially carried out on A–F₇. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Intra and intermolecular hydrogen bonding in formohydroxamic acid

The presence of hydrogen bonding interactions in several tautomeric forms of formohydroxamic acid (FHA) and 1:1 association among the tautomeric forms and water‐coordinated tautomeric forms of FHA is explored theoretically. Out of the seven equilibrium structures, four tautomeric forms have been selected for aggregation with single water molecule and dimer formation. Fifteen aggregates of FHA with H₂O have been optimized at MP2/AUG‐cc‐PVDZ level and analyzed for intramolecular and intermolecular H‐bond interactions. Twenty‐seven dimers of the four tautomeric forms have been obtained at MP2/6‐31+G* level. The stabilization energies associated with dimerization and adduct formation with water are the result of H‐bond interactions and range from very weak to medium. The atomic charges and NBO analysis indicate that the electrostatic and the charge transfer are the important components favoring H‐bond formation. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Density functional study of indole–pyrrole heterodimer potential energy hypersurface

A density‐functional study of indole–pyrrole heterodimer potential energy hypersurface (PES) was performed. Eight stationary points were located on the B3LYP/6‐31++G(d,p) PES, three of which correspond to real minima, all of them being characterized with an N? H … π type hydrogen bonding. In two of these minima (the local ones), pyrrole subunit acts as a hydrogen bond proton donor, while the global minimum corresponds to indole–H … π(‐pyrrole) arrangement. Besides the interaction and dissociation energies corrected for BSSE and the monomer relaxation energies and the relevant structural parameters, anharmonic N? H and N? H … π vibrational frequencies were calculated for various N? H oscillators involved in this interaction from the 1‐D DFT vibrational potentials. On the basis of anharmonic vibrational frequency analysis, it was concluded that the two types of N? H … π hydrogen bonded dimers (indole vs. pyrrole being a proton donor) should be distinguishable with spectroscopic methods. Various contributions to the overall anharmonic frequency shifts upon hydrogen bonding were calculated and discussed as well. The charge field perturbation (CFP) technique was employed to study the electrostatic + polarization influence of the proton accepting unit on the N? H(… π) vibrational potential. The second‐order perturbation theory analysis (SOPT) of the Fock matrix (i.e., its Kohn–Sham analog) within the natural bond orbital (NBO) basis, as well as various NBO deletion analyses revealed an essentially one‐directional charge transfer (CT) of a π(C? C) → σ*(N? H) type in the case of all three minima. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Fast,efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method

The AM1‐BCC method quickly and efficiently generates high‐quality atomic charges for use in condensed‐phase simulations. The underlying features of the electron distribution including formal charge and delocalization are first captured by AM1 atomic charges for the individual molecule. Bond charge corrections (BCCs), which have been parameterized against the HF/6‐31G* electrostatic potential (ESP) of a training set of compounds containing relevant functional groups, are then added using a formalism identical to the consensus BCI (bond charge increment) approach. As a proof of the concept, we fit BCCs simultaneously to 45 compounds including O‐, N‐, and S‐containing functionalities, aromatics, and heteroaromatics, using only 41 BCC parameters. AM1‐BCC yields charge sets of comparable quality to HF/6‐31G* ESP‐derived charges in a fraction of the time while reducing instabilities in the atomic charges compared to direct ESP‐fit methods. We then apply the BCC parameters to a small “test set” consisting of aspirin, d ‐glucose, and eryodictyol; the AM1‐BCC model again provides atomic charges of quality comparable with HF/6‐31G* RESP charges, as judged by an increase of only 0.01 to 0.02 atomic units in the root‐mean‐square (RMS) error in ESP. Based on these encouraging results, we intend to parameterize the AM1‐BCC model to provide a consistent charge model for any organic or biological molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 132–146, 2000     

溶剂分子对甲酰胺及其衍生物分子内氢转移催化作用的研究

在B3LYP/6-311+ +G* * 水平下,通过计算所形成二元团簇的能量来研究水、氨、甲醇、氟化氢等溶剂分子对甲酰胺及其衍生物分子内氢原子转移的催化作用.简单描绘了在有水、氨、甲醇和氟化氢等溶剂分子存在时,甲酰胺及其衍生物分子内氢原子转移的过程.结果表明,当有水、氨、甲醇、氟化氢等溶剂分子存在时,从甲酰胺甲酰胺酸转变的能垒会降低.而且不同的溶剂分子对甲酰胺(FA)、 N-甲基甲酰胺(MF)和N,N-二甲基甲酰胺(DMF)的催化能力各不相同.在这四种溶剂分子中,氟化氢的催化作用最强.     

Modeling the H bond donor strength of ?OH, ?NH,and ?CH sites by local molecular parameters

A quantum chemical model is introduced to predict the H‐bond donor strength of monofunctional organic compounds from their ground‐state electronic properties. The model covers ? OH, ? NH, and ? CH as H‐bond donor sites and was calibrated with experimental values for the Abraham H‐bond donor strength parameter A using the ab initio and density functional theory levels HF/6‐31G** and B3LYP/6‐31G**. Starting with the Morokuma analysis of hydrogen bonding, the electrostatic (ES), polarizability (PL), and charge transfer (CT) components were quantified employing local molecular parameters. With hydrogen net atomic charges calculated from both natural population analysis and the ES potential scheme, the ES term turned out to provide only marginal contributions to the Abraham parameter A, except for weak hydrogen bonds associated with acidic ? CH sites. Accordingly, A is governed by PL and CT contributions. The PL component was characterized through a new measure of the local molecular hardness at hydrogen, η(H), which in turn was quantified through empirically defined site‐specific effective donor and acceptor energies, EE_occ and EE_vac. The latter parameter was also used to address the CT contribution to A. With an initial training set of 77 compounds, HF/6‐31G** yielded a squared correlation coefficient, r², of 0.91. Essentially identical statistics were achieved for a separate test set of 429 compounds and for the recalibrated model when using all 506 compounds. B3LYP/6‐31G** yielded slightly inferior statistics. The discussion includes subset statistics for compounds containing ? OH, ? NH, and active ? CH sites and a nonlinear model extension with slightly improved statistics (r² = 0.92). © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

A Hybrid statistical mechanics—quantum chemical model for proton transfer in 5‐azauracil and 6‐azauracil in water solution

A hybrid statistical physics—quantum‐chemical methodology was implemented to study the water‐assisted intramolecular proton‐transfer processes in 5‐ and 6‐azauracils in aqueous solutions. The solvent effects were included in the model by explicit inclusion of two pairs of water molecules, which model the relevant part of the first hydration shell around the solute. The position of these water molecules was initially estimated by carrying out a classical Metropolis of dilute water solutions of the title compounds and subsequently analyzing solute–solvent intermolecular interactions in the Monte Carlo‐generated configurations. Sequentially to the statistical physics simulation, ab initio quantum mechanical (QM) level of theory was implemented. The effects of the water as solvent (at ab initio QM level) were introduced at two different levels—using solute–solvent clusters (four‐water molecules) and using the same clusters embedded in an external continuum. Full geometry optimizations of these complexes were carried out at MP2/6–31 + G(d, p) and conductor‐polarizable continuum model (C‐PCM)/MP2/6–31 + G(d, p). Single point calculations were performed at CCSD(T)/6–31 + G(d, p)//MP2/6–31 + G(d, p) computational level to obtain more accurate energies. According to our calculations hydrated azauracils should exist in three forms: mainly dioxo form and two hydroxy forms. The calculated proton transfer activation energies for tautomeric reactions of 5‐azauracil and 6‐azauracil show different pictures for these two compounds. According to C‐PCM/MP2/6–31 +G(d, p) data, water‐assisted proton transfer in 5‐azauracil realizes through two parallel reactions: 1,3,5‐triazine‐2,4(1H,3H)‐dione → 6‐hydroxy‐1,3,5‐triazin‐2(1H)‐one and 1,3,5‐triazine‐2,4(1H,3H)‐dione → 4‐hydroxy‐1,3,5‐triazin‐2(1H)‐one. Tautomeric equilibrium in 6‐azauracil in water could occur by two contiguous reactions: 1,2,4‐triazine‐3,5(2H,4H)‐dione → 5‐hydroxy‐1,2,4‐triazin‐3(2H)‐one and 5‐hydroxy‐1,2,4‐triazin‐3(2H)‐one → 3‐hydroxy‐1,2,4‐triazin‐5(2H)‐one. The proton transfer investigated reactions in 5‐ and 6‐azauracils involve concerted atomic movement. © 2015 Wiley Periodicals, Inc.     

Determination of net atomic charges using a modified partial equalization of orbital electronegativity method. IV. Application to hypervalent sulfur- and phosphorus-containing molecules

The parameters for an empirical point charge calculation method, Mulliken Population Constrained Potential Derived–Modified Partial Equalization of Orbital Electronegativity (MPCPD–MPEOE) method, for hypervalent sulfur, phosphorus, and the atoms bonded to these atoms, were determined. A new empirical method, MPCPD, which has inherent advantages over both the Mulliken population and potential derived methods, is proposed here. The MPCPD net atomic charges are transferable and reproduce the electrostatic potential. The electrical properties were calculated with the MPCPD and MPCPD–MPEOE charges. These properties agreed well with experimental results and with 6-31G** ab initio results. The MPEOE parameters for the sulfur atom in aromatic molecules and trivalent phosphorus in phosphite molecules were also determined. The dipole moments calculated with the MPCPD–MPEOE point charges agreed well with experimental results. All the parameters obtained in this work are consistent with those obtained in our previous work. © 1995 by John Wiley & Sons, Inc.     

The influence of formamide (model of protein unit) on the intramolecular proton transfer in the DNA simple base guanine: A density functional theory study

For the purpose of investigating the influence of protein unit on the intramolecular proton transfer (IPT) reactions in the simple base guanine, a simple model (formamide) of peptides is designed to biological system investigations, and five complexes of formamide–guanine (FG1, FG2, FG3, FG4, and FG5) are determined at the B3LYP/6‐311++G(d,p) level of theory. For comparison, HF and MP2 methods are also used in this paper. The proton transfer (PT) reaction processes of guanine and FGs have been investigated employing the B3LYP/6‐311++G(d,p) level of theory. The selected thermodynamic and kinetic parameters, such as the activation energies (E_a), changes of enthalpy (ΔH) and changes of free energies (ΔG), as well as the equilibrium constants (K_p) for those reaction processes, have also been obtained by calculational means. The calculated results indicate that the assisted and protected effects of formamide on IPT in guanine are site‐dependent. CH1 is the lowest activation energy needed PT process no matter where the formamide molecule is located in. The activation energy of CH1 with formamide in S2 is the lowest one (153.3 KJ/mol), whereas the one of CH5 with formamide in S5 is the highest (318.3 KJ/mol). © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Ab initio study of gas phase and water-assisted tautomerization of maleimide and formamide

Maleimide serves as an important starting material in the synthesis of drugs and enzyme inhibitors. In the present paper, knowing the importance of tautomerization in maleimide for its drug action, potential energy surface of maleimide is studied and its tautomerization has been discussed and compared with tautomerization of formamide. Gas phase tautomerization of maleimide requires large amount of energy (23·21 kcal/mol) in comparison to formamide (15·05 kcal/mol) at HF/6-31+G* level. Thus making the proton transfer reaction a difficult process in gas phase. Water molecule lowers the energy barrier of tautomerization thus facilitating the tautomerization of maleimide to 5-hydroxy-pyrrol-2-one. Water assisted tautomerization of maleimide requires 19·60 kcal/mol energy at HF/6-31+G* and 17·63 kcal/mol energy at B3LYP/6-31+G* level, a decrease of 3·61 and 5·96 kcal/mol over gas phase tautomerization. Whereas, tautomerization of formamide requires 14·16 and 12·84 kcal/mol energy, a decrease of 0·89 and 2·01 kcal/mol energy over gas phase tautomerization at HF/6-31+G* and B3LYP/6-31+G* level, respectively. Water-assisted tautomerization in maleimide and formamide showed that difference in energy barrier reduces to 2·83 kcal/mol from 10·41 kcal/mol (in gas phase) at B3LYP level, which resulted that maleimide readily undergoes tautomerization in water molecule.     

Electrostatic model for infrared intensities in a spectroscopically determined molecular mechanics force field

A new electrostatic model for the calculation of infrared intensities in molecular mechanics and molecular dynamics is presented. The model is based on atomic charges, atomic charge fluxes, and internal coordinate dipoles and their fluxes. The internal coordinate dipoles are used instead of atomic dipoles, thus simplifying the derivation of parameters. The model is designed to reproduce ab initio dipole derivatives, and the parameters can be obtained by (iterative) transformations from these, or by linear least squares fitting to them. A first application to linear alkanes has been made. For these molecules, the intensities can be predicted with an average accuracy of 30–40%. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 754–768, 1998     

Various atomic charge calculation schemes of CoMFA on HIF‐1 inhibitors of moracin analogs

Selection of appropriate partial charges in a molecule is crucial to derive good quantitative structure–activity relationship models. In this work, several partial atomic charges were assigned and tested in a comparative molecular field analysis (CoMFA) models. Many CoMFA models were generated for a series of hypoxia inducible factor 1 (HIF‐1) inhibitors using various partial atomic charges including charge equalization, Mülliken population analysis (MPA), natural population analysis, and electrostatic potential (ESP)‐derived charges. These atomic charges were investigated at various theoretical levels such as empirical, semiempirical, Hartree–Fock (HF), and density functional theory (DFT). Among them, Merz‐Singh‐Kollman (MK) ESP‐derived charges at the level of HF/6‐31G* gave the highest predictive q² with experimental pIC₅₀ values. With this charge scheme, a detailed analysis of CoMFA model was performed to understand the electrostatic interactions between ligand and receptor. More elaborate charge calculation schemes such as HF and DFT correlated more strongly with activity than empirical or semiempirical schemes. The choice of optimization methods was important. As geometries were fully optimized at the given levels of theory, the aligned structures were different. They differed considerably, especially for the flexible parts. This was likely the source of the substantial variation of q² values, even when the same steric factor was considered without electrostatic parameters. ESP‐derived charges were most appropriate to describe CoMFA electrostatic interactions among MPA, NBA, and ESP charges. Overall q² values vary considerably (0.8–0.5) depending on the charge schemes applied. The results demonstrate the need to consider more appropriate atomic charges rather than default CoMFA charges. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Experimental and Computational Studies on the 1:1 Complex of Anthranilic Acid with P-Toluenesulfonic Acid

Abstract

The proton-transfer compound (3) was synthesized in good yields by the reaction of anthranilic acid with p-toluenesulfonic acid. It was characterized by elemental analysis, infrared spectroscopy, and X-ray single-crystal determination. The crystal structure analysis of 2-carboxyanilinium p-toulenesulfonate (3) has revealed a one-dimensional hydrogen-bonded network structure, involving the tosylate anion, the carboxyl group, and the ammonium group. The H(N)···O distances range from 1.97 to 2.23 Å. The molecular geometry and vibrational frequencies of 3 were calculated using the ab-initio method (HF) with the 6–31G(d) and 6–31+G(d,p) basis sets. The computed results indicate that the optimized geometry reproduces the crystal structure well, and the assignments of fundamental vibrations also agree well with the theoretical frequencies. The intermolecular proton transfer process between the ionic (3) and nonionic (4) structures was also investigated with the theoretical computations. The nonionic form (4) is energetically more stable than the ionic form (3) and TS(3→4) by 9.76 and 7.01 kcal/mol, respectively, including the zero-point vibrational energy correction at the HF/6–31+G(d,p) level. In addition, the atomic charges, the molecular electrostatic potentials, the nucleus-independent chemical shifts, and the frontier molecular orbitals of 3 were carried out at the HF/6–31+G(d,p) level of theory.

Supplementary materials are available for this article. Go to the publisher's online edition of Phosphorus, Sulfer, and Silicon and the Related Elements for the following free supplemental files: Additional figures and tables     

Conformational preferences for hydroxyl groups in substituted tetrahydropyrans

Quantum mechanical (ab initio and semiempirical) and force field calculations are reported for representative torsion potentials in several tetrahydropyran derivatives. The overall agreement between the various methods is quite good except that the AMBER torsion profiles are sensitive to the choice of atomic point charges. Using electrostatic potential (ESP) derived atomic point charges determined with the STO-3G basis set we find that AMBER is able to match the best quantum mechanical results quite well. However, when the point charges are derived using the 6-31G* basis set we find that scaling the intramolecular electrostatic nonbond interactions is necessary. AM1 does not work very well for these compounds when compared to the ab initio methods and, therefore, should only be used in cases when ab initio calculations would be prohibitive. Based upon our results we feel that any force field that makes use of 6-31G* ESP derived atomic point charges will need to scale intramolecular interactions. Implications of scaling intramolecular interactions to the development of force fields based on 6-31G* ESP derived atomic point charges are discussed. © 1992 by John Wiley & Sons, Inc.     

Calculation of the electrostatic charges and energies for intercalation of aromatic drug molecules with DNA

In this work, we analyzed the influence of the charge model on the magnitudes of atomic charges and electrostatic energies for the binding of aromatic drug molecules with DNA. The dependence of the charge and energy on the level of theory (HF, DFT (B3LYP), MP2, semi‐empirical methods), basis set (STO‐3G, 3‐21G, 6‐31G, 6‐31G*, 6‐31G**), method of charge computation (Mulliken, Natural Population Analysis, CHelpG, Merz–Kollman), and force field charge (CHARMM27, AMBER99) has been tracked for typical aromatic drugs of different structure and charge state. Recommendations and restrictions have been formulated for the use of particular approaches in charge/electrostatic energy calculations. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Effective atomic charges in alanine dipeptide

A new set of effective atomic charges of different conformers of alanine dipeptide is presented. These charges are obtained by fitting the electrostatic potential resulting from the ab initio SCF wave function of the system obtained in a 6-31G basis set. A specific fit procedure is used providing charges weakly dependent on the fit points as well as on the geometry of the molecule. It is shown that these charges retain a reasonable chemical meaning. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 473–482, 1999     

4-羟甲基吡啶－水红移氢键的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-311++G**基函数对4-羟甲基吡啶与水形成的1:1和1:2（摩尔比）氢键复合物进行了理论计算研究,分别得到稳定的4-羟甲基吡啶－H₂O和4-羟甲基吡啶－(H₂O)₂氢键复合物3个和8个。经基组重叠误差和零点振动能校正后,最稳定的1:1和1:2氢键复合物的相互作用能分别为-20.536和-44.246 kJ/mol。振动分析显示O-H···N(O)氢键的形成使复合物中O-H键对称伸缩振动频率红移(减小)。自然键轨道分析表明,4-羟甲基吡啶与水形成最稳定的1:1和1:2氢键复合物时,分子间电荷转移分别为0.02642 e 和0.03813 e 。含时密度泛函理论TD-B3LYP/ 6-311++G**计算显示,相对于4-羟甲基吡啶单体分子,氢键H-OH···N和H-OH···OH的形成分别使最大吸收光谱波长兰移8~16纳米和红移4~11纳米。     

ChelpG and QTAIM atomic charge and dipole models for the infrared fundamental intensities of the fluorochloromethanes

ChelpG atomic charges and dipoles and the charge–charge flux–dipole flux (CCFDF) model have been used to quantitatively estimate the fundamental infrared intensities of the fluorochloromethanes. Since the ChelpG calculational procedure includes the constraint that the atomic charges and dipoles reproduce the equilibrium dipole moments the model results in accurate intensity values that have a root mean square error of 0.7 km mol^?1 compared to those determined directly from the MP2/6-311G++(3d,3p) electronic density and 23.1 km mol^?1 relative to the experimental intensities. Although these ChelpG results for total dipole moment derivatives are almost the same as those obtained previously using QTAIM (Quantum Theory of Atoms in Molecules) atomic charges and dipoles in the CCFDF model, their charge, charge flux and dipole flux contributions are completely different. Whereas the contributions calculated using the QTAIM parameters have values following expectations based on electronegativity concepts this is not true for those obtained from the ChelpG parameters. Mean dipole moment derivatives determined from experimental fundamental infrared intensities are compared with the ChelpG and QTAIM atomic charges. Furthermore, Generalized Atomic Polar Tensor Charges (GAPT) are found to need correction for their dynamic contributions if they are to be used as static atomic charges.     

Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations

Ab initio MP2/6-31G* interaction energies were calculated for more than 80 geometries of stacked cytosine dimer. Diffuse polarization functions were used to properly cover the dispersion energy. The results of ab initio calculations were compared with those obtained from three electrostatic empirical potential models, constructed as the sum of a Lennard-Jones potential (covering dispersion and repulsion contributions) and the electrostatic term. Point charges and point multipoles of the electrostatic term were also obtained at the MP2/6-31G* level of theory. The point charge MEP model (atomic charges derived from molecular electrostatic potential) satisfactorily reproduced the ab initio data. Addition of π-charges localized below and above the cytosine plane did not affect the calculated energies. The model employing the distributed multipole analysis gave worse agreement with the ab initio data than the MEP approach. The MP2 MEP charges were also derived using larger sets of atomic orbitals: cc-pVDZ, 6-311 + G(2d, p), and aug-cc-pVDZ. Differences between interaction energies calculated using these three sets of point charges and the MP2/6-31G* charges were smaller than 0.8 kcal/mol. The correlated ab initio calculations were also compared with the density functional theory (DFT) method. DFT calculations well reproduced the electrostatic part of interaction energy. They also covered some nonelectrostatic short-range effects which were not reproduced by the empirical potentials. The DFT method does not include the dispersion energy. This energy, approximated by an empirical term, was therefore added to the DFT interaction energy. The resulting interaction energy exhibited an artifact secondary minimum for a 3.9-4.0 vertical separation of bases. This defect is inherent in the DFT functionals, because it is not observed for the Hartree-Fock + dispersion interaction energy.© 1996 John Wiley & Sons, Inc.     

Potential-Derived point-charge model study of electrostatic interaction energies in some hydrogen-bonded systems

Mulliken's atomic charges (MC ) and potential derived (PD ) point charges obtained from STO -3G wave functions are used to study the electrostatic interaction energies for a series of representative hydrogenbonded complexes. The results of the above-mentioned models are compared with the more accurate results of segmental multipole moment (SMM ) expansion, and it is shown that the PD model is superior to the Mc model. The results of PD model are shown to be well correlated with the results of SMM expansion technique. Results of our calculations using 6-31G and 6-31G** PD charges are also reported here. Electrostatic interaction energies obtained using 6-31G** PD charges are compared with the 6-31G** SCF interaction energies available for the nine hydrogen-bonded dimers of ammonia, water, and hydrogen fluoride and a good con-elation between the two is shown. The interrelationship between the results of different basis sets are also examined for the PD point-charge model. The electrostatic interaction energies obtained using STO -3G PD model are shown to be well correlated to the results of 6-31G and 6-31G** PD models.