Biomedical subjects

Sason Shaik

Publications and source records attributed to Sason Shaik.

71 records · Page 4Linked to original sources

PubMed2003

An accurate barrier for the hydrogen exchange reaction from valence bond theory: is this theory coming of age?

One of the landmark achievements of quantum chemistry, specifically of MO-based methods that include electron correlation, was the precise calculation of the barrier for the hydrogen-exchange reaction (B. Liu, J. Chem. Phys. 1973, 58, 1925; P. Siegbahn, B. Liu, J. Chem. Phys. 1978, 68, 2457). This paper reports an accurate calculation of this barrier by two recently developed VB methods that use only the eight classical VB structures. To our knowledge, the present work is the first accurate ab initio VB barrier that matches an experimental value. Along with the accurate barrier, the VB method provides accurate bond energies and diabatic quantities that enable the barrier height to be analyzed by the VB state correlation diagram approach, VBSCD (S. Shaik, A. Shurki, Angew. Chem. 1999, 111, 616; Angew. Chem. Int. Ed. Engl. 1999, 38, 586). This is a proof of principal that VB theory with appropriate account of dynamic electron correlation can achieve quantitative accuracy of reaction barriers, and still retain a compact and interpretable wave function. A sample of S(N)2 barriers and dihalogen bonding energies, which are close to CCSD(T) and G2(+) values, show that the H(3) problem is not an isolated case, and while it is premature to conclude that VB theory has come of age, the occurrence of this event is clearly within sight.

Journal Article↗

PubMed2003

Can a single oxidant with two spin states masquerade as two different oxidants? A study of the sulfoxidation mechanism by cytochrome p450.

Density functional calculations were performed on the sulfoxidation reaction by a model compound I (Cpd I) of cytochrome P450. By contrast to previous alkane hydroxylation studies, which exhibit a dominant low-spin (LS) pathway, the sulfoxidation follows a dominant high-spin (HS) reaction. Thus, competing hydroxylation and sulfoxidation processes as observed for instance by Jones et al. (Volz, T. J.; Rock, D. A.; Jones, J. P. J. Am. Chem. Soc. 2002, 124, 9724) are the result of a two-state reactivity scenario, whereby the hydroxylation originates from the LS pathway and the sulfoxidation from the HS pathway. In this manner, two spin states of a single oxidant (Cpd I) can be disguised as two different oxidants. The calculations rule out the possibility that a second oxidant (the ferric peroxide, Cpd 0 species) interferes in the observed results of Jones et al.

Cytochrome P-450 Enzyme System↗

PubMed2003

Chemistry--a central pillar of human culture.

Explore the source record for details and available documents.

Chemistry↗

PubMed2003

A proton-shuttle mechanism mediated by the porphyrin in benzene hydroxylation by cytochrome p450 enzymes.

Benzene hydroxylation is a fundamental process in chemical catalysis. In nature, this reaction is catalyzed by the enzyme cytochrome P450 via oxygen transfer in a still debated mechanism of considerable complexity. The paper uses hybrid density functional calculations to elucidate the mechanisms by which benzene is converted to phenol, benzene oxide, and ketone, by the active species of the enzyme, the high-valent iron-oxo porphyrin species. The effects of the protein polarity and hydrogen-bonding donation to the active species are mimicked, as before (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893). It is verified that the reaction does not proceed either by hydrogen abstraction or by initial electron transfer (Ortiz de Montellano, P. R. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 8, pp 245-303). In accord with the latest experimental conclusions, the theoretical calculations show that the reactivity is an interplay of electrophilic and radicalar pathways, which involve an initial attack on the pi-system of the benzene to produce sigma-complexes (Korzekwa, K. R.; Swinney, D. C.; Trager, W. T. Biochemistry 1989, 28, 9019-9027). The dominant reaction channel is electrophilic and proceeds via the cationic sigma-complex,( 2)3, that involves an internal ion pair made from a cationic benzene moiety and an anionic iron porphyrin. The minor channel proceeds by intermediacy of the radical sigma-complex, (2)2, in which the benzene moiety is radicalar and the iron-porphyrin moiety is neutral. Ring closure in these intermediates produces the benzene oxide product ((2)4), which does not rearrange to phenol ((2)7) or cyclohexenone ((2)6). While such a rearrangement can occur post-enzymatically under physiological conditions by acid catalysis, the computations reveal a novel mechanism whereby the active species of the enzyme catalyzes directly the production of phenol and cyclohexenone. This enzymatic mechanism involves proton shuttles mediated by the porphyrin ring through the N-protonated intermediate, (2)5, which relays the proton either to the oxygen atom to form phenol ((2)7) or to the ortho-carbon atom to produce cyclohexenone product ((2)6). The formation of the phenol via this proton-shuttle mechanism will be competitive with the nonenzymatic conversion of benzene oxide to phenol by external acid catalysis. With the assumption that (2)5 is not fully thermalized, this novel mechanism would account also for the observation that there is a partial skeletal retention of the original hydrogen of the activated C-H bond, due to migration of the hydrogen from the site of hydroxylation to the adjacent carbon (so-called "NIH shift" (Jerina, D. M.; Daly, J. W. Science 1974, 185, 573-582)). Thus, in general, the computationally discovered mechanism of a porphyrin proton shuttle suggests thatthere is an enzymatic pathway that converts benzene directly to a phenol and ketone, in addition to nonenzymatic production of these species by conversion of arene oxide to phenol and ketone. The potential generality of protonated porphyrin intermediates in P450 chemistry is discussed in the light of the H/D exchange observed during some olefin epoxidation reactions (Groves, J. T.; Avaria-Neisser, G. E.; Fish, K. M.; Imachi, M.; Kuczkowski, R. J. Am. Chem. Soc. 1986, 108, 3837-3838) and the general observation of heme alkylation products (Kunze, K. L.; Mangold, B. L. K.; Wheeler, C.; Beilan, H. S.; Ortiz de Montellano, P. R. J. Biol. Chem. 1983, 258, 4202-4207). The competition, similarities, and differences between benzene oxidation viz. olefin epoxidation and alkanyl C-H hydroxylation are discussed, and comparison is made with relevant experimental and computational data. The dominance of low-spin reactivity in benzene hydroxylation viz. two-state reactivity (Shaik, S.; de Visser, S. P.; Ogliaro, F.; Schwarz, H.; Schröder, D. Curr. Opin. Chem. Biol. 2002, 6, 556-567) in olefin epoxidation and alkane hydroxylation is traced to the loss of benzene resonance energy during the bond activation step.

Benzene↗

PubMed2003

Fluorinated alcohols enable olefin epoxidation by H2O2: template catalysis.

Experimental observations show that direct olefin epoxidation by H(2)O(2), which is extremely sluggish otherwise, occurs in fluorinated alcohol (R(f)OH) solutions under mild conditions requiring no additional catalysts. Theoretical calculations of ethene and propene epoxidation by H(2)O(2) in the gas phase and in the presence of methanol and of two fluorinated alcohols, presented in this paper, show that the fluoro alcohol itself acts as a catalyst for the reaction by providing a template that stabilizes specifically the transition state (TS) of the reaction. Thus, much like an enzyme, the fluoro alcohol provides a complementary charge template that leads to the reduction of the barrier by 5-8 kcal mol(-)(1). Additionally, the fluoro alcohol template keeps the departing OH and hydroxyalkenyl moieties in close proximity and, by polarizing them, facilitates the hydrogen migration from the latter to form water and the epoxide product. The reduced activation energy and structural confinement of the TS over the fluoro alcohol template render the epoxidation reaction observable under mild synthetic conditions.

Journal Article↗

PubMed2003

Is the ruthenium analogue of compound I of cytochrome p450 an efficient oxidant? A theoretical investigation of the methane hydroxylation reaction.

High-valent metal-oxo complexes catalyze C-H bond activation by oxygen insertion, with an efficiency that depends on the identity of the transition metal and its oxidation state. Our study uses density functional calculations and theoretical analysis to derive fundamental factors of catalytic activity, by comparison of a ruthenium-oxo catalyst with its iron-oxo analogue toward methane hydroxylation. The study focuses on the ruthenium analogue of the active species of the enzyme cytochrome P450, which is known to be among the most potent catalysts for C-H activation. The computed reaction pathways reveal one high-spin (HS) and two low-spin (LS) mechanisms, all nascent from the low-lying states of the ruthenium-oxo catalyst (Ogliaro, F.; de Visser, S. P.; Groves, J. T.; Shaik, S. Angew. Chem. Int. Ed. 2001, 40, 2874-2878). These mechanisms involve a bond activation phase, in which the transition states (TS's) appear as hydrogen abstraction species, followed by a C-O bond making phase, through a rebound of the methyl radical on the metal-hydroxo complex. However, while the HS mechanism has a significant rebound barrier, and hence a long lifetime of the radical intermediate, by contrast, the LS ones are effectively concerted with small barriers to rebound, if at all. Unlike the iron catalyst, the hydroxylation reaction for the ruthenium analogue is expected to follow largely a single-state reactivity on the LS surface, due to a very large rebound barrier of the HS process and to the more efficient spin crossover expected for ruthenium. As such, ruthenium-oxo catalysts (Groves, J. T.; Shalyaev, K.; Lee, J. In The Porphyrin Handbook; Biochemistry and Binding: Activation of Small Molecules, Vol. 4; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Academic Press: New York, 2000; pp 17-40) are expected to lead to more stereoselective hydroxylations compared with the corresponding iron-oxo reactions. It is reasoned that the ruthenium-oxo catalyst should have larger turnover numbers compared with the iron-oxo analogue, due to lesser production of suicidal side products that destroy the catalyst (Ortiz de Montellano, P. R.; Beilan, H. S.; Kunze, K. L.; Mico, B. A. J. Biol. Chem. 1981, 256, 4395-4399). The computations reveal also that the ruthenium complex is more electrophilic than its iron analogue, having lower hydrogen abstraction barriers. These reactivity features of the ruthenium-oxo system are analyzed and shown to originate from a key fundamental factor, namely, the strong 4d(Ru)-2p(O,N) overlaps, which produce high-lying pi(Ru-O), sigma(Ru-O), and sigma(Ru-N) orbitals and thereby to lead to a preference of ruthenium for higher-valent oxidation states with higher electrophilicity, for the effectively concerted LS hydroxylation mechanism, and for less suicidal complexes. As such, the ruthenium-oxo species is predicted to be a more robust catalyst than its iron-oxo analogue.

Biomimetic Materials↗

PubMed2003

A conversation on VB vs MO theory: a never-ending rivalry?

Quantum mechanics has provided chemistry with two general theories, valence bond (VB) theory and molecular orbital (MO) theory. The two theories were developed at about the same time, but quickly diverged into rival schools that have competed, sometimes fervently, on charting the mental map and epistemology of chemistry. Three practitioners of MO and VB theory talk - fighting a little, trying to understand - of the past and present of these two approaches to describing bonding in molecules.

Journal Article↗

PubMed2002

What factors affect the regioselectivity of oxidation by cytochrome p450? A DFT study of allylic hydroxylation and double bond epoxidation in a model reaction.

Epoxidation (C=C) vis-à-vis allylic hydroxylation (C-H) reactions of propene with a model compound I (Cpd I) of the enzyme cytochrome P450 were studied using B3LYP density functional theory. Potential energy profiles and kinetic isotope effects (KIE) were calculated. The interactions in the protein pocket were mimicked by adding two external NH- - -S hydrogen bonds to the thiolate ligand and by introducing a nonpolar medium (with a dielectric constant, epsilon = 5.7) that can exert a polarization effect on the reacting species. A two-state reactivity (TSR) with high-spin (HS) and low-spin (LS) states were located for both processes (Ogliaro, F.; Harris, N.; Cohen, S.; Filatov, M.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 8977-8989. de Visser, S. P.; Ogilaro, F.; Harris, N.; Shaik, S. J. Am. Chem. Soc. 2001, 123, 3037-3047). The HS processes were found to be stepwise, whereas the LS processes were characterized as nonsynchronous but effectively concerted pathways. The computed KIE for C-H hydroxylation with and without tunneling corrections are large (>7), and they support the assignment of the corresponding transition states as hydrogen-abstraction species (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structures, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 1; pp 3-48). In the gas phase, epoxidation is energetically favorable by 3.4 kcal mol(-1). Inclusion of zero-point energy reduces this difference but still predicts C=C/C-H > 1. Environmental effects were found to have major impact on the C=C/C-H ratio as well as on the stereoselectivity of the processes. Thus, two NH- - -S hydrogen bonds away from the reaction center reverse the regioselectivity and prefer hydroxylation, namely, C=C/C-H <1. The polarity of the medium further accentuates the trend and leads to a change by 2 orders of magnitude in the regioselectivity, C=C/C-H << 1. Furthermore, since the environmental interactions prefer the LS over the HS reactions, both hydroxylation and epoxidation processes are rendered more stereoselective, again by 2 orders of magnitude. It follows, therefore, that Cpd I is a chameleon oxidant (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893; Ogliaro, F.; de Visser, S. P.; Cohen, S.; Kaneti, J.; Shaik, S. Chembiochem. 2001, 2, 848-851; Ogliaro, F.; de Visser, S. P.; Groves, J. T.; Shaik, S. Angew. Chem., Int. Ed. 2001, 40, 2874-2878) that tunes its reactivity and selectivity patterns in response to the protein environment in which it is accommodated. A valence bond (VB) model, akin to "redox mesomerism" (Bernadou, J.; Fabiano, A.-S.; Robert, A.; Meunier, B. J. Am. Chem. Soc. 1994, 116, 9375-9376), is constructed and enables the description of a chameleon transition state. It shows that the good donor ability of the thiolate ligand and the acceptor ability of the iron porphyrin create mixed-valent situations that endow the transition state with a great sensitivity to external perturbations as in the protein pocket. The model is used to discuss the computed results and to relate them to experimental findings.

Alkenes↗

PubMed2002

The 'push' effect of the thiolate ligand in cytochrome P450: a theoretical gauging.

The 'push' effect of the thiolate ligand in cytochrome P450 is investigated using density functional calculations. Theory supports Dawson's postulate that the 'push' effect is crucial for the heterolytic O-O bond cleavage of ferric-peroxide, as well as for controlling the Fe(III)/Fe(II) redox process and gating the catalytic cycle. Two energetic factors that contribute to the 'push' effect are revealed. The dominant one is the field factor (DeltaE(field)=54-103 kcal/mol) that accounts for the classical electrostatic repulsion with the negative charge of thiolate. The smaller factor is a quantum mechanical effect (DeltaE(QM)(sigma)=39 kcal/mol, DeltaE(QM)(pi)=4 kcal/mol), which is associated with the sigma- and pi-donor capabilities of thiolate. The effects of ligand replacement, changes in hydrogen bonding and dielectric screening are discussed in term of these quantities. In an environment with a dielectric constant of 5.7, the total 'push' effect is reduced to 29-33 kcal/mol. Manifestations of the 'push' effect on other properties of thiolate enzymes are discussed.

Cytochrome P-450 Enzyme System↗

PubMed2002

The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations.

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.

Cytochrome P-450 Enzyme System↗

PubMed2002

Searching for the second oxidant in the catalytic cycle of cytochrome P450: a theoretical investigation of the iron(III)-hydroperoxo species and its epoxidation pathways.

Iron(III)-hydroperoxo, [Por(CysS)Fe(III)-OOH](-), a key species in the catalytic cycle of cytochrome P450, was recently identified by EPR/ENDOR spectroscopies (Davydov, R.; Makris, T. M.; Kofman, V.; Werst, D. E.; Sligar, S. G.; Hoffman, B. M. J. Am. Chem. Soc. 2001, 123, 1403-1415). It constitutes the last station of the preparative steps of the enzyme before oxidation of an organic compound and is implicated as the second oxidant capable of olefin epoxidation (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555-3560), in addition to the penultimate active species, Compound I (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; pp 3-48). In response, we present a density functional study of a model species and its ethylene epoxidation pathways. The study characterizes a variety of properties of iron(III)-hydroperoxo, such as the O-O bonding, the Fe-S bonding, Fe-O and Fe-S stretching frequencies, its electron attachment, and ionization energies. Wherever possible these properties are compared with those of Compound I. The proton affinities for protonation on the proximal and distal oxygen atoms of iron(III)-hydroperoxo, and the effect of the thiolate ligand thereof, are determined. In accordance with previous results (Harris, D. L.; Loew, G. H. J. Am. Chem. Soc. 1998, 120, 8941-8948), iron(III)-hydroperoxo is a strong base (as compared with water), and its distal protonation leads to a barrier-free formation of Compound I. The origins of this barrier-free process are discussed using a valence bond approach. It is shown that the presence of the thiolate is essential for this process, in line with the "push effect" deduced by experimentalists (Sono, M.; Roach, M. P.; Coulter, E. D.; Dawson, J. H. Chem. Rev. 1996, 96, 2841-2887). Finally, four epoxidation pathways of iron(III)-hydroxperoxo are located, in which the species transfers oxygen to ethylene either from the proximal or from the distal sites, in both concerted and stepwise manners. The barriers for the four mechanisms are 37-53 kcal/mol, in comparison with 14 kcal/mol for epoxidation by Compound I. It is therefore concluded that iron(III)-hydroperoxo, as such, cannot be a second oxidant, in line with its significant basicity and poor electron-accepting capability. Possible versions of a second oxidant are discussed.

Catalysis↗

PubMed2002

Two-state reactivity mechanisms of hydroxylation and epoxidation by cytochrome P-450 revealed by theory.

Recent computational studies of alkane hydroxylation and alkene epoxidation by a model active species of the enzyme cytochrome P-450 reveal a two-state reactivity (TSR) scenario in which the information content of the product distribution is determined jointly by two states. TSR is used to reconcile the dilemma of the consensus 'rebound mechanism' of alkane hydroxylation, which emerged from experimental studies of ultra-fast radical clocks. The dilemma, stated succinctly as 'radicals are both present and absent and the rebound mechanism is both right and wrong', is simply understood once one is cognizant that the mechanism operates by two states, one low-spin (LS) the other high-spin (HS). In both states, bond activation proceeds in a manner akin to the rebound mechanism, but the LS mechanism is effectively concerted, whereas the HS is stepwise with incursion of radical intermediates.

Alkanes↗

PubMed2001