The interaction of formohydroxamic acid with carbon monoxide: FTIR matrix isolation and quantum chemistry studies

Argon matrix infrared spectra of the complexes formed between formohydroxamic acid (HCONHOH) and carbon monoxide have been recorded. The experimental results indicate formation of three isomeric complexes. In two complexes the NH or OH groups of formohydroxamic acid are attached to the carbon atom of the CO molecule and in the third complex the NH group interacts with the oxygen atom of CO. The formohydroxamic acid complex with two CO molecules is also trapped in the matrix. One of the two CO molecules interacts with the NH group and the second one with the OH group of HCONHOH, in both cases the site of interaction is the carbon atom of CO. Theoretical studies of the structure and spectral characteristics of the complexes were carried out on the DFT(B3LYP)/6-311++G(2d,2p) level. The calculated vibrational frequencies for the complexes present in the matrices are in good agreement with the experimental data. The calculations show also an additional potential energy minimum corresponding to the complex in which the OH group of formohydroxamic acid is attached to the oxygen atom of carbon monoxide.     

Theoretical investigation of hydrogen bonds between CO and HNF2, H2NF, and HNO

Ab initio quantum mechanics methods were applied to investigate the hydrogen bonds between CO and HNF2, H2NF, and HNO. We use the Hartree-Fock, MP2, and MP4(SDQ) theories with three basis sets 6-311++G(d,p), 6-311++G(2df,2p), and AUG-cc-pVDZ, and both the standard gradient and counterpoise-corrected gradient techniques to optimize the geometries in order to explore the effects of the theories, basis sets, and different optimization methods on this type of H bond. Eight complexes are obtained, including the two types of C...H-N and O...H-N hydrogen bonds: OC...HNF2(C(s)), OC...H2NF(C(s) and C1), and OC...HNO(C(s)), and CO...HNF2(C(s)), CO...H2NF(C(s) and C1), and CO...HNO(C(s)). The vibrational analysis shows that they have no imaginary frequencies and are minima in potential energy surfaces. The N-H bonds exhibit a small decrease with a concomitant blue shift of the N-H stretch frequency on complexation, except for OC...HNF2 and OC...H2NF(C1), which are red-shifting at high levels of theory and with large basis sets. The O...H-N hydrogen bonds are very weak, with 0 K dissociation energies of only 0.2-2.5 kJ/mol, but the C...H-N hydrogen bonds are stronger with dissociation energies of 2.7-7.0 kJ/mol at the MP2/AUG-cc-pVDZ level. It is notable that the IR intensity of the N-H stretch vibration decreases on complexation for the proton donor HNO but increases for HNF2 and H2NF. A calculation investigation of the dipole moment derivative leads to the conclusion that a negative permanent dipole moment derivative of the proton donor is not a necessary condition for the formation of the blue-shifting hydrogen bond. Natural bond orbital analysis shows that for the C...H-N hydrogen bonds a large electron density is transferred from CO to the donors, but for the O...H-N hydrogen bonds a small electron density transfer exists from the proton donor to the acceptor CO, which is unusual except for CO...H2NF(C(s)). From the fact that the bent hydrogen bonds in OC(CO)...H2NF(C(s)) are quite different from those in the others, we conclude that a greatly bent H-bond configuration shall inhibit both hyperconjugation and rehybridization.     

Theoretical study of the vibrational spectra of the hydrogen-bonded systems between pyridine-3-carboxamide (nicotinamide) and DMSO

The vibrational characteristics (vibrational frequencies and infrared intensities) for the hydrogen-bonded systems of nicotinamide (NA(Z) and NA(E)) with dimethyl sulfoxide (DMSO) have been predicted using ab initio SCF/6-31G(d,p) and DFT (BLYP/6-311++G(d,p)) calculations. The changes in the vibrational characteristics from free monomers to a complex have been calculated. The ab initio and BLYP calculations show that the complexation between nicotinamide (NA(Z) and NA(E)) and DMSO leads to large red shifts of the stretching vibrations for the hydrogen-bonded N-H bonds of nicotinamide and very strong increase in their IR intensity. The results from the BLYP/6-311++G(d,p) calculations show that the predicted red shifts of the nu(s)(NH) and nu(as)(NH) vibrations for the complex NA(E)-DMSO (1:2) (Deltanu(as)(NH)=-186 cm(-1) and Deltanu(s)(NH)=-198 cm(-1)) are in better agreement with the experimentally measured. The magnitudes of the wavenumber shifts are indicative of strong NH...O hydrogen-bonded interactions in both complexes. The calculations predict an increase of the IR intensity of nu(s)(NH) and nu(as)(NH) vibrations in the complexes up to 14 times. Having in mind that in more cases the predicted changes in the vibrational characteristics for the complexes studied are very near, it could be concluded that both conformers of nicotinamide, Z-conformer and E-conformer, are present in the solution forming the hydrogen-bonded complexes with DMSO.     

An ab initio correlated study of the potential energy surface for the HOBr.H2O complex

The potential energy surface (PES) for the HOBr.H(2)O complex has been investigated using second- and fourth-order M?ller-Plesset perturbation theory (MP2, MP4) and coupled cluster theory with single and doubles excitations (CCSD), and a perturbative approximation of triple excitations (CCSD-T), correlated ab initio levels of theory employing basis sets of triple zeta quality with polarization and diffuse functions up to the 6-311++G(3dp,3df ) standard Pople's basis set. Six stationary points being three minima, two first-order transition state (TS) structures and one second-order TS were located on the PES. The global minimum syn and the anti equilibrium structure are virtually degenerated [DeltaE(ele-nuc) approximately 0.3 kcal mol(-1), CCSD-T/6-311++G(3df,3pd) value], with the third minima being approximately 4 kcal mol(-1) away. IRC analysis was performed to confirm the correct connectivity of the two first-order TS structures. The CCSD-T/6-311++G(3df,3pd)//MP2/6-311G(d,p) barrier for the syn<-->anti interconversion is 0.3 kcal mol(-1), indicating that a mixture of the syn and anti forms of the HOBr.H(2)O complex is likely to exist.     

Searching blue-shifted O–H···Y hydrogen bond in CH<Subscript>3</Subscript>OH complexes with CF<Subscript>4</Subscript>, C<Subscript>2</Subscript>F<Subscript>2</Subscript>, OC,Ne, and He: a theoretical study

The hydrogen-bonded structures of the CH₃OH complexes with CF₄, C₂F₂, OC, Ne, and He are designated as the starting points for geometry optimizations without and with counterpoise (CP) correction at MP2 level of theory with the basis sets 6-31+G*, 6-31++G**, and 6-311++G**, respectively. Tight convergence criteria are applied throughout all geometry optimizations in order to reduce the computational errors. According to the optimizations without CP correction, a blue-shifted O–H···Y (where Y = F, O, Ne, or He) hydrogen bond exists in all these five complexes. The magnitudes of blue shifts of ν(O–H) of the former four complexes with respect to that of CH₃OH are reduced greatly when the polarization and diffuse functions of the hydrogen atoms are added (results from 6-31+G* versus those from 6-31++G**). However, for the complexes CH₃OH–CF₄ and CH₃OH–C₂F₂, our optimizations using the CP corrections did not find the hydrogen-bonded structure to be a stationary point. The energy minimum of both the complexes corresponds to a non-hydrogen-bonded structure.     

Hydrogen-bonded complexes of nitrous and nitric acids with ethene Ab initio and DFT studies

The complexes formed by ethene with nitric and nitrous (trans and cis) acids have been investigated by ab initio (SCF and MP2) and B3LYP calculations with 6-311++G(d,p) basis set. Full geometry optimisation has been performed for the complexes studied. The most stable structures of the complexes are established. Bearing in mind the corrected values of the dissociation energy the studied hydrogen-bonded complexes can be ordered as follows: C(2)H(4)...HONO(2)>C(2)H(4)...HONO-trans>C(2)H(4)...HONO-cis. In the complexes the acids act as proton donors forming the pi-type of hydrogen bond with ethene. The predicted changes in the vibrational characteristics (vibrational frequencies and infrared intensities) arising from the hydrogen bonding are in good agreement with the experimentally measured. The predicted frequency shift of the stretching OH vibration in the nitric acid is largest (-210 cm(-1)), followed by the shifts in the trans-HONO (-141 cm(-1)) and cis-HONO (-109 cm(-1)). The calculations predict an increase of the IR intensity of the stretching O-H vibration in the complexes from 6 to 10 times.     

Vibrational frequencies and infrared intensities of the hydrogen-bonded complexes of nitrous acid with ethers: ab initio and DFT studies

The vibrational spectra of the binary complexes formed by HONO-trans and HONO-cis with dimethyl and diethyl ethers have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimisation was made for the complexes studied. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The ab initio and DFT calculations show that the changes in the vibrational characteristics (vibrational frequencies and infrared intensities) upon hydrogen bonding for the hydrogen-bonded complex (CH3)2O...HONO-trans are larger than for the complex (CH3)2O...HONO-cis.     

Influence of H2 on the gas-phase decomposition of formic acid: a theoretical study

Gas-phase decomposition of formic acid results in final products CO + H2O and CO2 + H2. Experimentally, the CO/CO2 ratio tends to be large, in contradiction with mechanism studies, which show almost equal activation energies for dehydration and decarboxylation. In this work, the influence of H2 on the decomposition mechanism of HCOOH was explored using ab initio calculations at the CCSD(T)/6-311++G**//MP2/6-311++G** level. It was found that, in the presence of H2, the reaction channels leading to CO + H2O are more than those leading to CO2 + H2. With competitive energy, H2 addition to HCOOH can reduce the latter into HCHO, which then dissociates into CO + H2 catalyzed by H2O. Compared to trans-HCOOH, cis-HCOOH and cis-C(OH)2, conformers required for decarboxylation, are less populated due to interactions with H2.     

Structures of [(CO2)n(CH3OH)m](-) (n = 1-4, m = 1, 2) cluster anions

The infrared photodissociation spectra of [(CO 2) n (CH 3OH) m ] (-) ( n = 1-4, m = 1, 2) are measured in the 2700-3700 cm (-1) range. The observed spectra consist of an intense broad band characteristic of hydrogen-bonded OH stretching vibrations at approximately 3300 cm (-1) and congested vibrational bands around 2900 cm (-1). No photofragment signal is observed for [(CO 2) 1,2(CH 3OH) 1] (-) in the spectral range studied. Ab initio calculations are performed at the MP2/6-311++G** level to obtain structural information such as optimized structures, stabilization energies, and vibrational frequencies of [(CO 2) n (CH 3OH) m ] (-). Comparison between the experimental and the theoretical results reveals the structural properties of [(CO 2) n (CH 3OH) m ] (-): (1) the incorporated CH 3OH interacts directly with either CO 2 (-) or C 2O 4 (-) core by forming an O-HO linkage; (2) the introduction of CH 3OH promotes charge localization in the clusters via the hydrogen-bond formation, resulting in the predominance of CO 2 (-).(CH 3OH) m (CO 2) n-1 isomeric forms over C 2O 4 (-).(CH 3OH) m (CO 2) n-2 ; (3) the hydroxyl group of CH 3OH provides an additional solvation cite for neutral CO 2 molecules.     

Theoretical study of the structures, stability and vibrational spectra of the nitrous acid complexes with CH4

The structures, stability and vibrational spectra of the binary complexes CH4...HONO-trans and CH4...HONO-cis have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimization was made for the complexes studied. It was established that the complex CH4...HONO-trans is more stable by 0.41 kcal mol(-1) than the complex CH4...HONO-cis. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The changes in the vibrational characteristics of methane and trans-, cis-nitrous acid upon formation of the hydrogen bond show that the complexes CH4...HONO-trans and CH4...HONO-cis have geometry in which the OH group interacts with a methane molecule forming a single hydrogen bond. This fact is confirmed by relatively strong perturbation of the OH stretching vibration to lower frequencies and an increase of the infrared intensity of this vibration up to three times upon hydrogen bonding.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

Dimers of boroglycine and methylamine boronic acid: a computational comparison of the relative importance of dative versus hydrogen bonding

Boronic acids are widely used in materials science, pharmacology, and the synthesis of biologically active compounds. In this Article, geometrical structures and relative energies of dimers of boroglycine, H2N-CH2-B(OH)2, and its constitutional isomer H3C-NH-B(OH)2, were computed using second-order M?ller-Plesset perturbation theory and density functional theory; Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the MP2 calculations, and the Pople 6-311++G(d,p) basis set was employed for a majority of the DFT calculations. Effects of an aqueous environment were incorporated into the results using PCM and COSMO-RS methodology. The lowest-energy conformer of the H2N-CH2-B(OH)2 dimer was a six-membered ring structure (chair conformation; Ci symmetry) with two intermolecular B:N dative-bonds; it was 14.0 kcal/mol lower in energy at the MP2/aug-cc-pVDZ computational level than a conformer with the classic eight-centered ring structure (Ci symmetry) in which the boroglycine monomers are linked by a pair of H-O...H bonds. Compared to the results of MP2 calculations with correlation-consistent basis sets, DFT calculations using the PBE1PBE and TPSS functionals with the 6-311++G(d,p) basis set were significantly better at predicting relative conformational energies of the H2N-CH2-B(OH)2 and H3C-NH-B(OH)2 dimers than corresponding calculations using the BLYP, B3LYP, OLYP, and O3LYP functionals, particularly with respect to dative-bonded structures.     

Complexes of atmospheric alpha-dicarbonyls with water: FTIR matrix isolation and theoretical study

The complexes of glyoxal (Gly), methylglyoxal (MGly), and diacetyl (DAc) with water have been studied using Fourier transform infrared (FTIR) matrix isolation spectroscopy and MP2 calculations with 6-311++G(2d,2p) basis set. The analysis of the experimental spectra of the Gly(MGly,DAc)/H2O/Ar matrixes indicates formation of one Gly...H2O complex, three MGly...H2O complexes, and two DAc...H2O ones. All the complexes are stabilized by the O-H...O(C) hydrogen bond between the water molecule and carbonyl oxygen as evidenced by the strong perturbation of the O-H, C=O stretching vibrations. The blue shift of the CH stretching vibration in the Gly...H2O complex and in two MGly...H2O ones suggests that these complexes are additionally stabilized by the improper C-H...O(H2) hydrogen bonding. The theoretical calculations confirm the experimental findings. They evidence the stability of three hydrogen-bonded Gly...H2O and DAc...H2O complexes and six MGly...H2O ones stabilized by the O-H...O(C) hydrogen bond. The calculated vibrational frequencies and geometrical parameters indicate that one DAc..H2O complexes, two Gly...H2O, and three MGly...H2O ones are additionally stabilized by the improper hydrogen bonding between the C-H group and water oxygen. The comparison of the theoretical frequencies with the experimental ones allowed us to attribute the calculated structures to the complexes present in the matrixes.     

Ab initio comprehensive conformational analysis of 2'-deoxyuridine, the biologically significant DNA minor nucleoside, and reconstruction of its low-temperature matrix infrared spectrum

A comprehensive conformational analysis of isolated 2'-deoxyuridine (dU), a minor DNA nucleoside, has been performed by means of ab initio calculations at the MP2/6-311++G (d,p)//DFT B3LYP/6-31G (d,p) level of theory. At 298.15 and 420 K, all 94 allowed conformers of dU are within 8.96 and 7.91 kcal/mol Gibbs energy ranges, respectively. Syn orientation for the base and South (S) conformers for the sugar dominate at 298.15 K: syn/anti=62.3%:37.7% and S/N=77.2%:22.8%. At 420 K in the majority of conformers, the base is anti oriented and the population of North (N) sugars increases: syn/anti=39.3%:60.7% and S/N=63.0%:37.0%. Values of all conformational parameters and correlations between them, as well as their correlations with valence bonds, and also correlations between valence bonds and angles were estimated. In general, 14 types of intramolecular H-bonds were detected (1-3 H-bonds per conformer, the total number 175), namely, C1'H...O2 (16 H-bonds), C2'H1...O5' (9), C2'H2...O2 (21), C3'H...O2 (21), C5'H1...O2 (14), C5'H2...O2 (11), C6H...O4' (37), C6H...O5' (22), C3'H...HC6 (4), O5'H...HC6 (2), O3'H...O5' (5), O5'H...O4' (1), O5'H...O3' (4), and O5'H...O2 (8). Geometric, vibrational, structural-topological, and energetic features of the OH...O intramolecular H-bonds in dU conformers were determined. The close similarity between energetic and geometric characteristics of dU and thymidine DNA-like conformers in anti and relevant syn conformations and their transition states of the anti-->syn interconversion implies that mismatch DNA glycosylase discriminates between the two nucleosides, mainly because of the difference in the shapes of their bases. Convolution of calculated IR spectra of all the dU conformers within the limits 3400-3700 cm(-1) appears to be consistent with its low-temperature matrix IR spectrum (Ivanov et al. Spectrochim. Acta, Part A 2003, 59, 1959), wavenumber discrepancy not exceeding 1%. It was concluded that, for a reliable reproduction of the experimental spectrum, the whole set of conformers should be taken into consideration. The suggested method makes reconstruction of the isolated nucleoside IR spectrum at a physiological interval of temperature reasonably possible.     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Observation of cysteine thiolate and -S...H-O intermolecular hydrogen bond

The cysteine anion was produced in the gas phase by electrospray ionization and investigated by photoelectron spectroscopy at low temperature (70 K). The cysteine anion was found to exhibit the thiolate form [-SCH2CH(NH2)CO2H], rather than the expected carboxylate form [HSCH2CH(NH2)CO2-]. This observation was confirmed by two control experiments, that is, methyl cysteine [CH3SCH2CH(NH2)CO2-] and cysteine methyl ester [-SCH2CH(NH2)CO2CH3]. The electron binding energy of [-SCH2CH(NH2)CO2H] was measured to be about 0.7 eV blue-shifted relative to [-SCH2CH(NH2)CO2CH3] due to the formation of an intramolecular -S-...HO2C- hydrogen bond in the cysteine thiolate. Theoretical calculations at the CCSD(T)/6-311++G(2df,p) and B3LYP/6-311++G(2df,p) levels were carried out to estimate the strength of this intramolecular -S-...HO2C- hydrogen bond. Combining experimental measurements and theoretical calculations yielded an estimated value of 16.4 +/- 2.0 kcal/mol for the -S-...HO2C- intramolecular hydrogen-bond strength.     

Properties of the halogen-hydride interaction: an ab initio and "atoms in molecules" analysis

X-Cl...H-Y interactions are analyzed by applying ab initio methods as well as the Bader theory. All calculations were performed using Pople's basis sets (6-311++G(2df,2pd) and 6-311++G(3df,3pd)) as well as the Dunning-type bases (aug-cc-pVDZ and aug-cc-pVTZ) within the MP2 method. For the complexes analyzed here, X-Cl and H-Y may be treated as a Lewis acid and a Lewis base, respectively. The Cl...H interactions are rather weak or at most moderate since, for the strongest interaction of the F3...HLi complex, the binding energy calculated at the MP2/6-311++G(3df,3pd) level of approximation amounts to -3.4 kcal/mol, and the H...Cl distance is equal to 2.65 A, less than the corresponding sum of van der Waals radii. These interactions may be classified as halogen-hydride interactions. However, some of the complexes analyzed, especially F3SiCl...HBeF and F3SiCl...HBeF, are very weakly bound, probably by typical van der Waals interactions.     

AB INITIO SCF MO STUDIES ON THE COMPLEXES BETWEEN CH_3OH AND H_2CO

<正> The interaction between CH3OH and H2CO has been studied by ab ini-tio method at the level of STO-3G and 6-311G basis sets. It has been found that there are two possible complexes; a hydrogen bonded complex CH3OH...CH2O(Ⅰ) and an electron donor-acceptor complex CH3OH.....OCH2(Ⅱ).The stabilization energies of (Ⅰ) and (Ⅱ) are 14. 6 and 3. 6kJ/mol (STO-3G results) or 25. 1 and 17. 1kJ/mol (6-311G results) respectively. The nature of these complexes has been discussed by using the energy decomposition scheme.     

Using Raman Spectroscopy and ab initio Calculations to Investigate Intermolecular Hydrogen Bonds in Binary Mixture (Tetrahydrofuran+Water)

We analyzed the properties and structures of the hydrogen-bonded complexes of tetrahydrofuran(THF) and water by means of experimental Raman spectra and ab initio calculations.The optimized geometries and vibrational frequencies of the neat THF molecule and its hydrogen-bonded complexes with water(THF/H2O) were calculated at the MP2/6-311+G(d,p) level of theory.We found that the intermolecular hydrogen bonds which are formed from the binary mixtures of the neat THF and water with different molar ratios could...     

How many conformers determine the thymidine low-temperature matrix infrared spectrum? DFT and MP2 quantum chemical study

A comprehensive conformational analysis of isolated 2'-beta-deoxy-thymidine (T), canonical DNA nucleoside, has been performed by means of ab initio calculations at the MP2/6-311++G(d,p)//DFT B3LYP/6-31G(d,p) level of theory. At 298.15 K, all 92 conformers of isolated dT are within a 7.49 kcal/mol Gibbs energy range. Syn orientation for the base and South (S) conformers for the sugar dominate at this temperature: syn/anti = 61.6%:38.4% and S/N = 74.5%:25.5%. However, at 420 K, the majority of conformers contain anti base and the population of North (N) sugars increases: syn/anti = 38.0%:62.0% and S/N = 59.5%:40.5%. The whole conformational parameters (P, chi, gamma, delta, beta, epsilon, nu max) were analyzed as well as the energies of the OH...O intramolecular H-bonds on the basis of nu(OH) stretching vibrations. Convolution of calculated IR spectra of all of the T conformers appears consistent with its low-temperature matrix spectrum (Ivanov et al. Low Temp. Phys. 2003, 29, 809). The maximal discrepancy in frequencies between calculated and experimental spectra is less than 1%. A conclusion was made that for reliable reconstruction of the isolated nucleoside IR spectrum the quasi whole set of conformers should be taken into consideration. In essence, this result opens up a possibility to reconstruct IR spectra of isolated nucleosides at physiological temperatures with rather satisfactory probability.