Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

Molecular interaction between perthiolated [beta]-cyclodextrin (CD) and the guests molecules adamantaneacetic acid (AD) and ferroceneacetic acid (FC); and the effect of the interaction on the electron transition of CD anchored particles

The molecular interactions between the host molecule, perthiolated beta-cyclodextrin (CD), and the guest molecules, adamantaneacetic acid (AD) and ferroceneacetic acid (FC), have been inestigated theoretically in both the gas and aqueous phases. The major computations have been carried out at the theoretical levels, RHF/6-31G and B3LYP/6- 31G. MP2 electronic energies were also computed based at the geometries optimized by both the RHF and B3LYP methods in the gas phase to establish a better estimate of the correlation effect. The solvent phase computations were completed at the RHF/6-31G and B3LYP/6-31G levels using the PCM model. The most stable structures optimized in gas phase by both the RHF and B3LYP methods were used for the computations in solution. A method to systematically manipulate the relative position and orientation between the interacting molecules is proposed. In the gas phase, six trials with different host-guest relative positions and orientations were completed successfully with the B3LYP method for both the CD-AD and CD-FC complexes. Only four trials were completed with RHF method. In the gas phase, the best results from the RHF method gives for the association Gibbs free energy (ΔG°) values equal to -32.21kj/mol for CD-AD and -25.73kj/mol for CD-FC. And the best results from the B3LYP method have ΔG° equal to -47.57kj/mol for CD-AD and -41.09kj/mol for CD-FC. The MP2 correction significantly lowers ΔG° based on the geometries from both methods. For the RHF structure, the MP2 computations lowered ΔG° to -60.64kj/mol for CD-AD and -54.10 for CD-FC. For the structure from the B3LYP method, it was reduced to -59.87 kj/mol for CD-AD and -54.84 kj/mol for CDFC. The RHF solvent phase calculations yielded following results: ΔG°(aq) equals 107.2kj/mol for CD-AD and 111.4kj/mol for CD-FC. Compared with the results from the RHF method, the B3LYP method provided clearly better solvent phase results with ΔG° (aq) equal to 38.64kj/mol for CD-AD and 39.61kj/mol for CD-FC. These results qualitatively explain the experimental observations. However quantitatively they are in poor agreement with the experimental values available in the literature and those recently published by Liu et al. And the reason is believed to be omission of hydrophobic contribution to the association. Determining the global geometrical minima for these very large systems was very difficult and computationally time consuming, but after a very thorough search, these were identified. A relevant result of this search is that when the complexes, CD-AD and CD-FC, are formed, the AD and FC molecules are only partially embedded inside the CD cavity. The totally embedded complexes were found to have significantly higher energies. The semiempirical method, ZINDO, was employed to investigate the effect of complexation on the first electronic excitation of CD anchored to a metal nano-particle. The computational results revealed that after complexation to FC, the transition intensity declines to about 25% of the original value, and after complexation with AD, the intensity drops almost 50%. The tighter binding and transition intensity of CD-AD qualitatively agrees with the experimental result that the addition of AD to a solution of CD and FC restores the fluorescence of CD that was quenched by the addition of FC. A method to evaluate the “hydrophobic force” effect is proposed for future work.

An examination of the metabolic processes underpinning critical swimming in Atlantic cod (Gadus morhua L.) using in vivo 31P-NMR spectroscopy.

Traditionally, critical swimming speed has been defined as the speed when a fish can no longer propel itself forward, and is exhausted. To gain a better understanding of the metabolic processes at work during a U(crit) swim test, and that lead to fatigue, we developed a method using in vivo (31)P-NMR spectroscopy in combination with a Brett-type swim tunnel. Our data showed that a metabolic transition point is reached when the fish change from using steady state aerobic metabolism to non-steady state anaerobic metabolism, as indicated by a significant increase in inorganic phosphate levels from 0.3+/-0.3 to 9.5+/-3.4 mol g(-1), and a drop in intracellular pH from 7.48+/-0.03 to 6.81+/-0.05 in muscle. This coincides with the point when the fish change gait from subcarangiform swimming to kick-and-glide bursts. As the number of kicks increased, so too did the Pi concentration, and the pH(i) dropped. Both changes were maximal at U(crit). A significant drop in Gibbs free energy change of ATP hydrolysis from -55.6+/-1.4 to -49.8+/-0.7 kJ mol(-1) is argued to have been involved in fatigue. This confirms earlier findings that the traditional definition of U(crit), unlike other critical points that are typically marked by a transition from aerobic to anaerobic metabolism, is the point of complete exhaustion of both aerobic and anaerobic resources.

Deuterium isotope effect on the intramolecular electron transfer in Pseudomonas aeruginosa azurin

Intramolecular electron transfer in azurin in water and deuterium oxide has been studied over a broad temperature range. The kinetic deuterium isotope effect, kH/kD, is smaller than unity (0.7 at 298 K), primarily caused by the different activation entropies in water (−56.5 J K−1 mol−1) and in deuterium oxide (−35.7 J K−1 mol−1). This difference suggests a role for distinct protein solvation in the two media, which is supported by the results of voltammetric measurements: the reduction potential (E0′) of Cu2+/+ at 298 K is 10 mV more positive in D2O than in H2O. The temperature dependence of E0′ is also different, yielding entropy changes of −57 J K−1 mol−1 in water and −84 J K−1 mol−1 in deuterium oxide. The driving force difference of 10 mV is in keeping with the kinetic isotope effect, but the contribution to ΔS‡ from the temperature dependence of E0′ is positive rather than negative. Isotope effects are, however, also inherent in the nuclear reorganization Gibbs free energy and in the tunneling factor for the electron transfer process. A slightly larger thermal protein expansion in H2O than in D2O (0.001 nm K−1) is sufficient both to account for the activation entropy difference and to compensate for the different temperature dependencies of E0′. Thus, differences in driving force and thermal expansion appear as the most straightforward rationale for the observed isotope effect.

A thermodynamic coupling mechanism for GroEL-mediated unfolding.

Chaperonins prevent the aggregation of partially folded or misfolded forms of a protein and, thus, keep it competent for productive folding. It was suggested that GroEL, the chaperonin of Escherichia coli, exerts this function 1 unfolding such intermediates, presumably in a catalytic fashion. We investigated the kinetic mechanism of GroEL-induced protein unfolding by using a reduced and carbamidomethylated variant of RNase T1, RCAM-T1, as a substrate. RCAM-T1 cannot fold to completion, because the two disulfide bonds are missing, and it is, thus, a good model for long-lived folding intermediates. RCAM-T1 unfolds when GroEL is added, but GroEL does not change the microscopic rate constant of unfolding, ruling out that it catalyzes unfolding. GroEL unfolds RCAM-T1 because it binds with high affinity to the unfolded form of the protein and thereby shifts the overall equilibrium toward the unfolded state. GroEL can unfold a partially folded or misfolded intermediate by this thermodynamic coupling mechanism when the Gibbs free energy of the binding to GroEL is larger than the conformational stability of the intermediate and when the rate of its unfolding is high.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.

Estudo teórico de intermediários tetraédricos acidez / basicidade e estereosseletividade enzimáticos

The present work aimed first, the theoretical study of tetrahedral intermediate stability formed from carbonyl addition reactions using the second (MP2) and third (MP3) order Møller–Plesset perturbation theory. Linear correlations between electronic energy difference of reactions with Wiberg Indexes and C-O bond lengths were obtained, and was observed that the stability of adducts formed depends directly of electronic density involved between these atoms. The knowing of electronic parameters of these structures has an important hole due to the large use on reactions that in his course forms this tetrahedral intermediate. Employing the ONIOM (B3LYP:AMBER) methodology, was evaluated the stereoselectivity of a enzymatic reaction between CAL B enzyme and a long chain ester. In this study, were obtained the electronic energies of ground state and intermediate state of transesterification rate-determing step from two possible proquirals faces Re and Si. The objective was study the enantioselectivity of CAL B and rationalizes it using quantum theory of atoms in molecules (QTAIM). A theoretical study employing inorganic compounds was performed using ab initio CBS-QB3 method aiming to find a link between thermodynamic and equilibrium involving acids and bases. The results observed showed an excellent relationship between difference in Gibbs free energy, ΔG of acid dissociation reaction and ΔG of hydrolysis reaction of the corresponding conjugate base. It was also observed, a relationship between ΔG of hydrolysis reaction of conjugate acids and their corresponding atomic radius showing that stability plays an important role in hydrolysis reactions. The importance of solvation in acid/base behavior when compared to theoretical and experimental ΔG´s also was evaluated.

Síntese e caracterização de materiais mesoporosos para a captura de CO2: influência do óxido de níquel

Several materials are currently under study for the CO2 capture process, like the metal oxides and mixed metal oxides, zeolites, carbonaceous materials, metal-organic frameworks (MOF's) organosilica and modified silica surfaces. In this work, evaluated the adsorption capacity of CO2 in mesoporous materials of different structures, such as MCM-48 and SBA- 15 without impregnating and impregnated with nickel in the proportions 5 %, 10 % and 20 % (m/m), known as 5Ni-MCM-48, 10Ni-MCM-48, 20Ni-MCM-48 and 5Ni-SBA-15, 10NiSBA-15, 20Ni-SBA-15. The materials were characterized by means of X-ray diffraction (XRD), thermal analysis (TG and DTG), Fourier transform infrared spectroscopy (FT-IR), N2 adsorption and desorption (BET) and scanning electron microscopy (SEM) with EDS. The adsorption process was performed varying the pressure of 100 - 4000 kPa and keeping the temperature constant and equal to 298 K. At a pressure of 100 kPa, higher concentrations of adsorption occurred for the materials 5Ni-MCM-48 (0.795 mmol g-1 ) and SBA-15 (0.914 mmol g-1 ) is not impregnated, and at a pressure of 4000 kPa for MCM-48 materials (14.89 mmol g-1) and SBA-15 (9.97 mmol g-1) not impregnated. The results showed that the adsorption capacity varies positively with the specific area, however, has a direct dependency on the type and geometry of the porous structure of channels. The data were fitted using the Langmuir and Freundlich models and were evaluated thermodynamic parameters Gibbs free energy and entropy of the adsorption system

Adsorção de cromo (VI) em leito fixo utilizando partículas de areia recobertas com Quitosana

Neste trabalho foi avaliado o processo de recobrimento de partículas de areia com quitosana utilizando a técnica dip-coating, e analisado o emprego destas partículas como recheio de uma coluna de leito fixo no processo de adsorção de cromo (VI) em solução aquosa. A quitosana foi obtida a partir de resíduos de camarão e caracterizada. O estudo avaliou a influência do tamanho das partículas e da concentração da solução de quitosana no recobrimento das partículas de areia. Foram avaliados parâmetros termodinâmicos, isotermas de equilíbrio e parâmetros relacionados ao funcionamento do leito (vazão e pH da solução, diâmetro de partícula) para o processo de adsorção de cromo (VI) em solução aquosa. No recobrimento das partículas de areia, o tamanho não teve significância na resposta, enquanto a concentração da solução de recobrimento mostrou ter grande influência sobre o resultado, sendo que a menor concentração de quitosana dentro da faixa estudada (0,5% p/v) apresentou o melhor desempenho. O processo de cura física para o recobrimento das partículas de areia mostrou melhor desempenho para a adsorção de cromo (VI) em leito fixo frente ao processo físico-químico. As análises de superfície (MEV) e de difração de raio-X (EDX) comprovaram a mudança na superfície das partículas recobertas e a presença de cromo após a adsorção. O modelo de Sips foi o que melhor representou os dados experimentais de equilíbrio, com R2 >0,99% e EMR<3,5%, sendo que a capacidade máxima de adsorção foi de 46,93 mg g-1 obtida a 298 K. O processo se mostrou espontâneo, exotérmico e favorável, com valores de -4,49 a -4,66kJ mol-1 para e energia livre de Gibbs, -5,97kJ mol-1 para a variação de entalpia e -5,17x10-3 kJ mol-1K -1 para variação de entropia. O aumento do pH diminuiu a adsorção de cromo (VI), sendo que a melhor resposta foi obtida em pH 3, sendo que o diâmetro de partícula não teve efeito significante. O estudo da vazão da solução de cromo (VI) no desempenho do leito mostrou que no menor valor (2,5 mL min-1 ) a concentração de saída do leito foi próxima a zero, mantendo-se com concentrações de saída abaixo de 20% por pelo menos 20 min. O estudo da dessorção do leito mostrou que após cinco ciclos de trabalho o leito manteve 86% da capacidade de adsorção, com taxas de recuperação do cromo no processo de dessorção maiores que 95%.

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

Interaction of Multiwalled Carbon Nanotube Produces Structural Alteration of DNA

Multiwalled carbon nanotube (MWCNT) has been found to produce structural changes in Calf Thymus-DNA (CT-DNA). The interaction or binding of the multi-walled carbon nanotubes (MWCNT) was investigated in order to discover if it brings about any significant changes of the DNA double helix using CD spectra of the CT-DNA at two concentration levels of MWCNT representing an increasing MWCNT/DNA molar ratio. In addition, spectrophotometric titrations between MWCNT and CT-DNA were carried out in order to utilize spectral changes as a means of detecting specific binding modes of either intercalation or degradation of DNA. Interactions of MWCNT induced significant changes in the CD spectra of the B-form of natural DNA. The intensities of the positive CD band at 280 nm decreased significantly. This decrease was found to be concentration-dependent. Following spectrophotometric titrations; specific subtle conformational changes were observed with a molar ratio combination of 2:1 between MWCNT and CT-DNA and these were characterized by a formation constant of the order of 103 M-1 and a negative Gibbs free energy suggesting that MWCNT avidly binds to DNA. Thermodynamic considerations revealed that electrostatic interactions between the DNA base pairs and the MWCNT are taking place accounting for the negative free energy change, positive enthalpy change with a small entropy change. The results obtained in the study of the binding interactions of MWCNT with DNA confirm that a cytogenetic effect of MWCNT with DNA is a possibility in vivo.

Measurements of activity coefficients at infinite dilution of organic solutes and water on polar imidazolium-based ionic liquids

The activity coefficients at infinite dilution, gamma(infinity)(13), of 55 organic solutes and water in three ionic liquids with the common cation 1-butyl-3-methylimidazolium and the polar anions Cl--,Cl- [CH3SO3](-) and [(CH3)(2)PO4](-), were determined by (gas + liquid) chromatography at four temperatures in the range (358.15 to 388.15) K for alcohols and water, and T = (398.15 to 428.15) K for the other organic solutes including alkanes, cycloalkanes, alkenes, cycloalkenes, alkynes, ketones, ethers, cyclic ethers, aromatic hydrocarbons, esters, butyraldehyde, acetonitrile, pyridine, 1-nitropropane and thiophene. From the experimental gamma(infinity)(13) values, the partial molar excess Gibbs free energy, (G) over bar (E infinity)(m), enthalpy (H) over bar (E infinity)(m), and entropy (S) over bar (E infinity)(m), at infinite dilution, were estimated in order to provide more information about the interactions between the solutes and the ILs. Moreover, densities were measured and (gas + liquid) partition coefficients (KL) calculated. Selectivities at infinite dilution for some separation problems such as octane/benzene, cyclohexane/benzene and cyclohexane/thiophene were calculated using the measured gamma(infinity)(13), and compared with literature values for N-methyl-2-pyrrolidinone (NMP), sulfolane, and other ionic liquids with a common cation or anion of the ILs here studied. From the obtained infinite dilution selectivities and capacities, it can be concluded that the ILs studied may replace conventional entrainers applied for the separation processes of aliphatic/aromatic hydrocarbons.

Gibbs Energy of Formation of Cobalt Divanadium Tetroxide

The Gibbs energy of formation of V2O3-saturated spinel CoV2O4 has been measured in the temperature range 900–1700 K using a solid state galvanic cell, which can be represented as Pt, Co + CoV2O4 + V2O3/(CaO) ZrO2/Co + CoO, Pt. The standard free energy of formation of cobalt vanadite from component oxides can be represented as CoO (rs) + V2O3 (cor) → CoV2O4 (sp), ΔG° = −30,125 − 5.06T (± 150) J mole−1. Cation mixing on crystallographically nonequivalent sites of the spinel is responsible for the decrease in free energy with increasing temperature. A correlation between “second law” entropies of formation of cubic 2–3 spinels from component oxides with rock salt and corundum structures and cation distribution is presented. Based on the information obtained in this study and trends in the stability of aluminate and chromite spinels, it can be deduced that copper vanadite is unstable.

Phase equilibria in the system NiO-CaO-SiO2 and Gibbs energy of formation of CaNiSi2O6

Phase relations in the pseudoternary system NiO-CaO-SiO2 at 1373 K are established. The coexisting phases are identified by X-ray diffraction and energy-dispersive X-ray analysis of equilibrated samples. There is only one quaternary oxide CaNiSi2O6 with clinopyroxene structure. The Gibbs energy of formation of CaNiSi2O6 is measured using a solid state galvanic cell incorporating stabilized zirconia as the solid electrolyte in the temperature range of 1000 to 1400 K:Pt, Ni + SiO2 + CaSiO3 + CaNiSi2O6 \ (Y2O3)ZrO2 \ Ni + NiO, Pt From the electromotive force (emf) of the cell, the Gibbs energy of formation of CaNiSi2O6 from NiO, SiO2, and CaSiO3 is obtained. To derive the Gibbs energy of formation of the quaternary oxide from component binary oxides, the free energy of formation of CaSiO, is determined separately using a solid state cell based on single crystal CaF2 as the electrolyte: Pt, O-2, CaO + CaF2 \ CaF2 \ CaSiO3 + SiO2 + CaF2, O-2, Pt The results can be expressed by the following equations: NiO (r.s) + CaO (r.s) + 2SiO(2) (qz) --> CaNiSi2O6 (pyr) Delta G degrees = -115,700 + 10.63T (+/-100) J mol(-1) CaO (r.s) + SiO2 (qz) --> CaSiO3 (wol) Delta G degrees = -90,030 -0.61T (+/-60) J mol(-1).

Phase coexistence of cluster crystals: beyond the Gibbs phase rule

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants.