Theoretical study of the mechanism of acetaldehyde hydroxylation by compound I of CYP2E1

Recent experimental studies revealed that cytochrome P450 2E1 (CYP2E1) could metabolize not only ethanol but also its primary product, acetaldehyde, accompanying the well-known acetaldehyde dehydrogenases (ALDH) in the metabolism of acetaldehyde. Mechanistic aspects of acetaldehyde hydroxylation by Compound I model active species of CYP2E1 were investigated by means of B3LYP DFT calculations in the present paper. Our study results demonstrate that acetaldehyde hydroxylation by CYP2E1 is in accord with the effectively concerted mechanisms both on the high quartet spin state (HS) and on the low doublet spin state (LS). The rate-limiting step is H-abstraction, and the activation energy is about 11.7 approximately 14.0 kcal/mol on the quartet (doublet) reaction route, which is about one-half to one-third of that needed by methane hydroxylation. The phenomenon that the HS and LS reaction routes are both effectively concerted was shown for the first time to occur in trans-2-phenyl-iso-propylcyclopropane hydroxylation by Kumar et al. (see Figure 7 in the paper of Kumar, D.; de Visser, S. P.; Sharma, P. K.; Cohen, S.; Shaik, S. J. Am. Chem. Soc. 2004, 126, 1907) and was confirmed in our work of acetaldehyde hydroxylation by cytochrome P450. Theoretical exploration of the HS O-rebound barrier degradation is also presented in the present paper on the basis of Shaik's valence bond (VB) model.     

Synthesis and properties of a fluorene-capped isotruxene: a new unsymmetrical star-shaped pi-system

An unsymmetrical star-shaped pi-system (1) with an isotruxene core and three fluorene arms has been synthesized, and its photophysical, electrochemical, and thermochemical properties have been investigated and compared with the corresponding symmetrical truxene derivatives 2a-d (Kanibolotsky, A. L.; Berridge, R.; Skabara, P. J.; Perepichka, I. F.; Bradley, D. D. C.; Koeberg, M. J. Am. Chem. Soc. 2004, 126, 13695-13702). Stronger electronic couplings between the arms in 1 vs 2a-d lead to lower optical band gap, larger splitting in oxidation potentials, and superior thermostability. [structure: see text]     

Determinants of cyanuric acid and melamine assembly in water

While the recognition of cyanuric acid (CA) by melamine (M) and their derivatives has been known to occur in both water and organic solvents for some time, analysis of CA/M assembly in water has not been reported (Ranganathan, A.; Pedireddi, V. R.; Rao, C. N. R. J. Am. Chem. Soc.1999, 121, 1752-1753; Mathias, J. P.; Simanek, E. E.; Seto, C. T.; Whitesides, G. M. Macromol. Symp.1994, 77, 157-166; Zerkowski, J. A.; MacDonald, J. C.; Seto, C. T.; Wierda, D. A.; Whitesides, G. M. J. Am. Chem. Soc.1994, 116, 2382-2391; Mathias, J. P.; Seto, C. T.; Whitesides, G. M. Polym. Prepr.1993, 34, 92-93; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1993, 115, 905-916; Zerkowski, J. A.; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1992, 114, 5473-5475; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1990, 112, 6409-6411; Wang, Y.; Wei, B.; Wang, Q. J. Chem. Cryst.1990, 20, 79-84; ten Cate, M. G. J.; Huskens, J.; Crego-Calama, M.; Reinhoudt, D. N. Chem.-Eur. J.2004, 10, 3632-3639). We have examined assembly of CA/M, as well as assembly of soluble trivalent CA and M derivatives (TCA/TM), in aqueous solvent, using a combination of solution phase NMR, isothermal titration and differential scanning calorimetry (ITC/DSC), cryo-transmission electron microscopy (cryo-TEM), and synthetic chemistry. While the parent heterocycles coprecipitate in water, the trivalent system displays more controlled and cooperative assembly that occurs at lower concentrations than the parent and yields a stable nanoparticle suspension. The assembly of both parent and trivalent systems is rigorously 1:1 and proceeds as an exothermic, proton-transfer coupled process in neutral pH water. Though CA and M are considered canonical hydrogen-bonding motifs in organic solvents, we find that their assembly in water is driven in large part by enthalpically favorable surface-area burial, similar to what is observed with nucleic acid recognition. There are currently few synthetic systems capable of robust molecular recognition in water that do not rely on native recognition motifs, possibly due to an incomplete understanding of recognition processes in water. This study establishes a detailed conceptual framework for considering CA/M heterocycle recognition in water which enables the future design of molecular recognition systems that function in water.     

Competing mechanisms in the optically activated functionalization of the hydrogen-terminated Si(111) surface

Recent experiments have shown that organic monolayers on silicon surfaces can be formed through the optically activated surface reaction of H-terminated Si surfaces with terminally unsaturated organic molecules (Eves et al. J. Am. Chem. Soc. 2004, 126, 14318; Sun et al. J. Am. Chem. Soc. 2005, 127, 2514). Possible mechanisms for the formation of this monolayer involve the abstraction of a H atom either at the same attachment site of the molecule (Path A) or from a neighboring site (Path B). Using periodic Density Functional Theory calculations together with an efficient method for finding reaction pathways, we examine both optically activated reaction mechanisms for an alkene and an aldehyde reacting with H-Si(111). Our results show that while Path A is energetically more favorable its significant barrier is likely to limit its viability. Path B on the other hand encounters a much lower H atom abstraction barrier and appears to be more viable.     

New reaction simulator "LUMMOX" and its application for prediction of catalytic activities

We developed a new reaction simulator, "LUMMOX." It is an intermolecular interaction analyzer based on the theories of paired interacting orbitals (PIOs) and localized frontier orbitals (LFOs) that have been developed by Fujimoto et al. (Fukui, K.; Koga, N.; Fujimoto, H. J Am Chem Soc 1981, 103, 196; Fujimoto, H.; Koga, N.; Fukui, K. J Am Chem Soc 1981, 103, 7452; Fujimoto, H.; Satoh, S. J Phys Chem 1994, 98, 1436). LUMMOX runs on a Windows PC and displays graphic representation of orbital interactions. Prediction of activities, selectivities, and molecular weight of olefin polymerization catalysts are presented using PIO analysis and LFO calculation. Not only computational chemists but also experimental chemists can easily use this new system for catalyst design or molecular design from the point of view of orbital interaction.     

Modeling pi-pi interactions with the effective fragment potential method: the benzene dimer and substituents

This study compares the results of the general effective fragment potential (EFP2) method to the results of a previous combined coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] and symmetry-adapted perturbation theory (SAPT) study [Sinnokrot and Sherrill, J. Am. Chem. Soc., 2004, 126, 7690] on substituent effects in pi-pi interactions. EFP2 is found to accurately model the binding energies of the benzene-benzene, benzene-phenol, benzene-toluene, benzene-fluorobenzene, and benzene-benzonitrile dimers, as compared with high-level methods [Sinnokrot and Sherrill, J. Am. Chem. Soc., 2004, 126, 7690], but at a fraction of the computational cost of CCSD(T). In addition, an EFP-based Monte Carlo/simulated annealing study was undertaken to examine the potential energy surface of the substituted dimers.     

A combined ONIOM quantum chemical-molecular dynamics study of zinc-uracil bond breaking in yeast cytosine deaminase

A QM/MM method that combines ONIOM quantum chemistry and molecular dynamics is developed and applied to a step in the deamination of cytosine to uracil in yeast cytosine deaminase (yCD). A two-layer ONIOM calculation is used for the reaction complex, with an inner part treated at a high level for the chemical reaction (bond breaking) and a middle part treated at a lower level for relevant protein residues that are frozen in the quantum optimization. An outer layer (protein and solvent) is treated using MD. Configurations for the entire system are generated using MD and optimized with ONIOM. The method permits the use of high-level quantum calculations along with sufficient configurational sampling to approximate the potential of mean force for certain bond-breaking reactions. A previously proposed reaction mechanism for deamination (Sklenak, S.; Yao, L. S.; Cukier, R. I.; Yan, H. G. J. Am. Chem. Soc. 2004, 126, 14879) requires breaking the bond between a catalytic zinc and the O4 of uracil in order to permit product release. Using an ONIOM approach, direct bond cleavage was found to be energetically unfavorable. In the work presented here, the combined ONIOM MD method is used to show that the barrier for bond cleavage is small, approximately 3 kcal/mol, and, consequently, should not be the rate-limiting step in the reaction.     

Conserved mechanism of copper binding and transfer. A comparison of the copper-resistance proteins PcoC from Escherichia coli and CopC from Pseudomonas syringae

The copper-resistance proteins PcoC from Escherichia coli and CopC from Pseudomonas syringae exhibit 67% sequence identity, but the chemistry reported for PcoC (Peariso, K.; Huffman, D. L.; Penner-Hahn, J. E.; O'Halloran, T. V. J. Am. Chem. Soc. 2003, 125, 342-343) was distinctly different from that reported for CopC (Zhang, L.; Koay, M.; Maher, M. J.; Xiao, Z.; Wedd, A. G. J. Am. Chem. Soc. 2006, 128, 5834-5850). The source of the inconsistency has been identified, and His1 is confirmed as an unprecedented bidentate ligand in each protein. Access to a bona fide wild-type PcoC protein allowed unequivocal observation of intermediates involved in intermolecular redox copper transfer reactions.     

Structural studies of nonclassical cyclobutylmethyl cations by the ab initio method

Ab initio calculations at the MP2/cc-pVTZ level show that the cyclobutylmethyl cation is a nonclassical sigma-delocalized species, which is distinct from the global minimum C2-symmetric cyclopentyl cation (Schleyer, P. v. R.; Carneiro, J. W. de M.; Koch, W.; Raghavachari, K. J. Am. Chem. Soc. 1989, 111, 5475). Relatively lower level DFT calculations, on the other hand, show that the primary cyclobutylmethyl cation spontaneously collapses into the cyclopentyl cation (Prakash, G. K. S.; Reddy, V. P.; Rasul, G.; Casanova, J.; Olah, G. A. J. Am. Chem. Soc. 1998, 120, 13362). The secondary 1-cyclobutylethyl cation is also a nonclassical carbocation, as shown by calculations at the MP2/cc-pVTZ level. Two structures having energy minima are identified for the latter cation on the potential energy surface. The conformer in which the methyl group is in the exo orientation is a global minimum and is favored over the corresponding endo conformer by 1.2 kcal/mol at the MP2/cc-pVTZ//MP2/cc-pVTZ +ZPE level of calculations. The tertiary 1-cyclobutyl-1-methylethyl cation, at this level of calculations, also involves substantial nonclassical sigma-delocalization, showing that the nonclassical stabilization is more important for cyclobutylmethyl cations relative to the cyclopropylmethyl cations. The 13C NMR chemical shifts obtained from GIAO-CCSD(T)/tzp/dz calculations further substantiate the nonclassical structures for these carbocations.     

Do spin state changes matter in organometallic chemistry? A computational study

Spin changes occur often in organometallic chemistry, and their effect on kinetics is not well understood. We report computations on the singlet and triplet potential energy surfaces of several processes of this type and show that the topology of the individual surfaces, as well as of the crossing regions between them, can be used to rationalize the observed reactivity in all cases. In particular, the slow addition of dihydrogen to W[N(CH(2)CH(2)NSiMe(3))(3)]H (Schrock, R. R.; Shih, K. Y.; Dobbs, D. A.; Davis, W. M. J. Am. Chem. Soc. 1995, 117, 6609) is shown to be a "spin-blocked" reaction with a high barrier due to the crossing between reactant triplet and product singlet surfaces. In contrast, addition of CO to TpCo(CO) (Detrich, J. L.; Reinaud, O. M.; Rheingold, A. L.; Theopold, K. H. J. Am. Chem. Soc. 1995, 117, 11745) is fast because the triplet and singlet surfaces cross at low energy. Particular care is taken to use DFT methods which are in adequate agreement with experimental and high-level computational energetics for these systems.     

Two-state reactivity in alkane hydroxylation by non-heme iron-oxo complexes

Density functional theory is used to explore the mechanisms of alkane hydroxylation for four synthetic non-heme iron(IV)-oxo complexes with three target substrates (Kaizer, J.; Klinker, E. J.; Oh, N. Y.; Rohde; J.-U.; Song, W. J.; Stubna, A.; Kim, J.; Münck, E.; Nam, W.; Que, L., Jr. J. Am. Chem. Soc. 2004, 126, 472-473; Rohde, J.-U.; Que, L., Jr. Angew. Chem. Int. Ed. 2005, 44, 2255-2258.). The iron-oxo reagents possess triplet ground states and low-lying quintet excited states. The set of experimental and theoretical reactivity trends can be understood if the reactions proceed on the two spin states, namely two-state reactivity (TSR); an appropriate new model is presented. The TSR model makes testable predictions: (a) If crossing to the quintet state occurs, the hydroxylation will be effectively concerted; however, if the process transpires only on the triplet surface, stepwise hydroxylation will occur, and side products derived from radical intermediates would be observed (e.g., loss of stereochemistry). (b) In cases of crossing en route to the quintet transition state, one expects kinetic isotope effects (KIEs) typical of tunneling. (c) In situations where the two surfaces contribute to the rate, one expects intermediate KIEs and radical scrambling patterns that reflect the two processes. (d) Solvent effects on these reactions are expected to be very large.     

Characteristics of a miniature compartment-less glucose-O2 biofuel cell and its operation in a living plant

We report the temperature, pH, glucose concentration, NaCl concentration, and operating atmosphere dependence of the power output of a compartment-less miniature glucose-O(2) biofuel cell, comprised only of two bioelectrocatalyst-coated carbon fibers, each of 7 micro m diameter and 2 cm length (Mano, N.; Mao, F.; Heller, A. J. Am. Chem. Soc. 2002, 124, 12962). The bioelectrocatalyst of the anode consists of glucose oxidase from Aspergillus niger electrically "wired" by polymer I, having a redox potential of -0.19 V vs Ag/AgCl. That of the cathode consists of bilirubin oxidase from Trachyderma tsunodae "wired" by polymer II having a redox potential of +0.36 V vs Ag/AgCl (Mano, N.; Kim, H.-H.; Zhang, Y.; Heller, A. J. Am. Chem. Soc. 2002, 124, 6480. Mano, N.; Kim, H.-H.; Heller, A. J. Phys. Chem. B 2002, 106, 8842). Implantation of the fibers in the grape leads to an operating biofuel cell producing 2.4 micro W at 0.52 V.     

Dinuclear Zn(II) complex catalyzed phosphodiester cleavage proceeds via a concerted mechanism: a density functional theory study

Density functional theory (DFT) calculations were used to study the mechanism for the cleavage reaction of the RNA analogue HpPNP (HpPNP = 2-hydroxypropyl-4-nitrophenyl phosphate) catalyzed by the dinuclear Zn(II) complex of 1,3-bis(1,4,7-triazacyclonon-1-yl)-2-hydroxypropane (Zn(2)(L(2)O)). We present a binding mode in which each terminal phosphoryl oxygen atom binds to one zinc center, respectively, and the nucleophilic 2-hydroxypropyl group coordinates to one of the zinc ions, while the hydroxide from deprotonation of a water molecule coordinates to the other zinc ion. Our calculations found a concerted mechanism for the HpPNP cleavage with a 16.5 kcal/mol reaction barrier. An alternative proposed stepwise mechanism through a pentavalent oxyphosphorane dianion reaction intermediate for the HpPNP cleavage was found to be less feasible with a significantly higher energy barrier. In this stepwise mechanism, the deprotonation of the nucleophilic 2-hydroxypropyl group is accompanied with nucleophilic attack in the rate-determining step. Calculations of the nucleophile (18)O kinetic isotope effect (KIE) and leaving (18)O KIE for the concerted mechanism are in reasonably good agreement with the experimental values. Our results indicate a specific-base catalysis mechanism takes place in which the deprotonation of the nucleophilic 2-hydroxypropyl group occurs in a pre-equilibrium step followed by a nucleophilic attack on the phosphorus center. Detailed comparison of the geometric and electronic structure for the HpPNP cleavage reaction mechanisms in the presence/absence of catalyst revealed that the catalyst significantly altered the determining-step transition state to become far more associative or tight, that is, bond formation to the nucleophile was remarkably more advanced than leaving group bond fission in the catalyzed mechanism. Our results are consistent with and provide a reliable interpretation for the experimental observations that suggest the reaction occurs by a concerted mechanism (see Humphry, T.; Iyer, S.; Iranzo, O.; Morrow, J. R.; Richard, J. P.; Paneth, P.; Hengge, A. C. J. Am. Chem. Soc. 2008, 130, 17858-17866) and has a specific-base catalysis character (see Yang, M.-Y.; Iranzo, O.; Richard, J. P.; Morrow, J. R. J. Am. Chem. Soc. 2005, 127, 1064-1065).     

The phenylnitrene rearrangement in the inner phase of a hemicarcerand

The photochemistry of phenyl azide 1 and 13C-labeled phenyl azide 13C-1 incarcerated inside a hemicarcerand 4 was investigated. Low-temperature photolysis of hemicarceplex 41 and 413C-1 yields incarcerated 1-azacyclohepta-1,2,4,6-tetraene 42 and 413C-2 (18-50%), respectively, which were characterized by low-temperature FT-IR and 1H NMR and 13C NMR spectroscopy. After correction for the hemicarcerand-induced upfield shift, the 13C chemical shifts of incarcerated 13C-2 compare very well (Deltadelta 相似文献   

A coupled reference interaction site model/molecular dynamics study of the potential of mean force curve of the SN2 Cl- + CH3Cl reaction in water

An application of the coupled reference interaction site model (RISM)/simulation methodology to the calculation of the potential of mean force (PMF) curve in aqueous solution for the identity nucleophilic substitution reaction Cl(-) + CH(3)Cl is performed. The free energy of activation is calculated to be 27.1 kcal/mol which compares very well with the experimentally determined barrier height of 26.6 kcal/mol. Furthermore, the calculated PMF is almost superimposed with that previously calculated using the computationally rigorous Monte Carlo with importance sampling method (Chandrasekhar, J.; Smith, S. F.; Jorgensen, W. L. J. Am. Chem. Soc. 1985, 107, 154). Using the calculated PMF, a crude estimate of the solvated kinetic transmission coefficient also compares well with that of previous more accurate simulations. These results indicate that the coupled RISM/simulation method provides a cost-effective methodology for studying reactions in solution.     

Tension-induced morphological transition in mixed lipid bilayers

Recently, Rozovsky et al. reported on the morphology and dynamics of superstructures in three-component lipid bilayers containing saturated lipid, unsaturated lipid, and cholesterol (Rozovsky, S.; Kaizuka, Y.; Groves, J. T. J. Am. Chem. Soc. 2005, 127, 36). We suggest that the observed sequence of the striped-to-hexagonal morphological transition in mixed bilayers can be attributed to an enhanced membrane surface tension that is induced by the vesicle adhesion onto the solid surface.     

Chemistry of periodate-mediated cross-linking of 3,4-dihydroxylphenylalanine-containing molecules to proteins

Chemical cross-linking is an attractive approach to map peptide-protein and protein-protein complexes. Previously, we explored 3,4-dihydroxylphenylalanine (DOPA) as a protein cross-linking agent upon periodate oxidation (Burdine, L.; Gillette, T. G.; Lin, H.-J.; Kodadek, T. J. Am. Chem. Soc. 2004, 126, 11442-11443). We report here a study on the chemistry of DOPA-protein cross-linking. First, using a peptide nucleic acid templated system, we identified the alpha-amino, epsilon-amino of Lys, imidazole of His, and thiol of Cys as functional groups capable of attacking DOPA ortho-quinone. Second, we demonstrated that periodate-induced DOPA-protein cross-linking could be carried out efficiently at neutral pH in the presence of excess aliphatic 1,2-diols such as ethylene glycol, lactose, and adenosine triphosphate. This result indicated that DOPA-protein cross-linking and 1,2-diol oxidative cleavage proceed via different mechanisms and that carbohydrates will not interfere with this process when carried out in crude cell extracts or on intact cells.