Search

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.

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.

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.

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.

Calculations show that the transition structure for the synchronous oxygen transfer by Compound I is a second order saddle point. The process is unlikely.

Density functional calculations on oxo-manganese complexes of corrole (1) and porphyrin (2 and 3) show a fundamental difference. The ground state of I is the singlet manganese(V) state, 1A(MnV), in which corrole is a closed shell. In contrast, 2 and 3 have high-spin manganese(IV) states, 1A1u and 3A2u respectively. This difference and the state ordering for each system are rationalized based on the competition between the intrinsic tendency of manganese to prefer high-spin electronic configurations, vis-à-vis the general tendency to prefer double occupancy in the low-lying orbitals. The outcome of this competition is determined primarily by the identity of the macrocycle, corrole versus porphyrin. Corrole with a small cavity holds the MnO moiety with a high off-plane displacement, and thereby prefers the low-spin state. On the other hand, porphyrin with the wider cavity holds the MnO moiety closer to the plane, and thereby prefers high-spin states.

Bond stabilization through bending! Valence bond analysis shows that the σ frames of 1-3 (1: E = Si, E' = C; 2: E = E' = Si; 3: E = E' = C) are stabilized by trans bending (B), while π bonding is weakened. In acetylene (3) π bonding overrides the σ frame and establishes a linear molecule (3 L). In contrast, the σ frames dominate in silyne (1) and disilyne (2) and lead to trans-bent structures (1 B and 2 B)., (© 2001 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

Chameleon states: the ruthenium and iron metalloporphyrin analogues of compound I of cytochrome P450 (1; L = thiolate) possess low-lying states that change their electronic structure with solvent polarization. The ground state of the ruthenium complex is a low-spin electrophilic state, whereas the ground state of the iron complex is triradicaloid., (© 2001 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

The suicidal complex 4 2, which inactivates cytochrome P450 during olefin epoxidation, was shown by density functional calculations to be formed from the same high-spin intermediate ( 4 1-III) that leads to stereochemical scrambling., (© 2001 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

The epoxidation of ethene by a model for Compound I of cytochrome P450, studied by the use of density functional B3LYP calculations, involves two-state reactivity (TSR) with multiple electromer species, hence "multi-state epoxidation". The reaction is found to proceed in stepwise and effectively concerted manners. Several reactive states are involved; the reactant is an (oxo)iron(IV) porphyrin cation radical complex with two closely lying spin states (quartet and doublet), both of which react with ethene to form intermediate complexes with a covalent C-O bond and a carbon-centered radical (radical intermediates). The radical intermediates exist in two electromers that differ in the oxidation state of iron; Por(+)(*)Fe(III)OCH(2)CH(2)(*) and PorFe(IV)OCH(2)CH(2)(*) (Por = porphyrin). These radical intermediates exist in both the doublet- and quartet spin states. The quartet spin intermediates have substantial barriers for transformation to the quartet spin PorFe(III)-epoxide complex (2.3 kcal mol(-)(1) for PorFe(IV)OCH(2)CH(2)(*) and 7.2 kcal mol(-)(1) for Por(+)(*)Fe(III)OCH(2)CH(2)(*)). In contrast, the doublet spin radicals collapse to the corresponding PorFe(III)-epoxide complex with virtually no barriers. Consequently, the lifetimes of the radical intermediates are much longer on the quartet- than on the doublet spin surface. The loss of isomeric identity in the epoxide and rearrangements to other products arise therefore mostly, if not only, from the quartet process, while the doublet state epoxidation is effectively concerted (Scheme 7). Experimental trends are discussed in the light of the computed mechanistic scheme, and a comparison is made with closely related mechanistic schemes deduced from experiment.

Molecular dynamics has been used to investigate the reaction of a series of ketyl anion radicals and alkyl halides, CH2O(*)(-) + CH3X (X = F, Cl, Br) and NCCHO(*)(-) + CH3Cl. In addition to a floppy outer-sphere transition state which leads directly to ET products, there is a strongly bound transition state that yields both electron transfer (ET) and C-alkylated (SUB(C)) products. This common transition state has significant C-- C bonding and gives ET and SUB(C) products via a bifurcation on a single potential energy surface. Branching ratios have been estimated from ab initio classical trajectory calculations. The SUB(C) products are favored for transition states with short C--C bonds and ET for long C--C bonds. ET reactivity can be observed even at short distances of r(C)(-)(C) = ca. 2.4 A as in the transition state for the reaction NCCHO(*)(-) + CH3Cl. Therefore, the ET/SUB(C) reactivity is entangled over a significant range of the C--C distance. The mechanistic significance of the molecular dynamics study is discussed.

The flexibility of valence bond (VB) theory provides a new method of calculating pi-bond energies in the double-bonded species H(m)A=BH(n), where A, B = C, N, O, Si, P, S. This new method circumvents the problems usually associated with obtaining pi-bond strengths by targeting only the pi bond, while all other factors remain constant. In this manner, a clean separation between sigma- and pi effects can be achieved which highlights some expected trends in bond strength upon moving from left to right and up and down the Periodic Table. Intra-row pi bonds conform to the classic statement by Pauling [L. Pauling, The Natiure of the Chemical Bond, Cornell University Press, Ithaca, 1960, 3rd edition] regarding the relationship of heteronuclear bond strengths to their homonuclear constituents whereas inter-row pi bonds do not. This variance with Pauling's statement is shown to be due to the constraining effect of the underlying sigma bonds which prevents optimal p(pi)-p(pi) overlap. While Pauling's statement was based on the assumption that the resonance energy (RE) would be large for heteronuclear and small for homonuclear bonds, we have found large REs for all bonds studied herein; this leads to the conclusion that REs are dependent not only on the electronegativity difference but also the electronegativity sum of the constituent atoms. This situation where the bond is neither covalent nor ionic but originates in the covalent-ionic mixing has been termed charge shift (CS) bonding [S. Shaik, P. Maitre, G. Sini, P. C. Hiberty, J. Am. Chem. Soc. 1992, 114, 7861]. We have shown that CS bonding extends beyond single sigma bonds in first row molecules, thus supporting the idea that CS-bonding is a ubiquitous bonding form.

The Bergman cyclization of (Z)-hex-3-ene-1,5-diynes (1, enediynes), which produces pharmacologically important DNA-cleaving biradicals (1,4-benzyne, 2), was studied by using Hartree-Fock (HF) and density-functional theory (DFT) based valence bond (VB) methods (VB-HF and VB-DFT, respectively). We found that only three VB configurations are needed to arrive at results not too far from complete active space [CASSCF(6 x 6)] computations, while the quality of VB-DTF utilizing the same three configurations improves upon CASSCF(6 x 6) analogous to CASPT2. The dominant VB configuration in 1 contributes little to 2, while the most important biradical configuration in 2 plays a negligible role in 1. The avoided crossing of the energy curves of these two configurations along the reaction coordinate leads to the transition state (TS). As a consequence of the shape and position of the crossing section, the changes in geometry and in the electronic wavefunction along the reaction coordinate are non-synchronous; the TS is geometrically approximately 80% product-like and electronically approximately 70% reactant-like. While the pi resonance in the TS is very small, it is large (64.4 kcal mol(-1)) for 2 (cf. benzene=61.5 kcal mol(-1)). As a consequence, substituents operating on the sigma electrons should be much more effective in changing the Bergman reaction cyclization barrier. Furthermore, additional sigma resonance in 2 results in unusually high values for the nucleus-independent chemical shift (NICS, a direct measure for aromaticity). Similarly, the high NICS value of the TS is due mostly to sigma resonance to which the NICS procedure is relatively sensitive.

It is proposed that spin-crossing effects can dramatically affect reaction mechanisms, rate constants, branching ratios, and temperature behaviors of organometallic transformations. This phenomenon is termed two-state reactivity (TSR) and involves participation of spin inversion in the rate-determining step. While the present analysis is based on studies of transition metals under idealized conditions, several recent reports imply that TSR is by no means confined to the gas phase. In fact, participation of more than a single spin surface in the reaction pathways is proposed as a key feature in organometallic chemistry.

Two electromeric forms, a and b (a is the ground state in a solvent) exist for the hydroxo-iron complex 1, an intermediate in the rebound mechanism of alkane hydroxylation by cytochrome P450. Results of density functional and model solvent calculations of various species are in agreement with experimental findings, and imply the role of 1 a in the rebound mechanism.

A unified description of chemical reactivity is made possible with valence bond (VB) diagrams, such as those shown for the a simple case (only reactants and products must be considered) and a more complex system (an intermediate state also plays an important role; RC = reaction coordinate). Analysis of reactivity and mechanistic problems in organic and organometallic chemistry exemplifies the generality of the VB paradigm: in situ DNA repair, C-F and C-H bond activation, S RN 2 c mechanism, stepwise versus concerted cycloaddition, and a lot more., (© 1999 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

A twin of the transition state, which can be investigated spectroscopically and can thus supply information about the structure of the transition state, has now been characterized for the Cope rearrangement of semibullvalene (shown below). It involves an excited state with B 2 symmetry and results from a linear combination of the ground-state wave functions of (mirror-image) reactant and product., (© 1998 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

A prospective study was conducted on paediatric thermal injury patients admitted to the Burns Unit at King Fahd Hospital, Jeddah, Saudi Arabia over a 2-year period (December 1985 to December 1987). A total of 197 patients (out of 319) were paediatric, aged up to 18 years. Infants and toddlers accounted for 59.4 per cent and adolescents for 14.2 per cent. Scalding and flame injuries accounted for 98 per cent with most injuries occurring at home (97.5 per cent) and the overall paediatric mortality rate was 4.4 per cent.

Three hundred and nineteen patients with different types of burns were studied at King Fahd Hospital, Jeddah, Saudi Arabia over a 2-year period (December, 1985 to December, 1987). Scalding was the most common cause (56.4 per cent) of burn injuries compared with 41.4 per cent of patients who sustained flame injury; 84.6 per cent of the thermal injuries occurred at home, with children (less than or equal to 18 years of age) being affected most frequently (61.8 per cent). The overall mortality was 9.4 per cent.