An Ab Initio Investigation of the Structure and Alkaline Earth Divalent Cation Selectivity of 18-Crown-6

Infrared multiple photon dissociation action spectroscopy of protonated glycine, histidine, lysine, and arginine complexed with 18-crown-6 ether

Complexes of 18-crown-6 ether (18C6) with four protonated amino acids (AAs) are examined using infrared multiple photon dissociation (IRMPD) action spectroscopy utilizing light generated by the infrared free electron laser at the Centre Laser Infrarouge d'Orsay (CLIO). The AAs examined in this work include glycine (Gly) and the three basic AAs: histidine (His), lysine (Lys), and arginine (Arg). To identify the (AA)H⁺(18C6) conformations present in the experimental studies, the measured IRMPD spectra are compared to spectra calculated at the B3LYP/6-311+G(d,p) level of theory. Relative energies of various conformers and isomers are provided by single point energy calculations carried out at the B3LYP, B3P86, M06, and MP2(full) levels using the 6-311+G(2p,2d) basis set. The comparisons between the IRMPD and theoretical IR spectra indicate that 18C6 binds to Gly and His via the protonated backbone amino group, whereas protonated Lys prefers binding via the protonated side-chain amino group. Results for Arg are less definitive with strong evidence for binding to the protonated guanidino side chain (the calculated ground conformer at most levels of theory), but contributions from backbone binding to a zwitterionic structure are likely.

Aerosol Cross-Linked Crown Ether Diols Melded with Poly(vinyl alcohol) as Specialized Microfibrous Li+ Adsorbents

Crown ether (CE)-based Li⁺ adsorbent microfibers (MFs) were successfully fabricated through a combined use of CE diols, electrospinning, and aerosol cross-linking. The 14- to 16-membered CEs, with varied ring subunits and cavity dimensions, have two hydroxyl groups for covalent attachments to poly(vinyl alcohol) (PVA) as the chosen matrix. The CE diols were blended with PVA and transformed into microfibers via electrospinning, a highly effective technique in minimizing CE loss during MF fabrication. Subsequent aerosol glutaraldehyde (GA) cross-linking of the electrospun CE/PVA MFs stabilized the adsorbents in water. The aerosol technique is highly effective in cross-linking the MFs at short time (5 h) with minimal volume requirement of GA solution (2.4 mL g^-1 MF). GA cross-linking alleviated CE leakage from the fibers as the CEs were securely attached with PVA through covalent CE-GA-PVA linkages. Three types of CE/PVA MFs were fabricated and characterized through Fourier transform infrared-attenuated total reflection, ¹³C cross-polarization magic angle spinning NMR, field emission scanning electron microscope, N₂ adsorption/desorption, and universal testing machine. The MFs exhibited pseudo-second-order rate and Langmuir-type Li⁺ adsorption. At their saturated states, the MFs were able to use 90-99% CEs for 1:1 Li⁺ complexation, suggesting favorability of their microfibrous structures for CE accessibility to Li⁺. The MFs were highly Li⁺-selective in seawater. Neopentyl-bearing CE was most effective in blocking larger monovalents Na⁺ and K⁺, whereas the dibenzo CE was best in discriminating divalents Mg²⁺ and Ca²⁺. Experimental selectivity trends concur with the reaction enthalpies from density functional theory calculations, confirming the influence of CE structures and cavity dimensions in their "size-match" Li⁺ selectivity.

“Leaching or not leaching”: an alternative approach to antimicrobial materials via copolymers containing crown ethers as active groups

New copolymers containing MMA and 18C6 crown-ether pendants, with or without a PEG arm, were synthesized to check if sequestering structural alkali-earth ions from the bacterial outer membrane (OM) may lead to bacterial death.

What causes the weakest host to act as the strongest one? A theoretical study on the host-guest chemistry of five azacryptands and fluoride anions

In this work we have attempted to computationally analyze the important parameters affecting the selectivity in host-guest systems in order to show that the solvent effect in some host-guest systems is even more important than the hole-size fitting and/or host-guest interaction energy. For this purpose, the fluoride anion selectivities of the five most studied azacryptands with different affinities, moieties, cavity sizes and degrees of preorganization, in their fully protonated forms, were studied at PBE/TZVP and B3P86/TZVP levels of theory. The factors affecting the selectivity such as hole size fitting, steric effects, electronic properties, preorganization of hosts, as well as solvent effects were investigated. The results showed that among the five studied azacryptands, one of them, which ranks fifth for selectivity according to its weak interaction energy with fluoride anions, ranks fourth due to its good preorganization. However, a surprising improvement occurs when its selectivity for fluoride anions ranks first due to the existence of less solvent hindrance. The results of the calculations were further confirmed by a good correlation between the calculated and experimental formation constants.

Remarkable conformational flexibility of aqueous 18-crown-6 and its strontium(II) complex-ab initio molecular dynamics simulations

Ab initio QMCF-MD simulations of aqueous 18-crown-6 (18C6) and strontium(II)-18-crown-6 (18C6-Sr) were performed to gather insight into their hydration properties. Strongly different characteristics were found for the two solutes in terms of structure and dynamics such as H-bonding. They, however, have in common that their backbone shows high flexibility in aqueous medium, adopting structures significantly differing from idealized gas phase geometries. In particular, planar oxyethylene units stable in the picosecond range occurred in 18C6, while the strontium complex readily exhibits a bent structure. Detailed analysis of this high flexibility was done via two dimensional root mean square deviation plots as well as the evolution of dihedral angles and angles within the simulation trajectory. The vibrational spectra obtained from the QMCF-MD simulations are in excellent agreement with experimental data and show a pronounced blueshift upon complexation of the strontium(II) ion in 18C6.

Complexation of some alkali and alkaline earth metal cations by macrocyclic compounds containing four pyridine subunits – a DFT study

Tetradentate pyridine-based macrocyclic compounds offer useful ligands capable of efficient and selective complexation of M = Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺and Ca²⁺.

Conductometric Studies of Thermodynamics of 1,10-Didecyl-1,10-diaza-18-crown-6 Complexes with , , , , , , and Ions in Acetonitrile, Methanol, and Ethanol Solutions

The complexation reactions between 1,10-didecyl-1,10-diaza-18-crown-6 (DD18C6) and , , , , , and ions were studied conductometrically in acetonitrile, methanol, and ethanol at various temperatures. The formation constants of the resulting 1 : 1 complexes were calculated from the computer fitting of the molar conductance-mole ratio data at different temperatures. In the solvents studied, the stability of the resulting complexes varied in the order acetonitrile ≫ methanol > ethanol. The enthalpy and entropy changes of the complexation reactions were evaluated from the temperature dependence of formation constants. The versus plot of thermodynamic data obtained shows a linear correlation indicating the existence of enthalpy-entropy compensation in the complexation reactions. The ab initio studies calculated at B3LYP/6-31G level of theory indicate that the binding energy of complexes decreases with increasing cation size in the gas phase. The experimental data obtained in the solvents show different trends so that the stability constant of Li+and Mg2+complexes is lower than other alkali and alkaline earth ions, respectively due to the solvation of these small ions in solution.

Ab initio and density functional theoretical design and screening of model crown ether based ligand (host) for extraction of lithium metal ion (guest): effect of donor and electronic induction

The structures, energetic and thermodynamic parameters of model crown ethers with different donor, cavity and electron donating/ withdrawing functional group have been determined with ab initio MP2 and density functional theory in gas and solvent phase. The calculated values of binding energy/ enthalpy for lithium ion complexation are marginally higher for hard donor based aza and oxa crown compared to soft donor based thia and phospha crown. The calculated values of binding enthalpy for lithium metal ion with 12C4 at MP2 level of theory is in good agreement with the available experimental result. The binding energy is altered due to the inductive effect imparted by the electron donating/ withdrawing group in crown ether, which is well correlated with the values of electron transfer. The role of entropy for extraction of hydrated lithium metal ion by different donor and functional group based ligand has been demonstrated. The HOMO-LUMO gap is decreased and dipole moment of the ligand is increased from gas phase to organic phase because of the dielectric constant of the solvent. The gas phase binding energy is reduced in solvent phase as the solvent molecules weaken the metal-ligand binding. The theoretical values of extraction energy for LiCl salt from aqueous solution in different organic solvent is validated by the experimental trend. The study presented here should contribute to the design of model host ligand and screening of solvent for metal ion recognition and thus can contribute in planning the experiments.

Computational modeling of capillary electrophoretic behavior of primary amines using dual system of 18-crown-6 and β-cyclodextrin

Using capillary electrophoresis (CE) three chiral primary amine compounds 1-aminoindan (AI), 1-(1-naphthyl)ethylamine (NEA) and 1,2,3,4-tetrahydro-1-naphthylamine (THAN), exhibited only partial or no separation when β-cyclodextrin (βCD) was used as chiral selector. The use of 18-crown-6 (18C6) as a second additive with βCD resulted in an enhanced separation. A molecular modeling study, using molecular mechanics and the semiempirical PM6 calculations, was used to help explaining the mechanism of the enantiodifferentiation and to predict the separation process. Optimization of the structures of the complexes by the PM6 method indicate that the poor separation obtained in the presence of the βCD chiral selector alone is due to the small binding energy differences (ΔΔE) of 4.7, 1.1 and 1.2 kcal mol(-1) for AI, NEA and THAN, respectively. In the presence of 18C6 it was suggested that a sandwich compound between 18C6, amine and βCD is formed. Theoretical calculations show that a significant increase in the binding energy is obtained for the sandwich compounds indicating strong hydrophobic and van der Waals interactions that show enhanced enantiodifferentiation.

Ultraviolet photodepletion spectroscopy of dibenzo-18-crown-6-ether complexes with alkaline Earth metal divalent cations

Ultraviolet photodepletion spectra of dibenzo-18-crown-6-ether complexes with alkaline earth metal divalent cations (A(2+)-DB18C6, A = Ba, Sr, Ca, and Mg) were obtained in the gas phase using electrospray ionization quadrupole ion-trap reflectron time-of-flight mass spectrometry. Each spectrum exhibits the lowest energy absorption band in the wavenumber region of 35 400-37 800 cm(-1), which is tentatively assigned as the origin of the S(0)-S(1) transition of A(2+)-DB18C6. This origin band shows a red shift as the size of the metal dication increases from Mg(2+) to Ba(2+). The binding energies of the metal dications to DB18C6 at the S(0) state were calculated at the lowest energy structures optimized by the density functional theory and employed with the experimental energies of the origin bands to estimate the binding energies at the S(1) state. We suggest that the red shifts of the origin bands arise from the decrease in the binding energies of the metal dications at the S(1) state by nearly constant ratios with respect to the binding energies at the S(0) state, which decrease with increasing size of the metal dication. This unique relationship of the binding energies between the S(0) and S(1) states gives rise to a linear correlation between the relative shift of the origin band of A(2+)-DB18C6 and the binding energy of the metal dication at the S(0) state. The size effects of the metal cations on the properties of metal-DB18C6 complex ions are also manifested in the linear plot of the relative shift of the origin band as a function of the size to charge ratio of the metal cations, where the shifts of the origin bands for all DB18C6 complexes with alkali and alkaline earth metal cations are fit to the same line.

DFT modeling on the suitable crown ether architecture for complexation with Cs⁺ and Sr²⁺ metal ions

Crown ether architectures were explored for the inclusion of Cs(+) and Sr(2+) ions within nano-cavity of macrocyclic crown ethers using density functional theory (DFT) modeling. The modeling was undertaken to gain insight into the mechanism of the complexation of Cs(+) and Sr(2+) ion with this ligand experimentally. The selectivity of Cs(+) and Sr(2+) ions for a particular size of crown ether has been explained based on the fitting and binding interaction of the guest ions in the narrow cavity of crown ethers. Although, Di-Benzo-18-Crown-6 (DB18C6) and Di-Benzo-21-Crown-7 (DB21C7) provide suitable host architecture for Sr(2+) and Cs(+) ions respectively as the ion size match with the cavity of the host, but consideration of binding interaction along with the cavity matching both DB18C6 and DB21C7 prefers Sr(2+) ion. The calculated values of binding enthalpy of Cs metal ion with the crown ethers were found to be in good agreement with the experimental results. The gas phase binding enthalpy for Sr(2+) ion with crown ether was higher than Cs metal ion. The ion exchange reaction between Sr and Cs always favors the selection of Sr metal ion both in the gas and in micro-solvated systems. The gas phase selectivity remains unchanged in micro-solvated phase. We have demonstrated the effect of micro-solvation on the binding interaction between the metal ions (Cs(+) and Sr(2+)) and the macrocyclic crown ethers by considering micro-solvated metal ions up to eight water molecules directly attached to the metal ion and also by considering two water molecules attached to metal-ion-crown ether complexes. A metal ion exchange reaction involving the replacement of strontium ion in metal ion-crown ether complexes with cesium ion contained within a metal ion-water cluster serves as the basis for modeling binding preferences in solution. The calculated O-H stretching frequency of H(2)O molecule in micro-solvated metal ion-crown complexes is more red-shifted in comparison to hydrated metal ions. The calculated IR spectra can be compared with an experimental spectrum to determine the presence of micro-solvated metal ion-crown ether complexes in extractant phase.

Density functional theory study of the free and tetraprotonated spheroidal macrotricyclic ligands and the complexes with halide anions: F(-), Cl(-), Br(-)

Theoretical studies of the macrotricyclic tetramine hexaether (SC), its tetraprotonated form SC-4H(+), and the corresponding complexes X(-) subsetSC-4H(+) (This expression represents the structural properties of the halide inclusion complex formed though the free ligand SC-4H(+) and the halide anion X(-): the spherical halide anion X(-) is held by a tetrahedral array of (+)N--H...X(-) hydrogen bonds inside the intramolecular cavity of the tetraprotonated form SC-4H(+)) of SC-4H(+) with the halide anions: F(-), Cl(-), and Br(-) have been performed using density functional theory (DFT) with B3LYP/6-31G method implemented in the Gaussian 03 program package. The optimized geometric structures obtained from DFT calculations are used to perform Natural Bond Orbital (NBO) analysis. The three main types of hydrogen bonds (+)N--H...F(-), (+)N--H...Cl(-), and (+)N--H...Br(-) are investigated. The results indicate that hydrogen bonding interactions are dominant and the halide anions: F(-), Cl(-), and Br(-) offer lone pair electrons to the contacting sigma* (N--H) antibond orbital of SC-4H(+). For all the structures, the most pronounced changes in geometric parameters upon interaction are observed in the proton-donor molecule. The intermolecular interaction energies are predicted by using B3LYP/6-31G methods with basis set superposition error (BSSE) and zero-point energy (ZPE) correction.

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.

Conformational study of the structure of free 18-crown-6

A conformational search was performed for 18-crown-6 using the CONLEX method at the MM3 level. To have a more accurate energy order of the predicted conformations, the predicted conformations were geometry optimized at the HF/STO-3G level and the 198 lowest energy conformations, according to the HF/STO-3G energy order, were geometry optimized at the HF/6-31+G level. In addition, the 47 nonredundant lowest energy conformations, according to the MP2/6-31+G energy order at the HF/6-31+G optimized geometry, hereafter the MP2/6-31+G//HF/6-31+G energy order, were geometry optimized at the B3LYP/6-31+G level. According to the MP2/6-31+G//B3LYP/6-31+G energy order, three conformations had energies lower than the experimentally known Ci conformation of 18c6. At the MP2/6-31+G//B3LYP/6-31+G level, the S6 lowest energy conformation is more stable by 1.96 kcal/mol than this Ci conformation. This was confirmed by results at the MP2/6-31+G level with an energy difference of 1.84 kcal/mol. Comparison between the structure of the S6 conformation of 18c6 and the S4 lowest energy conformation of 12-crown-4, as well as other important conformations of both molecules, is made. It is concluded that the correlation energy is necessary to have an accurate energy order of the predicted conformations. A rationalization of the conformational energy order in terms of the hydrogen bonding and conformational dihedral angles is given. It is also suggested that to have a better energy order of the predicted conformations at the MM3 level, better empirical force fields corresponding to the hydrogen bond interactions are needed.

Experimental and Theoretical Study of the Vibrational Spectra of 12-Crown-4−Alkali Metal Cation Complexes

The vibrational, Raman, and IR, spectra of the five 12-crown-4 (12c4) complexes with Li+, Na+, K+, Rb+, and Cs+ alkali metal cations were measured. Except for a small shift of the position of some bands in the vibrational spectra of the Li+ complex, the vibrational spectra of the five complexes are so similar that it is concluded that the five complexes exist in the same conformation. B3LYP/6-31+G* force fields were calculated for six of the eight predicted conformations in a previous report (J. Phys. Chem. A 2005, 109, 8041) of the 12c4-Li+, Na+, and K+ complexes that are of symmetries higher than the C1 symmetry. These six conformations, in energy order, are of C4, Cs, Cs, C(2v), C(2v), and Cs symmetries. Comparison between the experimental and calculated vibrational frequencies assuming any of the above-mentioned six conformations shows that the five complexes exist in the C4 conformation. This agrees with the fact that the five alkali metal cations are larger than the 12c4 ring cavity. The B3LYP/6-31+G* force fields of the C4 conformation of the Li+, Na+ and K+ complexes were scaled using a set of eight scale factors and the scale factors were varied so as to minimize the difference between the calculated and experimental vibrational frequencies. The root-mean-square (rms) deviations of the calculated frequencies from the experimental frequencies were 7.7, 5.6, and 5.1 cm(-1) for the Li+, Na+, and K+ complexes, respectively. To account for the earlier results of the Li+ complex that the Cs conformation is more stable than the C4 conformation by 0.16 kcal/mol at the MP2/6-31+G* level, optimized geometries of the complex were calculated for the C4 and Cs conformations at the MP2/6-311++G** level. The C4 conformation was calculated to be more stable than the Cs conformation by 0.13 kcal/mol.

Determination of the Conformation of 2-Hydroxy- and 2-Aminobenzoic Acid Dimers Using 13C NMR and Density Functional Theory/Natural Bond Order Analysis: The Central Importance of the Carboxylic Acid Carbon

The 13C chemical shift for the carboxylic acid carbon provides a powerful diagnostic probe to determine the preferred isomeric dimer structures of benzoic acid derivatives undergoing intra- and intermolecular H-bonding in the gas, solution and crystalline phases. We have employed hybrid density functional calculations and natural bond orbital analysis to elucidate the electronic origins of the observed 13C shieldings and their relationship to isomeric stability. We find that delocalizing interactions from the carbonyl oxygen lone pairs (nO) into vicinal carbon-oxygen and carbon-carbon antibonds (sigmaCO*,sigmaCC*) make critical contributions to the 13C shieldings, and these nO --> sigmaCO*, nO --> sigmaCC* interactions are in turn sensitive to the intramolecular interactions that dictate dimer structure and stability. The carboxyl carbon atom can thus serve as a useful detector of subtle structural and conformational features in this pharmacologically important class of carboxylic acid interactions.

Conformational Stability of a Model Macrocycle Tetraamide: An ab Initio Study

Ab initio calculations are carried out to investigate the conformational stability of a model macrocyle tetraamide. The four amide groups in the selected model are present in the sequence: -(O=CNH)-Ph-(NHC=O)-CH=CH-(O=CNH)-Ph-(NHC=O)-CH=CH-. In this sequence, two phenyl rings and two ethene groups act as bridges between the amide units. Each amide motif bonds to a phenyl ring through its amide nitrogen and to an ethene group through its amide carbon. Four clearly distinct minimum-energy conformations are found upon full geometry optimization using the B3LYP/6-31+G(d) method. Frequency calculations using the same method confirm that the four conformations are indeed minima in the macrocycle potential energy surface. Relative to the most stable conformer, the other conformations are higher in energy by 0.86, 2.09, and 9.17 kcal/mol, respectively, at the MP2/6-31+G(d,p) level. The stability of the macrocycle conformations is correlated primarily to the existence and strength of intramolecular N-H...O=C hydrogen bonds. Additional stability to the conformations is found to come from weak Ph-H...O=C hydrogen bonding between a carbonyl oxygen and a hydrogen atom of a phenyl group. Solvent effects play an important role in the relative energies of the various conformations, as indicated by the simple SCRF = dipole model calculations for the case of aqueous solution.