A Selective and Functional Group-Tolerant Ruthenium-Catalyzed Olefin Metathesis/Transfer Hydrogenation Tandem Sequence Using Formic Acid as Hydrogen Source.

A ruthenium-catalyzed transfer hydrogenation of olefins utilizing formic acid as a hydrogen donor is described. The application of commercially available alkylidene ruthenium complexes opens access to attractive C(sp3)-C(sp3) bond formation in an olefin metathesis/transfer hydrogenation sequence under tandem catalysis conditions. High chemoselectivity of the developed methodology provides a remarkable synthetic tool for the reduction of various functionalized alkenes under mild reaction conditions. The developed methodology is applied for the formal synthesis of the drugs pentoxyverine and bencyclane.

Electrophilicity of oxalic acid monomer is enhanced in the dimer by intermolecular proton transfer.

We have analyzed the effect of excess electron attachment on the network of hydrogen bonds in the oxalic acid dimer (OA)2. The most stable anionic structures may be viewed as complexes of a neutral hydrogenated moiety HOA˙ coordinated to an anionic deprotonated moiety (OA-H)-. HOA˙ acts as a double proton donor and (OA-H)- as a double proton acceptor. Thus the excess electron attachment drives intermolecular proton transfer. We have identified several cyclic hydrogen bonded structures of (OA)2-. Their stability has been analyzed in terms of the stability of the involved conformers, the energetic penalty for deformation of these conformers to the geometry of the dimer, and the two-body interaction energy between the deformed HOA˙ and (OA-H)-. There are at least seven isomers of (OA)2- with stabilization energies in the range of 1.26-1.39 eV. These energies are dominated by attractive two-body interaction energies. The anions are vertically bound electronically by 3.0-3.4 eV and adiabatically bound by at least 1.6 eV. The computational predictions are consistent with the anion photoelectron spectrum of (OA)2-. The spectrum consists of a broad feature, with an onset of 2.5 eV and spanning to 4.3 eV. The electron vertical detachment energy (VDE) is assigned to be 3.3 eV.

Different Conformations of 2'-Deoxycytidine in the Gas and Solid Phases: Competition between Intra- and Intermolecular Hydrogen Bonds.

Computational results have been reported for 2'-deoxycytidine (dC), its gas phase isomers, tautomers, and their conformers, as well as for the crystalline phase. In addition to the neutral gas phase molecules, we have also considered associated radical anions and cations. The structural calculations were performed at the density functional and MP2 levels of theory. Vertical electron ionization energies and excess electron binding energies were determined using electron propagator theory. The α-anomer proved to be more stable by a fraction of kcal/mol than the biologically relevant canonical ß-anomer. The conformational space of canonical dC has been systematically probed. dC in the crystalline phase or DNA structures favors canonical anti conformations. These structures were used in past computational studies to model gas phase characteristics of dC. Our findings indicate, however, that the gas phase dC favors syn conformations. It has repercussions for earlier interpretations of gas phase experimental results based on these computational results. The thermodynamic dominance of syn conformations results from the formation of an intramolecular O5'-H13···O2 hydrogen bond. The IR spectra of the most stable syn and anti canonical conformers differ markedly in the region of frequencies corresponding to NH/OH stretching modes. The MP2 value of deprotonation enthalpy of dC of 1411.7 kJ/mol is in very good agreement with the experimental value of 1409 ± 2.5 kJ/mol. The most stable valence anions are characterized by electron vertical detachment energies (VDE) in the 0.8-1.0 eV range, in good agreement with the experimental VDE of 0.87 eV. The barrier for the glycosidic bond cleavage is significant in the neutral canonical dC, 40.0 kcal/mol, and it is reduced to 22 and 16 kcal/mol for the anionic and cationic radicals of dC, respectively. The cleavage reaction is exothermic by 4 kcal/mol for dC- and endothermic by 7 and 9 kcal/mol for dC+ and dC, respectively. We decomposed the crystal cohesive energy into repulsive one-body terms associated with the syn-anti conformational changes, and the attractive intermolecular interaction term. We exposed that the syn-anti conformational changes are very favorable for intermolecular interactions; in particular they make the imino-amino side of the cytosine residue accessible to intermolecular interactions.

Importance of Time Scale and Local Environment in Electron-Driven Proton Transfer. The Anion of Acetoacetic Acid.

Anion photoelectron spectroscopy (PES) and electron energy-loss spectroscopy (EELS) probe different regions of the anionic potential energy surface. These complementary techniques provided information about anionic states of acetoacetic acid (AA). Electronic structure calculations facilitated the identification of the most stable tautomers and conformers for both neutral and anionic AA and determined their relative stabilities and excess electron binding energies. The most stable conformers of the neutral keto and enol tautomers differ by less than 1 kcal/mol in terms of electronic energies corrected for zero-point vibrations. Thermal effects favor these conformers of the keto tautomer, which do not support an intramolecular hydrogen bond between the keto and the carboxylic groups. The valence anion displays a distinct minimum which results from proton transfer from the carboxylic to the keto group; thus, we name it an ol structure. The minimum is characterized by a short intramolecular hydrogen bond, a significant electron vertical detachment energy of 2.38 eV, but a modest adiabatic electron affinity of 0.33 eV. The valence anion was identified in the anion PES experiments, and the measured electron vertical detachment energy of 2.30 eV is in good agreement with our computational prediction. We conclude that binding an excess electron in a π* valence orbital changes the localization of a proton in the fully relaxed structure of the AA(-) anion. The results of EELS experiments do not provide evidence for an ultrarapid proton transfer in the lowest π* resonance of AA(-), which would be capable of competing with electron autodetachment. This observation is consistent with our computational results, indicating that major gas-phase conformers and tautomers of neutral AA do not support the intramolecular hydrogen bond that would facilitate ultrarapid proton transfer and formation of the ol valence anion. This is confirmed by our vibrational EELS spectrum. Anions formed by vertical electron attachment to dominant neutrals undergo electron autodetachment with or without vibrational excitations but are unable to relax to the ol structure on a time scale fast enough to compete with autodetachment.

Intermolecular interactions between molecules in various conformational states: the dimer of oxalic acid.

We considered stability of the dimer of oxalic acid. The global minimum energy structure identified by us is stabilized by two inter- and four intramolecular hydrogen bonds, whereas the most stable structure identified in previous studies is supported by two inter- and three intramolecular hydrogen bonds. The latter structure proves to be less stable by 25 meV than the former. The global minimum stability results from a balancing act between a moderately attractive two-body interaction energy and small repulsive one-body terms. We have analyzed zero-point vibrational corrections to the stability of various conformers of oxalic acid and their dimers. We have found that minimum energy structures with the most stabilizing sets of hydrogen bonds have the largest zero-point vibrational energy, contrary to a naive anticipation based on red shifts of OH stretching modes involved in hydrogen bonds.

Communication: Remarkable electrophilicity of the oxalic acid monomer: an anion photoelectron spectroscopy and theoretical study.

Our experimental and computational results demonstrate an unusual electrophilicity of oxalic acid, the simplest dicarboxylic acid. The monomer is characterized by an adiabatic electron affinity and electron vertical detachment energy of 0.72 and 1.08 eV (±0.05 eV), respectively. The electrophilicity results primarily from the bonding carbon-carbon interaction in the singly occupied molecular orbital of the anion, but it is further enhanced by intramolecular hydrogen bonds. The well-resolved structure in the photoelectron spectrum is reproduced theoretically, based on Franck-Condon factors for the vibronic anion → neutral transitions.

SSC: a tool for constructing libraries for systematic screening of conformers.

Even a relatively small molecule with 10-20 atoms might have a few local minima, which correspond to different conformers. The number of local minima quickly increases with molecular size and the most common algorithms, driven by calculated forces, frequently identify a minimum, which is closest to the initial structure, rather than the most stable conformer. Here we discuss how to perform a systematic search of the conformational space for a chain-like molecule. Our approach is fully automated and a user has control which chemical bonds will be probed and with which increments. Moreover, whole fragments of the molecule, which are adjacent to each selected rotational bond, are rotated in a properly selected cylindrical coordinate system and unchemical hybridizations and some "clashes" between neighboring groups, which are common when standard Z-matrices are used, are avoided. A library of potentially relevant conformers is created with a tool, which we call SSC, denoting Systematic Screening of Conformers. Each member of the library is prescreened at a predefined level of theory and the most promising conformers are identified. Finally, they are further evaluated at a higher level of theory to identify the most stable structures and their physicochemical properties. As an example, we demonstrate the results of this approach for 2'-deoxycytidine.

Is electronegativity a useful descriptor for the pseudo-alkali metal NH4?

Molecular ions in the form of "pseudo-atoms" are common structural motifs in chemistry, with properties that are transferrable between different compounds. We have determined one such property--the electronegativity--for the "pseudo-alkali metal" ammonium (NH(4)), and evaluated its reliability as a descriptor versus the electronegativities of the alkali metals. The computed properties of ammonium's binary complexes with astatine and of selected borohydrides confirm the similarity of NH(4) to the alkali metal atoms, although the electronegativity of NH(4) is relatively large in comparison to its cationic radius. We have paid particular attention to the molecular properties of ammonium (angular anisotropy, geometric relaxation and reactivity), which can cause deviations from the behaviour expected of a conceptual "true alkali metal" with this electronegativity. These deviations allow for the discrimination of effects associated with the molecular nature of NH(4).

Theoretical investigations on the formation and dehydrogenation reaction pathways of H(NH2BH2)(n)H (n = 1-4) oligomers: importance of dihydrogen interactions.

The H(NH(2)BH(2))(n)H oligomers are possible products from dehydrogenation of ammonia borane (NH(3)BH(3)) and ammonium borohydride (NH(4)BH(4)), which belong to a class of boron-nitrogen-hydrogen (BNH(x)) compounds that are promising materials for chemical hydrogen storage. Understanding the kinetics and reaction pathways of formation of these oligomers and their further dehydrogenation is essential for developing BNH(x)-based hydrogen storage materials. We have performed computational modeling using density functional theory (DFT), ab initio wave function theory, and Car-Parrinello molecular dynamics (CPMD) simulations on the energetics and formation pathways for the H(NH(2)BH(2))(n)H (n = 1-4) oligomers, polyaminoborane (PAB), from NH(3)BH(3) monomers and the subsequent dehydrogenation steps to form polyiminoborane (PIB). Through computational transition state searches and evaluation of the intrinsic reaction coordinates, we have investigated the B-N bond cleavage, the reactions of NH(3)BH(3) molecule with intermediates, dihydrogen release through intra- and intermolecular hydrogen transfer, dehydrocoupling/cyclization of the oligomers, and the dimerization of NH(3)BH(3) molecules. We find that the formation of H(NH(2)BH(2))(n+1)H oligomers occurs first through reactions of the H(NH(2)BH(2))(n)H oligomers with BH(3) followed by reactions with NH(3) and the release of H(2), where the BH(3) and NH(3) intermediates are formed through dissociation of NH(3)BH(3). We also find that the dimerization of the NH(3)BH(3) molecules to form cyclic c-(NH(2)BH(2))(2) is slightly exothermic, with an unexpected transition state that leads to the simultaneous release of two H(2) molecules. The dehydrogenations of the oligomers are also exothermic, typically by less than 10 kcal/(mol of H(2)), with the largest exothermicity for n = 3. The transition state search shows that the one-step direct dehydrocoupling cyclization of the oligomers is not a favored pathway because of high activation barriers. The dihydrogen bonding, in which protic (H(N)) hydrogens interact with hydridic (H(B)) hydrogens, plays a vital role in stabilizing different structures of the reactants, transition states, and products. The dihydrogen interaction (DHI) within the R-BH(2)(eta(2)-H(2)) moiety accounts for both the formation mechanisms of the oligomers and for the dehydrogenation of ammonia borane.

Combinatorial-computational-chemoinformatics (C3) approach to finding and analyzing low-energy tautomers.

Finding the most stable tautomer or a set of low-energy tautomers of molecules is critical in many aspects of molecular modelling or virtual screening experiments. Enumeration of low-energy tautomers of neutral molecules in the gas-phase or typical solvents can be performed by applying available organic chemistry knowledge. This kind of enumeration is implemented in a number of software packages and it is relatively reliable. However, in esoteric cases such as charged molecules in uncommon, non-aqueous solvents there is simply not enough available knowledge to make reliable predictions of low energy tautomers. Over the last few years we have been developing an approach to address the latter problem and we successfully applied it to discover the most stable anionic tautomers of nucleic acid bases that might be involved in the process of DNA damage by low-energy electrons and in charge transfer through DNA. The approach involves three steps: (1) combinatorial generation of a library of tautomers, (2) energy-based screening of the library using electronic structure methods, and (3) analysis of the information generated in step (2). In steps 1-3 we employ combinatorial, computational and chemoinformatics techniques, respectively. Therefore, this hybrid approach is named "Combinatorial*Computational*Chemoinformatics", or just abbreviated as C(3) (or C-cube) approach. This article summarizes our developments and most interesting methodological aspects of the C(3) approach. It can serve as an example how to identify the most stable tautomers of molecular systems for which common chemical knowledge had not been sufficient to make definite predictions.

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids.

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.

Structure and singly occupied molecular orbital analysis of anionic tautomers of guanine.

Recently, we reported the discovery of adiabatically bound anions of guanine that might be involved in the processes of DNA damage by low-energy electrons and in charge transfer through DNA. These anions correspond to some tautomers that have been ignored thus far. They were identified using a hybrid quantum mechanical-combinatorial approach in which an energy-based screening was performed on the library of 499 tautomers with their relative energies calculated with quantum chemistry methods. In the current study, we analyze the adiabatically bound anions of guanine in two aspects: (1) the geometries and excess electron distributions are analyzed and compared with anions of the most stable neutrals to identify the sources of stability; (2) the chemical space of guanine tautomers is explored to verify if these new tautomers are contained in a particular subspace of the tautomeric space. The first task involves the development of novel approaches-the quantum chemical data like electron density, orbital, and information on its bonding/antibonding character are coded into holograms and analyzed using chemoinformatics techniques. The second task is completed using substructure analysis and clustering techniques performed on molecules represented by 2D fingerprints. The major conclusion is that the high stability of adiabatically bound anions originates from the bonding character of the pi orbital occupied by the excess electron. This compensates for the antibonding character that usually causes significant buckling of the ring. Also, the excess electron is more homogenously distributed over both rings than in the case of anions of the most stable neutral species. In terms of 2D substructure, the most stable anionic tautomers generally have additional hydrogen atoms at C8 and/or C2 and they do not have hydrogen atoms attached to C4, C5, and C6. They also form an "island of stability" in the tautomeric space of guanine.

Cylindrical projection of electrostatic potential and image analysis tools for damaged DNA: the substitution of thymine with thymine glycol.

Changes of electrostatic potential around the DNA molecule resulting from chemical modifications of nucleotides may play a role in enzymatic recognition of damaged sites. Effects of chemical modifications of nucleotides on the structure of DNA have been characterized through electronic structure computations. Quantum mechanical structural optimizations of fragments of five pairs of nucleotides with thymine or thymine glycol were performed at the density functional level of theory with a B3LYP exchange-correlation functional and 6-31G(d,p) basis sets. The electrostatic potential (EP) around DNA fragments was projected on a cylindrical surface around the double helix. The 2D maps of EP of intact and damaged DNA fragments were compared using image analysis methods to identify and measure modifications of the EP that result from the occurrence of thymine glycol. It was found that distortions of phosphate groups and displacements of the accompanying countercations by up to approximately 0.5 angstroms along the axis of DNA are clearly reflected in the EP maps. Modifications of the EP in the major groove of DNA near the damaged site are also reported.

Photoelectron spectrum of valence anions of uracil and first-principles calculations of excess electron binding energies.

The photoelectron spectrum (PES) of the uracil anion is reported and discussed from the perspective of quantum chemical calculations of the vertical detachment energies (VDEs) of the anions of various tautomers of uracil. The PES peak maximum is found at an electron binding energy of 2.4 eV, and the width of the main feature suggests that the parent anions are in a valence rather than a dipole-bound state. The canonical tautomer as well as four tautomers that result from proton transfer from an NH group to a C atom were investigated computationally. At the Hartree-Fock and second-order Moller-Plesset perturbation theory levels, the adiabatic electron affinity (AEA) and the VDE have been converged to the limit of a complete basis set to within +/-1 meV. Post-MP2 electron-correlation effects have been determined at the coupled-cluster level of theory including single, double, and noniterative triple excitations. The quantum chemical calculations suggest that the most stable valence anion of uracil is the anion of a tautomer that results from a proton transfer from N1H to C5. It is characterized by an AEA of 135 meV and a VDE of 1.38 eV. The peak maximum is as much as 1 eV larger, however, and the photoelectron intensity is only very weak at 1.38 eV. The PES does not lend support either to the valence anion of the canonical tautomer, which is the second most stable anion, and whose VDE is computed at about 0.60 eV. Agreement between the peak maximum and the computed VDE is only found for the third most stable tautomer, which shows an AEA of approximately -0.1 eV and a VDE of 2.58 eV. This tautomer results from a proton transfer from N3H to C5. The results illustrate that the characteristics of biomolecular anions are highly dependent on their tautomeric form. If indeed the third most stable anion is observed in the experiment, then it remains an open question why and how this species is formed under the given conditions.

Adiabatically bound valence anions of guanine.

It is believed that guanine, a basic component of DNA and RNA, has the smallest affinity to an excess electron among all nucleic acid bases. Our experimental and computational findings indicate, however, that many so far neglected tautomers of guanine support adiabatically bound anionic states in the gas phase. The computed values of electron vertical detachment energy for the most stable anionic tautomers are within a broad range of the dominant feature of the photoelectron spectrum. We suggest that guanine might be the strongest excess electron acceptor among nucleic acid bases. Thus it might be critical to radiobiological damage of DNA and it might contribute to those chemical transformations of DNA that proceed through bound anionic states.

Differences in electrostatic potential around DNA fragments containing adenine and 8-oxo-adenine. An analysis based on regular cylindrical projection.

Changes of electrostatic potential (EP) around the DNA molecule resulting from chemical modifications of nucleotides may play a role in enzymatic recognition of damaged sites. Effects of chemical modifications of nucleotides on the structure of DNA have been characterized through large scale density functional theory computations. Quantum mechanical structural optimizations of DNA fragments with three pairs of nucleotides and accompanying counteractions were performed with a B3LYP exchange-correlation functional and 6-31G(d,p) basis sets. The "intact" DNA fragment contained adenine in the middle layer, while the "damaged" fragment had the adenine replaced with 8-oxo-adenine. The electrostatic potential around these DNA fragments was projected on a cylindrical surface around the double helix. The two-dimensional maps of EP of the intact and damaged DNA fragments were analyzed to identify these modifications of EP that result from the occurrence of 8-oxo-adenine (8oA). It was found that distortions of a phosphate group neighboring 8oA and displacements of the accompanying countercation are clearly reflected in the EP maps.

On the unusual stability of valence anions of thymine based on very rare tautomers: A computational study.

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).

Intermolecular proton transfer in anionic complexes of uracil with alcohols.

A series of 18 alcohols (ROH) has been designed with an enthalpy of deprotonation in the gas phase (H(DP)) in the range 13.8-16.3 eV. The effects of excess electron attachment to the binary alcohol-uracil (ROH...U) complexes have been studied at the density functional level with a B3LYP exchange-correlation functional and at the second-order Møller-Plesset perturbation theory level. The photoelectron spectra of anionic complexes of uracil with 3 alcohols (ethanol, 2,2,3,3,3-pentafluoropropanol, and 1,1,1,3,3,3-hexafluoro-2-propanol) have been measured with 2.54 eV photons. For ROHs with deprotonation enthalpies larger than 14.8 eV, only the ROH...U- minimum exists on the potential energy surface of the anionic complex. For alcohols with deprotonation enthalpies in the range 14.3-14.8 eV, two minima might exist on the anionic potential energy surface, which correspond to the RO-...HU* and ROH...U- structures. For ROHs with deprotonation enthalpies smaller than 14.3 eV, the excess electron attachment to the ROH...U complex always induces a barrier-free proton transfer from the hydroxyl group of ROH to the O8 atom of U, with the product being RO-...HU*.

Discovery of Most Stable Structures of Neutral and Anionic Phenylalanine through Automated Scanning of Tautomeric and Conformational Spaces.

We have developed a software tool for combinatorial generation of tautomers and conformers of small molecules. We have demonstrated it by performing a systematic search for the most stable structures of neutral and anionic phenylalanine (Phe) using electronic structure methods. For the neutral canonical tautomer we found out that the conformers with and without the intramolecular (O)H···NH2 hydrogen bond are similarly stable, within the error bars of our method. A unique IR signature of the conformer without the hydrogen bond has been identified. We also considered anions of Phe, both valence type and dipole-bound. We have found out that tautomers resulting from proton transfer from the carboxylic OH to the phenyl ring do support valence anions that are vertically strongly bound, with electron vertical detachment energies (VDE) in a range of 3.2-3.5 eV. The most stable conformer of these valence anions remains adiabatically unbound with respect to the canonical neutral by only 2.17 kcal/mol at the CCSD(T)/aug-cc-pVDZ level. On the basis of our past experience with valence anions of nucleic acid bases, we suggest that the valence anions of Phe identified in this report can be observed experimentally. The most stable conformer of canonical Phe is characterized by an adiabatic electron affinity of 53 meV (a dipole-bound state).