An ab initio study of amide proton shift tensor dependence on local protein structure.

Ab initio shielding tensor calculations were carried out on residues in human ubiquitin. Reported experimental data on isotropic and anisotropic components of the amide proton chemical shifts were used as benchmarks to test the validity of the chosen basis sets as well as methods in structure optimization and shielding calculations. The best agreement with the experimental values was observed when the 6-311**G and 6-311++G(2d,2p) basis sets were used to optimize the structure and to calculate the shielding tensor, respectively. The same method was employed in subsequent model calculations to characterize the dependence of amide proton shielding to the local structure. Both the isotropic and the anisotropic components of the symmetric tensor were found to depend very strongly on the hydrogen bond length. A weaker dependence can also be observed for the hydrogen bond angle. Antisymmetric tensor elements were found to be relatively small. This study permits separation of various local structure contributions to the amide proton shielding tensor that complements scarce experimental data.     

33S NMR spectroscopy 3. Substituent effects on 33S NMR parameters in 2-substituted ethanesulfonates

33S NMR parameters (chemical shifts and linewidths) in 2-substituted sodium ethanesulfonates, XCH2CH2SO3Na (X = H, CH3, OH, SH, NH2, Cl, Br, NH3+) depend upon the electronic properties of substituents. To explain experimental results and obtain additional information on the origin of the observed substituent effect (SE), sulfur isotropic absolute shielding constants have been calculated at DFT level of theory (B3LYP/6-311++G(2d,p)) by gauge-including atomic orbitals (GIAO) method. Data have been interpreted with the aid of natural bond orbital (NBO) method and natural chemical shielding (NCS) analysis. It has been demonstrated that in the class of compounds considered the diamagnetic contribution to sulfur-shielding constant is constant and the observed upfield shift of 33S resonance induced by electron-withdrawing substituents (reverse chemical shift effect) can be related to variations of the paramagnetic contribution. Substituents with different electronic properties cause variations in the polarization of S-C and S-O bonds of the -C-SO3- moiety thus determining changes of the electron density at sulfur nucleus and consequently the expansion or contraction of 3p sulfur orbitals. Also oxygen lone-pairs and sulfur core 2p electrons can play an active role in determining the paramagnetic contribution to sulfur shielding. With regard to linewidth variations, they can be ascribed primarily to changes in the nuclear quadrupole coupling constant values. B3LYP/6-311++G(2d,p) method allows obtaining a good reproducibility of SE on the electric field gradient (EFG) at sulfur, although its values tend to be underestimated significantly. Moreover, 17O shielding constants have been calculated.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Planar tetracoordinate carbon species involving beryllium substituents

Planar tetracoordinate carbon (ptC) arrangements can be achieved by employing multiple substituents based on beryllium, despite its rather weak pi-acceptor ability. A variety of ptC-containing examples, some with more than one ptC, have been designed computationally by elaborating the planar C(BeH) 4 (2-) prototype at B3LYP/6-311++G(3df,2p) and MP2/6-311++G(3df,2p) levels of theory for some small ptC representatives. The prototype prefers a D(2h) paramagnetic triplet ground state due to Hund's rule, rather than a singlet. The highly polarized C-Be bonding weakens the rigidity of the tetrahedral carbon in T(d)C(BeH) 4 enormously, and the enhancement of both C-Be and Be 4 peripheral covalent bonding exerted by the extra electrons stabilizes the ptC eventually. The delocalization of the two p pi electrons is only modest, but their density on the most electronegative carbon atom helps stabilize the ptC arrangement. This is in contrast to the conventional strategy to delocalize p(pi) lone pairs for stabilizing the ptC arrangement. Various strategies to achieve neutral derivatives with ptCs are demonstrated.     

Intramolecular hydrogen bond in the push–pull CF3-aminoenones: DFT and FTIR study,NBO analysis

Postulated conformers of trifluoromethylated β-aminoenones stabilized by intramolecular NH?O and N?HO bonds were studied by IR and NMR spectroscopy and evaluated with quantum chemical calculations (B3LYP/6-311+G(d,p), MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) and MP2/6-31G(d,p)) and NBO analysis. The influence of the nature of EWG, substituents at the nitrogen atom and double bond, and of orbital interactions of heteroatoms and double bonds in these structures on the proton affinity of basic and acid centers, strength of hydrogen bonds, and the energy of tautomeric transfers is discussed. The theoretical results agree satisfactorily with the experimental observations.     

Electrostatic potential minimum of the aromatic ring as a measure of substituent constant

The molecular electrostatic potential minimum (Vmin) observed for the arene pi system of a substituted benzene derivatives is found to correlate linearly with the substituent constant sigma(p) degrees . The use of Vmin as a measure of substituent effect is further confirmed by obtaining a linear correlation between Vmin and a thermodynamic measure of the substituent effect obtained from an isodesmic reaction scheme involving benzene derivatives. Vmin and the recently proposed electrostatic potential value at the nucleus of the para carbon atom (Vc) show a nearly identical trend toward quantification of substituent effects. Both quantities have been compared at three different density functional theory methods, viz. B3LYP/6-311+G(2d,2p), BPW91/6-311G(d,p), and B3LYP/aug-cc-pvtz, as well as the at the MP2/6-31+G(d,p) level of theory, showing remarkable consistency among them.     

Molecular structures of benzoic acid and 2-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of benzoic acid (C6H5COOH) and 2-hydroxybenzoic acid (C6H4OHCOOH) have been determined in the gas phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d,p)) predict two conformers for benzoic acid, one which is 25.0 kJ mol(-1) (MP2) lower in energy than the other. In the low-energy form, the carboxyl group is coplanar with the phenyl ring and the O-H group eclipses the C=O bond. Theoretical calculations (HF and MP2/6-311+G(d,p)) carried out for 2-hydroxybenzoic acid gave evidence for seven stable conformers but one low-energy form (11.7 kJ mol(-1) lower in energy (MP2)) which again has the carboxyl group coplanar with the phenyl ring, the O-H of the carboxyl group eclipsing the C=O bond and the C=O of the carboxyl group oriented toward the O-H group of the phenyl ring. The effects of internal hydrogen bonding in 2-hydroxybenzoic acid can be clearly observed by comparison of pertinent structural parameters between the two compounds. These differences for 2-hydroxybenzoic acid include a shorter exocyclic C-C bond, a lengthening of the ring C-C bond between the substituents, and a shortening of the carboxylic single C-O bond.     

Noncovalent interactions between cytosine and SWCNT: curvature dependence of complexes via pi...pi stacking and cooperative CH...pi/NH...pi

Fragments of C24H12, adapted from a variety of armchair [(n,n), (n = 5, 7, and 8)] and zigzag [(m,0) (m = 8, 10, and 12)] single-walled carbon nanotube (SWCNT), are used to model corresponding SWCNTs with different diameters and electronic structures. The parallel binding mainly through pi...pi stacking interaction, as well as the perpendicular binding via cooperative NH...pi and CH...pi between cytosine and the fragments of SWCNT have been extensively investigated with a GGA type of DFT, PW91LYP/6-311++G(d,p). The eclipsed tangential (ET) conformation with respect to the six-membered ring of cytosine and the central ring of SWCNT fragments is less stable than the slipped tangential (ST) conformation for the given fragment; perpendicular conformations with NH2 and CH ends have higher negative binding energy than those with NH and CH ends. At PW91LYP/6-311++G(d,p) level, two tangential complexes are less bound than perpendicular complexes. However, as electron correlation is treated with MP2/6-311G(d,p) for PW91LYP/6-311++G(d,p) optimized complexes, it turns out there is an opposite trend that two tangential complexes become more stable than three perpendicular complexes. This result implies that electron correlation, a primary source to dispersion energy, has more significant contributions to the pi...pi stacking complexes than to the complexes via cooperative NH...pi and CH...pi interactions. In addition, it was found for the first time that binding energies for two tangential complexes become more negative with increasing nanotube diameter, while those for three perpendicular complexes have a weaker dependence on the curvature; i.e., binding energies are slightly less and less negative. The performance of a novel hybrid DFT, MPWB1K, was also discussed.     

Arene-cation interactions of positive quadrupole moment aromatics and arene-anion interactions of negative quadrupole moment aromatics

Intermolecular interactions involving aromatic pi-electron density are widely believed to be governed by the aromatic molecular quadrupole moment, Theta(zz). Arene-cation binding is believed to occur primarily with negative Theta(zz) aromatics, and arene-anion binding is believed to occur largely with positive Theta(zz) aromatics. We have performed quantum mechanical computations that show the cation binding of positive Theta(zz) aromatics and the anion binding of negative Theta(zz) aromatics is quite common in the gas phase. The pi-electron density of hexafluorobenzene, the prototypical positive Theta(zz ) aromatic (experimental Theta(zz) = 9.5 +/- 0.5 DA), has a Li+ binding enthalpy of -4.37 kcal/mol at the MP2(full)/6-311G**level of theory. The RHF/6-311G** calculated Theta(zz) value of 1,4-dicyanobenzene is +11.81 DA, yet it has an MP2(full)/6-311G** Li+ binding enthalpy of -12.65 kcal/mol and a Na+ binding enthalpy of -3.72 kcal/mol. The pi-electron density of benzene, the prototypical negative Theta(zz) aromatic (experimental Theta(zz) = -8.7 +/- 0.5 DA), has a F- binding enthalpy of -5.51 kcal/mol. The RHF/6-311G** calculated Theta(zz) of C6H2I4 is -10.45 DA, yet it has an MP2(full)/6-311++G** calculated F- binding enthalpy of -20.13 kcal/mol. Our results show that as the aromatic Theta(zz) value increases the cation binding enthalpy decreases; a plot of cation binding enthalpies versus aromatic Theta(zz) gives a line of best of fit with R2 = 0.778. No such correlation exists between the aromatic Theta(zz) value and the anion binding enthalpy; the line of best fit has R2 = 0.297. These results are discussed in terms of electrostatic and polarizability contributions to the overall binding enthalpies.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1,S2,S3和S4).其中,在复合物S1和S3中,HOCl单体的5H原子作为质子供体,与HCOCl单体中作为质子受体的10原子相互作用,形成红移氢键复合物;在复合物S4中,HOCl单体的7Cl原子作为质子供体,与HCOCl单体中作为质子受体的IO原子相互作用,形成红移卤键复合物;而在复合物S2中,同时存在2C-3H…6O蓝移氢键和4Cl…5O相互作用.在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正,其值在-5.05与-14.76 kJ·mol-1之间.采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

Spectra and structure of silicon-containing compounds. Part XXXVIII: Infrared and Raman spectra,vibrational assignment,conformational stability,and ab initio calculations of vinyldifluorosilane

Infrared spectra (3500-50 cm(-1)) of gaseous and solid, and Raman spectrum (3500-30 cm(-1)) of liquid vinyldifluorosilane, CH(2)z.dbnd6;CHSiF(2)H, are reported. Both the cis and gauche rotamers have been identified in the fluid phases. From temperature-dependent FT-infrared spectra of krypton solutions, it is shown that the cis conformer is more stable than the gauche form by 119+/-12 cm(-1) (1.42+/-0.14 kJ mol(-1)). At ambient temperature there is 53+/-2% of the gauche conformer present. Complete vibrational assignments are provided for the cis conformer and several modes are identified for the gauche form. Harmonic force constants, fundamental frequencies, infrared intensities, and Raman activities have been obtained from MP2/6-31G(d) calculations with full electron correlation. The optimized geometries and conformational stabilities have also been obtained from ab initio MP2/6-31G(d), MP2/6-311+G(d,p), and MP2/6-311+G(2d,2p) calculations with full electron correlation as well as from density functional theory calculations (DFT) by the B3LYP method. The SiH bond distances (r(0)) of 1.472 and 1.471 A have been obtained for the cis and gauche conformers, respectively, from the silicon-hydrogen stretching frequencies. These results are compared to the corresponding quantities of the corresponding carbon analogue as well as with some similar molecules.     

The Cycloaddition Reactions of Angle Strained Cycloalkynes. A Theoretical Study

The potential energy surfaces for the cycloaddition reactions of angle strained cycloalkynes to ethylene have been studied using ab initio methods. All the stationary points were determined with the MP2/6–311G (d,p) method with some calculations performed at the CCSD(T)/6–311G++G(d,p)//MP2/6–311G (d,p) level. Three kinds of cycloalkyne species, including monocyclic alkynes, bridged bicyclic alkynes, and heterocyclic alkynes, have been chosen in this work as model reactants. Two different reaction pathways have been proposed: (A) 1,2‐carbon shift and (B) 1,2‐hydrogen shift. That is, reactants → [2+1]‐TS‐1 → spiro‐carbene intermediate → (A) TS‐A → Pro‐A or (B) TS‐B → Pro‐B. As a result, it is found that ground‐state cycloalkyne appears to react more like a monocarbene than like an alkyne or a vicinal dicarbene as conventionally proposed. Our theoretical investigations also suggest that a cycloalkyne with a small C‐C≡C bond angle should be a good candidate for cycloaddition to an olefin. Moreover, in the cycloaddition reaction of a small (≤ six‐membered ring six‐membered) ring cycloalkyne, both 1,2‐carbon and 1,2‐hydrogen migrations will compete with each other. On the other hand, reactions involving larger (≥ seven‐membered) ring cycloalkynes should proceed with a 1,2‐carbon shift, leading to the major [2+2] cycloadduct. Furthermore, a configuration mixing model has been used to rationalize the computational results and to develop an explanation for the barrier heights. The results obtained allow a number of predictions to be made.     

Infrared and Raman spectra, conformational stability, ab initio calculations and vibrational assignments for trans-3-chloropropenoyl chloride

The infrared (3500-30 cm(-1)) spectra of gaseous and solid and the Raman (3500-200 cm(-1)) spectra of the liquid with quantitative depolarization ratios and solid trans-3-chloropropenoyl chloride (trans-ClCHCHCClO) have been recorded. These data indicate that both the anti (carbonyl bond trans to the carbon-carbon double bond) and syn conformers are present in the fluid states but only the anti conformer is present in the crystalline state. The mid-infrared spectra of the sample dissolved in liquid xenon as a function of temperature (-55 to -100 degrees C) have been recorded. Utilizing conformer pairs at 870 and 725 cm(-1), 1215 and 1029 cm(-1), and 1215 and 1228 cm(-1), the enthalpy difference has been determined to be 136+/-5 cm(-1) (389+/-14 cal mol(-1)) with the anti conformer the more stable form. Optimized geometries and conformational stabilities were obtained from ab initio calculations at the levels of RHF/6-31G(d), MP2/6-31G(d), MP2/6-311 + + G(d,p), MP2/6-311 + + G(2d,2p) and MP2/6-311 + + G(2df,2pd) with only the latter two calculations predicting the anti rotamer to be the more stable form. The vibrational frequencies, harmonic force constants and infrared intensities were obtained from the MP2/6-31G(d) calculations, whereas the Raman activities and depolarization values were obtained from the RHF/6-31G(d) calculations. The spectra are interpreted in detail and the results are compared with those obtained for some related molecules.     

33S NMR shieldings and chemical bonding in compounds of sulfur

Sulfur nuclear magnetic resonance (NMR) chemical shieldings have been determined at the correlation-including density functional theory scaled B3LYP/6-311+G(nd,p)//B3LYP/6-311+G(d,p) and modified MP2/6-311+G(nd,p) estimated infinite order Møller-Plesset levels with n = 2 for sulfur. The calculations span the range of sulfur shieldings and show agreement with experiment of about 3% of the shielding range. The atoms-in-molecules delocalization index and a covalent bond order from specific localized orbitals in the DFT approach are used to characterize sulfur's bonding and to relate it, where possible, to the calculated shieldings. © 2004 Wiley Periodicals, Inc. Heteroatom Chem 15:216–224, 2004; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/hc.20000     

Structures and Electronic Properties of HOCl···HCOCl Complexes

B3LYP/6-311++G** and MP2/6-311++G** calculations were used to analyze the interaction between hypochlorous acid (HOCl) and formyl chloride (HCOCl). The results showed that there were four equilibrium geometries (S1, S2, S3, and S4) optimized at B3LYP/6-311++G** level, and all the equilibrium geometries were confirmed to be in stable states by analytical frequency calculations. Complexes S1 and S3 use the 5H atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift hydrogen bond systems. However, the blue-shifted hydrogen bond (2C-3H···6O) coexists with 4Cl···5O interaction in structures S2. As for S4, it uses the 7Cl atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift halogen bond system. Interaction energies between monomers in the four complexes corrected with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) lie in the range from −5.05 to −14.76 kJ·mol⁻¹ at MP2/6-311++G** level. The natural bond orbital (NBO) and atoms in molecules (AIM) theories have also been applied to explain the structures and the properties of the complexes.     

Saturated hydrocarbon-benzene complexes: theoretical study of cooperative CH/pi interactions

High-level ab initio calculations at the CCSD(T)/aug-cc-pVTZ//MP2/aug(d,p)-6-311G(d,p) level were employed to investigate the cooperative CH/pi effects between the pi face of benzene and several modeled saturated hydrocarbons, propane, isobutane, cyclopropane, cyclobutane, cyclopentane, cyclohexane, cyclopentane, cyclooctane, and bicyclo[2.2.2]octane. In all cases, multiple C-H groups (2-4) are found to interact with the pi face of benzene, with one C-H group pointing close to the center of the benzene ring. The geometries of these complexes are governed predominantly by electrostatic interaction between the interacting systems. The calculated interaction energies (10-14 kJ mol(-1)) are 2-3 times larger than that of the prototypical methane-benzene complex. The trends of geometries, interaction energies, binding properties, as well as electron-density topological properties were analyzed. The calculated interaction energies correlate well with the polarizabilities of the hydrocarbons. AIM analysis confirms the hydrogen-bonded nature of the CH/pi interactions. Significant changes in proton chemical shift and stretching frequency (blue shift) are predicted for the ring C-H bond in these complexes.     

Theoretical Study of Hydrogen Bonding Interaction of 1:1 Dimer of HNO with HArF

The hydrogen bonding interaction of 1:1 dimer formed between HNO and HArF molecule has been completely investigated in the present study using Second-order M?ller-Plesset Perturbation (MP2) method in conjunction with 6-311+G**, 6-311++G** and 6-311++G(2d,2p) basis sets. The standard and CP-corrected calculations have been employed to determine the equilibrium structures, the vibrational frequencies and interaction energies. The interaction energies of the dimers were also calculated at G2MP2 level. Two stable structures are found as the minima. Dimer I(H···F)is a five-membered cyclic hydrogen bonded structure and is more stable than the Dimer II(H···O). The blue-shifted N-H···F hydrogen bond is confirmed with standard and CP-corrected calculations by the MP2 and DFT methods in conjunction with different basis sets. The results obtained at MP2 in conjunction with different basis sets show there is a red-shifted hydrogen bond (Ar-H···O) in the Dimer II(H···O). The topological and electronic properties, the origin of red- and blue-shifted hydrogen bonds were investigated at MP2/6-311++G(2d,2p) with CP corrected calculations. From the NBO analysis, the reasonable explanations for the red- and blue-shifted hydrogen bonds were proposed.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

Theoretical studies of azapentalenes. Part 4: Theoretical study of the properties of 3a,6a-diazapentalene

A theoretical study of the properties of the isolated 3a,6a-diazapentalene by means of DFT, B3LYP/6-311++G(d,p) and ab initio methods, MP2/6-311++G(d,p), has been carried out. In addition, the complexes formed with hydrogen bond donor, acceptors, cations, and anions have been studied and analyzed. Ring opening into 1,5-diazocine as well as basicity and acidity properties of 3a,6a-diazapentalene have been explored. Their ability to form HB complexes and the complexes formed with anions and cations have been studied.     

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.