Biological Activity of Ionic Liquids and Their Application in Pharmaceutics and Medicine

Ionic liquids are remarkable chemical compounds, which find applications in many areas of modern science. Because of their highly tunable nature and exceptional properties, ionic liquids have become essential players in the fields of synthesis and catalysis, extraction, electrochemistry, analytics, biotechnology, etc. Apart from physical and chemical features of ionic liquids, their high biological activity has been attracting significant attention from biochemists, ecologists, and medical scientists. This Review is dedicated to biological activities of ionic liquids, with a special emphasis on their potential employment in pharmaceutics and medicine. The accumulated data on the biological activity of ionic liquids, including their antimicrobial and cytotoxic properties, are discussed in view of possible applications in drug synthesis and drug delivery systems. Dedicated attention is given to a novel active pharmaceutical ingredient-ionic liquid (API-IL) concept, which suggests using traditional drugs in the form of ionic liquid species. The main aim of this Review is to attract a broad audience of chemical, biological, and medical scientists to study advantages of ionic liquid pharmaceutics. Overall, the discussed data highlight the importance of the research direction defined as "Ioliomics", studies of ions in liquids in modern chemistry, biology, and medicine.

Density Functional Theory Analysis of Anthraquinone Derivative Hydrogenation over Palladium Catalyst

A density functional theory (DFT) analysis was conducted on the hydrogenation of 2-alkyl-anthraquinone (AQ), including 2-ethyl-9,10-anthraquinone (eAQ) and 2-ethyl-5,6,7,8-tetrahydro-9,10-anthraquinone (H₄ eAQ), to the corresponding anthrahydroquinone (AQH₂ ) over a Pd₆ H₂ cluster. Hydrogenation of H₄ eAQ is suggested to be more favorable than that of eAQ owing to a higher adsorption energy of the reactant (H₄ eAQ), lower barrier of activation energy, and smaller desorption energy of the target product (2-ethyl-5,6,7,8-tetrahydro-9,10-anthrahydroquinone, H₄ eAQH₂ ). For the most probable reaction routes, the energy barrier of the second hydrogenation step of AQ is circa 8 kcal mol^-1 higher than that of the first step. Electron transfer of these processes were systematically investigated. Facile electron transfer from Pd₆ H₂ cluster to AQ/AQH intermediate favors the hydrogenation of C=O. The electron delocalization over the boundary aromatic ring of AQ/AQH intermediate and the electron-withdrawing effect of C=O are responsible for the electron transfer. In addition, a pathway of the electron transfer is proposed for the adsorption and subsequent hydrogenation of AQ on the surface of Pd₆ H₂ cluster. The electron transfers from the abstracted H atom (reactive H) to a neighbor Pd atom (Pd_H ), and finally goes to the carbonyl group through the C₄ atom of AQ aromatic ring (C₄ ).

The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactions

The σ-hole and π-hole of the protonated 2-halogenated imidazolium cation (XC₃H₄N₂⁺; X = F, Cl, Br, I) were investigated and analyzed. The monomers of (CH₃)₃SiY(Y=F, Cl, Br, I), considered as the Lewis base, were combined with the σ-hole and π-hole of XC₃H₄N₂⁺ to form the σ-hole and π-hole interactions in the bimolecular complexes (CH₃)₃SiY · · · XC₃H₄N₂⁺ and (CH₃)₃SiY · · · C₃(X)H₄N₂⁺(X/Y=F, Cl, Br, I), respectively. For both the σ-hole and π-hole interactions, the equilibrium geometries of complexes show regular changes according to the sequence of heavy sequence of the noncovalent interaction acceptors and donors. The electrostatic energy is the main contribution in the formation of both kinds of interactions, it has linear relations with the V _S,max values of σ-hole and the V' _S,max values of π-hole. Both the σ-hole and π-hole interactions belong to the closed-shell and noncovalent interactions. The π-hole interactions are stronger than the σ-hole interactions. For the π-hole interactions, the contribution percents of the dispersion energies are somewhat greater than those of the σ-hole interactions, while it is contrary for the polarization energy. Graphical Abstract The protonated 2-halogenated imidazolium cation as the noncovalent interaction donor: the σ-hole and π-hole interactionsᅟ.

Incorporation of local structure into kriging models for the prediction of atomistic properties in the water decamer

Machine learning algorithms have been demonstrated to predict atomistic properties approaching the accuracy of quantum chemical calculations at significantly less computational cost. Difficulties arise, however, when attempting to apply these techniques to large systems, or systems possessing excessive conformational freedom. In this article, the machine learning method kriging is applied to predict both the intra-atomic and interatomic energies, as well as the electrostatic multipole moments, of the atoms of a water molecule at the center of a 10 water molecule (decamer) cluster. Unlike previous work, where the properties of small water clusters were predicted using a molecular local frame, and where training set inputs (features) were based on atomic index, a variety of feature definitions and coordinate frames are considered here to increase prediction accuracy. It is shown that, for a water molecule at the center of a decamer, no single method of defining features or coordinate schemes is optimal for every property. However, explicitly accounting for the structure of the first solvation shell in the definition of the features of the kriging training set, and centring the coordinate frame on the atom-of-interest will, in general, return better predictions than models that apply the standard methods of feature definition, or a molecular coordinate frame. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Association of symmetrical alkane diols with pyridine: DFT/GIAO calculation of 1 H NMR chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of the free diol and of its 1 : 1 and 1 : 2 hydrogen-bonded complexes with pyridine have been computed for five symmetrical alkane diols on the basis of density functional theory, by applying the gauge-including atomic orbital method to geometry-optimized conformers. For certain conformers, intramolecular OH···OH interactions, evidenced by high NMR OH proton shifts, are further enhanced on going from the free diol to the corresponding 1 : 1 diol/pyridine complex. This is confirmed by atoms-in-molecules and non-covalent interaction plots. The computed OH and CH proton shifts for the diol and the two complexes correlate well with values obtained by analysing data from the NMR titration of the diols in benzene against pyridine. Shift values for the diols in neat pyridine are calculated by weighting the shifts of the various protons in the three forms (free diol, 1 : 1 and 1 : 2 diol/pyridine complexes) according to the experimentally determined association constants. The results are in good agreement with those observed, and after empirical scaling, the root mean square difference is 0.18 ppm. Copyright © 2016 John Wiley & Sons, Ltd.

GenLocDip: A Generalized Program to Calculate and Visualize Local Electric Dipole Moments

Local dipole moments (i.e., dipole moments of atomic or molecular subsystems) are essential for understanding various phenomena in nanoscience, such as solvent effects on the conductance of single molecules in break junctions or the interaction between the tip and the adsorbate in atomic force microscopy. We introduce GenLocDip, a program for calculating and visualizing local dipole moments of molecular subsystems. GenLocDip currently uses the Atoms-In-Molecules (AIM) partitioning scheme and is interfaced to various AIM programs. This enables postprocessing of a variety of electronic structure output formats including cube and wavefunction files, and, in general, output from any other code capable of writing the electron density on a three-dimensional grid. It uses a modified version of Bader's and Laidig's approach for achieving origin-independence of local dipoles by referring to internal reference points which can (but do not need to be) bond critical points (BCPs). Furthermore, the code allows the export of critical points and local dipole moments into a POVray readable input format. It is particularly designed for fragments of large systems, for which no BCPs have been calculated for computational efficiency reasons, because large interfragment distances prevent their identification, or because a local partitioning scheme different from AIM was used. The program requires only minimal user input and is written in the Fortran90 programming language. To demonstrate the capabilities of the program, examples are given for covalently and non-covalently bound systems, in particular molecular adsorbates. © 2016 Wiley Periodicals, Inc.

Local electric dipole moments: A generalized approach

We present an approach for calculating local electric dipole moments for fragments of molecular or supramolecular systems. This is important for understanding chemical gating and solvent effects in nanoelectronics, atomic force microscopy, and intensities in infrared spectroscopy. Owing to the nonzero partial charge of most fragments, "naively" defined local dipole moments are origin-dependent. Inspired by previous work based on Bader's atoms-in-molecules (AIM) partitioning, we derive a definition of fragment dipole moments which achieves origin-independence by relying on internal reference points. Instead of bond critical points (BCPs) as in existing approaches, we use as few reference points as possible, which are located between the fragment and the remainder(s) of the system and may be chosen based on chemical intuition. This allows our approach to be used with AIM implementations that circumvent the calculation of critical points for reasons of computational efficiency, for cases where no BCPs are found due to large interfragment distances, and with local partitioning schemes other than AIM which do not provide BCPs. It is applicable to both covalently and noncovalently bound systems. © 2016 Wiley Periodicals, Inc.

Boron as an Electron-Pair Donor for B⋅⋅⋅Cl Halogen Bonds

MP2/aug'-cc-pVTZ calculations were performed to investigate boron as an electron-pair donor in halogen-bonded complexes (CO)₂ (HB):ClX and (N₂ )₂ (HB):ClX, for X=F, Cl, OH, NC, CN, CCH, CH₃ , and H. Equilibrium halogen-bonded complexes with boron as the electron-pair donor are found on all of the potential surfaces, except for (CO)₂ (HB):ClCH₃ and (N₂ )₂ (HB):ClF. The majority of these complexes are stabilized by traditional halogen bonds, except for (CO)₂ (HB):ClF, (CO)₂ (HB):ClCl, (N₂ )₂ (HB):ClCl, and (N₂ )₂ (HB):ClOH, which are stabilized by chlorine-shared halogen bonds. These complexes have increased binding energies and shorter B-Cl distances. Charge transfer stabilizes all complexes and occurs from the B lone pair to the σ* Cl-A orbital of ClX, in which A is the atom of X directly bonded to Cl. A second reduced charge-transfer interaction occurs in (CO)₂ (HB):ClX complexes from the Cl lone pair to the π* C≡O orbitals. Equation-of-motion coupled cluster singles and doubles (EOM-CCSD) spin-spin coupling constants, ^1x J(B-Cl), across the halogen bonds are also indicative of the changing nature of this bond. ^1x J(B-Cl) values for both series of complexes are positive at long distances, increase as the distance decreases, and then decrease as the halogen bonds change from traditional to chlorine-shared bonds, and begin to approach the values for the covalent bonds in the corresponding ions [(CO)₂ (HB)-Cl]⁺ and [(N₂ )₂ (HB)-Cl]⁺ . Changes in ¹¹ B chemical shieldings upon complexation correlate with changes in the charges on B.

Theoretical studies on structure, formation and nature of bond in a Disila-, Digerma- and distannacyclobutene ring

A theoretical study on the structure, formation and nature of E–E and C–E bonds (E[Formula: see text][Formula: see text][Formula: see text]Si, Ge, Sn) in catenated compounds of the group 14 elements including disila-, digerma- and distannacyclobutene ring formed from the reaction of 1-thiacyclohept-4-yne and E(NR)2SiR2 molecules [E[Formula: see text][Formula: see text][Formula: see text]Si, Ge, Sn, R[Formula: see text][Formula: see text][Formula: see text]t-Bu, H, F, Cl, Br] in 1:2 mole ratio has been investigated at the M06/def2-TZVPP level of theory. The results showed that the formation energy of the products of the reaction above with Si(NBr)2SiBr2 and Sn(NF)2SiF2 reactants has greatest and the smallest value, respectively. In agreement with the values of formation energies, both the calculated Wiberg bond indices (WBI) for E—E and C—E bonds and [Formula: see text](BCP) of E—E bond in the products of two above reactants have largest and smallest values, respectively. The nature of E—E bond in the products was also studied with atoms in molecules (AIM) and natural bond orbital (NBO) analyses. The data confirmed that the E—E bond is partly covalent. In addition, the nature of C—E bond was investigated with energy decomposition analysis (EDA) and it was shown that the covalent contribution is in the range 48–53% depending on the types of E atom and corresponding substituents.

Theoretical and Synthetic Study on the Existence, Structures, and Bonding of the Halide-Bridged [B2X7](-) (X = F, Cl, Br, I) Anions

While hydrogen bridging is very common in boron chemistry, halogen bridging is rather rare. The simplest halogen-bridged boron compounds are the [B2X7](-) anions (X = F, Cl, Br, I), of which only [B2F7](-) has been reported to exist experimentally. In this paper a detailed theoretical and synthetic study on the [B2X7](-) anions is presented. The structures of [B2X7](-) anions have been calculated at the MP2/def2-TZVPP level of theory, and their local minima have been shown to be of C2 symmetry in all cases. The bonding situation varies significantly between the different anions. While in [B2F7](-) the bonding is mainly governed by electrostatics, the charge is almost equally distributed over all atoms in [B2I7](-) and additional weak iodine···iodine interactions are observed. This was shown by an atoms in molecules (AIM) analysis. The thermodynamic stability of the [B2X7](-) anions was estimated in all phases (gas, solution, and solid state) based on quantum-chemical calculations and estimations of the lattice enthalpies using a volume-based approach. In the gas phase the formation of [B2X7](-) anions from [BX4](-) and BX3 is favored in accord with the high Lewis acidity of the BX3 molecules. In solution and in the solid state only [B2F7](-) is stable against dissociation. The other three anions are borderline cases, which might be detectable under favorable conditions. However, experimental attempts to identify [B2X7](-) (X = Cl, Br, I) anions in solution by (11)B NMR spectroscopy and to prepare stable [PNP][B2X7] salts failed.

Comparative studies on group III σ-hole and π-hole interactions

The σ-hole of M2 H6 (M = Al, Ga, In) and π-hole of MH3 (M = Al, Ga, In) were discovered and analyzed, the bimolecular complexes M2 H6 ···NH3 and MH3 ···N2 P2 F4 (M = Al, Ga, In) were constructed to carry out comparative studies on the group III σ-hole interactions and π-hole interactions. The two types of interactions are all partial-covalent interactions; the π-hole interactions are stronger than σ-hole interactions. The electrostatic energy is the largest contribution for forming the σ-hole and π-hole interaction, the polarization energy is also an important factor to form the M···N interaction. The electrostatic energy contributions to the interaction energy of the σ-hole interactions are somewhat greater than those of the π-hole interactions. However, the polarization contributions for the π-hole interactions are somewhat greater than those for the σ-hole interactions.

Electrostatic potentials and average electron densities of bioisosteres in methylsquarate and acetic acid

Background: The bioisosterism in −CO2H and −C4HO3 is exploited using the quantum theory of atoms in molecules and molecular electrostatic potentials (ESP). Results & discussion: Bioisosteres in methylsquarate and acetic acid, in the neutral/anionic forms, have average electron densities that differ by less than 2% (i.e., ∼0.01 atomic units) while irrespective of the capping group. The topography of the ESP reveals similarities in the case of the neutral species but not in the anionic forms. Conclusion: The nonclassical bioisosteres in methylsquarate and acetic acid have average electron densities that are similar and relatively insensitive to the ionization state (neutral or anionic) or its studied capping group (H, CH3, Cl or phenyl). The ESP reveals similarities in the topography of neutral molecules.

Muon-Substituted Malonaldehyde: Transforming a Transition State into a Stable Structure by Isotope Substitution

Isotope substitutions are usually conceived to play a marginal role on the structure and bonding pattern of molecules. However, a recent study [Angew. Chem. Int. Ed. 2014, 53, 13706-13709; Angew. Chem. 2014, 126, 13925-13929] further demonstrates that upon replacing a proton with a positively charged muon, as the lightest radioisotope of hydrogen, radical changes in the nature of the structure and bonding of certain species may take place. The present report is a primary attempt to introduce another example of structural transformation on the basis of the malonaldehyde system. Accordingly, upon replacing the proton between the two oxygen atoms of malonaldehyde with the positively charged muon a serious structural transformation is observed. By using the ab initio nuclear-electronic orbital non-Born-Oppenheimer procedure, the nuclear configuration of the muon-substituted species is derived. The resulting nuclear configuration is much more similar to the transition state of the proton transfer in malonaldehyde rather than to the stable configuration of malonaldehyde. The comparison of the "atoms in molecules" (AIM) structure of the muon-substituted malonaldehyde and the AIM structure of the stable and the transition-state configurations of malonaldehyde also unequivocally demonstrates substantial similarities of the muon-substituted malonaldehyde to the transition state.

Electron density, disorder and polymorphism: high-resolution diffraction studies of the highly polymorphic neuralgic drug carbamazepine

Analysis of neutron and high-resolution X-ray diffraction data on form (III) of carbamazepine at 100 K using the atoms in molecules (AIM) topological approach afforded excellent agreement between the experimental results and theoretical densities from the optimized gas-phase structure and from multipole modelling of static theoretical structure factors. The charge density analysis provides experimental confirmation of the partially localized π-bonding suggested by the conventional structural formula, but the evidence for any significant C-N π bonding is not strong. Hirshfeld atom refinement (HAR) gives H atom positional and anisotropic displacement parameters that agree very well with the neutron parameters. X-ray and neutron diffraction data on the dihydrate of carbemazepine strongly indicate a disordered orthorhombic crystal structure in the space group Cmca, rather than a monoclinic crystal structure in space group P2(1)/c. This disorder in the dihydrate structure has implications for both experimental and theoretical studies of polymorphism.

Ab initio and DFT studies on the structures, binding energies and nature of bonds in X2Y3 metal clusters (X+ = Li+, Na+ and K+; ${\rm Y}_{3}^{2-} = {\rm Zn}_{3}^{2-}, {\rm Cd}_{3}^{2-}$ and ${\rm Hg}_{3}^{2-}$)

A comprehensive theoretical study on the structure and stability of linear and triangular isomers of anionic clusters of zinc, cadmium and mercury [Formula: see text] and their binding with one and two alkali metal cations ( X + = Li +, Na +, K +) has been investigated at the density functional (BP86 and B3LYP) and ab initio (MP2, MP4 and coupled cluster single and double (CCSD)) methods. The results showed that in all cases, the triangular geometry with D3h symmetry is more stable than linear one. The calculated values of interaction energies (IE) between [Formula: see text] anions and two X+ cations, Wiberg bond indices (WBI) and the electron densities at bond critical points (BCP), ρ(BCP), for Y–X bonds show that among all complexes investigated here at all levels of theory Zn 3 Li 2 and Hg 3 K 2 have the largest and the smallest values of IE, WBI and ρ(BCP), respectively. The natural charges of the atoms and WBI involved in the bonding as well as the global value of the charge transfer ΔQ from [Formula: see text] to X+ cation in X2Y3 clusters, evaluated through natural population analysis, confirmed that covalent contribution in Y–X bond formation increases from K + to Li +. Also the energy decomposition analyses (EDA) were used to detect the nature of interaction in the complexes. The results confirmed that the contribution of electrostatic interactions in present complexes is almost more than 70%.

Unusual H-Bond Topology and Bifurcated H-bonds in the 2-Fluoroethanol Trimer

By using a combination of rotational spectroscopy and ab initio calculations, an unusual H-bond topology was revealed for the 2-fluoroethanol trimer. The trimer exhibits a strong heterochiral preference and adopts an open OH⋅⋅⋅OH H-bond topology while utilizing two types of bifurcated H-bonds involving organic fluorine. This is in stark contrast to the cyclic OH⋅⋅⋅OH H-bond topology adopted by trimers of water and other simple alcohols. The strengths of different H-bonds in the trimer were analyzed by using the quantum theory of atoms in molecules. The study showcases a remarkable example of a chirality-induced switch in H-bond topology in a simple transient chiral fluoroalcohol. It provides important insight into the H-bond topologies of small fluoroalcohol aggregates, which are proposed to play a key role in protein folding and in enantioselective reactions and separations where fluoroalcohols serve as a (co)solvent.