Computational studies on glyceraldehyde and glycine Maillard reaction-II

Density-functional theory (BLYP/6-31G*) and semiempirical theory was used to study compounds related to a strained cyclic hexaalkyne that was postulated by Staab et al. as a trimeric intermediate in the copper-mediated oxidative coupling of 2,2′-diethynylbiphenyl. Among the products of this reaction is an ortho-arene cyclyne. It was found that the hexaalkyne→ortho-arene cyclyne rearrangement is quite exothermic (ca. 130 kcal/mol). Of the semiempirical methods tested, the AM1 results were closer than PM3 or MNDO to BLYP/6-31G*. For the intermediate, two conformations were found: the expected D₃-symmetric conformation and an unexpected C₂-symmetric conformation. Structures and energies of hexamethyl- and dodecamethyl-substituted analogues of the intermediate and product were also calculated at the AM1 level, in order to test whether steric effects might stabilize the intermediate, relative to the product. It was found that even very sterically crowded products were much more stable than the relevant intermediates. This suggested that the large exothermicity of the hexaalkyne→ortho-arene cyclyne rearrangement might be used to drive formation of interesting strained products. An example used to illustrate this is the rearrangement of a binaphthyl analogue of Staab's biphenyl intermediate. BLYP and AM1 calculations suggest that the rearrangement of this intermediate would be exothermic by ca. 105–110 kcal/mol, despite strain induced by formation of three [5]helicene moieties in the product.     

Density functional computational studies on the glucose and glycine Maillard reaction: Formation of the Amadori rearrangement products

Theoretical energy changes of various intermediates leading to the formation of the Amadori rearrangement products (ARPs) under different mechanistic assumptions have been calculated, by using open chain glucose (O‐Glu)/closed chain glucose (A‐Glu and B‐Glu) and glycine (Gly) as a model for the Maillard reaction. Density functional theory (DFT) computations have been applied on the proposed mechanisms under different pH conditions. Thus, the possibility of the formation of different compounds and electronic energy changes for different steps in the proposed mechanisms has been evaluated. B‐Glu has been found to be more efficient than A‐Glu, and A‐Glu has been found more efficient than O‐Glu in the reaction. The reaction under basic condition is the most favorable for the formation of ARPs. Other reaction pathways have been computed and discussed in this work. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

酶解小麦蛋白产物-还原糖美拉德反应的光谱研究

采用紫外-可见吸收光谱和荧光光谱研究了酶解小麦蛋白产物与还原糖不同加热条件下的美拉德反应及其产物.美拉德反应在紫外区240和294 nm产生两个特征峰,荧光的最大激发和发射波长为347和450 nm;随反应进行,紫外光吸收和荧光强度迅速增加,表明美拉德反应进入高级阶段,产生的糠醛类、呋喃酮类、吡喃酮类、噻吩类及噻唑类等小分子物质表现较大的积累速率.温度升高,强度增加速率增大.在较高温度时,紫外光吸收出现最大平稳值;荧光强度则到达最大值后降低,表明小分子物质间或与肽聚合生成大分子黑素类物质,小分子物质的积累表现消除速率,反应进入终级阶段.     

Effects of Maillard Reaction on Volatile Compounds and Antioxidant Capacity of Cat Food Attractant

In this study, self-made cat food attractant was prepared by Maillard reaction using hydrolysate of grass carp waste as raw material and glucose and cysteine hydrochloride as substrate. Its volatile compounds, antioxidant capacity, and pet palatability were investigated. The volatile compounds of attractants were analyzed using gas chromatography–mass spectrometry (GC-MS) which showed that alcohols and aldehydes were the most volatile in self-made attractants, accounting for 34.29% and 33.52%, respectively. Furthermore, Maillard reaction could significantly increase the antioxidant activity of self-made attractant, including scavenging activity on OH and DPPH free radicals as well as the chelating ability of Fe²⁺. The acceptance and palatability of two kinds of cat food were studied by adding 3% self-made or commercial attractants. The results of this study also found that both attractants could remarkably improve the intake rate of cat food. However, the self-made group was significantly less than the commercial group in first smell, first bite, and feeding rate, which might be because of the absence of umami ingredients and spices in self-made attractants.     

Conformational preferences of the aldol adducts of oxadiazinones. H NMR spectroscopy and computational studies of N4-methyl and N4-isopropyloxadiazinones

The conformational properties of the aldol adducts of some N₄-isopropyl-oxadiazinones have been investigated by ¹H NMR spectroscopy and computational studies. An earlier study of the syn-aldol adducts of N₄-methyl-oxadiazinone 2 led to the conclusion that the solution and solid state conformation of these compounds involve syn-parallel arrangement of the C₂- and N₃-carbonyls of the oxadiazinones. However, the synthesis and asymmetric aldol reactions of an N₃-hydrocinnamoyl-N₄-isopropyl-oxadiazinone 4 has yielded aldol adducts 5a-e in which the orientation of the C₂- and N₃-carbonyls are most likely in the anti-parallel arrangement. These aldol adducts have been studied by ¹H NMR spectroscopy and the shielding aspect observed clearly suggests the presence of the anti-parallel arrangement. The installment of a N₄-d₆-isopropyl group further confirmed this assertion. Computational studies support the conclusion that solution state conformation of the N₄-methyl and N₄-isopropyl-oxadiazinones involves anti-parallel carbonyls in contrast to the solid state evidence of the X-ray crystallographic data of oxadiazinone 2.     

Computational studies on the mechanism and kinetics of Cl reaction with C2H5I

The dual‐level direct kinetics method has been used to investigate the multichannel reactions of C₂H₅I + Cl. Three hydrogen abstraction channels and one displacement process are found for the title reaction. The calculation indicates that the hydrogen abstraction from ? CH₂? group is the dominant reaction channel, and the displacement process may be negligible because of the high barrier. The rate constants for individual reaction channels are calculated by the improved canonical variational transition‐state theory with small‐curvature tunneling correction over the temperature range of 220–1500 K. Our results show that the tunneling correction plays an important role in the rate constant calculation in the low‐temperature range. Agreement between the calculated and experimental data available is good. The Arrhenius expression k(T) = 2.33 × 10^?16 T^1.83 exp(?185.01/T) over a wide temperature range is obtained. Furthermore, the kinetic isotope effects for the reaction C₂H₅I + Cl are estimated so as to provide theoretical estimation for future laboratory investigation. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

On Mercury(I) Oxo Compounds — Quasi‐Relativistic Computational and Experimental Studies

About 60 molecular species composed of up to 10 mercury atoms and of oxygen atoms and/or of some other elements or groups (such as halogen, OH₂, OH, H, alkali, NO₃) have been investigated quantum chemically. Different density functional approaches and the ab initio SCF‐MP2 method were applied, comparing different basis sets and different atomic core sizes. It is important not to treat the Hg 5s, p, d as inactive core shells, and to use sufficiently many polarization functions. The shape of the 〉O‐Hg‐Hg‐O〈 units is not favorable concerning the formation of lattices composed of Hg^I, O and OH only. Despite its bulkiness, the OHgHgO units can easily come into contact with each other and then disproportionate. This is prevented in the so‐called ternary M‐Hg^I oxides by the embedded oxometallate (oxoacidic) anions. Furthermore, the Hg^I and Hg^II oxide bond energies are less favorable towards the stability of Hg^I oxo compounds, as compared to Hg halidic or oxoacidic compounds. Both points are not promising concerning the search for Hg^I oxides/hydroxides, although the preparation of such compounds, including spacer groups, by topochemical reactions can still not be excluded. So far, experimental efforts towards the synthesis of such a new class of compounds have only demonstrated that Hg^II is strictly preferred over Hg^I in the formation of solids of binary Hg‐O or ternary A‐Hg‐O composition (A = electropositive metal such as alkali, in contrast to M = transition or semi‐metal). This is so even if compounds containing ‘electron rich Hg^δ— atoms’ (i.e. A‐Hg amalgams) are oxidized under mild conditions.     

Computational studies on density function theory for the bimolecular metathesis reaction CH3+HCl⇔︁CH4+Cl

Density function theory has been applied to alkyl radical reaction to get helpful data for its geometric parameters, energy, and vibrational frequency compared with results obtained by ab initio methods and experimental values. The geometry optimization of the transition state, the precursor complex and the successor complex were performed at the 6‐311G* basis set level. The transition state of the CH₄Cl system of the reaction computed was in agreement with the prediction of Benson. From analysis of the vibration frequency and the net charge on the atom of the precursor complex, transition state, successor complex and the isolated state, the reaction mechanism was derived which we complicated with the bond‐rupture electron‐transfer theory. The atom H in molecule HCl attacks the atom C, forming a transition state via the precursor complex and the electron transfer happens in the precursor complex. The reaction rate of the electron transfer determines the rate of the whole reaction to a certain extent, and active energy, electronic coupling matrix element, and reorganization energy were obtained. © 2001 Wiley Periodicals, Inc. Int J Quantum Chem, 2001     

Pincer‐Type Heck Catalysts and Mechanisms Based on PdIV Intermediates: A Computational Study

Pincer‐type palladium complexes are among the most active Heck catalysts. Due to their exceptionally high thermal stability and the fact that they contain Pd^II centers, controversial Pd^II/Pd^IV cycles have been often proposed as potential catalytic mechanisms. However, pincer‐type Pd^IV intermediates have never been experimentally observed, and computational studies to support the proposed Pd^II/Pd^IV mechanisms with pincer‐type catalysts have never been carried out. In this computational study the feasibility of potential catalytic cycles involving Pd^IV intermediates was explored. Density functional calculations were performed on experimentally applied aminophosphine‐, phosphine‐, and phosphite‐based pincer‐type Heck catalysts with styrene and phenyl bromide as substrates and (E)‐stilbene as coupling product. The potential‐energy surfaces were calculated in dimethylformamide (DMF) as solvent and demonstrate that Pd^II/Pd^IV mechanisms are thermally accessible and thus a true alternative to formation of palladium nanoparticles. Initial reaction steps of the lowest energy path of the catalytic cycle of the Heck reaction include dissociation of the chloride ligands from the neutral pincer complexes [{2,6‐C₆H₃(XPR₂)₂}Pd(Cl)] [X=NH, R=piperidinyl ( 1 a ); X=O, R=piperidinyl ( 1 b ); X=O, R=iPr ( 1 c ); X=CH₂, R=iPr ( 1 d )] to yield cationic, three‐coordinate, T‐shaped 14e^? palladium intermediates of type [{2,6‐C₆H₃(XPR₂)₂}Pd]⁺ ( 2 ). An alternative reaction path to generate complexes of type 2 (relevant for electron‐poor pincer complexes) includes initial coordination of styrene to 1 to yield styrene adducts [{2,6‐C₆H₃(XPR₂)₂}Pd(Cl)(CH₂?CHPh)] ( 4 ) and consecutive dissociation of the chloride ligand to yield cationic square‐planar styrene complexes [{2,6‐C₆H₃(XPR₂)₂}Pd(CH₂?CHPh)]⁺ ( 6 ) and styrene. Cationic styrene adducts of type 6 were additionally found to be the resting states of the catalytic reaction. However, oxidative addition of phenyl bromide to 2 result in pentacoordinate Pd^IV complexes of type [{2,6‐C₆H₃(XPR₂)₂}Pd(Br)(C₆H₅)]⁺ ( 11 ), which subsequently coordinate styrene (in trans position relative to the phenyl unit of the pincer cores) to yield hexacoordinate phenyl styrene complexes [{2,6‐C₆H₃(XPR₂)₂}Pd(Br)(C₆H₅)(CH₂?CHPh)]⁺ ( 12 ). Migration of the phenyl ligand to the olefinic bond gives cationic, pentacoordinate phenylethenyl complexes [{2,6‐C₆H₃(XPR₂)₂}Pd(Br)(CHPhCH₂Ph)]⁺ ( 13 ). Subsequent β‐hydride elimination induces direct HBr liberation to yield cationic, square‐planar (E)‐stilbene complexes with general formula [{2,6‐C₆H₃(XPR₂)₂}Pd(CHPh?CHPh)]⁺ ( 14 ). Subsequent liberation of (E)‐stilbene closes the catalytic cycle.     

Effect of glycine on the surface activity and micellar properties of <Emphasis Type="Italic">N</Emphasis>-decanoyl-<Emphasis Type="Italic">N</Emphasis>-methylglucamide

Surface activity, micelle formation and structure of N-decanoyl-N-methylglucamide in presence of increasing glycine concentrations have been investigated. Changes in the critical micelle concentration of the surfactant with the addition of the amino acid were examined by both surface tension and pyrene 1:3 ratio methods. The observed reduction in the critical micelle concentration was attributed to the structure making action of the additive. Micelle structure parameters were obtained as a function of the glycine concentration by using light scattering and fluorescence-quenching measurements. As regard to the effect of glycine of the micelle size, it was found that both the hydrodynamic radius and the micellar aggregation number increased with the glycine concentration. However, the observed micellar growth seems to be mainly determined by the increase in the aggregation number rather than by changes in the amount of water specifically associated to the micelle. This fact is also supported by the observed trend in the micellar partial specific volume as estimated by complementary density measurements.     

Computational and experimental studies on the mechanism of the photochemical carbonylation of group 6 Fischer carbene complexes

The photocarbonylation reaction of Group 6 Fischer carbene complexes has been studied by DFT and experimental procedures. The process occurs by intersystem crossing (ISC) from the lowest excited singlet state (S1) to the lowest triplet state (T1), the latter structure being decisive for the outcome of the reaction. Methylenepentacarbonylchromium(0) complexes, alkoxypentacarbonylchromium(0)carbene complexes, and alkoxyphosphinetetracarbonylchromium(0) carbene complexes have coordinatively unsaturated chromacyclopropanone T1 structures with a biradical character. The evolution of the metallacyclopropanone species occurs by a jump (spin inversion) to the S(0) hypersurface by coordination of a molecule of the solvent, leading to ketene-derived products in the presence of ketenophiles or reverting to the starting carbene complex in their absence. The T1 excited states obtained from methylenephosphinetetracarbonylchromium(0) complexes and pentacarbonyltungsten(0)carbene complexes are unable to produce the carbonylation. The reaction with ketenophiles is favored in coordinating solvents, which has been tested experimentally in the reaction of alkoxypentacarbonylchromium(0) complexes and imines.     

Computational studies on electron and proton transfer in phenol‐imidazole‐base triads

The electron and proton transfer in phenol‐imidazole‐base systems (base = NH₂^? or OH^?) were investigated by density‐functional theory calculations. In particular, the role of bridge imidazole on the electron and proton transfer was discussed in comparison with the phenol‐base systems (base = imidazole, H₂O, NH₃, OH^?, and NH₂^?). In the gas phase phenol‐imidazole‐base system, the hydrogen bonding between the phenol and the imidazole is classified as short strong hydrogen bonding, whereas that between the imidazole and the base is a conventional hydrogen bonding. The n value in spⁿ hybridization of the oxygen and carbon atoms of the phenolic CO sigma bond was found to be closely related to the CO bond length. From the potential energy surfaces without and with zero point energy correction, it can be concluded that the separated electron and proton transfer mechanism is suitable for the gas‐phase phenol‐imidazole‐base triads, in which the low‐barrier hydrogen bond is found and the delocalized phenolic proton can move freely in the single‐well potential. For the gas‐phase oxidized systems and all of the triads in water solvent, the homogeneous proton‐coupled electron transfer mechanism prevails. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Computational Study of the Substitution Effect on the Mechanism for the Aza‐ and Arsa‐Wittig Reactions

The aza‐ and arsa‐Wittig reactions HM=PH₃ + O=CHX → HM=CHX + O=PH₃ (M = N, As; X = H, F, Cl, Me, OMe, NMe₂, CMe₃) were examined using the density functional theory calculations. All of the structures were completely optimized at the B3LYP/6‐311++G** level of theory. The main finding of this work is that the difference between singlet‐triplet splitting of O=CHX and HM=PH₃ play an important role in determining the kinetic and thermodynamic stability of the aza‐ and arsa‐Wittig reactions. When HM=PH₃ with more ylidic character is utilized, the reaction has a smaller activation energy and a larger exothermicity.     

Ir‐Catalysed Asymmetric Allylic Substitutions with Cyclometalated (Phosphoramidite)Ir Complexes—Resting States,Catalytically Active (π‐Allyl)Ir Complexes and Computational Exploration

Mechanistic aspects of allylic substitutions with iridium catalysts derived from phosphoramidites by cyclometalation were investigated. The determination of resting states by ³¹P NMR spectroscopy led to the conclusion that the cyclometalation process is reversible. A novel, one‐pot procedure for the preparation of (π‐ allyl)Ir complexes was developed, and these complexes were characterised by X‐ray crystal structure analyses and spectral data. They are fully active catalysts of the allylic substitution reaction. DFT calculations on the allyl complexes, transition states of the allylic substitution and product olefin complexes gave further mechanistic insight.     

Real wave packet and flux analysis studies of the H + F2 → HF + F reaction

The H + F₂ → HF + F reaction on ground state potential energy surface is investigated using the quantum mechanical real wave packet and Flux analysis method based on centrifugal sudden approximation. The initial state selected reaction probabilities for total angular momentum J = 0 have been calculated by both methods while the probabilities for J > 0 have been calculated by Flux analysis method. The initial state selected reaction probabilities, integral cross sections and rate coefficients have been calculated for a broad range of collision energy. The results show a large rotational enhancement of the reaction probability. Some resonances were seen in the state‐to‐state reaction probabilities while state‐to‐all reaction probabilities and the reaction cross section do not manifest any oscillations and the initial state selected reaction rate constants are sensitive to the temperature. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Experimental and Computational Studies on the CH Amination Mechanism of Tetrahydrocarbazoles via Hydroperoxides

The acid‐catalyzed reactions of photochemically generated tetrahydrocarbazole peroxides with anilines have been studied experimentally and computationally to identify the underlying reaction mechanism. The kinetic data indicate a reaction order of one in the hydroperoxide and zero in the aniline. Computational investigations using density functional theory support the experimental findings and predict an initial tautomerization between an imine and enamine substructure of the primarily generated tetrahydrocarbazole peroxide to be the rate controlling step. The enamine tautomer then loses hydrogen peroxide upon protonation, generating a stabilized allylic carbocation that is reversibly trapped by solvent or aniline to form the isolated products.     

Configurational Simulations and Theoretical Calculations of Molecularly Imprinted Polymers of Histamine and 2‐(Trifluoromethyl)acrylic Acid Based on Computational Chemistry

In this article, we provide a theoretical discussion on the interactions between a template molecule and functional monomer(s) in the preparation of molecularly imprinted polymers (MIPs ). Density functional theory (DFT ) was used to compute the 3D structures, natural bond orbital, and binding energy in the template–monomer(s) complexes. Histamine (HA ) and 2‐(trifluoromethyl)acrylic acid (TFMAA ) were, respectively, selected as the template and the monomer. The computational process was performed according to B3LYP method with 6‐311 + (d,p) basis set under the different HA–TFMAA ratios from 1:1 to 1:10. The computational results show that the HA–TFMAA complex at the ratio of 1:5 yields no consequence. Furthermore, the HA–TFMAA complex at the ratio of 1:5 allowed the minimum binding energy and the steadiest condition, with four hydrogen bonds. Configurational simulations and theoretical calculations of the template–monomer complex can be used as an appropriate guiding tool for manufacturing MIPs with high specificity and selectivity, thereby avoiding repeated experiments and wastage of substantial time and investment.     

Evolution of a Multicomponent System: Computational and Mechanistic Studies on the Chemo‐ and Stereoselectivity of a Divergent Process

The evolution of a ternary molecular system (imine, diene and nitrile) is analyzed to disclose the pathways leading to a divergent synthetic outcome. The Lewis acid catalyzed reaction between cyclohexadiene, 2‐phenyl‐indol‐3‐one and acetonitrile yields the imino‐Diels–Alder adduct as the major product together with minor amounts of the Mannich–Ritter‐amidine product. The experimental and computational data show that the relative orientation of the initial reactants dictates the synthetic outcome. The exo approach between imine and diene leads to the Diels–Alder adduct in a concerted process, whereas the endo mode leads to a polarized intermediate, which is trapped by acetonitrile to yield the multicomponent adduct.     

Single crystal EPR and optical studies of paramagnetic ions doped zinc potassium phosphate hexahydrate—Part I: Cu(II)—a case of orthorhombic symmetry

Single crystal electron paramagnetic resonance (EPR) studies on Cu(II) doped zinc potassium phosphate hexahydrate (ZPPH) were carried out at room temperature. The angular variation spectra in the three orthogonal planes indicate that the paramagnetic impurity has entered the lattice substitutionally in place of Zn(II) and the spin Hamiltonian parameters calculated from these spectra are g(xx) = 2.188, g(yy) = 2.032, g(zz) = 2.373, Axx = 50 G, Ayy = 65.0 G and Azz = 80 G. The g and A tensors were coincident and these values matched fairly well with the values obtained from powder spectrum. The bonding parameters have also been calculated.