Theoretical studies of the tautomeric equilibria for five-member N-heterocycles in the gas phase and in solution

Tautomeric equilibria have been studied for five-member N-heterocycles and their methyl derivatives in the gas phase and in different solvents with dielectric constants of epsilon = 4.7-78.4. The free energy changes differently for tautomers upon solvation as compared to the gas phase, resulting in a shift of the equilibrium constant in solution. Solvents with increasing dielectric constant produce more negative solute-solvent interaction energies and increasing internal energies. The methyl-substituted imidazole and pyrrazole form delicate equilibria between two tautomeric forms. Depending on the solvent, the methyl-substituted triazoles and tetrazole have one or two major tautomers in solution. When estimating the relative solvation free energies by means of an explicit solvent model and using the FEP/MC method, one observes that the preferred tautomers differ in several cases from those predicted by the continuum solvent model. The 1,2-prototropic shift, as an intramolecular tautomerization path, requires about 50 kcal/mol activation energy for imidazole in the gas phase, and this route is also disfavored in a solution. The calculated activation free energy along the intramolecular path is 48-50 kcal/mol in chloroform and water as compared to a literature value of 13.6 kcal/mol for pyrrazole in DMSO. A molecular dynamics computer experiment favors the formation of an imidazole chain in chloroform, making the 1,3-tautomerization feasible along an intermolecular path in nonprotic solvents. In aqueous solution, one strong N-H...Ow hydrogen bond is formed for each species, whereas all other nitrogens in the ring form weaker, N...HwOw type hydrogen bonds. The tetrahydrofuran solvent acts as a hydrogen bond acceptor and forms N-H...Oether bonds. Molecules of the dichloromethane solvent are in favorable dipole-dipole interactions with the solute. The results obtained are useful in the design of N-heterocyclic ligands forming specified hydrogen bonds with protein side chains.     

Carbon Dioxide Capture and Release by Anions with Solvent‐Dependent Behaviour: A Theoretical Study

Recently, Clyburne and co‐workers [Science, 2014 , 344, 75–78] reported the novel synthesis of the elusive cyanoformate anion, NCCO₂^?. The stability of this anion is dependent on the dielectric constant of the local environment (polarity‐switchable solvent): it is stable in low‐polarity media and unstable in high‐polarity solvents; hence, capturing and then releasing CO₂. The possibility of extending such behaviour to other anions is explored herein. Specifically, the CO₂ capture process is studied for 26 anions in the gas phase and 3 distinct solvents (water, tetrahydrofuran, and toluene) by using the polarisable continuum model. Calculations are performed with the M06‐2X and B3LYP‐D3 density functionals and the aug‐cc‐pVTZ basis set. The design of new CO₂ complexes with the anion, which can be formed or destroyed on demand by changing the solvent, is possible; the results for the alkoxylate and thiolate anions are especially promising. The nature of the substituents connected to the atom that bonds to CO₂ in the anion is crucial in modulating the relative stability of the products—a key point for reversibility in the CO₂ capture process. A moderate interaction for the anion–CO₂ adduct—about 10 kcal mol^?1 relative free energy with respect to the isolated reactants in the gas phase—and a relevant effect in the dielectric constant of the local environment are also key ingredients to achieve solvent dependency.     

2,6-二巯基吡啶互变异构平衡体系溶剂效应的理论研究

在气相及甲苯、氯仿、乙腈和水等溶剂中对2,6-二巯基吡啶及其硫酮式互变异构体进行了HF/6-31G^**水平上的优化,其中溶液中的计算采用Onsager自洽反应场(SCRF)模型.探讨了溶剂对体系几何结构和能量的影响.结果表明:溶剂的存在与极性的增加有利于平衡体系中硫酮式异构体的存在.     

Effects of solvent on weak halogen bonds: Density functional theory calculations

In recent years, many applications of solution‐phase halogen bonding in anion recognition, catalysis, and pseudorotaxane formation have been reported. Moreover, a number of thermodynamic data of halogen bonding interactions in organic solution are now available. To obtain detailed information of the influence of the surrounding medium on weak halogen bonds, a series of dimeric complexes of halobenzene (PhX) with three electron donors (H₂O, HCHO, and NH₃) were investigated by means of DFT/PBE calculations in this work. The PCM implicit solvation approach was utilized to include the effects of three solvents (cyclohexane, chloroform, and water) as representatives for a wide range of dielectric constant. In some cases, halogen‐bond distances are shown to shorten in solution, accompanied by concomitant elongation of the C? X bonds. For the remaining systems, the intermolecular distances tend to increase or remain almost unchanged under solvent effects. In general, the solvent has a slight destabilizing effect on weak halogen bonds; the strength order of halogen bonds observed in vacuum remains unchanged in liquid phases. Particularly, the interaction strength attenuates in the order I > Br > Cl in solution, consistent with the experimental measurements of weak halogen bond door abilities. The similarities between halogen and hydrogen bonding in solution were also elucidated. The results presented herein would be very useful in future applications of halogen bonding in molecular recognition and medicinal chemistry. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Detection of Significant Aprotic Solvent Effects on the Conformational Distribution of Methyl 4‐Nitrophenyl Sulfoxide: From Gas‐Phase Rotational to Liquid‐Crystal NMR Spectroscopy

The conformational equilibrium of methyl 4‐nitrophenyl sulfoxide (MNPSO) was experimentally investigated in the gas phase by using microwave spectroscopy and in isotropic and nematic liquid‐crystal solutions, in which the solvents are nonaqueous and aprotic, by using NMR spectroscopy; moreover, it was theoretically studied in vacuo and in solution at different levels of theory. The overall set of results indicates a significant dependence of the solute conformational distribution on the solvent dielectric permittivity constant: when dissolved in low‐polarity media, the most stable conformation of MNPSO proved to be strongly twisted with respect to that in more polar solvents, in which the conformational distribution maximum essentially coincides with that obtained in the gas phase. We discuss a possible explanation of this behavior, which rests on electrostatic solute–solvent interactions and is supported by calculations of the solute electric dipole moment as a function of the torsional angle. This function shows that the least polar conformation of MNPSO is located at a twist angle close to that of the conformational distribution maximum found in less‐polar solvents. This fact, associated with a relatively flat torsional potential, can justify the stabilization of the twisted conformation by the less‐polar solvents.     

Origin of the solvent effects on the barrier to amide isomerization from the combined QM/MM Monte Carlo simulations

The origin of the solvent effects on the free energy of activation for the isomerization of N,N-dimethylformamide in water, CHCl₃ and CCl₄ has been investigated through statistical mechanical simulations using the combined quantum mechanical and molecular mechanical AM1/OPLS potential. The differential solvations between the ground state and transition state in various solvents can be attributed to the differences in molecular dipole moments in solution, and to the solvent polarization effects. In polar solvents, DMF is polarized more favorably in the ground state than in the rotated conformers, leading to greater solvent contributions. The modest solvent effects in CCl₄ are a reflection of its much smaller dielectric constant.     

The ab initio study and NBO analysis of the solvent dielectric constant effects and the implicit water molecules on the structural stability of the 1-(6-chloroquinoxalin-2-yl)-2-(4-(trifluoromethyl)-2,6-dinitrophenyl) hydrazine

Ab initio molecular orbital (MO) and density functional theory (DFT) methods were used to analyze the structure and the relative stability of 1-(6-chloroquinoxalin-2-yl)-2-(4-(trifluoromethyl)-2,6-dinitrophenyl) hydrazine in gas phase and the different solvent media. The effects of solvent dielectric constant and the implicit water molecules were investigated on the structural stability and intramolecular interactions. All used methods revealed that by the increase of solvent dielectric constant, the relative stability of the considered compound increase. Hence, the most stable structure is perceived in aqueous solution. Furthermore, natural bond orbital (NBO) analysis demonstrated that in the presence of implicit water molecules, the lone pair electrons of nitrogen have the most contribution in the resonance interactions of the aromatic rings and their stability. These facts may be the probable reasons behind the structural stability of the considered structure in the water solution based on energetic data and NBO analysis at the microscopic level.     

Solvent as a Probe of Active Site Motion and Chemistry during the Hydrogen Tunnelling Reaction in Morphinone Reductase

The reductive half‐reaction of morphinone reductase involves a hydride transfer from enzyme‐bound β‐nicotinamide adenine dinucleotide (NADH) to a flavin mononucleotide (FMN). We have previously demonstrated that this step proceeds via a quantum mechanical tunnelling mechanism. Herein, we probe the effect of the solvent on the active site chemistry. The pK_a of the reduced FMN N1 is 7.4±0.7, based on the pH‐dependence of the FMN midpoint potential. We rule out that protonation of the reduced FMN N1 is coupled to the preceding H‐transfer as both the rate and temperature‐dependence of the reaction are insensitive to changes in solution pH above and below this pK_a. Further, the solvent kinetic isotope effect is ～1.0 and both the 1° and 2° KIEs are insensitive to solution pH. The effect of the solvent’s dielectric constant is investigated and the rate of H‐transfer is found to be unaffected by changes in the dielectric constant between ～60 and 80. We suggest that, while there is crystallographic evidence for some water in the active site, the putative promoting motion involved in the H‐tunnelling reaction is insensitive to such changes.     

Glycine Dimers: Structure,Stability, and Medium Effects

We report a study on different ionization states and conformations of the bimolecular (Gly)₂ system by means of quantum mechanical calculations. Optimized geometries for energy minima of the glycine dimer, as well as relative energies and free energies were computed as a function of the medium: gas phase, nonpolar polarizable solvent, and aqueous solution. The polarizable continuum model was employed to account for solvation effects. Energy calculations were done using the MP2/aug‐cc‐pVTZ and B3LYP/6‐311+G(2df,2p) methods on B3LYP/6‐31+G(d,p) optimized structures (some single‐point energy calculations were also done using the B3PW91 and PBE1KCIS methods). Ionized forms of the glycine dimer (either zwitterion–zwitterion or neutral–zwitterion) are predicted to exist in all media, in contrast to amino acid monomers. In aqueous solution, dimerization is an exergonic process (?4 kcal mol^?1). Thus, according to our results, zwitterion–zwitterion Gly dimers might be abundant in supersaturated glycine aqueous solutions, a fact that has been connected with the structure of α‐glycine crystals but that remains controversial in the literature. Another noticeable result is that zwitterion–zwitterion interactions are substantially underestimated when computed using methods based on density functional theory. For comparison, some calculations for the dimer of the simplest chiral amino acid alanine were done as well and differences to the glycine dimer are discussed.     

Theoretical investigation of solvent effects on tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol

The density functional theory has been used to study the tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol(DDNP) in the gas phase and in 14 solvents at the B3LYP/6‐31G* level. The solvent effects on the tautomeric equilibria were investigated by the self‐consistent reaction field theory (SCRF) based on conductor polarized continuum model (CPCM) in apolar and polar solvents and by the hybrid continuum‐discrete model in protic solvent, respectively. Solvent effects on the computed molecular properties, such as molecular geometries, dipole moments, E_LUMO, E_HOMO, total energies for DDNP tautomers and transition state, tautomerization energies and solvation energies have been found to be evident. The tautomeric equilibrium of DDNP is solvent‐dependent to a certain extent. The tautomer I (cyclic azoxy form) is preferred in the gas phase, while in nonpolar solvents tautomer I and II (quinold form) exist in comparable amounts, and in highly polar solvents, the tautomeric equilibrium is shifted in favor of the more polar tautomer II . © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

DFT and Ab initio computational study on the reactivity sites of the GABA and its agonists,such as CACA,TACA, DABA,and muscimol: In the gas phase and dielectric media

The reactive behavior of GABA and its agonist molecules have been investigated using B3LYP hybrid density functional method at the 6‐311++G** basis set level, in the gas phase and dielectric media. The calculations have been performed to obtain optimized geometries, relative energies, net atomic charges, HOMO, and LUMO energies. Solvent effects have been analyzed using isodensity polarized continuum model (IPCM) for four different solvents, which are chloroform, ethanol, DMSO, and H₂O. The results have suggested that both the amino group and carboxyl group have been two of the most important active site for interactions between all compounds and their receptors. In addition they have suggested that the amino group reflects the difference in biological activity for each all molecules. The structures of these agonists have provided an essential foundation for subsequent structure‐activity analysis of ligand binding at GABA receptors, neuronal uptake inhibitors and transporters. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Systematic theoretical investigations on the tautomers of thymine in gas phase and solution

The tautomers of thymine are systematically calculated using various methods. The order of the relative stability, dipole moment as well as solvent effect of the 13 isomers are investigated. The one-step transition process by proton transfer in gas phase and in solution (H₂O or CH₃OH), which acts as both proton donor and acceptor, are exhaustively studied, meanwhile, the internal rotation between two related isomers in the gas phase is also investigated. Furthermore, the energy barrier of each transition is calculated and the results show that both H₂O and CH₃OH could reduce the energy barrier. Comparisons between the activation energies of all the producing processes of each tautomer as well as temperature effect on the producing processes are made. The calculation results also indicated that the diketo is the most stable isomer both in gas phase and in solution, additionally, it is thermodynamically as well as dynamically favored.     

Intramolecular proton transfer of serine in aqueous solution. Mechanism and energetics

Serine amino acid in aqueous solution is theoretically studied at the B3PW91/6-31+G^** level using a dielectric continuum solvent model. Neutral and zwitterionic structures in the gas phase and in solution are described and the proton-transfer mechanism is discussed. A neutral conformation in which the carboxyl hydrogen atom is already oriented toward the amino group seems to be the absolute energy minimum in the gas phase and the most stable neutral form in solution. The absolute energy minimum in solution is a zwitterionic form. The energy barrier for proton transfer is predicted to be very small, in particular when zero-point-energy contributions are added. Our calculations allow the dynamic aspects of the ionization mechanism to be discussed by incorporating nonequilibrium effects. Received: 28 June 1999 / Accepted: 13 October 1999 / Published online: 14 March 2000     

Solvent effect on analyte charge state,signal intensity,and stability in negative ion electrospray mass spectrometry; implications for the mechanism of negative ion formation

Department of Chemistry, University of New Orleans, New Orleans, Louisiana, USA The effect of solvent composition on negative ion electrospray ionization (ESI) mass spectrometry was examined. The onset potentials for ES1 of a series of chlorinated solvents and methanol were found to be within the range predicted by D. P. H. Smith, based on differences in the surface tension of the solvents used. The tendency toward electric discharge decreased with increasing percent weight of chlorine in the solvent. This effect has been attributed to an increasing propensity for electron capture for more highly chlorinated solvents. Addition of the electron scavenger gas SF, was even more effective at suppressing corona discharge phenomena. In a comparison of ultimate signal intensity obtainable for a test analyte in 10% methanol, the highest signal, which was stable over the widest range of temperatures, was exhibited by chloroform compared to dichloromethane, 1,2-dichloroethane, carbon tetrachloride, and methanol (100%). Chloroform, thus, is a recommended solvent for negative ion electrospray mass spectrometry (ES/MS) when solubility is not a limiting issue. Solvent polarity was shown to exhibit a profound influence on the distribution of charge states in negative ion ES/MS. For both chlorinated and nonchlorinated organic solvents, the higher the solution dielectric constant, the more the charge-state distribution is shifted toward higher charge states. These observations build on the “electrophoretic” mechanism of droplet charging. Solvents with high solution dielectric constants are considered to be most effective at stabilizing multiply charged ions (where charge separation is greatest), and they are likely to increase the level of droplet charging. Solvents with high basicities (gas phase and solution phase) and high proton affinities, yet low dielectric constants, favor lower charge states in ES mass spectra of lipid A and cardiolipin from Escherichia coli. This indicates that gas-phase processes and solvent basicity contribute much less toward ion formation than solution-phase solvation via preferred orientation of the solvent dipole.     

溶剂化效应对6-亚甲基环戊二烯酮与HCN的反应机理影响的理论研究

在B3LYP/6-311+G**计算水平上, 采用导体极化连续模型研究了溶剂化效应对6-亚甲基环戊二烯酮与HCN反应生成主要产物b类酸的反应机理的影响. 计算结果表明, 在溶剂中的反应机理与在气相中的反应机理一致. 溶剂化效应使反应路径中各驻点的自由能降低, 稳定化了各物质. 溶液中的活化自由能与气相相比也有所降低, 反应更容易发生, 其中CC进攻方式的活化自由能降低得更多.     

Conformational preferences of N-methoxycarbonyl proline dipeptide

The conformational study on N-methoxycarbonyl-L-proline-N'-methylamide (Moc-Pro-NHMe, prolylcarbamate) is carried out using ab initio HF and density functional B3LYP methods with the self-consistent reaction field method in the gas phase and in solution (chloroform, acetonitrile, and water). The replacement of the N-acetyl group by the N-methoxycarbonyl group results in the changes in conformational preferences, populations for backbone and prolyl puckering, and barriers to cis-trans isomerization of the prolyl residue in the gas phase and in solution, although there are small changes in the geometry of the prolyl peptide bond and the torsion angles of backbone and prolyl ring. The cis population increases with the increase of solvent polarity, as found for Ac-Pro-NHMe (prolylamide), but it is amplified by 9% in the gas phase and about 17% in solution for prolylcarbamate compared with those for prolylamide. It is found that the cis-trans isomerization for prolylcarbamate proceeds through the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond in the gas phase and in solution, as found for prolylamide. However, the rotational barriers to the cis-trans isomerization for prolylcarbamate are calculated to be 3.7-4.7 kcal/mol lower than those of prolylamide in the gas phase and in solution, and are found to be less sensitive to the solvent polarity. The calculated rotational barriers for prolylcarbamate in chloroform and water are in good agreement with the observed values. The shorter hydrogen-bond distance between the prolyl nitrogen and the amide H (H(NHMe)) of the NHMe group, the decrease in electron overlap of the prolyl C-N bond, and the favorable electrostatic interaction between the ester oxygen and the amide H(NHMe) for the transition state seem to play a role in lowering the rotational barrier of prolylcarbamate. The smaller molecular dipole moments of the ground- and transition-state structures for prolylcarbamate in the gas phase and in solution seem to be one of factors to make the rotational barrier less sensitive to the solvent polarity. As the solvent polarity increases (i.e., from the gas phase to chloroform to acetonitrile), the value of DeltaH(tc)(double dagger) decreases and the magnitude of DeltaS(tc)(double dagger) increases for prolylcarbamate, which results in a nearly constant value of the rotational barrier.     

How does halogen bonding behave in solution? A theoretical study using implicit solvation model

A systematic study of halogen bonding interactions in gas phase and in solution was carried out by means of quantum chemical DFT/B3LYP method. Three solvents with different polarities (chloroform, acetone, and water) were selected, and solvation effects were considered using the polarized continuum model (PCM). For charged halogen-bonded complexes, the strength of the interactions tends to significantly weaken in solution, with a concomitant elongation of intermolecular distances. For neutral systems, halogen bond distances are shown to shorten and the interaction energies change slightly. Computations also reveal that in the gas phase the binding affinities decrease in the order Cl(-) > Br(-) > I(-), while in solution the energy gaps of binding appear limited for the three halide anions. According to free energy results, many systems under investigation are stable in solution. Particularly, calculated free energies of formation of the complexes correlate well with halogen-bonding association constants determined experimentally. The differences of the effects of solvent upon halogen and hydrogen bonding were also elucidated. This study can establish fundamental characteristics of halogen bonding in media, which would be very helpful for applying this noncovalent interaction in medicinal chemistry and material design.     

Kinetics of the solvolysis of the cis-chloroazidobis-(1,2-diaminoethane)cobalt(III) ion in water+t-butanol mixtures

Summary Rate constants for the solvolysis of the complex ioncis-[Coen₂N₃Cl]⁺ have been determined for a range of temperatures in mixtures of water witht-butanol ranging up to 50% v/v of the latter. Linear plots of log (rate constant) against the reciprocal of the absolute temperature are obtained at all concentrations oft-butanol and the variations of the enthalpy and entropy of activation with solvent composition show extrema at compositions where extrema occur in the physical properties of the mixture which are influenced by solvent structure. The importance of solvent structure is shown by the curved plot found for log (rate constant) against reciprocal of dielectric constant at constant temperature. The application of a free energy cycle to the free energies of activation in water and in the mixture and the free energies of transfer of individual ionic species between water and the mixtures shows that the effect of changes in solvent structure on the reacting cation in the transition state dominates over that in the initial state.