Modelos contínuos do solvente: fundamentos

Continuum solvation models are nowadays widely used in the modeling of solvent effects and the range of applications goes from the calculation of partition coefficients to chemical reactions in solution. The present work presents a detailed explanation of the physical foundations of continuum models. We discuss the polarization of a dielectric and its representation through the volume and surface polarization charges. The Poisson equation for a dielectric was obtained and we have also derived and discuss the apparent surface charge method and its application for free energy of solvation calculations.

Estrutura e propriedades da difenil-4-amina sulfonato de sódio: um dopante de filmes condutores

This work contains the theoretical simulation of the conformation of diphenyl-4-amine sodium sulphonate (DASNa) and correlates its geometry with conductivity, showing that the conductivity increases as the molecule becomes more planar. The solvent effect was also evaluated, using water and dimethylsulfoxide. Some properties, such as bond distance, vibration al frequency and effective charge were calculated. The utilization of diphenyl-4-amine sodium sulphonate (DASNa) as dopant of aniline was investigated in view of the HOMO-LUMO energy gap.

Metodologia AGOA: a modelagem de clusters de hidratação no complexo aziridina···ácido fluorídrico

We present a theoretical study of solvent effect on C2H5N···HF hydrogen-bonded complex through the application of the AGOA methodology. By using the TIP4P model to orientate the configuration of water molecules, the hydration clusters generated by AGOA were obtained through the analysis of the molecular electrostatic potential (MEP) of solute (C2H5N···HF). Thereby, it was calculated the hydration energies on positive and negative MEP fields, which are maxima (PEMmax) and minima (PEMmin) when represent the -CH2- methylene groups and hydrofluoric acid, respectively. By taking into account the higher and lower hydration energy values of -370.6 kJ mol-1 and -74.3 kJ mol-1 for PEMmax and PEMmin of the C2H5N···HF, our analysis shows that these results corroborate the open ring reaction of aziridine, in which the preferential attack of water molecules occurs at the methylene groups of this heterocyclic.

Terpenoid transformations over gold catalysts

Nowadays biomass transformation has a great potential for the synthesis of value-added compounds with a wide range of applications. Terpenoids, extracted from biomass, are inexpensive and renewable raw materials which often have a biological activity and are widely used as important organic platform molecules in the development of new medicines as well as in the synthesis of fine chemicals and intermediates. At the same time, special attention is devoted to the application of gold catalysts to fine chemical synthesis due to their outstanding activity and/or selectivity for transformations of complex organic compounds. Conversion of renewable terpenoids in the presence of gold nanoparticles is one of the new and promising directions in the transformation of biomass to valuable chemicals. In the doctoral thesis, different kinds of natural terpenoids, such as α-pinene, myrtenol and carvone were selected as starting materials. Gold catalysts were utilized for the promising routes of these compounds transformation. Investigation of selective α-pinene isomerization to camphene, which is an important step in an industrial process towards the synthesis of camphor as well as other valuable substrates for the pharmaceutical industry, was performed. A high activity of heterogeneous gold catalysts in the Wagner-Meerwein rearrangement was demonstrated for the first time. Gold on alumina carrier was found to reach the α-pinene isomerization conversion up to 99.9% and the selectivity of 60-80%, thus making this catalyst very promising from an industrial viewpoint. A detailed investigation of kinetic regularities including catalyst deactivation during the reaction was performed. The one-pot terpene alcohol amination, which is a promising approach to the synthesis of valuable complex amines having specific physiological properties, was investigated. The general regularities of the one-pot natural myrtenol amination in the presence of gold catalysts as well as a correlation between catalytic activity, catalyst redox treatment and the support nature were obtained. Catalytic activity and product distribution were shown to be strongly dependent on the support properties, namely acidity and basicity. The gold-zirconia (Au/ZrO2) catalyst pretreated under oxidizing atmosphere was observed to be rather active, resulting in the total conversion of myrtenol and the selectivity to the corresponding amine of about 53%. The reaction kinetics was modelled based on the mechanistic considerations with the catalyst deactivation step incorporated in the mechanism. Carvone hydrogenation over a gold catalyst was studied with the general idea of investigating both the activity of gold catalysts in competitive hydrogenation of different functional groups and developing an approach to the synthesis of valuable carvone derivatives. Gold was found to promote stereo- and chemoselective carvone hydrogenation to dihydrocarvone with a predominant formation of the trans-isomer, which generally is a novel synthetic method for an industrially valuable dihydrocarvone. The solvent effect on the catalytic activity as well as on the ratio between trans- and cis-dihydrocarvone was evaluated.

Novel technology for preparation of optically active chemicals

Asymmetric synthesis using modified heterogeneous catalysts has gained lots of interest in the production of optically pure chemicals, such as pharmaceuticals, nutraceuticals, fragrances and agrochemicals. Heterogeneous modified catalysts capable of inducing high enantioselectivities are preferred in industrial scale due to their superior separation and handling properties. The topic has been intensively investigated both in industry and academia. The enantioselective hydrogenation of ethyl benzoylformate (EBF) to (R)-ethyl mandelate over (-)-cinchonidine (CD)-modified Pt/Al2O3 catalyst in a laboratory-scale semi-batch reactor was studied as a function of modifier concentration, reaction temperature, stirring rate and catalyst particle size. The main product was always (R)-ethyl mandelate while small amounts of (S)-ethyl mandelate were obtained as by product. The kinetic results showed higher enantioselectivity and lower initial rates approaching asymptotically to a constant value as the amount of modifier was increased. Additionally, catalyst deactivation due to presence of impurities in the feed was prominent in some cases; therefore activated carbon was used as a cleaning agent of the raw material to remove impurities prior to catalyst addition. Detailed characterizations methods (SEM, EDX, TPR, BET, chemisorption, particle size distribution) of the catalysts were carried out. Solvent effects were also studied in the semi-batch reactor. Solvents with dielectric constant (e) between 2 and 25 were applied. The enantiomeric excess (ee) increased with an increase of the dielectric coefficient up to a maximum followed by a nonlinear decrease. A kinetic model was proposed for the enantioselectivity dependence on the dielectric constant based on the Kirkwood treatment. The non-linear dependence of ee on (e) successfully described the variation of ee in different solvents. Systematic kinetic experiments were carried out in the semi-batch reactor. Toluene was used as a solvent. Based on these results, a kinetic model based on the assumption of different number of sites was developed. Density functional theory calculations were applied to study the energetics of the EBF adsorption on pure Pt(1 1 1). The hydrogenation rate constants were determined along with the adsorption parameters by non-linear regression analysis. A comparison between the model and the experimental data revealed a very good correspondence. Transient experiments in a fixed-bed reactor were also carried out in this work. The results demonstrated that continuous enantioselective hydrogenation of EBF in hexane/2-propanol 90/10 (v/v) is possible and that continuous feeding of (-)-cinchonidine is needed to maintain a high steady-state enantioselectivity. The catalyst showed a good stability and high enantioselectivity was achieved in the fixed-bed reactor. Chromatographic separation of (R)- and (S)-ethyl mandelate originating from the continuous reactor was investigated. A commercial column filled with a chiral resin was chosen as a perspective preparative-scale adsorbent. Since the adsorption equilibrium isotherms were linear within the entire investigated range of concentrations, they were determined by pulse experiments for the isomers present in a post-reaction mixture. Breakthrough curves were measured and described successfully by the dispersive plug flow model with a linear driving force approximation. The focus of this research project was the development of a new integrated production concept of optically active chemicals by combining heterogeneous catalysis and chromatographic separation technology. The proposed work is fundamental research in advanced process technology aiming to improve efficiency and enable clean and environmentally benign production of enantiomeric pure chemicals.

Investigations into the determination of polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) by capillary gas chromatography

Factors involved in the determination of PAHs (16 priority PAHs as an example) and PCBs (10 PCB congeners, representing 10 isomeric groups) by capillary gas chromatography coupled with mass spectrometry (GC/MS, for PAHs) and electron capture detection (GC/ECD , for PCBs) were studied, with emphasis on the effect of solvent. Having various volatilities and different polarities, solvent studied included dichloromethane, acetonitrile, hexan e, cyclohexane, isooctane, octane, nonane, dodecane, benzene, toluene, p-xylene, o-xylene, and mesitylene. Temperatures of the capillary column, the injection port, the GC/MS interface, the flow rates of carrier gas and make-up gas, and the injection volume were optimized by one factor at a time method or simplex optimization method. Under the optimized conditions, both peak height and peak area of 16 PAHs, especially the late-eluting PAHs, were significantly enhanced (1 to 500 times) by using relatively higher boiling point solvents such as p-xylene and nonane, compared with commonly used solvents like benzene and isooctane. With the improved sensitivity, detection limits of between 4.4 pg for naphthalene and 30.8 pg for benzo[g,h,i]perylene were obtained when p-xylene was used as an injection solvent. Effect of solvent on peak shape and peak intensity were found to be greatly dependent on temperature parameters, especially the initial temperature of the capillary column. The relationship between initial temperature and shape of peaks from 16 PAHs and 10 PCBs were studied and compared when toluene, p-xylene, isooctane, and nonane were used as injection solvents. If a too low initial temperature was used, fronting or split of peaks was observed. On the other hand, peak tailing occurred at a too high initial column temperature. The optimum initial temperature, at which both peak fronting and tailing were avoided and symmetrical peaks were obtained, depended on both solvents and the stationary phase of the column used. On a methyl silicone column, the alkane solvents provided wider optimum ranges of initial temperature than aromatic solvents did, for achieving well-shaped symmetrical GC peaks. On a 5% diphenyl: 1% vinyl: 94% dimethyl polysiloxane column, when the aromatic solvents were used, the optimum initial temperature ranges for solutes to form symmetrical peaks were improved to a similar degree as those when the alkanes were used as injection solvents. A mechanism, based on the properties of and possible interactions among the analyte, the injection solvent, and the stationary phase of the capillary column, was proposed to explain these observations. The effect of initial temperature on peak height and peak area of the 16 PAHs and the 10 PCBs was also studied. The optimum initial temperature was found to be dependent on the physical properties of the solvent used and the amount of the solvent injected. Generally, from the boiling point of the solvent to 10 0C above its boiling point was an optimum range of initial temperature at which cthe highest peak height and peak area were obtained.

Infrared study of intercomponent interactions in a switchable hydrogen-bonded rotaxane

The macrocycle in rotaxane 1 is preferentially hydrogen bonded to the succinamide station in the neutral form, but can be moved to the naphthalimide station by one-electron reduction of the latter. The hydrogen bonding between the amide NH groups of the macrocycle and the C=O groups in the binding stations in the thread was studied with IR spectroscopy in different solvents in both states. In addition, the solvent effect on the vibrational frequencies was analyzed; a correlation with the solvent acceptor number (AN) was observed. The conformational switching upon reduction could be detected by monitoring the hydrogen-bond-induced shifts of the v(CO) frequencies of the C=O groups of the succinamide and the reduced naphthalimide stations. The macrocycle was found to shield the encapsulated station from the solvent: wavenumbers of v(CO) bands of the C=O groups residing inside the macrocycle cavity remain unaffected by the solvent polarity.

Excited state electronic polarization and reappraisal of the n <- pi* emission of acetone in water

Electronic polarization of the acetone molecule in the excited n -> pi* state is considered and its influence on the solvent shift in the emission spectrum is analyzed. Using an iterative procedure the electronic polarizations of both the ground and the excited states are included and compared with previous results obtained with Car-Parrinello dynamics. Analysis of the emission transition obtained using CIS(D)/aug-cc-pVDZ on statistically uncorrelated solute-solvent structures, composed of acetone and twelve explicit water molecules embedded in the electrostatic field of remaining 263 water molecules, corroborates that the solvent effect is mild, calculated here between 80 and 380 cm (1). (c) 2010 Published by Elsevier B.V.

Electronic spectroscopy of biomolecules in solution: fluorescein dianion in water

A combined and sequential Monte Carlo-quantum mechanics methodology is used to describe the electronic absorption spectrum of the fluorescein dianion in water. Different sets of 100 statistically relevant configurations composed of the solute and several solvent molecules are sampled from the Monte Carlo simulation for a posteriori quantum mechanical calculations of the spectra. In the largest case the configurations are composed of fluorescein and 90 explicit water molecules embedded in the electrostatic field of all remaining water molecules within a distance of 11.3 angstrom. These configurations include 305 atoms and 842 valence electrons, justifying the use of a semi-empirical approach. The electronic spectrum is then calculated using the INDO/CIS method. The solvatochromic shift of fluorescein in water, compared with in isolation, is calculated using the discrete and explicit solvent models. The use of electrostatically embedded explicit water molecules, in INDO/CIS calculations, gives a good description of the spectral shift of the fluorescein dianion in aqueous environment. The results are verified to converge both statistically and with respect to the number of explicit solvent molecules used.

Spectroscopic and theoretical studies of some N-methoxy-N-methyl-2-[(4 `-substituted) phenylsulfonyl]propanamides

The analysis of the IR carbonyl band of the N-methoxy-N-methyl-2-[(4`-substituted)phenyisulfonyI]propanamides Y-PhSO(2)CH(Me)C(O)N(OMe)Me (Y = OMe 1, Me 2, H 3, Cl 4, NO(2) 5), supported by HF/6-31G(d,p) calculations of 3,, indicated the existence of two gauche conformers (g(1) and g(2)), the g, conformer being the most stable and the least polar one (in gas phase and in solution). Both conformers are present in solution of the non polar solvent (CCl(4)) for 1-5 and in solution of the more polar solvents (CHCl(3)) for 1. 4, 5 and (CH(2)Cl(2)) for 5, while only the g(1) conformer is present in solution of the most polar solvent (CH(3)CN) for 1-5. The g, and g2 conformers correspond to the enantiomeric pairs of diastereomers (diast(1) and diast(2)) whose relative configurations are [C(3)(R)N(R)]/[C(3)(s)N(s)] and [C(3)(R)N(s)]/[C(3)(s)N(R)], respectively. The computed carbonyl frequencies for g(1) (diast(1)) and g(2) (diast(2)) stereoisomers of3 match well the experimental values. The NBO analysis, for 3 shows the important role of the orbital interactions in conformer stabilization and the overall balance of these interactions corroborates that the g, conformer is more stable than the 92 one. The observed abnormal solvent effect on the relative intensities of the carbonyl doublet components is attributed to the molecular crowding in the g2 conformer which hinders its solvation in comparison to the g, conformer (diast(1)). X-ray single crystal analysis performed for 3 shows the existence Of two 92, and g(1b) conformers of diastereomers (diast2, and diast(1b)) whose absolute configurations are [C(3)(R)N(s)] and [C(3)(R)N(R)], respectively. The larger population and. thus, the larger stabilization of the g(2), conformer over the gib form in the crystals may be associated with a larger energy gain deriving from dipole moment coupling in the former conformer along with a series of C-H center dot center dot center dot O electrostatic and hydrogen bond interactions, (C) 2009 Elsevier B.V. All rights reserved.

Conformational preferences for some 2-substituted N-methoxy-N-methylacetamides through spectroscopic and theoretical studies

The analysis of the IR carbonyl band of the 2-substituted N-methoxy-N-methylacetamides Y-CH(2)C(O)-N(OMe)Me (Y = F1, OMe 2, OPh 3, Cl 4), supported by B3LYP/6-311++G(3df, 3pd) calculations along with the NBO analysis for 1-4, indicated the existence of cis-gauche conformers i.e. (c) and (g) for 1 and 3, (c(1), c(2)) and (g(1), g(2)) for 2, and (c) and (g(1), g(2)) for 4. In the gas phase, the g conformer population prevails over the c one, for 1 and 3, the (c(1) + c(2)) population prevails over the (g(1) + g(2)) one for 2, and the (g(1) + g(2)) conformer population is more abundant than (c) one for 4. In n-hexane solution, the cis conformer is more abundant for 1-3. The occurrence of Fermi resonance in the nu(CO) region, in n-hexane, precludes the estimative of relative populations of the (c, g(1), g(2)) conformers for 4. The SCI-PCM calculations agree with the solvent effect on the nu(CO) band component relative intensities for 1-3. NBO analysis showed that the n(N) -> pi.(CO), orbital interaction is the main factor which stabilizes the gauche (g, g(1), g(2)) conformers for 1-4 into a larger extent relative to the cis (c, c(1), c(2)) ones. The n(y) -> pi(.)(Co,) sigma(C-Y) -> pi.(CO,) pi(CO) -> sigma(C-Y) and 7co orbital interactions still contribute, but into a minor extent for the stabilization of the gauche conformers relative to the cis ones. The existence of some pyramidalization at the nitrogen atom of the Weinreb amides 1-4 is responsible for the occurrence of Y(delta)-(4)center dot center dot center dot O(delta)-(9) and Y(delta)-(4)center dot center dot center dot N(delta)-(7) short contacts in the gauche (g, g(1), g(2)) conformers, which originates strong repulsive Coulombic interactions, acting in opposition to the large orbital stabilization of the gauche conformer with respect to the cis one. Therefore, a delicate balance of the Coulombic and orbital interactions seems to be responsible for the observed stabilization of the gauche (g, g(1), g(2)) and cis (c, c(1), c(2)) conformers, both in the gas phase and in the solution for 1-4. However, the cis conformer predominance, in non polar solvents, for the 2-substituted N-methoxy-N-methyl acetamides 1-3, bearing in a first raw (fluorine and oxygen) atoms, is in the opposite direction to the gauche conformer preference for the corresponding 2-substituted N,N-dialkyl-acetamides. (C) 2010 Elsevier B.V. All rights reserved.

Density functional study of the 5-methylcytosine tautomers

B3LYP/6-31++G** calculations to study seven tautomers of 5-methylcytosine in aqueous media have been carried out. Optimized geometries and relative stabilities for the different tautomers have been calculated in the gas phase, including interaction with two discrete water molecules and taking into account the solvent effects by using the self-consistent reaction field theory. The role of specific and bulk contributions of solvent effect on the observable properties of the 5-methylcytosine is clarified. The amino-oxo form is the most abundant tautomer in aqueous media. A reaction pathway connecting amino-oxo and amino-hydroxy forms along the corresponding transition structures has been characterized. Good agreement between theoretical and available experimental results of harmonic vibration frequencies is found. (C) 2001 Elsevier B.V. B.V. All rights reserved.

The Mycobacterium tuberculosis Rv2540c DNA sequence encodes a bifunctional chorismate synthase

Background. The emergence of multi- and extensively-drug resistant Mycobacterium tuberculosis strains has created an urgent need for new agents to treat tuberculosis (TB). The enzymes of shikimate pathway are attractive targets to the development of antitubercular agents because it is essential for M. tuberculosis and is absent from humans. Chorismate synthase (CS) is the seventh enzyme of this route and catalyzes the NADH- and FMN-dependent synthesis of chorismate, a precursor of aromatic amino acids, naphthoquinones, menaquinones, and mycobactins. Although the M. tuberculosis Rv2540c (aroF) sequence has been annotated to encode a chorismate synthase, there has been no report on its correct assignment and functional characterization of its protein product. Results. In the present work, we describe DNA amplification of aroF-encoded CS from M. tuberculosis (MtCS), molecular cloning, protein expression, and purification to homogeneity. N-terminal amino acid sequencing, mass spectrometry and gel filtration chromatography were employed to determine identity, subunit molecular weight and oligomeric state in solution of homogeneous recombinant MtCS. The bifunctionality of MtCS was determined by measurements of both chorismate synthase and NADH:FMN oxidoreductase activities. The flavin reductase activity was characterized, showing the existence of a complex between FMN ox and MtCS. FMNox and NADH equilibrium binding was measured. Primary deuterium, solvent and multiple kinetic isotope effects are described and suggest distinct steps for hydride and proton transfers, with the former being more rate-limiting. Conclusion. This is the first report showing that a bacterial CS is bifunctional. Primary deuterium kinetic isotope effects show that C4-proS hydrogen is being transferred during the reduction of FMNox by NADH and that hydride transfer contributes significantly to the rate-limiting step of FMN reduction reaction. Solvent kinetic isotope effects and proton inventory results indicate that proton transfer from solvent partially limits the rate of FMN reduction and that a single proton transfer gives rise to the observed solvent isotope effect. Multiple isotope effects suggest a stepwise mechanism for the reduction of FMNox. The results on enzyme kinetics described here provide evidence for the mode of action of MtCS and should thus pave the way for the rational design of antitubercular agents. © 2008 Ely et al; licensee BioMed Central Ltd.

Crystallographic portrayal of different conformational states of a Lys49 phospholipase A2 homologue: Insights into structural determinants for myotoxicity and dimeric configuration

Catalytically inactive phospholipase A2 (PLA2) homologues play key roles in the pathogenesis induced by snake envenomation, causing extensive tissue damage via a mechanism still unknown. Although, the amino acid residues directly involved in catalysis are conserved, the substitution of Asp49 by Arg/Lys/Gln or Ser prevents the binding of the essential calcium ion and hence these proteins are incapable of hydrolyzing phospholipids. In this work, the crystal structure of a Lys49-PLA2 homologue from Bothrops brazili (MTX-II) was solved in two conformational states: (a) native, with Lys49 singly coordinated by the backbone oxygen atom of Val31 and (b) complexed with tetraethylene glycol (TTEG). Interestingly, the TTEG molecule was observed in two different coordination cages depending on the orientation of the nominal calcium-binding loop and of the residue Lys49. These structural observations indicate a direct role for the residue Lys49 in the functioning of a catalytically inactive PLA2 homologue suggesting a contribution of the active site-like region in the expression of pharmacological effects such as myotoxicity and edema formation. Despite the several crystal structures of Lys49-PLA2 homologues already determined, their biological assembly remains controversial with two possible conformations. The extended dimer with the hydrophobic channel exposed to the solvent and the compact dimer in which the active site-like region is occluded by the dimeric interface. In the MTX-II crystal packing analysis was found only the extended dimer as a possible stable quaternary arrangement. © 2012 Elsevier B.V.

Planejamento de novos derivados da Edaravona

Um estudo teórico detalhado da edaravona foi realizado usando o método TFD (Teoria do Funcional da Densidade) com o conjunto de base B3LYP/6-31G*, com o objetivo de esclarecer o mecanismo sequestraste de radical e a influência do tautomerismo da edaravona em condições ácidas, importantes para auxiliar na elucidação do mecanismo de sequestro de radical pela edaravona no processo isquêmico. Em estudos teóricos anteriores, a tautomerização, o efeito do solvente e a abstração do elétron não foram considerados. Em nosso estudo, a estabilidade e reatividade foram determinadas através de parâmetros geométricos e energéticos. Os cálculos foram realizados em fase gasosa e o método PCM em fase aquosa e metanólica. As condições ácidas ou básicas foram consideradas pelas energias de dissociação e protonação, com consequente formação das formas aniônicas e catiônicas, respectivamente. As propriedades antioxidantes foram calculadas com base nos valores de HOMO, potencial de ionização (PI) e energia de dissociação da ligação (EDL). Os valores de HOMO e PI mostraram que o tautomero N-H é melhor antioxidante pela saída do elétron e os valores de EDL mostraram que o tautomero O-H é melhor antioxidante pela saída de hidrogênio. A protonação é mais favorecida termodinamicamente que a desprotonação. Além disso, a energia de protonação explica, teoricamente, a baixa diferença entre a protonação do N e O. O efeito do solvente diminuiu as barreiras de energias para a isomerização nos tautomeros O-H ou N-H. Adicionalmente, foram avaliados três derivados de pirazolonas com suas capacidades antioxidantes comparadas com a edaravona, em um esforço para identificar o farmacóforo antioxidante. A capacidade antioxidante da antipirina, dipirona, fenilbutazona e edaravona foram determinadas pela inibição de dois radicais estáveis DPPH e ABTS. A edaravona e fenilbutazona foram as mais potente para a inibição de radical DPPH e ABTS que a dipirona, enquanto que a antipirina não mostrou atividade em todas as concentrações analisadas. Simultaneamente, o método TFD previu o valor do potencial antioxidante in silico para explicar a relação de estrutura-atividade (REA). Além disso, foram calculados seus valores de HOMO, EDL X-H, e energia de estabilização (DEiso). Todos os cálculos foram executados usando o programa Gaussian 03, Hyperchem 7.5, e ChemOffice 2005. Os resultados mostraram que o derivado com C-H na posição C-4 aumentou a abstração do elétron ou de hidrogênio. Finalmente, a estratégia geral empregada para planejar os novos derivados se baseou nos estudos de dezoito derivados da edaravona descritos na literatura. Os estudos de REA classificaram os compostos em três grupos, como mais ativo, ativo e menos ativo. Nove derivados foram planejados a partir do composto mais ativo selecionado pelo estudo teórico.