Gas-phase reactivity of protonated 2-oxazoline derivatives: mass spectrometry and computational studies

RATIONALE: Oxazolines have attracted the attention of researchers worldwide due to their versatility as carboxylic acid protecting groups, chiral auxiliaries, and ligands for asymmetric catalysis. Electrospray ionization tandem mass spectrometric (ESI-MS/MS) analysis of five 2-oxazoline derivatives has been conducted, in order to understand the influence of the side chain on the gas-phase dissociation of these protonated compounds under collision-induced dissociation (CID) conditions. METHODS: Mass spectrometric analyses were conducted in a quadrupole time-of-flight (Q-TOF) spectrometer fitted with electrospray ionization source. Protonation sites have been proposed on the basis of the gas-phase basicity, proton affinity, atomic charges, and a molecular electrostatic potential map obtained on the basis of the quantum chemistry calculations at the B3LYP/6-31 + G(d, p) and G2(MP2) levels. RESULTS: Analysis of the atomic charges, gas-phase basicity and proton affinities values indicates that the nitrogen atom is a possible proton acceptor site. On the basis of these results, two main fragmentation processes have been suggested: one taking place via neutral elimination of the oxazoline moiety (99 u) and another occurring by sequential elimination of neutral fragments with 72 u and 27 u. These processes should lead to formation of R+. CONCLUSIONS: The ESI-MS/MS experiments have shown that the side chain could affect the dissociation mechanism of protonated 2-oxazoline derivatives. For the compound that exhibits a hydroxyl at the lateral chain, water loss has been suggested to happen through an E2-type elimination, in an exothermic step. Copyright (C) 2012 John Wiley & Sons, Ltd.

Effect of hydrogen bond formation on the elastic molecular scattering: a case study with methanol

Rayleigh optical activities of small hydrogen-bonded methanol clusters containing two to five molecules are reported. For the methanol trimer, tetramer, and pentamer both cyclic and linear structures are considered. After the geometry optimizations, the dipole moments and the dipole polarizabilities (mean, interaction, and anisotropic components) are calculated using HF, MP2 and DFT (B3LYP, B3P86 and BH&HLYP) with aug-cc-pVDZ extended basis set. The polarizabilities are used to analyse the depolarization ratios and the Rayleigh scattering activities. The variations in the activity and in the depolarization for Rayleigh scattered radiation with the increase in the cluster size for both cyclic and linear structures are analysed.

Host-guest interactions between xanthones and water: the role of O-H center dot center dot center dot O, C-H center dot center dot center dot O, and pi center dot center dot center dot pi contacts in the channel- and cage-type frameworks

In this article were studied two xanthone derivatives known as 1,5-dihydroxy-8-methoxyxanthone (I) and 1,3,7-trihydroxy-8-methoxyxanthone (II), which show one water molecule into their crystal structures. In xanthone I, there are water wires contributing to build up channel-like cavities along the c axis, whereas in xanthone II the water is surrounded by three xanthone molecules forming a cage-type structure. The geometries of I and II were optimized using the density functional theory method with B3LYP functional, and the results were compared with crystal structure. Both theoretical and experimental investigations reveal a concordance between structural parameters, with the xanthone core presenting an almost flat conformation and substituents adopting the more stable orientations. In the two compounds, the hydroxyl group linked at position 1 is involved in a resonance-assisted hydrogen bond with the carbonyl group. Besides, the supramolecular arrangement of the host/guest systems are stabilized mainly by classical intermolecular hydrogen bonds (O-H center dot center dot center dot O) involving xanthone-to-water and xanthone-to-xanthone. In addition, C-H center dot center dot center dot O weak hydrogen bonds, as well as pi-pi interactions play an important role to stabilize the crystal self-assembly of xanthones I and II. The results reported here underline the role of inclusion of water molecules and their different arrangement into the crystal structure of two xanthone host/guest systems.

Theoretical models for the antitrypanosomal activity of thiosemicarbazone derivatives

Thiosemicarbazones are cruzain inhibitors which have been identified as potential antitrypanosomal agents. In this work, several molecular properties were calculated at the density functional theory (DFT)/B3LYP/6-311G* level for a set of 44 thiosemicarbazones. Unsupervised and supervised pattern recognition techniques (hierarchical cluster analysis, principal component analysis, kth-nearest neighbors, and soft independent modeling by class analogy) were used to obtain structureactivity relationship models, which are able to classify unknown compounds according to their activities. The chemometric analyses performed here revealed that 12 descriptors can be considered responsible for the discrimination between high and low activity compounds. Classification models were validated with an external test set, showing that predictive classifications were achieved with the selected variable set. The results obtained here are in good agreement with previous findings from the literature, suggesting that our models can be useful on further investigations on the molecular determinants for the antichagasic activity. (C) 2012 Wiley Periodicals, Inc.

Metal binding selectivity of oxa-aza macrocyclic ligand: a DFT study of first- and second-row transition metal for four coordination systems

A detailed theoretical study of the 1,7,1l,17-tetraoxa-2,6,12,16-tetraaza-cycloeicosane ligand ([20]AneN(4)O(4)) coordinated to Fe2+, Co2+, Ni2+, Ru2+, Rh2+, and Pd2+ transition metal ions was carried out with the B3LYP method. Two different cases were performed: when nitrogen is the donor atom (1a (q) ) and also with the oxygen as the donor atom (1b (q) ). For all the cases performed in this study 1a (q) structures were always more stable than the 1b (q) ones. Considering each row is possible to see that the energy increases with the increase of the atomic number. The M2+ cation binding energies for the 1a (q) complexes increase with the following order: Fe2+ < Ru2+ < Co2+ < Ni2+ < Rh2+ < Pd2+.

Electrostatic potentials and polarization effects in proton-molecule interactions by means of multipoles from the quantum theory of atoms in molecules

Some atomic multipoles (charges, dipoles and quadrupoles) from the Quantum Theory of Atoms in Molecules (QTAIM) and CHELPG charges are used to investigate interactions between a proton and a molecule (F2, Cl2, BF, AlF, BeO, MgO, LiH, H2CO, NH3, PH3, BF3, and CO2). Calculations were done at the B3LYP/6-311G(3d,3p) level. The main aspect of this work is the investigation of polarization effects over electrostatic potentials and atomic multipoles along a medium to long range of interaction distances. Large electronic charge fluxes and polarization changes are induced by a proton mainly when this positive particle approaches the least electronegative atom of diatomic heteronuclear molecules. The search for simple equations to describe polarization on electrostatic potentials from QTAIM quantities resulted in linear relations with r-4 (r is the interaction distance) for many cases. Moreover, the contribution from atomic dipoles to these potentials is usually the most affected contribution by polarization what reinforces the need for these dipoles to a minimal description of purely electrostatic interactions. Finally, CHELPG charges provide a description of polarization effects on electrostatic potentials that is in disagreement with physical arguments for certain of these molecules. (c) 2012 Wiley Periodicals, Inc.

Surface-enhanced Raman spectroscopy studies of organophosphorous model molecules and pesticides

This work reports the analytical application of surface-enhanced Raman spectroscopy (SERS) in the trace analysis of organophosphorous pesticides (trichlorfon and glyphosate) and model organophosphorous compounds (dimethyl methylphosphonate and o-ethyl methylphosphonothioate) bearing different functional groups. SERS measurements were carried out using Ag nanocubes with an edge square dimension of ca. 100 nm as substrates. Density functional theory (DFT) with the B3LYP functional was used for the optimization of ground state geometries and simulation of Raman spectra of the organophosphorous compounds and their silver complexes. Adsorption geometries and marker bands were identified for each of the investigated compound. Results indicate the usefulness of SERS methodology for the sensitive analyses of organophosphorous compounds through the use of vibrational spectroscopy.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

The basic antioxidant structure for flavonoid derivatives

An antioxidant structure-activity study is carried out in this work with ten flavonoid compounds using quantum chemistry calculations with the functional of density theory method. According to the geometry obtained by using the B3LYP/6-31G(d) method, the HOMO, ionization potential, stabilization energies, and spin density distribution showed that the flavonol is the more antioxidant nucleus. The spin density contribution is determinant for the stability of the free radical. The number of resonance structures is related to the pi-type electron system. 3-hydroxyflavone is the basic antioxidant structure for the simplified flavonoids studied here. The electron abstraction is more favored in the molecules where ether group and 3-hydroxyl are present, nonetheless 2,3-double bond and carbonyl moiety are facultative.

Corrole isomers: intrinsic gas-phase shapes via traveling wave ion mobility mass spectrometry and dissociation chemistries via tandem mass spectrometry

Corrole and four of its isomers with subtle structural changes promoted by exchange of nitrogen and carbon atoms in the corrole ring have been studied by traveling wave ion mobility mass spectrometry and collision induced dissociation experiments. Significant differences in shapes and charge distributions for their protonated molecules were found to lead to contrasting gas phase mobilities, most particularly for corrorin, the most "confused" isomer. Accordingly, corrorin was predicted by B3LYP/6-31g(d,p) and collisional cross section calculations to display the most compact tri-dimensional structure, whereas NCC4 and corrole were found to be the most planar isomers. Better resolution between the corrole isomers was achieved using the more polarizable and massive CO2 as the drift gas. Sequential losses of HF molecules were found to dominate the dissociation chemistry of the protonated molecules of these corrole isomers, but their unique structures caused contrasting labilities towards CID, whereas NCC4 showed a peculiar and structurally diagnostic loss of NH3, allowing its prompt differentiation from the other isomers.

First- and second-row transition metal oxa-aza macrocyclic complexes: a DFT study of an octahedral conformation

A theoretical study of structures of the 1,7,1 l,17-tetraoxa-2,6,12,16-tetraaza-cycloeicosane ligand ([20]AneN(4)O(4)) coordinated to Fe2+, Co2+, Ni2+, Ru2+, Rh2+, and Pd2+ transition metals ions was carried out with the DFT/B3LYP method. Complexes were fully optimized in C-s symmetry with the metal ions coordinated either to nitrogen (1a) or oxygen atoms (1b). For all the cases performed in this work, 1a was always more stable than 1b. Considering each row it is possible to see that the binding energy increases with the atomic number. The M2+ cation binding energies increase in the following order: Fe2+ < Ru2+ < Co2+ < Ni2+ < Rh2+ < Pd2+. In addition, it was observed the preference of Pd2+ and Rh2+ complexes for a tetrahedral arrangement, while Fe2+, Ru2+, Co2+, Ni2+ complexes had a preference for the octahedral arrangement. From the orbital representation results, it was seen that 1b unsymmetrical orbitals may influence the susceptibility over metal ions orientation toward heteroatoms orbitals.

Theoretical study of the XP3 (X = Al, B, Ga) clusters

The lowest singlet and triplet states of AlP3, GaP3 and BP3 molecules with C-s, C-2v and C-3v symmetries were characterized using the B3LYP functional and the aug-cc-pVTZ and aug-cc-pVQZ correlated consistent basis sets. Geometrical parameters and vibrational frequencies were calculated and compared to existent experimental and theoretical data. Relative energies were obtained with single point CCSD(T) calculations using the aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets, and then extrapolating to the complete basis set (CBS) limit. (C) 2011 Elsevier B.V. All rights reserved.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

A Combined Study Using Ligand-Based Design, Synthesis, and Pharmacological Evaluation of Analogues of the Acetaminophen Ortho-Regioisomer with Potent Analgesic Activity

A ligand-based drug design study was performed to acetaminophen regioisomers as analgesic candidates employing quantum chemical calculations at the DFT/B3LYP level of theory and the 6-31G* basis set. To do so, many molecular descriptors were used such as highest occupied molecular orbital, ionization potential, HO bond dissociation energies, and spin densities, which might be related to quench reactivity of the tyrosyl radical to give N-acetyl-p-benzosemiquinone-imine through an initial electron withdrawing or hydrogen atom abstraction. Based on this in silico work, the most promising molecule, orthobenzamol, was synthesized and tested. The results expected from the theoretical prediction were confirmed in vivo using mouse models of nociception such as writhing, paw licking, and hot plate tests. All biological results suggested an antinociceptive activity mediated by opioid receptors. Furthermore, at 90 and 120 min, this new compound had an effect that was comparable to morphine, the standard drug for this test. Finally, the pharmacophore model is discussed according to the electronic properties derived from quantum chemistry calculations.