How does product isotope effect prove the operation of a two-state "rebound" mechanism in C-H hydroxylation by cytochrome P450?

C-H hydroxylation is a fundamental process. In Nature it is catalyzed by the enzyme cytochrome P450, in a still-debated mechanism that poses a major intellectual challenge for both experiment and theory; currently, the opinions keep swaying between the original single-state rebound mechanism, a two-oxidant mechanism (where ferric peroxide participates as a second oxidant, in addition to the primary active species, the high-valent iron-oxo species), and two-state reactivity (TSR) mechanism (where two spin states are involved). Recent product isotope effect (PIE) measurements for the trans-2-phenyl-methyl cyclopropane probe (1), led Newcomb and co-workers (Newcomb, M.; Aebisher, D.; Shen, R.; Esala, R.; Chandrasena, P.; Hollenberg, P.; Coon, M. J. J. Am. Chem. Soc. 2003, 125, 6064-6065) to rule out TSR in favor of the two-oxidant scenario, since the direction of the PIE was at odds with the one predicted from calculations on methane hydroxylation. The present report describes a density functional theoretical study of C-H hydroxylation of the Newcomb probe, 1, leading to rearranged (3) and unrearranged (2) products. Our study shows that the reaction occurs via TSR in which the high-spin pathway gives dominant rearranged products, whereas the low-spin pathway favors unrearranged products. The calculated PIE(2/3) values based on TSR are found to be in excellent agreement with the experimental data of Newcomb and co-workers. This match between experiment and theory makes a strong case that the reaction occurs via TSR mechanism.     

Thermal equilibrium of high- and low-spin forms of cytochrome P450 BM-3: repositioning of the substrate?

We demonstrate that cytochrome P450 BM-3 in complex with N-palmitoylglycine undergoes a spin state change between room temperature, where optimal activity is seen, and low temperatures, where X-ray diffraction characterization has been carried out. On the basis of NMR measurements of the full-length protein, this spin state change is likely to be accompanied by a general structural rearrangement in the enzyme pocket. The substrate remains bound at all temperatures. We propose that the substrate may "slide" from a position directly atop the heme (thus displacing the ligating water) to the more distant position (thus restoring the ligating water) as the temperature is lowered. This proposal is evaluated on the basis of computational modeling of the protein-ligand complex, using a novel induced fit methodology. We thereby generate a structure with the ligand in close contact with the heme, similar in energy to the experimental structure. With this combination of theory and experiment we provide a specific proposal of how ligands may be positioned for chemistry for this enzyme.     

Radical clock substrates, their C-H hydroxylation mechanism by cytochrome P450, and other reactivity patterns: what does theory reveal about the clocks' behavior?

There is an ongoing and tantalizing controversy regarding the mechanism of a key process in nature, C-H hydroxylation, by the enzyme cytochrome P450 (Auclaire, K.; Hu, Z.; Little, D. M.; Ortiz de Montellano, P. R.; Groves, J. T. J. Am. Chem. Soc. 2002, 124, 6020-6027. Newcomb, M.; Aebisher, D.; Shen, R.; Esala, R.; Chandrasena, P.; Hollenberg, P. F.; Coon, M. J. J. Am. Chem. Soc. 2003, 125, 6064-6065). To definitely resolve this controversy, theory must first address the actual systems that have been used by experiment, and that generated the controversy. This is done in the present paper, which constitutes the first extensive theoretical study of such two experimental systems, trans-2-phenylmethyl-cyclopropane (1) and trans-2-phenyl-iso-propylcyclopropane (4). The theoretical study of these substrates reveals that the only low energy pathway for C-H hydroxylation is the two-state rebound mechanism described originally for methane hydroxylation (Ogliaro, F.; Harris, N.; Cohen, S.; Filatov, M.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 8977-8989). The paper shows that the scenario of a two-state rebound mechanism accommodates much of the experimental data. The computational results provide a good match to experimental results concerning the very different extents of rearrangement for 1 (20-30%) vs 4 (virtually none), lead to product isotope effect for the reaction of 1, in the direction of the experimental result, and predict as well the observed metabolic switching from methyl to phenyl hydroxylation, which occurs upon deuteration of the methyl group. Furthermore, the study reveals that an intimate ion pair species involving an alkyl carbocation derived from 4 gives no rearranged products, again in accord with experiment. This coherent match between theory and experiment cannot be merely accidental; it comes close to being aproof that the actual mechanism of C-H hydroxylation involves the two-state reactivity revealed by theory. Analysis of the rearrangement modes of the carbocations derived from 1 and 4 excludes the participation of free carbocations during the hydroxylation of these substrates. Finally, the mechanistic significance of product isotope effect (different isotope effects for the rearranged and unrearranged alcohol products) is analyzed. It is shown to be a sensitive probe of two-state reactivity; the size of the intrinsic product isotope effect and its direction reveal the structural differences of the hydrogen abstraction transition states in the low-spin vs high-spin reaction manifolds.     

Does compound I vary significantly between isoforms of cytochrome P450?

The cytochrome P450 (CYP) enzymes are important in many areas, including pharmaceutical development. Subtle changes in the electronic structure of the active species, Compound I, have been postulated previously to account partly for the experimentally observed differences in reactivity between isoforms. Current predictive models of CYP metabolism typically assume an identical Compound I in all isoforms. Here we present a method to calculate the electronic structure and to estimate the Fe-O bond enthalpy of Compound I, and apply it to several human and bacterial CYP isoforms. Conformational flexibility is accounted for by sampling large numbers of structures from molecular dynamics simulations, which are subsequently optimized with density functional theory (B3LYP) based quantum mechanics/molecular mechanics. The observed differences in Compound I between human isoforms are small: They are generally smaller than the spread of values obtained for the same isoform starting from different initial structures. Hence, it is unlikely that the variation in activity between human isoforms is due to differences in the electronic structure of Compound I. A larger difference in electronic structure is observed between the human isoforms and P450(cam) and may be explained by the slightly different hydrogen-bonding environment surrounding the cysteinyl sulfur. The presence of substrate in the active site of all isoforms studied appears to cause a slight decrease in the Fe-O bond enthalpy, apparently due to displacement of water out of the active site, suggesting that Compound I is less stable in the presence of substrate.     

New model for a theoretical density functional theory investigation of the mechanism of the carbonic anhydrase: how does the internal bicarbonate rearrangement occur?

A theoretical density functional theory (DFT, B3LYP) investigation has been carried out on the catalytic cycle of the carbonic anhydrase. A model system including the Glu106 and Thr199 residues and the "deep" water molecule has been used. It has been found that the nucleophilic attack of the zinc-bound OH on the CO(2) molecule has a negligible barrier (only 1.2 kcal mol(-1)). This small value is due to a hydrogen-bond network involving Glu106, Thr199, and the deep water molecule. The two usually proposed mechanisms for the internal bicarbonate rearrangement have been carefully examined. In the presence of the two Glu106 and Thr199 residues, the direct proton transfer (Lipscomb mechanism) is a two-step process, which proceeds via a proton relay network characterized by two activation barriers of 4.4 and 9.0 kcal mol(-1). This pathway can effectively compete with a rotational mechanism (Lindskog mechanism), which has a barrier of 13.2 kcal mol(-1). The fast proton transfer found here is basically due to the effect of the Glu106 residue, which stabilizes an intermediate situation where the Glu106 fragment is protonated. In the absence of Glu106, the barrier for the proton transfer is much larger (32.3 kcal mol(-1)) and the Lindskog mechanism becomes favored.     

Oxygen economy of cytochrome P450: what is the origin of the mixed functionality as a dehydrogenase-oxidase enzyme compared with its normal function?

The economy of dioxygen consumption by enzymes constitutes a fundamental problem in enzymatic chemistry (ref 1). Sometimes, the enzyme converts ALL the oxygen into water, without affecting the organic substrate, thereby acting as an "oxidase" (ref 1). Other times, the enzyme converts all the oxygen into water and causes desaturation in the substrate, thus exhibiting a mixed function as both "oxidase" and "dehydrogenase" (refs 2-5). The present paper describes density functional calculations demonstrating that the oxidase-dehydrogenase mixed activity occurs from the cationic intermediate species and requires electro-steric inhibition of the rebound process. Furthermore, the calculations reveal that the carbocation is formally nascent from an excited state of the active species of the enzyme (2Cpd I), in which the Fe=O moiety is singlet coupled as in the 1Deltag state of dioxygen! Thus, our results resolve an important mechanism and reveal the factors that underlie its observability.     

The migration of styrene butadiene latex during the drying of coating suspensions: when and how does migration of colloidal particles occur?

Surface elemental compositions of model latex clay coatings on an impervious substrate consolidated under various conditions were measured using the XPS technique, in order to clarify when and how colloidal latex particles migrate to the surface during drying. Under similar drying conditions, surface carbon content decreased with the addition of a water-soluble polymer to the coating colors, while remaining virtually unchanged for coatings of different coat weights made with a given color, indicating that surface carbon content variation is mainly caused by migration of latex rather than of water-soluble polymer. The results also showed that for coatings made with a given suspension, surface carbon content decreased with increasing delay time between coating and heating. For coatings frozen during consolidation and dried by sublimation, surface carbon content increased with increasing drying time before freezing. These results suggest that for the model coatings studied, latex migration mainly occurs after coating application before capillary formation during the initial drying stage when coatings are in the liquid phase, contradicting both the conventional capillary transport and boundary wall migration mechanisms. An alternative mechanism which attributes latex migration to surface trapping effect and to higher Brownian mobility of the smaller latex particles compared with pigment appears to provide a systematically consistent explanation to those phenomena. The new particle migration mechanism implies that segregation of colloidal particles is a ubiquitous phenomenon that would occur not only during the drying of paper coatings but also during consolidation of colloidal films containing particles of different sizes. This is of great importance in the control of surface compositions of nanocomposite coatings.     

CCQM primary methods symposium: how far does the light shine?

The Consultative Committee for Amount of Substance (CCQM) has the task of organizing a comprehensive set of international comparisons to establish the technical basis for the mutual recognition of measurement capabilities among the national measurement institutes (NMIs) in the field of chemical measurement. The challenge that the CCQM faces is to identify, design, and conduct a limited number of key comparisons to enable the assessment of measurement comparability among NMIs across the entirety of ’chemical measurement space’. This is no easy task because the field of chemical metrology is extremely diverse and multidimensional, owing to the number of measurand types, concentrations, and matrix types of importance. The ”CCQM primary methods symposium: how far does the light shine?” was organized to provide information and initiate discussions to assist in this challenge, and clarify how the concept of a primary method of measurement could be instrumental in achieving this goal.     

Metal bioavailability: how does its significance change in the time?

The significance of terms "metal bioavailability" and "bioavailable metal fraction" is evolved in the time, passing from a very simple concept to a complex concept bound to abiotic and biotic aspects. At the beginning metal toxicity was related to metal fraction present in water phase, than only free metal ion activity was considered and the free ion activity model (FIAM) was proposed. Successively, due to the exceptions observed and to the consciousness that metal bioavailability could be considered as dynamic characteristic the concept of metal bioavailability became very complex, depending on physical, chemical and biological factors.     

How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.     

What factors influence the ratio of C-H hydroxylation versus C=C epoxidation by a nonheme cytochrome P450 biomimetic?

Density functional calculations on a nonheme biomimetic (Fe=O(TMCS+) have been performed and its catalytic properties versus propene investigated. Our studies show that this catalyst is able to chemoselectively hydroxylate C=H bonds even in the presence of C=C double bonds. This phenomenon has been analyzed and found to occur due to Pauli repusions between protons on the TMCS ligand with protons attached to the approaching substrate. The geometries of the rate determining transition states indicate that the steric hindrance is larger in the epoxidation transition states than in the hydroxylation ones with much shorter distances; hence the hydroxylation pathway is favored over the epoxidation. Although, the reactant experiences close lying triplet and quintet spin states, the dominant reaction mechanism takes place on the quintet spin state surface; i.e., Fe=O(TMCS)+ reacts via single-state reactivity. Our calculations show that this spin state selectivity is the result of geometric orientation of the transition state structures, whereby the triplet ones are destabilized by electrostatic repulsions between the substrate and the ligand while the quintet spin transition states are aligned along the ideal axis. The reactivity patterns and geometries are compared with oxoiron species of dioxygenase and monoxygenase enzymes. Thus, Fe=O(TMCS)+ shows some similarities with P450 enzyme reactivity: it chemoselectively hydroxylates C=H bonds even in the presence of a C=C double bond and therefore is an acceptable P450 biomimetic. However, the absolute barriers of substrate oxidation by Fe=O(TMCS)+ are higher than the ones obtained with heme enzymes, but the chemoselectivity is lesser affected by external perturbations such as hydrogen bonding of a methanol molecule toward the thiolate sulfur or a dielectric constant. This is the first oxoiron complex whereby we calculated a chemoselective hydroxylation over epoxidation in the gas phase.     

Quality control of Photosystem II: where and how does the degradation of the D1 protein by FtsH proteases start under light stress?--Facts and hypotheses

Degradation of the reaction center-binding D1 protein of Photosystem II is central in photoinhibition of Photosystem II. In higher plant chloroplasts, Photosystem II complexes are abundant in the grana. It has been suggested that the Photosystem II complexes containing photodamaged D1 protein migrate for their repair from the grana to the non-appressed stroma thylakoids, where the photodamaged D1 protein is degraded by a specific protease(s) such as filamentation temperature sensitive H (FtsH) protease. There are several possible ways to activate the FtsH proteases. As FtsH is a membrane-bound ATP-dependent metalloprotease, it requires ATP and zinc as essential part of its catalytic mechanism. It is also suggested that a membrane protein(s) associated with FtsH is required for modulation of the FtsH activity. Here, we propose several possible mechanisms for activation of the proteases, which depend on oligomerization of the monomer subunits. In relation to the oligomerization of FtsH subunits, we also suggest unique distribution of active FtsH hexamers on the thylakoids: hexamers of the FtsH proteases are localized near the Photosystem II complexes at the grana. Degradation of the D1 protein probably takes place in the grana rather than in the stroma thylakoids to circumvent long-distance migration of both the Photosystem II complexes containing the photodamaged D1 protein and the proteases.     

AFM imaging of immobilized fibronectin: does the surface conjugation scheme affect the protein orientation/conformation?

Covalent grafting of biomolecules could potentially improve the biocompatibility of materials. However, these molecules have to be grafted in an active conformation to play their biological roles. The present work aims at verifying if the surface conjugation scheme of fibronectin (FN) affects the protein orientation/conformation and activity. FN was grafted onto plasma-treated fused silica using two different crosslinkers, glutaric anhydride (GA) or sulfosuccinimidyl 4-(p-maleimidophenyl)butyrate (SMPB). Fused silica was chosen as a model surface material because it presents a roughness well below the dimensions of FN, therefore allowing AFM analyses with appropriate depth resolution. Cell adhesion assays were performed to evaluate the bioactivity of grafted FN. Cell adhesion was found to be higher on GA-FN than on SMPB-FN. Since FN-radiolabeling assays allowed us to rule out a surface concentration effect (approximately 80 ng/cm2 of FN on both crosslinkers), it was hypothesized that FN adopted a more active conformation when grafted via GA. In this context, the FN conformation on both crosslinkers was investigated through AFM and contact angle analyses. Before FN grafting, GA- and SMPB-modified surfaces had a similar water contact angle, topography, and roughness. However, water contact angles of GA-FN and SMPB-FN surfaces clearly show differences in surface hydrophilicity, therefore indicating a dependence of protein organization toward the conjugation strategy. Furthermore, AFM results demonstrated that surface topography and roughness of both FN-conjugated surfaces were significantly different. Distribution analysis of FN height and diameter confirmed this observation as the protein dimensions were significantly larger on GA than SMPB. This study confirmed that the covalent immobilization scheme of biomolecules influences their conformation and, hence, their activity. Consequently, selecting the appropriate conjugation strategy is of paramount importance in retaining molecule bioactivity.     

Does substrate oxidation determine the regioselectivity of cyclohexene and propene oxidation by cytochrome p450?

DFT and QM/MM computations of allylic C-H hydroxylation versus C=C epoxidation in cyclohexene and propene by Compound I of P450cam demonstrate that the relative barriers for the oxidative processes themselves are not good predictors of the observed selectivity. However, a kinetic expression previously developed (Kozuch, S.; Shaik, S. J. Am. Chem. Soc. 2006, 128, 3355) for catalytic cycles under steady-state conditions, predicts, in accord with experiment, that propene will undergo exclusive C=C epoxidation, while cyclohexene will undergo both reactions with a small preference for epoxidation. The model expression for the effective barrier of the cycle forms a general basis for understanding and predicting the selectivity of P450 isozymes.     

The Role of Cytochrome P450 AbyV in the Final Stages of Abyssomicin C Biosynthesis

Abyssomicin C and its atropisomer are potent inhibitors of bacterial folate metabolism. They possess complex polycyclic structures, and their biosynthesis has been shown to involve several unusual enzymatic transformations. Using a combination of synthesis and in vitro assays we reveal that AbyV, a cytochrome P450 enzyme from the aby gene cluster, catalyses a key late-stage epoxidation required for the installation of the characteristic ether-bridged core of abyssomicin C. The X-ray crystal structure of AbyV has been determined, which in combination with molecular dynamics simulations provides a structural framework for our functional data. This work demonstrates the power of combining selective carbon-13 labelling with NMR spectroscopy as a sensitive tool to interrogate enzyme-catalysed reactions in vitro with no need for purification.     

The mechanism of the α-ketoacid-hydroxylamine amide-forming ligation

Three-ring circus! Surprisingly complex molecular acrobatics are observed in the mechanism of the α-ketoacid-hydroxylamine amide-forming ligation reaction. Although this remarkable reaction can already be used for the chemoselective union of large, unprotected peptide fragments the elucidated mechanism provides important clues to extending its application to larger and more complex biological targets.     

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.     

How does mechanism of biosorption determine the differences between the initial and equilibrium adsorption states?

The analysis of various models assumed to represent the influence of pH on heavy metals biosorption equilibrium is presented. It shows that all of them lead to the same mathematical expressions (e.g. the Langmuir or the Flory adsorption isotherm equations) when the pH effects are neglected. Even if considering the pH effects, some of them (competitive adsorption and ion-exchange models, for instance) still lead to analogical expressions for sorption isotherm equations. The accepted mechanism of biosorption may, however, influence strongly the differences between the initial and equilibrium states of biosorption system.     

How does the hydrophobic content of methacrylate ABA triblock copolymers affect polymersome formation?

Polymersomes are exciting self-assembled structures with great potential in pharmaceutical applications. A systematic investigation of a novel series of methacrylate-based polymersomes is reported in this study. Five well-defined ABA triblock copolymers with A being based on tri(ethylene glycol) methyl methacrylate and B being based on 2-(diethylamino)ethyl methacrylate (DMAEMA) were synthesized using a living polymerization method. The effect of the composition of the ABA triblock copolymers on the thickness of the hydrophobic membrane of the polymersomes and the release of a model drug is demonstrated.