key: cord-330025-6vql2u75 authors: Belova, Natalya V.; Girichev, Georgiy V.; Kotova, Vitaliya E.; Korolkova, Kseniya A.; Trang, Nguyen Hoang title: The molecular structure of 4-methylpyridine-N-oxide: Gas-phase electron diffraction and quantum chemical calculations date: 2018-03-15 journal: Journal of Molecular Structure DOI: 10.1016/j.molstruc.2017.11.070 sha: doc_id: 330025 cord_uid: 6vql2u75 Abstract The molecular structure of 4-methylpiridine-N-oxide, 4-MePyO, has been studied by gas-phase electron diffraction monitored by mass spectrometry (GED/MS) and quantum chemical (DFT) calculations. Both, quantum chemistry and GED analyses resulted in C S molecular symmetry with the planar pyridine ring. Obtained molecular parameters confirm the hyperconjugation in the pyridine ring and the sp2 hybridization concept of the nitrogen and carbon atoms in the ring. The experimental geometric parameters are in a good agreement with the parameters for non-substituted N-oxide and reproduced very closely by DFT calculations. The presence of the electron-donating CH3 substituent in 4-MePyO leads to a decrease of the ipso-angle and to an increase of r(N→O) in comparison with the non-substituted PyO. Electron density distribution analysis has been performed in terms of natural bond orbitals (NBO) scheme. The nature of the semipolar N→O bond is discussed. The considerable interest to the heterocyclic compounds containing N-oxide group is caused by their pronounced biological activity [1] . Among the representatives of this class, there are compounds with bactericidal, analgesic, anticonvulsant activity, as well as ones with apoptotic, antifungal and antimicrobial activity [2, 3] . These properties allow to use N-oxides as inhibitor of HIV-1 reverse transcriptase [4, 5] , as antiviral agent against various SARS coronavirus strains [6] , and as antiadhesive and quorum-sensitive inhibitor [7] . In the literature one can find the evidence, that some complexes with pyridine-N-oxides are used in agriculture to regulate plant growth rates [8] . It was shown that pyridine-N-oxide derivatives are genotoxic [9] . The biochemical activity of N-oxides is usually considered to be the result of complexation with the metalloporphyrins in living organisms [10e12] . The reactivity of Noxides may vary due to different substituents in the ring [13] . The variation of the substituents gives wide possibilities to chemical modifications of N-oxides and allows to influence their complexing properties, and, as the result, their biological activity. To understand the complexation ability of N-oxides with different substituents the detailed information about the properties of free molecules is needed. The available literature data about molecular structure of Noxides can not be considered comprehensive. Thus, only pyridine-N-oxide and three of its para-substituted compounds have been studied by gas-electron diffraction (GED) [14, 15] . Analyzing the average values of bond lengths in the pyridine ring the authors of [15] have concluded that the aromaticity increases in the series: 4nitropyridine-N-oxide (4-NO 2 PyO), 4-chloropyridine-N-oxide (4-ClPyO), pyridine-N-oxide (PyO), 4-methylpyridine-N-oxide(4-MePyO). The crystal structures of PyO, 4-NO 2 PyO and 4-MePyO have been studied by X-ray crystallography [16, 17] . The quantum chemical study of the substituent effect on the properties of pyridine-N-oxides has been performed for ten compounds with different substituents [18] . The comparative analysis of structural data obtained by different methods suggests that the molecular parameters of N-oxides obtained by GED [15] look rather strange. Thus, the NeC bond lengths in 4-NO 2 PyO, 4-MePyO, 4-ClPyO recommended in Ref. [15] are even longer than r(CeC) in benzene ring, and seems to be overestimated compare to the calculated values and X-ray data. Moreover, the bond distance r(N/O) in 4-MePyO molecule recommended by authors [15] appears to be significantly overestimated and conflicting with the calculated and X-ray values. Furthermore, the accuracy of GED data obtained in Refs. [14, 15] cannot be considered today as satisfactory. The discrepancies in molecular structure of substituted N-oxides according to GED and other methods motivated us to perform a new gas-electron diffraction study of 4-methylpyridine-N-oxide. This molecule has been chosen because of the largest structural contradictions. Now more sophisticated methods for the analysis of GED data are available, including the possibility of quantum chemical results to be used in the interpretation of experimental data. This allows us to determine the structural parameters more precisely compared to the studies performed in 1982. The electron diffraction patterns and the mass spectra were recorded simultaneously using the techniques described previously [19, 20] . Two series of GED/MS experiments at two different nozzleto-plate distances were performed. The conditions of GED/MS experiments and the relative abundance of the characteristic ions of 4-MePyO are shown in Tables 1 and 2, respectively. The temperature of stainless steel effusion cell was measured by a W/Re-5/20 thermocouple, which was standardised using the melting points of Sn and Al. Polycrystalline ZnO was used to establish the wavelength of the fast electrons. Microphotometric measurements were carried out by means of automatic techniques [32] . The molecular intensities were obtained in the ranges 2.5÷28.7 Å À1 (short camera) and 1.3÷15.5 Å À1 (long camera). The molecular intensities and the radial distributing curves are shown in Figs. 2 and 3 respectively. Quantum chemical calculations were carried out with GAUSSIAN 03 program set [21] . A hybrid DFT computational methods, namely B3LYP [22e26] and PBE0 [27] , as well as MP2 method, were used. Structure optimizations were followed by calculations of the vibrational frequencies in order to ensure that a minimum on the potential energy hyper-surface had been reached. Analysis of the potential functions of methyl group internal rotation shows that there is only one stable geometric configuration with H5C6C1C2 about 90 (see Fig. 1 for atom numbering). Although it should be noted that the internal rotation barrier is negligibly small, 0.007 kcal/mol (B3LYP/cc-pVTZ) or 0.03 kcal/mol (MP2/cc-pVTZ), thus the rotation of CH 3 group in 4-MePyO can be considered as free. The geometrical parameters of the calculated equilibrium structure are listed in Table 3 . Vibrational amplitudes and corrections, Dr ¼ r h1 À r a , were derived from theoretical force fields (B3LYP/cc-pVTZ) by Sipachev's method (approximation with taking into account the nonlinear kinematic effects at the level of the first order perturbation theory for the transformation of Cartesian coordinates into internal coordinates), using the program SHRINK [28e30] . Relevant values are listed in Table 4 (excluding nonbonded distances involving hydrogen). The NBO 5G program [31] , implemented for natural orbital analysis in PC GAMESS [32] , was used to obtain the net atomic charges, and to study the effect of hyperconjugation on the structure. B3LYP/aug-cc-pVTZ wave functions were used in the NBO analyses. Relevant values of net atomic charges, Wiberg bond indexes, and second order interaction energies (E (2) ) between donoracceptor orbitals are collected in Tables 6 and 7 . Visualization of the natural orbitals performed by the ChemCraft program [33] . The heaviest ion in the mass spectrum observed during the GED/ MS experiment was [Me-PyO] þ (Table 2) . No ions were detected, which could arise from impurities. This proves the monomeric vapour composition under the conditions of the GED experiment. The theoretical sM(s) functions were calculated with the following assumptions. Independent r h1 parameters were used to describe the molecular structure. Starting parameters from B3LYP/ cc-pVTZ calculations were refined. The differences between all CeC, as well as between all CeH bond lengths were constrained to calculated values (B3LYP/cc-pVTZ). In compliance with the quantum chemical calculation results, a planar skeleton with C 2v symmetry of the pyridine ring was assumed. The methyl group was considered in the framework of C 3v symmetry. Vibrational amplitudes were refined in groups with fixed differences. With the abovementioned assumptions, four bond distances and four bond angles (Table 3) were refined simultaneously with six groups of vibrational amplitudes. Vibrational amplitudes were refined in groups with fixed differences (Table 4 ). Only two correlation coefficients had absolute values larger than 0.7: (r(NeC3))/ (r(C2eC3)) ¼ À0.74, and (:NC3C2)/(:C3NC4) ¼ -0.84. The best agreement factor is R f ¼ 3.4%. Final results of the least squares analysis are given in Table 3 (geometric parameters) and Table 4 (vibrational amplitudes). Table 3 summarizes the structural parameters of 4-Me-PyO. We can note a perfect agreement between calculated and experimental molecular structure parameters. The molecule possesses almost C S symmetry (the torsional angle H5C6C1C2 is close to 90 ) with the planar pyridine ring. Table 5 compares gas-phase structures of PyO and 4-Me-PyO obtained in Refs. [14, 15] and in this work. The parameters obtained in our research is in a good agreement with the theoretical results (therein and [18] ) and, also, with the parameters for non-substituted N-oxide in contrast to the values recommended by Chiang J. F. et al. [15] . Furthermore, the parameters of free molecules of PyO [14] and 4-Me-PyO obtained in this work are in a good agreement with the crystal structures [34, 35] . Apparently, some increase of the semipolar N/O bond lengths in the crystal phase compared to free molecules is due to intermolecular hydrogen bonding. All parameters presented in Table 5 confirm the hyperconjugation in the pyridine ring and the sp2 hybridization concept of the nitrogen and carbon atoms in the ring. The exceptions are only the data by Chiang J. F. et al. [15] . which look rather strange. Unfortunately, the correct comparison of GED parameters of 4-MePyO and non-substituted pyridine-N-oxide looks to be impossible because of large uncertainties in GED values obtained for PyO [14] . According to quantum chemical calculations the insertion of the methyl group in para e position to the PyO has no significant effect on the structural parameters of the pyridine ring, excluding only the ipso-angle : C2C1C5, which changes from 120.6 in PyO to 115.9 in 4-MePyO (B3LYP/cc-pVTZ). It is interesting to note, that the substituent effects on the geometry of heterocyclic ring are rather similar to the tendencies for the benzene ring [36, 37] . Thus, the presence of the electron-donating CH 3 substituent leads to the decrease of the ipso-angle. Furthermore, according to the calculations the r(N/O) bond lengths appear to be also sensitive to a Fig. 1 . b Uncertainties in r h1 s¼(s sc 2 þ(2.5s LS ) 2 ) 1/2 (s sc ¼ 0,002r, s LS estandard deviation in least-squares refinement), for angles s ¼ 3s LS . c p i e parameter refined independently. Table 6 . The CeC and CeN bond orders in the ring for both molecules, as well as the values of the distances from Tables 3 and 5 show that these bonds do not correspond to either single nor double bonds. This confirms the existence of p-conjugation in the pyridine rings. The C1eC6 bond order for 4-MePyO corresponds to the single CeC bond. It is interesting to note that Q(NeO) in both molecules can not be interpreted as exact single or double bond. Moreover, the introduction of an electron-donating substituent increases the bond length N/O, slightly reduces the corresponding bond order and significantly increases the electron density on the oxygen. One can conclude that in this case the nucleophilic properties of the molecule increase. Furthermore, as it has been mentioned previously [18] the introduction of electron-donating substituents to the heterocyclic ring of N-oxides results in increasing ability for complex formation. Our experimental and calculated data (Tables 3, 6 and 7) reveal the incorrectness of the authors of [15] in their conclusion about more single NeO bond character in 4-MePyO. The delocalization of electron density between occupied Lewis type (bonds or lone pairs) NBO orbitals and formally unoccupied (antibonding) non-Lewis NBO orbitals corresponds to a stabilizing donor-acceptor interaction. The energy of these interactions can be estimated by the second order perturbation theory. Table 7 collects the largest values of calculated second order interaction energies (E 2 ) between donor-acceptor orbitals in PyO and 4-MePyO. In both cases several values correspond to the interactions related to the resonance in the pyridine ring, i.e. the interactions between p(CeC) and p*(CeC) and between p(CeC) and p*(CeN). The nature of semipolar NeO bond was explained in detail [18] . The largest second order interaction energies (E 2 ) in both molecules corre- To study the substituent effect on the aromaticity the nuclear independent chemical shifts NICS(1) have been calculated by using B3LYP/aug-cc-pVTZ approximation. NICS(1) values correspond to the negative isotropic shielding (in ppm) at 1 Å above the ring plane. The calculated values of NICS(1): À7.40 (PyO) and À7.36 (4-Me-PyO) are very close to each other. Thus, from the NICS(1) values one can conclude that the introduction of CH 3 group does not noticeably affect the aromaticity of the pyridine-N-oxide. Furthermore, the data of Table 5 do not confirm the conclusion of the authors of [15] about changing the aromaticity with substitution, which was made by analyzing the average values of bond length in the pyridine ring. Ab initio SCF and CI study of the electronic spectrum of pyridine Noxide 1-and 3-nitro-6-azabenzo[a]pyrens and their N-oxides: hifly mutagenic nitrated azaarenes Pridine oxide derivatives: structureactivity relationship for inhibition of human immunodeficiency virus and cytomegalovirus replication in cell culture Pyridine N-oxide derivatives are inhibitory to the human SARS and feline infectious peritonitis coronavirus in cell culture Aromatic Amine Oxides Genotoxic Impurities,strategies for Identification and Control Purification and properties of trimethylamine N-Oxide reductase from aerobic photosynthetic bacterium roseobacter denitrificans Non-enzymatic reduction of aliphatic tertiary amine N-oxides mediated by the haem moiety of cytochrome P450 Nonenzymatic reduction of brucine N-oxide by the heme group of ctocrome P450 Electronic factors effects on the reactivity of heteroaromatic Noxides Molecular structure of pyridine N oxide Molecular structures of 4-nitro-, 4-methyl-and 4-chloropyridine-N-oxides Substituent effect on the properties of pyridine-N-oxides Upgrading the EMR-100 Electrondiffraction Camera for Use with Gases Apparatus for Study of Molecular Structure of Valence-unsaturated Compounds, Instruments and Experimental Techniques Cr2 revisited Density-functional thermochemistry. III. The role of exact exchange Density-functional exchange-energy approximation with correct asymptotic behavior Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis Toward reliable density functional methods without adjustable parameters: the PBE0 model Local centrifugal distortions caused by internal motions of molecules Calculation of shrikage corrections in harmonic aproximation Advances in Molecular Structure Research ChemCraft version 1.6 (build 312) Crystal structure of 4-methylpyridine Noxide Isothermal and isochoric crystallization of highly hygroscopic pyridine N-oxide of aqueous solution Molecular geometry of substituted benzene derivatives. I. On the nature of the ring deformations indused by sublimations Group Electronegativities from benzene ring deformations: a quantum chemical study